Differential repression of alternative transcripts: a screen for miRNA targets

Generation of Transcript Length Variants and Reprogramming of mRNA Splicing During Atherosclerosis Progression in ApoE-Deficient Mice

Background. Variant 3'UTRs provide mRNAs with different binding sites for miRNAs or RNA-binding proteins (RBPs) allowing the establishment of new regulatory environments. Regulation of 3'UTR length impacts on the control of gene expression by regulating accessibility of miRNAs or RBPs to homologous sequences in mRNAs. Objective. Studying the dynamics of mRNA length variations in atherosclerosis (ATS) progression and reversion in ApoE-deficient mice exposed to a high-fat diet and treated with an αCD40-specific siRNA or with a sequence-scrambled siRNA as control. Methods. We gathered microarray mRNA expression data from the aortas of mice after 2 or 16 weeks of treatments, and used these data in a Bioinformatics analysis. Results. Here, we report the lengthening of the 5'UTR/3'UTRs and the shortening of the CDS in downregulated mRNAs during ATS progression. Furthermore, treatment with the αCD40-specific siRNA resulted in the partial reversion of the 3'UTR lengthening. Exon analysis showed that these length variations were actually due to changes in the number of exons embedded in mRNAs, and the further examination of transcripts co-expressed at weeks 2 and 16 in mice treated with the control siRNA revealed a process of mRNA isoform switching in which transcript variants differed in the patterns of alternative splicing or activated latent/cryptic splice sites. Conclusion. We document length variations in the 5'UTR/3'UTR and CDS of mRNAs downregulated during atherosclerosis progression and suggest a role for mRNA splicing reprogramming and transcript isoform switching in the generation of disease-related mRNA sequence diversity and variability.

Genomic analyses of intricate interaction of TE-lncRNA overlapping genes with miRNAs in human diseases

BACKGROUND

Transposable elements (TEs) are known to be inserted into genome to create transcript isoforms or to generate long non-coding RNA (lncRNA) sequences. The insertion of TEs generates a gene protein sequence within the genome, but also provides a microRNA (miRNA) regulatory region.

OBJECTIVE

To determine the effect of gene sequence changes caused by TE insertion on miRNA binding and to investigate the formation of an overlapping lncRNA that represses it.

METHODS

The distribution of overlapping regions between exons and TE regions with lncRNA was examined using the Bedtools. miRNAs that can bind to those overlapping regions were identified through the miRDB web program. For TE-lncRNA overlapping genes, bioinformatic analysis was conducted using DAVID web database. Differential expression analysis was conducted using data from the GEO dataset and TCGA.

RESULTS

Most TEs were distributed more frequently in untranslated regions than open reading frames. There were 30 annotated TE-lncRNA overlapping genes with same strand that could bind to the same miRNA. As a result of identifying the association between these 30 genes and diseases, TGFB2, FCGR2A, DCTN5, and IFI6 were associated with breast cancer, and HMGCS1, FRMD4A, EDNRB, and SNCA were associated with Alzheimer's disease. Analysis of the GEO and TCGA data showed that the relevant expression of miR-891a and miR-28, which bind to the TE overlapping region of DCTN5 and HMGCS1, decreased.

CONCLUSION

This study indicates that the interaction between TE-lncRNA overlapping genes and miRNAs can affect disease progression.

Determinants of Functional MicroRNA Targeting

MicroRNAs (miRNAs) play cardinal roles in regulating biological pathways and processes, resulting in significant physiological effects. To understand the complex regulatory network of miRNAs, previous studies have utilized massivescale datasets of miRNA targeting and attempted to computationally predict the functional targets of miRNAs. Many miRNA target prediction tools have been developed and are widely used by scientists from various fields of biology and medicine. Most of these tools consider seed pairing between miRNAs and their mRNA targets and additionally consider other determinants to improve prediction accuracy. However, these tools exhibit limited prediction accuracy and high false positive rates. The utilization of additional determinants, such as RNA modifications and RNA-binding protein binding sites, may further improve miRNA target prediction. In this review, we discuss the determinants of functional miRNA targeting that are currently used in miRNA target prediction and the potentially predictive but unappreciated determinants that may improve prediction accuracy.

Whole transcriptome sequencing reveals core genes related to spermatogenesis in Bactrian camels

Bactrian camels survive and reproduce better in extreme climatic conditions than other domestic animals can. However, the reproductive efficiency of camels under their natural pastoral conditions is low. Several factors affect mammalian reproductive performance, including testicular development, semen quality, libido, and mating ability. Testis is a main reproductive organ of the male and is responsible for producing spermatozoa and hormones. However, our understanding of the expression patterns of the genes in camel testis is minimal. Thus, we performed total RNA-sequencing to investigate the gene expression pattern. As a result, 1,538 differential expressed mRNAs (DEmRNAs), 702 differential expressed long non-coding RNAs (DElncRNAs), and 61 differential expressed microRNAs (DEmiRNAs) were identified between pubertal and adult Bactrian camel testes. Then the genomic features, length distribution, and other characteristics of the lncRNAs and mRNAs in the Bactrian camel testis were investigated. Target genes of the DEmiRNAs and DEmRNAs were further subjected to gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses. Genes, such as AMHR2, FGF1, ACTL7A, GATA4, WNT4, ID2, LAMA1, IGF1, INHBB, and TLR2, were mainly involved in the TGF-β, PI3K-AKT, Wnt, GnRH, and Hippo signaling pathways which relate to spermatogenesis. Some of the DEmiRNAs were predicted to be associated with numerous DElncRNAs and DEmRNAs through competing endogenous RNA (ceRNA) regulatory network. At last, the candidate genes were validated by RT-qPCR, dual fluorescent reporter gene, and a fluorescence in situ hybridization (FISH) assay. This research provides high-throughput RNA sequencing data of the testes of Bactrian camels across different developmental stages. It lays the foundation for further investigations on lncRNAs, miRNAs, and mRNAs that involved in Bactrian camel spermatogenesis.

Dynamic Variations of 3'UTR Length Reprogram the mRNA Regulatory Landscape

This paper concerns 3′-untranslated regions (3′UTRs) of mRNAs, which are non-coding regulatory platforms that control stability, fate and the correct spatiotemporal translation of mRNAs. Many mRNAs have polymorphic 3′UTR regions. Controlling 3′UTR length and sequence facilitates the regulation of the accessibility of functional effectors (RNA binding proteins, miRNAs or other ncRNAs) to 3′UTR functional boxes and motifs and the establishment of different regulatory landscapes for mRNA function. In this context, shortening of 3′UTRs would loosen miRNA or protein-based mechanisms of mRNA degradation, while 3′UTR lengthening would strengthen accessibility to these effectors. Alterations in the mechanisms regulating 3′UTR length would result in widespread deregulation of gene expression that could eventually lead to diseases likely linked to the loss (or acquisition) of specific miRNA binding sites. Here, we will review the mechanisms that control 3′UTR length dynamics and their alterations in human disorders. We will discuss, from a mechanistic point of view centered on the molecular machineries involved, the generation of 3′UTR variability by the use of alternative polyadenylation and cleavage sites, of mutually exclusive terminal alternative exons (exon skipping) as well as by the process of exonization of Alu cassettes to generate new 3′UTRs with differential functional features.

