Patterns, Mechanisms and Genetics of Speciation in Reptiles and Amphibians

Genetic basis and expression of ventral colour in polymorphic common lizards

Colour is an important visual cue that can correlate with sex, behaviour, life history or ecological strategies, and has evolved divergently and convergently across animal lineages. Its genetic basis in non-model organisms is rarely known, but such information is vital for determining the drivers and mechanisms of colour evolution. Leveraging genetic admixture in a rare contact zone between oviparous and viviparous common lizards (Zootoca vivipara), we show that females (N = 558) of the two otherwise morphologically indistinguishable reproductive modes differ in their ventral colouration (from pale to vibrant yellow) and intensity of melanic patterning. We find no association between female colouration and reproductive investment, and no evidence for selection on colour. Using a combination of genetic mapping and transcriptomic evidence, we identified two candidate genes associated with ventral colour differentiation, DGAT2 and PMEL. These are genes known to be involved in carotenoid metabolism and melanin synthesis respectively. Ventral melanic spots were associated with two genomic regions, including a SNP close to protein tyrosine phosphatase (PTP) genes. Using genome re-sequencing data, our results show that fixed coding mutations in the candidate genes cannot account for differences in colouration. Taken together, our findings show that the evolution of ventral colouration and its associations across common lizard lineages is variable. A potential genetic mechanism explaining the flexibility of ventral colouration may be that colouration in common lizards, but also across squamates, is predominantly driven by regulatory genetic variation.

Maternal transmission of bacterial microbiota during embryonic development in a viviparous lizard

IMPORTANCE

We investigated the presence and diversity of bacteria in the embryos of the viviparous lizard Sceloporus grammicus and their amniotic environment. We compared this diversity to that found in the maternal intestine, mouth, and cloaca. We detected bacterial DNA in the embryos, albeit with a lower bacterial species diversity than found in maternal tissues. Most of the bacterial species detected in the embryos were also found in the mother, although not all of them. Interestingly, we detected a high similarity in the composition of bacterial species among embryos from different mothers. These findings suggest that there may be a mechanism controlling the transmission of bacteria from the mother to the embryo. Our results highlight the possibility that the interaction between maternal bacteria and the embryo may affect the development of the lizards.

Repeated divergence of amphibians and reptiles across an elevational gradient in northern Madagascar

How environmental factors shape patterns of biotic diversity in tropical ecosystems is an active field of research, but studies examining the possibility of ecological speciation in terrestrial tropical ecosystems are scarce. We use the isolated rainforest herpetofauna on the Montagne d'Ambre (Amber Mountain) massif in northern Madagascar as a model to explore elevational divergence at the level of populations and communities. Based on intensive sampling and DNA barcoding of amphibians and reptiles along a transect ranging from ca. 470-1470 m above sea level (a.s.l.), we assessed a main peak in species richness at an elevation of ca. 1000 m a.s.l. with 41 species. The proportion of local endemics was highest (about 1/3) at elevations >1100 m a.s.l. Two species of chameleons (Brookesia tuberculata, Calumma linotum) and two species of frogs (Mantidactylus bellyi, M. ambony) studied in depth by newly developed microsatellite markers showed genetic divergence up the slope of the mountain, some quite strong, others very weak, but in each case with genetic breaks between 1100 and 1270 m a.s.l. Genetic clusters were found in transect sections significantly differing in bioclimate and herpetological community composition. A decrease in body size was detected in several species with increasing elevation. The studied rainforest amphibians and reptiles show concordant population genetic differentiation across elevation along with morphological and niche differentiation. Whether this parapatric or microallopatric differentiation will suffice for the completion of speciation is, however, unclear, and available phylogeographic evidence rather suggests that a complex interplay between ecological and allopatric divergence processes is involved in generating the extraordinary species diversity of Madagascar's biota. Our study reveals concordant patterns of diversification among main elevational bands, but suggests that these adaptational processes are only part of the complex of processes leading to species formation, among which geographical isolation is probably also important.

Allometric analysis of sexual dimorphism and morphological variation in two chromosome races of the Sceloporus grammicus complex (Squamata: Phrynosomatidae) from Mexico

