Phylogeography, historical demography and systematics of the world's smallest pythons (Pythonidae, Antaresia).

Advances from empirical studies in phylogeography, systematics and species delimitation highlight the importance of integrative approaches for quantifying taxonomic diversity. Genomic data have greatly improved our ability to discern both systematic diversity and evolutionary history. Here we combine analyses of mitochondrial DNA sequences, thousands of genome-wide SNPs and linear and geometric morphometrics on Antaresia, a clade of four currently recognised dwarf pythons from Australia and New Guinea (Antaresia childreni, A. stimsoni, A. maculosa and A. perthensis). Our integrative analyses of phylogenetics, population structure, species delimitation, historical demography and morphometrics revealed that the true evolutionary diversity is not well reflected in the current appraisal of the diversity of the group. We find that Antaresia childreni and A. stimsoni comprise a widespread network of populations connected by gene flow and without evidence of species-level divergence among them. However, A. maculosa shows considerable genetic structuring which leads us to recognise two subspecies in northeastern Australia and a new species in Torres Strait and New Guinea. These two contrasting cases of over and under estimation of diversity, respectively, illustrate the power of thorough integrative approaches into understanding evolution of biodiversity. Furthermore, our analyses of historical demographic patterns highlight the importance of the Kimberley, Pilbara and Cape York as origins of biodiversity in Australia.

Phylogenomics, biogeography and taxonomic revision of New Guinean pythons (Pythonidae, Leiopython) harvested for international trade.

The large and enigmatic New Guinean pythons in the genus Leiopython are harvested from the wild to supply the international trade in pets. Six species are currently recognized (albertisii, biakensis, fredparkeri, huonensis, meridionalis, montanus) but the taxonomy of this group has been controversial. We combined analysis of 421 nuclear loci and complete mitochondrial genomes with morphological data to construct a detailed phylogeny of this group, understand their biogeographic patterns and establish the systematic diversity of this genus. Our molecular genetic data support two major clades, corresponding to L. albertisii and L. fredparkeri, but offer no support for the other four species. Our morphological data also only support two species. We therefore recognize L. albertisii and L. fredparkeri as valid species and place L. biakensis, L. meridionalis, L. huonensis and L. montanus into synonymy. We found that L. albertisii and L. fredparkeri are sympatric in western New Guinea; an atypical pattern compared to other Papuan species complexes in which the distributions of sister taxa are partitioned to the north and south of the island's central mountain range. For the purpose of conservation management, overestimation of species diversity within Leiopython has resulted in the unnecessary allocation of resources that could have been expended elsewhere. We strongly caution against revising the taxonomy of geographically widespread species groups when little or no molecular genetic data and only small morphological samples are available.

Phylogenomics, Biogeography, and Morphometrics Reveal Rapid Phenotypic Evolution in Pythons After Crossing Wallace's Line.

Ecological opportunities can be provided to organisms that cross stringent biogeographic barriers towards environments with new ecological niches. Wallace's and Lyddeker's lines are arguably the most famous biogeographic barriers, separating the Asian and Australo-Papuan biotas. One of the most ecomorphologically diverse groups of reptiles, the pythons, is distributed across these lines, and are remarkably more diverse in phenotype and ecology east of Lydekker's line in Australo-Papua. We used an anchored hybrid enrichment approach, with near complete taxon sampling, to extract mitochondrial genomes and 376 nuclear loci to resolve and date their phylogenetic history. Biogeographic reconstruction demonstrates that they originated in Asia around 38-45 Ma and then invaded Australo-Papua around 23 Ma. Australo-Papuan pythons display a sizeable expansion in morphological space, with shifts towards numerous new adaptive optima in head and body shape, coupled with the evolution of new micro-habitat preferences. We provide an updated taxonomy of pythons and our study also demonstrates how ecological opportunity following colonization of novel environments can promote morphological diversification in a formerly ecomorphologically conservative group. [Adaptive radiation; anchored hybrid enrichment; biogeography; morphometrics; snakes.].

