COI Barcodes combined with multilocus data for representative Aporia taxa shed light on speciation in the high altitude Irano-Turanian mountain plateaus (Lepidoptera: Pieridae).

Even though the high plateaus of Qinghai-Tibet and Iran share many faunal elements, the historical biogeography of the species present in this area are not very well understood. We present a complete COI barcode library for Aporia Hübner and a first comprehensive phylogeny for the genus including all known species and majority of subspecies using ten available genes (COI-COII, ND1, ND5, Cytb, EF-1a, Wg, 16S, 28S-D2/D3 and 28S-D8). We then focus on A. leucodice (Eversmann, 1843) and related taxa in order to resolve some long-standing taxonomic issues in this species-group. Based on DNA sequence data as well as morphology, we raise Aporia illumina (Grum-Grshimailo 1890) stat. nov. (= pseudoillumina Tshikolovets 2021 syn. nov.) as a distinct species and designate a lectotype; synonymize Aporia leucodice leucodice Eversmann, 1843 (= A. l. morosevitshae Sheljuzhko, 1908 syn. nov.); and describe a new species, Aporia ahura sp. nov., from the Central Alborz Mountains in northern Iran.

Within-host adaptive speciation of commensal yoyo clams leads to ecological exclusion, not co-existence.

Symbionts dominate planetary diversity and three primary symbiont diversification processes have been proposed: co-speciation with hosts, speciation by host-switching, and within-host speciation. The last mechanism is prevalent among members of an extraordinary marine symbiosis in the Indian River Lagoon, Florida, composed of a host mantis shrimp, Lysiosquilla scabricauda, and seven host-specific commensal vasconielline "yoyo" clams (Galeommatoidea) that collectively occupy two distinct niches: burrow-wall-attached, and host-attached/ectocommensal. This within-host symbiont radiation provides a natural experiment to test how symbiont coexistence patterns are regulated in a common ancestral habitat. The competitive exclusion principle predicts that sister taxa produced by adaptive speciation (with distinct morphologies and within-burrow niches) are most likely to coexist whereas the neutral theory predicts no difference among adaptive and non-adaptive sister taxa co-occurrence. To test these predictions, we engaged in (1) field-censusing commensal species assemblages; (2) trophic niche analyses; (3) laboratory behavioral observations. Although predicted by both models, the field census found no mixed-niche commensal assemblages: multi-species burrows were exclusively composed of burrow-wall commensals. Their co-occurrence matched random assembly process expectations, but presence of the single ectocommensal species had a highly significant negative effect on recruitment of all burrow-wall commensal species (P < 0.001), including on its burrow-wall commensal sister species (P < 0.001). Our stable isotope data indicated that commensals are suspension feeders and that co-occurring burrow-wall commensals may exhibit trophic niche differentiation. The artificial burrow behavioral experiment yielded no evidence of spatial segregation among burrow-wall commensals, and it was terminated by a sudden breakdown of the host-commensal relationship resulting in a mass mortality of all commensals unattached to the host. This study system appears to contain two distinct, superimposed patterns of commensal distribution: (1) all burrow-wall commensal species; (2) the ectocommensal species. Burrow-wall commensals (the plesiomorphic condition) broadly adhere to neutral theory expectations of species assembly but the adaptive evolution of ectocommensalism has apparently led to ecological exclusion rather than coexistence, an inverse outcome of theoretical expectations. The ecological factors regulating the observed burrow-wall/ectocommensal exclusion are currently obscure but potentially include differential recruitment to host burrows and/or differential survival in "mixed" burrow assemblages, the latter potentially due to changes in host predatory behavior. Resampling host burrows during commensal recruitment peak periods and tracking burrow-wall commensal survival in host burrows with and without added ectocommensals could resolve this outstanding issue.

Phylogenomic analyses of all species of swordtail fishes (genus Xiphophorus) show that hybridization preceded speciation.

