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.

Facultative parthenogenesis discovered for the first time in Jamaican boa (Chilabothrus subflavus) using novel microsatellite markers.

The mode of reproduction most often seen in snakes is sexual, but studies have noted facultative parthenogenesis in at least six families. Here, we provide evidence for the first observed case of facultative parthenogenesis in a captive Jamaican boa (Chilabothrus subflavus). A 7-year-old female Jamaican boa, isolated since birth, was found to have produced a litter of 15 offspring. To provide molecular DNA evidence of parthenogenesis, 13 new microsatellite loci were isolated in the species. All offspring were found to be homozygous at each locus and only possess alleles found in the dam, implicating that they were born from asexual reproduction. Several developmental abnormalities, including stillbirths and spinal deformities, were noted in the litter which may be explained by their increased level of homozygosity. To preserve genetic diversity in the captive population, research should be conducted to understand the prevalence of this mode of reproduction and to guide future management decisions of this IUCN listed Vulnerable species.

Postprandial cardiac hypertrophy is sustained by mechanics, epigenetic, and metabolic reprogramming in pythons.

Constricting pythons, known for their ability to consume infrequent, massive meals, exhibit rapid and reversible cardiac hypertrophy following feeding. Our primary goal was to investigate how python hearts achieve this adaptive response after feeding. Isolated myofibrils increased force after feeding without changes in sarcomere ultrastructure and without increasing energy cost. Ca2+ transients were prolonged after feeding with no changes in myofibril Ca2+ sensitivity. Feeding reduced titin-based tension, resulting in decreased cardiac tissue stiffness. Feeding also reduced the activity of sirtuins, a metabolically linked class of histone deacetylases, and increased chromatin accessibility. Transcription factor enrichment analysis on transposase-accessible chromatin with sequencing revealed the prominent role of transcription factors Yin Yang1 and NRF1 in postfeeding cardiac adaptation. Gene expression also changed with the enrichment of translation and metabolism. Finally, metabolomics analysis and adenosine triphosphate production demonstrated that cardiac adaptation after feeding not only increased energy demand but also energy production. These findings have broad implications for our understanding of cardiac adaptation across species and hold promise for the development of innovative approaches to address cardiovascular diseases.

Comparative genomics sheds new light on the convergent evolution of infrared vision in snakes.

Infrared vision is a highly specialized sensory system that evolved independently in three clades of snakes. Apparently, convergent evolution occurred in the transient receptor potential ankyrin 1 (TRPA1) proteins of infrared-sensing snakes. However, this gene can only explain how infrared signals are received, and not the transduction and processing of those signals. We sequenced the genome of Xenopeltis unicolor, a key outgroup species of pythons, and performed a genome-wide analysis of convergence between two clades of infrared-sensing snakes. Our results revealed pervasive molecular adaptation in pathways associated with neural development and other functions, with parallel selection on loci associated with trigeminal nerve structural organization. In addition, we found evidence of convergent amino acid substitutions in a set of genes, including TRPA1 and TRPM2. The analysis also identified convergent accelerated evolution in non-coding elements near 12 genes involved in facial nerve structural organization and optic nerve development. Thus, convergent evolution occurred across multiple dimensions of infrared vision in vipers and pythons, as well as amino acid substitutions, non-coding elements, genes and functions. These changes enabled independent groups of snakes to develop and use infrared vision.

The good, the bad and the boa: An unexpected new species of a true boa revealed by morphological and molecular evidence.

