Patterns of Adaptive and Neutral Diversity Identify the Xiaoxiangling Mountains as a Refuge for the Giant Panda

Genetic variation plays a significant role in maintaining the evolutionary potential of a species. Comparing the patterns of adaptive and neutral diversity in extant populations is useful for understanding the local adaptations of a species. In this study, we determined the fine-scale genetic structure of 6 extant populations of the giant panda (Ailuropoda melanoleuca) using mtDNA and DNA fingerprints, and then overlaid adaptive variations in 6 functional Aime-MHC class II genes (DRA, DRB3, DQA1, DQA2, DQB1, and DQB2) on this framework. We found that: (1) analysis of the mtDNA and DNA fingerprint-based networks of the 6 populations identified the independent evolutionary histories of the 2 panda subspecies; (2) the basal (ancestral) branches of the fingerprint-based Sichuan-derived network all originated from the smallest Xiaoxiangling (XXL) population, suggesting the status of a glacial refuge in XXL; (3) the MHC variations among the tested populations showed that the XXL population exhibited extraordinary high levels of MHC diversity in allelic richness, which is consistent with the diversity characteristics of a glacial refuge; (4) the phylogenetic tree showed that the basal clades of giant panda DQB sequences were all occupied by XXL-specific sequences, providing evidence for the ancestor-resembling traits of XXL. Finally, we found that the giant panda had many more DQ alleles than DR alleles (33∶13), contrary to other mammals, and that the XXL refuge showed special characteristics in the DQB loci, with 7 DQB members of 9 XXL-unique alleles. Thus, this study identified XXL as a glacial refuge, specifically harboring the most number of primitive DQB alleles.     

一种从大熊猫粪便中提取DNA的改进方法

本研究描述一个改进的方法，使从大熊猫粪便中提取DNA用于PCR扩增变得更加容易。在粪便DNA的提取过程中采用一个新的预处理方法，将粪便用预冷的丙酮洗2～3次，除去粪便中含有的大量PCR抑制物，然后用蛋白酶K裂解、酚氯仿抽提，能提取到纯度很高的DNA供PCR扩增。本实验PCR扩增了大熊猫脑源性神经营养因子(BDNF)基因和线粒体细胞色素6基因片段，并进行测序分析，证实了提取的可靠性。对比本方法和未经丙酮预处理的方法提取的DNA进行PCR扩增，前者的扩增结果明显优于后者。     

粪便取样对亚洲黑熊脑源性神经营养因子基因全长序列分析

沿用本室改进的粪便提取方法,参照马来熊BDNF基因序列设计引物,首次从亚洲黑熊粪便DNA中扩增和克隆到包含完整核BDNF基因的753 bp片段,以毛发作阳性对照并进行重复实验,获得稳定一致结果。序列分析表明,亚洲黑熊的BDNF基因非常保守,与人相比,一致性达94.5%,与大熊猫比达98.9%。在推导的多肽序列中,其成熟区氨基酸序列与所有已报道哺乳动物的完全一致;对亚洲黑熊及其相关物种BDNF基因序列的比较分析,发现大熊猫与包括黑熊在内的熊科动物亲缘关系更近,而与小熊猫较远。文章首次采用非损伤性取样法在分子生物学水平对亚洲黑熊基因组核BDNF基因进行分析,不仅为亚洲黑熊的保护和繁育提供重要参考资料,为非损伤性取样在珍稀濒危野生动物研究中的应用拓宽了思路,也为亚洲黑熊及其近缘种的系统分类研究提供分子证据。     

Di‐, tri‐ and tetranucleotide microsatellite loci for the giant panda,Ailuropoda melanoleuca

We describe 10 polymorphic microsatellite loci for the giant panda, Ailuropoda melanoleuca. Microsatellite sequences were isolated from three partial genomic libraries of giant panda DNA that were enriched for (i) (GT), (ii) (GAA) & (CAA), and (iii) (GATA) repeat sequences. The markers were tested for polymorphism in up to 82 pandas. Number of alleles at each locus varied between four and 11, and the observed and expected heterozygosities varied between 0.267 and 0.732, and between 0.601 and 0.799, respectively.     

大熊猫和小熊猫粪便DNA提取的简易方法

