Comprehensive analysis of single nucleotide polymorphisms in human microRNAs

Phenome-wide association study on miRNA-related sequence variants: the UK Biobank

BACKGROUND

Genetic variants in the coding region could directly affect the structure and expression levels of genes and proteins. However, the importance of variants in the non-coding region, such as microRNAs (miRNAs), remain to be elucidated. Genetic variants in miRNA-related sequences could affect their biogenesis or functionality and ultimately affect disease risk. Yet, their implications and pleiotropic effects on many clinical conditions remain unknown.

METHODS

Here, we utilised genotyping and hospital records data in the UK Biobank (N = 423,419) to investigate associations between 346 genetic variants in miRNA-related sequences and a wide range of clinical diagnoses through phenome-wide association studies. Further, we tested whether changes in blood miRNA expression levels could affect disease risk through colocalisation and Mendelian randomisation analysis.

RESULTS

We identified 122 associations for six variants in the seed region of miRNAs, nine variants in the mature region of miRNAs, and 27 variants in the precursor miRNAs. These included associations with hypertension, dyslipidaemia, immune-related disorders, and others. Nineteen miRNAs were associated with multiple diagnoses, with six of them associated with multiple disease categories. The strongest association was reported between rs4285314 in the precursor of miR-3135b and celiac disease risk (odds ratio (OR) per effect allele increase = 0.37, P = 1.8 × 10-162). Colocalisation and Mendelian randomisation analysis highlighted potential causal role of miR-6891-3p in dyslipidaemia.

CONCLUSIONS

Our study demonstrates the pleiotropic effect of miRNAs and offers insights to their possible clinical importance.

Identification of microRNA editing sites in clear cell renal cell carcinoma

Clear cell renal cell carcinoma (ccRCC) is a malignant tumor originating from the renal tubular epithelium. Although the microRNAs (miRNAs) transcriptome of ccRCC has been extensively studied, the role of miRNAs editing in ccRCC is largely unknown. By analyzing small RNA sequencing profiles of renal tissues of 154 ccRCC patients and 22 normal controls, we identified 1025 miRNA editing sites from 246 pre-miRNAs. There were 122 editing events with significantly different editing levels in ccRCC compared to normal samples, which include two A-to-I editing events in the seed regions of hsa-mir-376a-3p and hsa-mir-376c-3p, respectively, and one C-to-U editing event in the seed region of hsa-mir-29c-3p. After comparing the targets of the original and edited miRNAs, we found that hsa-mir-376a-1_49g, hsa-mir-376c_48g and hsa-mir-29c_59u had many new targets, respectively. Many of these new targets were deregulated in ccRCC, which might be related to the different editing levels of hsa-mir-376a-3p, hsa-mir-376c-3p, hsa-mir-29c-3p in ccRCC compared to normal controls. Our study sheds new light on miRNA editing events and their potential biological functions in ccRCC.

Phylogenomic Analysis of micro-RNA Involved in Juvenile to Flowering-Stage Transition in Photophilic Rice and Its Sister Species

Vegetative to reproductive phase transition in phototropic plants is an important developmental process and is sequentially mediated by the expression of micro-RNA MIR172. To obtain insight into the evolution, adaptation, and function of MIR172 in photophilic rice and its wild relatives, we analyzed the genescape of a 100 kb segment harboring MIR172 homologs from 11 genomes. The expression analysis of MIR172 revealed its incremental accumulation from the 2-leaf to 10-leaf stage, with maximum expression coinciding with the flag-leaf stage in rice. Nonetheless, the microsynteny analysis of MIR172s revealed collinearity within the genus Oryza, but a loss of synteny was observed in (i) MIR172A in O. barthii (AA) and O. glaberima (AA); (ii) MIR172B in O. brachyantha (FF); and (iii) MIR172C in O. punctata (BB). Phylogenetic analysis of precursor sequences/region of MIR172 revealed a distinct tri-modal clade of evolution. The genomic information generated in this investigation through comparative analysis of MIRNA, suggests mature MIR172s to have evolved in a disruptive and conservative mode amongst all Oryza species with a common origin of descent. Further, the phylogenomic delineation provided an insight into the adaptation and molecular evolution of MIR172 to changing environmental conditions (biotic and abiotic) of phototropic rice through natural selection and the opportunity to harness untapped genomic regions from rice wild relatives (RWR).

Characterizing microRNA editing and mutation sites in Autism Spectrum Disorder

Autism Spectrum Disorder (ASD) is a neurodevelopmental disorder whose pathogenesis is still unclear. MicroRNAs (miRNAs) are a kind of endogenous small non-coding RNAs that play important roles in the post-transcriptional regulation of genes. Recent researches show that miRNAs are edited in multiple ways especially in central nervous systems. A-to-I editing of RNA catalyzed by Adenosine deaminases acting on RNA (ADARs) happens intensively in brain and is also noticed in other organs and tissues. Although miRNAs are widely edited in human brain, miRNA editing in ASD is still largely unexplored. In order to reveal the editing events of miRNAs in ASD, we analyzed 131 miRNA-seq samples from 8 different brain regions of ASD patients and normal controls. We identified 834 editing sites with significant editing levels, of which 70 sites showed significantly different editing levels in the superior frontal gyrus samples of ASD patients (ASD-SFG) when compared with those of control samples. The editing level of an A-to-I editing site in hsa-mir-376a-1 (hsa-mir-376a-1_9_A_g) in ASD-SFG is higher than that of normal controls, and the difference is exaggerated in individuals under 10 years. The increased expression of ADAR1 is consistent with the increased editing level of hsa-mir-376a-1_9_A_g in ASD-SFG samples compared to normal SFG samples. Furthermore, we verify that A-to-I edited hsa-mir-376a-5p directly represses GPR85 and NAPB, which may contribute to the abnormal neuronal development of ASD patients. These results provide new insights into the mechanism of ASD.

