Mapping epigenetic mutations in fission yeast using whole-genome next-generation sequencing

The TAC1 Gene in Candida albicans: Structure, Function, and Role in Azole Resistance: A Mini-Review

Candidiasis is a common fungal infection caused by Candida species, with Candida albicans being the most prevalent. Resistance to azole drugs, commonly used to treat Candida infections, poses a significant challenge. Transcriptional activator candidate 1 (TAC1) gene has emerged as a key player in regulating drug resistance in C. albicans. This review explores the structure and function of the TAC1 gene and its role in azole resistance. This gene encodes a transcription factor that controls the expression of genes involved in drug resistance, such as efflux pump genes (CDR1, CDR2, and MDR1) and ERG11. Mutations in TAC1 can increase these genes' expression and confer resistance to azoles. Various TAC1 gene mutations, mostly gain-of-function mutations, have been identified, which upregulate CDR1 and CDR2 expression, resulting in azole resistance. Understanding the mechanisms of azole resistance mediated by the TAC1 gene is crucial for the strategies in the effective antifungal development pipeline.

The leucine-NH4+ uptake regulator Any1 limits growth as part of a general amino acid control response to loss of La protein by fission yeast

The sla1⁺ gene of Schizosachharoymces pombe encodes La protein which promotes proper processing of precursor-tRNAs. Deletion of sla1 (sla1Δ) leads to disrupted tRNA processing and sensitivity to target of rapamycin (TOR) inhibition. Consistent with this, media containing NH4⁺ inhibits leucine uptake and growth of sla1Δ cells. Here, transcriptome analysis reveals that genes upregulated in sla1Δ cells exhibit highly significant overalp with general amino acid control (GAAC) genes in relevant transcriptomes from other studies. Growth in NH4⁺ media leads to additional induced genes that are part of a core environmental stress response (CESR). The sla1Δ GAAC response adds to evidence linking tRNA homeostasis and broad signaling in S. pombe. We provide evidence that deletion of the Rrp6 subunit of the nuclear exosome selectively dampens a subset of GAAC genes in sla1Δ cells suggesting that nuclear surveillance-mediated signaling occurs in S. pombe. To study the NH4⁺-effects, we isolated sla1Δ spontaneous revertants (SSR) of the slow growth phenotype and found that GAAC gene expression and rapamycin hypersensitivity were also reversed. Genome sequencing identified a F32V substitution in Any1, a known negative regulator of NH4⁺-sensitive leucine uptake linked to TOR. We show that ³H-leucine uptake by SSR-any1-F32V cells in NH4⁺-media is more robust than by sla1Δ cells. Moreover, F32V may alter any1⁺ function in sla1Δ vs. sla1⁺ cells in a distinctive way. Thus deletion of La, a tRNA processing factor leads to a GAAC response involving reprogramming of amino acid metabolism, and isolation of the any1-F32V rescuing mutant provides an additional specific link.

Multiple nutritional phenotypes of fission yeast mutants defective in genes encoding essential mitochondrial proteins

Mitochondria are essential for regulation of cellular respiration, energy production, small molecule metabolism, anti-oxidation and cell ageing, among other things. While the mitochondrial genome contains a small number of protein-coding genes, the great majority of mitochondrial proteins are encoded by chromosomal genes. In the fission yeast Schizosaccharomyces pombe, 770 proteins encoded by chromosomal genes are located in mitochondria. Of these, 195 proteins, many of which are implicated in translation and transport, are absolutely essential for viability. We isolated and characterized eight temperature-sensitive (ts) strains with mutations in essential mitochondrial proteins. Interestingly, they are also sensitive to limited nutrition (glucose and/or nitrogen), producing low-glucose-sensitive and ‘super-housekeeping' phenotypes. They fail to produce colonies under low-glucose conditions at the permissive temperature or lose cell viability under nitrogen starvation at the restrictive temperature. The majority of these ts mitochondrial mutations may cause defects of gene expression in the mitochondrial genome. mrp4 and mrp17 are defective in mitochondrial ribosomal proteins. ppr3 is defective in rRNA expression, and trz2 and vrs2 are defective in tRNA maturation. This study promises potentially large dividends because mitochondrial quiescent functions are vital for human brain and muscle, and also for longevity.

MutantHuntWGS: A Pipeline for Identifying Saccharomyces cerevisiae Mutations

MutantHuntWGS is a user-friendly pipeline for analyzing Saccharomyces cerevisiae whole-genome sequencing data. It uses available open-source programs to: (1) perform sequence alignments for paired and single-end reads, (2) call variants, and (3) predict variant effect and severity. MutantHuntWGS outputs a shortlist of variants while also enabling access to all intermediate files. To demonstrate its utility, we use MutantHuntWGS to assess multiple published datasets; in all cases, it detects the same causal variants reported in the literature. To encourage broad adoption and promote reproducibility, we distribute a containerized version of the MutantHuntWGS pipeline that allows users to install and analyze data with only two commands. The MutantHuntWGS software and documentation can be downloaded free of charge from https://github.com/mae92/MutantHuntWGS.

Nitrogen starvation reveals the mitotic potential of mutants in the S/MAPK pathways

The genetics of quiescence is an emerging field compared to that of growth, yet both states generate spontaneous mutations and genetic diversity fueling evolution. Reconciling mutation rates in dividing conditions and mutation accumulation as a function of time in non-dividing situations remains a challenge. Nitrogen-starved fission yeast cells reversibly arrest proliferation, are metabolically active and highly resistant to a variety of stresses. Here, we show that mutations in stress- and mitogen-activated protein kinase (S/MAPK) signaling pathways are enriched in aging cultures. Targeted resequencing and competition experiments indicate that these mutants arise in the first month of quiescence and expand clonally during the second month at the expense of the parental population. Reconstitution experiments show that S/MAPK modules mediate the sacrifice of many cells for the benefit of some mutants. These findings suggest that non-dividing conditions promote genetic diversity to generate a social cellular environment prone to kin selection.

Nitrogen-starved fission yeast cells survive for weeks without dividing. Here, the authors show that some of these surviving cells accumulate mutations in the stress- and mitogen-activated protein kinase pathways and outcompete their parental cells, which provide nutrients for the mutant cells.

