Replication of TCF4 through association and linkage studies in late-onset Fuchs endothelial corneal dystrophy

Fuchs endothelial corneal dystrophy (FECD) is a common, late-onset disorder of the corneal endothelium. Although progress has been made in understanding the genetic basis of FECD by studying large families in which the phenotype is transmitted in an autosomal dominant fashion, a recently reported genome-wide association study identified common alleles at a locus on chromosome 18 near TCF4 which confer susceptibility to FECD. Here, we report the findings of our independent validation study for TCF4 using the largest FECD dataset to date (450 FECD cases and 340 normal controls). Logistic regression with sex as a covariate was performed for three genetic models: dominant (DOM), additive (ADD), and recessive (REC). We found significant association with rs613872, the target marker reported by Baratz et al.(2010), for all three genetic models (DOM: P = 9.33×10⁻³⁵; ADD: P = 7.48×10⁻³⁰; REC: P = 5.27×10⁻⁶). To strengthen the association study, we also conducted a genome-wide linkage scan on 64 multiplex families, composed primarily of affected sibling pairs (ASPs), using both parametric and non-parametric two-point and multipoint analyses. The most significant linkage region localizes to chromosome 18 from 69.94cM to 85.29cM, with a peak multipoint HLOD = 2.5 at rs1145315 (75.58cM) under the DOM model, mapping 1.5 Mb proximal to rs613872. In summary, our study presents evidence to support the role of the intronic TCF4 single nucleotide polymorphism rs613872 in late-onset FECD through both association and linkage studies.     

RNA Toxicity and Missplicing in the Common Eye Disease Fuchs Endothelial Corneal Dystrophy

Fuchs endothelial corneal dystrophy (FECD) is an inherited degenerative disease that affects the internal endothelial cell monolayer of the cornea and can result in corneal edema and vision loss in severe cases. FECD affects ∼5% of middle-aged Caucasians in the United States and accounts for >14,000 corneal transplantations annually. Among the several genes and loci associated with FECD, the strongest association is with an intronic (CTG·CAG)_n trinucleotide repeat expansion in the TCF4 gene, which is found in the majority of affected patients. Corneal endothelial cells from FECD patients harbor a poly(CUG)_n RNA that can be visualized as RNA foci containing this condensed RNA and associated proteins. Similar to myotonic dystrophy type 1, the poly(CUG)_n RNA co-localizes with and sequesters the mRNA-splicing factor MBNL1, leading to missplicing of essential MBNL1-regulated mRNAs. Such foci and missplicing are not observed in similar cells from FECD patients who lack the repeat expansion. RNA-Seq splicing data from the corneal endothelia of FECD patients and controls reveal hundreds of differential alternative splicing events. These include events previously characterized in the context of myotonic dystrophy type 1 and epithelial-to-mesenchymal transition, as well as splicing changes in genes related to proposed mechanisms of FECD pathogenesis. We report the first instance of RNA toxicity and missplicing in a common non-neurological/neuromuscular disease associated with a repeat expansion. The FECD patient population with this (CTG·CAG)_n trinucleotide repeat expansion exceeds that of the combined number of patients in all other microsatellite expansion disorders.     

Reliability of Corneal Dynamic Scheimpflug Analyser Measurements in Virgin and Post-PRK Eyes

Purpose

To determine the measurement reliability of CorVis ST, a dynamic Scheimpflug analyser, in virgin and post-photorefractive keratectomy (PRK) eyes and compare the results between these two groups.

Methods

Forty virgin eyes and 42 post-PRK eyes underwent CorVis ST measurements performed by two technicians. Repeatability was evaluated by comparing three consecutive measurements by technician A. Reproducibility was determined by comparing the first measurement by technician A with one performed by technician B. Intraobserver and interobserver intraclass correlation coefficients (ICCs) were calculated. Univariate analysis of covariance (ANCOVA) was used to compare measured parameters between virgin and post-PRK eyes.

Results

The intraocular pressure (IOP), central corneal thickness (CCT) and 1_st applanation time demonstrated good intraobserver repeatability and interobserver reproducibility (ICC≧0.90) in virgin and post-PRK eyes. The deformation amplitude showed a good or close to good repeatability and reproducibility in both groups (ICC≧0.88). The CCT correlated positively with 1_st applanation time (r = 0.437 and 0.483, respectively, p<0.05) and negatively with deformation amplitude (r = −0.384 and −0.375, respectively, p<0.05) in both groups. Compared to post-PRK eyes, virgin eyes showed longer 1_st applanation time (7.29±0.21 vs. 6.96±0.17 ms, p<0.05) and lower deformation amplitude (1.06±0.07 vs. 1.17±0.08 mm, p<0.05).

