Bulk Segregant Analysis Reveals the Genetic Basis of a Natural Trait Variation in Fission Yeast

Engineering heterothallic strains in fission yeast

In poor nitrogen conditions, fission yeast cells mate, undergo meiosis and form spores that are resistant to deleterious environments. Natural isolates of Schizosaccharomyces pombe are homothallic. This allows them to naturally switch between the two h- and h+ mating types with a high frequency, thereby ensuring the presence of both mating partners in a population of cells. However, alteration of the mating type locus can abolish mating type switching or reduce it to a very low frequency. Such heterothallic strains have been isolated and are common in research laboratories due to the simplicity of their use for Mendelian genetics. In addition to the standard laboratory strains, a large collection of natural S. pombe isolates is now available, representing a powerful resource for investigating the genetic diversity and biology of fission yeast. However, most of these strains are homothallic, and only tedious or mutagenic strategies have been described to obtain heterothallic cells from a homothallic parent. Here, we describe a simple approach to generate heterothallic strains. It takes advantage of an alteration of the mating type locus that was previously identified in a mating type switching-deficient strain and the CRISPR-Cas9 editing tool, allowing for a one-step engineering of heterothallic cells with high efficiency.

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.

The ortholog of human REEP1-4 is required for autophagosomal enclosure of ER-phagy/nucleophagy cargos in fission yeast

Selective macroautophagy of the endoplasmic reticulum (ER) and the nucleus, known as ER-phagy and nucleophagy, respectively, are processes whose mechanisms remain inadequately understood. Through an imaging-based screen, we find that in the fission yeast Schizosaccharomyces pombe, Yep1 (also known as Hva22 or Rop1), the ortholog of human REEP1-4, is essential for ER-phagy and nucleophagy but not for bulk autophagy. In the absence of Yep1, the initial phase of ER-phagy and nucleophagy proceeds normally, with the ER-phagy/nucleophagy receptor Epr1 coassembling with Atg8. However, ER-phagy/nucleophagy cargos fail to reach the vacuole. Instead, nucleus- and cortical-ER-derived membrane structures not enclosed within autophagosomes accumulate in the cytoplasm. Intriguingly, the outer membranes of nucleus-derived structures remain continuous with the nuclear envelope-ER network, suggesting a possible outer membrane fission defect during cargo separation from source compartments. We find that the ER-phagy role of Yep1 relies on its abilities to self-interact and shape membranes and requires its C-terminal amphipathic helices. Moreover, we show that human REEP1-4 and budding yeast Atg40 can functionally substitute for Yep1 in ER-phagy, and Atg40 is a divergent ortholog of Yep1 and REEP1-4. Our findings uncover an unexpected mechanism governing the autophagosomal enclosure of ER-phagy/nucleophagy cargos and shed new light on the functions and evolution of REEP family proteins.

S. pombe wtf drivers use dual transcriptional regulation and selective protein exclusion from spores to cause meiotic drive

Meiotic drivers bias gametogenesis to ensure their transmission into more than half the offspring of a heterozygote. In Schizosaccharomyces pombe, wtf meiotic drivers destroy the meiotic products (spores) that do not inherit the driver from a heterozygote, thereby reducing fertility. wtf drivers encode both a Wtfpoison protein and a Wtfantidote protein using alternative transcriptional start sites. Here, we analyze how the expression and localization of the Wtf proteins are regulated to achieve drive. We show that transcriptional timing and selective protein exclusion from developing spores ensure that all spores are exposed to Wtf4poison, but only the spores that inherit wtf4 receive a dose of Wtf4antidote sufficient for survival. In addition, we show that the Mei4 transcription factor, a master regulator of meiosis, controls the expression of the wtf4poison transcript. This transcriptional regulation, which includes the use of a critical meiotic transcription factor, likely complicates the universal suppression of wtf genes without concomitantly disrupting spore viability. We propose that these features contribute to the evolutionary success of the wtf drivers.

