MicroRNAs targeting mutant K-ras by electrotransfer inhibit human colorectal adenocarcinoma cell growth in vitro and in vivo

Mutations of K-ras have been found in 30-60% of colorectal carcinomas and are believed to be associated with tumor initiation, tumor progression and metastasis formation. Therefore, silencing of mutant K-ras expression has become an attractive therapeutic strategy for colorectal cancer treatment. The aim of our study was to investigate the effect of microRNA (miRNA) molecules directed against K-ras (miRNA-K-ras) on K-ras expression level and the growth of colorectal carcinoma cell line LoVo in vitro and in vivo. In addition, we evaluated electroporation as a gene delivery method for transfection of LoVo cells and tumors with plasmid DNA encoding miRNA-K-ras (pmiRNA-K-ras). Results of our study indicated that miRNAs targeting K-ras efficiently reduced K-ras expression and cell survival after in vitro electrotransfection of LoVo cells with pmiRNA-K-ras. In vivo, electroporation has proven to be a simple and efficient delivery method for local administration of pmiRNA-K-ras molecules into LoVo tumors. This therapy shows pronounced antitumor effectiveness and has no side effects. The obtained results demonstrate that electrogene therapy with miRNA-K-ras molecules can be potential therapeutic strategy for treatment of colorectal cancers harboring K-ras mutations. © 2010 Nature Publishing Group All rights reserved.

No evidence for genetic association with the let-7 microRNA-binding site or other common KRAS variants in risk of endometriosis

STUDY QUESTION Is there a contribution of the minor allele at the KRAS single nucleotide polymorphism (SNP) rs61764370 in the let-7 microRNA-binding site to endometriosis risk? SUMMARY ANSWER We found no evidence for association between endometriosis risk and rs61764370 or any other SNPs in KRAS. WHAT IS KNOWN ALREADY The rs61764370 SNP in the 3' untranslated region of the KRAS gene is predicted to disrupt a complementary binding site (LCS6) for the let-7 microRNA, and was recently reported to be at a high frequency (31%) in 132 women of varying ancestry with endometriosis compared with frequencies in a database of population controls (up to 7.6% depending on ancestry), suggesting a strong effect of this KRAS SNP in the aetiology of endometriosis. STUDY DESIGN, SIZE AND DURATION This was a case-control study with a total of 11 206 subjects. The study was performed between February 2012 and July 2012. PARTICIPANTS/MATERIALS, SETTINGAND METHODS We first investigated a possible association between common markers in KRAS and endometriosis risk from our genome-wide association (GWA) data in 3194 surgically confirmed endometriosis cases and 7060 controls of European ancestry. Although rs61764370 was not genotyped on the GWA arrays, five SNPs typed in the study were highly correlated with this variant. The rs61764370 and two SNPs highly correlated with rs61764370 were then genotyped in 933 endometriosis cases and 952 controls using the Sequenom MassARRAY platform. MAIN RESULTS AND THE ROLE OF CHANCE There was no evidence for an association between rs61764370 and endometriosis risk P = 0.411 and odds ratio = 1.10 (95% confidence intervals: 0.88-1.36). We also found no evidence for an association between the highly correlated SNP rs17387019 and endometriosis. Their minor allele frequencies in cases and controls were of 0.087-0.091 similar to the population frequency reported previously for this variant in controls. Analyses of endometriosis cases with revised American Fertility Society stage III/IV disease also showed no evidence for an association between these SNPs and endometriosis risk. LIMITATIONS AND REASONS FOR CAUTION The GWA and genotyped data sets were not independent since individuals and cases from some families overlap. Controls in our GWA study were not screened for endometriosis. WIDER IMPLICATIONS OF THE FINDINGS The key SNP, rs61764370, was genotyped in a subset of samples. Our results do not support the suggestion that carrying the minor allele at rs61764370 contributes to a significant number of endometriosis cases and rs61764370 is, therefore, unlikely to be a useful marker of endometriosis risk. STUDY FUNDING/COMPETING INTEREST(S) The research was funded by grants from the Australian National Health and Medical Research Council and Wellcome Trust. None of the authors has competing interests for the study.

