Genetic diversity and population genetic structure of Python bivittatus in China

Abstract

The Burmese python (Python bivittatus) has recently suffered large population declines in the wild in China due to illegal capture, overexploitation, deforestation and the loss of its natural habitat. Greater knowledge of the genetic diversity and structure of wild P. bivittatus populations is needed to help ensure its effective management. In this study, we investigated the genetic diversity and population genetic structure of wild P. bivittatus in China in detail. 109 P. bivittatus individuals from five distribution areas in Guangdong (GD), Guangxi (GX), Hainan (HN), Fujian (FJ) and Yunnan (YN) province of China were collected, and their genetic structure and diversity were analyzed. Eight highly polymorphic microsatellite loci were utilized to reveal high levels of genetic diversity in the P. bivittatus population. Genetic diversity was highest in GX, and lowest in GD. All geographic populations demonstrated a bottleneck effect indicating recent population decline. F _st and AMOVA analyses revealed that there was moderate genetic differentiation among the five populations, and that only 10.59 % of total genetic diversity occurred among populations. F _st values between populations were positively correlated with their geographical distances. Genetic structure analyses revealed that the HN, GX and GD populations, which were geographically closest, were assigned to a genetic cluster, while the YN and FJ populations constituted a single cluster, respectively.

Keywords

Python bivittatus / Microsatellite / Genetic diversity / Population genetic

References

[1] Astarita G, Rourke BC, Andersen JB, Fu J, Kim JH, Bennett AF, Hicks JW, Piomelli D. Postprandial increase of oleoylethanolamide mobilization in small intestine of the Burmese python (Python molurus). Am J Phys, 2006, 290: 1407-1412.

[2] Barker DG, Barker TM. The distribution of the Burmese Python, Python molurus bivittatus. Bull Chicago Herp Soc, 2008, 43: 33-38.

[3] Barker DG, Barker TM. The distribution of the Burmese python, Python bivittatus. Bull Chicago Herp Soc, 2010, 45: 86-88.

[4] Botstein D, White RL, Skalnick MH, Davis RW. Construction of a genetic linkage map in man using restriction fragment lengthpolymorphism. Am J Hum Genet, 1980, 32: 314-331.

[5] Castoe TA, deKoning AP, Hall KT, Yokoyama KD, Gu W, Smith EN, Feschotte C, Uetz P, Ray DA, Dobry J, Bogden R, Mackessy SP, Bronikowski AM, Warren WC, Secor SM, Pollock DD. Sequencing the genome of the Burmese python (Python molurus bivittatus) as a model for studying extreme adaptations in snakes. Genome Biol, 2011, 12: 406.

[6] Castoe TA, Poole AW, deKoning AP, Jones KL, Tomback DF, Oyler-McCance SJ, Fike JA, Lance SL, Streicher JW, Smith EN, Pollock DD. Rapid microsatellite identification from illumina paired-end genomic sequencing in two birds and a snake. PLoS One, 2012, 7: e30953.

[7] Coatoe TA, Streicher JW, Meik JM, Ingrasci MJ, Poole AW, deKoning AP, Campbell JA, Parkinson CL, Smith EN, Pollock DD. Thousands of microsatellite loci from the venomous coralsnake Micrurus fulvius and variability of select loci across populations and related species. Mol Ecol Resour, 2012, 12: 1105-1113.

[8] Collins TM, Freeman B. Genetic characterization of populations of the nonindigenous Burmese python in Everglades National Park, 2008, Miami: Florida International University 6 8

[9] Cornuet JM, Luikart G. Description and power analysis of two tests for detecting recent population bottlenecks from allele frequency data. Genet, 1996, 144: 2001-2014.

[10] Dorcas ME, Willson JD, Gibbons JW. Can invasive Burmese pythons inhabit temperate regions of the southeastern United States?. Biol Invasions, 2011, 13: 793-802.

[11] Dorcas ME, Willson JD, Reed RN, Snow RW, Rochford MR, Miller MA, Meshaka WE, Andreadis PT, Mazzotti FJ, Romagosa CM, Hart KM. Severe mammal declines coincide with proliferation of invasive Burmese pythons in Everglades National Park. Proc Natl Acad Sci USA, 2012, 109: 2418-2422.

