Epistatic selection on a selfish Segregation Distorter supergene - drive, recombination, and genetic load

Supergenes in organismal and social development of insects: ideas and opportunities

Supergenes, or regions of the genome containing two or more linked functional mutations that control complex traits, are emerging as a common genetic basis for many striking phenotypic polymorphisms in insects. Now that we know that supergenes are common, we can seek common features of diverse supergene systems. Here, we lay out a framework of open questions (see graphical abstract) that can be addressed separately in each system and, ultimately, compared across systems to seek general patterns in supergene evolution. Few empirical studies have investigated what causes new supergene haplotypes to initially increase in frequency, but to not eventually fix in a population. Resolving the genotype-phenotype connection and isolating functional genes will provide more insight into the forms of selecting shaping supergene evolution. Ultimately, research on supergenes will help to broaden our understanding of how recombination rate variation influences the evolutionary trajectories of sexually reproducing organisms.

Out with the old, in with the new: Meiotic driving of sex chromosome evolution

Chromosomal regions with meiotic drivers exhibit biased transmission (> 50 %) over their competing homologous chromosomal region. These regions often have two prominent genetic features: suppressed meiotic crossing over and rapidly evolving multicopy gene families. Heteromorphic sex chromosomes (e.g., XY) often share these two genetic features with chromosomal regions exhibiting meiotic drive. Here, we discuss parallels between meiotic drive and sex chromosome evolution, how the divergence of heteromorphic sex chromosomes can be influenced by meiotic drive, experimental approaches to study meiotic drive on sex chromosomes, and meiotic drive in traditional and non-traditional model organisms with high-quality genome assemblies. The newly available diversity of high-quality sex chromosome sequences allows us to revisit conventional models of sex chromosome evolution through the lens of meiotic drive.

Long-Term Maintenance of Complex Chromosomal Inversion Polymorphism in Drosophila mediopunctata

Natural selection is known to favor specific gene combinations, thereby shaping the evolution of recombination rates, often through epistatic interactions. However, the dynamics of these interacting factors within natural populations remain poorly understood. In this study, we investigate the long-term maintenance of a complex polymorphism involving linked, nonoverlapping chromosomal inversions in a natural population of Drosophila mediopunctata. Remarkably, even after 30 years-equivalent to roughly 340 generations-two major features have remained unexpectedly stable: the linkage disequilibrium (LD) between inversions, which deviates significantly from the theoretical prediction of decay, and a consistent seasonal cycle pattern of heterozygous excess and homozygous deficiencies. We explored the roles of recombination suppression, epistatic selection, and overdominance in maintaining this stability, examining their alignment with previously described patterns. Our findings reveal that moderate selection coefficients, such as s = 0.0407, are sufficient to maintain the observed LD for the most common haplotypes, albeit leading to an unstable equilibrium. Simulations further reveal that the introduction of overdominance stabilizes the system, enabling the long-term persistence of this complex inversion polymorphism across various frequency scenarios. The stability of this system appears to hinge on a delicate balance between LD, recombination rates, and selective pressures, with overdominance playing a critical role. Our findings highlight the significance of epistatic interactions and selective pressures in shaping evolutionary pathways in natural populations and offer a compelling example of natural selection acting on a complex inversion polymorphism, providing valuable insights into the evolutionary dynamics governing inversion systems.

X centromeric drive may explain the prevalence of polycystic ovary syndrome and other conditions: Genomic structure of the human X chromosome pericentromeric region is consistent with meiotic drive associated with PCOS and other conditions

X chromosome centromeric drive may explain the prevalence of polycystic ovary syndrome and contribute to oocyte aneuploidy, menopause, and other conditions. The mammalian X chromosome may be vulnerable to meiotic drive because of X inactivation in the female germline. The human X pericentromeric region contains genes potentially involved in meiotic mechanisms, including multiple SPIN1 and ZXDC paralogs. This is consistent with a multigenic drive system comprising differential modification of the active and inactive X chromosome centromeres in female primordial germ cells and preferential segregation of the previously inactivated X chromosome centromere to the polar body at meiosis I. The drive mechanism may explain differences in X chromosome regulation in the female germlines of the human and mouse and, based on the functions encoded by the genes in the region, the transmission of X pericentromeric genetic or epigenetic variants to progeny could contribute to preeclampsia, autism, and differences in sexual differentiation.

