Reduced seed region-based off-target activity with lentivirus-mediated RNAi.

Along with silencing intended target genes, transfected siRNAs regulate numerous unintended transcripts through a mechanism in which the equivalent of a microRNA-like seed region in the siRNA recognizes complementary sequences in transcript 3' UTRs. Amelioration of this off-target silencing would lead to more accurate interpretation of RNA interference (RNAi) experiments and thus greatly enhance their value. We tested whether lentivirus-mediated delivery of shRNA is prone to the sequence-based off-target activity prevalent in siRNA experiments. We compared target gene silencing and overall impact on global gene expression caused by multiple sequences delivered as both transfected siRNAs and lentivirus vector-expressed shRNAs. At equivalent levels of target gene silencing, signatures induced by shRNAs were significantly smaller than those induced by cognate siRNAs and arose less frequently from seed region activity. Interestingly, the low level of seed region-based off-target activity exhibited by shRNAs resulted in down-regulation of transcripts that were largely distinct from those regulated by siRNAs. On the basis of these observations, we recommend lentivirus-mediated RNAi for pathway profiling experiments that measure whole genome transcriptional readouts as well as for large-scale screens when resources for extensive follow up are limited.

Microarray analysis shows that some microRNAs downregulate large numbers of target mRNAs.

MicroRNAs (miRNAs) are a class of noncoding RNAs that post-transcriptionally regulate gene expression in plants and animals. To investigate the influence of miRNAs on transcript levels, we transfected miRNAs into human cells and used microarrays to examine changes in the messenger RNA profile. Here we show that delivering miR-124 causes the expression profile to shift towards that of brain, the organ in which miR-124 is preferentially expressed, whereas delivering miR-1 shifts the profile towards that of muscle, where miR-1 is preferentially expressed. In each case, about 100 messages were downregulated after 12 h. The 3' untranslated regions of these messages had a significant propensity to pair to the 5' region of the miRNA, as expected if many of these messages are the direct targets of the miRNAs. Our results suggest that metazoan miRNAs can reduce the levels of many of their target transcripts, not just the amount of protein deriving from these transcripts. Moreover, miR-1 and miR-124, and presumably other tissue-specific miRNAs, seem to downregulate a far greater number of targets than previously appreciated, thereby helping to define tissue-specific gene expression in humans.

Transcripts targeted by the microRNA-16 family cooperatively regulate cell cycle progression.

microRNAs (miRNAs) are abundant, approximately 21-nucleotide, noncoding regulatory RNAs. Each miRNA may regulate hundreds of mRNA targets, but the identities of these targets and the processes they regulate are poorly understood. Here we have explored the use of microarray profiling and functional screening to identify targets and biological processes triggered by the transfection of human cells with miRNAs. We demonstrate that a family of miRNAs sharing sequence identity with miRNA-16 (miR-16) negatively regulates cellular growth and cell cycle progression. miR-16-down-regulated transcripts were enriched with genes whose silencing by small interfering RNAs causes an accumulation of cells in G(0)/G(1). Simultaneous silencing of these genes was more effective at blocking cell cycle progression than disruption of the individual genes. Thus, miR-16 coordinately regulates targets that may act in concert to control cell cycle progression.

Designing siRNA that distinguish between genes that differ by a single nucleotide.

