畜牧兽医学报  2021, Vol. 52 Issue (3): 662-675. DOI: 10.11843/j.issn.0366-6964.2021.03.010    PDF    
基于全基因组重测序技术的狮头鹅Indel标记分析
王慧芳1, 周光现1, 孙永峰2, 陈新企1, 雷栩斌1, 张锐毅1, 贾汝敏1, 赵志辉1     
1. 广东海洋大学滨海农业学院, 湛江 524088;
2. 吉林农业大学动物科学技术学院, 长春 130118
摘要:旨在利用重测序技术分析狮头鹅基因Indel差异并探讨其产蛋调控通路。本研究选用产蛋量有差异的成年狮头鹅和吉林白鹅母鹅各5只,采用全基因组重测序方法,对两个品种基因组水平的插入缺失突变情况进行比较分析。结果表明,测序共得到350.35 Gb数据,其中Q20与Q30平均值分别为96.74%、92.19%;通过注释分析发现了748 821.2个Indel变异位点。10只试验鹅的Indel位点在编码区(CDS)的分布趋势和数量与参考基因组相似,其中有11 737个Indel位点导致了插入或删除;在此基础上,进一步分析筛选产蛋相关通路,发现在GnRH信号通路、雌激素信号通路、催产素信号通路、催乳素信号通路以及孕酮介导的卵母细胞成熟通路中分别有15、1、1、1、18个变异基因参与狮头鹅产蛋调控。结果显示,通过Indel位点分析表明,GnRH信号通路、雌激素信号通路、催产素信号通路、催乳素信号通路和孕酮介导的卵母细胞成熟通路均与狮头鹅产蛋有关,并进一步筛选出可能的产蛋相关基因。这些数据将为我们研究狮头鹅繁殖性能提供有效基础和依据。
关键词狮头鹅    重测序    Indel    信号通路    产蛋    
Analysis of Indel Markers of Shitou Goose Based on Whole Genome Resequencing Technology
WANG Huifang1, ZHOU Guangxian1, SUN Yongfeng2, CHEN Xinqi1, LEI Xubin1, ZHANG Ruiyi1, JIA Rumin1, ZHAO Zhihui1     
1. College of Coastal Agricultural Sciences, Guangdong Ocean University, Zhanjiang 524088, China;
2. College of Animal Science and Technology, Jilin Agricultural University, Changchun 130118, China
Abstract: This study aimed to use resequencing technology to analyze the differences in Shitou goose gene Indel and to explore its egg production regulatory pathway. Five adult Shitou geese and five Jilin white geese with different egg production were selected, and the whole genome resequencing method was used to compare and analyze the indel mutations at the genome level of the two breeds. The results showed that 350.35 Gb data was obtained by sequencing, among which the average values of Q20 and Q30 reached 96.74% and 92.19%, respectively. Annotation analysis found 748 821.2 Indel mutation sites. The distribution trend and the number of Indel loci in the coding region (CDS) of 10 experimental geese were similar to those of the reference genome. Among them, 11 737 Indel loci resulted in insertion or deletion. On the basis, egg-laying related pathways were further analyzed and screened, and 15, 1, 1, 1, and 18 variant genes in the GnRH, estrogen, oxytocin, prolactin signaling pathways and progesterone-mediated oocyte maturation pathway, respectively involved in the regulation of goose egg-laying were found. The results suggest that the GnRH, estrogen, oxytocin, prolactin signaling pathways and progesterone-mediated oocyte maturation pathway are all related to Shitou goose egg-laying, and egg-laying related genes were further screened. These data will provide an effective basis for studying the reproductive performance of Shitou geese.
Key words: Shitou goose    resequencing    Indel    signaling pathway    egg-laying    

