Phillipe Loher

The complex transcriptional landscape of the anucleate human platelet

BackgroundHuman blood platelets are essential to maintaining normal hemostasis, and platelet dysfunction often causes bleeding or thrombosis. Estimates of genome-wide platelet RNA expression using microarrays have provided insights to the platelet transcriptome but were limited by the number of known transcripts. The goal of this effort was to deep-sequence RNA from leukocyte-depleted platelets to capture the complex profile of all expressed transcripts.ResultsFrom each of four healthy individuals we generated long RNA (≥40 nucleotides) profiles from total and ribosomal-RNA depleted RNA preparations, as well as short RNA (<40 nucleotides) profiles. Analysis of ~1 billion reads revealed that coding and non-coding platelet transcripts span a very wide dynamic range (≥16 PCR cycles beyond β-actin), a result we validated through qRT-PCR on many dozens of platelet messenger RNAs. Surprisingly, ribosomal-RNA depletion significantly and adversely affected estimates of the relative abundance of transcripts. Of the known protein-coding loci, ~9,500 are present in human platelets. We observed a strong correlation between mRNAs identified by RNA-seq and microarray for well-expressed mRNAs, but RNASeq identified many more transcripts of lower abundance and permitted discovery of novel transcripts.ConclusionsOur analyses revealed diverse classes of non-coding RNAs, including: pervasive antisense transcripts to protein-coding loci; numerous, previously unreported and abundant microRNAs; retrotransposons; and thousands of novel un-annotated long and short intronic transcripts, an intriguing finding considering the anucleate nature of platelets. The data are available through a local mirror of the UCSC genome browser and can be accessed at:http://cm.jefferson.edu/platelets_2012/.

Sex hormone-dependent tRNA halves enhance cell proliferation in breast and prostate cancers

Significance Although transfer RNAs (tRNAs) are best known as adapter molecules essential for translation, recent biochemical and computational evidence has led to a previously unexpected conceptual consensus that tRNAs are not always end products but can further serve as a source of small functional RNAs. Here we report that a novel type of tRNA-derived small RNA, termed SHOT-RNAs, are specifically and abundantly expressed in sex hormone-dependent breast and prostate cancers. SHOT-RNAs are produced from aminoacylated mature tRNAs by angiogenin-mediated cleavage of the anticodon loop, which is promoted by sex hormones and their receptors. We identified the complete repertoire of SHOT-RNAs, and also found their functional significance in cell proliferation. These results have unveiled a novel tRNA-engaged pathway in tumorigenesis. Sex hormones and their receptors play critical roles in the development and progression of the breast and prostate cancers. Here we report that a novel type of transfer RNA (tRNA)-derived small RNA, termed Sex HOrmone-dependent TRNA-derived RNAs (SHOT-RNAs), are specifically and abundantly expressed in estrogen receptor (ER)-positive breast cancer and androgen receptor (AR)-positive prostate cancer cell lines. SHOT-RNAs are not abundantly present in ER− breast cancer, AR− prostate cancer, or other examined cancer cell lines from other tissues. ER-dependent accumulation of SHOT-RNAs is not limited to a cell culture system, but it also occurs in luminal-type breast cancer patient tissues. SHOT-RNAs are produced from aminoacylated mature tRNAs by angiogenin-mediated anticodon cleavage, which is promoted by sex hormones and their receptors. Resultant 5′- and 3′-SHOT-RNAs, corresponding to 5′- and 3′-tRNA halves, bear a cyclic phosphate (cP) and an amino acid at the 3′-end, respectively. By devising a “cP-RNA-seq” method that is able to exclusively amplify and sequence cP-containing RNAs, we identified the complete repertoire of 5′-SHOT-RNAs. Furthermore, 5′-SHOT-RNA, but not 3′-SHOT-RNA, has significant functional involvement in cell proliferation. These results have unveiled a novel tRNA-engaged pathway in tumorigenesis of hormone-dependent cancers and implicate SHOT-RNAs as potential candidates for biomarkers and therapeutic targets.

Interactive exploration of RNA22 microRNA target predictions

MicroRNA (miRNA) target prediction is an important problem. Given an miRNA sequence the task is to determine the identity of the messenger RNAs targeted by it, the locations within them where the interactions happen and the specifics of the formed heteroduplexes. Here, we describe a web-based application, RNA22-GUI, which we have designed and implemented for the interactive exploration and in-context visualization of predictions by RNA22, one of the popular miRNA target prediction algorithms. Central to our design has been the requirement to provide informative and comprehensive visualization that is integrated with interactive search capabilities and permits one to selectively isolate and focus on relevant information that is distilled on-the-fly from a large repository of pre-compiled predictions. RNA22-GUI is currently available for Homo sapiens, Mus musculus, Drosophila melanogaster and Caenorhabditis elegans.

The human platelet: strong transcriptome correlations among individuals associate weakly with the platelet proteome

