Identification of RNA-binding protein hnRNP C targeting the 3'UTR of the TAP-associated glycoprotein tapasin in melanoma.

Deregulation or loss of the human leukocyte antigen class I (HLA-I) molecules on tumor cells leading to inhibition of CD8+ T cell recognition is an important tumor immune escape strategy, which could be caused by a posttranscriptional control of molecules in the HLA-I pathway mediated by RNA-binding proteins (RBPs). So far, there exists only limited information about the interaction of RBPs with HLA-I-associated molecules, but own work demonstrated a binding of the heterogeneous ribonucleoprotein C (hnRNP C) to the 3' untranslated region (UTR) of the TAP-associated glycoprotein tapasin (tpn). In this study, in silico analysis of pan-cancer TCGA datasets revealed that hnRNP C is higher expressed in tumor specimens compared to corresponding normal tissues, which is negatively correlated to tpn expression, T cell infiltration and the overall survival of tumor patients. Functional analysis demonstrated an upregulation of tpn expression upon siRNA-mediated downregulation of hnRNP C, which is accompanied by an increased HLA-I surface expression. Thus, hnRNP C has been identified to target tpn and its inhibition could improve the HLA-I surface expression on melanoma cells suggesting its use as a possible biomarker for T-cell-based tumor immunotherapies.

Ckip-1 3'UTR alleviates prolonged sleep deprivation induced cardiac dysfunction by activating CaMKK2/AMPK/cTNI pathway.

Sleep deprivation (SD) has emerged as a critical concern impacting human health, leading to significant damage to the cardiovascular system. However, the underlying mechanisms are still unclear, and the development of targeted drugs is lagging. Here, we used mice to explore the effects of prolonged SD on cardiac structure and function. Echocardiography analysis revealed that cardiac function was significantly decreased in mice after five weeks of SD. Real-time quantitative PCR (RT-q-PCR) and Masson staining analysis showed that cardiac remodeling marker gene Anp (atrial natriuretic peptide) and fibrosis were increased, Elisa assay of serum showed that the levels of creatine kinase (CK), creatine kinase-MB (CK-MB), ANP, brain natriuretic peptide (BNP) and cardiac troponin T (cTn-T) were increased after SD, suggesting that cardiac remodeling and injury occurred. Transcript sequencing analysis indicated that genes involved in the regulation of calcium signaling pathway, dilated cardiomyopathy, and cardiac muscle contraction were changed after SD. Accordingly, Western blotting analysis demonstrated that the cardiac-contraction associated CaMKK2/AMPK/cTNI pathway was inhibited. Since our preliminary research has confirmed the vital role of Casein Kinase-2 -Interacting Protein-1 (CKIP-1, also known as PLEKHO1) in cardiac remodeling regulation. Here, we found the levels of the 3' untranslated region of Ckip-1 (Ckip-1 3'UTR) decreased, while the coding sequence of Ckip-1 (Ckip-1 CDS) remained unchanged after SD. Significantly, adenovirus-mediated overexpression of Ckip-1 3'UTR alleviated SD-induced cardiac dysfunction and remodeling by activating CaMKK2/AMPK/cTNI pathway, which proposed the therapeutic potential of Ckip-1 3'UTR in treating SD-induced heart disease.

Quantitative profiling of m6A at single base resolution across the life cycle of rice and Arabidopsis.

