Association between genetic variants in DNA and histone methylation and telomere length

Telomere length, a biomarker of aging and age-related diseases, exhibits wide variation between individuals. Common genetic variation may explain some of the individual differences in telomere length. To date, however, only a few genetic variants have been identified in the previous genome-wide association studies. As emerging data suggest epigenetic regulation of telomere length, we investigated 72 single nucleotide polymorphisms (SNPs) in 46 genes that involve DNA and histone methylation as well as telomerase and telomere-binding proteins and DNA damage response. Genotyping and quantification of telomere length were performed in blood samples from 989 non-Hispanic white participants of the Sister Study, a prospective cohort of women aged 35–74 years. The association of each SNP with logarithmically-transformed relative telomere length was estimated using multivariate linear regression. Six SNPs were associated with relative telomere length in blood cells with p-values<0.05 (uncorrected for multiple comparisons). The minor alleles of BHMT rs3733890 G>A (p = 0.041), MTRR rs2966952 C>T (p = 0.002) and EHMT2 rs558702 G>A (p = 0.008) were associated with shorter telomeres, while minor alleles of ATM rs1801516 G>A (p = 0.031), MTR rs1805087 A>G (p = 0.038) and PRMT8 rs12299470 G>A (p = 0.019) were associated with longer telomeres. Five of these SNPs are located in genes coding for proteins involved in DNA and histone methylation. Our results are consistent with recent findings that chromatin structure is epigenetically regulated and may influence the genomic integrity of telomeric region and telomere length maintenance. Larger studies with greater coverage of the genes implicated in DNA methylation and histone modifications are warranted to replicate these findings.     

Pipeline for large-scale microdroplet bisulfite PCR-based sequencing allows the tracking of hepitype evolution in tumors

Cytosine methylation provides an epigenetic level of cellular plasticity that is important for development, differentiation and cancerogenesis. We adopted microdroplet PCR to bisulfite treated target DNA in combination with second generation sequencing to simultaneously assess DNA sequence and methylation. We show measurement of methylation status in a wide range of target sequences (total 34 kb) with an average coverage of 95% (median 100%) and good correlation to the opposite strand (rho = 0.96) and to pyrosequencing (rho = 0.87). Data from lymphoma and colorectal cancer samples for SNRPN (imprinted gene), FGF6 (demethylated in the cancer samples) and HS3ST2 (methylated in the cancer samples) serve as a proof of principle showing the integration of SNP data and phased DNA-methylation information into “hepitypes” and thus the analysis of DNA methylation phylogeny in the somatic evolution of cancer.     

Von Hippel-Lindau (VHL) inactivation in sporadic clear cell renal cancer: associations with germline VHL polymorphisms and etiologic risk factors

Renal tumor heterogeneity studies have utilized the von Hippel-Lindau VHL gene to classify disease into molecularly defined subtypes to examine associations with etiologic risk factors and prognosis. The aim of this study was to provide a comprehensive analysis of VHL inactivation in clear cell renal tumors (ccRCC) and to evaluate relationships between VHL inactivation subgroups with renal cancer risk factors and VHL germline single nucleotide polymorphisms (SNPs). VHL genetic and epigenetic inactivation was examined among 507 sporadic RCC/470 ccRCC cases using endonuclease scanning and using bisulfite treatment and Sanger sequencing across 11 CpG sites within the VHL promoter. Case-only multivariate analyses were conducted to identify associations between alteration subtypes and risk factors. VHL inactivation, either through sequence alterations or promoter methylation in tumor DNA, was observed among 86.6% of ccRCC cases. Germline VHL SNPs and a haplotype were associated with promoter hypermethylation in tumor tissue (OR = 6.10; 95% CI: 2.28–16.35, p = 3.76E-4, p-global = 8E-5). Risk of having genetic VHL inactivation was inversely associated with smoking due to a higher proportion of wild-type ccRCC tumors [former: OR = 0.70 (0.20–1.31) and current: OR = 0.56 (0.32–0.99); P-trend = 0.04]. Alteration prevalence did not differ by histopathologic characteristics or occupational exposure to trichloroethylene. ccRCC cases with particular VHL germline polymorphisms were more likely to have VHL inactivation through promoter hypermethylation than through sequence alterations in tumor DNA, suggesting that the presence of these SNPs may represent an example of facilitated epigenetic variation (an inherited propensity towards epigenetic variation) in renal tissue. A proportion of tumors from current smokers lacked VHL alterations and may represent a biologically distinct clinical entity from inactivated cases.     

