Evaluating the use of ABBA-BABA statistics to locate introgressed loci

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.

Geogenomic Predictors of Genetree Heterogeneity Explain Phylogeographic and Introgression History: A Case Study in an Amazonian Bird (Thamnophilus aethiops)

Can knowledge about genome architecture inform biogeographic and phylogenetic inference? Selection, drift, recombination, and gene flow interact to produce a genomic landscape of divergence wherein patterns of differentiation and genealogy vary nonrandomly across the genomes of diverging populations. For instance, genealogical patterns that arise due to gene flow should be more likely to occur on smaller chromosomes, which experience high recombination, whereas those tracking histories of geographic isolation (reduced gene flow caused by a barrier) and divergence should be more likely to occur on larger and sex chromosomes. In Amazonia, populations of many bird species diverge and introgress across rivers, resulting in reticulated genomic signals. Herein, we used reduced representation genomic data to disentangle the evolutionary history of 4 populations of an Amazonian antbird, Thamnophilus aethiops, whose biogeographic history was associated with the dynamic evolution of the Madeira River Basin. Specifically, we evaluate whether a large river capture event ca. 200 Ka, gave rise to reticulated genealogies in the genome by making spatially explicit predictions about isolation and gene flow based on knowledge about genomic processes. We first estimated chromosome-level phylogenies and recovered 2 primary topologies across the genome. The first topology (T1) was most consistent with predictions about population divergence and was recovered for the Z-chromosome. The second (T2), was consistent with predictions about gene flow upon secondary contact. To evaluate support for these topologies, we trained a convolutional neural network to classify our data into alternative diversification models and estimate demographic parameters. The best-fit model was concordant with T1 and included gene flow between non-sister taxa. Finally, we modeled levels of divergence and introgression as functions of chromosome length and found that smaller chromosomes experienced higher gene flow. Given that (1) genetrees supporting T2 were more likely to occur on smaller chromosomes and (2) we found lower levels of introgression on larger chromosomes (and especially the Z-chromosome), we argue that T1 represents the history of population divergence across rivers and T2 the history of secondary contact due to barrier loss. Our results suggest that a significant portion of genomic heterogeneity arises due to extrinsic biogeographic processes such as river capture interacting with intrinsic processes associated with genome architecture. Future phylogeographic studies would benefit from accounting for genomic processes, as different parts of the genome reveal contrasting, albeit complementary histories, all of which are relevant for disentangling the intricate geogenomic mechanisms of biotic diversification. [Amazonia; biogeography; demographic modeling; gene flow; gene tree; genome architecture; geogenomics; introgression; linked selection; neural network; phylogenomic; phylogeography; reproductive isolation; speciation; species tree.].

Conservation genomics provides insights into genetic resilience and adaptation of the endangered Chinese hazelnut, Corylus chinensis

Global climate change has increased concerns regarding biodiversity loss. However, many key conservation issues still required further research, including demographic history, deleterious mutation load, adaptive evolution, and putative introgression. Here we generated the first chromosome-level genome of the endangered Chinese hazelnut, Corylus chinensis, and compared the genomic signatures with its sympatric widespread C. kwechowensis-C. yunnanensis complex. We found large genome rearrangements across all Corylus species and identified species-specific expanded gene families that may be involved in adaptation. Population genomics revealed that both C. chinensis and the C. kwechowensis-C. yunnanensis complex had diverged into two genetic lineages, forming a consistent pattern of southwestern-northern differentiation. Population size of the narrow southwestern lineages of both species have decreased continuously since the late Miocene, whereas the widespread northern lineages have remained stable (C. chinensis) or have even recovered from population bottlenecks (C. kwechowensis-C. yunnanensis complex) during the Quaternary. Compared with C. kwechowensis-C. yunnanensis complex, C. chinensis showed significantly lower genomic diversity and higher inbreeding level. However, C. chinensis carried significantly fewer deleterious mutations than C. kwechowensis-C. yunnanensis complex, as more effective purging selection reduced the accumulation of homozygous variants. We also detected signals of positive selection and adaptive introgression in different lineages, which facilitated the accumulation of favorable variants and formation of local adaptation. Hence, both types of selection and exogenous introgression could have mitigated inbreeding and facilitated survival and persistence of C. chinensis. Overall, our study provides critical insights into lineage differentiation, local adaptation, and the potential for future recovery of endangered trees.

Genomic evidence for evolutionary history and local adaptation of two endemic apricots: Prunus hongpingensis and P. zhengheensis

Apricot, belonging to the Armeniaca section of Rosaceae, is one of the economically important crop fruits that has been extensively cultivated. The natural wild apricots offer valuable genetic resources for crop improvement. However, some of them are endemic, with small populations, and are even at risk of extinction. In this study we unveil chromosome-level genome assemblies for two southern China endemic apricots, Prunus hongpingensis (PHP) and P. zhengheensis (PZH). We also characterize their evolutionary history and the genomic basis of their local adaptation using whole-genome resequencing data. Our findings reveal that PHP and PZH are closely related to Prunus armeniaca and form a distinct lineage. Both species experienced a decline in effective population size following the Last Glacial Maximum (LGM), which likely contributed to their current small population sizes. Despite the observed decrease in genetic diversity and heterozygosity, we do not observe an increased accumulation of deleterious mutations in these two endemic apricots. This is likely due to the combined effects of a low inbreeding coefficient and strong purifying selection. Furthermore, we identify a set of genes that have undergone positive selection and are associated with local environmental adaptation in PHP and PZH, respectively. These candidate genes can serve as valuable genetic resources for targeted breeding and improvement of cultivated apricots. Overall, our study not only enriches our comprehension of the evolutionary history of apricot species but also offers crucial insights for the conservation and future breeding of other endemic species amidst rapid climate changes.

Into the Wild: A novel wild-derived inbred strain resource expands the genomic and phenotypic diversity of laboratory mouse models

The laboratory mouse has served as the premier animal model system for both basic and preclinical investigations for over a century. However, laboratory mice capture only a subset of the genetic variation found in wild mouse populations, ultimately limiting the potential of classical inbred strains to uncover phenotype-associated variants and pathways. Wild mouse populations are reservoirs of genetic diversity that could facilitate the discovery of new functional and disease-associated alleles, but the scarcity of commercially available, well-characterized wild mouse strains limits their broader adoption in biomedical research. To overcome this barrier, we have recently developed, sequenced, and phenotyped a set of 11 inbred strains derived from wild-caught Mus musculus domesticus. Each of these "Nachman strains" immortalizes a unique wild haplotype sampled from one of five environmentally distinct locations across North and South America. Whole genome sequence analysis reveals that each strain carries between 4.73-6.54 million single nucleotide differences relative to the GRCm39 mouse reference, with 42.5% of variants in the Nachman strain genomes absent from current classical inbred mouse strain panels. We phenotyped the Nachman strains on a customized pipeline to assess the scope of disease-relevant neurobehavioral, biochemical, physiological, metabolic, and morphological trait variation. The Nachman strains exhibit significant inter-strain variation in >90% of 1119 surveyed traits and expand the range of phenotypic diversity captured in classical inbred strain panels. These novel wild-derived inbred mouse strain resources are set to empower new discoveries in both basic and preclinical research.

Significant genomic introgression from grey junglefowl (Gallus sonneratii) to domestic chickens (Gallus gallus domesticus)

BACKGROUND

Chicken is one of the most numerous and widely distributed species around the world, and many studies support the multiple ancestral origins of domestic chickens. The research regarding the yellow skin phenotype in domestic chickens (regulated by BCO2) likely originating from the grey junglefowl serves as crucial evidence for demonstrating the multiple origins of chickens. However, beyond the BCO2 gene region, much remains unknown about the introgression from the grey junglefowl into domestic chickens. Therefore, in this study, based on whole-genome data of 149 samples including 4 species of wild junglefowls and 13 local domestic chicken breeds, we explored the introgression events from the grey junglefowl to domestic chickens.

