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.

Does selection favour the maintenance of porous species boundaries?

The endpoint of speciation has been viewed as complete isolation and the absence of gene flow between species. If the influx of genes from another species is maladaptive because species have different adaptations and genetic backgrounds, selection should favour the closing of species boundaries and zero gene flow, a process known as reinforcement. Recently, numerous cases of gene flow between species have been identified, many of which involved adaptive introgression of beneficial alleles. These cases could reflect transient states on the way to closed species boundaries or the result of declining strength or efficacy of selection for reinforcement as the level of gene flow approaches zero. An alternative hypothesis, however, is that selection favours porous species boundaries that allow beneficial alleles to cross, especially in changing environments. This perspective evaluates the conditions that would be needed for selection to favour porous species boundaries and the evidence for them. A contrast is made between hybridization in sexual eukaryotes and gene transfer via homologous recombination in bacteria. Current evidence is inconclusive on whether non-zero gene flow is favoured by selection. Studies are needed that quantify selection gradients on rates of gene flow and test for evolution towards intermediate values, especially experiments that manipulate conditions and track evolution for multiple generations. Estimation of gene flow networks for more clades and regional assemblages using emerging genome data will also allow the evolutionary determinants of interspecific gene flow to be better understood.

Museum genomics reveals the hybrid origin of an extinct crater lake endemic

Crater lake fishes are common evolutionary model systems, with recent studies suggesting a key role for gene flow in promoting rapid adaptation and speciation. However, the study of these young lakes can be complicated by human-mediated extinctions. Museum genomics approaches integrating genetic data from recently extinct species are therefore critical to understanding the complex evolutionary histories of these fragile systems. Here, we examine the evolutionary history of an extinct Southern Hemisphere crater lake endemic, the rainbowfish Melanotaenia eachamensis. We undertook comprehensive sampling of extant rainbowfish populations of the Atherton Tablelands of Australia alongside historical museum material to understand the evolutionary origins of the extinct crater lake population and the dynamics of gene flow across the ecoregion. The extinct crater lake species is genetically distinct from all other nearby populations due to historic introgression between two proximate riverine lineages, similar to other prominent crater lake speciation systems, but this historic gene flow has not been sufficient to induce a species flock. Our results suggest that museum genomics approaches can be successfully combined with extant sampling to unravel complex speciation dynamics involving recently extinct species.

Giraffe lineages are shaped by major ancient admixture events

Strong genetic structure has prompted discussion regarding giraffe taxonomy,1,2,3 including a suggestion to split the giraffe into four species: Northern (Giraffa c. camelopardalis), Reticulated (G. c. reticulata), Masai (G. c. tippelskirchi), and Southern giraffes (G. c. giraffa).4,5,6 However, their evolutionary history is not yet fully resolved, as previous studies used a simple bifurcating model and did not explore the presence or extent of gene flow between lineages. We therefore inferred a model that incorporates various evolutionary processes to assess the drivers of contemporary giraffe diversity. We analyzed whole-genome sequencing data from 90 wild giraffes from 29 localities across their current distribution. The most basal divergence was dated to 280 kya. Genetic differentiation, FST, among major lineages ranged between 0.28 and 0.62, and we found significant levels of ancient gene flow between them. In particular, several analyses suggested that the Reticulated lineage evolved through admixture, with almost equal contribution from the Northern lineage and an ancestral lineage related to Masai and Southern giraffes. These new results highlight a scenario of strong differentiation despite gene flow, providing further context for the interpretation of giraffe diversity and the process of speciation in general. They also illustrate that conservation measures need to target various lineages and sublineages and that separate management strategies are needed to conserve giraffe diversity effectively. Given local extinctions and recent dramatic declines in many giraffe populations, this improved understanding of giraffe evolutionary history is relevant for conservation interventions, including reintroductions and reinforcements of existing populations.

Ancestry dynamics and trait selection in a designer cat breed

The Bengal cat breed was developed from intercrosses between the Asian leopard cat, Prionailurus bengalensis, and the domestic cat, Felis catus, with a last common ancestor approximately 6 million years ago. Predicted to derive ∼94% of their genome from domestic cats, regions of the leopard cat genome are thought to account for the unique pelage traits and ornate color patterns of the Bengal breed, which are similar to those of ocelots and jaguars. We explore ancestry distribution and selection signatures in the Bengal breed by using reduced representation and whole-genome sequencing from 947 cats. The mean proportion of leopard cat DNA in the Bengal breed is 3.48%, lower than predicted from breed history, and is broadly distributed, covering 93% of the Bengal genome. Overall, leopard cat introgressions do not show strong signatures of selection across the Bengal breed. However, two popular color traits in Bengal cats, charcoal and pheomelanin intensity, are explained by selection of leopard cat genes whose expression is reduced in a domestic cat background, consistent with genetic incompatibility resulting from hybridization. We characterize several selective sweeps in the Bengal genome that harbor candidate genes for pelage and color pattern and that are associated with domestic, rather than leopard, cat haplotypes. We identify the molecular and phenotypic basis of one selective sweep as reduced expression of the Fgfr2 gene, which underlies glitter, a trait desired by breeders that affects hair texture and light reflectivity.

Museum skins enable identification of introgression associated with cytonuclear discordance

Increased sampling of genomes and populations across closely related species has revealed that levels of genetic exchange during and after speciation are higher than previously thought. One obvious manifestation of such exchange is strong cytonuclear discordance, where the divergence in mitochondrial DNA (mtDNA) differs from that for nuclear genes more (or less) than expected from differences between mtDNA and nuclear DNA (nDNA) in population size and mutation rate. Given genome-scale datasets and coalescent modelling, we can now confidently identify cases of strong discordance and test specifically for historical or recent introgression as the cause. Using population sampling, combining exon capture data from historical museum specimens and recently collected tissues we showcase how genomic tools can resolve complex evolutionary histories in the brachyotis group of rock-wallabies (Petrogale). In particular, applying population and phylogenomic approaches we can assess the role of demographic processes in driving complex evolutionary patterns and assess a role of ancient introgression and hybridisation. We find that described species are well supported as monophyletic taxa for nDNA genes, but not for mtDNA, with cytonuclear discordance involving at least four operational taxonomic units (OTUs) across four species which diverged 183-278 kya. ABC modelling of nDNA gene trees supports introgression during or after speciation for some taxon pairs with cytonuclear discordance. Given substantial differences in body size between the species involved, this evidence for gene flow is surprising. Heterogenous patterns of introgression were identified but do not appear to be associated with chromosome differences between species. These and previous results suggest that dynamic past climates across the monsoonal tropics could have promoted reticulation among related species.

The emerging importance of cross-ploidy hybridisation and introgression

Natural hybridisation is now recognised as pervasive in its occurrence across the Tree of Life. Resurgent interest in natural hybridisation fuelled by developments in genomics has led to an improved understanding of the genetic factors that promote or prevent species cross-mating. Despite this body of work overturning many widely held assumptions about the genetic barriers to hybridisation, it is still widely thought that ploidy differences between species will be an absolute barrier to hybridisation and introgression. Here, we revisit this assumption, reviewing findings from surveys of polyploidy and hybridisation in the wild. In a case study in the British flora, 203 hybrids representing 35% of hybrids with suitable data have formed via cross-ploidy matings, while a wider literature search revealed 59 studies (56 in plants and 3 in animals) in which cross-ploidy hybridisation has been confirmed with genetic data. These results show cross-ploidy hybridisation is readily overlooked, and potentially common in some groups. General findings from these studies include strong directionality of hybridisation, with introgression usually towards the higher ploidy parent, and cross-ploidy hybridisation being more likely to involve allopolyploids than autopolyploids. Evidence for adaptive introgression across a ploidy barrier and cases of cross-ploidy hybrid speciation shows the potential for important evolutionary outcomes.