Snakes of the genus Boa are outstanding elements of the New World biota with a broad sociological influence on pop culture. Historically, several taxa have been recognized in the past 300 years, being mostly described in the early days of binomial nomenclature. As a rule, these taxa were recognized based on a suite of phenotypic characters mainly those from the external morphology. However, there is a huge disagreement with respect to the current taxonomy and available molecular phylogenies. In order to reconcile both lines of evidence, we investigate the phylogenetic reconstruction (using mitochondrial and nuclear genes) of the genus in parallel to the detailed study of some phenotypic systems from a geographically representative sample of the cis-Andean mainland Boa constrictor. We used cyt-b only (744bp) from 73 samples, and cyt-b, ND4, NTF3, and ODC partial sequences (in a total of 2305 bp) from 35 samples, comprising nine currently recognized taxa (species or subspecies), to infer phylogenetic relationships of boas. Topologies recovered along all the analyses and genetic distances obtained allied to a unique combination of morphological traits (colouration, pholidosis, meristic, morphometric, and male genitalia features) allowed us to recognize B. constrictor lato sensu, B. nebulosa, B. occidentalis, B. orophias and a distinct lineage from the eastern coast of Brazil, which we describe here as a new species, diagnosing it from the previously recognized taxa. Finally, we discuss the minimally necessary changes in the taxonomy of Boa constrictor complex; the value of some usually disregarded phenotypic character system; and we highlight the urgency of continuing environmental policy to preserve one of the most impacted Brazilian hotspots, the Atlantic Forest, which represents an ecoregion full of endemism.

High-Voltage Toxin'Roll: Electrostatic Charge Repulsion as a Dynamic Venom Resistance Trait in Pythonid Snakes.

The evolutionary interplay between predator and prey has significantly shaped the development of snake venom, a critical adaptation for subduing prey. This arms race has spurred the diversification of the components of venom and the corresponding emergence of resistance mechanisms in the prey and predators of venomous snakes. Our study investigates the molecular basis of venom resistance in pythons, focusing on electrostatic charge repulsion as a defense against α-neurotoxins binding to the alpha-1 subunit of the postsynaptic nicotinic acetylcholine receptor. Through phylogenetic and bioactivity analyses of orthosteric site sequences from various python species, we explore the prevalence and evolution of amino acid substitutions that confer resistance by electrostatic repulsion, which initially evolved in response to predatory pressure by Naja (cobra) species (which occurs across Africa and Asia). The small African species Python regius retains the two resistance-conferring lysines (positions 189 and 191) of the ancestral Python genus, conferring resistance to sympatric Naja venoms. This differed from the giant African species Python sebae, which has secondarily lost one of these lysines, potentially due to its rapid growth out of the prey size range of sympatric Naja species. In contrast, the two Asian species Python brongersmai (small) and Python bivittatus (giant) share an identical orthosteric site, which exhibits the highest degree of resistance, attributed to three lysine residues in the orthosteric sites. One of these lysines (at orthosteric position 195) evolved in the last common ancestor of these two species, which may reflect an adaptive response to increased predation pressures from the sympatric α-neurotoxic snake-eating genus Ophiophagus (King Cobras) in Asia. All these terrestrial Python species, however, were less neurotoxin-susceptible than pythons in other genera which have evolved under different predatory pressure as: the Asian species Malayopython reticulatus which is arboreal as neonates and juveniles before rapidly reaching sizes as terrestrial adults too large for sympatric Ophiophagus species to consider as prey; and the terrestrial Australian species Aspidites melanocephalus which occupies a niche, devoid of selection pressure from α-neurotoxic predatory snakes. Our findings underline the importance of positive selection in the evolution of venom resistance and suggest a complex evolutionary history involving both conserved traits and secondary evolution. This study enhances our understanding of the molecular adaptations that enable pythons to survive in environments laden with venomous threats and offers insights into the ongoing co-evolution between venomous snakes and their prey.

Trypanosoma sp. infection in Boa constrictor snakes: morphological, hematological, clinical biochemistry, molecular, and phylogenetic characteristics.