Hybridization has been recognized to play important roles in evolution, however studies of the genetic consequence are still lagging behind in vertebrates due to the lack of appropriate experimental systems. Fish of the genus Xiphophorus are proposed to have evolved with multiple ancient and ongoing hybridization events. They have served as an informative research model in evolutionary biology and in biomedical research on human disease for more than a century. Here, we provide the complete genomic resource including annotations for all described 26 Xiphophorus species and three undescribed taxa and resolve all uncertain phylogenetic relationships. We investigate the molecular evolution of genes related to cancers such as melanoma and for the genetic control of puberty timing, focusing on genes that are predicted to be involved in pre-and postzygotic isolation and thus affect hybridization. We discovered dramatic size-variation of some gene families. These persisted despite reticulate evolution, rapid speciation and short divergence time. Finally, we clarify the hybridization history in the entire genus settling disputed hybridization history of two Southern swordtails. Our comparative genomic analyses revealed hybridization ancestries that are manifested in the mosaic fused genomes and show that hybridization often preceded speciation.

Beyond Haldane's rule: Sex-biased hybrid dysfunction for all modes of sex determination.

Haldane's rule occupies a special place in biology as one of the few 'rules' of speciation, with empirical support from hundreds of species. And yet, its classic purview is restricted taxonomically to the subset of organisms with heteromorphic sex chromosomes. I propose explicit acknowledgement of generalized hypotheses about Haldane's rule that frame sex bias in hybrid dysfunction broadly and irrespective of the sexual system. The consensus view of classic Haldane's rule holds that sex-biased hybrid dysfunction across taxa is a composite phenomenon that requires explanations from multiple causes. Testing of the multiple alternative hypotheses for Haldane's rule is, in many cases, applicable to taxa with homomorphic sex chromosomes, environmental sex determination, haplodiploidy, and hermaphroditism. Integration of a variety of biological phenomena about hybrids across diverse sexual systems, beyond classic Haldane's rule, will help to derive a more general understanding of the contributing forces and mechanisms that lead to predictable sex biases in evolutionary divergence and speciation.

A consequential one-night stand: Episodic historical hybridization leads to mitochondrial takeover in sympatric desert ant-eating spiders.

Disentangling the genomic intricacies underlying speciation and the causes of discordance between sources of evidence can offer remarkable insights into evolutionary dynamics. The ant-eating spider Zodarion nitidum, found across the Middle East and Egypt, displays yellowish and blackish morphs that co-occur sympatrically. These morphs additionally differ in behavioral and physiological features and show complete pre-mating reproductive isolation. In contrast, they possess similar sexual features and lack distinct differences in their mitochondrial DNA. We analyzed both Z. nitidum morphs and outgroups using genome-wide and additional mitochondrial DNA data. The genomic evidence indicated that Yellow and Black are reciprocally independent lineages without signs of recent admixture. Interestingly, the sister group of Yellow is not Black but Z. luctuosum, a morphologically distinct species. Genomic gene flow analyses pinpointed an asymmetric nuclear introgression event, with Yellow contributing nearly 5 % of its genome to Black roughly 320,000 years ago, intriguingly aligning with the independently estimated origin of the mitochondrial DNA of Black. We conclude that the blackish and yellowish morphs of Z. nitidum are long-diverged distinct species, and that the ancient and modest genomic introgression event registered resulted in a complete mitochondrial takeover of Black by Yellow. This investigation underscores the profound long-term effects that even modest hybridization events can have on the genome of organisms. It also exemplifies the utility of phylogenetic networks for estimating historical events and how integrating independent lines of evidence can increase the reliability of such estimations.

Establishing species boundaries in Bornean geckos.

Species delimitation using mitochondrial DNA (mtDNA) remains an important and accessible approach for discovering and delimiting species. However, delimiting species with a single locus (e.g. DNA barcoding) is biased towards overestimating species diversity. The highly diverse gecko genus Cyrtodactylus is one such group where delimitation using mtDNA remains the paradigm. In this study, we use genomic data to test putative species boundaries established using mtDNA within three recognized species of Cyrtodactylus on the island of Borneo. We predict that multi-locus genomic data will estimate fewer species than mtDNA, which could have important ramifications for the species diversity within the genus. We aim to (i) investigate the correspondence between species delimitations using mtDNA and genomic data, (ii) infer species trees for each target species, and (iii) quantify gene flow and identify migration patterns to assess population connectivity. We find that species diversity is overestimated and that species boundaries differ between mtDNA and nuclear data. This underscores the value of using genomic data to reassess mtDNA-based species delimitations for taxa lacking clear species boundaries. We expect the number of recognized species within Cyrtodactylus to continue increasing, but, when possible, genomic data should be included to inform more accurate species boundaries.

