The breakdown of genomic ancestry blocks in hybrid lineages given a finite number of recombination sites

simRestore: A decision-making tool for adaptive management of the native genetic status of wild populations

Anthropogenic hybridization, or higher and non-natural rates of gene flow directly and indirectly induced by human activities, is considered a significant threat to biodiversity. The primary concern for conservation is the potential for genomic extinction and loss of adaptiveness for native species due to the extensive introgression of non-native genes. To alleviate or reverse trends for such scenarios requires the direct integration of genomic data within a model framework for effective management. Towards this end, we developed the simRestore R program as a decision-making tool that integrates ecological and genomic information to simulate ancestry outcomes from optimized conservation strategies. In short, the program optimizes supplementation and removal strategies across generations until a set native genetic threshold is reached within the studied population. Importantly, in addition to helping with initial decision-making, simulations can be updated with the outcomes of ongoing efforts, allowing for the adaptive management of populations. After demonstrating functionality, we apply and optimize among actionable management strategies for the endangered Hawaiian duck for which the current primary threat is genetic extinction through ongoing anthropogenic hybridization with feral mallards. Simulations demonstrate that supplemental and removal efforts can be strategically tailored to move the genetic ancestry of Hawaii's hybrid populations towards Hawaiian duck without the need to completely start over. Further, we discuss ecological parameter sensitivity, including which factors are most important to ensure genetic outcomes (i.e. number of offspring). Finally, to facilitate use, the program is also available online as a Shiny Web application.

IBD sharing patterns as intra-breed admixture indicators in small ruminants

