SNPs in miRNAs and Target Sequences: Role in Cancer and Diabetes

The Role of Single Nucleotide Polymorphisms in MicroRNA Genes in Head and Neck Squamous Cell Carcinomas: Susceptibility and Prognosis

BACKGROUND

Head and neck squamous cell carcinoma (HNSCC) is one of the most prevalent cancers worldwide. The identification of molecular alterations adding to the individual risk of HNSCC development and progression is one of the most important challenges in studies on cancer genetics. MicroRNAs (miRNAs), which belong to the group of important post-transcriptional regulators of human gene expression, seem to be valuable options for consideration as key modifiers of individual cancer risk, and therefore may be helpful in predicting inter-individual differences in cancer risk, response to treatment and prognosis.

METHODS

There have not been many studies focused on the relationship between miRNA variants and HNSCC published in PubMed within the last 15 years. We found and analyzed 30 reviews, meta-analyses and research papers and revealed 14 SNPs which have been reported as significant in the context of HNSCC susceptibility and/or prognosis.

RESULTS

These 14 SNPs were located in 13 separate miRNAs. Among them, four were the most frequently studied (miRNA-146, -196, -149 and -499) and have been shown to have the greatest impact on the course of HNSCC. However, the presented results have been conflicting.

CONCLUSIONS

It must be concluded that, despite the years of studies, there are no conclusive reports demonstrating a significant role of SNPs in miRNAs in the context of the susceptibility to HNSCC or its prognosis.

Genetic polymorphism in untranslated regions of PRKCZ influences mRNA structure, stability and binding sites

BACKGROUND

Variations in untranslated regions (UTR) alter regulatory pathways impacting phenotype, disease onset, and course of disease. Protein kinase C Zeta (PRKCZ), a serine-threonine kinase, is implicated in cardiovascular, neurological and oncological disorders. Due to limited research on PRKCZ, this study aimed to investigate the impact of UTR genetic variants' on binding sites for transcription factors and miRNA. RNA secondary structure, eQTLs, and variation tolerance analysis were also part of the study.

METHODS

The data related to PRKCZ gene variants was downloaded from the Ensembl genome browser, COSMIC and gnomAD. The RegulomeDB database was used to assess the functional impact of 5' UTR and 3'UTR variants. The analysis of the transcription binding sites (TFBS) was done through the Alibaba tool, and the Kyoto Encyclopaedia of Genes and Genomes (KEGG) was employed to identify pathways associated with PRKCZ. To predict the effect of variants on microRNA binding sites, PolymiRTS was utilized for 3' UTR variants, and the SNPinfo tool was used for 5' UTR variants.

RESULTS

The results obtained indicated that a total of 24 variants present in the 3' UTR and 25 variants present in the 5' UTR were most detrimental. TFBS analysis revealed that 5' UTR variants added YY1, repressor, and Oct1, whereas 3' UTR variants added AP-2alpha, AhR, Da, GR, and USF binding sites. The study predicted TFs that influenced PRKCZ expression. RNA secondary structure analysis showed that eight 5' UTR and six 3' UTR altered the RNA structure by either removal or addition of the stem-loop. The microRNA binding site analysis highlighted that seven 3' UTR and one 5' UTR variant altered the conserved site and also created new binding sites. eQTLs analysis showed that one variant was associated with PRKCZ expression in the lung and thyroid. The variation tolerance analysis revealed that PRKCZ was an intolerant gene.

CONCLUSION

This study laid the groundwork for future studies aimed at targeting PRKCZ as a therapeutic target.

Amplification-Free Strategy for miRNA Quantification in Human Serum Using Single Particle ICP-MS and Gold Nanoparticles as Labels

MicroRNAs (miRNAs), which are short single-stranded RNA sequences between 18 and 24 nucleotides, are known to play a crucial role in gene expression. Changes in their expression are not only involved in many diseases but also as a response to physiological changes, such as physical exercise. In this work, a new analytical strategy for the sensitive and specific analysis of miRNA sequences in human plasma is presented. The developed strategy does not depend on any nucleic acid amplification process and can be obtained in direct correlation to the number of events obtained by using single-particle ICP-MS measurements. The high selectivity of the assay (up to single nucleotide polymorphisms) can be achieved by a double hybridization process of the target miRNA with a complementary capture oligonucleotide that is conjugated to a magnetic microparticle and simultaneously with a complementary reporter oligonucleotide conjugated to a gold nanoparticle. Thanks to the novel approach followed in this method, the stoichiometry of the oligonucleotide-nanoparticle conjugates does not need to be addressed for the quantification of the target miRNA, which also represents a big advantage over other similar methods. The optimized method is applied to the determination of a miRNA as a biomarker of physical exercise in non-spiked human serum samples, and the results are validated against rt-qPCR. The achieved sensitivity permits the direct differentiation among sedentary and sportive subjects. This general platform can be easily applied to any other sequence by only modifying the capture and reporter oligonucleotides, paving the way for multiple clinically interesting applications.