Cardiac myosin binding protein C and MAP-kinase activating death domain-containing gene polymorphisms and diastolic heart failure

Objective

Myosin binding protein C (MYBPC3) plays a role in ventricular relaxation. The aim of the study was to investigate the association between cardiac myosin binding protein C (MYBPC3) gene polymorphisms and diastolic heart failure (DHF) in a human case-control study.

Methods

A total of 352 participants of 1752 consecutive patients from the National Taiwan University Hospital and its affiliated hospital were enrolled. 176 patients diagnosed with DHF confirmed by echocardiography were recruited. Controls were matched 1-to-1 by age, sex, hypertension, diabetes, renal function and medication use. We genotyped 12 single nucleotide polymorphisms (SNPs) according to HapMap Han Chinese Beijing databank across a 40 kb genetic region containing the MYBPC3 gene and the neighboring DNA sequences to capture 100% of haplotype variance in all SNPs with minor allele frequencies ≧5%. We also analyzed associations of these tagging SNPs and haplotypes with DHF and linkage disequilibrium (LD) structure of the MYBPC3 gene.

Results

In a single locus analysis, SNP rs2290149 was associated with DHF (allele-specific p = 0.004; permuted p = 0.031). The SNP with a minor allele frequency of 9.4%, had an odds ratio 2.14 (95% CI 1.25–3.66; p = 0.004) for the additive model and 2.06 for the autosomal dominant model (GG+GA : AA, 95% CI 1.17–3.63; p = 0.013), corresponding to a population attributable risk fraction of 12.02%. The haplotypes in a LD block of rs2290149 (C-C-G-C) was also significantly associated with DHF (odds ratio 2.10 (1.53–2.89); permuted p = 0.029).

Conclusions

We identified a SNP (rs2290149) among the tagging SNP set that was significantly associated with early DHF in a Chinese population.     

Genetic association study identifies HSPB7 as a risk gene for idiopathic dilated cardiomyopathy

Dilated cardiomyopathy (DCM) is a structural heart disease with strong genetic background. Monogenic forms of DCM are observed in families with mutations located mostly in genes encoding structural and sarcomeric proteins. However, strong evidence suggests that genetic factors also affect the susceptibility to idiopathic DCM. To identify risk alleles for non-familial forms of DCM, we carried out a case-control association study, genotyping 664 DCM cases and 1,874 population-based healthy controls from Germany using a 50K human cardiovascular disease bead chip covering more than 2,000 genes pre-selected for cardiovascular relevance. After quality control, 30,920 single nucleotide polymorphisms (SNP) were tested for association with the disease by logistic regression adjusted for gender, and results were genomic-control corrected. The analysis revealed a significant association between a SNP in HSPB7 gene (rs1739843, minor allele frequency 39%) and idiopathic DCM (p = 1.06×10⁻⁶, OR = 0.67 [95% CI 0.57–0.79] for the minor allele T). Three more SNPs showed p < 2.21×10⁻⁵. De novo genotyping of these four SNPs was done in three independent case-control studies of idiopathic DCM. Association between SNP rs1739843 and DCM was significant in all replication samples: Germany (n = 564, n = 981 controls, p = 2.07×10⁻³, OR = 0.79 [95% CI 0.67–0.92]), France 1 (n = 433 cases, n = 395 controls, p = 3.73×10⁻³, OR = 0.74 [95% CI 0.60–0.91]), and France 2 (n = 249 cases, n = 380 controls, p = 2.26×10⁻⁴, OR = 0.63 [95% CI 0.50–0.81]). The combined analysis of all four studies including a total of n = 1,910 cases and n = 3,630 controls showed highly significant evidence for association between rs1739843 and idiopathic DCM (p = 5.28×10⁻¹³, OR = 0.72 [95% CI 0.65–0.78]). None of the other three SNPs showed significant results in the replication stage.This finding of the HSPB7 gene from a genetic search for idiopathic DCM using a large SNP panel underscores the influence of common polymorphisms on DCM susceptibility.     

