Genetic Variations in Key MicroRNA Processing Genes and Risk of Head and Neck Cancer: A Case-Control Study in Chinese Population

MicroRNAs (miRNAs) have been reported to play a key role in oncogenesis. Genetic variations in miRNA processing genes and miRNA binding sites may affect the biogenesis of miRNA and the miRNA-mRNA interactions, hence promoting tumorigenesis. In the present study, we hypothesized that potentially functional polymorphisms in miRNA processing genes may contribute to head and neck cancer (HNC) susceptibility. To test this hypothesis, we genotyped three SNPs at miRNA binding sites of miRNA processing genes (rs1057035 in 3′UTR of DICER, rs3803012 in 3′UTR of RAN and rs10773771 in 3′UTR of HIWI) with a case-control study including 397 HNC cases and 900 controls matched by age and sex in Chinese. Although none of three SNPs was significantly associated with overall risk of HNC, rs1057035 in 3′UTR of DICER was associated with a significantly decreased risk of oral cancer (TC/CC vs. TT: adjusted OR  = 0.65, 95% CI  = 0.46–0.92). Furthermore, luciferase activity assay showed that rs1057035 variant C allele led to significantly lower expression levels as compared to the T allele, which may be due to the relatively high inhibition of hsa-miR-574-3p on DICER mRNA. These findings indicated that rs1057035 located at 3′UTR of DICER may contribute to the risk of oral cancer by affecting the binding of miRNAs to DICER. Large-scale and well-designed studies are warranted to validate our findings.     

Molecular architecture of a miRNA-regulated 3' UTR

Animal genomes contain hundreds of microRNAs (miRNAs), small regulatory RNAs that control gene expression by binding to complementary sites in target mRNAs. Some rules that govern miRNA/target interaction have been elucidated but their general applicability awaits further experimentation on a case-by-case basis. We use here an assay system in transgenic nematodes to analyze the interaction of the Caenorhabditis elegans lsy-6 miRNA with 3' UTR sequences. In contrast to many previously described assay systems used to analyze miRNA/target interactions, our assay system operates within the cellular context in which lsy-6 normally functions, a single neuron in the nervous system of C. elegans. Through extensive mutational analysis, we define features in the known and experimentally validated target of lsy-6, the 3' UTR of the cog-1 homeobox gene, that are required for a functional miRNA/target interaction. We describe that both in the context of the cog-1 3' UTR and in the context of heterologous 3' UTRs, one or more seed matches are not a reliable predictor for a functional miRNA/target interaction. We rather find that two nonsequence specific contextual features beyond miRNA target sites are critical determinants of miRNA-mediated 3' UTR regulation. The contextual features reside 3' of lsy-6 binding sites in the 3' UTR and act in a combinatorial manner; mutation of each results in limited defects in 3' UTR regulation, but a combinatorial deletion results in complete loss of 3' UTR regulation. Together with two lsy-6 sites, these two contextual features are capable of imparting regulation on a heterologous 3' UTR. Moreover, the contextual features need to be present in a specific configuration relative to miRNA binding sites and could either represent protein binding sites or provide an appropriate structural context. We conclude that a given target site resides in a 3' UTR context that evolved beyond target site complementarity to support regulation by a specific miRNA. The large number of 3' UTRs that we analyzed in this study will also be useful to computational biologists in designing the next generation of miRNA/target prediction algorithms.     

ACACβ gene (rs2268388) and AGTR1 gene (rs5186) polymorphism and the risk of nephropathy in Asian Indian patients with type 2 diabetes