RESULTS

We successfully detected introgression regions besides BCO2, including two associated with growth trait (IGFBP2 and TKT), one associated with angiogenesis (TIMP3) and two members of the heat shock protein family (HSPB2 and CRYAB). Our findings suggest that the introgression from the grey junglefowl may impact the growth performance of chickens. Furthermore, we revealed introgression events from grey junglefowl at the BCO2 region in multiple domestic chicken breeds, indicating a phenomenon where the yellow skin phenotype likely underwent strong selection and was retained. Additionally, our haplotype analysis shed light on BCO2 introgression event from different sources of grey junglefowl into domestic chickens, possibly suggesting multiple genetic flows between the grey junglefowl and domestic chickens.

CONCLUSIONS

In summary, our findings provide evidences of the grey junglefowl contributing to the genetic diversity of domestic chickens, laying the foundation for a deeper understanding of the genetic composition within domestic chickens, and offering new perspectives on the impact of introgression on domestic chickens.

Adaptive introgression of a visual preference gene

Visual preferences are important drivers of mate choice and sexual selection, but little is known of how they evolve at the genetic level. In this study, we took advantage of the diversity of bright warning patterns displayed by Heliconius butterflies, which are also used during mate choice. Combining behavioral, population genomic, and expression analyses, we show that two Heliconius species have evolved the same preferences for red patterns by exchanging genetic material through hybridization. Neural expression of regucalcin1 correlates with visual preference across populations, and disruption of regucalcin1 with CRISPR-Cas9 impairs courtship toward conspecific females, providing a direct link between gene and behavior. Our results support a role for hybridization during behavioral evolution and show how visually guided behaviors contributing to adaptation and speciation are encoded within the genome.

The extent of introgression between incipient Clarkia species is determined by temporal environmental variation and mating system

Introgression is pervasive across the tree of life but varies across taxa, geography, and genomic regions. However, the factors modulating this variation and how they may be affected by global change are not well understood. Here, we used 200 genomes and a 15-y site-specific environmental dataset to investigate the effects of environmental variation and mating system divergence on the magnitude of introgression between a recently diverged outcrosser-selfer pair of annual plants in the genus Clarkia. These sister taxa diverged very recently and subsequently came into secondary sympatry where they form replicated contact zones. Consistent with observations of other outcrosser-selfer pairs, we found that introgression was asymmetric between taxa, with substantially more introgression from the selfer to the outcrosser. This asymmetry was caused by a bias in the direction of initial F1 hybrid formation and subsequent backcrossing. We also found extensive variation in the outcrosser's admixture proportion among contact zones, which was predicted nearly entirely by interannual variance in spring precipitation. Greater fluctuations in spring precipitation resulted in higher admixture proportions, likely mediated by the effects of spring precipitation on the expression of traits that determine premating reproductive isolation. Climate-driven hybridization dynamics may be particularly affected by global change, potentially reshaping species boundaries and adaptation to novel environments.

Recombination, admixture and genome instability shape the genomic landscape of Saccharomyces cerevisiae derived from spontaneous grape ferments

Cultural exchange of fermentation techniques has driven the spread of Saccharomyces cerevisiae across the globe, establishing natural populations in many countries. Despite this, Oceania is thought to lack native populations of S. cerevisiae, only being introduced after colonisation. Here we investigate the genomic landscape of 411 S. cerevisiae isolated from spontaneous grape fermentations in Australia across multiple locations, years, and grape cultivars. Spontaneous fermentations contained highly recombined mosaic strains that exhibited high levels of genome instability. Assigning genomic windows to putative ancestral origin revealed that few closely related starter lineages have come to dominate the genetic landscape, contributing most of the genetic variation. Fine-scale phylogenetic analysis of loci not observed in strains of commercial wine origin identified widespread admixture with European derived beer yeast along with three independent admixture events from potentially endemic Oceanic lineages that was associated with genome instability. Finally, we investigated Australian ecological niches for basal isolates, identifying phylogenetically distinct S. cerevisiae of non-European, non-domesticated origin associated with admixture loci. Our results illustrate the effect commercial use of microbes may have on local microorganism genetic diversity and demonstrates the presence of non-domesticated, potentially endemic lineages of S. cerevisiae in Australian niches that are actively admixing.

Gene flow and an anomaly zone complicate phylogenomic inference in a rapidly radiated avian family (Prunellidae)

BACKGROUND

Resolving the phylogeny of rapidly radiating lineages presents a challenge when building the Tree of Life. An Old World avian family Prunellidae (Accentors) comprises twelve species that rapidly diversified at the Pliocene-Pleistocene boundary.

RESULTS

Here we investigate the phylogenetic relationships of all species of Prunellidae using a chromosome-level de novo assembly of Prunella strophiata and 36 high-coverage resequenced genomes. We use homologous alignments of thousands of exonic and intronic loci to build the coalescent and concatenated phylogenies and recover four different species trees. Topology tests show a large degree of gene tree-species tree discordance but only 40-54% of intronic gene trees and 36-75% of exonic genic trees can be explained by incomplete lineage sorting and gene tree estimation errors. Estimated branch lengths for three successive internal branches in the inferred species trees suggest the existence of an empirical anomaly zone. The most common topology recovered for species in this anomaly zone was not similar to any coalescent or concatenated inference phylogenies, suggesting presence of anomalous gene trees. However, this interpretation is complicated by the presence of gene flow because extensive introgression was detected among these species. When exploring tree topology distributions, introgression, and regional variation in recombination rate, we find that many autosomal regions contain signatures of introgression and thus may mislead phylogenetic inference. Conversely, the phylogenetic signal is concentrated to regions with low-recombination rate, such as the Z chromosome, which are also more resistant to interspecific introgression.

CONCLUSIONS

Collectively, our results suggest that phylogenomic inference should consider the underlying genomic architecture to maximize the consistency of phylogenomic signal.

Karyotypic stasis and swarming influenced the evolution of viral tolerance in a species-rich bat radiation

The emergence of COVID-19 and severe acute respiratory syndrome (SARS) has prioritized understanding bats' viral tolerance. Myotis bats are exceptionally species rich and have evolved viral tolerance. They also exhibit swarming, a cryptic behavior where large, multi-species assemblages gather for mating, which has been hypothesized to promote interspecific hybridization. To resolve the coevolution of genome architecture and their unusual antiviral tolerance, we undertook a phylogenomic analysis of 60 Old World Myotis genomes. We demonstrate an extensive history of introgressive hybridization that has replaced the species phylogeny across 17%-93% of the genome except for pericentromeric regions of macrochromosomes. Introgression tracts were enriched on microchromosome regions containing key antiviral pathway genes overexpressed during viral challenge experiments. Together, these results suggest that the unusual Myotis karyotype may have evolved to selectively position immune-related genes in high recombining genomic regions prone to introgression of divergent alleles, including a diversity of interleukin loci responsible for the release of pro-inflammatory cytokines.

A graph-based pan-genome of Brassica oleracea provides new insights into its domestication and morphotype diversification

The domestication of Brassica oleracea has resulted in diverse morphological types with distinct patterns of organ development. Here we report a graph-based pan-genome of B. oleracea constructed from high-quality genome assemblies of different morphotypes. The pan-genome harbors over 200 structural variant hotspot regions enriched in auxin- and flowering-related genes. Population genomic analyses revealed that early domestication of B. oleracea focused on leaf or stem development. Gene flows resulting from agricultural practices and variety improvement were detected among different morphotypes. Selective-sweep and pan-genome analyses identified an auxin-responsive small auxin up-regulated RNA gene and a CLAVATA3/ESR-RELATED family gene as crucial players in leaf-stem differentiation during the early stage of B. oleracea domestication and the BoKAN1 gene as instrumental in shaping the leafy heads of cabbage and Brussels sprouts. Our pan-genome and functional analyses further revealed that variations in the BoFLC2 gene play key roles in the divergence of vernalization and flowering characteristics among different morphotypes, and variations in the first intron of BoFLC3 are involved in fine-tuning the flowering process in cauliflower. This study provides a comprehensive understanding of the pan-genome of B. oleracea and sheds light on the domestication and differential organ development of this globally important crop species.