N6-methyladenosine (m6A) plays critical roles in regulating mRNA metabolism. However, comprehensive m6A methylomes in different plant tissues with single-base precision have yet to be reported. Here, we present transcriptome-wide m6A maps at single-base resolution in different tissues of rice and Arabidopsis using m6A-SAC-seq. Our analysis uncovers a total of 205,691 m6A sites distributed across 22,574 genes in rice, and 188,282 m6A sites across 19,984 genes in Arabidopsis. The evolutionarily conserved m6A sites in rice and Arabidopsis ortholog gene pairs are involved in controlling tissue development, photosynthesis and stress response. We observe an overall mRNA stabilization effect by 3' UTR m6A sites in certain plant tissues. Like in mammals, a positive correlation between the m6A level and the length of internal exons is also observed in plant mRNA, except for the last exon. Our data suggest an active m6A deposition process occurring near the stop codon in plant mRNA. In addition, the MTA-installed plant mRNA m6A sites correlate with both translation promotion and translation suppression, depicting a more complicated regulatory picture. Our results therefore provide in-depth resources for relating single-base resolution m6A sites with functions in plants and uncover a suppression-activation model controlling m6A biogenesis across species.

Introduction of miR-3613-3p as a regulator of transforming growth factor-ß (TGF-ß) signaling pathway in colorectal cancer.

INTRODUCTION: Colorectal cancer (CRC) is the second common cancer and the fourth major reason of cancer death worldwide. Dysregulation of intracellular pathways, such as TGF-ß/SMAD signaling, contributes to CRC development. MicroRNAs (miRNAs) are post-transcriptional regulators that are involved in CRC pathogenesis. Here, we aimed to investigate the effect of miR-3613-3p on the TGF-ß /SMAD signaling pathway in CRC. METHODS & RESULTS: Bioinformatics analysis suggested that miR-3613-3p is a regulator of TGF-Β signaling downstream genes. Then, miR-3613-3p overexpression was followed by downregulation of TGF-ßR1, TGF-ßR2, and SMAD2 expression levels, detected by RT-qPCR. Additionally, dual luciferase assay supported the direct interaction of miR-3613-3p with 3'UTR sequences of TGF-ßR1 and TGF-ßR2 genes. Furthermore, reduced SMAD3 protein level following the miR-3613-3p overexpression verified its suppressive effect against TGF-ß signaling in HCT-116 cells, detected by western blot analysis. Finally, miR-3613-3p overexpression induced sub-G1 arrest in HCT116 cells, detected by flow cytometry, and promoted downregulation of cyclin D1 protein expression, which was detected by western blotting analysis. CONCLUSION: Our findings indicated that miR-3613-3p plays an important role in CRC by targeting the TGF-ß/SMAD signaling pathway and could be considered as a new candidate for further therapy investigations.

Deoxynivalenol induces m6A-mediated upregulation of p21 and growth arrest of mouse hippocampal neuron cells in vitro.

Hippocampal neurons maintain the ability of proliferation throughout life to support neurogenesis. Deoxynivalenol (DON) is a mycotoxin that exhibits brain toxicity, yet whether and how DON affects hippocampal neurogenesis remains unknown. Here, we use mouse hippocampal neuron cells (HT-22) as a model to illustrate the effects of DON on neuron proliferation and to explore underlying mechanisms. DON exposure significantly inhibits the proliferation of HT-22 cells, which is associated with an up-regulation of cell cycle inhibitor p21 at both mRNA and protein levels. Global and site-specific m6A methylation levels on the 3'UTR of p21 mRNA are significantly increased in response to DON treatment, whereas inhibition of m6A hypermethylation significantly alleviates DON-induced cell cycle arrest. Further mechanistic studies indicate that the m6A readers YTHDF1 and IGF2BP1 are responsible for m6A-mediated increase in p21 mRNA stability. Meanwhile, 3'UTR of E3 ubiquitin ligase TRIM21 mRNA is also m6A hypermethylated, and another m6A reader YTHDF2 binds to the m6A sites, leading to decreased TRIM21 mRNA stability. Consequently, TRIM21 suppression impairs ubiquitin-mediated p21 protein degradation. Taken together, m6A-mediated upregulation of p21, at both post-transcriptional and post-translational levels, contributes to DON-induced inhibition of hippocampal neuron proliferation. These results may provide new insights for epigenetic therapy of neurodegenerative diseases.