Sheep Experiment Station: a case study of introgression

Sheep breeders requested that the U.S. Sheep Experiment Station (USSES) to participate in national genetic evaluation through the National Sheep Improvement Program (NSIP). The reasons included the need for (1) a comparison of the productivity of industry and United States Department of Agriculture (USDA) lines, (2) transparency of USDA flocks, (3) genetic ties for NSIP by sampling of industry flocks, and (4) development of premium genetic lines for public release. In response, USSES began to incorporate external sires from NSIP participating flocks into the USSES Targhee flock. Our objective, based on a pedigree analysis, was to test if introgression of external genetics into the flock was achieved. The pedigree included 13,189 animals with mean maximum generations, mean complete generations, and mean equivalent complete generations of 4.2, 1.8, and 2.6, respectively. The mean generation interval was 3.1 yr. The reference population was defined as lambs born from 2021 to 2023 (n = 792). Two additional populations were defined as the current mature ewe flock (n = 123) and the current mature rams (n = 14). The Genetic Conservation Index averaged 7.7 for the full population and 25.7 for the reference population. Overall inbreeding was 0.003 for the full population and 0.006 for the reference population. The rate of inbreeding was 0.0003 per generation. Average relatedness was 0.015 for the full population and 0.018 for the reference population. The effective number of founders, effective number of ancestors, and founder genome equivalents contributing to the reference population were 60, 39, and 19.1, respectively. The ratio of the effective number of founders to the effective number of ancestors was 1.5, indicating the presence of genetic bottlenecks. Measures of effective population size ranged from 102 to 547. Of the 704 offspring produced by external sires, 17 ram lambs and 132 ewe lambs were retained for breeding. The USSES sires produced 299 offspring with 2 ram lambs and 51 ewe lambs retained. Incorporating external sires resulted in a cumulative percentage of genetic variance of 48.8, 49.1, and 44.2 of external genetics for the reference population, current mature ewe flock, and current mature rams, respectively. Stakeholder needs were addressed by introgression of external sires and participation in NSIP, but future selection practices need to be modified to maintain a minimum of 50% USSES core genetics in the flock.

Evolutionary radiation of the Eurasian Pinus species under pervasive gene flow

Evolutionary radiation, a pivotal aspect of macroevolution, offers valuable insights into evolutionary processes. The genus Pinus is the largest genus in conifers withc . $$ c. $$ 90% of the extant species emerged in the Miocene, which signifies a case of rapid diversification. Despite this remarkable history, our understanding of the mechanisms driving radiation within this expansive genus has remained limited. Using exome capture sequencing and a fossil-calibrated phylogeny, we investigated the divergence history, niche diversification, and introgression among 13 closely related Eurasian species spanning climate zones from the tropics to the boreal Arctic. We detected complex introgression among lineages in subsection Pinus at all stages of the phylogeny. Despite this widespread gene exchange, each species maintained its genetic identity and showed clear niche differentiation. Demographic analysis unveiled distinct population histories among these species, which further influenced the nucleotide diversity and efficacy of purifying and positive selection in each species. Our findings suggest that radiation in the Eurasian pines was likely fueled by interspecific recombination and further reinforced by their adaptation to distinct environments. Our study highlights the constraints and opportunities for evolutionary change, and the expectations of future adaptation in response to environmental changes in different lineages.

The temporal and genomic scale of selection following hybridization

Genomic evidence supports an important role for selection in shaping patterns of introgression along the genome, but frameworks for understanding the evolutionary dynamics within hybrid populations that underlie these patterns have been lacking. Due to the clock-like effect of recombination in hybrids breaking up parental haplotypes, drift and selection produce predictable patterns of ancestry variation at varying spatial genomic scales through time. Here, we develop methods based on the Discrete Wavelet Transform to study the genomic scale of local ancestry variation and its association with recombination rates and show that these methods capture temporal dynamics of drift and genome-wide selection after hybridization. We apply these methods to published datasets from hybrid populations of swordtail fish (Xiphophorus) and baboons (Papio) and to inferred Neanderthal introgression in modern humans. Across systems, upward of 20% of variation in local ancestry at the broadest genomic scales can be attributed to systematic selection against introgressed alleles, consistent with strong selection acting on early-generation hybrids. Signatures of selection at fine genomic scales suggest selection over longer time scales; however, we suggest that our ability to confidently infer selection at fine scales is likely limited by inherent biases in current methods for estimating local ancestry from contiguous segments of genomic similarity. Wavelet approaches will become widely applicable as genomic data from systems with introgression become increasingly available and can help shed light on generalities of the genomic consequences of interspecific hybridization.

Evaluating and Predicting the Performance of Sorghum Lines in an Elite by Exotic Backcross-Nested Association Mapping Population

Maintaining or introducing genetic diversity into plant breeding programs is necessary for continual genetic gain; however, diversity at the cost of reduced performance is not something sought by breeders. To this end, backcross-nested association mapping (BC-NAM) populations, in which the recurrent parent is an elite line, can be employed as a strategy to introgress diversity from unadapted accessions while maintaining agronomic performance. This study evaluates (i) the hybrid performance of sorghum lines from 18 BC1-NAM families and (ii) the potential of genomic prediction to screen lines from BC1-NAM families for hybrid performance prior to phenotypic evaluation. Despite the diverse geographical origins and agronomic performance of the unadapted parents for BC1-NAM families, many BC1-derived lines performed significantly better in the hybrid trials than the elite recurrent parent, R.Tx436. The genomic prediction accuracies for grain yield, plant height, and days to mid-anthesis were acceptable, but the prediction accuracies for plant height were lower than expected. While the prediction accuracies increased when including more individuals in the training set, improvements tended to plateau between two and five lines per family, with larger training sets being required for more complex traits such as grain yield. Therefore, genomic prediction models can be optimized in a large BC1-NAM population with a relatively low fraction of individuals needing to be evaluated. These results suggest that genomic prediction is an effective method of pre-screening lines within BC1-NAM families prior to evaluation in extensive hybrid field trials.

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.

Chromosomal Inversions and the Demography of Speciation in Drosophila montana and Drosophila flavomontana

Chromosomal inversions may play a central role in speciation given their ability to locally reduce recombination and therefore genetic exchange between diverging populations. We analyzed long- and short-read whole-genome data from sympatric and allopatric populations of 2 Drosophila virilis group species, Drosophila montana and Drosophila flavomontana, to understand if inversions have contributed to their divergence. We identified 3 large alternatively fixed inversions on the X chromosome and one on each of the autosomes 4 and 5. A comparison of demographic models estimated for inverted and noninverted (colinear) chromosomal regions suggests that these inversions arose before the time of the species split. We detected a low rate of interspecific gene flow (introgression) from D. montana to D. flavomontana, which was further reduced inside inversions and was lower in allopatric than in sympatric populations. Together, these results suggest that the inversions were already present in the common ancestral population and that gene exchange between the sister taxa was reduced within inversions both before and after the onset of species divergence. Such ancestrally polymorphic inversions may foster speciation by allowing the accumulation of genetic divergence in loci involved in adaptation and reproductive isolation inside inversions early in the speciation process, while gene exchange at colinear regions continues until the evolving reproductive barriers complete speciation. The overlapping X inversions are particularly good candidates for driving the speciation process of D. montana and D. flavomontana, since they harbor strong genetic incompatibilities that were detected in a recent study of experimental introgression.