[12] Evanno G, Regnaut S, Goudet J. Detecting the number of clusters of individuals using the software structure: a simulation study. Mol Ecol, 2005, 14: 2611-2620.

[13] Excoffier L, Lischer HE. Arlequin suite ver3. 5: a new series of programs to perform population genetics analyses under Linux and Windows. Mol Ecol Resour, 2010, 10: 564-567.

[14] Fontanella FM, Siddall ME. Isolation and characterization of 14 polymorphic microsatellite loci in the ringneck snake Diadophis punctatus (Colubridae: Dipsadinae). Conserv Genet, 2010, 11: 1193-1195.

[15] Frank JM, Michael SC, Kristen MH, Ray WS, Michael RR, Michael ED. Cold-induced mortality of invasive Burmese pythons in south Florida. Biol Invasions, 2011, 13: 143-151.

[16] Frankham R, Ballou JD, Briscoe DA. Introduction to conservation genetics, 2002, Oxford: Cambridge University 78 83

[17] Gibbs HL, Prior KA, Weatherhead PJ, Johnson G. Genetic structure of populations of the threatened eastern massasuaga rattlesnake, Sistrurus catenatus: evidence from microsatellite DNA markers. Mol Ecol, 1997, 6: 1123-1132.

[18] Groot TV, Bruins E, Breeuwer JA. Molecular genetic evidence for parthenogenesis in the Burmese python, Python molurus bivittatus. Heredity, 2003, 90: 130-135.

[19] Gong S, Shi H, Mo Y Phylogeography of the endangered black-breasted leaf turtle (Geoemyda spengleri) and conservation implications for other chelonians. Amphibia-Reptilia, 2009, 30: 57-62.

[20] Hunter ME, Hart KM. Rapid microsatellite marker development using next generation pyrosequencing to inform invasive Burmese python—Python molurus bivittatus—Management. Int J Mol Sci, 2013, 14: 4793-4804.

[21] Jensen JL, Bohonak AJ, Kelley ST (2005) Isolation by distance, web service. BMC Genetics 6:13. v.3.23. http://ibdws.sdsu.edu/

[22] Jordan PW, Goodman AE, Donnellan S. Microsatellite primers for Australian and New Guinean pythons isolated with an efficient marker development method for related species. Mol Ecol Notes, 2002, 2: 78-82.

[23] Kalinowski T, Taper ML, Marshall TC. Revising how the computer program cervus accommodates genotyping error increases success in paternity assignment. Mol Ecol, 2007, 16: 1099-1106.

[24] Ke WC, Wu LL, Wu WS, Chen YL. Identify of Python molurus based on Cyt b and CO I gene sequence. J Fujian Norm Univ, 2014, 30: 91-97.

[25] Li PP, Wang WS, Lü XP. Snake conservation and sustainable utilization in China: history, statues and future. J Shenyang Norm Univ, 2014, 31: 129-135.

[26] Lin LH, Ji X, Diong CH, Du Y, Lin CX. Phylogeography and population structure of the Reevese’s Butterfly Lizard (Leiolepis reevesii) inferred from mitochondrial DNA sequences. Mol Phylogenet Evol, 2010, 56: 601-607.

[27] Liu LJ, Song ZJ, Cao SY, Zhai YF, Wang AJ, Gao L, Xiu QY, Cao DH. The dynamic setting and processes of tectonic and magmatic evolution of the oblique collision zone between Indian and Eurasion plates: Exemplified by the tectonic evolution of the Three River region, eastern Tibet. Acta Petrol Sin, 2006, 22: 775-786.

[28] Liu P, Chen H, Zhao WG, Jia JB. Snakes species diversity and its conservation in China. Chin J Wildl, 2009, 30: 44-46.

[29] Liu YF, Cao T, Liu JC, Yu P, Zhao CP, Hou GY, Zhang LL. Complete mitochondrial genome sequence and analysis of the Burmese python Python bivittatus. Sichuan J Zool, 2014, 33: 23-30.

[30] Manjaka A, Carolyn AB, Gary DS, Olga R, Rick AB, Edward EL. Characterization of 20 microsatellite marker loci in the Malagasy tree boa (Sanzinia madagascariensis madagascariensis). Conserv Genet, 2009, 10: 1953-1956.

[31] Mao LQ. Phylogeography and population genetic structure of the Chinese Cobra, Naja Atra, 2011, Nanjing: Nanjing Noamal University 11 13

