The PNPLA3 rs738409 G-allele associates with reduced fasting serum triglyceride and serum cholesterol in Danes with impaired glucose regulation

Background and Aim

Non-alcoholic fatty liver disease (NAFLD) is a common condition, associated with hepatic insulin resistance and the metabolic syndrome including hyperglycaemia and dyslipidemia. We aimed at studying the potential impact of the NAFLD-associated PNPLA3 rs738409 G-allele on NAFLD-related metabolic traits in hyperglycaemic individuals.

Methods

The rs738409 variant was genotyped in the population-based Inter99 cohort examined by an oral glucose-tolerance test, and a combined study-sample consisting of 192 twins (96 twin pairs) and a sub-set of the Inter99 population (n = 63) examined by a hyperinsulinemic euglycemic clamp (n _total = 255). In Inter99, we analyzed associations of rs738409 with components of the WHO-defined metabolic syndrome (n = 5,847) and traits related to metabolic disease (n = 5,663). In the combined study sample we elucidated whether the rs738409 G-allele altered hepatic or peripheral insulin sensitivity. Study populations were divided into individuals with normal glucose-tolerance (NGT) and with impaired glucose regulation (IGR).

Results

The case-control study showed no associations with components of the metabolic syndrome or the metabolic syndrome. Among 1,357 IGR individuals, the rs738409 G-allele associated with decreased fasting serum triglyceride levels (per allele effect(β) = −9.9% [−14.4%;−4.0% (95% CI)], p = 5.1×10⁻⁵) and fasting total cholesterol (β = −0.2 mmol/l [−0.3;−0.01 mmol/l(95% CI)], p = 1.5×10⁻⁴). Meta-analyses showed no impact on hepatic or peripheral insulin resistance in carriers of the rs738409 G-allele.

Conclusion

Our findings suggest that the G-allele of PNPLA3 rs738409 associates with reduced fasting levels of cholesterol and triglyceride in individuals with IGR.     

Genetic variants on chromosome 1q41 influence ocular axial length and high myopia

As one of the leading causes of visual impairment and blindness, myopia poses a significant public health burden in Asia. The primary determinant of myopia is an elongated ocular axial length (AL). Here we report a meta-analysis of three genome-wide association studies on AL conducted in 1,860 Chinese adults, 929 Chinese children, and 2,155 Malay adults. We identified a genetic locus on chromosome 1q41 harboring the zinc-finger 11B pseudogene ZC3H11B showing genome-wide significant association with AL variation (rs4373767, β = −0.16 mm per minor allele, P_meta = 2.69×10⁻¹⁰). The minor C allele of rs4373767 was also observed to significantly associate with decreased susceptibility to high myopia (per-allele odds ratio (OR) = 0.75, 95% CI: 0.68–0.84, P_meta = 4.38×10⁻⁷) in 1,118 highly myopic cases and 5,433 controls. ZC3H11B and two neighboring genes SLC30A10 and LYPLAL1 were expressed in the human neural retina, retinal pigment epithelium, and sclera. In an experimental myopia mouse model, we observed significant alterations to gene and protein expression in the retina and sclera of the unilateral induced myopic eyes for the murine genes ZC3H11A, SLC30A10, and LYPLAL1. This supports the likely role of genetic variants at chromosome 1q41 in influencing AL variation and high myopia.     

Identification of novel susceptibility Loci for kawasaki disease in a Han chinese population by a genome-wide association study

