Genome-wide survey and analysis of microsatellites in giant panda (Ailuropoda melanoleuca), with a focus on the applications of a novel microsatellite marker system

Background

The giant panda (Ailuropoda melanoleuca) is a critically endangered species endemic to China. Microsatellites have been preferred as the most popular molecular markers and proven effective in estimating population size, paternity test, genetic diversity for the critically endangered species. The availability of the giant panda complete genome sequences provided the opportunity to carry out genome-wide scans for all types of microsatellites markers, which now opens the way for the analysis and development of microsatellites in giant panda.

Results

By screening the whole genome sequence of giant panda in silico mining, we identified microsatellites in the genome of giant panda and analyzed their frequency and distribution in different genomic regions. Based on our search criteria, a repertoire of 855,058 SSRs was detected, with mono-nucleotides being the most abundant. SSRs were found in all genomic regions and were more abundant in non-coding regions than coding regions. A total of 160 primer pairs were designed to screen for polymorphic microsatellites using the selected tetranucleotide microsatellite sequences. The 51 novel polymorphic tetranucleotide microsatellite loci were discovered based on genotyping blood DNA from 22 captive giant pandas in this study. Finally, a total of 15 markers, which showed good polymorphism, stability, and repetition in faecal samples, were used to establish the novel microsatellite marker system for giant panda. Meanwhile, a genotyping database for Chengdu captive giant pandas (n = 57) were set up using this standardized system. What’s more, a universal individual identification method was established and the genetic diversity were analysed in this study as the applications of this marker system.

Conclusion

The microsatellite abundance and diversity were characterized in giant panda genomes. A total of 154,677 tetranucleotide microsatellites were identified and 15 of them were discovered as the polymorphic and stable loci. The individual identification method and the genetic diversity analysis method in this study provided adequate material for the future study of giant panda.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1268-z) contains supplementary material, which is available to authorized users.     

卧龙圈养大熊猫遗传多样性现状及预测,

以中国最大的大熊猫圈养种群—四川卧龙中国大熊猫保护中心的圈养种群为对象,以８个大熊猫微卫星位点为分子标记, 探讨了大熊猫圈养种群的遗传多样性, 并与邛崃野生种群及其他７个濒危物种进行比较。微卫星数据表明, 圈养种群的遗传多样性水平(A＝5.5, He ＝0.620, Ho＝0.574) 低于邛崃野生种群(A＝9.8,He＝0.779,Ho＝0.581),但高于其他7 个濒危物种的种群(He＝0.13~0.46)。在此数据的基础上对未来100个世代内圈养种群遗传多样性的变化情况做出了预测。结果表明假设种群数量比现在扩大一倍, 经历100个世代后也只会使平均等位基因数少减少0.4。因此继续增加野生个体对保持遗传多样性的意义已经不大, 建议该圈养种群的保护策略应将重点放到制定更有效的繁殖计划以避免近交上。     

圈养大熊猫群体间的基因流状况分析

圈养群体的遗传管理一个非常重要的手段就是实现不同群体间的基因交流.为了全面评估大熊猫圈养群体间的基因流状况,本研究以卧龙中国大熊猫保护中心的31只圈养大熊猫(简称卧龙群体)和成都大熊猫繁育研究基地与楼观台陕西省珍稀野生动物抢救饲养研究中心的37只圈养大熊猫(简称成都群体,其中楼观台1只)为研究对象,以7个大熊猫微卫星位点为分子标记,发现卧龙群体和成都群体间的遗传分化水平很低(Fst=0.041,P=0.001);尽管整个圈养群体的近交程度较低(Fis=0.026),但是成都群体的近交系数(Fis=0.045)远高于卧龙群体的近交系数(Fis=0.002);谱系分析、贝叶斯分析和系统进化法分析均显示,这两个群体间存在着基因流,但是这种基因流是单向的.此结果提示各个大熊猫饲养单位之间必须实现更多的合作,将大熊猫群体作为一个管理单元进行管理,从而实现更多的基因流.     

Microsatellite variability reveals the necessity for genetic input from wild giant pandas (Ailuropoda melanoleuca) into the captive population

