Meiotic drive is associated with sexual incompatibility in Neurospora.

Evolution of Bateson-Dobzhansky-Muller (BDM) incompatibilities is thought to represent a key step in the formation of separate species. They are incompatible alleles that have evolved in separate populations and are exposed in hybrid offspring as hybrid sterility or lethality. In this study, we reveal a previously unconsidered mechanism promoting the formation of BDM incompatibilities, meiotic drive. Theoretical studies have evaluated the role that meiotic drive, the phenomenon whereby selfish elements bias their transmission to progeny at ratios above 50:50, plays in speciation, and have mostly concluded that drive could not result in speciation on its own. Using the model fungus Neurospora, we demonstrate that the large meiotic drive haplotypes, Sk-2 and Sk-3, contain putative sexual incompatibilities. Our experiments revealed that although crosses between Neurospora intermedia and Neurospora metzenbergii produce viable progeny at appreciable rates, when strains of N. intermedia carry Sk-2 or Sk-3 the proportion of viable progeny drops substantially. Additionally, it appears that Sk-2 and Sk-3 have accumulated different incompatibility phenotypes, consistent with their independent evolutionary history. This research illustrates how meiotic drive can contribute to reproductive isolation between populations, and thereby speciation.

The spore killers, fungal meiotic driver elements.

During meiosis, both alleles of any given gene should have equal chances of being inherited by the progeny. There are a number of reasons why, however, this is not the case, with one of the most intriguing instances presenting itself as the phenomenon of meiotic drive. Genes that are capable of driving can manipulate the ratio of alleles among viable meiotic products so that they are inherited in more than half of them. In many cases, this effect is achieved by direct antagonistic interactions, where the driving allele inhibits or otherwise eliminates the alternative allele. In ascomycete fungi, meiotic products are packaged directly into ascospores; thus, the effect of meiotic drive has been given the nefarious moniker, "spore killing." In recent years, many of the known spore killers have been elevated from mysterious phenotypes to well-described systems at genetic, genomic, and molecular levels. In this review, we describe the known diversity of spore killers and synthesize the varied pieces of data from each system into broader trends regarding genome architecture, mechanisms of resistance, the role of transposable elements, their effect on population dynamics, speciation and gene flow, and finally how they may be developed as synthetic drivers. We propose that spore killing is common, but that it is under-observed because of a lack of studies on natural populations. We encourage researchers to seek new spore killers to build on the knowledge that these remarkable genetic elements can teach us about meiotic drive, genomic conflict, and evolution more broadly.

An introgressed gene causes meiotic drive in Neurospora sitophila.

Meiotic drive elements cause their own preferential transmission following meiosis. In fungi, this phenomenon takes the shape of spore killing, and in the filamentous ascomycete Neurospora sitophila, the Sk-1 spore killer element is found in many natural populations. In this study, we identify the gene responsible for spore killing in Sk-1 by generating both long- and short-read genomic data and by using these data to perform a genome-wide association test. We name this gene Spk-1 Through molecular dissection, we show that a single 405-nt-long open reading frame generates a product that both acts as a poison capable of killing sibling spores and as an antidote that rescues spores that produce it. By phylogenetic analysis, we demonstrate that the gene has likely been introgressed from the closely related species Neurospora hispaniola, and we identify three subclades of N. sitophila, one where Sk-1 is fixed, another where Sk-1 is absent, and a third where both killer and sensitive strain are found. Finally, we show that spore killing can be suppressed through an RNA interference-based genome defense pathway known as meiotic silencing by unpaired DNA. Spk-1 is not related to other known meiotic drive genes, and similar sequences are only found within Neurospora These results shed light on the diversity of genes capable of causing meiotic drive, their origin and evolution, and their interaction with the host genome.

Inferring Adaptive Introgression Using Hidden Markov Models.

Adaptive introgression-the flow of adaptive genetic variation between species or populations-has attracted significant interest in recent years and it has been implicated in a number of cases of adaptation, from pesticide resistance and immunity, to local adaptation. Despite this, methods for identification of adaptive introgression from population genomic data are lacking. Here, we present Ancestry_HMM-S, a hidden Markov model-based method for identifying genes undergoing adaptive introgression and quantifying the strength of selection acting on them. Through extensive validation, we show that this method performs well on moderately sized data sets for realistic population and selection parameters. We apply Ancestry_HMM-S to a data set of an admixed Drosophila melanogaster population from South Africa and we identify 17 loci which show signatures of adaptive introgression, four of which have previously been shown to confer resistance to insecticides. Ancestry_HMM-S provides a powerful method for inferring adaptive introgression in data sets that are typically collected when studying admixed populations. This method will enable powerful insights into the genetic consequences of admixture across diverse populations. Ancestry_HMM-S can be downloaded from https://github.com/jesvedberg/Ancestry_HMM-S/.

