Comprehensive polyadenylation site maps in yeast and human reveal pervasive alternative polyadenylation.

The emerging discoveries on the link between polyadenylation and disease states underline the need to fully characterize genome-wide polyadenylation states. Here, we report comprehensive maps of global polyadenylation events in human and yeast generated using refinements to the Direct RNA Sequencing technology. This direct approach provides a quantitative view of genome-wide polyadenylation states in a strand-specific manner and requires only attomole RNA quantities. The polyadenylation profiles revealed an abundance of unannotated polyadenylation sites, alternative polyadenylation patterns, and regulatory element-associated poly(A)(+) RNAs. We observed differences in sequence composition surrounding canonical and noncanonical human polyadenylation sites, suggesting novel noncoding RNA-specific polyadenylation mechanisms in humans. Furthermore, we observed the correlation level between sense and antisense transcripts to depend on gene expression levels, supporting the view that overlapping transcription from opposite strands may play a regulatory role. Our data provide a comprehensive view of the polyadenylation state and overlapping transcription.

Noncoding RNAs that associate with YB-1 alter proliferation in prostate cancer cells.

The highly conserved, multifunctional YB-1 is a powerful breast cancer prognostic indicator. We report on a pervasive role for YB-1 in which it associates with thousands of nonpolyadenylated short RNAs (shyRNAs) that are further processed into small RNAs (smyRNAs). Many of these RNAs have previously been identified as functional noncoding RNAs (http://www.johnlab.org/YB1). We identified a novel, abundant, 3'-modified short RNA antisense to Dicer1 (Shad1) that colocalizes with YB-1 to P-bodies and stress granules. The expression of Shad1 was shown to correlate with that of YB-1 and whose inhibition leads to an increase in cell proliferation. Additionally, Shad1 influences the expression of additional prognostic markers of cancer progression such as DLX2 and IGFBP2. We propose that the examination of these noncoding RNAs could lead to better understanding of prostate cancer progression.

New class of gene-termini-associated human RNAs suggests a novel RNA copying mechanism.

Small (<200 nucleotide) RNA (sRNA) profiling of human cells using various technologies demonstrates unexpected complexity of sRNAs with hundreds of thousands of sRNA species present. Genetic and in vitro studies show that these RNAs are not merely degradation products of longer transcripts but could indeed have a function. Furthermore, profiling of RNAs, including the sRNAs, can reveal not only novel transcripts, but also make clear predictions about the existence and properties of novel biochemical pathways operating in a cell. For example, sRNA profiling in human cells indicated the existence of an unknown capping mechanism operating on cleaved RNA, a biochemical component of which was later identified. Here we show that human cells contain a novel type of sRNA that has non-genomically encoded 5' poly(U) tails. The presence of these RNAs at the termini of genes, specifically at the very 3' ends of known mRNAs, strongly argues for the presence of a yet uncharacterized endogenous biochemical pathway in cells that can copy RNA. We show that this pathway can operate on multiple genes, with specific enrichment towards transcript-encoding components of the translational machinery. Finally, we show that genes are also flanked by sense, 3' polyadenylated sRNAs that are likely to be capped.

miRNA-939 regulates human inducible nitric oxide synthase posttranscriptional gene expression in human hepatocytes.

Human inducible nitric oxide synthase (hiNOS) gene expression is regulated by transcriptional and posttranscriptional mechanisms. The purpose of this study was to determine whether specific microRNA (miRNA) directly regulate hiNOS gene expression. Sequence analysis of the 496-bp hiNOS 3'-untranslated region (3'-UTR) revealed five putative miR-939 binding sites. The hiNOS 3'-UTR conferred significant posttranscriptional blockade of luciferase activity in human A549, HCT8, and HeLa cells. The hiNOS 3'-UTR also exerted basal and cytokine-stimulated posttranscriptional repression in an orientation-dependent manner. Functional studies demonstrated that transfection of miR-939 into primary human hepatocytes (HCs) significantly inhibited cytokine-induced NO synthesis in a dose-dependent manner that was abrogated by a specific miR-939 inhibitor. MiR-939 (but not other miRNAs) abolished cytokine-stimulated hiNOS protein in human HC, but had no effect on hiNOS mRNA levels. Site-directed mutagenesis of miR-939 bindings sites at +99 or +112 bp in the hiNOS 3'-UTR increased reporter gene expression. Furthermore, intact miR-939 binding sites at +99 or +112 positions were required for posttranscriptional suppression by miR-939. Cytokine stimulation directly increased miR-939 levels in human HC. Transfection of miR-939 inhibitor (antisense miR-939) enhanced cytokine-induced hiNOS protein and increased NO synthesis in vitro in human HC. Finally, cytokine or LPS injection in vivo in mice increased hepatic miR-939 levels. Taken together, these data identify that miR-939 directly regulates hiNOS gene expression by binding in the 3'-UTR to produce a translational blockade. These findings suggest dual regulation of iNOS gene expression where cytokines induce iNOS transcription and also increase miR-939, leading to translational inhibition in a check-and-balance system.

