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.

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.

Widespread Co-translational RNA Decay Reveals Ribosome Dynamics.

It is generally assumed that mRNAs undergoing translation are protected from decay. Here, we show that mRNAs are, in fact, co-translationally degraded. This is a widespread and conserved process affecting most genes, where 5'-3' transcript degradation follows the last translating ribosome, producing an in vivo ribosomal footprint. By sequencing the ends of 5' phosphorylated mRNA degradation intermediates, we obtain a genome-wide drug-free measurement of ribosome dynamics. We identify general translation termination pauses in both normal and stress conditions. In addition, we describe novel codon-specific ribosomal pausing sites in response to oxidative stress that are dependent on the RNase Rny1. Our approach is simple and straightforward and does not require the use of translational inhibitors or in vitro RNA footprinting that can alter ribosome protection patterns.

Genetic Control of Chromatin States in Humans Involves Local and Distal Chromosomal Interactions.

Deciphering the impact of genetic variants on gene regulation is fundamental to understanding human disease. Although gene regulation often involves long-range interactions, it is unknown to what extent non-coding genetic variants influence distal molecular phenotypes. Here, we integrate chromatin profiling for three histone marks in lymphoblastoid cell lines (LCLs) from 75 sequenced individuals with LCL-specific Hi-C and ChIA-PET-based chromatin contact maps to uncover one of the largest collections of local and distal histone quantitative trait loci (hQTLs). Distal QTLs are enriched within topologically associated domains and exhibit largely concordant variation of chromatin state coordinated by proximal and distal non-coding genetic variants. Histone QTLs are enriched for common variants associated with autoimmune diseases and enable identification of putative target genes of disease-associated variants from genome-wide association studies. These analyses provide insights into how genetic variation can affect human disease phenotypes by coordinated changes in chromatin at interacting regulatory elements.

Single-cell transcriptome analysis reveals coordinated ectopic gene-expression patterns in medullary thymic epithelial cells.

Expression of tissue-restricted self antigens (TRAs) in medullary thymic epithelial cells (mTECs) is essential for the induction of self-tolerance and prevents autoimmunity, with each TRA being expressed in only a few mTECs. How this process is regulated in single mTECs and is coordinated at the population level, such that the varied single-cell patterns add up to faithfully represent TRAs, is poorly understood. Here we used single-cell RNA sequencing and obtained evidence of numerous recurring TRA-co-expression patterns, each present in only a subset of mTECs. Co-expressed genes clustered in the genome and showed enhanced chromatin accessibility. Our findings characterize TRA expression in mTECs as a coordinated process that might involve local remodeling of chromatin and thus ensures a comprehensive representation of the immunological self.

Set3 HDAC mediates effects of overlapping noncoding transcription on gene induction kinetics.

The Set3 histone deacetylase complex (Set3C) binds histone H3 dimethylated at lysine 4 (H3K4me2) to mediate deacetylation of histones in 5'-transcribed regions. To discern how Set3C affects gene expression, genome-wide transcription was analyzed in yeast undergoing a series of carbon source shifts. Deleting SET3 primarily caused changes during transition periods, as genes were induced or repressed. Surprisingly, a majority of Set3-affected genes are overlapped by noncoding RNA (ncRNA) transcription. Many Set3-repressed genes have H3K4me2 instead of me3 over promoter regions, due to either reduced H3K4me3 or ncRNA transcription from distal or antisense promoters. Set3C also represses internal cryptic promoters, but in different regions of genes than the Set2/Rpd3S pathway. Finally, Set3C stimulates some genes by repressing an overlapping antagonistic antisense transcript. These results show that overlapping noncoding transcription can fine-tune gene expression, not via the ncRNA but by depositing H3K4me2 to recruit the Set3C deacetylase.

Insights into RNA biology from an atlas of mammalian mRNA-binding proteins.

RNA-binding proteins (RBPs) determine RNA fate from synthesis to decay. Employing two complementary protocols for covalent UV crosslinking of RBPs to RNA, we describe a systematic, unbiased, and comprehensive approach, termed "interactome capture," to define the mRNA interactome of proliferating human HeLa cells. We identify 860 proteins that qualify as RBPs by biochemical and statistical criteria, adding more than 300 RBPs to those previously known and shedding light on RBPs in disease, RNA-binding enzymes of intermediary metabolism, RNA-binding kinases, and RNA-binding architectures. Unexpectedly, we find that many proteins of the HeLa mRNA interactome are highly intrinsically disordered and enriched in short repetitive amino acid motifs. Interactome capture is broadly applicable to study mRNA interactome composition and dynamics in varied biological settings.

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.

Opposing T cell responses in experimental autoimmune encephalomyelitis.

Experimental autoimmune encephalomyelitis is a model for multiple sclerosis. Here we show that induction generates successive waves of clonally expanded CD4+, CD8+ and γδ+ T cells in the blood and central nervous system, similar to gluten-challenge studies of patients with coeliac disease. We also find major expansions of CD8+ T cells in patients with multiple sclerosis. In autoimmune encephalomyelitis, we find that most expanded CD4+ T cells are specific for the inducing myelin peptide MOG35-55. By contrast, surrogate peptides derived from a yeast peptide major histocompatibility complex library of some of the clonally expanded CD8+ T cells inhibit disease by suppressing the proliferation of MOG-specific CD4+ T cells. These results suggest that the induction of autoreactive CD4+ T cells triggers an opposing mobilization of regulatory CD8+ T cells.

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.

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.

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.

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.

