Novel meta-analysis-derived type 2 diabetes risk loci do not determine prediabetic phenotypes

Background

Genome-wide association (GWA) studies identified a series of novel type 2 diabetes risk loci. Most of them were subsequently demonstrated to affect insulin secretion of pancreatic β-cells. Very recently, a meta-analysis of GWA data revealed nine additional risk loci with still undefined roles in the pathogenesis of type 2 diabetes. Using our thoroughly phenotyped cohort of subjects at an increased risk for type 2 diabetes, we assessed the association of the nine latest genetic variants with the predominant prediabetes traits, i.e., obesity, impaired insulin secretion, and insulin resistance.

Methodology/Principal Findings

One thousand five hundred and seventy-eight metabolically characterized non-diabetic German subjects were genotyped for the reported candidate single nucleotide polymorphisms (SNPs) JAZF1 rs864745, CDC123/CAMK1D rs12779790, TSPAN8/LGR5 rs7961581, THADA rs7578597, ADAMTS9 rs4607103, NOTCH2 rs10923931, DCD rs1153188, VEGFA rs9472138, and BCL11A rs10490072. Insulin sensitivity was derived from fasting glucose and insulin concentrations, oral glucose tolerance test (OGTT), and hyperinsulinemic-euglycemic clamp. Insulin secretion was estimated from OGTT data. After appropriate adjustment for confounding variables and Bonferroni correction for multiple comparisons (corrected α-level: p = 0.0014), none of the SNPs was reliably associated with adiposity, insulin sensitivity, or insulin secretion (all p≥0.0117, dominant inheritance model). The risk alleles of ADAMTS9 SNP rs4607103 and VEGFA SNP rs9472138 tended to associate with more than one measure of insulin sensitivity and insulin secretion, respectively, but did not reach formal statistical significance. The study was sufficiently powered (1-β = 0.8) to detect effect sizes of 0.19≤d≤0.25 (α = 0.0014) and 0.13≤d≤0.16 (α = 0.05).

Conclusions/Significance

In contrast to the first series of GWA-derived type 2 diabetes candidate SNPs, we could not detect reliable associations of the novel risk loci with prediabetic phenotypes. Possible weak effects of ADAMTS9 SNP rs4607103 and VEGFA SNP rs9472138 on insulin sensitivity and insulin secretion, respectively, await further confirmation by larger studies.     

Variance of the SGK1 gene is associated with insulin secretion in different European populations: results from the TUEF, EUGENE2, and METSIM studies

Hypothesis

Serum- and Glucocorticoid-inducible Kinase 1 (SGK1) is involved in the regulation of insulin secretion and may represent a candidate gene for the development of type 2 diabetes mellitus in humans.

Methods

Three independent European populations were analyzed for the association of SGK1 gene (SGK) variations and insulin secretion traits. The German TUEF project provided the screening population (N = 725), and four tagging SNPs (rs1763527, rs1743966, rs1057293, rs9402571) were investigated. EUGENE2 (N = 827) served as a replication cohort for the detected associations. Finally, the detected associations were validated in the METSIM study, providing 3798 non-diabetic and 659 diabetic (type 2) individuals.

Results

Carriers of the minor G allele in rs9402571 had significantly higher C-peptide levels in the 2 h OGTT (+10.8%, p = 0.04; dominant model) and higher AUC_C-Peptide/AUC_Glc ratios (+7.5%, p = 0.04) compared to homozygous wild type TT carriers in the screening population. As interaction analysis for BMI×rs9402571 was significant (p = 0.04) for the endpoint insulin secretion, we stratified the TUEF cohort for BMI, using a cut off point of BMI = 25. The effect on insulin secretion only remained significant in lean TUEF participants (BMI≤25). This finding was replicated in lean EUGENE2 rs9402571 minor allele carriers, who had a significantly higher AUC_Ins/AUC_Glc (TT: 226±7, XG: 246±9; p = 0.019). Accordingly, the METSIM trial revealed a lower prevalence of type 2 diabetes (OR: 0.85; 95%CI: 0.71–1.01; p = 0.065, dominant model) in rs9402571 minor allele carriers.

