Genome-scale analysis of interactions between genetic perturbations and natural variation.

Interactions between genetic perturbations and segregating loci can cause perturbations to show different phenotypic effects across genetically distinct individuals. To study these interactions on a genome scale in many individuals, we used combinatorial DNA barcode sequencing to measure the fitness effects of 8046 CRISPRi perturbations targeting 1721 distinct genes in 169 yeast cross progeny (or segregants). We identified 460 genes whose perturbation has different effects across segregants. Several factors caused perturbations to show variable effects, including baseline segregant fitness, the mean effect of a perturbation across segregants, and interacting loci. We mapped 234 interacting loci and found four hub loci that interact with many different perturbations. Perturbations that interact with a given hub exhibit similar epistatic relationships with the hub and show enrichment for cellular processes that may mediate these interactions. These results suggest that an individual's response to perturbations is shaped by a network of perturbation-locus interactions that cannot be measured by approaches that examine perturbations or natural variation alone.

Histone deacetylase Hos2 regulates protein expression noise by potentially modulating the protein translation machinery.

Non-genetic variations derived from expression noise at transcript or protein levels can result in cell-to-cell heterogeneity within an isogenic population. Although cells have developed strategies to reduce noise in some cellular functions, this heterogeneity can also facilitate varying levels of regulation and provide evolutionary benefits in specific environments. Despite several general characteristics of cellular noise having been revealed, the detailed molecular pathways underlying noise regulation remain elusive. Here, we established a dual-fluorescent reporter system in Saccharomyces cerevisiae and performed experimental evolution to search for mutations that increase expression noise. By analyzing evolved cells using bulk segregant analysis coupled with whole-genome sequencing, we identified the histone deacetylase Hos2 as a negative noise regulator. A hos2 mutant down-regulated multiple ribosomal protein genes and exhibited partially compromised protein translation, indicating that Hos2 may regulate protein expression noise by modulating the translation machinery. Treating cells with translation inhibitors or introducing mutations into several Hos2-regulated ribosomal protein genes-RPS9A, RPS28B and RPL42A-enhanced protein expression noise. Our study provides an effective strategy for identifying noise regulators and also sheds light on how cells regulate non-genetic variation through protein translation.

Dissecting quantitative trait nucleotides by saturation genome editing.

Genome editing technologies have the potential to transform our understanding of how genetic variation gives rise to complex traits through the systematic engineering and phenotypic characterization of genetic variants. However, there has yet to be a system with sufficient efficiency, fidelity, and throughput to comprehensively identify causal variants at the genome scale. Here we explored the ability of templated CRISPR editing systems to install natural variants genome-wide in budding yeast. We optimized several approaches to enhance homology-directed repair (HDR) with donor DNA templates, including donor recruitment to target sites, single-stranded donor production by bacterial retrons, and in vivo plasmid assembly. We uncovered unique advantages of each system that we integrated into a single superior system named MAGESTIC 3.0. We used MAGESTIC 3.0 to dissect causal variants residing in 112 quantitative trait loci across 32 environmental conditions, revealing an enrichment for missense variants and loci with multiple causal variants. MAGESTIC 3.0 will facilitate the functional analysis of the genome at single-nucleotide resolution and provides a roadmap for improving template-based genome editing systems in other organisms.

Size-exclusion chromatography combined with DIA-MS enables deep proteome profiling of extracellular vesicles from melanoma plasma and serum.

