A data analysis framework for combining multiple batches increases the power of isobaric proteomics experiments.

We present a framework for the analysis of multiplexed mass spectrometry proteomics data that reduces estimation error when combining multiple isobaric batches. Variations in the number and quality of observations have long complicated the analysis of isobaric proteomics data. Here we show that the power to detect statistical associations is substantially improved by utilizing models that directly account for known sources of variation in the number and quality of observations that occur across batches.In a multibatch benchmarking experiment, our open-source software (msTrawler) increases the power to detect changes, especially in the range of less than twofold changes, while simultaneously increasing quantitative proteome coverage by utilizing more low-signal observations. Further analyses of previously published multiplexed datasets of 4 and 23 batches highlight both increased power and the ability to navigate complex missing data patterns without relying on unverifiable imputations or discarding reliable measurements.

rDNA array length is a major determinant of replicative lifespan in budding yeast.

The complex processes and interactions that regulate aging and determine lifespan are not fully defined for any organism. Here, taking advantage of recent technological advances in studying aging in budding yeast, we discovered a previously unappreciated relationship between the number of copies of the ribosomal RNA gene present in its chromosomal array and replicative lifespan (RLS). Specifically, the chromosomal ribosomal DNA (rDNA) copy number (rDNA CN) positively correlated with RLS and this interaction explained over 70% of variability in RLS among a series of wild-type strains. In strains with low rDNA CN, SIR2 expression was attenuated and extrachromosomal rDNA circle (ERC) accumulation was increased, leading to shorter lifespan. Suppressing ERC formation by deletion of FOB1 eliminated the relationship between rDNA CN and RLS. These data suggest that previously identified rDNA CN regulatory mechanisms limit lifespan. Importantly, the RLSs of reported lifespan-enhancing mutations were significantly impacted by rDNA CN, suggesting that changes in rDNA CN might explain the magnitude of some of those reported effects. We propose that because rDNA CN is modulated by environmental, genetic, and stochastic factors, considering rDNA CN is a prerequisite for accurate interpretation of lifespan data.

Novel insights from a multiomics dissection of the Hayflick limit.

The process wherein dividing cells exhaust proliferative capacity and enter into replicative senescence has become a prominent model for cellular aging in vitro. Despite decades of study, this cellular state is not fully understood in culture and even much less so during aging. Here, we revisit Leonard Hayflick's original observation of replicative senescence in WI-38 human lung fibroblasts equipped with a battery of modern techniques including RNA-seq, single-cell RNA-seq, proteomics, metabolomics, and ATAC-seq. We find evidence that the transition to a senescent state manifests early, increases gradually, and corresponds to a concomitant global increase in DNA accessibility in nucleolar and lamin associated domains. Furthermore, we demonstrate that senescent WI-38 cells acquire a striking resemblance to myofibroblasts in a process similar to the epithelial to mesenchymal transition (EMT) that is regulated by t YAP1/TEAD1 and TGF-ß2. Lastly, we show that verteporfin inhibition of YAP1/TEAD1 activity in aged WI-38 cells robustly attenuates this gene expression program.

Loss of major nutrient sensing and signaling pathways suppresses starvation lethality in electron transport chain mutants.

The electron transport chain (ETC) is a well-studied and highly conserved metabolic pathway that produces ATP through generation of a proton gradient across the inner mitochondrial membrane coupled to oxidative phosphorylation. ETC mutations are associated with a wide array of human disease conditions and to aging-related phenotypes in a number of different organisms. In this study, we sought to better understand the role of the ETC in aging using a yeast model. A panel of ETC mutant strains that fail to survive starvation was used to isolate suppressor mutants that survive. These suppressors tend to fall into major nutrient sensing and signaling pathways, suggesting that the ETC is involved in proper starvation signaling to these pathways in yeast. These suppressors also partially restore ETC-associated gene expression and pH homeostasis defects, though it remains unclear whether these phenotypes directly cause the suppression or are simply effects. This work further highlights the complex cellular network connections between metabolic pathways and signaling events in the cell and their potential roles in aging and age-related diseases.

Differentiation reveals latent features of aging and an energy barrier in murine myogenesis.