Methods for Functional Characterization of Genetic Polymorphisms of Non-Coding Regulatory Regions of the Human Genome

Currently, numerous associations between genetic polymorphisms and various diseases have been characterized through the Genome-Wide Association Studies. Majority of the clinically significant polymorphisms are localized in non-coding regions of the genome. While modern bioinformatic resources make it possible to predict molecular mechanisms that explain influence of the non-coding polymorphisms on gene expression, such hypotheses require experimental verification. This review discusses the methods for elucidating molecular mechanisms underlying dependence of the disease pathogenesis on specific genetic variants within the non-coding sequences. A particular focus is on the methods for identification of transcription factors with binding efficiency dependent on polymorphic variations. Despite remarkable progress in bioinformatic resources enabling prediction of the impact of polymorphisms on the disease pathogenesis, there is still the need for experimental approaches to investigate this issue.

Whole transcriptome sequencing analyses of islets reveal ncRNA regulatory networks underlying impaired insulin secretion and increased β-cell mass in high fat diet-induced diabetes mellitus

AIM

Our study aims to identify novel non-coding RNA-mRNA regulatory networks associated with β-cell dysfunction and compensatory responses in obesity-related diabetes.

METHODS

Glucose metabolism, islet architecture and secretion, and insulin sensitivity were characterized in C57BL/6J mice fed on a 60% high-fat diet (HFD) or control for 24 weeks. Islets were isolated for whole transcriptome sequencing to identify differentially expressed (DE) mRNAs, miRNAs, IncRNAs, and circRNAs. Regulatory networks involving miRNA-mRNA, lncRNA-mRNA, and lncRNA-miRNA-mRNA were constructed and functions were assessed through Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses.

RESULTS

Despite compensatory hyperinsulinemia and a significant increase in β-cell mass with a slow rate of proliferation, HFD mice exhibited impaired glucose tolerance. In isolated islets, insulin secretion in response to glucose and palmitic acid deteriorated after 24 weeks of HFD. Whole transcriptomic sequencing identified a total of 1324 DE mRNAs, 14 DE miRNAs, 179 DE lncRNAs, and 680 DE circRNAs. Our transcriptomic dataset unveiled several core regulatory axes involved in the impaired insulin secretion in HFD mice, such as miR-6948-5p/Cacna1c, miR-6964-3p/Cacna1b, miR-3572-5p/Hk2, miR-3572-5p/Cckar and miR-677-5p/Camk2d. Additionally, proliferative and apoptotic targets, including miR-216a-3p/FKBP5, miR-670-3p/Foxo3, miR-677-5p/RIPK1, miR-802-3p/Smad2 and ENSMUST00000176781/Caspase9 possibly contribute to the increased β-cell mass in HFD islets. Furthermore, competing endogenous RNAs (ceRNA) regulatory network involving 7 DE miRNAs, 15 DE lncRNAs and 38 DE mRNAs might also participate in the development of HFD-induced diabetes.

CONCLUSIONS

The comprehensive whole transcriptomic sequencing revealed novel non-coding RNA-mRNA regulatory networks associated with impaired insulin secretion and increased β-cell mass in obesity-related diabetes.

The Search for and Functional Analysis of Genetic Variants in microRNA-Binding Sites using Massively Parallel Reporter Assay

A large-scale search for the genetic variants with a bias in the representation of alleles in transcriptome data (AE SNPs) and the binding sites in microRNA 3'-UTRs was performed and their functional significance was assessed using massively parallel reporter assay (MPRA). Of the 629,559 associated "SNP-gene" pairs (eQTLs) discovered in the human liver tissue according to the GTEx Analysis V8 data, 4394 polymorphic positions in the 3'-UTRs of the genes, which represent the eQTLs for these genes were selected. The TargetScanHuman 7.0 algorithm and PolymiRTS database were searched for the potential microRNA-binding sites. Of the predicted microRNA sites affected by eQTL-SNPs, we selected 51 sites with the best evidence of functionality according to Ago2-CLIP-seq, CLEAR-CLIP, and eCLIP-seq for RNA-binding proteins. For MPRA, a library of the plasmids carrying the main and alternative alleles for each AE SNP (in total, 102 constructs) was created. Allele-specific expression for 6 SNPs was detected by transfection of the HepG2 cell line with the constructed plasmid library and sequencing of target DNA and RNA sequences using the Illumina (MiSeq) platform.