中国是世界上家鹅品种最多、驯养历史最悠久的国家之一。由于我国幅员辽阔,各地自然生态条件复杂多样,不同时期的经济文化背景不同,对鹅的选择利用目的不同,逐步形成了具有不同遗传特性和生产性能的地方品种,分布遍及全国[1-4]。目前,随着鹅相关产品需求水平的增长,低繁殖力成为了制约鹅产业发展的主要因素[5-6]。影响鹅产蛋性能的因素是多面的,如遗传基因、生理状况、环境条件、饲料营养和使用年限等,各种应激(如停电停水、噪音等)都可显著降低鹅的产蛋量[7-10]。狮头鹅起源于鸿雁,是我国第一批受保护的畜禽品种资源之一,因长期受人类和自然环境的影响,有着高度的遗传多样性,其具有耐粗饲、生长速度快、抗逆性强、耗料少、肉质好等优良特性,被誉为“世界鹅王”[11-12]。狮头鹅成年种公鹅体重为10~12 kg,母鹅为9~10 kg,母鹅年产蛋25~32枚[13-14]。吉林白鹅是由吉林农业大学科技工作者历经20余年精心培育而成,其成年公鹅体重4.1 kg,母鹅3.2 kg,年产蛋数达100枚以上,具有产蛋多、耐粗饲、抗病力强等特点[15]。狮头鹅与吉林白鹅产蛋量相差悬殊,适合作为探索鹅繁殖性状相关研究的素材。

随着基因组学技术的发展,使得大规模、高通量检测狮头鹅基因组内的变异位点,筛选产蛋相关通路以及相关产蛋基因成为可能[16-18]。Indel (insertion deletion) 是影响基因组中生产繁殖性能的关键遗传标记,研究发现,生物学基因组中Indel起着关键性的作用[18-20]。Indel变异一般比SNP变异少,但同样反映样品与参考基因组之间的差异,并且编码区的Indel会引起移码突变,导致基因功能上的变化[21-22]。本研究采用第三代重测序技术,获得了狮头鹅插入/缺失(Indel)位点,经过数据的比较分析发现,大量基因与狮头鹅的产蛋性能相关。经过数据分析发现,GnRH信号通路、雌激素信号通路、催产素信号通路、催乳素信号通路以及孕酮介导的卵母细胞成熟通路与狮头鹅产蛋性状相关。利用比较基因组学筛选鉴定两种鹅群体中特异性Indel位点[23],为狮头鹅的群体遗传分析提供基础,并进一步筛选产蛋相关基因,为狮头鹅生产研究提供新的思路。

1 材料与方法 1.1 试验材料

本研究选取健康成年狮头鹅母鹅(年均产蛋25~32枚)和吉林白鹅母鹅(年均产蛋≥100枚)各5只(R1~R5:狮头鹅;R6~R10:吉林白鹅),静脉翅下采血3 mL,加入抗凝剂,置于-20 ℃保存备用。

1.2 血液DNA提取

血液样品DNA使用美基生物公司的HiPure Blood DNA Mini Kit试剂盒进行提取。经1%凝胶电泳检测DNA完整性,并使用Nanodrop检测其纯度。合格样品DNA送北京百迈客公司测序。

1.3 全基因组重测序

1.3.1 文库构建与测序   样品基因组DNA检测合格后,用机械打断的方法(超声波)将DNA片段化,然后进行片段纯化,末端修复,3′端加A以及连接测序接头。通过琼脂糖凝胶电泳选择片段的大小,PCR扩增目的片段形成测序文库,对质检合格的文库进行Illumina测序。

1.3.2 测序数据质控   测序数据通过公式Qphred =-10lg(e)(Qphred为碱基质量值;e为测序错误率)转化可以得到碱基测序的错误率。在碱基识别过程中计算Phred数值(一种预测碱基判别发生错误的概率模型),其对应关系如表 1所示。

表 1 预测碱基判别发生错误概率结果 Table 1 Prediction incorrect rate result on base recognition

1.3.3 低质量数据过滤   原始测序序列(Raw Reads)含有带接头的以及低质量的Reads,需对Raw Reads进行过滤得到Clean Reads后,方可用于后续的信息分析。数据过滤的主要步骤:1) 去除带接头(adapter)的reads;2) 过滤N含量超过10%的reads;3) 去除质量值低于10的碱基超过50%的reads。

1.4 Small Indel的检测与注释