[32] Marshall JC, Kingsbury BA, Minchella DJ. Microsatellite variation, population structure, and bottlenecks in the threatened copperbelly water snake. Conserv Genet, 2009, 10: 465-476.

[33] Maruyama T, Fuerst PA. Population bottlenecks and nonequilibrium models in population genetics. II. Number of alleles in a small population that was formed by a recent bottleneck. Genetics, 1985, 111: 675-689.

[34] Mauldin RE, Savarie PJ. Acetaminophen as an Oral Toxicant for Nile Monitor Lizards (Varanus niloticus) and Burmese pythons (Python molurus bivittatus). Wildl Res, 2010, 37: 215-222.

[35] Mccue MD, Bennett AF, Hicks JW. The effect of meal composition on specific dynamic action in Burmese pythons (Python molurus). Physiol Biochem Zool, 2005, 78: 182-192.

[36] Meshaka WE, Loftus WF, Steiner T. The herpetofauna of Everglades National Park. Fla Sci, 2000, 63: 84-103.

[37] Michael LA, Richard ME, Kandy LK, John SH, William EB, Tom CM, Richard EM. Cold weather and the potential range of invasive Burmese pythons. Biol Invasions, 2010, 12: 3649-3652.

[38] Nei M. Estimation of average heterozygosity and genetic distance from a small number of individuals. Genetics, 1978, 89: 583-590.

[39] Nei M, Maruyama T, Chakraborty R. The bottleneck effect and genetic variability in populations. Evolution, 1975, 29: 1-10.

[40] Paquin MM, Wylie GD, Routman EJ. Population structure of the giant garter snake, Thamnophis gigas. Conserv Genet, 2006, 7: 25-36.

[41] Pritchard JK, Stephens M, Donnelly P. Inference of population structure using multilocus genotype data. Genetics, 2000, 155: 945-959.

[42] Ratnoff OD, Rosenberg MJ, Everson B, Emanuelson M, Tulodziecki N. Notes on clotting in a Burmese python (Python molurus bivittatus). J Lab Clin Med, 1990, 115: 629-635.

[43] Saghai-Maroof MA, Soliman KM, Jorgensen RA, Allard RW. Ribosomal DNA spacer length in barley: Mendelian inheritance, chromosomal location, and population dynamics. Proc Natl Acad Sci USA, 1984, 81: 8014-8018.

[44] State Forestry Administration. The survey focused on key terrestrial wildlife in China, 2009, Beijing: Forestry 74 88

[45] Stuart B, Nguyen TQ, Thy N, Grismer L, Chan-Ard T, Iskandar D, Golynsky E, Lau MWN (2014) Python bivittatus. The IUCN Red List of Threatened Species. www.iucnredlist.org. Accessed on 23 Oct 2014

[46] Tzika AC, Koenig S, Miller R, Garcia G, Remy C, Milinkovitch MC. Population structure of an endemic vulnerable species, the Jamaican boa (Epicrates subflavus). Mol Ecol, 2008, 17: 533-544.

[47] Wang S, Xie Y. China red data book of endangered animals (amphibia and reptilia), 2004, Beijing: Science 208 211

[48] Wang T, Altimiras J, Axelsson M. Intracardiac flow separation in an in situperfused heart from Burmese python (Python molurus). J Exp Biol, 2002, 205: 2715-2723.

[49] Willson JD, Dorcas ME, Snow RW. Identifying plausible scenarios for the establishment of invasive Burmese pythons (Python molurus) in Southern Florida. Biol Invasions, 2011, 13: 1493-1504.

[50] Wu W, Chen ZW, Lu XF, Fu XJ, Long LN, Tan JM, Su ZR. Medicinal value of the Burmese python. J Snake, 2010, 22: 145-147.

[51] Yan J. The study on evolutional history of Hainan island’s ecological environment, 2006, Nanjing: Nanjing Agricultural University 40 41

[52] Zhang JJ, Zhong DL, Sang HQ, Zhou Y. Structural and geochronological evidence for multiple episodes of deformation since Paleocene along the Ailao Shan-Red River shear zone, Southeastern Asia. Chin J Geol, 2006, 41: 291-310.

[53] Zhang YZ. Problems and countermeasures on Snakes resources protection in China management. J Snake, 2010, 22: 343-348.