Patients with type 2 diabetes (T2DM) are usually obese and concurrent obesity results into activation of the renin–angiotensin-system (RAS) which is a risk factor for diabetic nephropathy (DN). Gene–gene interaction between acetyl-coenzymeA carboxylase beta (ACACβ) gene, which is involved in fatty acid metabolism and angiotensin II receptors (AGTR1) gene, which mediates RAS proteins actions on renal tissue, polymorphism with DN have not been studied earlier. The present study was designed with the aim to examine the association of an ACACβ (rs2268388) and AGTR1 (rs5186) gene polymorphism with the risk of DN in Asian Indians. 1,158 patients with T2DM belonging to two independently ascertained North Indian and one South Indian cohorts were genotyped for ACACβ (rs2268388) and AGTR1 (rs5186) polymorphism using real time PCR-based Taq-man assay and PCR–RFLP assays. In all the three cohorts, a significantly higher frequency of T allele and TT genotypes of ACACβ and C allele and CC genotypes of AGTR1 were found in patients with DN as compared to patients without nephropathy. Further, T allele of ACACβ and C allele of AGTR1 were found to be significantly associated with proteinuria, a hallmark of DN. We also found significant epistatic interactions between these two genes. TT genotypes of ACACβ gene and CC genotype of AGTR1 gene confers the risk of DN and both genes had significant epistatic interaction in Asian Indian patients with T2DM.     

ACE and AGTR1 genes polymorphisms in left ventricular hypertrophy pathogenesis in humans

The role of A2350G polymorphism in exon 17 of the ACE gene and A1166C - in 3'-UTR of the AGTR1 in the pathogenesis of left ventricular hypertrophy was studied in patients with essential hypertension (EH) and arterial hypertension combined with diabetes mellitus type 2 (AH + DM2). Patients with EH and AH + DM2 did not differ from the control sample of healthy individuals by allele or genotype frequencies. However, an association of both polymorphisms with LVH was detected in EH patients. The frequency of 1166C allele was higher in patients with LVH (33.6% vs 20.7% without LVH). A1166C polymorphism determined the magnitude of left ventricular mass index (LVMI) in EH patients as well (p = 0.007). 2350G allele frequency of the ACE gene was in 1.5, and GG genotype--in 3.5-fold higher in EH patients with LVH, as compared without LVH. LVMI was significantly higher in patients with GG genotype as compared with heterozygotes and AA homozygotes (p = 0.002). Thus the presence of 1166C allele of AGTR1 and 2350G allele of ACE can be considered as predisposing factors for LVH development in EH. In contrast, association of studied polymorphisms with presence or LVH degree was not detected in patients with arterial hypertension combined with DM2. This may indicate another structure of genetic component of predisposition to LVH in different causes.     

血管内皮生长因子基因3＇UTR不同基因型载体的构建与鉴定

目的：构建含SNP位点的血管内皮生长因子（VEGF）基因3＇UTR的荧光素酶报告基因载体,为进一步揭示VEGF基因3＇UTR的单核苷酸多态性（SNP）影响肺癌发病风险的分子机制奠定基础。方法：以rs3025039和rs3025040两个位点均为C纯合子的非癌症病人血液DNA为模板,扩增出两位点为C/C单体型、长度为1448 bp的VEGF基因3＇UTR目的片段,测序验证后将其克隆至pMIR-REPORT荧光素酶报告基因载体上,得到重组质粒pMIR-C/C。同时,我们以pMIR-C/C为模板定点突变两个SNP位点,得到具有T/T单体型的重组质粒pMIR-T/T。将各重组质粒转化大肠杆菌DH10B,筛选阳性克隆后提取质粒进行双酶切鉴定及DNA测序鉴定。结果：单菌落质粒测序验证显示带有C/C单体型的VEGF基因3＇UTR重组质粒pMIR-C/C构建成功;经两次定点突变,成功将pMIR-C/C质粒转变为pMIR-T/T,经测序验证未引入任何其他突变。同时生物信息学预测还显示rs3025040位点位于miR-199a/b与VEGF基因mRNA的结合位置,其改变可以影响miRNA与mRNA的结合效率。结论：本研究成功构建了含有两个连锁SNP的VEGF基因3＇UTR的荧光素酶报告基因载体,为今后VEGF基因3＇UTR的功能研究奠定基础。     

血管内皮生长因子基因3'UTR 不同基因型载体的构建与鉴定