Key methodological challenges in detecting circulating miRNAs in different biofluids

The technological advancement in diagnostic techniques has immensely improved the capability of predicting disease progression. Yet, there is a great interest in developing newer biomarkers that can enhance disease risk prediction thereby minimising the associated morbidity and mortality. Circulating miRNAs, a non-coding RNA molecule, are critical regulators in the pathophysiology of various complex multifactorial diseases. In recent years, circulating miRNAs have been enormously studied and are considered as an emerging biomarker due to their easy accessibility, stability, and detection by sequence-specific amplification methods. However, there is a distinct lack of consensus regarding the preanalytical factors such as preferred sample selection, methodological aspects, etc that may independently or together influence the detection of circulating miRNAs resulting in erroneous expression profiles. Therefore, the present review makes an attempt to highlight the various pre-analytical and analytical factors that can potentially influence the circulating miRNA levels. Literature on circulating miRNA's stability, processing and quantitation in different biofluids along with the effect of various controllable and uncontrollable factors influencing circulating miRNA expression have been summarised in the current review.

Characterizing Relevant MicroRNA Editing Sites in Parkinson's Disease

MicroRNAs (miRNAs) are extensively edited in human brains. However, the functional relevance of the miRNA editome is largely unknown in Parkinson's disease (PD). By analyzing small RNA sequencing profiles of brain tissues of 43 PD patients and 88 normal controls, we found that the editing levels of five A-to-I and two C-to-U editing sites are significantly correlated with the ages of normal controls, which is disrupted in PD patients. We totally identified 362 miRNA editing sites with significantly different editing levels in prefrontal cortices of PD patients (PD-PC) compared to results of normal controls. We experimentally validated that A-to-I edited miR-497-5p, with significantly higher expression levels in PD-PC compared to normal controls, directly represses OPA1 and VAPB. Furthermore, overexpression of A-to-I edited miR-497-5p downregulates OPA1 and VAPB in two cell lines, and inhibits proliferation of glioma cells. These results suggest that the hyperediting of miR-497-5p in PD contributes to enhanced progressive neurodegeneration of PD patients. Our results provide new insights into the mechanistic understanding, novel diagnostics, and therapeutic clues of PD.

A Single Variant in Pri-miRNA-155 Associated with Susceptibility to Hereditary Breast Cancer Promotes Aggressiveness in Breast Cancer Cells

Variants in genes encoding for microRNAs have been associated with their deregulation in breast cancer (BC). Sequencing of microRNAs deregulated in BC was performed using DNA from Chilean patients with a strong family history and negative for mutations in BRCA1/BRCA2. Seventeen variants were identified, three of which were selected for a case-control association study: rs376491654 (miR-335), rs755634302 (miR-497), and rs190708267 (miR-155). For rs190708267 C>T, the heterozygous T allele was detected in four BC cases and absent in controls, while homozygous TT cases were not detected. Variants were modelled in silico, cloned in a plasmid, expressed in BC cell lines, and functional in vitro assays were performed. Overexpression of the miR-155-T allele increased mature miR-155-5p levels in both BC cell lines, suggesting that its presence alters pre-miR-155 processing. Moreover, BC cells overexpressing the miR-155-T allele showed increased proliferation, migration, and resistance to cisplatin-induced death compared to miR-155-C overexpressing cells. Of note, the 3′UTR of APC, GSK3β, and PPP1CA genes, all into the canonical Wnt signaling pathway, were identified as direct targets. APC and GSK3β mRNA levels decreased while PP1 levels increased. These results suggest a pathogenic role of the variant rs190708267 (miR-155) in BRCA 1/2 negative BC, conferring susceptibility and promoting traits of aggressiveness.

Identification of microRNA editing sites in three subtypes of leukemia

Leukemia is an aberrant hyper-proliferation of immature blood cells that do not form solid tumors. The transcriptomes of microRNAs (miRNAs) of leukemia have been intensively explored. However, miRNA editing of leukemia has not been extensively studied. To identify miRNA editing patterns and explore their functional relevance in leukemia, we analyzed 200 small RNA sequencing profiles of three subtypes of leukemia and identified hundreds of miRNA editing sites in three subtypes of leukemia. Then, we compared the editing levels of identified miRNA editing sites in leukemia and normal controls. Many miRNAs were differential edited in different subtypes of leukemia. We also found the editing levels of 3'-A editing sites of hsa-mir-21-5p and hsa-mir-155-5p decreased in chronic lymphocytic leukemia patients with radiation treatments. By integrating PAR-CLIP sequencing profiles, we predicted the targets of original and edited miRNAs. One of the edited miRNA, hsa-let-7b_5c, with an additional cytosine at 5' end of hsa-let-7b-5p, potentially targeted VBP1 and CTDSP1. CTDSP1 was significantly downregulated in T-ALL compared to normal controls, which might be originated from the hyperediting of hsa-let-7b-5p in T-ALL. Our study provides a comprehensive view of miRNA editing in three different subtypes of leukemia.

Analysis of Homologous Regions of Small RNAs MIR397 and MIR408 Reveals the Conservation of Microsynteny among Rice Crop-Wild Relatives

