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.

Ehf and Fezf2 regulate late medullary thymic epithelial cell and thymic tuft cell development

Thymic epithelial cells are indispensable for T cell maturation and selection and the induction of central immune tolerance. The self-peptide repertoire expressed by medullary thymic epithelial cells is in part regulated by the transcriptional regulator Aire (Autoimmune regulator) and the transcription factor Fezf2. Due to the high complexity of mTEC maturation stages (i.e., post-Aire, Krt10+ mTECs, and Dclk1+ Tuft mTECs) and the heterogeneity in their gene expression profiles (i.e., mosaic expression patterns), it has been challenging to identify the additional factors complementing the transcriptional regulation. We aimed to identify the transcriptional regulators involved in the regulation of mTEC development and self-peptide expression in an unbiased and genome-wide manner. We used ATAC footprinting analysis as an indirect approach to identify transcription factors involved in the gene expression regulation in mTECs, which we validated by ChIP sequencing. This study identifies Fezf2 as a regulator of the recently described thymic Tuft cells (i.e., Tuft mTECs). Furthermore, we identify that transcriptional regulators of the ELF, ESE, ERF, and PEA3 subfamily of the ETS transcription factor family and members of the Krüppel-like family of transcription factors play a role in the transcriptional regulation of genes involved in late mTEC development and promiscuous gene expression.

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.

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.

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.

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.

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.

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 disease 1-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 large-scale genetic perturbations generated by the ENCODE Consortium 7 . We first create a systematic benchmarking pipeline to compare predictive models, assembling a dataset of 10,411 element-gene 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 enhancer-promoter 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.

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.

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.

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).

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.

Cell sorters see things more clearly now

Microscopy and fluorescence-activated cell sorting (FACS) are two of the most important tools for single-cell phenotyping in basic and biomedical research. Microscopy provides high-resolution snapshots of cell morphology and the inner workings of cells, while FACS isolates thousands of cells per second using simple parameters, such as the intensity of fluorescent protein labels. Recent technologies are now combining both methods to enable the fast isolation of cells with microscopic phenotypes of interest, thereby bridging a long-standing gap in the life sciences. In this Commentary, we discuss the technical advancements made by image-enabled cell sorting and highlight novel experimental strategies in functional genomics and single-cell research.

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.

Differential regulation of mRNA stability modulates transcriptional memory and facilitates environmental adaptation

Transcriptional memory, by which cells respond faster to repeated stimuli, is key for cellular adaptation and organism survival. Chromatin organization has been shown to play a role in the faster response of primed cells. However, the contribution of post-transcriptional regulation is not yet explored. Here we perform a genome-wide screen to identify novel factors modulating transcriptional memory in S. cerevisiae in response to galactose. We find that depletion of the nuclear RNA exosome increases GAL1 expression in primed cells. Our work shows that gene-specific differences in intrinsic nuclear surveillance factor association can enhance both gene induction and repression in primed cells. Finally, we show that primed cells present altered levels of RNA degradation machinery and that both nuclear and cytoplasmic mRNA decay modulate transcriptional memory. Our results demonstrate that mRNA post-transcriptional regulation, and not only transcription regulation, should be considered when investigating gene expression memory.

Deep phenotyping of two preclinical mouse models and a cohort of RBM20 mutation carriers reveals no sex-dependent disease severity in RBM20 cardiomyopathy

RBM20 cardiomyopathy is an arrhythmogenic form of dilated cardiomyopathy caused by mutations in the splicing factor RBM20. A recent study found a more severe phenotype in male patients with RBM20 cardiomyopathy patients than in female patients. Here, we aim to determine sex differences in an animal model of RBM20 cardiomyopathy and investigate potential underlying mechanisms. In addition, we aim to determine sex and gender differences in clinical parameters in a novel RBM20 cardiomyopathy patient cohort. We characterized an Rbm20 knockout (KO) mouse model, and show that splicing of key RBM20 targets, cardiac function, and arrhythmia susceptibility do not differ between sexes. Next, we performed deep phenotyping of these mice, and show that male and female Rbm20-KO mice possess transcriptomic and phosphoproteomic differences. Hypothesizing that these differences may influence the heart's ability to compensate for stress, we exposed Rbm20-KO mice to acute catecholaminergic stimulation and again found no functional differences. We also replicate the lack of functional differences in a mouse model with the Rbm20-R636Q mutation. Lastly, we present a patient cohort of 33 RBM20 cardiomyopathy patients and show that these patients do not possess sex and gender differences in disease severity. Current mouse models of RBM20 cardiomyopathy show more pronounced changes in gene expression and phosphorylation of cardiac proteins in male mice, but no sex differences in cardiac morphology and function. Moreover, other than reported before, male RBM20 cardiomyopathy patients do not present with worse cardiac function in a patient cohort from Germany and the Netherlands.NEW & NOTEWORTHY Optimal management of the cardiac disease is increasingly personalized, partly because of differences in outcomes between sexes. RBM20 cardiomyopathy has been described to be more severe in male patients, and this carries the risk that male patients are more scrutinized in the clinic than female patients. Our findings do not support this observation and suggest that treatment should not differ between male and female RBM20 cardiomyopathy patients, but instead should focus on the underlying disease mechanism.