Frequent and simultaneous epigenetic inactivation of TP53 pathway genes in acute lymphoblastic leukemia

Aberrant DNA methylation is one of the most frequent alterations in patients with Acute Lymphoblastic Leukemia (ALL). Using methylation bead arrays we analyzed the methylation status of 807 genes implicated in cancer in a group of ALL samples at diagnosis (n = 48). We found that 154 genes were methylated in more than 10% of ALL samples. Interestingly, the expression of 13 genes implicated in the TP53 pathway was downregulated by hypermethylation. Direct or indirect activation of TP53 pathway with 5-aza-2′-deoxycitidine, Curcumin or Nutlin-3 induced an increase in apoptosis of ALL cells. The results obtained with the initial group of 48 patients was validated retrospectively in a second cohort of 200 newly diagnosed ALL patients. Methylation of at least 1 of the 13 genes implicated in the TP53 pathway was observed in 78% of the patients, which significantly correlated with a higher relapse (p = 0.001) and mortality (p<0.001) rate being an independent prognostic factor for disease-free survival (DFS) (p = 0.006) and overall survival (OS) (p = 0.005) in the multivariate analysis. All these findings indicate that TP53 pathway is altered by epigenetic mechanisms in the majority of ALL patients and correlates with prognosis. Treatments with compounds that may reverse the epigenetic abnormalities or activate directly the p53 pathway represent a new therapeutic alternative for patients with ALL.     

Integrated genetic and epigenetic analysis identifies haplotype-specific methylation in the FTO type 2 diabetes and obesity susceptibility locus

Recent multi-dimensional approaches to the study of complex disease have revealed powerful insights into how genetic and epigenetic factors may underlie their aetiopathogenesis. We examined genotype-epigenotype interactions in the context of Type 2 Diabetes (T2D), focussing on known regions of genomic susceptibility. We assayed DNA methylation in 60 females, stratified according to disease susceptibility haplotype using previously identified association loci. CpG methylation was assessed using methylated DNA immunoprecipitation on a targeted array (MeDIP-chip) and absolute methylation values were estimated using a Bayesian algorithm (BATMAN). Absolute methylation levels were quantified across LD blocks, and we identified increased DNA methylation on the FTO obesity susceptibility haplotype, tagged by the rs8050136 risk allele A (p = 9.40×10⁻⁴, permutation p = 1.0×10⁻³). Further analysis across the 46 kb LD block using sliding windows localised the most significant difference to be within a 7.7 kb region (p = 1.13×10⁻⁷). Sequence level analysis, followed by pyrosequencing validation, revealed that the methylation difference was driven by the co-ordinated phase of CpG-creating SNPs across the risk haplotype. This 7.7 kb region of haplotype-specific methylation (HSM), encapsulates a Highly Conserved Non-Coding Element (HCNE) that has previously been validated as a long-range enhancer, supported by the histone H3K4me1 enhancer signature. This study demonstrates that integration of Genome-Wide Association (GWA) SNP and epigenomic DNA methylation data can identify potential novel genotype-epigenotype interactions within disease-associated loci, thus providing a novel route to aid unravelling common complex diseases.     

Correlation of global and gene-specific DNA methylation in maternal-infant pairs