MIRNAs are small non-coding RNAs that play important roles in a wide range of biological processes in plant growth and development. MIR397 (involved in drought, low temperature, and nitrogen and copper (Cu) starvation) and MIR408 (differentially expressed in response to environmental stresses such as copper, light, mechanical stress, dehydration, cold, reactive oxygen species, and drought) belong to conserved MIRNA families that either negatively or positively regulate their target genes. In the present study, we identified the homologs of MIR397 and MIR408 in Oryza sativa and its six wild progenitors, three non-Oryza species, and one dicot species. We analyzed the 100 kb segments harboring MIRNA homologs from 11 genomes to obtain a comprehensive view of their community evolution around these loci in the farthest (distant) relatives of rice. Our study showed that mature MIR397 and MIR408 were highly conserved among all Oryza species. Comparative genomics analyses also revealed that the microsynteny of the 100 kb region surrounding MIRNAs was only conserved in Oryza spp.; disrupted in Sorghum, maize, and wheat; and completely lost in Arabidopsis. There were deletions, rearrangements, and translocations within the 100 kb segments in Oryza spp., but the overall microsynteny of the region was maintained. The phylogenetic analyses of the precursor regions of all MIRNAs under study revealed a bimodal clade of common origin. This comparative analysis of miRNA involved in abiotic stress tolerance in plants provides a powerful tool for future Oryza research. Crop wild relatives (CWRs) offer multiple traits with potential to decrease the amount of yield loss owing to biotic and abiotic stresses. Using a comparative genomics approach, the exploration of CWRs as a source of tolerance to these stresses by understanding their evolution can be further used to leverage their yield potential.

A microRNA-microRNA crosstalk network inferred from genome-wide single nucleotide polymorphism variants in natural populations of Arabidopsis thaliana

To adapt to variable natural conditions, plants have evolved several strategies to respond to different environmental stresses. MicroRNA (miRNA)-mediated gene regulation is one of such strategies. Variants, e.g., single nucleotide polymorphisms (SNPs) within the mature miRNAs or their target sites may cause the alteration of regulatory networks and serious phenotype changes. In this study, we proposed a novel approach to construct a miRNA-miRNA crosstalk network in Arabidopsis thaliana based on the notion that two cooperative miRNAs toward common targets are under a strong pressure to be inherited together across ecotypes. By performing a genome-wide scan of the SNPs within the mature miRNAs and their target sites, we defined a "regulation fate profile" to describe a miRNA-target regulation being static (kept) or dynamic (gained or lost) across 1,135 ecotypes compared with the reference genome of Col-0. The cooperative miRNA pairs were identified by estimating the similarity of their regulation fate profiles toward the common targets. The reliability of the cooperative miRNA pairs was supported by solid expressional correlation, high PPImiRFS scores, and similar stress responses. Different combinations of static and dynamic miRNA-target regulations account for the cooperative miRNA pairs acting on various biological characteristics of miRNA conservation, expression, homology, and stress response. Interestingly, the targets that are co-regulated dynamically by both cooperative miRNAs are more likely to be responsive to stress. Hence, stress-related genes probably bear selective pressures in a certain group of ecotypes, in which miRNA regulations on the stress genes reprogram. Finally, three case studies showed that reprogramming miRNA-miRNA crosstalk toward the targets in specific ecotypes was associated with these ecotypes' climatic variables and geographical locations. Our study highlights the potential of miRNA-miRNA crosstalk as a genetic basis underlying environmental adaptation in natural populations.

MiRLog and dbmiR: prioritization and functional annotation tools to study human microRNA sequence variants

The recent identification of noncoding variants with pathogenic effects suggests that these variations could underlie a significant number of undiagnosed cases. Several computational methods have been developed to predict the functional impact of noncoding variants, but they exhibit only partial concordance and are not integrated with functional annotation resources, making the interpretation of these variants still challenging. MicroRNAs (miRNAs) are small noncoding RNA molecules that act as fine regulators of gene expression and play crucial functions in several biological processes, such as cell proliferation and differentiation. An increasing number of studies demonstrate a significant impact of miRNA single nucleotide variants (SNVs) both in Mendelian diseases and complex traits. To predict the functional effect of miRNA SNVs, we implemented a new meta‐predictor, MiRLog, and we integrated it into a comprehensive database, dbmiR, which includes a precompiled list of all possible miRNA allelic SNVs, providing their biological annotations at nucleotide and miRNA levels. MiRLog and dbmiR were used to explore the genetic variability of miRNAs in 15,708 human genomes included in the gnomAD project, finding several ultra‐rare SNVs with a potentially deleterious effect on miRNA biogenesis and function representing putative contributors to human phenotypes.

Interleukin 27 Signaling in Rheumatoid Arthritis Patients: Good or Evil?

The occurrence and development of rheumatoid arthritis (RA) is regulated by numerous cytokines. Interleukin 27 (IL-27) is a soluble cytokine that exerts biological effects by regulating the Janus tyrosine kinase (JAK)/signal transducer and activator of the transcription (STAT) signaling pathway via the IL-27 receptor. IL-27 is known for its pleiotropic roles in modulating inflammatory responses. Previous studies found that IL-27 levels are elevated in RA blood, synovial fluid, and rheumatoid nodules. Cellular and animal experiments indicated that IL-27 exerts multiple regulatory functions in RA patients via different mechanisms. IL-27 inhibits ectopic-like structure (ELS) formation and CD4⁺ T helper type 2 (Th2) cell, CD4⁺ T helper type 17 (Th17) cell, and osteoclast differentiation in RA, contributing to alleviating RA. However, IL-27 promotes Th1 cell differentiation, which may exacerbate RA synovitis. Moreover, IL-27 also acts on RA synovial fibroblasts (RA-FLSs) and regulatory T cells (Tregs), but some of its functions are unclear. There is currently insufficient evidence to determine whether IL-27 promotes or relieves RA. Targeting IL-27 signaling in RA treatment should be deliberate based on current knowledge.

MicroRNA editing patterns in Huntington's disease