Extensive terminal and asymmetric processing of small RNAs from rRNAs, snoRNAs, snRNAs, and tRNAs.

Deep sequencing studies frequently identify small RNA fragments of abundant RNAs. These fragments are thought to represent degradation products of their precursors. Using sequencing, computational analysis, and sensitive northern blot assays, we show that constitutively expressed non-coding RNAs such as tRNAs, snoRNAs, rRNAs and snRNAs preferentially produce small 5' and 3' end fragments. Similar to that of microRNA processing, these terminal fragments are generated in an asymmetric manner that predominantly favors either the 5' or 3' end. Terminal-specific and asymmetric processing of these small RNAs occurs in both mouse and human cells. In addition to the known processing of some 3' terminal tRNA-derived fragments (tRFs) by the RNase III endonuclease Dicer, we show that several RNase family members can produce tRFs, including Angiogenin that cleaves the TψC loop to generate 3' tRFs. The 3' terminal tRFs but not the 5' tRFs are highly complementary to human endogenous retroviral sequences in the genome. Despite their independence from Dicer processing, these tRFs associate with Ago2 and are capable of down regulating target genes by transcript cleavage in vitro. We suggest that endogenous 3' tRFs have a role in regulating the unwarranted expression of endogenous viruses through the RNA interference pathway.

An in-depth map of polyadenylation sites in cancer.

We present a comprehensive map of over 1 million polyadenylation sites and quantify their usage in major cancers and tumor cell lines using direct RNA sequencing. We built the Expression and Polyadenylation Database to enable the visualization of the polyadenylation maps in various cancers and to facilitate the discovery of novel genes and gene isoforms that are potentially important to tumorigenesis. Analyses of polyadenylation sites indicate that a large fraction (∼30%) of mRNAs contain alternative polyadenylation sites in their 3' untranslated regions, independent of the cell type. The shortest 3' untranslated region isoforms are preferentially upregulated in cancer tissues, genome-wide. Candidate targets of alternative polyadenylation-mediated upregulation of short isoforms include POLR2K, and signaling cascades of cell-cell and cell-extracellular matrix contact, particularly involving regulators of Rho GTPases. Polyadenylation maps also helped to improve 3' untranslated region annotations and identify candidate regulatory marks such as sequence motifs, H3K36Me3 and Pabpc1 that are isoform dependent and occur in a position-specific manner. In summary, these results highlight the need to go beyond monitoring only the cumulative transcript levels for a gene, to separately analysing the expression of its RNA isoforms.

An ultrastructural report of human articular cartilage resident cells in correlation with their phenotypic characteristics.