An indel introduced by Neanderthal introgression, rs3835124:ATTTATT > ATT, might contribute to prostate cancer risk by regulating PDK1 expression

INTRODUCTION

Prostate cancer is one of the most common cancer types in males and rs12621278:A > G has been suggested to be associated with this disease by previous genome-wide association studies. One thousand genomes project data analysis indicated that rs12621278:A > G is within two long-core haplotypes. However, the origin, causal variant(s), and molecular function of these haplotypes were remaining unclear.

MATERIALS AND METHODS

Population genetics analysis and functional genomics work was performed for this locus.

RESULTS

Phylogeny analysis verified that the rare haplotype is derived from Neanderthal introgression. Genome annotation suggested that three genetic variants in the core haplotypes, rs116108611:G > A, rs139972066:AAAAAAAA > AAAAAAAAA, and rs3835124:ATTTATT > ATT, are located in functional regions. Luciferase assay indicated that rs139972066:AAAAAAAA > AAAAAAAAA and rs116108611:G > A are not able to alter ITGA6 (integrin alpha 6) and ITGA6 antisense RNA 1 expression, respectively. In contrast, rs3835124:ATTTATT > ATT can significantly influence PDK1 (pyruvate dehydrogenase kinase 1) expression, which was verified by expression quantitative trait locus analysis. This genetic variant can alter transcription factor cut like homeobox 1 interaction efficiency. The introgressed haplotype was observed to be subject to positive selection in East Asian populations. The molecular function of the haplotype suggested that Neanderthal should be with lower PDK1 expression and further different energy homeostasis from modern human.

CONCLUSION

This study provided new insight into the contribution of Neanderthal introgression to human phenotypes.

Reproductive isolation in a three-way contact zone

Contact zones between divergent forms within a species provide insight into the role of gene flow in adaptation and speciation. Previous work has focused on contact zones involving only two divergent forms, but in nature, many more than two populations may overlap simultaneously and experience gene flow. Patterns of introgression in wild populations are, therefore, likely much more complicated than is often assumed. We begin to address this gap in current knowledge by investigating patterns of divergence and introgression across a complex natural contact zone. We use phenotypic and genomic data to confirm the existence of a three-way contact zone among divergent freshwater resident, saltwater resident and saltwater migratory three-spined stickleback (Gasterosteus aculeatus) on the island of North Uist, Scottish Western Isles. We find evidence for hybridization, mostly between saltwater resident and saltwater migratory forms. Despite hybridization, genomic analyses reveal pairwise islands of divergence between all forms that are maintained across the contact zone. Genomic cline analyses also provide evidence for selection and/or hybrid incompatibilities in divergent regions. Divergent genomic regions occur across multiple chromosomes and involve many known adaptive loci and several chromosomal inversions. We also identify distinct immune gene expression profiles between forms, but no evidence for transgressive expression in hybrids. Our results suggest that reproductive isolation is maintained in this three-way contact zone, despite some hybridization, and that reduced recombination in chromosomal inversions may play an important role in maintaining this isolation.

MAST: Phylogenetic Inference with Mixtures Across Sites and Trees

Hundreds or thousands of loci are now routinely used in modern phylogenomic studies. Concatenation approaches to tree inference assume that there is a single topology for the entire dataset, but different loci may have different evolutionary histories due to incomplete lineage sorting, introgression, and/or horizontal gene transfer; even single loci may not be treelike due to recombination. To overcome this shortcoming, we introduce an implementation of a multi-tree mixture model that we call MAST. This model extends a prior implementation by Boussau et al. (2009) by allowing users to estimate the weight of each of a set of pre-specified bifurcating trees in a single alignment. The MAST model allows each tree to have its own weight, topology, branch lengths, substitution model, nucleotide or amino acid frequencies, and model of rate heterogeneity across sites. We implemented the MAST model in a maximum-likelihood framework in the popular phylogenetic software, IQ-TREE. Simulations show that we can accurately recover the true model parameters, including branch lengths and tree weights for a given set of tree topologies, under a wide range of biologically realistic scenarios. We also show that we can use standard statistical inference approaches to reject a single-tree model when data are simulated under multiple trees (and vice versa). We applied the MAST model to multiple primate datasets and found that it can recover the signal of incomplete lineage sorting in the Great Apes, as well as the asymmetry in minor trees caused by introgression among several macaque species. When applied to a dataset of four Platyrrhine species for which standard concatenated maximum likelihood and gene tree approaches disagree, we observe that MAST gives the highest weight (i.e. the largest proportion of sites) to the tree also supported by gene tree approaches. These results suggest that the MAST model is able to analyse a concatenated alignment using maximum likelihood, while avoiding some of the biases that come with assuming there is only a single tree. We discuss how the MAST model can be extended in the future.

The interplay of local adaptation and gene flow may lead to the formation of supergenes

Supergenes are genetic architectures resulting in the segregation of alternative combinations of alleles underlying complex phenotypes. The co-segregation of alleles at linked loci is often facilitated by polymorphic chromosomal rearrangements suppressing recombination locally. Supergenes are involved in many complex polymorphisms, including sexual, colour or behavioural polymorphisms in numerous plants, fungi, mammals, fish, and insects. Despite a long history of empirical and theoretical research, the formation of supergenes remains poorly understood. Here, using a two-island population genetic model, we explore how gene flow and the evolution of overdominant chromosomal inversions may jointly lead to the formation of supergenes. We show that the evolution of inversions in differentiated populations, both under disruptive selection, leads to an increase in frequency of poorly adapted, immigrant haplotypes. Indeed, rare allelic combinations, such as immigrant haplotypes, are more frequently reshuffled by recombination than common allelic combinations, and therefore benefit from the recombination suppression generated by inversions. When an inversion capturing a locally adapted haplotype spreads but is associated with a fitness cost hampering its fixation (e.g. a recessive mutation load), the maintenance of a non-inverted haplotype in the population is enhanced; under certain conditions, the immigrant haplotype persists alongside the inverted local haplotype, while the standard local haplotype disappears. This establishes a stable, local polymorphism with two non-recombining haplotypes encoding alternative adaptive strategies, that is, a supergene. These results bring new light to the importance of local adaptation, overdominance, and gene flow in the formation of supergenes and inversion polymorphisms in general.

A CRISPR/Cas9-based vector system enables the fast breeding of selection-marker-free canola with Rcr1-rendered clubroot resistance

Breeding for disease resistance in major crops is of crucial importance for global food security and sustainability. However, common biotechnologies such as traditional transgenesis or genome editing do not provide an ideal solution, whereas transgenic crops free of selection markers such as cisgenic/intragenic crops might be suitable. In this study, after cloning and functional verification of the Rcr1 gene for resistance to clubroot (Plasmodiophora brassicae), we confirmed that the genes Rcr1, Rcr2, Rcr4, and CRa from Brassica rapa crops and the resistance gene from B. napus oilseed rape cv. 'Mendel' on chromosome A03 were identical in their coding regions. We also determined that Rcr1 has a wide distribution in Brassica breeding materials and renders potent resistance against multiple representative clubroot strains in Canada. We then modified a CRISPR/Cas9-based cisgenic vector system and found that it enabled the fast breeding of selection-marker-free transgenic crops with add-on traits, with selection-marker-free canola (B. napus) germplasms with Rcr1-rendered stable resistance to clubroot disease being successfully developed within 2 years. In the B. napus background, the intragenic vector system was able to remove unwanted residue sequences from the final product with high editing efficiency, and off-target mutations were not detected. Our study demonstrates the potential of applying this breeding strategy to other crops that can be transformed by Agrobacterium. Following the streamlined working procedure, intragenic germplasms can be developed within two generations, which could significantly reduce the breeding time and labor compared to traditional introgression whilst still achieving comparable or even better breeding results.

