Genomics and the origin of species

Ploidy as a leaky reproductive barrier: mechanisms, rates and evolutionary significance of interploidy gene flow

BACKGROUND

Whole-genome duplication (polyploidization) is a dominant force in sympatric speciation, particularly in plants. Genome doubling instantly poses a barrier to gene flow owing to the strong crossing incompatibilities between individuals differing in ploidy. The strength of the barrier, however, varies from species to species and recent genetic investigations revealed cases of rampant interploidy introgression in multiple ploidy-variable species.

SCOPE

Here, we review novel insights into the frequency of interploidy gene flow in natural systems and summarize the underlying mechanisms promoting interploidy gene flow. Field surveys, occasionally complemented by crossing experiments, suggest frequent opportunities for interploidy gene flow, particularly in the direction from diploid to tetraploid, and between (higher) polyploids. However, a scarcity of accompanying population genetic evidence and a virtual lack of integration of these approaches leave the underlying mechanisms and levels of realized interploidy gene flow in nature largely unknown. Finally, we discuss potential consequences of interploidy genome permeability on polyploid speciation and adaptation and highlight novel avenues that have just recently been opened by the very first genomic studies of ploidy-variable species. Standing in stark contrast with rapidly accumulating evidence for evolutionary importance of homoploid introgression, similar cases in ploidy-variable systems are yet to be documented.

CONCLUSIONS

The genomics era provides novel opportunity to re-evaluate the role of interploidy introgression in speciation and adaptation. To achieve this goal, interdisciplinary studies bordering ecology and population genetics and genomics are needed.

Whole-chromosome oligo-painting in licorice unveils interspecific chromosomal evolutionary relationships and possible origin of triploid genome species

Licorice is one of the most extensively studied medicinal plants in the world, whose roots and rhizomes have long been used as both a sweetener and an essential component in numerous herbal preparations. However, the genus Glycyrrhiza has a complex composition, and the interspecies chromosomal relationships, origin, and evolution are still largely unclear. Here, we develop a set of whole-chromosome painting probes that allowed identification of all eight chromosomes of licorice on same metaphase chromosomes. Comparative chromosome painting analyses in seven different Glycyrrhiza species revealed that the genus Glycyrrhiza maintained extraordinarily conserved chromosomal synteny after about 3-12 million years of divergence. No cytologically visible inter-chromosomal rearrangements were identified in any species. By comparative chromosomal karyotype analyses, we revealed interspecific chromosome evolutionary relationships and dramatic variable chromosomal karyotype after independent divergence and demonstrated that G. prostrate was the most closely related to the ancestral type among the seven Glycyrrhiza species. Furthermore, we also discovered a G. glandulosa seed with distinct triploid-genome for the first time in China, suggesting the existence of a polyploid evolutionary pathway in the genus Glycyrrhiza, which challenges the previous notion that only diploids of licorice existed in nature. This study expands our knowledge of the chromosome evolution of licorice and will lay an important foundation for the genome origin and evolution studies in the genus Glycyrrhiza.

Integrative analyses of convergent adaptation in sympatric extremophile fishes

The evolution of independent lineages along replicated environmental transitions frequently results in convergent adaptation, yet the degree to which convergence is present across multiple levels of biological organization is often unclear. Additionally, inherent biases associated with shared ancestry and variation in selective regimes across geographic replicates often pose challenges for confidently identifying patterns of convergence. We investigated a system in which three species of poeciliid fishes sympatrically occur in a toxic spring rich in hydrogen sulfide (H2S) and an adjacent nonsulfidic stream to examine patterns of adaptive evolution across levels of biological organization. We found convergence in morphological and physiological traits and genome-wide patterns of gene expression among all three species. In addition, there were shared signatures of selection on genes encoding H2S toxicity targets in the mitochondrial genomes of each species. However, analyses of nuclear genomes revealed neither evidence for substantial genomic islands of divergence around genes involved in H2S toxicity and detoxification nor substantial congruence of strongly differentiated regions across population pairs. These non-convergent, heterogeneous patterns of genomic divergence may indicate that sulfide tolerance is highly polygenic, with shared allele frequency shifts present at many loci with small effects along the genome. Alternatively, H2S tolerance may involve substantial genetic redundancy, with non-convergent, lineage-specific variation at multiple loci along the genome underpinning similar changes in phenotypes and gene expression. Overall, we demonstrate variability in the extent of convergence across organizational levels and highlight the challenges of linking patterns of convergence across scales.

Evolution of Sensory Receptors

Sensory receptors are at the interface between an organism and its environment and thus represent key sites for biological innovation. Here, we survey major sensory receptor families to uncover emerging evolutionary patterns. Receptors for touch, temperature, and light constitute part of the ancestral sensory toolkit of animals, often predating the evolution of multicellularity and the nervous system. In contrast, chemoreceptors exhibit a dynamic history of lineage-specific expansions and contractions correlated with the disparate complexity of chemical environments. A recurring theme includes independent transitions from neurotransmitter receptors to sensory receptors of diverse stimuli from the outside world. We then provide an overview of the evolutionary mechanisms underlying sensory receptor diversification and highlight examples where signatures of natural selection are used to identify novel sensory adaptations. Finally, we discuss sensory receptors as evolutionary hotspots driving reproductive isolation and speciation, thereby contributing to the stunning diversity of animals.

On the evolution and genetic diversity of the bread wheat D genome

Bread wheat (Triticum aestivum) became a globally dominant crop after incorporating the D genome from the donor species Aegilops tauschii, but the evolutionary history that shaped the D genome during this process remains to be clarified. Here, we propose a renewed evolutionary model linking Ae. tauschii and the hexaploid wheat D genome by constructing an ancestral haplotype map covering 762 Ae. tauschii and hexaploid wheat accessions. We dissected the evolutionary trajectories of Ae. tauschii lineages and reported a few independent intermediate accessions, demonstrating that low-frequency inter-sublineage gene flow had enriched the diversity of Ae. tauschii. We discovered that the D genome of hexaploid wheat was inherited from a unified ancestral template, but with a mosaic composition that was highly mixed and derived mainly from three Ae. tauschii L2 sublineages located in the Caspian coastal region. This result suggests that early agricultural activities facilitated innovations in D-genome composition and finalized the success of hexaploidization. We found that the majority (51.4%) of genetic diversity was attributed to novel mutations absent in Ae. tauschii, and we identified large Ae. tauschii introgressions from various lineages, which expanded the diversity of the wheat D genome and introduced beneficial alleles. This work sheds light on the process of wheat hexaploidization and highlights the evolutionary significance of the multi-layered genetic diversity of the bread wheat D genome.

Integrative taxonomy of Metarhizium anisopliae species complex, based on phylogenomics combined with morphometrics, metabolomics, and virulence data

Metarhizium anisopliae (Clavicipitaceae, Hypocreales) is a globally distributed entomopathogenic fungus, which has been largely studied and used in agriculture for its potent entomopathogenicity. Since its taxonomic establishment as a member of Metarhizium, many closely related taxa have been described with highly similar morphology (cryptic species). A holotype specimen of M. anisopliae is not extant, and the ex-neotype strain (CBS 130.71) does not form a monophyletic clade with other strains, up to now, recognized as M. anisopliae sensu stricto. In this study, we have conducted an integrative taxonomic treatment of M. anisopliae sensu lato by including the ex-neotype strain of M. anisopliae, other unknown strains from our collections identified as M. anisopliae s. lat., as well as other known species that have been previously delimited as closely related but distinct to M. anisopliae. By including whole-genome sequencing, morphometric analysis, LC-MS based metabolomics, and virulence assays, we have demonstrated that M. anisopliae s. str. should also include M. lepidiotae (synonym), and that M. anisopliae s. str. differentiates from the other species of the complex by its metabolome and less severe entomopathogenicity. New taxa, namely M. hybridum, M. neoanisopliae and M. parapingshaense spp. nov., are proposed. The novel taxa proposed here have strong phylogenomics support, corroborated by fine-scale differences in the length/width of conidia/phialides, while the metabolomics and virulence data still largely overlap. We have also demonstrated via population genomics data the existence of local clonal lineages, particularly the one corresponding to the persistence of a biocontrol candidate strain that has been used in the field application for three years. This study showcases the utility of combining various data sources for accurate delimitation of species within an important group of fungal biocontrol agents against pest insects.

