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

Evolution and genetic architecture of sex-limited polymorphism in cuckoos

Sex-limited polymorphism has evolved in many species including our own. Yet, we lack a detailed understanding of the underlying genetic variation and evolutionary processes at work. The brood parasitic common cuckoo (Cuculus canorus) is a prime example of female-limited color polymorphism, where adult males are monochromatic gray and females exhibit either gray or rufous plumage. This polymorphism has been hypothesized to be governed by negative frequency-dependent selection whereby the rarer female morph is protected against harassment by males or from mobbing by parasitized host species. Here, we show that female plumage dichromatism maps to the female-restricted genome. We further demonstrate that, consistent with balancing selection, ancestry of the rufous phenotype is shared with the likewise female dichromatic sister species, the oriental cuckoo (Cuculus optatus). This study shows that sex-specific polymorphism in trait variation can be resolved by genetic variation residing on a sex-limited chromosome and be maintained across species boundaries.

Transposon Removal Reveals Their Adaptive Fitness Contribution

Transposable elements are molecular parasites that persist in their host genome by generating new copies to outpace natural selection. Transposable elements exert a large influence on host genome evolution, in some cases providing adaptive changes. Here we measure the fitness effect of the transposable element insertions in the fission yeast Schizosaccharomyces pombe type strain by removing all insertions of its only native transposable element family, the long terminal repeat retrotransposon Tf2. We show that Tf2 elements provide a positive fitness contribution to its host. Tf2 ablation results in changes to the regulation of a mitochondrial gene and, consistently, the fitness effect are sensitive to growth conditions. We propose that Tf2 influences host fitness in a directed manner by dynamically rewiring the transcriptional response to metabolic stress.

Detection of mosaic and population-level structural variants with Sniffles2

Calling structural variations (SVs) is technically challenging, but using long reads remains the most accurate way to identify complex genomic alterations. Here we present Sniffles2, which improves over current methods by implementing a repeat aware clustering coupled with a fast consensus sequence and coverage-adaptive filtering. Sniffles2 is 11.8 times faster and 29% more accurate than state-of-the-art SV callers across different coverages (5-50×), sequencing technologies (ONT and HiFi) and SV types. Furthermore, Sniffles2 solves the problem of family-level to population-level SV calling to produce fully genotyped VCF files. Across 11 probands, we accurately identified causative SVs around MECP2, including highly complex alleles with three overlapping SVs. Sniffles2 also enables the detection of mosaic SVs in bulk long-read data. As a result, we identified multiple mosaic SVs in brain tissue from a patient with multiple system atrophy. The identified SV showed a remarkable diversity within the cingulate cortex, impacting both genes involved in neuron function and repetitive elements.

Subtelomeres: hotspots of genome variation

Eukaryotic cells contain multiple types of duplicated sequences. Typical examples are tandem repeat sequences including telomeres, centromeres, rDNA genes and transposable elements. Most of these sequences are unstable; thus, their copy numbers or sequences change rapidly in the course of evolution. In this review, I will describe roles of subtelomere regions, which are located adjacent to telomeres at chromosome ends, and recent discoveries about their sequence variation.

The Fission Yeast Mating-Type Switching Motto: "One-for-Two" and "Two-for-One"

Schizosaccharomyces pombe is an ascomycete fungus that divides by medial fission; it is thus commonly referred to as fission yeast, as opposed to the distantly related budding yeast Saccharomyces cerevisiae. The reproductive lifestyle of S. pombe relies on an efficient genetic sex determination system generating a 1:1 sex ratio and using alternating haploid/diploid phases in response to environmental conditions. In this review, we address how one haploid cell manages to generate two sister cells with opposite mating types, a prerequisite to conjugation and meiosis. This mating-type switching process depends on two highly efficient consecutive asymmetric cell divisions that rely on DNA replication, repair, and recombination as well as the structure and components of heterochromatin. We pay special attention to the intimate interplay between the genetic and epigenetic partners involved in this process to underscore the importance of basic research and its profound implication for a better understanding of chromatin biology.

Apparent Absence of Selective Pressure on Pneumocystis jirovecii Organisms in Patients with Prior Methotrexate Exposure

Pneumocystis jirovecii infections occur in patients treated with methotrexate (MTX) because of immunosuppressive effects of this highly potent dihydrofolate reductase (DHFR) inhibitor. Conversely, MTX may act as an anti-P. jirovecii drug and consequently may exert a selective pressure on this fungus. In this context, we compared the sequences of the dhfr gene of P. jirovecii isolates obtained from two groups of patients with P. jirovecii infections. The first group, with systemic diseases or malignancies, had prior exposure to MTX (21 patients), whereas the second group (22 patients), the control group, did not. Three single nucleotide polymorphisms (SNPs) were observed at positions 278, 312, and 381. The first one was located in the intronic region and the two others were synonymous. Based on these SNPs, three P. jirovecii dhfr alleles, named A, B, and C, were specified. Allele A was the most frequent, as it was observed in 18 patients (85.7%) and in 16 patients (72.7%) of the first and second groups, respectively. No significant difference in P. jirovecii dhfr gene diversity in the two patient groups was observed. In conclusion, these original results suggest that MTX does not exert an overt selective pressure on P. jirovecii organisms.

The evolutionary and ecological potential of yeast hybrids