BackgroundFor the anucleate platelet it has been unclear how well platelet transcriptomes correlate among different donors or across different RNA profiling platforms, and what the transcriptomes’ relationship is with the platelet proteome. We profiled the platelet transcriptome of 10 healthy young males (5 white and 5 black) with no notable clinical history using RNA sequencing and by Affymetrix microarray.ResultsWe found that the abundance of platelet mRNA transcripts was highly correlated across the 10 individuals, independently of race and of the employed technology. Our RNA-seq data showed that these high inter-individual correlations extend beyond mRNAs to several categories of non-coding RNAs. Pseudogenes represented a notable exception by exhibiting a difference in expression by race. Comparison of our mRNA signatures to a publicly available quantitative platelet proteome showed that most (87.5%) identified platelet proteins had a detectable corresponding mRNA. However, a high number of mRNAs that were present in the transcriptomes of all 10 individuals had no representation in the proteome. Spearman correlations of the relative abundances for those genes represented by both an mRNA and a protein showed a weak (~0.3) connection. Further analysis of the overlapping and non-overlapping platelet mRNAs and proteins identified gene groups corresponding to distinct cellular processes.ConclusionsThe results of our analyses provide novel insights for platelet biology, show only a weak connection between the platelet transcriptome and proteome, and indicate that it is feasible to assemble a platelet mRNA-ome that can serve as a reference for future platelet transcriptomic studies of human health and disease.Reviewed byThis article was reviewed by Dr Mikhail Dozmorov (nominated by Dr Yuri Gusev), Dr Neil Smalheiser and Dr Eugene Koonin.

Argonaute CLIP-Seq reveals miRNA targetome diversity across tissue types

To date, analyses of individual targets have provided evidence of a miRNA targetome that extends beyond the boundaries of messenger RNAs (mRNAs) and can involve non-Watson-Crick base pairing in the miRNA seed region. Here we report our findings from analyzing 34 Argonaute HITS-CLIP datasets from several human and mouse cell types. Investigation of the architectural (i.e. bulge vs. contiguous pairs) and sequence (Watson-Crick vs. G:U pairs) preferences for human and mouse miRNAs revealed that many heteroduplexes are “non-canonical” i.e. their seed region comprises G:U and bulge combinations. The genomic distribution of miRNA targets differed distinctly across cell types but remained congruent across biological replicates of the same cell type. For some cell types intergenic and intronic targets were more frequent whereas in other cell types mRNA targets prevailed. The findings suggest an expanded model of miRNA targeting that is more frequent than the standard model currently in use. Lastly, our analyses of data from different cell types and laboratories revealed consistent Ago-loaded miRNA profiles across replicates whereas, unexpectedly, the Ago-loaded targets exhibited a much more dynamic behavior across biological replicates.

Molecular dynamics simulations of Ago silencing complexes reveal a large repertoire of admissible ‘seed-less’ targets

To better understand the recognition mechanism of RISC and the repertoire of guide-target interactions we introduced G:U wobbles and mismatches at various positions of the microRNA (miRNA) ‘seed’ region and performed all-atom molecular dynamics simulations of the resulting Ago-miRNA:mRNA ternary complexes. Our simulations reveal that many modifications, including combinations of multiple G:U wobbles and mismatches in the seed region, are admissible and result in only minor structural fluctuations that do not affect overall complex stability. These results are further supported by analyses of HITS-CLIP data. Lastly, introduction of disruptive mutations revealed a bending motion of the PAZ domain along the L1/L2 ‘hinge’ and a subsequent opening of the nucleic-acid-binding channel. Our findings suggest that the spectrum of a miRNAs admissible targets is different from what is currently anticipated by the canonical seed-model. Moreover, they provide a likely explanation for the previously reported sequence-dependent regulation of unintended targeting by siRNAs.

Beyond the one-locus-one-miRNA paradigm: microRNA isoforms enable deeper insights into breast cancer heterogeneity

Here we describe our study of miRNA isoforms (isomiRs) in breast cancer (BRCA) and normal breast data sets from the Cancer Genome Atlas (TCGA) repository. We report that the full isomiR profiles, from both known and novel human-specific miRNA loci, are particularly rich in information and can distinguish tumor from normal tissue much better than the archetype miRNAs. IsomiR expression is also dependent on the patients race, exemplified by miR-183-5p, several isomiRs of which are upregulated in triple negative BRCA in white but not black women. Additionally, we find that an isomiRs 5′ endpoint and length, but not the genomic origin, are key determinants of the regulation of its expression. Overexpression of distinct miR-183-5p isomiRs in MDA-MB-231 cells followed by microarray analysis revealed that each isomiR has a distinct impact on the cellular transcriptome. Parallel integrative analysis of mRNA expression from BRCA data sets of the TCGA repository demonstrated that isomiRs can distinguish between the luminal A and luminal B subtypes and explain in more depth the molecular differences between them than the archetype molecules. In conclusion, our findings provide evidence that post-transcriptional studies of BRCA will benefit from transcending the one-locus-one-miRNA paradigm and taking into account all isoforms from each miRNA locus as well as the patients race.

Knowledge about the presence or absence of miRNA isoforms (isomiRs) can successfully discriminate amongst 32 TCGA cancer types.