On the Distribution of Tract Lengths During Adaptive Introgression.

Admixture is increasingly being recognized as an important factor in evolutionary genetics. The distribution of genomic admixture tracts, and the resulting effects on admixture linkage disequilibrium, can be used to date the timing of admixture between species or populations. However, the theory used for such prediction assumes selective neutrality despite the fact that many famous examples of admixture involve natural selection acting for or against admixture. In this paper, we investigate the effects of positive selection on the distribution of tract lengths. We develop a theoretical framework that relies on approximating the trajectory of the selected allele using a logistic function. By numerically calculating the expected allele trajectory, we also show that the approach can be extended to cases where the logistic approximation is poor due to the effects of genetic drift. Using simulations, we show that the model is highly accurate under most scenarios. We use the model to show that positive selection on average will tend to increase the admixture tract length. However, perhaps counter-intuitively, conditional on the allele frequency at the time of sampling, positive selection will actually produce shorter expected tract lengths. We discuss the consequences of our results in interpreting the timing of the introgression of EPAS1 from Denisovans into the ancestors of Tibetans.

Horizontal transmission and recombination maintain forever young bacterial symbiont genomes.

Bacterial symbionts bring a wealth of functions to the associations they participate in, but by doing so, they endanger the genes and genomes underlying these abilities. When bacterial symbionts become obligately associated with their hosts, their genomes are thought to decay towards an organelle-like fate due to decreased homologous recombination and inefficient selection. However, numerous associations exist that counter these expectations, especially in marine environments, possibly due to ongoing horizontal gene flow. Despite extensive theoretical treatment, no empirical study thus far has connected these underlying population genetic processes with long-term evolutionary outcomes. By sampling marine chemosynthetic bacterial-bivalve endosymbioses that range from primarily vertical to strictly horizontal transmission, we tested this canonical theory. We found that transmission mode strongly predicts homologous recombination rates, and that exceedingly low recombination rates are associated with moderate genome degradation in the marine symbionts with nearly strict vertical transmission. Nonetheless, even the most degraded marine endosymbiont genomes are occasionally horizontally transmitted and are much larger than their terrestrial insect symbiont counterparts. Therefore, horizontal transmission and recombination enable efficient natural selection to maintain intermediate symbiont genome sizes and substantial functional genetic variation.

Neurospora from Natural Populations: Population Genomics Insights into the Life History of a Model Microbial Eukaryote.

The ascomycete filamentous fungus Neurospora crassa played a historic role in experimental biology and became a model system for genetic research. Stimulated by a systematic effort to collect wild strains initiated by Stanford geneticist David Perkins, the genus Neurospora has also become a basic model for the study of evolutionary processes, speciation, and population biology. In this chapter, we will first trace the history that brought Neurospora into the era of population genomics. We will then cover the major contributions of population genomic investigations using Neurospora to our understanding of microbial biogeography and speciation, and review recent work using population genomics and genome-wide association mapping that illustrates the unique potential of Neurospora as a model for identifying the genetic basis of (potentially adaptive) phenotypes in filamentous fungi. The advent of population genomics has contributed to firmly establish Neurospora as a complete model system and we hope our review will entice biologists to include Neurospora in their research.

Combinations of Spok genes create multiple meiotic drivers in Podospora.

Meiotic drive is the preferential transmission of a particular allele during sexual reproduction. The phenomenon is observed as spore killing in multiple fungi. In natural populations of Podospora anserina, seven spore killer types (Psks) have been identified through classical genetic analyses. Here we show that the Spok gene family underlies the Psks. The combination of Spok genes at different chromosomal locations defines the spore killer types and creates a killing hierarchy within a population. We identify two novel Spok homologs located within a large (74-167 kbp) region (the Spok block) that resides in different chromosomal locations in different strains. We confirm that the SPOK protein performs both killing and resistance functions and show that these activities are dependent on distinct domains, a predicted nuclease and kinase domain. Genomic and phylogenetic analyses across ascomycetes suggest that the Spok genes disperse through cross-species transfer, and evolve by duplication and diversification within lineages.

Identification of rfk-1, a Meiotic Driver Undergoing RNA Editing in Neurospora.

