Engineering RNA export for measurement and manipulation of living cells.

A system for programmable export of RNA molecules from living cells would enable both non-destructive monitoring of cell dynamics and engineering of cells capable of delivering executable RNA programs to other cells. We developed genetically encoded cellular RNA exporters, inspired by viruses, that efficiently package and secrete cargo RNA molecules from mammalian cells within protective nanoparticles. Exporting and sequencing RNA barcodes enabled non-destructive monitoring of cell population dynamics with clonal resolution. Further, by incorporating fusogens into the nanoparticles, we demonstrated the delivery, expression, and functional activity of exported mRNA in recipient cells. We term these systems COURIER (controlled output and uptake of RNA for interrogation, expression, and regulation). COURIER enables measurement of cell dynamics and establishes a foundation for hybrid cell and gene therapies based on cell-to-cell delivery of RNA.

Cryo-EM structure of gas vesicles for buoyancy-controlled motility.

Gas vesicles are gas-filled nanocompartments that allow a diverse group of bacteria and archaea to control their buoyancy. The molecular basis of their properties and assembly remains unclear. Here, we report the 3.2 Å cryo-EM structure of the gas vesicle shell made from the structural protein GvpA that self-assembles into hollow helical cylinders closed off by cone-shaped tips. Two helical half shells connect through a characteristic arrangement of GvpA monomers, suggesting a mechanism of gas vesicle biogenesis. The fold of GvpA features a corrugated wall structure typical for force-bearing thin-walled cylinders. Small pores enable gas molecules to diffuse across the shell, while the exceptionally hydrophobic interior surface effectively repels water. Comparative structural analysis confirms the evolutionary conservation of gas vesicle assemblies and demonstrates molecular features of shell reinforcement by GvpC. Our findings will further research into gas vesicle biology and facilitate molecular engineering of gas vesicles for ultrasound imaging.

Massively parallel functional dissection of schizophrenia-associated noncoding genetic variants.

Schizophrenia (SCZ) is a highly heritable mental disorder with thousands of associated genetic variants located mostly in the noncoding space of the genome. Translating these associations into insights regarding the underlying pathomechanisms has been challenging because the causal variants, their mechanisms of action, and their target genes remain largely unknown. We implemented a massively parallel variant annotation pipeline (MVAP) to perform SCZ variant-to-function mapping at scale in disease-relevant neural cell types. This approach identified 620 functional variants (1.7%) that operate in a highly developmental context and neuronal-activity-dependent manner. Multimodal integration of epigenomic and CRISPRi screening data enabled us to link these functional variants to target genes, biological processes, and ultimately alterations of neuronal physiology. These results provide a multistage prioritization strategy to map functional single-nucleotide polymorphism (SNP)-to-gene-to-endophenotype relations and offer biological insights into the context-dependent molecular processes modulated by SCZ-associated genetic variation.

Functional viromic screens uncover regulatory RNA elements.

The number of sequenced viral genomes has surged recently, presenting an opportunity to understand viral diversity and uncover unknown regulatory mechanisms. Here, we conducted a screening of 30,367 viral segments from 143 species representing 96 genera and 37 families. Using a library of viral segments in 3' UTR, we identified hundreds of elements impacting RNA abundance, translation, and nucleocytoplasmic distribution. To illustrate the power of this approach, we investigated K5, an element conserved in kobuviruses, and found its potent ability to enhance mRNA stability and translation in various contexts, including adeno-associated viral vectors and synthetic mRNAs. Moreover, we identified a previously uncharacterized protein, ZCCHC2, as a critical host factor for K5. ZCCHC2 recruits the terminal nucleotidyl transferase TENT4 to elongate poly(A) tails with mixed sequences, delaying deadenylation. This study provides a unique resource for virus and RNA research and highlights the potential of the virosphere for biological discoveries.

Long Noncoding RNAs: Molecular Modalities to Organismal Functions.