Conclusions

The rs9402571 SGK genotype associates with increased insulin secretion in lean non-diabetic TUEF/EUGENE2 participants and with lower diabetes prevalence in METSIM. Our study in three independent European populations supports the conclusion that SGK variability affects diabetes risk.     

T regulatory cells in cord blood--FOXP3 demethylation as reliable quantitative marker

Background

Regulatory T-cells (Tregs), characterized as CD4+CD25^hi T-cells expressing FOXP3, play a crucial role in controlling healthy immune development during early immune maturation. Recently, FOXP3 demethylation was suggested to be a novel marker for natural Tregs in adults. In cord blood, the role and function of Tregs and its demethylation is poorly understood. We assessed FOXP3 demethylation in cord blood in relation to previously used Treg markers such as CD4⁺CD25^hi, FOXP3 mRNA, protein expression, and suppressive Treg function.

Methodology

Cord blood mononuclear cells (CBMC) were isolated from 70 healthy neonates, stimulated for 3 days with the microbial stimulus lipid A (LpA), and allergen Dermatophagoides pteronyssinus (Derp1). Tregs (CD4+CD25^hi, intracellular, mRNA FOXP3 expression, isolated cells), DNA methylation of the FOXP3-locus and suppressive Treg function were assessed.

Principal Findings

Demethylation of FOXP3 in whole blood was specific for isolated CD4+CD25^hi Tregs. Demethylation of FOXP3 was positively correlated with unstimulated and LpA-stimulated FOXP3 mRNA-expression (p≤0.05), and CD4⁺CD25^hi T-cells (p≤0.03). Importantly, increased FOXP3 demethylation correlated with more efficient suppressive capacity of Tregs (r = 0.72, p = 0.005). Furthermore, FOXP3 demethylation was positively correlated with Th2 cytokines (IL-5, IL-13) following LpA-stimulation (p = 0.006/0.04), with Th2 and IL-17 following Derp1+LpA-stimulations (p≤0.009), but not Th1 cytokines (IFN-γ).

Conclusions

FOXP3 demethylation reliable quantifies Tregs in cord blood. FOXP3 demethylation corresponds well with the suppressive potential of Tregs. The resulting strict correlation with functionally suppressive Tregs and the relative ease of measurement render it into a valuable novel marker for large field studies assessing Tregs as qualitative marker indicative of functional activity.     

The role of osteopontin (OPN/SPP1) haplotypes in the susceptibility to Crohn's disease

Background

Osteopontin represents a multifunctional molecule playing a pivotal role in chronic inflammatory and autoimmune diseases. Its expression is increased in inflammatory bowel disease (IBD). The aim of our study was to analyze the association of osteopontin (OPN/SPP1) gene variants in a large cohort of IBD patients.

Methodology/Principal Findings

Genomic DNA from 2819 Caucasian individuals (n = 841 patients with Crohn''s disease (CD), n = 473 patients with ulcerative colitis (UC), and n = 1505 healthy unrelated controls) was analyzed for nine OPN SNPs (rs2728127, rs2853744, rs11730582, rs11739060, rs28357094, rs4754 = p.Asp80Asp, rs1126616 = p.Ala236Ala, rs1126772 and rs9138). Considering the important role of osteopontin in Th17-mediated diseases, we performed analysis for epistasis with IBD-associated IL23R variants and analyzed serum levels of the Th17 cytokine IL-22. For four OPN SNPs (rs4754, rs1126616, rs1126772 and rs9138), we observed significantly different distributions between male and female CD patients. rs4754 was protective in male CD patients (p = 0.0004, OR = 0.69). None of the other investigated OPN SNPs was associated with CD or UC susceptibility. However, several OPN haplotypes showed significant associations with CD susceptibility. The strongest association was found for a haplotype consisting of the 8 OPN SNPs rs2728127-rs2853744-rs11730582-rs11439060-rs28357094-rs112661-rs1126772-rs9138 (omnibus p-value = 2.07×10⁻⁸). Overall, the mean IL-22 secretion in the combined group of OPN minor allele carriers with CD was significantly lower than that of CD patients with OPN wildtype alleles (p = 3.66×10⁻⁵). There was evidence for weak epistasis between the OPN SNP rs28357094 with the IL23R SNP rs10489629 (p = 4.18×10⁻²) and between OPN SNP rs1126616 and IL23R SNP rs2201841 (p = 4.18×10⁻²) but none of these associations remained significant after Bonferroni correction.