Extracellular vesicles (EVs) are important players in melanoma progression, but their use as clinical biomarkers has been limited by the difficulty of profiling blood-derived EV proteins with high depth of coverage, the requirement for large input amounts, and complex protocols. Here, we provide a streamlined and reproducible experimental workflow to identify plasma- and serum- derived EV proteins of healthy donors and melanoma patients using minimal amounts of sample input. SEC-DIA-MS couples size-exclusion chromatography to EV concentration and deep-proteomic profiling using data-independent acquisition. From as little as 200 µL of plasma per patient in a cohort of three healthy donors and six melanoma patients, we identified and quantified 2896 EV-associated proteins, achieving a 3.5-fold increase in depth compared to previously published melanoma studies. To compare the EV-proteome to unenriched blood, we employed an automated workflow to deplete the 14 most abundant proteins from plasma and serum and thereby approximately doubled protein group identifications versus native blood. The EV proteome diverged from corresponding unenriched plasma and serum, and unlike the latter, separated healthy donor and melanoma patient samples. Furthermore, known melanoma markers, such as MCAM, TNC, and TGFBI, were upregulated in melanoma EVs but not in depleted melanoma plasma, highlighting the specific information contained in EVs. Overall, EVs were significantly enriched in intact membrane proteins and proteins related to SNARE protein interactions and T-cell biology. Taken together, we demonstrated the increased sensitivity of an EV-based proteomic workflow that can be easily applied to larger melanoma cohorts and other indications.

Uterine injury during diestrus leads to placental and embryonic defects in future pregnancies in mice†.

Uterine injury from procedures such as Cesarean sections (C-sections) often have severe consequences on subsequent pregnancy outcomes, leading to disorders such as placenta previa, placenta accreta, and infertility. With rates of C-section at ~30% of deliveries in the USA and projected to continue to climb, a deeper understanding of the mechanisms by which these pregnancy disorders arise and opportunities for intervention are needed. Here we describe a rodent model of uterine injury on subsequent in utero outcomes. We observed three distinct phenotypes: increased rates of resorption and death, embryo spacing defects, and placenta accreta-like features of reduced decidua and expansion of invasive trophoblasts. We show that the appearance of embryo spacing defects depends entirely on the phase of estrous cycle at the time of injury. Using RNA-seq, we identified perturbations in the expression of components of the COX/prostaglandin pathway after recovery from injury, a pathway that has previously been demonstrated to play an important role in embryo spacing. Therefore, we demonstrate that uterine damage in this mouse model causes morphological and molecular changes that ultimately lead to placental and embryonic developmental defects.

Genome-scale analysis of interactions between genetic perturbations and natural variation.

Interactions between genetic perturbations and segregating loci can cause perturbations to show different phenotypic effects across genetically distinct individuals. To study these interactions on a genome scale in many individuals, we used combinatorial DNA barcode sequencing to measure the fitness effects of 7,700 CRISPRi perturbations targeting 1,712 distinct genes in 169 yeast cross progeny (or segregants). We identified 460 genes whose perturbation has different effects across segregants. Several factors caused perturbations to show variable effects, including baseline segregant fitness, the mean effect of a perturbation across segregants, and interacting loci. We mapped 234 interacting loci and found four hub loci that interact with many different perturbations. Perturbations that interact with a given hub exhibit similar epistatic relationships with the hub and show enrichment for cellular processes that may mediate these interactions. These results suggest that an individual's response to perturbations is shaped by a network of perturbation-locus interactions that cannot be measured by approaches that examine perturbations or natural variation alone.

Debugging and consolidating multiple synthetic chromosomes reveals combinatorial genetic interactions.

The Sc2.0 project is building a eukaryotic synthetic genome from scratch. A major milestone has been achieved with all individual Sc2.0 chromosomes assembled. Here, we describe the consolidation of multiple synthetic chromosomes using advanced endoreduplication intercrossing with tRNA expression cassettes to generate a strain with 6.5 synthetic chromosomes. The 3D chromosome organization and transcript isoform profiles were evaluated using Hi-C and long-read direct RNA sequencing. We developed CRISPR Directed Biallelic URA3-assisted Genome Scan, or "CRISPR D-BUGS," to map phenotypic variants caused by specific designer modifications, known as "bugs." We first fine-mapped a bug in synthetic chromosome II (synII) and then discovered a combinatorial interaction associated with synIII and synX, revealing an unexpected genetic interaction that links transcriptional regulation, inositol metabolism, and tRNASerCGA abundance. Finally, to expedite consolidation, we employed chromosome substitution to incorporate the largest chromosome (synIV), thereby consolidating >50% of the Sc2.0 genome in one strain.

