Biologically indeterminate yet ordered promiscuous gene expression in single medullary thymic epithelial cells.

To induce central T-cell tolerance, medullary thymic epithelial cells (mTEC) collectively express most protein-coding genes, thereby presenting an extensive library of tissue-restricted antigens (TRAs). To resolve mTEC diversity and whether promiscuous gene expression (PGE) is stochastic or coordinated, we sequenced transcriptomes of 6,894 single mTEC, enriching for 1,795 rare cells expressing either of two TRAs, TSPAN8 or GP2. Transcriptional heterogeneity allowed partitioning of mTEC into 15 reproducible subpopulations representing distinct maturational trajectories, stages and subtypes, including novel mTEC subsets, such as chemokine-expressing and ciliated TEC, which warrant further characterisation. Unexpectedly, 50 modules of genes were robustly defined each showing patterns of co-expression within individual cells, which were mainly not explicable by chromosomal location, biological pathway or tissue specificity. Further, TSPAN8+ and GP2+ mTEC were randomly dispersed within thymic medullary islands. Consequently, these data support observations that PGE exhibits ordered co-expression, although mechanisms underlying this instruction remain biologically indeterminate. Ordered co-expression and random spatial distribution of a diverse range of TRAs likely enhance their presentation and encounter with passing thymocytes, while maintaining mTEC identity.

An integrative transcriptomics approach identifies miR-503 as a candidate master regulator of the estrogen response in MCF-7 breast cancer cells.

Estrogen receptor α (ERα) is an important biomarker of breast cancer severity and a common therapeutic target. In response to estrogen, ERα stimulates a dynamic transcriptional program including both coding and noncoding RNAs. We generate a fine-scale map of expression dynamics by performing a temporal profiling of both messenger RNAs (mRNAs) and microRNAs (miRNAs) in MCF-7 cells (an ER+ model cell line for breast cancer) in response to estrogen stimulation. We identified three primary expression trends-transient, induced, and repressed-that were each enriched for genes with distinct cellular functions. Integrative analysis of mRNA and miRNA temporal expression profiles identified miR-503 as the strongest candidate master regulator of the estrogen response, in part through suppression of ZNF217-an oncogene that is frequently amplified in cancer. We confirmed experimentally that miR-503 directly targets ZNF217 and that overexpression of miR-503 suppresses MCF-7 cell proliferation. Moreover, the levels of ZNF217 and miR-503 are associated with opposite outcomes in breast cancer patient cohorts, with high expression of ZNF217 associated with poor survival and high expression of miR-503 associated with improved survival. Overall, these data indicate that miR-503 acts as a potent estrogen-induced candidate tumor suppressor miRNA that opposes cellular proliferation and has promise as a novel therapeutic for breast cancer. More generally, our work provides a systems-level framework for identifying functional interactions that shape the temporal dynamics of gene expression.

Pseudogenes transcribed in breast invasive carcinoma show subtype-specific expression and ceRNA potential.

BACKGROUND: Recent studies have shown that some pseudogenes are transcribed and contribute to cancer when dysregulated. In particular, pseudogene transcripts can function as competing endogenous RNAs (ceRNAs). The high similarity of gene and pseudogene nucleotide sequence has hindered experimental investigation of these mechanisms using RNA-seq. Furthermore, previous studies of pseudogenes in breast cancer have not integrated miRNA expression data in order to perform large-scale analysis of ceRNA potential. Thus, knowledge of both pseudogene ceRNA function and the role of pseudogene expression in cancer are restricted to isolated examples. RESULTS: To investigate whether transcribed pseudogenes play a pervasive regulatory role in cancer, we developed a novel bioinformatic method for measuring pseudogene transcription from RNA-seq data. We applied this method to 819 breast cancer samples from The Cancer Genome Atlas (TCGA) project. We then clustered the samples using pseudogene expression levels and integrated sample-paired pseudogene, gene and miRNA expression data with miRNA target prediction to determine whether more pseudogenes have ceRNA potential than expected by chance. CONCLUSIONS: Our analysis identifies with high confidence a set of 440 pseudogenes that are transcribed in breast cancer tissue. Of this set, 309 pseudogenes exhibit significant differential expression among breast cancer subtypes. Hierarchical clustering using only pseudogene expression levels accurately separates tumor samples from normal samples and discriminates the Basal subtype from the Luminal and Her2 subtypes. Correlation analysis shows more positively correlated pseudogene-parent gene pairs and negatively correlated pseudogene-miRNA pairs than expected by chance. Furthermore, 177 transcribed pseudogenes possess binding sites for co-expressed miRNAs that are also predicted to target their parent genes. Taken together, these results increase the catalog of putative pseudogene ceRNAs and suggest that pseudogene transcription in breast cancer may play a larger role than previously appreciated.

