Search

MicroRNAs (miRNAs) are short, noncoding RNAs that associate with Argonaute (AGO) to influence mRNA stability and translation, thereby regulating cellular determination and phenotype. While several individual miRNAs have been shown to control adipocyte function, including energy storage in white fat and energy dissipation in brown fat, a comprehensive analysis of miRNA activity in these tissues has not been performed. We used high-throughput sequencing of RNA isolated by cross-linking immunoprecipitation (HITS-CLIP) to comprehensively characterize the network of high-confidence, in vivo mRNA:miRNA interactions across white and brown fat, revealing >20,000 unique AGO binding sites. When coupled with miRNA and mRNA sequencing, we found an inverse correlation between depot-enriched miRNAs and their targets. To illustrate the functionality of our HITS-CLIP data set in identifying specific miRNA:mRNA interactions, we show that miR-29 is a novel regulator of leptin, an adipocyte-derived hormone that coordinates food intake and energy homeostasis. Two independent miR-29 binding sites in the leptin 3′ UTR were validated using luciferase assays, and miR-29 gain and loss of function modulated leptin mRNA and protein secretion in primary adipocytes. This work represents the only experimentally generated miRNA targetome in adipose tissue and identifies multiple regulatory pathways that may specify the unique identities of white and brown fat.

Although branchpoint recognition is an essential component of intron excision during the RNA splicing process, the branchpoint itself is frequently assumed to be a basal, rather than regulatory, sequence feature. However, this assumption has not been systematically tested due to the technical difficulty of identifying branchpoints and quantifying their usage. Here, we analyzed ∼1.31 trillion reads from 17,164 RNA sequencing data sets to demonstrate that almost all human introns contain multiple branchpoints. This complexity holds even for constitutive introns, 95% of which contain multiple branchpoints, with an estimated five to six branchpoints per intron. Introns upstream of the highly regulated ultraconserved poison exons of SR genes contain twice as many branchpoints as the genomic average. Approximately three-quarters of constitutive introns exhibit tissue-specific branchpoint usage. In an extreme example, we observed a complete switch in branchpoint usage in the well-studied first intron of HBB (β-globin) in normal bone marrow versus metastatic prostate cancer samples. Our results indicate that the recognition of most introns is unexpectedly complex and tissue-specific and suggest that alternative splicing catalysis typifies the majority of introns even in the absence of differences in the mature mRNA.

Cells communicate with each other to coordinate their gene activities at the population level through signaling pathways. It has been shown that many gene activities are oscillatory and that the frequency and phase of oscillatory gene expression encode various types of information. However, whether or how such oscillatory information is transmitted from cell to cell remains unknown. Here, we developed an integrated approach that combines optogenetic perturbations and single-cell bioluminescence imaging to visualize and reconstitute synchronized oscillatory gene expression in signal-sending and signal-receiving processes. We found that intracellular and intercellular periodic inputs of Notch signaling entrain intrinsic oscillations by frequency tuning and phase shifting at the single-cell level. In this way, the oscillation dynamics are transmitted through Notch signaling, thereby synchronizing the population of oscillators. Thus, this approach enabled us to control and monitor dynamic cell-to-cell transfer of oscillatory information to coordinate gene expression patterns at the population level.

To understand mammalian active DNA demethylation, various methods have been developed to map the genomic distribution of the demethylation intermediates 5-formylcysotine (5fC) and 5-carboxylcytosine (5caC). However, the majority of these methods requires a large number of cells to begin with. In this study, we describe low-input methylase-assisted bisulfite sequencing (liMAB-seq ) and single-cell MAB-seq (scMAB-seq), capable of profiling 5fC and 5caC at genome scale using ∼100 cells and single cells, respectively. liMAB-seq analysis of preimplantation embryos reveals the oxidation of 5mC to 5fC/5caC and the positive correlation between chromatin accessibility and processivity of ten-eleven translocation (TET) enzymes. scMAB-seq captures the cell-to-cell heterogeneity of 5fC and 5caC and reveals the strand-biased distribution of 5fC and 5caC. scMAB-seq also allows the simultaneous high-resolution mapping of sister chromatid exchange (SCE), facilitating the study of this type of genomic rearrangement. Therefore, our study not only establishes new methods for the genomic mapping of active DNA demethylation using limited numbers of cells or single cells but also demonstrates the utilities of the methods in different biological contexts.

Bromodomain proteins (BRD) are key chromatin regulators of genome function and stability as well as therapeutic targets in cancer. Here, we systematically delineate the contribution of human BRD proteins for genome stability and DNA double-strand break (DSB) repair using several cell-based assays and proteomic interaction network analysis. Applying these approaches, we identify 24 of the 42 BRD proteins as promoters of DNA repair and/or genome integrity. We identified a BRD-reader function of PCAF that bound TIP60-mediated histone acetylations at DSBs to recruit a DUB complex to deubiquitylate histone H2BK120, to allowing direct acetylation by PCAF, and repair of DSBs by homologous recombination. We also discovered the bromo-and-extra-terminal (BET) BRD proteins, BRD2 and BRD4, as negative regulators of transcription-associated RNA-DNA hybrids (R-loops) as inhibition of BRD2 or BRD4 increased R-loop formation, which generated DSBs. These breaks were reliant on topoisomerase II, and BRD2 directly bound and activated topoisomerase I, a known restrainer of R-loops. Thus, comprehensive interactome and functional profiling of BRD proteins revealed new homologous recombination and genome stability pathways, providing a framework to understand genome maintenance by BRD proteins and the effects of their pharmacological inhibition.

Addressing the complexity of organogenesis at a system-wide level requires a complete understanding of adult cell types, their origin, and precursor relationships. The Drosophila ovary has been a model to study how coordinated stem cell units, germline, and somatic follicle stem cells maintain and renew an organ. However, lack of cell type-specific tools have limited our ability to study the origin of individual cell types and stem cell units. Here, we used a single-cell RNA sequencing approach to uncover all known cell types of the developing ovary, reveal transcriptional signatures, and identify cell type-specific markers for lineage tracing. Our study identifies a novel cell type corresponding to the elusive follicle stem cell precursors and predicts subtypes of known cell types. Altogether, we reveal a previously unanticipated complexity of the developing ovary and provide a comprehensive resource for the systematic analysis of ovary morphogenesis.

