Genetic Architecture of Skin and Eye Color in an African-European Admixed Population

Variation in human skin and eye color is substantial and especially apparent in admixed populations, yet the underlying genetic architecture is poorly understood because most genome-wide studies are based on individuals of European ancestry. We study pigmentary variation in 699 individuals from Cape Verde, where extensive West African/European admixture has given rise to a broad range in trait values and genomic ancestry proportions. We develop and apply a new approach for measuring eye color, and identify two major loci (HERC2[OCA2] P = 2.3×10⁻⁶², SLC24A5 P = 9.6×10⁻⁹) that account for both blue versus brown eye color and varying intensities of brown eye color. We identify four major loci (SLC24A5 P = 5.4×10⁻²⁷, TYR P = 1.1×10⁻⁹, APBA2[OCA2] P = 1.5×10⁻⁸, SLC45A2 P = 6×10⁻⁹) for skin color that together account for 35% of the total variance, but the genetic component with the largest effect (∼44%) is average genomic ancestry. Our results suggest that adjacent cis-acting regulatory loci for OCA2 explain the relationship between skin and eye color, and point to an underlying genetic architecture in which several genes of moderate effect act together with many genes of small effect to explain ∼70% of the estimated heritability.     

A Large Scale Gene-Centric Association Study of Lung Function in Newly-Hired Female Cotton Textile Workers with Endotoxin Exposure

Background

Occupational exposure to endotoxin is associated with decrements in pulmonary function, but how much variation in this association is explained by genetic variants is not well understood.

Objective

We aimed to identify single nucleotide polymorphisms (SNPs) that are associated with the rate of forced expiratory volume in one second (FEV₁) decline by a large scale genetic association study in newly-hired healthy young female cotton textile workers.

Methods

DNA samples were genotyped using the Illumina Human CVD BeadChip. Change rate in FEV₁ was modeled as a function of each SNP genotype in linear regression model with covariate adjustment. We controlled the type 1 error in study-wide level by permutation method. The false discovery rate (FDR) and the family-wise error rate (FWER) were set to be 0.10 and 0.15 respectively.

Results

Two SNPs were found to be significant (P<6.29×10⁻⁵), including rs1910047 (P = 3.07×10⁻⁵, FDR = 0.0778) and rs9469089 (P = 6.19×10⁻⁵, FDR = 0.0967), as well as other eight suggestive (P<5×10⁻⁴) associated SNPs. Gene-gene and gene-environment interactions were also observed, such as rs1910047 and rs1049970 (P = 0.0418, FDR = 0.0895); rs9469089 and age (P = 0.0161, FDR = 0.0264). Genetic risk score analysis showed that the more risk loci the subjects carried, the larger the rate of FEV₁ decline occurred (P _trend = 3.01×10⁻¹⁸). However, the association was different among age subgroups (P = 7.11×10⁻⁶) and endotoxin subgroups (P = 1.08×10⁻²). Functional network analysis illustrates potential biological connections of all interacted genes.

Conclusions

Genetic variants together with environmental factors interact to affect the rate of FEV₁ decline in cotton textile workers.     

Insights into the Genetic Architecture of Early Stage Age-Related Macular Degeneration: A Genome-Wide Association Study Meta-Analysis

Genetic factors explain a majority of risk variance for age-related macular degeneration (AMD). While genome-wide association studies (GWAS) for late AMD implicate genes in complement, inflammatory and lipid pathways, the genetic architecture of early AMD has been relatively under studied. We conducted a GWAS meta-analysis of early AMD, including 4,089 individuals with prevalent signs of early AMD (soft drusen and/or retinal pigment epithelial changes) and 20,453 individuals without these signs. For various published late AMD risk loci, we also compared effect sizes between early and late AMD using an additional 484 individuals with prevalent late AMD. GWAS meta-analysis confirmed previously reported association of variants at the complement factor H (CFH) (peak P = 1.5×10⁻³¹) and age-related maculopathy susceptibility 2 (ARMS2) (P = 4.3×10⁻²⁴) loci, and suggested Apolipoprotein E (ApoE) polymorphisms (rs2075650; P = 1.1×10⁻⁶) associated with early AMD. Other possible loci that did not reach GWAS significance included variants in the zinc finger protein gene GLI3 (rs2049622; P = 8.9×10⁻⁶) and upstream of GLI2 (rs6721654; P = 6.5×10⁻⁶), encoding retinal Sonic hedgehog signalling regulators, and in the tyrosinase (TYR) gene (rs621313; P = 3.5×10⁻⁶), involved in melanin biosynthesis. For a range of published, late AMD risk loci, estimated effect sizes were significantly lower for early than late AMD. This study confirms the involvement of multiple established AMD risk variants in early AMD, but suggests weaker genetic effects on the risk of early AMD relative to late AMD. Several biological processes were suggested to be potentially specific for early AMD, including pathways regulating RPE cell melanin content and signalling pathways potentially involved in retinal regeneration, generating hypotheses for further investigation.     