Conclusions/Significance

Our study identified OPN haplotypes as modifiers of CD susceptibility, while the combined effects of certain OPN variants may modulate IL-22 secretion.     

Genes involved in the metabolism of poly-unsaturated fatty-acids (PUFA) and risk for Crohn's disease in children & young adults

Background and Objectives

Epidemiological evidence for the role of polyunsaturated fatty-acids (PUFA) in Crohn''s disease (CD) is unclear, although the key metabolite leucotriene B4 (LTB₄) is closely linked to the inflammatory process. We hypothesized that inherited variation in key PUFA metabolic enzymes may modify susceptibility for CD.

Methods and Principal Results

A case-control design was implemented at three pediatric gastroenterology clinics in Canada. Children ≤20 yrs diagnosed with CD and controls were recruited. 19 single nucleotide polymorphisms (SNPs) across the ALOX5 (4) CYP4F3 (5) and CYP4F2 (10) genes, were genotyped. Associations between SNPs/haplotypes and CD were examined. A total of 431 cases and 507 controls were studied. The mean (±SD) age of the cases was 12.4 (±3.3) years. Most cases were male (56.4%), had ileo-colonic disease (L3±L4, 52.7%) and inflammatory behavior (B1±p, 87%) at diagnosis. One genotyped CYP4F3 SNP (rs2683037) not in Hardy-Weinberg Equilibrium was excluded. No associations with the remaining 4 CYP4F3 SNPs with CD were evident. However haplotype analysis revealed associations with a two-marker haplotype (TG) (rs3794987 & rs1290617) (p = 0.02; permuted p = 0.08). CYP4F2 SNPs, rs3093158 (OR (recessive) = 0.56, 95% CI = 0.35–0.89; p = 0.01), rs2074902 (OR (trend) = 1.26, 95% CI = 1.00–1.60; p = 0.05), and rs2108622 (OR (recessive) = 1.6, 95% CI = 1.00–2.57; p = 0.05) were significantly associated whereas rs1272 (OR (recessive) = 0.58, 95% CI = 0.30–1.13; p = 0.10) showed suggestions for associations with CD. A haplotype comprising these 4 SNPs was significantly associated (p = 0.007, permuted p = 0.02) with CD. Associations with SNP rs3780901 in the ALOX5 gene were borderline non-significant (OR (dominant) = 1.29, 95% CI = 0.99–1.67; p = 0.056). A haplotype comprising the 4 ALOX5 SNPs (TCAA, p = 0.036) was associated with CD, but did not withstand corrections for multiple comparisons (permuted p = 0.14).

Conclusions

Inherited variation in enzymes involved in the synthesis/metabolism of LTB₄ may be associated with CD. These findings implicate PUFA metabolism as a important pathway in the CD pathogenesis.     

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.     

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.     

Caucasian and Asian specific rheumatoid arthritis risk loci reveal limited replication and apparent allelic heterogeneity in north Indians

