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.     

Association of Cytokines in Individuals Sensitive and Insensitive to Dust Mites in a Brazilian Population

Introduction

Allergic reaction to dust mites is a relatively common condition among children, triggering cutaneous and respiratory responses that have a great impact on the health of this population. Anaphylactic hypersensitivity is characterized by an exacerbated response involving the production of regulatory cytokines responsible for stimulating the production of IgE antibodies.

Objective

To investigate an association of variants in cytokine genes (IL1A ⁻⁸⁸⁹, IL1B ^{−511, +3962}, IL1R ¹⁹⁷⁰, IL1RA ¹¹¹⁰⁰, IL4RA ⁺¹⁹⁰², IL12 ⁻¹¹⁸⁸, IFNG ⁺⁸⁷⁴, TGFB1 ^{codon 10, codon 25}, TNFA ^{−308, −238}, IL2 ^{−330, +166}, IL4 ^{−1098, −590, −33}, IL6 ^{−174, nt565}, and IL10 ^{−1082, −819, −592}) between patients sensitive to dust mites and a control group.

Methods

A total of 254 patients were grouped as atopic and non-atopic according to sensitivity as evaluated by the Prick Test and to cytokine genotyping by the polymerase chain reaction-sequence specific primers (PCR-SSP) method using the Cytokine Genotyping Kit.

Results

A comparison between individuals allergic to Dermatophagoides farinae, Dermatophagoides pteronyssinus, and Blomia tropicalis and a non-atopic control group showed significant differences between allele and genotype frequencies in the regulatory regions of cytokine genes, with important evidence for IL4 ⁻⁵⁹⁰ in T/C (10.2% vs. 43.1%, odd ratio [OR] = 0.15, p = 5.2 10⁻⁸, pc = 0.0000011, and 95% confidence interval [95%CI] = 0.07–0.32) and T/T genotypes (42.9% vs. 13.8%, OR = 4.69, p = 2.5 10⁻⁶, pc = 0.000055, and 95%CI = 2.42–9.09). Other associations were observed in the pro-inflammatory cytokines IL1A ⁻⁸⁸⁹ (T/T, C, and T) and IL2 ⁻³³⁰ (G/T and T/T) and the anti-inflammatory cytokines IL4RA ⁺¹⁹⁰² (A and G), IL4 ⁻⁵⁹⁰ (T/C, T/T, C, and T), and IL10 ⁻⁵⁹² (A/A, C/A, A, and C).

Conclusion

Our results suggest a possible association between single nucleotide polymorphisms (SNPs) in cytokine genes and hypersensitivity to dust mites.     

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.     

Quantitative Assessment of the Association between rs2046210 at 6q25.1 and Breast Cancer Risk

Genome-wide association studies (GWAS) have identified several genetic susceptibility loci for breast cancer (BC). One of them, conducted among Chinese women, found an association of rs2046210 at 6q25.1 with the risk of BC recently. Since then, numerous association studies have been carried out to investigate the relationship between this polymorphism and BC risk in various populations. However, these have yielded contradictory results. We therefore performed a meta-analysis to clarify this inconsistency. Overall, a total of 235003 subjects based on 13 studies were included in our study. Significantly increased BC risk was detected in the pooled analysis [allele contrast: OR = 1.13, 95%CI = 1.10–1.17, P(Z) <10⁻⁵, P(Q) <10⁻⁴; dominant model: OR = 1.21, 95%CI = 1.14–1.27, P(Z) <10⁻⁵, P(Q) <10⁻⁴; recessive model: OR = 1.18, 95%CI = 1.12–1.24, P(Z) <10⁻⁵, P(Q) = 0.04]. In addition, our data revealed that rs2046210 conferred greater risk in estrogen receptor (ER)-negative tumors [OR = 1.27, 95%CI = 1.15–1.40, P(Z) <10⁻⁵, P(Q) <10⁻⁴] than in ER-positive ones [OR = 1.18, 95%CI = 1.09–1.28, P(Z) <10⁻⁴, P(Q) = 0.0003]. When stratified by ethnicity, significant associations were found in Caucasian and Asian populations, but not detected among Africans. There was evidence of heterogeneity (P<0.05), however, the heterogeneity largely disappeared after stratification by ethnicity. The present meta-analysis demonstrated that the rs2046210 polymorphism may be associated with increased BC susceptibility, but this association varies in different ethnicities.     

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.     

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.     

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 Genome-Wide Association Study of Total Serum and Mite-Specific IgEs in Asthma Patients