Harnessing the potential of bulk segregant analysis sequencing and its related approaches in crop breeding

Most plant traits are governed by polygenes including both major and minor genes. Linkage mapping and positional cloning have contributed greatly to mapping genomic loci controlling important traits in crop species. However, they are low-throughput, time-consuming, and have low resolution due to which their efficiency in crop breeding is reduced. In this regard, the bulk segregant analysis sequencing (BSA-seq) and its related approaches, viz., quantitative trait locus (QTL)-seq, bulk segregant RNA-Seq (BSR)-seq, and MutMap, have emerged as efficient methods to identify the genomic loci/QTLs controlling specific traits at high resolution, accuracy, reduced time span, and in a high-throughput manner. These approaches combine BSA with next-generation sequencing (NGS) and enable the rapid identification of genetic loci for qualitative and quantitative assessments. Many previous studies have shown the successful identification of the genetic loci for different plant traits using BSA-seq and its related approaches, as discussed in the text with details. However, the efficiency and accuracy of the BSA-seq depend upon factors like sequencing depth and coverage, which enhance the sequencing cost. Recently, the rapid reduction in the cost of NGS together with the expected cost reduction of third-generation sequencing in the future has further increased the accuracy and commercial applicability of these approaches in crop improvement programs. This review article provides an overview of BSA-seq and its related approaches in crop breeding together with their merits and challenges in trait mapping.

Two independent approaches converge to the cloning of a new Leptosphaeria maculans avirulence effector gene, AvrLmS-Lep2

Brassica napus (oilseed rape, canola) seedling resistance to Leptosphaeria maculans, the causal agent of blackleg (stem canker) disease, follows a gene-for-gene relationship. The avirulence genes AvrLmS and AvrLep2 were described to be perceived by the resistance genes RlmS and LepR2, respectively, present in B. napus 'Surpass 400'. Here we report cloning of AvrLmS and AvrLep2 using two independent methods. AvrLmS was cloned using combined in vitro crossing between avirulent and virulent isolates with sequencing of DNA bulks from avirulent or virulent progeny (bulked segregant sequencing). AvrLep2 was cloned using a biparental cross of avirulent and virulent L. maculans isolates and a classical map-based cloning approach. Taking these two approaches independently, we found that AvrLmS and AvrLep2 are the same gene. Complementation of virulent isolates with this gene confirmed its role in inducing resistance on Surpass 400, Topas-LepR2, and an RlmS-line. The gene, renamed AvrLmS-Lep2, encodes a small cysteine-rich protein of unknown function with an N-terminal secretory signal peptide, which is a common feature of the majority of effectors from extracellular fungal plant pathogens. The AvrLmS-Lep2/LepR2 interaction phenotype was found to vary from a typical hypersensitive response through intermediate resistance sometimes towards susceptibility, depending on the inoculation conditions. AvrLmS-Lep2 was nevertheless sufficient to significantly slow the systemic growth of the pathogen and reduce the stem lesion size on plant genotypes with LepR2, indicating the potential efficiency of this resistance to control the disease in the field.

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.

R-loops and regulatory changes in chronologically ageing fission yeast cells drive non-random patterns of genome rearrangements