Combining Molecular, Macroevolutionary, and Macroecological Perspectives on the Generation of Diversity

Charles Darwin presented a unified process of diversification driven by the gradual accumulation of heritable variation. The growth in DNA databases and the increase in genomic sequencing, combined with advances in molecular phylogenetic analyses, gives us an opportunity to realize Darwin's vision, connecting the generation of variation to the diversification of lineages. The rate of molecular evolution is correlated with the rate of diversification across animals and plants, but the relationship between genome change and speciation is complex: Mutation rates evolve in response to life history and niche; substitution rates are influenced by mutation, selection, and population size; rates of acquisition of reproductive isolation vary between populations; and traits, niches, and distribution can influence diversification rates. The connection between mutation rate and diversification rate is one part of the complex and varied story of speciation, which has theoretical importance for understanding the generation of biodiversity and also practical impacts on the use of DNA to understand the dynamics of speciation over macroevolutionary timescales.

Beyond gene flow: (non)-parallelism of secondary contact in a pair of highly differentiated sibling species

Replicated secondary contact zones can provide insights into the barriers to gene flow that are important during speciation and can reveal to which degree secondary contact may result in similar evolutionary outcomes. Here, we studied two secondary contact zones between highly differentiated Alpine butterflies of the genus Erebia using whole-genome resequencing data. We assessed the genomic relationships between populations and species and found hybridization to be rare, with no to little current or historical introgression in either contact zone. There are large similarities between contact zones, consistent with an allopatric origin of interspecific differentiation, with no indications for ongoing reinforcing selection. Consistent with expected reduced effective population size, we further find that scaffolds related to the Z-chromosome show increased differentiation compared to the already high levels across the entire genome, which could also hint towards a contribution of the Z chromosome to species divergence in this system. Finally, we detected the presence of the endosymbiont Wolbachia, which can cause reproductive isolation between its hosts, in all E. cassioides, while it appears to be fully or largely absent in contact zone populations of E. tyndarus. We discuss how this rare pattern may have arisen and how it may have affected the dynamics of speciation upon secondary contact.

Comparative genomics of dusky kob (Argyrosomus japonicus, Sciaenidae) conspecifics: Evidence for speciation and the genetic mechanisms underlying traits

Dusky kob (Argyrosomus japonicus) is a commercially important finfish, indigenous to South Africa, Australia, and China. Previous studies highlighted differences in genetic composition, life history, and morphology of the species across geographic regions. A draft genome sequence of 0.742 Gb (N50 = 5.49 Mb; BUSCO completeness = 97.8%) and 22,438 predicted protein-coding genes was generated for the South African (SA) conspecific. A comparison with the Chinese (CN) conspecific revealed a core set of 32,068 orthologous protein clusters across both genomes. The SA genome exhibited 440 unique clusters compared to 1928 unique clusters in the CN genome. Transportation and immune response processes were overrepresented among the SA accessory genome, whereas the CN accessory genome was enriched for immune response, DNA transposition, and sensory detection (FDR-adjusted p < 0.01). These unique clusters may represent an adaptive component of the species' pangenome that could explain population divergence due to differential environmental specialisation. Furthermore, 700 single-copy orthologues (SCOs) displayed evidence of positive selection between the SA and CN genomes, and globally these genomes shared only 92% similarity, suggesting they might be distinct species. These genes primarily play roles in metabolism and digestion, illustrating the evolutionary pathways that differentiate the species. Understanding these genomic mechanisms underlying adaptation and evolution within and between species provides valuable insights into growth and maturation of kob, traits that are particularly relevant to commercial aquaculture.

The importance of reproductive isolation in driving diversification and speciation within Peruvian mimetic poison frogs (Dendrobatidae)

To explain how populations with distinct warning signals coexist in close parapatry, we experimentally assessed intrinsic mechanisms acting as reproductive barriers within three poison-frog species from the Peruvian Amazon belonging to a Müllerian mimetic ring (Ranitomeya variabilis, Ranitomeya imitator and Ranitomeya fantastica). We tested the role of prezygotic and postzygotic isolation barriers between phenotypically different ecotypes of each species, using no-choice mating experiments and offspring survival analysis. Our results show that prezygotic mating preference did not occur except for one specific ecotype of R. imitator, and that all three species were able to produce viable inter-population F1 hybrids. However, while R. variabilis and R. imitator hybrids were able to produce viable F2 generations, we found that for R. fantastica, every F1 hybrid males were sterile while females remained fertile. This unexpected result, echoing with Haldane's rule of speciation, validated phylogenetic studies which tentatively diagnose these populations of R. fantastica as two different species. Our work suggests that postzygotic genetic barriers likely participate in the extraordinary phenotypic diversity observed within Müllerian mimetic Ranitomeya populations, by maintaining species boundaries.

Coupling of Barriers to Gene Exchange: Causes and Consequences

Coupling has emerged as a concept to describe the transition from differentiated populations to newly evolved species through the strengthening of reproductive isolation. However, the term has been used in multiple ways, and relevant processes have sometimes not been clearly distinguished. Here, we synthesize existing uses of the concept of coupling and find three main perspectives: (1) coupling as the build-up of linkage disequilibrium among loci underlying barriers to gene exchange, (2) coupling as the build-up of genome-wide linkage disequilibrium, and (3) coupling as the process generating a coincidence of distinct barrier effects. We compare and contrast these views, show the diverse processes involved and the complexity of the relationships among recombination, linkage disequilibrium, and reproductive isolation, and, finally, we emphasize how each perspective can guide new directions in speciation research. Although the importance of coupling for evolutionary divergence and speciation is well established, many theoretical and empirical questions remain unanswered.

Historical and contemporary processes driving the origin and structure of an admixed population within a contact zone between subspecies of a north temperate diadromous fish

Hybridization between divergent lineages can result in losses of distinct evolutionary taxa. Alternatively, hybridization can lead to increased genetic variability that may fuel local adaptation and the generation of novel traits and/or taxa. Here, we examined single-nucleotide polymorphisms generated using genotyping-by-sequencing in a population of Dolly Varden char (Pisces: Salmonidae) that is highly admixed within a contact zone between two subspecies (Salvelinus malma malma, Northern Dolly Varden [NDV] and S. m. lordi, Southern Dolly Varden [SDV]) in southwestern Alaska to assess the spatial distribution of hybrids and to test hypotheses on the origin of the admixed population. Ancestry analysis revealed that this admixed population is composed of advanced generation hybrids between NDV and SDV or advanced backcrosses to SDV; no F1 hybrids were detected. Coalescent-based demographic modelling supported the origin of this population about 55,000 years ago by secondary contact between NDV and SDV with low levels of contemporary gene flow. Ancestry in NDV and SDV varies within the watershed and ancestry in NDV was positively associated with distance upstream from the sea, contingent on habitat-type sampled, and negatively associated with the number of migrations that individual fish made to the sea. Our results suggest that divergence between subspecies over hundreds of thousands of years may not be associated with significant reproductive isolation, but that elevated diversity owing to hybridization may have contributed to adaptive divergence in habitat use and life history.

Phenotypic characterization and gene mapping of hybrid necrosis in Triticum durum-Haynaldia villosa amphiploids

KEY MESSAGE