Local environments, not invasive hybridization, influence cardiac performance of native trout under acute thermal stress

Climate-induced expansion of invasive hybridization (breeding between invasive and native species) poses a significant threat to the persistence of many native species worldwide. In the northern U.S. Rocky Mountains, hybridization between native cutthroat trout and non-native rainbow trout has increased in recent decades due, in part, to climate-driven increases in water temperature. It has been postulated that invasive hybridization may enhance physiological tolerance to climate-induced thermal stress because laboratory studies indicate that rainbow trout have a higher thermal tolerance than cutthroat trout. Here, we assessed whether invasive hybridization improves cardiac performance response to acute water temperature stress of native wild trout populations. We collected trout from four streams with a wide range of non-native admixture among individuals and with different temperature and streamflow regimes in the upper Flathead River drainage, USA. We measured individual cardiac performance (maximum heart rate, "MaxHR", and temperature at arrhythmia, "ArrTemp") during laboratory trials with increasing water temperatures (10-28°C). Across the study populations, we observed substantial variation in cardiac performance of individual trout when exposed to thermal stress. Notably, we found significant differences in the cardiac response to thermal regimes among native cutthroat trout populations, suggesting the importance of genotype-by-environment interactions in shaping the physiological performance of native cutthroat trout. However, rainbow trout admixture had no significant effect on cardiac performance (MaxHR and ArrTemp) within any of the three populations. Our results indicate that invasive hybridization with a warmer-adapted species does not enhance the cardiac performance of native trout under warming conditions. Maintaining numerous populations across thermally and hydrologically diverse stream environments will be crucial for native trout to adapt and persist in a warming climate.

Postzygotic barriers persist despite ongoing introgression in hybridizing Mimulus species

The evolution of postzygotic isolation is thought to be a key step in maintaining species boundaries upon secondary contact, yet the dynamics and persistence of hybrid incompatibilities in naturally hybridizing species are not well understood. Here, we explore these issues using genetic mapping in three independent populations of recombinant inbred lines between naturally hybridizing monkeyflowers, Mimulus guttatus and Mimulus nasutus, from the sympatric Catherine Creek population. We discover that the three M. guttatus founders differ dramatically in admixture history, with nearly a quarter of one founder's genome introgressed from M. nasutus. Comparative genetic mapping in the three RIL populations reveals three new putative inversions, each one segregating among the M. guttatus founders, two due to admixture. We find strong, genome-wide transmission ratio distortion in all RILs, but patterns are highly variable among the three populations. At least some of this distortion appears to be explained by epistatic selection favouring parental genotypes, but tests of inter-chromosomal linkage disequilibrium also reveal multiple candidate Dobzhansky-Muller incompatibilities. We also map several genetic loci for hybrid pollen viability, including two interacting pairs that coincide with peaks of distortion. Remarkably, even with this limited sample of three M. guttatus lines, we discover abundant segregating variation for hybrid incompatibilities with M. nasutus, suggesting this population harbours diverse contributors to postzygotic isolation. Moreover, even with substantial admixture, hybrid incompatibilities between Mimulus species persist, suggesting postzygotic isolation might be a potent force in maintaining species barriers in this system.

Morphological and genetic evidence suggest gene flow among native and naturalized mint species

PREMISE

Cultivation and naturalization of plants beyond their natural range can bring previously geographically isolated taxa together, increasing the opportunity for hybridization, the outcomes of which are not predictable. Here, we explored the phenotypic and genomic effects of interspecific gene flow following the widespread cultivation of Mentha spicata (spearmint), M. longifolia, and M. suaveolens.

METHODS

We morphologically evaluated 155 herbarium specimens of three Mentha species and sequenced the genomes of a subset of 93 specimens. We analyzed the whole genomes in a population and the phylogenetic framework and associated genomic classifications in conjunction with the morphological assessments.

RESULTS

The allopolyploid M. spicata, which likely evolved in cultivation, had altered trichome characters, that is possibly a product of human selection for a more palatable plant or a byproduct of selection for essential oils. There were signs of genetic admixture between mints, including allopolyploids, indicating that the reproductive barriers between Mentha species with differences in ploidy are likely incomplete. Still, despite gene flow between species, we found that genetic variants associated with the cultivated trichome morphology continue to segregate.

CONCLUSIONS

Although hybridization, allopolyploidization, and human selection during cultivation can increase species richness (e.g., by forming hybrid taxa), we showed that unless reproductive barriers are strong, these processes can also result in mixing of genes between species and the potential loss of natural biodiversity.

The absence of bumblebees on an oceanic island blurs the species boundary of two closely related orchids

Oceanic islands offer valuable natural laboratories for studying evolution. The Izu Islands, with their recent geological origin, provide an exceptional opportunity to explore the initial evolution on oceanic islands. Another noteworthy aspect is the absence of bumblebee species on most Izu Islands. We used ecological, morphological, and molecular data to investigate the impact of bumblebee absence on the evolution of two closely related orchid species, Goodyera henryi and Goodyera similis, focusing on Kozu Island, the Izu Islands. Our investigation revealed that while G. henryi exclusively relies on a bumblebee species for pollination on the mainland, G. similis is pollinated by scoliid wasps on both the mainland and the island. Intriguingly, all specimens initially categorized as G. henryi on Kozu Island are hybrids of G. henryi and G. similis, leading to the absence of pure G. henryi distribution on the island. These hybrids are pollinated by the scoliid wasp species that also pollinates G. similis on the island. The absence of bumblebees might result in sporadic and inefficient pollination of G. henryi by scoliid wasps, consequently promoting hybrid proliferation on the island. Our findings suggest that the absence of bumblebees can blur plant species boundaries.

Introgression and incomplete lineage sorting blurred phylogenetic relationships across the genomes of sclerophyllous oaks from southwest China

Resolving evolutionary relationships among closely related species with interspecific gene flow is challenging. Genome-scale data provide opportunities to clarify complex evolutionary relationships in closely related species and to observe variations in species relationships across the genomes of such species. The Himalayan-Hengduan subalpine oaks have a nearly completely sympatric distribution in southwest China and probably constitute a syngameon. In this study, we mapped resequencing data from different species in this group to the Quercus aquifolioides reference genome to obtain a high-quality filtered single nucleotide polymorphism (SNP) dataset. We also assembled their plastomes. We reconstructed their phylogenetic relationships, explored the level and pattern of introgression among these species and investigated gene tree variation in the genomes of these species using sliding windows. The same or closely related plastomes were found to be shared extensively among different species within a specific geographical area. Phylogenomic analyses of genome-wide SNP data found that most oaks in the Himalayan-Hengduan subalpine clade showed genetic coherence, but several species were found to be connected by introgression. The gene trees obtained using sliding windows showed that the phylogenetic relationships in the genomes of oaks are highly heterogeneous and therefore highly obscured. Our study found that all the oaks of the Himalayan-Hengduan subalpine clade from southwest China form a syngameon. The obscured phylogenetic relationships observed empirically across the genome are best explained by interspecific gene flow in conjunction with incomplete lineage sorting.