Small interfering RNAs (siRNAs), the guides that direct RNA interference (RNAi), provide a powerful tool to reduce the expression of a single gene in human cells. Ideally, dominant, gain-of-function human diseases could be treated using siRNAs that specifically silence the mutant disease allele, while leaving expression of the wild-type allele unperturbed. Previous reports suggest that siRNAs can be designed with single nucleotide specificity, but no rational basis for the design of siRNAs with single nucleotide discrimination has been proposed. We systematically identified siRNAs that discriminate between the wild-type and mutant alleles of two disease genes: the human Cu, Zn superoxide dismutase (SOD1) gene, which contributes to the progression of hereditary amyotrophic lateral sclerosis through the gain of a toxic property, and the huntingtin (HTT) gene, which causes Huntington disease when its CAG-repeat region expands beyond approximately 35 repeats. Using cell-free RNAi reactions in Drosophila embryo lysate and reporter assays and microarray analysis of off-target effects in cultured human cells, we identified positions within an siRNA that are most sensitive to mismatches. We also show that purine:purine mismatches imbue an siRNA with greater discriminatory power than other types of base mismatches. siRNAs in which either a G:U wobble or a mismatch is located in the "seed" sequence, the specialized siRNA guide region responsible for target binding, displayed lower levels of selectivity than those in which the mismatch was located 3' to the seed; this region of an siRNA is critical for target cleavage but not siRNA binding. Our data suggest that siRNAs can be designed to discriminate between the wild-type and mutant alleles of many genes that differ by just a single nucleotide.

Expression profiling reveals off-target gene regulation by RNAi.

RNA interference is thought to require near-identity between the small interfering RNA (siRNA) and its cognate mRNA. Here, we used gene expression profiling to characterize the specificity of gene silencing by siRNAs in cultured human cells. Transcript profiles revealed siRNA-specific rather than target-specific signatures, including direct silencing of nontargeted genes containing as few as eleven contiguous nucleotides of identity to the siRNA. These results demonstrate that siRNAs may cross-react with targets of limited sequence similarity.

A distinct Smoothened mutation causes severe cerebellar developmental defects and medulloblastoma in a novel transgenic mouse model.

Deregulated developmental processes in the cerebellum cause medulloblastoma, the most common pediatric brain malignancy. About 25 to 30% of cases are caused by mutations increasing the activity of the Sonic hedgehog (Shh) pathway, a critical mitogen in cerebellar development. The proto-oncogene Smoothened (Smo) is a key transducer of the Shh pathway. Activating mutations in Smo that lead to constitutive activity of the Shh pathway have been identified in human medulloblastoma. To understand the developmental and oncogenic effects of two closely positioned point mutations in Smo, we characterized NeuroD2-SmoA2 mice and compared them to NeuroD2-SmoA1 mice. While both SmoA1 and SmoA2 transgenes cause medulloblastoma with similar frequencies and timing, SmoA2 mice have severe aberrations in cerebellar development, whereas SmoA1 mice are largely normal during development. Intriguingly, neurologic function, as measured by specific tests, is normal in the SmoA2 mice despite extensive cerebellar dysplasia. We demonstrate how two nearly contiguous point mutations in the same domain of the encoded Smo protein can produce striking phenotypic differences in cerebellar development and organization in mice.

p53-responsive miR-194 inhibits thrombospondin-1 and promotes angiogenesis in colon cancers.

Thrombospondin-1 (TSP-1) is an endogenous inhibitor of angiogenesis encoded by the THBS1 gene, whose promoter is activated by p53. In advanced colorectal cancers (CRC), its expression is sustained or even slightly increased despite frequent loss of p53. Here, we determined that in HCT116 CRC cells, p53 activates the THBS1 primary transcript, but fails to boost THBS1 mRNA or protein levels, implying posttranscriptional regulation by microRNAs (miRNA). In a global miRNA gain-of-function screen done in the Dicer-deficient HCT116 variant, several miRNAs negatively regulated THBS1 mRNA and protein levels, one of them being miR-194. Notably, in agreement with published data, p53 upregulated miR-194 expression in THBS1 retrovirus-transduced HCT116 cells, leading to decreased TSP-1 levels. This negative effect was mediated by a single miR-194 complementary site in the THBS1 3'-untranslated region, and its elimination resulted in TSP-1 reactivation, impaired angiogenesis in Matrigel plugs, and reduced growth of HCT116 xenografts. Conversely, transient overexpression of miR-194 in HCT116/THBS1 cells boosted Matrigel angiogenesis, and its stable overexpression in Ras-induced murine colon carcinomas increased microvascular densities and vessel sizes. Although the overall contribution of miR-194 to neoplastic growth is context dependent, p53-induced activation of this GI tract-specific miRNA during ischemia could promote angiogenesis and facilitate tissue repair.