Identification of a novel FGFRL1 MicroRNA target site polymorphism for bone mineral density in meta-analyses of genome-wide association studies

MicroRNAs (miRNAs) are critical post-transcriptional regulators. Based on a previous genome-wide association (GWA) scan, we conducted a polymorphism in microRNAs' Target Sites (poly-miRTS)-centric multistage meta-analysis for lumbar spine (LS)-, total hip (HIP)-, and femoral neck (FN)-bone mineral density (BMD). In stage I, 41,102 poly-miRTSs were meta-analyzed in 7 cohorts with a genome-wide significance (GWS) α=0.05/41,102=1.22×10-6. By applying α=5×10-5 (suggestive significance), 11 poly-miRTSs were selected, with FGFRL1 rs4647940 and PRR5 rs3213550 as top signals for FN-BMD (P-value=7.67×10-6 and 1.58×10-5) in gender-combined sample. In stage II in silico replication (two cohorts), FGFRL1 rs4647940 was the only signal marginally replicated for FN-BMD (P-value=5.08×10-3) at α=0.10/11=9.09×10-3. PRR5 rs3213550 was also selected based on biological significance. In stage III de novo genotyping replication (two cohorts), FGFRL1 rs4647940 was the only signal significantly replicated for FN-BMD (P-value=7.55×10-6) at α=0.05/2=0.025 in gender-combined sample. Aggregating three stages, FGFRL1 rs4647940 was the single stage I-discovered and stages II- and III-replicated signal attaining GWS for FN-BMD (P-value=8.87×10-12). Dual-luciferase reporter assays demonstrated that FGFRL1 3' untranslated region harboring rs4647940 appears to be hsa-miR-140-5p's target site. In a zebrafish microinjection experiment, dre-miR-140-5p is shown to exert a dramatic impact on craniofacial skeleton formation. Taken together, we provided functional evidence for a novel FGFRL1 poly-miRTS rs4647940 in a previously known 4p16.3 locus, and experimental and clinical genetics studies have shown both FGFRL1 and hsa-miR-140-5p are important for bone formation. © The Author 2015. Published by Oxford University Press. All rights reserved.

Oncogenic MicroRNA-155 Down-regulates Tumor Suppressor CDC73 and Promotes Oral Squamous Cell Carcinoma Cell Proliferation IMPLICATIONS FOR CANCER THERAPEUTICS

The CDC73 gene is mutationally inactivated in hereditary and sporadic parathyroid tumors. It negatively regulates beta-catenin, cyclin D1, and c-MYC. Down-regulation of CDC73 has been reported in breast, renal, and gastric carcinomas. However, the reports regarding the role of CDC73 in oral squamous cell carcinoma (OSCC) are lacking. In this study we show that CDC73 is down-regulated in a majority of OSCC samples. We further show that oncogenic microRNA-155 (miR-155) negatively regulates CDC73 expression. Our experiments show that the dramatic up-regulation of miR-155 is an exclusive mechanism for down-regulation of CDC73 in a panel of human cell lines and a subset of OSCC patient samples in the absence of loss of heterozygosity, mutations, and promoter methylation. Ectopic expression of miR-155 in HEK293 cells dramatically reduced CDC73 levels, enhanced cell viability, and decreased apoptosis. Conversely, the delivery of a miR-155 antagonist (antagomir-155) to KB cells overexpressing miR-155 resulted in increased CDC73 levels, decreased cell viability, increased apoptosis, and marked regression of xenografts in nude mice. Cotransfection of miR-155 with CDC73 in HEK293 cells abrogated its pro-oncogenic effect. Reduced cell proliferation and increased apoptosis of KB cells were dependent on the presence or absence of the 3'-UTR in CDC73. In summary, knockdown of CDC73 expression due to overexpression of miR-155 not only adds a novelty to the list of mechanisms responsible for its down-regulation in different tumors, but the restoration of CDC73 levels by the use of antagomir-155 may also have an important role in therapeutic intervention of cancers, including OSCC.

Function of microRNA-146a and NF-κB in physiologic and pathologic hematopoiesis