Thymic tuft cells promote an IL-4-enriched medulla and shape thymocyte development.

The thymus is responsible for generating a diverse yet self-tolerant pool of T cells1. Although the thymic medulla consists mostly of developing and mature AIRE+ epithelial cells, recent evidence has suggested that there is far greater heterogeneity among medullary thymic epithelial cells than was previously thought2. Here we describe in detail an epithelial subset that is remarkably similar to peripheral tuft cells that are found at mucosal barriers3. Similar to the periphery, thymic tuft cells express the canonical taste transduction pathway and IL-25. However, they are unique in their spatial association with cornified aggregates, ability to present antigens and expression of a broad diversity of taste receptors. Some thymic tuft cells pass through an Aire-expressing stage and depend on a known AIRE-binding partner, HIPK2, for their development. Notably, the taste chemosensory protein TRPM5 is required for their thymic function through which they support the development and polarization of thymic invariant natural killer T cells and act to establish a medullary microenvironment that is enriched in the type 2 cytokine, IL-4. These findings indicate that there is a compartmentalized medullary environment in which differentiation of a minor and highly specialized epithelial subset has a non-redundant role in shaping thymic function.

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.

Assembly of functionally integrated human forebrain spheroids.

The development of the nervous system involves a coordinated succession of events including the migration of GABAergic (γ-aminobutyric-acid-releasing) neurons from ventral to dorsal forebrain and their integration into cortical circuits. However, these interregional interactions have not yet been modelled with human cells. Here we generate three-dimensional spheroids from human pluripotent stem cells that resemble either the dorsal or ventral forebrain and contain cortical glutamatergic or GABAergic neurons. These subdomain-specific forebrain spheroids can be assembled in vitro to recapitulate the saltatory migration of interneurons observed in the fetal forebrain. Using this system, we find that in Timothy syndrome-a neurodevelopmental disorder that is caused by mutations in the CaV1.2 calcium channel-interneurons display abnormal migratory saltations. We also show that after migration, interneurons functionally integrate with glutamatergic neurons to form a microphysiological system. We anticipate that this approach will be useful for studying neural development and disease, and for deriving spheroids that resemble other brain regions to assemble circuits in vitro.

Myt1l safeguards neuronal identity by actively repressing many non-neuronal fates.

Normal differentiation and induced reprogramming require the activation of target cell programs and silencing of donor cell programs. In reprogramming, the same factors are often used to reprogram many different donor cell types. As most developmental repressors, such as RE1-silencing transcription factor (REST) and Groucho (also known as TLE), are considered lineage-specific repressors, it remains unclear how identical combinations of transcription factors can silence so many different donor programs. Distinct lineage repressors would have to be induced in different donor cell types. Here, by studying the reprogramming of mouse fibroblasts to neurons, we found that the pan neuron-specific transcription factor Myt1-like (Myt1l) exerts its pro-neuronal function by direct repression of many different somatic lineage programs except the neuronal program. The repressive function of Myt1l is mediated via recruitment of a complex containing Sin3b by binding to a previously uncharacterized N-terminal domain. In agreement with its repressive function, the genomic binding sites of Myt1l are similar in neurons and fibroblasts and are preferentially in an open chromatin configuration. The Notch signalling pathway is repressed by Myt1l through silencing of several members, including Hes1. Acute knockdown of Myt1l in the developing mouse brain mimicked a Notch gain-of-function phenotype, suggesting that Myt1l allows newborn neurons to escape Notch activation during normal development. Depletion of Myt1l in primary postmitotic neurons de-repressed non-neuronal programs and impaired neuronal gene expression and function, indicating that many somatic lineage programs are actively and persistently repressed by Myt1l to maintain neuronal identity. It is now tempting to speculate that similar 'many-but-one' lineage repressors exist for other cell fates; such repressors, in combination with lineage-specific activators, would be prime candidates for use in reprogramming additional cell types.

A functional connection between translation elongation and protein folding at the ribosome exit tunnel in Saccharomyces cerevisiae.

Proteostasis needs to be tightly controlled to meet the cellular demand for correctly de novo folded proteins and to avoid protein aggregation. While a coupling between translation rate and co-translational folding, likely involving an interplay between the ribosome and its associated chaperones, clearly appears to exist, the underlying mechanisms and the contribution of ribosomal proteins remain to be explored. The ribosomal protein uL3 contains a long internal loop whose tip region is in close proximity to the ribosomal peptidyl transferase center. Intriguingly, the rpl3[W255C] allele, in which the residue making the closest contact to this catalytic site is mutated, affects diverse aspects of ribosome biogenesis and function. Here, we have uncovered, by performing a synthetic lethal screen with this allele, an unexpected link between translation and the folding of nascent proteins by the ribosome-associated Ssb-RAC chaperone system. Our results reveal that uL3 and Ssb-RAC cooperate to prevent 80S ribosomes from piling up within the 5' region of mRNAs early on during translation elongation. Together, our study provides compelling in vivo evidence for a functional connection between peptide bond formation at the peptidyl transferase center and chaperone-assisted de novo folding of nascent polypeptides at the solvent-side of the peptide exit tunnel.