Transcriptome-wide studies of merkel cell carcinoma and validation of intratumoral CD8+ lymphocyte invasion as an independent predictor of survival.

PURPOSE: Merkel cell carcinoma (MCC) is a polyomavirus-associated skin cancer that is frequently lethal and lacks established prognostic biomarkers. This study sought to identify biomarkers that improve prognostic accuracy and provide insight into MCC biology. PATIENTS AND METHODS: Gene expression profiles of 35 MCC tumors were clustered based on prognosis. The cluster of genes overexpressed in good-prognosis tumors was tested for biologic process enrichment. Relevant mRNA expression differences were confirmed by quantitative polymerase chain reaction and immunohistochemistry. An independent set of 146 nonoverlapping MCC tumors (median follow-up, 25 months among 116 living patients) was employed for biomarker validation. Univariate and multivariate Cox regression analyses were performed. RESULTS: Immune response gene signatures were prominent in patients with good prognoses. In particular, genes associated with cytotoxic CD8+ lymphocytes were overexpressed in tumors from patients with favorable prognoses. In the independent validation set, cases with robust intratumoral CD8+ lymphocyte infiltration had improved outcomes (100% MCC-specific survival, n = 26) compared with instances characterized by sparse infiltration (60% survival, n = 120). Only stage and intratumoral CD8 infiltration (but not age, sex, or CD8+ lymphocytes localized to the tumor-stroma interface) were significant in both univariate and multivariate Cox regression analyses. Notably, traditional histologic identification of tumor-infiltrating lymphocytes was not a significant independent predictor of survival. CONCLUSION: Intratumoral CD8+ lymphocyte infiltration can be readily assessed on paraffin-embedded tissue, is independently associated with improved MCC-specific survival, and therefore, may provide prognostic information that enhances established MCC staging protocols.

The myc-miR-17~92 axis blunts TGF{beta} signaling and production of multiple TGF{beta}-dependent antiangiogenic factors.

c-Myc stimulates angiogenesis in tumors through mechanisms that remain incompletely understood. Recent work indicates that c-Myc upregulates the miR-17∼92 microRNA cluster and downregulates the angiogenesis inhibitor thrombospondin-1, along with other members of the thrombospondin type 1 repeat superfamily. Here, we show that downregulation of the thrombospondin type 1 repeat protein clusterin in cells overexpressing c-Myc and miR-17∼92 promotes angiogenesis and tumor growth. However, clusterin downregulation by miR-17∼92 is indirect. It occurs as a result of reduced transforming growth factor-ß (TGFß) signaling caused by targeting of several regulatory components in this signaling pathway. Specifically, miR-17-5p and miR-20 reduce the expression of the type II TGFß receptor and miR-18 limits the expression of Smad4. Supporting these results, in human cancer cell lines, levels of the miR-17∼92 primary transcript MIR17HG negatively correlate with those of many TGFß-induced genes that are not direct targets of miR-17∼92 (e.g., clusterin and angiopoietin-like 4). Furthermore, enforced expression of miR-17∼92 in MIR17HG(low) cell lines (e.g., glioblastoma) results in impaired gene activation by TGFß. Together, our results define a pathway in which c-Myc activation of miR-17∼92 attenuates the TGFß signaling pathway to shut down clusterin expression, thereby stimulating angiogenesis and tumor cell growth.

MicroRNA miR-210 modulates cellular response to hypoxia through the MYC antagonist MNT.