Design, construction, and functional characterization of a tRNA neochromosome in yeast.

Here, we report the design, construction, and characterization of a tRNA neochromosome, a designer chromosome that functions as an additional, de novo counterpart to the native complement of Saccharomyces cerevisiae. Intending to address one of the central design principles of the Sc2.0 project, the ∼190-kb tRNA neochromosome houses all 275 relocated nuclear tRNA genes. To maximize stability, the design incorporates orthogonal genetic elements from non-S. cerevisiae yeast species. Furthermore, the presence of 283 rox recombination sites enables an orthogonal tRNA SCRaMbLE system. Following construction in yeast, we obtained evidence of a potent selective force, manifesting as a spontaneous doubling in cell ploidy. Furthermore, tRNA sequencing, transcriptomics, proteomics, nucleosome mapping, replication profiling, FISH, and Hi-C were undertaken to investigate questions of tRNA neochromosome behavior and function. Its construction demonstrates the remarkable tractability of the yeast model and opens up opportunities to directly test hypotheses surrounding these essential non-coding RNAs.

Consequences of a telomerase-related fitness defect and chromosome substitution technology in yeast synIX strains.

We describe the complete synthesis, assembly, debugging, and characterization of a synthetic 404,963 bp chromosome, synIX (synthetic chromosome IX). Combined chromosome construction methods were used to synthesize and integrate its left arm (synIXL) into a strain containing previously described synIXR. We identified and resolved a bug affecting expression of EST3, a crucial gene for telomerase function, producing a synIX strain with near wild-type fitness. To facilitate future synthetic chromosome consolidation and increase flexibility of chromosome transfer between distinct strains, we combined chromoduction, a method to transfer a whole chromosome between two strains, with conditional centromere destabilization to substitute a chromosome of interest for its native counterpart. Both steps of this chromosome substitution method were efficient. We observed that wild-type II tended to co-transfer with synIX and was co-destabilized with wild-type IX, suggesting a potential gene dosage compensation relationship between these chromosomes.

An encyclopedia of enhancer-gene regulatory interactions in the human genome.

Identifying transcriptional enhancers and their target genes is essential for understanding gene regulation and the impact of human genetic variation on disease1-6. Here we create and evaluate a resource of >13 million enhancer-gene regulatory interactions across 352 cell types and tissues, by integrating predictive models, measurements of chromatin state and 3D contacts, and largescale genetic perturbations generated by the ENCODE Consortium7. We first create a systematic benchmarking pipeline to compare predictive models, assembling a dataset of 10,411 elementgene pairs measured in CRISPR perturbation experiments, >30,000 fine-mapped eQTLs, and 569 fine-mapped GWAS variants linked to a likely causal gene. Using this framework, we develop a new predictive model, ENCODE-rE2G, that achieves state-of-the-art performance across multiple prediction tasks, demonstrating a strategy involving iterative perturbations and supervised machine learning to build increasingly accurate predictive models of enhancer regulation. Using the ENCODE-rE2G model, we build an encyclopedia of enhancer-gene regulatory interactions in the human genome, which reveals global properties of enhancer networks, identifies differences in the functions of genes that have more or less complex regulatory landscapes, and improves analyses to link noncoding variants to target genes and cell types for common, complex diseases. By interpreting the model, we find evidence that, beyond enhancer activity and 3D enhancer-promoter contacts, additional features guide enhancerpromoter communication including promoter class and enhancer-enhancer synergy. Altogether, these genome-wide maps of enhancer-gene regulatory interactions, benchmarking software, predictive models, and insights about enhancer function provide a valuable resource for future studies of gene regulation and human genetics.

Identification and quantification of small exon-containing isoforms in long-read RNA sequencing data.