Sexual dimorphism is a widespread feature in the Animal Kingdom. In lizards of the Sceloporus grammicus complex, studies of sexual dimorphism that analyze the allometric trajectories of body traits remain unexplored. Here we investigate sexual dimorphism in key phenotypic traits, including body size (snout-vent length, SVL) as well as head length (HL), head width (HW), and forearm length (FL). We use an allometric approach to detect differences in scale relationships among body parts in the S. grammicus complex in Mexico. We focus on two chromosomal races within this complex, F5 (2n = 34) and FM2 (2n = 46). In the complex, we found that males are larger than females in all morphological variables, and this pattern was confirmed in both races. We determined negative allometric trajectories (SVLvs.HL and HW), isometry (SVLvs.FL) and intersexual differences in the slopes of the SVLvs.HL and HW; the males showed steeper slopes. Thus, the growth of the head is more pronounced in males than females. Additionally, we found between-race differences in these trajectories (SVLvs.FL) and in all morphological variables (F5 lizards are larger than those of the FM2 race), which correlate with their chromosomal divergence. We discuss biological implications of our findings in relation to sexual selection and natural selection.

Slow and steady wins the race: Diversification rate is independent from body size and lifestyle in Malagasy skinks (Squamata: Scincidae: Scincinae)

Most of the unique and diverse vertebrate fauna that inhabits Madagascar derives from in situ diversification from colonisers that reached this continental island through overseas dispersal. The endemic Malagasy Scincinae lizards are amongst the most species-rich squamate groups on the island. They colonised all bioclimatic zones and display many ecomorphological adaptations to a fossorial (burrowing) lifestyle. Here we propose a new phylogenetic hypothesis for their diversification based on the largest taxon sampling so far compiled for this group. We estimated divergence times and investigated several aspects of their diversification (diversification rate, body size and fossorial lifestyle evolution, and biogeography). We found that diversification rate was constant throughout most of the evolutionary history of the group, but decreased over the last 6-4 million years and independently from body size and fossorial lifestyle evolution. Fossoriality has evolved from fully quadrupedal ancestors at least five times independently, which demonstrates that even complex morphological syndromes - in this case involving traits such as limb regression, body elongation, modification of cephalic scalation, depigmentation, and eyes and ear-opening regression - can evolve repeatedly and independently given enough time and eco-evolutionary advantages. Initial diversification of the group likely occurred in forests, and the divergence of sand-swimmer genera around 20 Ma appears linked to a period of aridification. Our results show that the large phenotypic variability of Malagasy Scincinae has not influenced diversification rate and that their rich species diversity results from a constant accumulation of lineages through time. By compiling large geographic and trait-related datasets together with the computation of a new time tree for the group, our study contributes important insights on the diversification of Malagasy vertebrates.

Genomic analyses of the scorpion mud turtle (Kinosternon scorpioides) (Linnaeus, 1766) in insular and continental Colombia: Evidence for multiple conservation and taxonomic units

The turtle genus Kinosternon is widespread with at least 25 species distributed from Mexico to northern Argentina. The taxonomy of this genus is controversial and requires a full revision using both morphological and molecular approaches. In this study, we did a genomic analysis on the species Kinosternon scorpioides distributed in insular and continental Colombia in order to define conservation units. Total DNA was extracted from 24 tissue samples and RADseq genotyping analysis was done. In addition, the intron R35 was amplified and sequenced for a subset of samples. A total of 35,507 SNPs combined with 1,047 bp of the intron were used for spatiotemporal colonization pattern reconstruction and phylogenetic analyses. In addition, SNPs were used for population structure inferences and allele frequency-based analyses. Reciprocal monophyly, significant differences in allele frequencies (Fst = 0.32 - 0.78), and evidence of reproductive isolation (no admixture/geneflow), indicate long-term divergence between groups (2-8 MYA), possibly due to geographical barriers. Four Evolutionarily Significant Units (ESUs) were defined within our samples. One ESU was represented by the insular subspecies K. scorpioides albogulare, found on San Andrés island, and three ESUs were defined for the subspecies K. s. scorpioides in continental Colombia: one trans-Andean, found in northwestern Colombia (Caribbean region) and two cis-Andean, found in eastern and southeastern Colombia in the Orinoco and Amazon regions, respectively. Colonization of this species occurred from an ancestral area on South of Central America region (~ 8.43 MYA), followed by an establishing of current populations on San Andrés Island and then, in the continent. First, in the Colombian Caribbean, next, in the Orinoco, and more recently, in the Amazon. We hypothesize that the emergence of the Panamá Isthmus, as well as the final uplift of the North Eastern Andes and Vaupes Arch, were key event leading to the differentiation of these ESUs. For management and conservation purposes, each of these ESUs should be considered as a separate management unit. A full revision of the taxonomy of the genus Kinosternon is warranted.