A Systemic Review and Meta-analysis on Natural Resistance-associated Macrophage Protein 1 (3'-Untranslated Region) and Nucleotide-binding Oligomerization Domain-2 (rs8057341) Polymorphisms and Leprosy Susceptibility in Asian and Caucasian Populations.

The current meta-analysis aims to explore the potential correlation between natural resistance-associated macrophage protein 1 (NRAMP1) (3'-Untranslated region [3'-UTR]) and nucleotide-binding oligomerization domain-2 (NOD2 [rs8057341]) gene polymorphisms and their association with leprosy susceptibility in both Asian and Caucasian populations. Datas were retrieved from case control studies with NOD 2 and NRAMP 1 gene polymorphism associated with leprosy disease. Leprosy emerges as a particularly distinctive ailment among women on a global scale. The NRAMP1 (3'-UTR) and NOD2 (rs8057341) genetic variations play a crucial role in the progression of leprosy. A systematic review of relevant case-control studies was conducted across several databases, including ScienceDirect, PubMed, Google Scholar, and Embase. Utilizing MetaGenyo and Review Manager 5.4 Version, statistical analyses were carried out. Nine case-control studies totaling 3281 controls and 3062 leprosy patients are included in the research, with the objective of examining the potential association between NRAMP1 (3'-UTR) and NOD2 (rs8057341) gene polymorphisms and leprosy risk. The review methodology was registered in PROSPERO (ID520883). The findings reveal a robust association between NRAMP1 (3'-UTR) and NOD2 (rs8057341) gene polymorphisms and leprosy risk across various genetic models. Although the funnel plot analysis did not identify publication bias, bolstering these findings and elucidating potential gene-gene and gene-environment interactions require further comprehensive epidemiological research. This study identified a strong correlation between polymorphisms in the NOD2 (rs8057341) genes and susceptibility to leprosy across two genetic models. Further comprehensive epidemiological investigations are warranted to validate these findings and explore potential interactions between these genes and environmental factors.

Additional prognostic value of polymorphisms within the 3'-untranslated region of programmed cell death pathway genes in early-stage breast cancer.

Introduction: The programmed cell death (PCD) pathway plays an important role in restricting cancer cell survival and proliferation. However, limited studies have investigated the association between genetic variants in the 3'-untranslated region of the PCD pathway genes and breast cancer outcomes. Methods: In this study, we genotyped 28 potentially functional single nucleotide polymorphisms (SNPs) in 23 PCD pathway genes in 1,177 patients with early-stage breast cancer (EBC) from a Han Chinese population. The median follow-up period was 174 months. Results: Among all the candidate SNPs, four independent SNPs (rs4900321 and rs7150025 in ATG2B, rs6753785 in BCL2L11, and rs2213181 in c-Kit) were associated with invasive disease-free survival (iDFS), distant disease-free survival (DDFS), breast cancer-specific survival (BCSS) and overall survival (OS), respectively. Further combined genotypes of these four SNPs revealed that the survival decreased as the number of unfavorable genotypes increased (Ptrend = 1.0 × 10-6, 8.5 × 10-8, 3.6 × 10-4, and 1.3 × 10-4 for iDFS, DDFS, BCSS, and OS, respectively). Receiver operating characteristic curve analysis demonstrated that incorporating unfavorable genotypes and clinicopathological variables improved the ability to predict EBC survival (P = 0.006, 0.004, 0.029, and 0.019 for iDFS, DDFS, BCSS, and OS, respectively). Additionally, rs6753785 and rs2213181 were associated with BCL2L11 and c-Kit mRNA expression, respectively. Conclusions: Our results suggest that these four SNPs may act as novel biomarkers for EBC survival, possibly by modulating the expression of the corresponding genes.

HuR facilitates miR-93-5p-induced activation of MAP3K2 translation via MAP3K2 3'UTR ARE2 in hepatocellular carcinoma.

MicroRNAs (miRNAs) can positively regulate gene expression through an unconventional RNA activation mechanism involving direct targeting 3' untranslated regions (UTRs). Our prior study found miR-93-5p activates mitogen-activated protein kinase kinase kinase 2 (MAP3K2) in hepatocellular carcinoma (HCC) via its 3'UTR. However, the underlying mechanism remains elusive. Here, we identified two candidate AU-rich element (ARE) motifs (ARE1 and ARE2) adjacent to the miR-93-5p binding site located within the MAP3K2 3'UTR using AREsite2. Luciferase reporter and translation assays validated that only ARE2 participated in MAP3K2 activation. Integrative analysis revealed that human antigen R (HuR), an ARE2-associated RNA-binding protein (RBP), physically and functionally interacted with the MAP3K2 3'UTR. Consequently, an HuR-ARE2 complex was shown to facilitate miR-93-5p-mediated upregulation of MAP3K2 expression. Furthermore, bioinformatics analysis and studies of HCC cells and specimens highlighted an oncogenic role for HuR and positive HuR-MAP3K2 expression correlation. HuR is also an enhancing factor in the positive feedback circuit comprising miR-93-5p, MAP3K2, and c-Jun demonstrated in our prior study. The newly identified HuR-ARE2 involvement enriches the mechanism of miR-93-5p-driven MAP3K2 activation and suggests new therapeutic strategies warranted for exploration in HCC.

hsa-miR-CHA2, a novel microRNA, exhibits anticancer effects by suppressing cyclin E1 in human non-small cell lung cancer cells.

Despite considerable therapeutic advancements, the global survival rate for lung cancer patients remains poor, posing challenges in developing an effective treatment strategy. In many cases, microRNAs (miRNAs) exhibit abnormal expression levels in cancers, including lung cancer. Dysregulated miRNAs often play a crucial role in the development and progression of cancer. Therefore, understanding the mechanisms underlying aberrant miRNA expression during carcinogenesis may provide crucial clues to develop novel therapeutics. In this study, we identified and cloned a novel miRNA, hsa-miR-CHA2, which is abnormally downregulated in non-small cell lung cancer (NSCLC)-derived cell lines and tissues of patients with NSCLC. Furthermore, we found that hsa-miR-CHA2 regulates the post-transcriptional levels of Cyclin E1 (CCNE1) by binding to the 3'-UTR of CCNE1 mRNA. CCNE1, a cell cycle regulator involved in the G1/S transition, is often amplified in various cancers. Notably, hsa-miR-CHA2 overexpression led to the alteration of the Rb-E2F pathway, a significant signaling pathway in the cell cycle, by targeting CCNE1 in A549 and SK-LU-1 cells. Subsequently, we confirmed that hsa-miR-CHA2 induced G1-phase arrest and exhibited an anti-proliferative effect by targeting CCNE1. Moreover, in subcutaneous xenograft mouse models, intra-tumoral injection of polyplexed hsa-miR-CHA2 mimic suppressed tumor growth and development. In conclusion, hsa-miR-CHA2 exhibited an anticancer effect by targeting CCNE1 both in vitro and in vivo. These findings suggest the potential role of hsa-miR-CHA2 as an important regulator of cell proliferation in molecular-targeted therapy for NSCLC.

Quantifying 3'UTR length from scRNA-seq data reveals changes independent of gene expression.

Although more than half of all genes generate transcripts that differ in 3'UTR length, current analysis pipelines only quantify the amount but not the length of mRNA transcripts. 3'UTR length is determined by 3' end cleavage sites (CS). We map CS in more than 200 primary human and mouse cell types and increase CS annotations relative to the GENCODE database by 40%. Approximately half of all CS are used in few cell types, revealing that most genes only have one or two major 3' ends. We incorporate the CS annotations into a computational pipeline, called scUTRquant, for rapid, accurate, and simultaneous quantification of gene and 3'UTR isoform expression from single-cell RNA sequencing (scRNA-seq) data. When applying scUTRquant to data from 474 cell types and 2134 perturbations, we discover extensive 3'UTR length changes across cell types that are as widespread and coordinately regulated as gene expression changes but affect mostly different genes. Our data indicate that mRNA abundance and mRNA length are two largely independent axes of gene regulation that together determine the amount and spatial organization of protein synthesis.

The 3'UTR Polymorphisms in the NLRP3 Gene Associated with the Risk of COPD and Their Putative Effects on the microRNA Mechanism.

Aims: Evaluating the association between a single nucleotide polymorphism in the 3' untranslated region (3'UTR) of the miRNA binding site of the NLRP3 gene and the occurrence and development of chronic obstructive pulmonary disease (COPD) and providing information to aid in the early detection and treatment of COPD. Materials and Methods: The regulatory single nuclear polymorphisms (SNPs) located in NLRP3 3'UTR were searched by using the dbSNP database and miRNA binding site prediction database. Meanwhile, samples from COPD patients and healthy controls in the same period were used for verification. The clinical baseline information of all subjects was collected, and the transcription level and protein expression level of NLRP3 and the expression level of inflammatory factors downstream of NLRP3 were detected. The effects of SNPs' single nucleotide changes on the transcription and expression of inflammatory factors were analyzed. Results: The study included 418 participants (249 in the COPD group and 169 in the control group). NLRP3 SNPs with miRNA binding sites include rs10754558 (G > C), rs1664774076 (ATAT > del), and rs1664775106 (C > G). Furthermore, two genotypes, GCG and GCA, were discovered to have a linkage mutation at 3'UTR 459-461. COPD susceptibility is tightly associated with the expression of the rs1664774076 del/del genotype, and the risk of COPD increased by 2.770 times (p = 0.003). Type 459-461 GCA was substantially related to the likelihood of developing COPD at various stages (p < 0.05). Except for rs10754558, all homozygous mutants increased NLRP3 mRNA and protein levels. NLRP3 had the greatest area under the receiver operating characteristic (ROC) curve for predicting the development and diagnosis of COPD when compared with its downstream inflammatory variables (AUC = 0.9291). Conclusions: The NLRP3 rs1664774076 del/del genotype is a COPD susceptibility gene, and the GCA genotype at 459-461 can be used as an early predictor of COPD exacerbation. The NLRP3 3'UTR polymorphism may alter the loss of miRNA binding sites, leading to an increase in NLRP3 expression. In the development of COPD, NLRP3 has a better diagnostic value than traditional inflammatory factors. The Clinical Trials Registration number Z: protocol KY01-2020-11-06.

