Exonuclease 1 (EXO1) gene variation and melanoma risk

ObjectiveDNA repair pathway genes play an important role in maintaining genomic integrity and protecting against cancer development. This study aimed to identify novel SNPs in the DNA repair-related genes associated with melanoma risk from a genome-wide association study (GWAS).MethodsA total of 8422 SNPs from the 165 DNA repair-related genes were extracted from a GWAS of melanoma risk, including 494 cases and 5628 controls from the Nurses’ Health Study (NHS) and the Health Professionals Follow-up Study (HPFS). We further replicated the top SNPs in a GWAS of melanoma risk from the MD Anderson Cancer Center (1804 cases and 1026 controls).ResultsA total of 3 SNPs with P value <0.001 were selected for in silico replication. One SNP was replicated: rs3902093 [A] in EXO1 promoter region (P_discovery = 6.6 × 10^?4, P_replication = 0.039, P_joint = 2.5 × 10^?4; OR_joint = 0.80, 95% CI: 0.71, 0.90). This SNP was associated with the expression of the EXO1; carriers of the A allele showed lower expression (P = 0.002).ConclusionOur study found that a promoter region SNP in the editing and processing nucleases gene EXO1 was associated with decreased expression of EXO1 and decreased melanoma risk. Further studies are warranted to validate this association and to investigate the potential mechanisms.     

Common genetic variants associated with lipid profiles in a Chinese pediatric population

Genome-wide association (GWA) studies have identified many candidate genes that are associated with blood lipid and lipoprotein concentrations. In this study, we want to know whether the results from European for lipid-related single-nucleotide polymorphisms (SNPs) are generalizable to Chinese children. We genotyped seven SNPs in Chinese school-age children (n = 3,503) and assessed the associations of these SNPs with lipids profiles and dyslipidemia. After false discovery rate correction, of the seven SNPs, six (rs2144300, p ~ 9.30 × 10^?3; rs1260333, p ~ 6.20 × 10^?11; rs1260326, p ~ 8.73 × 10^?11; rs10105606, p ~ 0.010; rs1748195, p ~ 0.016 and rs964184, p ~ 2.33 × 10^?13) showed strong association with triglycerides. Three SNPs (rs1260333, p ~ 3.30 × 10^?3; rs1260326, p ~ 4.39 × 10^?3 and rs2954029, p ~ 6.36 × 10^?4) showed strong association with total cholesterol. Two SNPs (rs10105606, p ~ 6.66 × 10^?4 and rs1748195, p ~ 2.55 × 10^?3) showed strong association with high density lipoprotein cholesterol. Four SNPs (rs1260333, p ~ 0.017; rs1260326, p ~ 0.013; rs2954029, p ~ 1.09 × 10^?3 and rs964184, p ~ 5.51 × 10^?3) showed strong association with low density lipoprotein cholesterol. There were significant associations between rs1260333 (OR is 0.82, 95 % CI 0.74–0.92, p ~ 3.96 × 10^?4), rs1260326 (OR is 0.82, 95 % CI 0.74–0.92, p ~ 5.31 × 10^?4), and rs964184 (OR is 1.36, 95 % CI 1.20–1.55, p ~ 1.89 × 10^?6) and dyslipidemia. These SNPs generated strong combined effects on lipid profiles and dyslipidemia. Our study demonstrates that SNPs associated with lipids from European GWA studies also play roles in Chinese children, which broadened the understanding of lipids metabolism.     

Genome-wide association study in Han Chinese identifies three novel loci for human height

Human height is a complex genetic trait with high heritability but discovery efforts in Asian populations are limited. We carried out a meta-analysis of genome-wide association studies (GWAS) for height in 6,534 subjects with in silico replication of 1,881 subjects in Han Chinese. We identified three novel loci reaching the genome-wide significance threshold (P < 5 × 10^?8), which mapped in or near ZNF638 (rs12612930, P = 2.02 × 10^?10), MAML2 (rs11021504, P = 7.81 × 10^?9), and C18orf12 (rs11082671, P = 1.87 × 10^?8). We also confirmed two loci previously reported in European populations including CS (rs3816804, P = 2.63 × 10^?9) and CYP19A1 (rs3751599, P = 4.80 × 10^?10). In addition, we provided evidence supporting 35 SNPs identified by previous GWAS (P < 0.05). Our study provides new insights into the genetic determination of biological regulation of human height.     