Spatiotemporal Regulation of Transcript Isoform Expression in the Hippocampus

Proper development and plasticity of hippocampal neurons require specific RNA isoforms to be expressed in the right place at the right time. Precise spatiotemporal transcript regulation requires the incorporation of essential regulatory RNA sequences into expressed isoforms. In this review, we describe several RNA processing strategies utilized by hippocampal neurons to regulate the spatiotemporal expression of genes critical to development and plasticity. The works described here demonstrate how the hippocampus is an ideal investigative model for uncovering alternate isoform-specific mechanisms that restrict the expression of transcripts in space and time.

MicroRNAs are indispensable for the proliferation and differentiation of adult neural progenitor cells in mice

More than two decades after the discovery of adult neurogenesis in humans, researchers still struggle to elucidate the underlying transcriptional and post-transcriptional mechanisms. RNA interference is a crucially important process in the central nervous system, and its role in adult neurogenesis is poorly understood. In this work, we address the role of Dicer-dependent microRNA biogenesis in neuronal differentiation of adult neural stem cells within the subventricular zone of the mouse brain. Loss of the Dicer1 gene in the tailless (Tlx)-positive cells did not cause the decline in their numbers, but severely affected differentiation. Thus, our findings identify yet another phenomenon associated with microRNA pathway deregulation in adult neural stem cells which might be of relevance both for neuroscience and clinical practice.

The Oncogenic Functions of Insulin-like Growth Factor 2 mRNA-Binding Protein 3 in Human Carcinomas

IGF2BP3 (also known as IMP3, KOC), a member of the insulin-like growth factor mRNA-binding protein family (IMPs), has been a research target in recent studies of promoting embryo development and exacerbating cancer. IGF2BP3 is ubiquitously expressed in early embryogenesis stages but limited in postembryonic stages, which is important in many physiological aspects such as stem cell renewal, morphological development and metabolism. A large number of studies show that IGF2BP3 interacts with many kinds of non-coding RNAs and proteins to promote cancer cell proliferation and metastasis and inhibit cancer cell apoptosis. As IGF2BP3 is highly expressed in advanced cancers and associated with poor overall survival rates of patients, it may be a potential molecular marker in cancer diagnosis for the detection of cancerous tissues and an indicator of cancer stages. Therefore, anti-IGF2BP3 drugs or monoclonal antibodies are expected as new therapeutic methods in cancer treatment. This review summarizes recent findings among IGF2BP3, RNA and proteins in cancer processes, with a focus on its cancer-promoting mechanisms and potential application as a new biomarker for cancer diagnosis and treatment.

The effects of alternative splicing on miRNA binding sites in bladder cancer

Eukaryotic organisms have developed a variety of mechanisms to regulate translation post-transcriptionally, including but not limited to the use of miRNA silencing in many species. One method of post-transcriptional regulation is through miRNAs that bind to the 3′ UTRs to regulate mRNA abundance and influence protein expression. Therefore, the diversity of mRNA 3′ UTRs mediating miRNA binding sites influence miRNA-mediated regulation. Alternative polyadenylation, by shortening mRNA isoforms, increases the diversity of 3′ UTRs; moreover, short mRNA isoforms elude miRNA-medicated repression. Because no current prediction methods for putative miRNA target sites consider whether or not 1) splicing-informed miRNA binding sites and/or 2) the use of 3′ UTRs provide higher resolution or functionality, we sought to identify not only the genome-wide impact of using exons in mRNA 3′ UTRs but also their functional connection to miRNA regulation and clinical outcomes in cancer. With a genome-wide expression of mRNA and miRNA quantified by 395 bladder cancer cases from The Cancer Genome Atlas (TCGA), we 1) demonstrate the diversity of 3′ UTRs affecting miRNA efficiency and 2) identify a set of genes clinically associated with mRNA expression in bladder cancer. Knowledge of 3′ UTR diversity will not only be a useful addition to current miRNA target prediction algorithms but also enhance the clinical utility of mRNA isoforms in the expression of mRNA in cancer. Thus, variability among cancer patient’s variability in molecular signatures based on these exon usage events in 3′ UTR along with miRNAs in bladder cancer may lead to better prognostic/treatment strategies for improved precision medicine.

A Fresh Look at Huntingtin mRNA Processing in Huntington's Disease

Huntington's disease (HD) is an inherited neurodegenerative disorder caused by a mutation that expands the polyglutamine (CAG) repeat in exon 1 of the huntingtin (HTT) gene. Wild-type HTT protein interacts with other proteins to protect cells against toxic stimuli, mediate vesicle transport and endocytosis, and modulate synaptic activity. Mutant HTT protein disrupts autophagy, vesicle transport, neurotransmitter signaling, and mitochondrial function. Although many of the activities of wild-type HTT protein and the toxicities of mutant HTT protein are characterized, less is known about the activities of HTT mRNA. Most putative HD therapies aim to target mutant HTT mRNA before it is translated into the protein. Therefore, it is imperative to learn as much as we can about how cells handle both wild-type and mutant HTT mRNA so that effective therapies can be designed. Here, we review the structure of wild-type and mutant HTT mRNA, with emphasis on their alternatively polyadenylated or spliced isoforms. We then consider the abundance of HTT mRNA isoforms in HD and discuss the potential implications of these findings. Evidence in the review should be used to guide future research aimed at developing mRNA-lowering therapies for HD.

Cleavage and polyadenylation: Ending the message expands gene regulation

Cleavage and polyadenylation (pA) is a fundamental step that is required for the maturation of primary protein encoding transcripts into functional mRNAs that can be exported from the nucleus and translated in the cytoplasm. 3'end processing is dependent on the assembly of a multiprotein processing complex on the pA signals that reside in the pre-mRNAs. Most eukaryotic genes have multiple pA signals, resulting in alternative cleavage and polyadenylation (APA), a widespread phenomenon that is important to establish cell state and cell type specific transcriptomes. Here, we review how pA sites are recognized and comprehensively summarize how APA is regulated and creates mRNA isoform profiles that are characteristic for cell types, tissues, cellular states and disease.

Distinct cellular toxicity of two mutant huntingtin mRNA variants due to translation regulation