Huntington's disease (HD) is a neurodegenerative disease. MicroRNAs (miRNAs) are small non-coding RNAs that mediate post-transcriptional regulation of target genes. Although miRNAs are extensively edited in human brains, the editome of miRNAs in brains of HD patients is largely unknown. By analyzing the small RNA sequencing profiles of brain tissues of 28 HD patients and 83 normal controls, 1182 miRNA editing sites with significant editing levels were identified. In addition to 27 A-to-I editing sites, we identified 3 conserved C-to-U editing sites in miRNAs of HD patients. 30 SNPs in the miRNAs of HD patients were also identified. Furthermore, 129 miRNA editing events demonstrated significantly different editing levels in prefrontal cortex samples of HD patients (HD-PC) when compared to those of healthy controls. We found that hsa-mir-10b-5p was edited to have an additional cytosine at 5'-end in HD-PC, and the edited hsa-mir-10b repressed GTPBP10 that was often downregulated in HD. The down-regulation of GTPBP10 might contribute to the progression of HD by causing gradual loss of function of mitochondrial. These results provide the first endeavor to characterize the miRNA editing events in HD and their potential functions.

Assessment of Association between miR-146a Polymorphisms and Expression of miR-146a, TRAF-6, and IRAK-1 Genes in Patients with Brucellosis

BACKGROUND

Brucellosis is a major zoonosis all over the world. MicroRNAs are significant gene expression regulators and could be involved during the infections and also genetic alterations in the miRNAs sequence can affect primary miRNAs and precursor miRNAs processing and thus alter miRNAs expression. Current research studied the impact of the miR-146a polymorphism on miR-146a, TRAF-6, and IRAK-1 genes expression in patients with brucellosis illness.

METHODS AND RESULTS

In this research, 25 patients with brucellosis and 25 healthy participants with determined genotypes for miR-SNP rs2910164 and miR-SNP rs57095329 were recruited. IRAK-1, TRAF-6, and miR-146a expressions in peripheral blood mononuclear cells (PBMCs) were specified by quantitative real- time PCR (qRT-PCR). Moreover, interleukin-1β (IL-1β) and tumor necrosis factor- alpha (TNF-α) serum levels were assessed by a sandwich enzyme-linked immunosorbent assay (ELISA) technique. There was no significant difference in the expression level of miR-146a, IRAK-1, and TRAF-6, among the patients with brucellosis and control group. TRAF-6 PBMCs expression levels in the distinctive genotypes of rs2910164 were significantly observed in patients (P = 0.048). No significant distinctions were found in miR-146a, IRAK-1, and TRAF-6 expression levels and among the rs57095329 different genotypes in brucellosis patients and controls. Meanwhile, no significant relationship was found between the rs2910164 and rs57095329 genotypes and the serum level of cytokines mentioned between the two groups. We did not find any association between expression of TRAF-6, miR-146a, and IRAK-1 in PBMCs, and cytokines serum levels with two single nucleotide polymorphisms (SNPs) in miR-146a.

CONCLUSIONS

To the best of writers' knowledge, this research is the first one evaluating the probable link between the miR-146a rs2910164 and rs57095329 variant with miRNAs, relevant cytokine levels, and target genes in brucellosis.

Association analysis of miRNA-related genetic polymorphisms in miR-143/145 and KRAS with colorectal cancer susceptibility and survival

Background: There is accumulating evidence of aberrant expression of miR-143 and miR-145 and their target gene KRAS in colorectal cancer (CRC). We hypothesize that single nucleotide polymorphisms (SNPs) within or near mRNA–microRNA (miRNA) binding sites may affect miRNA/target gene interaction, resulting in differential mRNA/protein expression and promoting the development and progression of CRC. Methods: We conducted a case–control study of 507 patients with CRC recruited from a tertiary hospital and 497 population-based controls to assess the association of genetic polymorphisms in miR-143/145 and the KRAS 3′ untranslated region (3′UTR) with susceptibility to CRC and patients’ survival. In addition, genetic variations of genomic regions located from 500 bp upstream to 500 bp downstream of the miR-143/miR-145 gene and the 3′UTR of KRAS were selected for analysis using the Haploview and HaploReg software. Results: Using publicly available expression profiling data, we found that miR-143/145 and KRAS expression were all reduced in rectal cancer tissue compared with adjacent non-neoplastic large intestinal mucosa. The rs74693964 C/T variant located 65 bp downstream of miR-145 genomic regions was observed to be associated with susceptibility to CRC (adjusted odds ratio (OR): 2.414, 95% CI: 1.385–4.206). Cumulative effects of miR-143 and miR-145 on CRC risk were observed (P_trend=0.03). Patients having CRC carrying variant genotype TT of KRAS rs712 had poorer survival (log-rank P=0.044, adjusted hazard ratio (HR): 4.328, 95% CI: 1.236–15.147). Conclusions: Our results indicate that miRNA-related polymorphisms in miR-143/145 and KRAS are likely to be deleterious and represent potential biomarkers for susceptibility to CRC and patients’ survival.

Potential functions of hsa-miR-155-5p and core genes in chronic myeloid leukemia and emerging role in human cancer: A joint bioinformatics analysis

Considering the critical roles of hsa-miR-155-5p participated in hematopoietic system, this study aims to clarify the possible pathogenesis of chronic myeloid leukemia (CML) induced by hsa-miR-155-5p.Three different strategies were employed, namely a network-based pipeline, a survival analysis and genetic screening method, and a simulation modeling approach, to assess the oncogenic role of hsa-miR-155-5p in CML. We identified new potential roles of hsa-miR-155-5p in CML, involving the BCR/ABL-mediated leukemogenesis through MAPK signaling. Several promising targets including E2F2, KRAS and FLI1 were screened as candidate diagnostic marker genes. The survival analysis revealed that mRNA expression of E2F2, KRAS and FLI1 was negatively correlated with hsa-miR-155-5p and these targets were significantly associated with poor overall survival. Furthermore, an overlap between CML-related genes and hsa-miR-155-5p target genes was revealed using competing endogenous RNA (ceRNA) networks analysis. Taken together, our results reveal the dynamic regulatory aspect of hsa-miR-155-5p as potential player in CML pathogenesis.