There are few reports of Trypanosoma in snakes, as well as little information about its pathogenicity in these animals. Thus, the present study aimed to characterize Trypanosoma found in Boa constrictor snakes, to verify the influence of the parasitism on hematological and clinical biochemistry parameters, and to perform a phylogenetic study of the isolates. Blood samples from sixty-one boas were analyzed for the presence of trypanosomatids and by hematological and clinical biochemistry assays. The flagellates that were found in this analysis were used for cell culture, morphometry, and molecular analysis. Later, molecular typing phylogenetic studies were performed. Nine positive animals (14.75%) were identified by microscopy analysis. The hematological results showed that parasitized animals presented significantly lower levels of packed cell volume, hemoglobin, mean corpuscular volume, and mean corpuscular hemoglobin. In the leukogram, eosinophils and heterophils counts were higher in parasitized animals. Considering the molecular analyses, the isolates presented a higher identity of the glyceraldehyde-3-phosphate dehydrogenase (GAPDH) and the 18S small subunit ribosomal RNA (SSU rRNA) gene fragments with Trypanosoma serpentis. The phylogenetic tree, using the GAPDH, clustered all isolates with T. serpentis and Trypanosoma cascavelli. This is the first description of T. serpentis parasitizing boas and of the clinical changes caused by trypanosomatid infection in snakes.

A Multiplex RT-PCR Method for the Detection of Reptarenavirus Infection.

Reptarenaviruses cause Boid Inclusion Body Disease (BIBD), a fatal disease of boid snakes with an economic and ecological impact, as it affects both captive and wild constrictor snakes. The clinical picture of BIBD is highly variable but often only limited. Intracytoplasmic inclusion bodies (IB), which develop in most cell types including blood cells, are the pathognomonic hallmark of BIBD; their detection represents the diagnostic gold standard of the disease. However, IBs are not consistently present in clinically healthy reptarenavirus carriers, which can, if undetected, lead to and maintain the spread of the disease within and between snake populations. Sensitive viral detection tools are required for screening and control purposes; however, the genetic diversity of reptarenaviruses hampers the reverse transcription (RT) PCR-based diagnostics. Here, we describe a multiplex RT-PCR approach for the molecular diagnosis of reptarenavirus infection in blood samples. The method allows the detection of a wide range of reptarenaviruses with the detection limit reaching 40 copies per microliter of blood. Using 245 blood samples with a reference RT-PCR result, we show that the technique performs as well as the segment-specific RT-PCRs in our earlier studies. It can identify virus carriers and serve to limit reptarenavirus spreading in captive snake collections.

Demonstration of Parthenogenetic Reproduction in a Pet Ball Python (Python regius) through Analysis of Early-Stage Embryos.

Parthenogenesis is an asexual form of reproduction, normally present in various animal and plant species, in which an embryo is generated from a single gamete. Currently, there are some species for which parthenogenesis is supposed but not confirmed, and the mechanisms that activate it are not well understood. A 10-year-old, wild-caught female ball python (Python regius) laid four eggs without any prior contact with a male. The eggs were not incubated and, after 3 days, were submitted to the University of Parma for analysis due to the suspicion of potential embryo presence. Examination of the egg content revealed residual blood vessels and a small red spot, indicative of an early-stage embryo. DNA was extracted from the three deceased embryos and from the mother's blood, five microsatellites were analyzed to ascertain the origin of the embryos. The captive history data, together with the genetic microsatellite analysis approach, demonstrated the parthenogenetic origin of all three embryos. The embryos were homozygous for each of the maternal microsatellites, suggesting a terminal fusion automixis mode of development.

Reptarenavirus S Segment RNA Levels Correlate with the Presence of Inclusion Bodies and the Number of L Segments in Snakes with Reptarenavirus Infection-Lessons Learned from a Large Breeding Colony.