Genotype Complements the Phenotype: Identification of the Pathogenicity of an LMNA Splice Variant by Nanopore Long-Read Sequencing in a Large DCM Family

Dilated cardiomyopathy (DCM) is a common cause of heart failure (HF) and is of familial origin in 20−40% of cases. Genetic testing by next-generation sequencing (NGS) has yielded a definite diagnosis in many cases; however, some remain elusive. In this study, we used a combination of NGS, human-induced pluripotent-stem-cell-derived cardiomyocytes (iPSC-CMs) and nanopore long-read sequencing to identify the causal variant in a multi-generational pedigree of DCM. A four-generation family with familial DCM was investigated. Next-generation sequencing (NGS) was performed on 22 family members. Skin biopsies from two affected family members were used to generate iPSCs, which were then differentiated into iPSC-CMs. Short-read RNA sequencing was used for the evaluation of the target gene expression, and long-read RNA nanopore sequencing was used to evaluate the relevance of the splice variants. The pedigree suggested a highly penetrant, autosomal dominant mode of inheritance. The phenotype of the family was suggestive of laminopathy, but previous genetic testing using both Sanger and panel sequencing only yielded conflicting evidence for LMNA p.R644C (rs142000963), which was not fully segregated. By re-sequencing four additional affected family members, further non-coding LMNA variants could be detected: rs149339264, rs199686967, rs201379016, and rs794728589. To explore the roles of these variants, iPSC-CMs were generated. RNA sequencing showed the LMNA expression levels to be significantly lower in the iPSC-CMs of the LMNA variant carriers. We demonstrated a dysregulated sarcomeric structure and altered calcium homeostasis in the iPSC-CMs of the LMNA variant carriers. Using targeted nanopore long-read sequencing, we revealed the biological significance of the variant c.356+1G>A, which generates a novel 5′ splice site in exon 1 of the cardiac isomer of LMNA, causing a nonsense mRNA product with almost complete RNA decay and haploinsufficiency. Using novel molecular analysis and nanopore technology, we demonstrated the pathogenesis of the rs794728589 (c.356+1G>A) splice variant in LMNA. This study highlights the importance of precise diagnostics in the clinical management and workup of cardiomyopathies.

Transcriptomic diversity in human medullary thymic epithelial cells

The induction of central T cell tolerance in the thymus depends on the presentation of peripheral self-epitopes by medullary thymic epithelial cells (mTECs). This promiscuous gene expression (pGE) drives mTEC transcriptomic diversity, with non-canonical transcript initiation, alternative splicing, and expression of endogenous retroelements (EREs) representing important but incompletely understood contributors. Here we map the expression of genome-wide transcripts in immature and mature human mTECs using high-throughput 5' cap and RNA sequencing. Both mTEC populations show high splicing entropy, potentially driven by the expression of peripheral splicing factors. During mTEC maturation, rates of global transcript mis-initiation increase and EREs enriched in long terminal repeat retrotransposons are up-regulated, the latter often found in proximity to differentially expressed genes. As a resource, we provide an interactive public interface for exploring mTEC transcriptomic diversity. Our findings therefore help construct a map of transcriptomic diversity in the healthy human thymus and may ultimately facilitate the identification of those epitopes which contribute to autoimmunity and immune recognition of tumor antigens.

KIR+CD8+ T cells suppress pathogenic T cells and are active in autoimmune diseases and COVID-19

In this work, we find that CD8+ T cells expressing inhibitory killer cell immunoglobulin-like receptors (KIRs) are the human equivalent of Ly49+CD8+ regulatory T cells in mice and are increased in the blood and inflamed tissues of patients with a variety of autoimmune diseases. Moreover, these CD8+ T cells efficiently eliminated pathogenic gliadin-specific CD4+ T cells from the leukocytes of celiac disease patients in vitro. We also find elevated levels of KIR+CD8+ T cells, but not CD4+ regulatory T cells, in COVID-19 patients, correlating with disease severity and vasculitis. Selective ablation of Ly49+CD8+ T cells in virus-infected mice led to autoimmunity after infection. Our results indicate that in both species, these regulatory CD8+ T cells act specifically to suppress pathogenic T cells in autoimmune and infectious diseases.

Transcriptional neighborhoods regulate transcript isoform lengths and expression levels