Kawasaki disease (KD) is an acute systemic vasculitis syndrome that primarily affects infants and young children. Its etiology is unknown; however, epidemiological findings suggest that genetic predisposition underlies disease susceptibility. Taiwan has the third-highest incidence of KD in the world, after Japan and Korea. To investigate novel mechanisms that might predispose individuals to KD, we conducted a genome-wide association study (GWAS) in 250 KD patients and 446 controls in a Han Chinese population residing in Taiwan, and further validated our findings in an independent Han Chinese cohort of 208 cases and 366 controls. The most strongly associated single-nucleotide polymorphisms (SNPs) detected in the joint analysis corresponded to three novel loci. Among these KD-associated SNPs three were close to the COPB2 (coatomer protein complex beta-2 subunit) gene: rs1873668 (p = 9.52×10⁻⁵), rs4243399 (p = 9.93×10⁻⁵), and rs16849083 (p = 9.93×10⁻⁵). We also identified a SNP in the intronic region of the ERAP1 (endoplasmic reticulum amino peptidase 1) gene (rs149481, p_best = 4.61×10⁻⁵). Six SNPs (rs17113284, rs8005468, rs10129255, rs2007467, rs10150241, and rs12590667) clustered in an area containing immunoglobulin heavy chain variable regions genes, with p_best-values between 2.08×10⁻⁵ and 8.93×10⁻⁶, were also identified. This is the first KD GWAS performed in a Han Chinese population. The novel KD candidates we identified have been implicated in T cell receptor signaling, regulation of proinflammatory cytokines, as well as antibody-mediated immune responses. These findings may lead to a better understanding of the underlying molecular pathogenesis of KD.     

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.     

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.     

Genetic loci associated with plasma phospholipid n-3 fatty acids: a meta-analysis of genome-wide association studies from the CHARGE Consortium

Long-chain n-3 polyunsaturated fatty acids (PUFAs) can derive from diet or from α-linolenic acid (ALA) by elongation and desaturation. We investigated the association of common genetic variation with plasma phospholipid levels of the four major n-3 PUFAs by performing genome-wide association studies in five population-based cohorts comprising 8,866 subjects of European ancestry. Minor alleles of SNPs in FADS1 and FADS2 (desaturases) were associated with higher levels of ALA (p = 3×10⁻⁶⁴) and lower levels of eicosapentaenoic acid (EPA, p = 5×10⁻⁵⁸) and docosapentaenoic acid (DPA, p = 4×10⁻¹⁵⁴). Minor alleles of SNPs in ELOVL2 (elongase) were associated with higher EPA (p = 2×10⁻¹²) and DPA (p = 1×10⁻⁴³) and lower docosahexaenoic acid (DHA, p = 1×10⁻¹⁵). In addition to genes in the n-3 pathway, we identified a novel association of DPA with several SNPs in GCKR (glucokinase regulator, p = 1×10⁻⁸). We observed a weaker association between ALA and EPA among carriers of the minor allele of a representative SNP in FADS2 (rs1535), suggesting a lower rate of ALA-to-EPA conversion in these subjects. In samples of African, Chinese, and Hispanic ancestry, associations of n-3 PUFAs were similar with a representative SNP in FADS1 but less consistent with a representative SNP in ELOVL2. Our findings show that common variation in n-3 metabolic pathway genes and in GCKR influences plasma phospholipid levels of n-3 PUFAs in populations of European ancestry and, for FADS1, in other ancestries.     

Association of genetic variants of melatonin receptor 1B with gestational plasma glucose level and risk of glucose intolerance in pregnant Chinese women

Background

This study aimed to explore the association of MTNR1B genetic variants with gestational plasma glucose homeostasis in pregnant Chinese women.

Methods

A total of 1,985 pregnant Han Chinese women were recruited and evaluated for gestational glucose tolerance status with a two-step approach. The four MTNR1B variants rs10830963, rs1387153, rs1447352, and rs2166706 which had been reported to associate with glucose levels in general non-pregnant populations, were genotyped in these women. Using an additive model adjusted for age and body mass index (BMI), association of these variants with gestational fasting and postprandial plasma glucose (FPG and PPG) levels were analyzed by multiple linear regression; relative risk of developing gestational glucose intolerance was calculated by logistic regression. Hardy-Weinberg Equilibrium was tested by Chi-square and linkage disequilibrium (LD) between these variants was estimated by measures of D′ and r².

Results

In the pregnant Chinese women, the MTNR1B variant rs10830963, rs1387153, rs2166706 and rs1447352 were shown to be associated with the increased 1 hour PPG level (p = 8.04×10⁻¹⁰, 5.49×10⁻⁶, 1.89×10⁻⁵ and 0.02, respectively). The alleles were also shown to be associated with gestational glucose intolerance with odds ratios (OR) of 1.64 (p = 8.03×10⁻¹¹), 1.43 (p = 1.94×10⁻⁶), 1.38 (p = 1.63×10⁻⁵) and 1.24 (p = 0.007), respectively. MTNR1B rs1387153, rs2166706 were shown to be associated with gestational FPG levels (p = 0.04). Our data also suggested that, the LD pattern of these variants in the studied women conformed to that in the general populations: rs1387153 and rs2166706 were in high LD, they linked moderately with rs10830963, but might not linked with rs1447352;rs10830963 might not link with rs1447352, either. In addition, the MTNR1B variants were not found to be associated with any other traits tested.