The Role of cis- and trans-Acting RNA Regulatory Elements in Leukemia

RNA molecules are a source of phenotypic diversity and an operating system that connects multiple genetic and metabolic processes in the cell. A dysregulated RNA network is a common feature of cancer. Aberrant expression of long non-coding RNA (lncRNA), micro RNA (miRNA), and circular RNA (circRNA) in tumors compared to their normal counterparts, as well as the recurrent mutations in functional regulatory cis-acting RNA motifs have emerged as biomarkers of disease development and progression, opening avenues for the design of novel therapeutic approaches. This review looks at the progress, challenges and future prospects of targeting cis-acting and trans-acting RNA elements for leukemia diagnosis and treatment.

OPTIMIR, a novel algorithm for integrating available genome-wide genotype data into miRNA sequence alignment analysis

Next-generation sequencing is an increasingly popular and efficient approach to characterize the full set of microRNAs (miRNAs) present in human biosamples. MiRNAs' detection and quantification still remain a challenge as they can undergo different posttranscriptional modifications and might harbor genetic variations (polymiRs) that may impact on the alignment step. We present a novel algorithm, OPTIMIR, that incorporates biological knowledge on miRNA editing and genome-wide genotype data available in the processed samples to improve alignment accuracy. OPTIMIR was applied to 391 human plasma samples that had been typed with genome-wide genotyping arrays. OPTIMIR was able to detect genotyping errors, suggested the existence of novel miRNAs and highlighted the allelic imbalance expression of polymiRs in heterozygous carriers. OPTIMIR is written in python, and freely available on the GENMED website (http://www.genmed.fr/index.php/fr/) and on Github (github.com/FlorianThibord/OptimiR).

miR-30a promoter variation contributes to the increased risk of colorectal cancer in an Iranian population

Cytochrome P450 family 24 subfamily A member 1 (CYP24A1) gene is overexpressed in many cancers including colorectal cancer (CRC) and correlated with tumor invasion, lymph node metastasis, and the reduced overall survival. We predicted that miR-30a and miR-125a regulate the CYp24A1 gene expression. Therefore, we performed a case-control study using 800 individuals, including 389 patients with CRC and 411 noncancer controls to evaluate the association between miR-30a rs2222722 and miR-125a rs12976445 polymorphisms, located at in the promoter region, and the risk of sporadic CRC in an Iranian population. The genotyping assay for both polymorphisms was performed using Tetra-primer amplification refractory mutation systems polymerase chain reaction. The results indicated that the frequency of the miR-30a rs2222722 CT genotype was significantly different in the studied groups ( P = 0.0001; odds ratio [OR] = 1.9; 95% confidence interval [CI], 1.39-2.60). Also, a significant difference was observed under the dominant inheritance model ( P = 0.0001; OR = 1.8; 95% CI, 1.33-2.43). The frequency of the miR-30a rs2222722 T allele was significantly associated with increased CRC risk in the studied population ( P = 0.0019; OR = 1.47; 95% CI, 1.15-1.89). Taken together, our study provides preliminary evidence that the rs2222722 polymorphism increases the susceptibility to CRC in an Iranian population. Therefore, the affecting factors on CYP24A1 gene expression such as microRNAs can be considered as risk factors for CRC.

Network Profiling of Brain-Expressed X-Chromosomal MicroRNA Genes Implicates Shared Key MicroRNAs in Intellectual Disability

MicroRNAs are endogenous non-protein-coding RNA molecules that regulate post-transcriptional gene expression. The majority of human miRNAs are brain-expressed and chromosome X is enriched in miRNA genes. We analyzed the genomic regions of 12 brain-expressed pre-miRNAs located on chromosome X coding for 18 mature miRNAs, aiming to investigate the involvement of miRNA sequence variants on X-linked intellectual disability (XLID). Genomic DNA samples from 135 unrelated Brazilian males with intellectual disability, suggestive of X-linked inheritance, were amplified through polymerase chain reaction and sequenced by Sanger sequencing. Although no sequence variations have been identified, suggesting an intense selective pressure, further computational analysis evidenced that eight mature miRNAs (miR-221-3p/222-3p, miR-223-3p, miR-361-5p, miR-362-5p, miR-504-5p.1, miR-505-3p.1, and miR-505-3p.2) act as critical regulators of X-linked and autosomal ID genes in a fully connected network. Enrichment approaches identify transcription regulation, nervous system development, and regulation of cell proliferation as the main common biological processes among the target ID genes. Besides, a clustered chromosomal coverage of the imputed miRNAs target genes and related regulators was found on X chromosome. Considering the role of miRNAs as fine-tuning regulators of gene expression, a systematic analysis of miRNAs' expression could uncover part of the genetic landscape subjacent to ID, being urgently necessary in patients with this condition, particularly XLID.

The rs2910164 variant is associated with reduced miR-146a expression but not cytokine levels in patients with type 2 diabetes

PURPOSE

Previous reports have demonstrated that genetic variations in microRNAs regulome could affect microRNAs-mediated regulation. Therefore, in the present study we were aimed at (1) comparison of microRNA 146-a (miR-146a) peripheral blood mononuclear cells (PBMCs) and plasma levels between diabetic patients and controls, and (2) investigating the possible association of rs2910164 with miR-146a and its related target genes expression and also serum cytokine levels.

METHODS

The study population consisted of 60 subjects including 30 type 2 diabetes (T2D) patients and 30 controls with determined genotypes for rs2910164. The RNA expression levels were determined by real-time PCR. Moreover, TNF-α, IL-6, IL-10 and IL-1β serum levels were measured using ELISA method.

RESULTS