Genetic variations in tumor-suppressor miRNA-encoding genes and their target genes: focus on breast cancer development and possible therapeutic strategies

MicroRNAs (miRNAs) negatively affect gene expression by binding to their specific mRNAs resulting in either mRNA destruction or translational repression. The aberrant expression of various miRNAs has been associated with a number of human cancer. Oncogenic or tumor-suppressor miRNAs regulate a variety of pathways involved in the development of breast cancer (BC), including cell proliferation, apoptosis, metastasis, cancer recurrence, and chemoresistance. Variations in miRNA-encoding genes and their target genes lead to dysregulated gene expression resulting in the development and progression of BC. The various therapeutic approaches to treat the disease include chemotherapy, radiation therapy, surgical removal, hormone therapy, chemotherapy, and targeted biological therapy. The purpose of the current review is to explore the genetic variations in tumor-suppressor miRNA-encoding genes and their target genes in association with the disease development and prognosis. The therapeutic interventions targeting the variants for better disease outcomes have also been discussed.

Increased hsa-miR-100-5p Expression Improves Hepatocellular Carcinoma Prognosis in the Asian Population with PLK1 Variant rs27770A>G

Hepatocellular carcinoma (HCC) has the highest incidence and mortality in the Asian population, and race is an independent risk factor affecting survival time in liver cancer. Micro RNAs (miRNAs) are remarkably dysregulated in HCC and closely associated with HCC prognosis. Recent studies show that genetic variability between ethnic groups may result in differences in the specificity of HCC miRNA biomarkers. Here, we reveal a high expression level of hsa-miR-100-5p, an HCC prognosis-related miRNA, which improves HCC prognosis in the Asian Population with Polo-like kinase 1 (PLK1) variant rs27770A>G. In this study, we discovered that hsa-miR-100-5p was downregulated in various HCC cell lines. While mimics transient transfection and mouse liver cancer model confirmed the interaction between hsa-miR-100-5p and PLK1, a stratified analysis based on the Cancer Genome Atlas Liver Hepatocellular Carcinoma (TCGA-LIHC) data suggest both low hsa-miR-100-5p expression level and high PLK1 expression level associated with poor HCC prognosis, especially in the Asian population. According to the 1000 Genomes Project database, the SNP rs27770 located in 3'UTR of PLK1 had a significantly higher G allele frequency in the East Asian population. Bioinformatics analysis suggested that rs27770 A>G affects PLK1 mRNA secondary structure and alters the hsa-miR-100-5p/PLK1 interaction by forming an additional seedless binding site. This racial variation caused PLK1 to be more vulnerable to hsa-miR-100-5p inhibition, resulting in hsa-miR-100-5p being more favorable for HCC prognosis in the Asian population.

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.

Circ_0004535/miR-1827/CASP8 network involved in type 2 diabetes mellitus with nonalcoholic fatty liver disease

Diagnostic delay in type 2 diabetes mellitus (T2DM) with nonalcoholic fatty liver disease (NAFLD) patients often leads to a serious public health problem. Understanding the pathophysiological mechanisms of disease will help develop more effective treatments. High-throughput sequencing was used to determine the expression levels of circRNAs, and mRNAs in health controls, T2DM patients, and T2DM with NAFLD patients. Differentially expressed genes (DEcircRs, DEmRs) in T2DM with NAFLD were identified by differential analysis. The miRNAs with targeted relationship with the DEcircRs and DEmRs were respectively predicted to construct a ceRNA regulatory network. In addition, enrichment analysis of DEmRs in the ceRNA network was performed. The expression of important DEcircRs was further validated by quantitative real-time PCR (qRT-PCR). The steatosis was detected in glucose treated LO2 cells by overexpressing circ_0004535, and CASP8. There were 586 DEmRs, and 10 DEcircRs in both T2DM and T2DM with NAFLD patients. Combined with predicted results and differential analysis, the ceRNA networks were constructed. The DEmRs in the ceRNA networks were mainly enriched in Toll-like receptor signaling pathway, and apoptosis. Importantly, dual luciferase experiments validated the targeted binding of hsa_circ_0004535 and hsa-miR-1827 or hsa-miR-1827 and CASP8. qRT-PCR experiments validated that hsa_circ_0004535, and CASP8 was downregulated and hsa-miR-1827 was upregulated expression in peripheral blood of T2DM with NAFLD patients. Abnormal cell morphology, and increased lipid droplet fusion were observed in the glucose treated LO2 cells, overexpression of circ_0004535 and CASP8 ameliorated these changes. Our work provides a deeper understanding of ceRNA mediated pathogenesis of T2DM with NAFLD and provides a novel strategy for treatment.

