A lethal mitonuclear incompatibility in complex I of natural hybrids.

The evolution of reproductive barriers is the first step in the formation of new species and can help us understand the diversification of life on Earth. These reproductive barriers often take the form of hybrid incompatibilities, in which alleles derived from two different species no longer interact properly in hybrids1-3. Theory predicts that hybrid incompatibilities may be more likely to arise at rapidly evolving genes4-6 and that incompatibilities involving multiple genes should be common7,8, but there has been sparse empirical data to evaluate these predictions. Here we describe a mitonuclear incompatibility involving three genes whose protein products are in physical contact within respiratory complex I of naturally hybridizing swordtail fish species. Individuals homozygous for mismatched protein combinations do not complete embryonic development or die as juveniles, whereas those heterozygous for the incompatibility have reduced complex I function and unbalanced representation of parental alleles in the mitochondrial proteome. We find that the effects of different genetic interactions on survival are non-additive, highlighting subtle complexity in the genetic architecture of hybrid incompatibilities. Finally, we document the evolutionary history of the genes involved, showing signals of accelerated evolution and evidence that an incompatibility has been transferred between species via hybridization.

Analysis of ancestry heterozygosity suggests that hybrid incompatibilities in threespine stickleback are environment dependent.

Hybrid incompatibilities occur when interactions between opposite ancestry alleles at different loci reduce the fitness of hybrids. Most work on incompatibilities has focused on those that are "intrinsic," meaning they affect viability and sterility in the laboratory. Theory predicts that ecological selection can also underlie hybrid incompatibilities, but tests of this hypothesis using sequence data are scarce. In this article, we compiled genetic data for F2 hybrid crosses between divergent populations of threespine stickleback fish (Gasterosteus aculeatus L.) that were born and raised in either the field (seminatural experimental ponds) or the laboratory (aquaria). Because selection against incompatibilities results in elevated ancestry heterozygosity, we tested the prediction that ancestry heterozygosity will be higher in pond-raised fish compared to those raised in aquaria. We found that ancestry heterozygosity was elevated by approximately 3% in crosses raised in ponds compared to those raised in aquaria. Additional analyses support a phenotypic basis for incompatibility and suggest that environment-specific single-locus heterozygote advantage is not the cause of selection on ancestry heterozygosity. Our study provides evidence that, in stickleback, a coarse-albeit indirect-signal of environment-dependent hybrid incompatibility is reliably detectable and suggests that extrinsic incompatibilities can evolve before intrinsic incompatibilities.

Imbalanced segregation of recombinant haplotypes in hybrid populations reveals inter- and intrachromosomal Dobzhansky-Muller incompatibilities.

Dobzhansky-Muller incompatibilities (DMIs) are a major component of reproductive isolation between species. DMIs imply negative epistasis and are exposed when two diverged populations hybridize. Mapping the locations of DMIs has largely relied on classical genetic mapping. Approaches to date are hampered by low power and the challenge of identifying DMI loci on the same chromosome, because strong initial linkage of parental haplotypes weakens statistical tests. Here, we propose new statistics to infer negative epistasis from haplotype frequencies in hybrid populations. When two divergent populations hybridize, the variance in heterozygosity at two loci decreases faster with time at DMI loci than at random pairs of loci. When two populations hybridize at near-even admixture proportions, the deviation of the observed variance from its expectation becomes negative for the DMI pair. This negative deviation enables us to detect intermediate to strong negative epistasis both within and between chromosomes. In practice, the detection window in hybrid populations depends on the demographic scenario, the recombination rate, and the strength of epistasis. When the initial proportion of the two parental populations is uneven, only strong DMIs can be detected with our method unless migration prevents parental haplotypes from being lost. We use the new statistics to infer candidate DMIs from three hybrid populations of swordtail fish. We identify numerous new DMI candidates, some of which are inferred to interact with several loci within and between chromosomes. Moreover, we discuss our results in the context of an expected enrichment in intrachromosomal over interchromosomal DMIs.

PRDM9 losses in vertebrates are coupled to those of paralogs ZCWPW1 and ZCWPW2.