Phenotypical, physiological and genetic characterization was carried out on the hybrid necrosis gene from Haynaldia villosa, and the related gene Ne-V was mapped to chromosome arm 2VL. Introducing genetic variation from wild relatives into common wheat through wide crosses is a vital strategy for enriching genetic diversity and promoting wheat breeding. However, hybrid necrosis, a genetic autoimmunity syndrome, often occurs in the offspring of interspecific or intraspecific crosses, restricting both the selection of hybrid parents and the pyramiding of beneficial genes. To utilize the germplasms of Haynaldia villosa (2n = 2x = 14, VV), we conducted wide hybridization between durum wheat (2n = 4x = 28, AABB) and multiple H. villosa accessions to synthesize the amphiploids (2n = 6x = 42, AABBVV). This study revealed that 61.5% of amphiploids derived from the above crosses exhibited hybrid necrosis, with some amphiploids even dying before reaching maturity. However, the initiation time and severity of necrosis varied dramatically among the progenies, suggesting that there were multiple genetic loci or multiple alleles in the same genetic locus conferring to hybrid necrosis in H. villosa accessions. Genetic analysis was performed on the F2 and derived F2:3 populations, which were constructed between amphiploid STH59-1 with normal leaves and amphiploid STH59-2 with necrotic leaves. A semidominant hybrid necrosis-related gene, Ne-V, was mapped to an 11.8-cM genetic interval on the long arm of chromosome 2V, representing a novel genetic locus identified in Triticum-related species. In addition, the hybrid necrosis was correlated with enhanced H2O2 accumulation and cell death, and it was influenced by the temperature and light. Our findings provide a foundation for cloning the Ne-V gene and exploring its molecular mechanism.

Sexual selection and speciation in the Anthropocene

Anthropogenic change threatens global biodiversity by causing severe ecological disturbance and extinction. Here, we consider the effects of anthropogenic change on one process that generates biodiversity. Sexual selection (a potent evolutionary force and driver of speciation) is highly sensitive to the environment and, thus, vulnerable to anthropogenic ecological change. Anthropogenic alterations to sexual display and mate preference can make it harder to distinguish between conspecific and heterospecific mates or can weaken divergence via sexual selection, leading to higher rates of hybridization and biodiversity loss. Occasionally, anthropogenically altered sexual selection can abet diversification, but this appears less likely than biodiversity loss. In our rapidly changing world, a full understanding of sexual selection and speciation requires a global change perspective.

Molecular and agro-morphological characterization of new barley genotypes in arid environments

BACKGROUND

Genetic diversity, population structure, agro-morphological traits, and molecular characteristics, are crucial for either preserving genetic resources or developing new cultivars. Due to climate change, water availability for agricultural use is progressively diminishing. This study used 100 molecular markers (25 TRAP, 22 SRAP, 23 ISTR, and 30 SSR). Additionally, 15 morphological characteristics were utilized to evaluate the optimal agronomic traits of 12 different barley genotypes under arid conditions.

RESULTS

Substantial variations, ranging from significant to highly significant, were observed in the 15 agromorphological parameters evaluated among the 12 genotypes. The KSU-B101 barley genotype demonstrated superior performance in five specific traits: spike number per plant, 100-grain weight, spike number per square meter, harvest index, and grain yield. These results indicate its potential for achieving high yields in arid regions. The Sahrawy barley genotype exhibited the highest values across five parameters, namely leaf area, spike weight per plant, spike length, spike weight per square meter, and biological yield, making it a promising candidate for animal feed. The KSU-B105 genotype exhibited early maturity and a high grain count per spike, which reflects its early maturity and ability to produce a high number of grains per spike. This suggests its suitability for both animal feed and human food in arid areas. Based on marker data, the molecular study found that the similarity coefficients between the barley genotypes ranged from 0.48 to 0.80, with an average of 0.64. The dendrogram constructed from these data revealed three distinct clusters with a similarity coefficient of 0.80. Notably, the correlation between the dendrogram and its similarity matrix was high (0.903), indicating its accuracy in depicting the genetic relationships. The combined analysis revealed a moderate correlation between the morphological and molecular analysis, suggesting alignment between the two characterization methods.

CONCLUSIONS

The morphological and molecular analyses of the 12 barley genotypes in this study effectively revealed the varied genetic characteristics of their agro-performance in arid conditions. KSU-B101, Sahrawy, and KSU-B105 have emerged as promising candidates for different agricultural applications in arid regions. Further research on these genotypes could reveal their full potential for breeding programs.

Genomic landscapes of divergence among island bird populations: Evidence of parallel adaptation but at different loci?

When populations colonise new environments, they may be exposed to novel selection pressures but also suffer from extensive genetic drift due to founder effects, small population sizes and limited interpopulation gene flow. Genomic approaches enable us to study how these factors drive divergence, and disentangle neutral effects from differentiation at specific loci due to selection. Here, we investigate patterns of genetic diversity and divergence using whole-genome resequencing (>22× coverage) in Berthelot's pipit (Anthus berthelotii), a passerine endemic to the islands of three north Atlantic archipelagos. Strong environmental gradients, including in pathogen pressure, across populations in the species range, make it an excellent system in which to explore traits important in adaptation and/or incipient speciation. First, we quantify how genomic divergence accumulates across the speciation continuum, that is, among Berthelot's pipit populations, between sub species across archipelagos, and between Berthelot's pipit and its mainland ancestor, the tawny pipit (Anthus campestris). Across these colonisation timeframes (2.1 million-ca. 8000 years ago), we identify highly differentiated loci within genomic islands of divergence and conclude that the observed distributions align with expectations for non-neutral divergence. Characteristic signatures of selection are identified in loci associated with craniofacial/bone and eye development, metabolism and immune response between population comparisons. Interestingly, we find limited evidence for repeated divergence of the same loci across the colonisation range but do identify different loci putatively associated with the same biological traits in different populations, likely due to parallel adaptation. Incipient speciation across these island populations, in which founder effects and selective pressures are strong, may therefore be repeatedly associated with morphology, metabolism and immune defence.

Recurrent gene flow events occurred during the diversification of clownfishes of the skunk complex

Clownfish (subfamily Amphiprioninae) are an iconic group of coral reef fish that evolved a mutualistic interaction with sea anemones, which triggered the adaptive radiation of the clade. Within clownfishes, the "skunk complex" is particularly interesting. Besides ecological speciation, interspecific gene flow and hybrid speciation are thought to have shaped the evolution of the group. We investigated the mechanisms characterizing the diversification of this complex. By taking advantage of their disjunct geographical distribution, we obtained whole-genome data of sympatric and allopatric populations of the three main species of the complex (Amphiprion akallopisos, A. perideraion and A. sandaracinos). We examined population structure, genomic divergence and introgression signals and performed demographic modelling to identify the most realistic diversification scenario. We excluded scenarios of strict isolation or hybrid origin of A. sandaracinos. We discovered moderate gene flow from A. perideraion to the ancestor of A. akallopisos + A. sandaracinos and weak gene flow between the species in the Indo-Australian Archipelago throughout the diversification of the group. We identified introgressed regions in A. sandaracinos and detected in A. perideraion two large regions of high divergence from the two other species. While we found that gene flow has occurred throughout the species' diversification, we also observed that recent admixture was less pervasive than initially thought, suggesting a role of host repartition or behavioural barriers in maintaining the genetic identity of the species in sympatry.

Chromosomal inversions from an initial ecotypic divergence drive a gradual repeated radiation of Galápagos beetles

Island faunas exhibit some of the most iconic examples where similar forms repeatedly evolve within different islands. Yet, whether these deterministic evolutionary trajectories within islands are driven by an initial, singular divergence and the subsequent exchange of individuals and adaptive genetic variation between islands remains unclear. Here, we study a gradual, repeated evolution of low-dispersive highland ecotypes from a dispersive lowland ecotype of Calosoma beetles along the island progression of the Galápagos. We show that repeated highland adaptation involved selection on multiple shared alleles within extensive chromosomal inversions that originated from an initial adaptation event on the oldest island. These highland inversions first spread through dispersal of highland individuals. Subsequent admixture with the lowland ecotype resulted in polymorphic dispersive populations from which the highland populations evolved on the youngest islands. Our findings emphasize the significance of an ancient divergence in driving repeated evolution and highlight how a mixed contribution of inter-island colonization and within-island evolution can shape parallel species communities.