CPEB3 Maintains Developmental Competence of the Oocyte.

Mammalian oocyte development depends on the temporally controlled translation of maternal transcripts, particularly in the coordination of meiotic and early embryonic development when transcription has ceased. The translation of mRNA is regulated by various RNA-binding proteins. We show that the absence of cytoplasmic polyadenylation element-binding protein 3 (CPEB3) negatively affects female reproductive fitness. CPEB3-depleted oocytes undergo meiosis normally but experience early embryonic arrest due to a disrupted transcriptome, leading to aberrant protein expression and the subsequent failure of embryonic transcription initiation. We found that CPEB3 stabilizes a subset of mRNAs with a significantly longer 3'UTR that is enriched in its distal region with cytoplasmic polyadenylation elements. Overall, our results suggest that CPEB3 is an important maternal factor that regulates the stability and translation of a subclass of mRNAs that are essential for the initiation of embryonic transcription and thus for embryonic development.

Determinants of gastric cancer immune escape identified from non-coding immune-landscape quantitative trait loci.

The landscape of non-coding mutations in cancer progression and immune evasion is largely unexplored. Here, we identify transcrptome-wide somatic and germline 3' untranslated region (3'-UTR) variants from 375 gastric cancer patients from The Cancer Genome Atlas. By performing gene expression quantitative trait loci (eQTL) and immune landscape QTL (ilQTL) analysis, we discover 3'-UTR variants with cis effects on expression and immune landscape phenotypes, such as immune cell infiltration and T cell receptor diversity. Using a massively parallel reporter assay, we distinguish between causal and correlative effects of 3'-UTR eQTLs in immune-related genes. Our approach identifies numerous 3'-UTR eQTLs and ilQTLs, providing a unique resource for the identification of immunotherapeutic targets and biomarkers. A prioritized ilQTL variant signature predicts response to immunotherapy better than standard-of-care PD-L1 expression in independent patient cohorts, showcasing the untapped potential of non-coding mutations in cancer.