Origin of the Laurentian Great Lakes fish fauna through upward adaptive radiation cascade prior to the Last Glacial Maximum.

The evolutionary histories of adaptive radiations can be marked by dramatic demographic fluctuations. However, the demographic histories of ecologically-linked co-diversifying lineages remain understudied. The Laurentian Great Lakes provide a unique system of two such radiations that are dispersed across depth gradients with a predator-prey relationship. We show that the North American Coregonus species complex ("ciscoes") radiated rapidly prior to the Last Glacial Maximum (80-90 kya), a globally warm period, followed by rapid expansion in population size. Similar patterns of demographic expansion were observed in the predator species, Lake Charr (Salvelinus namaycush), following a brief time lag, which we hypothesize was driven by predator-prey dynamics. Diversification of prey into deep water created ecological opportunities for the predators, facilitating their demographic expansion, which is consistent with an upward adaptive radiation cascade. This study provides a new timeline and environmental context for the origin of the Laurentian Great Lakes fish fauna, and firmly establishes this system as drivers of ecological diversification and rapid speciation through cyclical glaciation.

Trait-mediated speciation and human-driven extinctions in proboscideans revealed by unsupervised Bayesian neural networks.

Species life-history traits, paleoenvironment, and biotic interactions likely influence speciation and extinction rates, affecting species richness over time. Birth-death models inferring the impact of these factors typically assume monotonic relationships between single predictors and rates, limiting our ability to assess more complex effects and their relative importance and interaction. We introduce a Bayesian birth-death model using unsupervised neural networks to explore multifactorial and nonlinear effects on speciation and extinction rates using fossil data. It infers lineage- and time-specific rates and disentangles predictor effects and importance through explainable artificial intelligence techniques. Analysis of the proboscidean fossil record revealed speciation rates shaped by dietary flexibility and biogeographic events. The emergence of modern humans escalated extinction rates, causing recent diversity decline, while regional climate had a lesser impact. Our model paves the way for an improved understanding of the intricate dynamics shaping clade diversification.

The genetic basis of incipient sexual isolation in Drosophila melanogaster.

Speciation is a fundamental evolutionary process but the genetic changes accompanying speciation are difficult to determine since true species do not produce viable and fertile offspring. Partially reproductively isolated incipient species are useful for assessing genetic changes that occur prior to speciation. Drosophila melanogaster from Zimbabwe, Africa are partially sexually isolated from other D. melanogaster populations whose males have poor mating success with Zimbabwe females. We used the North American D. melanogaster Genetic Reference Panel (DGRP) to show that there is significant genetic variation in mating success of DGRP males with Zimbabwe females, to map genetic variants and genes associated with variation in mating success and to determine whether mating success to Zimbabwe females is associated with other quantitative traits previously measured in the DGRP. Incipient sexual isolation is highly polygenic and associated with the common African inversion In(3R)K and the amount of the sex pheromone 5,9-heptacosadiene in DGRP females. We functionally validated the effect of eight candidate genes using RNA interference to provide testable hypotheses for future studies investigating the molecular genetic basis of incipient sexual isolation in D. melanogaster.

Evolution of secondary metabolites, morphological structures and associated gene expression patterns in galls induced by four closely related aphid species on a host plant species.