Computational analysis of five neurodegenerative diseases reveals shared and specific genetic loci

Neurodegenerative diseases (ND) are heterogeneous disorders of the central nervous system that share a chronic and selective process of neuronal cell death. A computational approach to investigate shared genetic and specific loci was applied to 5 different ND: Amyotrophic lateral sclerosis (ALS), Alzheimer's disease (AD), Parkinson's disease (PD), Multiple sclerosis (MS), and Lewy body dementia (LBD). The datasets were analyzed separately, and then we compared the obtained results. For this purpose, we applied a genetic correlation analysis to genome-wide association datasets and revealed different genetic correlations with several human traits and diseases. In addition, a clumping analysis was carried out to identify SNPs genetically associated with each disease. We found 27 SNPs in AD, 6 SNPs in ALS, 10 SNPs in PD, 17 SNPs in MS, and 3 SNPs in LBD. Most of them are located in non-coding regions, with the exception of 5 SNPs on which a protein structure and stability prediction was performed to verify their impact on disease. Furthermore, an analysis of the differentially expressed miRNAs of the 5 examined pathologies was performed to reveal regulatory mechanisms that could involve genes associated with selected SNPs. In conclusion, the results obtained constitute an important step toward the discovery of diagnostic biomarkers and a better understanding of the diseases.

miR-499a rs3746444 A>G Polymorphism Is Correlated with Type 2 Diabetes Mellitus and Diabetic Polyneuropathy in a Romanian Cohort: A Preliminary Study

Type 2 diabetes mellitus (T2DM) is a common metabolic disorder that results from complex interactions of both environmental and genetic factors. Many single nucleotide polymorphisms (SNPs), including noncoding RNA genes, have been investigated for their association with susceptibility to T2DM and its complications, with little evidence available regarding Caucasians. The aim of the present study was to establish whether four miRNA SNPs (miR-27a rs895819 T>C, miR-146a rs2910164 G>C, miR-196a2 rs11614913 C>T, and miR-499a rs3746444 A>G) are correlated with susceptibility to T2DM and/or diabetic polyneuropathy (DPN) in a Romanian population. A total of 167 adult T2DM patients and 324 age- and sex-matched healthy controls were included in our study. miRNA SNPs were detected by real-time PCR using a TaqMan genotyping assay. A significant association with T2DM was observed only for the miR-499a rs3746444 A>G SNP in all the tested models, and the frequencies of both the miR-499a rs3746444 AG and the GG genotypes were higher in the T2DM patients compared to the controls. No correlation was observed for the miR-27a rs895819 T>C, miR-146a rs2910164 G>C, or miR-196a2 rs11614913 C>T SNPs in any genetic model. When we assessed the association of these SNPs with DPN separately, we found a positive association for the miR-499a rs3746444 SNP in both codominant and dominant models (OR 6.47, 95% CI: 1.71-24.47; OR 2.30, 95% CI: 1.23-4.29, respectively). In conclusion, this study shows that miR-499a rs3746444 A>G may influence both T2DM and DPN susceptibility, with carriers of the GG genotype and the G allele being at an increased risk in the Romanian population.

SNPs in 3'UTR miRNA Target Sequences Associated with Individual Drug Susceptibility

The complementary interaction of microRNAs (miRNAs) with their binding sites in the 3'untranslated regions (3'UTRs) of target gene mRNAs represses translation, playing a leading role in gene expression control. MiRNA recognition elements (MREs) in the 3'UTRs of genes often contain single nucleotide polymorphisms (SNPs), which can change the binding affinity for target miRNAs leading to dysregulated gene expression. Accumulated data suggest that these SNPs can be associated with various human pathologies (cancer, diabetes, neuropsychiatric disorders, and cardiovascular diseases) by disturbing the interaction of miRNAs with their MREs located in mRNA 3'UTRs. Numerous data show the role of SNPs in 3'UTR MREs in individual drug susceptibility and drug resistance mechanisms. In this review, we brief the data on such SNPs focusing on the most rigorously proven cases. Some SNPs belong to conventional genes from the drug-metabolizing system (in particular, the genes coding for cytochromes P450 (CYP 450), phase II enzymes (SULT1A1 and UGT1A), and ABCB3 transporter and their expression regulators (PXR and GATA4)). Other examples of SNPs are related to the genes involved in DNA repair, RNA editing, and specific drug metabolisms. We discuss the gene-by-gene studies and genome-wide approaches utilized or potentially utilizable to detect the MRE SNPs associated with individual response to drugs.