Species delimitation and systematics of the green pythons (Morelia viridis complex) of melanesia and Australia.

Molecular data sets and the increasing use of integrative systematics is revealing cryptic diversity in a range of taxa - particularly in remote and poorly sampled landscapes like the island of New Guinea. Green pythons (Morelia viridis complex) are one of the most conspicuous elements of this island's fauna, with large numbers taken from the wild to supply international demand for exotic pets. We test hypotheses about species boundaries in green pythons from across New Guinea and Australia with mitochondrial genomes, 389 nuclear exons, and comprehensive assessment of morphological variation. Strong genetic structuring of green python populations and species delimitation methods confirm the presence of two species, broadly occurring north and south of New Guinea's central mountains. Our data also support three subspecies within the northern species. Subtle but consistent morphological divergence among the putative taxa is concordant with patterns of molecular divergence. Our extensive sampling identifies several zones of hitherto unknown biogeographical significance on the island of New Guinea. We revise the taxonomy of the group, discuss the relevance of our findings in the context of Papuan biogeography and the implications of our systematic changes for the conservation management of these taxa.

Parallel selective pressures drive convergent diversification of phenotypes in pythons and boas.

Pythons and boas are globally distributed and distantly related radiations with remarkable phenotypic and ecological diversity. We tested whether pythons, boas and their relatives have evolved convergent phenotypes when they display similar ecology. We collected geometric morphometric data on head shape for 1073 specimens representing over 80% of species. We show that these two groups display strong and widespread convergence when they occupy equivalent ecological niches and that the history of phenotypic evolution strongly matches the history of ecological diversification, suggesting that both processes are strongly coupled. These results are consistent with replicated adaptive radiation in both groups. We argue that strong selective pressures related to habitat-use have driven this convergence. Pythons and boas provide a new model system for the study of macro-evolutionary patterns of morphological and ecological evolution and they do so at a deeper level of divergence and global scale than any well-established adaptive radiation model systems.

Phylogenomic analyses reveal novel relationships among snake families.

Snakes are a diverse and important group of vertebrates. However, relationships among the major groups of snakes have remained highly uncertain, with recent studies hypothesizing very different (and typically weakly supported) relationships. Here, we address family-level snake relationships with new phylogenomic data from 3776 nuclear loci from ultraconserved elements (1.40million aligned base pairs, 52% missing data overall) sampled from 29 snake species that together represent almost all families, a dataset ∼100 times larger than used in previous studies. We found relatively strong support from species-tree analyses (NJst) for most relationships, including three largely novel clades: (1) a clade uniting the boas, pythons and their relatives, (2) a clade placing cylindrophiids and uropeltids with this clade, and (3) a clade uniting bolyeriids (Round Island boas) with pythonids and their relatives (xenopeltids and loxocemids). Relationships among families of advanced snakes (caenophidians) were also strongly supported. The results show the potential for phylogenomic analyses to resolve difficult groups, but also show a surprising sensitivity of the analyses to the inclusion or exclusion of outgroups.

Phylogeographic and population genetic analyses reveal multiple species of Boa and independent origins of insular dwarfism.

Boa is a Neotropical genus of snakes historically recognized as monotypic despite its expansive distribution. The distinct morphological traits and color patterns exhibited by these snakes, together with the wide diversity of ecosystems they inhabit, collectively suggest that the genus may represent multiple species. Morphological variation within Boa also includes instances of dwarfism observed in multiple offshore island populations. Despite this substantial diversity, the systematics of the genus Boa has received little attention until very recently. In this study we examined the genetic structure and phylogenetic relationships of Boa populations using mitochondrial sequences and genome-wide SNP data obtained from RADseq. We analyzed these data at multiple geographic scales using a combination of phylogenetic inference (including coalescent-based species delimitation) and population genetic analyses. We identified extensive population structure across the range of the genus Boa and multiple lines of evidence for three widely-distributed clades roughly corresponding with the three primary land masses of the Western Hemisphere. We also find both mitochondrial and nuclear support for independent origins and parallel evolution of dwarfism on offshore island clusters in Belize and Cayos Cochinos Menor, Honduras.