Conclusions

The MTNR1B is likely to be involved in the regulation of glucose homeostasis during pregnancy.     

Genome-wide gene-environment study identifies glutamate receptor gene GRIN2A as a Parkinson's disease modifier gene via interaction with coffee

Our aim was to identify genes that influence the inverse association of coffee with the risk of developing Parkinson''s disease (PD). We used genome-wide genotype data and lifetime caffeinated-coffee-consumption data on 1,458 persons with PD and 931 without PD from the NeuroGenetics Research Consortium (NGRC), and we performed a genome-wide association and interaction study (GWAIS), testing each SNP''s main-effect plus its interaction with coffee, adjusting for sex, age, and two principal components. We then stratified subjects as heavy or light coffee-drinkers and performed genome-wide association study (GWAS) in each group. We replicated the most significant SNP. Finally, we imputed the NGRC dataset, increasing genomic coverage to examine the region of interest in detail. The primary analyses (GWAIS, GWAS, Replication) were performed using genotyped data. In GWAIS, the most significant signal came from rs4998386 and the neighboring SNPs in GRIN2A. GRIN2A encodes an NMDA-glutamate-receptor subunit and regulates excitatory neurotransmission in the brain. Achieving P_2df = 10⁻⁶, GRIN2A surpassed all known PD susceptibility genes in significance in the GWAIS. In stratified GWAS, the GRIN2A signal was present in heavy coffee-drinkers (OR = 0.43; P = 6×10⁻⁷) but not in light coffee-drinkers. The a priori Replication hypothesis that “Among heavy coffee-drinkers, rs4998386_T carriers have lower PD risk than rs4998386_CC carriers” was confirmed: OR_Replication = 0.59, P_Replication = 10⁻³; OR_Pooled = 0.51, P_Pooled = 7×10⁻⁸. Compared to light coffee-drinkers with rs4998386_CC genotype, heavy coffee-drinkers with rs4998386_CC genotype had 18% lower risk (P = 3×10⁻³), whereas heavy coffee-drinkers with rs4998386_TC genotype had 59% lower risk (P = 6×10⁻¹³). Imputation revealed a block of SNPs that achieved P_2df<5×10⁻⁸ in GWAIS, and OR = 0.41, P = 3×10⁻⁸ in heavy coffee-drinkers. This study is proof of concept that inclusion of environmental factors can help identify genes that are missed in GWAS. Both adenosine antagonists (caffeine-like) and glutamate antagonists (GRIN2A-related) are being tested in clinical trials for treatment of PD. GRIN2A may be a useful pharmacogenetic marker for subdividing individuals in clinical trials to determine which medications might work best for which patients.     

The D299G/T399I Toll-like receptor 4 variant associates with body and liver fat: results from the TULIP and METSIM Studies

Background

Toll-like-receptor 4 (TLR) is discussed to provide a molecular link between obesity, inflammation and insulin resistance. Genetic studies with replications in non-diabetic individuals in regard to their fat distribution or insulin resistance according to their carrier status of a common toll-like receptor 4 (TLR4) variant (TLR4^D299G/T399I) are still lacking.

Methodology/Principal Findings