Gall-forming insects induce various types of galls on their host plants by altering gene expression in host plant organs, and recent studies have been conducted for gene expression in galls. However, the evolutionary trajectories of gene expression patterns and the resulting phenotypes have not yet been studied using multiple related species. We investigated the speciation and the diversification process of galls induced by four closely related aphid species (Hormaphidini) on a host plant species (Hamamelis japonica) by examining the phylogenetic congruence between the geographical divergences of aphids and the host plant, and by comparing their gene expression patterns and resulting phenotypes. Phylogenetic analysis of aphids and the host plant showed that geographical isolation among host plant populations has interrupted gene flow in aphids and accelerated the speciation process. The concentration of phenolics and the complexity of the internal structure of galls were correlated with the expression levels of genes for the biosynthesis of phenolics and morphogenesis respectively. These results suggest that the expression levels of genes for the biosynthesis of phenolics and morphogenesis have evolutionarily increased in galls accelerated by the speciation process of aphids due to the distribution change of the host plant, leading to the related phenotypic evolution. Our study showed the evolutionary process of phenotypic traits in galls in the wild from both gene expression and actual phenotype levels.

A Diffusion-Based Approach for Simulating Forward-in-Time State-Dependent Speciation and Extinction Dynamics.

We establish a general framework using a diffusion approximation to simulate forward-in-time state counts or frequencies for cladogenetic state-dependent speciation-extinction (ClaSSE) models. We apply the framework to various two- and three-region geographic-state speciation-extinction (GeoSSE) models. We show that the species range state dynamics simulated under tree-based and diffusion-based processes are comparable. We derive a method to infer rate parameters that are compatible with given observed stationary state frequencies and obtain an analytical result to compute stationary state frequencies for a given set of rate parameters. We also describe a procedure to find the time to reach the stationary frequencies of a ClaSSE model using our diffusion-based approach, which we demonstrate using a worked example for a two-region GeoSSE model. Finally, we discuss how the diffusion framework can be applied to formalize relationships between evolutionary patterns and processes under state-dependent diversification scenarios.

Phylogeography and paleoclimatic range dynamics explain variable outcomes to contact across a species' range.

Replicability of divergence after contact is a poorly characterized process, particularly in the contexts of phylogeography and postglacial range dynamics within species. Using contact zones located at the leading-, mid- and rear-edges of a species' range, we examined variation in outcomes to contact between divergent lineages of Campanula americana. We investigated whether contact zones vary in quantity and directionality of gene flow, how phylogeographic structure differs between contact zones, and how historic range dynamics may affect outcomes to contact. We found that all contact zones formed at similar times via primary contact yet detected significant admixture in only the rear-edge (RE) contact zone. In the northern leading-edge contact zone and the mid-range Virginia contact zone, gene flow was minimal and asymmetric. In the southern RE contact zone, gene flow was strong and symmetric. Asymmetric admixture in the leading-edge and Virginia contact zones matches the directionality of a known cosmopolitan cytonuclear incompatibility between lineages of C. americana. Our results emphasize the dependence of speciation processes on phylogeographic structure, evolutionary history and range dynamics.

The biogeographic and evolutionary processes shaping population divergence in Laupala.

Speciation generates biodiversity and the mechanisms involved are thought to vary across the tree of life and across environments. For example, well-studied adaptive radiations are thought to be fuelled by divergent ecological selection, but additionally are influenced heavily by biogeographic, genomic and demographic factors. Mechanisms of non-adaptive radiations, producing ecologically cryptic taxa, have been less well-studied but should likewise be influenced by these latter factors. Comparing among contexts can help pinpoint universal mechanisms and outcomes, especially if we integrate biogeographic, ecological and evolutionary processes. We investigate population divergence in the swordtail cricket Laupala cerasina, a wide-spread endemic on Hawai'i Island and one of 38 ecologically cryptic Laupala species. The nine sampled populations show striking population genetic structure at small spatio-temporal scales. The rapid differentiation among populations and species of Laupala shows that neither a specific geographical context nor ecological opportunity are pre-requisites for rapid divergence. Spatio-temporal patterns in population divergence, population size change, and gene flow are aligned with the chronosequence of the four volcanoes on which L. cerasina occurs and reveal the composite effects of geological dynamics and Quaternary climate change on population dynamics. Spatio-temporal patterns in genetic variation along the genome reveal the interplay of genetic and genomic architecture in shaping population divergence. In early phases of divergence, we find elevated differentiation in genomic regions harbouring mating song loci. In later stages of divergence, we find a signature of linked selection that interacts with recombination rate variation. Comparing our findings with recent work on complementary systems supports the conclusion that mostly universal factors influence the speciation process.