Bayesian logistic regression in detection of gene–steroid interaction for cancer at <Emphasis Type="Italic">PDLIM5</Emphasis> locus

The PDZ and LIM domain 5 (PDLIM5) gene may play a role in cancer, bipolar disorder, major depression, alcohol dependence and schizophrenia; however, little is known about the interaction effect of steroid and PDLIM5 gene on cancer. This study examined 47 single-nucleotide polymorphisms (SNPs) within the PDLIM5 gene in the Marshfield sample with 716 cancer patients (any diagnosed cancer, excluding minor skin cancer) and 2848 noncancer controls. Multiple logistic regression model in PLINK software was used to examine the association of each SNP with cancer. Bayesian logistic regression in PROC GENMOD in SAS statistical software, ver. 9.4 was used to detect gene–steroid interactions influencing cancer. Single marker analysis using PLINK identified 12 SNPs associated with cancer (P < 0.05); especially, SNP rs6532496 revealed the strongest association with cancer (P = 6.84 × 10^?3); while the next best signal was rs951613 (P = 7.46 × 10^?3). Classic logistic regression in PROC GENMOD showed that both rs6532496 and rs951613 revealed strong gene–steroid interaction effects (OR = 2.18, 95% CI = 1.31?3.63 with P = 2.9 × 10^?3 for rs6532496 and OR = 2.07, 95% CI = 1.24 ?3.45 with P = 5.43 × 10^?3 for rs951613, respectively). Results from Bayesian logistic regression showed stronger interaction effects (OR = 2.26, 95% CI = 1.2 ?3.38 for rs6532496 and OR = 2.14, 95% CI = 1.14 ?3.2 for rs951613, respectively). All the 12 SNPs associated with cancer revealed significant gene–steroid interaction effects (P < 0.05); whereas 13 SNPs showed gene–steroid interaction effects without main effect on cancer. SNP rs4634230 revealed the strongest gene–steroid interaction effect (OR = 2.49, 95% CI = 1.5 ?4.13 with P = 4.0 × 10^?4 based on the classic logistic regression and OR = 2.59, 95% CI = 1.4 ?3.97 from Bayesian logistic regression; respectively). This study provides evidence of common genetic variants within the PDLIM5 gene and interactions between PLDIM5 gene polymorphisms and steroid use influencing cancer.     

Genetic variants in the JAK1 gene confer higher risk of Behcet’s disease with ocular involvement in Han Chinese

Recent surveys have identified SLC22A4, SLC22A5, RUNX1, JAK1 as susceptibility genes for various immune-related diseases. An association study was performed in 738 Behcet’s patients with ocular involvement and 1,873 controls using the iPLEX system method. The first-stage study for 30 SNPs showed that SNPs rs2780815, rs310241, rs3790532 in JAK1 were associated with Behcet’s disease in Han Chinese (Pc_{(Bonferroni correction)} = 0.022–7.7 × 10^?3). The G allele and AA genotype of SNP rs2834643 in RUNX1 (Pc = 0.041–1.75 × 10^?3), but none of the other SNPs, were associated with Behcet’s disease. Haplotype analysis for the SLC22A4, SLC22A5 genes showed an increased tendency for AGTCTGCCGC frequency in patients compared with controls; however, the significance was lost after Bonferroni correction (P = 0.004, Pc > 0.05). Subsequently, we further replicated the significantly associated SNPs using another independent cohort. Replication and combining studies showed that three SNPs rs2780815, rs310241, rs3790532 in JAK1, but not SNP rs2834643 in RUNX1, were consistently associated with Behcet’s disease (replication: Pc = 0.012–9.60 × 10^?4; combining: Pc = 0.030–1.90 × 10^?4). SNPs rs2780815, rs310241, rs3790532 were estimated to confer a population attributable risk of 35.0, 28.0, 27.0 %, respectively. We found a strong association between HLA-B51 with Behcet’s disease in Chinese Han population (P = 1.35 × 10^?73; OR = 5.15; 95 % CI 4.28–6.19). GMDR analysis showed that no gene–gene interaction was detectable between JAK1 and HLA-B51. Logistic analysis indicated that the JAK1 gene was an independent risk factor for Behcet’s disease (P > 0.05). Real-time PCR analysis showed that no difference on the expression of JAK1 in PBMCs or LPS-stimulated PBMCs between individuals with the different rs1762780815 genotypes studied (P > 0.05). In conclusion, this study suggests that JAK1, but not SLC22A4, SLC22A5 and RUNX1, contributes to the genetic susceptibility to Behcet’s disease with ocular involvement.     