The recent discovery of progenitors based on their differential fibronectin-adhesion (FAA-CPs) and migratory-based (MCPs) assay has evoked interest due to their superiority in terms of their efficient chondrogenesis and reduced hypertrophic propensity. This study aims to isolate and enrich three articular cartilage subsets, chondrocytes, FAA-CPs, and MCPs, and compare their undifferentiated and chondrogenic differentiated status, using in-vitro phenotypical characterization in correlation with ultrastructural analysis using Transmission Electron Microscopy (TEM). Following informed consent, cartilage shavings were procured from a non-diseased human ankle joint and cultured to obtain the three subsets. Chondrocytes exhibited higher CD106 and lower CD49b and CD146 levels. Following chondrogenic differentiation, corroborative results were seen, with the MCP group showing the highest GAG/DNA ratio levels and uptake of extracellular matrix stain as compared to the FAA-CP group. TEM analysis of the chondrocytes revealed the presence of more autolytic cells with disintegrated cytoplasm and plasma membrane. The differentiated FAA-CPs and MCPs displayed higher collagen and rough endoplasmic reticulum. The results presented in this study provide novel information on the ultrastructural characteristics of cartilage resident cells, with the chondrocyte group displaying features of terminal differentiation. Both progenitor subtypes showed superiority in varied contexts, with greater collagen fibrils and greater GAG content in MCPs. The display of preferential and differentiation traits sheds insight on the necessity to enrich progenitors and coculturing them with the general pool of constituent cells to combine their advantages and reduce their drawbacks to achieve a regenerative tissue displaying genuine hyaline-like repair while limiting their terminal differentiation.

A sensitive non-radioactive northern blot method to detect small RNAs.

The continuing discoveries of potentially active small RNAs at an unprecedented rate using high-throughput sequencing have raised the need for methods that can reliably detect and quantitate the expression levels of small RNAs. Currently, northern blot is the most widely used method for validating small RNAs that are identified by methods such as high-throughput sequencing. We describe a new northern blot-based protocol (LED) for small RNA (approximately 15-40 bases) detection using digoxigenin (DIG)-labeled oligonucleotide probes containing locked nucleic acids (LNA) and 1-ethyl-3-(3-dimethylaminopropyl) carbodiimide for cross-linking the RNA to the membrane. LED generates clearly visible signals for RNA amounts as low as 0.05 fmol. This method requires as little as a few seconds of membrane exposure to outperform the signal intensity using overnight exposure of isotope-based methods, corresponding to approximately 1000-fold improvement in exposure-time. In contrast to commonly used radioisotope-based methods, which require freshly prepared and hazardous probes, LED probes can be stored for at least 6 months, facilitate faster and more cost-effective experiments, and are more environmentally friendly. A detailed protocol of LED is provided in the Supplementary Data.

Characterization of viral and human RNAs smaller than canonical MicroRNAs.

Recently identified small (20 to 40 bases) RNAs, such as microRNAs (miRNAs) and endogenous small interfering RNAs (siRNAs) participate in important cellular pathways. In this report, we systematically characterized several novel features of human and viral RNA products smaller than miRNAs. We found that Kaposi sarcoma-associated herpesvirus K12-1 miRNA (23 bases) associates with a distinct, unusually small (17-base) RNA (usRNA) that can effectively downregulate a K12-1 miRNA target, human RAD21, suggesting that stable degradation-like products may also contribute to gene regulation. High-throughput sequencing reveals a diverse set of human miRNA-derived usRNAs and other non-miRNA-derived usRNAs. Human miRNA-derived usRNAs preferentially match to 5' ends of miRNAs and are also more likely to associate with the siRNA effector protein Ago2 than with Ago1. Many non-miRNA-derived usRNAs associate with Ago proteins and also frequently contain C-rich 3'-specific motifs that are overrepresented in comparison to Piwi-interacting RNAs and transcription start site-associated RNAs. We postulate that approximately 30% of usRNAs could have evolved to participate in biological processes, including gene silencing.

Regulatory circuit of human microRNA biogenesis.

miRNAs (microRNAs) are a class of endogenous small RNAs that are thought to negatively regulate protein production. Aberrant expression of many miRNAs is linked to cancer and other diseases. Little is known about the factors that regulate the expression of miRNAs. We have identified numerous regulatory elements upstream of miRNA genes that are likely to be essential to the transcriptional and posttranscriptional regulation of miRNAs. Newly identified regulatory motifs occur frequently and in multiple copies upstream of miRNAs. The motifs are highly enriched in G and C nucleotides, in comparison with the nucleotide composition of miRNA upstream sequences. Although the motifs were predicted using sequences that are upstream of miRNAs, we find that 99% of the top-predicted motifs preferentially occur within the first 500 nucleotides upstream of the transcription start sites of protein-coding genes; the observed preference in location underscores the validity and importance of the motifs identified in this study. Our study also raises the possibility that a considerable number of well-characterized, disease-associated transcription factors (TFs) of protein-coding genes contribute to the abnormal miRNA expression in diseases such as cancer. Further analysis of predicted miRNA-protein interactions lead us to hypothesize that TFs that include c-Myb, NF-Y, Sp-1, MTF-1, and AP-2alpha are master-regulators of miRNA expression. Our predictions are a solid starting point for the systematic elucidation of the causative basis for aberrant expression patterns of disease-related (e.g., cancer) miRNAs. Thus, we point out that focused studies of the TFs that regulate miRNAs will be paramount in developing cures for miRNA-related diseases. The identification of the miRNA regulatory motifs was facilitated by a new computational method, K-Factor. K-Factor predicts regulatory motifs in a set of functionally related sequences, without relying on evolutionary conservation.