Reptarenaviruses cause boid inclusion body disease (BIBD), a fatal disease particularly impacting captive boa constrictor collections. The development of cytoplasmic inclusion bodies (IBs) comprising reptarenavirus nucleoprotein (NP) in many cell types of affected snakes is characteristic of BIBD. However, snakes can harbor reptarenaviruses without showing IBs, hence representing carriers and a potential source of transmission. The RNA genome of reptarenaviruses comprises a small (S) and a large (L) segment, and the snakes with BIBD commonly carry a swarm of reptarenavirus segments. To design sensitive and reliable tools for the diagnosis of reptarenavirus infection in snake colonies, we used metatranscriptomics to determine the reptarenavirus segments present in a large boa constrictor breeding colony. The analysis identified one reptarenavirus S segment and three L segments in the colony. The sequence data served to design real-time reverse transcription-PCR (RT-PCR) targeting the found S segment. This allowed us to identify all infected animals and to quantify the S segment RNA levels, which we found to correlate with the presence of IBs. We further found a positive correlation between the number of L segments and the S segment RNA level, which could suggest that L segment excess also contributes to the IB formation. Information on cohousing of the snakes showed a clear association of reptarenavirus infection with cohousing in general and cohousing with infected animals. Information on breeding and offspring confirmed that vertical transmission occurred. Furthermore, our data suggest that some animals might be able to clear the infection or at least exhibit transient or intermittent viremia. IMPORTANCE Boid inclusion body disease (BIBD) is caused by reptarenavirus infection, and while reptarenavirus nucleoprotein is the main component of the inclusion bodies (IBs) characteristic of BIBD, not all reptarenavirus-infected snakes demonstrate IBs in their cells. Identification of infected individuals is critical for controlling the spread of the disease; however, the genetic divergence of reptarenaviruses complicates reverse transcription-PCR (RT-PCR)-based diagnostics. Here, we tested a next-generation-sequencing-based approach to establish a tailored "colony-specific" set of diagnostic tools for the detection of reptarenavirus small (S) and large (L) genome segments. With this approach, we could demonstrate that an S-segment-specific RT-PCR is highly effective in identifying the infected individuals. We further found the S segment RNA level to positively correlate with the presence of IBs and the number of L segments, which could direct future studies to identify the BIBD pathogenetic mechanisms.

Utility of the burmese Python as a model for studying plasticity of extreme physiological systems.

Non-traditional animal models present an opportunity to discover novel biology that has evolved to allow such animals to survive in extreme environments. One striking example is the Burmese python (Python molurus bivittatus), which exhibits extreme physiological adaptation in various metabolic organs after consuming a large meal following long periods of fasting. The response to such a large meal in pythons involves a dramatic surge in metabolic rate, lipid overload in plasma, and massive but reversible organ growth through the course of digestion. Multiple studies have reported the physiological responses in post-prandial pythons, while the specific molecular control of these processes is less well-studied. Investigating the mechanisms that coordinate organ growth and adaptive responses offers the opportunity to gain novel insight that may be able to treat various pathologies in humans. Here, we summarize past research on the post-prandial physiological changes in the Burmese python with a focus on the gastrointestinal tract, heart, and liver. Specifically, we address our recent molecular discoveries in the post-prandial python liver which demonstrate transient adaptations that may reveal new therapeutic targets. Lastly, we explore new biology of the aquaporin 7 gene that is potently upregulated in mammalian cardiac myocytes by circulating factors in post-prandial python plasma.

A community-science approach identifies genetic variants associated with three color morphs in ball pythons (Python regius).

Color morphs in ball pythons (Python regius) provide a unique and largely untapped resource for understanding the genetics of coloration in reptiles. Here we use a community-science approach to investigate the genetics of three color morphs affecting production of the pigment melanin. These morphs-Albino, Lavender Albino, and Ultramel-show a loss of melanin in the skin and eyes, ranging from severe (Albino) to moderate (Lavender Albino) to mild (Ultramel). To identify genetic variants causing each morph, we recruited shed skins of pet ball pythons via social media, extracted DNA from the skins, and searched for putative loss-of-function variants in homologs of genes controlling melanin production in other vertebrates. We report that the Albino morph is associated with missense and non-coding variants in the gene TYR. The Lavender Albino morph is associated with a deletion in the gene OCA2. The Ultramel morph is associated with a missense variant and a putative deletion in the gene TYRP1. Our study is one of the first to identify genetic variants associated with color morphs in ball pythons and shows that pet samples recruited from the community can provide a resource for genetic studies in this species.

Boid Inclusion Body Disease Is Also a Disease of Wild Boa Constrictors.

