Genotype patterns and characteristics of PRNP in the Korean population

Creutzfeldt-Jakob disease (CJD), included in the human transmissible spongiform encephalopathies (TSE), is widely known to be caused by an abnormal accumulation of misfolding prion protein in the brain. Human prion protein gene (PRNP) is mapped in chromosome 20p13 and many single nucleotide polymorphisms (SNPs) in PRNP have been discovered. However, the functionality of SNPs in PRNP is yet unclear, though several SNPs have been known as important mutation related with susceptibility human prion diseases. Our aim is to identify specific genotype patterns and characteristics in the PRNP genomic region and to understand susceptibility among Korean discriminated prion disease patients, suspected CJD patients and the KARE data group. Here, we have researched genotypes and SNPs allele frequencies in PRNP in discriminated prion disease patients group (n = 22), suspected prion diseases patients group (n = 163) and the Korea Association REsource (KARE) data group (n = 296) in Korea. The sequencing regions were promoter region, exon1 and exon2 with their junction parts among 481 samples. A total of 25 SNPs were shown in this study. Nucleotide frequencies of all SNPs are exceedingly tended to bias toward dominant homozygote types except in rs2756271. Genotype frequencies at codon 129 and 219 coding region were similar with previous studies in Korea and Japan. Pathogenic mutations such as 102P/L, 200E/K and 203V/I were observed in discriminated CJD patients group, and 180V/I and 232M/R were shown in suspected prion disease patients group and the KARE data group. A total of 10 SNPs were newly identified, six in the promoter region, one in exon 2 and three in the 3′ UTR. The strong and unique linkage disequilibrium (D' = 0.94, r² = 0.89) was observed between rs57633656 and rs1800014 which is located in codon 219 coding region. We expect that these data can be provided to determine specific susceptibility and a protective factor of prion diseases not only in Koreans but also in East Asians.     

Genotyping of five single nucleotide polymorphisms in the OCA2 and HERC2 genes associated with blue‐brown eye color in the Japanese population

Human eye color is a polymorphic phenotype influenced by multiple genes. It has recently been reported that three single nucleotide polymorphisms (SNPs) within intron 1 of the OCA2 gene (rs7495174, rs4778241, rs4778138) and two SNPs in intron 86 (rs12913832) and the 3′ UTR region (rs1129038) of the HERC2 gene—located in the upstream of the OCA2 locus —have a high statistical association with human eye color. The present study is the first to examine in detail the genotype and haplotype frequencies for these five SNPs in an Asian (Japanese) population (n = 523) comprising solely brown‐eyed individuals. Comparison of the genotype and haplotype distributions in Japanese with those in African and European subjects revealed significant differences between Japanese and other populations. Analysis of haplotypes consisting of four SNPs at the HERC2‐OCA2 locus (rs12913832/rs7495174/rs4778241/rs4778138) showed that the most frequent haplotype in the Japanese population is A‐GAG (0.568), while the frequency of this haplotype is rather low in the European population, even in the brown‐eyed group (0.167). The haplotype distribution in the Japanese population was significantly different from that in the brown‐eyed European group (F_ST = 0.18915). Copyright © 2009 John Wiley & Sons, Ltd.     

A functional variant in ST2 gene is associated with risk of hypertension via interfering MiR‐202‐3p

Recent studies have suggested that interleukin 1 receptor‐like 1 (ST2) plays a critical role in pathogenesis of several cardiovascular disease conditions. In this study, we examined association of 13 single nucleotide polymorphisms (SNPs) of ST2 gene with essential hypertension (EH) risk in 1151 patients with EH and 1135 controls. Our study showed that variants rs11685424, rs12999364 and rs3821204 are highly associated with an increase in risk of EH, while rs6543116 is associated with a decrease risk of EH. Notably, in silico analyses suggested the G>C change of rs3821204, which located within the 3′UTR of soluble ST2 mRNA, disrupted a putative binding site for miR202‐3p. Functional analyses suggested that miR‐202‐3p significantly decreased soluble ST2‐G mRNA stability and inhibited its endogenous expression. Furthermore, we found increased plasma‐soluble ST2 (sST2) level was highly associated with CC genotype of rs3821204 in vivo. Taken together, our findings provide the first evidence that genetic variants in ST2 gene are associated with EH risk and variant rs3821204 may influence the development of EH by controlling sST2 expression.     

Genotype patterns and characteristics of PRNP in the Korean population

Creutzfeldt-Jakob disease (CJD), included in the human transmissible spongiform encephalopathies (TSE), is widely known to be caused by an abnormal accumulation of misfolding prion protein in the brain. Human prion protein gene (PRNP) is mapped in chromosome 20p13 and many single nucleotide polymorphisms (SNPs) in PRNP have been discovered. However, the functionality of SNPs in PRNP is yet unclear, though several SNPs have been known as important mutation related with susceptibility human prion diseases. Our aim is to identify specific genotype patterns and characteristics in the PRNP genomic region and to understand susceptibility among Korean discriminated prion disease patients, suspected CJD patients and the KARE data group. Here, we have researched genotypes and SNPs allele frequencies in PRNP in discriminated prion disease patients group (n = 22), suspected prion diseases patients group (n = 163) and the Korea Association REsource (KARE) data group (n = 296) in Korea. The sequencing regions were promoter region, exon1 and exon2 with their junction parts among 481 samples. A total of 25 SNPs were shown in this study. Nucleotide frequencies of all SNPs are exceedingly tended to bias toward dominant homozygote types except in rs2756271. Genotype frequencies at codon 129 and 219 coding region were similar with previous studies in Korea and Japan. Pathogenic mutations such as 102P/L, 200E/K and 203V/I were observed in discriminated CJD patients group, and 180V/I and 232M/R were shown in suspected prion disease patients group and the KARE data group. A total of 10 SNPs were newly identified, six in the promoter region, one in exon 2 and three in the 3′ UTR. The strong and unique linkage disequilibrium (D' = 0.94, r² = 0.89) was observed between rs57633656 and rs1800014 which is located in codon 219 coding region. We expect that these data can be provided to determine specific susceptibility and a protective factor of prion diseases not only in Koreans but also in East Asians.     

Deep sequencing study of the MTHFR gene to identify variants associated with myelomeningocele

INTRODUCTION: Neural tube defects (NTDs) are congenital anomalies caused by a combination of genetic and environmental influences. A defect below the head region resulting in protuberance of meninges and nervous tissue is termed myelomeningocele (MM). MM, the most common NTD compatible with survival, occurs in approximately 1 in 1000 births worldwide. Maternal preconceptional and periconceptional folate supplementation reduces the risk of NTDs by up to 70%. A key enzyme in folate metabolism is 5, 10‐methylene‐tetrahydrofolate reductase (MTHFR). OBJECTIVES: Sequence the 12 exons of the MTHFR gene among 96 subjects with MM to identify variants potentially contributing to the disease trait. METHODS: Exons were amplified by polymerase chain reaction, and the products were sequenced with the Sanger method to reveal sequence variants compared to MTHFR reference sequences. Association of variants was examined by Fisher's test. RESULTS: A novel variant c.171+3G>T was identified in intron 1 in one affected subject. The variant was not found in the subject's unaffected mother's DNA, and the unaffected father's DNA was unavailable. We found significant differences in allele frequencies for seven SNPs in MM subjects compared with ethnically matched reference populations reported in the single nucleotide polymorphism database. CONCLUSION: We identified a novel variant c.171+3G>T in the MTHFR gene that potentially affects splicing in an affected subject. In addition, we observed five SNPs (rs13306561, rs2274976, rs2066462, rs12121543, and rs1476413) in the MTHFR gene not previously shown to associate with MM. The current study provides additional evidence that multiple variations in the MTHFR gene are associated with MM. Birth Defects Research (Part A) 2012. © 2012 Wiley Periodicals, Inc.     

Evaluation of polymorphisms in microRNA‐binding sites and pancreatic cancer risk in Chinese population

As promising biomarkers and therapy targets, microRNAs (miRNAs) are involved in various physiological and tumorigenic processes. Genetic variants in miRNA‐binding sites can lead to dysfunction of miRNAs and contribute to disease. However, systematic investigation of the miRNA‐related single nucleotide polymorphisms (SNPs) for pancreatic cancer (PC) risk remains elusive. We performed integrative bioinformatics analyses to select 31 SNPs located in miRNA‐target binding sites using the miRNASNP v2.0, a solid database providing miRNA‐related SNPs for genetic research, and investigated their associations with risk of PC in two large case‐control studies totally including 1847 cases and 5713 controls. We observed that the SNP rs3802266 is significantly associated with increased risk of PC (odds ratio (OR) = 1.21, 95% confidence intervals (CI) = 1.11‐1.31, P = 1.29E‐05). Following luciferase reporter gene assays show that rs3802266‐G creates a stronger binding site for miR‐181a‐2‐3p in 3′ untranslated region (3′UTR) of the gene ZHX2. Expression quantitative trait loci (eQTL) analysis suggests that ZHX2 expression is lower in individuals carrying rs3802266‐G with increased PC risk. In conclusion, our findings highlight the involvement of miRNA‐binding SNPs in PC susceptibility and provide new clues for PC carcinogenesis.     

Remarkable genetic diversity detected at river buffalo prolactin receptor (PRLR) gene and association studies with milk fatty acid composition

Prolactin is an anterior pituitary peptide hormone involved in many different endocrine activities and is essential for reproductive performance. This action is mediated by its receptor, the prolactin receptor, encoded by the PRLR gene. In this study, we sequenced and characterized the Mediterranean river buffalo PRLR gene (from exon 3 to 10), and we found remarkable genetic diversity. In particular, we found 24 intronic polymorphisms and 13 exonic SNPs, seven of which were non‐synonymous. Furthermore, the polymorphisms identified in the 3′‐UTR were investigated to establish their possible influence on microRNA binding sites. Considering all the amino acid changes and the observed allelic combinations, it is possible to deduce at least six different translations of the buffalo prolactin receptor and, consequently, the presence at the PRLR gene of at least six alleles. Furthermore, we identified a deletion of a CACTACC heptamer between nucleotides 1102 and 1103 of exon 10 (3′‐UTR), and we developed an allele‐specific PCR to identify the carriers of this genetic marker. Finally, the SNP g.11188A>G, detected in exon 10 and responsible for the amino acid replacement p.His328Arg, was genotyped in 308 Italian Mediterranean river buffaloes, and an association study with milk fat traits was carried out. The statistical analysis showed a tendency that approached significance for the AA genotype with higher contents of odd branched‐chain fatty acids. Thus, our results suggest that the PRLR gene is a good candidate for gene association studies with qualitative traits related to buffalo milk production.     

A program for annotating and predicting the effects of single nucleotide polymorphisms,SnpEff

We describe a new computer program, SnpEff, for rapidly categorizing the effects of variants in genome sequences. Once a genome is sequenced, SnpEff annotates variants based on their genomic locations and predicts coding effects. Annotated genomic locations include intronic, untranslated region, upstream, downstream, splice site, or intergenic regions. Coding effects such as synonymous or non-synonymous amino acid replacement, start codon gains or losses, stop codon gains or losses, or frame shifts can be predicted.

Here the use of SnpEff is illustrated by annotating ~356,660 candidate SNPs in ~117 Mb unique sequences, representing a substitution rate of ~1/305 nucleotides, between the Drosophila melanogaster w¹¹¹⁸; iso-2; iso-3 strain and the reference y¹; cn¹ bw¹ sp¹ strain. We show that ~15,842 SNPs are synonymous and ~4,467 SNPs are non-synonymous (N/S ~0.28). The remaining SNPs are in other categories, such as stop codon gains (38 SNPs), stop codon losses (8 SNPs), and start codon gains (297 SNPs) in the 5′UTR. We found, as expected, that the SNP frequency is proportional to the recombination frequency (i.e., highest in the middle of chromosome arms). We also found that start-gain or stop-lost SNPs in Drosophila melanogaster often result in additions of N-terminal or C-terminal amino acids that are conserved in other Drosophila species. It appears that the 5′ and 3′ UTRs are reservoirs for genetic variations that changes the termini of proteins during evolution of the Drosophila genus.

As genome sequencing is becoming inexpensive and routine, SnpEff enables rapid analyses of whole-genome sequencing data to be performed by an individual laboratory.     

Influence of CTNNB1 rs2953 polymorphism on schizophrenia susceptibility in Chinese Han population through modifying miR‐485 binding to CTNNB1

Schizophrenia (SCZ) and bipolar disorder (BD) are two major neuropsychiatric diseases that are the most substantial causes of disability and mortality worldwide. CTNNB1 encodes beta‐catenin, an important protein in canonical Wnt signaling. We aimed to investigate the association between the rs2953 of CTNNB1 and the risk of SCZ and BD and to further explore the function of rs2953. A total of 1658 samples (548 SCZ cases, 512 BD cases, and 598 controls) were examined in terms of the genotype of CTNNB1 rs2953. The mRNA expression level of CTNNB1 significantly increased in the SCZ and BD groups compared with that in the control group. Significant association remained between CTNNB1 3′‐untranslated region (UTR) variant rs2953 and SCZ susceptibility (additive and dominant model) after gender and age were adjusted. rs2953 disrupted the binding of CTNNB1 and miR‐485. miR‐485 significantly suppressed the luciferase activity of CTNNB1‐T vector by binding to the CTNNB1 3′‐UTR containing the T allele of rs2953. The mRNA expression of CTNNB1 can be used as a biomarker for the diagnosis of SCZ and BD. The 3′‐UTR variant rs2953 in CTNNB1 influences the risk of SCZ in the Han Chinese population and modifies the binding of miR‐485 to CTNNB1.     

Rs884225 polymorphism is associated with primary hypertension by compromising interaction between epithelial growth factor receptor (EGFR) and miR-214

Genetic variations in the 3′UTR of mRNAs as well as sequences of microRNAs (miRNAs) and long non-coding RNAs (lncRNAs) can affect gene expression by interfering with the binding between them. In this study, we investigated the role of the following polymorphisms in the risk of hypertension: the 774T > C (rs17337023) polymorphism located in the EGFR 3’ untranslated region (3’UTR), the rs884225 polymorphism located in the sequence of miR-214, and the single nucleotide polymorphisms (SNPs) rs325797437, rs344501106, rs81286029 and rs318656749 located in the promoter of lncRNA MEG3. Taqman genotyping assays and haplotype analysis tools were used to measure the MEG3 haplotypes and the rs17337023 and rs884225 polymorphisms genotypes. The relationship between MEG3, miR-214 and EGFR was validated using computational analysis and luciferase assays. Unlike other polymorphisms, only patients grouped according to their rs884225 genotypes exhibited varied EGFR mRNA and protein levels, which indicated that the rs884225 genotype is associated with the expression of EGFR mRNA and protein levels. MiR-214 was confirmed to bind to MEG3 and 3’UTR of EGFR by showing that the transfection of exogenous miR-214 significantly down-regulated the luciferase activity of A549 and H460 cells transfected with wild-type MEG3 or wild-type EGFR 3’ UTR. Additionally, MEG3 overexpression inhibited miR-214 expression while elevating the EGFR mRNA and protein expressions. Meanwhile, MEG3 down-regulation demonstrated an opposite result, thus establishing the MEG3/miR-214/EGRF signalling pathway. Our study confirmed that the T > C substitution of rs884225 polymorphism located in miR-214 binding site in the 3’UTR of EGFR is associated with increased risk of primary hypertension.     

DNA repair gene polymorphisms, pre-natal factors and the frequency of somatic mutations in the glycophorin-A gene among healthy newborns

This study investigates variation in somatic mutation frequency, as measured by the glycophorin-A (GPA) somatic mutation assay, in relation to polymorphic variation among 435 newborn babies in DNA repair genes XRCC1, XRCC3 and XRCC4 and gender, parental age, social class and smoking habits. The three polymorphisms under investigation were an Arg --> Gln substitution at codon 399 in exon 10 of XRCC1, a Thr --> Met substitution at codon 241 in exon 7 of XRCC3 and an Ile --> Thr substitution at codon 401 in exon 4 of XRCC4. The study population is an extension of that previously analysed for GPA mutations and XRCC1 polymorphisms. A significant difference was seen in the earlier work in the genotype distribution for the XRCC1 Arg399Gln polymorphism between the main population and the small number with extreme values for NN variant frequency and this was maintained in the larger study group (OR 3.20 [95% CI: 1.16, 8.81]) P = 0.043). No such association was seen for XRCC3 or XRCC4 polymorphisms. When adjustments were made for multiple testing, neither N0 nor NN variant frequencies in the main study population were found to be influenced by the polymorphisms in XRCC1, XRCC3, or XRCC4. In addition, neither maternal or paternal smoking, age or social class nor the gender of the offspring were found to affect variant frequencies nor were variant frequencies influenced by any interaction between any of these factors and genotype. It is concluded that the genotypic variation in DNA repair genes examined in this study has no discernable effect on the genesis of the somatic mutations observed at birth.     

Analysis of the SIM1 Contribution to Polygenic Obesity in the French Population

SIM1 (single‐minded 1) haploinsufficiency is responsible for obesity in both humans and mice, but the contribution of frequent DNA variation to polygenic obesity is unknown. Sequencing of all exons, exon/intron boundaries, 870 base pairs (bp) of the putative promoter, and 1,095 bp of the 3′UTR of SIM1 gene in 143 obese children and 24 control adults identified 13 common variants. After analysis of the linkage disequilibrium (LD) structure, association study of eight variants was performed in 1,275 obese children and severely obese adults, in 1,395 control subjects, and in 578 obesity‐selected pedigrees. A nominal evidence of association was found for the nonsynonymous variant P352T C/A (rs3734354) (P = 0.01, OR = 0.81 (0.70–0.95)), the +2,004 TGA ?/insT SNP (rs35180395) (P = 0.02, OR = 1.21 (1.02–1.43)), the +2,215A/G TGA SNP (rs9386126) (P = 0.002, OR = 0.81 (0.71–0.93)), and pooled childhood/adult obesity. Even though transmission disequilibrium test (TDT) further supported the association of P352T and +2,004 ?/inst T with obesity, none of these nominal associations remained significant after a multiple testing Bonferroni correction. Therefore, our study excludes a major contribution of SIM1 common variants in exons, 5′ and 3′ UTR regions in polygenic obesity susceptibility in French Europeans.     

Rare BRCA1 haplotypes including 3’UTR SNPs associated with breast cancer risk

Genetic markers identifying women at an increased risk of developing breast cancer exist, yet the majority of inherited risk remains elusive. While numerous BRCA1 coding sequence mutations are associated with breast cancer risk, BRCA1 mutations account for less then 5% of breast cancer risk. Since 3' untranslated region (3'UTR) polymorphisms disrupting microRNA (miRNA) binding can be functional and can act as genetic markers of cancer risk, we tested the hypothesis that such polymorphisms in the 3'UTR of BRCA1 and haplotypes containing these functional polymorphisms may be associated with breast cancer risk. We sequenced the BRCA1 3'UTR from breast cancer patients to identify miRNA disrupting polymorphisms. We further evaluated haplotypes of this region including the identified 3'UTR variants in a large population of controls and breast cancer patients (n=221) with known breast cancer subtypes and ethnicities. We identified three 3'UTR variants in BRCA1 that are polymorphic in breast cancer populations, and haplotype analysis including these variants revealed that breast cancer patients harbor five rare haplotypes not generally found among controls (9.50% for breast cancer chromosomes, 0.11% for control chromosomes, p=0.0001). Three of these rare haplotypes contain the rs8176318 BRCA1 3'UTR functional variant. These haplotypes are not biomarkers for BRCA1 coding mutations, as they are found rarely in BRCA1 mutant breast cancer patients (1/129 patients= 0.78%). These rare BRCA1 haplotypes and 3'UTR SNPs may represent new genetic markers of breast cancer risk.