The extent of sequence complementarity correlates with the potency of cellular miRNA-mediated restriction of HIV-1.

MicroRNAs (miRNAs) are 22-nt non-coding RNAs involved in the regulation of cellular gene expression and potential cellular defense against viral infection. Using in silico analyses, we predicted target sites for 22 human miRNAs in the HIV genome. Transfection experiments using synthetic miRNAs showed that five of these miRNAs capably decreased HIV replication. Using one of these five miRNAs, human miR-326 as an example, we demonstrated that the degree of complementarity between the predicted viral sequence and cellular miR-326 correlates, in a Dicer-dependent manner, with the potency of miRNA-mediated restriction of viral replication. Antagomirs to miR-326 that knocked down this cell endogenous miRNA increased HIV-1 replication in cells, suggesting that miR-326 is physiologically functional in moderating HIV-1 replication in human cells.

Two-tiered approach identifies a network of cancer and liver disease-related genes regulated by miR-122.

MicroRNAs function as important regulators of gene expression and are commonly linked to development, differentiation, and diseases such as cancer. To better understand their roles in various biological processes, identification of genes targeted by microRNAs is necessary. Although prediction tools have significantly helped with this task, experimental approaches are ultimately required for extensive target search and validation. We employed two independent yet complementary high throughput approaches to map a large set of mRNAs regulated by miR-122, a liver-specific microRNA implicated in regulation of fatty acid and cholesterol metabolism, hepatitis C infection, and hepatocellular carcinoma. The combination of luciferase reporter-based screening and shotgun proteomics resulted in the identification of 260 proteins significantly down-regulated in response to miR-122 in at least one method, 113 of which contain predicted miR-122 target sites. These proteins are enriched for functions associated with the cell cycle, differentiation, proliferation, and apoptosis. Among these miR-122-sensitive proteins, we identified a large group with strong connections to liver metabolism, diseases, and hepatocellular carcinoma. Additional analyses, including examination of consensus binding motifs for both miR-122 and target sequences, provide further insight into miR-122 function.

Kaposi's sarcoma-associated herpesvirus ORF57 promotes escape of viral and human interleukin-6 from microRNA-mediated suppression.

Kaposi's sarcoma-associated herpesvirus (KSHV) lytic infection increases the expression of viral and human interleukin-6 (vIL-6 and hIL-6, respectively), an important factor for cell growth and pathogenesis. Here, we report genome-wide analysis of viral RNA targets of KSHV ORF57 by a novel UV-cross-linking and immunoprecipitation (CLIP) assay. We identified 11 viral transcripts as putative ORF57 targets and demonstrate that vIL-6 mRNA is an authentic target of ORF57. Disrupting the ORF57 gene in the KSHV genome leads to inefficient expression of vIL-6. With transient transfection, the expression of vIL-6 could be enhanced greatly in the presence of ORF57 in a dose-dependent manner. We found that the open reading frame (ORF) region of vIL-6 RNA contains an MRE (MTA [ORF57]-responsive element) composed of two motifs, MRE-A and MRE-B, and binding of ORF57 to these two motifs stabilizes vIL-6 RNA and promotes vIL-6 translation. We demonstrate that vIL-6 MRE-B bears an miR-1293 binding site and that, mechanistically, ORF57 competes with miR-1293 for the same binding site to interact with vIL-6 RNA, thereby preventing vIL-6 RNA from association with the miR-1293-specified RNA-induced silencing complex (RISC). Consistent with this, ORF57 also interacts with an miR-608 binding site in the hIL-6 ORF and prevents miR-608 repression of hIL-6. Collectively, our results identify a novel function of ORF57 in being responsible for stabilization of viral and human IL-6 RNAs and the corresponding enhancement of RNA translation. In addition, our data provide the first evidence that a tumor virus may use a viral protein to interfere with microRNA (miRNA)-mediated repression of an miRNA target to induce cell proliferation and tumorigenesis during virus infection.

Sequence signatures and mRNA concentration can explain two-thirds of protein abundance variation in a human cell line.