采集了大熊猫和小熊猫的新鲜粪便样品 ,使用 1 0 0 %乙醇保存。通过重复离心富集研究动物的肠道脱落细胞 ,并使用乙醇和双蒸水洗涤以除去抑制物。用 1 %的SDS快速裂解细胞 ,离心除去残渣后 ,向裂解液中加入蛋白酶进行消化。消化结束后使用等体积的酚 /氯仿抽提 ,乙醇沉淀DNA。用双蒸水溶解粪便DNA后 ,使用PCR产物纯化试剂盒对粪便DNA进行纯化。电泳检测结果显示 ,从乙醇保存的大、小熊猫粪便样品中抽提到高质量的粪便DNA。对线粒体控制区、细胞色素b基因、 1 2SrRNA基因的PCR扩增反应以及测序结果也证实了样品保存方法和DNA抽提方法可靠而高效。此方法使用实验室内常用的分子生物学试剂 ,不仅克服了分子粪便学研究中常见的抑制物粪便DNA微量降解严重等障碍 ,与商业化的粪便抽提试剂盒 (QIAampDNAStoolMiniKit,Qiagen)相比还是一种经济的试验方法 (抽提反应成本为试剂盒的 1 / 5 )。文中还对粪便DNA内细菌基因组等背景DNA可能对分子粪便学试验结果的影响进行了探讨。在基于PCR技术的遗传学研究中 ,对于植食性动物而言 ,粪便内的背景DNA对目标动物DNA片断的扩增和序列测定未见影响 ;但对于肉食性动物 ,则必须考虑被捕食者基因组对试验可能产生的影响 ,应谨慎对待     

Characterization of a single nucleotide polymorphism in the coding sequence of the bovine transferrin gene.

A single nucleotide polymorphism was identified in the coding sequence of the bovine transferrin gene. Two alleles (SSCP1 and SSCP2) were detected by SSCP analysis and the mutation point was identified and confirmed by direct sequencing of the PCR products. The relationship between protein and DNA polymorphism was established. Protein variants A, D1 and E correspond to SSCP allele 1 and variant D2 corresponds to SSCP allele 2. DNA sequences from genotypes AA, AE, AD2, D1E, D2E and D2D2 reveal an A/G substitution at position 1455 of the cDNA which causes a Gly/Glu substitution which could be responsible for the mobility difference between D1 and D2 variants. Because of the number of variants, this suggests that other SNPs exist in the bovine transferrin gene. A linkage analysis between the SSCPs and two microsatellites (UWCA46 and CSSM019) mapped the transferrin gene to BTA1. Two-point analysis revealed a tight linkage within the transferrin protein variants and the SSCPs.     

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