Skeletal muscle experiences a decline in lean mass and regenerative potential with age, in part due to intrinsic changes in progenitor cells. However, it remains unclear how age-related changes in progenitors manifest across a differentiation trajectory. Here, we perform single-cell RNA sequencing (RNA-seq) on muscle mononuclear cells from young and aged mice and profile muscle stem cells (MuSCs) and fibro-adipose progenitors (FAPs) after differentiation. Differentiation increases the magnitude of age-related change in MuSCs and FAPs, but it also masks a subset of age-related changes present in progenitors. Using a dynamical systems approach and RNA velocity, we find that aged MuSCs follow the same differentiation trajectory as young cells but stall in differentiation near a commitment decision. Our results suggest that differentiation reveals latent features of aging and that fate commitment decisions are delayed in aged myogenic cells in vitro.

Solo: Doublet Identification in Single-Cell RNA-Seq via Semi-Supervised Deep Learning.

Single-cell RNA sequencing (scRNA-seq) measurements of gene expression enable an unprecedented high-resolution view into cellular state. However, current methods often result in two or more cells that share the same cell-identifying barcode; these "doublets" violate the fundamental premise of single-cell technology and can lead to incorrect inferences. Here, we describe Solo, a semi-supervised deep learning approach that identifies doublets with greater accuracy than existing methods. Solo embeds cells unsupervised using a variational autoencoder and then appends a feed-forward neural network layer to the encoder to form a supervised classifier. We train this classifier to distinguish simulated doublets from the observed data. Solo can be applied in combination with experimental doublet detection methods to further purify scRNA-seq data to true single cells. It is freely available from https://github.com/calico/solo. A record of this paper's transparent peer review process is included in the Supplemental Information.

A tps1Δ persister-like state in Saccharomyces cerevisiae is regulated by MKT1.

Trehalose metabolism in yeast has been linked to a variety of phenotypes, including heat resistance, desiccation tolerance, carbon-source utilization, and sporulation. The relationships among the several phenotypes of mutants unable to synthesize trehalose are not understood, even though the pathway is highly conserved. One of these phenotypes is that tps1Δ strains cannot reportedly grow on media containing glucose or fructose, even when another carbon source they can use (e.g. galactose) is present. Here we corroborate the recent observation that a small fraction of yeast tps1Δ cells do grow on glucose, unlike the majority of the population. This is not due to a genetic alteration, but instead resembles the persister phenotype documented in many microorganisms and cancer cells undergoing lethal stress. We extend these observations to show that this phenomenon is glucose-specific, as it does not occur on another highly fermented carbon source, fructose. We further demonstrate that this phenomenon appears to be related to mitochondrial complex III function, but unrelated to inorganic phosphate levels in the cell, as had previously been suggested. Finally, we found that this phenomenon is specific to S288C-derived strains, and is the consequence of a variant in the MKT1 gene.

Learning causal networks using inducible transcription factors and transcriptome-wide time series.

We present IDEA (the Induction Dynamics gene Expression Atlas), a dataset constructed by independently inducing hundreds of transcription factors (TFs) and measuring timecourses of the resulting gene expression responses in budding yeast. Each experiment captures a regulatory cascade connecting a single induced regulator to the genes it causally regulates. We discuss the regulatory cascade of a single TF, Aft1, in detail; however, IDEA contains > 200 TF induction experiments with 20 million individual observations and 100,000 signal-containing dynamic responses. As an application of IDEA, we integrate all timecourses into a whole-cell transcriptional model, which is used to predict and validate multiple new and underappreciated transcriptional regulators. We also find that the magnitudes of coefficients in this model are predictive of genetic interaction profile similarities. In addition to being a resource for exploring regulatory connectivity between TFs and their target genes, our modeling approach shows that combining rapid perturbations of individual genes with genome-scale time-series measurements is an effective strategy for elucidating gene regulatory networks.

Murine single-cell RNA-seq reveals cell-identity- and tissue-specific trajectories of aging.