Genetic Determinants of Long-Term Changes in Blood Lipid Concentrations: 10-Year Follow-Up of the GLACIER Study

Recent genome-wide meta-analyses identified 157 loci associated with cross-sectional lipid traits. Here we tested whether these loci associate (singly and in trait-specific genetic risk scores [GRS]) with longitudinal changes in total cholesterol (TC) and triglyceride (TG) levels in a population-based prospective cohort from Northern Sweden (the GLACIER Study). We sought replication in a southern Swedish cohort (the MDC Study; N = 2,943). GLACIER Study participants (N = 6,064) were genotyped with the MetaboChip array. Up to 3,495 participants had 10-yr follow-up data available in the GLACIER Study. The TC- and TG-specific GRSs were strongly associated with change in lipid levels (β = 0.02 mmol/l per effect allele per decade follow-up, P = 2.0×10⁻¹¹ for TC; β = 0.02 mmol/l per effect allele per decade follow-up, P = 5.0×10⁻⁵ for TG). In individual SNP analysis, one TC locus, apolipoprotein E (APOE) rs4420638 (β = 0.12 mmol/l per effect allele per decade follow-up, P = 2.0×10⁻⁵), and two TG loci, tribbles pseudokinase 1 (TRIB1) rs2954029 (β = 0.09 mmol/l per effect allele per decade follow-up, P = 5.1×10⁻⁴) and apolipoprotein A-I (APOA1) rs6589564 (β = 0.31 mmol/l per effect allele per decade follow-up, P = 1.4×10⁻⁸), remained significantly associated with longitudinal changes for the respective traits after correction for multiple testing. An additional 12 loci were nominally associated with TC or TG changes. In replication analyses, the APOE rs4420638, TRIB1 rs2954029, and APOA1 rs6589564 associations were confirmed (P≤0.001). In summary, trait-specific GRSs are robustly associated with 10-yr changes in lipid levels and three individual SNPs were strongly associated with 10-yr changes in lipid levels.     

Genome-Wide Joint Meta-Analysis of SNP and SNP-by-Smoking Interaction Identifies Novel Loci for Pulmonary Function

Genome-wide association studies have identified numerous genetic loci for spirometic measures of pulmonary function, forced expiratory volume in one second (FEV₁), and its ratio to forced vital capacity (FEV₁/FVC). Given that cigarette smoking adversely affects pulmonary function, we conducted genome-wide joint meta-analyses (JMA) of single nucleotide polymorphism (SNP) and SNP-by-smoking (ever-smoking or pack-years) associations on FEV₁ and FEV₁/FVC across 19 studies (total N = 50,047). We identified three novel loci not previously associated with pulmonary function. SNPs in or near DNER (smallest P _JMA =5.00×10⁻¹¹), HLA-DQB1 and HLA-DQA2 (smallest P _JMA =4.35×10⁻⁹), and KCNJ2 and SOX9 (smallest P _JMA =1.28×10⁻⁸) were associated with FEV₁/FVC or FEV₁ in meta-analysis models including SNP main effects, smoking main effects, and SNP-by-smoking (ever-smoking or pack-years) interaction. The HLA region has been widely implicated for autoimmune and lung phenotypes, unlike the other novel loci, which have not been widely implicated. We evaluated DNER, KCNJ2, and SOX9 and found them to be expressed in human lung tissue. DNER and SOX9 further showed evidence of differential expression in human airway epithelium in smokers compared to non-smokers. Our findings demonstrated that joint testing of SNP and SNP-by-environment interaction identified novel loci associated with complex traits that are missed when considering only the genetic main effects.     

Modeling of Environmental Effects in Genome-Wide Association Studies Identifies SLC2A2 and HP as Novel Loci Influencing Serum Cholesterol Levels

