A scalable and cGMP-compatible autologous organotypic cell therapy for Dystrophic Epidermolysis Bullosa.

We present Dystrophic Epidermolysis Bullosa Cell Therapy (DEBCT), a scalable platform producing autologous organotypic iPS cell-derived induced skin composite (iSC) grafts for definitive treatment. Clinical-grade manufacturing integrates CRISPR-mediated genetic correction with reprogramming into one step, accelerating derivation of COL7A1-edited iPS cells from patients. Differentiation into epidermal, dermal and melanocyte progenitors is followed by CD49f-enrichment, minimizing maturation heterogeneity. Mouse xenografting of iSCs from four patients with different mutations demonstrates disease modifying activity at 1 month. Next-generation sequencing, biodistribution and tumorigenicity assays establish a favorable safety profile at 1-9 months. Single cell transcriptomics reveals that iSCs are composed of the major skin cell lineages and include prominent holoclone stem cell-like signatures of keratinocytes, and the recently described Gibbin-dependent signature of fibroblasts. The latter correlates with enhanced graftability of iSCs. In conclusion, DEBCT overcomes manufacturing and safety roadblocks and establishes a reproducible, safe, and cGMP-compatible therapeutic approach to heal lesions of DEB patients.

Cohesin-mediated 3D contacts tune enhancer-promoter regulation.

Enhancers are key drivers of gene regulation thought to act via 3D physical interactions with the promoters of their target genes. However, genome-wide depletions of architectural proteins such as cohesin result in only limited changes in gene expression, despite a loss of contact domains and loops. Consequently, the role of cohesin and 3D contacts in enhancer function remains debated. Here, we developed CRISPRi of regulatory elements upon degron operation (CRUDO), a novel approach to measure how changes in contact frequency impact enhancer effects on target genes by perturbing enhancers with CRISPRi and measuring gene expression in the presence or absence of cohesin. We systematically perturbed all 1,039 candidate enhancers near five cohesin-dependent genes and identified 34 enhancer-gene regulatory interactions. Of 26 regulatory interactions with sufficient statistical power to evaluate cohesin dependence, 18 show cohesin-dependent effects. A decrease in enhancer-promoter contact frequency upon removal of cohesin is frequently accompanied by a decrease in the regulatory effect of the enhancer on gene expression, consistent with a contact-based model for enhancer function. However, changes in contact frequency and regulatory effects on gene expression vary as a function of distance, with distal enhancers (e.g., >50Kb) experiencing much larger changes than proximal ones (e.g., <50Kb). Because most enhancers are located close to their target genes, these observations can explain how only a small subset of genes - those with strong distal enhancers - are sensitive to cohesin. Together, our results illuminate how 3D contacts, influenced by both cohesin and genomic distance, tune enhancer effects on gene expression.

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.

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.

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.

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.

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.

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.

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.

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.