The evolutionary network of whiptail lizards reveals predictable outcomes of hybridization

Speciation by hybridization: the mind-boggling nature, educational, and research value of the largest group of unisexual vertebrates

Thirteen species of North American lizards are remarkable because only females exist, which reproduce by cloning unfertilized eggs. Their closest relatives reproduce sexually, with eggs fertilized by sperm from males, as in most vertebrates. The unisexual species originated through hybridization, dispensing with sex and males in a single generation. These lizards hold tremendous potential in science education as a fascinating model for learning about fundamental biological concepts, and in research for developing knowledge with medical applications for reproductive biology, embryonic development, and genetic interactions. These lizards maintain genome integrity in a hybrid state in which recombination is absent, but do not suffer from conditions or disorders such as Down's Syndrome or cancer that are caused by aneuploidy in humans. The multifarious impacts of hybridization on the diversity of species in this group present an exceptional opportunity to deepen understanding of the complicated process of evolutionary diversification.

Allopatric Speciation and Interspecific Gene Flow Driven by Niche Conservatism of Diploderma Tree Lizards in Taiwan

Allopatric speciation is a widely accepted hypothesis for species distributed across geographic barriers. Meanwhile, niche conservatism, the tendency of species to retain their ancestral ecological traits, helps reinforce genetic differentiation by stabilising species distributions over time and reducing the role of competition in shaping range boundaries. In contrast, hybridisation can occur at the edges of distribution after secondary contact following climatic or geological events, leading to a reduction in genetic divergence between divergent lineages. In this study, we investigated the role of geographic barriers, niche conservatism and gene flow in the speciation history of Diploderma species in Taiwan, where geographically distinct taxa share similar environmental preferences. By using ddRAD-seq data, seven distinct genetic clusters were identified with two putatively new cryptic species in D. brevipes and D. polygonatum. Most sister species pairs share similar climatic niches based on niche equivalency and similarity tests. We further detected significant historical gene flow between lineages of D. brevipes and D. polygonatum, where secondary contact might have occurred because of palaeoclimate changes and historical demographic expansion. Our results demonstrate that niche conservatism does not always act in concert to strengthen the result of allopatric speciation; instead, it may also lead to gene flow between divergent lineages following secondary contact. On the other hand, postdivergence gene flow may be a creating force generating phenotypic diversity in sexually selected traits in our study system. The underestimated species diversity of Diploderma in Taiwan requires further taxonomic work in the future.

Genomic data support reticulate evolution in whiptail lizards from the Brazilian Caatinga

Species relationships have traditionally been represented by phylogenetic trees, but not all evolutionary histories fit into bifurcating divergence models. Introgressive hybridization challenges this assumption by sometimes [or maybe often] leading to mitochondrial introgression, wherein one species' mitochondrial genome is entirely replaced by another's (mitochondrial capture). Such processes result in mitonuclear discrepancies, complicating species delimitation and phylogenetic inference. In our study, we used ultraconserved elements (UCE) and mitogenomic data to investigate the evolutionary history of the Ameivula ocellifera complex, a group of South American whiptail lizards widely distributed in semiarid environments of the Caatinga Domain in Brazil. We examine mitonuclear discordances, assessing reticulate evolution, evaluating species limits, and testing for adaptive mitochondrial capture that could explain higher introgression in the mitochondrial genome compared to nuclear DNA. Our findings support the occurrence of an ancient reticulation event during the diversification of these lizards, driven by introgressive hybridization, leading to mitochondrial capture, and explaining mitonuclear discrepancies. Overall, we did not find clear evidence of positive selection across mitochondrial protein-coding genes suggesting adaptive mitochondrial capture of individuals with introgressed mtDNA. Thus, the genetic diversification and mitogenome evolution could be neutral, with selection against hybridization in the autosomal loci only, or even mediated by mitonuclear incompatibilities. Analyses of mtDNA genomes alongside network and species delimitation methods were crucial for identifying and validating individuals with introgressed mtDNA as a distinct species, demonstrating the potential of genome sampling, and using innovative analytical techniques for elucidating speciation processes in the presence of introgressive hybridization.

Natural repeated backcrosses lead to triploidy and tetraploidy in parthenogenetic butterfly lizards (Leiolepis: Agamidae)