Secondary Contact, Introgressive Hybridization, and Genome Stabilization in Sticklebacks

Advances in genomic studies have revealed that hybridization in nature is pervasive and raised questions about the dynamics of different genetic and evolutionary factors following the initial hybridization event. While recent research has proposed that the genomic outcomes of hybridization might be predictable to some extent, many uncertainties remain. With comprehensive whole-genome sequence data, we investigated the genetic introgression between 2 divergent lineages of 9-spined sticklebacks (Pungitius pungitius) in the Baltic Sea. We found that the intensity and direction of selection on the introgressed variation has varied across different genomic elements: while functionally important regions displayed reduced rates of introgression, promoter regions showed enrichment. Despite the general trend of negative selection, we identified specific genomic regions that were enriched for introgressed variants, and within these regions, we detected footprints of selection, indicating adaptive introgression. Geographically, we found the selection against the functional changes to be strongest in the vicinity of the secondary contact zone and weaken as a function of distance from the initial contact. Altogether, the results suggest that the stabilization of introgressed variation in the genomes is a complex, multistage process involving both negative and positive selection. In spite of the predominance of negative selection against introgressed variants, we also found evidence for adaptive introgression variants likely associated with adaptation to Baltic Sea environmental conditions.

Phylogenomics and species delimitation of an abundant and little-studied Amazonian forest spiny rat

Species delimitation studies based on integrating different datasets such as genomic, morphometric, and cytogenetics data are rare in studies focused on Neotropical rodents. As a consequence, the evolutionary history of most of these genera remains poorly understood. Proechimys is a highly diverse and widely distributed genus of Neotropical spiny rats with unique traits like multiple sympatry, micro-habitat segregation, and fuzzy species limits. Here, we applied RAD-Seq to infer the phylogenetic relationships, estimate the species boundaries, and estimate the divergence times for Proechimys, one of the most common and least studied small mammals in the Amazon. We tested whether inferred lineages in the phylogenetic trees could be considered distinct species based on the genomic dataset and morphometric data. Analyses revealed the genus is not monophyletic, with Proechimys hoplomyoides sister to a group of Hoplomys gymnurus + all other Proechimys species, contesting the generic status of Hoplomys. There are five main clades in Proechimys stricto sensu (excluding H. gymnurus and P. hoplomyoides). Species delimitation analyses supported 25 species within the genus Proechimys. The five main clades in Proechimys stricto sensu also showed similar ages for their origins, and two rapid diversification events were identified in the Early Pliocene and in the Early Pleistocene. Most cases of sympatry in Proechimys occur among species from the different main clades, and although Proechimys is an inhabitant of the Amazon, three species occupied the Cerrado biome during the Pleistocene. We could associate available nominal taxon, cytogenetics information, and DNA sequences in Genbank to most of the 25 species we hypothesized from our delimitation analyses. Based on our analyses, we estimate that eight forms represent putative new species that need a taxonomic revision.

Comparative genomic analyses provide new insights into evolutionary history and conservation genomics of gorillas

Genome sequencing is a powerful tool to understand species evolutionary history, uncover genes under selection, which could be informative of local adaptation, and infer measures of genetic diversity, inbreeding and mutational load that could be used to inform conservation efforts. Gorillas, critically endangered primates, have received considerable attention and with the recently sequenced Bwindi mountain gorilla population, genomic data is now available from all gorilla subspecies and both mountain gorilla populations. Here, we reanalysed this rich dataset with a focus on evolutionary history, local adaptation and genomic parameters relevant for conservation. We estimate a recent split between western and eastern gorillas of 150,000-180,000 years ago, with gene flow around 20,000 years ago, primarily between the Cross River and Grauer's gorilla subspecies. This gene flow event likely obscures evolutionary relationships within eastern gorillas: after excluding putatively introgressed genomic regions, we uncover a sister relationship between Virunga mountain gorillas and Grauer's gorillas to the exclusion of Bwindi mountain gorillas. This makes mountain gorillas paraphyletic. Eastern gorillas are less genetically diverse and more inbred than western gorillas, yet we detected lower genetic load in the eastern species. Analyses of indels fit remarkably well with differences in genetic diversity across gorilla taxa as recovered with nucleotide diversity measures. We also identified genes under selection and unique gene variants specific for each gorilla subspecies, encoding, among others, traits involved in immunity, diet, muscular development, hair morphology and behavior. The presence of this functional variation suggests that the subspecies may be locally adapted. In conclusion, using extensive genomic resources we provide a comprehensive overview of gorilla genomic diversity, including a so-far understudied Bwindi mountain gorilla population, identify putative genes involved in local adaptation, and detect population-specific gene flow across gorilla species.

Phylogenomics reveal Populusgonggaensis as a hybrid between P.lasiocarpa and P.cathayana (Salicaceae)

High levels of intra-specific polymorphism and frequent hybridisation make it difficult to define species and correctly apply their scientific names. Populus L. is a challenging genus with plentiful natural and artificial hybrids. This study is a part of the project 'Flora of Pan-Himalaya' and aims to determine the taxonomic identity of P.gonggaensis N. Chao & J.R. He and to find out whether it is of hybrid origin. Whole-genome sequencing data were obtained from 57 samples. The SNP matrix was developed for phylogenetic reconstruction, ABBA-BABA statistics, PCA and ADMIXTURE analysis. The results indicate that P.gonggaensis is a spontaneous hybrid between P.lasiocarpa and P.cathayana. This study points out the importance of SNP data and comprehensive analyses for discovering the potential interspecific hybridisation and clarifies the usage of the name. In addition, the lectotype of P.gonggaensis was designated.

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.

The genetic basis of a recent transition to live-bearing in marine snails

Key innovations are fundamental to biological diversification, but their genetic basis is poorly understood. A recent transition from egg-laying to live-bearing in marine snails (Littorina spp.) provides the opportunity to study the genetic architecture of an innovation that has evolved repeatedly across animals. Individuals do not cluster by reproductive mode in a genome-wide phylogeny, but local genealogical analysis revealed numerous small genomic regions where all live-bearers carry the same core haplotype. Candidate regions show evidence for live-bearer-specific positive selection and are enriched for genes that are differentially expressed between egg-laying and live-bearing reproductive systems. Ages of selective sweeps suggest that live-bearer-specific alleles accumulated over more than 200,000 generations. Our results suggest that new functions evolve through the recruitment of many alleles rather than in a single evolutionary step.

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.

Ecological Diversification in an Adaptive Radiation of Plants: The Role of De Novo Mutation and Introgression

Adaptive radiations are characterized by rapid ecological diversification and speciation events, leading to fuzzy species boundaries between ecologically differentiated species. Adaptive radiations are therefore key systems for understanding how species are formed and maintained, including the role of de novo mutations versus preexisting variation in ecological adaptation and the genome-wide consequences of hybridization events. For example, adaptive introgression, where beneficial alleles are transferred between lineages through hybridization, may fuel diversification in adaptive radiations and facilitate adaptation to new environments. In this study, we employed whole-genome resequencing data to investigate the evolutionary origin of hummingbird-pollinated flowers and to characterize genome-wide patterns of phylogenetic discordance and introgression in Penstemon subgenus Dasanthera, a small and diverse adaptive radiation of plants. We found that magenta hummingbird-adapted flowers have apparently evolved twice from ancestral blue-violet bee-pollinated flowers within this radiation. These shifts in flower color are accompanied by a variety of inactivating mutations to a key anthocyanin pathway enzyme, suggesting that independent de novo loss-of-function mutations underlie the parallel evolution of this trait. Although patterns of introgression and phylogenetic discordance were heterogenous across the genome, a strong effect of gene density suggests that, in general, natural selection opposes introgression and maintains genetic differentiation in gene-rich genomic regions. Our results highlight the importance of both de novo mutation and introgression as sources of evolutionary change and indicate a role for de novo mutation in driving parallel evolution in adaptive radiations.