In most mammals and likely throughout vertebrates, the gene PRDM9 specifies the locations of meiotic double strand breaks; in mice and humans at least, it also aids in their repair. For both roles, many of the molecular partners remain unknown. Here, we take a phylogenetic approach to identify genes that may be interacting with PRDM9 by leveraging the fact that PRDM9 arose before the origin of vertebrates but was lost many times, either partially or entirely-and with it, its role in recombination. As a first step, we characterize PRDM9 domain composition across 446 vertebrate species, inferring at least 13 independent losses. We then use the interdigitation of PRDM9 orthologs across vertebrates to test whether it coevolved with any of 241 candidate genes coexpressed with PRDM9 in mice or associated with recombination phenotypes in mammals. Accounting for the phylogenetic relationship among a subsample of 189 species, we find two genes whose presence and absence is unexpectedly coincident with that of PRDM9: ZCWPW1, which was recently shown to facilitate double strand break repair, and its paralog ZCWPW2, as well as, more tentatively, TEX15 and FBXO47ZCWPW2 is expected to be recruited to sites of PRDM9 binding; its tight coevolution with PRDM9 across vertebrates suggests that it is a key interactor within mammals and beyond, with a role either in recruiting the recombination machinery or in double strand break repair.

Predictability and parallelism in the contemporary evolution of hybrid genomes.

Hybridization between species is widespread across the tree of life. As a result, many species, including our own, harbor regions of their genome derived from hybridization. Despite the recognition that this process is widespread, we understand little about how the genome stabilizes following hybridization, and whether the mechanisms driving this stabilization tend to be shared across species. Here, we dissect the drivers of variation in local ancestry across the genome in replicated hybridization events between two species pairs of swordtail fish: Xiphophorus birchmanni × X. cortezi and X. birchmanni × X. malinche. We find unexpectedly high levels of repeatability in local ancestry across the two types of hybrid populations. This repeatability is attributable in part to the fact that the recombination landscape and locations of functionally important elements play a major role in driving variation in local ancestry in both types of hybrid populations. Beyond these broad scale patterns, we identify dozens of regions of the genome where minor parent ancestry is unusually low or high across species pairs. Analysis of these regions points to shared sites under selection across species pairs, and in some cases, shared mechanisms of selection. We show that one such region is a previously unknown hybrid incompatibility that is shared across X. birchmanni × X. cortezi and X. birchmanni × X. malinche hybrid populations.

Genomic insights into variation in thermotolerance between hybridizing swordtail fishes.

Understanding how organisms adapt to changing environments is a core focus of research in evolutionary biology. One common mechanism is adaptive introgression, which has received increasing attention as a potential route to rapid adaptation in populations struggling in the face of ecological change, particularly global climate change. However, hybridization can also result in deleterious genetic interactions that may limit the benefits of adaptive introgression. Here, we used a combination of genome-wide quantitative trait locus mapping and differential gene expression analyses between the swordtail fish species Xiphophorus malinche and X. birchmanni to study the consequences of hybridization on thermotolerance. While these two species are adapted to different thermal environments, we document a complicated architecture of thermotolerance in hybrids. We identify a region of the genome that contributes to reduced thermotolerance in individuals heterozygous for X. malinche and X. birchmanni ancestry, as well as widespread misexpression in hybrids of genes that respond to thermal stress in the parental species, particularly in the circadian clock pathway. We also show that a previously mapped hybrid incompatibility between X. malinche and X. birchmanni contributes to reduced thermotolerance in hybrids. Together, our results highlight the challenges of understanding the impact of hybridization on complex ecological traits and its potential impact on adaptive introgression.

Assortative mating and persistent reproductive isolation in hybrids.

The emergence of new species is driven by the establishment of mechanisms that limit gene flow between populations. A major challenge is reconciling the theoretical and empirical importance of assortative mating in speciation with the ease with which it can fail. Swordtail fish have an evolutionary history of hybridization and fragile prezygotic isolating mechanisms. Hybridization between two swordtail species likely arose via pollution-mediated breakdown of assortative mating in the 1990s. Here we track unusual genetic patterns in one hybrid population over the past decade using whole-genome sequencing. Hybrids in this population formed separate genetic clusters by 2003, and maintained near-perfect isolation over 25 generations through strong ancestry-assortative mating. However, we also find that assortative mating was plastic, varying in strength over time and disappearing under manipulated conditions. In addition, a nearby population did not show evidence of assortative mating. Thus, our findings suggest that assortative mating may constitute an intermittent and unpredictable barrier to gene flow, but that variation in its strength can have a major effect on how hybrid populations evolve. Understanding how reproductive isolation varies across populations and through time is critical to understanding speciation and hybridization, as well as their dependence on disturbance.