Our results showed that the miR-146a expression levels were significantly decreased in PBMCs (P = 0.004) and plasma (P = 0.008) samples of patients with T2D compared to healthy participants. In addition, we observed that IRAK1 mRNA expression-but not TLR4, TRAF6 and NFĸB-was significantly increased in patients with T2D compared to controls (P = 0.028). The relative expression levels of miR-146a in plasma and PBMCs samples of diabetic patients with the rs2910164 GG genotypes were significantly higher than that in CC (P < 0.05). Moreover, no significant differences were found in miR-146a targets and cytokine levels between the rs2910164 different genotypes.

CONCLUSION

Our study demonstrated that miR-146a circulating levels were significantly elevated in controls compared with T2D patients. In addition, we identified that rs2910164-C allele is associated with reduced expression levels of the miR-146a but not its mRNAs targets and cytokine levels in diabetic patients.

Up-regulation of ribosome biogenesis by MIR196A2 genetic variation promotes endometriosis development and progression

Aberrant miRNA expression has been reported in endometriosis and miRNA gene polymorphisms have been linked to cancer. Because certain ovarian cancers arise from endometriosis, we genotyped seven cancer-related miRNA single nucleotide polymorphisms (MiRSNPs) to investigate their possible roles in endometriosis. Genetic variants in MIR196A2 (rs11614913) and MIR100 (rs1834306) were found to be associated with endometriosis development and related clinical phenotypes, such as infertility and pain. Downstream analysis of the MIR196A2 risk allele revealed upregulation of rRNA editing and protein synthesis genes, suggesting hyper-activation of ribosome biogenesis as a driving force for endometriosis progression. Clinical studies confirmed higher levels of small nucleolar RNAs and ribosomal proteins in atypical endometriosis lesions, and this was more pronounced in the associated ovarian clear cell carcinomas. Treating ovarian clear cells with CX5461, an RNA polymerase I inhibitor, suppressed cell growth and mobility followed by cell cycle arrest at G2/M stage and apoptosis. Our study thus uncovered a novel tumorigenesis pathway triggered by the cancer-related MIR196A2 risk allele during endometriosis development and progression. We suggest that anti-RNA polymerase I therapy may be efficacious for treating endometriosis and associated malignancies.

Comparing miRNA structure of mirtrons and non-mirtrons

BACKGROUND

MicroRNAs proceeds through the different canonical and non-canonical pathways; the most frequent of the non-canonical ones is the splicing-dependent biogenesis of mirtrons. We compare the mirtrons and non-mirtrons of human and mouse to explore how their maturation appears in the precursor structure around the miRNA.

RESULTS

We found the coherence of the overhang lengths what indicates the dependence between the cleavage sites. To explain this dependence we suggest the 2-lever model of the Dicer structure that couples the imprecisions in Drosha and Dicer. Considering the secondary structure of all animal pre-miRNAs we confirmed that single-stranded nucleotides tend to be located near the miRNA boundaries and in its center and are characterized by a higher mutation rate. The 5' end of the canonical 5' miRNA approaches the nearest single-stranded nucleotides what suggests the extension of the loop-counting rule from the Dicer to the Drosha cleavage site. A typical structure of the annotated mirtron pre-miRNAs differs from the canonical pre-miRNA structure and possesses the 1- and 2 nt hanging ends at the hairpin base. Together with the excessive variability of the mirtron Dicer cleavage site (that could be partially explained by guanine at its ends inherited from splicing) this is one more evidence for the 2-lever model. In contrast with the canonical miRNAs the mirtrons have higher snp densities and their pre-miRNAs are inversely associated with diseases. Therefore we supported the view that mirtrons are under positive selection while canonical miRNAs are under negative one and we suggested that mirtrons are an intrinsic source of silencing variability which produces the disease-promoting variants. Finally, we considered the interference of the pre-miRNA structure and the U2snRNA:pre-mRNA basepairing. We analyzed the location of the branchpoints and found that mirtron structure tends to expose the branchpoint site what suggests that the mirtrons can readily evolve from occasional hairpins in the immediate neighbourhood of the 3' splice site.

CONCLUSION

The miRNA biogenesis manifests itself in the footprints of the secondary structure. Close inspection of these structural properties can help to uncover new pathways of miRNA biogenesis and to refine the known miRNA data, in particular, new non-canonical miRNAs may be predicted or the known miRNAs can be re-classified.

High-Quality Overlapping Paired-End Reads for the Detection of A-to-I Editing on Small RNA

Paired-end RNA sequencing (RNA-seq) is usually applied to the quantification of long transcripts such as messenger or long non-coding RNAs, in which case overlapping pairs are discarded. In contrast, RNA-seq on short RNAs (≤ 200 nt) is typically carried out in single-end mode, as the additional cost associated with paired-end would only translate into redundant sequence information. Here, we exploit paired-end sequencing of short RNAs as a strategy to filter out sequencing errors and apply this method to the identification of adenosine-to-inosine (A-to-I) RNA editing events on human precursor microRNA (pre-miRNA) and mature miRNA. Combined with RNA immunoprecipitation sequencing (RIP-seq) of A-to-I RNA editing enzymes, this method takes full advantage of deep sequencing technology to identify RNA editing sites with unprecedented resolution in terms of editing efficiency.

A comprehensive review of web-based non-coding RNA resources for cancer research

Non-coding RNAs include many kinds of RNAs that did not encode proteins. Recent evidences reveal that ncRNAs play critical roles in initiation and progression of cancers. But it is not easy for cancer biologists and medical doctors to easily know the potential roles of ncRNAs in cancer and retrieve the information of ncRNAs under their investigations. To make the available web-based resources more accessible and understandable, we made a comprehensive review for 49 web-based resources of three types of ncRNAs, i.e., microRNAs (miRNAs), long non-coding RNAs (lncRNAs) and circular RNAs (circRNAs). We also listed some preferred resources for 6 different types of analyses related to ncRNAs.