Recent findings in yeast genetics and genomics have advanced our understanding of the evolutionary potential unlocked by hybridization, especially in the genus Saccharomyces. We now have a clearer picture of the prevalence of yeast hybrids in the environment, their ecological and evolutionary history, and the genetic mechanisms driving (and constraining) their adaptation. Here, we describe how the instability of hybrid genomes determines fitness across large evolutionary scales, highlight new hybrid strain engineering techniques, and review tools for comparative hybrid genome analysis. The recent push to take yeast research back 'into the wild' has resulted in new genomic and ecological resources. These provide an arena for quantitative genetics and allow us to investigate the architecture of complex traits and mechanisms of adaptation to rapidly changing environments. The vast genetic diversity of hybrid populations can yield insights beyond those possible with isogenic lines. Hybrids offer a limitless supply of genetic variation that can be tapped for industrial strain improvement but also, combined with experimental evolution, can be used to predict population responses to future climate change - a fundamental task for biologists.

Identification of an integrase-independent pathway of retrotransposition

Retroviruses and long terminal repeat retrotransposons rely on integrase (IN) to insert their complementary DNA (cDNA) into the genome of host cells. Nevertheless, in the absence of IN, retroelements can retain notable levels of insertion activity. We have characterized the IN-independent pathway of Tf1 and found that insertion sites had homology to the primers of reverse transcription, which form single-stranded DNAs at the termini of the cDNA. In the absence of IN activity, a similar bias was observed with HIV-1. Our studies showed that the Tf1 insertions result from single-strand annealing, a noncanonical form of homologous recombination mediated by Rad52. By expanding our analysis of insertions to include repeat sequences, we found most formed tandem elements by inserting at preexisting copies of a related transposable element. Unexpectedly, we found that wild-type Tf1 uses the IN-independent pathway as an alternative mode of insertion.

The teenage years of yeast population genomics trace history, admixing and getting wilder

Population genomics studies the evolutionary processes that shape intraspecies genetic variations. In this review, I explore the insights into yeast-population genomics that have emerged from recent advances in sequencing. Genomes of the model Saccharomyces cerevisiae and many new yeast species from around the world are being used to address various aspects of population biology, including geographical origin, the level of introgression, domestication signatures, and outcrossing frequency. New long-read sequencing has enabled a greater capacity to quantify these variations at a finer resolution from complete de novo genomes at the population scale to phasing subgenomes of different origins. These resources provide a platform to dissect the relationship between phenotypes across environmental niches.

Allorecognition genes drive reproductive isolation in Podospora anserina

Allorecognition, the capacity to discriminate self from conspecific non-self, is a ubiquitous organismal feature typically governed by genes evolving under balancing selection. Here, we show that in the fungus Podospora anserina, allorecognition loci controlling vegetative incompatibility (het genes), define two reproductively isolated groups through pleiotropic effects on sexual compatibility. These two groups emerge from the antagonistic interactions of the unlinked loci het-r (encoding a NOD-like receptor) and het-v (encoding a methyltransferase and an MLKL/HeLo domain protein). Using a combination of genetic and ecological data, supported by simulations, we provide a concrete and molecularly defined example whereby the origin and coexistence of reproductively isolated groups in sympatry is driven by pleiotropic genes under balancing selection.

Insights into the ecology of Schizosaccharomyces species in natural and artificial habitats

The fission yeast genus Schizosaccharomyces contains important model organisms for biological research. In particular, S. pombe is a widely used model eukaryote. So far little is known about the natural and artificial habitats of species in this genus. Finding out where S. pombe and other fission yeast species occur and how they live in their habitats can promote better understanding of their biology. Here we investigate in which substrates S. pombe, S. octosporus, S. osmophilus and S. japonicus are present. To this end about 2100 samples consisting of soil, tree sap fluxes, fresh fruit, dried fruit, honey, cacao beans, molasses and other substrates were analyzed. Effective isolation methods that allow efficient isolation of the above mentioned species were developed. Based on the frequency of isolating different fission yeast species in various substrates and on extensive literature survey, conclusions are drawn on their ecology. The results suggest that the primary habitat of S. pombe and S. octosporus is honeybee honey. Both species were also frequently detected on certain dried fruit like raisins, mango or pineapple to which they could be brought by the honey bees during ripening or during drying. While S. pombe was regularly isolated from grape mash and from fermented raw cacao beans S. octosporus was never isolated from fresh fruit. The main habitat of S. osmophilus seems to be solitary bee beebread. It was rarely isolated from raisins. S. japonicus was mainly found in forest substrates although it occurs on fruit and in fruit fermentations, too.

Ineffective silencing of transposable elements on an avian W Chromosome

One of the defining features of transposable elements (TEs) is their ability to move to new locations in the host genome. To minimise the potentially deleterious effects of de novo TE insertions, hosts have evolved several mechanisms to control TE activity, including recombination-mediated removal and epigenetic silencing; however, increasing evidence suggests that silencing of TEs is often incomplete. The crow family experienced a recent radiation of LTR retrotransposons (LTRs), offering an opportunity to gain insight into the regulatory control of young, potentially still active TEs. We quantified the abundance of TE-derived transcripts across several tissues in 15 Eurasian crows (Corvus (corone) spp.) raised under common garden conditions and find evidence for ineffective TE suppression on the female-specific W Chromosome. Using RNA-seq data, we show that ~ 9.5% of all transcribed TEs had considerably greater (average: 16-fold) transcript abundance in female crows, and that more than 85% of these female-biased TEs originated on the W Chromosome. After accounting for differences in TE density among chromosomal classes, W-linked TEs were significantly more highly expressed than TEs residing on other chromosomes, consistent with ineffective silencing on the former. Together, our results suggest that the crow W Chromosome acts as a source of transcriptionally active TEs, with possible negative fitness consequences for female birds analogous to Drosophila (an X/Y system), where overexpression of Y-linked TEs is associated with male-specific aging and fitness loss ('toxic Y').

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.