During inflammation and infection, hematopoietic stem and progenitor cells (HSPCs) are stimulated to proliferate and differentiate into mature immune cells, especially of the myeloid lineage. MicroRNA-146a (miR-146a) is a critical negative regulator of inflammation. Deletion of the gene encoding miR-146a—expressed in all blood cell types—produces effects that appear as dysregulated inflammatory hematopoiesis, leading to a decline in the number and quality of hematopoietic stem cells (HSCs), excessive myeloproliferation, and, ultimately, to exhaustion of the HSCs and hematopoietic neoplasms. Six-week-old deleted mice are normal, with no effect on cell numbers, but by 4 months bone marrow hypercellularity can be seen, and by 8 months marrow exhaustion is becoming evident. The ability of HSCs to replenish the entire hematopoietic repertoire in a myelo-ablated mouse also declines precipitously as miR-146a-deficient mice age. In the absence of miR-146a, LPS-mediated serial inflammatory stimulation accelerates the effects of aging. This chronic inflammatory stress on HSCs in deleted mice involves a molecular axis consisting of upregulation of the signaling protein TRAF6 leading to excessive activity of the transcription factor NF-κB and overproduction of the cytokine IL-6. At the cellular level, transplant studies show that the defects are attributable to both an intrinsic problem in the miR-146a-deficient HSCs and extrinsic effects of miR-146a-deficient lymphocytes and non-hematopoietic cells. This study has identified a microRNA, miR-146a, to be a critical regulator of HSC homeostasis during chronic inflammatory challenge in mice and has provided a molecular connection between chronic inflammation and the development of bone marrow failure and myeloproliferative neoplasms. This may have implications for human hematopoietic malignancies, such as myelodysplastic syndrome, which frequently displays downregulated miR-146a expression.

Revealing, illuminating, and modifying proteins in human diseases using noncanonical amino acids

To better understand human diseases, much recent work has focused on proteins to either identify disease targets through proteomics or produce therapeutics via protein engineering. Noncanonical amino acids (ncAAs) are tools for altering the chemical and physical properties of proteins, providing a facile strategy not only to label proteins but also to engineer proteins with novel properties. My thesis research has focused on the development and applications of noncanonical amino acids in identifying, imaging, and engineering proteins for studying human diseases. Chapter 1 introduces the concept of ncAAs and reveals insights to how I chose my thesis projects.

ncAAs have been incorporated to tag and enrich newly synthesized proteins for mass spectrometry through a method termed BONCAT, or bioorthogonal noncanonical amino acid tagging. Chapter 2 describes the investigation of the proteomic response of human breast cancer cells to induced expression of tumor suppressor microRNA miR-126 by combining BONCAT with another proteomic method, SILAC or stable isotope labeling by amino acids in cell culture. This proteomic analysis led to the discovery of a direct target of miR-126, shedding new light on its role in suppressing cancer metastasis.

In addition to mass spectrometry, ncAAs can also be utilized to fluorescently label proteins. Chapter 3 details the synthesis of a set of cell-permeant cyclooctyne probes and demonstration of selective labeling of newly synthesized proteins in live mammalian cells using azidohomoalanine. Similar to live cell imaging, the ability to selectively label a particular cell type within a mixed cell population is important to interrogating many biological systems, such as tumor microenvironments. By taking advantage of the metabolic differences between cancer and normal cells, Chapter 5 discusses efforts to develop selective labeling of cancer cells using a glutamine analogue.

Furthermore, Chapter 4 describes the first demonstration of global replacement at polar amino acid positions and its application in developing an alternative PEGylation strategy for therapeutic proteins. Polar amino acids typically occupy solvent-exposed positions on the protein surface, and incorporation of noncanonical amino acids at these positions should allow easier modification and cause less perturbation compared to replacements at the interior positions of proteins.

The genomic analysis of erythrocyte microRNA expression in sickle cell diseases.