ERAP1 genetic variations associated with HLA-B27 interaction and disease severity of syndesmophytes formation in Taiwanese ankylosing spondylitis

IntroductionAnkylosing spondylitis (AS) is a familial, heritable disease specified by syndesmophyte formation leading to an ankylosed spine. Endoplasmic reticulum aminopeptidase 1 (ERAP1) genetic variations have been widely proved to be associated with AS in several ethnic populations. The aim of this study was to investigate whether ERAP1 single nucleotide polymorphisms (SNPs) are associated with AS susceptibility and disease severity in Taiwanese.MethodsFour ERAP1 SNPs (rs27037, rs27980, rs27044 and rs30187) were genotyped in 797 Taiwanese AS patients and 1,150 healthy controls. Distributions of genotype and alleles were compared between AS patients and healthy controls, and among AS patients stratified by clinical parameters.ResultsThe SNP rs27037T allele appeared to be a risk factor for AS susceptibility (P = 5.5 × 10^-5, OR 1.30, 95% CI: 1.15 to 1.48; GT+TT vs. GG P = 9.3 × 10^-5, OR 1.49, 95% CI: 1.22 to 1.82). In addition, the coding SNP (cSNP) rs27044G allele (P = 1.5 × 10^-4, OR 1.28, 95% CI: 1.13 to 1.46; CG+GG vs. CC, P = 1.7 × 10^-3, OR 1.44, 95% CI: 1.15 to 1.81) and the cSNP rs30187T allele (P = 1.7 × 10^-3, OR 1.23, 95% CI: 1.08 to 1.40; CT+TT vs. CC P = 6.1 × 10^-3, OR 1.38, 95% CI: 1.10 to 1.74) were predisposing factors for AS. Notably, the rs27044G allele carriers (CG+GG vs. CC, P = 0.015, OR 1.59, 95% CI: 1.33 to 2.30) and rs30187T allele carriers (CT+TT vs. CC, P = 0.011, OR 1.63, 95% CI: 1.12 to 2.38) were susceptible to syndesmophyte formation in AS patients. Furthermore, two cSNPs (rs27044 and rs30187) strongly associated with HLA-B27 positivity in AS patients. Finally, the ERAP1 SNP haplotype TCG (rs27037T/rs27980C/rs27044G) is a major risk factor for AS (adjusted P <0.00001, OR 1.38, 95% CI: 1.12 to 1.58) in Taiwanese.ConclusionsThis study provides the first evidence of ERAP1 SNPs involving syndesmophyte formation. The interactions between ERAP1 SNPs and HLA-B27 play critical roles in pMHC I pathway processing contributing to the pathogenesis of AS in multiple populations.     

Suggestion of GLYAT gene underlying variation of bone size and body lean mass as revealed by a bivariate genome-wide association study