Widespread introgression across a phylogeny of 155 Drosophila genomes

Genome-scale sequence data have invigorated the study of hybridization and introgression, particularly in animals. However, outside of a few notable cases, we lack systematic tests for introgression at a larger phylogenetic scale across entire clades. Here, we leverage 155 genome assemblies from 149 species to generate a fossil-calibrated phylogeny and conduct multilocus tests for introgression across 9 monophyletic radiations within the genus Drosophila. Using complementary phylogenomic approaches, we identify widespread introgression across the evolutionary history of Drosophila. Mapping gene-tree discordance onto the phylogeny revealed that both ancient and recent introgression has occurred across most of the 9 clades that we examined. Our results provide the first evidence of introgression occurring across the evolutionary history of Drosophila and highlight the need to continue to study the evolutionary consequences of hybridization and introgression in this genus and across the tree of life.

The search for Schizosaccharomyces fission yeasts in environmental meta-transcriptomes

Fission yeast is an important model organism in evolutionary genetics and cell biology research. Nevertheless, most research is limited to a single laboratory strain and knowledge of its natural occurrence is limited, which reduces our understanding of its life history and hinders isolation of new strains from nature. Understanding the natural diversity of fission yeast can provide insight into its genetic and phenotypic diversity and the evolutionary processes that shaped these. Here we aimed to identify candidate natural habitats of fission yeasts by searching through a large collection of publicly available environmental metatranscriptomic datasets. Using a custom pipeline, we processed over 13,000 NCBI SRA accessions, from a wide range of 34 different environmental categories. Overall, we found a very low abundance of putative yeast transcripts, with most fission yeast signatures coming from the categories of 'food' and 'terrestrial arthropods'. Additionally, a signal could be found in a variety of marine and fresh aquatic habitats. Our results do not provide a conclusive answer on the natural habitat of fission yeasts, but our analysis further narrows the range of locations where fission yeasts naturally occur.

Molecular Mechanisms and Evolutionary Consequences of Spore Killers in Ascomycetes

In this review, we examine the fungal spore killers. These are meiotic drive elements that cheat during sexual reproduction to increase their transmission into the next generation. Spore killing has been detected in a number of ascomycete genera, including Podospora, Neurospora, Schizosaccharomyces, Bipolaris, and Fusarium. There have been major recent advances in spore killer research that have increased our understanding of the molecular identity, function, and evolutionary history of the known killers. The spore killers vary in the mechanism by which they kill and are divided into killer-target and poison-antidote drivers. In killer-target systems, the drive locus encodes an element that can be described as a killer, while the target is an allele found tightly linked to the drive locus but on the nondriving haplotype. The poison-antidote drive systems encode both a poison and an antidote element within the drive locus. The key to drive in this system is the restricted distribution of the antidote: only the spores that inherit the drive locus receive the antidote and are rescued from the toxicity of the poison. Spore killers also vary in their genome architecture and can consist of a single gene or multiple linked genes. Due to their ability to distort meiosis, spore killers gain a selective advantage at the gene level that allows them to increase in frequency in a population over time, even if they reduce host fitness, and they may have significant impact on genome architecture and macroevolutionary processes such as speciation.

Reactivation of transposable elements following hybridization in fission yeast

Hybridization is thought to reactivate transposable elements (TEs) that were efficiently suppressed in the genomes of the parental hosts. Here, we provide evidence for this “genomic shock hypothesis” in the fission yeast Schizosaccharomyces pombe. In this species, two divergent lineages (Sp and Sk) have experienced recent, likely human-induced, hybridization. We used long-read sequencing data to assemble genomes of 37 samples derived from 31 S. pombe strains spanning a wide range of ancestral admixture proportions. A comprehensive TE inventory revealed exclusive presence of long terminal repeat (LTR) retrotransposons. Sequence analysis of active full-length elements, as well as solo LTRs, revealed a complex history of homologous recombination. Population genetic analyses of syntenic sequences placed insertion of many solo LTRs before the split of the Sp and Sk lineages. Most full-length elements were inserted more recently, after hybridization. With the exception of a single full-length element with signs of positive selection, both solo LTRs and, in particular, full-length elements carry signatures of purifying selection indicating effective removal by the host. Consistent with reactivation upon hybridization, the number of full-length LTR retrotransposons, varying extensively from zero to 87 among strains, significantly increases with the degree of genomic admixture. This study gives a detailed account of global TE diversity in S. pombe, documents complex recombination histories within TE elements, and provides evidence for the “genomic shock hypothesis.”

Diverse mating phenotypes impact the spread of wtf meiotic drivers in Schizosaccharomyces pombe