Obligatory parthenogenesis in vertebrates is restricted to squamate reptiles and evolved through hybridisation. Parthenogens can hybridise with sexual species, resulting in individuals with increased ploidy levels. We describe two successive hybridisations of the parthenogenetic butterfly lizards (genus Leiolepis) in Vietnam with a parental sexual species. Contrary to previous proposals, we document that parthenogenetic L. guentherpetersi has mitochondrial DNA and two haploid sets from L. guttata and one from L. reevesii, suggesting that it is the result of a backcross of a parthenogenetic L. guttata × L. reevesii hybrid with a L. guttata male increasing ploidy from 2n to 3n. Within the range of L. guentherpetersi, we found an adult tetraploid male with three L. guttata and one L. reevesii haploid genomes. It probably originated from fertilisation of an unreduced triploid L. guentherpetersi egg by a L. guttata sperm. Although its external morphology resembles that of the maternal species, it possessed exceptionally large erythrocytes and was likely sterile. As increased ploidy level above triploidy or tetraploidy appears to be harmful for amniotes, all-female asexual lineages should evolve a strategy to prevent incorporation of other haploid genomes from a sexual species by avoiding fertilisation by sexual males.

Understanding Species Boundaries that Arise from Complex Histories: Gene Flow Across the Speciation Continuum in the Spotted Whiptail Lizards

-Gene flow between diverging lineages challenges the resolution of species boundaries and the understanding of evolutionary history in recent radiations. Here, we integrate phylogenetic and coalescent tools to resolve reticulate patterns of diversification and use a perspective focused on evolutionary mechanisms to distinguish interspecific and intraspecific taxonomic variation. We use this approach to resolve the systematics for one of the most intensively studied but difficult to understand groups of reptiles: the spotted whiptail lizards of the genus Aspidoscelis (A. gularis complex). Whiptails contain the largest number of unisexual species known within any vertebrate group and the spotted whiptail complex has played a key role in the generation of this diversity through hybrid speciation. Understanding lineage boundaries and the evolutionary history of divergence and reticulation within this group is therefore key to understanding the generation of unisexual diversity in whiptails. Despite this importance, long-standing confusion about their systematics has impeded understanding of which gonochoristic species have contributed to the formation of unisexual lineages. Using reduced representation genomic data, we resolve patterns of divergence and gene flow within the spotted whiptails and clarify patterns of hybrid speciation. We find evidence that biogeographically structured ecological and environmental variation has been important in morphological and genetic diversification, as well as the maintenance of species boundaries in this system. Our study elucidates how gene flow among lineages and the continuous nature of speciation can bias the practice of species delimitation and lead taxonomists operating under different frameworks to different conclusions (here we propose that a 2 species arrangement best reflects our current understanding). In doing so, this study provides conceptual and methodological insights into approaches to resolving diversification patterns and species boundaries in rapid radiations with complex histories, as well as long-standing taxonomic challenges in the field of systematic biology.

Notes from the taxonomic disaster zone: Evolutionary drivers of intractable species boundaries in an Australian lizard clade (Scincidae: Ctenotus)

Genomic-scale datasets, sophisticated analytical techniques, and conceptual advances have disproportionately failed to resolve species boundaries in some groups relative to others. To understand the processes that underlie taxonomic intractability, we dissect the speciation history of an Australian lizard clade that arguably represents a "worst-case" scenario for species delimitation within vertebrates: the Ctenotus inornatus species group, a clade beset with decoupled genetic and phenotypic breaks, uncertain geographic ranges, and parallelism in purportedly diagnostic morphological characters. We sampled hundreds of localities to generate a genomic perspective on population divergence, structure, and admixture. Our results revealed rampant paraphyly of nominate taxa in the group, with lineages that are either morphologically cryptic or polytypic. Isolation-by-distance patterns reflect spatially continuous differentiation among certain pairs of putative species, yet genetic and geographic distances are decoupled in other pairs. Comparisons of mitochondrial and nuclear gene trees, tests of nuclear introgression, and historical demographic modelling identified gene flow between divergent candidate species. Levels of admixture are decoupled from phylogenetic relatedness; gene flow is often higher between sympatric species than between parapatric populations of the same species. Such idiosyncratic patterns of introgression contribute to species boundaries that are fuzzy while also varying in fuzziness. Our results suggest that "taxonomic disaster zones" like the C. inornatus species group result from spatial variation in the porosity of species boundaries and the resulting patterns of genetic and phenotypic variation. This study raises questions about the origin and persistence of hybridizing species and highlights the unique insights provided by taxa that have long eluded straightforward taxonomic categorization.

