miRDeep*: an integrated application tool for miRNA identification from RNA sequencing data

miRDeep and its varieties are widely used to quantify known and novel micro RNA (miRNA) from small RNA sequencing (RNAseq). This article describes miRDeep*, our integrated miRNA identification tool, which is modeled off miRDeep, but the precision of detecting novel miRNAs is improved by introducing new strategies to identify precursor miRNAs. miRDeep* has a user-friendly graphic interface and accepts raw data in FastQ and Sequence Alignment Map (SAM) or the binary equivalent (BAM) format. Known and novel miRNA expression levels, as measured by the number of reads, are displayed in an interface, which shows each RNAseq read relative to the pre-miRNA hairpin. The secondary pre-miRNA structure and read locations for each predicted miRNA are shown and kept in a separate figure file. Moreover, the target genes of known and novel miRNAs are predicted using the TargetScan algorithm, and the targets are ranked according to the confidence score. miRDeep* is an integrated standalone application where sequence alignment, pre-miRNA secondary structure calculation and graphical display are purely Java coded. This application tool can be executed using a normal personal computer with 1.5 GB of memory. Further, we show that miRDeep* outperformed existing miRNA prediction tools using our LNCaP and other small RNAseq datasets. miRDeep* is freely available online at http://www.australianprostatecentre.org/research/software/mirdeep-star

miRPlant : an integrated tool for identification of plant miRNA from RNA sequencing data

Background Small RNA sequencing is commonly used to identify novel miRNAs and to determine their expression levels in plants. There are several miRNA identification tools for animals such as miRDeep, miRDeep2 and miRDeep*. miRDeep-P was developed to identify plant miRNA using miRDeep’s probabilistic model of miRNA biogenesis, but it depends on several third party tools and lacks a user-friendly interface. The objective of our miRPlant program is to predict novel plant miRNA, while providing a user-friendly interface with improved accuracy of prediction. Result We have developed a user-friendly plant miRNA prediction tool called miRPlant. We show using 16 plant miRNA datasets from four different plant species that miRPlant has at least a 10% improvement in accuracy compared to miRDeep-P, which is the most popular plant miRNA prediction tool. Furthermore, miRPlant uses a Graphical User Interface for data input and output, and identified miRNA are shown with all RNAseq reads in a hairpin diagram. Conclusions We have developed miRPlant which extends miRDeep* to various plant species by adopting suitable strategies to identify hairpin excision regions and hairpin structure filtering for plants. miRPlant does not require any third party tools such as mapping or RNA secondary structure prediction tools. miRPlant is also the first plant miRNA prediction tool that dynamically plots miRNA hairpin structure with small reads for identified novel miRNAs. This feature will enable biologists to visualize novel pre-miRNA structure and the location of small RNA reads relative to the hairpin. Moreover, miRPlant can be easily used by biologists with limited bioinformatics skills.

Mapping the regulatory sequences controlling 93 breast cancer-associated miRNA genes leads to the identification of two functional promoters of the Hsa-mir-200b cluster, methylation of which is associated with metastasis or hormone receptor status in advanced breast cancer

MicroRNAs (miRNAs) are small non-coding RNAs of 20 nt in length that are capable of modulating gene expression post-transcriptionally. Although miRNAs have been implicated in cancer, including breast cancer, the regulation of miRNA transcription and the role of defects in this process in cancer is not well understood. In this study we have mapped the promoters of 93 breast cancer-associated miRNAs, and then looked for associations between DNA methylation of 15 of these promoters and miRNA expression in breast cancer cells. The miRNA promoters with clearest association between DNA methylation and expression included a previously described and a novel promoter of the Hsa-mir-200b cluster. The novel promoter of the Hsa-mir-200b cluster, denoted P2, is located 2 kb upstream of the 5′ stemloop and maps within a CpG island. P2 has comparable promoter activity to the previously reported promoter (P1), and is able to drive the expression of miR-200b in its endogenous genomic context. DNA methylation of both P1 and P2 was inversely associated with miR-200b expression in eight out of nine breast cancer cell lines, and in vitro methylation of both promoters repressed their activity in reporter assays. In clinical samples, P1 and P2 were differentially methylated with methylation inversely associated with miR-200b expression. P1 was hypermethylated in metastatic lymph nodes compared with matched primary breast tumours whereas P2 hypermethylation was associated with loss of either oestrogen receptor or progesterone receptor. Hypomethylation of P2 was associated with gain of HER2 and androgen receptor expression. These data suggest an association between miR-200b regulation and breast cancer subtype and a potential use of DNA methylation of miRNA promoters as a component of a suite of breast cancer biomarkers.