Causes and consequences of a complex recombinational landscape in the ant Cardiocondyla obscurior

Eusocial Hymenoptera have the highest recombination rates among all multicellular animals studied so far, but it is unclear why this is and how this affects the biology of individual species. A high-resolution linkage map for the ant Cardiocondyla obscurior corroborates genome-wide high recombination rates reported for ants (8.1 cM/Mb). However, recombination is locally suppressed in regions that are enriched with TEs, that have strong haplotype divergence, or that show signatures of epistatic selection in C. obscurior The results do not support the hypotheses that high recombination rates are linked to phenotypic plasticity or to modulating selection efficiency. Instead, genetic diversity and the frequency of structural variants correlate positively with local recombination rates, potentially compensating for the low levels of genetic variation expected in haplodiploid social Hymenoptera with low effective population size. Ultimately, the data show that recombination contributes to within-population polymorphism and to the divergence of the lineages within C. obscurior.

Divergence and conservation of the meiotic recombination machinery

Sexually reproducing eukaryotes use recombination between homologous chromosomes to promote chromosome segregation during meiosis. Meiotic recombination is almost universally conserved in its broad strokes, but specific molecular details often differ considerably between taxa, and the proteins that constitute the recombination machinery show substantial sequence variability. The extent of this variation is becoming increasingly clear because of recent increases in genomic resources and advances in protein structure prediction. We discuss the tension between functional conservation and rapid evolutionary change with a focus on the proteins that are required for the formation and repair of meiotic DNA double-strand breaks. We highlight phylogenetic relationships on different time scales and propose that this remarkable evolutionary plasticity is a fundamental property of meiotic recombination that shapes our understanding of molecular mechanisms in reproductive biology.

How chromosomal inversions reorient the evolutionary process

Inversions are structural mutations that reverse the sequence of a chromosome segment and reduce the effective rate of recombination in the heterozygous state. They play a major role in adaptation, as well as in other evolutionary processes such as speciation. Although inversions have been studied since the 1920s, they remain difficult to investigate because the reduced recombination conferred by them strengthens the effects of drift and hitchhiking, which in turn can obscure signatures of selection. Nonetheless, numerous inversions have been found to be under selection. Given recent advances in population genetic theory and empirical study, here we review how different mechanisms of selection affect the evolution of inversions. A key difference between inversions and other mutations, such as single nucleotide variants, is that the fitness of an inversion may be affected by a larger number of frequently interacting processes. This considerably complicates the analysis of the causes underlying the evolution of inversions. We discuss the extent to which these mechanisms can be disentangled, and by which approach.

New perspectives on the causes and consequences of male meiotic drive

Gametogenesis is vulnerable to selfish genetic elements that bias their transmission to the next generation by cheating meiosis. These so-called meiotic drivers are widespread in plants, animals, and fungi and can impact genome evolution. Here, we summarize recent progress on the causes and consequences of meiotic drive in males, where selfish elements attack vulnerabilities in spermatogenesis. Advances in genomics provide new insights into the organization and dynamics of driving chromosomes in natural populations. Common themes, including small RNAs, gene duplications, and heterochromatin, emerged from these studies. Interdisciplinary approaches combining evolutionary genomics with molecular and cell biology are beginning to unravel the mysteries of drive and suppression mechanisms. These approaches also provide insights into fundamental processes in spermatogenesis and chromatin regulation.

Meiotic drive of noncentromeric loci in mammalian meiosis II eggs