Evolution: Spectral speciation.

Speciation is a complex process sparked by multitudes of environmental stressors and culminating in adaptive, and perhaps novel, phenotypic traits. A new study presents evidence supporting spectral niche-partitioning in a cyanobacterial clade specializing in far-red photosynthesis.

Underground speciation: Unraveling the systematics and evolution of the highly diverse tuco-tucos (genus Ctenomys) with genomic data.

Subterranean rodents of the genus Ctenomys (tuco-tucos) are endemic to South America and have experienced relatively recent radiation. There are about 67 recognized species that originated in approximately 1-2 MY. They stand out for their species richness, extraordinary chromosomal diversity, and wide range of habitat they occupy in the continent. Phylogenetic relationships among species of tuco-tucos have been challenging to resolve. Groups of closely-related species have been suggested, but their relationships must be resolved. This study estimates the phylogeny of the genus using massive sequencing, generating thousands of independent molecular markers obtained by RADseq, with a taxonomic sampling that includes 66% of the recognized species. The sequences obtained were mapped against the C. sociabilis genome, recovering up to 1,215 widely shared RAD loci with more than 19,000 polymorphic sites. Our new phylogenetic hypothesis corroborated the species groups previously proposed with cytochrome b gene sequences and provided a much greater resolution of the relationships among species groups. The frater group is sister to all other tuco-tucos, whereas some of the earlierliest proposals placed the sociabilis group as sister to all other tuco-tucos. Ctenomys leucodon, previously proposed as an independent lineage, is associated with the frater group with moderate statistical support. The magellanicus and mendocinus are sister groups in a major clade formed by the boliviensis, talarum, tucumanus, torquatus, and opimus groups. Ctenomys viperinus, included in the phylogeny for the first time, belongs to the tucumanus group. This multi-locus phylogenetic hypothesis provides insights into the historical biogeography of understanding this highly diverse genus.

Identifiability of Level-1 Species Networks from Gene Tree Quartets.

When hybridization or other forms of lateral gene transfer have occurred, evolutionary relationships of species are better represented by phylogenetic networks than by trees. While inference of such networks remains challenging, several recently proposed methods are based on quartet concordance factors-the probabilities that a tree relating a gene sampled from the species displays the possible 4-taxon relationships. Building on earlier results, we investigate what level-1 network features are identifiable from concordance factors under the network multispecies coalescent model. We obtain results on both topological features of the network, and numerical parameters, uncovering a number of failures of identifiability related to 3-cycles in the network. Addressing these identifiability issues is essential for designing statistically consistent inference methods.

Phylogenomic insights into the historical biogeography, character-state evolution, and species diversification rates of Cypripedioideae (Orchidaceae).

Cypripedioideae (slipper orchids; Orchidaceae) currently consist of âˆ¼200 herbaceous species with a strikingly disjunctive distribution in tropical and temperate regions of both hemispheres. In this study, an updated phylogeny with representatives from all five cypripedioid genera was presented based on maximum likelihood and Bayesian inference of plastome and low-copy nuclear genes. Phylogenomic analyses indicated that each genus is monophyletic, but some relationships (e.g., those among Cypripedium sects. Acaulia, Arietinum, Bifolia, Flabellinervia, Obtusipetala and Palangshanensia) conflict with those in previous studies based on Sanger data. Cypripedioideae appeared to have arisen in South America and/or the adjacent Qinghai-Tibet Plateau and Hengduan Mountains âˆ¼35 Mya. We inferred multiple dispersal events between East Asia and North America in Cypripedium, and between mainland Southeast Asia and the Malay Archipelago in Paphiopedilum. In the Americas, divergences among four genera (except Cypripedium) occurred around 31-20 Mya, long before the closure of the Isthmus of Panama, indicating the importance of long-distance dispersal. Evolutionary patterns between morphological and plastome character evolution suggested several traits, genome size and NDH genes, which are likely to have contributed to the success of slipper orchids in alpine floras and low-elevation forests. Species diversification rates were notably higher in epiphytic clades of Paphiopedilum than in other, terrestrial cypripedioids, paralleling similar accelerations associated with epiphytism in other groups. This study also suggested that sea-level fluctuations and mountain-building processes promoted the diversification of the largest genera, Paphiopedilum and Cypripedium.