Genome-wide association studies (GWAS) have identified 38 larger genetic regions affecting classical blood lipid levels without adjusting for important environmental influences. We modeled diet and physical activity in a GWAS in order to identify novel loci affecting total cholesterol, LDL cholesterol, HDL cholesterol, and triglyceride levels. The Swedish (SE) EUROSPAN cohort (N _SE = 656) was screened for candidate genes and the non-Swedish (NS) EUROSPAN cohorts (N _NS = 3,282) were used for replication. In total, 3 SNPs were associated in the Swedish sample and were replicated in the non-Swedish cohorts. While SNP rs1532624 was a replication of the previously published association between CETP and HDL cholesterol, the other two were novel findings. For the latter SNPs, the p-value for association was substantially improved by inclusion of environmental covariates: SNP rs5400 (p _{SE,unadjusted} = 3.6×10⁻⁵, p _SE,adjusted = 2.2×10⁻⁶, p _{NS,unadjusted} = 0.047) in the SLC2A2 (Glucose transporter type 2) and rs2000999 (p _{SE,unadjusted} = 1.1×10⁻³, p _SE,adjusted = 3.8×10⁻⁴, p _{NS,unadjusted} = 0.035) in the HP gene (Haptoglobin-related protein precursor). Both showed evidence of association with total cholesterol. These results demonstrate that inclusion of important environmental factors in the analysis model can reveal new genetic susceptibility loci.     

Obesity-Related Genomic Loci Are Associated with Type 2 Diabetes in a Han Chinese Population

Background and Aims

Obesity is a well-known risk factor for type 2 diabetes. Genome-wide association studies have identified a number of genetic loci associated with obesity. The aim of this study is to examine the contribution of obesity-related genomic loci to type 2 diabetes in a Chinese population.

Methods

We successfully genotyped 18 obesity-related single nucleotide polymorphisms among 5338 type 2 diabetic patients and 4663 controls. Both individual and joint effects of these single nucleotide polymorphisms on type 2 diabetes and quantitative glycemic traits (assessing β-cell function and insulin resistance) were analyzed using logistic and linear regression models, respectively.

Results

Two single nucleotide polymorphisms near MC4R and GNPDA2 genes were significantly associated with type 2 diabetes before adjusting for body mass index and waist circumference (OR (95% CI) = 1.14 (1.06, 1.22) for the A allele of rs12970134, P = 4.75×10⁻⁴; OR (95% CI) = 1.10 (1.03, 1.17) for the G allele of rs10938397, P = 4.54×10⁻³). When body mass index and waist circumference were further adjusted, the association of MC4R with type 2 diabetes remained significant (P = 1.81×10⁻²) and that of GNPDA2 was attenuated (P = 1.26×10⁻¹), suggesting the effect of the locus including GNPDA2 on type 2 diabetes may be mediated through obesity. Single nucleotide polymorphism rs2260000 within BAT2 was significantly associated with type 2 diabetes after adjusting for body mass index and waist circumference (P = 1.04×10⁻²). In addition, four single nucleotide polymorphisms (near or within SEC16B, BDNF, MAF and PRL genes) showed significant associations with quantitative glycemic traits in controls even after adjusting for body mass index and waist circumference (all P values<0.05).

Conclusions

This study indicates that obesity-related genomic loci were associated with type 2 diabetes and glycemic traits in the Han Chinese population.     

New Variants Including ARG1 Polymorphisms Associated with C-Reactive Protein Levels Identified by Genome-Wide Association and Pathway Analysis

C-reactive protein (CRP) is a general marker of systemic inflammation and cardiovascular disease (CVD). The genetic contribution to differences in CRP levels remains to be explained, especially in non-European populations. Thus, the aim of this study was to identify genetic loci associated with CRP levels in Korean population. We performed genome-wide association studies (GWAS) using SNPs from 8,529 Korean individuals (7,626 for stage 1 and 903 for stage 2). We also performed pathway analysis. We identified a new genetic locus associated with CRP levels upstream of ARG1 gene (top significant SNP: rs9375813, P_meta = 2.85×10⁻⁸), which encodes a key enzyme of the urea cycle counteract the effects of nitric oxide, in addition to known CRP (rs7553007, P_meta = 1.72×10⁻¹⁶) and HNF1A loci (rs2259816, P_meta = 2.90×10⁻¹⁰). When we evaluated the associations between the CRP-related SNPs with cardiovascular disease phenotypes, rs9375813 (ARG1) showed a marginal association with hypertension (P = 0.0440). To identify more variants and pathways, we performed pathway analysis and identified six candidate pathways comprised of genes related to inflammatory processes and CVDs (CRP, HNF1A, PCSK6, CD36, and ABCA1). In addition to the previously reported loci (CRP, HNF1A, and IL6) in diverse ethnic groups, we identified novel variants in the ARG1 locus associated with CRP levels in Korean population and a number of interesting genes related to inflammatory processes and CVD through pathway analysis.     