Aging is a pleiotropic process affecting many aspects of mammalian physiology. Mammals are composed of distinct cell type identities and tissue environments, but the influence of these cell identities and environments on the trajectory of aging in individual cells remains unclear. Here, we performed single-cell RNA-seq on >50,000 individual cells across three tissues in young and old mice to allow for direct comparison of aging phenotypes across cell types. We found transcriptional features of aging common across many cell types, as well as features of aging unique to each type. Leveraging matrix factorization and optimal transport methods, we found that both cell identities and tissue environments exert influence on the trajectory and magnitude of aging, with cell identity influence predominating. These results suggest that aging manifests with unique directionality and magnitude across the diverse cell identities in mammals.

Simultaneous Profiling of DNA Accessibility and Gene Expression Dynamics with ATAC-Seq and RNA-Seq.

This chapter describes sequencing-based methods for profiling dynamic changes in DNA accessibility and gene expression in Saccharomyces cerevisiae. Assay for transposase-accessible chromatin with high-throughput sequencing (ATAC-Seq) is a powerful technique for identifying nucleosome-free regions of the genome. Combining ATAC-Seq with RNA Sequencing (RNA-Seq) is a rapid approach for studying the relationship between genome structure and changes in global patterns of gene expression from a single experiment. A laboratory protocol is presented for these methods as well as examples of typical results and visualizations.

A new experimental platform facilitates assessment of the transcriptional and chromatin landscapes of aging yeast.

Replicative aging of Saccharomyces cerevisiae is an established model system for eukaryotic cellular aging. A limitation in yeast lifespan studies has been the difficulty of separating old cells from young cells in large quantities. We engineered a new platform, the Miniature-chemostat Aging Device (MAD), that enables purification of aged cells at sufficient quantities for genomic and biochemical characterization of aging yeast populations. Using MAD, we measured DNA accessibility and gene expression changes in aging cells. Our data highlight an intimate connection between aging, growth rate, and stress. Stress-independent genes that change with age are highly enriched for targets of the signal recognition particle (SRP). Combining MAD with an improved ATAC-seq method, we find that increasing proteasome activity reduces rDNA instability usually observed in aging cells and, contrary to published findings, provide evidence that global nucleosome occupancy does not change significantly with age.

An estradiol-inducible promoter enables fast, graduated control of gene expression in fission yeast.

The fission yeast Schizosaccharomyces pombe lacks a diverse toolkit of inducible promoters for experimental manipulation. Available inducible promoters suffer from slow induction kinetics, limited control of expression levels and/or a requirement for defined growth medium. In particular, no S. pombe inducible promoter systems exhibit a linear dose-response, which would allow expression to be tuned to specific levels. We have adapted a fast, orthogonal promoter system with a large dynamic range and a linear dose response, based on ß-estradiol-regulated function of the human oestrogen receptor, for use in S. pombe. We show that this promoter system, termed Z3 EV, turns on quickly, can reach a maximal induction of 20-fold, and exhibits a linear dose response over its entire induction range, with few off-target effects. We demonstrate the utility of this system by regulating the mitotic inhibitor Wee1 to create a strain in which cell size is regulated by ß-estradiol concentration. This promoter system will be of great utility for experimentally regulating gene expression in fission yeast. Copyright © 2017 John Wiley & Sons, Ltd.

LINC00520 is induced by Src, STAT3, and PI3K and plays a functional role in breast cancer.

Long non-coding RNAs (lncRNAs) have been implicated in normal cellular homeostasis as well as pathophysiological conditions, including cancer. Here we performed global gene expression profiling of mammary epithelial cells transformed by oncogenic v-Src, and identified a large subset of uncharacterized lncRNAs potentially involved in breast cancer development. Specifically, our analysis revealed a novel lncRNA, LINC00520 that is upregulated upon ectopic expression of oncogenic v-Src, in a manner that is dependent on the transcription factor STAT3. Similarly, LINC00520 is also increased in mammary epithelial cells transformed by oncogenic PI3K and its expression is decreased upon knockdown of mutant PIK3CA. Additional expression profiling highlight that LINC00520 is elevated in a subset of human breast carcinomas, with preferential enrichment in the basal-like molecular subtype. ShRNA-mediated depletion of LINC00520 results in decreased cell migration and loss of invasive structures in 3D. RNA sequencing analysis uncovers several genes that are differentially expressed upon ectopic expression of LINC00520, a significant subset of which are also induced in v-Src-transformed MCF10A cells. Together, these findings characterize LINC00520 as a lncRNA that is regulated by oncogenic Src, PIK3CA and STAT3, and which may contribute to the molecular etiology of breast cancer.