Meiotic drivers are genetic elements that break Mendel's law of segregation to be transmitted into more than half of the offspring produced by a heterozygote. The success of a driver relies on outcrossing (mating between individuals from distinct lineages) because drivers gain their advantage in heterozygotes. It is, therefore, curious that Schizosaccharomyces pombe, a species reported to rarely outcross, harbors many meiotic drivers. To address this paradox, we measured mating phenotypes in S. pombe natural isolates. We found that the propensity for cells from distinct clonal lineages to mate varies between natural isolates and can be affected both by cell density and by the available sexual partners. Additionally, we found that the observed levels of preferential mating between cells from the same clonal lineage can slow, but not prevent, the spread of a wtf meiotic driver in the absence of additional fitness costs linked to the driver. These analyses reveal parameters critical to understanding the evolution of S. pombe and help explain the success of meiotic drivers in this species.

Trends in yeast diversity discovery

Yeasts, usually defined as unicellular fungi, occur in various fungal lineages. Hence, they are not a taxonomic unit, but rather represent a fungal lifestyle shared by several unrelated lineages. Although the discovery of new yeast species occurs at an increasing speed, at the current rate it will likely take hundreds of years, if ever, before they will all be documented. Many parts of the earth, including many threatened habitats, remain unsampled for yeasts and many others are only superficially studied. Cold habitats, such as glaciers, are home to a specific community of cold-adapted yeasts, and, hence, there is some urgency to study such environments at locations where they might disappear soon due to anthropogenic climate change. The same is true for yeast communities in various natural forests that are impacted by deforestation and forest conversion. Many countries of the so-called Global South have not been sampled for yeasts, despite their economic promise. However, extensive research activity in Asia, especially China, has yielded many taxonomic novelties. Comparative genomics studies have demonstrated the presence of yeast species with a hybrid origin, many of them isolated from clinical or industrial environments. DNA-metabarcoding studies have demonstrated the prevalence, and in some cases dominance, of yeast species in soils and marine waters worldwide, including some surprising distributions, such as the unexpected and likely common presence of Malassezia yeasts in marine habitats.

Segregating Complete Tf2 Elements Are Largely Neutral in Fission Yeast

Transposable elements (TEs) comprise a large proportion of the eukaryote genomes. Yet it remains poorly understood how TEs influence the fitness of the hosts carrying them. Here, we empirically test the impact of TEs on the host fitness in the fission yeast Schizosaccharomyces pombe. We find that two families of TEs (Tf1 and Tf2 elements), both of which belong to long terminal repeat retrotransposons, are highly polymorphic among individual S. pombe strains. Only 13 complete Tf2 elements are identified in S. pombe laboratory strain 972. These 13 Tf2 elements integrated into host genomes in very recent time and are segregating within the S. pombe population. Through knocking out each of the 13 Tf2 elements in S. pombe strain 972, we find Tf2 knockout does not affect the host fitness, and Tf2 elements do not alter the expression of nearby genes. Challenged by diverse forms of stress, the Tf2 knockout strains do not exhibit different growth rates from wild-type strain. Together, we conclude that segregating complete Tf2 elements insertions are largely neutral to host fitness in the fission yeast. Our study provides genome-wide empirical support for the selfish nature of TEs in fission yeast.

Genome Comparisons of the Fission Yeasts Reveal Ancient Collinear Loci Maintained by Natural Selection

Fission yeasts have a unique life history and exhibit distinct evolutionary patterns from other yeasts. Besides, the species demonstrate stable genome structures despite the relatively fast evolution of their genomic sequences. To reveal what could be the reason for that, comparative genomic analyses were carried out. Our results provided evidence that the structural and sequence evolution of the fission yeasts were correlated. Moreover, we revealed ancestral locally collinear blocks (aLCBs), which could have been inherited from their last common ancestor. These aLCBs proved to be the most conserved regions of the genomes as the aLCBs contain almost eight genes/blocks on average in the same orientation and order across the species. Gene order of the aLCBs is mainly fission-yeast-specific but supports the idea of filamentous ancestors. Nevertheless, the sequences and gene structures within the aLCBs are as mutable as any sequences in other parts of the genomes. Although genes of certain Gene Ontology (GO) categories tend to cluster at the aLCBs, those GO enrichments are not related to biological functions or high co-expression rates, they are, rather, determined by the density of essential genes and Rec12 cleavage sites. These data and our simulations indicated that aLCBs might not only be remnants of ancestral gene order but are also maintained by natural selection.

Hi-C scaffolded short- and long-read genome assemblies of the California sea lion are broadly consistent for syntenic inference across 45 million years of evolution