We performed a cross-sectional analysis in individuals phenotyped for prediabetic traits as body fat composition (including magnetic resonance imaging), blood glucose levels and insulin resistance (oral glucose tolerance testing, euglycemic hyperinsulinemic clamp), according to TLR4 genotype determined by candidate SNP analyses (rs4986790). We analyzed N = 1482 non-diabetic individuals from the TÜF/TULIP cohort (South Germany, aged 39±13 y, BMI 28.5±7.9, mean±SD) and N = 5327 non-diabetic participants of the METSIM study (Finland, males aged 58±6 y, BMI 26.8±3.8) for replication purposes. German TLR4^D299G/T399I carriers had a significantly increased body fat (XG in rs4986790: +6.98%, p = 0.03, dominant model, adjusted for age, gender) and decreased insulin sensitivity (XG: −15.3%, Matsuda model, p = 0.04; XG: −20.6%, p = 0.016, clamp; both dominant models adjusted for age, gender, body fat). In addition, both liver fat (AG: +49.7%; p = 0.002) and visceral adipose tissue (AG: +8.2%; p = 0.047, both adjusted for age, gender, body fat) were significantly increased in rs4986790 minor allele carriers, and the effect on liver fat remained significant also after additional adjustment for visceral fat (p = 0.014). The analysis in METSIM confirmed increased body fat content in association with the rare G allele in rs4986790 (AG: +1.26%, GG: +11.0%; p = 0.010, additive model, adjusted for age) and showed a non-significant trend towards decreased insulin sensitivity (AG: −0.99%, GG: −10.62%).

Conclusions/Significance

TLR4^D299G/T399I associates with increased total body fat, visceral fat, liver fat and decreased insulin sensitivity in non-diabetic Caucasians and may contribute to diabetes risk. This finding supports the role of TLR4 as a molecular link between obesity and insulin resistance.     

Variation in MSRA modifies risk of neonatal intestinal obstruction in cystic fibrosis

Meconium ileus (MI), a life-threatening intestinal obstruction due to meconium with abnormal protein content, occurs in approximately 15 percent of neonates with cystic fibrosis (CF). Analysis of twins with CF demonstrates that MI is a highly heritable trait, indicating that genetic modifiers are largely responsible for this complication. Here, we performed regional family-based association analysis of a locus that had previously been linked to MI and found that SNP haplotypes 5′ to and within the MSRA gene were associated with MI (P = 1.99×10⁻⁵ to 1.08×10⁻⁶; Bonferroni P = 0.057 to 3.1×10⁻³). The haplotype with the lowest P value showed association with MI in an independent sample of 1,335 unrelated CF patients (OR = 0.72, 95% CI [0.53–0.98], P = 0.04). Intestinal obstruction at the time of weaning was decreased in CF mice with Msra null alleles compared to those with wild-type Msra resulting in significant improvement in survival (P = 1.2×10⁻⁴). Similar levels of goblet cell hyperplasia were observed in the ilea of the Cftr ^−/− and Cftr ^−/− Msra ^−/− mice. Modulation of MSRA, an antioxidant shown to preserve the activity of enzymes, may influence proteolysis in the developing intestine of the CF fetus, thereby altering the incidence of obstruction in the newborn period. Identification of MSRA as a modifier of MI provides new insight into the biologic mechanism of neonatal intestinal obstruction caused by loss of CFTR function.     

Dietary Factors Impact on the Association between CTSS Variants and Obesity Related Traits

Background/Aims

Cathepsin S, a protein coded by the CTSS gene, is implicated in adipose tissue biology–this protein enhances adipose tissue development. Our hypothesis is that common variants in CTSS play a role in body weight regulation and in the development of obesity and that these effects are influenced by dietary factors–increased by high protein, glycemic index and energy diets.

Methods

Four tag SNPs (rs7511673, rs11576175, rs10888390 and rs1136774) were selected to capture all common variation in the CTSS region. Association between these four SNPs and several adiposity measurements (BMI, waist circumference, waist for given BMI and being a weight gainer–experiencing the greatest degree of unexplained annual weight gain during follow-up or not) given, where applicable, both as baseline values and gain during the study period (6–8 years) were tested in 11,091 European individuals (linear or logistic regression models). We also examined the interaction between the CTSS variants and dietary factors–energy density, protein content (in grams or in % of total energy intake) and glycemic index–on these four adiposity phenotypes.

Results

We found several associations between CTSS polymorphisms and anthropometric traits including baseline BMI (rs11576175 (SNP N°2), p = 0.02, β = −0.2446), and waist change over time (rs7511673 (SNP N°1), p = 0.01, β = −0.0433 and rs10888390 (SNP N°3), p = 0.04, β = −0.0342). In interaction with the percentage of proteins contained in the diet, rs11576175 (SNP N°2) was also associated with the risk of being a weight gainer (p_interaction = 0.01, OR = 1.0526)–the risk of being a weight gainer increased with the percentage of proteins contained in the diet.

Conclusion