The impact of geographic isolation and host shifts on population divergence of the rare cicada Subpsaltria yangi.

The contributions of divergent selection and spatial isolation to population divergence are among the main focuses of evolutionary biology. Here we employed integrated methods to explore genomic divergence, demographic history and calling-song differentiation in the cicada Subpsaltria yangi, and compared the genotype and calling-song phenotype of different populations occurring in distinct habitats. Our results indicate that this species comprises four main lineages with unique sets of haplotypes and calling-song structure, which are distinctly associated with geographic isolation and habitats. The populations occurring on the Loess Plateau underwent substantial expansion at âˆ¼0.130-0.115 Ma during the Last Interglacial. Geographic distance and host shift between pairs of populations predict genomic divergence, with geographic distance and acoustical signal together explaining > 60% of the divergence among populations. Differences in calling songs could reflect adaptation of populations to novel environments with different host plants, habitats and predators, which may have resulted from neutral divergence at the molecular level followed by natural selection. Geomorphic barriers and climate oscillations associated with Pleistocene glaciation may have been primary factors in shaping the population genetic structure of this species. Ultimately this may couple with a host shift in leading toward allopatric speciation in S. yangi, i.e., isolation by distance. Our findings improve understanding of divergence in allopatry of herbivorous insects, and may inform future studies on the molecular mechanisms underlying the association between genetic/phenotypic changes and adaptation of insects to novel niches and host plants.

The centromere landscapes of four karyotypically diverse Papaver species provide insights into chromosome evolution and speciation.

Understanding the roles played by centromeres in chromosome evolution and speciation is complicated by the fact that centromeres comprise large arrays of tandemly repeated satellite DNA, which hinders high-quality assembly. Here, we used long-read sequencing to generate nearly complete genome assemblies for four karyotypically diverse Papaver species, P. setigerum (2n = 44), P. somniferum (2n = 22), P. rhoeas (2n = 14), and P. bracteatum (2n = 14), collectively representing 45 gapless centromeres. We identified four centromere satellite (cenSat) families and experimentally validated two representatives. For the two allopolyploid genomes (P. somniferum and P. setigerum), we characterized the subgenomic distribution of each satellite and identified a "homogenizing" phase of centromere evolution in the aftermath of hybridization. An interspecies comparison of the peri-centromeric regions further revealed extensive centromere-mediated chromosome rearrangements. Taking these results together, we propose a model for studying cenSat competition after hybridization and shed further light on the complex role of the centromere in speciation.

Evidence for gene flow and trait reversal during radiation of Mexican Goodeid fish.

Understanding the phylogeographic history of a group and identifying the factors contributing to speciation is an important challenge in evolutionary biology. The Goodeinae are a group of live-bearing fishes endemic to Mexico. Here, we develop genomic resources for species within the Goodeinae and use phylogenomic approaches to characterise their evolutionary history. We sequenced, assembled and annotated the genomes of four Goodeinae species, including Ataeniobius toweri, the only matrotrophic live-bearing fish without a trophotaenia in the group. We estimated timings of species divergence and examined the extent and timing of introgression between the species to assess if this may have occurred during an early radiation, or in more recent episodes of secondary contact. We used branch-site models to detect genome-wide positive selection across Goodeinae, and we specifically asked whether this differs in A. toweri, where loss of placental viviparity has recently occurred. We found evidence of gene flow between geographically isolated species, suggesting vicariant speciation was supplemented by limited post-speciation gene flow, and gene flow may explain previous uncertainties about Goodeid phylogeny. Genes under positive selection in the group are likely to be associated with the switch to live-bearing. Overall, our studies suggest that both volcanism-driven vicariance and changes in reproductive mode influenced radiation in the Goodeinae.