Transporting Antitumor Drug Tamoxifen and Its Metabolites, 4-Hydroxytamoxifen and Endoxifen by Chitosan Nanoparticles

Synthetic and natural polymers are often used as drug delivery systems in vitro and in vivo. Biodegradable chitosan of different sizes were used to encapsulate antitumor drug tamoxifen (Tam) and its metabolites 4-hydroxytamoxifen (4-Hydroxytam) and endoxifen (Endox). The interactions of tamoxifen and its metabolites with chitosan 15, 100 and 200 KD were investigated in aqueous solution, using FTIR, fluorescence spectroscopic methods and molecular modeling. The structural analysis showed that tamoxifen and its metabolites bind chitosan via both hydrophilic and hydrophobic contacts with overall binding constants of K _tam-ch-15  = 8.7 (±0.5)×10³ M⁻¹, K _tam-ch-100  = 5.9 (±0.4)×10⁵ M⁻¹, K _tam-ch-200  = 2.4 (±0.4)×10⁵ M⁻¹ and K _{hydroxytam-ch-15}  = 2.6(±0.3)×10⁴ M⁻¹, K _{hydroxytam – ch-100}  = 5.2 (±0.7)×10⁶ M⁻¹ and K _{hydroxytam-ch-200}  = 5.1 (±0.5)×10⁵ M⁻¹, K _endox-ch-15  = 4.1 (±0.4)×10³ M⁻¹, K _endox-ch-100  = 1.2 (±0.3)×10⁶ M⁻¹ and K _endox-ch-200  = 4.7 (±0.5)×10⁵ M⁻¹ with the number of drug molecules bound per chitosan (n) 2.8 to 0.5. The order of binding is ch-100>200>15 KD with stronger complexes formed with 4-hydroxytamoxifen than tamoxifen and endoxifen. The molecular modeling showed the participation of polymer charged NH₂ residues with drug OH and NH₂ groups in the drug-polymer adducts. The free binding energies of −3.46 kcal/mol for tamoxifen, −3.54 kcal/mol for 4-hydroxytamoxifen and −3.47 kcal/mol for endoxifen were estimated for these drug-polymer complexes. The results show chitosan 100 KD is stronger carrier for drug delivery than chitosan-15 and chitosan-200 KD.     

Characterization of the major histocompatibility complex locus association with Beh?et’s disease in Iran

IntroductionThe aim of this study was to characterize the association of human leukocyte antigen (HLA) B alleles and major histocompatibility complex (MHC) single nucleotide polymorphisms (SNPs) with Behçet’s disease (BD) in an Iranian dataset.MethodsThe association of three SNPs in the MHC region previously identified as the most associated in high-density genotyping studies was tested in a case–control study on 973 BD patients and 825 controls from Iran, and the association of HLA-B alleles was tested in a subset of 681 patients and 414 controls.ResultsWe found that HLA-B*51 (P = 4.11 × 10⁻⁴¹, OR [95% CI] = 4.63[3.66-5.85]) and B*15 confer risk for BD (P = 2.83 × 10⁻², OR [95% CI] = 1.75[1.08-2.84]) in Iranian, and in B*51 negative individuals, only the B*15 allele is significantly associated with BD (P = 2.51 × 10⁻³, OR [95% CI] = 2.40[1.37-4.20]). rs76546355, formerly known as rs116799036, located between HLA-B and MICA (MHC class I polypeptide-related sequence A), demonstrated the same level of association with BD as HLA-B*51 (P_adj = 1.78 × 10⁻⁴⁶, OR [95% CI] = 5.46[4.21-7.09], and P_adj = 8.34 × 10⁻⁴⁸, OR [95% CI] = 5.44[4.20-7.05], respectively) in the HLA-B allelotyped subset, while rs2848713 was less associated (P_adj = 7.14 × 10⁻³⁵, OR [95% CI] = 3.73[2.97-4.69]) and rs9260997 was not associated (P_adj = 1.00 × 10⁻¹). Additionally, we found that B*51 genotype-phenotype correlations do not survive Bonferroni correction, while carriers of the rs76546355 risk allele predominate in BD cases with genital ulcers, positive pathergy test and positive BD family history (2.31 × 10⁻⁴ ≤ P ≤ 1.59 × 10⁻³).ConclusionsWe found that the HLA-B*51 allele and the rs76546355/rs116799036 MHC SNP are independent genetic risk factors for BD in Iranian, and that positivity for the rs76546355/rs116799036 risk allele, but not for B*51, does correlate with specific demographic characteristics or clinical manifestations in BD patients.