Ecological diversification in an adaptive radiation of plants: the role of de novo mutation and introgression

Adaptive radiations are characterized by rapid ecological diversification and speciation events, leading to fuzzy species boundaries between ecologically differentiated species. Adaptive radiations are therefore key systems for understanding how species are formed and maintained, including the role of de novo mutations vs. pre-existing variation in ecological adaptation and the genome-wide consequences of hybridization events. For example, adaptive introgression, where beneficial alleles are transferred between lineages through hybridization, may fuel diversification in adaptive radiations and facilitate adaptation to new environments. In this study, we employed whole-genome resequencing data to investigate the evolutionary origin of hummingbird-pollinated flowers and to characterize genome-wide patterns of phylogenetic discordance and introgression in Penstemon subgenus Dasanthera, a small and diverse adaptive radiation of plants. We found that magenta hummingbird-adapted flowers have apparently evolved twice from ancestral blue-violet bee-pollinated flowers within this radiation. These shifts in flower color are accompanied by a variety of inactivating mutations to a key anthocyanin pathway enzyme, suggesting that independent de novo loss-of-function mutations underlie parallel evolution of this trait. Although patterns of introgression and phylogenetic discordance were heterogenous across the genome, a strong effect of gene density suggests that, in general, natural selection opposes introgression and maintains genetic differentiation in gene-rich genomic regions. Our results highlight the importance of both de novo mutation and introgression as sources of evolutionary change and indicate a role for de novo mutation in driving parallel evolution in adaptive radiations.

Whole genome sequencing of Crassostrea ariakensis (Mollusca: Ostreidae) and C. hongkongensis expands understandings of stress resistance in sessile oysters

To understand the environmental adaptations among sessile bivalves lacking adaptive immunity, a series of analyses were conducted, with special emphasis on the widely distributed C. ariakensis. Employing Pacbio sequencing and Hi-C technologies, whole genome for each of a C. ariakensis (southern China) and C. hongkongensis individual was generated, with the contig N50 reaching 6.2 and 13.0 Mb, respectively. Each genome harbored over 30,000 protein-coding genes, with approximately half of each genome consisting of repeats. Genome alignment suggested possible introgression between C. gigas and C. ariakensis (northern China), and re-sequencing data corroborated this result and indicated significant gene flow between C. gigas and C. ariakensis. These introgressed candidates, well-represented by genes related to immunity and osmotic pressure, may be associated with environmental stresses. Gene family dynamics modeling suggested immune-related genes were well represented among the expanded genes in C. ariakensis. These outcomes could be attributed to the spread of C. ariakensis.

The genomes of Darwin's primroses reveal chromosome-scale adaptive introgression and differential permeability of species boundaries

Introgression is an important source of genetic variation that can determine species adaptation to environmental conditions. Yet, definitive evidence of the genomic and adaptive implications of introgression in nature remains scarce. The widespread hybrid zones of Darwin's primroses (Primula elatior, Primula veris, and Primula vulgaris) provide a unique natural laboratory for studying introgression in flowering plants and the varying permeability of species boundaries. Through analysis of 650 genomes, we provide evidence of an introgressed genomic region likely to confer adaptive advantage in conditions of soil toxicity. We also document unequivocal evidence of chloroplast introgression, an important precursor to species-wide chloroplast capture. Finally, we provide the first evidence that the S-locus supergene, which controls heterostyly in primroses, does not introgress in this clade. Our results contribute novel insights into the adaptive role of introgression and demonstrate the importance of extensive genomic and geographical sampling for illuminating the complex nature of species boundaries.

Leveraging shared ancestral variation to detect local introgression

Introgression is a common evolutionary phenomenon that results in shared genetic material across non-sister taxa. Existing statistical methods such as Patterson's D statistic can detect introgression by measuring an excess of shared derived alleles between populations. The D statistic is effective to detect genome-wide patterns of introgression but can give spurious inferences of introgression when applied to local regions. We propose a new statistic, D+, that leverages both shared ancestral and derived alleles to infer local introgressed regions. Incorporating both shared derived and ancestral alleles increases the number of informative sites per region, improving our ability to identify local introgression. We use a coalescent framework to derive the expected value of this statistic as a function of different demographic parameters under an instantaneous admixture model and use coalescent simulations to compute the power and precision of D+. While the power of D and D+ is comparable, D+ has better precision than D. We apply D+ to empirical data from the 1000 Genome Project and Heliconius butterflies to infer local targets of introgression in humans and in butterflies.

Genome evolution is surprisingly predictable after initial hybridization

Over the past two decades, evolutionary biologists have come to appreciate that hybridization, or genetic exchange between distinct lineages, is remarkably common - not just in particular lineages but in taxonomic groups across the tree of life. As a result, the genomes of many modern species harbor regions inherited from related species. This observation has raised fundamental questions about the degree to which the genomic outcomes of hybridization are repeatable and the degree to which natural selection drives such repeatability. However, a lack of appropriate systems to answer these questions has limited empirical progress in this area. Here, we leverage independently formed hybrid populations between the swordtail fish Xiphophorus birchmanni and X. cortezi to address this fundamental question. We find that local ancestry in one hybrid population is remarkably predictive of local ancestry in another, demographically independent hybrid population. Applying newly developed methods, we can attribute much of this repeatability to strong selection in the earliest generations after initial hybridization. We complement these analyses with time-series data that demonstrates that ancestry at regions under selection has remained stable over the past ~40 generations of evolution. Finally, we compare our results to the well-studied X. birchmanni×X. malinche hybrid populations and conclude that deeper evolutionary divergence has resulted in stronger selection and higher repeatability in patterns of local ancestry in hybrids between X. birchmanni and X. cortezi.

Genomic landscape of introgression from the ghost lineage in a gobiid fish uncovers the generality of forces shaping hybrid genomes

Extinct lineages can leave legacies in the genomes of extant lineages through ancient introgressive hybridization. The patterns of genomic survival of these extinct lineages provide insight into the role of extinct lineages in current biodiversity. However, our understanding on the genomic landscape of introgression from extinct lineages remains limited due to challenges associated with locating the traces of unsampled 'ghost' extinct lineages without ancient genomes. Herein, we conducted population genomic analyses on the East China Sea (ECS) lineage of Chaenogobius annularis, which was suspected to have originated from ghost introgression, with the aim of elucidating its genomic origins and characterizing its landscape of introgression. By combining phylogeographic analysis and demographic modelling, we demonstrated that the ECS lineage originated from ancient hybridization with an extinct ghost lineage. Forward simulations based on the estimated demography indicated that the statistic γ of the HyDe analysis can be used to distinguish the differences in local introgression rates in our data. Consistent with introgression between extant organisms, we found reduced introgression from extinct lineage in regions with low recombination rates and with functional importance, thereby suggesting a role of linked selection that has eliminated the extinct lineage in shaping the hybrid genome. Moreover, we identified enrichment of repetitive elements in regions associated with ghost introgression, which was hitherto little known but was also observed in the re-analysis of published data on introgression between extant organisms. Overall, our findings underscore the unexpected similarities in the characteristics of introgression landscapes across different taxa, even in cases of ghost introgression.

Complex Evolutionary History With Extensive Ancestral Gene Flow in an African Primate Radiation