BACKGROUND: Since mature erythrocytes are terminally differentiated cells without nuclei and organelles, it is commonly thought that they do not contain nucleic acids. In this study, we have re-examined this issue by analyzing the transcriptome of a purified population of human mature erythrocytes from individuals with normal hemoglobin (HbAA) and homozygous sickle cell disease (HbSS). METHODS AND FINDINGS: Using a combination of microarray analysis, real-time RT-PCR and Northern blots, we found that mature erythrocytes, while lacking ribosomal and large-sized RNAs, contain abundant and diverse microRNAs. MicroRNA expression of erythrocytes was different from that of reticulocytes and leukocytes, and contributed the majority of the microRNA expression in whole blood. When we used microRNA microarrays to analyze erythrocytes from HbAA and HbSS individuals, we noted a dramatic difference in their microRNA expression pattern. We found that miR-320 played an important role for the down-regulation of its target gene, CD71 during reticulocyte terminal differentiation. Further investigation revealed that poor expression of miR-320 in HbSS cells was associated with their defective downregulation CD71 during terminal differentiation. CONCLUSIONS: In summary, we have discovered significant microRNA expression in human mature erythrocytes, which is dramatically altered in HbSS erythrocytes and their defect in terminal differentiation. Thus, the global analysis of microRNA expression in circulating erythrocytes can provide mechanistic insights into the disease phenotypes of erythrocyte diseases.

MicroRNA 203 Expression in Keratinocytes Is Dependent on Regulation of p53 Levels by E6.

A screen of microRNA (miRNA) expression following differentiation in human foreskin keratinocytes (HFKs) identified changes in several miRNAs, including miRNA 203 (miR-203), which has previously been shown to play an important role in epithelial cell biology by regulating p63 levels. We investigated how expression of human papillomavirus type 16 (HPV16) oncoproteins E6 and E7 affected miR-203 expression during proliferation and differentiation of HFKs. We demonstrated that miR-203 expression is reduced in HFKs where p53 function is compromised, either by the viral oncoprotein E6 or by knockout of p53 using short hairpin RNAs (p53i). We show that the induction of miR-203 observed during calcium-induced differentiation of HFKs is significantly reduced in HFKs expressing E6 and in p53i HFKs. Induction of miR-203 in response to DNA damage is also reduced in the absence of p53. We report that proliferation of HFKs is dependent on the level of miR-203 expression and that overexpression of miR-203 can reduce overproliferation in E6/E7-expressing and p53i HFKs. In summary, these results indicate that expression of miR-203 is dependent on p53, which may explain how expression of HPV16 E6 can disrupt the balance between proliferation and differentiation, as well as the response to DNA damage, in keratinocytes.

MicroRNA-34c inversely couples the biological functions of the runt-related transcription factor RUNX2 and the tumor suppressor p53 in osteosarcoma

Osteosarcoma (OS) is a primary bone tumor that is most prevalent during adolescence. RUNX2, which stimulates differentiation and suppresses proliferation of osteoblasts, is deregulated in OS. Here, we define pathological roles of RUNX2 in the etiology of OS and mechanisms by which RUNX2 expression is stimulated. RUNX2 is often highly expressed in human OS biopsies and cell lines. Small interference RNA (siRNA)-mediated depletion of RUNX2 inhibits growth of U2OS OS cells. RUNX2 levels are inversely linked to loss of p53 (which predisposes to OS) in distinct OS cell lines and osteoblasts. RUNX2 protein levels decrease upon stabilization of p53 with the MDM2 inhibitor Nutlin-3. Elevated RUNX2 protein expression is post-transcriptionally regulated and directly linked to diminished expression of several validated RUNX2 targeting microRNAs (miRNAs) in human OS cells compared to mesenchymal progenitor cells. The p53-dependent miR-34c is the most significantly down-regulated RUNX2 targeting miRNA in OS. Exogenous supplementation of miR-34c markedly decreases RUNX2 protein levels, while 3UTR reporter assays establish RUNX2 as a direct target of miR-34c in OS cells. Importantly, Nutlin-3 mediated stabilization of p53 increases expression of miR-34c and decreases RUNX2. Thus, a novel RUNX2-p53-miR34 network controls cell growth of osseous cells and is compromised in OS.

Human factors required for mycobacterium tuberculosis macrophage infection

Tese de doutoramento, Farmácia (Biologia Celular e Molecular), Universidade de Lisboa, Faculdade de Farmácia, 2014