Aberrant repair of DNA double-strand breaks can recombine distant chromosomal breakpoints. Chromosomal rearrangements compromise genome function and are a hallmark of ageing. Rearrangements are challenging to detect in non-dividing cell populations, because they reflect individually rare, heterogeneous events. The genomic distribution of de novo rearrangements in non-dividing cells, and their dynamics during ageing, remain therefore poorly characterized. Studies of genomic instability during ageing have focussed on mitochondrial DNA, small genetic variants, or proliferating cells. To characterize genome rearrangements during cellular ageing in non-dividing cells, we interrogated a single diagnostic measure, DNA breakpoint junctions, using Schizosaccharomyces pombe as a model system. Aberrant DNA junctions that accumulated with age were associated with microhomology sequences and R-loops. Global hotspots for age-associated breakpoint formation were evident near telomeric genes and linked to remote breakpoints elsewhere in the genome, including the mitochondrial chromosome. Formation of breakpoint junctions at global hotspots was inhibited by the Sir2 histone deacetylase and might be triggered by an age-dependent de-repression of chromatin silencing. An unexpected mechanism of genomic instability may cause more local hotspots: age-associated reduction in an RNA-binding protein triggering R-loops at target loci. This result suggests that biological processes other than transcription or replication can drive genome rearrangements. Notably, we detected similar signatures of genome rearrangements that accumulated in old brain cells of humans. These findings provide insights into the unique patterns and possible mechanisms of genome rearrangements in non-dividing cells, which can be promoted by ageing-related changes in gene-regulatory proteins.

Recombination and biased segregation of mitochondrial genomes during crossing and meiosis of different Schizosaccharomyces pombe strains

During meiosis, tethering of parental mitochondria to opposite cell poles inhibits the mixing of mitochondria with different genomes and ensures uniparental inheritance in thestandard laboratory strain of fission yeast. We here investigate mitochondrial inheritance in crosses between natural isolates using tetrad dissection and next-generation sequencing. We find that colonies grown from single spores can sometimes carry a mix of mitochondrial genotypes, that mitochondrial genomes can recombine during meiosis, that in some cases tetrads do not follow the 2:2 segregation pattern, and that certain crosses may feature a weak bias towards one of the parents. Together, these findings paint a more nuanced picture of mitochondrial inheritance in the wild.

Proteomic analysis of meiosis and characterization of novel short open reading frames in the fission yeast Schizosaccharomyces pombe

Meiosis is the process by which haploid gametes are produced from diploid precursor cells. We used stable isotope labeling by amino acids in cell culture (SILAC) to characterize the meiotic proteome in the fission yeast Schizosaccharomyces pombe. We compared relative levels of proteins extracted from cells harvested around meiosis I with those of meiosis II, and proteins from premeiotic S phase with the interval between meiotic divisions, when S phase is absent. Our proteome datasets revealed peptides corresponding to short open reading frames (sORFs) that have been previously identified by ribosome profiling as new translated regions. We verified expression of selected sORFs by Western blotting and analyzed the phenotype of deletion mutants. Our data provide a resource for studying meiosis that may help understand differences between meiosis I and meiosis II and how S phase is suppressed between the two meiotic divisions.

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.

Killer Meiotic Drive and Dynamic Evolution of the wtf Gene Family

Natural selection works best when the two alleles in a diploid organism are transmitted to offspring at equal frequencies. Despite this, selfish loci known as meiotic drivers that bias their own transmission into gametes are found throughout eukaryotes. Drive is thought to be a powerful evolutionary force, but empirical evolutionary analyses of drive systems are limited by low numbers of identified meiotic drive genes. Here, we analyze the evolution of the wtf gene family of Schizosaccharomyces pombe that contains both killer meiotic drive genes and suppressors of drive. We completed assemblies of all wtf genes for two S. pombe isolates, as well as a subset of wtf genes from over 50 isolates. We find that wtf copy number can vary greatly between isolates and that amino acid substitutions, expansions and contractions of DNA sequence repeats, and nonallelic gene conversion between family members all contribute to dynamic wtf gene evolution. This work demonstrates the power of meiotic drive to foster rapid evolution and identifies a recombination mechanism through which transposons can indirectly mobilize meiotic drivers.

Estimation of QTL heritability based on pooled sequencing data

Motivation

Bulked segregant analysis combined with next generation sequencing has proven to be a simple and efficient approach for fast mapping of quantitative trait loci (QTLs). However, how to estimate the proportion of phenotypic variance explained by a QTL (or termed QTL heritability) in such pooled QTL mapping is an unsolved problem.