Understanding the drivers of speciation is fundamental in evolutionary biology, and recent studies highlight hybridization as an important evolutionary force. Using whole-genome sequencing data from 22 species of guenons (tribe Cercopithecini), one of the world's largest primate radiations, we show that rampant gene flow characterizes their evolutionary history and identify ancient hybridization across deeply divergent lineages that differ in ecology, morphology, and karyotypes. Some hybridization events resulted in mitochondrial introgression between distant lineages, likely facilitated by cointrogression of coadapted nuclear variants. Although the genomic landscapes of introgression were largely lineage specific, we found that genes with immune functions were overrepresented in introgressing regions, in line with adaptive introgression, whereas genes involved in pigmentation and morphology may contribute to reproductive isolation. In line with reports from other systems that hybridization might facilitate diversification, we find that some of the most species-rich guenon clades are of admixed origin. This study provides important insights into the prevalence, role, and outcomes of ancestral hybridization in a large mammalian radiation.

Genomic divergence and introgression between cryptic species of a widespread North American songbird

Analysis of genomic variation among related populations can sometimes reveal distinct species that were previously undescribed due to similar morphological appearances, and close examination of such cases can provide much insight regarding speciation. Genomic data can also reveal the role of reticulate evolution in differentiation and speciation. White-breasted nuthatches (Sitta carolinensis) are widely distributed North American songbirds that are currently classified as a single species but have been suspected to represent a case of cryptic speciation. Previous genetic analyses suggested four divergent groups, but it was unclear whether these represented multiple reproductively isolated species. Using extensive genomic sampling of over 350 white-breasted nuthatches from across North America and a new chromosome-level reference genome, we asked if white-breasted nuthatches are comprised of multiple species and whether introgression has occurred between divergent populations. Genomic variation of over 300,000 loci revealed four highly differentiated populations (Pacific, n = 45; Eastern, n = 23; Rocky Mountains North, n = 138; and Rocky Mountains South, n = 150) with geographic ranges that are adjacent. We observed a moderate degree of admixture between Rocky Mountain populations but only a small number of hybrids between the Rockies and the Eastern population. The rarity of hybrids together with high levels of differentiation between populations is supportive of populations having some level of reproductive isolation. Between populations, we show evidence for introgression from a divergent ghost lineage of white-breasted nuthatches into the Rocky Mountains South population, which is otherwise closely related to Rocky Mountains North. We conclude that white-breasted nuthatches are best considered at least three species and that ghost lineage introgression has contributed to differentiation between the two Rocky Mountain populations. White-breasted nuthatches provide a dramatic case of morphological similarity despite high genomic differentiation, and the varying levels of reproductive isolation among the four groups provide an example of the speciation continuum.

Testing Times: Challenges in Disentangling Admixture Histories in Recent and Complex Demographies

Paleogenomics has expanded our knowledge of human evolutionary history. Since the 2020s, the study of ancient DNA has increased its focus on reconstructing the recent past. However, the accuracy of paleogenomic methods in answering questions of historical and archaeological importance amidst the increased demographic complexity and decreased genetic differentiation within the historical period remains an open question. We used two simulation approaches to evaluate the limitations and behavior of commonly used methods, qpAdm and the f 3 -statistic, on admixture inference. The first is based on branch-length data simulated from four simple demographic models of varying complexities and configurations. The second, an analysis of Eurasian history composed of 59 populations using whole-genome data modified with ancient DNA conditions such as SNP ascertainment, data missingness, and pseudo-haploidization. We show that under conditions resembling historical populations, qpAdm can identify a small candidate set of true sources and populations closely related to them. However, in typical ancient DNA conditions, qpAdm is unable to further distinguish between them, limiting its utility for resolving fine-scaled hypotheses. Notably, we find that complex gene-flow histories generally lead to improvements in the performance of qpAdm and observe no bias in the estimation of admixture weights. We offer a heuristic for admixture inference that incorporates admixture weight estimate and P -values of qpAdm models, and f 3 -statistics to enhance the power to distinguish between multiple plausible candidates. Finally, we highlight the future potential of qpAdm through whole-genome branch-length f 2 -statistics, demonstrating the improved demographic inference that could be achieved with advancements in f -statistic estimations.

Phylogenomics Reveals High Levels of Incomplete Lineage Sorting at the Ancestral Nodes of the Macaque Radiation

The genus Macaca includes 23 species assigned into 4 to 7 groups. It exhibits the largest geographic range and represents the most successful example of adaptive radiation of nonhuman primates. However, intrageneric phylogenetic relationships among species remain controversial and have not been resolved so far. In this study, we conducted a phylogenomic analysis on 16 newly generated and 8 published macaque genomes. We found strong evidence supporting the division of this genus into 7 species groups. Incomplete lineage sorting (ILS) was the primary factor contributing to the discordance observed among gene trees; however, we also found evidence of hybridization events, specifically between the ancestral arctoides/sinica and silenus/nigra lineages that resulted in the hybrid formation of the fascicularis/mulatta group. Combined with fossil data, our phylogenomic data were used to establish a scenario for macaque radiation. These findings provide insights into ILS and potential ancient introgression events that were involved in the radiation of macaques, which will lead to a better understanding of the rapid speciation occurring in nonhuman primates.

Inferring Historical Introgression with Deep Learning

Resolving phylogenetic relationships among taxa remains a challenge in the era of big data due to the presence of genetic admixture in a wide range of organisms. Rapidly developing sequencing technologies and statistical tests enable evolutionary relationships to be disentangled at a genome-wide level, yet many of these tests are computationally intensive and rely on phased genotypes, large sample sizes, restricted phylogenetic topologies, or hypothesis testing. To overcome these difficulties, we developed a deep learning-based approach, named ERICA, for inferring genome-wide evolutionary relationships and local introgressed regions from sequence data. ERICA accepts sequence alignments of both population genomic data and multiple genome assemblies, and efficiently identifies discordant genealogy patterns and exchanged regions across genomes when compared with other methods. We further tested ERICA using real population genomic data from Heliconius butterflies that have undergone adaptive radiation and frequent hybridization. Finally, we applied ERICA to characterize hybridization and introgression in wild and cultivated rice, revealing the important role of introgression in rice domestication and adaptation. Taken together, our findings demonstrate that ERICA provides an effective method for teasing apart evolutionary relationships using whole genome data, which can ultimately facilitate evolutionary studies on hybridization and introgression.

Indirect genetic effects are shaped by demographic history and ecology in Arabidopsis thaliana

The phenotype of an individual can be affected by the genes of its conspecifics through indirect genetic effects (IGEs). IGEs have been studied across different organisms including wild and domesticated animals and plants, but little is known about their genetic architecture. Here, in a large-scale intraspecific interaction experiment, we show that the contribution of IGEs to the biomass variation of Arabidopsis thaliana is comparable to values classically reported in animals. Moreover, we identify 11 loci explaining 85.1% of the variability in IGEs. We find that positive IGE alleles (that is, those with positive effects on neighbour biomass) occur both in relict accessions from southern Eurasia and in post-glacial colonizers from northern Scandinavia, and that they are likely to have two divergent origins: for nine loci, they evolved in the post-glacial colonizers independently from the relicts, while the two others were introgressed in the post-glacial colonizer from the relicts. Finally, we find that variation in IGEs probably reflects divergent adaptations to the contrasting environments of the edges and the centre of the native range of the species. These findings reveal a surprisingly tractable genetic basis of IGEs in A. thaliana that is shaped by the ecology and the demographic history of the species.

False discovery rates of qpAdm-based screens for genetic admixture