We have known for decades that long noncoding RNAs (lncRNAs) can play essential functions across most forms of life. The maintenance of chromosome length requires an lncRNA (e.g., hTERC) and two lncRNAs in the ribosome that are required for protein synthesis. Thus, lncRNAs can represent powerful RNA machines. More recently, it has become clear that mammalian genomes encode thousands more lncRNAs. Thus, we raise the question: Which, if any, of these lncRNAs could also represent RNA-based machines? Here we synthesize studies that are beginning to address this question by investigating fundamental properties of lncRNA genes, revealing new insights into the RNA structure-function relationship, determining cis- and trans-acting lncRNAs in vivo, and generating new developments in high-throughput screening used to identify functional lncRNAs. Overall, these findings provide a context toward understanding the molecular grammar underlying lncRNA biology.

Snake Venom Gland Organoids.

Wnt dependency and Lgr5 expression define multiple mammalian epithelial stem cell types. Under defined growth factor conditions, such adult stem cells (ASCs) grow as 3D organoids that recapitulate essential features of the pertinent epithelium. Here, we establish long-term expanding venom gland organoids from several snake species. The newly assembled transcriptome of the Cape coral snake reveals that organoids express high levels of toxin transcripts. Single-cell RNA sequencing of both organoids and primary tissue identifies distinct venom-expressing cell types as well as proliferative cells expressing homologs of known mammalian stem cell markers. A hard-wired regional heterogeneity in the expression of individual venom components is maintained in organoid cultures. Harvested venom peptides reflect crude venom composition and display biological activity. This study extends organoid technology to reptilian tissues and describes an experimentally tractable model system representing the snake venom gland.

Spatial Multiplexing of Fluorescent Reporters for Imaging Signaling Network Dynamics.

In order to analyze how a signal transduction network converts cellular inputs into cellular outputs, ideally one would measure the dynamics of many signals within the network simultaneously. We found that, by fusing a fluorescent reporter to a pair of self-assembling peptides, it could be stably clustered within cells at random points, distant enough to be resolved by a microscope but close enough to spatially sample the relevant biology. Because such clusters, which we call signaling reporter islands (SiRIs), can be modularly designed, they permit a set of fluorescent reporters to be efficiently adapted for simultaneous measurement of multiple nodes of a signal transduction network within single cells. We created SiRIs for indicators of second messengers and kinases and used them, in hippocampal neurons in culture and intact brain slices, to discover relationships between the speed of calcium signaling, and the amplitude of PKA signaling, upon receiving a cAMP-driving stimulus.

SARS-CoV-2 Reverse Genetics Reveals a Variable Infection Gradient in the Respiratory Tract.

The mode of acquisition and causes for the variable clinical spectrum of coronavirus disease 2019 (COVID-19) remain unknown. We utilized a reverse genetics system to generate a GFP reporter virus to explore severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pathogenesis and a luciferase reporter virus to demonstrate sera collected from SARS and COVID-19 patients exhibited limited cross-CoV neutralization. High-sensitivity RNA in situ mapping revealed the highest angiotensin-converting enzyme 2 (ACE2) expression in the nose with decreasing expression throughout the lower respiratory tract, paralleled by a striking gradient of SARS-CoV-2 infection in proximal (high) versus distal (low) pulmonary epithelial cultures. COVID-19 autopsied lung studies identified focal disease and, congruent with culture data, SARS-CoV-2-infected ciliated and type 2 pneumocyte cells in airway and alveolar regions, respectively. These findings highlight the nasal susceptibility to SARS-CoV-2 with likely subsequent aspiration-mediated virus seeding to the lung in SARS-CoV-2 pathogenesis. These reagents provide a foundation for investigations into virus-host interactions in protective immunity, host susceptibility, and virus pathogenesis.

Promoter-Intrinsic and Local Chromatin Features Determine Gene Repression in LADs.

It is largely unclear whether genes that are naturally embedded in lamina-associated domains (LADs) are inactive due to their chromatin environment or whether LADs are merely secondary to the lack of transcription. We show that hundreds of human promoters become active when moved from their native LAD position to a neutral context in the same cells, indicating that LADs form a repressive environment. Another set of promoters inside LADs is able to "escape" repression, although their transcription elongation is attenuated. By inserting reporters into thousands of genomic locations, we demonstrate that escaper promoters are intrinsically less sensitive to LAD repression. This is not simply explained by promoter strength but by the interplay between promoter sequence and local chromatin features that vary strongly across LADs. Enhancers also differ in their sensitivity to LAD chromatin. This work provides a general framework for the systematic understanding of gene regulation by repressive chromatin.