In this study, we investigated how IBD patterns shared between individuals of the same breed could be informative of its admixture level, with the underlying assumption that the most admixed breeds, i.e. the least genetically isolated, should have a much more fragmented genome. We considered 111 goat breeds (i.e. 2501 individuals) and 156 sheep breeds (i.e. 3304 individuals) from Europe, Africa and Asia, for which beadchip SNP genotypes had been performed. We inferred the breed's level of admixture from: (i) the proportion of the genome shared by breed's members (i.e. "genetic integrity level" assessed from ADMIXTURE software analyses), and (ii) the "AV index" (calculated from Reynolds' genetic distances), used as a proxy for the "genetic distinctiveness". In both goat and sheep datasets, the statistical analyses (comparison of means, Spearman correlations, LM and GAM models) revealed that the most genetically isolated breeds, also showed IBD profiles made up of more shared IBD segments, which were also longer. These results pave the way for further research that could lead to the development of admixture indicators, based on the characterization of intra-breed shared IBD segments, particularly effective as they would be independent of the knowledge of the whole genetic landscape in which the breeds evolve. Finally, by highlighting the fragmentation experienced by the genomes subjected to crossbreeding carried out over the last few generations, the study reminds us of the need to preserve local breeds and the integrity of their adaptive architectures that have been shaped over the centuries.

SARS-CoV-2 Recombinants: Genomic Comparison between XBF and Its Parental Lineages

Recombination events are very common and represent one of the primary drivers of RNA virus evolution. The XBF SARS-CoV-2 lineage is one of the most recently generated recombinants during the COVID-19 pandemic. It is a recombinant of BA.5.2.3 and BA.2.75.3, both descendants of lineages that caused many concerns (BA.5 and BA.2.75, respectively). Here, we performed a genomic survey focused on comparing the recombinant XBF with its parental lineages to provide a comprehensive assessment of the evolutionary potential, epidemiological trajectory, and potential risks. Genetic analyses indicated that although XBF initially showed the typical expansion depicted by a steep curve, causing several concerns, currently there is no indication of significant expansion potential or a contagion rate surpassing that of other currently active or previously prevalent lineages. BSP indicated that the peak has been reached around 19 October 2022 and then the genetic variability suffered slight oscillations until early 5 March 2023 when the population size reduced for the last time starting its last plateau that is still lasting. Structural analyses confirmed its reduced potential, also indicating that properties of NTDs and RBDs of XBF and its parental lineages present no significant difference. Of course, cautionary measures must still be taken and genome-based monitoring remains the best tool for detecting any important changes in viral genome composition.

Admixture of evolutionary rates across a butterfly hybrid zone

Hybridization is a major evolutionary force that can erode genetic differentiation between species, whereas reproductive isolation maintains such differentiation. In studying a hybrid zone between the swallowtail butterflies Papilio syfanius and Papilio maackii (Lepidoptera: Papilionidae), we made the unexpected discovery that genomic substitution rates are unequal between the parental species. This phenomenon creates a novel process in hybridization, where genomic regions most affected by gene flow evolve at similar rates between species, while genomic regions with strong reproductive isolation evolve at species-specific rates. Thus, hybridization mixes evolutionary rates in a way similar to its effect on genetic ancestry. Using coalescent theory, we show that the rate-mixing process provides distinct information about levels of gene flow across different parts of genomes, and the degree of rate-mixing can be predicted quantitatively from relative sequence divergence ([Formula: see text]) between the hybridizing species at equilibrium. Overall, we demonstrate that reproductive isolation maintains not only genomic differentiation, but also the rate at which differentiation accumulates. Thus, asymmetric rates of evolution provide an additional signature of loci involved in reproductive isolation.

Hybridization underlies localized trait evolution in cavefish

Introgressive hybridization may play an integral role in local adaptation and speciation (Taylor and Larson, 2019). In the Mexican tetra Astyanax mexicanus, cave populations have repeatedly evolved traits including eye loss, sleep loss, and albinism. Of the 30 caves inhabited by A. mexicanus, Chica cave is unique because it contains multiple pools inhabited by putative hybrids between surface and cave populations (Mitchell et al., 1977), providing an opportunity to investigate the impact of hybridization on complex trait evolution. We show that hybridization between cave and surface populations may contribute to localized variation in traits associated with cave evolution, including pigmentation, eye development, and sleep. We also uncover an example of convergent evolution in a circadian clock gene in multiple cavefish lineages and burrowing mammals, suggesting a shared genetic mechanism underlying circadian disruption in subterranean vertebrates. Our results provide insight into the role of hybridization in facilitating phenotypic evolution.

•

Hybridization leads to a localized difference in sleep duration within a single cave

•

Genomic analysis identifies coding differences in Cry1A across cave pools

•

Changes in Cry1A appear to be conserved in cavefish and burrowing mammals

OUP accepted manuscript

STUDY QUESTION

What information and support should be offered to donors, intended parents and donor-conceived people, in general and in consideration of the availability of direct-to-consumer genetic testing and matching services?

SUMMARY ANSWER

For donors, intended parents and donor-conceived offspring, recommendations are made that cover information needs and informed consent, psychosocial implications and disclosure.

WHAT IS KNOWN ALREADY

Trends indicate that the use of donor-assisted conception is growing and guidance is needed to help these recipients/intended parents, the donors and offspring, navigate the rapidly changing environment in which donor-assisted conception takes place.

STUDY DESIGN, SIZE, DURATION

A working group (WG) collaborated on writing recommendations based, where available, on evidence collected from a literature search and expert opinion. Draft recommendations were published for stakeholder review and adapted where relevant based on the comments received.

PARTICIPANTS/MATERIALS, SETTING, METHODS

Papers retrieved from PUBMED were included from 1 January 2014 up to 31 August 2020, focusing on studies published since direct-to-consumer genetic testing has become more widespread and accessible. The current paper is limited to reproductive donation performed in medically assisted reproduction (MAR) centres (and gamete banks): donation outside the medical context was not considered.

MAIN RESULTS AND THE ROLE OF CHANCE

In total, 32 recommendations were made for information provision and support to donors, 32 for intended parents and 27 for donor-conceived offspring requesting information/support.

LIMITATIONS, REASONS FOR CAUTION

The available evidence in the area of reproductive donation is limited and diverse with regards to the context and types of donation. General conclusions and recommendations are largely based on expert opinion and may need to be adapted in light of future research.

WIDER IMPLICATIONS OF THE FINDINGS

These recommendations provide guidance to MAR centres and gamete banks on good practice in information provision and support but should also be considered by regulatory bodies and policymakers at a national and international level to guide regulatory and legislative efforts towards the protection of donors and donor-conceived offspring.

STUDY FUNDING/COMPETING INTEREST(S)

The development of this good practice paper was funded by European Society of Human Reproduction and Embryology (ESHRE), covering expenses associated with the WG meetings, the literature searches and dissemination. The WG members did not receive any payment. The authors have no conflicts of interest to declare.

DISCLAIMER

This document represents the views of ESHRE, which are the result of consensus between the relevant ESHRE stakeholders and where relevant based on the scientific evidence available at the time of preparation.

The recommendations should be used for informational and educational purposes. They should not be interpreted as setting a standard of care, or be deemed inclusive of all proper methods of care nor exclusive of other methods of care reasonably directed to obtaining the same results. They do not replace the need for application of clinical judgement to each individual presentation, nor variations based on locality and facility type.

^†ESHRE pages content is not externally peer reviewed. The manuscript has been approved by the Executive Committee of ESHRE.

A dense linkage map for a large repetitive genome: discovery of the sex-determining region in hybridizing fire-bellied toads (Bombina bombina and Bombina variegata)

Genomic analysis of hybrid zones offers unique insights into emerging reproductive isolation and the dynamics of introgression. Because hybrid genomes consist of blocks inherited from one or the other parental taxon, linkage information is essential. In most cases, the spectrum of local ancestry tracts can be efficiently uncovered from dense linkage maps. Here, we report the development of such a map for the hybridizing toads, Bombina bombina and Bombina variegata (Anura: Bombinatoridae). Faced with the challenge of a large (7–10 Gb), repetitive genome, we set out to identify a large number of Mendelian markers in the nonrepetitive portion of the genome that report B. bombina vs B. variegata ancestry with appropriately quantified statistical support. Bait sequences for targeted enrichment were selected from a draft genome assembly, after filtering highly repetitive sequences. We developed a novel approach to infer the most likely diplotype per sample and locus from the raw read mapping data, which is robust to over-merging and obviates arbitrary filtering thresholds. Validation of the resulting map with 4755 markers underscored the large-scale synteny between Bombina and Xenopus tropicalis. By assessing the sex of late-stage F2 tadpoles from histological sections, we identified the sex-determining region in the Bombina genome to 7 cM on LG5, which is homologous to X. tropicalis chromosome 5, and inferred male heterogamety. Interestingly, chromosome 5 has been repeatedly recruited as a sex chromosome in anurans with XY sex determination.

Estimating the time since admixture from phased and unphased molecular data

After admixture, recombination breaks down genomic blocks of contiguous ancestry. The breakdown of these blocks forms a new “molecular clock” that ticks at a much faster rate than the mutation clock, enabling accurate dating of admixture events in the recent past. However, existing theory on the breakdown of these blocks, or the accumulation of delineations between blocks, so‐called “junctions”, has mostly been limited to using regularly spaced markers on phased data. Here, we present an extension to the theory of junctions using the ancestral recombination graph that describes the expected number of junctions for any distribution of markers along the genome. Furthermore, we provide a new framework to infer the time since admixture using unphased data. We demonstrate both the phased and unphased methods on simulated data and show that our new extensions have improved accuracy with respect to previous methods, especially for smaller population sizes and more ancient admixture times. Lastly, we demonstrate the applicability of our method on three empirical data sets, including labcrosses of yeast (Saccharomyces cerevisae) and two case studies of hybridization in swordtail fish and Populus trees.

Demographic history shapes genomic ancestry in hybrid zones

Demographic factors such as migration rate and population size can impede or facilitate speciation. In hybrid zones, reproductive boundaries between species are tested and demography mediates the opportunity for admixture between lineages that are partially isolated. Genomic ancestry is a powerful tool for revealing the history of admixed populations, but models and methods based on local ancestry are rarely applied to structured hybrid zones. To understand the effects of demography on ancestry in hybrids zones, we performed individual-based simulations under a stepping-stone model, treating migration rate, deme size, and hybrid zone age as parameters. We find that the number of ancestry junctions (the transition points between genomic regions with different ancestries) and heterogenicity (the genomic proportion heterozygous for ancestry) are often closely connected to demographic history. Reducing deme size reduces junction number and heterogenicity. Elevating migration rate increases heterogenicity, but migration affects junction number in more complex ways. We highlight the junction frequency spectrum as a novel and informative summary of ancestry that responds to demographic history. A substantial proportion of junctions are expected to fix when migration is limited or deme size is small, changing the shape of the spectrum. Our findings suggest that genomic patterns of ancestry could be used to infer demographic history in hybrid zones.

Ancestral Admixture Is the Main Determinant of Global Biodiversity in Fission Yeast

Mutation and recombination are key evolutionary processes governing phenotypic variation and reproductive isolation. We here demonstrate that biodiversity within all globally known strains of Schizosaccharomyces pombe arose through admixture between two divergent ancestral lineages. Initial hybridization was inferred to have occurred ∼20-60 sexual outcrossing generations ago consistent with recent, human-induced migration at the onset of intensified transcontinental trade. Species-wide heritable phenotypic variation was explained near-exclusively by strain-specific arrangements of alternating ancestry components with evidence for transgressive segregation. Reproductive compatibility between strains was likewise predicted by the degree of shared ancestry. To assess the genetic determinants of ancestry block distribution across the genome, we characterized the type, frequency, and position of structural genomic variation using nanopore and single-molecule real-time sequencing. Despite being associated with double-strand break initiation points, over 800 segregating structural variants exerted overall little influence on the introgression landscape or on reproductive compatibility between strains. In contrast, we found strong ancestry disequilibrium consistent with negative epistatic selection shaping genomic ancestry combinations during the course of hybridization. This study provides a detailed, experimentally tractable example that genomes of natural populations are mosaics reflecting different evolutionary histories. Exploiting genome-wide heterogeneity in the history of ancestral recombination and lineage-specific mutations sheds new light on the population history of S. pombe and highlights the importance of hybridization as a creative force in generating biodiversity.

Identifying hybrids & the genomics of hybridization: Mallards & American black ducks of Eastern North America

Resolving evolutionary relationships and establishing population structure depends on molecular diagnosability that is often limited for closely related taxa. Here, we use 3,200 ddRAD-seq loci across 290 mallards, American black ducks, and putative hybrids to establish population structure and estimate hybridization rates. We test between traditional assignment probability and accumulated recombination events based analyses to assign hybrids to generational classes. For hybrid identification, we report the distribution of recombination events complements ADMIXTURE simulation by extending resolution past F4 hybrid status; however, caution against hybrid assignment based on accumulated recombination events due to an inability to resolve F1 hybrids. Nevertheless, both analyses suggest that there are relatively few backcrossed stages before a lineage's hybrid ancestry is lost and the offspring are effectively parental again. We conclude that despite high rates of observed interspecific hybridization between mallards and black ducks in the middle part of the 20th century, our results do not support the predicted hybrid swarm. Conversely, we report that mallard samples genetically assigned to western and non-western clusters. We indicate that these non-western mallards likely originated from game-farm stock, suggesting landscape level gene flow between domestic and wild conspecifics.

Hybridization is a recurrent evolutionary stimulus in wild yeast speciation

Hybridization can result in reproductively isolated and phenotypically distinct lineages that evolve as independent hybrid species. How frequently hybridization leads to speciation remains largely unknown. Here we examine the potential recurrence of hybrid speciation in the wild yeast Saccharomyces paradoxus in North America, which comprises two endemic lineages SpB and SpC, and an incipient hybrid species, SpC*. Using whole-genome sequences from more than 300 strains, we uncover the hybrid origin of another group, SpD, that emerged from hybridization between SpC* and one of its parental species, the widespread SpB. We show that SpD has the potential to evolve as a novel hybrid species, because it displays phenotypic novelties that include an intermediate transcriptome profile, and partial reproductive isolation with its most abundant sympatric parental species, SpB. Our findings show that repetitive cycles of divergence and hybridization quickly generate diversity and reproductive isolation, providing the raw material for speciation by hybridization.

Hybridization can contribute to diversity from the genomic to the species level. Here, Eberlein, Hénault et al. investigate genomic, transcriptomic and phenotypic variation among wild lineages of the yeast Saccharomyces paradoxus and suggest that an incipient species has formed by recurrent hybridization.

Genomic consequences of a recent three-way admixture in supplemented wild brown trout populations revealed by local ancestry tracts

Understanding the evolutionary consequences of human-mediated introductions of domesticated strains into the wild and their subsequent admixture with natural populations is of major concern in conservation biology. However, the genomic impacts of stocking from distinct sources (locally derived vs. divergent) on the genetic integrity of wild populations remain poorly understood. We designed an approach based on estimating local ancestry along individual chromosomes to provide a detailed picture of genomic admixture in supplemented populations. We used this approach to document admixture consequences in the brown trout Salmo trutta, for which decades of stocking practices have profoundly impacted the genetic make-up of wild populations. In southern France, small local Mediterranean populations have been subject to successive introductions of domestic strains derived from the Atlantic and Mediterranean lineages. To address the impact of stocking, we evaluate the extent of admixture from both domestic strains within populations, using 75,684 mapped SNPs obtained from double-digested restriction site-associated DNA sequencing. Then, the chromosomal ancestry profiles of admixed individuals reveal a wider diversity of hybrid and introgressed genotypes than estimated using classical methods for inferring ancestry and hybrid pedigrees. In addition, the length distribution of introgressed tracts retained different timings of introgression between the two domestic strains. We finally reveal opposite consequences of admixture on the level of polymorphism of the recipient populations between domestic strains. Our study illustrates the potential of using the information contained in the genomic mosaic of ancestry tracts in combination with classical methods based on allele frequencies for analysing multiple-way admixture with population genomic data.

Signatures of hybridization and speciation in genomic patterns of ancestry

Genomes sampled from hybrid zones between nascent species provide important clues into the speciation process. With advances in genome sequencing and single nucleotide polymorphism (SNP) genotyping, it is now feasible to measure variation in gene flow with high genomic resolution. This progress motivates the development of conceptual and analytical frameworks for hybrid zones that complement well-established cline approaches. We extend the perspective that genomic distributions of ancestry are sensitive indicators of hybridization history. We use simulations to examine the behavior of the number of ancestry junctions-a simple summary of genomic patterns-in hybrid zones under increasingly realistic scenarios. Neutral simulations revealed that ancestry junction number is shaped by population structure, migration rate, and population size. Modeling multiple genetic architectures of hybrid dysfunction, with an emphasis on epistatic hybrid incompatibilities, showed that selection reduces junction number near loci that confer reproductive barriers. The magnitude of this signature was affected by the form of selection, dominance, and genomic location (autosome vs. sex chromosome) of incompatible loci. Our results suggest that researchers can identify loci involved in reproductive isolation by scanning hybrid genomes for local reductions in junction number. We outline necessary directions for future theory and method development to realize this goal.