The inheritance of DNA methylation patterns is a popular theory to explain the influence of parental genetic and environmental factors on the phenotype of their offspring but few studies have examined this relationship in humans. Using 120 paired maternal-umbilical cord blood samples randomly selected from a prospective birth cohort in Bangladesh, we quantified DNA methylation by pyrosequencing seven CpG positions in the promoter region of p16, four CpG positions in the promoter region of p53, LINE-1 and Alu. Positive correlations were observed between maternal and umbilical cord blood at p16, LINE-1, and Alu but not p53. Multiple linear regression models observed a significant association between maternal and umbilical cord blood at LINE-1 and Alu (LINE-1: β = 0.63, p<0.0001; Alu: β = 0.28, p = 0.009). After adjusting for multiple comparisons, maternal methylation of p16 at position 4 significantly predicted methylation at the same position in umbilical cord blood (β = 0.43, p = <0.0001). These models explained 48%, 5% and 16% of the observed variability in umbilical cord %5mC for LINE-1, Alu and p16 at position 4, respectively. These results suggest that DNA methylation in maternal blood was correlated with her offspring at LINE-1, Alu, and p16 but not p53. Additional studies are needed to confirm whether these observed associations were due to the inheritance of epigenetic events or the shared environment between mother and fetus. Future studies should also use a multi-generational family-based design that would quantify both maternal and paternal contributions to DNA methylation in offspring across more than one generation.     

Case-control study of vitamin D, dickkopf homolog 1 (DKK1) gene methylation, VDR gene polymorphism and the risk of colon adenoma in African Americans

Background

There are sparse data on genetic, epigenetic and vitamin D exposure in African Americans (AA) with colon polyp. Consequently, we evaluated serum 25(OH) D levels, vitamin D receptor (VDR) polymorphisms and the methylation status of the tumor suppressor gene dickkopf homolog 1 (DKK1) as risk factors for colon polyp in this population.

Methods

The case-control study consisted of 93 patients with colon polyp (cases) and 187 healthy individuals (controls) at Howard University Hospital. Serum levels of 25(OH)D (including D3, D2, and total) were measured by liquid chromatography-mass spectrometry. DNA analysis focused on 49 single nucleotide polymorphisms (SNPs) in the VDR gene. Promoter methylation analysis of DKK1 was also performed. The resulting data were processed in unadjusted and multivariable logistic regression analyses.

Results

Cases and controls differed in vitamin D status (D₃<50 nmol/L: Median of 35.5 in cases vs. 36.8 in controls nmol/L; P = 0.05). Low levels of 25(OH)D₃ (<50 nmol/L) were observed in 86% of cases and 68% of controls and it was associated with higher risks of colon polyp (odds ratio of 2.7, 95% confidence interval 1.3–3.4). The SNP analysis showed no association between 46 VDR polymorphisms and colon polyp. The promoter of the DKK1 gene was unmethylated in 96% of the samples.

Conclusion

We found an inverse association between serum 25(OH)D₃ and colon polyp in AAs. VDR SNPs and DKK1 methylation were not associated with colon polyp. Vitamin D levels may in part explain the higher incidence of polyp in AAs.     

Association between Common Germline Genetic Variation in 94 Candidate Genes or Regions and Risks of Invasive Epithelial Ovarian Cancer

Background

Recent studies have identified several single nucleotide polymorphisms (SNPs) in the population that are associated with variations in the risks of many different diseases including cancers such as breast, prostate and colorectal. For ovarian cancer, the known highly penetrant susceptibility genes (BRCA1 and BRCA2) are probably responsible for only 40% of the excess familial ovarian cancer risks, suggesting that other susceptibility genes of lower penetrance exist.

Methods

We have taken a candidate approach to identifying moderate risk susceptibility alleles for ovarian cancer. To date, we have genotyped 340 SNPs from 94 candidate genes or regions, in up to 1,491 invasive epithelial ovarian cancer cases and 3,145 unaffected controls from three different population based studies from the UK, Denmark and USA.

Results