Sk-2 is a meiotic drive element that was discovered in wild populations of Neurospora fungi over 40 years ago. While early studies quickly determined that Sk-2 transmits itself through sexual reproduction in a biased manner via spore killing, the genetic factors responsible for this phenomenon have remained mostly unknown. Here, we identify and characterize rfk-1, a gene required for Sk-2-based spore killing. The rfk-1 gene contains four exons, three introns, and two stop codons, the first of which undergoes RNA editing to a tryptophan codon during sexual development. Translation of an unedited rfk-1 transcript in vegetative tissue is expected to produce a 102-amino acid protein, whereas translation of an edited rfk-1 transcript in sexual tissue is expected to produce a protein with 130 amino acids. These findings indicate that unedited and edited rfk-1 transcripts exist and that these transcripts could have different roles with respect to the mechanism of meiotic drive by spore killing. Regardless of RNA editing, spore killing only succeeds if rfk-1 transcripts avoid silencing caused by a genome defense process called meiotic silencing by unpaired DNA (MSUD). We show that rfk-1's MSUD avoidance mechanism is linked to the genomic landscape surrounding the rfk-1 gene, which is located near the Sk-2 border on the right arm of chromosome III. In addition to demonstrating that the location of rfk-1 is critical to spore-killing success, our results add to accumulating evidence that MSUD helps protect Neurospora genomes from complex meiotic drive elements.

Convergent evolution of complex genomic rearrangements in two fungal meiotic drive elements.

Meiotic drive is widespread in nature. The conflict it generates is expected to be an important motor for evolutionary change and innovation. In this study, we investigated the genomic consequences of two large multi-gene meiotic drive elements, Sk-2 and Sk-3, found in the filamentous ascomycete Neurospora intermedia. Using long-read sequencing, we generated the first complete and well-annotated genome assemblies of large, highly diverged, non-recombining regions associated with meiotic drive elements. Phylogenetic analysis shows that, even though Sk-2 and Sk-3 are located in the same chromosomal region, they do not form sister clades, suggesting independent origins or at least a long evolutionary separation. We conclude that they have in a convergent manner accumulated similar patterns of tandem inversions and dense repeat clusters, presumably in response to similar needs to create linkage between genes causing drive and resistance.

Large-scale suppression of recombination predates genomic rearrangements in Neurospora tetrasperma.

A common feature of eukaryote genomes is large chromosomal regions where recombination is absent or strongly reduced, but the factors that cause this reduction are not well understood. Genomic rearrangements have often been implicated, but they may also be a consequence of recombination suppression rather than a cause. In this study, we generate eight high-quality genomic data sets of the filamentous ascomycete Neurospora tetrasperma, a fungus that lacks recombination over most of its largest chromosome. The genomes surprisingly reveal collinearity of the non-recombining regions and although large inversions are enriched in these regions, we conclude these inversions to be derived and not the cause of the suppression. To our knowledge, this is the first time that non-recombining, genic regions as large as 86% of a full chromosome (or 8 Mbp), are shown to be collinear. These findings are of significant interest for our understanding of the evolution of sex chromosomes and other supergene complexes.

Structural, biochemical, and computational characterization of the glycoside hydrolase family 7 cellobiohydrolase of the tree-killing fungus Heterobasidion irregulare.

Root rot fungi of the Heterobasidion annosum complex are the most damaging pathogens in temperate forests, and the recently sequenced Heterobasidion irregulare genome revealed over 280 carbohydrate-active enzymes. Here, H. irregulare was grown on biomass, and the most abundant protein in the culture filtrate was identified as the only family 7 glycoside hydrolase in the genome, which consists of a single catalytic domain, lacking a linker and carbohydrate-binding module. The enzyme, HirCel7A, was characterized biochemically to determine the optimal conditions for activity. HirCel7A was crystallized and the structure, refined at 1.7 Å resolution, confirms that HirCel7A is a cellobiohydrolase rather than an endoglucanase, with a cellulose-binding tunnel that is more closed than Phanerochaete chrysosporium Cel7D and more open than Hypocrea jecorina Cel7A, suggesting intermediate enzyme properties. Molecular simulations were conducted to ascertain differences in enzyme-ligand interactions, ligand solvation, and loop flexibility between the family 7 glycoside hydrolase cellobiohydrolases from H. irregulare, H. jecorina, and P. chrysosporium. The structural comparisons and simulations suggest significant differences in enzyme-ligand interactions at the tunnel entrance in the -7 to -4 binding sites and suggest that a tyrosine residue at the tunnel entrance of HirCel7A may serve as an additional ligand-binding site. Additionally, the loops over the active site in H. jecorina Cel7A are more closed than loops in the other two enzymes, which has implications for the degree of processivity, endo-initiation, and substrate dissociation. Overall, this study highlights molecular level features important to understanding this biologically and industrially important family of glycoside hydrolases.