Identification of microRNAs associated with allergic airway disease using a genetically diverse mouse population.

BACKGROUND: Allergic airway diseases (AADs) such as asthma are characterized in part by granulocytic airway inflammation. The gene regulatory networks that govern granulocyte recruitment are poorly understood, but evidence is accruing that microRNAs (miRNAs) play an important role. To identify miRNAs that may underlie AADs, we used two complementary approaches that leveraged the genotypic and phenotypic diversity of the Collaborative Cross (CC) mouse population. In the first approach, we sought to identify miRNA expression quantitative trait loci (eQTL) that overlap QTL for AAD-related phenotypes. Specifically, CC founder strains and incipient lines of the CC were sensitized and challenged with house dust mite allergen followed by measurement of granulocyte recruitment to the lung. Total lung RNA was isolated and miRNA was measured using arrays for CC founders and qRT-PCR for incipient CC lines. RESULTS: Among CC founders, 92 miRNAs were differentially expressed. We measured the expression of 40 of the most highly expressed of these 92 miRNAs in the incipient lines of the CC and identified 18 eQTL corresponding to 14 different miRNAs. Surprisingly, half of these eQTL were distal to the corresponding miRNAs, and even on different chromosomes. One of the largest-effect local miRNA eQTL was for miR-342-3p, for which we identified putative causal variants by bioinformatic analysis of the effects of single nucleotide polymorphisms on RNA structure. None of the miRNA eQTL co-localized with QTL for eosinophil or neutrophil recruitment. In the second approach, we constructed putative miRNA/mRNA regulatory networks and identified three miRNAs (miR-497, miR-351 and miR-31) as candidate master regulators of genes associated with neutrophil recruitment. Analysis of a dataset from human keratinocytes transfected with a miR-31 inhibitor revealed two target genes in common with miR-31 targets correlated with neutrophils, namely Oxsr1 and Nsf. CONCLUSIONS: miRNA expression in the allergically inflamed murine lung is regulated by genetic loci that are smaller in effect size compared to mRNA eQTL and often act in trans. Thus our results indicate that the genetic architecture of miRNA expression is different from mRNA expression. We identified three miRNAs, miR-497, miR-351 and miR-31, that are candidate master regulators of genes associated with neutrophil recruitment. Because miR-31 is expressed in airway epithelia and is predicted to target genes with known links to neutrophilic inflammation, we suggest that miR-31 is a potentially novel regulator of airway inflammation.

miR-29 is an important driver of aging-related phenotypes.

Aging is a consequence of complex molecular changes, but whether a single microRNA (miRNA) can drive aging remains unclear. A miRNA known to be upregulated during both normal and premature aging is miR-29. We find miR-29 to also be among the top miRNAs predicted to drive aging-related gene expression changes. We show that partial loss of miR-29 extends the lifespan of Zmpste24-/- mice, an established model of progeria, indicating that miR-29 is functionally important in this accelerated aging model. To examine whether miR-29 alone is sufficient to promote aging-related phenotypes, we generated mice in which miR-29 can be conditionally overexpressed (miR-29TG). miR-29 overexpression is sufficient to drive many aging-related phenotypes and led to early lethality. Transcriptomic analysis of both young miR-29TG and old WT mice reveals shared downregulation of genes associated with extracellular matrix organization and fatty acid metabolism, and shared upregulation of genes in pathways linked to inflammation. These results highlight the functional importance of miR-29 in controlling a gene expression program that drives aging-related phenotypes.

Complexity of microRNA function and the role of isomiRs in lipid homeostasis.