Long-read sequencing for 29 immune cell subsets reveals disease-linked isoforms.

Alternative splicing events are a major causal mechanism for complex traits, but they have been understudied due to the limitation of short-read sequencing. Here, we generate a full-length isoform annotation of human immune cells from an individual by long-read sequencing for 29 cell subsets. This contains a number of unannotated transcripts and isoforms such as a read-through transcript of TOMM40-APOE in the Alzheimer's disease locus. We profile characteristics of isoforms and show that repetitive elements significantly explain the diversity of unannotated isoforms, providing insight into the human genome evolution. In addition, some of the isoforms are expressed in a cell-type specific manner, whose alternative 3'-UTRs usage contributes to their specificity. Further, we identify disease-associated isoforms by isoform switch analysis and by integration of several quantitative trait loci analyses with genome-wide association study data. Our findings will promote the elucidation of the mechanism of complex diseases via alternative splicing.

Variant-to-function analysis of the childhood obesity chr12q13 locus implicates rs7132908 as a causal variant within the 3' UTR of FAIM2.

The ch12q13 locus is among the most significant childhood obesity loci identified in genome-wide association studies. This locus resides in a non-coding region within FAIM2; thus, the underlying causal variant(s) presumably influence disease susceptibility via cis-regulation. We implicated rs7132908 as a putative causal variant by leveraging our in-house 3D genomic data and public domain datasets. Using a luciferase reporter assay, we observed allele-specific cis-regulatory activity of the immediate region harboring rs7132908. We generated isogenic human embryonic stem cell lines homozygous for either rs7132908 allele to assess changes in gene expression and chromatin accessibility throughout a differentiation to hypothalamic neurons, a key cell type known to regulate feeding behavior. The rs7132908 obesity risk allele influenced expression of FAIM2 and other genes and decreased the proportion of neurons produced by differentiation. We have functionally validated rs7132908 as a causal obesity variant that temporally regulates nearby effector genes and influences neurodevelopment and survival.