After adjusting for population stratification by genomic control, 18 SNPs (5.3%) were significant at the 5% level, and 5 SNPs (1.5%) were significant at the 1% level. The most significant association was for the SNP rs2107425, located on chromosome 11p15.5, which has previously been identified as a susceptibility allele for breast cancer from a genome wide association study (P-trend = 0.0012). When SNPs/genes were stratified into 7 different pathways or groups of validation SNPs, the breast cancer associated SNPs were the only group of SNPs that were significantly associated with ovarian cancer risk (P-heterogeneity = 0.0003; P-trend = 0.0028; adjusted (for population stratification) P-trend = 0.006). We did not find statistically significant associations when the combined data for all SNPs were analysed using an admixture maximum likelihood (AML) experiment-wise test for association (P-heterogeneity = 0.051; P-trend = 0.068).

Conclusion

These data suggest that a proportion of the SNPs we evaluated were associated with ovarian cancer risk, but that the effect sizes were too small to detect associations with individual SNPs.     

Frequent and distinct aberrations of DNA methylation patterns in fibrolamellar carcinoma of the liver

Background

Gene silencing due to aberrant DNA methylation is a frequent event in hepatocellular carcinoma (HCC) and also in hepatocellular adenoma (HCA). However, very little is known about epigenetic defects in fibrolamellar carcinoma (FLC), a rare variant of hepatocellular carcinoma that displays distinct clinical and morphological features.

Methodology/Principal Findings

We analyzed the methylation status of the APC, CDH1, cyclinD2, GSTπ1, hsa-mir-9-1, hsa-mir-9-2, and RASSF1A gene in a series of 15 FLC and paired normal liver tissue specimens by quantitative high-resolution pyrosequencing. Results were compared with common HCC arising in non-cirrhotic liver (n = 10). Frequent aberrant hypermethylation was found for the cyclinD2 (19%) and the RASSF1A (38%) gene as well as for the microRNA genes mir-9-1 (13%) and mir-9-2 (33%). In contrast to common HCC the APC and CDH1 (E-cadherin) genes were found devoid of any DNA methylation in FLC, whereas the GSTπ1 gene showed comparable DNA methylation in tumor and surrounding tissue at a moderate level. Changes in global DNA methylation level were measured by analyzing methylation status of the highly repetitive LINE-1 sequences. No evidence of global hypomethylation could be found in FLCs, whereas HCCs without cirrhosis showed a significant reduction in global methylation level as described previously.

Conclusions

FLCs display frequent and distinct gene-specific hypermethylation in the absence of significant global hypomethylation indicating that these two epigenetic aberrations are induced by different pathways and that full-blown malignancy can develop in the absence of global loss of DNA methylation. Only quantitative DNA methylation detection methodology was able to identify these differences.     

Genetic variants in TGF-β pathway are associated with ovarian cancer risk

The transforming growth factor-β (TGF-β) signaling pathway is involved in a diverse array of cellular processes responsible for tumorigenesis. In this case-control study, we applied a pathway-based approach to evaluate single-nucleotide polymorphisms (SNPs) in the TGF-β signaling pathway as predictors of ovarian cancer risk. We systematically genotyped 218 SNPs from 21 genes in the TGF-β signaling pathway in 417 ovarian cancer cases and 417 matched control subjects. We analyzed the associations of these SNPs with ovarian cancer risk, performed haplotype analysis and identified potential cumulative effects of genetic variants. We also performed analysis to identify higher-order gene-gene interactions influencing ovarian cancer risk. Individual SNP analysis showed that the most significant SNP was SMAD6: rs4147407, with an adjusted odds ratio (OR) of 1.60 (95% confidence interval [CI], 1.14–2.24, P = 0.0066). Cumulative genotype analysis of 13 SNPs with significant main effects exhibited a clear dose-response trend of escalating risk with increasing number of unfavorable genotypes. In gene-based analysis, SMAD6 was identified as the most significant gene associated with ovarian cancer risk. Haplotype analysis further revealed that two haplotype blocks within SMAD6 were significantly associated with decreased ovarian cancer risk, as compared to the most common haplotype. Gene-gene interaction analysis further categorized the study population into subgroups with different ovarian cancer risk. Our findings suggest that genetic variants in the TGF-β signaling pathway are associated with ovarian cancer risk and may facilitate the identification of high-risk subgroups in the general population.     