Reptarenaviruses cause boid inclusion body disease (BIBD), a potentially fatal disease, occurring in captive constrictor snakes boas and pythons worldwide. Classical BIBD, characterized by the formation of pathognomonic cytoplasmic inclusion bodies (IBs), occurs mainly in boas, whereas in pythons, for example, reptarenavirus infection most often manifests as central nervous system signs with limited IB formation. The natural hosts of reptarenaviruses are unknown, although free-ranging/wild constrictor snakes are among the suspects. Here, we report BIBD with reptarenavirus infection in indigenous captive and wild boid snakes in Costa Rica using histology, immunohistology, transmission electron microscopy, and next-generation sequencing (NGS). The snakes studied represented diagnostic postmortem cases of captive and wild-caught snakes since 1989. The results from NGS on archival paraffin blocks confirm that reptarenaviruses were already present in wild boa constrictors in Costa Rica in the 1980s. Continuous sequences that were de novo assembled from the low-quality RNA obtained from paraffin-embedded tissue allowed the identification of a distinct pair of reptarenavirus S and L segments in all studied animals; in most cases, reference assembly could recover almost complete segments. Sampling of three prospective cases in 2018 allowed an examination of fresh blood or tissues and resulted in the identification of additional reptarenavirus segments and hartmanivirus coinfection. Our results show that BIBD is not only a disease of captive snakes but also occurs in indigenous wild constrictor snakes in Costa Rica, suggesting boa constrictors to play a role in natural reptarenavirus circulation. IMPORTANCE The literature describes cases of boid inclusion body disease (BIBD) in captive snakes since the 1970s, and in the 2010s, others and ourselves identified reptarenaviruses as the causative agent. BIBD affects captive snakes globally, but the origin and the natural host of reptarenaviruses remain unknown. In this report, we show BIBD and reptarenavirus infections in two native Costa Rican constrictor snake species, and by studying archival samples, we show that both the viruses and the disease have been present in free-ranging/wild snakes in Costa Rica at least since the 1980s. The diagnosis of BIBD in wild boa constrictors suggests that this species plays a role in the circulation of reptarenaviruses. Additional sample collection and analysis would help to clarify this role further and the possibility of, e.g., vector transmission from an arthropod host.

Reference genome of the rubber boa, Charina bottae (Serpentes: Boidae).

The rubber boa, Charina bottae is a semi-fossorial, cold-temperature adapted snake that ranges across the wetter and cooler ecoregions of the California Floristic Province. The rubber boa is 1 of 2 species in the family Boidae native to California and currently has 2 recognized subspecies, the Northern rubber boa C. bottae bottae and the Southern rubber boa C. bottae umbratica. Recent genomic work on C. bottae indicates that these 2 subspecies are collectively composed of 4 divergent lineages that separated during the late Miocene. Analysis of habitat suitability indicates that C. bottae umbratica montane sky-island populations from southern California will lose the majority of their habit over the next 70 yr, and is listed as Threatened under the California Endangered Species Act. Here, we report a new, chromosome-level assembly of C. bottae bottae as part of the California Conservation Genomics Project (CCGP). Consistent with the reference genome strategy of the CCGP, we used Pacific Biosciences HiFi long reads and Hi-C chromatin-proximity sequencing technology to produce a de novo assembled genome. The assembly comprises 289 scaffolds covering 1,804,944,895 bp, has a contig N50 of 37.3 Mb, a scaffold N50 of 97 Mb, and BUSCO completeness score of 96.3%, and represents the first reference genome for the Boidae snake family. This genome will enable studies of genetic differentiation and connectivity among C. bottae bottae and C. bottae umbratica populations across California and help manage locally endemic lineages as they confront challenges from human-induced climate warming, droughts, and wildfires across California.

Cytogenetic Analysis of the Members of the Snake Genera Cylindrophis, Eryx, Python, and Tropidophis.