Conclusions

CorVis ST demonstrated reliable measurements for CCT, IOP, and 1_st applanation time, as well as relatively reliable measurement for deformation amplitude in both virgin and post-PRK eyes. There were differences in 1_st applanation time and deformation amplitude between virgin and post-PRK eyes, which may reflect corneal biomechanical changes occurring after the surgery in the latter.     

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.     

Ocular Hypertension: General Characteristics and Estimated Cerebrospinal Fluid Pressure. The Beijing Eye Study 2011

Purpose

To examine characteristics of ocular hypertensive subjects and potential associations with estimated cerebrospinal fluid pressure (estCSFP).

Methods

The population-based Beijing Eye Study 2011 included 3468 individuals with a mean age of 64.6±9.8 years. Ocular hypertension was defined as intraocular pressure (IOP) >21 mmHg, normal optic nerve head appearance and normal retinal nerve fiber layer thickness. IOP was corrected for its dependence on central corneal thickness (CCT) and corneal curvature radius. Estimated CSFP was calculated as CSFP [mmHg] = 0.44×Body Mass Index [kg/m²]+0.16×Diastolic Blood Pressure [mmHg]−0.18×Age [Years]−1.91. Estimated trans-lamina cribrosa pressure difference (estTLCPD) was IOP–estCSFP.

Results

EstCSFP (10.5±3.6 mmHg versus 9.0±3.7 mmHg; P = 0.003) and estTLCPD (12.0±4.4 mmHg versus 5.4±3.8 mmHg; P<0.001) were higher in the ocular hypertensive group than in the normotensive group. In binary regression analysis, ocular hypertension was associated with increased estCSFP (P = 0.03; odds ratio (OR): 1.08; 95% confidence interval (CI): 1.01, 1.17) after adjusting for prevalence of arterial hypertension (P = 0.07; OR: 1.79; 95%CI: 0.96, 3.34), retinal nerve fiber layer thickness (P = 0.03; OR: 0.97; 95%CI: 0.95, 0.997) and blood glucose concentration (P = 0.006; OR: 1.17; 95%CI: 1.04, 1.30).

Conclusions

Ocular hypertensive subjects (with IOP correction for CCT and corneal curvature) as compared to ocular normotensive subjects had a significantly higher estCSFP in univariate analysis and in multivariate analysis. Despite of a higher estCSFP, estTLCPD was still markedly higher in ocular hypertensive eyes than in ocular normotensive eyes.     

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.     

Determinants of Posterior Corneal Biometric Measurements in a Multi-Ethnic Asian Population

Purpose

To describe the corneal and anterior segment determinants of posterior corneal arc length (PCAL) and posterior corneal curvature (PCC).

Methods

Cross-sectional, population-based study of 1069 subjects (1069 eyes) aged 40–80 years, from three major Asian ethnic groups. All underwent anterior segment optical coherence tomography imaging and analysis with Zhongshan Angle Assessment Program. Our main outcome measures were determinants of PCAL and PCC using adjusted, multivariate linear regression analysis, adjusted for confounders to obtain the estimated marginal means (EMM) with standard error (SE).

Results

The overall mean (± SD) of PCC was: 6.51±0.39 mm; and PCAL was: 12.52±0.59 mm. Malays had a relatively longer PCAL (EMM = 12.74 mm, SE = 0.04 mm) than Chinese (EMM = 12.48 mm, SE = 0.03 mm, P<0.001), and Indians (EMM = 12.42 mm, SE = 0.03 mm, P<0.001). Anterior segment parameters had weak-moderate correlations with PCAL, which included: anterior chamber depth (ACD) (r = 0.55, P<0.001), PCC (r = 0.27, P<0.001), anterior corneal curvature (ACC) (r = 0.14, P<0.001) and central corneal thickness (CCT) (r = −0.07, P = 0.023). In multivariate analysis, anterior segment parameters explained only 37.6% of the variance of PCAL, with ACD being the most important determinant (partial R²  = 0.300; P<0.001). The determinants of PCC included ACC, PCAL and CCT (explaining 72.1% variation of PCC), with ACC being the most important determinant (partial R²  = 0.683; P<0.001).

Conclusion

There was moderate correlation of PCAL with ACD, but anterior segment parameters accounted for only a small proportion of the variation in PCAL. The significant differences in PCAL and PCC amongst different Asian ethnic groups suggests that there is a need to consider this factor when planning for anterior segment surgeries such as endothelial keratoplasty.     

Novel loci for adiponectin levels and their influence on type 2 diabetes and metabolic traits: a multi-ethnic meta-analysis of 45,891 individuals