A Deep Neural Network for Predicting and Engineering Alternative Polyadenylation.

Alternative polyadenylation (APA) is a major driver of transcriptome diversity in human cells. Here, we use deep learning to predict APA from DNA sequence alone. We trained our model (APARENT, APA REgression NeT) on isoform expression data from over 3 million APA reporters. APARENT's predictions are highly accurate when tasked with inferring APA in synthetic and human 3'UTRs. Visualizing features learned across all network layers reveals that APARENT recognizes sequence motifs known to recruit APA regulators, discovers previously unknown sequence determinants of 3' end processing, and integrates these features into a comprehensive, interpretable, cis-regulatory code. We apply APARENT to forward engineer functional polyadenylation signals with precisely defined cleavage position and isoform usage and validate predictions experimentally. Finally, we use APARENT to quantify the impact of genetic variants on APA. Our approach detects pathogenic variants in a wide range of disease contexts, expanding our understanding of the genetic origins of disease.

A Suite of Transgenic Driver and Reporter Mouse Lines with Enhanced Brain-Cell-Type Targeting and Functionality.

Modern genetic approaches are powerful in providing access to diverse cell types in the brain and facilitating the study of their function. Here, we report a large set of driver and reporter transgenic mouse lines, including 23 new driver lines targeting a variety of cortical and subcortical cell populations and 26 new reporter lines expressing an array of molecular tools. In particular, we describe the TIGRE2.0 transgenic platform and introduce Cre-dependent reporter lines that enable optical physiology, optogenetics, and sparse labeling of genetically defined cell populations. TIGRE2.0 reporters broke the barrier in transgene expression level of single-copy targeted-insertion transgenesis in a wide range of neuronal types, along with additional advantage of a simplified breeding strategy compared to our first-generation TIGRE lines. These novel transgenic lines greatly expand the repertoire of high-precision genetic tools available to effectively identify, monitor, and manipulate distinct cell types in the mouse brain.

Vitamin A-Retinoic Acid Signaling Regulates Hematopoietic Stem Cell Dormancy.

Dormant hematopoietic stem cells (dHSCs) are atop the hematopoietic hierarchy. The molecular identity of dHSCs and the mechanisms regulating their maintenance or exit from dormancy remain uncertain. Here, we use single-cell RNA sequencing (RNA-seq) analysis to show that the transition from dormancy toward cell-cycle entry is a continuous developmental path associated with upregulation of biosynthetic processes rather than a stepwise progression. In addition, low Myc levels and high expression of a retinoic acid program are characteristic for dHSCs. To follow the behavior of dHSCs in situ, a Gprc5c-controlled reporter mouse was established. Treatment with all-trans retinoic acid antagonizes stress-induced activation of dHSCs by restricting protein translation and levels of reactive oxygen species (ROS) and Myc. Mice maintained on a vitamin A-free diet lose HSCs and show a disrupted re-entry into dormancy after exposure to inflammatory stress stimuli. Our results highlight the impact of dietary vitamin A on the regulation of cell-cycle-mediated stem cell plasticity. VIDEO ABSTRACT.

Enhancer Function and Evolutionary Roles of Human Accelerated Regions.

Human accelerated regions (HARs) are the fastest-evolving sequences in the human genome. When HARs were discovered in 2006, their function was mysterious due to scant annotation of the noncoding genome. Diverse technologies, from transgenic animals to machine learning, have consistently shown that HARs function as gene regulatory enhancers with significant enrichment in neurodevelopment. It is now possible to quantitatively measure the enhancer activity of thousands of HARs in parallel and model how each nucleotide contributes to gene expression. These strategies have revealed that many human HAR sequences function differently than their chimpanzee orthologs, though individual nucleotide changes in the same HAR may have opposite effects, consistent with compensatory substitutions. To fully evaluate the role of HARs in human evolution, it will be necessary to experimentally and computationally dissect them across more cell types and developmental stages.

TRACE generates fluorescent human reporter cell lines to characterize epigenetic pathways.