Study of miRNA polymorphisms and their potential association with breast cancer risk in Australian Caucasian populations

This project established a large and well characterised prospective breast cancer DNA biobank and used this biobank to conduct genetic studies in breast cancer. The thesis presented the results of these high-throughput genotyping studies in two separate Australian Caucasian case-control populations and identified association between three novel genetic variants in microRNA genes and breast cancer risk.

A large-scale analysis of genetic variants within putative miRNA binding sites in prostate cancer

Prostate cancer is the second most common malignancy among men worldwide. Genome-wide association studies have identified 100 risk variants for prostate cancer, which can explain approximately 33% of the familial risk of the disease. We hypothesized that a comprehensive analysis of genetic variations found within the 3' untranslated region of genes predicted to affect miRNA binding (miRSNP) can identify additional prostate cancer risk variants. We investigated the association between 2,169 miRSNPs and prostate cancer risk in a large-scale analysis of 22,301 cases and 22,320 controls of European ancestry from 23 participating studies. Twenty-two miRSNPs were associated (P<2.3×10(-5)) with risk of prostate cancer, 10 of which were within 7 genes previously not mapped by GWAS studies. Further, using miRNA mimics and reporter gene assays, we showed that miR-3162-5p has specific affinity for the KLK3 rs1058205 miRSNP T-allele, whereas miR-370 has greater affinity for the VAMP8 rs1010 miRSNP A-allele, validating their functional role. SIGNIFICANCE Findings from this large association study suggest that a focus on miRSNPs, including functional evaluation, can identify candidate risk loci below currently accepted statistical levels of genome-wide significance. Studies of miRNAs and their interactions with SNPs could provide further insights into the mechanisms of prostate cancer risk.

Genetic association analysis of miRNA SNPs implicates MIR145 in breast cancer susceptibility

Background MicroRNAs (miRNAs) are important small non-coding RNA molecules that regulate gene expression in cellular processes related to the pathogenesis of cancer. Genetic variation in miRNA genes could impact their synthesis and cellular effects and single nucleotide polymorphisms (SNPs) are one example of genetic variants studied in relation to breast cancer. Studies aimed at identifying miRNA SNPs (miR-SNPs) associated with breast malignancies could lead towards further understanding of the disease and to develop clinical applications for early diagnosis and treatment. Methods We genotyped a panel of 24 miR-SNPs using multiplex PCR and chip-based matrix assisted laser desorption ionization time-of-flight (MALDI-TOF) mass spectrometry (MS) analysis in two Caucasian breast cancer case control populations (Primary population: 173 cases and 187 controls and secondary population: 679 cases and 301 controls). Association to breast cancer susceptibility was determined using chi-square (X 2 ) and odds ratio (OR) analysis. Results Statistical analysis showed six miR-SNPs to be non-polymorphic and twelve of our selected miR-SNPs to have no association with breast cancer risk. However, we were able to show association between rs353291 (located in MIR145) and the risk of developing breast cancer in two independent case control cohorts (p = 0.041 and p = 0.023). Conclusions Our study is the first to report an association between a miR-SNP in MIR145 and breast cancer risk in individuals of Caucasian background. This finding requires further validation through genotyping of larger cohorts or in individuals of different ethnicities to determine the potential significance of this finding as well as studies aimed to determine functional significance. Keywords: Association analysis; Breast cancer; microRNA; miR-SNPs; MIR145

Integration of miRNA and mRNA expression with DNA copy number of Ewing sarcoma cell lines