Immunoglobulin E (IgE) is one of the central players in asthma and allergic diseases. Although the serum IgE level, a useful endophenotype, is generally increased in patients with asthma, genetic factors influencing IgE regulation in asthma are still not fully understood. To identify the genetic variations associated with total serum and mite-specific IgEs in asthmatics, a genome-wide association study (GWAS) of 657,366 single nucleotide polymorphisms (SNPs) was performed in 877 Korean asthmatics. This study found that several new genes might be associated with total IgE in asthmatics, such as CRIM1 (rs848512, P = 1.18×10⁻⁶; rs711254, P = 6.73×10⁻⁶), ZNF71 (rs10404342, P = 7.60×10⁻⁶), TLN1 (rs4879926, P = 7.74×10⁻⁶), and SYNPO2 (rs1472066, P = 8.36×10⁻⁶; rs1038770, P = 8.66×10⁻⁶). Regarding the association of specific IgE to house dust mites, it was observed that intergenic SNPs nearby to OPRK1 and LOC730217 might be associated with Dermatophagoides pteronyssinus (D.p.) and Dermatophagoides farinae (D.f.) in asthmatics, respectively. In further pathway analysis, the phosphatidylinositol signaling system and adherens junction pathways were estimated to play a role in the regulation of total IgE levels in asthma. Although functional evaluations and replications of these results in other populations are needed, this GWAS of serum IgE in asthmatics could facilitate improved understanding of the role of the newly identified genetic variants in asthma and its related phenotypes.     

Genome-Wide Association Studies Identify Two Novel BMP15 Mutations Responsible for an Atypical Hyperprolificacy Phenotype in Sheep

Some sheep breeds are naturally prolific, and they are very informative for the studies of reproductive genetics and physiology. Major genes increasing litter size (LS) and ovulation rate (OR) were suspected in the French Grivette and the Polish Olkuska sheep populations, respectively. To identify genetic variants responsible for the highly prolific phenotype in these two breeds, genome-wide association studies (GWAS) followed by complementary genetic and functional analyses were performed. Highly prolific ewes (cases) and normal prolific ewes (controls) from each breed were genotyped using the Illumina OvineSNP50 Genotyping Beadchip. In both populations, an X chromosome region, close to the BMP15 gene, harbored clusters of markers with suggestive evidence of association at significance levels between 1E⁻⁰⁵ and 1E⁻⁰⁷. The BMP15 candidate gene was then sequenced, and two novel non-conservative mutations called FecX^Gr and FecX^O were identified in the Grivette and Olkuska breeds, respectively. The two mutations were associated with the highly prolific phenotype (p_FecX^Gr = 5.98E⁻⁰⁶ and p_FecX^O = 2.55E⁻⁰⁸). Homozygous ewes for the mutated allele showed a significantly increased prolificacy (FecX^Gr/FecX^Gr, LS = 2.50±0.65 versus FecX⁺/FecX^Gr, LS = 1.93±0.42, p<1E⁻⁰³ and FecX^O/FecX^O, OR = 3.28±0.85 versus FecX⁺/FecX^O, OR = 2.02±0.47, p<1E⁻⁰³). Both mutations are located in very well conserved motifs of the protein and altered the BMP15 signaling activity in vitro using a BMP-responsive luciferase test in COV434 granulosa cells. Thus, we have identified two novel mutations in the BMP15 gene associated with increased LS and OR. Notably, homozygous FecX^Gr/FecX^Gr Grivette and homozygous FecX^O/FecX^O Olkuska ewes are hyperprolific in striking contrast with the sterility exhibited by all other known homozygous BMP15 mutations. Our results bring new insights into the key role played by the BMP15 protein in ovarian function and could contribute to a better understanding of the pathogenesis of women′s fertility disorders.     

