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.

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.

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.

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.

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.

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.

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.

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.

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.

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).

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.

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.

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.

The parathyroid hormone mRNA 3'-untranslated region AU-rich element is an unstructured functional element.

Parathyroid hormone (PTH) gene expression is regulated post-transcriptionally by hypocalcemia and hypophosphatemia. This regulation is dependent upon binding of protective trans-acting factors to a specific element in the PTH mRNA 3'-untranslated region (UTR). We have previously demonstrated that a 63-nucleotide (nt) AU-rich PTH mRNA element is sufficient to confer regulation of RNA stability by calcium and phosphate in an in vitro degradation assay (IVDA). The 63-nt element consists of a core 26-nt minimal binding sequence and flanking regions. We have now studied the functionality of this element in HEK293 cells using reporter genes and showed that it destabilizes mRNAs for green fluorescent protein (GFP) and growth hormone, similar to its effect in the IVDA. To understand how the cis-element functions as an instability element, we have analyzed its structure by RNase H, primer extension, and computer modeling. The results indicate that the PTH mRNA 3'-UTR and in particular the region of the cis-element are dominated by significant open regions with little folded base pairing. Mutation analysis of the 26-nt core element demonstrated the importance of defined nucleotides for protein-RNA binding. In the GFP reporter system, the same mutations that prevented binding were also ineffective in destabilizing GFP mRNA in HEK293 cells. This is the first study of an AU-rich element that relates function to structure. The PTH mRNA 3'-UTR cis-acting element is an open region that utilizes the distinct sequence pattern to determine mRNA stability by its interaction with trans-acting factors.

An algorithm for detecting homologues of known structured RNAs in genomes.

Distinct RNA structures are frequently involved in a wide-range of functions in various biological mechanisms. The three dimensional RNA structures solved by X-ray crystallography and various well-established RNA phylogenetic structures indicate that functional RNAs have characteristic RNA structural motifs represented by specific combinations of base pairings and conserved nucleotides in the loop region. Discovery of well-ordered RNA structures and their homologues in genome-wide searches will enhance our ability to detect the RNA structural motifs and help us to highlight their association with functional and regulatory RNA elements. We present here a novel computer algorithm, HomoStRscan, that takes a single RNA sequence with its secondary structure to search for homologous-RNAs in complete genomes. This novel algorithm completely differs from other currently used search algorithms of homologous structures or structural motifs. For an arbitrary segment (or window) given in the target sequence, that has similar size to the query sequence, HomoStRscan finds the most similar structure to the input query structure and computes the maximal similarity score (MSS) between the two structures. The homologousRNA structures are then statistically inferred from the MSS distribution computed in the target genome. The method provides a flexible, robust and fine search tool for any homologous structural RNAs.