Ewing sarcoma is an aggressive and poorly differentiated malignancy of bone and soft tissue. It primarily affects children, adolescents, and young adults, with a slight male predominance. It is characterized by a translocation between chromosomes 11 and 22 resulting in the EWSR1-FLI1fusion transcription factor. The aim of this study is to identify putative Ewing sarcoma target genes through an integrative analysis of three microarray data sets. Array comparative genomic hybridization is used to measure changes in DNA copy number, and analyzed to detect common chromosomal aberrations. mRNA and miRNA microarrays are used to measure expression of protein-coding and miRNA genes, and these results integrated with the copy number data. Chromosomal aberrations typically contain also bystanders in addition to the driving tumor suppressor and oncogenes, and integration with expression helps to identify the true targets. Correlation between expression of miRNAs and their predicted target mRNAs is also evaluated to assess the results of post-transcriptional miRNA regulation on mRNA levels. The highest frequencies of copy number gains were identified in chromosome 8, 1q, and X. Losses were most frequent in 9p21.3, which also showed an enrichment of copy number breakpoints relative to the rest of the genome. Copy number losses in 9p21.3 were found have a statistically significant effect on the expression of MTAP, but not on CDKN2A, which is a known tumor-suppressor in the same locus. MTAP was also down-regulated in the Ewing sarcoma cell lines compared to mesenchymal stem cells. Genes exhibiting elevated expression in association with copy number gains and up-regulation compared to the reference samples included DCAF7, ENO2, MTCP1, andSTK40. Differentially expressed miRNAs were detected by comparing Ewing sarcoma cell lines against mesenchymal stem cells. 21 up-regulated and 32 down-regulated miRNAs were identified, includingmiR-145, which has been previously linked to Ewing sarcoma. The EWSR1-FLI1 fusion gene represses miR-145, which in turn targets FLI1 forming a mutually repressive feedback loop. In addition higher expression linked to copy number gains and compared to mesenchymal stem cells, STK40 was also found to be a target of four different miRNAs that were all down-regulated in Ewing sarcoma cell lines compared to the reference samples. SLCO5A1 was identified as the only up-regulated gene within a frequently gained region in chromosome 8. This region was gained in over 90 % of the cell lines, and also with a higher frequency than the neighboring regions. In addition, SLCO5A1 was found to be a target of three miRNAs that were down-regulated compared to the mesenchymal stem cells.

Circulating miRNA profile in HCV infected serum: novel insight into pathogenesis

Changes in circulating miRNA profiles have been associated with different diseases. Here we demonstrate the circulating miRNA profile in serum of HCV infected individuals using a microRNA array that profiles the expression of 940 miRNAs. Serum samples from two HCV genotype -1 and two HCV genotype -3 infected individuals were compared with healthy controls. Expression levels of miR-134, miR-198, miR-320c and miR-483-5p that were commonly upregulated in case of both genotypes were validated in 36 individual patient serum samples. Serum miR-134, miR-320c and miR-483-5p were significantly upregulated during HCV infection. miR-320c and miR-483-5p were also upregulated in HCV-JFH1 infected cells and cell culture supernatant. Pathway analysis of putative target genes of these miRNAs indicated involvement of PI3K-Akt, NFKB and MAPK signaling pathways. Results revealed novel insights on the role of circulating miRNAs in mediating pathogenesis in HCV-infected cells.

Epigenetics and miRNA during bacteria-induced host immune responses

Various cellular processes including the pathogen-specific immune responses, host-pathogen interactions and the related evasion mechanisms rely on the ability of the immune cells to be reprogrammed accurately and in many cases instantaneously. In this context, the exact functions of epigenetic and miRNA-mediated regulation of genes, coupled with recent advent in techniques that aid such studies, make it an attractive field for research. Here, we review examples that involve the epigenetic and miRNA control of the host immune system during infection with bacteria. Interestingly, many pathogens utilize the epigenetic and miRNA machinery to modify and evade the host immune responses. Thus, we believe that global epigenetic and miRNA mapping of such host-pathogen interactions would provide key insights into their cellular functions and help to identify various determinants for therapeutic value.