Small exons are pervasive in transcriptomes across organisms, and their quantification in RNA isoforms is crucial for understanding gene functions. Although long-read RNA-seq based on Oxford Nanopore Technologies (ONT) offers the advantage of covering transcripts in full length, its lower base accuracy poses challenges for identifying individual exons, particularly microexons (≤ 30 nucleotides). Here, we systematically assess small exons quantification in synthetic and human ONT RNA-seq datasets. We demonstrate that reads containing small exons are often not properly aligned, affecting the quantification of relevant transcripts. Thus, we develop a local-realignment method for misaligned exons (MisER), which remaps reads with misaligned exons to the transcript references. Using synthetic and simulated datasets, we demonstrate the high sensitivity and specificity of MisER for the quantification of transcripts containing small exons. Moreover, MisER enabled us to identify small exons with a higher percent spliced-in index (PSI) in neural, particularly neural-regulated microexons, when comparing 14 neural to 16 non-neural tissues in humans. Our work introduces an improved quantification method for long-read RNA-seq and especially facilitates studies using ONT long-reads to elucidate the regulation of genes involving small exons.

Digital assay for rapid electronic quantification of clinical pathogens using DNA nanoballs.

Fast and accurate detection of nucleic acids is key for pathogen identification. Methods for DNA detection generally rely on fluorescent or colorimetric readout. The development of label-free assays decreases costs and test complexity. We present a novel method combining a one-pot isothermal generation of DNA nanoballs with their detection by electrical impedance. We modified loop-mediated isothermal amplification by using compaction oligonucleotides that self-assemble the amplified target into nanoballs. Next, we use capillary-driven flow to passively pass these nanoballs through a microfluidic impedance cytometer, thus enabling a fully compact system with no moving parts. The movement of individual nanoballs is detected by a change in impedance providing a quantized readout. This approach is flexible for the detection of DNA/RNA of numerous targets (severe acute respiratory syndrome coronavirus 2, HIV, ß-lactamase gene, etc.), and we anticipate that its integration into a standalone device would provide an inexpensive (<$5), sensitive (10 target copies), and rapid test (<1 hour).

Gut barrier defects, intestinal immune hyperactivation and enhanced lipid catabolism drive lethality in NGLY1-deficient Drosophila.

Intestinal barrier dysfunction leads to inflammation and associated metabolic changes. However, the relative impact of gut bacteria versus non-bacterial insults on animal health in the context of barrier dysfunction is not well understood. Here, we establish that loss of Drosophila N-glycanase 1 (Pngl) in a specific intestinal cell type leads to gut barrier defects, causing starvation and JNK overactivation. These abnormalities, along with loss of Pngl in enterocytes and fat body, result in Foxo overactivation, leading to hyperactive innate immune response and lipid catabolism and thereby contributing to lethality. Germ-free rearing of Pngl mutants rescued their developmental delay but not lethality. However, raising Pngl mutants on isocaloric, fat-rich diets partially rescued lethality. Our data indicate that Pngl functions in Drosophila larvae to establish the gut barrier, and that the lethality caused by loss of Pngl is primarily mediated through non-bacterial induction of immune and metabolic abnormalities.

Pooled Genome-Scale CRISPR Screens in Single Cells.

Assigning functions to genes and learning how to control their expression are part of the foundation of cell biology and therapeutic development. An efficient and unbiased method to accomplish this is genetic screening, which historically required laborious clone generation and phenotyping and is still limited by scale today. The rapid technological progress on modulating gene function with CRISPR-Cas and measuring it in individual cells has now relaxed the major experimental constraints and enabled pooled screening with complex readouts from single cells. Here, we review the principles and practical considerations for pooled single-cell CRISPR screening. We discuss perturbation strategies, experimental model systems, matching the perturbation to the individual cells, reading out cell phenotypes, and data analysis. Our focus is on single-cell RNA sequencing and cell sorting-based readouts, including image-enabled cell sorting. We expect this transformative approach to fuel biomedical research for the next several decades.

Large scale microfluidic CRISPR screening for increased amylase secretion in yeast.