Combining next-generation sequencing and progeny testing for rapid identification of induced recessive and dominant mutations in maize M2 individuals

Molecular identification of mutant alleles responsible for certain phenotypic alterations is a central goal of genetic analyses. In this study we describe a rapid procedure suitable for the identification of induced recessive and dominant mutations applied to two Zea mays mutants expressing a dwarf and a pale green phenotype, respectively, which were obtained through pollen ethyl methanesulfonate (EMS) mutagenesis. First, without prior backcrossing, induced mutations (single nucleotide polymorphisms, SNPs) segregating in a (M₂ ) family derived from a heterozygous (M₁ ) parent were identified using whole-genome shotgun (WGS) sequencing of a small number of (M₂ ) individuals with mutant and wild-type phenotypes. Second, the state of zygosity of the mutation causing the phenotype was determined for each sequenced individual by phenotypic segregation analysis of the self-pollinated (M₃ ) offspring. Finally, we filtered for segregating EMS-induced SNPs whose state of zygosity matched the determined state of zygosity of the mutant locus in each sequenced (M₂ ) individuals. Through this procedure, combining sequencing of individuals and Mendelian inheritance, three and four SNPs in linkage passed our zygosity filter for the homozygous dwarf and heterozygous pale green mutation, respectively. The dwarf mutation was found to be allelic to the an1 locus and caused by an insertion in the largest exon of the AN1 gene. The pale green mutation affected the nuclear W2 gene and was caused by a non-synonymous amino acid exchange in encoded chloroplast DNA polymerase with a predicted deleterious effect. This coincided with lower cpDNA levels in pale green plants.

Whole-Genome Sequencing Identifies a Rice Grain Shape Mutant, gs9-1

BACKGROUND

Breeding for genes controlling key agronomic traits is an important goal of rice genetic improvement. To gain insight into genes controlling grain morphology, we screened M3 plants derived from 1,000 whole-genome sequenced (WGS) M2 Kitaake mutants to identify lines with altered grain size.

RESULTS

In this study, we isolated a mutant, named fast-neutron (FN) 60-4, which exhibits a significant reduction in grain size. We crossed FN60-4 with the parental line Kitaake and analyzed the resulting backcross population. Segregation analysis of 113 lines from the BC2F2 population revealed that the mutant phenotype is controlled by a single semi-dominant locus. Mutant FN60-4 is reduced 20% in plant height and 8.8% in 1000-grain weight compared with Kitaake. FN60-4 also exhibits an 8% reduction in cell number and a 9% reduction in cell length along the vertical axis of the glume. We carried out whole-genome sequencing of DNA pools extracted from segregants with long grains or short grains, and revealed that one gene, LOC_Os09g02650, cosegregated with the grain size phenotype in the BC1F2 and BC2F2 populations. This mutant allele was named grain shape 9-1 (gs9-1). gs9-1 carries a 3-bp deletion that affects two amino acids. This locus is a new allele of the BC12/GDD1/MTD1 gene that encodes a kinesin-like protein involved in cell-cycle progression, cellulose microfibril deposition and gibberellic acid (GA) biosynthesis. The GA biosynthesis-related gene KO2 is down-regulated in gs9-1. The dwarf phenotype of gs9-1 can be rescued by adding exogenous GA3. In contrast to the phenotypes for the other alleles, the gs9-1 is less severe, consistent with the nature of the mutation, which does not disrupt the open reading frame as observed for the other alleles.

CONCLUSIONS

In this study, we isolated a mutant, which exhibits altered grain shape and identified the mutated gene, gs9-1. Our study reveals that gs9-1 is a semi-dominant gene that carries a two-amino acid mutation. gs9-1 is allelic to the BC12/GDD1/MTD1 gene involved in GA biosynthesis. These results demonstrate the efficiency and convenience of cloning genes from the whole-genome sequenced Kitaake mutant population to advance investigations into genes controlling key agronomic traits in rice.

Dependency of Heterochromatin Domains on Replication Factors

Chromatin structure regulates both genome expression and dynamics in eukaryotes, where large heterochromatic regions are epigenetically silenced through the methylation of histone H3K9, histone deacetylation, and the assembly of repressive complexes. Previous genetic screens with the fission yeast Schizosaccharomyces pombe have led to the identification of key enzymatic activities and structural constituents of heterochromatin. We report here on additional factors discovered by screening a library of deletion mutants for silencing defects at the edge of a heterochromatic domain bound by its natural boundary-the IR-R+ element-or by ectopic boundaries. We found that several components of the DNA replication progression complex (RPC), including Mrc1/Claspin, Mcl1/Ctf4, Swi1/Timeless, Swi3/Tipin, and the FACT subunit Pob3, are essential for robust heterochromatic silencing, as are the ubiquitin ligase components Pof3 and Def1, which have been implicated in the removal of stalled DNA and RNA polymerases from chromatin. Moreover, the search identified the cohesin release factor Wpl1 and the forkhead protein Fkh2, both likely to function through genome organization, the Ssz1 chaperone, the Fkbp39 proline cis-trans isomerase, which acts on histone H3P30 and P38 in Saccharomyces cerevisiae, and the chromatin remodeler Fft3. In addition to their effects in the mating-type region, to varying extents, these factors take part in heterochromatic silencing in pericentromeric regions and telomeres, revealing for many a general effect in heterochromatin. This list of factors provides precious new clues with which to study the spatiotemporal organization and dynamics of heterochromatic regions in connection with DNA replication.

Revitalization of a Forward Genetic Screen Identifies Three New Regulators of Fungal Secondary Metabolism in the Genus Aspergillus