SNPs in bone-related miRNAs are associated with the osteoporotic phenotype

Biogenesis and function of microRNAs can be influenced by genetic variants in the pri-miRNA sequences leading to phenotypic variability. This study aims to identify single nucleotide polymorphisms (SNPs) affecting the expression levels of bone-related mature microRNAs and thus, triggering an osteoporotic phenotype. An association analysis of SNPs located in pri-miRNA sequences with bone mineral density (BMD) was performed in the OSTEOMED2 cohort (n = 2183). Functional studies were performed for assessing the role of BMD-associated miRNAs in bone cells. Two SNPs, rs6430498 in the miR-3679 and rs12512664 in the miR-4274, were significantly associated with femoral neck BMD. Further, we measured these BMD-associated microRNAs in trabecular bone from osteoporotic hip fractures comparing to non-osteoporotic bone by qPCR. Both microRNAs were found overexpressed in fractured bone. Increased matrix mineralization was observed after miR-3679-3p inhibition in human osteoblastic cells. Finally, genotypes of rs6430498 and rs12512664 were correlated with expression levels of miR-3679 and miR-4274, respectively, in osteoblasts. In both cases, the allele that generated higher microRNA expression levels was associated with lower BMD values. In conclusion, two osteoblast-expressed microRNAs, miR-3679 and miR-4274, were associated with BMD; their overexpression could contribute to the osteoporotic phenotype. These findings open new areas for the study of bone disorders.

MiRNAs and miRNA Polymorphisms Modify Drug Response

Differences in expression of drug response-related genes contribute to inter-individual variation in drugs’ biological effects. MicroRNAs (miRNAs) are small noncoding RNAs emerging as new players in epigenetic regulation of gene expression at post-transcriptional level. MiRNAs regulate the expression of genes involved in drug metabolism, drug transportation, drug targets and downstream signal molecules directly or indirectly. MiRNA polymorphisms, the genetic variations affecting miRNA expression and/or miRNA-mRNA interaction, provide a new insight into the understanding of inter-individual difference in drug response. Here, we provide an overview of the recent progress in miRNAs mediated regulation of biotransformation enzymes, drug transporters, and nuclear receptors. We also describe the implications of miRNA polymorphisms in cancer chemotherapy response.

Accurate detection for a wide range of mutation and editing sites of microRNAs from small RNA high-throughput sequencing profiles

Various types of mutation and editing (M/E) events in microRNAs (miRNAs) can change the stabilities of pre-miRNAs and/or complementarities between miRNAs and their targets. Small RNA (sRNA) high-throughput sequencing (HTS) profiles can contain many mutated and edited miRNAs. Systematic detection of miRNA mutation and editing sites from the huge volume of sRNA HTS profiles is computationally difficult, as high sensitivity and low false positive rate (FPR) are both required. We propose a novel method (named MiRME) for an accurate and fast detection of miRNA M/E sites using a progressive sequence alignment approach which refines sensitivity and improves FPR step-by-step. From 70 sRNA HTS profiles with over 1.3 billion reads, MiRME has detected thousands of statistically significant M/E sites, including 3'-editing sites, 57 A-to-I editing sites (of which 32 are novel), as well as some putative non-canonical editing sites. We demonstrated that a few non-canonical editing sites were not resulted from mutations in genome by integrating the analysis of genome HTS profiles of two human cell lines, suggesting the existence of new editing types to further diversify the functions of miRNAs. Compared with six existing studies or methods, MiRME has shown much superior performance for the identification and visualization of the M/E sites of miRNAs from the ever-increasing sRNA HTS profiles.

MicroRNA hsa-miR-4717-5p regulates RGS2 and may be a risk factor for anxiety-related traits

Regulator of G-protein Signaling 2 (RGS2) is a key regulator of G-protein-coupled signaling pathways involved in fear and anxiety. Data from rodent models and genetic analysis of anxiety-related traits and disorders in humans suggest down-regulation of RGS2 expression to be a risk factor for anxiety. Here we investigated, whether genetic variation in microRNAs mediating posttranscriptional down-regulation of RGS2 may be a risk factor for anxiety as well. 75 microRNAs predicted to regulate RGS2 were identified by four bioinformatic algorithms and validated experimentally by luciferase reporter gene assays. Specificity was confirmed for six microRNAs (hsa-miR-1271-5p, hsa-miR-22-3p, hsa-miR-3591-3p, hsa-miR-377-3p, hsa-miR-4717-5p, hsa-miR-96-5p) by disrupting their seed sequence at the 3' untranslated region of RGS2. Hsa-miR-4717-5p showed the most robust effect on RGS2 and regulated two other candidate genes of anxiety disorders (CNR1 and IKBKE) as well. Two SNPs (rs150925, rs161427) within and 1,000 bp upstream of the hostgene of hsa-miR-4717-5p (MIR4717) show a minor allele frequency greater than 0.05. Both were in high linkage disequilibrium (r(2) = 1, D' = 1) and both major (G) alleles showed a trend for association with panic disorder with comorbid agoraphobia in one of two patient/control samples (combined n(patients) = 497). Dimensional anxiety traits, as described by Anxiety Sensitivity Index (ASI) and Agoraphobic Cognitions Questionnaire (ACQ) were significantly higher among carriers of both major (G) alleles in a combined patient/control sample (n(combined) = 831). Taken together, data indicate that MIR4717 regulates human RGS2 and contributes to the genetic risk towards anxiety-related traits.

Genetic variants in microRNA genes: impact on microRNA expression, function, and disease