Polymorphisms in the TNFA and IL6 Genes Represent Risk Factors for Autoimmune Thyroid Disease

Background

Autoimmune thyroid disease (AITD) comprises diseases including Hashimoto''s thyroiditis and Graves'' disease, both characterized by reactivity to autoantigens causing, respectively, inflammatory destruction and autoimmune stimulation of the thyroid-stimulating hormone receptor. AITD is the most common thyroid disease and the leading form of autoimmune disease in women. Cytokines are key regulators of the immune and inflammatory responses; therefore, genetic variants at cytokine-encoding genes are potential risk factors for AITD.

Methods

Polymorphisms in the IL6-174 G/C (rs1800795), TNFA-308 G/A (rs1800629), IL1B-511 C/T (rs16944), and IFNGR1-56 T/C (rs2234711) genes were assessed in a case-control study comprising 420 Hashimoto''s thyroiditis patients, 111 Graves'' disease patients and 735 unrelated controls from Portugal. Genetic variants were discriminated by real-time PCR using TaqMan SNP genotyping assays.

Results

A significant association was found between the allele A in TNFA-308 G/A and Hashimoto''s thyroiditis, both in the dominant (OR = 1.82, CI = 1.37–2.43, p-value = 4.4×10⁻⁵) and log-additive (OR = 1.64, CI = 1.28–2.10, p-value = 8.2×10⁻⁵) models. The allele C in IL6-174 G/C is also associated with Hashimoto''s thyroiditis, however, only retained significance after multiple testing correction in the log-additive model (OR = 1.28, CI = 1.06–1.54, p-value = 8.9×10⁻³). The group with Graves'' disease also registered a higher frequency of the allele A in TNFA-308 G/A compared with controls both in the dominant (OR = 1.85, CI = 1.19–2.87, p-value = 7.0×10⁻³) and log-additive (OR = 1.69, CI = 1.17–2.44, p-value = 6.6×10⁻³) models. The risk for Hashimoto''s thyroiditis and Graves'' disease increases with the number of risk alleles (OR for two risk alleles is, respectively, 2.27 and 2.59).

Conclusions

This study reports significant associations of genetic variants in TNFA and IL6 with the risk for AITD, highlighting the relevance of polymorphisms in inflammation-related genes in the etiopathogenesis of AITD.     

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.     

CX3CR1 Is a Modifying Gene of Survival and Progression in Amyotrophic Lateral Sclerosis

The objective of this study was to investigate the association of functional variants of the human CX3CR1 gene (Fractalkine receptor) with the risk of Amyotrophic Lateral Sclerosis (ALS), the survival and the progression rate of the disease symptoms in a Spanish ALS cohort. 187 ALS patients (142 sporadic [sALS] and 45 familial) and 378 controls were recruited. We investigated CX3CR1 V249I (rs3732379) and T280M (rs3732378) genotypes and their haplotypes as predictors of survival, the progression rate of the symptoms (as measured by ALSFRS-R and FVC decline) and the risk of suffering ALS disease. The results indicated that sALS patients with CX3CR1 249^I/I or 249^V/I genotypes presented a shorter survival time (42.27±4.90) than patients with 249^V/V genotype (67.65±7.42; diff −25.49 months 95%CI [−42.79,−8.18]; p = 0.004; adj-p = 0.018). The survival time was shorter in sALS patients with spinal topography and CX3CR1 249^I alleles (diff = −29.78 months; 95%CI [−49.42,−10.14]; p = 0.003). The same effects were also observed in the spinal sALS patients with 249^I–280^M haplotype (diff = −27.02 months; 95%CI [−49.57, −4.48]; p = 0.019). In the sALS group, the CX3CR1 249^I variant was associated with a faster progression of the disease symptoms (OR = 2.58; 95IC% [1.32, 5.07]; p = 0.006; adj-p = 0.027). There was no evidence for association of these two CX3CR1 variants with ALS disease risk. The association evidenced herein is clinically relevant and indicates that CX3CR1 could be a disease-modifying gene in sALS. The progression rate of the disease''s symptoms and the survival time is affected in patients with one or two copies of the CX3CR1 249I allele. The CX3CR1 is the most potent ALS survival genetic factor reported to date. These results reinforce the role of the immune system in ALS pathogenesis.     