Huntington's disease (HD) is a neurodegenerative disorder caused by CAG repeat expansion within exon1 of the HTT gene. The gene generates two mRNA variants that carry either a short or long 3' untranslated region (3'UTR) while encoding the same protein. It remains unknown whether the two mRNA variants play distinct roles in HD pathogenesis. We found that the long HTT 3'UTR was capable of guiding mRNA to neuronal dendrites, suggesting that some long-form HTT mRNA is transported to dendrites for local protein synthesis. To assay roles of two HTT mRNA variants in cell bodies, we expressed mRNA harboring HTT exon1 containing 23x or 145x CAGs with the short or long 3'UTR. We found that mutant mRNA containing the short 3'UTR produced more protein aggregates and caused more apoptosis in both cultured neurons and HEK293 cells, compared with mutant mRNA containing the long 3'UTR. Although the two 3'UTRs did not affect mRNA stability, we detected higher levels of protein synthesis from mRNA containing the short 3'UTR than from mRNA containing the long 3'UTR. These results indicate that the long HTT 3'UTR suppresses translation. Thus, short-form mutant HTT mRNA will be more efficient in producing toxic protein than its long-form counterpart.

Computational and Experimental Identification of Tissue-Specific MicroRNA Targets

In this chapter we discuss computational methods for the prediction of microRNA (miRNA) targets. More specifically, we consider machine learning-based approaches and explain why these methods have been relatively unsuccessful in reducing the number of false positive predictions. Further we suggest approaches designed to improve their performance by considering tissue-specific target regulation. We argue that the miRNA targetome differs depending on the tissue type and introduce a novel algorithm that predicts miRNA targets specifically for colorectal cancer. We discuss features of miRNAs and target sites that affect target recognition, and how next-generation sequencing data can support the identification of novel miRNAs, differentially expressed miRNAs and their tissue-specific mRNA targets. In addition, we introduce some experimental approaches for the validation of miRNA targets as well as web-based resources sharing predicted and validated miRNA target interactions.

RNA binding proteins: a common denominator of neuronal function and dysfunction

In eukaryotic cells, gene activity is not directly reflected by protein levels because mRNA processing, transport, stability, and translation are co- and post-transcriptionally regulated. These processes, collectively known as the ribonome, are tightly controlled and carried out by a plethora of trans-acting RNA-binding proteins (RBPs) that bind to specific cis elements throughout the RNA sequence. Within the nervous system, the role of RBPs in brain function turns out to be essential due to the architectural complexity of neurons exemplified by a relatively small somal size and an extensive network of projections and connections. Thus far, RBPs have been shown to be indispensable for several aspects of neurogenesis, neurite outgrowth, synapse formation, and plasticity. Consequently, perturbation of their function is central in the etiology of an ever-growing spectrum of neurological diseases, including fragile X syndrome and the neurodegenerative disorders frontotemporal lobar degeneration and amyotrophic lateral sclerosis.

Regulation of genes in the arachidonic acid metabolic pathway by RNA processing and RNA-mediated mechanisms

Arachidonic acid (AA) is converted by enzymes in an important metabolic pathway to produce molecules known collectively as eicosanoids, 20 carbon molecules with significant physiological and pathological functions in the human body. Cyclooxygenase (COX) enzymes work in one arm of the pathway to produce prostaglandins (PGs) and thromboxanes (TXs), while the actions of 5-lipoxygenase (ALOX5 or 5LO) and its associated protein (ALOX5AP or FLAP) work in the other arm of the metabolic pathway to produce leukotrienes (LTs). The expression of the COX and ALOX5 enzymes that convert AA to eicosanoids is highly regulated at the post- or co-transcriptional level by alternative mRNA splicing, alternative mRNA polyadenylation, mRNA stability, and microRNA (miRNA) regulation. This review article will highlight these mechanisms of mRNA modulation.

Mir-20a regulates in vitro mineralization and BMP signaling pathway by targeting BMP-2 transcript in fish

MicroRNAs (miRNAs) are important regulators of vertebrate development but their role during skeletogenesis remains unknown. In this regard, we investigated the mineralogenic activity of miR-20a, a miRNA associated with osteogenesis, in fish bone-derived cells. Expression of miR-20a was up-regulated during differentiation and its overexpression inhibited mineralization, suggesting a role in fish tissue calcification. In this regard, a conserved miR-20a binding site was identified in bone morphogenetic protein 2 (BMP-2) 3'UTR and its functionality was evidenced through luciferase assays, and further confirmed by western-blot and qPCR. Type II BMP receptor (BMPR2) is also targeted by miR-20a in mammalian systems and evidence was collected for the presence of a binding site in fish sequences. We propose that miR-20a is a regulator of BMP pathway through specific action on BMP-2 and possibly BMPR2. Overexpression of miR-20a was also shown to up-regulate matrix Gla protein (MGP) transcript, a physiological inhibitor of calcification previously found to form a complex with BMP-2. We propose that MGP may play a role in the anti-mineralogenic effect promoted by miR-20a by decreasing availability of BMP-2. This study gives new insights into miRNA-mediated regulation of BMP-2, and sheds light into the potential role of miR-20a as a regulator of skeletogenesis.

MicroRNA: a prognostic biomarker and a possible druggable target for circumventing multidrug resistance in cancer chemotherapy

Multidrug resistance (MDR) is a major obstacle to successful cancer treatment. It is often associated with an increased efflux of a variety of structurally unrelated anticancer drugs by ATP-binding cassette (ABC) transporters including P-gp, ABCG2 and MRP1. MicroRNAs (miRNAs) are small non-coding RNAs that govern posttranscriptional regulation of target genes by interacting with specific sequences in their 3′ untranslated region (3′UTR), thereby promoting mRNA degradation or suppressing translation. Accumulating evidence suggests that alterations in miRNAs contribute to resistance to anticancer drugs. While miRNAs are well-known to be dysregulated in cancer, recent literature revealed that miRNA levels in biological samples may be correlated with chemotherapy response. This review summarized the coordinated network by which miRNA regulated MDR transporters. The usefulness of miRNAs as prognostic biomarkers for predicting chemotherapeutic outcome is discussed. MiRNAs may also represent druggable targets for circumvention of MDR.

A dynamic interplay between alternative polyadenylation and microRNA regulation: implications for cancer (Review)

Alternative polyadenylation and microRNA regulation are both mechanisms of post-transcriptional regulation of gene expression. Alternative polyadenylation often results in mRNA isoforms with the same coding sequence but different lengths of 3' UTRs, while microRNAs regulate gene expression by binding to specific mRNA 3' UTRs. In this sense, different isoforms of an mRNA may be differentially regulated by microRNAs, sometimes resulting in cellular proliferation and this mechanism is being speculated on as a potential cause for cancer development.

Fine-tuning of microRNA-mediated repression of mRNA by splicing-regulated and highly repressive microRNA recognition element

BACKGROUND

MicroRNAs are very small non-coding RNAs that interact with microRNA recognition elements (MREs) on their target messenger RNAs. Varying the concentration of a given microRNA may influence the expression of many target proteins. Yet, the expression of a specific target protein can be fine-tuned by alternative cleavage and polyadenylation to the corresponding mRNA.

RESULTS

This study showed that alternative splicing of mRNA is a fine-tuning mechanism in the cellular regulatory network. The splicing-regulated MREs are often highly repressive MREs. This phenomenon was observed not only in the hsa-miR-148a-regulated DNMT3B gene, but also in many target genes regulated by hsa-miR-124, hsa-miR-1, and hsa-miR-181a. When a gene contains multiple MREs in transcripts, such as the VEGF gene, the splicing-regulated MREs are again the highly repressive MREs. Approximately one-third of the analysable human MREs in MiRTarBase and TarBase can potentially perform the splicing-regulated fine-tuning. Interestingly, the high (+30%) repression ratios observed in most of these splicing-regulated MREs indicate associations with functions. For example, the MRE-free transcripts of many oncogenes, such as N-RAS and others may escape microRNA-mediated suppression in cancer tissues.

CONCLUSIONS

This fine-tuning mechanism revealed associations with highly repressive MRE. Since high-repression MREs are involved in many important biological phenomena, the described association implies that splicing-regulated MREs are functional. A possible application of this observed association is in distinguishing functionally relevant MREs from predicted MREs.

DNA damage induces targeted, genome-wide variation of poly(A) sites in budding yeast

Systemic response to DNA damage and other stresses is a complex process that includes changes in the regulation and activity of nearly all stages of gene expression. One gene regulatory mechanism used by eukaryotes is selection among alternative transcript isoforms that differ in polyadenylation [poly(A)] sites, resulting in changes either to the coding sequence or to portions of the 3' UTR that govern translation, stability, and localization. To determine the extent to which this means of regulation is used in response to DNA damage, we conducted a global analysis of poly(A) site usage in Saccharomyces cerevisiae after exposure to the UV mimetic, 4-nitroquinoline 1-oxide (4NQO). Two thousand thirty-one genes were found to have significant variation in poly(A) site distributions following 4NQO treatment, with a strong bias toward loss of short transcripts, including many with poly(A) sites located within the protein coding sequence (CDS). We further explored one possible mechanism that could contribute to the widespread differences in mRNA isoforms. The change in poly(A) site profile was associated with an inhibition of cleavage and polyadenylation in cell extract and a decrease in the levels of several key subunits in the mRNA 3'-end processing complex. Sequence analysis identified differences in the cis-acting elements that flank putatively suppressed and enhanced poly(A) sites, suggesting a mechanism that could discriminate between variable and constitutive poly(A) sites. Our analysis indicates that variation in mRNA length is an important part of the regulatory response to DNA damage.

Stability, delivery and functions of human sperm RNAs at fertilization

Increasing attention has focused on the significance of RNA in sperm, in light of its contribution to the birth and long-term health of a child, role in sperm function and diagnostic potential. As the composition of sperm RNA is in flux, assigning specific roles to individual RNAs presents a significant challenge. For the first time RNA-seq was used to characterize the population of coding and non-coding transcripts in human sperm. Examining RNA representation as a function of multiple methods of library preparation revealed unique features indicative of very specific and stage-dependent maturation and regulation of sperm RNA, illuminating their various transitional roles. Correlation of sperm transcript abundance with epigenetic marks suggested roles for these elements in the pre- and post-fertilization genome. Several classes of non-coding RNAs including lncRNAs, CARs, pri-miRNAs, novel elements and mRNAs have been identified which, based on factors including relative abundance, integrity in sperm, available knockout data of embryonic effect and presence or absence in the unfertilized human oocyte, are likely to be essential male factors critical to early post-fertilization development. The diverse and unique attributes of sperm transcripts that were revealed provides the first detailed analysis of the biology and anticipated clinical significance of spermatozoal RNAs.

Differential microRNA regulation correlates with alternative polyadenylation pattern between breast cancer and normal cells

Alternative polyadenylation (APA) could result in mRNA isoforms with variable lengths of 3′ UTRs. Gain of microRNA target sites in the 3′ UTR of a long mRNA isoform may cause different regulation from the corresponding short isoform. It has been known that cancer cells globally exhibit a lower ratio of long and short isoforms (LSR); that is, they tend to express larger amounts of short isoforms. The objective of this study is to illustrate the relationship between microRNA differential regulation and LSR. We retrieved public APA annotations and isoform expression profiles of breast cancer and normal cells from a high-throughput sequencing method study specific for the mRNA 3′ end. Combining microRNA expression profiles, we performed statistical analysis to reveal and estimate microRNA regulation on APA patterns in a global scale. First, we found that the amount of microRNA target sites in the alternative UTR (aUTR), the region only present in long isoforms, could affect the LSR of the target genes. Second, we observed that the genes whose aUTRs were targeted by up-regulated microRNAs in cancer cells had an overall lower LSR. Furthermore, the target sites of up-regulated microRNAs tended to appear in aUTRs. Finally, we demonstrated that the amount of target sites for up-regulated microRNAs in aUTRs correlated with the LSR change between cancer and normal cells. The results indicate that up-regulation of microRNAs might cause lower LSRs of target genes in cancer cells through degradation of their long isoforms. Our findings provide evidence of how microRNAs might play a crucial role in APA pattern shifts from normal to cancerous or proliferative states.

Alterations in polyadenylation and its implications for endocrine disease

INTRODUCTION

Polyadenylation is the process in which the pre-mRNA is cleaved at the poly(A) site and a poly(A) tail is added - a process necessary for normal mRNA formation. Genes with multiple poly(A) sites can undergo alternative polyadenylation (APA), producing distinct mRNA isoforms with different 3' untranslated regions (3' UTRs) and in some cases different coding regions. Two thirds of all human genes undergo APA. The efficiency of the polyadenylation process regulates gene expression and APA plays an important part in post-transcriptional regulation, as the 3' UTR contains various cis-elements associated with post-transcriptional regulation, such as target sites for micro-RNAs and RNA-binding proteins. Implications of alterations in polyadenylation for endocrine disease: Alterations in polyadenylation have been found to be causative of neonatal diabetes and IPEX (immune dysfunction, polyendocrinopathy, enteropathy, X-linked) and to be associated with type I and II diabetes, pre-eclampsia, fragile X-associated premature ovarian insufficiency, ectopic Cushing syndrome, and many cancer diseases, including several types of endocrine tumor diseases.