Transcription, mRNA decay, translation and protein degradation are essential processes during eukaryotic gene expression, but their relative global contributions to steady-state protein concentrations in multi-cellular eukaryotes are largely unknown. Using measurements of absolute protein and mRNA abundances in cellular lysate from the human Daoy medulloblastoma cell line, we quantitatively evaluate the impact of mRNA concentration and sequence features implicated in translation and protein degradation on protein expression. Sequence features related to translation and protein degradation have an impact similar to that of mRNA abundance, and their combined contribution explains two-thirds of protein abundance variation. mRNA sequence lengths, amino-acid properties, upstream open reading frames and secondary structures in the 5' untranslated region (UTR) were the strongest individual correlates of protein concentrations. In a combined model, characteristics of the coding region and the 3'UTR explained a larger proportion of protein abundance variation than characteristics of the 5'UTR. The absolute protein and mRNA concentration measurements for >1000 human genes described here represent one of the largest datasets currently available, and reveal both general trends and specific examples of post-transcriptional regulation.

Pyrosequencing of small non-coding RNAs in HIV-1 infected cells: evidence for the processing of a viral-cellular double-stranded RNA hybrid.

Small non-coding RNAs of 18-25 nt in length can regulate gene expression through the RNA interference (RNAi) pathway. To characterize small RNAs in HIV-1-infected cells, we performed linker-ligated cloning followed by high-throughput pyrosequencing. Here, we report the composition of small RNAs in HIV-1 productively infected MT4 T-cells. We identified several HIV-1 small RNA clones and a highly abundant small 18-nt RNA that is antisense to the HIV-1 primer-binding site (PBS). This 18-nt RNA apparently originated from the dsRNA hybrid formed by the HIV-1 PBS and the 3' end of the human cellular tRNAlys3. It was found to associate with the Ago2 protein, suggesting its possible function in the cellular RNAi machinery for targeting HIV-1.

The A-rich RNA sequences of HIV-1 pol are important for the synthesis of viral cDNA.

The bias of A-rich codons in HIV-1 pol is thought to be a record of hypermutations in viral genomes that lack biological functions. Bioinformatic analysis predicted that A-rich sequences are generally associated with minimal local RNA structures. Using codon modifications to reduce the amount of A-rich sequences within HIV-1 genomes, we have reduced the flexibility of RNA sequences in pol to analyze the functional significance of these A-rich 'structurally poor' RNA elements in HIV-1 pol. Our data showed that codon modification of HIV-1 sequences led to a suppression of virus infectivity by 5-100-fold, and this defect does not correlate with, viral entry, viral protein expression levels, viral protein profiles or virion packaging of genomic RNA. Codon modification of HIV-1 pol correlated with an enhanced dimer stability of the viral RNA genome, which was associated with a reduction of viral cDNA synthesis both during HIV-1 infection and in a cell free reverse transcription assay. Our data provided direct evidence that the HIV-1 A-rich pol sequence is not merely an evolutionary artifact of enzyme-induced hypermutations, and that HIV-1 has adapted to rely on A-rich RNA sequences to support the synthesis of viral cDNA during reverse transcription, highlighting the utility of using 'structurally poor' RNA domains in regulating biological process.

Data mining of imperfect double-stranded RNA in 3' untranslated regions of eukaryotic mRNAs.

Recent developments in the study of RNA silencing indicate that double-stranded RNA (dsRNA) can be used in eukaryotes to block expression of a corresponding cellular gene. There is also a large class of small non-coding RNAs having potential to form a distinct, stable stem-loop in numbers of eukaryotic genomes. We had reported that a large imperfect dsRNA structure with hundreds of base-pairs (bp) in the 3' untranslated region (3' UTR) of cytotoxic ribonuclease was correlated with the translation suppression. In this study, we search for such dsRNAs in a 3' UTR database. The occurrence rate of large dsRNA in 3' UTRs ranges from 0.01% in plant to 0.30% in vertebrate mRNAs. However, small imperfect dsRNAs of ~ 30 bp are much more prevalent than large ones. The small dsRNAs are statistically very significant and uniquely well-ordered. Most of them have the conserved structural features of pre-miRNAs. Our data mining of the dsRNAs in the 3' UTR database can be used to explore RNA-based regulation of gene expression.

MicroRNAs in human immunodeficiency virus-1 infection.

Initially reported for Caenorhabditis elegans, microRNA (miRNA) has been shown to regulate gene expression in plants, flies, and mammals . Here, we report on our approaches to investigate the role of miRNAs in human immunodeficiency virus (HIV)-1 infection. Using computer-directed foldings, we first identify potential sequences in HIV-1 that putatively encode miRNAs. Subsequently, we use Northern blotting of RNAs isolated from HIV-infected cells to confirm expression of predicted miRNA sequences. Finally, we use a scanning algorithm to search 3' untranslated regions (UTRs) of human messenger RNAs (mRNAs) in the attempt to predict potential sites targeted by HIV-1 miRNAs.