Isoniazid Inhibits the Heme-Based Reactivity of Mycobacterium tuberculosis Truncated Hemoglobin N

Isoniazid represents a first-line anti-tuberculosis medication in prevention and treatment. This prodrug is activated by a mycobacterial catalase-peroxidase enzyme called KatG in Mycobacterium tuberculosis), thereby inhibiting the synthesis of mycolic acid, required for the mycobacterial cell wall. Moreover, isoniazid activation by KatG produces some radical species (e.g., nitrogen monoxide), that display anti-mycobacterial activity. Remarkably, the ability of mycobacteria to persist in vivo in the presence of reactive nitrogen and oxygen species implies the presence in these bacteria of (pseudo-)enzymatic detoxification systems, including truncated hemoglobins (trHbs). Here, we report that isoniazid binds reversibly to ferric and ferrous M. tuberculosis trHb type N (or group I; Mt-trHbN(III) and Mt-trHbN(II), respectively) with a simple bimolecular process, which perturbs the heme-based spectroscopic properties. Values of thermodynamic and kinetic parameters for isoniazid binding to Mt-trHbN(III) and Mt-trHbN(II) are K = (1.1±0.1)×10⁻⁴ M, k _on = (5.3±0.6)×10³ M⁻¹ s⁻¹ and k _off = (4.6±0.5)×10⁻¹ s⁻¹; and D = (1.2±0.2)×10⁻³ M, d _on = (1.3±0.4)×10³ M⁻¹ s⁻¹, and d _off = 1.5±0.4 s⁻¹, respectively, at pH 7.0 and 20.0°C. Accordingly, isoniazid inhibits competitively azide binding to Mt-trHbN(III) and Mt-trHbN(III)-catalyzed peroxynitrite isomerization. Moreover, isoniazid inhibits Mt-trHbN(II) oxygenation and carbonylation. Although the structure of the Mt-trHbN-isoniazid complex is not available, here we show by docking simulation that isoniazid binding to the heme-Fe atom indeed may take place. These data suggest a direct role of isoniazid to impair fundamental functions of mycobacteria, e.g. scavenging of reactive nitrogen and oxygen species, and metabolism.     

Targeting 160 Candidate Genes for Blood Pressure Regulation with a Genome-Wide Genotyping Array

The outcome of Genome-Wide Association Studies (GWAS) has challenged the field of blood pressure (BP) genetics as previous candidate genes have not been among the top loci in these scans. We used Affymetrix 500K genotyping data of KORA S3 cohort (n = 1,644; Southern-Germany) to address (i) SNP coverage in 160 BP candidate genes; (ii) the evidence for associations with BP traits in genome-wide and replication data, and haplotype analysis. In total, 160 gene regions (genic region±10 kb) covered 2,411 SNPs across 11.4 Mb. Marker densities in genes varied from 0 (n = 11) to 0.6 SNPs/kb. On average 52.5% of the HAPMAP SNPs per gene were captured. No evidence for association with BP was obtained for 1,449 tested SNPs. Considerable associations (P<10⁻³) were detected for the genes, where >50% of HAPMAP SNPs were tagged. In general, genes with higher marker density (>0.2 SNPs/kb) revealed a better chance to reach close to significance associations. Although, none of the detected P-values remained significant after Bonferroni correction (P<0.05/2319, P<2.15×10⁻⁵), the strength of some detected associations was close to this level: rs10889553 (LEPR) and systolic BP (SBP) (P = 4.5×10⁻⁵) as well as rs10954174 (LEP) and diastolic BP (DBP) (P = 5.20×10⁻⁵). In total, 12 markers in 7 genes (ADRA2A, LEP, LEPR, PTGER3, SLC2A1, SLC4A2, SLC8A1) revealed considerable association (P<10⁻³) either with SBP, DBP, and/or hypertension (HYP). None of these were confirmed in replication samples (KORA S4, HYPEST, BRIGHT). However, supportive evidence for the association of rs10889553 (LEPR) and rs11195419 (ADRA2A) with BP was obtained in meta-analysis across samples stratified either by body mass index, smoking or alcohol consumption. Haplotype analysis highlighted LEPR and PTGER3. In conclusion, the lack of associations in BP candidate genes may be attributed to inadequate marker coverage on the genome-wide arrays, small phenotypic effects of the loci and/or complex interaction with life-style and metabolic parameters.     