Sex chromosomes as supergenes of speciation: why amphibians defy the rules?

As reflected by the two rules of speciation (Haldane's rule and the large X-/Z-effect), sex chromosomes are expected to behave like supergenes of speciation: they recombine only in one sex (XX females or ZZ males), supposedly recruit sexually antagonistic genes and evolve faster than autosomes, which can all contribute to pre-zygotic and post-zygotic isolation. While this has been mainly studied in organisms with conserved sex-determining systems and highly differentiated (heteromorphic) sex chromosomes like mammals, birds and some insects, these expectations are less clear in organismal groups where sex chromosomes repeatedly change and remain mostly homomorphic, like amphibians. In this article, we review the proposed roles of sex-linked genes in isolating nascent lineages throughout the speciation continuum and discuss their support in amphibians given current knowledge of sex chromosome evolution and speciation modes. Given their frequent recombination and lack of differentiation, we argue that amphibian sex chromosomes are not expected to become supergenes of speciation, which is reflected by the rarity of empirical studies consistent with a 'large sex chromosome effect' in frogs and toads. The diversity of sex chromosome systems in amphibians has a high potential to disentangle the evolutionary mechanisms responsible for the emergence of sex-linked speciation genes in other organisms. This article is part of the theme issue 'Genomic architecture of supergenes: causes and evolutionary consequences'.

Composition of a chemical signalling trait varies with phylogeny and precipitation across an Australian lizard radiation

The environment presents challenges to the transmission and detection of animal signalling systems, resulting in selective pressures that can drive signal divergence amongst populations in disparate environments. For chemical signals, climate is a potentially important selective force because factors such as temperature and moisture influence the persistence and detection of chemicals. We investigated an Australian lizard radiation (Heteronotia) to explore relationships between a sexually dimorphic chemical signalling trait (epidermal pore secretions) and two key climate variables: temperature and precipitation. We reconstructed the phylogeny of Heteronotia with exon capture phylogenomics, estimated phylogenetic signal in amongst-lineage chemical variation and assessed how chemical composition relates to temperature and precipitation using multivariate phylogenetic regressions. High estimates of phylogenetic signal indicate that the composition of epidermal pore secretions varies amongst lineages in a manner consistent with Brownian motion, although there are deviations to this, with stark divergences coinciding with two phylogenetic splits. Accounting for phylogenetic non-independence, we found that amongst-lineage chemical variation is associated with geographic variation in precipitation but not temperature. This contrasts somewhat with previous lizard studies, which have generally found an association between temperature and chemical composition. Our results suggest that geographic variation in precipitation can affect the evolution of chemical signalling traits, possibly influencing patterns of divergence amongst lineages and species.

Population diversification in the frog Mantidactylus bellyi on an isolated massif in northern Madagascar based on genetic, morphological, bioacoustic and ecological evidence

In the processes that give rise to new species, changes first occur at the population level. But with the continuous nature of the divergence process, change in biological properties delimiting the shift from “individuals of divergent populations” towards “individuals of distinct species”, as well as abiotic factors driving the change, remain largely ambivalent. Here we study diversification processes at the population level in a semi-aquatic frog, Mantidactylus (Brygoomantis) bellyi, across the diverse vegetation types of Montagne d’Ambre National Park (MANP), Madagascar. Genetic diversity was assessed with seven newly developed microsatellite markers as well as mitochondrial DNA sequences and concordance with patterns of ecological, morphological, and bioacoustic divergence evaluated. We found M. bellyi lacking mitochondrial differentiation within MANP, while microsatellite datasets partitioned them into three highly differentiated, geographically separated subpopulations (with indications for up to five subpopulations). The molecular grouping–primarily clustering individuals by geographic proximity–was coincident with differences in mean depth and width of waters, suggesting a possible role of fluvial characteristics in genetic exchange in this stream-breeding species. Genetic clustering not consistent with differences in call properties, except for dominant call frequencies under the two-subpopulations model. Morphological divergence was mostly consistent with the genetic clustering; subpopulations strongly differed by their snout-vent length, with individuals from high-elevation subpopulations smaller than those from populations below 1000 m above sea level. These results exemplify how mountains and environmental conditions might primarily shape genetic and morphological divergence in frog populations, without strongly affecting their calls.

Genome and cuticular hydrocarbon-based species delimitation shed light on potential drivers of speciation in a Neotropical ant species complex

Geographic separation that leads to the evolution of reproductive isolation between populations generally is considered the most common form of speciation. However, speciation may also occur in the absence of geographic barriers due to phenotypic and genotypic factors such as chemical cue divergence, mating signal divergence, and mitonuclear conflict. Here, we performed an integrative study based on two genome-wide techniques (3RAD and ultraconserved elements) coupled with cuticular hydrocarbon (CHC) and mitochondrial (mt) DNA sequence data, to assess the species limits within the Ectatomma ruidum species complex, a widespread and conspicuous group of Neotropical ants for which heteroplasmy (i.e., presence of multiple mtDNA variants in an individual) has been recently discovered in some populations from southeast Mexico. Our analyses indicate the existence of at least five distinct species in this complex: two widely distributed across the Neotropics, and three that are restricted to southeast Mexico and that apparently have high levels of heteroplasmy. We found that species boundaries in the complex did not coincide with geographic barriers. We therefore consider possible roles of alternative drivers that may have promoted the observed patterns of speciation, including mitonuclear incompatibility, CHC differentiation, and colony structure. Our study highlights the importance of simultaneously assessing different sources of evidence to disentangle the species limits of taxa with complicated evolutionary histories.

Evidence for ecological processes driving speciation among endemic lizards of Madagascar

Although genetic patterns produced by population isolation during speciation are well documented, the biogeographic and ecological processes that trigger speciation remain poorly understood. Alternative hypotheses for the biogeography and ecology of speciation include geographic isolation combined with niche conservation (soft allopatry) or parapatric distribution on an environmental gradient with niche divergence (ecological speciation). Here, we use species' distributions, environmental data, and two null models (the Random Translation and Rotation and the Background Similarity Test) to test these alternative hypotheses among 28 sister pairs of microendemic lizards in Madagascar. Our results demonstrate strong bimodal peaks along a niche divergence-conservation spectrum, with at least 25 out of 28 sister pairs exhibiting either niche conservation or divergence, and the remaining pairs showing weak ecological signals. Yet despite these significant results, we do not find strong associations of niche conservation with allopatric distributions or niche divergence with parapatric distributions. Our findings thus provide strong evidence of a role for ecological processes driving speciation, rather than the classic expectation of speciation through geographic isolation, but demonstrate that the link between ecological speciation and parapatry is complex and requires further analysis of a broader taxonomic sample to fully resolve.

Ecological Divergence and the History of Gene Flow in the Nearctic Milksnakes (Lampropeltis triangulum Complex)

Many phylogeographic studies on species with large ranges have found genetic-geographic structure associated with changes in habitat and physical barriers preventing or reducing gene flow. These interactions with geographic space, contemporary and historical climate, and biogeographic barriers have complex effects on contemporary population genetic structure and processes of speciation. While allopatric speciation at biogeographic barriers is considered the primary mechanism for generating species, more recently it has been shown that parapatric modes of divergence may be equally or even more common. With genomic data and better modeling capabilities, we can more clearly define causes of speciation in relation to biogeography and migration between lineages, the location of hybrid zones with respect to the ecology of parental lineages, and differential introgression of genes between taxa. Here, we examine the origins of three Nearctic milksnakes (Lampropeltis elapsoides, Lampropeltis triangulum and Lampropeltis gentilis) using genome-scale data to better understand species diversification. Results from artificial neural networks show that a mix of a strong biogeographic barrier, environmental changes, and physical space has affected genetic structure in these taxa. These results underscore conspicuous environmental changes that occur as the sister taxa L. triangulum and L. gentilis diverged near the Great Plains into the forested regions of the Eastern Nearctic. This area has been recognized as a region for turnover for many vertebrate species, but as we show here the contemporary boundary does not isolate these sister species. These two species likely formed in the mid-Pleistocene and have remained partially reproductively isolated over much of this time, showing differential introgression of loci. We also demonstrate that when L. triangulum and L. gentilis are each in contact with the much older L. elapsoides, some limited gene flow has occurred. Given the strong agreement between nuclear and mtDNA genomes, along with estimates of ecological niche, we suggest that all three lineages should continue to be recognized as unique species. Furthermore, this work emphasizes the importance of considering complex modes of divergence and differential allelic introgression over a complex landscape when testing mechanisms of speciation. [Cline; delimitation; Eastern Nearctic; Great Plains; hybrids; introgression; speciation.].

De novo venom gland transcriptomics of Calliophis bivirgata flaviceps: uncovering the complexity of toxins from the Malayan blue coral snake

Background:

The Malayan blue coral snake, Calliophis bivirgata flaviceps, is a medically important venomous snake in Southeast Asia. However, the complexity and diversity of its venom genes remain little explored.

Methods:

To address this, we applied high-throughput next-generation sequencing to profile the venom gland cDNA libraries of C. bivirgata flaviceps. The transcriptome was de novo assembled, followed by gene annotation, multiple sequence alignment and analyses of the transcripts.

Results:

A total of 74 non-redundant toxin-encoding genes from 16 protein families were identified, with 31 full-length toxin transcripts. Three-finger toxins (3FTx), primarily delta-neurotoxins and cardiotoxin-like/cytotoxin-like proteins, were the most diverse and abundantly expressed. The major 3FTx (Cb_FTX01 and Cb_FTX02) are highly similar to calliotoxin, a delta-neurotoxin previously reported in the venom of C. bivirgata. This study also revealed a conserved tyrosine residue at position 4 of the cardiotoxin-like/cytotoxin-like protein genes in the species. These variants, proposed as Y-type CTX-like proteins, are similar to the H-type CTX from cobras. The substitution is conservative though, preserving a less toxic form of elapid CTX-like protein, as indicated by the lack of venom cytotoxicity in previous laboratory and clinical findings. The ecological role of these toxins, however, remains unclear. The study also uncovered unique transcripts that belong to phospholipase A₂ of Groups IA and IB, and snake venom metalloproteinases of PIII subclass, which show sequence variations from those of Asiatic elapids.

Conclusion:

The venom gland transcriptome of C. bivirgata flaviceps from Malaysia was de novo assembled and annotated. The diversity and expression profile of toxin genes provide insights into the biological and medical importance of the species.

Mass of genes rather than master genes underlie the genomic architecture of amphibian speciation

The genetic architecture of speciation, i.e., how intrinsic genomic incompatibilities promote reproductive isolation (RI) between diverging lineages, is one of the best-kept secrets of evolution. To directly assess whether incompatibilities arise in a limited set of large-effect speciation genes, or in a multitude of loci, we examined the geographic and genomic landscapes of introgression across the hybrid zones of 41 pairs of frog and toad lineages in the Western Palearctic region. As the divergence between lineages increases, phylogeographic transitions progressively become narrower, and larger parts of the genome resist introgression. This suggests that anuran speciation proceeds through a gradual accumulation of multiple barrier loci scattered across the genome, which ultimately deplete hybrid fitness by intrinsic postzygotic isolation, with behavioral isolation being achieved only at later stages. Moreover, these loci were disproportionately sex linked in one group (Hyla) but not in others (Rana and Bufotes), implying that large X-effects are not necessarily a rule of speciation with undifferentiated sex chromosomes. The highly polygenic nature of RI and the lack of hemizygous X/Z chromosomes could explain why the speciation clock ticks slower in amphibians compared to other vertebrates. The clock-like dynamics of speciation combined with the analytical focus on hybrid zones offer perspectives for more standardized practices of species delimitation.

A new species of wood lizard (Hoplocercinae, Enyalioides) from the Río Huallaga Basin in Central Peru

We report the discovery of a new species of Enyalioides from the premontane forest of the Río Huallaga basin in central Peru. The most similar and phylogenetically related species are E. binzayedi and E. rudolfarndti. However, the new species differs from E. binzayedi (state of character in parentheses) by having dorsal scales strongly keeled on paravertebral region and feebly keeled or smooth elsewhere (prominent medial keel on each dorsal scale), more dorsals in transverse row between dorsolateral crests at midbody 26–39, x̄ = 30.44 (22–31, x̄ = 27.57), and a conspicuous posteromedial black patch in the gular region of males (absent). Contrarily, adult males of the new species and E. rudolfarndti are readily distinguished by having a conspicuous posteromedial black patch in the gular region, absent in E. rudolfarndti, and by lacking a conspicuous orange blotch (faint if present) on the antehumeral region, as in E. rudolfarndti. We also present an updated molecular phylogenetic tree of hoplocercines, which strongly supports both referral of the newly discovered species to Enyalioides and its status as distinct from those recognized previously.

Heterogeneous relationships between rates of speciation and body size evolution across vertebrate clades

Several theories predict that rates of phenotypic evolution should be related to the rate at which new lineages arise. However, drawing general conclusions regarding the coupling between these fundamental evolutionary rates has been difficult due to the inconsistent nature of previous results combined with uncertainty over the most appropriate methodology with which to investigate such relationships. Here we propose and compare the performance of several different approaches for testing associations between lineage-specific rates of speciation and phenotypic evolution using phylogenetic data. We then use the best-performing method to test relationships between rates of speciation and body size evolution in five major vertebrate clades (amphibians, birds, mammals, ray-finned fish and squamate reptiles) at two phylogenetic scales. Our results provide support for the long-standing view that rates of speciation and morphological evolution are generally positively related at broad macroevolutionary scales, but they also reveal a substantial degree of heterogeneity in the strength and direction of these associations at finer scales across the vertebrate tree of life.

A new species of frog (Terrarana, Strabomantidae, Phrynopus) from the Peruvian Andean grasslands

We describe a new, medium-sized species of terrestrial frog of the genus Phrynopus from a single locality in the central Andes of Peru (Departamento de Huánuco) at 3,730 meters of elevation. Phylogenetic analyses supported Phrynopus remotum sp. nov. as an independent lineage, sister to most of its congeners. The new species is morphologically distinguishable by the presence of small tubercles on upper eyelids and heels, an areolate venter, and the absence of dorsolateral folds or ridges. This species inhabits the highlands adjacent to the Marañón Dry valley. The only sympatric amphibian species recorded is the marsupial frog Gastrotheca peruana.

Integrating hybrid zone analyses in species delimitation: lessons from two anuran radiations of the Western Mediterranean

Molecular ecologists often rely on phylogenetic evidence for assessing the species-level systematics of newly discovered lineages. Alternatively, the extent of introgression at phylogeographic transitions can provide a more direct test to assign candidate taxa into subspecies or species categories. Here, we compared phylogenetic versus hybrid zone approaches of species delimitation in two groups of frogs from the Western Mediterranean region (Discoglossus and Pelodytes), by using genomic data (ddRAD). In both genera, coalescent analyses recovered almost all nominal taxa as "species". However, the least-diverged pairs D. g. galganoi/jeanneae and P. punctatus/hespericus admix over hundreds of kilometers, suggesting that they have not yet developed strong reproductive isolation and should be treated as conspecifics. In contrast, the comparatively older D. scovazzi/pictus and P. atlanticus/ibericus form narrow contact zones, consistent with species distinctiveness. Due to their complementarity, we recommend taxonomists to combine phylogenomics with hybrid zone analyses to scale the gray zone of speciation, i.e., the evolutionary window separating widely admixing lineages versus nascent reproductively isolated species. The radically different transitions documented here conform to the view that genetic incompatibilities accumulating with divergence generate a weak barrier to gene flow for long periods of time, until their effects multiply and the speciation process then advances rapidly. Given the variability of the gray zone among taxonomic groups, at least from our current abilities to measure it, we recommend to customize divergence thresholds within radiations to categorize lineages for which no direct test of speciation is possible.