MicroRNA-34 induces cardiomyocyte apoptosis and accounts for the anti-apoptotic effect of Tanshinone IIA in myocardial infarction

MicroARN (miARN) ont récemment émergé comme un acteur central du gène réseau de régulation impliqués dans la prise du destin cellulaire. L'apoptose, un actif processus, par lequel des cellules déclenchent leur auto-destruction en réponse à un signal, peut être contrôlé par les miARN. Il a également été impliqué dans une variété de maladies humaines, comme les maladies du cœur, et a été pensé comme une cible pour le traitement de la maladie. Tanshinone IIA (TIIA), un monomère de phenanthrenequinones utilisé pour traiter maladies cardiovasculaires, est connu pour exercer des effets cardioprotecteurs de l'infarctus du myocarde en ciblant l'apoptose par le renforcement de Bcl-2 expression. Pour explorer les liens potentiels entre le miARN et l'action anti-apoptotique de TIIA, nous étudié l'implication possible des miARN. Nous avons constaté que l'expression de tous les trois membres de la famille miR-34, miR-34a, miR-34b et miR-34c ont été fortement régulée à la hausse après l'exposition soit à la doxorubicine, un agent endommageant l'ADN ou de pro-oxydant H2O2 pendant 24 heures. Cette régulation à la hausse causé significativement la mort cellulaire par apoptose, comme déterminé par fragmentation de l'ADN, et les effets ont été renversés par les ARNs antisens de ces miARN. Le prétraitement des cellules avec TIIA avant l'incubation avec la doxorubicine ou H2O2 a empêché surexpression de miR-34 et a réduit des apoptose. Nous avons ensuite établi BCL2L2, API5 et TCL1, en plus de BCL2, comme les gènes nouveaux cibles pour miR-34. Nous avons également élucidé que la répression des ces gènes par MiR-34 explique l'effet proapoptotique dans les cardiomyocytes. Ce que la régulation positive de ces gènes par TIIA realisée par la répression de l'expression de miR-34 est probable le mécanisme moléculaire de son effet bénéfique contre ischémique lésions cardiaques.

The Role of MicroRNA Regulation of Cardiac Ion Channel in Arrhythmia