Genetic divergence and diversity in the Mona and Virgin Islands Boas, Chilabothrus monensis (Epicrates monensis) (Serpentes: Boidae), West Indian snakes of special conservation concern.

Habitat fragmentation reduces the extent and connectivity of suitable habitats, and can lead to changes in population genetic structure. Limited gene flow among isolated demes can result in increased genetic divergence among populations, and decreased genetic diversity within demes. We assessed patterns of genetic variation in the Caribbean boa Chilabothrus monensis (Epicrates monensis) using two mitochondrial and seven nuclear markers, and relying on the largest number of specimens of these snakes examined to date. Two disjunct subspecies of C. monensis are recognized: the threatened C. m. monensis, endemic to Mona Island, and the rare and endangered C. m. granti, which occurs on various islands of the Puerto Rican Bank. Mitochondrial and nuclear markers revealed unambiguous genetic differences between the taxa, and coalescent species delimitation methods indicated that these snakes likely are different evolutionary lineages, which we recognize at the species level, C. monensis and C. granti. All examined loci in C. monensis (sensu stricto) are monomorphic, which may indicate a recent bottleneck event. Each population of C. granti exclusively contains private mtDNA haplotypes, but five of the seven nuclear genes assayed are monomorphic, and nucleotide diversity is low in the two remaining markers. The faster pace of evolution of mtDNA possibly reflects the present-day isolation of populations of C. granti, whereas the slower substitution rate of nuDNA may instead mirror the relatively recent episodes of connectivity among the populations facilitated by the lower sea level during the Pleistocene. The small degree of overall genetic variation in C. granti suggests that demes of this snake could be managed as a single unit, a practice that would significantly increase their effective population size.

Toward a Tree-of-Life for the boas and pythons: multilocus species-level phylogeny with unprecedented taxon sampling.

Snakes in the families Boidae and Pythonidae constitute some of the most spectacular reptiles and comprise an enormous diversity of morphology, behavior, and ecology. While many species of boas and pythons are familiar, taxonomy and evolutionary relationships within these families remain contentious and fluid. A major effort in evolutionary and conservation biology is to assemble a comprehensive Tree-of-Life, or a macro-scale phylogenetic hypothesis, for all known life on Earth. No previously published study has produced a species-level molecular phylogeny for more than 61% of boa species or 65% of python species. Using both novel and previously published sequence data, we have produced a species-level phylogeny for 84.5% of boid species and 82.5% of pythonid species, contextualized within a larger phylogeny of henophidian snakes. We obtained new sequence data for three boid, one pythonid, and two tropidophiid taxa which have never previously been included in a molecular study, in addition to generating novel sequences for seven genes across an additional 12 taxa. We compiled an 11-gene dataset for 127 taxa, consisting of the mitochondrial genes CYTB, 12S, and 16S, and the nuclear genes bdnf, bmp2, c-mos, gpr35, rag1, ntf3, odc, and slc30a1, totaling up to 7561 base pairs per taxon. We analyzed this dataset using both maximum likelihood and Bayesian inference and recovered a well-supported phylogeny for these species. We found significant evidence of discordance between taxonomy and evolutionary relationships in the genera Tropidophis, Morelia, Liasis, and Leiopython, and we found support for elevating two previously suggested boid species. We suggest a revised taxonomy for the boas (13 genera, 58 species) and pythons (8 genera, 40 species), review relationships between our study and the many other molecular phylogenetic studies of henophidian snakes, and present a taxonomic database and alignment which may be easily used and built upon by other researchers.

Species identification of protected carpet pythons suitable for degraded forensic samples.