The germline produces haploid gametes through a specialized cell division called meiosis. In general, homologous chromosomes from each parent segregate randomly to the daughter cells during meiosis, providing parental alleles with an equal chance of transmission. Meiotic drivers are selfish elements who cheat this process to increase their transmission rate. In female meiosis, selfish centromeres and noncentromeric drivers cheat by preferentially segregating to the egg cell. Selfish centromeres cheat in meiosis I (MI), while noncentromeric drivers can cheat in both meiosis I and meiosis II (MII). Here, we highlight recent advances on our understanding of the molecular mechanisms underlying these genetic cheating strategies, especially focusing on mammalian systems, and discuss new models of how noncentromeric selfish drivers can cheat in MII eggs.

Reproducible evaluation of transposable element detectors with McClintock 2 guides accurate inference of Ty insertion patterns in yeast

BACKGROUND

Many computational methods have been developed to detect non-reference transposable element (TE) insertions using short-read whole genome sequencing data. The diversity and complexity of such methods often present challenges to new users seeking to reproducibly install, execute, or evaluate multiple TE insertion detectors.

RESULTS

We previously developed the McClintock meta-pipeline to facilitate the installation, execution, and evaluation of six first-generation short-read TE detectors. Here, we report a completely re-implemented version of McClintock written in Python using Snakemake and Conda that improves its installation, error handling, speed, stability, and extensibility. McClintock 2 now includes 12 short-read TE detectors, auxiliary pre-processing and analysis modules, interactive HTML reports, and a simulation framework to reproducibly evaluate the accuracy of component TE detectors. When applied to the model microbial eukaryote Saccharomyces cerevisiae, we find substantial variation in the ability of McClintock 2 components to identify the precise locations of non-reference TE insertions, with RelocaTE2 showing the highest recall and precision in simulated data. We find that RelocaTE2, TEMP, TEMP2 and TEBreak provide consistent estimates of [Formula: see text]50 non-reference TE insertions per strain and that Ty2 has the highest number of non-reference TE insertions in a species-wide panel of [Formula: see text]1000 yeast genomes. Finally, we show that best-in-class predictors for yeast applied to resequencing data have sufficient resolution to reveal a dyad pattern of integration in nucleosome-bound regions upstream of yeast tRNA genes for Ty1, Ty2, and Ty4, allowing us to extend knowledge about fine-scale target preferences revealed previously for experimentally-induced Ty1 insertions to spontaneous insertions for other copia-superfamily retrotransposons in yeast.

CONCLUSION