Circulating levels of adiponectin, a hormone produced predominantly by adipocytes, are highly heritable and are inversely associated with type 2 diabetes mellitus (T2D) and other metabolic traits. We conducted a meta-analysis of genome-wide association studies in 39,883 individuals of European ancestry to identify genes associated with metabolic disease. We identified 8 novel loci associated with adiponectin levels and confirmed 2 previously reported loci (P = 4.5×10⁻⁸–1.2×10⁻⁴³). Using a novel method to combine data across ethnicities (N = 4,232 African Americans, N = 1,776 Asians, and N = 29,347 Europeans), we identified two additional novel loci. Expression analyses of 436 human adipocyte samples revealed that mRNA levels of 18 genes at candidate regions were associated with adiponectin concentrations after accounting for multiple testing (p<3×10⁻⁴). We next developed a multi-SNP genotypic risk score to test the association of adiponectin decreasing risk alleles on metabolic traits and diseases using consortia-level meta-analytic data. This risk score was associated with increased risk of T2D (p = 4.3×10⁻³, n = 22,044), increased triglycerides (p = 2.6×10⁻¹⁴, n = 93,440), increased waist-to-hip ratio (p = 1.8×10⁻⁵, n = 77,167), increased glucose two hours post oral glucose tolerance testing (p = 4.4×10⁻³, n = 15,234), increased fasting insulin (p = 0.015, n = 48,238), but with lower in HDL-cholesterol concentrations (p = 4.5×10⁻¹³, n = 96,748) and decreased BMI (p = 1.4×10⁻⁴, n = 121,335). These findings identify novel genetic determinants of adiponectin levels, which, taken together, influence risk of T2D and markers of insulin resistance.     

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.     

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.     

Association of Common Variants in TCF4 and PTPRG with Fuchs' Corneal Dystrophy: A Systematic Review and Meta-Analysis

Topic

A meta-analysis of TCF4 and PTPRG gene variants in Fuchs'' corneal dystrophy (FCD).

Clinical relevance

To identify novel genetic markers in patients with FCD in different ethnic populations.

Methods

MEDLINE and EMBASE were searched for eligible genetic studies on TCF4 and PTPRG in FCD. Odds ratios (OR) and 95% confidence intervals (CI) of each single-nucleotide polymorphism (SNP) in allelic, dominant and recessive models were estimated using fixed-effect model if I²<50% in the test for heterogeneity, otherwise the random effects model was used.

Results

Thirty-three records were obtained, with 8 being suitable for meta-analysis, which included five SNPs in TCF4 and two in PTPRG. There were 1610 FCD patients and 1565 controls tested for TCF4 rs613872. This SNP was strongly associated with FCD in Caucasians (P = 5.0×10⁻¹⁰⁶), with the risk allele G conferring an OR of 3.95 (95% CI: 3.49–4.46). A further 4 TCF4 SNPs (rs17595731, rs2286812, rs618869 and rs9954153) were also significantly associated with FCD in Caucasians (P<10⁻⁸). However, we found no SNP associated with FCD in Chinese. In addition, there was no significant association between FCD and PTPRG.

Conclusion

TCF4 rs613872 is strongly associated with FCD in Caucasians but not in Chinese, which may suggest ethnic diversity in FCD SNP associations. SNPs in PTPRG were not associated with FCD in Caucasians or Chinese populations. Results of this meta-analysis indicate the need for larger-scale and multi-ethnic genetic studies on FCD to further explore the associated gene variants and their roles on the mechanism and genetic basis of FCD.     

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.     

Physical Activity Attenuates the Genetic Predisposition to Obesity in 20,000 Men and Women from EPIC-Norfolk Prospective Population Study

Background

We have previously shown that multiple genetic loci identified by genome-wide association studies (GWAS) increase the susceptibility to obesity in a cumulative manner. It is, however, not known whether and to what extent this genetic susceptibility may be attenuated by a physically active lifestyle. We aimed to assess the influence of a physically active lifestyle on the genetic predisposition to obesity in a large population-based study.

Methods and Findings