Genetically encoded biosensors are powerful tools to monitor cellular behavior, but the difficulty in generating appropriate reporters for chromatin factors hampers our ability to dissect epigenetic pathways. Here, we present TRACE (transgene reporters across chromatin environments), a high-throughput, genome-wide technique to generate fluorescent human reporter cell lines responsive to manipulation of epigenetic factors. By profiling GFP expression from a large pool of individually barcoded lentiviral integrants in the presence and absence of a perturbation, we identify reporters responsive to pharmacological inhibition of the histone lysine demethylase LSD1 and genetic ablation of the PRC2 subunit SUZ12. Furthermore, by manipulating the HIV-1 host factor LEDGF through targeted deletion or fusion to chromatin reader domains, we alter lentiviral integration site preferences, thus broadening the types of chromatin examined by TRACE. The phenotypic reporters generated through TRACE will allow the genetic interrogation of a broad range of epigenetic pathways, furthering our mechanistic understanding of chromatin biology.

Impact of chromatin context on Cas9-induced DNA double-strand break repair pathway balance.

DNA double-strand break (DSB) repair is mediated by multiple pathways. It is thought that the local chromatin context affects the pathway choice, but the underlying principles are poorly understood. Using a multiplexed reporter assay in combination with Cas9 cutting, we systematically measure the relative activities of three DSB repair pathways as a function of chromatin context in >1,000 genomic locations. This reveals that non-homologous end-joining (NHEJ) is broadly biased toward euchromatin, while the contribution of microhomology-mediated end-joining (MMEJ) is higher in specific heterochromatin contexts. In H3K27me3-marked heterochromatin, inhibition of the H3K27 methyltransferase EZH2 reverts the balance toward NHEJ. Single-stranded template repair (SSTR), often used for precise CRISPR editing, competes with MMEJ and is moderately linked to chromatin context. These results provide insight into the impact of chromatin on DSB repair pathway balance and guidance for the design of Cas9-mediated genome editing experiments.

Mouse models to investigate in situ cell fate decisions induced by p53.

Investigating how transcription factors control complex cellular processes requires tools that enable responses to be visualised at the single-cell level and their cell fate to be followed over time. For example, the tumour suppressor p53 (also called TP53 in humans and TRP53 in mice) can initiate diverse cellular responses by transcriptional activation of its target genes: Puma to induce apoptotic cell death and p21 to induce cell cycle arrest/cell senescence. However, it is not known how these processes are regulated and initiated in different cell types. Also, the context-dependent interaction partners and binding loci of p53 remain largely elusive. To be able to examine these questions, we here developed knock-in mice expressing triple-FLAG-tagged p53 to facilitate p53 pull-down and two p53 response reporter mice, knocking tdTomato and GFP into the Puma/Bbc3 and p21 gene loci, respectively. By crossing these reporter mice into a p53-deficient background, we show that the new reporters reliably inform on p53-dependent and p53-independent initiation of both apoptotic or cell cycle arrest/senescence programs, respectively, in vitro and in vivo.

Fate-Mapping of GM-CSF Expression Identifies a Discrete Subset of Inflammation-Driving T Helper Cells Regulated by Cytokines IL-23 and IL-1ß.

Pathogenic lymphocytes initiate the development of chronic inflammatory diseases. The cytokine granulocyte-macrophage colony-stimulating factor (GM-CSF) (encoded by Csf2) is a key communicator between pathogenic lymphocytes and tissue-invading inflammatory phagocytes. However, the molecular properties of GM-CSF-producing cells and the mode of Csf2 regulation in vivo remain unclear. To systematically study and manipulate GM-CSF+ cells and their progeny in vivo, we generated a fate-map and reporter of GM-CSF expression mouse strain (FROG). We mapped the phenotypic and functional profile of auto-aggressive T helper (Th) cells during neuroinflammation and identified the signature and pathogenic memory of a discrete encephalitogenic Th subset. These cells required interleukin-23 receptor (IL-23R) and IL-1R but not IL-6R signaling for their maintenance and pathogenicity. Specific ablation of this subset interrupted the inflammatory cascade, despite the unperturbed tissue accumulation of other Th subsets (e.g., Th1 and Th17), highlighting that GM-CSF expression not only marks pathogenic Th cells, but that this subset mediates immunopathology and tissue destruction.