The study of aflatoxin in Aspergillus spp. has garnered the attention of many researchers due to aflatoxin's carcinogenic properties and frequency as a food and feed contaminant. Significant progress has been made by utilizing the model organism Aspergillus nidulans to characterize the regulation of sterigmatocystin (ST), the penultimate precursor of aflatoxin. A previous forward genetic screen identified 23 A. nidulans mutants involved in regulating ST production. Six mutants were characterized from this screen using classical mapping (five mutations in mcsA) and complementation with a cosmid library (one mutation in laeA). The remaining mutants were backcrossed and sequenced using Illumina and Ion Torrent sequencing platforms. All but one mutant contained one or more sequence variants in predicted open reading frames. Deletion of these genes resulted in identification of mutant alleles responsible for the loss of ST production in 12 of the 17 remaining mutants. Eight of these mutations were in genes already known to affect ST synthesis (laeA, mcsA, fluG, and stcA), while the remaining four mutations (in laeB, sntB, and hamI) were in previously uncharacterized genes not known to be involved in ST production. Deletion of laeB, sntB, and hamI in A. flavus results in loss of aflatoxin production, confirming that these regulators are conserved in the aflatoxigenic aspergilli. This report highlights the multifaceted regulatory mechanisms governing secondary metabolism in Aspergillus Additionally, these data contribute to the increasing number of studies showing that forward genetic screens of fungi coupled with whole-genome resequencing is a robust and cost-effective technique.IMPORTANCE In a postgenomic world, reverse genetic approaches have displaced their forward genetic counterparts. The techniques used in forward genetics to identify loci of interest were typically very cumbersome and time-consuming, relying on Mendelian traits in model organisms. The current work was pursued not only to identify alleles involved in regulation of secondary metabolism but also to demonstrate a return to forward genetics to track phenotypes and to discover genetic pathways that could not be predicted through a reverse genetics approach. While identification of mutant alleles from whole-genome sequencing has been done before, here we illustrate the possibility of coupling this strategy with a genetic screen to identify multiple alleles of interest. Sequencing of classically derived mutants revealed several uncharacterized genes, which represent novel pathways to regulate and control the biosynthesis of sterigmatocystin and of aflatoxin, a societally and medically important mycotoxin.

Genetic Analysis of Schizosaccharomyces pombe

In this introduction we discuss some basic genetic tools and techniques that are used with the fission yeast Schizosaccharomyces pombe Genes commonly used for selection or as reporters are discussed, with an emphasis on genes that permit counterselection, intragenic complementation, or colony-color assays. S. pombe is most stable as a haploid organism. We describe its mating-type system, how to perform genetic crosses and methods for selecting and propagating diploids. We discuss the relative merits of tetrad dissection and random spore preparation in strain construction and genetic analyses. Finally, we present several types of mutant screens, with an evaluation of their respective strengths and limitations in the light of emerging technologies such as next-generation sequencing.

RNA interference is essential for cellular quiescence

Quiescent cells play a predominant role in most organisms. Here we identify RNA interference (RNAi) as a major requirement for quiescence (G₀ phase of the cell cycle) in Schizosaccharomyces pombe RNAi mutants lose viability at G₀ entry and are unable to maintain long-term quiescence. We identified suppressors of G₀ defects in cells lacking Dicer (dcr1Δ), which mapped to genes involved in chromosome segregation, RNA polymerase-associated factors, and heterochromatin formation. We propose a model in which RNAi promotes the release of RNA polymerase in cycling and quiescent cells: (i) RNA polymerase II release mediates heterochromatin formation at centromeres, allowing proper chromosome segregation during mitotic growth and G₀ entry, and (ii) RNA polymerase I release prevents heterochromatin formation at ribosomal DNA during quiescence maintenance. Our model may account for the codependency of RNAi and histone H3 lysine 9 methylation throughout eukaryotic evolution.

Mapping Challenging Mutations by Whole-Genome Sequencing

Whole-genome sequencing provides a rapid and powerful method for identifying mutations on a global scale, and has spurred a renewed enthusiasm for classical genetic screens in model organisms. The most commonly characterized category of mutation consists of monogenic, recessive traits, due to their genetic tractability. Therefore, most of the mapping methods for mutation identification by whole-genome sequencing are directed toward alleles that fulfill those criteria (i.e., single-gene, homozygous variants). However, such approaches are not entirely suitable for the characterization of a variety of more challenging mutations, such as dominant and semidominant alleles or multigenic traits. Therefore, we have developed strategies for the identification of those classes of mutations, using polymorphism mapping in Caenorhabditis elegans as our model for validation. We also report an alternative approach for mutation identification from traditional recombinant crosses, and a solution to the technical challenge of sequencing sterile or terminally arrested strains where population size is limiting. The methods described herein extend the applicability of whole-genome sequencing to a broader spectrum of mutations, including classes that are difficult to map by traditional means.

Genome sequencing of the Trichoderma reesei QM9136 mutant identifies a truncation of the transcriptional regulator XYR1 as the cause for its cellulase-negative phenotype

Background

Trichoderma reesei is the main industrial source of cellulases and hemicellulases required for the hydrolysis of biomass to simple sugars, which can then be used in the production of biofuels and biorefineries. The highly productive strains in use today were generated by classical mutagenesis. As byproducts of this procedure, mutants were generated that turned out to be unable to produce cellulases. In order to identify the mutations responsible for this inability, we sequenced the genome of one of these strains, QM9136, and compared it to that of its progenitor T. reesei QM6a.

Results

In QM9136, we detected a surprisingly low number of mutagenic events in the promoter and coding regions of genes, i.e. only eight indels and six single nucleotide variants. One of these indels led to a frame-shift in the Zn₂Cys₆ transcription factor XYR1, the general regulator of cellulase and xylanase expression, and resulted in its C-terminal truncation by 140 amino acids. Retransformation of strain QM9136 with the wild-type xyr1 allele fully recovered the ability to produce cellulases, and is thus the reason for the cellulase-negative phenotype. Introduction of an engineered xyr1 allele containing the truncating point mutation into the moderate producer T. reesei QM9414 rendered this strain also cellulase-negative. The correspondingly truncated XYR1 protein was still able to enter the nucleus, but failed to be expressed over the basal constitutive level.

Conclusion

The missing 140 C-terminal amino acids of XYR1 are therefore responsible for its previously observed auto-regulation which is essential for cellulases to be expressed. Our data present a working example of the use of genome sequencing leading to a functional explanation of the QM9136 cellulase-negative phenotype.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1526-0) contains supplementary material, which is available to authorized users.