MicroRNAs (miRNAs) are important regulators of gene expression and like any other gene, their coding sequences are subject to genetic variation. Variants in miRNA genes can have profound effects on miRNA functionality at all levels, including miRNA transcription, maturation, and target specificity, and as such they can also contribute to disease. The impact of variants in miRNA genes is the focus of the present review. To put these effects into context, we first discuss the requirements of miRNA transcripts for maturation. In the last part an overview of available databases and tools and experimental approaches to investigate miRNA variants related to human disease is presented.

New insights about miRNAs in cystic fibrosis

The molecular basis of cystic fibrosis (CF) is a mutation-related defect in the epithelial-cell chloride channel called CF transmembrane conductance regulator (CFTR). This defect alters chloride ion transport and impairs water transport across the cell membrane. Marked clinical heterogeneity occurs even among patients carrying the same mutation in the CFTR gene. Recent studies suggest that such heterogeneity could be related to epigenetic factors and/or miRNAs, which are small noncoding RNAs that modulate the expression of various proteins via post-transcriptional inhibition of gene expression. In the respiratory system, it has been shown that the dysregulation of miRNAs could participate in and lead to pathogenicity in several diseases. In CF airways, recent studies have proposed that miRNAs may modulate disease progression by affecting the production of either CFTR or various proteins that are dysregulated in the CF lung. Herein, we provide an overview of studies showing how miRNAs may modulate CF pathology and the efforts to develop miRNA-based treatments and/or to consider miRNAs as biomarkers. The identification of miRNAs involved in CF disease progression opens up new avenues toward treatments targeting selected clinical components of CF, independently from the CFTR mutation.

Revealing editing and SNPs of microRNAs in colon tissues by analyzing high-throughput sequencing profiles of small RNAs

Background

Editing and mutations in microRNAs (miRNAs) can change the stability of pre-miRNAs and/or complementarities between miRNAs and their targets. Small RNA (sRNA) high-throughput sequencing (HTS) profiles contain miRNAs that are originated from mutated DNAs or are edited during their biogenesis procedures. It is largely unknown whether miRNAs are edited in colon tissues since existing studies mainly focused their attention on the editing of miRNAs in brain tissues.

Results

Through comprehensive analysis of four high-throughput sequencing profiles of normal and cancerous colon tissues, we identified 548 editing and/or SNPs in miRNAs that are significant in at least one of the sequencing profiles used. Our results show that the most abundant editing events of miRNAs in colon tissues are 3'-A and 3'-U. In addition to four known A-to-I editing sites previously reported in brain tissues, four novel A-to-I editing sites are also identified in colon tissues.

Conclusions

This suggests that A-to-I editing of miRNAs potentially is a commonly existing mechanism in different tissues to diversify the possible functional roles of miRNAs, but only a small portion of different miRNAs are edited by the A-to-I mechanism at a significant level. Our results suggest that there are other types of editing in miRNAs through unknown mechanisms. Furthermore, several SNPs in miRNAs are also identified.

Heterogeneity and individuality: microRNAs in mental disorders

MicroRNAs are about 22 nucleotide long single-stranded RNA molecules, negatively regulating gene expression of a single gene or a gene network. In neural tissues, they have been implicated in developmental and neuroplasticity-related processes, such as neurogenesis, differentiation, apoptosis and long-term potentiation. Their molecular mode of action is reminiscent of findings of genome-wide association studies in mental disorders, unable to attribute the risk of disease to a specific gene, but rather to multiple genes, gene-networks and gene-environment interaction. As such, microRNAs are an attractive target for research. Here, we review clinical studies conducted in humans on microRNAs in mental disorders with a particular focus on schizophrenia, bipolar disorder, major depressive disorder and anxiety disorders. The majority of clinical studies have focused on schizophrenia. The most robust finding has been reported for rs1625579 located in MIR137HG, which was associated with schizophrenia on a genome-wide level. Concerning bipolar disorder, major depression and anxiety disorders, promising results have been published, but only a considerably smaller number of clinical studies is available and genome-wide association studies did not suggest a direct link to microRNAs so far. Expression of microRNAs as biomarkers of mental disorders and treatment response is currently emerging with preliminary results. Larger-scaled genetic and functional studies along with translational research are needed to enhance our understanding of microRNAs in mental disorders. These studies will aid in disentangling the complex genetic nature of these disorders and possibly contribute to the development of novel, individualized diagnostic and therapeutic approaches.

Micro spies from the brain to the periphery: new clues from studies on microRNAs in neuropsychiatric disorders

microRNAs (miRNAs) are small non-coding RNAs (20-22 nucleotides) playing a major role in post-transcriptional regulation of gene expression. miRNAs are predicted to regulate more than 50% of all the protein-coding genes. Increasing evidence indicates that they may play key roles in the biological pathways that regulate neurogenesis and synaptic plasticity, as well as in neurotransmitter homeostasis in the adult brain. In this article we review recent studies suggesting that miRNAs may be involved in the pathophysiology of neuropsychiatric disorders and in the action of psychotropic drugs, in particular by analyzing the contribution of genomic studies in patients' peripheral tissues. Alterations in miRNA expression have been observed in schizophrenia, bipolar disorder, major depression, Parkinson's disease, Alzheimer's disease and other neuropsychiatric conditions. In particular, intriguing findings concern the identification of disease-associated miRNA signatures in peripheral tissues, or modifications in miRNA profiles induced by drug treatments. Furthermore, genetic variations in miRNA sequences and miRNA-related genes have been described in neuropsychiatric diseases. Overall, though still at a preliminary stage, several lines of evidence indicate an involvement of miRNAs in both the pathophysiology and pharmacotherapy of neuropsychiatric disorders. In this regard, the data obtained in peripheral tissues may provide further insights into the etiopathogenesis of several brain diseases and contribute to identify new biomarkers for diagnostic assessment improvement and treatment personalization.