Key to our ability to increase recombinant protein production through secretion is a better understanding of the pathways that interact to translate, process and export mature proteins to the surrounding environment, including the supporting cellular machinery that supplies necessary energy and building blocks. By combining droplet microfluidic screening with large-scale CRISPR libraries that perturb the expression of the majority of coding and non-coding genes in S. cerevisiae, we identified 345 genes for which an increase or decrease in gene expression resulted in increased secretion of α-amylase. Our results show that modulating the expression of genes involved in the trafficking of vesicles, endosome to Golgi transport, the phagophore assembly site, the cell cycle and energy supply improve α-amylase secretion. Besides protein-coding genes, we also find multiple long non-coding RNAs enriched in the vicinity of genes associated with endosomal, Golgi and vacuolar processes. We validated our results by overexpressing or deleting selected genes, which resulted in significant improvements in α-amylase secretion. The advantages, in terms of precision and speed, inherent to CRISPR based perturbations, enables iterative testing of new strains for increased protein secretion.

Mislocalization of pathogenic RBM20 variants in dilated cardiomyopathy is caused by loss-of-interaction with Transportin-3.

Severe forms of dilated cardiomyopathy (DCM) are associated with point mutations in the alternative splicing regulator RBM20 that are frequently located in the arginine/serine-rich domain (RS-domain). Such mutations can cause defective splicing and cytoplasmic mislocalization, which leads to the formation of detrimental cytoplasmic granules. Successful development of personalized therapies requires identifying the direct mechanisms of pathogenic RBM20 variants. Here, we decipher the molecular mechanism of RBM20 mislocalization and its specific role in DCM pathogenesis. We demonstrate that mislocalized RBM20 RS-domain variants retain their splice regulatory activity, which reveals that aberrant cellular localization is the main driver of their pathological phenotype. A genome-wide CRISPR knockout screen combined with image-enabled cell sorting identified Transportin-3 (TNPO3) as the main nuclear importer of RBM20. We show that the direct RBM20-TNPO3 interaction involves the RS-domain, and is disrupted by pathogenic variants. Relocalization of pathogenic RBM20 variants to the nucleus restores alternative splicing and dissolves cytoplasmic granules in cell culture and animal models. These findings provide proof-of-principle for developing therapeutic strategies to restore RBM20's nuclear localization in RBM20-DCM patients.

Striated muscle-specific base editing enables correction of mutations causing dilated cardiomyopathy.

Dilated cardiomyopathy is the second most common cause for heart failure with no cure except a high-risk heart transplantation. Approximately 30% of patients harbor heritable mutations which are amenable to CRISPR-based gene therapy. However, challenges related to delivery of the editing complex and off-target concerns hamper the broad applicability of CRISPR agents in the heart. We employ a combination of the viral vector AAVMYO with superior targeting specificity of heart muscle tissue and CRISPR base editors to repair patient mutations in the cardiac splice factor Rbm20, which cause aggressive dilated cardiomyopathy. Using optimized conditions, we repair >70% of cardiomyocytes in two Rbm20 knock-in mouse models that we have generated to serve as an in vivo platform of our editing strategy. Treatment of juvenile mice restores the localization defect of RBM20 in 75% of cells and splicing of RBM20 targets including TTN. Three months after injection, cardiac dilation and ejection fraction reach wild-type levels. Single-nuclei RNA sequencing uncovers restoration of the transcriptional profile across all major cardiac cell types and whole-genome sequencing reveals no evidence for aberrant off-target editing. Our study highlights the potential of base editors combined with AAVMYO to achieve gene repair for treatment of hereditary cardiac diseases.

Atlas of mRNA translation and decay for bacteria.