The genome of the black-footed cat: Revealing a rich natural history and urgent conservation priorities for small felids

Habitat degradation and loss of genetic diversity are common threats faced by almost all of today's wild cats. Big cats, such as tigers and lions, are of great concern and have received considerable conservation attention through policies and international actions. However, knowledge of and conservation actions for small wild cats are lagging considerably behind. The black-footed cat, Felis nigripes, one of the smallest felid species, is experiencing increasing threats with a rapid reduction in population size. However, there is a lack of genetic information to assist in developing effective conservation actions. A de novo assembly of a high-quality chromosome-level reference genome of the black-footed cat was made, and comparative genomics and population genomics analyses were carried out. These analyses revealed that the most significant genetic changes in the evolution of the black-footed cat are the rapid evolution of sensory and metabolic-related genes, reflecting genetic adaptations to its characteristic nocturnal hunting and a high metabolic rate. Genomes of the black-footed cat exhibit a high level of inbreeding, especially for signals of recent inbreeding events, which suggest that they may have experienced severe genetic isolation caused by habitat fragmentation. More importantly, inbreeding associated with two deleterious mutated genes may exacerbate the risk of amyloidosis, the dominant disease that causes mortality of about 70% of captive individuals. Our research provides comprehensive documentation of the evolutionary history of the black-footed cat and suggests that there is an urgent need to investigate genomic variations of small felids worldwide to support effective conservation actions.

Artifactual Orthologs and the Need for Diligent Data Exploration in Complex Phylogenomic Datasets: A museomic case study from the Andean flora

The Andes mountains of western South America are a globally important biodiversity hotspot, yet there is a paucity of resolved phylogenies for plant clades from this region. Filling an important gap to our understanding of the World's richest flora, we present the first phylogeny of Freziera (Pentaphylacaceae), an Andean-centered, cloud forest radiation. Our dataset was obtained via yrid-enriched target sequence capture of Angiosperms353 universal loci for 50 of the ca. 75 spp., obtained almost entirely from herbarium specimens. We identify high phylogenomic complexity in Freziera, including the presence of data artifacts. Via by-eye observation of gene trees, detailed examination of warnings from recently improved assembly pipelines, and gene tree filtering, we identified that artifactual orthologs (i.e., the presence of only one copy of a multi-copy gene due to differential assembly) were an important source of gene tree heterogeneity that had a negative impact on phylogenetic inference and support. These artifactual orthologs may be common in plant phylogenomic datasets, where multiple instances of genome duplication are common. After accounting for artifactual orthologs as source of gene tree error, we identified a significant, but non-specific signal of introgression using Patterson's D and f4 statistics. Despite phylogenomic complexity, we were able to resolve Freziera into nine well-supported subclades whose evolution has been shaped by multiple evolutionary processes, including incomplete lineage sorting, historical gene flow, and gene duplication. Our results highlight the complexities of plant phylogenomics, which are heightened in Andean radiations, and show the impact of filtering data processing artifacts and standard filtering approaches on phylogenetic inference.

Reticulate evolution of the tertiary relict Osmanthus

When interspecific gene flow is common, species relationships are more accurately represented by a phylogenetic network than by a bifurcating tree. This study aimed to uncover the role of introgression in the evolution of Osmanthus, the only genus of the subtribe Oleinae (Oleaceae) with its distribution center in East Asia. We built species trees, detected introgression, and constructed networks using multiple kinds of sequencing data (whole genome resequencing, transcriptome sequencing, and Sanger sequencing of nrDNA) combined with concatenation and coalescence approaches. Then, based on well-understood species relationships, historical biogeographic analyses and diversification rate estimates were employed to reveal the history of Osmanthus. Osmanthus originated in mid-Miocene Europe and dispersed to the eastern Tibetan Plateau in the late Miocene. Thereafter, it continued to spread eastwards. Phylogenetic conflict is common within the 'Core Osmanthus' clade and is seen at both early and late stages of diversification, leading to hypotheses of net-like species relationships. Incomplete lineage sorting proved ineffective in explaining phylogenetic conflicts and thus supported introgression as the main cause of conflicts. This study elucidates the diversification history of a relict genus in the subtropical regions of eastern Asia and reveals that introgression had profound effects on its evolutionary history.

Unidirectional trans-Atlantic gene flow and a mixed spawning area shape the genetic connectivity of Atlantic bluefin tuna

The commercially important Atlantic bluefin tuna (Thunnus thynnus), a large migratory fish, has experienced notable recovery aided by accurate resource assessment and effective fisheries management efforts. Traditionally, this species has been perceived as consisting of eastern and western populations, spawning respectively in the Mediterranean Sea and the Gulf of Mexico, with mixing occurring throughout the Atlantic. However, recent studies have challenged this assumption by revealing weak genetic differentiation and identifying a previously unknown spawning ground in the Slope Sea used by Atlantic bluefin tuna of uncertain origin. To further understand the current and past population structure and connectivity of Atlantic bluefin tuna, we have assembled a unique dataset including thousands of genome-wide single-nucleotide polymorphisms (SNPs) from 500 larvae, young of the year and spawning adult samples covering the three spawning grounds and including individuals of other Thunnus species. Our analyses support two weakly differentiated but demographically connected ancestral populations that interbreed in the Slope Sea. Moreover, we also identified signatures of introgression from albacore (Thunnus alalunga) into the Atlantic bluefin tuna genome, exhibiting varied frequencies across spawning areas, indicating strong gene flow from the Mediterranean Sea towards the Slope Sea. We hypothesize that the observed genetic differentiation may be attributed to increased gene flow caused by a recent intensification of westward migration by the eastern population, which could have implications for the genetic diversity and conservation of western populations. Future conservation efforts should consider these findings to address potential genetic homogenization in the species.

The Artefactual Branch Effect and Phylogenetic Conflict: Species Delimitation with Gene Flow in Mangrove Pit Vipers (Trimeresurus purpureomaculatus-erythrurus Complex)

Mangrove pit vipers of the Trimeresurus purpureomaculatus-erythrurus complex are the only species of viper known to naturally inhabit mangroves. Despite serving integral ecological functions in mangrove ecosystems, the evolutionary history, distribution, and species boundaries of mangrove pit vipers remain poorly understood, partly due to overlapping distributions, confusing phenotypic variations, and the lack of focused studies. Here, we present the first genomic study on mangrove pit vipers and introduce a robust hypothesis-driven species delimitation framework that considers gene flow and phylogenetic uncertainty in conjunction with a novel application of a new class of speciation-based delimitation model implemented through the program Delineate. Our results showed that gene flow produced phylogenetic conflict in our focal species and substantiates the artefactual branch effect where highly admixed populations appear as divergent nonmonophyletic lineages arranged in a stepwise manner at the basal position of clades. Despite the confounding effects of gene flow, we were able to obtain unequivocal support for the recognition of a new species based on the intersection and congruence of multiple lines of evidence. This study demonstrates that an integrative hypothesis-driven approach predicated on the consideration of multiple plausible evolutionary histories, population structure/differentiation, gene flow, and the implementation of a speciation-based delimitation model can effectively delimit species in the presence of gene flow and phylogenetic conflict.

Major patterns in the introgression history of Heliconius butterflies

