skip to main content


Title: A multiplexed, three-dimensional pooling and next-generation sequencing strategy for creating barcoded mutant arrays: construction of a Schizosaccharomyces pombe transposon insertion library
Abstract Arrayed libraries of defined mutants have been used to elucidate gene function in the post-genomic era. Yeast haploid gene deletion libraries have pioneered this effort, but are costly to construct, do not reveal phenotypes that may occur with partial gene function and lack essential genes required for growth. We therefore devised an efficient method to construct a library of barcoded insertion mutants with a wider range of phenotypes that can be generalized to other organisms or collections of DNA samples. We developed a novel but simple three-dimensional pooling and multiplexed sequencing approach that leveraged sequence information to reduce the number of required sequencing reactions by orders of magnitude, and were able to identify the barcode sequences and DNA insertion sites of 4391 Schizosaccharomyces pombe insertion mutations with only 40 sequencing preparations. The insertion mutations are in the genes and untranslated regions of nonessential, essential and noncoding RNA genes, and produced a wider range of phenotypes compared to the cognate deletion mutants, including novel phenotypes. This mutant library represents both a proof of principle for an efficient method to produce novel mutant libraries and a valuable resource for the S. pombe research community.  more » « less
Award ID(s):
1908875
NSF-PAR ID:
10345124
Author(s) / Creator(s):
; ; ; ; ; ; ;
Date Published:
Journal Name:
Nucleic Acids Research
ISSN:
0305-1048
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Auxin is a hormone that is required for hypocotyl elongation during seedling development. In response to auxin, rapid changes in transcript and protein abundance occur in hypocotyls, and some auxin responsive gene expression is linked to hypocotyl growth. To functionally validate proteomic studies, a reverse genetics screen was performed on mutants in auxin‐regulated proteins to identify novel regulators of plant growth. This uncovered a long hypocotyl mutant, which we calledslim shady, in an annotated insertion line inIMMUNOREGULATORY RNA‐BINDING PROTEIN(IRR). Overexpression of theIRRgene failed to rescue theslim shadyphenotype and characterization of a second T‐DNA allele of IRR found that it had a wild‐type (WT) hypocotyl length. Theslim shadymutant has an elevated expression of numerous genes associated with the brassinosteroid‐auxin‐phytochrome (BAP) regulatory module compared to WT, including transcription factors that regulate brassinosteroid, auxin, and phytochrome pathways. Additionally,slim shadyseedlings fail to exhibit a strong transcriptional response to auxin. Using whole genome sequence data and genetic complementation analysis with SALK_015201C, we determined that a novel single nucleotide polymorphism inPHYTOCHROME Bwas responsible for theslim shadyphenotype. This is predicted to induce a frameshift and premature stop codon at leucine 1125, within the histidine kinase‐related domain of the carboxy terminus of PHYB, which is required for phytochrome signaling and function. Genetic complementation analyses withphyb‐9confirmed thatslim shadyis a mutant allele ofPHYB. This study advances our understanding of the molecular mechanisms in seedling development, by furthering our understanding of how light signaling is linked to auxin‐dependent cell elongation. Furthermore, this study highlights the importance of confirming the genetic identity of research material before attributing phenotypes to known mutations sourced from T‐DNA stocks.

     
    more » « less
  2. SUMMARY

    Flowering of the reference legumeMedicago truncatulais promoted by winter cold (vernalization) followed by long‐day photoperiods (VLD) similar to winter annual Arabidopsis. However, Medicago lacksFLCandCO, key regulators of Arabidopsis VLD flowering.Most plants have twoINHIBITOR OF GROWTH(ING) genes (ING1andING2), encoding proteins with an ING domain with two anti‐parallel alpha‐helices and a plant homeodomain (PHD) finger, but their genetic role has not been previously described.In Medicago,Mting1gene‐edited mutants developed and flowered normally, but anMting2‐1 Tnt1insertion mutant and gene‐editedMting2mutants had developmental abnormalities including delayed flowering particularly in VLD, compact architecture, abnormal leaves with extra leaflets but no trichomes, and smaller seeds and barrels.Mting2mutants had reduced expression of activators of flowering, including theFT‐like geneMtFTa1, and increased expression of the candidate repressorMtTFL1c, consistent with the delayed flowering of the mutant.MtING2overexpression complementedMting2‐1, but did not accelerate flowering in wild type. The MtING2 PHD finger bound H3K4me2/3 peptides weaklyin vitro, but analysis of gene‐edited mutants indicated that it was dispensable to MtING2 function in wild‐type plants. RNA sequencing experiments indicated that >7000 genes are mis‐expressed in theMting2‐1mutant, consistent with its strong mutant phenotypes. Interestingly, ChIP‐seq analysis identified >5000 novel H3K4me3 locations in the genome ofMting2‐1mutants compared to wild type R108. Overall, our mutant study has uncovered an important physiological role of a plantING2gene in development, flowering, and gene expression, which likely involves an epigenetic mechanism.

     
    more » « less
  3. Transcription factor (TF)–promoter pairs have been repurposed from native hosts to provide tools to measure intracellular biochemical production titer and dynamically control gene expression. Most often, native TF–promoter systems require rigorous screening to obtain desirable characteristics optimized for biotechnological applications. High-throughput techniques may provide a rational and less labor-intensive strategy to engineer user-defined TF–promoter pairs using fluorescence-activated cell sorting and deep sequencing methods (sort-seq). Based on the designed promoter library’s distribution characteristics, we elucidate sequence–function interactions between the TF and DNA. In this work, we use the sort-seq method to study the sequence–function relationship of a σ 54 -dependent, butanol-responsive TF–promoter pair, BmoR-P BMO derived from Thauera butanivorans , at the nucleotide level to improve biosensor characteristics, specifically an improved dynamic range. Activities of promoters from a mutagenized P BMO library were sorted based on gfp expression and subsequently deep sequenced to correlate site-specific sequences with changes in dynamic range. We identified site-specific mutations that increase the sensor output. Double mutant and a single mutant, CA(129,130)TC and G(205)A, in P BMO promoter increased dynamic ranges of 4-fold and 1.65-fold compared with the native system, respectively. In addition, sort-seq identified essential sites required for the proper function of the σ 54 -dependent promoter biosensor in the context of the host. This work can enable high-throughput screening methods for strain development. 
    more » « less
  4. Transcription factor (TF)–promoter pairs have been repurposed from native hosts to provide tools to measure intracellular biochemical production titer and dynamically control gene expression. Most often, native TF–promoter systems require rigorous screening to obtain desirable characteristics optimized for biotechnological applications. High-throughput techniques may provide a rational and less labor-intensive strategy to engineer user-defined TF–promoter pairs using fluorescence-activated cell sorting and deep sequencing methods (sort-seq). Based on the designed promoter library’s distribution characteristics, we elucidate sequence–function interactions between the TF and DNA. In this work, we use the sort-seq method to study the sequence–function relationship of a σ54-dependent, butanol-responsive TF–promoter pair, BmoR-PBMO derived from Thauera butanivorans, at the nucleotide level to improve biosensor characteristics, specifically an improved dynamic range. Activities of promoters from a mutagenized PBMO library were sorted based on gfp expression and subsequently deep sequenced to correlate site-specific sequences with changes in dynamic range. We identified site-specific mutations that increase the sensor output. Double mutant and a single mutant, CA(129,130)TC and G(205)A, in PBMO promoter increased dynamic ranges of 4-fold and 1.65-fold compared with the native system, respectively. In addition, sort-seq identified essential sites required for the proper function of the σ54-dependent promoter biosensor in the context of the host. This work can enable high-throughput screening methods for strain development. 
    more » « less
  5. As transposon sequencing (TnSeq) assays have become prolific in the microbiology field, it is of interest to scrutinize their potential drawbacks. TnSeq data consist of millions of nucleotide sequence reads that are generated by PCR amplification of transposon-genomic junctions. Reads mapping to the junctions are enumerated thus providing information on the number of transposon insertion mutations in each individual gene. Here we explore the possibility that PCR amplification of transposon insertions in a TnSeq library skews the results by introducing bias into the detection and/or enumeration of insertions. We compared the detection and frequency of mapped insertions when altering the number of PCR cycles, and when including a nested PCR, in the enrichment step. Additionally, we present nCATRAs - a novel, amplification-free TnSeq method where the insertions are enriched via CRISPR/Cas9-targeted transposon cleavage and subsequent Oxford Nanopore MinION sequencing. nCATRAs achieved 54 and 23% enrichment of the transposons and transposon-genomic junctions, respectively, over background genomic DNA. These PCR-based and PCR-free experiments demonstrate that, overall, PCR amplification does not significantly bias the results of TnSeq insofar as insertions in the majority of genes represented in our library were similarly detected regardless of PCR cycle number and whether or not PCR amplification was employed. However, the detection of a small subset of genes which had been previously described as essential is sensitive to the number of PCR cycles. We conclude that PCR-based enrichment of transposon insertions in a TnSeq assay is reliable, but researchers interested in profiling putative essential genes should carefully weigh the number of amplification cycles employed in their library preparation protocols. In addition, nCATRAs is comparable to traditional PCR-based methods (Kendall’s correlation=0.896–0.897) although the latter remain superior owing to their accessibility and high sequencing depth. 
    more » « less