以中国最大的大熊猫圈养种群—四川卧龙中国大熊猫保护中心的圈养种群为对象,以８个大熊猫微卫星位点为分子标记, 探讨了大熊猫圈养种群的遗传多样性, 并与邛崃野生种群及其他７个濒危物种进行比较。微卫星数据表明, 圈养种群的遗传多样性水平(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。因此继续增加野生个体对保持遗传多样性的意义已经不大, 建议该圈养种群的保护策略应将重点放到制定更有效的繁殖计划以避免近交上。     

The transcription of MGAT4A glycosyl transferase is increased in white cells of peripheral blood of Type 2 Diabetes patients

Background

The giant panda (Ailuropoda melanoleuca) is one of the most endangered animals due to habitat fragmentation and loss. Although the captive breeding program for this species is now nearly two decades old, researches on the genetic background of such captive populations, especially on adaptive molecular polymorphism of major histocompatibility complex (MHC), are still limited. In this study, we characterized adaptive variation of the giant panda's MHC DQA gene by PCR amplification of its antigen-recognizing region (i.e. the exon 2) and subsequent single-strand conformational polymorphism (SSCP) and sequence analyses.

Results

The results revealed a low level of DQA exon 2 diversity in this rare animal, presenting 6 alleles from 61 giant panda individuals. The observed polymorphism was restricted to 9 amino acid substitutions, all of which occurred at and adjacent to positions forming the functionally important antigen-binding sites. All the samples were in Hardy-Weinberg proportions. A significantly higher rate of non-synonymous than synonymous substitutions at the antigen-binding sites indicated positive selection for diversity in the locus.

Conclusion

The DQA allelic diversity of giant pandas was low relative to other vertebrates. Nonetheless, the pandas exhibited more alleles in DQA than those in DRB, suggesting the alpha chain genes would play a leading role when coping with certain pathogens and thus should be included in conservation genetic investigation. The microsatellite and MHC loci might predict long-term persistence potential and short-term survival ability, respectively. Consequently, it is recommended to utilize multiple suites of microsatellite markers and multiple MHC loci to detect overall genetic variation in order to design unbiased conservation strategies.     

The faecal flora of the giant panda (Ailuropoda melanoleuca)

The faecal floras of two adult (male and female) and one infant (male) giant panda kept at the Ueno Zoo, Tokyo, Japan were examined and shown to be quite different from those of other animals. The predominant bacteria in the adults were Streptococcus (including Enterococcus) and Enterobacteriaceae, while obligate anaerobes had minor populations. Fastidious anaerobes were not detected. The predominant bacteria in the suckling infant were Lactobacillus and Streptococcus, followed by Bifidobacterium. After the infant began to eat bamboo leaves the number of Lactobacillus decreased and Bifidobacterium became undetectable, whereas Enterobacteriaceae became one of the most predominant flora. The most dominant streptococcus isolated from the female panda was identified as Streptococcus bovis, but those from the male adult and the weaned infant were not identified as any known species.     

The faecal flora of the giant panda (Ailuropoda melanoleuca)

The faecal floras of two adult (male and female) and one infant (male) giant panda kept at the Ueno Zoo, Tokyo, Japan were examined and shown to be quite different from those of other animals. The predominant bacteria in the adults were Streptococcus (including Enterococcus ) and Enterobacteriaceae, while obligate anaerobes had minor populations. Fastidious anaerobes were not detected. The predominant bacteria in the suckling infant were Lactobacillus and Streptococcus , followed by Bifidobacterium. After the infant began to eat bamboo leaves the number of Lactobacillus decreased and Bifidobacterium became undetectable, whereas Enterobacteriaceae became one of the most predominant flora. The most dominant streptococcus isolated from the female panda was identified as Streptococcus bovis , but those from the male adult and the weaned infant were not identified as any known species.     

A reliable, non-invasive PCR method for giant panda (Ailuropoda melanoleuca) sex identification

We developed an inexpensive, fast and reliable PCR method for sex identification of giant panda (Ailuropoda melanoleuca) by using one pair of primers to co-amplify homologous fragments with size polymorphism that located at amelogenin (AMEL) exon 5. In giant panda, a 63 bp deletion in exon 5 of Y-linked allele provides a significant discrimination between AMELX and AMELY, thus the amplification products can be distinguished simply by agarose gel electrophoresis, exhibiting sex-specific banding patterns (male: 237 bp, 174 bp; female: 237 bp). Both blood and feces samples from known-sex giant pandas were successfully amplified. Cross species test also revealed that this method could be applied to other Ursidae species. These authors contributed equally to this work.     

同域分布的珍稀野生动物对放牧的行为响应策略

保护区内放牧活动对野生动物保护存在负面影响,明确不同物种对放牧干扰的行为响应对制定更有针对性的保护管理政策具有重要意义。使用红外相机研究卧龙自然保护区放牧活动对多种珍稀野生动物的影响,分析放牧激励政策实施前后大熊猫(Ailuropoda melanoleuca)及其同域分布的小熊猫(Ailurus fulgens)、川金丝猴(Rhinopithecus roxellana)、水鹿(Rusa unicolor) 4种珍稀野生动物的照片数、空间分布以及活动模式的变化,探讨这4种动物对放牧的行为响应策略。结果表明:(1)一期(2012—2013),2012年10月实施了禁马政策,同年12月实施放牧(牛羊)激励政策)家畜照片数量很少,4种野生动物照片数相对较多。二期(2014—2015)家畜的照片数显著增加(P0.01),小熊猫(P0.05)与川金丝猴(P0.01)的照片数均显著减少,大熊猫、水鹿的照片数也呈减少趋势;到三期(2016—2017),大熊猫、小熊猫及水鹿3种关注野生动物的照片数基本回升到激励政策实施前的水平,无川金丝猴照片记录。(2)一期,4种野生动物在研究区域有较广的分布;二期,大熊猫、小熊猫的空间分布范围均缩小,无川金丝猴空间分布信息,而家畜、水鹿的空间分布范围有所增加;到三期,大熊猫、小熊猫的空间分布基本恢复到放牧激励政策实施前的区域,无川金丝猴的空间分布信息。(3)放牧激励政策实施前后,大熊猫、小熊猫及川金丝猴活动模式无明显变化,但水鹿的活动更加集中于傍晚,以避开人类与家畜的活动高峰。同域分布的不同的野生动物对人类活动(如放牧)的行为响应策略不同,各保护区在制定相关保护政策时应综合考虑人类干扰对多个物种的影响,增加决策的科学性与合理性。     

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.     

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.     