Although a broad range of methods exists for reconstructing population history from genome-wide single nucleotide polymorphism data, just a few methods gained popularity in archaeogenetics: principal component analysis (PCA); ADMIXTURE, an algorithm that models individuals as mixtures of multiple ancestral sources represented by actual or inferred populations; formal tests for admixture such as f3-statistics and D/f4-statistics; and qpAdm, a tool for fitting two-component and more complex admixture models to groups or individuals. Despite their popularity in archaeogenetics, which is explained by modest computational requirements and ability to analyze data of various types and qualities, protocols relying on qpAdm that screen numerous alternative models of varying complexity and find "fitting" models (often considering both estimated admixture proportions and p-values as a composite criterion of model fit) remain untested on complex simulated population histories in the form of admixture graphs of random topology. We analyzed genotype data extracted from such simulations and tested various types of high-throughput qpAdm protocols ("rotating" and "non-rotating", with or without temporal stratification of target groups and proxy ancestry sources, and with or without a "model competition" step). We caution that high-throughput qpAdm protocols may be inappropriate for exploratory analyses in poorly studied regions/periods since their false discovery rates varied between 12% and 68% depending on the details of the protocol and on the amount and quality of simulated data (i.e., >12% of fitting two-way admixture models imply gene flows that were not simulated). We demonstrate that for reducing false discovery rates of qpAdm protocols to nearly 0% it is advisable to use large SNP sets with low missing data rates, the rotating qpAdm protocol with a strictly enforced rule that target groups do not pre-date their proxy sources, and an unsupervised ADMIXTURE analysis as a way to verify feasible qpAdm models. Our study has a number of limitations: for instance, these recommendations depend on the assumption that the underlying genetic history is a complex admixture graph and not a stepping-stone model.

Genetic and cultural adaptations underlie the establishment of dairy pastoralism in the Tibetan Plateau

BACKGROUND

Domestication and introduction of dairy animals facilitated the permanent human occupation of the Tibetan Plateau. Yet the history of dairy pastoralism in the Tibetan Plateau remains poorly understood. Little is known how Tibetans adapted to milk and dairy products.

RESULTS

We integrated archeological evidence and genetic analysis to show the picture that the dairy ruminants, together with dogs, were introduced from West Eurasia into the Tibetan Plateau since ~ 3600 years ago. The genetic admixture between the exotic and indigenous dogs enriched the candidate lactase persistence (LP) allele 10974A > G of West Eurasian origin in Tibetan dogs. In vitro experiments demonstrate that - 13838G > A functions as a LP allele in Tibetans. Unlike multiple LP alleles presenting selective signatures in West Eurasians and South Asians, the de novo origin of Tibetan-specific LP allele - 13838G > A with low frequency (~ 6-7%) and absence of selection corresponds - 13910C > T in pastoralists across eastern Eurasia steppe.

CONCLUSIONS

Results depict a novel scenario of genetic and cultural adaptations to diet and expand current understanding of the establishment of dairy pastoralism in the Tibetan Plateau.

Understanding natural selection and similarity: Convergent, parallel and repeated evolution

Parallel and convergent evolution offer some of the most compelling evidence for the significance of natural selection in evolution, as the emergence of similar adaptive solutions is unlikely to occur by random chance alone. However, these terms are often employed inconsistently, leading to misinterpretation and confusion, and recently proposed definitions have unintentionally diminished the emphasis on the evolution of similar adaptive solutions. Here, I examine various conceptual frameworks and definitions related to parallel and convergent evolution and propose a consolidated framework that enhances our comprehension of these evolutionary patterns. The primary aim of this framework is to harmonize the concepts of parallel and convergent evolution together with natural selection and the idea of similarity. Both concepts involve the evolution of similar adaptive solutions as a result of environmental challenges. The distinction lies in ancestral phenotypes. Parallel evolution takes place when the ancestral phenotypes (before selection) of the lineages are similar. Convergent evolution happens when the lineages have distinct ancestral phenotypes (before selection). Because an ancestral-based distinction will inevitably lead to cases where uncertainty in the distinction may arise, the framework includes a general term, repeated evolution, which can be used as a term applying to the evolution of similar phenotypes and genotypes as well as similar responses to environmental pressures. Based on the argument that genetic similarity may frequently arise without selection, the framework posits that the similarity of genetic sequences is not of great interest unless linked to the actions of natural selection or to the origins (mutation, standing genetic variation, gene flow) and locations of the similar sequences.

Ancient hybridization leads to the repeated evolution of red flowers across a monkeyflower radiation

The reuse of old genetic variation can promote rapid diversification in evolutionary radiations, but in most cases, the historical events underlying this divergence are not known. For example, ancient hybridization can generate new combinations of alleles that sort into descendant lineages, potentially providing the raw material to initiate divergence. In the Mimulus aurantiacus species complex, there is evidence for widespread gene flow among members of this radiation. In addition, allelic variation in the MaMyb2 gene is responsible for differences in flower color between the closely related ecotypes of subspecies puniceus, contributing to reproductive isolation by pollinators. Previous work suggested that MaMyb2 was introgressed into the red-flowered ecotype of puniceus. However, additional taxa within the radiation have independently evolved red flowers from their yellow-flowered ancestors, raising the possibility that this introgression had a more ancient origin. In this study, we used repeated tests of admixture from whole-genome sequence data across this diverse radiation to demonstrate that there has been both ancient and recurrent hybridization in this group. However, most of the signal of this ancient introgression has been removed due to selection, suggesting that widespread barriers to gene flow are in place between taxa. Yet, a roughly 30 kb region that contains the MaMyb2 gene is currently shared only among the red-flowered taxa. Patterns of admixture, sequence divergence, and extended haplotype homozygosity across this region confirm a history of ancient hybridization, where functional variants have been preserved due to positive selection in red-flowered taxa but lost in their yellow-flowered counterparts. The results of this study reveal that selection against gene flow can reduce genomic signatures of ancient hybridization, but that historical introgression can provide essential genetic variation that facilitates the repeated evolution of phenotypic traits between lineages.

The complete and fully-phased diploid genome of a male Han Chinese

Since the release of the complete human genome, the priority of human genomic study has now been shifting towards closing gaps in ethnic diversity. Here, we present a fully phased and well-annotated diploid human genome from a Han Chinese male individual (CN1), in which the assemblies of both haploids achieve the telomere-to-telomere (T2T) level. Comparison of this diploid genome with the CHM13 haploid T2T genome revealed significant variations in the centromere. Outside the centromere, we discovered 11,413 structural variations, including numerous novel ones. We also detected thousands of CN1 alleles that have accumulated high substitution rates and a few that have been under positive selection in the East Asian population. Further, we found that CN1 outperforms CHM13 as a reference genome in mapping and variant calling for the East Asian population owing to the distinct structural variants of the two references. Comparison of SNP calling for a large cohort of 8869 Chinese genomes using CN1 and CHM13 as reference respectively showed that the reference bias profoundly impacts rare SNP calling, with nearly 2 million rare SNPs miss-called with different reference genomes. Finally, applying the CN1 as a reference, we discovered 5.80 Mb and 4.21 Mb putative introgression sequences from Neanderthal and Denisovan, respectively, including many East Asian specific ones undetected using CHM13 as the reference. Our analyses reveal the advances of using CN1 as a reference for population genomic studies and paleo-genomic studies. This complete genome will serve as an alternative reference for future genomic studies on the East Asian population.

Cycles of fusion and fission enabled rapid parallel adaptive radiations in African cichlids

Although some lineages of animals and plants have made impressive adaptive radiations when provided with ecological opportunity, the propensities to radiate vary profoundly among lineages for unknown reasons. In Africa's Lake Victoria region, one cichlid lineage radiated in every lake, with the largest radiation taking place in a lake less than 16,000 years old. We show that all of its ecological guilds evolved in situ. Cycles of lineage fusion through admixture and lineage fission through speciation characterize the history of the radiation. It was jump-started when several swamp-dwelling refugial populations, each of which were of older hybrid descent, met in the newly forming lake, where they fused into a single population, resuspending old admixture variation. Each population contributed a different set of ancient alleles from which a new adaptive radiation assembled in record time, involving additional fusion-fission cycles. We argue that repeated fusion-fission cycles in the history of a lineage make adaptive radiation fast and predictable.

Into the Wild: A novel wild-derived inbred strain resource expands the genomic and phenotypic diversity of laboratory mouse models