With the advent of chromatin-interaction maps, chromosome-level genome assemblies have become a reality for a wide range of organisms. Scaffolding quality is, however, difficult to judge. To explore this gap, we generated multiple chromosome-scale genome assemblies of an emerging wild animal model for carcinogenesis, the California sea lion (Zalophus californianus). Short-read assemblies were scaffolded with two independent chromatin interaction mapping data sets (Hi-C and Chicago), and long-read assemblies with three data types (Hi-C, optical maps and 10X linked reads) following the "Vertebrate Genomes Project (VGP)" pipeline. In both approaches, 18 major scaffolds recovered the karyotype (2n = 36), with scaffold N50s of 138 and 147 Mb, respectively. Synteny relationships at the chromosome level with other pinniped genomes (2n = 32-36), ferret (2n = 34), red panda (2n = 36) and domestic dog (2n = 78) were consistent across approaches and recovered known fissions and fusions. Comparative chromosome painting and multicolour chromosome tiling with a panel of 264 genome-integrated single-locus canine bacterial artificial chromosome probes provided independent evaluation of genome organization. Broad-scale discrepancies between the approaches were observed within chromosomes, most commonly in translocations centred around centromeres and telomeres, which were better resolved in the VGP assembly. Genomic and cytological approaches agreed on near-perfect synteny of the X chromosome, and in combination allowed detailed investigation of autosomal rearrangements between dog and sea lion. This study presents high-quality genomes of an emerging cancer model and highlights that even highly fragmented short-read assemblies scaffolded with Hi-C can yield reliable chromosome-level scaffolds suitable for comparative genomic analyses.

Selective sorting of ancestral introgression in maize and teosinte along an elevational cline

While often deleterious, hybridization can also be a key source of genetic variation and pre-adapted haplotypes, enabling rapid evolution and niche expansion. Here we evaluate these opposing selection forces on introgressed ancestry between maize (Zea mays ssp. mays) and its wild teosinte relative, mexicana (Zea mays ssp. mexicana). Introgression from ecologically diverse teosinte may have facilitated maize's global range expansion, in particular to challenging high elevation regions (> 1500 m). We generated low-coverage genome sequencing data for 348 maize and mexicana individuals to evaluate patterns of introgression in 14 sympatric population pairs, spanning the elevational range of mexicana, a teosinte endemic to the mountains of Mexico. While recent hybrids are commonly observed in sympatric populations and mexicana demonstrates fine-scale local adaptation, we find that the majority of mexicana ancestry tracts introgressed into maize over 1000 generations ago. This mexicana ancestry seems to have maintained much of its diversity and likely came from a common ancestral source, rather than contemporary sympatric populations, resulting in relatively low FST between mexicana ancestry tracts sampled from geographically distant maize populations. Introgressed mexicana ancestry in maize is reduced in lower-recombination rate quintiles of the genome and around domestication genes, consistent with pervasive selection against introgression. However, we also find mexicana ancestry increases across the sampled elevational gradient and that high introgression peaks are most commonly shared among high-elevation maize populations, consistent with introgression from mexicana facilitating adaptation to the highland environment. In the other direction, we find patterns consistent with adaptive and clinal introgression of maize ancestry into sympatric mexicana at many loci across the genome, suggesting that maize also contributes to adaptation in mexicana, especially at the lower end of its elevational range. In sympatric maize, in addition to high introgression regions we find many genomic regions where selection for local adaptation maintains steep gradients in introgressed mexicana ancestry across elevation, including at least two inversions: the well-characterized 14 Mb Inv4m on chromosome 4 and a novel 3 Mb inversion Inv9f surrounding the macrohairless1 locus on chromosome 9. Most outlier loci with high mexicana introgression show no signals of sweeps or local sourcing from sympatric populations and so likely represent ancestral introgression sorted by selection, resulting in correlated but distinct outcomes of introgression in different contemporary maize populations.

PRINCESS: comprehensive detection of haplotype resolved SNVs, SVs, and methylation

Long-read sequencing has been shown to have advantages in structural variation (SV) detection and methylation calling. Many studies focus either on SV, methylation, or phasing of SNV; however, only the combination of variants provides a comprehensive insight into the sample and thus enables novel findings in biology or medicine. PRINCESS is a structured workflow that takes raw sequence reads and generates a fully phased SNV, SV, and methylation call set within a few hours. PRINCESS achieves high accuracy and long phasing even on low coverage datasets and can resolve repetitive, complex medical relevant genes that often escape detection. PRINCESS is publicly available at https://github.com/MeHelmy/princess under the MIT license.

Towards population-scale long-read sequencing

Long-read sequencing technologies have now reached a level of accuracy and yield that allows their application to variant detection at a scale of tens to thousands of samples. Concomitant with the development of new computational tools, the first population-scale studies involving long-read sequencing have emerged over the past 2 years and, given the continuous advancement of the field, many more are likely to follow. In this Review, we survey recent developments in population-scale long-read sequencing, highlight potential challenges of a scaled-up approach and provide guidance regarding experimental design. We provide an overview of current long-read sequencing platforms, variant calling methodologies and approaches for de novo assemblies and reference-based mapping approaches. Furthermore, we summarize strategies for variant validation, genotyping and predicting functional impact and emphasize challenges remaining in achieving long-read sequencing at a population scale.

Invasion and maintenance of meiotic drivers in populations of ascomycete fungi

Meiotic drivers (MDs) are selfish genetic elements that are able to become overrepresented among the products of meiosis. This transmission advantage makes it possible for them to spread in a population even when they impose fitness costs on their host organisms. Whether an MD can invade a population, and subsequently reach fixation or coexist in a stable polymorphism, depends on the one hand on the biology of the host organism, including its life cycle, mating system, and population structure, and on the other hand on the specific fitness effects of the driving allele on the host. Here, we present a population genetic model for spore killing, a type of drive specific to fungi. We show how ploidy level, rate of selfing, and efficiency of spore killing affect the invasion probability of a driving allele and the conditions for its stable coexistence with a nondriving allele. Our model can be adapted to different fungal life cycles, and is applied here to two well-studied genera of filamentous ascomycetes known to harbor spore-killing elements, Neurospora and Podospora. We discuss our results in the light of recent empirical findings for these two systems.