Gene flow between species, although usually deleterious, is an important evolutionary process that can facilitate adaptation and lead to species diversification. It also makes estimation of species relationships difficult. Here, we use the full-likelihood multispecies coalescent (MSC) approach to estimate species phylogeny and major introgression events in Heliconius butterflies from whole-genome sequence data. We obtain a robust estimate of species branching order among major clades in the genus, including the 'melpomene-silvaniform' group, which shows extensive historical and ongoing gene flow. We obtain chromosome-level estimates of key parameters in the species phylogeny, including species divergence times, present-day and ancestral population sizes, as well as the direction, timing, and intensity of gene flow. Our analysis leads to a phylogeny with introgression events that differ from those obtained in previous studies. We find that Heliconius aoede most likely represents the earliest-branching lineage of the genus and that 'silvaniform' species are paraphyletic within the melpomene-silvaniform group. Our phylogeny provides new, parsimonious histories for the origins of key traits in Heliconius, including pollen feeding and an inversion involved in wing pattern mimicry. Our results demonstrate the power and feasibility of the full-likelihood MSC approach for estimating species phylogeny and key population parameters despite extensive gene flow. The methods used here should be useful for analysis of other difficult species groups with high rates of introgression.

Phylogenomic analyses reveal reticulate evolution between Neomicrocalamus and Temochloa (Poaceae: Bambusoideae)

Neomicrocalamus and Temochloa are closely related to bamboo genera. However, when considered with newly discovered and morphologically similar material from China and Vietnam, the phylogenetic relationship among these three groups was ambiguous in the analyses based on DNA regions. Here, as a means of investigating the relationships among the three bamboo groups and exploring potential sources of genomic conflicts, we present a phylogenomic examination based on the whole plastome, single-nucleotide polymorphism (SNP), and single-copy nuclear (SCN) gene datasets. Three different phylogenetic hypotheses were found. The inconsistency is attributed to the combination of incomplete lineage sorting and introgression. The origin of newly discovered bamboos is from introgressive hybridization between Temochloa liliana (which contributed 80.7% of the genome) and Neomicrocalamus prainii (19.3%), indicating that the newly discovered bamboos are closer to T. liliana in genetics. The more similar morphology and closer distribution elevation also imply a closer relationship between Temochloa and newly discovered bamboos.

Evolutionarily labile dispersal behavior and discontinuous habitats enhance population differentiation in island versus continentally distributed swallows

The causes of population divergence in vagile groups remain a paradox in evolutionary biology: dispersive species should be able to colonize new areas, a prerequisite for allopatric speciation, but dispersal also facilitates gene flow, which erodes population differentiation. Strong dispersal ability has been suggested to enhance divergence in patchy habitats and inhibit divergence in continuous landscapes, but empirical support for this hypothesis is lacking. Here we compared patterns of population divergence in a dispersive clade of swallows distributed across both patchy and continuous habitats. The Pacific Swallow (Hirundo tahitica) has an insular distribution throughout Southeast Asia and the Pacific, while its sister species, the Welcome Swallow (H. neoxena), has a continental distribution in Australia. We used whole-genome data to demonstrate strong genetic structure and limited introgression among insular populations, but not among continental populations. Demographic models show that historic changes in habitat connectivity have contributed to population structure within the clade. Swallows appear to exhibit evolutionarily labile dispersal behavior in which they reduce dispersal propensity after island colonization despite retaining strong flight ability. Our data support the hypothesis that fragmented habitats enhance population differentiation in vagile groups, and suggest that labile dispersal behavior is a key mechanism underlying this pattern.

Inference of Admixture Origins in Indigenous African Cattle

Present-day African cattle retain a unique genetic profile composed of a mixture of the Bos taurus and Bos indicus populations introduced into the continent at different time periods. However, details of the admixture history and the exact origins of the source populations remain obscure. Here, we infer the source of admixture in the earliest domestic cattle in Africa, African taurine. We detect a significant contribution (up to ∼20%) from a basal taurine lineage, which might represent the now-extinct African aurochs. In addition, we show that the indicine ancestry of African cattle, although most closely related to so-far sampled North Indian indicine breeds, has a small amount of additional genetic affinity to Southeast Asian indicine breeds. Our findings support the hypothesis of aurochs introgression into African taurine and generate a novel hypothesis that the origin of indicine ancestry in Africa might be different indicine populations than the ones found in North India today.

The role of conflict in shaping plant biodiversity

Although intrinsic postzygotic reproductive barriers can play a fundamental role in speciation, their underlying evolutionary causes are widely debated. One hypothesis is that incompatibilities result from genomic conflicts. Here, I synthesize the evidence that conflict generates incompatibilities in plants, thus playing a creative role in plant biodiversity. While much evidence supports a role for conflict in several classes of incompatibility, integrating knowledge of incompatibility alleles with natural history can provide further essential tests. Moreover, comparative work can shed light on the relative importance of conflict in causing incompatibilities, including the extent to which their evolution is repeatable. Together, these approaches can provide independent lines of evidence that conflict causes incompatibilities, cementing its role in plant speciation.

Hybridization in agricultural weeds: A review from ecological, evolutionary, and management perspectives

Agricultural weeds frequently hybridize with each other or with related crop species. Some hybrid weeds exhibit heterosis (hybrid vigor), which may be stabilized through mechanisms like genome duplication or vegetative reproduction. Even when heterosis is not stabilized, hybridization events diversify weed gene pools and often enable adaptive introgression. Consequently, hybridization may promote weed evolution and exacerbate weed-crop competition. However, hybridization does not always increase weediness. Even when viable and fertile, hybrid weeds sometimes prove unsuccessful in crop fields. This review provides an overview of weed hybridization and its management implications. We describe intrinsic and extrinsic factors that influence hybrid fitness in agroecosystems. We also survey the rapidly growing literature on crop-weed hybridization and the link between hybridization and invasiveness. These topics are increasingly relevant in this era of genetic tools for crop improvement, intensive and simplified cropping systems, and globalized trade. The review concludes with suggested research priorities, including hybridization in the context of climate change, plant-insect interactions, and redesigned weed management programs. From a weed management perspective, hybridization is one of many reasons that researchers and land managers must diversify their weed control toolkits.

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.

Extensive Phylogenomic Discordance and the Complex Evolutionary History of the Neotropical Cat Genus Leopardus

Even in the genomics era, the phylogeny of Neotropical small felids comprised in the genus Leopardus remains contentious. We used whole-genome resequencing data to construct a time-calibrated consensus phylogeny of this group, quantify phylogenomic discordance, test for interspecies introgression, and assess patterns of genetic diversity and demographic history. We infer that the Leopardus radiation started in the Early Pliocene as an initial speciation burst, followed by another in its subgenus Oncifelis during the Early Pleistocene. Our findings challenge the long-held notion that ocelot (Leopardus pardalis) and margay (L. wiedii) are sister species and instead indicate that margay is most closely related to the enigmatic Andean cat (L. jacobita), whose whole-genome data are reported here for the first time. In addition, we found that the newly sampled Andean tiger cat (L. tigrinus pardinoides) population from Colombia associates closely with Central American tiger cats (L. tigrinus oncilla). Genealogical discordance was largely attributable to incomplete lineage sorting, yet was augmented by strong gene flow between ocelot and the ancestral branch of Oncifelis, as well as between Geoffroy's cat (L. geoffroyi) and southern tiger cat (L. guttulus). Contrasting demographic trajectories have led to disparate levels of current genomic diversity, with a nearly tenfold difference in heterozygosity between Andean cat and ocelot, spanning the entire range of variability found in extant felids. Our analyses improved our understanding of the speciation history and diversity patterns in this felid radiation, and highlight the benefits to phylogenomic inference of embracing the many heterogeneous signals scattered across the genome.