We genotyped 12 SNPs in obesity-susceptibility loci in a population-based sample of 20,430 individuals (aged 39–79 y) from the European Prospective Investigation of Cancer (EPIC)-Norfolk cohort with an average follow-up period of 3.6 y. A genetic predisposition score was calculated for each individual by adding the body mass index (BMI)-increasing alleles across the 12 SNPs. Physical activity was assessed using a self-administered questionnaire. Linear and logistic regression models were used to examine main effects of the genetic predisposition score and its interaction with physical activity on BMI/obesity risk and BMI change over time, assuming an additive effect for each additional BMI-increasing allele carried. Each additional BMI-increasing allele was associated with 0.154 (standard error [SE] 0.012) kg/m² (p = 6.73×10⁻³⁷) increase in BMI (equivalent to 445 g in body weight for a person 1.70 m tall). This association was significantly (p _interaction = 0.005) more pronounced in inactive people (0.205 [SE 0.024] kg/m² [p = 3.62×10⁻¹⁸; 592 g in weight]) than in active people (0.131 [SE 0.014] kg/m² [p = 7.97×10⁻²¹; 379 g in weight]). Similarly, each additional BMI-increasing allele increased the risk of obesity 1.116-fold (95% confidence interval [CI] 1.093–1.139, p = 3.37×10⁻²⁶) in the whole population, but significantly (p _interaction = 0.015) more in inactive individuals (odds ratio [OR] = 1.158 [95% CI 1.118–1.199; p = 1.93×10⁻¹⁶]) than in active individuals (OR = 1.095 (95% CI 1.068–1.123; p = 1.15×10⁻¹²]). Consistent with the cross-sectional observations, physical activity modified the association between the genetic predisposition score and change in BMI during follow-up (p _interaction = 0.028).

Conclusions

Our study shows that living a physically active lifestyle is associated with a 40% reduction in the genetic predisposition to common obesity, as estimated by the number of risk alleles carried for any of the 12 recently GWAS-identified loci. Please see later in the article for the Editors'' Summary     

Highly Significant Association between Two Common Single Nucleotide Polymorphisms in CORIN Gene and Preeclampsia in Caucasian Women

Preeclampsia is a frequent medical complication during pregnancy. Corin, a serine protease which activates pro-atrial natriuretic peptide, has recently been shown to be involved in the pathophysiology of preeclampsia. The aim of this study was to search for CORIN gene variations and their association to preeclampsia in Caucasian and African women. Our study population was composed of 571 pregnant women (295 with preeclampsia and 276 normotensive controls) matched for maternal and gestational age, and ethnic origin. The 22 exons of the CORIN gene were sequenced in a discovery sample (n = 260), where 31 single nucleotide polymorphisms were identified. In a replication sample (n = 311), 4 single nucleotide polymorphisms were tested. Two minor alleles (C for rs2271036 and G for rs2271037) were significantly associated to preeclampsia. Adjusted odds ratios [95% confidence interval] were 2.5 [1.2–3.8] (p = 0.007) and 2.3 [1.5–3.5] (p = 1.3×10⁻⁴), respectively. These associations were ethnic-specific, as only found in the Caucasian of subjects (odds ratio = 3.5 [1.8–6.6], p = 1.1×10⁻⁴; odds ratio = 3.1 [1.7–5.8], p = 2.1×10⁻⁴, for each single nucleotide polymorphism, respectively). The two single nucleotide polymorphisms are in almost perfect linkage disequilibrium (r² = 0.93). No specific association was found with severe preeclampsia, early-onset preeclampsia nor fetal growth retardation. In conclusion, this is the first report of a highly significant association between these two single nucleotide polymorphisms in CORIN gene and preeclampsia. Our findings further support the probability of a critical role of corin in preeclamspia pathophysiology at the uteroplacental interface.     

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 SNPs of CD40 Gene with Multiple Sclerosis in Russians

Multiple sclerosis (MS) is a serious, incurable neurological disease. In 2009, the ANZgene studies detected the suggestive association of located upstream of CD40 gene in chromosome 20q13 (p = 1.3×10⁻⁷). Identification of the causal variant(s) in the CD40 locus leads to a better understanding of the mechanism underlying the development of autoimmune pathologies. We determined the genotypes of rs6074022, rs1883832, rs1535045, and rs11086996 in patients with MS (n = 1684) and in the control group (n = 879). Two SNPs were significantly associated with MS: rs6074022 (additive model C allele OR = 1.27, 95% CI = [1.12–1.45], p = 3×10⁻⁴) and rs1883832 (additive model T allele OR = 1.20, 95% CI = [1.05–1.38], p = 7×10⁻³). In the meta-analysis of our results and the results of four previous studies, we obtain the association p-value of 2.34×10⁻¹², which confirmed the association between MS and rs6074022 at a genome-wide significant level. Next, we demonstrated that the model including rs6074022 only sufficiently described the association. From our analysis, we can speculate that the association between rs1883832 and MS was induced by LD, whereas rs6074022 was a marker in stronger LD with the functional variant or was the functional variant itself. Our results indicated that the functional variants were located in the upstream region of the gene CD40 and were in higher LD with rs6074022 than LD with rs1883832.     

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.