miRNA polymorphisms are associated with risk of developing chronic lymphocytic leukemia

Chronic Lymphocytic Leukemia (CLL) is the most frequent leukemia of adults in Western countries and shows a ~8.5-fold increased relative risk in first-degree relatives. Up to date several studies have identified low-penetrance susceptibility alleles in CLL. Nevertheless, these studies scarcely study regions that do not encode proteins such as microRNAs (miRNAs). Abnormalities in miRNAs, as altered expression patterns and mutations, have been described in CLL, suggesting their implication in the development of the disease. Polymorphisms in these miRNAs may deregulate miRNAs expression levels and affect to the miRNA function. However, despite accumulating evidence that inherited genetic variation in miRNA genes can contribute to the predisposition for CLL, the role of these in the risk of CLL has not been extensively studied. Therefore, the aim of this study was to find new genetic markers of risk to CLL. To that end, we made a systematic search for SNPs in miRNAs and miRNAs deregulated in CLL and genotyped 213 polymorphisms in 401 samples of Spanish individuals. The literature search resulted in more than 100 miRNAs deregulated in CLL and 43 polymorphisms studied in the disease. Out of 213 genotyped SNPs, 13 showed to be significantly associated with CLL risk. rs2682818 in pre-mature miR618 was the most significant result, with 0.49 fold decreased risk to CLL. Interestingly, a previous study associated this SNP with an increased risk of developing follicular lymphoma. Secondly, rs10173558 SNP in mir- 1302-4 showed the highest risk association, with a 5.24 fold increased risk, but there were no previous works studying it. Finally, rs61992671 in miR412, previously associated with CLL risk, showed also association in our sample. In conclusion, we find 13 alleles which could contribute to the risk of CLL. However, new large-scale studies including functional analyses will be needed to validate our findings.

Genetic Variants in MiRNA Processing Genes and Pre-MiRNAs Are Associated with the Risk of Chronic Lymphocytic Leukemia

Genome wide association studies (GWAS) have identified several low-penetrance susceptibility alleles in chronic lymphocytic leukemia (CLL). Nevertheless, these studies scarcely study regions that are implicated in non-coding molecules such as microRNAs (miRNAs). Abnormalities in miRNAs, as altered expression patterns and mutations, have been described in CLL, suggesting their implication in the development of the disease. Genetic variations in miRNAs can affect levels of miRNA expression if present in pre-miRNAs and in miRNA biogenesis genes or alter miRNA function if present in both target mRNA and miRNA sequences. Therefore, the present study aimed to evaluate whether polymorphisms in pre-miRNAs, and/or miRNA processing genes contribute to predisposition for CLL. A total of 91 SNPs in 107 CLL patients and 350 cancer-free controls were successfully analyzed using TaqMan Open Array technology. We found nine statistically significant associations with CLL risk after FDR correction, seven in miRNA processing genes (rs3805500 and rs6877842 in DROSHA, rs1057035 in DICER1, rs17676986 in SND1, rs9611280 in TNRC6B, rs784567 in TRBP and rs11866002 in CNOT1) and two in pre-miRNAs (rs11614913 in miR196a2 and rs2114358 in miR1206). These findings suggest that polymorphisms in genes involved in miRNAs biogenesis pathway as well as in pre-miRNAs contribute to the risk of CLL. Large-scale studies are needed to validate the current findings.

单细胞真核生物的miRNA系统及其进化意义

miRNA系统在高等多细胞真核生物中得到了广泛深入的研究。近年来,人们在单细胞真核生物上的miRNA研究也取得了重要进展。这不仅丰富了人们对miRNA在整个生物界中的认识,更重要的是对于揭示miRNA这一表达调节系统是如何在生物界中起源进化的问题具有重要意义。该文结合作者在最低等单细胞真核生物——贾第虫上的研究结果,对该领域的研究进展作一概述,并对有关miRNA这一系统的起源进化问题进行了探讨。

大脑表达 miRNA 在大脑进化中的作用及其遗传机制初步研究