RNA molecules with structure dependent functions are uniquely folded.

Cis-acting elements in post-transcriptional regulation of gene expression are often correlated with distinct local RNA secondary structure. These structures are expected to be significantly more ordered than those anticipated at random because of evolutionary constraints and intrinsic structural properties. In this study, we introduce a computing method to calculate two quantitative measures, NRd and Stscr, for estimating the uniqueness of an RNA secondary structure. NRd is a normalized score based on evaluating how different a natural RNA structure is from those predicted for its randomly shuffled variants. The lower the score NRd the more well ordered is the natural RNA structure. The statistical significance of NRd compared with that computed from structural comparisons among large numbers of randomly permuted sequences is represented by a standardized score, STSCR: We tested the method on the trans-activation response element and Rev response element of HIV-1 mRNA, internal ribosome entry sequence of hepatitis C virus, Tetrahymena thermophila rRNA intron, 100 tRNAs and 14 RNase P RNAs. Our data indicate that functional RNA structures have high Stscr, while other structures have low Stscr. We conclude that RNA functional molecules and/or cis-acting elements with structure dependent functions possess well ordered conformations and they are uniquely folded as measured by this technique.

Mapping and characterization of the minimal internal ribosome entry segment in the human c-myc mRNA 5' untranslated region.

The human c-myc proto-oncogene is transcribed from four alternative promoters generating transcripts with 5' untranslated regions of various lengths. These transcripts encode two proteins, c-Myc1 and c-Myc2, from two initiation codons, CUG and AUG, respectively. We and others have previously demonstrated that the region of c-myc transcripts between nucleotides (nt) -363 and -94 upstream from the CUG start codon contained an internal ribosome entry site leading to the cap-independent translation of c-myc open reading frames (ORFs). Here, we mapped a 50-nt sequence (-143 -94), which is sufficient to promote internal translation initiation of c-myc ORFs. Interestingly, this 50-nt element can be further dissected into two segments of 14 nt, each capable of activating internal translation initiation. We also demonstrate that this 50-nt element acts as the ribosome landing site from which the preinitiation ribosomal complex scans the mRNA until the CUG or AUG start codons.

siRNA, miRNA and HIV: promises and challenges.

Small interfering RNA (siRNA) and microRNA (miRNA) are small RNAs of 18-25 nucleotides (nt) in length that play important roles in regulating gene expression. They are incorporated into an RNA-induced silencing complex (RISC) and serve as guides for silencing their corresponding target mRNAs based on complementary base-pairing. The promise of gene silencing has led many researchers to consider siRNA as an anti-viral tool. However, in long-term settings, many viruses appear to escape from this therapeutical strategy. An example of this may be seen in the case of human immunodeficiency virus type-1 (HIV-1) which is able to evade RNA silencing by either mutating the siRNA-targeted sequence or by encoding for a partial suppressor of RNAi (RNA interference). On the other hand, because miRNA targeting does not require absolute complementarity of base-pairing, mutational escape by viruses from miRNA-specified silencing may be more difficult to achieve. In this review, we discuss stratagems used by various viruses to avoid the cells' antiviral si/mi-RNA defenses and notions of how viruses might control and regulate host cell genes by encoding viral miRNAs (vmiRNAs).

HIV-1 encoded candidate micro-RNAs and their cellular targets.

MicroRNAs (miRNAs) are small RNAs of 21-25 nucleotides that specifically regulate cellular gene expression at the post-transcriptional level. miRNAs are derived from the maturation by cellular RNases III of imperfect stem loop structures of ~ 70 nucleotides. Evidence for hundreds of miRNAs and their corresponding targets has been reported in the literature for plants, insects, invertebrate animals, and mammals. While not all of these miRNA/target pairs have been functionally verified, some clearly serve roles in regulating normal development and physiology. Recently, it has been queried whether the genome of human viruses like their cellular counterpart also encode miRNA. To date, there has been only one report pertaining to this question. The Epstein-Barr virus (EBV) has been shown to encode five miRNAs. Here, we extend the analysis of miRNA-encoding potential to the human immunodeficiency virus (HIV). Using computer-directed analyses, we found that HIV putatively encodes five candidate pre-miRNAs. We then matched deduced mature miRNA sequences from these 5 pre-miRNAs against a database of 3' untranslated sequences (UTR) from the human genome. These searches revealed a large number of cellular transcripts that could potentially be targeted by these viral miRNA (vmiRNA) sequences. We propose that HIV has evolved to use vmiRNAs as a means to regulate cellular milieu for its benefit.