In this paper we report on the identification of a section of mitochondrial DNA that can be used to identify the species of protected and illegally traded pythons of the genus Morelia. Successful enforcement of wildlife laws requires forensic tests that can identify the species nominated in the relevant legislation. The potentially degraded state of evidentiary samples requires that forensic investigation using molecular genetic species identification is optimized to interrogate small fragments of DNA. DNA was isolated from 35 samples of Morelia spilota from which the complete cytochrome b was sequenced. The ND6 gene was also sequenced in 32 of these samples. Additional DNA sequences were generated from 9 additional species of Morelia. The sequences were aligned by Geneious and imported into MEGA to create phylogenetic trees based on the entire complex of approximately 1,706 base pairs (bp). To mimic degraded DNA, which is usually found in forensic cases, short sub-sections of the full alignment were used to generate phylogenetic trees. The sub-sections that had the greatest DNA sequence information were in parts of the cytochrome b gene. Our results highlight that legislation is presently informed by inadequate taxonomy. We demonstrated that a 278 bp region of the cytochrome b gene recovered the topology of the phylogenetic tree found with the entire gene sequence and correctly identified species of Morelia with a high degree of confidence. The locus described in this report will assist in the successful prosecution of alleged illegal trade in python species.

Differentiation between wild type and heterozygous albino ball pythons (Python regius) by PCR and qPCR.

Python regius or ball pythons are the famous exotic pets because of their beautiful color and pattern. The albino ball python is one type of ball python, but it is very difficult to determine the difference of phenotype between wildtype and heterozygous genotype of albino (het albino). In this study, PCR and qPCR can distinguish between wildtype and het albino. The PCR product size of wildtype and het albino was 415 bp, but the intensity of PCR product of wildtype was more intense than that of het albinos. No PCR amplicon was found in albinos and the Ct value of wildtype was lower than Ct of het albinos. The molecular detection technique, especially PCR and qPCR, can determine the difference between wildtype and het albinos of ball pythons.

Molecular phylogeny and historical biogeography of West Indian boid snakes (Chilabothrus).

The evolutionary and biogeographic history of West Indian boid snakes (Epicrates), a group of nine species and 14 subspecies, was once thought to be well understood; however, new research has indicated that we are missing a clear understanding of the evolutionary relationships of this group. Here, we present the first multilocus, species-tree based analyses of the evolutionary relationships, divergence times, and historical biogeography of this clade with data from 10 genes and 6256 bp. We find evidence for a single colonization of the Caribbean from mainland South America in the Oligocene or early Miocene, followed by a radiation throughout the Greater Antilles and Bahamas. These findings support the previous suggestion that Epicrates sensu lato Wagler is paraphyletic with respect to the anacondas (Eunectes Wagler), and hence we restrict Epicrates to the mainland clade and use the available name Chilabothrus Duméril and Bibron for the West Indian clade. Our results suggest some diversification occurred within island banks, though most species divergence events seem to have occurred in allopatry. We also find evidence for a remarkable diversification within the Bahamian archipelago suggesting that the recognition of another Bahamian endemic species C. strigilatus is warranted.

Molecular systematics and historical biogeography of tree boas (Corallus spp.).

Inferring the evolutionary and biogeographic history of taxa occurring in a particular region is one way to determine the processes by which the biodiversity of that region originated. Tree boas of the genus Corallus are an ancient clade and occur throughout Central and South America and the Lesser Antilles, making it an excellent group for investigating Neotropical biogeography. Using sequenced portions of two mitochondrial and three nuclear loci for individuals of all recognized species of Corallus, we infer phylogenetic relationships, present the first molecular analysis of the phylogenetic placement of the enigmatic C. cropanii, develop a time-calibrated phylogeny, and explore the biogeographic history of the genus. We found that Corallus diversified within mainland South America, via over-water dispersals to the Lesser Antilles and Central America, and via the traditionally recognized Panamanian land bridge. Divergence time estimates reject the South American Caribbean-Track as a general biogeographic model for Corallus and implicate a role for events during the Oligocene and Miocene in diversification such as marine incursions and the uplift of the Andes. Our findings also suggest that recognition of the island endemic species, C. grenadensis and C. cookii, is questionable as they are nested within the widely distributed species, C. hortulanus. Our results highlight the importance of using widespread taxa when forming and testing biogeographic hypotheses in complex regions and further illustrate the difficulty of forming broadly applicable hypotheses regarding patterns of diversification in the Neotropical region.

