Search

RNA-sequencing (RNA-Seq) is a powerful transcriptome profiling technology enabling transcript discovery and quantification. RNA-Seq data are large, and most commonly used as a source of genelevel quantification measurements, whilst the underlying assemblies of reads, if inspected, are usually viewed as sequence reads mapped on to a reference genome. Whilst sufficient for many needs, when the underlying transcript assemblies are complex, this visualisation approach can be limiting; errors in assembly can be difficult to spot and interpretation of splicing events is challenging.Here we report on the development of a graph-based visualisation method as a complementary approach to understanding transcript diversity and read assembly from short-read RNA-Seq data. Following the mapping of reads to the reference genome, read-to-read comparison is performed on all reads mapping to a given gene, producing a matrix of weighted similarity scores between reads. This is used to produce an RNA assembly graph where nodes represent reads derived from a cDNA and edges similarity scores between reads, above a defined threshold. Visualisation of resulting graphs is performed using Graphia Professional. This tool can render the often large and complex graph topologies that result from DNA/RNA sequence assembly in 3D space and supports info rmatio no verlay on to nodes, e.g. transcript models. We have also implemented an analysis pipeline for the creation of RNA assembly graphs with both a command-line and web-based interface that allows users to create and visualise these data. Here we demonstrate the utility of this approach on RNA-Seq data, including the unusual structure of these graphs and how they can be used to identify issues in assembly, repetitive sequences within transcripts and splice variants. We believe this approach has the potential to significantly improve our understanding of transcript complexity.

Proteins need to be tightly regulated as they control biological processes in most normal cellular functions. The precise mechanisms of regulation are rarely completely understood but can involve binding of endogenous ligands and/or partner proteins at specific locations on a protein that can modulate function. Often, these additional secondary binding sites appear separate to the primary binding site, which, for example for an enzyme, may bind a substrate. In previous work, we have uncovered several examples in which secondary binding sites were discovered on proteins using fragment screening approaches. In each case, we were able to establish that the newly identified secondary binding site was biologically relevant as it was able to modulate function by the binding of a small molecule. In this study, we investigate how often secondary binding sites are located on proteins by analyzing 24 protein targets for which we have performed a fragment screen using X-ray crystallography. Our analysis shows that, surprisingly, the majority of proteins contain secondary binding sites based on their ability to bind fragments. Furthermore, sequence analysis of these previously unknown sites indicate high conservation, which suggests that they may have a biological function, perhaps via an allosteric mechanism. Comparing the physicochemical properties of the secondary sites with known primary ligand binding sites also shows broad similarities indicating that many of the secondary sites may be druggable in nature with small molecules that could provide new opportunities to modulate potential therapeutic targets.

Despite their clinicopathologic heterogeneity, malignant germ cell tumors (GCT) share molecular abnormalities that are likely to be functionally important. In this study, we investigated the potential significance of downregulation of the let-7 family of tumor suppressor microRNAs in malignant GCTs. Microarray results from pediatric and adult samples (n = 45) showed that LIN28, the negative regulator of let-7 biogenesis, was abundant in malignant GCTs, regardless of patient age, tumor site, or histologic subtype. Indeed, a strong negative correlation existed between LIN28 and let-7 levels in specimens with matched datasets. Low let-7 levels were biologically significant, as the sequence complementary to the 2 to 7 nt common let-7 seed “GAGGUA” was enriched in the 3′ untranslated regions of mRNAs upregulated in pediatric and adult malignant GCTs, compared with normal gonads (a mixture of germ cells and somatic cells). We identified 27 mRNA targets of let-7 that were upregulated in malignant GCT cells, confirming significant negative correlations with let-7 levels. Among 16 mRNAs examined in a largely independent set of specimens by quantitative reverse transcription PCR, we defined negative-associations with let-7e levels for six oncogenes, including MYCN, AURKB, CCNF, RRM2, MKI67, and C12orf5 (when including normal control tissues). Importantly, LIN28 depletion in malignant GCT cells restored let-7 levels and repressed all of these oncogenic let-7 mRNA targets, with LIN28 levels correlating with cell proliferation and MYCN levels. Conversely, ectopic expression of let-7e was sufficient to reduce proliferation and downregulate MYCN, AURKB, and LIN28, the latter via a double-negative feedback loop. We conclude that the LIN28/let-7 pathway has a critical pathobiologic role in malignant GCTs and therefore offers a promising target for therapeutic intervention. Cancer Res; 73(15); 4872–84. ©2013 AACR.