PERSPECTIVES

Recent developments in high-throughput sequencing have made it possible to characterize polyadenylation genome-wide. Antisense elements inhibiting or enhancing specific poly(A) site usage can induce desired alterations in polyadenylation, and thus hold the promise of new therapeutic approaches.

SUMMARY

This review gives a detailed description of alterations in polyadenylation in endocrine disease, an overview of the current literature on polyadenylation and summarizes the clinical implications of the current state of research in this field.

Alternative polyadenylation: new insights from global analyses

Recent studies have revealed widespread mRNA alternative polyadenylation (APA) in eukaryotes and its dynamic spatial and temporal regulation. APA not only generates proteomic and functional diversity, but also plays important roles in regulating gene expression. Global deregulation of APA has been demonstrated in a variety of human diseases. Recent exciting advances in the field have been made possible in a large part by high throughput analyses using newly developed experimental tools. Here I review the recent progress in global studies of APA and the insights that have emerged from these and other studies that use more conventional methods.

Genetics and regulatory impact of alternative polyadenylation in human B-lymphoblastoid cells

Gene expression varies widely between individuals of a population, and regulatory change can underlie phenotypes of evolutionary and biomedical relevance. A key question in the field is how DNA sequence variants impact gene expression, with most mechanistic studies to date focused on the effects of genetic change on regulatory regions upstream of protein-coding sequence. By contrast, the role of RNA 3'-end processing in regulatory variation remains largely unknown, owing in part to the challenge of identifying functional elements in 3' untranslated regions. In this work, we conducted a genomic survey of transcript ends in lymphoblastoid cells from genetically distinct human individuals. Our analysis mapped the cis-regulatory architecture of 3' gene ends, finding that transcript end positions did not fall randomly in untranslated regions, but rather preferentially flanked the locations of 3' regulatory elements, including miRNA sites. The usage of these transcript length forms and motifs varied across human individuals, and polymorphisms in polyadenylation signals and other 3' motifs were significant predictors of expression levels of the genes in which they lay. Independent single-gene experiments confirmed the effects of polyadenylation variants on steady-state expression of their respective genes, and validated the regulatory function of 3' cis-regulatory sequence elements that mediated expression of these distinct RNA length forms. Focusing on the immune regulator IRF5, we established the effect of natural variation in RNA 3'-end processing on regulatory response to antigen stimulation. Our results underscore the importance of two mechanisms at play in the genetics of 3'-end variation: the usage of distinct 3'-end processing signals and the effects of 3' sequence elements that determine transcript fate. Our findings suggest that the strategy of integrating observed 3'-end positions with inferred 3' regulatory motifs will prove to be a critical tool in continued efforts to interpret human genome variation.