Regulation of messenger RNA stability is pivotal for programmed gene expression in bacteria and is achieved by a myriad of molecular mechanisms. By bulk sequencing of 5' monophosphorylated mRNA decay intermediates (5'P), we show that cotranslational mRNA degradation is conserved among both Gram-positive and -negative bacteria. We demonstrate that, in species with 5'-3' exonucleases, the exoribonuclease RNase J tracks the trailing ribosome to produce an in vivo single-nucleotide toeprint of the 5' position of the ribosome. In other species lacking 5'-3' exonucleases, ribosome positioning alters endonucleolytic cleavage sites. Using our metadegradome (5'P degradome) sequencing approach, we characterize 5'P mRNA decay intermediates in 96 species including Bacillus subtilis, Escherichia coli, Synechocystis spp. and Prevotella copri and identify codon- and gene-level ribosome stalling responses to stress and drug treatment. We also apply 5'P sequencing to complex clinical and environmental microbiomes and demonstrate that metadegradome sequencing provides fast, species-specific posttranscriptional characterization of responses to drug or environmental perturbations. Finally we produce a degradome atlas for 96 species to enable analysis of mechanisms of RNA degradation in bacteria. Our work paves the way for the application of metadegradome sequencing to investigation of posttranscriptional regulation in unculturable species and complex microbial communities.

Gut barrier defects, increased intestinal innate immune response, and enhanced lipid catabolism drive lethality in N -glycanase 1 deficient Drosophila.

Intestinal barrier dysfunction leads to inflammation and associated metabolic changes. However, the relative impact of infectious versus non-infectious mechanisms on animal health in the context of barrier dysfunction is not well understood. Here, we establish that loss of Drosophila N -glycanase 1 (Pngl) leads to gut barrier defects, which cause starvation and increased JNK activity. These defects result in Foxo overactivation, which induces a hyperactive innate immune response and lipid catabolism, thereby contributing to lethality associated with loss of Pngl . Notably, germ-free rearing of Pngl mutants did not rescue lethality. In contrast, raising Pngl mutants on isocaloric, fat-rich diets improved animal survival in a dosage-dependent manner. Our data indicate that Pngl functions in Drosophila larvae to establish the gut barrier, and that the immune and metabolic consequences of loss of Pngl are primarily mediated through non-infectious mechanisms.

A scalable, GMP-compatible, autologous organotypic cell therapy for Dystrophic Epidermolysis Bullosa.

Background: Gene editing in induced pluripotent stem (iPS) cells has been hailed to enable new cell therapies for various monogenetic diseases including dystrophic epidermolysis bullosa (DEB). However, manufacturing, efficacy and safety roadblocks have limited the development of genetically corrected, autologous iPS cell-based therapies. Methods: We developed Dystrophic Epidermolysis Bullosa Cell Therapy (DEBCT), a new generation GMP-compatible (cGMP), reproducible, and scalable platform to produce autologous clinical-grade iPS cell-derived organotypic induced skin composite (iSC) grafts to treat incurable wounds of patients lacking type VII collagen (C7). DEBCT uses a combined high-efficiency reprogramming and CRISPR-based genetic correction single step to generate genome scar-free, COL7A1 corrected clonal iPS cells from primary patient fibroblasts. Validated iPS cells are converted into epidermal, dermal and melanocyte progenitors with a novel 2D organoid differentiation protocol, followed by CD49f enrichment and expansion to minimize maturation heterogeneity. iSC product characterization by single cell transcriptomics was followed by mouse xenografting for disease correcting activity at 1 month and toxicology analysis at 1-6 months. Culture-acquired mutations, potential CRISPR-off targets, and cancer-driver variants were evaluated by targeted and whole genome sequencing. Findings: iPS cell-derived iSC grafts were reproducibly generated from four recessive DEB patients with different pathogenic mutations. Organotypic iSC grafts onto immune-compromised mice developed into stable stratified skin with functional C7 restoration. Single cell transcriptomic characterization of iSCs revealed prominent holoclone stem cell signatures in keratinocytes and the recently described Gibbin-dependent signature in dermal fibroblasts. The latter correlated with enhanced graftability. Multiple orthogonal sequencing and subsequent computational approaches identified random and non-oncogenic mutations introduced by the manufacturing process. Toxicology revealed no detectable tumors after 3-6 months in DEBCT-treated mice. Interpretation: DEBCT successfully overcomes previous roadblocks and represents a robust, scalable, and safe cGMP manufacturing platform for production of a CRISPR-corrected autologous organotypic skin graft to heal DEB patient wounds.