Association of the sirtuin and mitochondrial uncoupling protein genes with carotid plaque

Objective

Sirtuins (SIRTs) and mitochondrial uncoupling proteins (UCPs) have been implicated in cardiovascular diseases through the control of reactive oxygen species production. This study sought to investigate the association between genetic variants in the SIRT and UCP genes and carotid plaque.

Methods

In a group of 1018 stroke-free subjects from the Northern Manhattan Study with high-definition carotid ultrasonography and genotyping, we investigated the associations of 85 single nucleotide polymorphisms (SNPs) in the 11 SIRT and UCP genes with the presence and number of carotid plaques, and evaluated interactions of SNPs with sex, smoking, diabetes and hypertension as well as interactions between SNPs significantly associated with carotid plaque.

Results

Overall, 60% of subjects had carotid plaques. After adjustment for demographic and vascular risk factors, T-carriers of the SIRT6 SNP rs107251 had an increased risk for carotid plaque (odds ratio, OR = 1.71, 95% CI = 1.23–2.37, Bonferroni-corrected p = 0.03) and for a number of plaques (rate ratio, RR = 1.31, 1.18–1.45, Bonferroni-corrected p = 1.4×10⁻⁵), whereas T-carriers of the UCP5 SNP rs5977238 had an decreased risk for carotid plaque (OR = 0.49, 95% CI = 0.32–0.74, Bonferroni-corrected p = 0.02) and plaque number (RR = 0.64, 95% CI = 0.52–0.78, Bonferroni-corrected p = 4.9×10⁻⁴). Some interactions with a nominal p≤0.01 were found between sex and SNPs in the UCP1 and UCP3 gene; between smoking, diabetes, hypertension and SNPs in UCP5 and SIRT5; and between SNPs in the UCP5 gene and the UCP1, SIRT1, SIRT3, SIRT5, and SIRT6 genes in association with plaque phenotypes.

Conclusion

We observed significant associations between genetic variants in the SIRT6 and UCP5 genes and atherosclerotic plaque. We also found potential effect modifications by sex, smoking and vascular risk factors of the SIRT/UCP genes in the associations with atherosclerotic plaque. Further studies are needed to validate our observations.     

A genome-wide association search for type 2 diabetes genes in African Americans

African Americans are disproportionately affected by type 2 diabetes (T2DM) yet few studies have examined T2DM using genome-wide association approaches in this ethnicity. The aim of this study was to identify genes associated with T2DM in the African American population. We performed a Genome Wide Association Study (GWAS) using the Affymetrix 6.0 array in 965 African-American cases with T2DM and end-stage renal disease (T2DM-ESRD) and 1029 population-based controls. The most significant SNPs (n = 550 independent loci) were genotyped in a replication cohort and 122 SNPs (n = 98 independent loci) were further tested through genotyping three additional validation cohorts followed by meta-analysis in all five cohorts totaling 3,132 cases and 3,317 controls. Twelve SNPs had evidence of association in the GWAS (P<0.0071), were directionally consistent in the Replication cohort and were associated with T2DM in subjects without nephropathy (P<0.05). Meta-analysis in all cases and controls revealed a single SNP reaching genome-wide significance (P<2.5×10⁻⁸). SNP rs7560163 (P = 7.0×10⁻⁹, OR (95% CI) = 0.75 (0.67–0.84)) is located intergenically between RND3 and RBM43. Four additional loci (rs7542900, rs4659485, rs2722769 and rs7107217) were associated with T2DM (P<0.05) and reached more nominal levels of significance (P<2.5×10⁻⁵) in the overall analysis and may represent novel loci that contribute to T2DM. We have identified novel T2DM-susceptibility variants in the African-American population. Notably, T2DM risk was associated with the major allele and implies an interesting genetic architecture in this population. These results suggest that multiple loci underlie T2DM susceptibility in the African-American population and that these loci are distinct from those identified in other ethnic populations.     