Single nucleotide polymorphisms can create alternative polyadenylation signals and affect gene expression through loss of microRNA-regulation

Alternative polyadenylation (APA) can for example occur when a protein-coding gene has several polyadenylation (polyA) signals in its last exon, resulting in messenger RNAs (mRNAs) with different 3′ untranslated region (UTR) lengths. Different 3′UTR lengths can give different microRNA (miRNA) regulation such that shortened transcripts have increased expression. The APA process is part of human cells' natural regulatory processes, but APA also seems to play an important role in many human diseases. Although altered APA in disease can have many causes, we reasoned that mutations in DNA elements that are important for the polyA process, such as the polyA signal and the downstream GU-rich region, can be one important mechanism. To test this hypothesis, we identified single nucleotide polymorphisms (SNPs) that can create or disrupt APA signals (APA-SNPs). By using a data-integrative approach, we show that APA-SNPs can affect 3′UTR length, miRNA regulation, and mRNA expression—both between homozygote individuals and within heterozygote individuals. Furthermore, we show that a significant fraction of the alleles that cause APA are strongly and positively linked with alleles found by genome-wide studies to be associated with disease. Our results confirm that APA-SNPs can give altered gene regulation and that APA alleles that give shortened transcripts and increased gene expression can be important hereditary causes for disease.

Variants in DNA that affect gene expression—so-called regulatory variants—are thought to play important roles in common complex diseases, such as cancer. In contrast to variants in protein-coding regions, regulatory variants do not affect protein sequence and function. Instead, regulatory variants affect the amount of protein produced. The 3′ untranslated region (UTR) is one gene region that is critically important for gene regulation; cancers for example, often express genes with shortened 3′UTRs that, compared with full-length 3′UTRs, have higher and more stable expression levels. We have investigated one kind of regulatory variant that can affect the 3′UTR length and thereby cause disease. We identified several such variants in different genes and found that these variants affected the genes' expression. Some of these variants were also strongly linked with known markers for disease, suggesting that these regulatory variants are important hereditary causes for disease.

Mechanisms and consequences of alternative polyadenylation

Alternative polyadenylation (APA) is emerging as a widespread mechanism used to control gene expression. Like alternative splicing, usage of alternative poly(A) sites allows a single gene to encode multiple mRNA transcripts. In some cases, this changes the mRNA coding potential; in other cases, the code remains unchanged but the 3' UTR length is altered, influencing the fate of mRNAs in several ways, for example, by altering the availability of RNA binding protein sites and microRNA binding sites. The mechanisms governing both global and gene-specific APA are only starting to be deciphered. Here we review what is known about these mechanisms and the functional consequences of alternative polyadenylation.

Voltage-dependant anion channels: novel insights into isoform function through genetic models

Voltage-dependant Anion Channels, also known as mitochondrial porins, are pore-forming proteins located in the mitochondrial outer membrane (MOM) that, in addition to forming complexes with other proteins that localize to the MOM, also function as the main conduit for transporting metabolites between the cytoplasm and mitochondria. VDACs are encoded by a multi-member gene family, and the number of isoforms and specific functions of VDACs varies between species. Translating the well-described in vitro characteristics of the VDAC isoforms into in vivo functions has been a challenge, with the generation of animal models of VDAC deficiency providing much of the available information about isoform-specific roles in biology. Here, we review the approaches used to create these insect and mammalian animal models, and the conclusions reached by studying the consequences of loss of function mutations on the genetic, physiologic, and biochemical properties of the resulting models. This article is part of a Special Issue entitled: VDAC structure, function, and regulation of mitochondrial metabolism.

Can the microRNA signature distinguish between thyroid tumors of uncertain malignant potential and other well-differentiated tumors of the thyroid gland?