Genome-wide association studies and meta-analysis indicate that several genes/loci are consistently associated with rheumatoid arthritis (RA) in European and Asian populations. To evaluate the transferability status of these findings to an ethnically diverse north Indian population, we performed a replication analysis. We investigated the association of 47 single-nucleotide polymorphisms (SNPs) at 43 of these genes/loci with RA in a north Indian cohort comprising 983 RA cases and 1007 age and gender matched controls. Genotyping was done using Infinium human 660w-quad. Association analysis by chi-square test implemented in plink was carried out in two steps. Firstly, association of the index or surrogate SNP (r2>0.8, calculated from reference GIH Hap-Map population) was tested. In the second step, evidence for allelic/locus heterogeneity at aforementioned genes/loci was assessed for by testing additional flanking SNPs in linkage equilibrium with index/surrogate marker.Of the 44 European specific index SNPs, neither index nor surrogate SNPs were present for nine SNPs in the genotyping array. Of the remaining 35, associations were replicated at seven genes namely PTPN22 (rs1217407, p = 3×10⁻³); IL2–21 (rs13119723, p = 0.008); HLA-DRB1 (rs660895, p = 2.56×10⁻⁵; rs6457617, p = 1.6×10⁻⁰⁹; rs13192471, p = 6.7×10⁻¹⁶); TNFA1P3 (rs9321637, p = 0.03); CCL21 (rs13293020, p = 0.01); IL2RA (rs2104286, p = 1.9×10⁻⁴) and ZEB1 (rs2793108, p = 0.006). Of the three Asian specific loci tested, rs2977227 in PADI4 showed modest association (p<0.02). Further, of the 140 SNPs (in LE with index/surrogate variant) tested, association was observed at 11 additional genes: PTPRC, AFF3, CD28, CTLA4, PXK, ANKRD55, TAGAP, CCR6, BLK, CD40 and IL2RB. This study indicates limited replication of European and Asian index SNPs and apparent allelic heterogeneity in RA etiology among north Indians warranting independent GWAS in this population. However, replicated associations of HLA-DRB1, PTPN22 (which confer ∼50% of the heritable risk to RA) and IL2RA suggest that cross-ethnicity fine mapping of such loci is apposite for identification of causal variants.     

Genome-wide association meta-analysis of cortical bone mineral density unravels allelic heterogeneity at the RANKL locus and potential pleiotropic effects on bone

Previous genome-wide association (GWA) studies have identified SNPs associated with areal bone mineral density (aBMD). However, this measure is influenced by several different skeletal parameters, such as periosteal expansion, cortical bone mineral density (BMD_C) cortical thickness, trabecular number, and trabecular thickness, which may be under distinct biological and genetic control. We have carried out a GWA and replication study of BMD_C, as measured by peripheral quantitative computed tomography (pQCT), a more homogenous and valid measure of actual volumetric bone density. After initial GWA meta-analysis of two cohorts (ALSPAC n = 999, aged ∼15 years and GOOD n = 935, aged ∼19 years), we attempted to replicate the BMD_C associations that had p<1×10⁻⁵ in an independent sample of ALSPAC children (n = 2803) and in a cohort of elderly men (MrOS Sweden, n = 1052). The rs1021188 SNP (near RANKL) was associated with BMD_C in all cohorts (overall p = 2×10⁻¹⁴, n = 5739). Each minor allele was associated with a decrease in BMD_C of ∼0.14SD. There was also evidence for an interaction between this variant and sex (p = 0.01), with a stronger effect in males than females (at age 15, males −6.77mg/cm³ per C allele, p = 2×10⁻⁶; females −2.79 mg/cm³ per C allele, p = 0.004). Furthermore, in a preliminary analysis, the rs1021188 minor C allele was associated with higher circulating levels of sRANKL (p<0.005). We show this variant to be independent from the previously aBMD associated SNP (rs9594738) and possibly from a third variant in the same RANKL region, which demonstrates important allelic heterogeneity at this locus. Associations with skeletal parameters reflecting bone dimensions were either not found or were much less pronounced. This finding implicates RANKL as a locus containing variation associated with volumetric bone density and provides further insight into the mechanism by which the RANK/RANKL/OPG pathway may be involved in skeletal development.     

Association of six single nucleotide polymorphisms with gestational diabetes mellitus in a Chinese population

Background

To investigate whether the candidate genes that confer susceptibility to type 2 diabetes mellitus are also correlated with gestational diabetes mellitus (GDM) in pregnant Chinese women.

Methodology/Principal Findings