Tangled Banks, Braided Rivers and Complex Hierarchies: Beyond Microevolution and Macroevolution

Ever since the Modern Synthesis, a debate about the relationship between microevolution and macroevolution has persisted - specifically, whether they are equivalent, distinct, or explain one another. How one answers this has become shorthand for a much broader set of theoretical debates in evolutionary biology. Here, we examine microevolution and macroevolution in the context of the vast proliferation of data, knowledge, and theory since the advent of the Modern Synthesis. We suggest that traditional views on microevolution and macroevolution are too binary and reductive. For example, patterns and processes are not confined to micro- and macro- domains; they are interconnected at various temporal and spatial scales and across hierarchical entities. Further, biological entities have variably fuzzy boundaries, and evolutionary processes that influence macroevolution occur at micro- and macro- levels. In addition, these conceptual advances in phylodynamics have yet to be fully integrated with contemporary macroevolutionary approaches. Finally, holding microevolution and macroevolution as distinct domains thwarts synthesis and collaboration on important research questions. We propose that the focal entities and processes considered by evolutionary studies be contextualized within the newfound complexity of the multidimensional, multi-modal, multi-level phylogenetic system.

Genome sequencing and comparative genomics reveal insights into pathogenicity and evolution of Fusarium zanthoxyli, the causal agent of stem canker in prickly ash

BACKGROUND

Fusarium zanthoxyli is a destructive pathogen causing stem canker in prickly ash, an ecologically and economically important forest tree. However, the genome lack of F. zanthoxyli has hindered research on its interaction with prickly ash and the development of precise control strategies for stem canker.

RESULTS

In this study, we sequenced and annotated a relatively high-quality genome of F. zanthoxyli with a size of 43.39 Mb, encoding 11,316 putative genes. Pathogenicity-related factors are predicted, comprising 495 CAZymes, 217 effectors, 156 CYP450s, and 202 enzymes associated with secondary metabolism. Besides, a comparative genomics analysis revealed Fusarium and Colletotrichum diverged from a shared ancestor approximately 141.1 ~ 88.4 million years ago (MYA). Additionally, a phylogenomic investigation of 12 different phytopathogens within Fusarium indicated that F. zanthoxyli originated approximately 34.6 ~ 26.9 MYA, and events of gene expansion and contraction within them were also unveiled. Finally, utilizing conserved domain prediction, the results revealed that among the 59 unique genes, the most enriched domains were PnbA and ULP1. Among the 783 expanded genes, the most enriched domains were PKc_like kinases and those belonging to the APH_ChoK_Like family.

CONCLUSION

This study sheds light on the genetic basis of F. zanthoxyli's pathogenicity and evolution which provides valuable information for future research on its molecular interactions with prickly ash and the development of effective strategies to combat stem canker.

Comparative genomics of macaques and integrated insights into genetic variation and population history

The crab-eating macaques ( Macaca fascicularis ) and rhesus macaques ( M. mulatta ) are widely studied nonhuman primates in biomedical and evolutionary research. Despite their significance, the current understanding of the complex genomic structure in macaques and the differences between species requires substantial improvement. Here, we present a complete genome assembly of a crab-eating macaque and 20 haplotype-resolved macaque assemblies to investigate the complex regions and major genomic differences between species. Segmental duplication in macaques is ∼42% lower, while centromeres are ∼3.7 times longer than those in humans. The characterization of ∼2 Mbp fixed genetic variants and ∼240 Mbp complex loci highlights potential associations with metabolic differences between the two macaque species (e.g., CYP2C76 and EHBP1L1 ). Additionally, hundreds of alternative splicing differences show post-transcriptional regulation divergence between these two species (e.g., PNPO ). We also characterize 91 large-scale genomic differences between macaques and humans at a single-base-pair resolution and highlight their impact on gene regulation in primate evolution (e.g., FOLH1 and PIEZO2 ). Finally, population genetics recapitulates macaque speciation and selective sweeps, highlighting potential genetic basis of reproduction and tail phenotype differences (e.g., STAB1 , SEMA3F , and HOXD13 ). In summary, the integrated analysis of genetic variation and population genetics in macaques greatly enhances our comprehension of lineage-specific phenotypes, adaptation, and primate evolution, thereby improving their biomedical applications in human diseases.

Diverse pathways to speciation revealed by marine snails

Speciation is a key evolutionary process that is not yet fully understood. Combining population genomic and ecological data from multiple diverging pairs of marine snails (Littorina) supports the search for speciation mechanisms. Placing pairs on a one-dimensional speciation continuum, from undifferentiated populations to species, obscured the complexity of speciation. Adding multiple axes helped to describe either speciation routes or reproductive isolation in the snails. Divergent ecological selection repeatedly generated barriers between ecotypes, but appeared less important in completing speciation while genetic incompatibilities played a key role. Chromosomal inversions contributed to genomic barriers, but with variable impact. A multidimensional (hypercube) approach supported framing of questions and identification of knowledge gaps and can be useful to understand speciation in many other systems.

Divergent dynamics of sexual and habitat isolation at the transition between stick insect populations and species

Speciation is often viewed as a continuum along which populations diverge until they become reproductively-isolated species. However, such divergence may be heterogeneous, proceeding in fits and bursts, rather than being uniform and gradual. We show in Timema stick insects that one component of reproductive isolation evolves non-uniformly across this continuum, whereas another does not. Specifically, we use thousands of host-preference and mating trials to study habitat and sexual isolation among 42 pairs of taxa spanning a range of genomic differentiation and divergence time. We find that habitat isolation is uncoupled from genomic differentiation within species, but accumulates linearly with it between species. In contrast, sexual isolation accumulates linearly across the speciation continuum, and thus exhibits similar dynamics to morphological traits not implicated in reproductive isolation. The results show different evolutionary dynamics for different components of reproductive isolation and highlight a special relevance for species status in the process of speciation.

Genetic architecture of ecological divergence between Oryza rufipogon and Oryza nivara

Ecological divergence due to habitat difference plays a prominent role in the formation of new species, but the genetic architecture during ecological speciation and the mechanism underlying phenotypic divergence remain less understood. Two wild ancestors of rice (Oryza rufipogon and Oryza nivara) are a progenitor-derivative species pair with ecological divergence and provide a unique system for studying ecological adaptation/speciation. Here, we constructed a high-resolution linkage map and conducted a quantitative trait locus (QTL) analysis of 19 phenotypic traits using an F2 population generated from a cross between the two Oryza species. We identified 113 QTLs associated with interspecific divergence of 16 quantitative traits, with effect sizes ranging from 1.61% to 34.1% in terms of the percentage of variation explained (PVE). The distribution of effect sizes of QTLs followed a negative exponential, suggesting that a few genes of large effect and many genes of small effect were responsible for the phenotypic divergence. We observed 18 clusters of QTLs (QTL hotspots) on 11 chromosomes, significantly more than that expected by chance, demonstrating the importance of coinheritance of loci/genes in ecological adaptation/speciation. Analysis of effect direction and v-test statistics revealed that interspecific differentiation of most traits was driven by divergent natural selection, supporting the argument that ecological adaptation/speciation would proceed rapidly under coordinated selection on multiple traits. Our findings provide new insights into the understanding of genetic architecture of ecological adaptation and speciation in plants and help effective manipulation of specific genes or gene cluster in rice breeding.

The effect of divergent and parallel selection on the genomic landscape of divergence