Electronic supplementary material

The online version of this article (doi:10.1186/s13075-015-0585-6) contains supplementary material, which is available to authorized users.     

Genome-Wide Association Study of Breast Cancer in the Japanese Population

Breast cancer is the most common malignancy among women in worldwide including Japan. Several studies have identified common genetic variants to be associated with the risk of breast cancer. Due to the complex linkage disequilibrium structure and various environmental exposures in different populations, it is essential to identify variants associated with breast cancer in each population, which subsequently facilitate the better understanding of mammary carcinogenesis. In this study, we conducted a genome-wide association study (GWAS) as well as whole-genome imputation with 2,642 cases and 2,099 unaffected female controls. We further examined 13 suggestive loci (P<1.0×10⁻⁵) using an independent sample set of 2,885 cases and 3,395 controls and successfully validated two previously-reported loci, rs2981578 (combined P-value of 1.31×10⁻¹², OR = 1.23; 95% CI = 1.16–.30) on chromosome 10q26 (FGFR2), rs3803662 (combined P-value of 2.79×10⁻¹¹, OR = 1.21; 95% CI = 1.15–.28) and rs12922061 (combined P-value of 3.97×10⁻¹⁰, OR = 1.23; 95% CI = 1.15–.31) on chromosome 16q12 (TOX3-LOC643714). Weighted genetic risk score on the basis of three significantly associated variants and two previously reported breast cancer associated loci in East Asian population revealed that individuals who carry the most risk alleles in category 5 have 2.2 times higher risk of developing breast cancer in the Japanese population than those who carry the least risk alleles in reference category 1. Although we could not identify additional loci associated with breast cancer, our study utilized one of the largest sample sizes reported to date, and provided genetic status that represent the Japanese population. Further local and international collaborative study is essential to identify additional genetic variants that could lead to a better, accurate prediction for breast cancer.     

Identification of Novel Genetic Loci Associated with Thyroid Peroxidase Antibodies and Clinical Thyroid Disease

Autoimmune thyroid diseases (AITD) are common, affecting 2-5% of the general population. Individuals with positive thyroid peroxidase antibodies (TPOAbs) have an increased risk of autoimmune hypothyroidism (Hashimoto''s thyroiditis), as well as autoimmune hyperthyroidism (Graves'' disease). As the possible causative genes of TPOAbs and AITD remain largely unknown, we performed GWAS meta-analyses in 18,297 individuals for TPOAb-positivity (1769 TPOAb-positives and 16,528 TPOAb-negatives) and in 12,353 individuals for TPOAb serum levels, with replication in 8,990 individuals. Significant associations (P<5×10⁻⁸) were detected at TPO-rs11675434, ATXN2-rs653178, and BACH2-rs10944479 for TPOAb-positivity, and at TPO-rs11675434, MAGI3-rs1230666, and KALRN-rs2010099 for TPOAb levels. Individual and combined effects (genetic risk scores) of these variants on (subclinical) hypo- and hyperthyroidism, goiter and thyroid cancer were studied. Individuals with a high genetic risk score had, besides an increased risk of TPOAb-positivity (OR: 2.18, 95% CI 1.68–2.81, P = 8.1×10⁻⁸), a higher risk of increased thyroid-stimulating hormone levels (OR: 1.51, 95% CI 1.26–1.82, P = 2.9×10⁻⁶), as well as a decreased risk of goiter (OR: 0.77, 95% CI 0.66–0.89, P = 6.5×10⁻⁴). The MAGI3 and BACH2 variants were associated with an increased risk of hyperthyroidism, which was replicated in an independent cohort of patients with Graves'' disease (OR: 1.37, 95% CI 1.22–1.54, P = 1.2×10⁻⁷ and OR: 1.25, 95% CI 1.12–1.39, P = 6.2×10⁻⁵). The MAGI3 variant was also associated with an increased risk of hypothyroidism (OR: 1.57, 95% CI 1.18–2.10, P = 1.9×10⁻³). This first GWAS meta-analysis for TPOAbs identified five newly associated loci, three of which were also associated with clinical thyroid disease. With these markers we identified a large subgroup in the general population with a substantially increased risk of TPOAbs. The results provide insight into why individuals with thyroid autoimmunity do or do not eventually develop thyroid disease, and these markers may therefore predict which TPOAb-positives are particularly at risk of developing clinical thyroid dysfunction.     