Reproductive isolation of hybrid populations driven by genetic incompatibilities.

Despite its role in homogenizing populations, hybridization has also been proposed as a means to generate new species. The conceptual basis for this idea is that hybridization can result in novel phenotypes through recombination between the parental genomes, allowing a hybrid population to occupy ecological niches unavailable to parental species. Here we present an alternative model of the evolution of reproductive isolation in hybrid populations that occurs as a simple consequence of selection against genetic incompatibilities. Unlike previous models of hybrid speciation, our model does not incorporate inbreeding, or assume that hybrids have an ecological or reproductive fitness advantage relative to parental populations. We show that reproductive isolation between hybrids and parental species can evolve frequently and rapidly under this model, even in the presence of substantial ongoing immigration from parental species and strong selection against hybrids. An interesting prediction of our model is that replicate hybrid populations formed from the same pair of parental species can evolve reproductive isolation from each other. This non-adaptive process can therefore generate patterns of species diversity and relatedness that resemble an adaptive radiation. Intriguingly, several known hybrid species exhibit patterns of reproductive isolation consistent with the predictions of our model.

Early social learning triggers neurogenomic expression changes in a swordtail fish.

Mate choice can play a pivotal role in the nature and extent of reproductive isolation between species. Mating preferences are often dependent on an individual's social experience with adult phenotypes throughout development. We show that olfactory preference in a swordtail fish (Xiphophorus malinche) is affected by previous experience with adult olfactory signals. We compare transcriptome-wide gene expression levels of pooled sensory and brain tissues between three treatment groups that differ by social experience: females with no adult exposure, females exposed to conspecifics and females exposed to heterospecifics. We identify potential functionally relevant genes and biological pathways differentially expressed not only between control and exposure groups, but also between groups exposed to conspecifics and heterospecifics. Based on our results, we speculate that vomeronasal receptor type 2 paralogs may detect species-specific pheromone components and thus play an important role in reproductive isolation between species.

Admix'em: a flexible framework for forward-time simulations of hybrid populations with selection and mate choice.

UNLABELLED: We introduce a new forward-time simulator, Admix'em, that allows for rapid and realistic simulations of admixed populations with selection. Complex selection can be achieved through user-defined fitness and mating-preference probability functions. Users can specify realistic genomic landscapes and model neutral SNPs in addition to sites under selection. Admix'em is designed to simulate selection in admixed populations but can also be used as a general population simulator. Usage and examples are in the supplement. AVAILABILITY AND IMPLEMENTATION: C ++ and OpenMP, supports 64-bit Linux/Unix-like platforms. https://github.com/melop/admixem CONTACT: rcui@age.mpg.de SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Determining epistatic selection in admixed populations.

When two diverging species begin hybridizing, selection against hybridization is likely driven not by single substitutions, but by interactions between incompatible mutations. To identify these incompatibilities in natural populations, researchers examine the extent of nonrandom associations between ancestry at physically unlinked loci in admixed populations. In this approach, which we call 'AD scans', locus pairs with significantly positive 'ancestry disequilibrium' (AD, i.e. locus pairs that positively covary by ancestry) represent incompatible alleles. Past research has uniformly revealed an excess of locus pairs with significantly positive AD, suggesting that dozens to hundreds of incompatibilities separate species. With forward simulations, we show that many realistic demographic scenarios, including recent and/or ongoing hybridization, generate a bias towards positive ancestry disequilibrium. We suggest steps that researchers can take to avoid pitfalls in interpreting AD scans, and present a novel measure of AD, which minimizes but does not fully eliminate bias in the AD distribution. We also show, by simulation, that the tail of the AD distribution is enriched for true incompatibilities. To illustrate the potential power and appropriate caution in interpretation of AD scans, we reanalyse previously published data from two admixed populations of Xiphophorus fishes. Our results imply that the prevalence of positive LD in admixed populations does not in itself support the idea that two-locus incompatibilities are widespread, but the co-enrichment of top AD hits across the two Xiphophorus populations supports the idea that AD scans can identify candidate interspecific incompatibilities.