The term 'thyroid tumors of uncertain malignant potential' (TT-UMP) was coined by surgical pathologists to define well-differentiated tumors (WDT) showing inconclusive morphological evidence of malignancy or benignity. We have analyzed the expression of microRNA (miRNA) in a training set of 42 WDT of different histological subtypes: seven follicular tumors of UMP (FT-UMP), six WDT-UMP, seven follicular thyroid adenomas (FTA), 11 conventional papillary thyroid carcinomas (C-PTC), five follicular variants of PTC (FV-PTC), and six follicular thyroid carcinomas (FTC), which led to the identification of about 40 deregulated miRNAs. A subset of these altered miRNAs was independently validated by qRT-PCR, which included 18 supplementary TT-UMP (eight WDT-UMP and ten FT-UMP). Supervised clustering techniques were used to predict the first 42 samples. Based on the four possible outcomes (FTA, C-PTC, FV-PTC, and FTC), about 80% of FTA and C-PTC and 50% of FV-PTC and FTC samples were correctly assigned. Analysis of the independent set of 18 WDT-UMP by quantitative RT-PCR for the selection of the six most discriminating miRNAs was unable to separate FT-UMP from WDT-UMP, suggesting that the miRNA signature is insufficient in characterizing these two clinical entities. We conclude that considering FT-UMP and WDT-UMP as distinct and specific clinical entities may improve the diagnosis of WDT of the thyroid gland. In this context, a small set of miRNAs (i.e. miR-7, miR-146a, miR-146b, miR-200b, miR-221, and miR-222) appears to be useful, though not sufficient per se, in distinguishing TT-UMP from other WDT of the thyroid gland.

MicroRNA regulation by RNA-binding proteins and its implications for cancer

Non-protein-coding transcripts have been conserved throughout evolution, indicating that crucial functions exist for these RNAs. For example, microRNAs (miRNAs) have been found to modulate most cellular processes. The protein classes of RNA-binding proteins include essential regulators of miRNA biogenesis, turnover and activity. RNA-RNA and protein-RNA interactions are essential for post-transcriptional regulation in normal development and may be deregulated in disease. In reviewing emerging concepts of the interplay between miRNAs and RNA-binding proteins, we highlight the implications of these complex layers of regulation in cancer initiation and progression.

Integration of expressed sequence tag data flanking predicted RNA secondary structures facilitates novel non-coding RNA discovery

Many computational methods have been used to predict novel non-coding RNAs (ncRNAs), but none, to our knowledge, have explicitly investigated the impact of integrating existing cDNA-based Expressed Sequence Tag (EST) data that flank structural RNA predictions. To determine whether flanking EST data can assist in microRNA (miRNA) prediction, we identified genomic sites encoding putative miRNAs by combining functional RNA predictions with flanking ESTs data in a model consistent with miRNAs undergoing cleavage during maturation. In both human and mouse genomes, we observed that the inclusion of flanking ESTs adjacent to and not overlapping predicted miRNAs significantly improved the performance of various methods of miRNA prediction, including direct high-throughput sequencing of small RNA libraries. We analyzed the expression of hundreds of miRNAs predicted to be expressed during myogenic differentiation using a customized microarray and identified several known and predicted myogenic miRNA hairpins. Our results indicate that integrating ESTs flanking structural RNA predictions improves the quality of cleaved miRNA predictions and suggest that this strategy can be used to predict other non-coding RNAs undergoing cleavage during maturation.

MicroRNAs in Drosophila development

Micro-ribonucleic acids (miRNAs) are small (21-24 nucleotide), endogenously expressed, noncoding RNAs that have emerged as important posttranscriptional regulators of gene expression. MiRNAs have been identified and cloned from diverse eukaryotic organisms where they have been shown to control important physiological and developmental processes such as apoptosis, cell division, and differentiation. A high level of conservation of some miRNAs across phyla further emphasizes their importance as posttranscriptional regulators. Research in a variety of model systems has been instrumental in dissecting the biological functions of miRNAs. In this chapter, we discuss the current literature on the role of miRNAs as developmental regulators in Drosophila.

Alternative mRNA polyadenylation in eukaryotes: an effective regulator of gene expression

Alternative RNA processing mechanisms, including alternative splicing and alternative polyadenylation, are increasingly recognized as important regulators of gene expression. This article will focus on what has recently been described about alternative polyadenylation in development, differentiation, and disease in higher eukaryotes. We will also describe how the evolving global methodologies for examining the cellular transcriptome, both experimental and bioinformatic, are revealing new details about the complex nature of alternative 3' end formation, as well as interactions with other RNA-mediated and RNA processing mechanisms.

Impact of probe annotation on the integration of miRNA-mRNA expression profiles for miRNA target detection

MicroRNAs (miRNAs) are small non-coding RNAs that mediate gene expression at the post-transcriptional and translational levels by an imperfect binding to target mRNA 3′UTR regions. While the ab-initio computational prediction of miRNA–mRNA interactions still poses significant challenges, it is possible to overcome some of its limitations by carefully integrating into the analysis the paired expression profiles of miRNAs and mRNAs. In this work, we show how the choice of a proper probe annotation for microarray platforms is an essential requirement to achieve good sensitivity in the identification of miRNA–mRNA interactions. We compare the results obtained from the analysis of the same expression profiles using both gene and transcript based custom CDFs that we have developed for a number of different annotations (ENSEMBL, RefSeq, AceView). In all cases, transcript-based annotations clearly improve the effectiveness of data integration and thus provide a more reliable confirmation of computationally predicted miRNA–mRNA interactions.

Reprogramming of 3' untranslated regions of mRNAs by alternative polyadenylation in generation of pluripotent stem cells from different cell types

BACKGROUND

The 3' untranslated regions (3'UTRs) of mRNAs contain cis elements involved in post-transcriptional regulation of gene expression. Over half of all mammalian genes contain multiple polyadenylation sites that lead to different 3'UTRs for a gene. Studies have shown that the alternative polyadenylation (APA) pattern varies across tissues, and is dynamically regulated in proliferating or differentiating cells. Generation of induced pluripotent stem (iPS) cells, in which differentiated cells are reprogrammed to an embryonic stem (ES) cell-like state, has been intensively studied in recent years. However, it is not known how 3'UTRs are regulated during cell reprogramming.

METHODS/MAIN FINDINGS

Using a computational method that robustly examines APA across DNA microarray data sets, we analyzed 3'UTR dynamics in generation of iPS cells from different cell types. We found that 3'UTRs shorten during reprogramming of somatic cells, the extent of which depends on the type of source cell. By contrast, reprogramming of spermatogonial cells involves 3'UTR lengthening. The alternative polyadenylation sites that are highly responsive to change of cell state in generation of iPS cells are also highly regulated during embryonic development in opposite directions. Compared with other sites, they are more conserved, can lead to longer alternative 3'UTRs, and are associated with more cis elements for polyadenylation. Consistently, reprogramming of somatic cells and germ cells involves significant upregulation and downregulation, respectively, of mRNAs encoding polyadenylation factors, and RNA processing is one of the most significantly regulated biological processes during cell reprogramming. Furthermore, genes containing target sites of ES cell-specific microRNAs (miRNAs) in different portions of 3'UTR are distinctively regulated during cell reprogramming, suggesting impact of APA on miRNA targeting.