Sampling affects population genetic inference: A case study of the Allen's (Selasphorus sasin) and rufous hummingbird (Selasphorus rufus)

Gene flow can affect evolutionary inference when species are undersampled. Here, we evaluate the effects of gene flow and geographic sampling on demographic inference of 2 hummingbirds that hybridize, Allen's hummingbird (Selasphorus sasin) and rufous hummingbird (Selasphorus rufus). Using whole-genome data and extensive geographic sampling, we find widespread connectivity, with introgression far beyond the Allen's × rufous hybrid zone, although the Z chromosome resists introgression beyond the hybrid zone. We test alternative hypotheses of speciation history of Allen's, rufous, and Calliope (S. calliope) hummingbird and find that rufous hummingbird is the sister taxon to Allen's hummingbird, and Calliope hummingbird is the outgroup. A model treating the 2 subspecies of Allen's hummingbird as a single panmictic population fit observed genetic data better than models treating the subspecies as distinct populations, in contrast to morphological and behavioral differences and analyses of spatial population structure. With additional sampling, our study builds upon recent studies that came to conflicting conclusions regarding the evolutionary histories of these 2 species. Our results stress the importance of thorough geographic sampling when assessing demographic history in the presence of gene flow.

Genetic swamping of the critically endangered Scottish wildcat was recent and accelerated by disease

The European wildcat population in Scotland is considered critically endangered as a result of hybridization with introduced domestic cats,1,2 though the time frame over which this gene flow has taken place is unknown. Here, using genome data from modern, museum, and ancient samples, we reconstructed the trajectory and dated the decline of the local wildcat population from viable to severely hybridized. We demonstrate that although domestic cats have been present in Britain for over 2,000 years,3 the onset of hybridization was only within the last 70 years. Our analyses reveal that the domestic ancestry present in modern wildcats is markedly over-represented in many parts of the genome, including the major histocompatibility complex (MHC). We hypothesize that introgression provides wildcats with protection against diseases harbored and introduced by domestic cats, and that this selection contributes to maladaptive genetic swamping through linkage drag. Using the case of the Scottish wildcat, we demonstrate the importance of local ancestry estimates to both understand the impacts of hybridization in wild populations and support conservation efforts to mitigate the consequences of anthropogenic and environmental change.

Fish biogeography and hybridization: do contemporary distributions predict introgression history?

Freshwater ecosystems frequently house diverse assemblages of closely related fish taxa, which can be particularly prone to hybridization and introgression. While extensive introgression may be expected among biogeographically proximate lineages, recent analyses imply that contemporary distributions do not always accurately predict hybridization history. Here, we use the ABBA-BABA approach to test biogeographic hypotheses regarding the extent of hybridization in the recent evolution of New Zealand's species-rich freshwater Galaxias vulgaris fish complex. Genome-wide comparisons reveal significant increases in introgression associated with increasing geographic overlap of taxa. The estimator DP, which assesses the net proportion of a genome originating from introgression, shows a particularly strong relationship with biogeographic overlap (R2 = .43; p = .005). Our analyses nevertheless reveal surprisingly substantial signatures of introgression among taxa that currently have disjunct distributions within drainages (e.g., separate subcatchments). These "anomalies" imply that current biogeography is not always an accurate predictor of introgression history. Our study suggests that both modern and ancient biogeographic shifts, including recent anthropogenic range fragmentation and tectonically driven riven capture events, have influenced introgression histories in this dynamic freshwater fish radiation.

Signals of differential introgression in the genome of natural hybrids of Caribbean anoles

Hybridization facilitates recombination between divergent genetic lineages and can be shaped by both neutral and selective processes. Upon hybridization, loci with no net fitness effects introgress randomly from parental species into the genomes of hybrid individuals. Conversely, alleles from one parental species at some loci may provide a selective advantage to hybrids, resulting in patterns of introgression that do not conform to random expectations. We investigated genomic patterns of differential introgression in natural hybrids of two species of Caribbean anoles, Anolis pulchellus and A. krugi in Puerto Rico. Hybrids exhibit A. pulchellus phenotypes but possess A. krugi mitochondrial DNA, originated from multiple, independent hybridization events, and appear to have replaced pure A. pulchellus across a large area in western Puerto Rico. Combining genome-wide SNP datasets with bioinformatic methods to identify signals of differential introgression in hybrids, we demonstrate that the genomes of hybrids are dominated by pulchellus-derived alleles and show only 10%-20% A. krugi ancestry. The majority of A. krugi loci in hybrids exhibit a signal of non-random differential introgression and include loci linked to genes involved in development and immune function. Three of these genes (delta like canonical notch ligand 1, jagged1 and notch receptor 1) affect cell differentiation and growth and interact with mitochondrial function. Our results suggest that differential non-random introgression for a subset of loci may be driven by selection favouring the inheritance of compatible mitochondrial and nuclear-encoded genes in hybrids.

The ecological importance of hybridization

Hybridization as an evolutionary process has been studied in depth over the past few decades. Research has focused on its role in shaping reproductive barriers, its adaptive value, and its genomic consequences. In contrast, our knowledge of ecological dimensions of hybridization is still in its infancy, despite hybridization being an inherently ecological interaction. Using examples from various organisms, we show that hybridization can affect and be affected by non-reproductive interactions, including predation, competition, parasitism, mutualism, commensalism, and organism-environment interactions, with significant implications for community structure and ecosystem functioning. However, since these dimensions of hybridization have mostly been revealed from studies designed to decipher other evolutionary processes, we argue that much of the eco-evolutionary importance of hybridization is yet to be discovered.

Integrating morphological, molecular and cytogenetic data for F2 sea turtle hybrids diagnosis revealed balanced chromosomal sets

Hybridization could be considered part of the evolutionary history of many species. The hybridization among sea turtle species on the Brazilian coast is atypical and occurs where nesting areas and reproductive seasons overlap. Integrated analysis of morphology and genetics is still scarce, and there is no evidence of the parental chromosome set distribution in sea turtle interspecific hybrids. In this study, chromosome markers previously established for pure sea turtle species were combined with morphological and molecular analyses aiming to recognize genetic composition and chromosome sets in possible interspecific hybrids initially identified by mixed morphology. The data showed that one hybrid could be an F2 individual among Caretta caretta × Eretmochelys imbricata × Chelonia mydas, and another is resulting from backcross between C. caretta × Lepidochelys olivacea. Native alleles of different parental lineages were reported in the hybrids, and, despite this, it was verified that the hybrid chromosome sets were still balanced. Thus, how sea turtle hybridism can affect genetic features in the long term is a concern, as the implications of the crossing-over in hybrid chromosomal sets and the effects on genetic function are still unpredictable.

PhyloCoalSimulations: A Simulator for Network Multispecies Coalescent Models, Including a New Extension for the Inheritance of Gene Flow