MicroRNAs (miRNAs) are key posttranscriptional regulators of biological pathways that govern lipid metabolic phenotypes. Recent advances in high-throughput small RNA sequencing technology have revealed the complex and dynamic repertoire of miRNAs. Specifically, it has been demonstrated that a single genomic locus can give rise to multiple, functionally distinct miRNA isoforms (isomiR). There are several mechanisms by which isomiRs can be generated, including processing heterogeneity and posttranscriptional modifications, such as RNA editing, exonuclease-mediated nucleotide trimming, and/or nontemplated nucleotide addition (NTA). NTAs are dominant at the 3'-end of a miRNA, are most commonly uridylation or adenlyation events, and are catalyzed by one or more of several nucleotidyl transferase enzymes. 3' NTAs can affect miRNA stability and/or activity and are physiologically regulated, whereas modifications to the 5'-ends of miRNAs likely alter miRNA targeting activity. Recent evidence also suggests that the biogenesis of specific miRNAs, or small RNAs that act as miRNAs, can occur through unconventional mechanisms that circumvent key canonical miRNA processing steps. The unveiling of miRNA diversity has significantly added to our view of the complexity of miRNA function. In this review we present the current understanding of the biological relevance of isomiRs and their potential role in regulating lipid metabolism.

Ageing compromises mouse thymus function and remodels epithelial cell differentiation.

Ageing is characterised by cellular senescence, leading to imbalanced tissue maintenance, cell death and compromised organ function. This is first observed in the thymus, the primary lymphoid organ that generates and selects T cells. However, the molecular and cellular mechanisms underpinning these ageing processes remain unclear. Here, we show that mouse ageing leads to less efficient T cell selection, decreased self-antigen representation and increased T cell receptor repertoire diversity. Using a combination of single-cell RNA-seq and lineage-tracing, we find that progenitor cells are the principal targets of ageing, whereas the function of individual mature thymic epithelial cells is compromised only modestly. Specifically, an early-life precursor cell population, retained in the mouse cortex postnatally, is virtually extinguished at puberty. Concomitantly, a medullary precursor cell quiesces, thereby impairing maintenance of the medullary epithelium. Thus, ageing disrupts thymic progenitor differentiation and impairs the core immunological functions of the thymus.

Experimental design for single-cell RNA sequencing.

Single-cell RNA sequencing (scRNA-seq) has opened new avenues for the characterization of heterogeneity in a large variety of cellular systems. As this is a relatively new technique, the field is fast evolving. Here, we discuss general considerations in experimental design and the two most popular approaches, plate-based Smart-Seq2 and microdroplet-based scRNA-seq at the example of 10x Chromium. We discuss advantages and disadvantages of both methods and point out major factors to consider in designing successful experiments.

miRquant 2.0: an Expanded Tool for Accurate Annotation and Quantification of MicroRNAs and their isomiRs from Small RNA-Sequencing Data.

Small non-coding RNAs, in particular microRNAs, are critical for normal physiology and are candidate biomarkers, regulators, and therapeutic targets for a wide variety of diseases. There is an ever-growing interest in the comprehensive and accurate annotation of microRNAs across diverse cell types, conditions, species, and disease states. Highthroughput sequencing technology has emerged as the method of choice for profiling microRNAs. Specialized bioinformatic strategies are required to mine as much meaningful information as possible from the sequencing data to provide a comprehensive view of the microRNA landscape. Here we present miRquant 2.0, an expanded bioinformatics tool for accurate annotation and quantification of microRNAs and their isoforms (termed isomiRs) from small RNA-sequencing data. We anticipate that miRquant 2.0 will be useful for researchers interested not only in quantifying known microRNAs but also mining the rich well of additional information embedded in small RNA-sequencing data.

Essential Function of Dicer in Resolving DNA Damage in the Rapidly Dividing Cells of the Developing and Malignant Cerebellum.

Maintenance of genomic integrity is critical during neurodevelopment, particularly in rapidly dividing cerebellar granule neuronal precursors that experience constitutive replication-associated DNA damage. As Dicer was recently recognized to have an unexpected function in the DNA damage response, we examined whether Dicer was important for preserving genomic integrity in the developing brain. We report that deletion of Dicer in the developing mouse cerebellum resulted in the accumulation of DNA damage leading to cerebellar progenitor degeneration, which was rescued with p53 deficiency; deletion of DGCR8 also resulted in similar DNA damage and cerebellar degeneration. Dicer deficiency also resulted in DNA damage and death in other rapidly dividing cells including embryonic stem cells and the malignant cerebellar progenitors in a mouse model of medulloblastoma. Together, these results identify an essential function of Dicer in resolving the spontaneous DNA damage that occurs during the rapid proliferation of developmental progenitors and malignant cells.

Addressing Bias in Small RNA Library Preparation for Sequencing: A New Protocol Recovers MicroRNAs that Evade Capture by Current Methods.