Fast forward genetics to identify mutations causing a high light tolerant phenotype in Chlamydomonas reinhardtii by whole-genome-sequencing

BACKGROUND

High light tolerance of microalgae is a desired phenotype for efficient cultivation in large scale production systems under fluctuating outdoor conditions. Outdoor cultivation requires the use of either wild-type or non-GMO derived mutant strains due to safety concerns. The identification and molecular characterization of such mutants derived from untagged forward genetics approaches was limited previously by the tedious and time-consuming methods involving techniques such as classical meiotic mapping. The combination of mapping with next generation sequencing technologies offers alternative strategies to identify genes involved in high light adaptation in untagged mutants.

RESULTS

We used the model alga Chlamydomonas reinhardtii in a non-GMO mutation strategy without any preceding crossing step or pooled progeny to identify genes involved in the regulatory processes of high light adaptation. To generate high light tolerant mutants, wildtype cells were mutagenized only to a low extent, followed by a stringent selection. We performed whole-genome sequencing of two independent mutants hit1 and hit2 and the parental wildtype. The availability of a reference genome sequence and the removal of shared bakground variants between the wildtype strain and each mutant, enabled us to identify two single nucleotide polymorphisms within the same gene Cre02.g085050, hereafter called LRS1 (putative Light Response Signaling protein 1). These two independent single amino acid exchanges are both located in the putative WD40 propeller domain of the corresponding protein LRS1. Both mutants exhibited an increased rate of non-photochemical-quenching (NPQ) and an improved resistance against chemically induced reactive oxygen species. In silico analyses revealed homology of LRS1 to the photoregulatory protein COP1 in plants.

CONCLUSIONS

In this work we identified the nuclear encoded gene LRS1 as an essential factor for high light adaptation in C. reinhardtii. The causative random mutation within this gene was identified by a rapid and efficient method, avoiding any preceding crossing step, meiotic mapping, or pooled progeny. Our results open up new insights into mechanisms of high light adaptation in microalgae and at the same time provide a simplified strategy for non-GMO forward genetics, a crucial precondition that could result in the identification of key factors for economically relevant biological processes within algae.

Dicer promotes transcription termination at sites of replication stress to maintain genome stability

Nuclear RNAi is an important regulator of transcription and epigenetic modification, but the underlying mechanisms remain elusive. Using a genome-wide approach in the fission yeast S. pombe, we have found that Dcr1, but not other components of the canonical RNAi pathway, promotes the release of Pol II from the 3? end of highly transcribed genes, and, surprisingly, from antisense transcription of rRNA and tRNA genes, which are normally transcribed by Pol I and Pol III. These Dcr1-terminated loci correspond to sites of replication stress and DNA damage, likely resulting from transcription-replication collisions. At the rDNA loci, release of Pol II facilitates DNA replication and prevents homologous recombination, which would otherwise lead to loss of rDNA repeats especially during meiosis. Our results reveal a novel role for Dcr1-mediated transcription termination in genome maintenance and may account for widespread regulation of genome stability by nuclear RNAi in higher eukaryotes.

Using next-generation sequencing to isolate mutant genes from forward genetic screens

The long-lasting success of forward genetic screens relies on the simple molecular basis of the characterized phenotypes, which are typically caused by mutations in single genes. Mapping the location of causal mutations using genetic crosses has traditionally been a complex, multistep procedure, but next-generation sequencing now allows the rapid identification of causal mutations at single-nucleotide resolution even in complex genetic backgrounds. Recent advances of this mapping-by-sequencing approach include methods that are independent of reference genome sequences, genetic crosses and any kind of linkage information, which make forward genetics amenable for species that have not been considered for forward genetic screens so far.

Mudi, a web tool for identifying mutations by bioinformatics analysis of whole-genome sequence