Colonization of North America Boosted the Diversification of Whiptail Lizards

Diversification is frequently associated with change-anything from colonizing a new area to evolving a new trait. Once a lineage changes, the organisms may be able to exploit previously unavailable ecological opportunities and release pressures from predators, parasites, and competitors, which may increase the speciation rate. Modern teiid lizards originated in South America but managed to colonize and diversify in North America. We assessed whether geographic distribution, body size, and body temperatures are associated with teiid diversification using GeoHiSSE and inverse equal-splits statistics with simulation tests. We also estimated speciation rates with MiSSE to account for the effect of unmeasured variables. Moreover, we assessed the ecological niche overlap between North American (including Caribbean) teiids and their sister clade in South America. Our results indicate that only distribution range affected diversification, but we discuss that the available data might not have been enough to assess the effect of body temperatures. We also show that North American teiids have a broader ecological niche encompassing almost all environmental conditions used by their sister clade in South America but expanding mainly toward arid areas. Our results suggest that this expansion significantly impacted teiid diversification due to the seizing of ecological opportunities or ecological release, but we do not discard possible effects of phenotypic evolution.

Post-meiotic mechanism of facultative parthenogenesis in gonochoristic whiptail lizard species

Facultative parthenogenesis (FP) has historically been regarded as rare in vertebrates, but in recent years incidences have been reported in a growing list of fish, reptile, and bird species. Despite the increasing interest in the phenomenon, the underlying mechanism and evolutionary implications have remained unclear. A common finding across many incidences of FP is either a high degree of homozygosity at microsatellite loci or low levels of heterozygosity detected in next-generation sequencing data. This has led to the proposal that second polar body fusion following the meiotic divisions restores diploidy and thereby mimics fertilization. Here, we show that FP occurring in the gonochoristic Aspidoscelis species A. marmoratus and A. arizonae results in genome-wide homozygosity, an observation inconsistent with polar body fusion as the underlying mechanism of restoration. Instead, a high-quality reference genome for A. marmoratus and analysis of whole-genome sequencing from multiple FP and control animals reveals that a post-meiotic mechanism gives rise to homozygous animals from haploid, unfertilized oocytes. Contrary to the widely held belief that females need to be isolated from males to undergo FP, females housed with conspecific and heterospecific males produced unfertilized eggs that underwent spontaneous development. In addition, offspring arising from both fertilized eggs and parthenogenetic development were observed to arise from a single clutch. Strikingly, our data support a mechanism for facultative parthenogenesis that removes all heterozygosity in a single generation. Complete homozygosity exposes the genetic load and explains the high rate of congenital malformations and embryonic mortality associated with FP in many species. Conversely, for animals that develop normally, FP could potentially exert strong purifying selection as all lethal recessive alleles are purged in a single generation.

Phylogenomics of Psammodynastes and Buhoma (Elapoidea: Serpentes), with the description of a new Asian snake family

Asian mock vipers of the genus Psammodynastes and African forest snakes of the genus Buhoma are two genera belonging to the snake superfamily Elapoidea. The phylogenetic placements of Psammodynastes and Buhoma within Elapoidea has been extremely unstable which has resulted in their uncertain and debated taxonomy. We used ultraconserved elements and traditional nuclear and mitochondrial markers to infer the phylogenetic relationships of these two genera with other elapoids. Psammodynastes, for which a reference genome has been sequenced, were found, with strong branch support, to be a relatively early diverging split within Elapoidea that is sister to a clade consisting of Elapidae, Micrelapidae and Lamprophiidae. Hence, we allocate Psammodynastes to its own family, Psammodynastidae new family. However, the phylogenetic position of Buhoma could not be resolved with a high degree of confidence. Attempts to identify the possible sources of conflict in the rapid radiation of elapoid snakes suggest that both hybridisation/introgression during the rapid diversification, including possible ghost introgression, as well as incomplete lineage sorting likely have had a confounding role. The usual practice of combining mitochondrial loci with nuclear genomic data appears to mislead phylogeny reconstructions in rapid radiation scenarios, especially in the absence of genome scale data.

Interspecific behavioural interference and range dynamics: current insights and future directions