In this study, 1764 unrelated pregnant women were recruited, of which 725 women had GDM and 1039 served as controls. Six single nucleotide polymorphisms (rs7754840 in CDKAL1, rs391300 in SRR, rs2383208 in CDKN2A/2B, rs4402960 in IGF2BP2, rs10830963 in MTNR1B, rs4607517 in GCK) were genotyped using TaqMan allelic discrimination assays. The genotype and allele distributions of each SNP between the GDM cases and controls and the combined effects of alleles for the risk of developing GDM were analyzed. We found that the rs4402960, rs2383208 and rs391300 were statistically associated with GDM (OR = 1.207, 95%CI = 1.029–1.417, p = 0.021; OR = 1.242, 95%CI = 1.077–1.432, p = 0.003; OR = 1.202, 95%CI = 1.020–1.416, P = 0.028, respectively). In addition, the effect was greater under a recessive model in rs391300 (OR = 1.820, 95%CI = 1.226–2.701, p = 0.003). Meanwhile, the joint effect of these three loci indicated an additive effect of multiple alleles on the risk of developing GDM with an OR of 1.196 per allele (p = 1.08×10⁻⁴). We also found that the risk alleles of rs2383208 (b = −0.085, p = 0.003), rs4402960 (b = −0.057, p = 0.046) and rs10830963 (b = −0.096, p = 0.001) were associated with HOMA-B, while rs7754840 was associated with decrease in insulin AUC during a 100 g OGTT given at the time of GDM diagnosis (b = −0.080, p = 0.007).

Conclusions/Significance

Several risk alleles of type 2 diabetes were associated with GDM in pregnant Chinese women. The effects of these SNPs on GDM might be through the impairment of beta cell function and these risk loci contributed additively to the disease.     

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.     

A genome-wide association study of upper aerodigestive tract cancers conducted within the INHANCE consortium

Genome-wide association studies (GWAS) have been successful in identifying common genetic variation involved in susceptibility to etiologically complex disease. We conducted a GWAS to identify common genetic variation involved in susceptibility to upper aero-digestive tract (UADT) cancers. Genome-wide genotyping was carried out using the Illumina HumanHap300 beadchips in 2,091 UADT cancer cases and 3,513 controls from two large European multi-centre UADT cancer studies, as well as 4,821 generic controls. The 19 top-ranked variants were investigated further in an additional 6,514 UADT cancer cases and 7,892 controls of European descent from an additional 13 UADT cancer studies participating in the INHANCE consortium. Five common variants presented evidence for significant association in the combined analysis (p≤5×10⁻⁷). Two novel variants were identified, a 4q21 variant (rs1494961, p = 1×10⁻⁸) located near DNA repair related genes HEL308 and FAM175A (or Abraxas) and a 12q24 variant (rs4767364, p = 2×10⁻⁸) located in an extended linkage disequilibrium region that contains multiple genes including the aldehyde dehydrogenase 2 (ALDH2) gene. Three remaining variants are located in the ADH gene cluster and were identified previously in a candidate gene study involving some of these samples. The association between these three variants and UADT cancers was independently replicated in 5,092 UADT cancer cases and 6,794 controls non-overlapping samples presented here (rs1573496-ADH7, p = 5×10⁻⁸; rs1229984-ADH1B, p = 7×10⁻⁹; and rs698-ADH1C, p = 0.02). These results implicate two variants at 4q21 and 12q24 and further highlight three ADH variants in UADT cancer susceptibility.     

Impact of IL28B-related single nucleotide polymorphisms on liver histopathology in chronic hepatitis C genotype 2 and 3

Background and Aims

Recently, several genome-wide association studies have revealed that single nucleotide polymorphisms (SNPs) in proximity to IL28B predict spontaneous clearance of HCV infection as well as outcome following peginterferon and ribavirin therapy among HCV genotype 1 infected patients. The present study aimed to evaluate the impact of IL28B SNP variability on liver histology in the context of a phase III treatment trial (NORDynamIC) for treatment-naïve patients with chronic HCV genotype 2 or 3 infection, where pretreatment liver biopsies were mandatory.