In forward genetics, identification of mutations is a time-consuming and laborious process. Modern whole-genome sequencing, coupled with bioinformatics analysis, has enabled fast and cost-effective mutation identification. However, for many experimental researchers, bioinformatics analysis is still a difficult aspect of whole-genome sequencing. To address this issue, we developed a browser-accessible and easy-to-use bioinformatics tool called Mutation discovery (Mudi; http://naoii.nig.ac.jp/mudi_top.html), which enables 'one-click' identification of causative mutations from whole-genome sequence data. In this study, we optimized Mudi for pooled-linkage analysis aimed at identifying mutants in yeast model systems. After raw sequencing data are uploaded, Mudi performs sequential analysis, including mapping, detection of variant alleles, filtering and removal of background polymorphisms, prioritization, and annotation. In an example study of suppressor mutants of ptr1-1 in the fission yeast Schizosaccharomyces pombe, pooled-linkage analysis with Mudi identified mip1(+) , a component of Target of Rapamycin Complex 1 (TORC1), as a novel component involved in RNA interference (RNAi)-related cell-cycle control. The accessibility of Mudi will accelerate systematic mutation analysis in forward genetics.

A genetic screen based on in vivo RNA imaging reveals centrosome-independent mechanisms for localizing gurken transcripts in Drosophila

We have screened chromosome arm 3L for ethyl methanesulfonate−induced mutations that disrupt localization of fluorescently labeled gurken (grk) messenger (m)RNA, whose transport along microtubules establishes both major body axes of the developing Drosophila oocyte. Rapid identification of causative mutations by single-nucleotide polymorphism recombinational mapping and whole-genomic sequencing allowed us to define nine complementation groups affecting grk mRNA localization and other aspects of oogenesis, including alleles of elg1, scaf6, quemao, nudE, Tsc2/gigas, rasp, and Chd5/Wrb, and several null alleles of the armitage Piwi-pathway gene. Analysis of a newly induced kinesin light chain allele shows that kinesin motor activity is required for both efficient grk mRNA localization and oocyte centrosome integrity. We also show that initiation of the dorsoanterior localization of grk mRNA precedes centrosome localization, suggesting that microtubule self-organization contributes to breaking axial symmetry to generate a unique dorsoventral axis.

Next-generation mapping of genetic mutations using bulk population sequencing

Next-generation sequencing platforms have made it possible to very rapidly map genetic mutations in Arabidopsis using whole-genome resequencing against pooled members of an F2 mapping population. In the case of recessive mutations, all individuals expressing the phenotype will be homozygous for the mutant genome at the locus responsible for the phenotype, while all other loci segregate roughly equally for both parental lines due to recombination. Importantly, genomic regions flanking the recessive mutation will be in linkage disequilibrium and therefore also be homozygous due to genetic hitchhiking. This information can be exploited to quickly and effectively identify the causal mutation. To this end, sequence data generated from members of the pooled population exhibiting the mutant phenotype are first aligned to the reference genome. Polymorphisms between the mutant and mapping line are then identified and used to determine the homozygous, nonrecombinant region harboring the mutation. Polymorphisms in the identified region are filtered to provide a short list of markers potentially responsible for the phenotype of interest, which is followed by validation at the bench. Although the focus of recent studies has been on the mapping of point mutations exhibiting recessive phenotypes, the techniques employed can be extended to incorporate more complicated scenarios such as dominant mutations and those caused by insertions or deletions in genomic sequence. This chapter describes detailed procedures for performing next-generation mapping against an Arabidopsis mutant and discusses how different mutations might be approached.

Forward genetics by genome sequencing reveals that rapid cyanide release deters insect herbivory of Sorghum bicolor

Whole genome sequencing has allowed rapid progress in the application of forward genetics in model species. In this study, we demonstrated an application of next-generation sequencing for forward genetics in a complex crop genome. We sequenced an ethyl methanesulfonate-induced mutant of Sorghum bicolor defective in hydrogen cyanide release and identified the causal mutation. A workflow identified the causal polymorphism relative to the reference BTx623 genome by integrating data from single nucleotide polymorphism identification, prior information about candidate gene(s) implicated in cyanogenesis, mutation spectra, and polymorphisms likely to affect phenotypic changes. A point mutation resulting in a premature stop codon in the coding sequence of dhurrinase2, which encodes a protein involved in the dhurrin catabolic pathway, was responsible for the acyanogenic phenotype. Cyanogenic glucosides are not cyanogenic compounds but their cyanohydrins derivatives do release cyanide. The mutant accumulated the glucoside, dhurrin, but failed to efficiently release cyanide upon tissue disruption. Thus, we tested the effects of cyanide release on insect herbivory in a genetic background in which accumulation of cyanogenic glucoside is unchanged. Insect preference choice experiments and herbivory measurements demonstrate a deterrent effect of cyanide release capacity, even in the presence of wild-type levels of cyanogenic glucoside accumulation. Our gene cloning method substantiates the value of (1) a sequenced genome, (2) a strongly penetrant and easily measurable phenotype, and (3) a workflow to pinpoint a causal mutation in crop genomes and accelerate in the discovery of gene function in the postgenomic era.

Partitioning of the nuclear and mitochondrial tRNA 3'-end processing activities between two different proteins in Schizosaccharomyces pombe

tRNase Z is an essential endonuclease responsible for tRNA 3'-end maturation. tRNase Z exists in a short form (tRNase Z(S)) and a long form (tRNase Z(L)). Prokaryotes have only tRNase Z(S), whereas eukaryotes can have both forms of tRNase Z. Most eukaryotes characterized thus far, including Saccharomyces cerevisiae, Caenorhabditis elegans, Drosophila melanogaster, and humans, contain only one tRNase Z(L) gene encoding both nuclear and mitochondrial forms of tRNase Z(L). In contrast, Schizosaccharomyces pombe contains two essential tRNase Z(L) genes (trz1 and trz2) encoding two tRNase Z(L) proteins, which are targeted to the nucleus and mitochondria, respectively. Trz1 protein levels are notably higher than Trz2 protein levels. Here, using temperature-sensitive mutants of trz1 and trz2, we provide in vivo evidence that trz1 and trz2 are involved in nuclear and mitochondrial tRNA 3'-end processing, respectively. In addition, trz2 is also involved in generation of the 5'-ends of other mitochondrial RNAs, whose 5'-ends coincide with the 3'-end of tRNA. Thus, our results provide a rare example showing partitioning of the nuclear and mitochondrial tRNase Z(L) activities between two different proteins in S. pombe. The evolution of two tRNase Z(L) genes and their differential expression in fission yeast may avoid toxic off-target effects.

Impaired coenzyme A synthesis in fission yeast causes defective mitosis, quiescence-exit failure, histone hypoacetylation and fragile DNA

Biosynthesis of coenzyme A (CoA) requires a five-step process using pantothenate and cysteine in the fission yeast Schizosaccharomyces pombe. CoA contains a thiol (SH) group, which reacts with carboxylic acid to form thioesters, giving rise to acyl-activated CoAs such as acetyl-CoA. Acetyl-CoA is essential for energy metabolism and protein acetylation, and, in higher eukaryotes, for the production of neurotransmitters. We isolated a novel S. pombe temperature-sensitive strain ppc1-537 mutated in the catalytic region of phosphopantothenoylcysteine synthetase (designated Ppc1), which is essential for CoA synthesis. The mutant becomes auxotrophic to pantothenate at permissive temperature, displaying greatly decreased levels of CoA, acetyl-CoA and histone acetylation. Moreover, ppc1-537 mutant cells failed to restore proliferation from quiescence. Ppc1 is thus the product of a super-housekeeping gene. The ppc1-537 mutant showed combined synthetic lethal defects with five of six histone deacetylase mutants, whereas sir2 deletion exceptionally rescued the ppc1-537 phenotype. In synchronous cultures, ppc1-537 cells can proceed to the S phase, but lose viability during mitosis failing in sister centromere/kinetochore segregation and nuclear division. Additionally, double-strand break repair is defective in the ppc1-537 mutant, producing fragile broken DNA, probably owing to diminished histone acetylation. The CoA-supported metabolism thus controls the state of chromosome DNA.

Rapid positional cloning of zebrafish mutations by linkage and homozygosity mapping using whole-genome sequencing

Forward genetic screens in zebrafish have identified >9000 mutants, many of which are potential disease models. Most mutants remain molecularly uncharacterized because of the high cost, time and labor investment required for positional cloning. These costs limit the benefit of previous genetic screens and discourage future screens. Drastic improvements in DNA sequencing technology could dramatically improve the efficiency of positional cloning in zebrafish and other model organisms, but the best strategy for cloning by sequencing has yet to be established. Using four zebrafish inner ear mutants, we developed and compared two approaches for 'cloning by sequencing': one based on bulk segregant linkage (BSFseq) and one based on homozygosity mapping (HMFseq). Using BSFseq we discovered that mutations in lmx1b and jagged1b cause abnormal ear morphogenesis. With HMFseq we validated that the disruption of cdh23 abolishes the ear's sensory functions and identified a candidate lesion in lhfpl5a predicted to cause nonsyndromic deafness. The success of HMFseq shows that the high intrastrain polymorphism rate in zebrafish eliminates the need for time-consuming map crosses. Additionally, we analyzed diversity in zebrafish laboratory strains to find areas of elevated diversity and areas of fixed homozygosity, reinforcing recent findings that genome diversity is clustered. We present a database of >15 million sequence variants that provides much of this approach's power. In our four test cases, only a single candidate single nucleotide polymorphism (SNP) remained after subtracting all database SNPs from a mutant's critical region. The saturation of the common SNP database and our open source analysis pipeline MegaMapper will improve the pace at which the zebrafish community makes unique discoveries relevant to human health.

Accelerating forward genetics for cell wall deconstruction

The elucidation of the genes involved in cell wall synthesis and assembly remains one of the biggest challenges of cell wall biology. Although traditional genetic approaches, using simple yet elegant screens, have identified components of the cell wall, many unknowns remain. Exhausting the genetic toolbox by performing sensitized screens, adopting chemical genetics or combining these with improved cell wall imaging, hold the promise of new gene discovery and function. With the recent introduction of next-generation sequencing technologies, it is now possible to quickly and efficiently map and clone genes of interest in record time. The combination of a classical genetics approach and cutting edge technology will propel cell wall biology in plants forward into the future.

SNP-Ratio Mapping (SRM): identifying lethal alleles and mutations in complex genetic backgrounds by next-generation sequencing

We present a generally applicable method allowing rapid identification of causal alleles in mutagenized genomes by next-generation sequencing. Currently used approaches rely on recovering homozygotes or extensive backcrossing. In contrast, SNP-ratio mapping allows rapid cloning of lethal and/or poorly transmitted mutations and second-site modifiers, which are often in complex genetic/transgenic backgrounds.

Whole-Genome Sequencing of Sordaria macrospora Mutants Identifies Developmental Genes

The study of mutants to elucidate gene functions has a long and successful history; however, to discover causative mutations in mutants that were generated by random mutagenesis often takes years of laboratory work and requires previously generated genetic and/or physical markers, or resources like DNA libraries for complementation. Here, we present an alternative method to identify defective genes in developmental mutants of the filamentous fungus Sordaria macrospora through Illumina/Solexa whole-genome sequencing. We sequenced pooled DNA from progeny of crosses of three mutants and the wild type and were able to pinpoint the causative mutations in the mutant strains through bioinformatics analysis. One mutant is a spore color mutant, and the mutated gene encodes a melanin biosynthesis enzyme. The causative mutation is a G to A change in the first base of an intron, leading to a splice defect. The second mutant carries an allelic mutation in the pro41 gene encoding a protein essential for sexual development. In the mutant, we detected a complex pattern of deletion/rearrangements at the pro41 locus. In the third mutant, a point mutation in the stop codon of a transcription factor-encoding gene leads to the production of immature fruiting bodies. For all mutants, transformation with a wild type-copy of the affected gene restored the wild-type phenotype. Our data demonstrate that whole-genome sequencing of mutant strains is a rapid method to identify developmental genes in an organism that can be genetically crossed and where a reference genome sequence is available, even without prior mapping information.

Whole-Genome Sequencing to Identify Mutants and Polymorphisms in Chlamydomonas reinhardtii

Whole-genome sequencing (WGS) provides a new platform for the identification of mutations that produce a mutant phenotype. We used Illumina sequencing to identify the mutational profile of three Chlamydomonas reinhardtii mutant strains. The three strains have more than 38,000 changes from the reference genome. NG6 is aflagellate and maps to 269 kb with only one nonsynonymous change; the V(12)E mutation falls in the FLA8 gene. Evidence that NG6 is a fla8 allele comes from swimming revertants that are either true or pseudorevertants. NG30 is aflagellate and maps to 458 kb that has six nonsynonomous changes. Evidence that NG30 has a causative nonsense allele in IFT80 comes from rescue of the nonswimming phenotype with a fragment bearing only this gene. This gene has been implicated in Jeune asphyxiating thoracic dystrophy. Electron microscopy of ift80-1 (NG30) shows a novel basal body phenotype. A bar or cap is observed over the distal end of the transition zone, which may be an intermediate in preparing the basal body for flagellar assembly. In the acetate-requiring mutant ac17, we failed to find a nonsynonymous change in the 676 kb mapped region, which is incompletely assembled. In these strains, 43% of the changes occur on two of the 17 chromosomes. The excess on chromosome 6 surrounds the mating-type locus, which has numerous rearrangements and suppressed recombination, and the changes extend beyond the mating-type locus. Unexpectedly, chromosome 16 shows an unexplained excess of single nucleotide polymorphisms and indels. Overall, WGS in combination with limited mapping allows fast and accurate identification of point mutations in Chlamydomonas.

Altered nuclear tRNA metabolism in La-deleted Schizosaccharomyces pombe is accompanied by a nutritional stress response involving Atf1p and Pcr1p that is suppressible by Xpo-t/Los1p

Deletion of the sla1(+) gene, which encodes a homologue of the human RNA-binding protein La in Schizosaccharomyces pombe, causes irregularities in tRNA processing, with altered distribution of pre-tRNA intermediates. We show, using mRNA profiling, that cells lacking sla1(+) have increased mRNAs from amino acid metabolism (AAM) genes and, furthermore, exhibit slow growth in Edinburgh minimal medium. A subset of these AAM genes is under control of the AP-1-like, stress-responsive transcription factors Atf1p and Pcr1p. Although S. pombe growth is resistant to rapamycin, sla1-Δ cells are sensitive, consistent with deficiency of leucine uptake, hypersensitivity to NH4, and genetic links to the target of rapamycin (TOR) pathway. Considering that perturbed intranuclear pre-tRNA metabolism and apparent deficiency in tRNA nuclear export in sla1-Δ cells may trigger the AAM response, we show that modest overexpression of S. pombe los1(+) (also known as Xpo-t), encoding the nuclear exportin for tRNA, suppresses the reduction in pre-tRNA levels, AAM gene up-regulation, and slow growth of sla1-Δ cells. The conclusion that emerges is that sla1(+) regulates AAM mRNA production in S. pombe through its effects on nuclear tRNA processing and probably nuclear export. Finally, the results are discussed in the context of stress response programs in Saccharomyces cerevisiae.

Rapid mapping and identification of mutations in Caenorhabditis elegans by restriction site-associated DNA mapping and genomic interval pull-down sequencing

Forward genetic screens provide a powerful approach for inferring gene function on the basis of the phenotypes associated with mutated genes. However, determining the causal mutation by traditional mapping and candidate gene sequencing is often the rate-limiting step, especially when analyzing many mutants. We report two genomic approaches for more rapidly determining the identity of the affected genes in Caenorhabditis elegans mutants. First, we report our use of restriction site-associated DNA (RAD) polymorphism markers for rapidly mapping mutations after chemical mutagenesis and mutant isolation. Second, we describe our use of genomic interval pull-down sequencing (GIPS) to selectively capture and sequence megabase-sized portions of a mutant genome. Together, these two methods provide a rapid and cost-effective approach for positional cloning of C. elegans mutant loci, and are also applicable to other genetic model systems.

Next-generation mapping of Arabidopsis genes

Next-generation genomic sequencing technologies have made it possible to directly map mutations responsible for phenotypes of interest via direct sequencing. However, most mapping strategies proposed to date require some prior genetic analysis, which can be very time-consuming even in genetically tractable organisms. Here we present a de novo method for rapidly and robustly mapping the physical location of EMS mutations by sequencing a small pooled F₂ population. This method, called Next Generation Mapping (NGM), uses a chastity statistic to quantify the relative contribution of the parental mutant and mapping lines to each SNP in the pooled F₂ population. It then uses this information to objectively localize the candidate mutation based on its exclusive segregation with the mutant parental line. A user-friendly, web-based tool for performing NGM analysis is available at http://bar.utoronto.ca/NGM. We used NGM to identify three genes involved in cell-wall biology in Arabidopsis thaliana, and, in a power analysis, demonstrate success in test mappings using as few as ten F₂ lines and a single channel of Illumina Genome Analyzer data. This strategy can easily be applied to other model organisms, and we expect that it will also have utility in crops and any other eukaryote with a completed genome sequence.

Bulk segregant analysis followed by high-throughput sequencing reveals the Neurospora cell cycle gene, ndc-1, to be allelic with the gene for ornithine decarboxylase, spe-1

With the advent of high-throughput DNA sequencing, it is now straightforward and inexpensive to generate high-density small nucleotide polymorphism (SNP) maps. Here we combined high-throughput sequencing with bulk segregant analysis to expedite mutation mapping. The general map location of a mutation can be identified by a single backcross to a strain enriched in SNPs compared to a standard wild-type strain. Bulk segregant analysis simultaneously increases the likelihood of determining the precise nature of the mutation. We present here a high-density SNP map between Neurospora crassa Mauriceville-1-c (FGSC2225) and OR74A (FGSC2489), the strains most typically used by Neurospora researchers to carry out mapping crosses. We further have demonstrated the utility of the Mauriceville sequence and our approach by mapping the mutation responsible for the only existing temperature-sensitive (ts) cell cycle mutation in Neurospora, nuclear division cycle-1 (ndc-1). The single T-to-C point mutation maps to the gene encoding ornithine decarboxylase (ODC), spe-1 (NCU01271), and changes a Phe to a Ser residue within a highly conserved motif next to the catalytic site of the enzyme. By growth on spermidine and complementation with a wild-type spe-1 gene, we showed that the defect in spe-1 is responsible for the ts ndc-1 mutation. Based on our results, we propose changing ndc-1 to spe-1(ndc), which reflects that this mutation results in an ODC with a specific nuclear division defect.

Comparative whole genome sequencing reveals phenotypic tRNA gene duplication in spontaneous Schizosaccharomyces pombe La mutants

We used a genetic screen based on tRNA-mediated suppression (TMS) in a Schizosaccharomyces pombe La protein (Sla1p) mutant. Suppressor pre-tRNA^SerUCA-C47:6U with a debilitating substitution in its variable arm fails to produce tRNA in a sla1-rrm mutant deficient for RNA chaperone-like activity. The parent strain and spontaneous mutant were analyzed using Solexa sequencing. One synonymous single-nucleotide polymorphism (SNP), unrelated to the phenotype, was identified. Further sequence analyses found a duplication of the tRNA^SerUCA-C47:6U gene, which was shown to cause the phenotype. Ninety percent of 28 isolated mutants contain duplicated tRNA^SerUCA-C47:6U genes. The tRNA gene duplication led to a disproportionately large increase in tRNA^SerUCA-C47:6U levels in sla1-rrm but not sla1-null cells, consistent with non-specific low-affinity interactions contributing to the RNA chaperone-like activity of La, similar to other RNA chaperones. Our analysis also identified 24 SNPs between ours and S. pombe 972h- strain yFS101 that was recently sequenced using Solexa. By including mitochondrial (mt) DNA in our analysis, overall coverage increased from 52% to 96%. mtDNA from our strain and yFS101 shared 14 mtSNPs relative to a ‘reference’ mtDNA, providing the first identification of these S. pombe mtDNA discrepancies. Thus, strain-specific and spontaneous phenotypic mutations can be mapped in S. pombe by Solexa sequencing.

A new dominant peroxiredoxin allele identified by whole-genome re-sequencing of random mutagenized yeast causes oxidant-resistance and premature aging

The combination of functional genomics with next generation sequencing facilitates new experimental strategies for addressing complex biological phenomena. Here, we report the identification of a gain-of-function allele of peroxiredoxin (thioredoxin peroxidase, Tsa1p) via whole-genome re-sequencing of a dominantSaccharomyces cerevisiae mutant obtained by chemical mutagenesis. Yeast strain K6001, a screening system for lifespan phenotypes, was treated with ethyl methanesulfonate (EMS). We isolated an oxidative stress-resistant mutant (B7) which transmitted this phenotype in a background-independent, monogenic and dominant way. By massive parallel pyrosequencing, we generated an 38.8 fold whole-genome coverage of the strains, which differed in 12,482 positions from the reference (S288c) genome. Via a subtraction strategy, we could narrow this number to 13 total and 4 missense nucleotide variations that were specific for the mutant. Via expression in wild type backgrounds, we show that one of these mutations, exchanging a residue in the peroxiredoxin Tsa1p, was responsible for the mutant phenotype causing background-independent dominant oxidative stress-resistance. These effects were not provoked by altered Tsa1p levels, nor could they be simulated by deletion, haploinsufficiency or over-expression of the wild-type allele. Furthermore, via both a mother-enrichment technique and a micromanipulation assay, we found a robust premature aging phenotype of this oxidant-resistant strain. Thus, TSA1-B7 encodes for a novel dominant form of peroxiredoxin, and establishes a new connection between oxidative stress and aging. In addition, this study shows that the re-sequencing of entire genomes is becoming a promising alternative for the identification of functional alleles in approaches of classic molecular genetics.

A piggyBac transposon-based mutagenesis system for the fission yeast Schizosaccharomyces pombe

The TTAA-specific transposon piggyBac (PB), originally isolated from the cabbage looper moth, Trichoplusia ni, has been utilized as an insertional mutagenesis tool in various eukaryotic organisms. Here, we show that PB transposes in the fission yeast Schizosaccharomyces pombe and leaves almost no footprints. We developed a PB-based mutagenesis system for S. pombe by constructing a strain with a selectable transposon excision marker and an integrated transposase gene. PB transposition in this strain has low chromosomal distribution bias as shown by deep sequencing-based insertion site mapping. Using this system, we obtained loss-of-function alleles of klp5 and klp6, and a gain-of-function allele of dam1 from a screen for mutants resistant to the microtubule-destabilizing drug thiabendazole. From another screen for cdc25-22 suppressors, we obtained multiple alleles of wee1 as expected. The success of these two screens demonstrated the usefulness of this PB-mediated mutagenesis tool for fission yeast.

CENP-B preserves genome integrity at replication forks paused by retrotransposon LTR

Centromere-binding protein B (CENP-B) is a widely conserved DNA binding factor associated with heterochromatin and centromeric satellite repeats1. In fission yeast, CENP-B homologs have been shown to silence Long Terminal Repeat (LTR) retrotransposons by recruiting histone deacetylases2. However, CENP-B factors also have unexplained roles in DNA replication3, 4. Here, we show that a molecular function of CENP-B is to promote replication fork progression through the LTR. Mutants have increased genomic instability caused by replication fork blockage that depends on the DNA binding factor Switch Activating Protein 1 (Sap1), which is directly recruited by the LTR. The loss of Sap1-dependent barrier activity allows the unhindered progression of the replication fork, but results in rearrangements deleterious to the retrotransposon. We conclude that retrotransposons influence replication polarity through recruitment of Sap1 and transposition near replication fork blocks, while CENP-B counteracts this activity and promotes fork stability. Our results may account for the role of LTR in fragile sites, and for the association of CENP-B with pericentromeric heterochromatin and tandem satellite repeats.

A strategy for direct mapping and identification of mutations by whole-genome sequencing

Mutant screens have proven powerful for genetic dissection of a myriad of biological processes, but subsequent identification and isolation of the causative mutations are usually complex and time consuming. We have made the process easier by establishing a novel strategy that employs whole-genome sequencing to simultaneously map and identify mutations without the need for any prior genetic mapping.

Analysis of multiple ethyl methanesulfonate-mutagenized Caenorhabditis elegans strains by whole-genome sequencing

Whole-genome sequencing (WGS) of organisms displaying a specific mutant phenotype is a powerful approach to identify the genetic determinants of a plethora of biological processes. We have previously validated the feasibility of this approach by identifying a point-mutated locus responsible for a specific phenotype, observed in an ethyl methanesulfonate (EMS)-mutagenized Caenorhabditis elegans strain. Here we describe the genome-wide mutational profile of 17 EMS-mutagenized genomes as assessed with a bioinformatic pipeline, called MAQGene. Surprisingly, we find that while outcrossing mutagenized strains does reduce the total number of mutations, a striking mutational load is still observed even in outcrossed strains. Such genetic complexity has to be taken into account when establishing a causative relationship between genotype and phenotype. Even though unintentional, the 17 sequenced strains described here provide a resource of allelic variants in almost 1000 genes, including 62 premature stop codons, which represent candidate knockout alleles that will be of further use for the C. elegans community to study gene function.

The impact of whole genome sequencing on model system genetics: get ready for the ride

Much of our understanding of how organisms develop and function is derived from the extraordinarily powerful, classic approach of screening for mutant organisms in which a specific biological process is disrupted. Reaping the fruits of such forward genetic screens in metazoan model systems like Drosophila, Caenorhabditis elegans, or zebrafish traditionally involves time-consuming positional cloning strategies that result in the identification of the mutant locus. Whole genome sequencing (WGS) has begun to provide an effective alternative to this approach through direct pinpointing of the molecular lesion in a mutated strain isolated from a genetic screen. Apart from significantly altering the pace and costs of genetic analysis, WGS also provides new perspectives on solving genetic problems that are difficult to tackle with conventional approaches, such as identifying the molecular basis of multigenic and complex traits.