Bone and muscle, two major tissue types of musculoskeletal system, have strong genetic determination. Abnormality in bone and/or muscle may cause musculoskeletal diseases such as osteoporosis and sarcopenia. Bone size phenotypes (BSPs), such as hip bone size (HBS), appendicular bone size (ABS), are genetically correlated with body lean mass (mainly muscle mass). However, the specific genes shared by these phenotypes are largely unknown. In this study, we aimed to identify the specific genes with pleiotropic effects on BSPs and appendicular lean mass (ALM). We performed a bivariate genome-wide association study (GWAS) by analyzing ~690,000 SNPs in 1,627 unrelated Han Chinese adults (802 males and 825 females) followed by a replication study in 2,286 unrelated US Caucasians (558 males and 1,728 females). We identified 14 interesting single nucleotide polymorphisms (SNPs) that may contribute to variation of both BSPs and ALM, with p values <10^?6 in discovery stage. Among them, the association of three SNPs (rs2507838, rs7116722, and rs11826261) in/near GLYAT (glycine-N-acyltransferase) gene was replicated in US Caucasians, with p values ranging from 1.89 × 10^?3 to 3.71 × 10^?4 for ALM–ABS, from 5.14 × 10^?3 to 1.11 × 10^?2 for ALM–HBS, respectively. Meta-analyses yielded stronger association signals for rs2507838, rs7116722, and rs11826261, with pooled p values of 1.68 × 10^?8, 7.94 × 10^?8, 6.80 × 10^?8 for ALB–ABS and 1.22 × 10^?4, 9.85 × 10^?5, 3.96 × 10^?4 for ALM–HBS, respectively. Haplotype allele ATA based on these three SNPs was also associated with ALM–HBS and ALM–ABS in both discovery and replication samples. Interestingly, GLYAT was previously found to be essential to glucose metabolism and energy metabolism, suggesting the gene’s dual role in both bone development and muscle growth. Our findings, together with the prior biological evidence, suggest the importance of GLYAT gene in co-regulation of bone phenotypes and body lean mass.     

Six Novel Loci Associated with Circulating VEGF Levels Identified by a Meta-analysis of Genome-Wide Association Studies

Vascular endothelial growth factor (VEGF) is an angiogenic and neurotrophic factor, secreted by endothelial cells, known to impact various physiological and disease processes from cancer to cardiovascular disease and to be pharmacologically modifiable. We sought to identify novel loci associated with circulating VEGF levels through a genome-wide association meta-analysis combining data from European-ancestry individuals and using a dense variant map from 1000 genomes imputation panel. Six discovery cohorts including 13,312 samples were analyzed, followed by in-silico and de-novo replication studies including an additional 2,800 individuals. A total of 10 genome-wide significant variants were identified at 7 loci. Four were novel loci (5q14.3, 10q21.3, 16q24.2 and 18q22.3) and the leading variants at these loci were rs114694170 (MEF2C, P = 6.79x10^-13), rs74506613 (JMJD1C, P = 1.17x10^-19), rs4782371 (ZFPM1, P = 1.59x10^-9) and rs2639990 (ZADH2, P = 1.72x10^-8), respectively. We also identified two new independent variants (rs34528081, VEGFA, P = 1.52x10^-18; rs7043199, VLDLR-AS1, P = 5.12x10^-14) at the 3 previously identified loci and strengthened the evidence for the four previously identified SNPs (rs6921438, LOC100132354, P = 7.39x10^-1467; rs1740073, C6orf223, P = 2.34x10^-17; rs6993770, ZFPM2, P = 2.44x10^-60; rs2375981, KCNV2, P = 1.48x10^-100). These variants collectively explained up to 52% of the VEGF phenotypic variance. We explored biological links between genes in the associated loci using Ingenuity Pathway Analysis that emphasized their roles in embryonic development and function. Gene set enrichment analysis identified the ERK5 pathway as enriched in genes containing VEGF associated variants. eQTL analysis showed, in three of the identified regions, variants acting as both cis and trans eQTLs for multiple genes. Most of these genes, as well as some of those in the associated loci, were involved in platelet biogenesis and functionality, suggesting the importance of this process in regulation of VEGF levels. This work also provided new insights into the involvement of genes implicated in various angiogenesis related pathologies in determining circulating VEGF levels. The understanding of the molecular mechanisms by which the identified genes affect circulating VEGF levels could be important in the development of novel VEGF-related therapies for such diseases.     

Genetic susceptibility to obesity and diet intakes: association and interaction analyses in the Malmö Diet and Cancer Study