Data mining of functional RNA structures in genomic sequences.

The normal functions of genomes depend on the precise expression of messenger RNAs and noncoding RNAs (ncRNAs) such as transfer RNAs and microRNAs in eukaryotes. These ncRNAs and functional RNA structures (FRSs) act as regulators or response elements for cellular factors and participate in transcription, posttranscriptional processing, and translation. Knowledge discovery of these FRSs in huge DNA/RNA sequence databases is a very important step to reach our goal of going from genomic sequence data to biological knowledge for understanding RNA-based regulation. Analyses of a large number of FRSs have indicated that the FRS can be well characterized by some quantitative measures such as significance and well-ordered scores of the local segment. Various data mining tools have been developed and successfully applied to FRS discovery in genomic sequence databases. Here, we summarize our efforts in the computational discovery of structured features of ncRNAs and FRSs within complex genomes by EDscan and SigED.

Stability of a long noncoding viral RNA depends on a 9-nt core element at the RNA 5' end to interact with viral ORF57 and cellular PABPC1.

Kaposi sarcoma-associated herpesvirus (KSHV) ORF57, also known as Mta (mRNA transcript accumulation), enhances viral intron-less transcript accumulation and promotes splicing of intron-containing viral RNA transcripts. In this study, we identified KSHV PAN, a long non-coding polyadenylated nuclear RNA as a main target of ORF57 by a genome-wide CLIP (cross-linking and immunoprecipitation) approach. KSHV genome lacking ORF57 expresses only a minimal amount of PAN. In cotransfection experiments, ORF57 alone increased PAN expression by 20-30-fold when compared to vector control. This accumulation function of ORF57 was dependent on a structured RNA element in the 5' PAN, named MRE (Mta responsive element), but not much so on an ENE (expression and nuclear retention element) in the 3' PAN previously reported by other studies. We showed that the major function of the 5' PAN MRE is increasing the RNA half-life of PAN in the presence of ORF57. Further mutational analyses revealed a core motif consisting of 9 nucleotides in the MRE-II , which is responsible for ORF57 interaction and function. The 9-nt core in the MRE-II also binds cellular PABPC1, but not the E1B-AP5 which binds another region of the MRE-II. In addition, we found that PAN RNA is partially exportable in the presence of ORF57. Together, our data provide compelling evidence as to how ORF57 functions to accumulate a non-coding viral RNA in the course of virus lytic infection.

Genomic analyses of musashi1 downstream targets show a strong association with cancer-related processes.

Musashi1 (Msi1) is a highly conserved RNA-binding protein with pivotal functions in stem cell maintenance, nervous system development, and tumorigenesis. Despite its importance, only three direct mRNA targets have been characterized so far: m-numb, CDKN1A, and c-mos. Msi1 has been shown to affect their translation by binding to short elements located in the 3'-untranslated region. To better understand Msi1 functions, we initially performed an RIP-Chip analysis in HEK293T cells; this method consists of isolation of specific RNA-protein complexes followed by identification of the RNA component via microarrays. A group of 64 mRNAs was found to be enriched in the Msi1-associated population compared with controls. These genes belong to two main functional categories pertinent to tumorigenesis: 1) cell cycle, cell proliferation, cell differentiation, and apoptosis and 2) protein modification (including ubiquitination and ubiquitin cycle). To corroborate our findings, we examined the impact of Msi1 expression on both mRNA (transcriptomic) and protein (proteomic) expression levels. Genes whose mRNA levels were affected by Msi1 expression have a Gene Ontology distribution similar to RIP-Chip results, reinforcing Msi1 participation in cancer-related processes. The proteomics study revealed that Msi1 can have either positive or negative effects on gene expression of its direct targets. In summary, our results indicate that Msi1 affects a network of genes and could function as a master regulator during development and tumor formation.

Aberrant expression of oncogenic and tumor-suppressive microRNAs in cervical cancer is required for cancer cell growth.