Replication of 6 Obesity Genes in a Meta-Analysis of Genome-Wide Association Studies from Diverse Ancestries

Obesity is a major public health problem with a significant genetic component. Multiple DNA polymorphisms/genes have been shown to be strongly associated with obesity, typically in populations of European descent. The aim of this study was to verify the extent to which 6 confirmed obesity genes (FTO, CTNNBL1, ADRB2, LEPR, PPARG and UCP2 genes) could be replicated in 8 different samples (n = 11,161) and to explore whether the same genes contribute to obesity-susceptibility in populations of different ancestries (five Caucasian, one Chinese, one African-American and one Hispanic population). GWAS-based data sets with 1000 G imputed variants were tested for association with obesity phenotypes individually in each population, and subsequently combined in a meta-analysis. Multiple variants at the FTO locus showed significant associations with BMI, fat mass (FM) and percentage of body fat (PBF) in meta-analysis. The strongest association was detected at rs7185735 (P-value = 1.01×10⁻⁷ for BMI, 1.80×10⁻⁶ for FM, and 5.29×10⁻⁴ for PBF). Variants at the CTNNBL1, LEPR and PPARG loci demonstrated nominal association with obesity phenotypes (meta-analysis P-values ranging from 1.15×10⁻³ to 4.94×10⁻²). There was no evidence of association with variants at ADRB2 and UCP2 genes. When stratified by sex and ethnicity, FTO variants showed sex-specific and ethnic-specific effects on obesity traits. Thus, it is likely that FTO has an important role in the sex- and ethnic-specific risk of obesity. Our data confirmed the role of FTO, CTNNBL1, LEPR and PPARG in obesity predisposition. These findings enhanced our knowledge of genetic associations between these genes and obesity-related phenotypes, and provided further justification for pursuing functional studies of these genes in the pathophysiology of obesity. Sex and ethnic differences in genetic susceptibility across populations of diverse ancestries may contribute to a more targeted prevention and customized treatment of obesity.     

A Genome-Wide Association Study Identifies Potential Susceptibility Loci for Hirschsprung Disease

Hirschsprung disease (HSCR) is a congenital and heterogeneous disorder characterized by the absence of intramural nervous plexuses along variable lengths of the hindgut. Although RET is a well-established risk factor, a recent genome-wide association study (GWAS) of HSCR has identified NRG1 as an additional susceptibility locus. To discover additional risk loci, we performed a GWAS of 123 sporadic HSCR patients and 432 unaffected controls using a large-scale platform with coverage of over 1 million polymorphic markers. The result was that our study replicated the findings of RET-CSGALNACT2-RASGEF1A genomic region (_rawP = 5.69×10⁻¹⁹ before a Bonferroni correction; _corrP = 4.31×10⁻¹³ after a Bonferroni correction) and NRG1 as susceptibility loci. In addition, this study identified SLC6A20 (_adjP = 2.71×10⁻⁶), RORA (_adjP = 1.26×10⁻⁵), and ABCC9 (_adjP = 1.86×10⁻⁵) as new potential susceptibility loci under adjusting the already known loci on the RET-CSGALNACT2-RASGEF1A and NRG1 regions, although none of the SNPs in these genes passed the Bonferroni correction. In further subgroup analysis, the RET-CSGALNACT2-RASGEF1A genomic region was observed to have different significance levels among subgroups: short-segment (S-HSCR, _corrP = 1.71×10⁻⁵), long-segment (L-HSCR, _corrP = 6.66×10⁻⁴), and total colonic aganglionosis (TCA, _corrP>0.05). This differential pattern in the significance level suggests that other genomic loci or mechanisms may affect the length of aganglionosis in HSCR subgroups during enteric nervous system (ENS) development. Although functional evaluations are needed, our findings might facilitate improved understanding of the mechanisms of HSCR pathogenesis.     