Sequence features of genes and their flanking regulatory regions are determinants of RNA transcript isoform expression and have been used as context-independent plug-and-play modules in synthetic biology. However, genetic context-including the adjacent transcriptional environment-also influences transcript isoform expression levels and boundaries. We used synthetic yeast strains with stochastically repositioned genes to systematically disentangle the effects of sequence and context. Profiling 120 million full-length transcript molecules across 612 genomic perturbations, we observed sequence-independent alterations to gene expression levels and transcript isoform boundaries that were influenced by neighboring transcription. We identified features of transcriptional context that could predict these alterations and used these features to engineer a synthetic circuit where transcript length was controlled by neighboring transcription. This demonstrates how positional context can be leveraged in synthetic genome engineering.

Assembly-dependent translation of subunits 6 (Atp6) and 9 (Atp9) of ATP synthase in yeast mitochondria

The yeast mitochondrial ATP synthase is an assembly of 28 subunits of 17 types of which 3 (subunits 6, 8, and 9) are encoded by mitochondrial genes, while the 14 others have a nuclear genetic origin. Within the membrane domain (FO) of this enzyme, the subunit 6 and a ring of 10 identical subunits 9 transport protons across the mitochondrial inner membrane coupled to ATP synthesis in the extra-membrane structure (F1) of ATP synthase. As a result of their dual genetic origin, the ATP synthase subunits are synthesized in the cytosol and inside the mitochondrion. How they are produced in the proper stoichiometry from two different cellular compartments is still poorly understood. The experiments herein reported show that the rate of translation of the subunits 9 and 6 is enhanced in strains with mutations leading to specific defects in the assembly of these proteins. These translation modifications involve assembly intermediates interacting with subunits 6 and 9 within the final enzyme and cis-regulatory sequences that control gene expression in the organelle. In addition to enabling a balanced output of the ATP synthase subunits, these assembly-dependent feedback loops are presumably important to limit the accumulation of harmful assembly intermediates that have the potential to dissipate the mitochondrial membrane electrical potential and the main source of chemical energy of the cell.

High-speed fluorescence image-enabled cell sorting

Fast and selective isolation of single cells with unique spatial and morphological traits remains a technical challenge. Here, we address this by establishing high-speed image-enabled cell sorting (ICS), which records multicolor fluorescence images and sorts cells based on measurements from image data at speeds up to 15,000 events per second. We show that ICS quantifies cell morphology and localization of labeled proteins and increases the resolution of cell cycle analyses by separating mitotic stages. We combine ICS with CRISPR-pooled screens to identify regulators of the nuclear factor κB (NF-κB) pathway, enabling the completion of genome-wide image-based screens in about 9 hours of run time. By assessing complex cellular phenotypes, ICS substantially expands the phenotypic space accessible to cell-sorting applications and pooled genetic screening.

Human KIR + CD8 + T cells target pathogenic T cells in Celiac disease and are active in autoimmune diseases and COVID-19

Previous reports show that Ly49 ⁺ CD8 ⁺ T cells can suppress autoimmunity in mouse models of autoimmune diseases. Here we find a markedly increased frequency of CD8 ⁺ T cells expressing inhibitory Killer cell Immunoglobulin like Receptors (KIR), the human equivalent of the Ly49 family, in the blood and inflamed tissues of various autoimmune diseases. Moreover, KIR ⁺ CD8 ⁺ T cells can efficiently eliminate pathogenic gliadin-specific CD4 ⁺ T cells from Celiac disease (CeD) patients' leukocytes in vitro . Furthermore, we observe elevated levels of KIR ⁺ CD8 ⁺ T cells, but not CD4 ⁺ regulatory T cells, in COVID-19 and influenza-infected patients, and this correlates with disease severity and vasculitis in COVID-19. Expanded KIR ⁺ CD8 ⁺ T cells from these different diseases display shared phenotypes and similar T cell receptor sequences. These results characterize a regulatory CD8 ⁺ T cell subset in humans, broadly active in both autoimmune and infectious diseases, which we hypothesize functions to control self-reactive or otherwise pathogenic T cells.

ONE-SENTENCE SUMMARY

Here we identified KIR ⁺ CD8 ⁺ T cells as a regulatory CD8 ⁺ T cell subset in humans that suppresses self-reactive or otherwise pathogenic CD4 ⁺ T cells.

Patient-derived gene and protein expression signatures of NGLY1 deficiency

N-Glycanase 1 (NGLY1) deficiency is a rare and complex genetic disorder. Although recent studies have shed light on the molecular underpinnings of NGLY1 deficiency, a systematic characterization of gene and protein expression changes in patient-derived cells has been lacking. Here, we performed RNA-sequencing and mass spectrometry to determine the transcriptomes and proteomes of 66 cell lines representing 4 different cell types derived from 14 NGLY1 deficient patients and 17 controls. Although NGLY1 protein levels were up to 9.5-fold downregulated in patients compared to parents, residual and likely non-functional NGLY1 protein was detectable in all patient-derived lymphoblastoid cell lines. Consistent with the role of NGLY1 as a regulator of the transcription factor Nrf1, we observed a cell type-independent downregulation of proteasomal genes in NGLY1 deficient cells. In contrast, genes involved in ribosome biogenesis and mRNA processing were upregulated in multiple cell types. In addition, we observed cell type-specific effects. For example, genes and proteins involved in glutathione synthesis, such as the glutamate-cysteine ligase subunits GCLC and GCLM, were downregulated specifically in lymphoblastoid cells. We provide a web application that enables access to all results generated in this study at https://apps.embl.de/ngly1browser. This resource will guide future studies of NGLY1 deficiency in directions that are most relevant to patients.