The recent discovery of two independently evolved XX/XY sex determination systems in the snake genera Python and Boa sparked a new drive to study the evolution of sex chromosomes in poorly studied lineages of snakes, where female heterogamety was previously assumed. Therefore, we examined seven species from the genera Eryx, Cylindrophis, Python, and Tropidophis by conventional and molecular cytogenetic methods. Despite the fact that these species have similar karyotypes in terms of chromosome number and morphology, we detected variability in the distribution of heterochromatin, telomeric repeats, and rDNA loci. Heterochromatic blocks were mainly detected in the centromeric regions in all species, although accumulations were detected in pericentromeric and telomeric regions in a few macrochromosomes in several of the studied species. All species show the expected topology of telomeric repeats at the edge of all chromosomes, with the exception of Eryx muelleri, where additional accumulations were detected in the centromeres of three pairs of macrochromosomes. The rDNA loci accumulate in one pair of microchromosomes in all Eryx species and in Cylindrophis ruffus, in one macrochromosome pair in Tropidophis melanurus and in two pairs of microchromosomes in Python regius. Sex-specific differences were not detected, suggesting that these species likely have homomorphic, poorly differentiated sex chromosomes.

Animal Models of Exercise From Rodents to Pythons.

Acute and chronic animal models of exercise are commonly used in research. Acute exercise testing is used, often in combination with genetic, pharmacological, or other manipulations, to study the impact of these manipulations on the cardiovascular response to exercise and to detect impairments or improvements in cardiovascular function that may not be evident at rest. Chronic exercise conditioning models are used to study the cardiac phenotypic response to regular exercise training and as a platform for discovery of novel pathways mediating cardiovascular benefits conferred by exercise conditioning that could be exploited therapeutically. The cardiovascular benefits of exercise are well established, and, frequently, molecular manipulations that mimic the pathway changes induced by exercise recapitulate at least some of its benefits. This review discusses approaches for assessing cardiovascular function during an acute exercise challenge in rodents, as well as practical and conceptual considerations in the use of common rodent exercise conditioning models. The case for studying feeding in the Burmese python as a model for exercise-like physiological adaptation is also explored.

Short '1.2× Genome' Infectious Clone Initiates Kolmiovirid Replication in Boa constrictor Cells.

Human hepatitis D virus (HDV) depends on hepatitis B virus co-infection and its glycoproteins for infectious particle formation. HDV was the sole known deltavirus for decades and believed to be a human-only pathogen. However, since 2018, several groups reported finding HDV-like agents from various hosts but without co-infecting hepadnaviruses. In vitro systems enabling helper virus-independent replication are key for studying the newly discovered deltaviruses. Others and we have successfully used constructs containing multimers of the deltavirus genome for the replication of various deltaviruses via transfection in cell culture. Here, we report the establishment of deltavirus infectious clones with 1.2× genome inserts bearing two copies of the genomic and antigenomic ribozymes. We used Swiss snake colony virus 1 as the model to compare the ability of the previously reported "2× genome" and the "1.2× genome" infectious clones to initiate replication in cell culture. Using immunofluorescence, qRT-PCR, immuno- and northern blotting, we found the 2× and 1.2× genome clones to similarly initiate deltavirus replication in vitro and both induced a persistent infection of snake cells. The 1.2× genome constructs enable easier introduction of modifications required for studying deltavirus replication and cellular interactions.

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.

Revisiting the Karyotype Evolution of Neotropical Boid Snakes: A Puzzle Mediated by Chromosomal Fissions.

The Boidae family is an ancient group of snakes widely distributed across the Neotropical region, where several biogeographic events contributed towards shaping their evolution and diversification. Most species of this family have a diploid number composed of 2n = 36; however, among Booidea families, the Boidae stands out by presenting the greatest chromosomal diversity, with 2n ranging between 36 and 44 chromosomes and an undifferentiated XY sex chromosome system. Here, we applied a comparative chromosome analysis using cross-species chromosome paintings in five species representing four Boidae genera, to decipher the evolutionary dynamics of some chromosomes in these Neotropical snakes. Our study included all diploid numbers (2n = 36, 40, and 44) known for this family and our comparative chromosomal mappings point to a strong evolutionary relationship among the genera Boa, Corallus, Eunectes, and Epicrates. The results also allowed us to propose the cytogenomic diversification that had occurred in this family: a process mediated by centric fissions, including fission events of the putative and undifferentiated XY sex chromosome system in the 2n = 44 karyotype, which is critical in solving the puzzle of the karyotype evolution of boid snakes.