Rapid perturbation of protein function permits the ability to define primary molecular responses while avoiding down-stream cumulative effects of protein dysregulation. The auxin-inducible degron (AID) system was developed as a tool to achieve rapid and inducible protein degradation in non-plant systems. However, tagging proteins at their endogenous loci results in chronic, auxin-independent degradation by the proteasome. To correct this deficiency, we expressed the Auxin Response Transcription Factor (ARF) in an improved inducible degron system. ARF is absent from previously engineered AID systems, but ARF is a critical component of native auxin signaling. In plants, ARF directly interacts with AID in the absence of auxin and we found that expression of the ARF Phox and Bem1 (PB1) domain suppresses constitutive degradation of AID-tagged proteins. Moreover, the rate of auxin-induced AID degradation is substantially faster in the ARF-AID system. To test the ARF-AID system in a quantitative and sensitive manner, we measured genome-wide changes in nascent transcription after rapidly depleting the ZNF143 transcription factor. Transciptional profiling indicates that ZNF143 activates transcription in cis and ZNF143 regulates promoter-proximal paused RNA Polymerase density. Rapidly inducible degradation systems that preserve the target protein’s native expression levels and patterns will revolutionize the study of biological systems by enabling specific and temporally defined protein dysregulation.

Eukaryotic mRNAs are subject to multiple types of tailing that critically influence mRNA stability and translatability. To investigate RNA tails at the genomic scale, we previously developed TAIL-seq, but its low sensitivity precluded its application to biological materials of minute quantity. In this study, we report a new version of TAIL-seq (mRNA TAIL-seq [mTAIL-seq]) with enhanced sequencing depth for mRNAs (by ∼1000-fold compared with the previous version). The improved method allows us to investigate the regulation of poly(A) tails in Drosophila oocytes and embryos. We found that maternal mRNAs are polyadenylated mainly during late oogenesis, prior to fertilization, and that further modulation occurs upon egg activation. Wispy, a noncanonical poly(A) polymerase, adenylates the vast majority of maternal mRNAs, with a few intriguing exceptions such as ribosomal protein transcripts. By comparing mTAIL-seq data with ribosome profiling data, we found a strong coupling between poly(A) tail length and translational efficiency during egg activation. Our data suggest that regulation of poly(A) tails in oocytes shapes the translatomic landscape of embryos, thereby directing the onset of animal development. By virtue of the high sensitivity, low cost, technical robustness, and broad accessibility, mTAIL-seq will be a potent tool to improve our understanding of mRNA tailing in diverse biological systems.

A fundamental challenge in understanding cardiac biology and disease is that the remarkable heterogeneity in cell type composition and functional states have not been well characterized at single-cell resolution in maturing and diseased mammalian hearts. Massively parallel single-nucleus RNA sequencing (snRNA-seq) has emerged as a powerful tool to address these questions by interrogating the transcriptome of tens of thousands of nuclei isolated from fresh or frozen tissues. snRNA-seq overcomes the technical challenge of isolating intact single cells from complex tissues, including the maturing mammalian hearts; reduces biased recovery of easily dissociated cell types; and minimizes aberrant gene expression during the whole-cell dissociation. Here we applied sNucDrop-seq, a droplet microfluidics-based massively parallel snRNA-seq method, to investigate the transcriptional landscape of postnatal maturing mouse hearts in both healthy and disease states. By profiling the transcriptome of nearly 20,000 nuclei, we identified major and rare cardiac cell types and revealed significant heterogeneity of cardiomyocytes, fibroblasts, and endothelial cells in postnatal developing hearts. When applied to a mouse model of pediatric mitochondrial cardiomyopathy, we uncovered profound cell type-specific modifications of the cardiac transcriptional landscape at single-nucleus resolution, including changes of subtype composition, maturation states, and functional remodeling of each cell type. Furthermore, we employed sNucDrop-seq to decipher the cardiac cell type-specific gene regulatory network (GRN) of GDF15, a heart-derived hormone and clinically important diagnostic biomarker of heart disease. Together, our results present a rich resource for studying cardiac biology and provide new insights into heart disease using an approach broadly applicable to many fields of biomedicine.

Postnatal spermatogonial stem cells (SSCs) progress through proliferative and developmental stages to populate the testicular niche prior to productive spermatogenesis. To better understand, we conducted extensive genomic profiling at multiple postnatal stages on subpopulations enriched for particular markers (THY1, KIT, OCT4, ID4, or GFRa1). Overall, our profiles suggest three broad populations of spermatogonia in juveniles: (1) epithelial-like spermatogonia (THY1+; high OCT4, ID4, and GFRa1), (2) more abundant mesenchymal-like spermatogonia (THY1+; moderate OCT4 and ID4; high mesenchymal markers), and (3) (in older juveniles) abundant spermatogonia committing to gametogenesis (high KIT+). Epithelial-like spermatogonia displayed the expected imprinting patterns, but, surprisingly, mesenchymal-like spermatogonia lacked imprinting specifically at paternally imprinted loci but fully restored imprinting prior to puberty. Furthermore, mesenchymal-like spermatogonia also displayed developmentally linked DNA demethylation at meiotic genes and also at certain monoallelic neural genes (e.g., protocadherins and olfactory receptors). We also reveal novel candidate receptor–ligand networks involving SSCs and the developing niche. Taken together, neonates/juveniles contain heterogeneous epithelial-like or mesenchymal-like spermatogonial populations, with the latter displaying extensive DNA methylation/chromatin dynamics. We speculate that this plasticity helps SSCs proliferate and migrate within the developing seminiferous tubule, with proper niche interaction and membrane attachment reverting mesenchymal-like spermatogonial subtype cells back to an epithelial-like state with normal imprinting profiles.

Maturation of the large ribosomal subunit (LSU) in eukaryotes is a complex and highly coordinated process that requires the concerted action of a large, dynamic, ribonucleoprotein complex, the LSU processome. While we know that >80 ribosome biogenesis factors are required throughout the course of LSU assembly, little is known about how these factors interact with each other within the LSU processome. To interrogate its organization and architecture, we took a systems biology approach and performed a semi-high-throughput, array-based, directed yeast two-hybrid assay. Assaying 4800 protein–protein interactions, we identified 232 high-confidence, binary-interacting protein pairs, representing a fourfold increase from current knowledge. The resulting LSU processome interactome map has enhanced our understanding of the organization and function of the biogenesis factors within the LSU processome, revealing both novel and previously identified subcomplexes and hub proteins, including Nop4.