Phylogenetic study of Baylisascaris schroederi isolated from Qinling subspecies of giant panda in China based on combined nuclear 5.8S and the second internal transcribed spacer (ITS-2) ribosomal DNA sequences

The nuclear ribosomal DNA (rDNA) region spanning 5.8S rDNA and the second internal transcribed spacer (ITS-2) of Baylisascaris schroederi isolated from the Qinling subspecies of giant panda in Shaanxi Province, China were amplified and sequenced. Sequence variations in the two rDNA regions within B. schroederi and among species in the family Ascarididae were examined. The lengths of B. schroederi 5.8S and ITS-2 rDNA sequences were 156 bp and 327 bp, respectively, and no nucleotide variation was found in these two rDNA regions among the 20 B. schroederi samples examined, and these ITS-2 sequences were identical to that of B. schroederi isolated from giant panda in Sichuan province, China. The inter-species differences in 5.8S and ITS-2 rDNA sequences among members of the family Ascarididae were 0-1.3% and 0-17.7%, respectively. Phylogenetic relationships among species in the Ascarididae were re-constructed by Bayesian inference (Bayes), maximum parsimony (MP), and maximum likelihood (ML) analyses, based on combined sequences of 5.8S and ITS-2 rDNA. All B. schroederi samples clustered together and sistered to B. transfuga with high posterior probabilities/bootstrap values, which further confirmed that nematodes isolated from the Qinling subspecies of giant panda in Shaanxi Province, China represent B. schroederi. Because of the large number of ambiguously aligned sequence positions (difficulty of inferring homology by positions), ITS-2 sequence alone is likely unsuitable for phylogenetic analyses at the family level, but the combined 5.8S and ITS-2 rDNA sequences provide alternative genetic markers for the identification of B. schroederi and for phylogenetic analysis of parasites in the family Ascarididae.     

Extensive polymorphism of the BOLA-DRB3 gene distinguished by PCR-RFLP

A polymerase chain reaction (PCR)-based method is described for typing of alleles of the bovine lymphocyte antigen (BoLA)-DRB3 gene. A total of 30 DRB3 alleles were distinguished by digestion of PCR amplification products of BoLA-DRB3 exon 2 with RsaI, BstYI and HaeIII (PCR-RFLP). All restriction fragment patterns, with the exception of one HaeIII pattern, were consistent with restriction sites that were found among 14 previously sequenced DRB3 alleles. The PCR-RFLP typing method was evaluated on 168 genomic DNA samples collected from animals of 10 cattle breeds, 48 of which were typed in the Fourth International BoLA Workshop for BoLA-DRB and -DQ by conventional restriction fragment length polymorphism (RFLP) analysis using heterologous and homologous DNA probes. Thirty-one DRB/DQ haplotypes containing 23 DRB3 alleles were identified among the 48 workshop animals analysed. Using PCR-RFLP, 11 DRB3 alleles were identified in 18 workshop animals for which DRB RFLPs were not informative. PCR-RFLP typing of additional animals revealed five new DRB3 alleles, of which three contained a putatively located three basepair deletion in the identical position as found for the sequenced allele DRB*2A. PCR-RFLP was shown to be a rapid and sensitive method for the detection of polymorphism in a functionally relevant domain of the BoLA-DRB3 gene and should be useful for studying the evolution of DRB polymorphism in cattle and other Bovidae.     

黄麂Mhc-DRB 基因多态性及其维持机制