IsomiR-eQTL: A Cancer-Specific Expression Quantitative Trait Loci Database of miRNAs and Their Isoforms

The identification of expression quantitative trait loci (eQTL) is an important component in efforts to understand how genetic variants influence disease risk. MicroRNAs (miRNAs) are short noncoding RNA molecules capable of regulating the expression of several genes simultaneously. Recently, several novel isomers of miRNAs (isomiRs) that differ slightly in length and sequence composition compared to their canonical miRNAs have been reported. Here we present isomiR-eQTL, a user-friendly database designed to help researchers find single nucleotide polymorphisms (SNPs) that can impact miRNA (miR-eQTL) and isomiR expression (isomiR-eQTL) in 30 cancer types. The isomiR-eQTL includes a total of 152,671 miR-eQTLs and 2,390,805 isomiR-eQTLs at a false discovery rate (FDR) of 0.05. It also includes 65,733 miR-eQTLs overlapping known cancer-associated loci identified through genome-wide association studies (GWAS). To the best of our knowledge, this is the first study investigating the impact of SNPs on isomiR expression at the genome-wide level. This database may pave the way for researchers toward finding a model for personalised medicine in which miRNAs, isomiRs, and genotypes are utilised.

Mutational Signature and Integrative Genomic Analysis of Human Papillomavirus-Associated Penile Squamous Cell Carcinomas from Latin American Patients

Simple Summary

DNA sequencing has been crucial to comprehending cancer mutational patterns, leading to the identification of driver genes and altered signaling pathways. Thus, identifying new pathogenic variants and their impact on tumor onset, progression, and treatment response has fueled tumor biology research. Here, we present novel findings addressing the first whole-exome sequencing (WES) of human papillomavirus (HPV)-associated penile squamous cell carcinoma (PSCC) from Latin Americans and its association with pathogenesis. We also compared the molecular profile of the tumors to that of three previous studies from populations with different genetic and socioeconomic backgrounds, the majority of which was HPV-negative. We describe the most altered genes and the main pathogenic variants found in the Latin Americans, ten of which are exclusive to our study sample. The data allowed us to identify molecular pathways and druggable targets with potential treatment value for this still-neglected HPV-associated carcinoma.

Abstract

High-throughput DNA sequencing has allowed for the identification of genomic alterations and their impact on tumor development, progression, and therapeutic responses. In PSCC, for which the incidence has progressively increased worldwide, there are still limited data on the molecular mechanisms involved in the disease pathogenesis. In this study, we characterized the mutational signature of 30 human papillomavirus (HPV)-associated PSCC cases from Latin Americans, using whole-exome sequencing. Copy number variations (CNVs) were also identified and compared to previous array-generated data. Enrichment analyses were performed to reveal disrupted pathways and to identify alterations mapped to HPV integration sites (HPVis) and miRNA–mRNA hybridization regions. Among the most frequently mutated genes were NOTCH1, TERT, TTN, FAT1, TP53, CDKN2A, RYR2, CASP8, FBXW7, HMCN2, and ITGA8. Of note, 92% of these altered genes were localized at HPVis. We also found mutations in ten novel genes (KMT2C, SMARCA4, PTPRB, AJUBA, CR1, KMT2D, NBEA, FAM135B, GTF2I, and CIC), thus increasing our understanding of the potential HPV-disrupted pathways. Therefore, our study reveals innovative targets with potential therapeutic benefits for HPV-associated PSCCs. The CNV analysis by sequencing (CNV-seq) revealed five cancer-associated genes as the most frequent with gains (NOTCH1, MYC, NUMA1, PLAG1, and RAD21), while 30% of the tumors showed SMARCA4 with loss. Additionally, four cancer-associated genes (CARD11, CSMD3, KDR, and TLX3) carried untranslated regions (UTRs) variants, which may impact gene regulation by affecting the miRNAs hybridization regions. Altogether, these data contribute to the characterization of the mutational spectrum and its impact on cellular signaling pathways in PSCC, thus reinforcing the pivotal role of HPV infection in the molecular pathogenesis of these tumors.