Chemosensory age discrimination in the snake Boa constrictor (Serpentes: Boidae).

Many snakes are able to use their chemosensory system to detect scent of conspecifics, which is important in many social contexts. Age discrimination based on chemical cues may be especially important to ensure access to sexually mature potential partners. In this study, we used 24 individual Boa constrictor snakes (12 adults mature and 12 non-mature individuals) that had been captured in different areas of Ecuador, and were maintained in captivity at the Vivarium of Quito. We used tongue-flick experiments to examine whether these snakes were able to discriminate between scents from mature and non-mature individuals. Results showed that B. constrictor snakes used chemical cues to recognize conspecifics and that the scent of individuals of different ages elicited chemosensory responses of different magnitudes. The scents from adult conspecifics elicited the quickest and highest chemosensory responses (i.e., short latency times and high tongue-flick rates), although we did not find differential responses to scent of males and females. The magnitude of the responses was lower to scent of sub adult individuals, and then even lower to scent of juvenile snakes, but in all cases the scent of snakes was discriminated from a blank control. We discuss the potential chemical mechanisms that may allow age recognition and its implications for social and sexual behavior of this snake species.

Haemogregarine Genetic Diversity in Captive African Rock Pythons from Nigeria Suggests a Geographical Pattern.

OBJECTIVE: Understanding the processes that drive parasite diversification, distribution, and abundance is central to disentangle the dynamics and evolution of diseases. In this study, we screened African rock pythons from Nigeria for the presence of blood parasites to assess their distribution, diversity, and phylogenetic relationships. METHODS: A total of 21 captive African rock pythons collected from across 11 locations in Nigeria were sampled between August 2016 and January 2017. Samples were microscopically and genetically analyzed. RESULTS: From the blood smears analyzed, 10 (47.6%) snakes were found to be infected with haemogregarines. Eight of the infected samples were genetically assessed and confirmed as haemogregarines of the recently described Bartazoon group. Two haplotypes were retrieved, of which one was distributed in the northern-central sampled localities and the other in the southern localities. The two haplotypes were clustered in a clade of haemogregarines from snake, gecko, and rodent hosts, and among them, the haemogregarine species Hepatozoon ayorgbor was described from the ball python Python regius. CONCLUSIONS: Two haemogregarine haplotypes in Nigerian rock pythons, which appear to have a geographical pattern across the country, were detected in this study.

Is the Karyotype of Neotropical Boid Snakes Really Conserved? Cytotaxonomy, Chromosomal Rearrangements and Karyotype Organization in the Boidae Family.

Boids are primitive snakes from a basal lineage that is widely distributed in Neotropical region. Many of these species are both morphologically and biogeographically divergent, and the relationship among some species remains uncertain even with evolutionary and phylogenetic studies being proposed for the group. For a better understanding of the evolutionary relationship between these snakes, we cytogenetically analysed 7 species and 3 subspecies of Neotropical snakes from the Boidae family using different chromosomal markers. The karyotypes of Boa constrictor occidentalis, Corallus hortulanus, Eunectes notaeus, Epicrates cenchria and Epicrates assisi are presented here for the first time with the redescriptions of the karyotypes of Boa constrictor constrictor, B. c. amarali, Eunectes murinus and Epicrates crassus. The three subspecies of Boa, two species of Eunectes and three species of Epicrates exhibit 2n = 36 chromosomes. In contrast, C. hortulanus presented a totally different karyotype composition for the Boidae family, showing 2n = 40 chromosomes with a greater number of macrochromosomes. Furthermore, chromosomal mapping of telomeric sequences revealed the presence of interstitial telomeric sites (ITSs) on many chromosomes in addition to the terminal markings on all chromosomes of all taxa analysed, with the exception of E. notaeus. Thus, we demonstrate that the karyotypes of these snakes are not as highly conserved as previously thought. Moreover, we provide an overview of the current cytotaxonomy of the group.