Insight in Genome-Wide Association of Metabolite Quantitative Traits by Exome Sequence Analyses

Metabolite quantitative traits carry great promise for epidemiological studies, and their genetic background has been addressed using Genome-Wide Association Studies (GWAS). Thus far, the role of less common variants has not been exhaustively studied. Here, we set out a GWAS for metabolite quantitative traits in serum, followed by exome sequence analysis to zoom in on putative causal variants in the associated genes. ¹H Nuclear Magnetic Resonance (¹H-NMR) spectroscopy experiments yielded successful quantification of 42 unique metabolites in 2,482 individuals from The Erasmus Rucphen Family (ERF) study. Heritability of metabolites were estimated by SOLAR. GWAS was performed by linear mixed models, using HapMap imputations. Based on physical vicinity and pathway analyses, candidate genes were screened for coding region variation using exome sequence data. Heritability estimates for metabolites ranged between 10% and 52%. GWAS replicated three known loci in the metabolome wide significance: CPS1 with glycine (P-value  = 1.27×10⁻³²), PRODH with proline (P-value  = 1.11×10⁻¹⁹), SLC16A9 with carnitine level (P-value  = 4.81×10⁻¹⁴) and uncovered a novel association between DMGDH and dimethyl-glycine (P-value  = 1.65×10⁻¹⁹) level. In addition, we found three novel, suggestively significant loci: TNP1 with pyruvate (P-value  = 1.26×10⁻⁸), KCNJ16 with 3-hydroxybutyrate (P-value  = 1.65×10⁻⁸) and 2p12 locus with valine (P-value  = 3.49×10⁻⁸). Exome sequence analysis identified potentially causal coding and regulatory variants located in the genes CPS1, KCNJ2 and PRODH, and revealed allelic heterogeneity for CPS1 and PRODH. Combined GWAS and exome analyses of metabolites detected by high-resolution ¹H-NMR is a robust approach to uncover metabolite quantitative trait loci (mQTL), and the likely causative variants in these loci. It is anticipated that insight in the genetics of intermediate phenotypes will provide additional insight into the genetics of complex traits.     

Polymorphism Located between CPT1B and CHKB,and HLA-DRB1*1501-DQB1*0602 Haplotype Confer Susceptibility to CNS Hypersomnias (Essential Hypersomnia)

Background

SNP rs5770917 located between CPT1B and CHKB, and HLA-DRB1*1501-DQB1*0602 haplotype were previously identified as susceptibility loci for narcolepsy with cataplexy. This study was conducted in order to investigate whether these genetic markers are associated with Japanese CNS hypersomnias (essential hypersomnia: EHS) other than narcolepsy with cataplexy.

Principal Findings

EHS was significantly associated with SNP rs5770917 (P_allele = 3.6×10⁻³; OR = 1.56; 95% c.i.: 1.12–2.15) and HLA-DRB1*1501-DQB1*0602 haplotype (P _positivity = 9.2×10⁻¹¹; OR = 3.97; 95% c.i.: 2.55–6.19). No interaction between the two markers (SNP rs5770917 and HLA-DRB1*1501-DQB1*0602 haplotype) was observed in EHS.

Conclusion

CPT1B, CHKB and HLA are candidates for susceptibility to CNS hypersomnias (EHS), as well as narcolepsy with cataplexy.     

SORL1 Is Genetically Associated with Late-Onset Alzheimer’s Disease in Japanese,Koreans and Caucasians