MicroRNAs (miRNAs) play important roles in cancer development. By cloning and sequencing of a HPV16(+) CaSki cell small RNA library, we isolated 174 miRNAs (including the novel miR-193c) which could be grouped into 46 different miRNA species, with miR-21, miR-24, miR-27a, and miR-205 being most abundant. We chose for further study 10 miRNAs according to their cloning frequency and associated their levels in 10 cervical cancer- or cervical intraepithelial neoplasia-derived cell lines. No correlation was observed between their expression with the presence or absence of an integrated or episomal HPV genome. All cell lines examined contained no detectable miR-143 and miR-145. HPV-infected cell lines expressed a different set of miRNAs when grown in organotypic raft cultured as compared to monolayer cell culture, including expression of miR-143 and miR-145. This suggests a correlation between miRNA expression and tissue differentiation. Using miRNA array analyses for age-matched normal cervix and cervical cancer tissues, in combination with northern blot verification, we identified significantly deregulated miRNAs in cervical cancer tissues, with miR-126, miR-143, and miR-145 downregulation and miR-15b, miR-16, miR-146a, and miR-155 upregulation. Functional studies showed that both miR-143 and miR-145 are suppressive to cell growth. When introduced into cell lines, miR-146a was found to promote cell proliferation. Collectively, our data indicate that downregulation of miR-143 and miR-145 and upregulation of miR-146a play a role in cervical carcinogenesis.

RNA-binding protein HuR interacts with thrombomodulin 5'untranslated region and represses internal ribosome entry site-mediated translation under IL-1 beta treatment.

Reduction in host-activated protein C levels and resultant microvascular thrombosis highlight the important functional role of protein C anticoagulant system in the pathogenesis of sepsis and septic shock. Thrombomodulin (TM) is a critical factor to activate protein C in mediating the anticoagulation and anti-inflammation effects. However, TM protein content is decreased in inflammation and sepsis, and the mechanism is still not well defined. In this report, we identified that the TM 5' untranslated region (UTR) bearing the internal ribosome entry site (IRES) element controls TM protein expression. Using RNA probe pulldown assay, HuR was demonstrated to interact with the TM 5'UTR. Overexpression of HuR protein inhibited the activity of TM IRES, whereas on the other hand, reducing the HuR protein level reversed this effect. When cells were treated with IL-1beta, the IRES activity was suppressed and accompanied by an increased interaction between HuR and TM 5'UTR. In the animal model of sepsis, we found the TM protein expression level to be decreased while concurrently observing the increased interaction between HuR and TM mRNA in liver tissue. In summary, HuR plays an important role in suppression of TM protein synthesis in IL-1beta treatment and sepsis.

Evidence that HIV-1 encodes an siRNA and a suppressor of RNA silencing.

In plants and invertebrate animals, RNA silencing is a form of nucleic acid-based adaptive immunity. By contrast, jawed vertebrates have evolved complex protein-based adaptive immunity. Although short interfering RNAs (siRNAs) have been used as artificial tools to silence viral infection in human cells, it remains unknown whether mammalian viruses naturally elicit such immunity in vertebral cells. Here, we report the evidence that HIV-1 encodes viral siRNA precursors in its genome and that natural HIV-1 infection provokes nucleic acid-based immunity in human cells. To combat this cellular defense, HIV-1 has evolved in its Tat protein a suppressor of RNA silencing (SRS) function. Tat abrogates the cell's RNA-silencing defense by subverting the ability of Dicer to process precursor double-stranded RNAs into siRNAs.

Statistical inference for well-ordered structures in nucleotide sequences.

Distinct, local structures are frequently correlated with functional RNA elements involved in post-transcriptional regulation of gene expression. Discovery of microRNAs (miRNAs) suggests that there are a large class of small non-coding RNAs in eukaryotic genomes. These miRNAs have the potential to form distinct fold-back stem-loop structures. The prediction of those well-ordered folding sequences (WFS) in genomic sequences is very helpful for our understanding of RNA-based gene regulation and the determination of local RNA elements with structure-dependent functions. In this study, we describe a novel method for discovering the local WFS in a nucleotide sequence by Monte Carlo simulation and RNA folding. In the approach the quality of a local WFS is assessed by the energy difference (E(diff)) between the optimal structure folded in the local segment and its corresponding optimal, restrained structure where all the previous base pairings formed in the optimal structure are prohibited. Distinct WFS can be discovered by scanning successive segments along a sequence for evaluating the difference between E(diff) of the natural sequence and those computed from randomly shuffled sequences. Our results indicate that the statistically significant WFS detected in the genomic sequences of Caenorhabditis elegans (C.elegans) F49E12, T07C5, T07D1, T10H9, Y56A3A and Y71G12B are coincident with known fold-back stem-loops found in miRNA precursors. The potential and implications of our method in searching for miRNAs in genomes is discussed.