We consider the evolution of phylogenetic gene trees along phylogenetic species networks, according to the network multispecies coalescent process, and introduce a new network coalescent model with correlated inheritance of gene flow. This model generalizes two traditional versions of the network coalescent: with independent or common inheritance. At each reticulation, multiple lineages of a given locus are inherited from parental populations chosen at random, either independently across lineages or with positive correlation according to a Dirichlet process. This process may account for locus-specific probabilities of inheritance, for example. We implemented the simulation of gene trees under these network coalescent models in the Julia package PhyloCoalSimulations, which depends on PhyloNetworks and its powerful network manipulation tools. Input species phylogenies can be read in extended Newick format, either in numbers of generations or in coalescent units. Simulated gene trees can be written in Newick format, and in a way that preserves information about their embedding within the species network. This embedding can be used for downstream purposes, such as to simulate species-specific processes like rate variation across species, or for other scenarios as illustrated in this note. This package should be useful for simulation studies and simulation-based inference methods. The software is available open source with documentation and a tutorial at https://github.com/cecileane/PhyloCoalSimulations.jl.

Limited introgression from non-native commercial strains and signatures of adaptation in the key pollinator Bombus terrestris

Insect pollination is fundamental for natural ecosystems and agricultural crops. The bumblebee species Bombus terrestris has become a popular choice for commercial crop pollination worldwide due to its effectiveness and ease of mass rearing. Bumblebee colonies are mass produced for the pollination of more than 20 crops and imported into over 50 countries including countries outside their native ranges, and the risk of invasion by commercial non-native bumblebees is considered an emerging issue for global conservation and biological diversity. Here, we use genome-wide data from seven wild populations close to and far from farms using commercial colonies, as well as commercial populations, to investigate the implications of utilizing commercial bumblebee subspecies in the UK. We find evidence for generally low levels of introgression between commercial and wild bees, with higher admixture proportions in the bees occurring close to farms. We identify genomic regions putatively involved in local and global adaptation, and genes in locally adaptive regions were found to be enriched for functions related to taste receptor activity, oxidoreductase activity, fatty acid and lipid biosynthetic processes. Despite more than 30 years of bumblebee colony importation into the UK, we observe low impact on the genetic integrity of local B. terrestris populations, but we highlight that even limited introgression might negatively affect locally adapted populations.

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.

Phylogeography of a widespread Australian freshwater fish, western carp gudgeon (Eleotridae: Hypseleotris klunzingeri): Cryptic species, hybrid zones, and strong intra-specific divergences

Despite belonging to the most abundant and widespread genus of freshwater fishes in the region, the carp gudgeons of eastern Australia (genus Hypseleotris) have proved taxonomically and ecologically problematic to science since the 19th century. Several molecular studies and a recent taxonomic revision have now shed light on the complex biology and evolutionary history that underlies this group. These studies have demonstrated that carp gudgeons include a sexual/unisexual complex (five sexual species plus an assortment of hemiclonal lineages), many members of which also co-occur with an independent sexual relative, the western carp gudgeon (H. klunzingeri). Here, we fill yet another knowledge gap for this important group by presenting a detailed molecular phylogeographic assessment of the western carp gudgeon across its entire and extensive geographic range. We use a suite of nuclear genetic markers (SNPs and allozymes) plus a matrilineal genealogy (cytb) to demonstrate that H. klunzingeri s.l. also displays considerable taxonomic and phylogeographic complexity. All molecular datasets concur in recognizing the presence of multiple candidate species, two instances of historic between-species admixture, and the existence of a natural hybrid zone between two of the three candidate species found in the Murray-Darling Basin. We also discuss the major phylogeographic patterns evident within each taxon. Together, these analyses provide a robust molecular, taxonomic, and distributional framework to underpin future morphological and ecological investigations on this prominent member of regional freshwater ecosystems in eastern Australia.

Phylogenetic analysis of endogenous viral elements in the rice genome reveals local chromosomal evolution in Oryza AA-genome species

Introduction

Rice genomes contain endogenous viral elements homologous to rice tungro bacilliform virus (RTBV) from the pararetrovirus family Caulimoviridae. These viral elements, known as endogenous RTBV-like sequences (eRTBVLs), comprise five subfamilies, eRTBVL-A, -B, -C, -D, and -X. Four subfamilies (A, B, C, and X) are present to a limited degree in the genomes of the Asian cultivated rice Oryza sativa (spp. japonica and indica) and the closely related wild species Oryza rufipogon.

Methods

The eRTBVL-D sequences are widely distributed within these and other Oryza AA-genome species. Fifteen eRTBVL-D segments identified in the japonica (Nipponbare) genome occur mostly at orthologous chromosomal positions in other AA-genome species. The eRTBVL-D sequences were inserted into the genomes just before speciation of the AA-genome species.

Results and discussion

Ten eRTBVL-D segments are located at six loci, which were used for our evolutionary analyses during the speciation of the AA-genome species. The degree of genetic differentiation varied among the eRTBVL-D segments. Of the six loci, three showed phylogenetic trees consistent with the standard speciation pattern (SSP) of the AA-genome species (Type A), and the other three represented phylogenies different from the SSP (Type B). The atypical phylogenetic trees for the Type B loci revealed chromosome region-specific evolution among the AA-genome species that is associated with phylogenetic incongruences: complex genome rearrangements between eRTBVL-D segments, an introgression between the distant species, and low genetic diversity of a shared eRTBVL-D segment. Using eRTBVL-D as an indicator, this study revealed the phylogenetic incongruence of local chromosomal regions with different topologies that developed during speciation.

River drainage reorganization and reticulate evolution in the Two-Lined Salamander (Eurycea bislineata) species complex

The origin and eventual loss of biogeographic barriers can create alternating periods of allopatry and secondary contact, facilitating gene flow among distinct metapopulations and generating reticulate evolutionary histories that are not adequately described by a bifurcating evolutionary tree. One such example may exist in the two-lined salamander (Eurycea bislineata) species complex, where discordance among morphological and molecular datasets has created a "vexing taxonomic challenge". Previous phylogeographic analyses of mitochondrial DNA (mtDNA) suggested that the reorganization of Miocene paleodrainages drove vicariance and dispersal, but the inherent limitations of a single-locus dataset precluded the evaluation of subsequent gene flow. Here, we generate triple-enzyme restriction site-associated DNA sequencing (3RAD) data for >100 individuals representing all major mtDNA lineages and use a suite of complementary methods to demonstrate that discordance among earlier datasets is best explained by a reticulate evolutionary history influenced by river drainage reorganization. Systematics of such groups should acknowledge these complex histories and relationships that are not strictly hierarchical.

Diphyllon (Berberidaceae) in southern East Asia

East Asia is the richest region of plant biodiversity in the northern temperate zone, and its radiation provides key insights for understanding rapid speciation, including evolutionary patterns and processes. However, it is challenging to investigate the recent evolutionary radiation among plants because of the lack of genetic divergence, phenotypic convergence, and interspecific gene flow. Epimedium sect. Diphyllon is a rarely studied plant lineage endemic to East Asia, especially highly diversified in its southern part. In this study, we report a robust phylogenomic analysis based on genotyping-by-sequencing data of this lineage. The results revealed a clear biogeographic pattern for Epimedium sect. Diphyllon with recognition into two major clades corresponding to the Sino-Himalayan and Sino-Japanese subkingdoms of East Asian Flora and rapid diversification of the extant species dated to the Pleistocene. Evolutionary radiation of Epimedium sect. Diphyllon is characterized by recent and predominant parallel evolution and atavism between the two subkingdom regions, with extensive reticulating hybridization within each region during the course of diversification in southern East Asia. A parallel-atavism-introgression hypothesis is referred to in explaining the radiation of plant diversity in southern East Asia, which represents a potential model for the rapid diversification of plants under global climate cooling in the late Tertiary. Our study advances our understanding of the evolutionary processes of plant radiation in East Asia as well as in other biodiversity hotspot regions.

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 four 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 two 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) gene-trees 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.