Recent advances in sequencing technology have helped unveil the unexpected complexity and diversity of small RNAs. A critical step in small RNA library preparation for sequencing is the ligation of adapter sequences to both the 5' and 3' ends of small RNAs. Studies have shown that adapter ligation introduces a significant but widely unappreciated bias in the results of high-throughput small RNA sequencing. We show that due to this bias the two widely used Illumina library preparation protocols produce strikingly different microRNA (miRNA) expression profiles in the same batch of cells. There are 102 highly expressed miRNAs that are >5-fold differentially detected and some miRNAs, such as miR-24-3p, are over 30-fold differentially detected. While some level of bias in library preparation is not surprising, the apparent massive differential bias between these two widely used adapter sets is not well appreciated. In an attempt to mitigate this bias, the new Bioo Scientific NEXTflex V2 protocol utilizes a pool of adapters with random nucleotides at the ligation boundary. We show that this protocol is able to detect robustly several miRNAs that evade capture by the Illumina-based methods. While these analyses do not indicate a definitive gold standard for small RNA library preparation, the results of the NEXTflex protocol do correlate best with RT-qPCR. As increasingly more laboratories seek to study small RNAs, researchers should be aware of the extent to which the results may differ with different protocols, and should make an informed decision about the protocol that best fits their study.

Small tRNA-derived RNAs are increased and more abundant than microRNAs in chronic hepatitis B and C.

Persistent infections with hepatitis B virus (HBV) or hepatitis C virus (HCV) account for the majority of cases of hepatic cirrhosis and hepatocellular carcinoma (HCC) worldwide. Small, non-coding RNAs play important roles in virus-host interactions. We used high throughput sequencing to conduct an unbiased profiling of small (14-40 nts) RNAs in liver from Japanese subjects with advanced hepatitis B or C and hepatocellular carcinoma (HCC). Small RNAs derived from tRNAs, specifically 30-35 nucleotide-long 5' tRNA-halves (5' tRHs), were abundant in non-malignant liver and significantly increased in humans and chimpanzees with chronic viral hepatitis. 5' tRH abundance exceeded microRNA abundance in most infected non-cancerous tissues. In contrast, in matched cancer tissue, 5' tRH abundance was reduced, and relative abundance of individual 5' tRHs was altered. In hepatitis B-associated HCC, 5' tRH abundance correlated with expression of the tRNA-cleaving ribonuclease, angiogenin. These results demonstrate that tRHs are the most abundant small RNAs in chronically infected liver and that their abundance is altered in liver cancer.

Beta cell 5'-shifted isomiRs are candidate regulatory hubs in type 2 diabetes.

Next-generation deep sequencing of small RNAs has unveiled the complexity of the microRNA (miRNA) transcriptome, which is in large part due to the diversity of miRNA sequence variants ("isomiRs"). Changes to a miRNA's seed sequence (nucleotides 2-8), including shifted start positions, can redirect targeting to a dramatically different set of RNAs and alter biological function. We performed deep sequencing of small RNA from mouse insulinoma (MIN6) cells (widely used as a surrogate for the study of pancreatic beta cells) and developed a bioinformatic analysis pipeline to profile isomiR diversity. Additionally, we applied the pipeline to recently published small RNA-seq data from primary human beta cells and whole islets and compared the miRNA profiles with that of MIN6. We found that: (1) the miRNA expression profile in MIN6 cells is highly correlated with those of primary human beta cells and whole islets; (2) miRNA loci can generate multiple highly expressed isomiRs with different 5'-start positions (5'-isomiRs); (3) isomiRs with shifted start positions (5'-shifted isomiRs) are highly expressed, and can be as abundant as their unshifted counterparts (5'-reference miRNAs). Finally, we identified 10 beta cell miRNA families as candidate regulatory hubs in a type 2 diabetes (T2D) gene network. The most significant candidate hub was miR-29, which we demonstrated regulates the mRNA levels of several genes critical to beta cell function and implicated in T2D. Three of the candidate miRNA hubs were novel 5'-shifted isomiRs: miR-375+1, miR-375-1 and miR-183-5p+1. We showed by in silico target prediction and in vitro transfection studies that both miR-375+1 and miR-375-1 are likely to target an overlapping, but distinct suite of beta cell genes compared to canonical miR-375. In summary, this study characterizes the isomiR profile in beta cells for the first time, and also highlights the potential functional relevance of 5'-shifted isomiRs to T2D.