Abstract Isoforms of human miRNAs (isomiRs) are constitutively expressed with tissue- and disease-subtype-dependencies. We studied 10 271 tumor datasets from The Cancer Genome Atlas (TCGA) to evaluate whether isomiRs can distinguish amongst 32 TCGA cancers. Unlike previous approaches, we built a classifier that relied solely on ‘binarized’ isomiR profiles: each isomiR is simply labeled as ‘present’ or ‘absent’. The resulting classifier successfully labeled tumor datasets with an average sensitivity of 90% and a false discovery rate (FDR) of 3%, surpassing the performance of expression-based classification. The classifier maintained its power even after a 15× reduction in the number of isomiRs that were used for training. Notably, the classifier could correctly predict the cancer type in non-TCGA datasets from diverse platforms. Our analysis revealed that the most discriminatory isomiRs happen to also be differentially expressed between normal tissue and cancer. Even so, we find that these highly discriminating isomiRs have not been attracting the most research attention in the literature. Given their ability to successfully classify datasets from 32 cancers, isomiRs and our resulting ‘Pan-cancer Atlas’ of isomiR expression could serve as a suitable framework to explore novel cancer biomarkers.

Nuclear and mitochondrial tRNA-lookalikes in the human genome

We are interested in identifying and characterizing loci of the human genome that harbor sequences resembling known mitochondrial and nuclear tRNAs. To this end, we used the known nuclear and mitochondrial tRNA genes (the “tRNA-Reference” set) to search for “tRNA-lookalikes” and found many such loci at different levels of sequence conservation. We find that the large majority of these tRNA-lookalikes resemble mitochondrial tRNAs and exhibit a skewed over-representation in favor of some mitochondrial anticodons. Our analysis shows that the tRNA-lookalikes have infiltrated specific chromosomes and are preferentially located in close proximity to known nuclear tRNAs (z-score ≤ −2.54, P-value ≤ 0.00394). Examination of the transcriptional potential of these tRNA-lookalike loci using public transcript annotations revealed that more than 20% of the lookalikes are transcribed as part of either known protein-coding pre-mRNAs, known lncRNAs, or known non-protein-coding RNAs, while public RNA-seq data perfectly agreed with the endpoints of tRNA-lookalikes. Interestingly, we found that tRNA-lookalikes are significantly depleted in known genetic variations associated with human health and disease whereas the known tRNAs are enriched in such variations. Lastly, a manual comparative analysis of the cloverleaf structure of several of the transcribed tRNA-lookalikes revealed no disruptive mutations suggesting the possibility that these loci give rise to functioning tRNA molecules.

MINTbase: a framework for the interactive exploration of mitochondrial and nuclear tRNA fragments

Motivation: It has been known that mature transfer RNAs (tRNAs) that are encoded in the nuclear genome give rise to short molecules, collectively known as tRNA fragments or tRFs. Recently, we reported that, in healthy individuals and in patients, tRFs are constitutive, arise from mitochondrial as well as from nuclear tRNAs, and have composition and abundances that depend on a person’s sex, population origin and race as well as on tissue, disease and disease subtype. Our findings as well as similar work by other groups highlight the importance of tRFs and presage an increase in the community’s interest in elucidating the roles of tRFs in health and disease. Results: We created MINTbase, a web-based framework that serves the dual-purpose of being a content repository for tRFs and a tool for the interactive exploration of these newly discovered molecules. A key feature of MINTbase is that it deterministically and exhaustively enumerates all possible genomic locations where a sequence fragment can be found and indicates which fragments are exclusive to tRNA space, and thus can be considered as tRFs: this is a very important consideration given that the genomes of higher organisms are riddled with partial tRNA sequences and with tRNA-lookalikes whose aberrant transcripts can be mistaken for tRFs. MINTbase is extremely flexible and integrates and presents tRF information from multiple yet interconnected vantage points (‘vistas’). Vistas permit the user to interactively personalize the information that is returned and the manner in which it is displayed. MINTbase can report comparative information on how a tRF is distributed across all anticodon/amino acid combinations, provides alignments between a tRNA and multiple tRFs with which the user can interact, provides details on published studies that reported a tRF as expressed, etc. Importantly, we designed MINTbase to contain all possible tRFs that could ever be produced by mature tRNAs: this allows us to report on their genomic distributions, anticodon/amino acid properties, alignments, etc. while giving users the ability to at-will investigate candidate tRF molecules before embarking on focused experimental explorations. Lastly, we also introduce a new labeling scheme that is tRF-sequence-based and allows users to associate a tRF with a universally unique label (‘tRF-license plate’) that is independent of a genome assembly and does not require any brokering mechanism. Availability and Implementation: MINTbase is freely accessible at http://cm.jefferson.edu/MINTbase/. Dataset submissions to MINTbase can be initiated at http://cm.jefferson.edu/MINTsubmit/. Contact: isidore.rigoutsos@jefferson.edu Supplementary information: Supplementary data are available at Bioinformatics online.