While the role of selection in divergence along the speciation continuum is theoretically well understood, defining specific signatures of selection in the genomic landscape of divergence is empirically challenging. Modelling approaches can provide insight into the potential role of selection on the emergence of a heterogenous genomic landscape of divergence. Here, we extend and apply an individual-based approach that simulates the phenotypic and genotypic distributions of two populations under a variety of selection regimes, genotype-phenotype maps, modes of migration, and genotype-environment interactions. We show that genomic islands of high differentiation and genomic valleys of similarity may respectively form under divergent and parallel selection between populations. For both types of between-population selection, negative and positive frequency-dependent selection within populations generated genomic islands of higher magnitude and genomic valleys of similarity, respectively. Divergence rates decreased under strong dominance with divergent selection, as well as in models including genotype-environment interactions under parallel selection. For both divergent and parallel selection models, divergence rate was higher under an intermittent migration regime between populations, in contrast to a constant level of migration across generations, despite an equal number of total migrants. We highlight that interpreting a particular evolutionary history from an observed genomic pattern must be done cautiously, as similar patterns may be obtained from different combinations of evolutionary processes. Modelling approaches such as ours provide an opportunity to narrow the potential routes that generate the genomic patterns of specific evolutionary histories.

Epigenetic and Genetic Differentiation Between Coregonus Species Pairs

Phenotypic diversification is classically associated with genetic differentiation and gene expression variation. However, increasing evidence suggests that DNA methylation is involved in evolutionary processes due to its phenotypic and transcriptional effects. Methylation can increase mutagenesis and could lead to increased genetic divergence between populations experiencing different environmental conditions for many generations, though there has been minimal empirical research on epigenetically induced mutagenesis in diversification and speciation. Whitefish, freshwater members of the salmonid family, are excellent systems to study phenotypic diversification and speciation due to the repeated divergence of benthic-limnetic species pairs serving as natural replicates. Here we investigate whole genome genetic and epigenetic differentiation between sympatric benthic-limnetic species pairs in lake and European whitefish (Coregonus clupeaformis and Coregonus lavaretus) from four lakes (N = 64). We found considerable, albeit variable, genetic and epigenetic differences between species pairs. All SNP types were enriched at CpG sites supporting the mutagenic nature of DNA methylation, though C>T SNPs were most common. We also found an enrichment of overlaps between outlier SNPs with the 5% highest FST between species and differentially methylated loci. This could possibly represent differentially methylated sites that have caused divergent genetic mutations between species, or divergent selection leading to both genetic and epigenetic variation at these sites. Our results support the hypothesis that DNA methylation contributes to phenotypic divergence and mutagenesis during whitefish speciation.

Bayesian hybrid index and genomic cline estimation with the R package gghybrid

Admixture, the interbreeding of individuals from differentiated source populations, is now known to be a widespread phenomenon. Genomic studies of natural hybridisation can help to answer many questions on the impacts of admixture on adaptive evolution, reproductive isolation, and speciation. When a large variety of admixture proportions between two source populations exist, both geographic and genomic cline analysis are suitable methods for inferring biased, restricted or excessive gene flow at individual loci into the foreign genomic background, providing evidence for reproductive isolation, selection across an environmental transition, balancing selection, and adaptive introgression. Genomic cline analysis replaces geographic location with genome-wide hybrid index and is therefore useable in circumstances that violate geographic cline assumptions. Here, I introduce gghybrid, an R package for simple and flexible Bayesian estimation of Buerkle's hybrid index and Fitzpatrick's logit-logistic genomic clines using bi-allelic data, suitable for both small and large datasets. gghybrid allows any ploidy and uses Structure input file format. It has separate functions for hybrid index and cline estimation, treating each individual and locus respectively as an independent analysis, making it highly parallelisable. Admixture proportions from other software can alternatively be used in cline analysis, alongside parental allele frequencies. Parameters can be fixed and samples pooled for statistical model comparison with AIC or waic. Here, I describe the functions, pipeline, and statistical properties of gghybrid. Simulations reveal that model comparison with waic is preferred, and use of Bayesian posterior distributions and p values to select candidate non-null loci is problematic and should be avoided.

Patterns of admixture and introgression in a mosaic Mytilus galloprovincialis and Mytilus edulis hybrid zone in SW England

The study of hybrid zones offers important insights into speciation. Earlier studies on hybrid populations of the marine mussel species Mytilus edulis and Mytilus galloprovincialis in SW England provided evidence of admixture but were constrained by the limited number of molecular markers available. We use 57 ancestry-informative SNPs, most of which have been mapped genetically, to provide evidence of distinctive differences between admixed populations in SW England and asymmetrical introgression from M. edulis to M. galloprovincialis. We combine the genetic study with analysis of phenotypic traits of potential ecological and adaptive significance. We demonstrate that hybrid individuals have brown mantle edges unlike the white or purple in the parental species, suggesting allelic or non-allelic genomic interactions. We report differences in gonad development stage between the species consistent with a prezygotic barrier between the species. By incorporating results from publications dating back to 1980, we confirm the long-term stability of the hybrid zone despite higher viability of M. galloprovincialis. This stability coincides with a dramatic change in temperature of UK coastal waters and suggests that these hybrid populations might be resisting the effects of global warming. However, a single SNP locus associated with the Notch transmembrane signalling protein shows a markedly different pattern of variation to the others and might be associated with adaptation of M. galloprovincialis to colder northern temperatures.

Multiple contact zones and karyotypic evolution in a neotropical frog species complex

Previous studies of DNA sequence and karyotypic data have revealed high genetic diversity in the Physalaemus cuvieri - Physalaemus ephippifer species complex-a group of small leptodactylid frogs in South America. To date, seven major genetic lineages have been recognized in this group, with species delimitation tests supporting four to seven of them as valid species. Among these, only P. ephippifer shows heteromorphic sex chromosomes, but the implications of cytogenetic divergence for the evolution of this group are unknown. We analyzed karyotypic, mitochondrial DNA, and 3RAD genomic data to characterize a putative contact zone between P. ephippifer and P. cuvieri Lineage 1, finding evidence for admixture and karyotypic evolution. We also describe preliminary evidence for admixture between two other members of this species complex-Lineage 1 and Lineage 3 of P. cuvieri. Our study sheds new light on evolutionary relationships in the P. cuvieri - P. ephippifer species complex, suggesting an important role of karyotypic divergence in its evolutionary history and underscoring the importance of hybridization as a mechanism of sex chromosome evolution in amphibians.

Evidence that genetic drift not adaptation drives fast-Z and large-Z effects in Ficedula flycatchers

The sex chromosomes have been hypothesized to play a key role in driving adaptation and speciation across many taxa. The reason for this is thought to be the hemizygosity of the heteromorphic part of sex chromosomes in the heterogametic sex, which exposes recessive mutations to natural and sexual selection. The exposure of recessive beneficial mutations increases their rate of fixation on the sex chromosomes, which results in a faster rate of evolution. In addition, genetic incompatibilities between sex-linked loci are exposed faster in the genomic background of hybrids of divergent lineages, which makes sex chromosomes contribute disproportionately to reproductive isolation. However, in birds, which show a Z/W sex determination system, the role of adaptation versus genetic drift as the driving force of the faster differentiation of the Z chromosome (fast-Z effect) and the disproportionate role of the Z chromosome in reproductive isolation (large-Z effect) are still debated. Here, we address this debate in the bird genus Ficedula flycatchers based on population-level whole-genome sequencing data of six species. Our analysis provides evidence for both faster lineage sorting and reduced gene flow on the Z chromosome than the autosomes. However, these patterns appear to be driven primarily by the increased role of genetic drift on the Z chromosome, rather than an increased rate of adaptive evolution. Genomic scans of selective sweeps and fixed differences in fact suggest a reduced action of positive selection on the Z chromosome.

Genomic Architecture Predicts Tree Topology, Population Structuring, and Demographic History in Amazonian Birds