Sequence variation in tRNA genes influences the structure, modification, and stability of tRNA; affects translation fidelity; impacts the activity of numerous isodecoders in metazoans; and leads to human diseases. To comprehensively define the effects of sequence variation on tRNA function, we developed a high-throughput in vivo screen to quantify the activity of a model tRNA, the nonsense suppressor SUP4oc of Saccharomyces cerevisiae. Using a highly sensitive fluorescent reporter gene with an ochre mutation, fluorescence-activated cell sorting of a library of SUP4oc mutant yeast strains, and deep sequencing, we scored 25,491 variants. Unexpectedly, SUP4oc tolerates numerous sequence variations, accommodates slippage in tertiary and secondary interactions, and exhibits genetic interactions that suggest an alternative functional tRNA conformation. Furthermore, we used this methodology to define tRNA variants subject to rapid tRNA decay (RTD). Even though RTD normally degrades tRNAs with exposed 5′ ends, mutations that sensitize SUP4oc to RTD were found to be located throughout the sequence, including the anti-codon stem. Thus, the integrity of the entire tRNA molecule is under surveillance by cellular quality control machinery. This approach to assess activity at high throughput is widely applicable to many problems in tRNA biology.

The mammalian intestinal epithelium has a unique organization in which crypts harboring stem cells produce progenitors and finally clonal populations of differentiated cells. Remarkably, the epithelium is replaced every 3–5 d throughout adult life. Disrupted maintenance of the intricate balance of proliferation and differentiation leads to loss of epithelial integrity or barrier function or to cancer. There is a tight correlation between the epigenetic status of genes and expression changes during differentiation; however, the mechanism of how changes in DNA methylation direct gene expression and the progression from stem cells to their differentiated descendants is unclear. Using conditional gene ablation of the maintenance methyltransferase Dnmt1, we demonstrate that reducing DNA methylation causes intestinal crypt expansion in vivo. Determination of the base-resolution DNA methylome in intestinal stem cells and their differentiated descendants shows that DNA methylation is dynamic at enhancers, which are often associated with genes important for both stem cell maintenance and differentiation. We establish that the loss of DNA methylation at intestinal stem cell gene enhancers causes inappropriate gene expression and delayed differentiation.

More than half of human genes use alternative cleavage and polyadenylation (ApA) to generate mRNA transcripts that differ in the lengths of their 3′ untranslated regions (UTRs), thus altering the post-transcriptional fate of the message and likely the protein output. The extent of 3′ UTR variation across tissues and the functional role of ApA remain poorly understood. We developed a sequencing method called 3′-seq to quantitatively map the 3′ ends of the transcriptome of diverse human tissues and isogenic transformation systems. We found that cell type-specific gene expression is accomplished by two complementary programs. Tissue-restricted genes tend to have single 3′ UTRs, whereas a majority of ubiquitously transcribed genes generate multiple 3′ UTRs. During transformation and differentiation, single-UTR genes change their mRNA abundance levels, while multi-UTR genes mostly change 3′ UTR isoform ratios to achieve tissue specificity. However, both regulation programs target genes that function in the same pathways and processes that characterize the new cell type. Instead of finding global shifts in 3′ UTR length during transformation and differentiation, we identify tissue-specific groups of multi-UTR genes that change their 3′ UTR ratios; these changes in 3′ UTR length are largely independent from changes in mRNA abundance. Finally, tissue-specific usage of ApA sites appears to be a mechanism for changing the landscape targetable by ubiquitously expressed microRNAs.

Understanding how distinct cell types arise from multipotent progenitor cells is a major quest in stem cell biology. The liver and pancreas share many aspects of their early development and possibly originate from a common progenitor. However, how liver and pancreas cells diverge from a common endoderm progenitor population and adopt specific fates remains elusive. Using RNA sequencing (RNA-seq), we defined the molecular identity of liver and pancreas progenitors that were isolated from the mouse embryo at two time points, spanning the period when the lineage decision is made. The integration of temporal and spatial gene expression profiles unveiled mutually exclusive signaling signatures in hepatic and pancreatic progenitors. Importantly, we identified the noncanonical Wnt pathway as a potential developmental regulator of this fate decision and capable of inducing the pancreas program in endoderm and liver cells. Our study offers an unprecedented view of gene expression programs in liver and pancreas progenitors and forms the basis for formulating lineage-reprogramming strategies to convert adult hepatic cells into pancreatic cells.

Monitoring protein synthesis is essential to our understanding of gene expression regulation, as protein abundance is thought to be predominantly controlled at the level of translation. Mass-spectrometric and RNA sequencing methods have been recently developed for investigating mRNA translation at a global level, but these still involve technical limitations and are not widely applicable. In this study, we describe a novel system-wide proteomic approach for direct monitoring of translation, termed puromycin-associated nascent chain proteomics (PUNCH-P), which is based on incorporation of biotinylated puromycin into newly synthesized proteins under cell-free conditions followed by streptavidin affinity purification and liquid chromatography-tandem mass spectrometry analysis. Using PUNCH-P, we measured cell cycle-specific fluctuations in synthesis for >5000 proteins in mammalian cells, identified proteins not previously implicated in cell cycle processes, and generated the first translational profile of a whole mouse brain. This simple and economical technique is broadly applicable to any cell type and tissue, enabling the identification and quantification of rapid proteome responses under various biological conditions.

Cytosine methylation is a key epigenetic mark in many organisms, important for both transcriptional control and genome integrity. While relatively stable during somatic growth, DNA methylation is reprogrammed genome-wide during mammalian reproduction. Reprogramming is essential for zygotic totipotency and to prevent transgenerational inheritance of epimutations. However, the extent of DNA methylation reprogramming in plants remains unclear. Here, we developed sensors reporting with single-cell resolution CG and non-CG methylation in Arabidopsis. Live imaging during reproduction revealed distinct and sex-specific dynamics for both contexts. We found that CHH methylation in the egg cell depends on DOMAINS REARRANGED METHYLASE 2 (DRM2) and RNA polymerase V (Pol V), two main actors of RNA-directed DNA methylation, but does not depend on Pol IV. Our sensors provide insight into global DNA methylation dynamics at the single-cell level with high temporal resolution and offer a powerful tool to track CG and non-CG methylation both during development and in response to environmental cues in all organisms with methylated DNA, as we illustrate in mouse embryonic stem cells.