目的:构建含SNP位点的血管内皮生长因子(VEGF)基因3’UTR的荧光素酶报告基因载体,为进一步揭示VEGF基因3’UTR的单核苷酸多态性(SNP)影响肺癌发病风险的分子机制奠定基础。方法:以rs3025039和rs3025040两个位点均为C纯合子的非癌症病人血液DNA为模板,扩增出两位点为C/C单体型、长度为1448 bp的VEGF基因3’UTR目的片段,测序验证后将其克隆至pMIR-REPORT荧光素酶报告基因载体上,得到重组质粒pMIR-C/C。同时,我们以pMIR-C/C为模板定点突变两个SNP位点,得到具有T/T单体型的重组质粒pMIR-T/T。将各重组质粒转化大肠杆菌DH10B,筛选阳性克隆后提取质粒进行双酶切鉴定及DNA测序鉴定。结果:单菌落质粒测序验证显示带有C/C单体型的VEGF基因3’UTR重组质粒pMIR-C/C构建成功;经两次定点突变,成功将pMIR-C/C质粒转变为pMIR-T/T,经测序验证未引入任何其他突变。同时生物信息学预测还显示rs3025040位点位于miR-199a/b与VEGF基因mRNA的结合位置,其改变可以影响miRNA与mRNA的结合效率。结论:本研究成功构建了含有两个连锁SNP的VEGF基因3’UTR的荧光素酶报告基因载体,为今后VEGF基因3’UTR的功能研究奠定基础。     

Association of AGTR1 (A1166C) and ACE (I/D) Polymorphisms with Breast Cancer Risk in North Indian Population

Renin angiotensin system (RAS) comprising Angiotensin converting enzyme (ACE), Angiotensin II (Ang II) and its receptor Angiotensin II receptor type I (AGTR1), plays a critical role in several diseases including cancer. A single nucleotide polymorphism (SNP) A1166C located in 3′ untranslated region (UTR) of AGTR1 and an insertion/deletion (I/D) polymorphism present in intron 16 of ACE gene have been associated with many diseases, but their association with Breast cancer (BCa) is still debatable. Here, we for the first time investigated the association of these polymorphisms in a North Indian BCa cohort including 161 patients and 152 healthy women. The polymorphisms were evaluated by polymerase chain reaction (PCR) and restriction fragment length polymorphism (RFLP) respectively. The association between these polymorphisms and BCa risk was estimated by calculating Odds Ratio (OR) and chi-square (χ²) test. The DD genotype/D allele of ACE (I/D) polymorphism and “AC and CC” genotype/C allele of AGTR1 (A1166C) polymorphism were associated with higher risk of BCa when evaluated independently. Furthermore, interaction analysis of “AC and CC” and DD genotype and combination of “C and D” alleles of both polymorphisms revealed significantly greater BCa risk than that observed independently. Conclusively, women harboring “AC or CC” genotype/C allele for AGTR1 (A1166C) polymorphism and DD genotype/D allele for ACE (I/D) polymorphisms have a predisposition to develop more aggressive disease with advanced staging and larger tumor size. Our study indicates importance of genetic screening based on these polymorphisms for women, who may have higher risk of BCa.     

Host microRNAs exhibit differential propensity to interact with SARS-CoV-2 and variants of concern