Gene–environment interactions need to be studied to better understand the obesity. We aimed at determining whether genetic susceptibility to obesity associates with diet intake levels and whether diet intakes modify the genetic susceptibility. In 29,480 subjects of the population-based Malmö Diet and Cancer Study (MDCS), we first assessed association between 16 genome-wide association studies identified obesity-related single-nucleotide polymorphisms (SNPs) with body mass index (BMI) and associated traits. We then conducted association analyses between a genetic risk score (GRS) comprising of 13 replicated SNPs and the individual SNPs, and relative dietary intakes of fat, carbohydrates, protein, fiber and total energy intake, as well as interaction analyses on BMI and associated traits among 26,107 nondiabetic MDCS participants. GRS associated strongly with increased BMI (P = 3.6 × 10^?34), fat mass (P = 6.3 × 10^?28) and fat-free mass (P = 1.3 × 10^?24). Higher GRS associated with lower total energy intake (P = 0.001) and higher intake of fiber (P = 2.3 × 10^?4). No significant interactions were observed between GRS and the studied dietary intakes on BMI or related traits. Of the individual SNPs, after correcting for multiple comparisons, NEGR1 rs2815752 associated with diet intakes and BDNF rs4923461 showed interaction with protein intake on BMI. In conclusion, our study does not provide evidence for a major role for macronutrient-, fiber- or total energy intake levels in modifying genetic susceptibility to obesity measured as GRS. However, our data suggest that the number of risk alleles as well as some of the individual obesity loci may have a role in regulation of food and energy intake and that some individual loci may interact with diet.     

Monoamine oxidase A gene polymorphisms and enzyme activity associated with risk of gout in Taiwan aborigines

Taiwanese aborigines have a high prevalence of hyperuricemia and gout. Uric acid levels and urate excretion have correlated with dopamine-induced glomerular filtration response. MAOs represent one of the major renal dopamine metabolic pathways. We aimed to identify the monoamine oxidase A (MAOA, Xp11.3) gene variants and MAO-A enzyme activity associated with gout risk. This study was to investigate the association between gout and the MAOA single-nucleotide polymorphisms (SNPs) rs5953210, rs2283725, and rs1137070 as well as between gout and the COMT SNPs rs4680 Val158Met for 374 gout cases and 604 controls. MAO-A activity was also measured. All three MAOA SNPs were significantly associated with gout. A synonymous MAOA SNP, rs1137070 Asp470Asp, located in exon 14, was associated with the risk of having gout (P = 4.0 × 10^?5, adjusted odds ratio 1.46, 95% confidence intervals [CI]: 1.11–1.91). We also showed that, when compared to individuals with the MAOA GAT haplotype, carriers of the AGC haplotype had a 1.67-fold (95% CI: 1.28–2.17) higher risk of gout. Moreover, we found that MAOA enzyme activity correlated positively with hyperuricemia and gout (P for trend = 2.00 × 10^?3 vs. normal control). We also found that MAOA enzyme activity by rs1137070 allele was associated with hyperuricemia and gout (P for trend = 1.53 × 10^?6 vs. wild-type allele). Thus, our results show that some MAOA alleles, which have a higher enzyme activity, predispose to the development of gout.     

Gene-environment interactions and obesity traits among postmenopausal African-American and Hispanic women in the Women’s Health Initiative SHARe Study

Genome-wide association studies (GWAS) of obesity measures have identified associations with single nucleotide polymorphisms (SNPs). However, no large-scale evaluation of gene-environment interactions has been performed. We conducted a search of gene-environment (G × E) interactions in post-menopausal African-American and Hispanic women from the Women’s Health Initiative SNP Health Association Resource GWAS study. Single SNP linear regression on body mass index (BMI) and waist-to-hip circumference ratio (WHR) adjusted for multidimensional-scaling-derived axes of ancestry and age was run in race-stratified data with 871,512 SNPs available from African-Americans (N = 8,203) and 786,776 SNPs from Hispanics (N = 3,484). Tests of G × E interaction at all SNPs for recreational physical activity (m h/week), dietary energy intake (kcal/day), alcohol intake (categorical), cigarette smoking years, and cigarette smoking (ever vs. never) were run in African-Americans and Hispanics adjusted for ancestry and age at interview, followed by meta-analysis of G × E interaction terms. The strongest evidence for concordant G × E interactions in African-Americans and Hispanics was for smoking and marker rs10133840 (Q statistic P = 0.70, beta = ?0.01, P = 3.81 × 10^?7) with BMI as the outcome. The strongest evidence for G × E interaction within a cohort was in African-Americans with WHR as outcome for dietary energy intake and rs9557704 (SNP × kcal = ?0.04, P = 2.17 × 10^?7). No results exceeded the Bonferroni-corrected statistical significance threshold.     

Genetic contribution to variation in DNA methylation at maternal smoking-sensitive loci in exposed neonates