Human MicroRNA targets.

MicroRNAs (miRNAs) interact with target mRNAs at specific sites to induce cleavage of the message or inhibit translation. The specific function of most mammalian miRNAs is unknown. We have predicted target sites on the 3' untranslated regions of human gene transcripts for all currently known 218 mammalian miRNAs to facilitate focused experiments. We report about 2,000 human genes with miRNA target sites conserved in mammals and about 250 human genes conserved as targets between mammals and fish. The prediction algorithm optimizes sequence complementarity using position-specific rules and relies on strict requirements of interspecies conservation. Experimental support for the validity of the method comes from known targets and from strong enrichment of predicted targets in mRNAs associated with the fragile X mental retardation protein in mammals. This is consistent with the hypothesis that miRNAs act as sequence-specific adaptors in the interaction of ribonuclear particles with translationally regulated messages. Overrepresented groups of targets include mRNAs coding for transcription factors, components of the miRNA machinery, and other proteins involved in translational regulation, as well as components of the ubiquitin machinery, representing novel feedback loops in gene regulation. Detailed information about target genes, target processes, and open-source software for target prediction (miRanda) is available at http://www.microrna.org. Our analysis suggests that miRNA genes, which are about 1% of all human genes, regulate protein production for 10% or more of all human genes.

Prediction of human microRNA targets.

MicroRNAs (miRNAs) are small, nonprotein-coding RNAs that regulate gene expression. Although hundreds of human miRNA genes have been discovered, the functions of most of these are unknown. Computational predictions indicate that miRNAs, which account for at least 1% of human protein-coding genes, regulate protein production for thousands of or possibly all of human genes. We discuss the functions of mammalian miRNAs and the experimental and computational methods used to detect and predict human miRNA target genes. Anticipating their impact on genome-wide discovery of miRNA targets, we describe the various computational tools and web-based resources available to predict miRNA targets.

Comparative protein structure modeling by iterative alignment, model building and model assessment.

Comparative or homology protein structure modeling is severely limited by errors in the alignment of a modeled sequence with related proteins of known three-dimensional structure. To ameliorate this problem, we have developed an automated method that optimizes both the alignment and the model implied by it. This task is achieved by a genetic algorithm protocol that starts with a set of initial alignments and then iterates through re-alignment, model building and model assessment to optimize a model assessment score. During this iterative process: (i) new alignments are constructed by application of a number of operators, such as alignment mutations and cross-overs; (ii) comparative models corresponding to these alignments are built by satisfaction of spatial restraints, as implemented in our program MODELLER; (iii) the models are assessed by a variety of criteria, partly depending on an atomic statistical potential. When testing the procedure on a very difficult set of 19 modeling targets sharing only 4-27% sequence identity with their template structures, the average final alignment accuracy increased from 37 to 45% relative to the initial alignment (the alignment accuracy was measured as the percentage of positions in the tested alignment that were identical to the reference structure-based alignment). Correspondingly, the average model accuracy increased from 43 to 54% (the model accuracy was measured as the percentage of the C(alpha) atoms of the model that were within 5 A of the corresponding C(alpha) atoms in the superposed native structure). The present method also compares favorably with two of the most successful previously described methods, PSI-BLAST and SAM. The accuracy of the final models would be increased further if a better method for ranking of the models were available.

Tools for comparative protein structure modeling and analysis.