Complete sequences of Schizosaccharomyces pombe subtelomeres reveal multiple patterns of genome variation

Genome sequences have been determined for many model organisms; however, repetitive regions such as centromeres, telomeres, and subtelomeres have not yet been sequenced completely. Here, we report the complete sequences of subtelomeric homologous (SH) regions of the fission yeast Schizosaccharomyces pombe. We overcame technical difficulties to obtain subtelomeric repetitive sequences by constructing strains that possess single SH regions of a standard laboratory strain. In addition, some natural isolates of S. pombe were analyzed using previous sequencing data. Whole sequences of SH regions revealed that each SH region consists of two distinct parts with mosaics of multiple common segments or blocks showing high variation among subtelomeres and strains. Subtelomere regions show relatively high frequency of nucleotide variations among strains compared with the other chromosomal regions. Furthermore, we identified subtelomeric RecQ-type helicase genes, tlh3 and tlh4, which add to the already known tlh1 and tlh2, and found that the tlh1-4 genes show high sequence variation with missense mutations, insertions, and deletions but no severe effects on their RNA expression. Our results indicate that SH sequences are highly polymorphic and hot spots for genome variation. These features of subtelomeres may have contributed to genome diversity and, conversely, various diseases.

Vicariance followed by secondary gene flow in a young gazelle species complex

Grant's gazelles have recently been proposed to be a species complex comprising three highly divergent mtDNA lineages (Nanger granti, N. notata and N. petersii). The three lineages have nonoverlapping distributions in East Africa, but without any obvious geographical divisions, making them an interesting model for studying the early‐stage evolutionary dynamics of allopatric speciation in detail. Here, we use genomic data obtained by restriction site‐associated (RAD) sequencing of 106 gazelle individuals to shed light on the evolutionary processes underlying Grant's gazelle divergence, to characterize their genetic structure and to assess the presence of gene flow between the main lineages in the species complex. We date the species divergence to 134,000 years ago, which is recent in evolutionary terms. We find population subdivision within N. granti, which coincides with the previously suggested two subspecies, N. g. granti and N. g. robertsii. Moreover, these two lineages seem to have hybridized in Masai Mara. Perhaps more surprisingly given their extreme genetic differentiation, N. granti and N. petersii also show signs of prolonged admixture in Mkomazi, which we identified as a hybrid population most likely founded by allopatric lineages coming into secondary contact. Despite the admixed composition of this population, elevated X chromosomal differentiation suggests that selection may be shaping the outcome of hybridization in this population. Our results therefore provide detailed insights into the processes of allopatric speciation and secondary contact in a recently radiated species complex.

Hybridization and cryptic speciation in the Iberian endemic plant genus Phalacrocarpum (Asteraceae-Anthemideae)

Understanding the role and impact of reticulation in phylogenetic inquiry has improved with extended use of high throughput sequencing data. Yet, due to the dynamism of genomes over evolutionary time, disentangling old hybridization events remains a serious challenge. Phalacrocarpum (DC.) Willk. is one of the 27 Iberian endemic plant genera, currently considered monotypic but including three subspecies. Its uncertain phylogenetic relationships within tribe Anthemideae (Asteraceae) point to an Early Miocene divergence from its sister group, and its persistent taxonomic instability has been proposed to be due to hybridization. We aim at understanding the evolutionary history of this genus using SNPs called from a genotyping-by-sequencing (GBS) analysis, Sanger sequences-from three plastid DNA regions (psbJ-petA, petB-petD, trnH-psbA) and the nuclear ribosomal ITS regions (cloned)-as well as leaf morphometric multivariate analysis. SNP data and Sanger sequences strongly support the unforeseen existence of a cryptic species in the eastern populations of P. oppositifolium subsp. anomalum. Broad molecular and morphometric patterns of variation found in conflictive populations from the Sanabria Valley region convincingly identify a recent previously undocumented hybrid zone. By contrast, evidence is less conclusive on relationships between subspecies hoffmannseggii, oppositifolium and a second conflictive group distributed along the Galician-Portuguese border (Orense massifs). Although genetic clustering analysis of SNP data suggests that the former subspecies was the maternal progenitor in hybridization events that gave rise to the other two groups, we found considerable uniqueness of ITS ribotypes and plastid haplotypes in them. This result, in the context of Pleistocene climatically-driven range shifts in NW Iberian Peninsula, can be due to periods of isolation, genetic bottlenecks and drift superimposed on old hybridization events. Our study confirms the idea that unravelling old hybridization events may be compromised by the suite of evolutionary processes accumulated subsequently, particularly in areas with a history of climatic instability.

Higher Gene Flow in Sex-Related Chromosomes than in Autosomes during Fungal Divergence