Epigenome-wide DNA methylation association studies have identified highly replicable genomic loci sensitive to maternal smoking during gestation. The role of inter-individual genetic variation in influencing DNA methylation, leading to the possibility of confounding or bias of such associations, has not been assessed. We investigated whether the DNA methylation levels at the top 10 CpG sites previously associated with exposure to maternal smoking during gestation were associated with individual genetic variation at the genome-wide level. Genome-wide association tests between DNA methylation at the top 10 candidate CpG and genome-wide SNPs were performed in 736 case and control participants of the California Childhood Leukemia Study. Three of the strongest maternal-smoking sensitive CpG sites in newborns were significantly associated with SNPs located proximal to each gene: cg18146737 in the GFI1 gene with rs141819830 (P = 8.2×10^?44), cg05575921 in the AHRR gene with rs148405299 (P = 5.3×10^?10), and cg12803068 in the MYO1G gene with rs61087368 (P = 1.3×10^?18). For the GFI1 CpG cg18146737, the underlying genetic variation at rs141819830 confounded the association between maternal smoking and DNA methylation in our data (the regression coefficient changed from ?0.02 [P = 0.139] to ?0.03 [P = 0.015] after including the genotype). Our results suggest that further studies using DNA methylation at cg18146737, cg05575921, or cg12803068 that aim to assess exposure to maternal smoking during gestation should include genotype at the corresponding SNP. New methods are required for adequate and routine inclusion of genotypic influence on DNA methylation in epigenome-wide association studies to control for potential confounding.     

Identification of BACE2 as an avid ?-amyloid-degrading protease

Background

Glutathione S-transferase omega-1 and 2 genes (GSTO1, GSTO2), residing within an Alzheimer and Parkinson disease (AD and PD) linkage region, have diverse functions including mitigation of oxidative stress and may underlie the pathophysiology of both diseases. GSTO polymorphisms were previously reported to associate with risk and age-at-onset of these diseases, although inconsistent follow-up study designs make interpretation of results difficult. We assessed two previously reported SNPs, GSTO1 rs4925 and GSTO2 rs156697, in AD (3,493 ADs vs. 4,617 controls) and PD (678 PDs vs. 712 controls) for association with disease risk (case-controls), age-at-diagnosis (cases) and brain gene expression levels (autopsied subjects).

Results

We found that rs156697 minor allele associates with significantly increased risk (odds ratio = 1.14, p = 0.038) in the older ADs with age-at-diagnosis > 80 years. The minor allele of GSTO1 rs4925 associates with decreased risk in familial PD (odds ratio = 0.78, p = 0.034). There was no other association with disease risk or age-at-diagnosis. The minor alleles of both GSTO SNPs associate with lower brain levels of GSTO2 (p = 4.7 × 10^-11-1.9 × 10^-27), but not GSTO1. Pathway analysis of significant genes in our brain expression GWAS, identified significant enrichment for glutathione metabolism genes (p = 0.003).

Conclusion

These results suggest that GSTO locus variants may lower brain GSTO2 levels and consequently confer AD risk in older age. Other glutathione metabolism genes should be assessed for their effects on AD and other chronic, neurologic diseases.     

Paternal genetic variants and risk of obstructive heart defects: A parent-of-origin approach

Previous research on risk factors for obstructive heart defects (OHDs) focused on maternal and infant genetic variants, prenatal environmental exposures, and their potential interaction effects. Less is known about the role of paternal genetic variants or environmental exposures and risk of OHDs. We examined parent-of-origin effects in transmission of alleles in the folate, homocysteine, or transsulfuration pathway genes on OHD occurrence in offspring. We used data on 569 families of liveborn infants with OHDs born between October 1997 and August 2008 from the National Birth Defects Prevention Study to conduct a family-based case-only study. Maternal, paternal, and infant DNA were genotyped using an Illumina Golden Gate custom single nucleotide polymorphism (SNP) panel. Relative risks (RR), 95% confidence interval (CI), and likelihood ratio tests from log-linear models were used to estimate the parent-of-origin effect of 877 SNPs in 60 candidate genes in the folate, homocysteine, and transsulfuration pathways on the risk of OHDs. Bonferroni correction was applied for multiple testing. We identified 3 SNPs in the transsulfuration pathway and 1 SNP in the folate pathway that were statistically significant after Bonferroni correction. Among infants who inherited paternally-derived copies of the G allele for rs6812588 in the RFC1 gene, the G allele for rs1762430 in the MGMT gene, and the A allele for rs9296695 and rs4712023 in the GSTA3 gene, RRs for OHD were 0.11 (95% CI: 0.04, 0.29, P = 9.16x10^-7), 0.30 (95% CI: 0.17, 0.53, P = 9.80x10^-6), 0.34 (95% CI: 0.20, 0.57, P = 2.28x10^-5), and 0.34 (95% CI: 0.20, 0.58, P = 3.77x10^-5), respectively, compared to infants who inherited maternally-derived copies of the same alleles. We observed statistically significant decreased risk of OHDs among infants who inherited paternal gene variants involved in folate and transsulfuration pathways.     