HuR promotes triglyceride synthesis and intestinal fat absorption.

Triacylglyceride (TAG) synthesis in the small intestine determines the absorption of dietary fat, but the underlying mechanisms remain to be further studied. Here, we report that the RNA-binding protein HuR (ELAVL1) promotes TAG synthesis in the small intestine. HuR associates with the 3' UTR of Dgat2 mRNA and intron 1 of Mgat2 pre-mRNA. Association of HuR with Dgat2 3' UTR stabilizes Dgat2 mRNA, while association of HuR with intron 1 of Mgat2 pre-mRNA promotes the processing of Mgat2 pre-mRNA. Intestinal epithelium-specific HuR knockout reduces the expression of DGAT2 and MGAT2, thereby reducing the dietary fat absorption through TAG synthesis and mitigating high-fat-diet (HFD)-induced non-alcoholic fatty liver disease (NAFLD) and obesity. Our findings highlight a critical role of HuR in promoting dietary fat absorption.

Nuclear export is a limiting factor in eukaryotic mRNA metabolism.

The eukaryotic mRNA life cycle includes transcription, nuclear mRNA export and degradation. To quantify all these processes simultaneously, we perform thiol-linked alkylation after metabolic labeling of RNA with 4-thiouridine (4sU), followed by sequencing of RNA (SLAM-seq) in the nuclear and cytosolic compartments of human cancer cells. We develop a model that reliably quantifies mRNA-specific synthesis, nuclear export, and nuclear and cytosolic degradation rates on a genome-wide scale. We find that nuclear degradation of polyadenylated mRNA is negligible and nuclear mRNA export is slow, while cytosolic mRNA degradation is comparatively fast. Consequently, an mRNA molecule generally spends most of its life in the nucleus. We also observe large differences in the nuclear export rates of different 3'UTR transcript isoforms. Furthermore, we identify genes whose expression is abruptly induced upon metabolic labeling. These transcripts are exported substantially faster than average mRNAs, suggesting the existence of alternative export pathways. Our results highlight nuclear mRNA export as a limiting factor in mRNA metabolism and gene regulation.