Nonrecombining sex chromosomes are widely found to be more differentiated than autosomes among closely related species, due to smaller effective population size and/or to a disproportionally large-X effect in reproductive isolation. Although fungal mating-type chromosomes can also display large nonrecombining regions, their levels of differentiation compared with autosomes have been little studied. Anther-smut fungi from the Microbotryum genus are castrating pathogens of Caryophyllaceae plants with largely nonrecombining mating-type chromosomes. Using whole genome sequences of 40 fungal strains, we quantified genetic differentiation among strains isolated from the geographically overlapping North American species and subspecies of Silene virginica and S. caroliniana. We inferred that gene flow likely occurred at the early stages of divergence and then completely stopped. We identified large autosomal genomic regions with chromosomal inversions, with higher genetic divergence than the rest of the genomes and highly enriched in selective sweeps, supporting a role of rearrangements in preventing gene flow in genomic regions involved in ecological divergence. Unexpectedly, the nonrecombining mating-type chromosomes showed lower divergence than autosomes due to higher gene flow, which may be promoted by adaptive introgressions of less degenerated mating-type chromosomes. The fact that both mating-type chromosomes are always heterozygous and nonrecombining may explain such patterns that oppose to those found for XY or ZW sex chromosomes. The specific features of mating-type chromosomes may also apply to the UV sex chromosomes determining sexes at the haploid stage in algae and bryophytes and may help test general hypotheses on the evolutionary specificities of sex-related chromosomes.

Atypical meiosis can be adaptive in outcrossed Schizosaccharomyces pombe due to wtf meiotic drivers

Killer meiotic drivers are genetic parasites that destroy ‘sibling’ gametes lacking the driver allele. The fitness costs of drive can lead to selection of unlinked suppressors. This suppression could involve evolutionary tradeoffs that compromise gametogenesis and contribute to infertility. Schizosaccharomyces pombe, an organism containing numerous gamete (spore)-killing wtf drivers, offers a tractable system to test this hypothesis. Here, we demonstrate that in scenarios analogous to outcrossing, wtf drivers generate a fitness landscape in which atypical spores, such as aneuploids and diploids, are advantageous. In this context, wtf drivers can decrease the fitness costs of mutations that disrupt meiotic fidelity and, in some circumstances, can even make such mutations beneficial. Moreover, we find that S. pombe isolates vary greatly in their ability to make haploid spores, with some isolates generating up to 46% aneuploid or diploid spores. This work empirically demonstrates the potential for meiotic drivers to shape the evolution of gametogenesis.

Advances in understanding the evolution of fungal genome architecture

Diversity within the fungal kingdom is evident from the wide range of morphologies fungi display as well as the various ecological roles and industrial purposes they serve. Technological advances, particularly in long-read sequencing, coupled with the increasing efficiency and decreasing costs across sequencing platforms have enabled robust characterization of fungal genomes. These sequencing efforts continue to reveal the rampant diversity in fungi at the genome level. Here, we discuss studies that have furthered our understanding of fungal genetic diversity and genomic evolution. These studies revealed the presence of both small-scale and large-scale genomic changes. In fungi, research has recently focused on many small-scale changes, such as how hypermutation and allelic transmission impact genome evolution as well as how and why a few specific genomic regions are more susceptible to rapid evolution than others. High-throughput sequencing of a diverse set of fungal genomes has also illuminated the frequency, mechanisms, and impacts of large-scale changes, which include chromosome structural variation and changes in chromosome number, such as aneuploidy, polyploidy, and the presence of supernumerary chromosomes. The studies discussed herein have provided great insight into how the architecture of the fungal genome varies within species and across the kingdom and how modern fungi may have evolved from the last common fungal ancestor and might also pave the way for understanding how genomic diversity has evolved in all domains of life.

Near-chromosome level genome assembly of the fruit pest Drosophila suzukii using long-read sequencing

Over the past decade, the spotted wing Drosophila, Drosophila suzukii, has invaded Europe and America and has become a major agricultural pest in these areas, thereby prompting intense research activities to better understand its biology. Two draft genome assemblies already exist for this species but contain pervasive assembly errors and are highly fragmented, which limits their values. Our purpose here was to improve the assembly of the D. suzukii genome and to annotate it in a way that facilitates comparisons with D. melanogaster. For this, we generated PacBio long-read sequencing data and assembled a novel, high-quality D. suzukii genome assembly. It is one of the largest Drosophila genomes, notably because of the expansion of its repeatome. We found that despite 16 rounds of full-sib crossings the D. suzukii strain that we sequenced has maintained high levels of polymorphism in some regions of its genome. As a consequence, the quality of the assembly of these regions was reduced. We explored possible origins of this high residual diversity, including the presence of structural variants and a possible heterogeneous admixture pattern of North American and Asian ancestry. Overall, our assembly and annotation constitute a high-quality genomic resource that can be used for both high-throughput sequencing approaches, as well as manipulative genetic technologies to study D. suzukii.

Discovery and population genomics of structural variation in a songbird genus

Structural variation (SV) constitutes an important type of genetic mutations providing the raw material for evolution. Here, we uncover the genome-wide spectrum of intra- and interspecific SV segregating in natural populations of seven songbird species in the genus Corvus. Combining short-read (N = 127) and long-read re-sequencing (N = 31), as well as optical mapping (N = 16), we apply both assembly- and read mapping approaches to detect SV and characterize a total of 220,452 insertions, deletions and inversions. We exploit sampling across wide phylogenetic timescales to validate SV genotypes and assess the contribution of SV to evolutionary processes in an avian model of incipient speciation. We reveal an evolutionary young (~530,000 years) cis-acting 2.25-kb LTR retrotransposon insertion reducing expression of the NDP gene with consequences for premating isolation. Our results attest to the wealth and evolutionary significance of SV segregating in natural populations and highlight the need for reliable SV genotyping.