The hypoxia-inducible factor (HIF) pathway is essential for cell survival under low oxygen and plays an important role in tumor cell homeostasis. We investigated the function of miR-210, the most prominent microRNA upregulated by hypoxia and a direct transcriptional target of HIFs. miR-210 expression was elevated in multiple cancer types and correlated with metastasis of breast and melanoma tumors. miR-210 overexpression in cancer cell lines bypassed hypoxia-induced cell cycle arrest and partially reversed the hypoxic gene expression signature. We identified MNT, a known MYC antagonist, as a miR-210 target. MNT mRNA contains multiple miR-210 binding sites in the 3' UTR and its knockdown phenocopied miR-210 overexpression. Furthermore, loss of MYC abolished miR-210-mediated override of hypoxia-induced cell cycle arrest. Comparison of miR-210 and MYC overexpression with MNT knockdown signatures also indicated that miR-210 triggered a "MYC-like" transcriptional response. Thus, miR-210 influences the hypoxia response in tumor cells through targeting a key transcriptional repressor of the MYC-MAX network.

Coordinated regulation of cell cycle transcripts by p53-Inducible microRNAs, miR-192 and miR-215.

Cell cycle arrest in response to DNA damage is an important antitumorigenic mechanism. MicroRNAs (miRNAs) were recently shown to play key regulatory roles in cell cycle progression. For example, miR-34a is induced in response to p53 activation and mediates G(1) arrest by down-regulating multiple cell cycle-related transcripts. Here we show that genotoxic stress promotes the p53-dependent up-regulation of the homologous miRNAs miR-192 and miR-215. Like miR-34a, activation of miR-192/215 induces cell cycle arrest, suggesting that multiple miRNA families operate in the p53 network. Furthermore, we define a downstream gene expression signature for miR-192/215 expression, which includes a number of transcripts that regulate G(1) and G(2) checkpoints. Of these transcripts, 18 transcripts are direct targets of miR-192/215, and the observed cell cycle arrest likely results from a cooperative effect among the modulations of these genes by the miRNAs. Our results showing a role for miR-192/215 in cell proliferation combined with recent observations that these miRNAs are underexpressed in primary cancers support the idea that miR-192 and miR-215 function as tumor suppressors.

c-Myb oncoprotein is an essential target of the dleu2 tumor suppressor microRNA cluster.

The dleu2 tumor suppressor locus encodes two microRNAs, miR-15a and miR-16, which are thought to play an important role in B-cell neoplasms. However, relatively little is known about proteins that regulate or are regulated by this microRNA cluster. Here we demonstrate that the Pax5 oncoprotein downregulates the dleu2 gene and at the same time boosts expression of its own heterodimeric partner c-Myb. Interestingly, c-Myb upregulation occurs primarily at a post-transcriptional level, suggesting that it might be a target for microRNAs such as miR-15a/16. Indeed, miR-15a/16 have predicted binding sites in the c-Myb 3'-UTR and through them diminish protein output in luciferase sensor assays. Moreover, forced overexpression of miR-15a/16 reduces endogenous c-Myb levels and compromises Pax5 function. Conversely, restoration of c-Myb levels partly alleviates tumors suppressive effects of miR-15a/16, suggesting that c-Myb is a key downstream target of this microRNA cluster.

MicroRNAs in the miR-106b family regulate p21/CDKN1A and promote cell cycle progression.

microRNAs in the miR-106b family are overexpressed in multiple tumor types and are correlated with the expression of genes that regulate the cell cycle. Consistent with these observations, miR-106b family gain of function promotes cell cycle progression, whereas loss of function reverses this phenotype. Microarray profiling uncovers multiple targets of the family, including the cyclin-dependent kinase inhibitor p21/CDKN1A. We show that p21 is a direct target of miR-106b and that its silencing plays a key role in miR-106b-induced cell cycle phenotypes. We also show that miR-106b overrides a doxorubicin-induced DNA damage checkpoint. Thus, miR-106b family members contribute to tumor cell proliferation in part by regulating cell cycle progression and by modulating checkpoint functions.

Genome-wide resources of endoribonuclease-prepared short interfering RNAs for specific loss-of-function studies.