人类区别于其它动物的最本质特征是其拥有其他动物所无法比拟的大脑容量及高 级的认知能力。即使与其近亲－非人灵长类相比，人也拥有比非人灵长类大好几倍的脑 容量和更为发达的认知能力。现在一般认为，人类大脑的形成是适应性选择(达尔文正 选择)的结果。但是到目前为止，对人类起源过程中大脑容量增大及认知能力提高的遗 传学机制却知之甚少。以前的研究表明，几个与大脑发育相关的蛋白质编码基因在人类 起源中受到了正向选择。同时，也有证据表明人类大脑的进化也可能是基因表达调控变 化的结果。因此寻找人与非人灵长类大脑表达基因的调控差异或许能够进一步为人与非 人灵长类为何有如此巨大的差异提供分子生物学水平上的解释。 MicroRNA(miRNA)是一类在转录后水平调控基因表达的不编码蛋白质的小 RNA(长度20－24个碱基)。通过调控靶基因的表达，miRNA参与了众多的生理过程。而 很多大脑表达的miRNA的表达量在大脑发育过程中呈显著的变化，这表明miRNA参与 了大脑的发育调控。由此可知，miRNA对其靶基因调控效率的改变很可能引起大脑发育 调控的改变。 本文通过寻找大脑表达的保守microRNA 及对其靶基因的预测，同时用比较基因组 学的方法发现：有很多microRNA 的靶基因3’UTR 的靶位点中存在人类特异突变位点。 我们推测这些人类特异突变位点可能改变microRNA 对其靶基因的表达调控效应，而调 控效应的改变可能在进化过程中对认知能力的提高发挥重要作用。利用体外报告基因系 统，我们发现有几个预测的靶基因是对应miRNA 的真实靶基因。其中，miR-127 的靶 基因（SEMA3F）3’UTR 的靶位点中所含的一个人类特有的突变增强了miR-127 对 SEMA3F 的调控效率。将该位点突变回复为黑猩猩的位点使得miR-127 对SEMA3F 的 调控效率降低到黑猩猩的水平。这表明miR-127 对SEMA3F 的调控效率的改变确实是 由该人类特异突变位点引起的。我们提供了人类特异突变位点能够引起miR-127 对SEMA3F 的调控效率的改变的体外证据，但是体内的调控模式是否如此尚需进一步的工 作。 总之，本文通过体外试验表明，miRNA靶基因3’UTR的序列变异具有功能效应，它 有可能是人类中枢神经系统在起源和演化中发挥关键作用的重要遗传机制之一。

Udział miRNA w regulacji alternatywnego splicingu u roślin

Wydział Biologii

Prediction and Verification of miRNA Expression in Human and Rat Retinas.

PURPOSE: MicroRNAs (miRNAs) play a global role in regulating gene expression and have important tissue-specific functions. Little is known about their role in the retina. The purpose of this study was to establish the retinal expression of those miRNAs predicted to target genes involved in vision. METHODS: miRNAs potentially targeting important "retinal" genes, as defined by expression pattern and implication in disease, were predicted using a published algorithm (TargetScan; Envisioneering Medical Technologies, St. Louis, MO). The presence of candidate miRNAs in human and rat retinal RNA was assessed by RT-PCR. cDNA levels for each miRNA were determined by quantitative PCR. The ability to discriminate between miRNAs varying by a single nucleotide was assessed. The activity of miR-124 and miR-29 against predicted target sites in Rdh10 and Impdh1 was tested by cotransfection of miRNA mimics and luciferase reporter plasmids. RESULTS: Sixty-seven miRNAs were predicted to target one or more of the 320 retinal genes listed herein. All 11 candidate miRNAs tested were expressed in the retina, including miR-7, miR-124, miR135a, and miR135b. Relative levels of individual miRNAs were similar between rats and humans. The Rdh10 3'UTR, which contains a predicted miR-124 target site, mediated the inhibition of luciferase activity by miR-124 mimics in cell culture. CONCLUSIONS: Many miRNAs likely to regulate genes important for retinal function are present in the retina. Conservation of miRNA retinal expression patterns from rats to humans supports evidence from other tissues that disruption of miRNAs is a likely cause of a range of visual abnormalities.