Transcriptional profiling is a powerful approach for understanding development and disease. Current cell type-specific RNA purification methods have limitations, including cell dissociation trauma or inability to identify all RNA species. Here, we describe “mouse thiouracil (TU) tagging,” a genetic and chemical intersectional method for covalent labeling and purification of cell type-specific RNA in vivo. Cre-induced expression of uracil phosphoribosyltransferase (UPRT) provides spatial specificity; injection of 4-thiouracil (4TU) provides temporal specificity. Only UPRT+ cells exposed to 4TU produce thio-RNA, which is then purified for RNA sequencing (RNA-seq). This method can purify transcripts from spatially complex and rare (Tie2:Cre+ brain endothelia/microglia (76% validated by expression pattern), or temporally dynamic transcripts, such as those acutely induced by lipopolysaccharide (LPS) injection. Moreover, generating chimeric mice via UPRT+ bone marrow transplants identifies immune versus niche spleen RNA. TU tagging provides a novel method for identifying actively transcribed genes in specific cells at specific times within intact mice.

Large-scale sequencing studies are rapidly identifying putative oncogenic mutations in human tumors. However, discrimination between passenger and driver events in tumorigenesis remains challenging and requires in vivo validation studies in reliable animal models of human cancer. In this study, we describe a novel strategy for in vivo validation of candidate tumor suppressors implicated in invasive lobular breast carcinoma (ILC), which is hallmarked by loss of the cell–cell adhesion molecule E-cadherin. We describe an approach to model ILC by intraductal injection of lentiviral vectors encoding Cre recombinase, the CRISPR/Cas9 system, or both in female mice carrying conditional alleles of the Cdh1 gene, encoding for E-cadherin. Using this approach, we were able to target ILC-initiating cells and induce specific gene disruption of Pten by CRISPR/Cas9-mediated somatic gene editing. Whereas intraductal injection of Cas9-encoding lentiviruses induced Cas9-specific immune responses and development of tumors that did not resemble ILC, lentiviral delivery of a Pten targeting single-guide RNA (sgRNA) in mice with mammary gland-specific loss of E-cadherin and expression of Cas9 efficiently induced ILC development. This versatile platform can be used for rapid in vivo testing of putative tumor suppressor genes implicated in ILC, providing new opportunities for modeling invasive lobular breast carcinoma in mice.

The discovery of transcription factors (TFs) controlling pathways in health and disease is of paramount interest. We designed a widely applicable method, dubbed barcorded synthetic tandem repeat promoter screening (BC-STAR-PROM), to identify signal-activated TFs without any a priori knowledge about their properties. The BC-STAR-PROM library consists of ∼3000 luciferase expression vectors, each harboring a promoter (composed of six tandem repeats of synthetic random DNA) and an associated barcode of 20 base pairs (bp) within the 3′ untranslated mRNA region. Together, the promoter sequences encompass >400,000 bp of random DNA, a sequence complexity sufficient to capture most TFs. Cells transfected with the library are exposed to a signal, and the mRNAs that it encodes are counted by next-generation sequencing of the barcodes. This allows the simultaneous activity tracking of each of the ∼3000 synthetic promoters in a single experiment. Here we establish proof of concept for BC-STAR-PROM by applying it to the identification of TFs induced by drugs affecting actin and tubulin cytoskeleton dynamics. BC-STAR-PROM revealed that serum response factor (SRF) is the only immediate early TF induced by both actin polymerization and microtubule depolymerization. Such changes in cytoskeleton dynamics are known to occur during the cell division cycle, and real-time bioluminescence microscopy indeed revealed cell-autonomous SRF–myocardin-related TF (MRTF) activity bouts in proliferating cells.

In a previous analysis of 2300 mRNAs via whole-mount fluorescent in situ hybridization in cellularizing Drosophila embryos, we found that 70% of the transcripts exhibited some form of subcellular localization. To see whether this prevalence is unique to early Drosophila embryos, we examined ∼8000 transcripts over the full course of embryogenesis and ∼800 transcripts in late third instar larval tissues. The numbers and varieties of new subcellular localization patterns are both striking and revealing. In the much larger cells of the third instar larva, virtually all transcripts observed showed subcellular localization in at least one tissue. We also examined the prevalence and variety of localization mechanisms for >100 long noncoding RNAs. All of these were also found to be expressed and subcellularly localized. Thus, subcellular RNA localization appears to be the norm rather than the exception for both coding and noncoding RNAs. These results, which have been annotated and made available on a recompiled database, provide a rich and unique resource for functional gene analyses, some examples of which are provided.

Diverse proteins are known to be capable of forming amyloid aggregates, self-seeding fibrillar assemblies that may be biologically functional or pathological. Well-known examples include neurodegenerative disease-associated proteins that misfold as amyloid, fungal prion proteins that can transition to a self-propagating amyloid form and certain bacterial proteins that fold as amyloid at the cell surface and promote biofilm formation. To further explore the diversity of amyloidogenic proteins, generally applicable methods for identifying them are critical. Here we describe a cell-based method for generating amyloid aggregates that relies on the natural ability of Escherichia coli cells to elaborate amyloid fibrils at the cell surface. We use several different yeast prion proteins and the human huntingtin protein to show that protein secretion via this specialized export pathway promotes acquisition of the amyloid fold specifically for proteins that have an inherent amyloid-forming propensity. Furthermore, our findings establish the potential of this E. coli-based system to facilitate the implementation of high-throughput screens for identifying amyloidogenic proteins and modulators of amyloid aggregation.

We report a multifunctional gene-trapping approach, which generates full-length Citrine fusions with endogenous proteins and conditional mutants from a single integration event of the FlipTrap vector. We identified 170 FlipTrap zebrafish lines with diverse tissue-specific expression patterns and distinct subcellular localizations of fusion proteins generated by the integration of an internal citrine exon. Cre-mediated conditional mutagenesis is enabled by heterotypic lox sites that delete Citrine and “flip” in its place mCherry with a polyadenylation signal, resulting in a truncated fusion protein. Inducing recombination with Cerulean-Cre results in fusion proteins that often mislocalize, exhibit mutant phenotypes, and dramatically knock down wild-type transcript levels. FRT sites in the vector enable targeted genetic manipulation of the trapped loci in the presence of Flp recombinase. Thus, the FlipTrap captures the functional proteome, enabling the visualization of full-length fluorescent fusion proteins and interrogation of function by conditional mutagenesis and targeted genetic manipulation.

Transfer ribonucleic acids (tRNAs) are essential for protein synthesis. However, key gene products involved in tRNA biogenesis and subcellular movement remain to be discovered. We conducted the first comprehensive unbiased analysis of the role of nearly an entire proteome in tRNA biology and describe 162 novel and 12 previously known Saccharomyces cerevisiae gene products that function in tRNA processing, turnover, and subcellular movement. tRNA nuclear export is of particular interest because it is essential, but the known tRNA exporters (Los1 [exportin-t] and Msn5 [exportin-5]) are unessential. We report that mutations of CRM1 (Exportin-1), MEX67/MTR2 (TAP/p15), and five nucleoporins cause accumulation of unspliced tRNA, a hallmark of defective tRNA nuclear export. CRM1 mutation genetically interacts with los1Δ and causes altered tRNA nuclear–cytoplasmic distribution. The data implicate roles for the protein and mRNA nuclear export machineries in tRNA nuclear export. Mutations of genes encoding actin cytoskeleton components and mitochondrial outer membrane proteins also cause accumulation of unspliced tRNA, likely due to defective splicing on mitochondria. Additional gene products, such as chromatin modification enzymes, have unanticipated effects on pre-tRNA end processing. Thus, this genome-wide screen uncovered putative novel pathways for tRNA nuclear export and extensive links between tRNA biology and other aspects of cell physiology.

In the mammalian pancreas, endocrine cells undergo lineage allocation upon emergence from a bipotent duct/endocrine progenitor pool, which resides in the “trunk epithelium.” Major questions remain regarding how niche environments are organized within this epithelium to coordinate endocrine differentiation with programs of epithelial growth, maturation, and morphogenesis. We used EdU pulse-chase and tissue-reconstruction approaches to analyze how endocrine progenitors and their differentiating progeny are assembled within the trunk as it undergoes remodeling from an irregular plexus of tubules to form the eventual mature, branched ductal arbor. The bulk of endocrine progenitors is maintained in an epithelial “plexus state,” which is a transient intermediate during epithelial maturation within which endocrine cell differentiation is continually robust and surprisingly long-lived. Within the plexus, local feedback effects derived from the differentiating and delaminating endocrine cells nonautonomously regulate the flux of endocrine cell birth as well as proliferative growth of the bipotent cell population using Notch-dependent and Notch-independent influences, respectively. These feedback effects in turn maintain the plexus state to ensure prolonged allocation of endocrine cells late into gestation. These findings begin to define a niche-like environment guiding the genesis of the endocrine pancreas and advance current models for how differentiation is coordinated with the growth and morphogenesis of the developing pancreatic epithelium.

Pancreatic ductal adenocarcinoma (PDAC) is a genomically diverse, prevalent, and almost invariably fatal malignancy. Although conventional genetically engineered mouse models of human PDAC have been instrumental in understanding pancreatic cancer development, these models are much too labor-intensive, expensive, and slow to perform the extensive molecular analyses needed to adequately understand this disease. Here we demonstrate that retrograde pancreatic ductal injection of either adenoviral-Cre or lentiviral-Cre vectors allows titratable initiation of pancreatic neoplasias that progress into invasive and metastatic PDAC. To enable in vivo CRISPR/Cas9-mediated gene inactivation in the pancreas, we generated a Cre-regulated Cas9 allele and lentiviral vectors that express Cre and a single-guide RNA. CRISPR-mediated targeting of Lkb1 in combination with oncogenic Kras expression led to selection for inactivating genomic alterations, absence of Lkb1 protein, and rapid tumor growth that phenocopied Cre-mediated genetic deletion of Lkb1. This method will transform our ability to rapidly interrogate gene function during the development of this recalcitrant cancer.

The retinoblastoma tumor suppressor (pRb) protein associates with chromatin and regulates gene expression. Numerous studies have identified Rb-dependent RNA signatures, but the proteomic effects of Rb loss are largely unexplored. We acutely ablated Rb in adult mice and conducted a quantitative analysis of RNA and proteomic changes in the colon and lungs, where Rb[superscript KO] was sufficient or insufficient to induce ectopic proliferation, respectively. As expected, Rb[superscript KO] caused similar increases in classic pRb/E2F-regulated transcripts in both tissues, but, unexpectedly, their protein products increased only in the colon, consistent with its increased proliferative index. Thus, these protein changes induced by Rb loss are coupled with proliferation but uncoupled from transcription. The proteomic changes in common between Rb[superscript KO] tissues showed a striking decrease in proteins with mitochondrial functions. Accordingly, RB1 inactivation in human cells decreased both mitochondrial mass and oxidative phosphorylation (OXPHOS) function. Rb[superscript KO] cells showed decreased mitochondrial respiratory capacity and the accumulation of hypopolarized mitochondria. Additionally, RB/Rb loss altered mitochondrial pyruvate oxidation from [superscript 13]C-glucose through the TCA cycle in mouse tissues and cultured cells. Consequently, Rb[superscript KO] cells have an enhanced sensitivity to mitochondrial stress conditions. In summary, proteomic analyses provide a new perspective on Rb/RB1 mutation, highlighting the importance of pRb for mitochondrial function and suggesting vulnerabilities for treatment., National Institutes of Health (U.S.) (2-P01-CA42063)

Repetitive nucleotide or amino acid sequences are often engineered into probes and biosensors to achieve functional readouts and robust signal amplification. However, these repeated sequences are notoriously prone to aberrant deletion and degradation, impacting the ability to correctly detect and interpret biological functions. Here, we introduce a facile and generalizable approach to solve this often unappreciated problem by modifying the nucleotide sequences of the target mRNA to make them nonrepetitive but still functional (“synonymous”). We first demonstrated the procedure by designing a cassette of synonymous MS2 RNA motifs and tandem coat proteins for RNA imaging and showed a dramatic improvement in signal and reproducibility in single-RNA detection in live cells. The same approach was extended to enhancing the stability of engineered fluorescent biosensors containing a fluorescent resonance energy transfer (FRET) pair of fluorescent proteins on which a great majority of systems thus far in the field are based. Using the synonymous modification to FRET biosensors, we achieved correct expression of full-length sensors, eliminating the aberrant truncation products that often were assumed to be due to nonspecific proteolytic cleavages. Importantly, the biological interpretations of the sensor are significantly different when a correct, full-length biosensor is expressed. Thus, we show here a useful and generally applicable method to maintain the integrity of expressed genes, critical for the correct interpretation of probe readouts.

Although important for gene regulation, most studies of genome organization use either fluorescence in situ hybridization (FISH) or chromosome conformation capture (3C) methods. FISH directly visualizes the spatial relationship of sequences but is usually applied to a few loci at a time. The frequency at which sequences are ligated together by formaldehyde cross-linking can be measured genome-wide by 3C methods, with higher frequencies thought to reflect shorter distances. FISH and 3C should therefore give the same views of genome organization, but this has not been tested extensively. We investigated the murine HoxD locus with 3C carbon copy (5C) and FISH in different developmental and activity states and in the presence or absence of epigenetic regulators. We identified situations in which the two data sets are concordant but found other conditions under which chromatin topographies extrapolated from 5C or FISH data are not compatible. We suggest that products captured by 3C do not always reflect spatial proximity, with ligation occurring between sequences located hundreds of nanometers apart, influenced by nuclear environment and chromatin composition. We conclude that results obtained at high resolution with either 3C methods or FISH alone must be interpreted with caution and that views about genome organization should be validated by independent methods.

BRCA1 is a breast and ovarian tumor suppressor. Given its numerous incompletely understood functions and the possibility that more exist, we performed complementary systematic screens in search of new BRCA1 protein-interacting partners. New BRCA1 functions and/or a better understanding of existing ones were sought. Among the new interacting proteins identified, genetic interactions were detected between BRCA1 and four of the interactors: TONSL, SETX, TCEANC, and TCEA2. Genetic interactions were also detected between BRCA1 and certain interactors of TONSL, including both members of the FACT complex. From these results, a new BRCA1 function in the response to transcription-associated DNA damage was detected. Specifically, new roles for BRCA1 in the restart of transcription after UV damage and in preventing or repairing damage caused by stabilized R loops were identified. These roles are likely carried out together with some of the newly identified interactors. This new function may be important in BRCA1 tumor suppression, since the expression of several interactors, including some of the above-noted transcription proteins, is repeatedly aberrant in both breast and ovarian cancers.

DNA replication remains unfinished in many Drosophila polyploid cells, which harbor disproportionately fewer copies of late-replicating chromosomal regions. By analyzing paired-end high-throughput sequence data from polytene larval salivary gland cells, we define 112 underreplicated (UR) euchromatic regions 60–480 kb in size. To determine the effects of underreplication on genome integrity, we analyzed anomalous read pairs and breakpoint reads throughout the euchromatic genome. Each UR euchromatic region contains many different deletions 10–500 kb in size, while very few deletions are present in fully replicated chromosome regions or UR zones from embryo DNA. Thus, during endocycles, stalled forks within UR regions break and undergo local repair instead of remaining stable and generating nested forks. As a result, each salivary gland cell contains hundreds of unique deletions that account for their copy number reductions. Similar UR regions and deletions were observed in ovarian DNA, suggesting that incomplete replication, fork breakage, and repair occur widely in polytene cells. UR regions are enriched in genes encoding immunoglobulin superfamily proteins and contain many neurally expressed and homeotic genes. We suggest that the extensive somatic DNA instability described here underlies position effect variegation, molds the structure of polytene chromosomes, and should be investigated for possible functions.

Tissues may adopt diverse strategies to establish specific transcriptional programs in daughter lineages. In intestinal crypts, enhancers for genes expressed in both major cell types appear broadly permissive in stem and specified progenitor cells. In blood, another self-renewing tissue, it is unclear when chromatin becomes permissive for transcription of genes expressed in distinct terminal lineages. Using chromatin immunoprecipitation (ChIP) combined with deep sequencing (ChIP-seq) to profile activating histone marks, we studied enhancer dynamics in primary mouse blood stem, progenitor, and specified cells. Stem and multipotent progenitor cells show scant H3K4me2 marking at enhancers bound by specific transcription factors in their committed progeny. Rather, enhancers are modulated dynamically and serially, with substantial loss and gain of H3K4me2, at each cellular transition. Quantitative analysis of these dynamics accurately modeled hematopoiesis according to Waddington’s notion of epigenotypes. Delineation of enhancers in terminal blood lineages coincides with cell specification, and enhancers active in single lineages show well-positioned H3K4me2- and H3K27ac-marked nucleosomes and DNaseI hypersensitivity in other cell types, revealing limited lineage fidelity. These findings demonstrate that enhancer chronology in blood cells differs markedly from that in intestinal crypts. Chromatin dynamics in hematopoiesis provide a useful foundation to consider classical observations such as cellular reprogramming and multilineage locus priming.

Long noncoding RNAs (lncRNAs) are thought to play important roles in regulating gene transcription, but few have well-defined expression patterns or known biological functions during mammalian development. Using a conservative pipeline to identify lncRNAs that have important biological functions, we identified 363 lncRNAs in the lung and foregut endoderm. Importantly, we show that these lncRNAs are spatially correlated with transcription factors across the genome. In-depth expression analyses of lncRNAs with genomic loci adjacent to the critical transcription factors Nkx2.1, Gata6, Foxa2 (forkhead box a2), and Foxf1 mimic the expression patterns of their protein-coding neighbor. Loss-of-function analysis demonstrates that two lncRNAs, LL18/NANCI (Nkx2.1-associated noncoding intergenic RNA) and LL34, play distinct roles in endoderm development by controlling expression of critical developmental transcription factors and pathways, including retinoic acid signaling. In particular, we show that LL18/NANCI acts upstream of Nkx2.1 and downstream from Wnt signaling to regulate lung endoderm gene expression. These studies reveal that lncRNAs play an important role in foregut and lung endoderm development by regulating multiple aspects of gene transcription, often through regulation of transcription factor expression.

Despite having high-resolution structures for eukaryotic large ribosomal subunits, it remained unclear how these ribonucleoprotein complexes are constructed in living cells. Nevertheless, knowing where ribosomal proteins interact with ribosomal RNA (rRNA) provides a strategic platform to investigate the connection between spatial and temporal aspects of 60S subunit biogenesis. We previously found that the function of individual yeast large subunit ribosomal proteins (RPLs) in precursor rRNA (pre-rRNA) processing correlates with their location in the structure of mature 60S subunits. This observation suggested that there is an order by which 60S subunits are formed. To test this model, we used proteomic approaches to assay changes in the levels of ribosomal proteins and assembly factors in preribosomes when RPLs functioning in early, middle, and late steps of pre-60S assembly are depleted. Our results demonstrate that structural domains of eukaryotic 60S ribosomal subunits are formed in a hierarchical fashion. Assembly begins at the convex solvent side, followed by the polypeptide exit tunnel, the intersubunit side, and finally the central protuberance. This model provides an initial paradigm for the sequential assembly of eukaryotic 60S subunits. Our results reveal striking differences and similarities between assembly of bacterial and eukaryotic large ribosomal subunits, providing insights into how these RNA–protein particles evolved.

Genetically engineered mouse models (GEMMs) have greatly expanded our knowledge of pancreatic ductal adenocarcinoma (PDAC) and serve as a critical tool to identify and evaluate new treatment strategies. However, the cost and time required to generate conventional pancreatic cancer GEMMs limits their use for investigating novel genetic interactions in tumor development and maintenance. To address this problem, we developed flexible embryonic stem cell (ESC)-based GEMMs that facilitate the rapid generation of genetically defined multiallelic chimeric mice without further strain intercrossing. The ESCs harbor a latent Kras mutant (a nearly ubiquitous feature of pancreatic cancer), a homing cassette, and other genetic elements needed for rapid insertion and conditional expression of tetracycline-controlled transgenes, including fluorescence-coupled shRNAs capable of efficiently silencing gene function by RNAi. This system produces a disease that recapitulates the progression of pancreatic cancer in human patients and enables the study and visualization of the impact of gene perturbation at any stage of pancreas cancer progression. We describe the use of this approach to dissect temporal roles for the tumor suppressor Pten and the oncogene c-Myc in pancreatic cancer development and maintenance.

Fertilization precisely choreographs parental genomes by using gamete-derived cellular factors and activating genome regulatory programs. However, the mechanism remains elusive owing to the technical difficulties of preparing large numbers of high-quality preimplantation cells. Here, we collected >14 × 104 high-quality mouse metaphase II oocytes and used these to establish detailed transcriptional profiles for four early embryo stages and parthenogenetic development. By combining these profiles with other public resources, we found evidence that gene silencing appeared to be mediated in part by noncoding RNAs and that this was a prerequisite for post-fertilization development. Notably, we identified 817 genes that were differentially expressed in embryos after fertilization compared with parthenotes. The regulation of these genes was distinctly different from those expressed in parthenotes, suggesting functional specialization of particular transcription factors prior to first cell cleavage. We identified five transcription factors that were potentially necessary for developmental progression: Foxd1, Nkx2-5, Sox18, Myod1, and Runx1. Our very large-scale whole-transcriptome profile of early mouse embryos yielded a novel and valuable resource for studies in developmental biology and stem cell research. The database is available at http://dbtmee.hgc.jp.

RNAi combined with next-generation sequencing has proven to be a powerful and cost-effective genetic screening platform in mammalian cells. Still, this technology has its limitations and is incompatible with in situ mutagenesis screens on a genome-wide scale. Using p53 as a proof-of-principle target, we readapted the CRISPR (clustered regularly interspaced short palindromic repeats)/Cas9 (CRISPR associated 9) genome-editing system to demonstrate the feasibility of this methodology for targeted gene disruption positive selection assays. By using novel “all-in-one” lentiviral and retroviral delivery vectors heterologously expressing both a codon-optimized Cas9 and its synthetic guide RNA (sgRNA), we show robust selection for the CRISPR-modified Trp53 locus following drug treatment. Furthermore, by linking Cas9 expression to GFP fluorescence, we use an “all-in-one” system to track disrupted Trp53 in chemoresistant lymphomas in the Eμ-myc mouse model. Deep sequencing analysis of the tumor-derived endogenous Cas9-modified Trp53 locus revealed a wide spectrum of mutants that were enriched with seemingly limited off-target effects. Taken together, these results establish Cas9 genome editing as a powerful and practical approach for positive in situ genetic screens.

When adapting to environmental stress, cells attenuate and reprogram their translational output. In part, these altered translation profiles are established through changes in the interactions between RNA-binding proteins and mRNAs. The Argonaute 2 (Ago2)/microRNA (miRNA) machinery has been shown to participate in stress-induced translational up-regulation of a particular mRNA, CAT-1; however, a detailed, transcriptome-wide understanding of the involvement of Ago2 in the process has been lacking. Here, we profiled the overall changes in Ago2–mRNA interactions upon arsenite stress by cross-linking immunoprecipitation (CLIP) followed by high-throughput sequencing (CLIP-seq). Ago2 displayed a significant remodeling of its transcript occupancy, with the majority of 3′ untranslated region (UTR) and coding sequence (CDS) sites exhibiting stronger interaction. Interestingly, target sites that were destined for release from Ago2 upon stress were depleted in miRNA complementarity signatures, suggesting an alternative mode of interaction. To compare the changes in Ago2-binding patterns across transcripts with changes in their translational states, we measured mRNA profiles on ribosome/polysome gradients by RNA sequencing (RNA-seq). Increased Ago2 occupancy correlated with stronger repression of translation for those mRNAs, as evidenced by a shift toward lighter gradient fractions upon stress, while release of Ago2 was associated with the limited number of transcripts that remained translated. Taken together, these data point to a role for Ago2 and the mammalian miRNAs in mediating the translational component of the stress response.

The mammalian circadian timing system consists of a master pacemaker in the suprachiasmatic nucleus (SCN) in the hypothalamus, which is thought to set the phase of slave oscillators in virtually all body cells. However, due to the lack of appropriate in vivo recording technologies, it has been difficult to study how the SCN synchronizes oscillators in peripheral tissues. Here we describe the real-time recording of bioluminescence emitted by hepatocytes expressing circadian luciferase reporter genes in freely moving mice. The technology employs a device dubbed RT-Biolumicorder, which consists of a cylindrical cage with reflecting conical walls that channel photons toward a photomultiplier tube. The monitoring of circadian liver gene expression revealed that hepatocyte oscillators of SCN-lesioned mice synchronized more rapidly to feeding cycles than hepatocyte clocks of intact mice. Hence, the SCN uses signaling pathways that counteract those of feeding rhythms when their phase is in conflict with its own phase.

Fertilization is a crucial yet poorly characterized event in eukaryotes. Our previous discovery that the broadly conserved protein HAP2 (GCS1) functioned in gamete membrane fusion in the unicellular green alga Chlamydomonas and the malaria pathogen Plasmodium led us to exploit the rare biological phenomenon of isogamy in Chlamydomonas in a comparative transcriptomics strategy to uncover additional conserved sexual reproduction genes. All previously identified Chlamydomonas fertilization-essential genes fell into related clusters based on their expression patterns. Out of several conserved genes in a minus gamete cluster, we focused on Cre06.g280600, an ortholog of the fertilization-related Arabidopsis GEX1. Gene disruption, cell biological, and immunolocalization studies show that CrGEX1 functions in nuclear fusion in Chlamydomonas. Moreover, CrGEX1 and its Plasmodium ortholog, PBANKA_113980, are essential for production of viable meiotic progeny in both organisms and thus for mosquito transmission of malaria. Remarkably, we discovered that the genes are members of a large, previously unrecognized family whose first-characterized member, KAR5, is essential for nuclear fusion during yeast sexual reproduction. Our comparative transcriptomics approach provides a new resource for studying sexual development and demonstrates that exploiting the data can lead to the discovery of novel biology that is conserved across distant taxa.

Hypocretin (orexin; Hcrt)-containing neurons of the hypothalamus are essential for the normal regulation of sleep and wake behaviors and have been implicated in feeding, anxiety, depression, and reward. The absence of these neurons causes narcolepsy in humans and model organisms. However, little is known about the molecular phenotype of these cells; previous attempts at comprehensive profiling had only limited sensitivity or were inaccurate. We generated a Hcrt translating ribosome affinity purification (bacTRAP) line for comprehensive translational profiling of all ribosome-bound transcripts in these neurons in vivo. From this profile, we identified >6000 transcripts detectably expressed above background and 188 transcripts that are highly enriched in these neurons, including all known markers of the cells. Blinded analysis of in situ hybridization databases suggests that ∼60% of these are expressed in a Hcrt marker-like pattern. Fifteen of these were confirmed with double labeling and microscopy, including the transcription factor Lhx9. Ablation of this gene results in a >30% loss specifically of Hcrt neurons, without a general disruption of hypothalamic development. Polysomnography and activity monitoring revealed a profound hypersomnolence in these mice. These data provide an in-depth and accurate profile of Hcrt neuron gene expression and suggest that Lhx9 may be important for specification or survival of a subset of these cells.

Antisense oligonucleotides (ASOs) are versatile molecules that can be designed to specifically alter splicing patterns of target pre-mRNAs. Here we exploit this feature to phenocopy a genetic disease. Spinal muscular atrophy (SMA) is a motor neuron disease caused by loss-of-function mutations in the SMN1 gene. The related SMN2 gene expresses suboptimal levels of functional SMN protein due to alternative splicing that skips exon 7; correcting this defect—e.g., with ASOs—is a promising therapeutic approach. We describe the use of ASOs that exacerbate SMN2 missplicing and phenocopy SMA in a dose-dependent manner when administered to transgenic Smn−/− mice. Intracerebroventricular ASO injection in neonatal mice recapitulates SMA-like progressive motor dysfunction, growth impairment, and shortened life span, with α-motor neuron loss and abnormal neuromuscular junctions. These SMA-like phenotypes are prevented by a therapeutic ASO that restores correct SMN2 splicing. We uncovered starvation-induced splicing changes, particularly in SMN2, which likely accelerate disease progression. These results constitute proof of principle that ASOs designed to cause sustained splicing defects can be used to induce pathogenesis and rapidly and accurately model splicing-associated diseases in animals. This approach allows the dissection of pathogenesis mechanisms, including spatial and temporal features of disease onset and progression, as well as testing of candidate therapeutics.

Epigenetic memory mediated by Polycomb group (PcG) proteins must be maintained during cell division, but must also be flexible to allow cell fate transitions. Here we quantify dynamic chromatin-binding properties of PH::GFP and PC::GFP in living Drosophila in two cell types that undergo defined differentiation and mitosis events. Quantitative fluorescence recovery after photobleaching (FRAP) analysis demonstrates that PcG binding has a higher plasticity in stem cells than in more determined cells and identifies a fraction of PcG proteins that binds mitotic chromatin with up to 300-fold longer residence times than in interphase. Mathematical modeling examines which parameters best distinguish stem cells from differentiated cells. We identify phosphorylation of histone H3 at Ser 28 as a potential mechanism governing the extent and rate of mitotic PC dissociation in different lineages. We propose that regulation of the kinetic properties of PcG–chromatin binding is an essential factor in the choice between stability and flexibility in the establishment of cell identities.

Legumes and many nonleguminous plants enter symbiotic interactions with microbes, and it is poorly understood how host plants respond to promote beneficial, symbiotic microbial interactions while suppressing those that are deleterious or pathogenic. Trans-acting siRNAs (tasiRNAs) negatively regulate target transcripts and are characterized by siRNAs spaced in 21-nucleotide (nt) “phased” intervals, a pattern formed by DICER-LIKE 4 (DCL4) processing. A search for phased siRNAs (phasiRNAs) found at least 114 Medicago loci, the majority of which were defense-related NB-LRR-encoding genes. We identified three highly abundant 22-nt microRNA (miRNA) families that target conserved domains in these NB-LRRs and trigger the production of trans-acting siRNAs. High levels of small RNAs were matched to >60% of all ∼540 encoded Medicago NB-LRRs; in the potato, a model for mycorrhizal interactions, phasiRNAs were also produced from NB-LRRs. DCL2 and SGS3 transcripts were also cleaved by these 22-nt miRNAs, generating phasiRNAs, suggesting synchronization between silencing and pathogen defense pathways. In addition, a new example of apparent “two-hit” phasiRNA processing was identified. Our data reveal complex tasiRNA-based regulation of NB-LRRs that potentially evolved to facilitate symbiotic interactions and demonstrate miRNAs as master regulators of a large gene family via the targeting of highly conserved, protein-coding motifs, a new paradigm for miRNA function.