Results

In this paper, we propose a method called PQHE to estimate QTL heritability using pooled sequencing data obtained under different experimental designs. Simulation studies indicated that our method is correct and feasible. Four practical examples from rice and yeast are demonstrated, each representing a different situation.

Availability and implementation

The R scripts of our method are open source under GPLv3 license at http://genetics.fafu.edu.cn/PQHE or https://github.com/biotangweiqi/PQHE. The R scripts require the R package rootSolve.

Contact

wuwr@fafu.edu.cn.

Supplementary information

Supplementary data are available at Bioinformatics online.

Uncovering Natural Longevity Alleles from Intercrossed Pools of Aging Fission Yeast Cells

Chronological lifespan of non-dividing yeast cells is a quantitative trait that reflects cellular aging. By monitoring allele frequencies in aging segregant pools, Ellis et al. uncover regulatory variants in the 5'-untranslated regions of two genes...

Quantitative traits often show large variation caused by multiple genetic factors . One such trait is the chronological lifespan of non-dividing yeast cells, serving as a model for cellular aging. Screens for genetic factors involved in aging typically assay mutants of protein-coding genes. To identify natural genetic variants contributing to cellular aging, we exploited two strains of the fission yeast, Schizosaccharomyces pombe, that differ in chronological lifespan. We generated segregant pools from these strains and subjected them to advanced intercrossing over multiple generations to break up linkage groups. We chronologically aged the intercrossed segregant pool, followed by genome sequencing at different times to detect genetic variants that became reproducibly enriched as a function of age. A region on Chromosome II showed strong positive selection during aging. Based on expected functions, two candidate variants from this region in the long-lived strain were most promising to be causal: small insertions and deletions in the 5′-untranslated regions of ppk31 and SPBC409.08. Ppk31 is an ortholog of Rim15, a conserved kinase controlling cell proliferation in response to nutrients, while SPBC409.08 is a predicted spermine transmembrane transporter. Both Rim15 and the spermine-precursor, spermidine, are implicated in aging as they are involved in autophagy-dependent lifespan extension. Single and double allele replacement suggests that both variants, alone or combined, have subtle effects on cellular longevity. Furthermore, deletion mutants of both ppk31 and SPBC409.08 rescued growth defects caused by spermidine. We propose that Ppk31 and SPBC409.08 may function together to modulate lifespan, thus linking Rim15/Ppk31 with spermidine metabolism.

A Cloning-Free Method for CRISPR/Cas9-Mediated Genome Editing in Fission Yeast

The CRISPR/Cas9 system, which relies on RNA‐guided DNA cleavage to induce site-specific DNA double-strand breaks, is a powerful tool for genome editing. This system has been successfully adapted for the fission yeast Schizosaccharomyces pombe by expressing Cas9 and the single-guide RNA (sgRNA) from a plasmid. In the procedures published to date, the cloning step that introduces a specific sgRNA target sequence into the plasmid is the most tedious and time-consuming. To increase the efficiency of applying the CRISPR/Cas9 system in fission yeast, we here developed a cloning-free procedure that uses gap repair in fission yeast cells to assemble two linear DNA fragments, a gapped Cas9-encoding plasmid and a PCR-amplified sgRNA insert, into a circular plasmid. Both fragments contain only a portion of the ura4 or bsdMX marker so that only the correctly assembled plasmid can confer uracil prototrophy or blasticidin resistance. We show that this gap-repair-based and cloning-free CRISPR/Cas9 procedure permits rapid and efficient point mutation knock-in, endogenous N-terminal tagging, and genomic sequence deletion in fission yeast.

Veni, vidi, vici: the success of wtf meiotic drivers in fission yeast