Geographic barriers are frequently invoked to explain genetic structuring across the landscape. However, inferences on the spatial and temporal origins of population variation have been largely limited to evolutionary neutral models, ignoring the potential role of natural selection and intrinsic genomic processes known as genomic architecture in producing heterogeneity in differentiation across the genome. To test how variation in genomic characteristics (e.g. recombination rate) impacts our ability to reconstruct general patterns of differentiation between species that cooccur across geographic barriers, we sequenced the whole genomes of multiple bird populations that are distributed across rivers in southeastern Amazonia. We found that phylogenetic relationships within species and demographic parameters varied across the genome in predictable ways. Genetic diversity was positively associated with recombination rate and negatively associated with species tree support. Gene flow was less pervasive in genomic regions of low recombination, making these windows more likely to retain patterns of population structuring that matched the species tree. We further found that approximately a third of the genome showed evidence of selective sweeps and linked selection, skewing genome-wide estimates of effective population sizes and gene flow between populations toward lower values. In sum, we showed that the effects of intrinsic genomic characteristics and selection can be disentangled from neutral processes to elucidate spatial patterns of population differentiation.

Peripatric speciation within Torreya fargesii (Taxaceae) in the Hengduan Mountains inferred from multi-loci phylogeography

BACKGROUND

The Hengduan Mountains (HDM) are one of the major global biodiversity hotspots in the world. Several evolutionary scenarios, especially in-situ diversification, have been proposed to account for the high species richness of temperate plants. However, peripatric speciation, an important mode of allopatric speciation, has seldom been reported in this region.

RESULTS

Here, two chloroplast DNA regions and 14 nuclear loci were sequenced for 112 individuals from 10 populations of Torreya fargesii var. fargesii and 63 individuals from 6 populations of T. fargesii var. yunnanensis. Population genetic analyses revealed that the two varieties are well differentiated genetically (FST, 0.5765) and have uneven genetic diversity (π, 0.00221 vs. 0.00073 on an average of nuclear loci). The gene genealogical relationship showed that T. fargesii var. yunnanensis is inferred as derived from T. fargesii var. fargesii, which was further supported by the coalescent simulations (DIYABC, fastsimcoal2 and IMa2). By the coalescent simulations, the divergence time (~ 2.50-3.65 Ma) and the weak gene flow between the two varieties were detected. The gene flow was asymmetrical and only occurred in later stages of divergence, which is caused by second contact due to the population expansion (~ 0.61 Ma) in T. fargesii var. fargesii. In addition, niche modeling indicated that the two varieties are differentiated geographically and ecologically and have unbalanced distribution range.

CONCLUSIONS

Overall, T. fargesii var. fargesii is always parapatric with respect to T. fargesii var. yunnanensis, and the latter derived from the former in peripatry of the HDM following a colonization from central China during the late Pliocene. Our findings demonstrate that peripatric speciation following dispersal events may be an important evolutionary scenario for the formation of biodiversity hotspot of the HDM.

New chromosome-scale genomes provide insights into marine adaptations of sea snakes (Hydrophis: Elapidae)

BACKGROUND

Sea snakes underwent a complete transition from land to sea within the last ~ 15 million years, yet they remain a conspicuous gap in molecular studies of marine adaptation in vertebrates.

RESULTS

Here, we generate four new annotated sea snake genomes, three of these at chromosome-scale (Hydrophis major, H. ornatus and H. curtus), and perform detailed comparative genomic analyses of sea snakes and their closest terrestrial relatives. Phylogenomic analyses highlight the possibility of near-simultaneous speciation at the root of Hydrophis, and synteny maps show intra-chromosomal variations that will be important targets for future adaptation and speciation genomic studies of this system. We then used a strict screen for positive selection in sea snakes (against a background of seven terrestrial snake genomes) to identify genes over-represented in hypoxia adaptation, sensory perception, immune response and morphological development.

CONCLUSIONS

We provide the best reference genomes currently available for the prolific and medically important elapid snake radiation. Our analyses highlight the phylogenetic complexity and conserved genome structure within Hydrophis. Positively selected marine-associated genes provide promising candidates for future, functional studies linking genetic signatures to the marine phenotypes of sea snakes and other vertebrates.

Diversity across organisational scale emerges through dispersal ability and speciation dynamics in tropical fish

BACKGROUND

Biodiversity exists at different levels of organisation: e.g. genetic, individual, population, species, and community. These levels of organisation all exist within the same system, with diversity patterns emerging across organisational scales through several key processes. Despite this inherent interconnectivity, observational studies reveal that diversity patterns across levels are not consistent and the underlying mechanisms for variable continuity in diversity across levels remain elusive. To investigate these mechanisms, we apply a spatially explicit simulation model to simulate the global diversification of tropical reef fishes at both the population and species levels through emergent population-level processes.

RESULTS

We find significant relationships between the population and species levels of diversity which vary depending on both the measure of diversity and the spatial partitioning considered. In turn, these population-species relationships are driven by modelled biological trait parameters, especially the divergence threshold at which populations speciate.

CONCLUSIONS

To explain variation in multi-level diversity patterns, we propose a simple, yet novel, population-to-species diversity partitioning mechanism through speciation which disrupts continuous diversity patterns across organisational levels. We expect that in real-world systems this mechanism is driven by the molecular dynamics that determine genetic incompatibility, and therefore reproductive isolation between individuals. We put forward a framework in which the mechanisms underlying patterns of diversity across organisational levels are universal, and through this show how variable patterns of diversity can emerge through organisational scale.

Quantitative trait loci underlying a speciation phenotype

Sexual signalling traits and their associated genetic components play a crucial role in the speciation process, as divergence in these traits can contribute to sexual isolation. Despite their importance, our understanding of the genetic basis of variable sexual signalling traits linked to speciation remains limited. In this study, we present new genetic evidence of Quantitative Trait Loci (QTL) underlying divergent sexual signalling behaviour, specifically pulse rate, in the Hawaiian cricket Laupala. By performing RNA sequencing on the brain and central nervous system of the parental species, we annotate these QTL regions and identify candidate genes associated with pulse rate. Our findings provide insights into the genetic processes driving reproductive isolation during speciation, with implications for understanding the mechanisms underlying species diversity.

Partial and asymmetrical reproductive isolation between two sympatric tropical shrub species: Cnidoscolus aconitifolius and C. souzae (Euphorbiaceae)

Reproductive isolation is conferred by several barriers that occur at different stages of reproduction. Comprehensive reviews on the topic have identified that barriers occurring prior to zygote formation are often stronger than those that occur afterward. However, the overrepresentation of temperate perennial herbs in the current literature precludes any generalization of this pattern to plants that present other life forms and patterns of distribution. Here, we assessed reproductive isolation barriers and their absolute contribution to reproductive isolation and asymmetry in Cnidoscolus aconitifolius and C. souzae, two closely related tropical shrub species that co-occur on the Yucatan peninsula. The reproductive barriers assessed were phenological mismatch, pollinator differentiation, pollen-pistil incompatibility (three pre-zygotic barriers), fruit set failure, and seed unviability (post-zygotic barriers). Reproductive isolation between the study species was found to be complete in the direction C. aconitifolius to C. souzae, but only partial in the opposite direction. One post-zygotic barrier was the strongest example. Most barriers, particularly the pre-zygotic examples, were asymmetrical and predicted the direction of heterospecific pollen flow and hybrid formation from C. souzae to C. aconitifolius. Both parental species, as well as the hybrids, were diploid and had a chromosome number 2n = 36. More studies with tropical woody perennials are required to fully determine whether this group of plants consistently shows stronger post-zygotic barriers.

Darwin's 'mystery of mysteries': the role of sexual selection in plant speciation

Sexual selection is considered one of the key processes that contribute to the emergence of new species. While the connection between sexual selection and speciation has been supported by comparative studies, the mechanisms that mediate this connection remain unresolved, especially in plants. Similarly, it is not clear how speciation processes within plant populations translate into large-scale speciation dynamics. Here, we review the mechanisms through which sexual selection, pollination, and mate choice unfold and interact, and how they may ultimately produce reproductive isolation in plants. We also overview reproductive strategies that might influence sexual selection in plants and illustrate how functional traits might connect speciation at the population level (population differentiation, evolution of reproductive barriers; i.e. microevolution) with evolution above the species level (macroevolution). We also identify outstanding questions in the field, and suitable data and tools for their resolution. Altogether, this effort motivates further research focused on plants, which might potentially broaden our general understanding of speciation by sexual selection, a major concept in evolutionary biology.