Male-specific genetic effect on hypertension and metabolic disorders

Genetic risk factors for hypertension may have age or gender specificity and pleiotropic effects. This study aims to measure the risk of genetic and non-genetic factors in the occurrence of hypertension and related diseases, with consideration of potential confounding factors and age-gender stratification. A discovery set of 352,228 genotyped plus 1.8 million imputed single-nucleotide polymorphisms were analyzed for 2,886 hypertensive cases and 3,440 healthy controls obtained from two community-based cohorts in Korea, and selected gene variants were replicated in the Health Examinee cohort (665 cases and 1,285 controls). Genome-wide association analyses were conducted in 12 groups stratified by age and gender after adjusting for potential covariates under three genetic models. Age, rural area residence, body mass index, family history of hypertension, male gender, current alcohol drinking status, and current smoking status were significantly associated with hypertension (P = 4 × 10^?151 to 0.011). Five gene variants, rs11066280 (C12orf51), rs12229654 and rs3782889 (MYL2), rs2072134 (OAS3), rs2093395 (TREML2), and rs17249754 (ATP2B1), were found to be associated with hypertension mostly in men (P = 4.76 × 10^?14 to 4.46 × 10^?7 in the joint analysis); three SNPs (rs11066280, rs12229654, and rs3782889) remained significant after Bonferroni correction in an independent population. Three gene variants, rs12229654, rs17249754, and rs11066280, were significantly associated with metabolic disorders such as hyperlipidemia and diabetes (P = 0.00071 to 0.0097, respectively). Careful consideration of the potential confounding effects in future genome-wide association studies is necessary to uncover the genetic underpinnings of complex diseases.     

Genome-wide scan revealed that polymorphisms in the PNPLA3, SAMM50, and PARVB genes are associated with development and progression of nonalcoholic fatty liver disease in Japan

We examined the genetic background of nonalcoholic fatty liver disease (NAFLD) in the Japanese population, by performing a genome-wide association study (GWAS). For GWAS, 392 Japanese NAFLD subjects and 934 control individuals were analyzed. For replication studies, 172 NAFLD and 1,012 control subjects were monitored. After quality control, 261,540 single-nucleotide polymorphisms (SNPs) in autosomal chromosomes were analyzed using a trend test. Association analysis was also performed using multiple logistic regression analysis using genotypes, age, gender and body mass index (BMI) as independent variables. Multiple linear regression analyses were performed to evaluate allelic effect of significant SNPs on biochemical traits and histological parameters adjusted by age, gender, and BMI. Rs738409 in the PNPLA3 gene was most strongly associated with NAFLD after adjustment (P = 6.8 × 10^?14, OR = 2.05). Rs2896019, and rs381062 in the PNPLA3 gene, rs738491, rs3761472, and rs2143571 in the SAMM50 gene, rs6006473, rs5764455, and rs6006611 in the PARVB gene had also significant P values (<2.0 × 10^?10) and high odds ratios (1.84–2.02). These SNPs were found to be in the same linkage disequilibrium block and were associated with decreased serum triglycerides and increased aspartate aminotransferase (AST) and alanine aminotransferase (ALT) in NAFLD patients. These SNPs were associated with steatosis grade and NAFLD activity score (NAS). Rs738409, rs2896019, rs738491, rs6006473, rs5764455, and rs6006611 were associated with fibrosis. Polymorphisms in the SAMM50 and PARVB genes in addition to those in the PNPLA3 gene were observed to be associated with the development and progression of NAFLD.     