McClintock ( https://github.com/bergmanlab/mcclintock/ ) provides a user-friendly pipeline for the identification of TEs in short-read WGS data using multiple TE detectors, which should benefit researchers studying TE insertion variation in a wide range of different organisms. Application of the improved McClintock system to simulated and empirical yeast genome data reveals best-in-class methods and novel biological insights for one of the most widely-studied model eukaryotes and provides a paradigm for evaluating and selecting non-reference TE detectors in other species.

Evidence of epistasis in regions of long-range linkage disequilibrium across five complex diseases in the UK Biobank and eMERGE datasets

Leveraging linkage disequilibrium (LD) patterns as representative of population substructure enables the discovery of additive association signals in genome-wide association studies (GWASs). Standard GWASs are well-powered to interrogate additive models; however, new approaches are required for invesigating other modes of inheritance such as dominance and epistasis. Epistasis, or non-additive interaction between genes, exists across the genome but often goes undetected because of a lack of statistical power. Furthermore, the adoption of LD pruning as customary in standard GWASs excludes detection of sites that are in LD but might underlie the genetic architecture of complex traits. We hypothesize that uncovering long-range interactions between loci with strong LD due to epistatic selection can elucidate genetic mechanisms underlying common diseases. To investigate this hypothesis, we tested for associations between 23 common diseases and 5,625,845 epistatic SNP-SNP pairs (determined by Ohta's D statistics) in long-range LD (>0.25 cM). Across five disease phenotypes, we identified one significant and four near-significant associations that replicated in two large genotype-phenotype datasets (UK Biobank and eMERGE). The genes that were most likely involved in the replicated associations were (1) members of highly conserved gene families with complex roles in multiple pathways, (2) essential genes, and/or (3) genes that were associated in the literature with complex traits that display variable expressivity. These results support the highly pleiotropic and conserved nature of variants in long-range LD under epistatic selection. Our work supports the hypothesis that epistatic interactions regulate diverse clinical mechanisms and might especially be driving factors in conditions with a wide range of phenotypic outcomes.

Reproducible evaluation of transposable element detectors with McClintock 2 guides accurate inference of Ty insertion patterns in yeast

BACKGROUND

Many computational methods have been developed to detect non-reference transposable element (TE) insertions using short-read whole genome sequencing data. The diversity and complexity of such methods often present challenges to new users seeking to reproducibly install, execute, or evaluate multiple TE insertion detectors.

RESULTS

We previously developed the McClintock meta-pipeline to facilitate the installation, execution, and evaluation of six first-generation short-read TE detectors. Here, we report a completely re-implemented version of McClintock written in Python using Snakemake and Conda that improves its installation, error handling, speed, stability, and extensibility. McClintock 2 now includes 12 short-read TE detectors, auxiliary pre-processing and analysis modules, interactive HTML reports, and a simulation framework to reproducibly evaluate the accuracy of component TE detectors. When applied to the model microbial eukaryote Saccharomyces cerevisiae, we find substantial variation in the ability of McClintock 2 components to identify the precise locations of non-reference TE insertions, with RelocaTE2 showing the highest recall and precision in simulated data. We find that RelocaTE2, TEMP, TEMP2 and TEBreak provide a consistent and biologically meaningful view of non-reference TE insertions in a species-wide panel of ∼1000 yeast genomes, as evaluated by coverage-based abundance estimates and expected patterns of tRNA promoter targeting. Finally, we show that best-in-class predictors for yeast have sufficient resolution to reveal a dyad pattern of integration in nucleosome-bound regions upstream of yeast tRNA genes for Ty1, Ty2, and Ty4, allowing us to extend knowledge about fine-scale target preferences first revealed experimentally for Ty1 to natural insertions and related copia-superfamily retrotransposons in yeast.

CONCLUSION

McClintock (https://github.com/bergmanlab/mcclintock/) provides a user-friendly pipeline for the identification of TEs in short-read WGS data using multiple TE detectors, which should benefit researchers studying TE insertion variation in a wide range of different organisms. Application of the improved McClintock system to simulated and empirical yeast genome data reveals best-in-class methods and novel biological insights for one of the most widely-studied model eukaryotes and provides a paradigm for evaluating and selecting non-reference TE detectors for other species.

A method for low-coverage single-gamete sequence analysis demonstrates adherence to Mendel's first law across a large sample of human sperm

Recently published single-cell sequencing data from individual human sperm (n=41,189; 969-3377 cells from each of 25 donors) offer an opportunity to investigate questions of inheritance with improved statistical power, but require new methods tailored to these extremely low-coverage data (∼0.01× per cell). To this end, we developed a method, named rhapsodi, that leverages sparse gamete genotype data to phase the diploid genomes of the donor individuals, impute missing gamete genotypes, and discover meiotic recombination breakpoints, benchmarking its performance across a wide range of study designs. We then applied rhapsodi to the sperm sequencing data to investigate adherence to Mendel's Law of Segregation, which states that the offspring of a diploid, heterozygous parent will inherit either allele with equal probability. While the vast majority of loci adhere to this rule, research in model and non-model organisms has uncovered numerous exceptions whereby 'selfish' alleles are disproportionately transmitted to the next generation. Evidence of such 'transmission distortion' (TD) in humans remains equivocal in part because scans of human pedigrees have been under-powered to detect small effects. After applying rhapsodi to the sperm data and scanning for evidence of TD, our results exhibited close concordance with binomial expectations under balanced transmission. Together, our work demonstrates that rhapsodi can facilitate novel uses of inferred genotype data and meiotic recombination events, while offering a powerful quantitative framework for testing for TD in other cohorts and study systems.