扩展功能
文章信息
- 涂洪梅, 周闯, 王冠楠, 成美玲, 岳碧松, 孟杨
- TU Hongmei, ZHOU Chuang, WANG Guannan, CHENG Meiling, YUE Bisong, MENG Yang
- 基于转录组测序的大熊猫多态性微卫星标记筛选
- Development of Polymorphic Microsatellite Markers for Ailuropoda melanoleuca Based on RNA-Sequencing
- 四川动物, 2020, 39(1): 15-22
- Sichuan Journal of Zoology, 2020, 39(1): 15-22
- 10.11984/j.issn.1000-7083.20190276
-
文章历史
- 收稿日期: 2019-08-07
- 接受日期: 2019-11-12
大熊猫Ailuropoda melanoleuca是全球生物多样性保护的旗舰物种(Wei et al., 2012)。由于受到人类人口扩张和栖息地丧失、破碎化等的影响,大熊猫的生存面临严重威胁(Huang et al., 2015)。近几十年来,我国为保护这一珍稀物种采取了一系列的保护措施,使得全国大熊猫种群得以逐渐复壮。全国第四次大熊猫数量调查显示,截止2014年底,野生大熊猫数量为1 864只,与第三次调查相比增长了16.8%(唐小平等,2015)。如今大熊猫栖息地存在不同程度的破碎化,野生大熊猫种群被进一步分割,部分小种群还面临着灭绝的风险,大熊猫的生存状况不容乐观,因此对大熊猫的保护及遗传多样性研究仍是重点(青菁,2016)。
微卫星标记称为简单重复序列(simple sequence repeats,SSRs)或短串联重复序列,由核心序列和侧翼序列组成,其中核心序列由1~6 bp核苷酸序列组成,在真核生物基因组中分布广泛。微卫星标记多态性含量丰富、易检测、小片段或者降解的DNA都能因含有足够的微卫星而进行PCR扩增,同时微卫星的检测在种间、种外都能进行(张正义等,2017)。作为分子生物学中最常用的一种遗传标记(Ma et al., 2004;Selkoe & Toonen,2006),微卫星分子标记被广泛应用于物种鉴定、群体遗传学和遗传图谱的研究(Huang et al., 2015;修云芳等,2015;Wang et al., 2016)。二代测序(next-generation sequencing,NGS)技术又被称为高通量测序,近年来应用十分广泛,使用二代测序技术开发微卫星标记的研究日趋增多。传统的微卫星分子标记开发方法,如DNA文库筛选和利用亲缘物种跨种扩增等往往耗时长、数量少、工作量大。与传统方法相比,使用NGS筛选微卫星分子标记成本更低、效益更高(Zheng et al., 2013;Liu et al., 2017)。因此,高通量测序由于其高效性迅速成为开发DNA分子标记的首选(岳华梅等,2016)。
从GenBank下载大熊猫完整基因组,使用本实验室自主开发的Krait(Du et al., 2017)对基因组中微卫星序列进行搜索和统计。结合本实验室已测得的6只大熊猫转录组数据,使用CandiSSR(Xia et al., 2016)对基因组和转录组中多态性微卫星进行对比与筛选。CandiSSR是一种能快速有效地进行多态性微卫星筛选的软件,能够基于多个组装序列鉴定出具有多态性的微卫星。本研究根据CandiSSR鉴定出的大熊猫基因组中的多态性微卫星,并挑选其中评分较为理想的20个,提取25只大熊猫血液DNA为模板进行PCR扩增,将扩增产物送至成都擎科梓熙生物科技有限公司分型并分析结果,验证所选位点的多态性及其他遗传参数。本研究结果开发了更多优良微卫星位点,为大熊猫种群遗传学的进一步研究提供了资源。
1 研究方法 1.1 基因组中的微卫星统计分析GenBank下载大熊猫全基因组(AilMel_1.0)。微卫星的检索标准分别限定为单碱基重复12次以上,二碱基重复7次以上,三碱基重复5次以上,四、五和六碱基重复4次以上,这与之前使用Krait进行研究的方法一致(Qi et al., 2015;Xu et al., 2016)。在本研究中,重复单位模式为循环排列和/或反向互补的被归类为一种SSR类型。例如,AGC模式由AGC、GCA、CAG、GCT、TGC和CTG组成(Jurka & Pethiyagoda,1995)。
1.2 RNA-seq数据分析采集6只大熊猫个体的RNA-seq数据。使用FastQC(https://www.bioinformatics.babraham.ac.uk/projects/fastqc)对reads进行质控,使用Cutadapt(Martin,2015)去除接头序列和较低质量reads(Phred score < 20)。使用Trinity(Haas et al., 2013)在默认参数下对clean reads进行组装,保留长度超过200 bp的转录本进行分析。采用CD-HIT-EST(Li & Godzik,2006)删除带有默认参数的冗余序列,并使用UniGenes发现微卫星标记。
1.3 多态性微卫星检测及引物设计基于本实验室已组装好的6只大熊猫的转录组序列文件和从GenBank下载的大熊猫全基因组序列,使用CandiSSR(Xia et al., 2016)检索多态性微卫星并由此设计引物。从筛选结果中选择20个多态性二碱基微卫星位点,用于25个大熊猫血液DNA样品的多态性验证。每个引物对(表 4)的正向引物在5’-末端用FAM荧光染料标记进行基因分型分析。
1.4 样品采集和DNA提取25只大熊猫个体的血液样本来自成都大熊猫繁育研究基地。使用DNA Extraction Kit(Tsingke,北京)试剂盒,并按照说明书从血液样品中提取基因组总DNA。以提取的基因组DNA为模板进行PCR扩增。
1.5 PCR条件和扩增PCR反应体系为20 μL:2.5 μL 10×Taq缓冲液,1.0~3.0 μL dNTP(各2.5 mmol·L-1),各引物1 μL(10 mmol·L-1),0.2 μL ExTaq聚合酶(5 U·mL-1;TaKaRa,日本),基因组总DNA约200 ng作为模板。PCR反应程序为:95 ℃ 5 min;94 ℃ 30 s,55~63 ℃ 30 s,72 ℃ 12 s,35个循环;72 ℃ 10 min,10 ℃保存。
1.6 基因分型所有样品的基因分型检测在成都擎科梓熙生物科技有限公司进行。基因分型时将每个荧光引物扩增的产物单独放在一条泳道里电泳分型。所有样品的基因分型在ABI 3730 DNA Analyzer上进行,使用GeneMapper 4.0检测各样本等位基因数,等位基因大小相对于分子内标GS500LIZ决定。
1.7 微卫星数据分析使用CERVUS(Marshall et al., 1998)计算等位基因数、观察杂合度、期望杂合度和多态信息含量(polymorphic information content,PIC)。使用GENEPOP(Raymond & Rousset,1995)计算哈迪-温伯格平衡和连锁不平衡的偏差。作为衡量多态信息含量的指标,PIC>0.50为高度多态性位点,0.25 < PIC < 0.50为中度多态性位点,PIC < 0.25为低度多态性位点(Botstein et al., 1980;Shete et al., 2000;Nagy,2012)。
2 结果 2.1 微卫星类型和丰度Krait搜索和统计结果显示,在大熊猫的基因组中,不完美的微卫星是最常见的类型,共有3 177 509个,相对丰度为1 415.17 loci/Mb,相对密度为40 727.66 bp/Mb。其次是完美微卫星,数量最少的是复合型微卫星,只有53 787个,相对丰度为23.96 loci/Mb,相对密度为980.90 bp/Mb(表 1)。完美微卫星的数量、长度、相对丰度和相对密度如表 2所示。6种完美微卫星总数为879 113个,总长度为15 386 054 bp。其中最常见的类型是单碱基微卫星,有415 846个,且长度最长,达到了6 176 980 bp,最高丰度为185.21 loci/Mb,最高密度为2 751.06 bp/Mb,占基因组微卫星总数的47.30%。其次是二碱基微卫星,占微卫星总数的25.74%。而三碱基、四碱基和五碱基微卫星较少,六碱基微卫星则最少,仅占微卫星总数的0.34%。在大熊猫基因组中前10个最常见的重复单元是(A)n、(AC)n、(AG)n、(AAAT)n、(AAAG)n、(AT)n、(C)n、(AAAC)n、(AAGG)n和(AAT)n。
完美微卫星 | 复合微卫星 | 不完美微卫星 | |
微卫星数量 | 879 113 | 53 787 | 3 177 509 |
相对丰度/(loci/Mb) | 391.53 | 23.96 | 1 415.17 |
相对密度/(bp/Mb) | 6 852.52 | 980.90 | 40 727.66 |
类型 | 数量 | 长度/bp | 占比/% | 平均长度/bp | 相对丰度/(loci/Mb) | 相对密度/(bp/Mb) |
单碱基 | 415 846 | 6 176 980 | 47.30 | 14.85 | 185.21 | 2 751.06 |
二碱基 | 226 315 | 4 269 300 | 25.74 | 18.86 | 100.79 | 1 901.43 |
三碱基 | 38 434 | 679 929 | 4.37 | 17.69 | 17.12 | 302.82 |
四碱基 | 167 994 | 3 458 168 | 19.11 | 20.59 | 74.82 | 1 540.17 |
五碱基 | 27 570 | 724 445 | 3.14 | 26.28 | 12.28 | 322.65 |
六碱基 | 2 954 | 77 232 | 0.34 | 26.14 | 1.32 | 34.40 |
CandiSSR分析结果显示,从大熊猫转录组与基因组中共筛选出326个多态性微卫星,其中二碱基最多,占69.94%,共228个。三、四、五、六碱基所占比例分别为9.51%、14.11%、5.21%、1.22%。本研究根据标准差、缺失率以及序列长短等因素,选出20个二碱基微卫星用于后续的多态性评估验证(表 3)。
重复单元 | 标准差 | 缺失率 |
(AC)15 | 3.80 | 0.25 |
(CT)8 | 2.75 | 0.25 |
(AT)10 | 2.10 | 0.12 |
(TG)10 | 1.96 | 0.12 |
(CA)15 | 1.89 | 0.25 |
(GT)15 | 1.57 | 0.25 |
(GT)8 | 1.36 | 0.12 |
(CT)8 | 1.34 | 0.25 |
(GT)13 | 1.25 | 0.25 |
(TG)15 | 1.16 | 0.12 |
(AT)12 | 1.12 | 0.25 |
(CA)7 | 1.11 | 0.25 |
(AC)7 | 1.00 | 0.25 |
(CA)9 | 0.99 | 0.12 |
(CA)7 | 0.76 | 0.25 |
(AT)7 | 0.70 | 0.12 |
(CA)8 | 0.7 | 0.12 |
(AT)10 | 0.69 | 0.25 |
(GT)7 | 0.69 | 0.25 |
(CA)8 | 0.64 | 0.12 |
选出的20个微卫星位点能够在25只大熊猫个体的血液DNA中成功扩增。每个位点的等位基因数为2~8。观察杂合度为0~1.00,平均值为0.472,期望杂合度为0.280~0.784,平均值为0.532。使用Bonferroni校正后,有4个微卫星位点显著偏离哈迪-温伯格平衡(P < 0.01),20个微卫星位点未观察到显著的连锁不平衡(P>0.01)。20个位点多态信息含量为0.246~0.734,其中具有高度多态性的位点9个(PIC>0.50),11个位点呈中度多态性(0.25 < PIC < 0.50)(表 4)。
位点 | 引物序列(5’-3’) | 重复单元 | 等位基因数量 | 长度/bp | 观察杂合度 | 期望杂合度 | P | 多态信息含量 | Fis estimates | |
W & C | R & H | |||||||||
GP1 | F:CCCCACTACGCTCTTACACC R:TATTTGTGTCCCGGCCAGTC |
(AC)15 | 7 | 172 192 | 0.52 | 0.784 | 0.002 0 | 0.734 | 0.341 1 | 0.269 0 |
GP2 | F:GGTTGAGTGAGTTCTGGTTTGC R:CACTGTGCCCACATCAGAGT |
(AC)7 | 3 | 144 148 | 0.32 | 0.523 | 0.016 9 | 0.448 | 0.393 4 | 0.239 2 |
GP3 | F:GCAGGTTCCCTCGACAACAA R:TCCCCACCTGTTTTGCAGTC |
(AT)10 | 4 | 178 184 | 0.64 | 0.607 | 0.678 9 | 0.547 | -0.054 9 | -0.056 1 |
GP4 | F:TGTGCGTATATCTATGTGTG R:TTCTGTTGCATCTGTTAGCA |
(AT)10 | 5 | 152 160 | 0.64 | 0.656 | 0.312 5 | 0.576 | 0.024 1 | -0.016 2 |
GP5 | F:CCTGAATCTGTCCTGGGGAAG R:GGTGAGGGAGGGATCCAAGA |
(AT)12 | 3 | 178 182 | 0.32 | 0.280 | 1.000 0 | 0.246 | -0.146 3 | -0.076 4 |
GP6 | F:TTTCTCTCGTGCAGCCAGTT R:TGACTGCACATACCAAGTGG |
(AT)7 | 2 | 138 142 | 0.32 | 0.372 | 0.589 0 | 0.298 | 0.142 9 | 0.146 2 |
GP7 | F:TCATTCCCAGAGCCAGCTTG R:TGGGTGTGTGTGATTGCACT |
(CA)10 | 3 | 194 198 | 0.28 | 0.520 | 0.016 4 | 0.399 | 0.466 7 | 0.251 7 |
GP8 | F:TGTCACCAGTCAGGCTAACG R:GTCTAGGCTCTCCCCTTCCA |
(CA)15 | 4 | 174 182 | 0.80 | 0.665 | 0.580 7 | 0.606 | -0.207 5 | -0.188 6 |
GP9 | F:GGCGAGGGTTCCGTATTTCT R:GAGGGCTGGGAAGAGAGGAA |
(CA)7 | 2 | 186 190 | 0 | 0.327 | 0 | 0.269 | 1.000 0 | 1.041 7 |
GP10 | F:GTCCCTTTGGCTGTGTGAGA R:GGAACAGGGGTGTTTGGGAA |
(CA)8 | 3 | 182 196 | 0.44 | 0.482 | 0.026 2 | 0.393 | 0.089 7 | 0.497 2 |
GP11 | F:AGGAGTCACGAGGGTAGCTG R:CCACGGTTAGAGCAGGTAGG |
(CA)8 | 4 | 194 200 | 0.24 | 0.675 | 0 | 0.595 | 0.649 2 | 0.729 0 |
GP12 | F:ATGAGTTCGAGCCCCATGTG R:ATTGCTGGCCCCTCTGACT |
(CA)9 | 3 | 192 200 | 0.48 | 0.545 | 0.039 4 | 0.428 | 0.122 0 | 0.527 7 |
GP13 | F:CCTGTGTTGGGCTCCATTCT R:CGTTCTGACGTAAAACCAGGA |
(CT)8 | 4 | 180 194 | 0.52 | 0.621 | 0.074 5 | 0.535 | 0.165 8 | 0.060 2 |
GP14 | F:GCTCTATTCCCTTGGCAGCA R:TACCATCTGGCAGGGTAGCT |
(CT)8 | 2 | 122 128 | 0.28 | 0.301 | 1.000 0 | 0.252 | 0.071 8 | 0.073 4 |
GP16 | F:AAGCATTACCTTCCCCGGAA R:TTGTGTGTATGTGTTCACGT |
(GT)13 | 3 | 154 164 | 0.36 | 0.402 | 0.702 5 | 0.331 | 0.107 4 | 0.066 8 |
GP17 | F:CTCTGACCAGCGAGTCCTTG R:TGCGAAAGGAGAAAACTCCAGT |
(GT)15 | 8 | 122 144 | 1.00 | 0.773 | 0.000 8 | 0.722 | -0.301 5 | -0.125 0 |
GP18 | F:ACTGCTCAAATCACACTTCT R:TGGCTTCTGTTGTAGACTCAGA |
(GT)7 | 3 | 162 166 | 0.48 | 0.431 | 1.000 0 | 0.361 | -0.116 3 | -0.064 2 |
GP20 | F:GTCATGGTGTGCGTGCATG R:ACTGCTCTGTGCTTGACTCC |
(GT)8 | 2 | 192 194 | 0.48 | 0.372 | 0.275 8 | 0.298 | -0.297 3 | -0.301 5 |
GP25 | F:TGACAGCCAAGTGGAAAGAGT R:CTTGCTGCCACACTCTGGTA |
(TG)10 | 5 | 194 204 | 0.56 | 0.616 | 0.130 7 | 0.562 | 0.093 1 | 0.297 7 |
GP26 | F:TTGGGGTCAGAGGGAAGTGA R:AGTGCTGACCGTGAGTTCTG |
(TG)15 | 4 | 184 192 | 0.76 | 0.680 | 0.161 3 | 0.602 | -0.120 4 | 0.082 9 |
转录组平台的构建在很大程度上推动了DNA分子标记的开发。近年来利用转录组数据获得含有微卫星的序列,并对其进行遗传多样性的研究在国际上已有很多成功报道(Tian et al., 2014;Xiao et al., 2015),而利用微卫星标记进行野生大熊猫遗传多样性分析从2001年就已开始(Lu et al., 2001)。在二代测序技术广泛应用之前,开发微卫星标记的方法较多,比较常用的方法包括构建DNA文库筛选法、小片段DNA克隆法、富集法、利用亲缘物种筛选等。Shen等(2005)利用Dynal磁珠法筛选了37对大熊猫微卫星引物,Wu等(2008)使用富集法筛选出33个微卫星位点,Li等(2010)通过构建DNA文库法筛选出7个稳定的四碱基微卫星位点。这些位点的筛选与标记开发在一定程度上可以用于大熊猫种群遗传多样性评估、遗传结构分析等相关研究。由于部分微卫星位点扩增不稳定,在使用粪便或者毛发等样本提取的DNA扩增效果不好时就限制了这些位点的应用。随着测序成本的降低及测序平台通量的不断扩大,利用二代测序技术对基因组进行随机测序即可获得海量的微卫星序列。与传统微卫星标记开发相比,二代测序技术对样品的处理和使用简单、处理速度快、通量高,能在较短时间内开发出数量多且类型丰富的微卫星。目前利用二代测序技术对微卫星位点进行筛选与开发也越来越多(Ritchie et al., 2016;李薇等,2017;宋琪等,2019)。
本研究对大熊猫基因组中的微卫星序列进行搜索和统计,使用的是本实验室开发的Krait:一种用于全基因组微卫星调查和引物设计的超快工具(Du et al., 2017)。该软件是在2013年本实验室开发的微卫星搜索及统计软件MSDB(Du et al., 2013)基础之上进行开发的,有着更强大的功能。Krait可以很容易地用于筛选微卫星标记,用户可以高效地设计用于微卫星放大的引物。李午佼等(2014)曾使用MSDB分析了大熊猫基因组微卫星的分布特征。研究表明,Krait比MSDB搜索检测速度更快、更高效,能够快速研究DNA序列中的微卫星(Du et al., 2017)。对比2个软件的分析结果,发现Krait不仅能快速筛选微卫星标记,还能有效地检测到更多微卫星重复序列,在种类及丰度等分布特征的识别上更为准确。因此本研究使用Krait对大熊猫基因组中微卫星序列进行检测的结果是可信、可靠的。
到目前为止,也有较多的微卫星标记成功用于野生大熊猫的种群遗传多样性分析(He et al., 2008;Hu et al., 2010)。据统计,目前公布的大熊猫的微卫星标记共109种,其中超过60%是二碱基微卫星(乔麦菊等,2019)。根据目前研究进展,用于野生大熊猫遗传多样性调查研究的微卫星标记多,且存在差异,故无法进行统一的比较分析,未来关于大熊猫微卫星标记的开发在稳定实用的基础上,应尽量统一标准。本研究中使用6只大熊猫的转录组数据与大熊猫基因组相结合,利用软件筛选多态性微卫星,共筛选出了326个多态性微卫星,处理过程耗时短,减少了很多工作量。而从中挑选的20个微卫星位点通过实验验证,显示部分位点多态性较高,证明该方法是可行的,可以对后续位点筛选提供资源。
群体的遗传多样性主要表现在等位基因数、杂合度和多态信息含量3个方面(Senanan et al., 2004)。本研究选取的20个位点中,高度多态性位点有9个,中度多态性位点有11个,其中,有4个显著偏离哈迪-温伯格平衡(P < 0.01)。后续实验可将这20个位点中多态性较高的位点选用粪便DNA进行扩增和其他参数的评估,筛选出能够应用于野生大熊猫遗传多样性研究的微卫星标记,为野生大熊猫种群的群体结构及遗传研究提供一定的理论指导。
根据分析结果,本研究认为NGS是一种有效开发大熊猫微卫星标记的方法。利用大熊猫的NGS序列数据可以开发大量的微卫星标记。微卫星分子标记的开发将有助于评估大熊猫的遗传多样性和种群结构,并制定有效的保护和管理策略。
李薇, 李久煊, 荆慧芳, 等. 2017. 基于高通量测序的达氏鲟微卫星标记筛选[J]. 动物学杂志, 52(3): 449-457. |
李午佼, 李玉芝, 杜联明, 等. 2014. 大熊猫和北极熊基因组微卫星分布特征比较分析[J]. 四川动物, 33(6): 874-878. |
乔麦菊, 冉江洪, 张和民, 等. 2019. 微卫星标记在大熊猫研究中的应用进展[J]. 兽类学报, 39(1): 103-110. |
青菁. 2016.四川省大熊猫栖息地破碎化现状研究及廊道规划[D].南充: 西华师范大学. http://cdmd.cnki.com.cn/Article/CDMD-10638-1016748842.htm
|
宋琪, 刘金龙, 郭宪光, 等. 2019. 基于Roche 454 GS FLX高通量测序的叶城沙蜥基因组微卫星特征分析[J]. 四川动物, 38(1): 62-67. |
唐小平, 贾建生, 王志臣, 等. 2015. 全国第四次大熊猫调查方案设计及主要结果分析[J]. 林业资源管理, (1): 11-16. |
修云芳, 李碧春, 陈玉村, 等. 2015. 应用微卫星分型进行小熊猫亲子鉴定(英文)[J]. 兽类学报, 35(1): 55-64. |
岳华梅, 翟晴, 宋明月, 等. 2016. 基于转录组测序的兴国红鲤微卫星标记筛选[J]. 淡水渔业, 46(1): 24-28. |
张正义, 邢秀梅, 胡鹏飞, 等. 2017. 微卫星标记在动物遗传多样性分析的研究进展[J]. 经济动物学报, 21(3): 164-168. |
Botstein D, White RL, Skolnick M, et al. 1980. Construction of a genetic linkage map in man using restriction fragment length polymorphisms[J]. American Journal of Human Genetics, 32: 314-324. |
Du L, Zhang C, Liu Q, et al. 2017. Krait:an ultrafast tool for genome-wide survey of microsatellites and primer design[J]. Bioinformatics, 34(4): 681-683. |
Du LM, Li YZ, Zhang XY, et al. 2013. MSDB:a user-friendly program for reporting distribution and building databases of microsatellites from genome sequences[J]. Journal of Heredity, 104(1): 154-157. |
Haas BJ, Papanicolaou A, Yassour M, et al. 2013. De novo transcript sequence reconstruction from RNA-seq using the Trinity platform for reference generation and analysis[J]. Nature Protocols, 8(8): 1494-1512. |
He W, Lin L, Shen FJ, et al. 2008. Genetic diversities of the giant panda (Ailuropoda melanoleuca) in Wanglang and Baoxing Nature Reserve[J]. Conservation Genetics, 9: 1541-1546. |
Hu Y, Zhan X, Qi D, et al. 2010. Spatial genetic structure and dispersal of giant pandas on a mountain-range scale[J]. Conservation Genetics, 11(6): 2145-2155. |
Huang J, Li YZ, Du LM, et al. 2015. Genome-wide survey and analysis of microsatellites in giant panda (Ailuropoda melanoleuca), with a focus on the applications of a novel microsatellite marker system[J]. BMC Genomics, 16(1): 61. DOI:10.1186/s12864-015-1268-z |
Jurka J, Pethiyagoda C. 1995. Simple repetitive DNA sequences from primates:compilation and analysis[J]. Journal of Molecular Evolution, 40(2): 120-126. |
Li WZ, Godzik A. 2006. Cd-hit:a fast program for clustering and comparing large sets of protein or nucleotide sequences[J]. Bioinformatics, 22(13): 1658-1659. |
Li YZ, Xu X, Shen FJ, et al. 2010. Development of new tetranucleotide microsatellite loci and assessment of genetic variation of giant panda in two largest giant panda captive breeding populations[J]. Journal of Zoology (London), 282(1): 39-46. |
Liu HG, Yang Z, Tang HY, et al. 2017. Microsatellite development and characterization for Saurogobio dabryi Bleeker, 1871 in a Yangtze river-connected lake, China[J]. Journal of Genetics, 96: 1-4. |
Lu Z, Johnson WE, Menotti-Raymond M, et al. 2001. Patterns of genetic diversity in remaining giant panda populations[J]. Conservation Biology, 15(6): 1596-1607. |
Ma HY, Yue YS, Lu Y. 2004. Microsatellite DNA molecular marker and its application in animal breeding and genetics[J]. Journal of Animal Science & Veterinary Medicine, 23(5): 16-19. |
Marshall TC, Slate JB, Kruuk LE, et al. 1998. Statistical confidence for likelihood-based paternity inference in natural populations[J]. Molecular Ecology, 7(5): 639-655. |
Martin M. 2015. Cutadapt removes adapter sequences from high-throughput sequencing reads[J/OL]. EMBnet.jour-nal, 17(1): 10-12[2019-05-05]. https://journal.embnet.org/index.php/embnetjournal/article/view/200/479. DOI: 10.14806/ej.17.1.200.
|
Nagy S. 2012. PICcalc:an online program to calculate polymorphic information content for molecular genetic studies[J]. Biochemical Genetics, 50(9-10): 670-672. |
Qi WH, Jiang XM, Du LM, et al. 2015. Genome-wide survey and analysis of microsatellite sequences in bovid species[J/OL]. PLoS ONE, 10: e0133667[2019-05-05]. https://doi.org/10.1371/journal.pone.0133667.
|
Raymond M, Rousset F. 1995. GENEPOP (version 1.2):population genetics software for exact tests and ecumenicism[J]. Heredity, 86(3): 248-249. |
Ritchie H, Jamieson AJ, Piertney SB. 2016. Isolation and characterisation of microsatellite DNA markers in the deep-sea amphipod Paralicella tenuipes by Illumina Miseq sequencing[J]. Journal of Heredity, 107(4): 367-371. |
Selkoe KA, Toonen RJ. 2006. Microsatellites for ecologists:a practical guide to using and evaluating microsatellite markers[J]. Ecology Letters, 9(5): 615-629. |
Senanan W, Kapuscunski AR, Na-Nakorn U, et al. 2004. Genetic impacts of hybrid catfish faming (Clarias macrocephalus×C. gariepinus) on native catfish populations in central Thailand[J]. Aquaculture, 235(1-4): 167-184. |
Shen FJ, Watts P, He ZZ, et al. 2005. Enrichment of giant panda microsatellite markers using dynal magnet beads[J]. Acta Genetica Sinica, 32(5): 457-462. |
Shete S, Tiwari H, Elston RC. 2000. On estimating the heterozygosity and polymorphic information content value[J]. Theoretical Population Biology, 57(3): 265-271. |
Tian CX, Liang XF, Yang M, et al. 2014. New microsatellite loci for the mandarin fish Siniperca chuatsi and their application in population genetic analysis[J]. Genetics and Molecular Research, 13(1): 546-558. |
Wang D, Hu Y, Ma T, et al. 2016. Noninvasive genetics provides insights into the population size and genetic diversity of an Amur tiger population in China[J]. Integrative Zoology, 11(1): 16-24. |
Wei F, Hu Y, Zhu L, et al. 2012. Black and white and read all over:the past, present and future of giant panda genetics[J]. Molecular Ecology, 21(23): 5660-5674. |
Wu H, Zhan XJ, Zhang ZJ, et al. 2008. Thirty-three microsatellite loci for noninvasive genetic studies of the giant panda (Ailuropoda melanoleuca)[J]. Conservation Genetics, 10(3): 649-652. |
Xia EH, Yao QY, Zhang HB, et al. 2016. CandiSSR: an efficient pipeline used for identifying candidate polymorphic SSRs based on multiple assembled sequences[J/OL]. Frontiers in Plant Science, 6: 1171[2019-05-05]. https://doi.org/10.3389/fpls.2015.01171.
|
Xiao TQ, Lu CY, Xu YL, et al. 2015. Screening of SSR markers as sociated with scale cover pattern and mapped to a genetic linkage map of common carp (Cyprinus carpio L.)[J]. Journal of Applied Genetics, 56(2): 261-269. |
Xu Y, Hu Z, Wang C, et al. 2016. Characterization of perfect microsatellite based on genome-wide and chromosome level in rhesus monkey (Macaca mulatta)[J]. Gene, 592: 269-275. |
Zheng XF, Pan C, Diao Y, et al. 2013. Development of microsatellite markers by transcriptome sequencing in two species of Amorphophallus (Araceae)[J/OL]. BMC Genomics, 14(1): 490[2019-05-05]. https://bmcgenomics.biomedcentral.com/articles/10.1186/1471-2164-14-490.
|