Effects of MTNR1B variants on fasting glucose levels in a Korean population

Fasting glucose level is the most basic and widely used indicator of diabetes. Several genome wide association studies have reported that the gene encoding melatonin receptor 1B (MTNR1B) exerts a major effect on serum fasting glucose levels. We tested for the association between single nucleotide polymorphisms (SNPs) in the MTNR1B gene and fasting glucose levels in a Korean population consisting of 8,229 subjects taken from two community-based cohorts. The mean age of the subjects in the study population was 51.9 years. For this study, we selected 363 SNPs located in the MTNR1B gene, which is located on chromosome 11. Multivariate linear regression models were used to test for genotypic effects on fasting glucose levels while adjusting for age and sex under an additive model. The MTNR1B SNP most highly associated with fasting glucose levels was rs10830962 (p=1.95×10^?5), followed by rs3847554 (p=3.16×10^?4). Replication of these initial findings is important to better understand the correlation between MTNR1B variations and their effects, especially in Asian populations.     

Metabolic heritability at birth: implications for chronic disease research

Recent genome-wide association studies of the adult human metabolome have identified genetic variants associated with relative levels of several acylcarnitines, which are important clinical correlates for chronic conditions such as type 2 diabetes and obesity. We have previously shown that these same metabolite levels are highly heritable at birth; however, no studies to our knowledge have examined genetic associations with these metabolites measured at birth. Here, we examine, in 743 newborns, 58 single nucleotide polymorphisms (SNPs) in 11 candidate genes previously associated with differing relative levels of short-chain acylcarnitines in adults. Six SNPs (rs2066938, rs3916, rs3794215, rs555404, rs558314, rs1799958) in the short-chain acyl-CoA dehydrogenase gene (ACADS) were associated with neonatal C4 levels. Most significant was the G allele of rs2066938, which was associated with significantly higher levels of C4 (P = 1.5 × 10^?29). This SNP explains 25 % of the variation in neonatal C4 levels, which is similar to the variation previously reported in adult C4 levels. There were also significant (P < 1 × 10^?4) associations between neonatal levels of C5-OH and SNPs in the solute carrier family 22 genes (SLC22A4 and SLC22A5) and the 3-methylcrotonyl-CoA carboxylase 1 gene (MCCC1). We have replicated, in newborns, SNP associations between metabolic traits and the ACADS and SLC22A4 genes observed in adults. This research has important implications not only for the identification of rare inborn errors of metabolism but also for personalized medicine and early detection of later life risks for chronic conditions.     

Gene-Gene Associations with the Susceptibility of Kawasaki Disease and Coronary Artery Lesions

Kawasaki disease (KD) is a systemic vasculitis primarily affecting children < 5 years old. Genes significantly associated with KD mostly involve cardiovascular, immune, and inflammatory responses. Recent studies have observed stronger associations for KD risk with multiple genes compared to individual genes. Therefore, we investigated whether gene combinations influenced KD susceptibility or coronary artery lesion (CAL) formation. We examined 384 single-nucleotide polymorphisms (SNPs) for 159 immune-related candidate genes in DNA samples from KD patients with CAL (n = 73), KD patients without CAL (n = 153), and cohort controls (n = 575). Individual SNPs were first assessed by univariate analysis (UVA) and multivariate analysis (MVA). We used multifactor dimensionality reduction (MDR) to examine individual SNPs in one-, two-, and three-locus best fit models. UVA identified 53 individual SNPs that were significantly associated with KD risk or CAL formation (p < 0.10), while 35 individual SNPs were significantly associated using MVA (p ≤ 0.05). Significant associations in MDR analysis were only observed for the two-locus models after permutation testing (p ≤ 0.05). In logistic regression, combined possession of PDE2A (rs341058) and CYFIP2 (rs767007) significantly increased KD susceptibility (OR = 3.54; p = 4.14 x 10⁻⁷), while combinations of LOC100133214 (rs2517892) and IL2RA (rs3118470) significantly increased the risk of CAL in KD patients (OR = 5.35; p = 7.46 x 10⁻⁵). Our results suggest varying gene-gene associations respectively predispose individuals to KD risk or its complications of CAL.