Meiotic drivers are selfish DNA loci that can bias their own transmission into gametes. Owing to their transmission advantages, meiotic drivers can spread in populations even if the drivers or linked variants decrease organismal fitness. Meiotic drive was first formally described in the 1950s and is thought to be a powerful force shaping eukaryotic genomes. Classic genetic analyses have detected the action of meiotic drivers in plants, filamentous fungi, insects and vertebrates. Several of these drive systems have limited experimental tractability and relatively little is known about the molecular mechanisms of meiotic drive. Recently, however, meiotic drivers were discovered in a yeast species. The Schizosaccharomyces pombe wtf gene family contains several active meiotic drive genes. This review summarizes what is known about the wtf family and highlights its potential as a highly tractable experimental model for molecular and evolutionary characterization of meiotic drive.

The natural diversity and ecology of fission yeast

While the fission yeast is a powerful model of eukaryote biology, there have been few studies of quantitative genetics, phenotypic or genetic diversity. Here I survey the small collection of fission yeast diversity research. I discuss what we can infer about the ecology and origins of Schizosaccharomyces pombe from microbiology field studies and the few strains that have been collected.

A large gene family in fission yeast encodes spore killers that subvert Mendel's law

Spore killers in fungi are selfish genetic elements that distort Mendelian segregation in their favor. It remains unclear how many species harbor them and how diverse their mechanisms are. Here, we discover two spore killers from a natural isolate of the fission yeast Schizosaccharomyces pombe. Both killers belong to the previously uncharacterized wtf gene family with 25 members in the reference genome. These two killers act in strain-background-independent and genome-location-independent manners to perturb the maturation of spores not inheriting them. Spores carrying one killer are protected from its killing effect but not that of the other killer. The killing and protecting activities can be uncoupled by mutation. The numbers and sequences of wtf genes vary considerably between S. pombe isolates, indicating rapid divergence. We propose that wtf genes contribute to the extensive intraspecific reproductive isolation in S. pombe, and represent ideal models for understanding how segregation-distorting elements act and evolve.

The repeatability of genome-wide mutation rate and spectrum estimates

Over the last decade, mutation studies have grown in popularity due to the affordability and accessibility of whole genome sequencing. As the number of species in which spontaneous mutation has been directly estimated approaches 20 across two domains of life, questions arise over the repeatability of results in such experiments. Five species were identified in which duplicate mutation studies have been performed. Across these studies the difference in estimated spontaneous mutation rate is at most, weakly significant (p < 0.01). However, a highly significant (p < 10(-5)), threefold difference in the rate of insertions/deletions (indels) exists between two recent studies in Schizosaccharomyces pombe. Upon investigation of the ancestral genome sequence for both studies, a possible anti-mutator allele was identified. The observed variation in indel rate may imply that the use of indel markers, such as microsatellites, for the investigation of genetic diversity within and among populations may be inappropriate because of the assumption of uniform mutation rate within a species.

Experimental evolution in fungi: An untapped resource

Historically, evolutionary biology has been considered an observational science. Examining populations and inferring evolutionary histories mold evolutionary theories. In contrast, laboratory evolution experiments make use of the amenability of traditional model organisms to study fundamental processes underlying evolution in real time in simple, but well-controlled, environments. With advances in high-throughput biology and next generation sequencing, it is now possible to propagate hundreds of parallel populations over thousands of generations and to quantify precisely the frequencies of various mutations over time. Experimental evolution combines the ability to simultaneously monitor replicate populations with the power to vary individual parameters to test specific evolutionary hypotheses, something that is impractical or infeasible in natural populations. Many labs are now conducting laboratory evolution experiments in nearly all model systems including viruses, bacteria, yeast, nematodes, and fruit flies. Among these systems, fungi occupy a unique niche: with a short generation time, small compact genomes, and sexual cycles, fungi are a particularly valuable and largely untapped resource for propelling future growth in the field of experimental evolution. Here, we describe the current state of fungal experimental evolution and why fungi are uniquely positioned to answer many of the outstanding questions in the field. We also review which fungal species are most well suited for experimental evolution.