The following resources for comparative protein structure modeling and analysis are described (http://salilab.org): MODELLER, a program for comparative modeling by satisfaction of spatial restraints; MODWEB, a web server for automated comparative modeling that relies on PSI-BLAST, IMPALA and MODELLER; MODLOOP, a web server for automated loop modeling that relies on MODELLER; MOULDER, a CPU intensive protocol of MODWEB for building comparative models based on distant known structures; MODBASE, a comprehensive database of annotated comparative models for all sequences detectably related to a known structure; MODVIEW, a Netscape plugin for Linux that integrates viewing of multiple sequences and structures; and SNPWEB, a web server for structure-based prediction of the functional impact of a single amino acid substitution.

Genome-wide transcript profiling reveals novel breast cancer-associated intronic sense RNAs.

Non-coding RNAs (ncRNAs) play major roles in development and cancer progression. To identify novel ncRNAs that may identify key pathways in breast cancer development, we performed high-throughput transcript profiling of tumor and normal matched-pair tissue samples. Initial transcriptome profiling using high-density genome-wide tiling arrays revealed changes in over 200 novel candidate genomic regions that map to intronic regions. Sixteen genomic loci were identified that map to the long introns of five key protein-coding genes, CRIM1, EPAS1, ZEB2, RBMS1, and RFX2. Consistent with the known role of the tumor suppressor ZEB2 in the cancer-associated epithelial to mesenchymal transition (EMT), in situ hybridization reveals that the intronic regions deriving from ZEB2 as well as those from RFX2 and EPAS1 are down-regulated in cells of epithelial morphology, suggesting that these regions may be important for maintaining normal epithelial cell morphology. Paired-end deep sequencing analysis reveals a large number of distinct genomic clusters with no coding potential within the introns of these genes. These novel transcripts are only transcribed from the coding strand. A comprehensive search for breast cancer associated genes reveals enrichment for transcribed intronic regions from these loci, pointing to an underappreciated role of introns or mechanisms relating to their biology in EMT and breast cancer.

Detection of homologous proteins by an intermediate sequence search.

We developed a variant of the intermediate sequence search method (ISS(new)) for detection and alignment of weakly similar pairs of protein sequences. ISS(new) relates two query sequences by an intermediate sequence that is potentially homologous to both queries. The improvement was achieved by a more robust overlap score for a match between the queries through an intermediate. The approach was benchmarked on a data set of 2369 sequences of known structure with insignificant sequence similarity to each other (BLAST E-value larger than 0.001); 2050 of these sequences had a related structure in the set. ISS(new) performed significantly better than both PSI-BLAST and a previously described intermediate sequence search method. PSI-BLAST could not detect correct homologs for 1619 of the 2369 sequences. In contrast, ISS(new) assigned a correct homolog as the top hit for 121 of these 1619 sequences, while incorrectly assigning homologs for only nine targets; it did not assign homologs for the remainder of the sequences. By estimate, ISS(new) may be able to assign the folds of domains in approximately 29,000 of the approximately 500,000 sequences unassigned by PSI-BLAST, with 90% specificity (1 - false positives fraction). In addition, we show that the 15 alignments with the most significant BLAST E-values include the nearly best alignments constructed by ISS(new).

Enhanced detection of small RNAs using a nonradioactive approach.

Recent advancements in high-throughput sequencing have led to the identification of many new classes of small noncoding RNAs such as endo-siRNAs. Unfortunately, reliable quantification of RNAs by sequencing is difficult due to artifacts arising from various factors involved in cDNA library preparation. Northern blot is one of the leading methods used to confirm the presence of a given RNA sequence because it can accurately quantify the cellular abundance, the size of the small RNA and reveal the presence of potential precursors and RNA isoforms. Here, we present a comprehensive description of LNA probe design along with a recently developed highly sensitive and cost-effective nonradioactive northern blot approach termed LED. LED combines a cross-linking method (EDC) and digoxigenin (DIG) labeling, and it can detect small RNAs with concentrations as low as 0.05 fmol and requires as little as a few seconds of membrane exposure for signal generation.

Structural characterization of components of protein assemblies by comparative modeling and electron cryo-microscopy.

We explore structural characterization of protein assemblies by a combination of electron cryo-microscopy (cryoEM) and comparative protein structure modeling. Specifically, our method finds an optimal atomic model of a given assembly subunit and its position within an assembly by fitting alternative comparative models into a cryoEM map. The alternative models are calculated by MODELLER [J. Mol. Biol. 234 (1993) 313] from different sequence alignments between the modeled protein and its template structures. The fitting of these models into a cryoEM density map is performed either by FOLDHUNTER [J. Mol. Biol. 308 (2001) 1033] or by a new density fitting module of MODELLER (Mod-EM). Identification of the most accurate model is based on the correlation between the model accuracy and the quality of fit into the cryoEM density map. To quantify this correlation, we created a benchmark consisting of eight proteins of different structural folds with corresponding density maps simulated at five resolutions from 5 to 15 angstroms, with three noise levels each. Each of the proteins in the set was modeled based on 300 different alignments to their remotely related templates (12-32% sequence identity), spanning the range from entirely inaccurate to essentially accurate alignments. The benchmark revealed that one of the most accurate models can usually be identified by the quality of its fit into the cryoEM density map, even for noisy maps at 15 angstroms resolution. Therefore, a cryoEM density map can be helpful in improving the accuracy of a comparative model. Moreover, a pseudo-atomic model of a component in an assembly may be built better with comparative models of the native subunit sequences than with experimentally determined structures of their homologs.

The developmental miRNA profiles of zebrafish as determined by small RNA cloning.

MicroRNAs (miRNAs) represent a family of small, regulatory, noncoding RNAs that are found in plants and animals. Here, we describe the miRNA profile of the zebrafish Danio rerio resolved in a developmental and cell-type-specific manner. The profiles were obtained from larger-scale sequencing of small RNA libraries prepared from developmentally staged zebrafish, and two adult fibroblast cell lines derived from the caudal fin (ZFL) and the liver epithelium (SJD). We identified a total of 154 distinct miRNAs expressed from 343 miRNA genes. Other experimental/computational sources support an additional 10 miRNAs encoded by 19 genes. The miRNAs can be classified into 87 distinct families. Cross-species comparison indicates that 81 families are conserved in mammals, 17 of which also have at least one member conserved in an invertebrate. Our analysis reveals that the zygotes are essentially devoid of miRNAs and that their expression begins during the blastula period with a zebrafish-specific family of miRNAs encoded by closely spaced multicopy genes. Computational predictions of zebrafish miRNA targets are provided that take into account the depth of evolutionary conservation. Besides miRNAs, we identified a prominent class of repeat-associated small interfering RNAs (rasiRNAs).

MicroRNA targets in Drosophila.

BACKGROUND: The recent discoveries of microRNA (miRNA) genes and characterization of the first few target genes regulated by miRNAs in Caenorhabditis elegans and Drosophila melanogaster have set the stage for elucidation of a novel network of regulatory control. We present a computational method for whole-genome prediction of miRNA target genes. The method is validated using known examples. For each miRNA, target genes are selected on the basis of three properties: sequence complementarity using a position-weighted local alignment algorithm, free energies of RNA-RNA duplexes, and conservation of target sites in related genomes. Application to the D. melanogaster, Drosophila pseudoobscura and Anopheles gambiae genomes identifies several hundred target genes potentially regulated by one or more known miRNAs. RESULTS: These potential targets are rich in genes that are expressed at specific developmental stages and that are involved in cell fate specification, morphogenesis and the coordination of developmental processes, as well as genes that are active in the mature nervous system. High-ranking target genes are enriched in transcription factors two-fold and include genes already known to be under translational regulation. Our results reaffirm the thesis that miRNAs have an important role in establishing the complex spatial and temporal patterns of gene activity necessary for the orderly progression of development and suggest additional roles in the function of the mature organism. In addition the results point the way to directed experiments to determine miRNA functions. CONCLUSIONS: The emerging combinatorics of miRNA target sites in the 3' untranslated regions of messenger RNAs are reminiscent of transcriptional regulation in promoter regions of DNA, with both one-to-many and many-to-one relationships between regulator and target. Typically, more than one miRNA regulates one message, indicative of cooperative translational control. Conversely, one miRNA may have several target genes, reflecting target multiplicity. As a guide to focused experiments, we provide detailed online information about likely target genes and binding sites in their untranslated regions, organized by miRNA or by gene and ranked by likelihood of match. The target prediction algorithm is freely available and can be applied to whole genome sequences using identified miRNA sequences.