CONCLUSIONS/SIGNIFICANCE

Taken together, these findings indicate that reprogramming of 3'UTRs by APA, which result from regulation of both general polyadenylation activity and cell type-specific factors and can reset post-transcriptional gene regulatory programs in the cell, is an integral part of iPS cell generation, and the APA pattern can be a good biomarker for cell type and state, useful for sample classification. The results also suggest that perturbation of the mRNA polyadenylation machinery or RNA processing activity may facilitate generation of iPS cells.

Characterization of full-length sequenced cDNA inserts (FLIcs) from Atlantic salmon (Salmo salar)

Background

Sequencing of the Atlantic salmon genome is now being planned by an international research consortium. Full-length sequenced inserts from cDNAs (FLIcs) are an important tool for correct annotation and clustering of the genomic sequence in any species. The large amount of highly similar duplicate sequences caused by the relatively recent genome duplication in the salmonid ancestor represents a particular challenge for the genome project. FLIcs will therefore be an extremely useful resource for the Atlantic salmon sequencing project. In addition to be helpful in order to distinguish between duplicate genome regions and in determining correct gene structures, FLIcs are an important resource for functional genomic studies and for investigation of regulatory elements controlling gene expression. In contrast to the large number of ESTs available, including the ESTs from 23 developmental and tissue specific cDNA libraries contributed by the Salmon Genome Project (SGP), the number of sequences where the full-length of the cDNA insert has been determined has been small.

Results

High quality full-length insert sequences from 560 pre-smolt white muscle tissue specific cDNAs were generated, accession numbers [GenBank: BT043497 - BT044056]. Five hundred and ten (91%) of the transcripts were annotated using Gene Ontology (GO) terms and 440 of the FLIcs are likely to contain a complete coding sequence (cCDS). The sequence information was used to identify putative paralogs, characterize salmon Kozak motifs, polyadenylation signal variation and to identify motifs likely to be involved in the regulation of particular genes. Finally, conserved 7-mers in the 3'UTRs were identified, of which some were identical to miRNA target sequences.

Conclusion

This paper describes the first Atlantic salmon FLIcs from a tissue and developmental stage specific cDNA library. We have demonstrated that many FLIcs contained a complete coding sequence (cCDS). This suggests that the remaining cDNA libraries generated by SGP represent a valuable cCDS FLIc source. The conservation of 7-mers in 3'UTRs indicates that these motifs are functionally important. Identity between some of these 7-mers and miRNA target sequences suggests that they are miRNA targets in Salmo salar transcripts as well.

Progressive lengthening of 3' untranslated regions of mRNAs by alternative polyadenylation during mouse embryonic development

The 3' untranslated regions (3' UTRs) of mRNAs contain cis-acting elements for posttranscriptional regulation of gene expression. Here, we report that mouse genes tend to express mRNAs with longer 3' UTRs as embryonic development progresses. This global regulation is controlled by alternative polyadenylation and coordinates with initiation of organogenesis and aspects of embryonic development, including morphogenesis, differentiation, and proliferation. Using myogenesis of C2C12 myoblast cells as a model, we recapitulated this process in vitro and found that 3' UTR lengthening is likely caused by weakening of mRNA polyadenylation activity. Because alternative 3' UTR sequences are typically longer and have higher AU content than constitutive ones, our results suggest that lengthening of 3' UTR can significantly augment posttranscriptional control of gene expression during embryonic development, such as microRNA-mediated regulation.

Genomic structure and expression analysis of the RNase kappa family ortholog gene in the insect Ceratitis capitata

Cc RNase is the founding member of the recently identified RNase kappa family, which is represented by a single ortholog in a wide range of animal taxonomic groups. Although the precise biological role of this protein is still unknown, it has been shown that the recombinant proteins isolated so far from the insect Ceratitis capitata and from human exhibit ribonucleolytic activity. In this work, we report the genomic organization and molecular evolution of the RNase kappa gene from various animal species, as well as expression analysis of the ortholog gene in C. capitata. The high degree of amino acid sequence similarity, in combination with the fact that exon sizes and intronic positions are extremely conserved among RNase kappa orthologs in 15 diverse genomes from sea anemone to human, imply a very significant biological function for this enzyme. In C. capitata, two forms of RNase kappa mRNA (0.9 and 1.5 kb) with various lengths of 3' UTR were identified as alternative products of a single gene, resulting from the use of different polyadenylation signals. Both transcripts are expressed in all insect tissues and developmental stages. Sequence analysis of the extended region of the longer transcript revealed the existence of three mRNA instability motifs (AUUUA) and five poly(U) tracts, whose functional importance in RNase kappa mRNA decay remains to be explored.

A PASport to cellular proliferation

The use of alternative polyadenylation sites produces mRNA isoforms with different 3' untranslated regions. A recent report in Science (Sandberg et al., 2008) suggests that alternative polyadenylation is connected to microRNA-mediated regulation of gene expression as part of a global program for cellular proliferation.

Posttranscriptional regulation of BK channel splice variant stability by miR-9 underlies neuroadaptation to alcohol

Tolerance represents a critical component of addiction. The large-conductance calcium- and voltage-activated potassium channel (BK) is a well-established alcohol target, and an important element in behavioral and molecular alcohol tolerance. We tested whether microRNA, a newly discovered class of gene expression regulators, plays a role in the development of tolerance. We show that in adult mammalian brain, alcohol upregulates microRNA miR-9 and mediates posttranscriptional reorganization in BK mRNA splice variants by miR-9-dependent destabilization of BK mRNAs containing 3'UTRs with a miR-9 Recognition Element (MRE). Different splice variants encode BK isoforms with different alcohol sensitivities. Computational modeling indicates that this miR-9-dependent mechanism contributes to alcohol tolerance. Moreover, this mechanism can be extended to include regulation of additional miR-9 targets relevant to alcohol abuse. Our results describe a mechanism of multiplex regulation of stability of alternatively spliced mRNA by microRNA in drug adaptation and neuronal plasticity.

Regulation of ABCG2 expression at the 3' untranslated region of its mRNA through modulation of transcript stability and protein translation by a putative microRNA in the S1 colon cancer cell line