CTSS variants seem to be nominally associated to obesity related traits and this association may be modified by dietary protein intake.     

The role of inflammatory pathway genetic variation on maternal metabolic phenotypes during pregnancy

Background

Since mediators of inflammation are associated with insulin resistance, and the risk of developing diabetes mellitus and gestational diabetes, we hypothesized that genetic variation in members of the inflammatory gene pathway impact glucose levels and related phenotypes in pregnancy. We evaluated this hypothesis by testing for association between genetic variants in 31 inflammatory pathway genes in the Hyperglycemia and Adverse Pregnancy Outcome (HAPO) cohort, a large multiethnic multicenter study designed to address the impact of glycemia less than overt diabetes on pregnancy outcome.

Results

Fasting, 1-hour, and 2-hour glucose, fasting and 1-hour C-peptide, and HbA1c levels were measured in blood samples obtained from HAPO participants during an oral glucose tolerance test at 24-32 weeks gestation. We tested for association between 458 SNPs mapping to 31 genes in the inflammatory pathway and metabolic phenotypes in 3836 European ancestry and 1713 Thai pregnant women. The strongest evidence for association was observed with TNF alpha and HbA1c (rs1052248; 0.04% increase per allele C; p-value = 4.4×10⁻⁵), RETN and fasting plasma glucose (rs1423096; 0.7 mg/dl decrease per allele A; p-value = 1.1×10⁻⁴), IL8 and 1 hr plasma glucose (rs2886920; 2.6 mg/dl decrease per allele T; p-value = 1.3×10⁻⁴), ADIPOR2 and fasting C-peptide (rs2041139; 0.55 ug/L decrease per allele A; p-value = 1.4×10⁻⁴), LEPR and 1-hour C-peptide (rs1171278; 0.62 ug/L decrease per allele T; p-value = 2.4×10⁻⁴), and IL6 and 1-hour plasma glucose (rs6954897; −2.29 mg/dl decrease per allele G, p-value = 4.3×10⁻⁴).

Conclusions

Based on the genes surveyed in this study the inflammatory pathway is unlikely to have a strong impact on maternal metabolic phenotypes in pregnancy although variation in individual members of the pathway (e.g. RETN, IL8, ADIPOR2, LEPR, IL6, and TNF alpha,) may contribute to metabolic phenotypes in pregnant women.     

A genome-wide association scan on the levels of markers of inflammation in Sardinians reveals associations that underpin its complex regulation

Identifying the genes that influence levels of pro-inflammatory molecules can help to elucidate the mechanisms underlying this process. We first conducted a two-stage genome-wide association scan (GWAS) for the key inflammatory biomarkers Interleukin-6 (IL-6), the general measure of inflammation erythrocyte sedimentation rate (ESR), monocyte chemotactic protein-1 (MCP-1), and high-sensitivity C-reactive protein (hsCRP) in a large cohort of individuals from the founder population of Sardinia. By analysing 731,213 autosomal or X chromosome SNPs and an additional ∼1.9 million imputed variants in 4,694 individuals, we identified several SNPs associated with the selected quantitative trait loci (QTLs) and replicated all the top signals in an independent sample of 1,392 individuals from the same population. Next, to increase power to detect and resolve associations, we further genotyped the whole cohort (6,145 individuals) for 293,875 variants included on the ImmunoChip and MetaboChip custom arrays. Overall, our combined approach led to the identification of 9 genome-wide significant novel independent signals—5 of which were identified only with the custom arrays—and provided confirmatory evidence for an additional 7. Novel signals include: for IL-6, in the ABO gene (rs657152, p = 2.13×10⁻²⁹); for ESR, at the HBB (rs4910472, p = 2.31×10⁻¹¹) and UCN119B/SPPL3 (rs11829037, p = 8.91×10⁻¹⁰) loci; for MCP-1, near its receptor CCR2 (rs17141006, p = 7.53×10⁻¹³) and in CADM3 (rs3026968, p = 7.63×10⁻¹³); for hsCRP, within the CRP gene (rs3093077, p = 5.73×10⁻²¹), near DARC (rs3845624, p = 1.43×10⁻¹⁰), UNC119B/SPPL3 (rs11829037, p = 1.50×10⁻¹⁴), and ICOSLG/AIRE (rs113459440, p = 1.54×10⁻⁰⁸) loci. Confirmatory evidence was found for IL-6 in the IL-6R gene (rs4129267); for ESR at CR1 (rs12567990) and TMEM57 (rs10903129); for MCP-1 at DARC (rs12075); and for hsCRP at CRP (rs1205), HNF1A (rs225918), and APOC-I (rs4420638). Our results improve the current knowledge of genetic variants underlying inflammation and provide novel clues for the understanding of the molecular mechanisms regulating this complex process.     