RNA interference (RNAi) has become an important technique for loss-of-gene-function studies in mammalian cells. To achieve reliable results in an RNAi experiment, efficient and specific silencing triggers are required. Here we present genome-wide data sets for the production of endoribonuclease-prepared short interfering RNAs (esiRNAs) for human, mouse and rat. We used an algorithm to predict the optimal region for esiRNA synthesis for every protein-coding gene of these three species. We created a database, RiDDLE, for retrieval of target sequences and primer information. To test this in silico resource experimentally, we generated 16,242 esiRNAs that can be used for RNAi screening in human cells. Comparative analyses with chemically synthesized siRNAs demonstrated a high silencing efficacy of esiRNAs and a 12-fold reduction of downregulated off-target transcripts as detected by microarray analysis. Hence, the presented esiRNA libraries offer an efficient, cost-effective and specific alternative to presently available mammalian RNAi resources.

Widespread siRNA "off-target" transcript silencing mediated by seed region sequence complementarity.

Transfected siRNAs and miRNAs regulate numerous transcripts that have only limited complementarity to the active strand of the RNA duplex. This process reflects natural target regulation by miRNAs, but is an unintended ("off-target") consequence of siRNA-mediated silencing. Here we demonstrate that this unintended off-target silencing is widespread, and occurs in a manner reminiscent of target silencing by miRNAs. A high proportion of unintended transcripts silenced by siRNAs showed 3' UTR sequence complementarity to the seed region of the siRNA. Base mismatches within the siRNA seed region reduced the set of original off-target transcripts but generated new sets of silenced transcripts with sequence complementarity to the mismatched seed sequence. An inducible shRNA silenced a subset of transcripts that were silenced by an siRNA of the same sequence, demonstrating that unintended silencing is sequence mediated and is independent of delivery method. In all cases, off-target transcript silencing was accompanied by loss of the corresponding protein and occurred with dependence on siRNA concentration similar to that of silencing of the target transcript. Thus, short stretches of sequence complementarity to the siRNA or shRNA seed region are key to the silencing of unintended transcripts.

Position-specific chemical modification of siRNAs reduces "off-target" transcript silencing.

Transfected siRNAs regulate numerous transcripts sharing limited complementarity to the RNA duplex. This unintended ("off-target") silencing can hinder the use of RNAi to define gene function. Here we describe position-specific, sequence-independent chemical modifications that reduced silencing of partially complementary transcripts by all siRNAs tested. Silencing of perfectly matched targets was unaffected by these modifications. The chemical modification also reduced off-target phenotypes in growth inhibition studies. Key to the modification was 2'-O-methyl ribosyl substitution at position 2 in the guide strand, which reduced silencing of most off-target transcripts with complementarity to the seed region of the siRNA guide strand. The sharp position dependence of 2'-O-methyl ribosyl modification contrasts with the broader position dependence of base-pair substitutions within the seed region, suggesting a role for position 2 of the guide strand distinct from its effects on pairing to target transcripts.

T lymphocyte activation gene identification by coregulated expression on DNA microarrays.

High-capacity methods for assessing gene function have become increasingly important because of the increasing number of newly identified genes emerging from large-scale genome sequencing and cDNA cloning efforts. We investigated the use of DNA microarrays to identify uncharacterized genes specifically involved in human T cell activation. Activation of human peripheral blood T lymphocytes induced significant changes in hundreds of transcripts, but most of these were not unique to T cell activation. Variation of experimental parameters and analysis techniques allowed better enrichment for gene expression changes unique to T cell activation. Best results were achieved by identification of genes that were most highly coregulated with the T-cell-specific transcript interleukin 2 (IL2) in a "compendium" of experiments involving both T cells and other cell types. Among the genes most highly coregulated with IL2 were many genes known to function during T cell activation, together with ESTs of unknown function. Four of these ESTs were extended to novel full-length clones encoding T-cell-regulated proteins with predicted functions in GTP metabolism, cell organization, and signal transduction.