A significant number of SARS-CoV-2-infected individuals naturally overcome viral infection, suggesting the existence of a potent endogenous antiviral mechanism. As an innate defense mechanism, microRNA (miRNA) pathways in mammals have evolved to restrict viruses, besides regulating endogenous mRNAs. In this study, we systematically examined the complete repertoire of human miRNAs for potential binding sites on SARS-CoV-2 Wuhan-Hu-1, Beta, Delta, and Omicron. Human miRNA and viral genome interaction were analyzed using RNAhybrid 2.2 with stringent parameters to identify highly bonafide miRNA targets. Using publicly available data, we filtered for miRNAs expressed in lung epithelial cells/tissue and oral keratinocytes, concentrating on the miRNAs that target SARS-CoV-2 S protein mRNAs. Our results show a significant loss of human miRNA and SARS-CoV-2 interactions in Omicron (130 miRNAs) compared to Wuhan-Hu-1 (271 miRNAs), Beta (279 miRNAs), and Delta (275 miRNAs). In particular, hsa-miR-3150b-3p and hsa-miR-4784 show binding affinity for S protein of Wuhan strain but not Beta, Delta, and Omicron. Loss of miRNA binding sites on N protein was also observed for Omicron. Through Ingenuity Pathway Analysis (IPA), we examined the experimentally validated and highly predicted functional role of these miRNAs. We found that hsa-miR-3150b-3p and hsa-miR-4784 have several experimentally validated or highly predicted target genes in the Toll-like receptor, IL-17, Th1, Th2, interferon, and coronavirus pathogenesis pathways. Focusing on the coronavirus pathogenesis pathway, we found that hsa-miR-3150b-3p and hsa-miR-4784 are highly predicted to target MAPK13. Exploring miRNAs to manipulate viral genome/gene expression can provide a promising strategy with successful outcomes by targeting specific VOCs.     

ACE and AGTR1 Polymorphisms in Pathogenesis of Human Left Ventricular Hypertrophy

The A2350G polymorphism of exon 17 of the angiotensin I-converting enzyme gene (ACE) and the A1166C polymorphism of the 3-untranslated region (3-UTR) of the angiotensin II type 1 receptor gene (AGTR1) were tested for association with left ventricular hypertrophy (LVH) in patients with essential hypertension (EH) or arterial hypertension (AH) combined with diabetes mellitus type 2 (DM2). The patients with EH or AH + DM2 did not differ significantly in ACE or AGTR1 allele or genotype frequencies from healthy subjects. Both polymorphisms were associated with LVH in EH. AGTR1 allele 1166C was more frequent in patients with LVH than without (33.6 vs. 20.7%) and affected the left ventricular mass index (LVMI) in patients with EH (p = 0.007). The frequency of ACE allele 2350G in EH patients with LVH was 1.5 times higher, and that of genotype GG was 3.5 times higher, than in patients without LVH. LVMI differed significantly (p = 0.002) between patients with different ACE genotypes, being maximum in homozygotes GG and minimum in homozygotes AA. Thus, AGTR1 allele 1166C and ACE allele 2350G were identified as predisposing to LVH in EH. The two polymorphisms were not associated with the incidence or severity of LVH in patients with AH and DM2.     

Human-chimpanzee DNA sequence variation in the four major genes of the renin angiotensin system

The renin angiotensin system (RAS) is involved in blood pressure control and water/sodium metabolism. The genes encoding the proteins of this system are candidate genes for essential hypertension. The RAS involves four main molecules: angiotensinogen, renin, angiotensin I-converting enzyme, and the angiotensin II type 1 receptor (encoded by the genes AGT, REN, DCP1, and AGTR1, respectively). We performed a molecular screening over 17,037 bp of the coding and 5' and 3' untranslated regions of these genes, from three to six common chimpanzees. We identified 44 single-nucleotide polymorphisms (SNPs) in chimpanzee samples, including 18 coding-region SNPs, 5 of which led to an amino acid replacement. We observed common and different features at various sites (synonymous, nonsynonymous, and noncoding) within and between the four chimpanzee genes: (1) the nucleotide diversity at noncoding sites was similar; (2) the nucleotide diversity at nonsynonymous sites was low, probably reflecting purifying selection, except for the AGT gene; (3) the nucleotide diversity at synonymous sites, which was dependent on the G+C content at the third position of the codon, was high, except for the AGTR1 gene. Comparison of the chimpanzee SNPs with those previously reported for humans identified 119 sites with fixed differences (including 62 coding sites, 17 of which resulted in amino acid differences between the species). Analysis of polymorphism within species and divergence between species shed light on the evolutionary constraints on these genes. In particular, comparison of the pattern of mutation at polymorphic and fixed sites between humans and chimpanzees suggested that the high G+C content of the DCP1 gene was maintained by positive selection at its silent sites. Finally, we propose 68 ancestral alleles for the human RAS genes and discuss the implications for their use in future hypertension-susceptibility association studies.     