Polymorphisms of homologous recombination genes and clinical outcomes of non-small cell lung cancer patients treated with definitive radiotherapy

The repair of DNA double-strand breaks (DSBs) is the major mechanism to maintain genomic stability in response to irradiation. We hypothesized that genetic polymorphisms in DSB repair genes may affect clinical outcomes among non-small cell lung cancer (NSCLC) patients treated with definitive radio(chemo)therapy. We genotyped six potentially functional single nucleotide polymorphisms (SNPs) (i.e., RAD51 −135G>C/rs1801320 and −172G>T/rs1801321, XRCC2 4234G>C/rs3218384 and R188H/rs3218536 G>A, XRCC3 T241M/rs861539 and NBN E185Q/rs1805794) and estimated their associations with overall survival (OS) and radiation pneumonitis (RP) in 228 NSCLC patients. We found a predictive role of RAD51 −135G>C SNP in RP development (adjusted hazard ratio [HR] = 0.52, 95% confidence interval [CI], 0.31–0.86, P = 0.010 for CG/CC vs. GG). We also found that RAD51 −135G>C and XRCC2 R188H SNPs were independent prognostic factors for overall survival (adjusted HR = 1.70, 95% CI, 1.14–2.62, P = 0.009 for CG/CC vs. GG; and adjusted HR = 1.70; 95% CI, 1.02–2.85, P = 0.043 for AG vs. GG, respectively) and that the SNP-survival association was most pronounced in the presence of RP. Our study suggests that HR genetic polymorphisms, particularly RAD51 −135G>C, may influence overall survival and radiation pneumonitis in NSCLC patients treated with definitive radio(chemo)therapy. Large studies are needed to confirm our findings.     

Genetic and non-genetic influences during pregnancy on infant global and site specific DNA methylation: role for folate gene variants and vitamin B12

Inter-individual variation in patterns of DNA methylation at birth can be explained by the influence of environmental, genetic and stochastic factors. This study investigates the genetic and non-genetic determinants of variation in DNA methylation in human infants. Given its central role in provision of methyl groups for DNA methylation, this study focuses on aspects of folate metabolism. Global (LUMA) and gene specific (IGF2, ZNT5, IGFBP3) DNA methylation were quantified in 430 infants by Pyrosequencing®. Seven polymorphisms in 6 genes (MTHFR, MTRR, FOLH1, CβS, RFC1, SHMT) involved in folate absorption and metabolism were analysed in DNA from both infants and mothers. Red blood cell folate and serum vitamin B₁₂ concentrations were measured as indices of vitamin status. Relationships between DNA methylation patterns and several covariates viz. sex, gestation length, maternal and infant red cell folate, maternal and infant serum vitamin B₁₂, maternal age, smoking and genotype were tested. Length of gestation correlated positively with IGF2 methylation (rho = 0.11, p = 0.032) and inversely with ZNT5 methylation (rho = −0.13, p = 0.017). Methylation of the IGFBP3 locus correlated inversely with infant vitamin B₁₂ concentration (rho = −0.16, p = 0.007), whilst global DNA methylation correlated inversely with maternal vitamin B₁₂ concentrations (rho = 0.18, p = 0.044). Analysis of common genetic variants in folate pathway genes highlighted several associations including infant MTRR 66G>A genotype with DNA methylation (χ² = 8.82, p = 0.003) and maternal MTHFR 677C>T genotype with IGF2 methylation (χ² = 2.77, p = 0.006). These data support the hypothesis that both environmental and genetic factors involved in one-carbon metabolism influence DNA methylation in infants. Specifically, the findings highlight the importance of vitamin B₁₂ status, infant MTRR genotype and maternal MTHFR genotype, all of which may influence the supply of methyl groups for DNA methylation. In addition, gestational length appears to be an important determinant of infant DNA methylation patterns.     