Recent success in breeding giant pandas in captivity has encouraged panda conservationists to believe that the ex situ population is ready to serve as a source for supporting the wild population. In this study, we used 11 microsatellite DNA markers to assess the amount and distribution of genetic variability present in the two largest captive populations (Chengdu Research Base of Giant Panda Breeding, Sichuan Province and the China Research and Conservation Center for the Giant Panda at Wolong, Sichuan Province). The data were compared with those samples from wild pandas living in two key giant panda nature reserves (Baoxing Nature Reserve and Wanglang Nature Reserve). The results show that the captive populations have retained lower levels of allelic diversity and heterozygosity compared to isolated wild populations. However, low inbreeding coefficients indicate that captive populations are under careful genetic management. Excessive heterozygosity suggests that the two captive populations have experienced a genetic bottleneck, presumably caused by founder effects. Moreover, evidence of increased genetic divergence demonstrates restricted breeding options within facilities. Based on these results, we conclude that the genetic diversity in the captive populations is not optimal. Introduction of genetic materials from wild pandas and improved exchange of genetic materials among institutions will be necessary for the captive pandas to be representative of the wild populations.     

Major histocompatibility complex class II variation in the giant panda (Ailuropoda melanoleuca)

Habitat destruction and human activity have greatly impacted the natural history of the giant panda (Ailuropoda melanoleuca). Although the genetic diversity of neutral markers has been examined in this endangered species, no previous work has examined adaptive molecular polymorphisms in the giant panda. Here, the major histocompatibility complex (MHC) class II DRB locus was investigated in the giant panda, using single-strand conformation polymorphism (SSCP) and sequence analysis. Comparisons of DNA samples extracted from faecal and blood samples from the same individual revealed that the two materials yielded similar quantities and qualities of DNA, as well as identical SSCP patterns and allelic sequences, demonstrating the reliability of DNA isolation from panda faeces. Analysis of faecal samples from 60 giant pandas revealed relatively low number of alleles: seven alleles. However, the alleles were quite divergent, varying from each other by a range of 7-47 nucleotide substitutions (4-25 amino acid substitutions). Construction of a neighbour-joining tree and comparisons among DRB alleles from other species revealed that both similar and highly divergent alleles survived in the bottlenecked panda populations. Despite species-specific primers used and excellent faecal DNA isolated, a lower level of heterozygosity than expected was still observed due to inbreeding. There were three types of evidence supporting the presence of balancing selection in the giant panda: (i) an obvious excess of nonsynonymous substitutions over synonymous at the antigen-binding positions; (ii) trans-species evolution of two alleles between the giant panda and other felids; and (iii) a more even distribution of alleles than expected from neutrality.     

Genetic variation in the major histocompatibility complex of the European brown hare (Lepus europaeus) across distinct phylogeographic areas

The major histocompatibility complex is one of the best studied systems in vertebrates providing evidence for the long-term action of selection. Here, we examined the intra- and inter-population genetic diversity of the MHC class II DRB locus in European brown hare (Lepus europaeus) and correlated the results with genetic variability already estimated from the MHC DQA locus and from maternally (mitochondrial DNA (mtDNA)) and biparentally (allozymes, microsatellites) inherited loci. L. europaeus showed remarkable genetic polymorphism in both DQA and DRB1 loci. The Anatolian populations exhibited the highest genetic polymorphism for both loci. Balancing selection has established increased variability in the European populations despite the founder effects after the last glaciation. Different evolutionary rates were traced for DRB1 and DQA loci, as evidenced by the higher number of common DRB1 than DQA alleles and the greater differences between DRB1 alleles with common origin in comparison with DQA alleles. The high number of rare alleles with low frequencies detected implies that frequency-dependent selection drives MHC evolution in the brown hare through the advantage of rare alleles. Both loci were under the influence of positive selection within the peptide-binding region. The functional polymorphism, recorded as amino acid substitutions within the binding pockets, fell also within distinct geographic patterns, yet it was much narrower than the genetic polymorphism. We hypothesize that certain structural and functional characteristics of the binding pockets set limitations to the actual shape of genetic polymorphism in MHC.     

Decreased microbial diversity and <Emphasis Type="Italic">Lactobacillus</Emphasis> group in the intestine of geriatric giant pandas (<Emphasis Type="Italic">Ailuropoda melanoleuca</Emphasis>)