Gene Silencing Efficiency and INF-β Induction Effects of Splicing miRNA 155-Based Artificial miRNA with Pre-miRNA Stem-Loop Structures

Artificial microRNA (miRNA) expression vectors have been developed and used for RNA interference. The secondary structure of artificial miRNA is important for RNA interference efficacy. We designed two groups of six artificial splicing miRNA 155-based miRNAs (SM155-based miRNAs) with the same target in the coding region or 3' UTR of a target gene and studied their RNA silencing efficiency and interferon β (IFN-β) induction effects. SM155-based miRNA with a mismatch at the +1 position and a bulge at the +11, +12 positions in a miRNA precursor stem-loop structure showed the highest gene silencing efficiency and lowest IFN-β induction effect (increased IFN-β mRNA level by 10% in both target cases), regardless of the specificity of the target sequence, suggesting that pSM155-based miRNA with this design could be a valuable miRNA expression vector.     

Single Nucleotide Polymorphisms in the 3′UTR of VPAC-1 Cooperate in Modulating Gene Expression and Impact Differently on the Interaction with miR525-5p

Complex immune and neurodegenerative disorders are the result of multiple interactions between common genetic variations having, individually, a weak effect on the disease susceptibility or resistance. Interestingly, some genes have been found to be associated with more than one disease although not necessarily the same SNPs are involved. In this context, single nucleotide polymorphisms in the 3′UTR region of type 1 receptor (VPAC-1) for vasoactive intestinal peptide (VIP) have been reported to be associated with some immune-mediated as well as with neurodegenerative diseases such as Alzheimer''s Disease (AD). Here, we demonstrate that variations at the 3′UTR of the VPAC-1 gene act synergistically to affect the expression of the luciferase as well as of the GFP reporter genes expressed in HEK293T cells. Moreover, the miRNA 525-5p, previously shown by us to target the 3′UTR of VPAC-1, is more efficient in decreasing GFP expression when co-expressed with constructs carrying the allele C at rs896 (p<10^-3) suggesting that this miRNA regulates VPAC-1 expression at different levels depending on rs896 polymorphism and thus adding complexity to the network of disease susceptibility.     

SNPs occur in regions with less genomic sequence conservation

Rates of SNPs (single nucleotide polymorphisms) and cross-species genomic sequence conservation reflect intra- and inter-species variation, respectively. Here, I report SNP rates and genomic sequence conservation adjacent to mRNA processing regions and show that, as expected, more SNPs occur in less conserved regions and that functional regions have fewer SNPs. Results are confirmed using both mouse and human data. Regions include protein start codons, 3' splice sites, 5' splice sites, protein stop codons, predicted miRNA binding sites, and polyadenylation sites. Throughout, SNP rates are lower and conservation is higher at regulatory sites. Within coding regions, SNP rates are highest and conservation is lowest at codon position three and the fewest SNPs are found at codon position two, reflecting codon degeneracy for amino acid encoding. Exon splice sites show high conservation and very low SNP rates, reflecting both splicing signals and protein coding. Relaxed constraint on the codon third position is dramatically seen when separating exonic SNP rates based on intron phase. At polyadenylation sites, a peak of conservation and low SNP rate occurs from 30 to 17 nt preceding the site. This region is highly enriched for the sequence AAUAAA, reflecting the location of the conserved polyA signal. miRNA 3' UTR target sites are predicted incorporating interspecies genomic sequence conservation; SNP rates are low in these sites, again showing fewer SNPs in conserved regions. Together, these results confirm that SNPs, reflecting recent genetic variation, occur more frequently in regions with less evolutionarily conservation.