The laboratory mouse has served as the premier animal model system for both basic and preclinical investigations for a century. However, laboratory mice capture a narrow subset of the genetic variation found in wild mouse populations. This consideration inherently restricts the scope of potential discovery in laboratory models and narrows the pool of potentially identified phenotype-associated variants and pathways. Wild mouse populations are reservoirs of predicted functional and disease-associated alleles, but the sparsity of commercially available, well-characterized wild mouse strains limits their broader adoption in biomedical research. To overcome this barrier, we have recently imported, sequenced, and phenotyped a set of 11 wild-derived inbred strains developed from wild-caught Mus musculus domesticus. Each of these "Nachman strains" immortalizes a unique wild haplotype sampled from five environmentally diverse locations across North and South America: Saratoga Springs, New York, USA; Gainesville, Florida, USA; Manaus, Brazil; Tucson, Arizona, USA; and Edmonton, Alberta, Canada. Whole genome sequence analysis reveals that each strain carries between 4.73-6.54 million single nucleotide differences relative to the mouse reference assembly, with 42.5% of variants in the Nachman strain genomes absent from classical inbred mouse strains. We phenotyped the Nachman strains on a customized pipeline to assess the scope of disease-relevant neurobehavioral, biochemical, physiological, metabolic, and morphological trait variation. The Nachman strains exhibit significant inter-strain variation in >90% of 1119 surveyed traits and expand the range of phenotypic diversity captured in classical inbred strain panels alone. Taken together, our work introduces a novel wild-derived inbred mouse strain resource that will enable new discoveries in basic and preclinical research. These strains are currently available through The Jackson Laboratory Repository under laboratory code NachJ.

Genomes of cultivated and wild Capsicum species provide insights into pepper domestication and population differentiation

Pepper (Capsicum spp.) is one of the earliest cultivated crops and includes five domesticated species, C. annuum var. annuum, C. chinense, C. frutescens, C. baccatum var. pendulum and C. pubescens. Here, we report a pepper graph pan-genome and a genome variation map of 500 accessions from the five domesticated Capsicum species and close wild relatives. We identify highly differentiated genomic regions among the domesticated peppers that underlie their natural variations in flowering time, characteristic flavors, and unique resistances to biotic and abiotic stresses. Domestication sweeps detected in C. annuum var. annuum and C. baccatum var. pendulum are mostly different, and the common domestication traits, including fruit size, shape and pungency, are achieved mainly through the selection of distinct genomic regions between these two cultivated species. Introgressions from C. baccatum into C. chinense and C. frutescens are detected, including those providing genetic sources for various biotic and abiotic stress tolerances.

Post-hybridization introgression and natural selection promoted genomic divergence of Aegilops speltoides and the four S*-genome diploid species

Understanding how different driving forces have promoted biological divergence and speciation is one of the central issues in evolutionary biology. The Triticum/Aegilops species complex contains 13 diploid species belonging to the A-, B- and D-lineages and offers an ideal system to address the evolutionary dynamics of lineage fusion and splitting. Here, we sequenced the whole genomes of one S-genome species (Aegilops speltoides) of the B-lineage and four S*-genome diploid species (Aegilops bicornis, Aegilops longissima, Aegilops sharonensis and Aegilops searsii) of the D-lineage at the population level. We performed detailed comparisons of the five species and with the other four representative A-, B- and D-lineage species. Our estimates identified frequent genetic introgressions from A- and B-lineages to the D-lineage species. A remarkable observation is the contrasting distributions of putative introgressed loci by the A- and B-lineages along all the seven chromosomes to the extant D-lineage species. These genetic introgressions resulted in high levels of genetic divergence at centromeric regions between Ae. speltoides (B-lineage) and the other four S*-genome diploid species (D-lineage), while natural selection is a potential contributor to divergence among the four S*-genome species at telomeric regions. Our study provides a genome-wide view on how genetic introgression and natural selection acted together yet chromosome-regionally divided to promote genomic divergence among the five S- and S*-genome diploid species, which provides new and nuanced insights into the evolutionary history of the Triticum/Aegilops species complex.

A genomic assessment of the marine-speciation paradox within the toothed whale superfamily Delphinoidea

The impact of post-divergence gene flow in speciation has been documented across a range of taxa in recent years, and may have been especially widespread in highly mobile, wide-ranging marine species, such as cetaceans. Here, we studied individual genomes from nine species across the three families of the toothed whale superfamily Delphinoidea (Delphinidae, Phocoenidae and Monodontidae). To investigate the role of post-divergence gene flow in the speciation process, we used a multifaceted approach, including (i) phylogenomics, (ii) the distribution of shared derived alleles and (iii) demographic inference. We found the divergence of lineages within Delphinoidea did not follow a process of pure bifurcation, but was much more complex. Sliding-window phylogenomics reveal a high prevalence of discordant topologies within the superfamily, with further analyses indicating these discordances arose due to both incomplete lineage sorting and gene flow. D-statistics and f-branch analyses supported gene flow between members of Delphinoidea, with the vast majority of gene flow occurring as ancient interfamilial events. Demographic analyses provided evidence that introgressive gene flow has likely ceased between all species pairs tested, despite reports of contemporary interspecific hybrids. Our study provides the first steps towards resolving the large complexity of speciation within Delphinoidea; we reveal the prevalence of ancient interfamilial gene flow events prior to the diversification of each family, and suggest that contemporary hybridisation events may be disadvantageous, as hybrid individuals do not appear to contribute to the parental species' gene pools.

Range-wide evolutionary relationships and historical demography of brown bears (Ursus arctos) revealed by whole-genome sequencing of isolated central Asian populations

Phylogeographic studies uncover hidden pathways of divergence and inform conservation. Brown bears (Ursus arctos) have one of the broadest distributions of all land mammals, ranging from Eurasia to North America, and are an important model for evolutionary studies. Although several whole genomes were available for individuals from North America, Europe and Asia, limited whole-genome data were available from Central Asia, including the highly imperilled brown bears in the Gobi Desert. To fill this knowledge gap, we sequenced whole genomes from nine Asian brown bears from the Gobi Desert of Mongolia, Northern Mongolia and the Himalayas of Pakistan. We combined these data with published brown bear sequences from Europe, Asia and North America, as well as other bear species. Our goals were to determine the evolutionary relationships among brown bear populations worldwide, their genetic diversity and their historical demography. Our analyses revealed five major lineages of brown bears based on a filtered set of 684,081 single nucleotide polymorphisms. We found distinct evolutionary lineages of brown bears in the Gobi, Himalayas, northern Mongolia, Europe and North America. The lowest level of genetic diversity and the highest level of inbreeding were found in Pakistan, the Gobi Desert and Central Italy. Furthermore, the effective population size (Ne ) for all brown bears decreased over the last 70,000 years. Our results confirm the genetic distinctiveness and ancient lineage of brown bear subspecies in the Gobi Desert of Mongolia and the Himalayas of Pakistan and highlight their importance for conservation.

Interploidy Introgression Shaped Adaptation during the Origin and Domestication History of Brassica napus

Polyploidy is recurrent across the tree of life and known as an evolutionary driving force in plant diversification and crop domestication. How polyploid plants adapt to various habitats has been a fundamental question that remained largely unanswered. Brassica napus is a major crop cultivated worldwide, resulting from allopolyploidy between unknown accessions of diploid B. rapa and B. oleracea. Here, we used whole-genome resequencing data of accessions representing the majority of morphotypes and ecotypes from the species B. rapa, B. oleracea, and B. napus to investigate the role of polyploidy during domestication. To do so, we first reconstructed the phylogenetic history of B. napus, which supported the hypothesis that the emergence of B. napus derived from the hybridization of European turnip of B. rapa and wild B. oleracea. These analyses also showed that morphotypes of swede and Siberian kale (used as vegetable and fodder) were domesticated before rapeseed (oil crop). We next observed that frequent interploidy introgressions from sympatric diploids were prominent throughout the domestication history of B. napus. Introgressed genomic regions were shown to increase the overall genetic diversity and tend to be localized in regions of high recombination. We detected numerous candidate adaptive introgressed regions and found evidence that some of the genes in these regions contributed to phenotypic diversification and adaptation of different morphotypes. Overall, our results shed light on the origin and domestication of B. napus and demonstrate interploidy introgression as an important mechanism that fuels rapid diversification in polyploid species.