It has been established beyond doubt that giant panda genome lacks lignin-degrading related enzyme, gastrointestinal microbes may play a vital role in digestion of highly fibrous bamboo diet. However, there is not much information available about the intestinal bacteria composition in captive giant pandas with different ages. In this study, we compared the intestinal bacterial community of 12 captive giant pandas from three different age groups (subadults, adults, and geriatrics) through PCR-denaturing gradient gel electrophoresis (DGGE) and real-time PCR analysis. Results indicated that microbial diversity in the intestine of adults was significantly higher than that of the geriatrics (p < 0.05), but not significant compared to the subadults (p > 0.05). The predominant bands in DGGE patterns shared by the twelve pandas were related to Firmicutes and Proteobacteria. Additionally, in comparison to healthy individuals, antibiotic-treated animals showed partial microbial dysbiosis. Real-time PCR analyses confirmed a significantly higher abundance of the Lactobacillus in the fecal microbiota of adults (p < 0.05), while other bacterial groups and species detected did not significantly differ among the three age groups (p > 0.05). This study revealed that captive giant pandas with different ages showed different intestinal bacteria composition.     

基于转录组测序的大熊猫多态性微卫星标记筛选

结合已公布的大熊猫Ailuropoda melanoleuca基因组和本实验室所测6只大熊猫的转录组数据,筛选多态性微卫星位点并分析其组成及特征。结果显示:共获得326个多态性微卫星位点,其中二碱基多态性微卫星最多,共228个,占69.93%;三、四、五、六碱基所占比例分别为9.51%、14.11%、5.21%、1.22%。根据分析结果中缺失率与标准差2项指标以及位点序列长度,选取20个多态性二碱基微卫星位点,用于25只大熊猫个体血液DNA进行PCR验证并做后续分析。结果表明:不同位点的等位基因数为2～8,平均等位基因数为3.70,观测杂合度、期望杂合度分别为0～1.000、0.280～0.784,平均值分别为0.472和0.532。在Bonferroni校正后,证实4个位点显著偏离哈迪-温伯格平衡,所有位点未观察到显著连锁不平衡(P>0.01)。20个位点多态信息含量(PIC)在0.246～0.734,其中具有高度多态性的位点9个(PIC>0.50),11个位点呈中度多态性(0.25相似文献   

Major histocompatibility complex alleles associated with parasite susceptibility in wild giant pandas

Major histocompatibility complex (MHC) polymorphism is thought to be driven by antagonistic coevolution between pathogens and hosts, mediated through either overdominance or frequency-dependent selection. However, investigations under natural conditions are still rare for endangered mammals which often exhibit depleted variation, and the mechanism of selection underlying the maintenance of characteristics remains a considerable debate. In this study, 87 wild giant pandas were used to investigate MHC variation associated with parasite load. With the knowledge of the MHC profile provided by the genomic data of the giant panda, seven DRB1, seven DQA1 and eight DQA2 alleles were identified at each single locus. Positive selection evidenced by a significantly higher number of non-synonymous substitutions per non-synonymous codon site relative to synonymous substitutions per synonymous codon site could only be detected at the DRB1 locus, which leads to the speculation that DRB1 may have a more important role in dealing with parasite infection for pandas. Coprological analyses revealed that 55.17% of individuals exhibited infection with 1–2 helminthes and 95.3% of infected pandas carried Baylisascaris shroederi. Using a generalized linear model, we found that Aime-DRB1*10 was significantly associated with parasite infection, but no resistant alleles could be detected. MHC heterozygosity of the pandas was found to be uncorrelated with the infection status or the infection intensity. These results suggested that the possible selection mechanisms in extant wild pandas may be frequency dependent rather than being determined by overdominance selection. Our findings could guide the candidate selection for the ongoing reintroduction or translocation of pandas.     

圈养大熊猫冷冻精液对其种群遗传多样性的作用

在过去34年的圈养大熊猫种群保护工作中,我们成功建立了全球最大的大熊猫精子库,目前已保存50只大熊猫个体总计7 000余支细管冷冻精液(冻精)。冷冻精液一方面可以使物种的遗传资源得到长久保存,另一方面可以通过人工授精的方式促进种群繁育。但是,圈养大熊猫冷冻精液对其种群遗传多样性的作用尚未有明确报道。本研究首先根据成都大熊猫繁育研究基地2000—2014年冷冻精液人工授精数据,对比分析了冻精人工授精个体和圈养种群的遗传多样性。结果显示,冻精人工授精个体遗传多样性均高于同年圈养种群的平均遗传多样性,表明在繁殖年份中冻精人工授精可以显著提高圈养大熊猫种群的遗传多样性。统计精子库中所有冻精个体的平均血缘系数并与圈养种群进行对比分析,探究冷冻精液对圈养种群遗传多样性的潜在作用。结果显示,精子库中有21只已死亡个体的精液,其中有66.67%的个体平均血缘系数低于圈养种群;有14只20岁以上个体的精液,其中有50.00%的个体平均血缘系数低于圈养种群;另有15只20岁以下个体的精液,其中有53.33%的个体平均血缘系数低于圈养种群,表明冷冻精液对圈养种群遗传多样性的保护具有重要价值。综上所述,冷冻精液不但有效保存了大熊猫遗传资源,而且在保护圈养种群遗传多样性方面具有积极的促进作用。     

Preliminary investigation on iodine nutrition in captive giant pandas

Background/ObjectiveThe giant panda belongs to the family Ursidae and, as a species of bear, still retains the simple digestive system of a Carnivoran. However, under the pressure of a specific habitat they had to adapt to a plant mono-diet consisting of bamboo with different species and growth stages around the year. A plant-based diet has relatively low iodine content with risk of iodine deficiency. Furthermore, bamboo contains cyanogenic glycosides releasing cyanide whose detoxification metabolite the thiocyanate acts as antagonist against iodine uptake and storage in the thyroid. To date very little is known about the iodine nutritional status of the giant panda, thus this study was conducted to receive the first information about the iodine nutrition of captive giant panda.Subjects/MethodsHere we investigated the iodine content of bamboo with different plant parts/vegetation stage and species and further compounds of the captive giant panda diet. Next, the urinary iodine (UI) and urinary thiocyanate (UT) levels of infant, sub-adult, adult and geriatric captive giant pandas was measured during the periods when the pandas consume both bamboo leaves- and culm (bamboo leaf-culm stage). Afterwards, the UI of 19 adult giant pandas was measured again for the different iodine intake during bamboo shoot stage. Finally, in this study part also the fecal iodine concentration was analyzed for calculation of total iodine excretion in relation to the iodine intake.ResultsBamboo leaves had the highest iodine content (453 μg/kg dry matter (DM)), followed by the shoots (84 μg/kg DM, p < 0.05), while bamboo culm had the lowest value (12 μg/kg DM, p < 0.05). During bamboo leaf-culm stage, giant pandas of different age groups had different UI and UT levels (p < 0.05). Furthermore, UI and UT were positively correlated among sub-adult, adult and geriatric giant pandas (p < 0.05). In adult giant pandas during bamboo shoot stage, the iodine excretion in feces was not different from that in urine while their total iodine excretion was less than their iodine intake (p < 0.05). Moreover, during bamboo shoot stage, the UI level of adult giant pandas was much lower than noted during bamboo leaf-culm stage (p < 0.05).ConclusionsOur results indicate that UI of captive giant pandas was related to their age as well as to the vegetation stage/part of bamboo they consumed reflecting a different periodic iodine supply. Thiocyanate and fecal excretion should be emphasized when considering the iodine nutrition of giant pandas.     

基于线粒体12S对中国两大山系大熊猫西氏贝蛔虫种群遗传多样性的分析

大熊猫是中国特有的珍稀濒危物种,而西氏贝蛔虫是大熊猫体内最为常见的肠道寄生虫,对野生和圈养大熊猫危害极大。考虑到大熊猫因分布区域的差异而形成了不同的亚种以及寄生虫与宿主间广泛的协同进化关系,西氏贝蛔虫是否也存在与大熊猫相适应的亚种分化一直是野生动物学家极其关注和热议的话题。为此,本文选择中国两大山系（岷山和邛崃）大熊猫种群体内共计34株西氏贝蛔虫虫体样本进行种群遗传多态性研究。利用PCR技术扩增出了岷山（14株）和邛崃（20株）西氏贝蛔虫的线粒体12S基因全序列并对其做了遗传多样性分析。结果表明：（1）34个样本包含9个单倍型,呈现出一个高单倍性多样性和低核苷酸多样性的特点;（2）负的Tajima's D和Fu's Fs中性检验值及“多峰型”的种群歧点分布图暗示种群不久前曾经历过突增长的现象;（3）低的种群间的分化系数和高的基因流表明两个地理种群间未形成显著的遗传分化;（4）系统发育树和单倍型网络图表明两山系种群分布无区域特异性。因此,岷山和邛崃山系的大熊猫体内的西氏贝蛔虫种群遗传变异性较低,分化不明显。该发现不仅暗示了西氏贝蛔虫与其宿主（大熊猫）的进化不同步,而且还为不同区域大熊猫西氏贝蛔虫病的监控提供了理论依据。     

No evidence for MHC‐based mate choice in wild giant pandas

Major histocompatibility complex genes (MHC), a gene cluster that controls the immune response to parasites, are regarded as an important determinant of mate choice. However, MHC‐based mate choice studies are especially rare for endangered animals. The giant panda (Ailuropoda melanoleuca), a flagship species, has suffered habitat loss and fragmentation. We investigated the genetic variation of three MHC class II loci, including DRB1, DQA1, and DQA2, for 19 mating‐pairs and 11 parent‐pairs of wild giant pandas based on long‐term field behavior observations and genetic samples. We tested four hypotheses of mate choice based on this MHC variation. We found no supporting evidence for the MHC‐based heterosis, genetic diversity, genetic compatibility and “good gene” hypotheses. These results suggest that giant pandas may not use MHC‐based signals to select mating partners, probably because limited mating opportunities or female‐biased natal dispersal restricts selection for MHC‐based mate choice, acknowledging the caveat of the small sample size often encountered in endangered animal studies. Our study provides insight into the mate choice mechanisms of wild giant pandas and highlights the need to increase the connectivity and facilitate dispersal among fragmented populations and habitats.     

Sixteen novel microsatellite loci developed for the giant panda (Ailuropoda melanoleuca)

Sixteen novel microsatellite DNA loci were developed from the giant panda (Ailuropoda melanoleuca) using a magnetic-bead capture method. A total of 115 alleles were obtained for these markers, ranging from 4 to 12 alleles per locus (average 7.188). These loci exhibited high levels of polymorphic information content and expected heterozygosity, 0.558–0.855 (average 0.729) and 0.628–0.885 (average 0.778), respectively. Therefore, the allelic polymorphism and heterozygosity show that the giant pandas raised in China Research and Conservation Center possess abundant genetic variation. In addition, if the three markers showing null alleles were excluded, the remaining 13 microsatellite loci still presented extremely low non-exclusion probabilities of parentage (0.002), paternity (0.000) and identity (0.000). As a result, this new suit of microsatellite markers would be a very informative tool for the genetic and conservation studies of giant pandas.     

微卫星标记在大熊猫研究中的应用进展

大熊猫是我国保护最为成功、研究最为深入的珍稀动物之一,可以为其它珍稀濒危物种的保护研究工作提供参考。20世纪70年代末期借助无线电颈圈,大熊猫的生态学研究工作取得了突破性进展,近20年来微卫星标记和非损伤性遗传取样技术的联合使用,将大熊猫的保护研究工作提升到一个崭新的高度。本文在综合所有已发表大熊猫微卫星标记的基础上,梳理了微卫星标记在圈养大熊猫亲子鉴定与遗传管理,野生大熊猫个体识别与种群数量调查、遗传多样性评估、扩散和种群遗传结构研究中的应用情况,并着重介绍了其中29个重要的微卫星标记。同时指出目前微卫星标记的使用存在标记选择不统一、等位基因读数无统一规程等问题,并对应用前景进行了前瞻。     

Domestication does not narrow MHC diversity in Sus scrofa

The Major Histocompatibility Complex (MHC) is a multigene family of outstanding polymorphism. MHC molecules bind antigenic peptides in the peptide-binding region (PBR) that consists of five binding pockets (P). In this study, we compared the genetic diversity of domestic pigs to that of the modern representatives of their wild ancestors, the wild boar, in two MHC loci, the oligomorphic DQA and the polymorphic DRB1. MHC nucleotide polymorphism was compared with the actual functional polymorphism in the PBR and the binding pockets P1, P4, P6, P7, and P9. The analysis of approximately 200 wild boars collected throughout Europe and 120 domestic pigs from four breeds (three pureblood, Pietrain, Leicoma, and Landrace, and one mixed Danbred) revealed that wild boars and domestic pigs share the same levels of nucleotide and amino acid polymorphism, allelic richness, and heterozygosity. Domestication did not appear to act as a bottleneck that would narrow MHC diversity. Although the pattern of polymorphism was uniform between the two loci, the magnitude of polymorphism was different. For both loci, most of the polymorphism was located in the PBR region and the presence of positive selection was supported by a statistically significant excess of nonsynonymous substitutions over synonymous substitutions in the PBR. P4 and P6 were the most polymorphic binding pockets. Functional polymorphism, i.e., the number and the distribution of pocket variants within and among populations, was significantly narrower than genetic polymorphism, indicative of a hierarchical action of selection pressures on MHC loci.     

圈养小熊猫遗传多样性与种群遗传结构分析

小熊猫是亚洲特有的珍稀濒危动物,目前受到栖息地减少、片断化和人类活动干扰等威胁。中国圈养小熊猫已经有60 多年历史,约55 个机构曾经饲养过小熊猫,现今圈养数量有400 多只,评估小熊猫圈养种群的遗传多样性和遗传结构对科学维持圈养种群和保存遗传种质资源意义重大。本研究利用19 个微卫星座位,对中国境内11 个小熊猫圈养种群的116 只个体进行了遗传多样性评估及遗传结构分析。结果显示11 个种群都具有较高的遗传多样性,平均基因丰富度3.505 ± 1.033 (北京)至4.026 ± 1.219 (冕宁),期望杂合度0.631 ± 0.225(黄山)至0.782 ±0.171 (温岭)。其中福州和无锡种群极显著偏离Hardy-Weinberg 平衡。整个圈养群体内各个种群遗传分化系数为0.055,呈显著分化,表明11 个种群遗传分化水平较高。Bayesian 遗传聚类分析将11 个种群聚为三个遗传簇,与野生种群的遗传聚类结果一致。结论:小熊猫圈养种群与野生种群相比,同样具有较高的遗传多样性。因此,圈养小熊猫遗传管理的重点不再是引进野生个体充实圈养种群,应制订科学的繁殖计划,避免近交,从而维持圈养种群的遗传多样性。     

圈养繁殖在大熊猫保护生物学中的作用

栖息地的日益恶化、片段化和种群孤岛化已成为保护大熊猫的严重障碍。尽管就地保护濒危物种是最有希望的措施,但对于大熊猫而言,研究表明圈养繁殖也是保护这一极危物和中的有效手段之一,对保护遗传多样性起以了重要作用,是延缓绝灭速率的保证。本文阐述了大圈养繁殖的必要性及其在保护生物学中的作用。     

应用微卫星分型方法进行大熊猫父亲鉴定

1979年我们利用水稻花粉单倍体植株的 茎节、幼穗、叶鞘和叶片等组织进行离体培养, 诱导出二倍体植株。     

Non-invasive giant panda paternity exclusion

Giant panda hair samples obtained by noninvasive methods served as a source of DNA for amplification of seven giant panda microsatellite loci utilizing the polymerase chain reaction. Thirteen giant pandas held in Chinese zoos were tested for identification of paternity. Some males listed as sires have been excluded as the biological father of captive-born giant pandas. Because of the death of some potential sires, paternity is still not assigned for some giant pandas, although there is a high likelihood that paternity assignment could be made if postmortem samples are available for genetic analysis. The DNA microsatellite variation assayed by the test we have developed provides a rapid, highly informative, and noninvasive method for paternity identification in giant pandas. © 1994 Wiley-Liss, Inc.