Ancient hybridization and genomic stabilization in a swordtail fish.

A rapidly increasing body of work is revealing that the genomes of distinct species often exhibit hybrid ancestry, presumably due to postspeciation hybridization between closely related species. Despite the growing number of documented cases, we still know relatively little about how genomes evolve and stabilize following hybridization, and to what extent hybridization is functionally relevant. Here, we examine the case of Xiphophorus nezahualcoyotl, a teleost fish whose genome exhibits significant hybrid ancestry. We show that hybridization was relatively ancient and is unlikely to be ongoing. Strikingly, the genome of X. nezahualcoyotl has largely stabilized following hybridization, distinguishing it from examples such as human-Neanderthal hybridization. Hybridization-derived regions are remarkably distinct from other regions of the genome, tending to be enriched in genomic regions with reduced constraint. These results suggest that selection has played a role in removing hybrid ancestry from certain functionally important regions. Combined with findings in other systems, our results raise many questions about the process of genomic stabilization and the role of selection in shaping patterns of hybrid ancestry in the genome.

Phylogenomics analyses of all species of Swordtails (Genus Xiphophorus ) highlights hybridization precedes speciation.

Hybridization has been recognized as an important driving force for evolution, however studies of the genetic consequence and its cause are still lagging behind in vertebrates due to the lack of appropriate experimental systems. Fish of the central American genus Xiphophorus were proposed to have evolved with multiple ancient and ongoing hybridization events, and served as a valuable research model in evolutionary biology and in biomedical research on human disease for more than a century. Here, we provide the complete genome resource and its annotation of all 26 Xiphophorus species. On this dataset we resolved the so far conflicting phylogeny. Through comparative genomic analyses we investigated the molecular evolution of genes related to melanoma, for a main sexually selected trait and for the genetic control of puberty timing, which are predicted to be involved in pre-and postzygotic isolation and thus to influence the probability of interspecific hybridization in Xiphophorus . We demonstrate dramatic size-variation of some gene families across species, despite the reticulate evolution and short divergence time. Finally, we clarify the hybridization history in the genus Xiphophorus genus, settle the long dispute on the hybridization origin of two Southern swordtails, highlight hybridizations precedes speciation, and reveal the distribution of hybridization ancestry remaining in the fused genome.

Patterns of recombination in snakes reveal a tug-of-war between PRDM9 and promoter-like features.

In some mammals, notably humans, recombination occurs almost exclusively where the protein PRDM9 binds, whereas in vertebrates lacking an intact PRDM9, such as birds and canids, recombination rates are elevated near promoter-like features. To determine whether PRDM9 directs recombination in nonmammalian vertebrates, we focused on an exemplar species with a single, intact PRDM9 ortholog, the corn snake (Pantherophis guttatus). Analyzing historical recombination rates along the genome and crossovers in pedigrees, we found evidence that PRDM9 specifies the location of recombination events, but we also detected a separable effect of promoter-like features. These findings reveal that the uses of PRDM9 and promoter-like features need not be mutually exclusive and instead reflect a tug-of-war that is more even in some species than others.

Pervasive gene flow despite strong and varied reproductive barriers in swordtails.

One of the mechanisms that can lead to the formation of new species occurs through the evolution of reproductive barriers. However, recent research has demonstrated that hybridization has been pervasive across the tree of life even in the presence of strong barriers. Swordtail fishes (genus Xiphophorus) are an emerging model system for studying the interface between these barriers and hybridization. We document overlapping mechanisms that act as barriers between closely related species, X. birchmanni and X. cortezi, by combining genomic sequencing from natural hybrid populations, artificial crosses, behavioral assays, sperm performance, and developmental studies. We show that strong assortative mating plays a key role in maintaining subpopulations with distinct ancestry in natural hybrid populations. Lab experiments demonstrate that artificial F1 crosses experience dysfunction: crosses with X. birchmanni females were largely inviable and crosses with X. cortezi females had a heavily skewed sex ratio. Using F2 hybrids we identify several genomic regions that strongly impact hybrid viability. Strikingly, two of these regions underlie genetic incompatibilities in hybrids between X. birchmanni and its sister species X. malinche. Our results demonstrate that ancient hybridization has played a role in the origin of this shared genetic incompatibility. Moreover, ancestry mismatch at these incompatible regions has remarkably similar consequences for phenotypes and hybrid survival in X. cortezi × X. birchmanni hybrids as in X. malinche × X. birchmanni hybrids. Our findings identify varied reproductive barriers that shape genetic exchange between naturally hybridizing species and highlight the complex evolutionary outcomes of hybridization.