Molecular identification of python species: development and validation of a novel assay for forensic investigations.

Python snake species are often encountered in illegal activities and the question of species identity can be pertinent to such criminal investigations. Morphological identification of species of pythons can be confounded by many issues and molecular examination by DNA analysis can provide an alternative and objective means of identification. Our paper reports on the development and validation of a PCR primer pair that amplifies a segment of the mitochondrial cytochrome b gene that has been suggested previously as a good candidate locus for differentiating python species. We used this DNA region to perform species identification of pythons, even when the template DNA was of poor quality, as might be the case with forensic evidentiary items. Validation tests are presented to demonstrate the characteristics of the assay. Tests involved the cross-species amplification of this marker in non-target species, minimum amount of DNA template required, effects of degradation on product amplification and a blind trial to simulate a casework scenario that provided 100% correct identity. Our results demonstrate that this assay performs reliably and robustly on pythons and can be applied directly to forensic investigations where the presence of a species of python is in question.

[The snakes of Senegal: an annotated species list]. / Les serpents du Sénégal: liste commentée des espèces.

Between 1990 and 2001, the laboratory of Paludologie of IRD at Dakar collected over 5,500 snakes from all over Senegal. By studying this collection, an entirely new species was discovered for science and eight new species for Senegal. The presence of many specimens of rare species and the great number of localities from which snakes were collected also allowed us to solve several delicate taxonomic problems and better to specify the biogeography of Senegalese snakes. Currently, the presence in Senegal of at least 55 different species has been established with certainty. The data collected simultaneously on the burden of snakebites for public health shows considerable differences according to geographic area, with populations from south-eastern Senegal being more exposed to the risk of death by snakebite.

[Transfer of exotic tick Aponomma latum (Koch, 1844) (Acari: Ixodida: Ixodidae) on ball pythons (Python regius Shaw, 1802) brought to Poland]. / Zawlekanie egzotycznego kleszcza Aponomma latum (Koch, 1844) (Acari: Ixodida: Ixodidae) na sprowadzanych do Polski pytonach królewskich (Python regius Shaw, 1802).

103 specimens of Python regius brought to Poland between October 2002 and March 2004 were examined. Occurrence of tick Aponomma latum was reported from 80.6% of the examined reptiles. 549 specimens of A. latum were collected including 341 males, 149 females and 59 nymphs at the various stage of engorgement. Tick A. latum is frequently transferred beyond its natural range of occurrence--Afrotropical region.

Detecting an elusive invasive species: a diagnostic PCR to detect Burmese python in Florida waters and an assessment of persistence of environmental DNA.

Recent studies have demonstrated that detection of environmental DNA (eDNA) from aquatic vertebrates in water bodies is possible. The Burmese python, Python bivittatus, is a semi-aquatic, invasive species in Florida where its elusive nature and cryptic coloration make its detection difficult. Our goal was to develop a diagnostic PCR to detect P. bivittatus from water-borne eDNA, which could assist managers in monitoring this invasive species. First, we used captive P. bivittatus to determine whether reptilian DNA could be isolated and amplified from water samples. We also evaluated the efficacy of two DNA isolation methods and two DNA extraction kits commonly used in eDNA preparation. A fragment of the mitochondrial cytochrome b gene from P. bivittatus was detected in all water samples isolated with the sodium acetate precipitate and the QIAamp DNA Micro Kit. Next, we designed P. bivittatus-specific primers and assessed the degradation rate of eDNA in water. Our primers did not amplify DNA from closely related species, and we found that P. bivittatus DNA was consistently detectable up to 96 h. Finally, we sampled water from six field sites in south Florida. Samples from five sites, where P. bivittatus has been observed, tested positive for eDNA. The final site was negative and had no prior documented evidence of P. bivittatus. This study shows P. bivittatus eDNA can be isolated from water samples; thus, this method is a new and promising technique for the management of invasive reptiles.