Genetic divergences and hybridization within the Sebastes inermis complex

The Sebastes inermis complex includes three sympatric species (Sebastes cheni, viz Sebastes inermis, and Sebastes ventricosus) with clear ecomorphological differences, albeit incomplete reproductive isolation. The presence of putative morphological hybrids (PMH) with plausibly higher fitness than the parent species indicates the need to confirm whether hybridization occurs within the complex. In this sense, we assessed the dynamics of genetic divergence and hybridization within the species complex using a panel of 10 microsatellite loci, and sequences of the mitochondrial control region (D-loop) and the intron-free rhodopsin (RH1) gene. The analyses revealed the presence of three distinct genetic clusters, large genetic distances using D-loop sequences, and distinctive mutations within the RH1 gene. These results are consistent with the descriptions of the three species. Two microsatellite loci had signatures of divergent selection, indicating that they are linked to genomic regions that are crucial for speciation. Furthermore, nonsynonymous mutations within the RH1 gene detected in S. cheni and "Kumano" (a PMH) suggest dissimilar adaptations related to visual perception in dim-light environments. The presence of individuals with admixed ancestry between two species confirmed hybridization. The presence of nonsynonymous mutations within the RH1 gene and the admixed ancestry of the "Kumano" morphotype highlight the potential role of hybridization in generating novelties within the species complex. We discuss possible outcomes of hybridization within the species complex, considering hybrid fitness and assortative mating. Overall, our findings indicate that the genetic divergence of each species is maintained in the presence of hybridization, as expected in a scenario of speciation-with-gene-flow.

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.

Sperm-dependent asexual species and their role in ecology and evolution

Sexual reproduction is the primary mode of reproduction in eukaryotes, but some organisms have evolved deviations from classical sex and switched to asexuality. These asexual lineages have sometimes been viewed as evolutionary dead ends, but recent research has revealed their importance in many areas of general biology. Our review explores the understudied, yet important mechanisms by which sperm-dependent asexuals that produce non-recombined gametes but rely on their fertilization, can have a significant impact on the evolution of coexisting sexual species and ecosystems. These impacts are concentrated around three major fields. Firstly, sperm-dependent asexuals can potentially impact the gene pool of coexisting sexual species by either restricting their population sizes or by providing bridges for interspecific gene flow whose type and consequences substantially differ from gene flow mechanisms expected under sexual reproduction. Secondly, they may impact on sexuals' diversification rates either directly, by serving as stepping-stones in speciation, or indirectly, by promoting the formation of pre- and postzygotic reproduction barriers among nascent species. Thirdly, they can potentially impact on spatial distribution of species, via direct or indirect (apparent) types of competition and Allee effects. For each such mechanism, we provide empirical examples of how natural sperm-dependent asexuals impact the evolution of their sexual counterparts. In particular, we highlight that these broad effects may last beyond the tenure of the individual asexual lineages causing them, which challenges the traditional perception that asexual lineages are short-lived evolutionary dead ends and minor sideshows. Our review also proposes new research directions to incorporate the aforementioned impacts of sperm-dependent asexuals. These research directions will ultimately enhance our understanding of the evolution of genomes and biological interactions in general.

Can demographic histories explain long-term isolation and recent pulses of asymmetric gene flow between highly divergent grey fox lineages?

Secondary contact zones between deeply divergent, yet interfertile, lineages provide windows into the speciation process. North American grey foxes (Urocyon cinereoargenteus) are divided into western and eastern lineages that diverged approximately 1 million years ago. These ancient lineages currently hybridize in a relatively narrow zone of contact in the southern Great Plains, a pattern more commonly observed in smaller-bodied taxa, which suggests relatively recent contact after a long period of allopatry. Based on local ancestry inference with whole-genome sequencing (n = 43), we identified two distinct Holocene pulses of admixture. The older pulse (500-3500 YBP) reflected unidirectional gene flow from east to west, whereas the more recent pulse (70-200 YBP) of admixture was bi-directional. Augmented with genotyping-by-sequencing data from 216 additional foxes, demographic analyses indicated that the eastern lineage declined precipitously after divergence, remaining small throughout most of the late Pleistocene, and expanding only during the Holocene. Genetic diversity in the eastern lineage was highest in the southeast and lowest near the contact zone, consistent with a westward expansion. Concordantly, distribution modelling indicated that during their isolation, the most suitable habitat occurred far east of today's contact zone or west of the Great Plains. Thus, long-term isolation was likely caused by the small, distant location of the eastern refugium, with recent contact reflecting a large increase in suitable habitat and corresponding demographic expansion from the eastern refugium. Ultimately, long-term isolation in grey foxes may reflect their specialized bio-climatic niche. This system presents an opportunity for future investigation of potential pre- and post-zygotic isolating mechanisms.

Temporal Variation in Selection Influences Microgeographic Local Adaptation

AbstractEcological heterogeneity can lead to local adaptation when populations exhibit fitness trade-offs among habitats. However, the degree to which local adaptation is affected by the spatial and temporal scale of environmental variation is poorly understood. A multiyear reciprocal transplant experiment was performed with populations of the annual plant Leptosiphon parviflorus living on adjacent serpentine and nonserpentine soil. Local adaptation over this small geographic scale was observed, but there were differences in the temporal variability of selection across habitats. On serpentine soil, the local population had a consistently large survival advantage, presumably as a result of the temporal stability in selection imposed by soil cation content. In contrast, a fecundity advantage was observed for the sandstone population on its native soil type but only in the two study years with the highest rainfall. A manipulative greenhouse experiment demonstrated that the fitness advantage of the sandstone population in its native soil type depends critically on water availability. The temporal variability in local adaptation driven by variation in precipitation suggests that continued drought conditions have the potential to erode local adaptation in these populations. These results show how different selective factors can influence spatial and temporal patterns of variation in fitness trade-offs.

Community-wide genome sequencing reveals 30 years of Darwin's finch evolution

A fundamental goal in evolutionary biology is to understand the genetic architecture of adaptive traits. Using whole-genome data of 3955 of Darwin's finches on the Galápagos Island of Daphne Major, we identified six loci of large effect that explain 45% of the variation in the highly heritable beak size of Geospiza fortis, a key ecological trait. The major locus is a supergene comprising four genes. Abrupt changes in allele frequencies at the loci accompanied a strong change in beak size caused by natural selection during a drought. A gradual change in Geospiza scandens occurred across 30 years as a result of introgressive hybridization with G. fortis. This study shows how a few loci with large effect on a fitness-related trait contribute to the genetic potential for rapid adaptive radiation.

Genomics of plant speciation

Studies of plants have been instrumental for revealing how new species originate. For several decades, botanical research has complemented and, in some cases, challenged concepts on speciation developed via the study of other organisms while also revealing additional ways in which species can form. Now, the ability to sequence genomes at an unprecedented pace and scale has allowed biologists to settle decades-long debates and tackle other emerging challenges in speciation research. Here, we review these recent genome-enabled developments in plant speciation. We discuss complications related to identification of reproductive isolation (RI) loci using analyses of the landscape of genomic divergence and highlight the important role that structural variants have in speciation, as increasingly revealed by new sequencing technologies. Further, we review how genomics has advanced what we know of some routes to new species formation, like hybridization or whole-genome duplication, while casting doubt on others, like population bottlenecks and genetic drift. While genomics can fast-track identification of genes and mutations that confer RI, we emphasize that follow-up molecular and field experiments remain critical. Nonetheless, genomics has clarified the outsized role of ancient variants rather than new mutations, particularly early during speciation. We conclude by highlighting promising avenues of future study. These include expanding what we know so far about the role of epigenetic and structural changes during speciation, broadening the scope and taxonomic breadth of plant speciation genomics studies, and synthesizing information from extensive genomic data that have already been generated by the plant speciation community.