To discover susceptibility genes of late-onset Alzheimer’s disease (LOAD), we conducted a 3-stage genome-wide association study (GWAS) using three populations: Japanese from the Japanese Genetic Consortium for Alzheimer Disease (JGSCAD), Koreans, and Caucasians from the Alzheimer Disease Genetic Consortium (ADGC). In Stage 1, we evaluated data for 5,877,918 genotyped and imputed SNPs in Japanese cases (n = 1,008) and controls (n = 1,016). Genome-wide significance was observed with 12 SNPs in the APOE region. Seven SNPs from other distinct regions with p-values <2×10⁻⁵ were genotyped in a second Japanese sample (885 cases, 985 controls), and evidence of association was confirmed for one SORL1 SNP (rs3781834, P = 7.33×10⁻⁷ in the combined sample). Subsequent analysis combining results for several SORL1 SNPs in the Japanese, Korean (339 cases, 1,129 controls) and Caucasians (11,840 AD cases, 10,931 controls) revealed genome wide significance with rs11218343 (P = 1.77×10⁻⁹) and rs3781834 (P = 1.04×10⁻⁸). SNPs in previously established AD loci in Caucasians showed strong evidence of association in Japanese including rs3851179 near PICALM (P = 1.71×10⁻⁵) and rs744373 near BIN1 (P = 1.39×10⁻⁴). The associated allele for each of these SNPs was the same as in Caucasians. These data demonstrate for the first time genome-wide significance of LOAD with SORL1 and confirm the role of other known loci for LOAD in Japanese. Our study highlights the importance of examining associations in multiple ethnic populations.     

Association of Aminoacyl-tRNA Synthetases Gene Polymorphisms with the Risk of Congenital Heart Disease in the Chinese Han Population

Aminoacyl-tRNA synthetases (ARSs) are in charge of cellular protein synthesis and have additional domains that function in a versatile manner beyond translation. Eight core ARSs (EPRS, MRS, QRS, RRS, IRS, LRS, KRS, DRS) combined with three nonenzymatic components form a complex known as multisynthetase complex (MSC).We hypothesize that the single-nucleotide polymorphisms (SNPs) of the eight core ARS coding genes might influence the susceptibility of sporadic congenital heart disease (CHD). Thus, we conducted a case-control study of 984 CHD cases and 2953 non-CHD controls in the Chinese Han population to evaluate the associations of 16 potentially functional SNPs within the eight ARS coding genes with the risk of CHD. We observed significant associations with the risk of CHD for rs1061248 [G/A; odds ratio (OR) = 0.90, 95% confidence interval (CI) = 0.81–0.99; P = 3.81×10⁻²], rs2230301 [A/C; OR = 0.73, 95%CI = 0.60–0.90, P = 3.81×10⁻²], rs1061160 [G/A; OR = 1.18, 95%CI = 1.06–1.31; P = 3.53×10⁻³] and rs5030754 [G/A; OR = 1.39, 95%CI = 1.11–1.75; P = 4.47×10⁻³] of EPRS gene. After multiple comparisons, rs1061248 conferred no predisposition to CHD. Additionally, a combined analysis showed a significant dosage-response effect of CHD risk among individuals carrying the different number of risk alleles (P _trend = 5.00×10⁻⁴). Compared with individuals with “0–2” risk allele, those carrying “3”, “4” or “5 or more” risk alleles had a 0.97-, 1.25- or 1.38-fold increased risk of CHD, respectively. These findings indicate that genetic variants of the EPRS gene may influence the individual susceptibility to CHD in the Chinese Han population.     

Interactions between Urinary 4-tert-Octylphenol Levels and Metabolism Enzyme Gene Variants on Idiopathic Male Infertility

Octylphenol (OP) and Trichlorophenol (TCP) act as endocrine disruptors and have effects on male reproductive function. We studied the interactions between 4-tert-Octylphenol (4-t-OP), 4-n- Octylphenol (4-n-OP), 2,3,4-Trichlorophenol (2,3,4-TCP), 2,4,5-Trichlorophenol (2,4,5-TCP) urinary exposure levels and polymorphisms in selected xenobiotic metabolism enzyme genes among 589 idiopathic male infertile patients and 396 controls in a Han-Chinese population. Ultra high performance liquid chromatography-tandem mass spectrometry (UPLC-MS/MS) was used to measure alkylphenols and chlorophenols in urine. Polymorphisms were genotyped using the SNPstream platform and the Taqman method. Among four phenols that were detected, we found that only exposure to 4-t-OP increased the risk of male infertility (P _trend = 1.70×10⁻⁷). The strongest interaction was between 4-t-OP and rs4918758 in CYP2C9 (P _inter = 6.05×10⁻⁷). It presented a significant monotonic increase in risk estimates for male infertility with increasing 4-t-OP exposure levels among men with TC/CC genotype (low level compared with non-exposed, odds ratio (OR) = 2.26, 95% confidence intervals (CI) = 1.06, 4.83; high level compared with non-exposed, OR = 9.22, 95% CI = 2.78, 30.59), but no associations observed among men with TT genotype. We also found interactions between 4-t-OP and rs4986894 in CYP2C19, and between rs1048943 in CYP1A1, on male infertile risk (P _inter = 8.09×10⁻⁷, P _inter = 3.73×10⁻⁴, respectively).We observed notable interactions between 4-t-OP exposure and metabolism enzyme gene polymorphisms on idiopathic infertility in Han-Chinese men.     