Individual variation in preference behavior in sailfin fish refines the neurotranscriptomic pathway for mate preference.

Social interactions can drive distinct gene expression profiles which may vary by social context. Here we use female sailfin molly fish (Poecilia latipinna) to identify genomic profiles associated with preference behavior in distinct social contexts: male interactions (mate choice) versus female interactions (shoaling partner preference). We measured the behavior of 15 females interacting in a non-contact environment with either two males or two females for 30 min followed by whole-brain transcriptomic profiling by RNA sequencing. We profiled females that exhibited high levels of social affiliation and great variation in preference behavior to identify an order of magnitude more differentially expressed genes associated with behavioral variation than by differences in social context. Using a linear model (limma), we took advantage of the individual variation in preference behavior to identify unique gene sets that exhibited distinct correlational patterns of expression with preference behavior in each social context. By combining limma and weighted gene co-expression network analyses (WGCNA) approaches we identified a refined set of 401 genes robustly associated with mate preference that is independent of shoaling partner preference or general social affiliation. While our refined gene set confirmed neural plasticity pathways involvement in moderating female preference behavior, we also identified a significant proportion of discovered that our preference-associated genes were enriched for 'immune system' gene ontology categories. We hypothesize that the association between mate preference and transcriptomic immune function is driven by the less well-known role of these genes in neural plasticity which is likely involved in higher-order learning and processing during mate choice decisions.

Genomes of two Extinct-in-the-Wild reptiles from Christmas Island reveal distinct evolutionary histories and conservation insights.

Genomics can play important roles in biodiversity conservation, especially for Extinct-in-the-Wild species where genetic factors greatly influence risk of total extinction and probability of successful reintroductions. The Christmas Island blue-tailed skink (Cryptoblepharus egeriae) and Lister's gecko (Lepidodactylus listeri) are two endemic reptile species that went extinct in the wild shortly after the introduction of a predatory snake. After a decade of management, captive populations have expanded from 66 skinks and 43 geckos to several thousand individuals; however, little is known about patterns of genetic variation in these species. Here, we use PacBio HiFi long-read and Hi-C sequencing to generate highly contiguous reference genomes for both reptiles, including the XY chromosome pair in the skink. We then analyse patterns of genetic diversity to infer ancient demography and more recent histories of inbreeding. We observe high genome-wide heterozygosity in the skink (0.007 heterozygous sites per base-pair) and gecko (0.005), consistent with large historical population sizes. However, nearly 10% of the blue-tailed skink reference genome falls within long (>1 Mb) runs of homozygosity (ROH), resulting in homozygosity at all major histocompatibility complex (MHC) loci. In contrast, we detect a single ROH in Lister's gecko. We infer from the ROH lengths that related skinks may have established the captive populations. Despite a shared recent extinction in the wild, our results suggest important differences in these species' histories and implications for management. We show how reference genomes can contribute evolutionary and conservation insights, and we provide resources for future population-level and comparative genomic studies in reptiles.

Complex hybridization between deeply diverged fish species in a disturbed ecosystem.

Over the past two decades researchers have documented the extent of natural hybridization between closely related species using genomic tools. Many species across the tree of life show evidence of past hybridization with their evolutionary relatives. In some cases, this hybridization is complex-involving gene flow between more than two species. While hybridization is common over evolutionary timescales, some researchers have proposed that it may be even more common in contemporary populations where anthropogenic disturbance has modified a myriad of aspects of the environments in which organisms live and reproduce. Here, we develop a flexible tool for local ancestry inference in hybrids derived from three source populations and describe a complex, recent hybridization event between distantly related swordtail fish lineages (Xiphophorus) and its potential links to anthropogenic disturbance.