La fibrillation auriculaire (FA) est le trouble du rythme le plus fréquemment observé en pratique clinique. Elle constitue un risque important de morbi-mortalité. Le traitement de la FA reste un défi majeur en lien avec les nombreux effets secondaires associés aux approches thérapeutiques actuelles. Dans ce contexte, une meilleure compréhension des mécanismes sous-jacents à la FA est essentielle pour le développement de nouvelles thérapies offrant un meilleur rapport bénéfice/risque pour les patients. La FA est caractérisée par i) un remodelage électrique délétère associé le plus souvent ii) à un remodelage structurel du myocarde favorisant la récurrence et le maintien de l’arythmie. La diminution de la période réfractaire effective au sein du tissu auriculaire est un élément clef du remodelage électrique. Le remodelage structurel, quant à lui, se manifeste principalement par une fibrose tissulaire qui altère la propagation de l’influx électrique dans les oreillettes. Les mécanismes moléculaires impliqués dans la mise en place de ces deux substrats restent mal connus. Récemment, le rôle des microARNs (miARNs) a été pointé du doigt dans de nombreuses pathologies notamment cardiaques. Dans ce contexte les objectifs principaux de ce travail ont été i) d'acquérir une compréhension approfondie du rôle des miARNs dans la régulation de l’expression des canaux ioniques et ii) de mieux comprendre le rôle de ces molécules dans l’installation d’un substrat favorable a la FA. Nous avons, dans un premier temps, effectué une analyse bio-informatique combinée à des approches expérimentales spécifiques afin d’identifier clairement les miARNs démontrant un fort potentiel de régulation des gènes codant pour l’expression des canaux ioniques cardiaques humains. Nous avons identifié un nombre limité de miARNs cardiaques qui possédaient ces propriétés. Sur la base de ces résultats, nous avons démontré que l’altération de l'expression des canaux ioniques, observée dans diverse maladies cardiaques (par exemple, les cardiomyopathies, l’ischémie myocardique, et la fibrillation auriculaire), peut être soumise à ces miARNs suggérant leur implication dans l’arythmogénèse. La régulation du courant potassique IK1 est un facteur déterminant du remodelage électrique auriculaire associée à la FA. Les mécanismes moléculaires sous-jacents sont peu connus. Nous avons émis l’hypothèse que l'altération de l’expression des miARNs soit corrélée à l’augmentation de l’expression d’IK1 dans la FA. Nous avons constaté que l’expression de miR-26 est réduite dans la FA et qu’elle régule IK1 en modulant l’expression de sa sous-unité Kir2.1. Nous avons démontré que miR-26 est sous la répression transcriptionnelle du facteur nucléaire des lymphocytes T activés (NFAT) et que l’activité accrue de NFATc3/c4, aboutit à une expression réduite de miR-26. En conséquence IK1 augmente lors de la FA. Nous avons enfin démontré que l’interférence in vivo de miR-26 influence la susceptibilité à la FA en régulant IK1, confirmant le rôle prépondérant de miR-26 dans le remodelage auriculaire électrique. La fibrose auriculaire est un constituant majeur du remodelage structurel associé à la FA, impliquant l'activation des fibroblastes et l’influx cellulaire du Ca2 +. Nous avons cherché à déterminer i) si le canal perméable au Ca2+, TRPC3, jouait un rôle dans la fibrose auriculaire en favorisant l'activation des fibroblastes et ii) étudié le rôle potentiel des miARNs dans ce contexte. Nous avons démontré que les canaux TRPC3 favorisent l’influx du Ca2 +, activant la signalisation Ca2 +-dépendante ERK et en conséquence activent la prolifération des fibroblastes. Nous avons également démontré que l’expression du TRPC3 est augmentée dans la FA et que le blocage in vivo de TRPC3 empêche le développement de substrats reliés à la FA. Nous avons par ailleurs validé que miR-26 régule les canaux TRPC3 en diminuant leur expression dans les fibroblastes. Enfin, nous avons montré que l'expression réduite du miR-26 est également due à l’activité augmentée de NFATc3/c4 dans les fibroblastes, expliquant ainsi l’augmentation de TRPC3 lors de la FA, confirmant la contribution de miR-26 dans le processus de remodelage structurel lié à la FA. En conclusion, nos résultats mettent en évidence l'importance des miARNs dans la régulation des canaux ioniques cardiaques. Notamment, miR-26 joue un rôle important dans le remodelage électrique et structurel associé à la FA et ce, en régulant IK1 et l’expression du canal TRPC3. Notre étude démasque ainsi un mécanisme moléculaire de contrôle de la FA innovateur associant des miARNs. miR-26 en particulier représente apres ces travaux une nouvelle cible thérapeutique prometteuse pour traiter la FA.

Systems biology of the human MHC class I immunopeptidome

Le système de différenciation entre le « soi » et le « non-soi » des vertébrés permet la détection et le rejet de pathogènes et de cellules allogéniques. Il requiert la surveillance de petits peptides présentés à la surface cellulaire par les molécules du complexe majeur d’histocompatibilité de classe I (CMH I). Les molécules du CMH I sont des hétérodimères composés par une chaîne lourde encodée par des gènes du CMH et une chaîne légère encodée par le gène β2-microglobuline. L’ensemble des peptides est appelé l’immunopeptidome du CMH I. Nous avons utilisé des approches en biologie de systèmes pour définir la composition et l’origine cellulaire de l’immunopeptidome du CMH I présenté par des cellules B lymphoblastoïdes dérivés de deux pairs de fratries avec un CMH I identique. Nous avons découvert que l’immunopeptidome du CMH I est spécifique à l’individu et au type cellulaire, qu’il dérive préférentiellement de transcrits abondants, est enrichi en transcrits possédant d’éléments de reconnaissance par les petits ARNs, mais qu’il ne montre aucun biais ni vers les régions génétiques invariables ni vers les régions polymorphiques. Nous avons également développé une nouvelle méthode qui combine la spectrométrie de masse, le séquençage de nouvelle génération et la bioinformatique pour l’identification à grand échelle de peptides du CMH I, dont ceux résultants de polymorphismes nucléotidiques simples non-synonymes (PNS-ns), appelés antigènes mineurs d’histocompatibilité (AMHs), qui sont les cibles de réponses allo-immunitaires. La comparaison de l’origine génomique de l’immunopeptidome de soeurs avec un CMH I identique a révélé que 0,5% des PNS-ns étaient représentés dans l’immunopeptidome et que 0,3% des peptides du CMH I seraient immunogéniques envers une des deux soeurs. En résumé, nous avons découvert des nouveaux facteurs qui modèlent l’immunopeptidome du CMH I et nous présentons une nouvelle stratégie pour l’indentification de ces peptides, laquelle pourrait accélérer énormément le développement d’immunothérapies ciblant les AMHs.