将在参考基因组上定位后的Clean Reads结果进行分析,检测是否存在Small Indel变异。使用GATK软件工具包检测Small Indel插入与缺失;使用Picard(http://sourceforge.net/projects/picard/)过滤冗余reads(MarkDuplicates),以保证检测结果的准确性;然后使用GATK的HaplotypeCaller(局部单体型组装)算法进行Indel的变异检测;通过SnpEff软件实现Small Indel的注释。

1.5 DNA水平的变异基因挖掘

比较筛选狮头鹅与吉林白鹅两个群体特有的差异Indel位点,通过BLAST变异基因与GO、COG、KEGG等功能数据库比对,筛选产蛋相关通路。通过KEGG注释筛选产蛋相关候选基因,用以分析基因功能。

2 结果 2.1 数据质量控制

由电泳检测结果(图 1)可以看出,DNA纯度以及完整性良好。另外检测得出OD260 nm/OD280 nm值均在1.88左右,其浓度、体积和质量平均值分别为:173.4 ng·μL-1、99.1 μL和18.7 μg。

R1~R5. 狮头鹅;R6~R10. 吉林白鹅。下同 R1-R5. Shitou goose; R6-R10. Jilin white goose. The same as below 图 1 琼脂糖凝胶电泳检测结果 Fig. 1 The detection results of agarose gel electrophoresis

通过对原始数据各碱基比例分布和碱基质量分布的分析来进行质量控制,得到图 2。从图中可以看出,AT曲线以及GC曲线均重合,证明AT碱基以及GC碱基比例均相等,无明显偏向性,曲线较平缓,未发生分离,测序结果正常,该数据可用于下一步分析。

N. 测序中识别不出的碱基 N. The unrecognized bases in sequencing 图 2 样品各碱基比例分布 Fig. 2 The proportion of each base of the samples

通过对狮头鹅和吉林白鹅全基因组重测序,原始数据过滤处理后共得到350.35 Gb的Clean Data,各样本的Raw data在71 553 908~208 410 322 Mb之间,Q20与Q30平均值分别达到96.74%、92.19%。GC含量在43.67%~45.64%之间(表 2),结果表明,测序质量良好,可进行后续分析。

表 2 测序数据质量表 Table 2 Summary of sequence data quality
2.2 与参考基因组的比对分析

参考基因组为Anser cygnoides (Swan goose),其版本为AnsCyg_PRJNA183603_v1.0;下载网址:https://www.ncbi.nlm.nih.gov/genome/?term=Anser+cygnoides,参考基因组的大小为1 134 Mb(表 3)。由表 4可知,所有样本数据与参考基因组的比对率在93.18%~94.96%之间,平均比对率为94.15%,平均覆盖深度为26×,1×覆盖度(至少有1个碱基的覆盖)在99.47%以上,5×覆盖度(至少有5个碱基的覆盖)在96.84%以上。结果表明,测序物种同参考基因组相似度高,证明测序准确,可进行后续的相关分析。

表 3 参考基因组信息 Table 3 Reference genome information
表 4 测序深度及覆盖度统计 Table 4 Sequencing depth and coverage statistics
2.3 Indel序列长度及分布分析

利用GATK从10只试验鹅与参考基因组比对结果中检测Indel变异(小于50 bp碱基片段的插入与缺失)。然后对Indel序列的分布位置进行注释及分析(表 5)。

表 5 Indel位点的检测及注释 Table 5 Detection and annotation of Indel loci

通过SnpEff软件注释分析Indel位点,分别从10只试验鹅中发现了786 695、782 331、766 853、775 202、718 792、717 274、708 825、738 934、753 315、739 991个Indel位点(平均有748 821.2个),其数量和分布与参考基因组相似。在这10只试验鹅中,Indel位点主要分布在内含子,吉林白鹅与狮头鹅Indel位点在内含子中分别占有54.18%和52.12%;其次是基因间区,两种鹅分别占有36.05%与33.91%(图 3)。大部分的Indel位点分布在内含子与基因间区中,CDS区中分布的Indel位点最少。

图 3 吉林白鹅和狮头鹅相同Indels在基因组上的分布 Fig. 3 Genomic distribution of the same Indels of Jilin white goose and Shitou goose

在这10只试验鹅全基因组中,Indel位点的长度分布趋势也基本一致。对群体Indel突变位点的碱基长度进行统计并分析,得到图 4,发现全基因组中长度为1 bp的Indel位点数量超过310 000个(约占全基因组所有Indel的43%);长度为2~4 bp的Indel位点数量为2 529 498个;长度为5~9 bp的Indel位点数量为860 657个;长度大于等于10 bp的Indel位点数量共808 872个。并且在编码区中,长度为1 bp的Indel位点数量共22 754个(约占编码区所有Indel的33%);长度为2~4 bp的Indel位点数量为22 426个;长度为5~9 bp的Indel位点数量为10 080个;长度大于等于10 bp的Indel位点数量共13 491个。另外,图 4说明了在片段变异10 bp长度的范围内,随着碱基片段长度增加,发生Indel变异的数量随之减少,且数量波动较大。

图 4 全基因组和编码区的Indel长度分布情况 Fig. 4 Distribution map of Indel length in the population's entire genome and coding region

10只试验鹅的Indel位点在CDS区的分布趋势和数量也与参考基因组相似,分别有7 090、7 026、7 022、7 060、6 696、6 649、6 651、6 862、6 868和6 827个Indel位点(表 6)。其中共有11 737个Indel位点导致了插入或删除。

表 6 Indel变异在全基因组和编码区的统计情况 Table 6 Statistics of Indel mutations in the whole genome and coding regions
2.4 变异基因注释分析

通过GO注释分析富集了到细胞组分(cellular component)17个、分子功能(molecular function)20个和生物过程(biological process)23个,详见表 7

表 7 变异基因功能富集分析的相关注释 Table 7 Annotation of function enrich analysis of variant genes

通过KEGG注释分析富集到253条通路,其中筛选出与产蛋相关的通路,分别为GnRH信号通路、雌激素信号通路、催乳素信号通路、催产素信号通路以及孕酮介导的卵母细胞成熟通路。通过比对测序得到的群体基因注释列表发现了36个产蛋相关变异基因,详见表 8。这为狮头鹅后期高产基因的筛选奠定了基础。

表 8 产蛋相关通路及变异基因 Table 8 Egg-laying related pathways and variant genes
3 讨论

狮头鹅是我国乃至世界著名的大型水禽品种,也是我国南方养殖范围较广的经济家禽,其产蛋期一般为9月至次年4月[24-25]。吉林白鹅是我国北方养殖较广的经济家禽,其产蛋旺季在2~6月,研究表明,两种鹅日照长短不一样[26]。甄霆等[27]研究发现,在浙东白鹅的繁殖季节通过调节光照可影响FSH和INH的表达,从而延长其产蛋时间,缩短抱窝和就巢时间,增加鹅的产蛋量,说明鹅的繁殖性能受光照的影响,合理地控制光照可延长鹅的产蛋期,缩短休产期,提高鹅的产蛋性能。姚颖等[28-31]研究发现,浙东白鹅约83.04%的产蛋是在有光照的情况下,仅有16.96%的产蛋是在无光照的情况下,说明产蛋行为与光照密切相关[28]。产蛋性状的遗传力较低,其由微效多基因控制,三代测序技术的发展为研究狮头鹅产蛋性状提供了新的思路。由于狮头鹅较强的就巢性和产蛋季节性导致其繁殖力较低,从而造成狮头鹅产蛋性能低下。提高狮头鹅的产蛋性能一直是研究的主要目标之一。

本研究利用重测序技术对5只狮头鹅和5只吉林白鹅全基因组测序后共得到350.35 Gb数据,其中Q20与Q30平均值分别达到96.74%、92.19%,测序数据结果可靠。在此基础上,通过注释分析得到Indel位点平均有748 821.2个变异。大部分(53%左右)的Indel位点分布在内含子中,蛋白质编码区中分布的Indel位点最少(11 737个)。研究还发现,Indel位点的长度分布趋势也基本一致,10只试验鹅的Indel位点在CDS区的分布趋势和数量也与参考基因组相似。在此基础上,通过plink软件对SNPs进行主成分分析(PCA)发现(数据未显示),吉林白鹅与狮头鹅两个群体能够有效地分离。从SNP层次观测,两种鹅具有较大的差异性,这种差异性包含潜在主导其不同性状的SNP位点。吉林白鹅体型小,年产蛋量高,而狮头鹅体型大,年产蛋量低。张泽宾[32]对22只野鸭和56只家鸭全基因组测序得到310万个Indel变异位点,对Indel注释发现,绝大多数的变异位于基因间区及内含子等非编码区域,此研究发现,编码区中的Indel突变比率与其它区间相比较低。刘杰[33]分别选取北京油鸡和科宝肉鸡中高采食量个体48只和低采食量个体48只进行全基因组重测序,得到的Indel变异位点在内含子和基因间区最多,编码区较少。章双杰等[34]对4种鸭群进行全基因组重测序分析,通过检测到的Indel位点发现,其Indel变异大部分存在于基因间或是基因内的内含子区,证明编码区的变异相对较少。以上研究结果均与本试验Indel变异分布结果一致,非编码区的Indel变异多于编码区,这一结果符合生物学特性。对差异分析筛选的候选Indel标记的作用可进一步研究,为鹅生长性状遗传机理的相关研究奠定理论基础。

在此基础上,通过数据库比对分析进一步筛选产蛋相关通路,发现在促性腺激素释放激素(GnRH) 信号通路、催乳素(PRL)信号通路、催产素(OT)信号通路和孕酮介导的卵母细胞成熟通路均有Indel位点变异。GnRH作为HPG轴的关键信号分子与繁殖周期的变化呈相关性。产蛋性状是由多个相关基因共同表达和相互作用来调控的,并受多种外源因素影响[35-37]。常见的产蛋相关主要基因有促性腺激素释放激素(GnRH)[38]、促性腺激素抑制激素(GnIH)[39]、催乳素(PRL)[40]、促卵泡激素(FSH)[41]、雌激素受体(ESR)[42]等。张克山等[43]研究表明,在产蛋前期或高峰期,四川白鹅的GnRH mRNA表达量高于休产期,并指出在调控鹅产蛋性状中GnRH信号通路发挥重要作用。催产素通路在鹅的繁殖、就巢和产蛋等方面有重要作用。本课题组成功克隆了狮头鹅PRL基因并进行了原核表达,研究发现,狮头鹅与其他禽类PRL具有高度保守性[44]。催乳素通路是细胞生长和分化过程中调控鹅性成熟的主要通路,同时也是鹅产蛋性状的候选遗传标记基因之一[45]PRL基因的多态性主要与产卵性状相关。催乳素通路在鹅卵巢发育和产卵过程中也发挥重要作用。孕酮介导的卵母细胞成熟通路对卵泡发育起调节作用。卵泡和排卵取决于促性腺激素的释放。在家禽生产中,孕酮的分泌与卵巢的生长和状态紧密相关。在排卵过程中,孕酮作为正反馈信号在排卵前触发GnRH/LH通路[46]。总之,本研究筛选出的产蛋相关变异通路和调节途径,以及进一步筛选出的与产蛋相关的候选基因,为研究狮头鹅繁殖性能提供有效基础和依据。

4 结论

通过全基因组测序后注释分析Indel变异位点,发现大部分的Indel位点分布在内含子中,蛋白质编码区中分布的Indel位点最少;另外发现,Indel位点长度分布趋势基本一致,在CDS区的分布趋势和数量与参考基因组相似,非编码区的Indel变异数量多于编码区,这一结果符合生物学特性。比对测序得到的群体基因注释列表筛选出了产蛋相关变异基因,本研究突破了传统育种的局限性,有效地为低遗传力产蛋性状的遗传进展加快了步伐,提供研究鹅繁殖性能重要的分子标记理论数据,可以在实际选育工作中提高鹅养殖产业的经济效益。

参考文献
[1] PROSSER D J, WU J X, ELLIS E C, et al. Modelling the distribution of chickens, ducks, and geese in China[J]. Agric Ecosyst Environ, 2011, 141(3-4): 381–389. DOI: 10.1016/j.agee.2011.04.002
[2] ZHENG R B, SMITH L M, PROSSER D J, et al. Investigating home range, movement pattern, and habitat selection of bar-headed geese during breeding season at Qinghai Lake, China[J]. Animals (Basel), 2018, 8(10): 182.
[3] ZHANG P Y, ZOU Y A, XIE Y H, et al. Shifts in distribution of herbivorous geese relative to hydro-logical variation in East Dongting Lake wetland, China[J]. Sci Total Environ, 2018, 636: 30–38. DOI: 10.1016/j.scitotenv.2018.04.247
[4] 刘思扬, 刘秋翔. 饶平狮头鹅[J]. 广东畜牧兽医科技, 2019, 44(6): 10–13.
LIU S Y, LIU Q X. Lion-head goose from Raoping[J]. Guangdong Journal of Animal and Veterinary Science, 2019, 44(6): 10–13. (in Chinese)
[5] YANG Y Z, YAO Y, CAO Z F, et al. Histological characteristics of follicles and reproductive hormone secretion during ovarian follicle development in laying geese[J]. Poult Sci, 2019, 98(11): 6063–6070. DOI: 10.3382/ps/pez278
[6] YAO Y, YANG Y Z, GU T T, et al. Comparison of the broody behavior characteristics of different breeds of geese[J]. Poult Sci, 2019, 98(11): 5226–5233. DOI: 10.3382/ps/pez366
[7] LIU G J, CHEN Z F, ZHAO X H, et al. Meta-analysis: Supplementary artificial light and goose reproduction[J]. Anim Reprod Sci, 2020, 214: 106278. DOI: 10.1016/j.anireprosci.2020.106278
[8] MA W Q, CHENG H Z, ZHAO D H, et al. Effects of dietary Enteromorpha powder supplementation on productive performance, egg quality, and antioxidant performance during the late laying period in Zi geese[J]. Poult Sci, 2020, 99(2): 1062–1068. DOI: 10.1016/j.psj.2019.10.003
[9] LAMERIS T K, DE JONG M E, BOOM M P, et al. Climate warming may affect the optimal timing of reproduction for migratory geese differently in the low and high Arctic[J]. Oecologia, 2019, 191(4): 1003–1014. DOI: 10.1007/s00442-019-04533-7
[10] 罗艺, 刘作兰, 王启贵, 等. 饲粮能量水平对初产期四川白鹅生产性能、蛋品质和孵化性能的影响[J]. 畜牧兽医学报, 2018, 49(2): 439–448.
LUO Y, LIU Z L, WANG Q G, et al. Effects of dietary energy concentrations on growth perfor-mance, egg quality and hatching performance in primiparous laying Sichuan white geese[J]. Acta Veterinaria et Zootechnica Sinica, 2018, 49(2): 439–448. (in Chinese)
[11] 林潮佳. 世界鹅王——狮头鹅[J]. 生物学通报, 1986(11): 42.
LIN C J. King of world goose-Lion head goose[J]. Bulletin of Biology, 1986(11): 42. (in Chinese)
[12] 包海缘, 欧阳宏佳, 唐军, 等. 不同生长期狮头鹅体尺性状、屠宰性能和种蛋品质的测定与分析[J]. 养禽与禽病防治, 2019(6): 36–39.
BAO H Y, OUYANG H J, TANG J, et al. Determination and analysis of traits, slaughter performance and egg quality of lion head goose in different growth stages[J]. Poultry Husbandry and Disease Control, 2019(6): 36–39. (in Chinese)
[13] 林树欣, 潘育璇, 黄靓, 等. 狮头鹅种鹅的饲养管理[J]. 广东畜牧兽医科技, 2019, 44(5): 26–27, 35.
LIN S X, PAN Y X, HUANG L, et al. Feeding and management of Shitou Goose for breeding[J]. Guangdong Journal of Animal and Veterinary Science, 2019, 44(5): 26–27, 35. (in Chinese)
[14] 孙永泰. 狮头鹅的饲养技术[J]. 水禽世界, 2016(2): 59.
SUN Y T. Breeding technology of lion head goose[J]. Waterfowl World, 2016(2): 59. (in Chinese)
[15] 王雪丽, 孙姝芬. 吉林白鹅简介[J]. 农村科学实验, 2006(12): 28.
WANG X L, SUN S F. Introduction of Jilin white goose[J]. Rural Science Experiment, 2006(12): 28. (in Chinese)
[16] GAO G, ZHANG K, ZHAO X, et al. Molecular cloning of the goose GnRH gene and identification of GnRH polymorphisms associated with laying traits[J]. Br Poult Sci, 2020, 61(5): 502–507. DOI: 10.1080/00071668.2020.1758298
[17] WU Y P, ZHAO X Y, CHEN L, et al. Transcriptomic analyses of the hypothalamic-pituitary-gonadal axis identify candidate genes related to egg production in Xinjiang Yili geese[J]. Animals (Basel), 2020, 10(1): 90.
[18] ZHAO Q Q, CHEN J P, ZHANG X H, et al. Genome-wide association analysis reveals key genes responsible for egg production of lion head goose[J]. Front Genet, 2020, 10(1): 1391.
[19] LI X P, LU Y L, LIU X F, et al. Identification of chicken FSHR gene promoter and the correlations between polymorphisms and egg production in Chinese native hens[J]. Reprod Domest Anim, 2019, 54(4): 702–711. DOI: 10.1111/rda.13412
[20] CABANILLAS R, DIÑEIRO M, CIFUENTES G A, et al. Comprehensive genomic diagnosis of non-syndromic and syndromic hereditary hearing loss in Spanish patients[J]. BMC Med Genomics, 2018, 11(1): 58. DOI: 10.1186/s12920-018-0375-5
[21] LI X J, YE J W, HAN X L, et al. Whole-genome sequencing identifies potential candidate genes for reproductive traits in pigs[J]. Genomics, 2020, 112(1): 199–206. DOI: 10.1016/j.ygeno.2019.01.014
[22] 唐莹, 王金玉, 张跟喜, 等. Myf5基因8个SNPs位点与京海黄鸡生长和繁殖性状的关联分析[J]. 畜牧兽医学报, 2014, 45(6): 863–870.
TANG Y, WANG J Y, ZHANG G X, et al. Eight SNPs of the Myf5 gene and its relationship with growth and reproductive traits in Jinghai yellow chicken[J]. Acta Veterinaria et Zootechnica Sinica, 2014, 45(6): 863–870. (in Chinese)
[23] HARDISON R C. Comparative genomics[J]. PLoS Biol, 2003, 1(2): E58. DOI: 10.1371/journal.pbio.0000058
[24] 吕敏娜, 陈伟峰, 廖申权. 狮头鹅之健康养殖[J]. 畜牧兽医科技信息, 2019(9): 143–145.
LV M N, CHEN W F, LIAO S Q. Healthy breeding of lion head goose[J]. Chinese Journal of Animal Husbandry and Veterinary Medicine, 2019(9): 143–145. (in Chinese)
[25] TANG J, SHEN X, OUYANG H J, et al. Transcriptome analysis of pituitary gland revealed candidate genes and gene networks regulating the growth and development in goose[J/OL]. Anim Biotechnol, 2020, doi: 10.1080/10495398.2020.1801457.
[26] 隋玉健, 杨薇, 许云鹏, 等. 复合中草药添加剂对吉林白鹅生产性能的影响[J]. 黑龙江畜牧兽医, 2020(20): 106–108, 113.
SUI Y J, YANG W, XU Y P, et al. Effect of compound Chinese herbal medicine additive on the production performance of Jilin white goose[J]. Heilongjiang Animal Science and Veterinary Medicine, 2020(20): 106–108, 113. (in Chinese)
[27] 甄霆, 张响英, 张海波, 等. 光照对浙东白鹅产蛋和相关基本表达量的影响[J]. 黑龙江畜牧兽医, 2018(8): 67–68.
ZHEN T, ZHANG X Y, ZHANG H B, et al. Effects of light treatments on egg production and related gene expression level in Eastern Zhejiang White Goose[J]. Heilongjiang Animal Science and Veterinary Medicine, 2018(8): 67–68. (in Chinese)
[28] 姚颖, 曹正锋, 杨耀宗, 等. 浙东白鹅产蛋与就巢行为学研究[J]. 中国畜牧杂志, 2018, 54(4): 113–116.
YAO Y, CAO Z F, YANG Y Z, et al. Study on egg laying and nesting behavior of eastern Zhejiang White Goose[J]. Chinese Journal of Animal Science, 2018, 54(4): 113–116. (in Chinese)
[29] OUYANG Q Y, HU S Q, WANG G S, et al. Comparative transcriptome analysis suggests key roles for 5-hydroxytryptamlne receptors in control of goose egg production[J]. Genes (Basel), 2020, 11(4): 455.
[30] LIU H H, WANG J W, LI L, et al. Transcriptome analysis revealed the possible regulatory pathways initiating female geese broodiness within the hypothalamic-pituitary-gonadal axis[J]. PLoS One, 2018, 13(2): e0191213.
[31] CAO Z, MENG B, FAN R, et al. Comparative proteomic analysis of ovaries from Huoyan geese between pre-laying and laying periods using an iTRAQ-based approach[J]. Poult Sci, 2018, 97(6): 2170–2182.
[32] 张泽宾. 利用全基因组重测序研究家鸭驯化的选择信号及驯化时间[D]. 北京: 中国农业大学, 2018.
ZHANG Z B. Whole-genome resequencing reveals signatures of selection and timing of duck domestication[D]. Beijing: China Agricultural University, 2018. (in Chinese)
[33] 刘杰. 利用高通量测序技术挖掘影响鸡剩余采食量的遗传变异和功能基因[D]. 北京: 中国农业科学院, 2017.
LIU J. Exploring genomic variants and genes related to residual feed intake in chickens by high-throughput sequencing[D]. Beijing: Chinese Academy of Agricultural Sciences, 2017. (in Chinese)
[34] 章双杰, 王靖, 杨建生, 等. 长江下游地区地方鸭种全基因组序列测定及遗传多样性分析[J]. 中国家禽, 2016, 38(23): 5–9.
ZHANG S J, WANG J, YANG J S, et al. Analysis on complete genome sequencing and genetic diversity of native duck breeds in the lower reaches of Yangtze River[J]. China Poultry, 2016, 38(23): 5–9. (in Chinese)
[35] RADWAN L M. Genetic improvement of egg laying traits in Fayoumi chickens bred under conditions of heat stress through selection and gene expression studies[J]. J Therm Biol, 2020, 89: 102546.
[36] WOLC A, JANKOWSKI T, ARANGO J, et al. Investigating the genetic determination of clutch traits in laying hens[J]. Poult Sci, 2019, 98(1): 39–45.
[37] ROSI N ′ SKI A, NOWACZEWSKI S, KONTECKA H, et al. Analysis of the laying rhythm and reproductive traits of geese[J]. Folia Biol (Krakow), 2006, 54(3-4): 145–152.
[38] DUBOIS E A, ZANDBERGEN M A, PEUTE J, et al. Evolutionary development of three gonadotropin-releasing hormone (GnRH) systems in vertebrates[J]. Brain Res Bull, 2002, 57(3-4): 413–418.
[39] 胡彦竞科, 张克山, 钟航, 等. 鹅GnIH基因多态性及其与产蛋量的关联性研究[J]. 中国畜牧杂志, 2016, 52(21): 17–19, 104.
HUYAN J K, ZHANG K S, ZHONG H, et al. Study on GnIH gene polymorphism of goose and its association with egg production[J]. Chinese Journal of Animal Science, 2016, 52(21): 17–19, 104. (in Chinese)
[40] GOFFIN V, BINART N, TOURAINE P, et al. Prolactin: the new biology of an old hormone[J]. Annu Rev Physiol, 2002, 64: 47–67.
[41] 郭军, 汤青萍, 章双杰, 等. 利用抑制消减杂交技术筛选鹅就巢行为相关基因[J]. 畜牧兽医学报, 2011, 42(10): 1477–1484.
GUO J, TANG Q P, ZHANG S J, et al. Identification of broodiness-related geese genes by suppression subtractive hybridization[J]. Acta Veterinaria et Zootechnica Sinica, 2011, 42(10): 1477–1484. (in Chinese)
[42] 张响英, 唐现文, 王利刚, 等. 实时荧光定量RT-PCR检测狮头鹅繁殖周期中PRL mRNA的表达[J]. 江苏农业科学, 2012, 40(9): 50–52.
ZHANG X Y, TANG X W, WANG L G, et al. Real-time fluorescence quantitative RT-PCR detection of PRL mRNA expression during the reproductive cycle of lion head goose[J]. Jiangsu Agricultural Sciences, 2012, 40(9): 50–52. (in Chinese)
[43] 张克山, 胡彦竞科, 韩笑哲, 等. 鹅不同繁殖时期GnRHGnIH基因表达和激素浓度变化分析[J]. 畜牧兽医学报, 2016, 47(8): 1720–1726.
ZHANG K S, HUYAN J K, HAN X Z, et al. Analysis of the serum concentrations and mRNA expression levels of GnRH and GnIH in geese during different reproductive periods[J]. Acta Veterinaria et Zootechnica Sinica, 2016, 47(8): 1720–1726. (in Chinese)
[44] 贾汝敏, 吴慧英, 刘铀, 等. 狮头鹅PRL基因的克隆及原核表达[J]. 畜牧兽医学报, 2011, 42(3): 343–348.
JIA R M, WU H Y, LIU Y, et al. Cloning and prokaryotic expression of PRL gene in Shitou Goose[J]. Acta Veterinaria et Zootechnica Sinica, 2011, 42(3): 343–348. (in Chinese)
[45] CHEN J, LIU H L, CAI Y F, et al. Mutations in the exon 10 of prolactin receptor gene change the egg production performance in Wanjiang white goose[J]. Mol Biol Rep, 2012, 39(1): 475–483. DOI: 10.1007/s11033-011-0761-y
[46] XIAO Y Q, NI Y D, HUANG Y B, et al. Effects of kisspeptin-10 on progesterone secretion in cultured chicken ovarian granulosa cells from preovulatory (F1-F3) follicles[J]. Peptides, 2011, 32(10): 2091–2097.