Association of Genetic Variants Influencing Lipid Levels with Coronary Artery Disease in Japanese Individuals

Background/Objective

In Japanese populations, we performed a replication study of genetic loci previously identified in European-descent populations as being associated with lipid levels and risk of coronary artery disease (CAD).

Methods

We genotyped 48 single nucleotide polymorphisms (SNPs) from 22 candidate loci that had previously been identified by genome-wide association (GWA) meta-analyses for low-density lipoprotein cholesterol (LDL-C), high-density lipoprotein cholesterol (HDL-C), and/or triglycerides in Europeans. We selected 22 loci with 2 parallel tracks from 95 reported loci: 16 significant loci (p<1×10⁻³⁰ in Europeans) and 6 other loci including those with suggestive evidence of lipid associations in 1292 GWA-scanned Japanese samples. Genotyping was done in 4990 general population samples, and 1347 CAD cases and 1337 controls. For 9 SNPs, we further examined CAD associations in an additional panel of 3052 CAD cases and 6335 controls.

Principal Findings

Significant lipid associations (one-tailed p<0.05) were replicated for 18 of 22 loci in Japanese samples, with significant inter-ethnic heterogeneity at 4 loci–APOB, APOE-C1, CETP, and APOA5–and allelic heterogeneity. The strongest association was detected at APOE rs7412 for LDL-C (p = 1.3×10⁻⁴¹), CETP rs3764261 for HDL-C (p = 5.2×10⁻²⁴), and APOA5 rs662799 for triglycerides (p = 5.8×10⁻⁵⁴). CAD association was replicated and/or verified for 4 loci: SORT1 rs611917 (p = 1.7×10⁻⁸), APOA5 rs662799 (p = 0.0014), LDLR rs1433099 (p = 2.1×10⁻⁷), and APOE rs7412 (p = 6.1×10⁻¹³).

Conclusions

Our results confirm that most of the tested lipid loci are associated with lipid traits in the Japanese, further indicating that in genetic susceptibility to lipid levels and CAD, the related metabolic pathways are largely common across the populations, while causal variants at individual loci can be population-specific.     

Common variants in CRP and LEPR influence high sensitivity C-reactive protein levels in North Indians

Background

High sensitivity C-reactive protein (hsCRP) levels are shown to be influenced by genetic variants in Europeans; however, little is explored in Indian population.

Methods

Herein, we comprehensively evaluated association of all previously reported genetic determinants of hsCRP levels, including 18 cis (proximal to CRP gene) and 73 trans-acting (distal to CRP gene) variants in 4,200 North Indians of Indo-European ethnicity. First, we evaluated association of 91 variants from 12 candidate loci with hsCRP levels in 2,115 North Indians (1,042 non-diabetic subjects and 1,073 patients with type 2 diabetes). Then, cis and trans-acting variants contributing maximally to hsCRP level variation were further replicated in an independent 2,085 North Indians (1,047 patients with type 2 diabetes and 1,038 non-diabetic subjects).

Results

We found association of 12 variants from CRP, LEPR, IL1A, IL6, and IL6R with hsCRP levels in non-diabetic subjects. However, only rs3093059-CRP [β = 0.33, P = 9.6×10⁻⁵] and the haplotype harboring rs3093059 risk allele [β = 0.32 µg/mL, P = 1.4×10⁻⁴/P_perm = 9.0×10⁻⁴] retained significance after correcting for multiple testing. The cis-acting variant rs3093059-CRP had maximum contribution to the variance in hsCRP levels (1.14%). Among, trans-acting variants, rs1892534-LEPR was observed to contribute maximally to hsCRP level variance (0.59%). Associations of rs3093059-CRP and rs1892534-LEPR were confirmed by replication and attained higher significance after meta-analysis [β_meta = 0.26/0.22; P_meta = 4.3×10⁻⁷/7.4×10⁻³ and β_meta = −0.15/−0.12; P_meta = 2.0×10⁻⁶/1.6×10⁻⁶ for rs3093059 and rs1892534, respectively in non-diabetic subjects and all subjects taken together].

Conclusion

In conclusion, we identified rs3093059 in CRP and rs1892534 in LEPR as major cis and trans-acting contributor respectively, to the variance in hsCRP levels in North Indian population.