Hybridization and Transgressive Evolution Generate Diversity in an Adaptive Radiation of Anolis Lizards

Interspecific hybridization may act as a major force contributing to the evolution of biodiversity. Although generally thought to reduce or constrain divergence between 2 species, hybridization can, paradoxically, promote divergence by increasing genetic variation or providing novel combinations of alleles that selection can act upon to move lineages toward new adaptive peaks. Hybridization may, then, play a key role in adaptive radiation by allowing lineages to diversify into new ecological space. Here, we test for signatures of historical hybridization in the Anolis lizards of Puerto Rico and evaluate 2 hypotheses for the role of hybridization in facilitating adaptive radiation-the hybrid swarm origins hypothesis and the syngameon hypothesis. Using whole genome sequences from all 10 species of Puerto Rican anoles, we calculated D and f-statistics (from ABBA-BABA tests) to test for introgression across the radiation and employed multispecies network coalescent methods to reconstruct phylogenetic networks that allow for hybridization. We then analyzed morphological data for these species to test for patterns consistent with transgressive evolution, a phenomenon in which the trait of a hybrid lineage is found outside of the range of its 2 parents. Our analyses uncovered strong evidence for introgression at multiple stages of the radiation, including support for an ancient hybrid origin of a clade comprising half of the extant Puerto Rican anole species. Moreover, we detected significant signals of transgressive evolution for 2 ecologically important traits, head length and toepad width, the latter of which has been described as a key innovation in Anolis. [Adaptive radiation; introgression; multispecies network coalescent; phenotypic evolution; phylogenetic network; reticulation; syngameon; transgressive segregation.].

A syntelog-based pan-genome provides insights into rice domestication and de-domestication

BACKGROUND

Asian rice is one of the world's most widely cultivated crops. Large-scale resequencing analyses have been undertaken to explore the domestication and de-domestication genomic history of Asian rice, but the evolution of rice is still under debate.

RESULTS

Here, we construct a syntelog-based rice pan-genome by integrating and merging 74 high-accuracy genomes based on long-read sequencing, encompassing all ecotypes and taxa of Oryza sativa and Oryza rufipogon. Analyses of syntelog groups illustrate subspecies divergence in gene presence-and-absence and haplotype composition and identify massive genomic regions putatively introgressed from ancient Geng/japonica to ancient Xian/indica or its wild ancestor, including almost all well-known domestication genes and a 4.5-Mbp centromere-spanning block, supporting a single domestication event in main rice subspecies. Genomic comparisons between weedy and cultivated rice highlight the contribution from wild introgression to the emergence of de-domestication syndromes in weedy rice.

CONCLUSIONS

This work highlights the significance of inter-taxa introgression in shaping diversification and divergence in rice evolution and provides an exploratory attempt by utilizing the advantages of pan-genomes in evolutionary studies.

North African fox genomes show signatures of repeated introgression and adaptation to life in deserts

Elucidating the evolutionary process of animal adaptation to deserts is key to understanding adaptive responses to climate change. Here we generated 82 individual whole genomes of four fox species (genus Vulpes) inhabiting the Sahara Desert at different evolutionary times. We show that adaptation of new colonizing species to a hot arid environment has probably been facilitated by introgression and trans-species polymorphisms shared with older desert resident species, including a putatively adaptive 25 Mb genomic region. Scans for signatures of selection implicated genes affecting temperature perception, non-renal water loss and heat production in the recent adaptation of North African red foxes (Vulpes vulpes), after divergence from Eurasian populations approximately 78 thousand years ago. In the extreme desert specialists, Rueppell's fox (V. rueppellii) and fennec (V. zerda), we identified repeated signatures of selection in genes affecting renal water homeostasis supported by gene expression and physiological differences. Our study provides insights into the mechanisms and genetic underpinnings of a natural experiment of repeated adaptation to extreme conditions.

Introgression between highly divergent fungal sister species

To understand how species evolve and adapt to changing environments, it is important to study gene flow and introgression due to their influence on speciation and radiation events. Here, we apply a novel experimental system for investigating these mechanisms using natural populations. The system is based on two fungal sister species with morphological and ecological similarities occurring in overlapping habitats. We examined introgression between these species by conducting whole genome sequencing of individuals from populations in North America and Europe. We assessed genome-wide nucleotide divergence and performed crossing experiments to study reproductive barriers. We further used ABBA-BABA statistics together with a network analysis to investigate introgression, and conducted demographic modelling to gain insight into divergence times and introgression events. The results revealed that the species are highly divergent and incompatible in vitro. Despite this, small regions of introgression were scattered throughout the genomes and one introgression event likely involves a ghost population (extant or extinct). This study demonstrates that introgression can be found among divergent species and that population histories can be studied without collections of all the populations involved. Moreover, the experimental system is shown to be a useful tool for research on reproductive isolation in natural populations.

Suppressing a phosphohydrolase of cytokinin nucleotide enhances grain yield in rice

One-step and two-step pathways are proposed to synthesize cytokinin in plants. The one-step pathway is mediated by LONELY GUY (LOG) proteins. However, the enzyme for the two-step pathway remains to be identified. Here, we show that quantitative trait locus GY3 may boost grain yield by more than 20% through manipulating a two-step pathway. Locus GY3 encodes a LOG protein that acts as a 5'-ribonucleotide phosphohydrolase by excessively consuming the cytokinin precursors, which contrasts with the activity of canonical LOG members as phosphoribohydrolases in a one-step pathway. The residue S41 of GY3 is crucial for the dephosphorylation of iPRMP to produce iPR. A solo-LTR insertion within the promoter of GY3 suppressed its expression and resulted in a higher content of active cytokinins in young panicles. Introgression of GY302428 increased grain yield per plot by 7.4% to 16.3% in all investigated indica backgrounds, which demonstrates the great value of GY302428 in indica rice production.

Genomic conservation of crop wild relatives: A case study of citrus

Conservation of crop wild relatives is critical for plant breeding and food security. The lack of clarity on the genetic factors that lead to endangered status or extinction create difficulties when attempting to develop concrete recommendations for conserving a citrus wild relative: the wild relatives of crops. Here, we evaluate the conservation of wild kumquat (Fortunella hindsii) using genomic, geographical, environmental, and phenotypic data, and forward simulations. Genome resequencing data from 73 accessions from the Fortunella genus were combined to investigate population structure, demography, inbreeding, introgression, and genetic load. Population structure was correlated with reproductive type (i.e., sexual and apomictic) and with a significant differentiation within the sexually reproducing population. The effective population size for one of the sexually reproducing subpopulations has recently declined to ~1,000, resulting in high levels of inbreeding. In particular, we found that 58% of the ecological niche overlapped between wild and cultivated populations and that there was extensive introgression into wild samples from cultivated populations. Interestingly, the introgression pattern and accumulation of genetic load may be influenced by the type of reproduction. In wild apomictic samples, the introgressed regions were primarily heterozygous, and genome-wide deleterious variants were hidden in the heterozygous state. In contrast, wild sexually reproducing samples carried a higher recessive deleterious burden. Furthermore, we also found that sexually reproducing samples were self-incompatible, which prevented the reduction of genetic diversity by selfing. Our population genomic analyses provide specific recommendations for distinct reproductive types and monitoring during conservation. This study highlights the genomic landscape of a wild relative of citrus and provides recommendations for the conservation of crop wild relatives.