Case-control study evaluating risk factors for SARS-CoV-2 outbreak amongst healthcare personnel at a tertiary care center

BACKGROUND

Despite several outbreaks of SARS-CoV-2 amongst healthcare personnel (HCP) exposed to COVID-19 patients globally, risk factors for transmission remain poorly understood.

METHODS

We conducted an outbreak investigation and case-control study to evaluate SARS-CoV-2 transmission risk in an outbreak among HCP at an academic medical center in California that was confirmed by whole genome sequencing.

RESULTS

A total of 7/9 cases and 93/182 controls completed a voluntary survey about risk factors. Compared to controls, cases reported significantly more patient contact time. Cases were also significantly more likely to have performed airway procedures on the index patient, particularly placing the patient on high flow nasal cannula, continuous positive airway pressure (CPAP), or bilevel positive airway pressure (BiPAP) (OR = 11.6; 95% CI = 1.7 -132.1).

DISCUSSION

This study highlights the risk of nosocomial infection of SARS-CoV-2 from patients who become infectious midway into their hospitalization. Our findings also reinforce the importance of patient contact time and aerosol-generating procedures as key risk factors for HCP infection with SARS-CoV-2.

CONCLUSIONS

Re-testing patients for SARS-CoV-2 after admission in suspicious cases and using N95 masks for all aerosol-generating procedures regardless of initial patient SARS-CoV-2 test results can help reduce the risk of SARS-COV-2 transmission to HCP.

High-Throughput Nucleotide Resolution Predictions of Assay Limitations Increase the Reliability and Concordance of Clinical Tests

PURPOSE

The ability of next-generation sequencing (NGS) assays to interrogate thousands of genomic loci has revolutionized genetic testing. However, translation to the clinic is impeded by false-negative results that pose a risk to patients. In response, regulatory bodies are calling for reliability measures to be reported alongside NGS results. Existing methods to estimate reliability do not account for sample- and position-specific variability, which can be significant. Here, we report an approach that computes reliability metrics for every genomic position and sample interrogated by an NGS assay.

METHODS

Our approach predicts the limit of detection (LOD), the lowest reliably detectable variant fraction, by taking technical factors into account. We initially explored how LOD is affected by input material amount, library conversion rate, sequencing coverage, and sequencing error rate. This revealed that LOD depends heavily on genomic context and sample properties. Using these insights, we developed a computational approach to predict LOD on the basis of a biophysical model of the NGS workflow. We focused on targeted assays for cell-free DNA, but, in principle, this approach applies to any NGS assay.

RESULTS

We validated our approach by showing that it accurately predicts LOD and distinguishes reliable from unreliable results when screening 580 lung cancer samples for actionable mutations. Compared with a standard variant calling workflow, our approach avoided most false negatives and improved interassay concordance from 94% to 99%.

CONCLUSION

Our approach, which we name LAVA (LOD-aware variant analysis), reports the LOD for every position and sample interrogated by an NGS assay. This enables reliable results to be identified and improves the transparency and safety of genetic tests.

Integrated single-cell transcriptomics and epigenomics reveals strong germinal center-associated etiology of autoimmune risk loci

[Figure: see text].

Regulation of Msh4-Msh5 association with meiotic chromosomes in budding yeast

In the baker's yeast Saccharomyces cerevisiae, most of the meiotic crossovers are generated through a pathway involving the highly conserved mismatch repair related Msh4-Msh5 complex. To understand the role of Msh4-Msh5 in meiotic crossing over, we determined its genome wide in vivo binding sites in meiotic cells. We show that Msh5 specifically associates with DSB hotspots, chromosome axes, and centromeres on chromosomes. A basal level of Msh5 association with these chromosomal features is observed even in the absence of DSB formation (spo11Δ mutant) at the early stages of meiosis. But efficient binding to DSB hotspots and chromosome axes requires DSB formation and resection and is enhanced by double Holliday junction structures. Msh5 binding is also correlated to DSB frequency and enhanced on small chromosomes with higher DSB and crossover density. The axis protein Red1 is required for Msh5 association with the chromosome axes and DSB hotspots but not centromeres. Although binding sites of Msh5 and other pro-crossover factors like Zip3 show extensive overlap, Msh5 associates with centromeres independent of Zip3. These results on Msh5 localization in wild type and meiotic mutants have implications for how Msh4-Msh5 works with other pro-crossover factors to ensure crossover formation.

Single-molecule, full-length transcript isoform sequencing reveals disease-associated RNA isoforms in cardiomyocytes