Identifying co-regulating microrna groups

Background: Current miRNA target prediction tools have the common problem that their false positive rate is high. This renders identification of co-regulating groups of miRNAs and target genes unreliable. In this study, we describe a procedure to identify highly probable co-regulating miRNAs and the corresponding co-regulated gene groups. Our procedure involves a sequence of statistical tests: (1) identify genes that are highly probable miRNA targets; (2) determine for each such gene, the minimum number of miRNAs that co-regulate it with high probability; (3) find, for each such gene, the combination of the determined minimum size of miRNAs that co-regulate it with the lowest p-value; and (4) discover for each such combination of miRNAs, the group of genes that are co-regulated by these miRNAs with the lowest p-value computed based on GO term annotations of the genes.
Results: Our method identifies 4, 3 and 2-term miRNA groups that co-regulate gene groups of size at least 3 in human. Our result suggests some interesting hypothesis on the functional role of several miRNAs through a "guilt by association" reasoning. For example, miR-130, miR-19 and miR-101 are known neurodegenerative diseases associated miRNAs. Our 3-term miRNA table shows that miR-130/19/101 form a co-regulating group of rank 22 (p-value =1.16 × 10^-2). Since miR-144 is co-regulating with miR-130, miR-19 and miR-101 of rank 4 (p-value = 1.16 × 10^-2) in our 4-term miRNA table, this suggests hsa-miR-144 may be neurodegenerative diseases related miRNA. Conclusions: This work identifies highly probable co-regulating miRNAs, which are refined from the prediction by computational tools using (1) signal-to-noise ratio to get high accurate regulating miRNAs for every gene, and (2) Gene Ontology to obtain functional related co-regulating miRNA groups. Our result has partly been supported by biological experiments. Based on prediction by TargetScanS, we found highly probable target gene groups in the Supplementary Information. This result might help biologists to find small set of miRNAs for genes of interest rather than huge amount of miRNA set.

Promotion of Hendra virus replication by microRNA 146a

Hendra virus is a highly pathogenic zoonotic paramyxovirus in the genus Henipavirus. Thirty-nine outbreaks of Hendra virus have been reported since its initial identification in Queensland, Australia, resulting in seven human infections and four fatalities. Little is known about cellular host factors impacting Hendra virus replication. In this work, we demonstrate that Hendra virus makes use of a microRNA (miRNA) designated miR-146a, an NF-κB-responsive miRNA upregulated by several innate immune ligands, to favor its replication. miR-146a is elevated in the blood of ferrets and horses infected with Hendra virus and is upregulated by Hendra virus in human cells in vitro. Blocking miR-146a reduces Hendra virus replication in vitro, suggesting a role for this miRNA in Hendra virus replication. In silico analysis of miR-146a targets identified ring finger protein (RNF)11, a member of the A20 ubiquitin editing complex that negatively regulates NF-κB activity, as a novel component of Hendra virus replication. RNA interference-mediated silencing of RNF11 promotes Hendra virus replication in vitro, suggesting that increased NF-κB activity aids Hendra virus replication. Furthermore, overexpression of the IκB superrepressor inhibits Hendra virus replication. These studies are the first to demonstrate a host miRNA response to Hendra virus infection and suggest an important role for host miRNAs in Hendra virus disease.