XACT-Seq Comprehensively Defines the Promoter-Position and Promoter-Sequence Determinants for Initial-Transcription Pausing.

Pausing by RNA polymerase (RNAP) during transcription elongation, in which a translocating RNAP uses a "stepping" mechanism, has been studied extensively, but pausing by RNAP during initial transcription, in which a promoter-anchored RNAP uses a "scrunching" mechanism, has not. We report a method that directly defines the RNAP-active-center position relative to DNA with single-nucleotide resolution (XACT-seq; "crosslink-between-active-center-and-template sequencing"). We apply this method to detect and quantify pausing in initial transcription at 411 (∼4,000,000) promoter sequences in vivo in Escherichia coli. The results show initial-transcription pausing can occur in each nucleotide addition during initial transcription, particularly the first 4 to 5 nucleotide additions. The results further show initial-transcription pausing occurs at sequences that resemble the consensus sequence element for transcription-elongation pausing. Our findings define the positional and sequence determinants for initial-transcription pausing and establish initial-transcription pausing is hard coded by sequence elements similar to those for transcription-elongation pausing.

Systematic identification of genotype-dependent enhancer variants in eosinophilic esophagitis.

Eosinophilic esophagitis (EoE) is a rare atopic disorder associated with esophageal dysfunction, including difficulty swallowing, food impaction, and inflammation, that develops in a small subset of people with food allergies. Genome-wide association studies (GWASs) have identified 9 independent EoE risk loci reaching genome-wide significance (p < 5 × 10-8) and 27 additional loci of suggestive significance (5 × 10-8 < p < 1 × 10-5). In the current study, we perform linkage disequilibrium (LD) expansion of these loci to nominate a set of 531 variants that are potentially causal. To systematically interrogate the gene regulatory activity of these variants, we designed a massively parallel reporter assay (MPRA) containing the alleles of each variant within their genomic sequence context cloned into a GFP reporter library. Analysis of reporter gene expression in TE-7, HaCaT, and Jurkat cells revealed cell-type-specific gene regulation. We identify 32 allelic enhancer variants, representing 6 genome-wide significant EoE loci and 7 suggestive EoE loci, that regulate reporter gene expression in a genotype-dependent manner in at least one cellular context. By annotating these variants with expression quantitative trait loci (eQTL) and chromatin looping data in related tissues and cell types, we identify putative target genes affected by genetic variation in individuals with EoE. Transcription factor enrichment analyses reveal possible roles for cell-type-specific regulators, including GATA3. Our approach reduces the large set of EoE-associated variants to a set of 32 with allelic regulatory activity, providing functional insights into the effects of genetic variation in this disease.

High-throughput characterization of functional variants highlights heterogeneity and polygenicity underlying lung cancer susceptibility.

Genome-wide association studies (GWASs) have identified numerous lung cancer risk-associated loci. However, decoding molecular mechanisms of these associations is challenging since most of these genetic variants are non-protein-coding with unknown function. Here, we implemented massively parallel reporter assays (MPRAs) to simultaneously measure the allelic transcriptional activity of risk-associated variants. We tested 2,245 variants at 42 loci from 3 recent GWASs in East Asian and European populations in the context of two major lung cancer histological types and exposure to benzo(a)pyrene. This MPRA approach identified one or more variants (median 11 variants) with significant effects on transcriptional activity at 88% of GWAS loci. Multimodal integration of lung-specific epigenomic data demonstrated that 63% of the loci harbored multiple potentially functional variants in linkage disequilibrium. While 22% of the significant variants showed allelic effects in both A549 (adenocarcinoma) and H520 (squamous cell carcinoma) cell lines, a subset of the functional variants displayed a significant cell-type interaction. Transcription factor analyses nominated potential regulators of the functional variants, including those with cell-type-specific expression and those predicted to bind multiple potentially functional variants across the GWAS loci. Linking functional variants to target genes based on four complementary approaches identified candidate susceptibility genes, including those affecting lung cancer cell growth. CRISPR interference of the top functional variant at 20q13.33 validated variant-to-gene connections, including RTEL1, SOX18, and ARFRP1. Our data provide a comprehensive functional analysis of lung cancer GWAS loci and help elucidate the molecular basis of heterogeneity and polygenicity underlying lung cancer susceptibility.