New species based on the biological species concept within the complex of Lariophagus distinguendus (Hymenoptera, Chalcidoidea, Pteromalidae), a parasitoid of household pests

The pteromalid parasitoid Lariophagus distinguendus (Foerster) belongs to the Hymenoptera, a megadiverse insect order with high cryptic diversity. It attacks stored product pest beetles in human storage facilities. Recently, it has been shown to consist of two separate species. To further study its cryptic diversity, strains were collected to compare their relatedness using barcoding and nuclear genes. Nuclear genes identified two clusters which agree with the known two species, whereas the barcode fragment determined an additional third Clade. Total reproductive isolation (RI) according to the biological species concept (BSC) was investigated in crossing experiments within and between clusters using representative strains. Sexual isolation exists between all studied pairs, increasing from slight to strong with genetic distance. Postzygotic barriers mostly affected hybrid males, pointing to Haldane's rule. Hybrid females were only affected by unidirectional Spiroplasma-induced cytoplasmic incompatibility and behavioural sterility, each in one specific strain combination. RI was virtually absent between strains separated by up to 2.8% COI difference, but strong or complete in three pairs from one Clade each, separated by at least 7.2%. Apparently, each of these clusters represents one separate species according to the BSC, highlighting cryptic diversity in direct vicinity to humans. In addition, these results challenge the recent 'turbo-taxonomy' practice of using 2% COI differences to delimitate species, especially within parasitic Hymenoptera. The gradual increase in number and strength of reproductive barriers between strains with increasing genetic distance also sheds light on the emergence of barriers during the speciation process in L. distinguendus.

Different patterns of introgression in a three species hybrid zone among European cave salamanders

Hybrid zones occur where genetically distinct populations meet, mate and produce offspring with mixed ancestry. In Plethodontid salamanders, introgressive hybridization is a common phenomenon, where hybrids backcross with parental populations leading to the spread of new alleles into the parental genomes. Whereas many hybrid zones have been reported in American Plethodontid salamanders, only a single hybrid zone has been documented in European plethodontids so far, which is located at the Apuan Alps in the Italian Peninsula. Here, we describe a previously unreported hybrid zone in the Northern Apennines involving all the three Plethodontid salamander species inhabiting the Italian Peninsula. We found 21 new Speleomantes sites of occurrence, from a hitherto unexplored area located at the boundaries between three Speleomantes species ranges. Using mitochondrial (Cytb and ND2 genes) and nuclear markers (two diagnostic SNPs at the NCX1 gene), we revealed a three-way contact zone where all the three mainland species hybridize: S. strinatii, S. ambrosii and S. italicus. We observed a strong mitonuclear discordance, with mitochondrial markers showing a conspicuous geographic pattern, while diagnostic nuclear SNPs coexisted in both the same populations and individuals, providing evidence of hybridization in many possible combinations. The introgression is asymmetric, with S. italicus mitogenome usually associated with S. a. ambrosii and, to a lesser extent, to S. strinatii nuclear alleles. This finding confirms that Plethodontid are a group of choice to investigate hybridization mechanisms and suggests that behavioural, genetic and ecological components may concur in determining the direction and extent of introgression.

Selection on the Fly: Short-Term Adaptation to an Altered Sexual Selection Regime in Drosophila pseudoobscura

Experimental evolution studies are powerful approaches to examine the evolutionary history of lab populations. Such studies have shed light on how selection changes phenotypes and genotypes. Most of these studies have not examined the time course of adaptation under sexual selection manipulation, by resequencing the populations' genomes at multiple time points. Here, we analyze allele frequency trajectories in Drosophila pseudoobscura where we altered their sexual selection regime for 200 generations and sequenced pooled populations at 5 time points. The intensity of sexual selection was either relaxed in monogamous populations (M) or elevated in polyandrous lines (E). We present a comprehensive study of how selection alters population genetics parameters at the chromosome and gene level. We investigate differences in the effective population size-Ne-between the treatments, and perform a genome-wide scan to identify signatures of selection from the time-series data. We found genomic signatures of adaptation to both regimes in D. pseudoobscura. There are more significant variants in E lines as expected from stronger sexual selection. However, we found that the response on the X chromosome was substantial in both treatments, more pronounced in E and restricted to the more recently sex-linked chromosome arm XR in M. In the first generations of experimental evolution, we estimate Ne to be lower on the X in E lines, which might indicate a swift adaptive response at the onset of selection. Additionally, the third chromosome was affected by elevated polyandry whereby its distal end harbors a region showing a strong signal of adaptive evolution especially in E lines.

Speciation and development

Understanding general principles about the origin of species remains one of the foundational challenges in evolutionary biology. The genomic divergence between groups of individuals can spawn hybrid inviability and hybrid sterility, which presents a tantalizing developmental problem. Divergent developmental programs may yield either conserved or divergent phenotypes relative to ancestral traits, both of which can be responsible for reproductive isolation during the speciation process. The genetic mechanisms of developmental evolution involve cis- and trans-acting gene regulatory change, protein-protein interactions, genetic network structures, dosage, and epigenetic regulation, all of which also have roots in population genetic and molecular evolutionary processes. Toward the goal of demystifying Darwin's "mystery of mysteries," this review integrates microevolutionary concepts of genetic change with principles of organismal development, establishing explicit links between population genetic process and developmental mechanisms in the production of macroevolutionary pattern. This integration aims to establish a more unified view of speciation that binds process and mechanism.

Comparative genomics reveals the hybrid origin of a macaque group

Although species can arise through hybridization, compelling evidence for hybrid speciation has been reported only rarely in animals. Here, we present phylogenomic analyses on genomes from 12 macaque species and show that the fascicularis group originated from an ancient hybridization between the sinica and silenus groups ~3.45 to 3.56 million years ago. The X chromosomes and low-recombination regions exhibited equal contributions from each parental lineage, suggesting that they were less affected by subsequent backcrossing and hence could have played an important role in maintaining hybrid integrity. We identified many reproduction-associated genes that could have contributed to the development of the mixed sexual phenotypes characteristic of the fascicularis group. The phylogeny within the silenus group was also resolved, and functional experimentation confirmed that all extant Western silenus species are susceptible to HIV-1 infection. Our study provides novel insights into macaque evolution and reveals a hybrid speciation event that has occurred only very rarely in primates.

Target enrichment of long open reading frames and ultraconserved elements to link microevolution and macroevolution in non-model organisms

Despite the increasing accessibility of high-throughput sequencing, obtaining high-quality genomic data on non-model organisms without proximate well-assembled and annotated genomes remains challenging. Here, we describe a workflow that takes advantage of distant genomic resources and ingroup transcriptomes to select and jointly enrich long open reading frames (ORFs) and ultraconserved elements (UCEs) from genomic samples for integrative studies of microevolutionary and macroevolutionary dynamics. This workflow is applied to samples of the African unionid bivalve tribe Coelaturini (Parreysiinae) at basin and continent-wide scales. Our results indicate that ORFs are efficiently captured without prior identification of intron-exon boundaries. The enrichment of UCEs was less successful, but nevertheless produced substantial data sets. Exploratory continent-wide phylogenetic analyses with ORF supercontigs (>515,000 parsimony informative sites) resulted in a fully resolved phylogeny, the backbone of which was also retrieved with UCEs (>11,000 informative sites). Variant calling on ORFs and UCEs of Coelaturini from the Malawi Basin produced ~2000 SNPs per population pair. Estimates of nucleotide diversity and population differentiation were similar for ORFs and UCEs. They were low compared to previous estimates in molluscs, but comparable to those in recently diversifying Malawi cichlids and other taxa at an early stage of speciation. Skimming off-target sequence data from the same enriched libraries of Coelaturini from the Malawi Basin, we reconstructed the maternally-inherited mitogenome, which displays the gene order inferred for the most recent common ancestor of Unionidae. Overall, our workflow and results provide exciting perspectives for integrative genomic studies of microevolutionary and macroevolutionary dynamics in non-model organisms.