Novel biotic interactions in shifting communities play a key role in determining the ability of species' ranges to track suitable habitat. To date, the impact of biotic interactions on range dynamics have predominantly been studied in the context of interactions between different trophic levels or, to a lesser extent, exploitative competition between species of the same trophic level. Yet, both theory and a growing number of empirical studies show that interspecific behavioural interference, such as interspecific territorial and mating interactions, can slow down range expansions, preclude coexistence, or drive local extinction, even in the absence of resource competition. We conducted a systematic review of the current empirical research into the consequences of interspecific behavioural interference on range dynamics. Our findings demonstrate there is abundant evidence that behavioural interference by one species can impact the spatial distribution of another. Furthermore, we identify several gaps where more empirical work is needed to test predictions from theory robustly. Finally, we outline several avenues for future research, providing suggestions for how interspecific behavioural interference could be incorporated into existing scientific frameworks for understanding how biotic interactions influence range expansions, such as species distribution models, to build a stronger understanding of the potential consequences of behavioural interference on the outcome of future range dynamics.

Interspecific gene flow obscures phylogenetic relationships in an important insect pest species complex

As genomic data proliferates, the prevalence of post-speciation gene flow is making species boundaries and relationships increasingly ambiguous. Although current approaches inferring fully bifurcating phylogenies based on concatenated datasets provide simple and robust answers to many species relationships, they may be inaccurate because the models ignore inter-specific gene flow and incomplete lineage sorting. To examine the potential error resulting from ignoring gene flow, we generated both a RAD-seq and a 500 protein-coding loci highly multiplexed amplicon (HiMAP) dataset for a monophyletic group of 12 species defined as the Bactrocera dorsalis sensu lato clade. With some of the world's worst agricultural pests, the taxonomy of the B. dorsalis s.l. clade is important for trade and quarantines. However, taxonomic confusion confounds resolution due to intra- and interspecific phenotypic variation and convergence, mitochondrial introgression across half of the species, and viable hybrids. We compared the topological convergence of our datasets using concatenated phylogenetic and various multispecies coalescent approaches, some of which account for gene flow. All analyses agreed on species delimitation, but there was incongruence between species relationships. Under concatenation, both datasets suggest identical species relationships with mostly high statistical support. However, multispecies coalescent and multispecies network approaches suggest markedly different hypotheses and detected significant gene flow. We suggest that the network approaches are likely more accurate because gene flow violates the assumptions of the concatenated phylogenetic analyses, but the data-reductive requirements of network approaches resulted in reduced statistical support and could not unambiguously resolve gene flow directions. Our study highlights the importance of testing for gene flow, particularly with phylogenomic datasets, even when concatenated approaches receive high statistical support.

Mechanisms of Intrinsic Postzygotic Isolation: From Traditional Genic and Chromosomal Views to Genomic and Epigenetic Perspectives

Intrinsic postzygotic isolation typically appears as reduced viability or fertility of interspecific hybrids caused by genetic incompatibilities between diverged parental genomes. Dobzhansky-Muller interactions among individual genes, and chromosomal rearrangements causing problems with chromosome synapsis and recombination in meiosis, have both long been considered as major mechanisms behind intrinsic postzygotic isolation. Recent research has, however, suggested that the genetic basis of intrinsic postzygotic isolation can be more complex and involves, for example, overall divergence of the DNA sequence or epigenetic changes. Here, we review the mechanisms of intrinsic postzygotic isolation from genic, chromosomal, genomic, and epigenetic perspectives across diverse taxa. We provide empirical evidence for these mechanisms, discuss their importance in the speciation process, and highlight questions that remain unanswered.

The evolutionary outcomes of climate-change-induced hybridization in insect populations

Rapid range shifts are one of the most frequent responses to climate change in insect populations. Climate-induced range shifts can lead to the breakdown of isolation barriers, and thus, to an increase in hybridization and introgression. Long-term evolutionary consequences such as the formation of hybrid zones, introgression, speciation, and extinction have been predicted as a result of climate-induced hybridization. Our review shows that there has been an increase in the number of published cases of climate-induced hybridization in insects, and that the formation of hybrid zones and introgression seems to be, at the moment, the most frequent outcomes. Although introgression is considered positive, since it increases species' genetic diversity, in the long term, it could lead to negative outcomes such as species fusion or genetic swamping.