A genome-wide association study of nephrolithiasis in the Japanese population identifies novel susceptible Loci at 5q35.3, 7p14.3, and 13q14.1

Nephrolithiasis is a common nephrologic disorder with complex etiology. To identify the genetic factor(s) for nephrolithiasis, we conducted a three-stage genome-wide association study (GWAS) using a total of 5,892 nephrolithiasis cases and 17,809 controls of Japanese origin. Here we found three novel loci for nephrolithiasis: RGS14-SLC34A1-PFN3-F12 on 5q35.3 (rs11746443; P = 8.51×10⁻¹², odds ratio (OR) = 1.19), INMT-FAM188B-AQP1 on 7p14.3 (rs1000597; P = 2.16×10⁻¹⁴, OR = 1.22), and DGKH on 13q14.1 (rs4142110; P = 4.62×10⁻⁹, OR = 1.14). Subsequent analyses in 21,842 Japanese subjects revealed the association of SNP rs11746443 with the reduction of estimated glomerular filtration rate (eGFR) (P = 6.54×10⁻⁸), suggesting a crucial role for this variation in renal function. Our findings elucidated the significance of genetic variations for the pathogenesis of nephrolithiasis.     

The effect of chromosome 9p21 variants on cardiovascular disease may be modified by dietary intake: evidence from a case/control and a prospective study

Background

One of the most robust genetic associations for cardiovascular disease (CVD) is the Chromosome 9p21 region. However, the interaction of this locus with environmental factors has not been extensively explored. We investigated the association of 9p21 with myocardial infarction (MI) in individuals of different ethnicities, and tested for an interaction with environmental factors.

Methods and Findings

We genotyped four 9p21 SNPs in 8,114 individuals from the global INTERHEART study. All four variants were associated with MI, with odds ratios (ORs) of 1.18 to 1.20 (1.85×10⁻⁸≤p≤5.21×10⁻⁷). A significant interaction (p = 4.0×10⁻⁴) was observed between rs2383206 and a factor-analysis-derived “prudent” diet pattern score, for which a major component was raw vegetables. An effect of 9p21 on MI was observed in the group with a low prudent diet score (OR = 1.32, p = 6.82×10⁻⁷), but the effect was diminished in a step-wise fashion in the medium (OR = 1.17, p = 4.9×10⁻³) and high prudent diet scoring groups (OR = 1.02, p = 0.68) (p = 0.014 for difference). We also analyzed data from 19,129 individuals (including 1,014 incident cases of CVD) from the prospective FINRISK study, which used a closely related dietary variable. In this analysis, the 9p21 risk allele demonstrated a larger effect on CVD risk in the groups with diets low or average for fresh vegetables, fruits, and berries (hazard ratio [HR] = 1.22, p = 3.0×10⁻⁴, and HR = 1.35, p = 4.1×10⁻³, respectively) compared to the group with high consumption of these foods (HR = 0.96, p = 0.73) (p = 0.0011 for difference). The combination of the least prudent diet and two copies of the risk allele was associated with a 2-fold increase in risk for MI (OR = 1.98, p = 2.11×10⁻⁹) in the INTERHEART study and a 1.66-fold increase in risk for CVD in the FINRISK study (HR = 1.66, p = 0.0026).

Conclusions

The risk of MI and CVD conferred by Chromosome 9p21 SNPs appears to be modified by a prudent diet high in raw vegetables and fruits. Please see later in the article for the Editors'' Summary     

Genetic determinants of serum testosterone concentrations in men