利用牛DRB3 特异性引物(LA31 和LA32),通过聚合酶链式反应(PCR)、单链构象多态性(SSCP)以及克隆测序技术,从12 只黄麂个体中共获得20 个DRB 第二外显子等位基因,其中6 个个体具有3 ～ 4 个等位基因,提示利用该引物从黄麂中至少可以扩增出2 个DRB 位点。所有序列均无插入、缺失和终止密码子。基于序列比对(与牛DRB3 和鹿科DRB 基因同源性非常高),以及所检测到的氨基酸变异位点主要位于抗原结合区,推测本文所获得的黄麂序列为表达的、且具有重要功能的DRB 位点。抗原结合区氨基酸位点的非同义替换(d_N )显著大于同义替换(d_S )(P < 0.01),说明历史上黄麂DRB 基因经历过正选择作用。CODMEL 程序中的模型M7 和M8 似然比检测(Likelihood ratio test,LRT)结果同样支持上述推论。进一步利用经验贝叶斯法准确地检测出6 个受正选择作用的氨基酸位点(位点11、37、61、67、71、86),其中的5 个位点位于PBR 区。因此,正选择作用可能是维持黄麂DRB 基因多态性的主要机制之一。基于DRB 外显子2 序列利用邻接法(NJ)
构建了部分偶蹄动物系统发生关系,在NJ 树上,黄麂DRB 基因与其它鹿科动物DRB 基因呈镶嵌式分布,提示跨物种进化是维持黄麂DRB 基因多态性的另一重要机制。此外,黄麂两个等位基因(Mure-DRB1 和Mure-DRB11)和马鹿的两个等位基因(Ceel-DRB34 和Ceel-DRB46)与牛科的等位基因构成一个独立的进化枝,说明黄麂和马鹿的某些DRB 基因具有非常古老的谱系。     

大熊猫及其近缘种rDNA序列变异和系统进化关系

应用ｒＤＮＡ间隔区Ｓｏｕｔｈｅｒｎ转换技术研究大熊猫及其近缘种的分子系统关系。通过比较大熊猫、小熊猫、黑熊、马来熊、浣熊和猞猁的ｒＤＮＡ间隔区限制性内切酶图谱,用最大似然法和简约法构建它们的分子系统树。结果表明大熊猫与熊具有较近的亲缘关系,与小熊猫和浣熊的亲缘关系较远。     

大熊猫生长激素受体(GHR) cDNA 的克隆与序列分析

根据已报道的若干物种GHR 基因cDNA 序列设计引物, 利用RT- PCR 技术首次从大熊猫肝脏组织总RNA中扩增出GHR 基因编码区全长cDNA 序列, 克隆于pGEM®-T 载体后进行测序和序列分析。结果表明,大熊猫GHR 的ORF为1 917 bp , 编码638 个氨基酸的前体蛋白, 由18 个氨基酸的信号肽和620 个氨基酸的成熟肽组成,与人、狗、猪GHR 结构相似, 大熊猫GHR 成熟肽由246 个氨基酸的胞外区、24 个氨基酸的跨膜区和350 个氨基酸的胞内区组成, 并具GHR 的特征性结构。序列相似性比较显示, 大熊猫GHR 与哺乳类GHR 具有69 %～93 %的高序列相似性, 与爬行类和鸟类的序列相似性也达到60 % , 而与鱼类的序列相似性较低, 仅为30 %左右。与其它哺乳动物GHR 相比, 大熊猫GHR 在氨基酸序列上也存在明显的特异性。     

Comparison of two selective media for the detection and enumeration of Lactobacilli in human faeces

The enumeration of faecal bacteria is an important requirement for many studies of bowel health. One approach is the use of selective culture media for the culture and identification of genera or species from faeces. This study compares the culture of Lactobacilli from dilution series of faecal samples from six healthy human volunteers on two commonly used media, LAMVAB and Rogosa agar. Colonies were counted after a 72-h anaerobic incubation at 37 degrees C, and colony morphology recorded by a single observer. DNA was isolated from a representative number of colonies and genus-specific PCR, single-stranded conformation polymorphism (SSCP) and DNA sequencing performed. Total colony counts ranged from <3.00 to 7.48 log(10) cfu/g of faeces for LAMVAB and 5.09 to 7.66 log(10) cfu/g for Rogosa. For each subject, the total colony count was higher on Rogosa than that obtained with LAMVAB agar. SSCP analysis and DNA sequencing indicated that colony morphology was not an accurate predictor of genus identity. Growth of two species, Lactobacillus acidophilus and Lactobacillus gasseri, was not supported on LAMVAB medium. Rogosa agar was more likely to support growth of non-Lactobacillus species. Therefore, neither medium gave a fully accurate representation of the Lactobacilli species present in human faecal samples.