Genome-wide RNA-Seq of Human Motor Neurons Implicates Selective ER Stress Activation in Spinal Muscular Atrophy.

Spinal muscular atrophy (SMA) is caused by mutations in the SMN1 gene. Because this gene is expressed ubiquitously, it remains poorly understood why motor neurons (MNs) are one of the most affected cell types. To address this question, we carried out RNA sequencing studies using fixed, antibody-labeled, and purified MNs produced from control and SMA patient-derived induced pluripotent stem cells (iPSCs). We found SMA-specific changes in MNs, including hyper-activation of the ER stress pathway. Functional studies demonstrated that inhibition of ER stress improves MN survival in vitro even in MNs expressing low SMN. In SMA mice, systemic delivery of an ER stress inhibitor that crosses the blood-brain barrier led to the preservation of spinal cord MNs. Therefore, our study implies that selective activation of ER stress underlies MN death in SMA. Moreover, the approach we have taken would be broadly applicable to the study of disease-prone human cells in heterogeneous cultures.

miR-142 regulates the tumorigenicity of human breast cancer stem cells through the canonical WNT signaling pathway.

MicroRNAs (miRNAs) are important regulators of stem and progenitor cell functions. We previously reported that miR-142 and miR-150 are upregulated in human breast cancer stem cells (BCSCs) as compared to the non-tumorigenic breast cancer cells. In this study, we report that miR-142 efficiently recruits the APC mRNA to an RNA-induced silencing complex, activates the canonical WNT signaling pathway in an APC-suppression dependent manner, and activates the expression of miR-150. Enforced expression of miR-142 or miR-150 in normal mouse mammary stem cells resulted in the regeneration of hyperproliferative mammary glands in vivo. Knockdown of endogenous miR-142 effectively suppressed organoid formation by BCSCs and slowed tumor growth initiated by human BCSCs in vivo. These results suggest that in some tumors, miR-142 regulates the properties of BCSCs at least in part by activating the WNT signaling pathway and miR-150 expression.

Topological organization of multichromosomal regions by the long intergenic noncoding RNA Firre.

RNA, including long noncoding RNA (lncRNA), is known to be an abundant and important structural component of the nuclear matrix. However, the molecular identities, functional roles and localization dynamics of lncRNAs that influence nuclear architecture remain poorly understood. Here, we describe one lncRNA, Firre, that interacts with the nuclear-matrix factor hnRNPU through a 156-bp repeating sequence and localizes across an ~5-Mb domain on the X chromosome. We further observed Firre localization across five distinct trans-chromosomal loci, which reside in spatial proximity to the Firre genomic locus on the X chromosome. Both genetic deletion of the Firre locus and knockdown of hnRNPU resulted in loss of colocalization of these trans-chromosomal interacting loci. Thus, our data suggest a model in which lncRNAs such as Firre can interface with and modulate nuclear architecture across chromosomes.

Transposable elements modulate human RNA abundance and splicing via specific RNA-protein interactions.

BACKGROUND: Transposable elements (TEs) have significantly influenced the evolution of transcriptional regulatory networks in the human genome. Post-transcriptional regulation of human genes by TE-derived sequences has been observed in specific contexts, but has yet to be systematically and comprehensively investigated. Here, we study a collection of 75 CLIP-Seq experiments mapping the RNA binding sites for a diverse set of 51 human proteins to explore the role of TEs in post-transcriptional regulation of human mRNAs and lncRNAs via RNA-protein interactions. RESULTS: We detect widespread interactions between RNA binding proteins (RBPs) and many families of TE-derived sequence in the CLIP-Seq data. Further, alignment coverage peaks on specific positions of the TE consensus sequences, illuminating a diversity of TE-specific RBP binding motifs. Evidence of binding and conservation of these motifs in the nonrepetitive transcriptome suggests that TEs have generally appropriated existing sequence preferences of the RBPs. Depletion assays for numerous RBPs show that TE-derived binding sites affect transcript abundance and splicing similarly to nonrepetitive sites. However, in a few cases the effect of RBP binding depends on the specific TE family bound; for example, the ubiquitously expressed RBP HuR confers transcript stability unless bound to an Alu element. CONCLUSIONS: Our meta-analysis suggests a widespread role for TEs in shaping RNA-protein regulatory networks in the human genome.