Methods

Three hundred and thirty-nine Caucasian patients had samples available for IL28B genotyping (rs12979860) of whom 314 had pretreatment liver biopsies that were evaluated using the Ishak protocol, allowing for detailed grading and staging of liver histopathology.

Results

IL28B CC_rs12979860 genotype in HCV genotype 3 infected patients was associated with higher ALT levels (p<0.0001), higher AST to platelet ratio index (APRI; p = 0.001), and higher baseline viral load (p<0.0001) as compared to patients with the CT or TT genotypes. Additionally the CC_rs12979860 genotype entailed more pronounced portal inflammation (p = 0.02) and steatosis (p = 0.03). None of these associations were noted among HCV genotype 2 infected patients.

Conclusion

This study shows that the CC_rs12979860 SNP is associated with more pronounced liver histopathology in patients chronically infected with HCV genotype 3, which may be secondary to higher viral load. The finding that IL28B variability did not impact on liver pathology or viral load among genotype 2 infected patients implies that IL28B may differentially regulate the course of genotype 2 and 3 infection.     

Evidence for more than one Parkinson's disease-associated variant within the HLA region

Parkinson''s disease (PD) was recently found to be associated with HLA in a genome-wide association study (GWAS). Follow-up GWAS''s replicated the PD-HLA association but their top hits differ. Do the different hits tag the same locus or is there more than one PD-associated variant within HLA? We show that the top GWAS hits are not correlated with each other (0.00≤r²≤0.15). Using our GWAS (2000 cases, 1986 controls) we conducted step-wise conditional analysis on 107 SNPs with P<10⁻³ for PD-association; 103 dropped-out, four remained significant. Each SNP, when conditioned on the other three, yielded P_SNP1 = 5×10⁻⁴, P_SNP2 = 5×10⁻⁴, P_SNP3 = 4×10⁻³ and P_SNP4 = 0.025. The four SNPs were not correlated (0.01≤r²≤0.20). Haplotype analysis (excluding rare SNP2) revealed increasing PD risk with increasing risk alleles from OR = 1.27, P = 5×10⁻³ for one risk allele to OR = 1.65, P = 4×10⁻⁸ for three. Using additional 843 cases and 856 controls we replicated the independent effects of SNP1 (P_{conditioned-on-SNP4} = 0.04) and SNP4 (P_{conditioned-on-SNP1} = 0.04); SNP2 and SNP3 could not be replicated. In pooled GWAS and replication, SNP1 had OR_{conditioned-on-SNP4} = 1.23, P_{conditioned-on-SNP4} = 6×10⁻⁷; SNP4 had OR_{conditioned-on-SNP1} = 1.18, P_{conditioned-on-SNP1} = 3×10⁻³; and the haplotype with both risk alleles had OR = 1.48, P = 2×10⁻¹². Genotypic OR increased with the number of risk alleles an individual possessed up to OR = 1.94, P = 2×10⁻¹¹ for individuals who were homozygous for the risk allele at both SNP1 and SNP4. SNP1 is a variant in HLA-DRA and is associated with HLA-DRA, DRB5 and DQA2 gene expression. SNP4 is correlated (r² = 0.95) with variants that are associated with HLA-DQA2 expression, and with the top HLA SNP from the IPDGC GWAS (r² = 0.60). Our findings suggest more than one PD-HLA association; either different alleles of the same gene, or separate loci.     

Genome-wide association studies of the PR interval in African Americans