Epigenome-Wide DNA Methylation in Hearing Ability: New Mechanisms for an Old Problem

Epigenetic regulation of gene expression has been shown to change over time and may be associated with environmental exposures in common complex traits. Age-related hearing impairment is a complex disorder, known to be heritable, with heritability estimates of 57–70%. Epigenetic regulation might explain the observed difference in age of onset and magnitude of hearing impairment with age. Epigenetic epidemiology studies using unrelated samples can be limited in their ability to detect small effects, and recent epigenetic findings in twins underscore the power of this well matched study design. We investigated the association between venous blood DNA methylation epigenome-wide and hearing ability. Pure-tone audiometry (PTA) and Illumina HumanMethylation array data were obtained from female twin volunteers enrolled in the TwinsUK register. Two study groups were explored: first, an epigenome-wide association scan (EWAS) was performed in a discovery sample (n = 115 subjects, age range: 47–83 years, Illumina 27 k array), then replication of the top ten associated probes from the discovery EWAS was attempted in a second unrelated sample (n = 203, age range: 41–86 years, Illumina 450 k array). Finally, a set of monozygotic (MZ) twin pairs (n = 21 pairs) within the discovery sample (Illumina 27 k array) was investigated in more detail in an MZ discordance analysis. Hearing ability was strongly associated with DNA methylation levels in the promoter regions of several genes, including TCF25 (cg01161216, p = 6.6×10⁻⁶), FGFR1 (cg15791248, p = 5.7×10⁻⁵) and POLE (cg18877514, p = 6.3×10⁻⁵). Replication of these results in a second sample confirmed the presence of differential methylation at TCF25 (p(replication) = 6×10⁻⁵) and POLE (p(replication) = 0.016). In the MZ discordance analysis, twins'' intrapair difference in hearing ability correlated with DNA methylation differences at ACP6 (cg01377755, r = −0.75, p = 1.2×10⁻⁴) and MEF2D (cg08156349, r = −0.75, p = 1.4×10⁻⁴). Examination of gene expression in skin, suggests an influence of differential methylation on expression, which may account for the variation in hearing ability with age.     

Differing Roles for TCF4 and COL8A2 in Central Corneal Thickness and Fuchs Endothelial Corneal Dystrophy

Fuchs endothelial corneal dystrophy (FECD) is the most common late-onset, vision-threatening corneal dystrophy in the United States, affecting about 4% of the population. Advanced FECD involves a thickening of the cornea from stromal edema and changes in Descemet membrane. To understand the relationship between FECD and central corneal thickness (CCT), we characterized common genetic variation in COL8A2 and TCF4, genes previously implicated in CCT and/or FECD. Other genes previously associated with FECD (PITX2, ZEB1, SLC4A11), and genes only known to affect CCT (COL5A1, FOXO1, AVGR8, ZNF469) were also interrogated. FECD probands, relatives and controls were recruited from 32 clinical sites; a total of 532 cases and 204 controls were genotyped and tested for association of FECD case/control status, a 7-step FECD severity scale and CCT, adjusting for age and sex. Association of FECD grade with TCF4 was highly significant (OR  = 6.01 at rs613872; p = 4.8×10⁻²⁵), and remained significant when adjusted for changes in CCT (OR  = 4.84; p = 2.2×10⁻¹⁶). Association of CCT with TCF4 was also significant (p = 6.1×10⁻⁷), but was abolished with adjustment for FECD grade (p = 0.92). After adjusting for FECD grade, markers in other genes examined were modestly associated (p ∼ 0.001) with FECD and/or CCT. Thus, common variants in TCF4 appear to influence FECD directly, and CCT secondarily via FECD. Additionally, changes in corneal thickness due to the effect of other loci may modify disease severity, age-at-onset, or other biomechanical characteristics.