Long noncoding RNAs regulate adipogenesis.

The prevalence of obesity has led to a surge of interest in understanding the detailed mechanisms underlying adipocyte development. Many protein-coding genes, mRNAs, and microRNAs have been implicated in adipocyte development, but the global expression patterns and functional contributions of long noncoding RNA (lncRNA) during adipogenesis have not been explored. Here we profiled the transcriptome of primary brown and white adipocytes, preadipocytes, and cultured adipocytes and identified 175 lncRNAs that are specifically regulated during adipogenesis. Many lncRNAs are adipose-enriched, strongly induced during adipogenesis, and bound at their promoters by key transcription factors such as peroxisome proliferator-activated receptor γ (PPARγ) and CCAAT/enhancer-binding protein α (CEBPα). RNAi-mediated loss of function screens identified functional lncRNAs with varying impact on adipogenesis. Collectively, we have identified numerous lncRNAs that are functionally required for proper adipogenesis.

Differential analysis of gene regulation at transcript resolution with RNA-seq.

Differential analysis of gene and transcript expression using high-throughput RNA sequencing (RNA-seq) is complicated by several sources of measurement variability and poses numerous statistical challenges. We present Cuffdiff 2, an algorithm that estimates expression at transcript-level resolution and controls for variability evident across replicate libraries. Cuffdiff 2 robustly identifies differentially expressed transcripts and genes and reveals differential splicing and promoter-preference changes. We demonstrate the accuracy of our approach through differential analysis of lung fibroblasts in response to loss of the developmental transcription factor HOXA1, which we show is required for lung fibroblast and HeLa cell cycle progression. Loss of HOXA1 results in significant expression level changes in thousands of individual transcripts, along with isoform switching events in key regulators of the cell cycle. Cuffdiff 2 performs robust differential analysis in RNA-seq experiments at transcript resolution, revealing a layer of regulation not readily observable with other high-throughput technologies.

Concordant regulation of translation and mRNA abundance for hundreds of targets of a human microRNA.

MicroRNAs (miRNAs) regulate gene expression posttranscriptionally by interfering with a target mRNA's translation, stability, or both. We sought to dissect the respective contributions of translational inhibition and mRNA decay to microRNA regulation. We identified direct targets of a specific miRNA, miR-124, by virtue of their association with Argonaute proteins, core components of miRNA effector complexes, in response to miR-124 transfection in human tissue culture cells. In parallel, we assessed mRNA levels and obtained translation profiles using a novel global approach to analyze polysomes separated on sucrose gradients. Analysis of translation profiles for approximately 8,000 genes in these proliferative human cells revealed that basic features of translation are similar to those previously observed in rapidly growing Saccharomyces cerevisiae. For approximately 600 mRNAs specifically recruited to Argonaute proteins by miR-124, we found reductions in both the mRNA abundance and inferred translation rate spanning a large dynamic range. The changes in mRNA levels of these miR-124 targets were larger than the changes in translation, with average decreases of 35% and 12%, respectively. Further, there was no identifiable subgroup of mRNA targets for which the translational response was dominant. Both ribosome occupancy (the fraction of a given gene's transcripts associated with ribosomes) and ribosome density (the average number of ribosomes bound per unit length of coding sequence) were selectively reduced for hundreds of miR-124 targets by the presence of miR-124. Changes in protein abundance inferred from the observed changes in mRNA abundance and translation profiles closely matched changes directly determined by Western analysis for 11 of 12 proteins, suggesting that our assays captured most of miR-124-mediated regulation. These results suggest that miRNAs inhibit translation initiation or stimulate ribosome drop-off preferentially near the start site and are not consistent with inhibition of polypeptide elongation, or nascent polypeptide degradation contributing significantly to miRNA-mediated regulation in proliferating HEK293T cells. The observation of concordant changes in mRNA abundance and translational rate for hundreds of miR-124 targets is consistent with a functional link between these two regulatory outcomes of miRNA targeting, and the well-documented interrelationship between translation and mRNA decay.