Testosterone concentrations in men are associated with cardiovascular morbidity, osteoporosis, and mortality and are affected by age, smoking, and obesity. Because of serum testosterone''s high heritability, we performed a meta-analysis of genome-wide association data in 8,938 men from seven cohorts and followed up the genome-wide significant findings in one in silico (n = 871) and two de novo replication cohorts (n = 4,620) to identify genetic loci significantly associated with serum testosterone concentration in men. All these loci were also associated with low serum testosterone concentration defined as <300 ng/dl. Two single-nucleotide polymorphisms at the sex hormone-binding globulin (SHBG) locus (17p13-p12) were identified as independently associated with serum testosterone concentration (rs12150660, p = 1.2×10⁻⁴¹ and rs6258, p = 2.3×10⁻²²). Subjects with ≥3 risk alleles of these variants had 6.5-fold higher risk of having low serum testosterone than subjects with no risk allele. The rs5934505 polymorphism near FAM9B on the X chromosome was also associated with testosterone concentrations (p = 5.6×10⁻¹⁶). The rs6258 polymorphism in exon 4 of SHBG affected SHBG''s affinity for binding testosterone and the measured free testosterone fraction (p<0.01). Genetic variants in the SHBG locus and on the X chromosome are associated with a substantial variation in testosterone concentrations and increased risk of low testosterone. rs6258 is the first reported SHBG polymorphism, which affects testosterone binding to SHBG and the free testosterone fraction and could therefore influence the calculation of free testosterone using law-of-mass-action equation.     

A genome-wide scan of Ashkenazi Jewish Crohn's disease suggests novel susceptibility loci

Crohn''s disease (CD) is a complex disorder resulting from the interaction of intestinal microbiota with the host immune system in genetically susceptible individuals. The largest meta-analysis of genome-wide association to date identified 71 CD–susceptibility loci in individuals of European ancestry. An important epidemiological feature of CD is that it is 2–4 times more prevalent among individuals of Ashkenazi Jewish (AJ) descent compared to non-Jewish Europeans (NJ). To explore genetic variation associated with CD in AJs, we conducted a genome-wide association study (GWAS) by combining raw genotype data across 10 AJ cohorts consisting of 907 cases and 2,345 controls in the discovery stage, followed up by a replication study in 971 cases and 2,124 controls. We confirmed genome-wide significant associations of 9 known CD loci in AJs and replicated 3 additional loci with strong signal (p<5×10⁻⁶). Novel signals detected among AJs were mapped to chromosomes 5q21.1 (rs7705924, combined p = 2×10⁻⁸; combined odds ratio OR = 1.48), 2p15 (rs6545946, p = 7×10⁻⁹; OR = 1.16), 8q21.11 (rs12677663, p = 2×10⁻⁸; OR = 1.15), 10q26.3 (rs10734105, p = 3×10⁻⁸; OR = 1.27), and 11q12.1 (rs11229030, p = 8×10⁻⁹; OR = 1.15), implicating biologically plausible candidate genes, including RPL7, CPAMD8, PRG2, and PRG3. In all, the 16 replicated and newly discovered loci, in addition to the three coding NOD2 variants, accounted for 11.2% of the total genetic variance for CD risk in the AJ population. This study demonstrates the complementary value of genetic studies in the Ashkenazim.     

Associations of HLA-DP variants with hepatitis B virus infection in southern and northern Han Chinese populations: a multicenter case-control study

Background

Human leukocyte antigen DP (HLA-DP) locus has been reported to be associated with hepatitis B virus (HBV) infection in populations of Japan and Thailand. We aimed to examine whether the association can be replicated in Han Chinese populations.

Methodology/Principal Findings

Two HLA-DP variants rs2395309 and rs9277535 (the most strongly associated SNPs from each HLA-DP locus) were genotyped in three independent Han cohorts consisting of 2 805 cases and 1 796 controls. By using logistic regression analysis, these two SNPs in the HLA-DPA1 and HLA-DPB1 genes were significantly associated with HBV infection in Han Chinese populations (P = 0.021∼3.36×10⁻⁸ at rs2395309; P = 8.37×10⁻³∼2.68×10⁻¹⁰ at rs9277535). In addition, the genotype distributions of both sites (rs2395309 and rs9277535) were clearly different between southern and northern Chinese population (P = 8.95×10⁻⁵ at rs2395309; P = 1.64×10⁻⁹ at rs9277535). By using asymptomatic HBV carrier as control group, our study showed that there were no associations of two HLA-DP variants with HBV progression (P = 0.305∼0.822 and 0.163∼0.881 in southern Chinese population, respectively; P = 0.097∼0.697 and 0.198∼0.615 in northern Chinese population, respectively).