Alternative splicing generates differing RNA isoforms that govern phenotypic complexity of eukaryotes. Its malfunction underlies many diseases, including cancer and cardiovascular diseases. Comparative analysis of RNA isoforms at the genome-wide scale has been difficult. Here, we establish an experimental and computational pipeline that performs de novo transcript annotation and accurately quantifies transcript isoforms from cDNA sequences with a full-length isoform detection accuracy of 97.6%. We generate a searchable, quantitative human transcriptome annotation with 31,025 known and 5,740 novel transcript isoforms ( http://steinmetzlab.embl.de/iBrowser/ ). By analyzing the isoforms in the presence of RNA Binding Motif Protein 20 (RBM20) mutations associated with aggressive dilated cardiomyopathy (DCM), we identify 121 differentially expressed transcript isoforms in 107 cardiac genes. Our approach enables quantitative dissection of complex transcript architecture instead of mere identification of inclusion or exclusion of individual exons, as exemplified by the discovery of IMMT isoforms mis-spliced by RBM20 mutations. Thereby we achieve a path to direct differential expression testing independent of an existing annotation of transcript isoforms, providing more immediate biological interpretation and higher resolution transcriptome comparisons.

Fast and inexpensive whole-genome sequencing library preparation from intact yeast cells

Through the increase in the capacity of sequencing machines massively parallel sequencing of thousands of samples in a single run is now possible. With the improved throughput and resulting drop in the price of sequencing, the cost and time for preparation of sequencing libraries have become the major bottleneck in large-scale experiments. Methods using a hyperactive variant of the Tn5 transposase efficiently generate libraries starting from cDNA or genomic DNA in a few hours and are highly scalable. For genome sequencing, however, the time and effort spent on genomic DNA isolation limit the practicability of sequencing large numbers of samples. Here, we describe a highly scalable method for preparing high-quality whole-genome sequencing libraries directly from Saccharomyces cerevisiae cultures in less than 3 h at 34 cents per sample. We skip the rate-limiting step of genomic DNA extraction by directly tagmenting lysed yeast spheroplasts and add a nucleosome release step prior to enrichment PCR to improve the evenness of genomic coverage. Resulting libraries do not show any GC bias and are comparable in quality to libraries processed from genomic DNA with a commercially available Tn5-based kit. We use our protocol to investigate CRISPR/Cas9 on- and off-target edits and reliably detect edited variants and shared polymorphisms between strains. Our protocol enables rapid preparation of unbiased and high-quality, sequencing-ready indexed libraries for hundreds of yeast strains in a single day at a low price. By adjusting individual steps of our workflow, we expect that our protocol can be adapted to other organisms.

Recommendations for accurate genotyping of SARS-CoV-2 using amplicon-based sequencing of clinical samples

OBJECTIVES

Genotyping of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has been instrumental in monitoring viral evolution and transmission during the pandemic. The quality of the sequence data obtained from these genotyping efforts depends on several factors, including the quantity/integrity of the input material, the technology, and laboratory-specific implementation. The current lack of guidelines for SARS-CoV-2 genotyping leads to inclusion of error-containing genome sequences in genomic epidemiology studies. We aimed to establish clear and broadly applicable recommendations for reliable virus genotyping.

METHODS

We established and used a sequencing data analysis workflow that reliably identifies and removes technical artefacts; such artefacts can result in miscalls when using alternative pipelines to process clinical samples and synthetic viral genomes with an amplicon-based genotyping approach. We evaluated the impact of experimental factors, including viral load and sequencing depth, on correct sequence determination.

RESULTS

We found that at least 1000 viral genomes are necessary to confidently detect variants in the SARS-CoV-2 genome at frequencies of ≥10%. The broad applicability of our recommendations was validated in over 200 clinical samples from six independent laboratories. The genotypes we determined for clinical isolates with sufficient quality cluster by sampling location and period. Our analysis also supports the rise in frequencies of 20A.EU1 and 20A.EU2, two recently reported European strains whose dissemination was facilitated by travel during the summer of 2020.

CONCLUSIONS

We present much-needed recommendations for the reliable determination of SARS-CoV-2 genome sequences and demonstrate their broad applicability in a large cohort of clinical samples.

iPSC Modeling of RBM20-Deficient DCM Identifies Upregulation of RBM20 as a Therapeutic Strategy

Recent advances in induced pluripotent stem cell (iPSC) technology and directed differentiation of iPSCs into cardiomyocytes (iPSC-CMs) make it possible to model genetic heart disease in vitro. We apply CRISPR/Cas9 genome editing technology to introduce three RBM20 mutations in iPSCs and differentiate them into iPSC-CMs to establish an in vitro model of RBM20 mutant dilated cardiomyopathy (DCM). In iPSC-CMs harboring a known causal RBM20 variant, the splicing of RBM20 target genes, calcium handling, and contractility are impaired consistent with the disease manifestation in patients. A variant (Pro633Leu) identified by exome sequencing of patient genomes displays the same disease phenotypes, thus establishing this variant as disease causing. We find that all-trans retinoic acid upregulates RBM20 expression and reverts the splicing, calcium handling, and contractility defects in iPSC-CMs with different causal RBM20 mutations. These results suggest that pharmacological upregulation of RBM20 expression is a promising therapeutic strategy for DCM patients with a heterozygous mutation in RBM20.

Briganti et al. use iPSC and CRISPR/Cas9 to create a model of RBM20-deficient dilated cardiomyopathy (DCM) that recapitulates mRNA splicing and contractile defects of the disease. They evaluate pharmacological upregulation of RBM20 as a therapeutic strategy. All-trans retinoic acid upregulates RBM20 expression and ameliorates the in vitro hallmarks of disease.

The chaperone-binding activity of the mitochondrial surface receptor Tom70 protects the cytosol against mitoprotein-induced stress

Most mitochondrial proteins are synthesized as precursors in the cytosol and post-translationally transported into mitochondria. The mitochondrial surface protein Tom70 acts at the interface of the cytosol and mitochondria. In vitro import experiments identified Tom70 as targeting receptor, particularly for hydrophobic carriers. Using in vivo methods and high-content screens, we revisit the question of Tom70 function and considerably expand the set of Tom70-dependent mitochondrial proteins. We demonstrate that the crucial activity of Tom70 is its ability to recruit cytosolic chaperones to the outer membrane. Indeed, tethering an unrelated chaperone-binding domain onto the mitochondrial surface complements most of the defects caused by Tom70 deletion. Tom70-mediated chaperone recruitment reduces the proteotoxicity of mitochondrial precursor proteins, particularly of hydrophobic inner membrane proteins. Thus, our work suggests that the predominant function of Tom70 is to tether cytosolic chaperones to the outer mitochondrial membrane, rather than to serve as a mitochondrion-specifying targeting receptor.

Single-cell analyses reveal SARS-CoV-2 interference with intrinsic immune response in the human gut

Exacerbated pro-inflammatory immune response contributes to COVID-19 pathology. However, despite the mounting evidence about SARS-CoV-2 infecting the human gut, little is known about the antiviral programs triggered in this organ. To address this gap, we performed single-cell transcriptomics of SARS-CoV-2-infected intestinal organoids. We identified a subpopulation of enterocytes as the prime target of SARS-CoV-2 and, interestingly, found the lack of positive correlation between susceptibility to infection and the expression of ACE2. Infected cells activated strong pro-inflammatory programs and produced interferon, while expression of interferon-stimulated genes was limited to bystander cells due to SARS-CoV-2 suppressing the autocrine action of interferon. These findings reveal that SARS-CoV-2 curtails the immune response and highlights the gut as a pro-inflammatory reservoir that should be considered to fully understand SARS-CoV-2 pathogenesis.

Identification of leukemic and pre-leukemic stem cells by clonal tracking from single-cell transcriptomics

Cancer stem cells drive disease progression and relapse in many types of cancer. Despite this, a thorough characterization of these cells remains elusive and with it the ability to eradicate cancer at its source. In acute myeloid leukemia (AML), leukemic stem cells (LSCs) underlie mortality but are difficult to isolate due to their low abundance and high similarity to healthy hematopoietic stem cells (HSCs). Here, we demonstrate that LSCs, HSCs, and pre-leukemic stem cells can be identified and molecularly profiled by combining single-cell transcriptomics with lineage tracing using both nuclear and mitochondrial somatic variants. While mutational status discriminates between healthy and cancerous cells, gene expression distinguishes stem cells and progenitor cell populations. Our approach enables the identification of LSC-specific gene expression programs and the characterization of differentiation blocks induced by leukemic mutations. Taken together, we demonstrate the power of single-cell multi-omic approaches in characterizing cancer stem cells.

CRISPRi screens reveal genes modulating yeast growth in lignocellulose hydrolysate

BACKGROUND

Baker's yeast is a widely used eukaryotic cell factory, producing a diverse range of compounds including biofuels and fine chemicals. The use of lignocellulose as feedstock offers the opportunity to run these processes in an environmentally sustainable way. However, the required hydrolysis pretreatment of lignocellulosic material releases toxic compounds that hamper yeast growth and consequently productivity.

RESULTS

Here, we employ CRISPR interference in S. cerevisiae to identify genes modulating fermentative growth in plant hydrolysate and in presence of lignocellulosic toxins. We find that at least one-third of hydrolysate-associated gene functions are explained by effects of known toxic compounds, such as the decreased growth of YAP1 or HAA1, or increased growth of DOT6 knock-down strains in hydrolysate.

CONCLUSION

Our study confirms previously known genetic elements and uncovers new targets towards designing more robust yeast strains for the utilization of lignocellulose hydrolysate as sustainable feedstock, and, more broadly, paves the way for applying CRISPRi screens to improve industrial fermentation processes.

Combined transient ablation and single-cell RNA-sequencing reveals the development of medullary thymic epithelial cells

Medullary thymic epithelial cells (mTECs) play a critical role in central immune tolerance by mediating negative selection of autoreactive T cells through the collective expression of the peripheral self-antigen compartment, including tissue-specific antigens (TSAs). Recent work has shown that gene-expression patterns within the mTEC compartment are heterogenous and include multiple differentiated cell states. To further define mTEC development and medullary epithelial lineage relationships, we combined lineage tracing and recovery from transient in vivo mTEC ablation with single-cell RNA-sequencing in Mus musculus. The combination of bioinformatic and experimental approaches revealed a non-stem transit-amplifying population of cycling mTECs that preceded Aire expression. We propose a branching model of mTEC development wherein a heterogeneous pool of transit-amplifying cells gives rise to Aire- and Ccl21a-expressing mTEC subsets. We further use experimental techniques to show that within the Aire-expressing developmental branch, TSA expression peaked as Aire expression decreased, implying Aire expression must be established before TSA expression can occur. Collectively, these data provide a roadmap of mTEC development and demonstrate the power of combinatorial approaches leveraging both in vivo models and high-dimensional datasets.

Specialized cells in the immune system known as T cells protect the body from infection by destroying disease-causing microbes, such as bacteria or viruses. T cells use proteins on their surface called receptors to stick to infectious microbes and remove them from the body. Some newly developed T-cells, however, contain receptors that recognize and bind to cells that belong in the body. If these faulty T cells are released, they can attack healthy tissues and cause an autoimmune disease.

After a new T cell is developed, it gets carried to a gland in the chest known as the thymus. Cells in the thymus called mTECs screen T cells for receptors that may bind to the body’s tissues. mTECs do this by presenting T cells with proteins that are commonly found on the surface of healthy cells in the body. If a T cell recognizes any of these ‘tissue specific proteins’, it is destroyed or given a new role in the body. Some faulty T cells, however, still manage to evade detection. One way to uncover why this might happen is to investigate how mTECs develop. Previous work showed that mTECs transition through various stages before reaching their final form. However, the order in which these events occur remained unclear.

To gain a better understanding of these developmental steps, Wells, Miller et al. extracted mTECs from the thymus of mice and analyzed the genetic make-up of individual cells. This uncovered a missing link in mTEC development: a new type of cell that is the immediate predecessor of the final mTEC. These ‘predecessor’ cells were actively growing, highlighting that mTECs can be constantly generated in the body. By probing the genes that generate tissue-specific proteins in mTECs, Wells, Miller et al. revealed that these proteins were only produced for short periods and in the late stages of mTEC development.

These findings contribute to our understanding of how mTECs develop to screen T cells. Mapping these developmental stages will make it easier to identify when faulty T cells are able to evade mTECs. This will lead to earlier detection of autoimmune diseases which could result in better treatments.

Rpb4 and Puf3 imprint and post-transcriptionally control the stability of a common set of mRNAs in yeast

Gene expression involving RNA polymerase II is regulated by the concerted interplay between mRNA synthesis and degradation, crosstalk in which mRNA decay machinery and transcription machinery respectively impact transcription and mRNA stability. Rpb4, and likely dimer Rpb4/7, seem the central components of the RNA pol II governing these processes. In this work we unravel the molecular mechanisms participated by Rpb4 that mediate the posttranscriptional events regulating mRNA imprinting and stability. By RIP-Seq, we analysed genome-wide the association of Rpb4 with mRNAs and demonstrated that it targeted a large population of more than 1400 transcripts. A group of these mRNAs was also the target of the RNA binding protein, Puf3. We demonstrated that Rpb4 and Puf3 physically, genetically, and functionally interact and also affect mRNA stability, and likely the imprinting, of a common group of mRNAs. Furthermore, the Rpb4 and Puf3 association with mRNAs depends on one another. We also demonstrated, for the first time, that Puf3 associates with chromatin in an Rpb4-dependent manner. Our data also suggest that Rpb4 could be a key element of the RNA pol II that coordinates mRNA synthesis, imprinting and stability in cooperation with RBPs.

TIF-Seq2 disentangles overlapping isoforms in complex human transcriptomes

Eukaryotic transcriptomes are complex, involving thousands of overlapping transcripts. The interleaved nature of the transcriptomes limits our ability to identify regulatory regions, and in some cases can lead to misinterpretation of gene expression. To improve the understanding of the overlapping transcriptomes, we have developed an optimized method, TIF-Seq2, able to sequence simultaneously the 5' and 3' ends of individual RNA molecules at single-nucleotide resolution. We investigated the transcriptome of a well characterized human cell line (K562) and identified thousands of unannotated transcript isoforms. By focusing on transcripts which are challenging to be investigated with RNA-Seq, we accurately defined boundaries of lowly expressed unannotated and read-through transcripts putatively encoding fusion genes. We validated our results by targeted long-read sequencing and standard RNA-Seq for chronic myeloid leukaemia patient samples. Taking the advantage of TIF-Seq2, we explored transcription regulation among overlapping units and investigated their crosstalk. We show that most overlapping upstream transcripts use poly(A) sites within the first 2 kb of the downstream transcription units. Our work shows that, by paring the 5' and 3' end of each RNA, TIF-Seq2 can improve the annotation of complex genomes, facilitate accurate assignment of promoters to genes and easily identify transcriptionally fused genes.

GOTI, a method to identify genome-wide off-target effects of genome editing in mouse embryos

Genome editing holds great potential for correcting pathogenic mutations. We developed a method called GOTI (genome-wide off-target analysis by two-cell embryo injection) to detect off-target mutations by editing one blastomere of two-cell mouse embryos using either CRISPR-Cas9 or base editors. GOTI directly compares edited and non-edited cells without the interference of genetic background and thus could detect potential off-target variants with high sensitivity. Notably, the GOTI method was designed to detect potential off-target variants of any genome editing tools by the combination of experimental and computational approaches, which is critical for accurate evaluation of the safety of genome editing tools. Here we provide a detailed protocol for GOTI, including mice mating, two-cell embryo injection, embryonic day 14.5 embryo digestion, fluorescence-activated cell sorting, whole-genome sequencing and data analysis. To enhance the utility of GOTI, we also include a computational workflow called GOTI-seq (https://github.com/sydaileen/GOTI-seq) for the sequencing data analysis, which can generate the final genome-wide off-target variants from raw sequencing data directly. The protocol typically takes 20 d from the mice mating to sequencing and 7 d for sequencing data analysis.

Targeted Perturb-seq enables genome-scale genetic screens in single cells

The transcriptome contains rich information on molecular, cellular and organismal phenotypes. However, experimental and statistical limitations constrain sensitivity and throughput of genetic screening with single-cell transcriptomics readout. To overcome these limitations, we introduce targeted Perturb-seq (TAP-seq), a sensitive, inexpensive and platform-independent method focusing single-cell RNA-seq coverage on genes of interest, thereby increasing the sensitivity and scale of genetic screens by orders of magnitude. TAP-seq permits routine analysis of thousands of CRISPR-mediated perturbations within a single experiment, detects weak effects and lowly expressed genes, and decreases sequencing requirements by up to 50-fold. We apply TAP-seq to generate perturbation-based enhancer-target gene maps for 1,778 enhancers within 2.5% of the human genome. We thereby show that enhancer-target association is jointly determined by three-dimensional contact frequency and epigenetic states, allowing accurate prediction of enhancer targets throughout the genome. In addition, we demonstrate that TAP-seq can identify cell subtypes with only 100 sequencing reads per cell.

Targeted Perturb-seq enables genome-scale genetic screens in single cells

The transcriptome contains rich information on molecular, cellular and organismal phenotypes. However, experimental and statistical limitations constrain sensitivity and throughput of genetic screening with single-cell transcriptomics readout. To overcome these limitations, we introduce targeted Perturb-seq (TAP-seq), a sensitive, inexpensive and platform-independent method focusing single-cell RNA-seq coverage on genes of interest, thereby increasing the sensitivity and scale of genetic screens by orders of magnitude. TAP-seq permits routine analysis of thousands of CRISPR-mediated perturbations within a single experiment, detects weak effects and lowly expressed genes, and decreases sequencing requirements by up to 50-fold. We apply TAP-seq to generate perturbation-based enhancer-target gene maps for 1,778 enhancers within 2.5% of the human genome. We thereby show that enhancer-target association is jointly determined by three-dimensional contact frequency and epigenetic states, allowing accurate prediction of enhancer targets throughout the genome. In addition, we demonstrate that TAP-seq can identify cell subtypes with only 100 sequencing reads per cell.

Loss of N-Glycanase 1 Alters Transcriptional and Translational Regulation in K562 Cell Lines

N-Glycanase 1 (NGLY1) deficiency is an ultra-rare, complex and devastating neuromuscular disease. Patients display multi-organ symptoms including developmental delays, movement disorders, seizures, constipation and lack of tear production. NGLY1 is a deglycosylating protein involved in the degradation of misfolded proteins retrotranslocated from the endoplasmic reticulum (ER). NGLY1-deficient cells have been reported to exhibit decreased deglycosylation activity and an increased sensitivity to proteasome inhibitors. We show that the loss of NGLY1 causes substantial changes in the RNA and protein landscape of K562 cells and results in downregulation of proteasomal subunits, consistent with its processing of the transcription factor NFE2L1. We employed the CMap database to predict compounds that can modulate NGLY1 activity. Utilizing our robust K562 screening system, we demonstrate that the compound NVP-BEZ235 (Dactosilib) promotes degradation of NGLY1-dependent substrates, concurrent with increased autophagic flux, suggesting that stimulating autophagy may assist in clearing aberrant substrates during NGLY1 deficiency.