The PR interval on the electrocardiogram reflects atrial and atrioventricular nodal conduction time. The PR interval is heritable, provides important information about arrhythmia risk, and has been suggested to differ among human races. Genome-wide association (GWA) studies have identified common genetic determinants of the PR interval in individuals of European and Asian ancestry, but there is a general paucity of GWA studies in individuals of African ancestry. We performed GWA studies in African American individuals from four cohorts (n = 6,247) to identify genetic variants associated with PR interval duration. Genotyping was performed using the Affymetrix 6.0 microarray. Imputation was performed for 2.8 million single nucleotide polymorphisms (SNPs) using combined YRI and CEU HapMap phase II panels. We observed a strong signal (rs3922844) within the gene encoding the cardiac sodium channel (SCN5A) with genome-wide significant association (p<2.5×10⁻⁸) in two of the four cohorts and in the meta-analysis. The signal explained 2% of PR interval variability in African Americans (beta  = 5.1 msec per minor allele, 95% CI  = 4.1–6.1, p = 3×10⁻²³). This SNP was also associated with PR interval (beta = 2.4 msec per minor allele, 95% CI = 1.8–3.0, p = 3×10⁻¹⁶) in individuals of European ancestry (n = 14,042), but with a smaller effect size (p for heterogeneity <0.001) and variability explained (0.5%). Further meta-analysis of the four cohorts identified genome-wide significant associations with SNPs in SCN10A (rs6798015), MEIS1 (rs10865355), and TBX5 (rs7312625) that were highly correlated with SNPs identified in European and Asian GWA studies. African ancestry was associated with increased PR duration (13.3 msec, p = 0.009) in one but not the other three cohorts. Our findings demonstrate the relevance of common variants to African Americans at four loci previously associated with PR interval in European and Asian samples and identify an association signal at one of these loci that is more strongly associated with PR interval in African Americans than in Europeans.     

FTO gene associated fatness in relation to body fat distribution and metabolic traits throughout a broad range of fatness

Background

A common single nucleotide polymorphism (SNP) of FTO (rs9939609, T/A) is associated with total body fatness. We investigated the association of this SNP with abdominal and peripheral fatness and obesity-related metabolic traits in middle-aged men through a broad range of fatness present already in adolescence.

Methodology/Principal Findings

Obese young Danish men (n = 753, BMI≥31.0 kg/m²) and a randomly selected group (n = 879) from the same population were examined in three surveys (mean age 35, 46 and 49 years, respectively). The traits included anthropometrics, body composition, oral glucose tolerance test, blood lipids, blood pressure, fibrinogen and aspartate aminotransferase. Logistic regression analysis was used to assess the age-adjusted association between the phenotypes and the odds ratios for the FTO rs9939609 (TT and TA genotype versus the AA genotype), for anthropometrics and body composition estimated per unit z-score. BMI was strongly associated with the AA genotype in all three surveys: OR = 1.17, p = 1.1*10⁻⁶, OR = 1.20, p = 1.7*10⁻⁷, OR = 1.17, p = 3.4*10⁻³, respectively. Fat body mass index was also associated with the AA genotype (OR = 1.21, p = 4.6*10⁻⁷ and OR = 1.21, p = 1.0*10⁻³). Increased abdominal fatness was associated with the AA genotype when measured as waist circumference (OR = 1.21, p = 2.2*10⁻⁶ and OR = 1.19, p = 5.9*10⁻³), sagittal abdominal diameter (OR = 1.17, p = 1.3*10⁻⁴ and OR = 1.18, p = 0.011) and intra-abdominal adipose tissue (OR = 1.21, p = 0.005). Increased peripheral fatness measured as hip circumference (OR = 1.19, p = 1.3*10⁻⁵ and OR = 1.18, p = 0.004) and lower body fat mass (OR = 1.26, p = 0.002) was associated with the AA genotype. The AA genotype was significantly associated with decreased Stumvoll insulin sensitivity index (OR = 0.93, p = 0.02) and with decreased non-fasting plasma HDL-cholesterol (OR = 0.57, p = 0.037), but not with any other of the metabolic traits. However, all significant results for both body fat distribution and metabolic traits were explained by a mediating effect of total fat mass.

Conclusion

The association of the examined FTO SNP to general fatness throughout the range of fatness was confirmed, and this association explains the relation between the SNP and body fat distribution and decreased insulin sensitivity and HDL-cholesterol. The SNP was not significantly associated with other metabolic traits suggesting that they are not derived from the general accumulation of body fat.