ABCG2 is recognized as an important efflux transporter in clinical pharmacology and is potentially important in resistance to chemotherapeutic drugs. To identify epigenetic mechanisms regulating ABCG2 mRNA expression at its 3' untranslated region (3'UTR), we performed 3' rapid amplification of cDNA ends with the S1 parental colon cancer cell line and its drug-resistant ABCG2-overexpressing counterpart. We found that the 3'UTR is >1,500 bp longer in parental cells and, using the miRBase TARGETs database, identified a putative microRNA (miRNA) binding site, distinct from the recently reported hsa-miR520h site, in the portion of the 3'UTR missing from ABCG2 mRNA in the resistant cells. We hypothesized that the binding of a putative miRNA at the 3'UTR of ABCG2 suppresses the expression of ABCG2. In resistant S1MI80 cells, the miRNA cannot bind to ABCG2 mRNA because of the shorter 3'UTR, and thus, mRNA degradation and/or repression on protein translation is relieved, contributing to overexpression of ABCG2. This hypothesis was rigorously tested by reporter gene assays, mutational analysis at the miRNA binding sites, and forced expression of miRNA inhibitors or mimics. The removal of this epigenetic regulation by miRNA could be involved in the overexpression of ABCG2 in drug-resistant cancer cells.

Proliferating cells express mRNAs with shortened 3' untranslated regions and fewer microRNA target sites

Messenger RNA (mRNA) stability, localization, and translation are largely determined by sequences in the 3' untranslated region (3'UTR). We found a conserved increase in expression of mRNAs terminating at upstream polyadenylation sites after activation of primary murine CD4+ T lymphocytes. This program, resulting in shorter 3'UTRs, is a characteristic of gene expression during immune cell activation and correlates with proliferation across diverse cell types and tissues. Forced expression of full-length 3'UTRs conferred reduced protein expression. In some cases the reduction in protein expression could be reversed by deletion of predicted microRNA target sites in the variably included region. Our data indicate that gene expression is coordinately regulated, such that states of increased proliferation are associated with widespread reductions in the 3'UTR-based regulatory capacity of mRNAs.

Proliferating cells express mRNAs with shortened 3' untranslated regions and fewer microRNA target sites

Messenger RNA (mRNA) stability, localization, and translation are largely determined by sequences in the 3' untranslated region (3'UTR). We found a conserved increase in expression of mRNAs terminating at upstream polyadenylation sites after activation of primary murine CD4+ T lymphocytes. This program, resulting in shorter 3'UTRs, is a characteristic of gene expression during immune cell activation and correlates with proliferation across diverse cell types and tissues. Forced expression of full-length 3'UTRs conferred reduced protein expression. In some cases the reduction in protein expression could be reversed by deletion of predicted microRNA target sites in the variably included region. Our data indicate that gene expression is coordinately regulated, such that states of increased proliferation are associated with widespread reductions in the 3'UTR-based regulatory capacity of mRNAs.

3' end mRNA processing: molecular mechanisms and implications for health and disease

Recent advances in the understanding of the molecular mechanism of mRNA 3' end processing have uncovered a previously unanticipated integrated network of transcriptional and RNA-processing mechanisms. A variety of human diseases impressively reflect the importance of the precision of the complex 3' end-processing machinery and gene specific deregulation of 3' end processing can result from mutations of RNA sequence elements that bind key specific processing factors. Interestingly, more general deregulation of 3' end processing can be caused either by mutations of these processing factors or by the disturbance of the well-coordinated equilibrium between these factors. From a medical perspective, both loss of function and gain of function can be functionally relevant, and an increasing number of different disease entities exemplifies that inappropriate 3' end formation of human mRNAs can have a tremendous impact on health and disease. Here, we review the mechanistic hallmarks of mRNA 3' end processing, highlight the medical relevance of deregulation of this important step of mRNA maturation and illustrate the implications for diagnostic and therapeutic strategies.

MicroRNAs in carcinogenesis

MicroRNAs are an abundant class of noncoding RNAs, typically 20-23 nucleotides in length that are often evolutionarily conserved in metazoans and expressed in a cell and tissue specific manner. MicroRNAs exert their gene regulatory activity primarily by imperfectly base pairing to the 3' UTR of their target mRNAs, leading to mRNA degradation or translational inhibition. In cancer, microRNAs are often dysregulated with their expression patterns being correlated with clinically relevant tumor characteristics. Recently, microRNAs were shown to be directly involved in cancer initiation and progression. This review focuses primarily on emerging developments in the microRNA field that impact our understanding of how these molecules contribute to carcinogenesis.

Gene expression trees in lymphoid development

BACKGROUND

The regulatory processes that govern cell proliferation and differentiation are central to developmental biology. Particularly well studied in this respect is the lymphoid system due to its importance for basic biology and for clinical applications. Gene expression measured in lymphoid cells in several distinguishable developmental stages helps in the elucidation of underlying molecular processes, which change gradually over time and lock cells in either the B cell, T cell or Natural Killer cell lineages. Large-scale analysis of these gene expression trees requires computational support for tasks ranging from visualization, querying, and finding clusters of similar genes, to answering detailed questions about the functional roles of individual genes.

RESULTS

We present the first statistical framework designed to analyze gene expression data as it is collected in the course of lymphoid development through clusters of co-expressed genes and additional heterogeneous data. We introduce dependence trees for continuous variates, which model the inherent dependencies during the differentiation process naturally as gene expression trees. Several trees are combined in a mixture model to allow inference of potentially overlapping clusters of co-expressed genes. Additionally, we predict microRNA targets.

CONCLUSION

Computational results for several data sets from the lymphoid system demonstrate the relevance of our framework. We recover well-known biological facts and identify promising novel regulatory elements of genes and their functional assignments. The implementation of our method (licensed under the GPL) is available at http://algorithmics.molgen.mpg.de/Supplements/ExpLym/.

Spatial preferences of microRNA targets in 3' untranslated regions

BACKGROUND

MicroRNAs are an important class of regulatory RNAs which repress animal genes by preferentially interacting with complementary sequence motifs in the 3' untranslated region (UTR) of target mRNAs. Computational methods have been developed which can successfully predict which microRNA may target which mRNA on a genome-wide scale.

RESULTS

We address how predicted target sites may be affected by alternative polyadenylation events changing the 3'UTR sequence. We find that two thirds of targeted genes have alternative 3'UTRs, with 40% of predicted target sites located in alternative UTR segments. We propose three classes based on whether the target sites fall within constitutive and/or alternative UTR segments, and examine the spatial distribution of predicted targets in alternative UTRs. In particular, there is a strong preference for targets to be located in close vicinity of the stop codon and the polyadenylation sites.

CONCLUSION

The transcript diversity seen in non-coding regions, as well as the relative location of miRNA target sites defined by it, has a potentially large impact on gene regulation by miRNAs and should be taken into account when defining, predicting or validating miRNA targets.