Mast cells have essential effector and immunoregulatory functions in IgE-associated allergic disorders and certain innate and adaptive immune responses, but the role of miRNAs in regulating mast cell functions is almost completely unexplored. To examine the role of the activation-induced miRNA miR-221 in mouse mast cells, we developed robust lentiviral systems for miRNA overexpression and depletion. While miR-221 favored mast cell adhesion and migration towards SCF or antigen in trans-well migration assays, as well as cytokine production and degranulation in response to IgE-antigen complexes, neither miR-221 overexpression, nor its ablation, interfered with mast cell differentiation. Transcriptional profiling of miR-221-overexpressing mast cells revealed modulation of many transcripts, including several associated with the cytoskeleton; indeed, miR-221 overexpression was associated with reproducible increases in cortical actin in mast cells, and with altered cellular shape and cell cycle in murine fibroblasts. Our bioinformatics analysis showed that this effect was likely mediated by the composite effect of miR-221 on many primary and secondary targets in resting cells. Indeed, miR- 221-induced cellular alterations could not be recapitulated by knockdown of one of the major targets of miR-221. We propose a model in which miR-221 has two different roles in mast cells: in resting cells, basal levels of miR-221 contribute to the regulation of the cell cycle and cytoskeleton, a general mechanism probably common to other miR-221- expressing cell types, such as fibroblasts. Vice versa, upon induction in response to mast cell stimulation, miR-221 effects are mast cell-specific and activation-dependent, contributing to the regulation of degranulation, cytokine production and cell adherence. Our studies provide new insights into the roles of miR-221 in mast cell biology, and identify novel mechanisms that may contribute to mast cell-related pathological conditions, such as asthma, allergy and mastocytosis.

Background We hypothesised that differences in microRNA expression profiles contribute to the contrasting natural history and clinical outcome of the two most common types of malignant germ cell tumour (GCT), yolk sac tumours (YSTs) and germinomas. Results By direct comparison, using microarray data for paediatric GCT samples and published qRT-PCR data for adult samples, we identified microRNAs significantly up-regulated in YSTs (n = 29 paediatric, 26 adult, 11 overlapping) or germinomas (n = 37 paediatric). By Taqman qRT-PCR we confirmed differential expression of 15 of 16 selected microRNAs and further validated six of these (miR-302b, miR-375, miR-200b, miR-200c, miR-122, miR-205) in an independent sample set. Interestingly, the miR-302 cluster, which is over-expressed in all malignant GCTs, showed further over-expression in YSTs versus germinomas, representing six of the top eight microRNAs over-expressed in paediatric YSTs and seven of the top 11 in adult YSTs. To explain this observation, we used mRNA expression profiles of paediatric and adult malignant GCTs to identify 10 transcription factors (TFs) consistently over-expressed in YSTs versus germinomas, followed by linear regression to confirm associations between TF and miR-302 cluster expression levels. Using the sequence motif analysis environment iMotifs, we identified predicted binding sites for four of the 10 TFs (GATA6, GATA3, TCF7L2 and MAF) in the miR-302 cluster promoter region. Finally, we showed that miR-302 family over-expression in YST is likely to be functionally significant, as mRNAs down-regulated in YSTs were enriched for 3' untranslated region sequences complementary to the common seed of miR-302a~miR-302d. Such mRNAs included mediators of key cancer-associated processes, including tumour suppressor genes, apoptosis regulators and TFs. Conclusions Differential microRNA expression is likely to contribute to the relatively aggressive behaviour of YSTs and may enable future improvements in clinical diagnosis and/or treatment.

MicroRNAs (miRNAs) are important genetic regulators of development, differentiation, growth, and metabolism. The mammalian genome encodes ≈500 known miRNA genes. Approximately 50% are expressed from non-protein-coding transcripts, whereas the rest are located mostly in the introns of coding genes. Intronic miRNAs are generally transcribed coincidentally with their host genes. However, the nature of the primary transcript of intergenic miRNAs is largely unknown. We have performed a large-scale analysis of transcription start sites, polyadenylation signals, CpG islands, EST data, transcription factor-binding sites, and expression ditag data surrounding intergenic miRNAs in the human genome to improve our understanding of the structure of their primary transcripts. We show that a significant fraction of primary transcripts of intergenic miRNAs are 3–4 kb in length, with clearly defined 5′ and 3′ boundaries. We provide strong evidence for the complete transcript structure of a small number of human miRNAs.