Genetic variants of TSLP and asthma in an admixed urban population

Background

Thymic stromal lymphopoietin (TSLP), an IL7-like cytokine produced by bronchial epithelial cells is upregulated in asthma and induces dendritic cell maturation supporting a Th2 response. Environmental pollutants, including tobacco smoke and diesel exhaust particles upregulate TSLP suggesting that TSLP may be an interface between environmental pollution and immune responses in asthma. Since asthma is prevalent in urban communities, variants in the TSLP gene may be important in asthma susceptibility in these populations.

Objectives

To determine whether genetic variants in TSLP are associated with asthma in an urban admixed population.

Methodology and Main Results

Ten tag-SNPs in the TSLP gene were analyzed for association with asthma using 387 clinically diagnosed asthmatic cases and 212 healthy controls from an urban admixed population. One SNP (rs1898671) showed nominally significant association with asthma (odds ratio (OR) = 1.50; 95% confidence interval (95% CI): 1.09–2.05, p = 0.01) after adjusting for age, BMI, income, education and population stratification. Association results were consistent using two different approaches to adjust for population stratification. When stratified by smoking status, the same SNP showed a significantly increased risk associated with asthma in ex-smokers (OR = 2.00, 95% CI: 1.04–3.83, p = 0.04) but not significant in never-smokers (OR = 1.34; 95% CI: 0.93–1.94, p = 0.11). Haplotype-specific score test indicated that an elevated risk for asthma was associated with a specific haplotype of TSLP involving SNP rs1898671 (OR = 1.58, 95% CI: 1.10–2.27, p = 0.01). Association of this SNP with asthma was confirmed in an independent large population-based cohort consortium study (OR = 1.15, 95% CI: 1.07–1.23, p = 0.0003) and the results stratified by smoking status were also validated (ex-smokers: OR = 1.21, 95% CI: 1.08–1.34, p = 0.003; never-smokers: OR = 1.06, 95% CI: 0.94–1.17, p = 0.33).

Conclusions

Genetic variants in TSLP may contribute to asthma susceptibility in admixed urban populations with a gene and environment interaction.     

Genetic determinants for body iron store and type 2 diabetes risk in US men and women

Background

High body iron store has been associated with an increased risk of type 2 diabetes (T2D); it remains unknown whether the genetic variants related to body iron status affect T2D risk. We aimed at comprehensively investigating the associations between the genetic variants related to body iron status and the T2D risk.

Methodology/Principal Findings

Six common SNPs related to body iron status from recent genome-wide association (GWA) studies were determined in the Nurses’ Health Study (NHS; 1,467 diabetic cases and 1,754 controls) and the Health Professionals Follow-up Study (HPFS; 1,124, diabetic cases and 1,298 controls). Plasma levels of ferritin, soluble transferrin receptor (sTfR), and transferrin were measured in NHS. Significant associations were observed for loci in TPMRSS6 with sTfR (P = 3.47×10⁻⁶), TF with transferrin (P = 0.0002 to 1.72×10⁻¹⁰); and HFE with ferritin (P = 0.017 to 1.6×10⁻⁸), sTfR (P = 0.007 to 7.9×10⁻⁶), and transferrin (P = 0.006 to 0.0007). The six SNPs together explained 5.7%, 2.7%, and 13.3% of the variation in plasma levels of ferritin, sTfR, and transferrin. After adjustment for the conventional risk factors, the T allele of SNP rs855791 in the TPMRSS6 gene was significantly associated with a 19% decreased risk of T2D (OR = 0.81; 95% CI = 0.66–0.98; P = 0.03) in men. Multiple tests attenuated this significant association to null. No associations were observed in women. SNPs at HFE and TF were not associated with diabetes risk in either sex. Dietary iron intake did not modify the associations of the newly identified loci with diabetes risk.

Conclusions/Significance

The newly identified iron-related SNP rs855791 in TPMRSS6 was nominally associated with a decreased risk of T2D in men but not in women. The apparent differences by gender warrant further study.