The effect of miR-23b-3p on regulating GH by targeting POU1F1 in Yanbian yellow cattle.

This study aimed to evaluate the effect of miR-23b-3p on growth hormone (GH) in pituitary cells of Yanbian yellow cattle. The mRNA and protein levels of GH and miR-23b-3p target genes were measured by real time fluorescence quantitative PCR (qPCR) and Western blot, respectively. The target relationship of miR-23b-3p was validated by double luciferase reporter gene system. The results showed that GH mRNA and protein levels in pituitary cells of Yanbian yellow cattle were significantly lower in the miR-23b-3p-mi group than in the NC group (P<0.01), while GH mRNA and protein levels were higher in the miR-23b-3p-in group than in the iNC group (P<0.05). The result of bioinformatics analysis and double luciferase reporter gene system validation proved that miR-23b-3p targeted 3'UTR of pituitary specific transcription factor 1 (POU1F1). POU1F1 mRNA and protein levels were lower miR-23b-3p-mi group than in the NC group (P<0.01), while POU1F1 mRNA and protein levels were higher in the miR-23b-3p-in group than in the iNC group (P<0.01). These results demonstrated that miR-23b-3p could regulate GH expression in pituitary cells by regulating POU1F1 gene.

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.

Association of deletion polymorphism rs10573247 in the HMGA2 gene with the risk of breast cancer: bioinformatic and experimental analyses.

BACKGROUND: The high mobility group A2 (HMGA2) gene is expressed extensively during early embryonic development but is inactivated in adulthood, and it is also reactivated in various benign and malignant tumors, including breast cancer. We first assessed the potential functional significance of the unstudied deletion polymorphism rs10573247 at the 3'UTR of HMGA2 on miRNA binding using bioinformatic tools, and subsequently, the association between this polymorphism and breast cancer susceptibility was investigated. MATERIALS AND METHODS: We applied the RNAhybrid tool to predict the functional effects of polymorphism rs10573247 located within the 3' UTR of the HMGA2 gene on miRNA binding. Then, following DNA extraction, 141 breast cancer patients and 123 healthy controls were genotyped for polymorphism rs10573247 using RFLP-PCR with the restriction enzyme Eam1104I. RESULTS: Our bioinformatic data have shown that polymorphism rs10573247 is located in the region that serves as a potential target site for eight miRNAs binding. Among them, miR-3125 exhibited decreased binding affinity for the allele delTT (MFE = -21.8) when compared to the allele TT (MFE = -23.9), but miR-4476 increased binding affinity for the allele delTT (MFE = -22.4) compared to the allele TT (MFE = -22.2). In addition, our results showed that the genotype TT/delTT (p = 0.005) and the genotype delTT/delTT (p = 0.029) were significantly associated with an increased risk of developing breast cancer compared to the genotype TT/TT using RFLP-PCR. DISCUSSION AND CONCLUSION: Our findings suggest that polymorphism rs10573247 may contribute to the risk of breast cancer through the functional effect of this polymorphism on miRNA binding.