Genomic evidence for a hybrid origin of the yeast opportunistic pathogen Candida albicans

Background

Opportunistic yeast pathogens of the genus Candida are an important medical problem. Candida albicans, the most prevalent Candida species, is a natural commensal of humans that can adopt a pathogenic behavior. This species is highly heterozygous and cannot undergo meiosis, adopting instead a parasexual cycle that increases genetic variability and potentially leads to advantages under stress conditions. However, the origin of C. albicans heterozygosity is unknown, and we hypothesize that it could result from ancestral hybridization. We tested this idea by analyzing available genomes of C. albicans isolates and comparing them to those of hybrid and non-hybrid strains of other Candida species.

Results

Our results show compelling evidence that C. albicans is an evolved hybrid. The genomic patterns observed in C. albicans are similar to those of other hybrids such as Candida orthopsilosis MCO456 and Candida inconspicua, suggesting that it also descends from a hybrid of two divergent lineages. Our analysis indicates that most of the divergence between haplotypes in C. albicans heterozygous blocks was already present in a putative heterozygous ancestor, with an estimated 2.8% divergence between homeologous chromosomes. The levels and patterns of ancestral heterozygosity found cannot be fully explained under the paradigm of vertical evolution and are not consistent with continuous gene flux arising from lineage-specific events of admixture.

Conclusions

Although the inferred level of sequence divergence between the putative parental lineages (2.8%) is not clearly beyond current species boundaries in Saccharomycotina, we show here that all analyzed C. albicans strains derive from a single hybrid ancestor and diverged by extensive loss of heterozygosity. This finding has important implications for our understanding of C. albicans evolution, including the loss of the sexual cycle, the origin of the association with humans, and the evolution of virulence traits.

Intraspecific Diversity of Fission Yeast Mitochondrial Genomes

The fission yeast Schizosaccharomyces pombe is an important model organism, but its natural diversity and evolutionary history remain under-studied. In particular, the population genomics of the S. pombe mitochondrial genome (mitogenome) has not been thoroughly investigated. Here, we assembled the complete circular-mapping mitogenomes of 192 S. pombe isolates de novo, and found that these mitogenomes belong to 69 nonidentical sequence types ranging from 17,618 to 26,910 bp in length. Using the assembled mitogenomes, we identified 20 errors in the reference mitogenome and discovered two previously unknown mitochondrial introns. Analyzing sequence diversity of these 69 types of mitogenomes revealed two highly distinct clades, with only three mitogenomes exhibiting signs of inter-clade recombination. This diversity pattern suggests that currently available S. pombe isolates descend from two long-separated ancestral lineages. This conclusion is corroborated by the diversity pattern of the recombination-repressed K-region located between donor mating-type loci mat2 and mat3 in the nuclear genome. We estimated that the two ancestral S. pombe lineages diverged about 31 million generations ago. These findings shed new light on the evolution of S. pombe and the data sets generated in this study will facilitate future research on genome evolution.

Dramatically diverse Schizosaccharomyces pombe wtf meiotic drivers all display high gamete-killing efficiency

Meiotic drivers are selfish alleles that can force their transmission into more than 50% of the viable gametes made by heterozygotes. Meiotic drivers are known to cause infertility in a diverse range of eukaryotes and are predicted to affect the evolution of genome structure and meiosis. The wtf gene family of Schizosaccharomyces pombe includes both meiotic drivers and drive suppressors and thus offers a tractable model organism to study drive systems. Currently, only a handful of wtf genes have been functionally characterized and those genes only partially reflect the diversity of the wtf gene family. In this work, we functionally test 22 additional wtf genes for meiotic drive phenotypes. We identify eight new drivers that share between 30–90% amino acid identity with previously characterized drivers. Despite the vast divergence between these genes, they generally drive into >85% of gametes when heterozygous. We also identify three wtf genes that suppress other wtf drivers, including two that also act as autonomous drivers. Additionally, we find that wtf genes do not underlie a weak (64% allele transmission) meiotic driver on chromosome 1. Finally, we find that some Wtf proteins have expression or localization patterns that are distinct from the poison and antidote proteins encoded by drivers and suppressors, suggesting some wtf genes may have non-meiotic drive functions. Overall, this work expands our understanding of the wtf gene family and the burden selfish driver genes impose on S. pombe.

During gametogenesis, the two gene copies at a given locus, known as alleles, are each transmitted to 50% of the gametes (e.g. sperm). However, some alleles cheat so that they are found in more than the expected 50% of gametes, often at the expense of fertility. This selfish behavior is known as meiotic drive. Some members of the wtf gene family in the fission yeast Schizosaccharomyces pombe kill the gametes (spores) that do not inherit them, resulting in meiotic drive favoring the wtf allele. Other wtf genes act as suppressors of drive. However, the wtf gene family is diverse and only a small subset of the genes has been characterized. Here we analyze the functions of other members of this gene family and found eight new drivers as well as three new suppressors of drive. Surprisingly, we find that drive is relatively insensitive to changes in wtf gene sequence as highly diverged wtf genes execute gamete killing with similar efficiency. Finally, we also find that the expression and localization of some Wtf proteins are distinct from those of known drivers and suppressors, suggesting that these proteins may have non-meiotic drive functions.