Conclusions

Our results confirmed that two SNPs (rs2395309 and rs9277535) in the HLA-DP loci were strongly associated with HBV infection in southern and northern Han Chinese populations, but not with HBV progression.     

Comparison of 6q25 Breast Cancer Hits from Asian and European Genome Wide Association Studies in the Breast Cancer Association Consortium (BCAC)

The 6q25.1 locus was first identified via a genome-wide association study (GWAS) in Chinese women and marked by single nucleotide polymorphism (SNP) rs2046210, approximately 180 Kb upstream of ESR1. There have been conflicting reports about the association of this locus with breast cancer in Europeans, and a GWAS in Europeans identified a different SNP, tagged here by rs12662670. We examined the associations of both SNPs in up to 61,689 cases and 58,822 controls from forty-four studies collaborating in the Breast Cancer Association Consortium, of which four studies were of Asian and 39 of European descent. Logistic regression was used to estimate odds ratios (OR) and 95% confidence intervals (CI). Case-only analyses were used to compare SNP effects in Estrogen Receptor positive (ER+) versus negative (ER−) tumours. Models including both SNPs were fitted to investigate whether the SNP effects were independent. Both SNPs are significantly associated with breast cancer risk in both ethnic groups. Per-allele ORs are higher in Asian than in European studies [rs2046210: OR (A/G) = 1.36 (95% CI 1.26–1.48), p = 7.6×10⁻¹⁴ in Asians and 1.09 (95% CI 1.07–1.11), p = 6.8×10⁻¹⁸ in Europeans. rs12662670: OR (G/T) = 1.29 (95% CI 1.19–1.41), p = 1.2×10⁻⁹ in Asians and 1.12 (95% CI 1.08–1.17), p = 3.8×10⁻⁹ in Europeans]. SNP rs2046210 is associated with a significantly greater risk of ER− than ER+ tumours in Europeans [OR (ER−) = 1.20 (95% CI 1.15–1.25), p = 1.8×10⁻¹⁷ versus OR (ER+) = 1.07 (95% CI 1.04–1.1), p = 1.3×10⁻⁷, p_{heterogeneity} = 5.1×10⁻⁶]. In these Asian studies, by contrast, there is no clear evidence of a differential association by tumour receptor status. Each SNP is associated with risk after adjustment for the other SNP. These results suggest the presence of two variants at 6q25.1 each independently associated with breast cancer risk in Asians and in Europeans. Of these two, the one tagged by rs2046210 is associated with a greater risk of ER− tumours.     

Genetic variants at chromosomes 2q35, 5p12, 6q25.1, 10q26.13, and 16q12.1 influence the risk of breast cancer in men

Male breast cancer accounts for approximately 1% of all breast cancer. To date, risk factors for male breast cancer are poorly defined, but certain risk factors and genetic features appear common to both male and female breast cancer. Genome-wide association studies (GWAS) have recently identified common single nucleotide polymorphisms (SNPs) that influence female breast cancer risk; 12 of these have been independently replicated. To examine if these variants contribute to male breast cancer risk, we genotyped 433 male breast cancer cases and 1,569 controls. Five SNPs showed a statistically significant association with male breast cancer: rs13387042 (2q35) (odds ratio (OR)  = 1.30, p = 7.98×10⁻⁴), rs10941679 (5p12) (OR = 1.26, p = 0.007), rs9383938 (6q25.1) (OR = 1.39, p = 0.004), rs2981579 (FGFR2) (OR = 1.18, p = 0.03), and rs3803662 (TOX3) (OR = 1.48, p = 4.04×10⁻⁶). Comparing the ORs for male breast cancer with the published ORs for female breast cancer, three SNPs—rs13387042 (2q35), rs3803662 (TOX3), and rs6504950 (COX11)—showed significant differences in ORs (p<0.05) between sexes. Breast cancer is a heterogeneous disease; the relative risks associated with loci identified to date show subtype and, based on these data, gender specificity. Additional studies of well-defined patient subgroups could provide further insight into the biological basis of breast cancer development.