Summary

In the presence of epistasis multilocus association tests of human complex traits can provide powerful methods to detect susceptibility variants. We undertook multilocus analyses in 1924 type 2 diabetes cases and 2938 controls from the Wellcome Trust Case Control Consortium (WTCCC). We performed a two-dimensional genome-wide association (GWA) scan using joint two-locus tests of association including main and epistatic effects in 70,236 markers tagging common variants. We found two-locus association at 79 SNP-pairs at a Bonferroni-corrected P-value = 0.05 (uncorrected P-value = 2.14 × 10−11). The 79 pair-wise results always contained rs11196205 in TCF7L2 paired with 79 variants including confirmed variants in FTO, TSPAN8, and CDKAL1, which are associated in the absence of epistasis. However, the majority (82%) of the 79 variants did not have compelling single-locus association signals (P-value = 5 × 10−4). Analyses conditional on the single-locus effects at TCF7L2 established that the joint two-locus results could be attributed to single-locus association at TCF7L2 alone. Interaction analyses among the peak 80 regions and among 23 previously established diabetes candidate genes identified five SNP-pairs with case-control and case-only epistatic signals. Our results demonstrate the feasibility of systematic scans in GWA data, but confirm that single-locus association can underlie and obscure multilocus findings.

OBJECTIVE

Despite the clinical importance of an accurate diagnosis in individuals with monogenic forms of diabetes, restricted access to genetic testing leaves many patients with undiagnosed diabetes. Recently, common variation near the HNF1 homeobox A (HNF1A) gene was shown to influence C-reactive protein levels in healthy adults. We hypothesized that serum levels of high-sensitivity C-reactive protein (hs-CRP) could represent a clinically useful biomarker for the identification of HNF1A mutations causing maturity-onset diabetes of the young (MODY).

RESEARCH DESIGN AND METHODS

Serum hs-CRP was measured in subjects with HNF1A-MODY (n = 31), autoimmune diabetes (n = 316), type 2 diabetes (n = 240), and glucokinase (GCK) MODY (n = 24) and in nondiabetic individuals (n = 198). The discriminative accuracy of hs-CRP was evaluated through receiver operating characteristic (ROC) curve analysis, and performance was compared with standard diagnostic criteria. Our primary analyses excluded ∼11% of subjects in whom the single available hs-CRP measurement was >10 mg/l.

RESULTS

Geometric mean (SD range) hs-CRP levels were significantly lower (P ≤ 0.009) for HNF1A-MODY individuals, 0.20 (0.03–1.14) mg/l, than for any other group: autoimmune diabetes 0.58 (0.10–2.75) mg/l, type 2 diabetes 1.33 (0.28–6.14) mg/l, GCK-MODY 1.01 (0.19–5.33) mg/l, and nondiabetic 0.48 (0.10–2.42) mg/l. The ROC-derived C-statistic for discriminating HNF1A-MODY and type 2 diabetes was 0.8. Measurement of hs-CRP, either alone or in combination with current diagnostic criteria, was superior to current diagnostic criteria alone. Sensitivity and specificity for the combined criteria approached 80%.

CONCLUSIONS

Serum hs-CRP levels are markedly lower in HNF1A-MODY than in other forms of diabetes. hs-CRP has potential as a widely available, cost-effective screening test to support more precise targeting of MODY diagnostic testing.

OBJECTIVE

Assignment of the correct molecular diagnosis in diabetes is necessary for informed decisions regarding treatment and prognosis. Better clinical markers would facilitate discrimination and prioritization for genetic testing between diabetes subtypes. Serum 1,5 anhydroglucitol (1,5AG) levels were reported to differentiate maturity-onset diabetes of the young due to HNF1A mutations (HNF1A-MODY) from type 2 diabetes, but this requires further validation. We evaluated serum 1,5AG in a range of diabetes subtypes as an adjunct for defining diabetes etiology.

RESEARCH DESIGN AND METHODS

1,5AG was measured in U.K. subjects with: HNF1A-MODY (n = 23), MODY due to glucokinase mutations (GCK-MODY, n = 23), type 1 diabetes (n = 29), latent autoimmune diabetes in adults (LADA, n = 42), and type 2 diabetes (n = 206). Receiver operating characteristic curve analysis was performed to assess discriminative accuracy of 1,5AG for diabetes etiology.

RESULTS

Mean (SD range) 1,5AG levels were: GCK-MODY 13.06 μg/ml (5.74–29.74), HNF1A-MODY 4.23 μg/ml (2.12–8.44), type 1 diabetes 3.09 μg/ml (1.45–6.57), LADA 3.46 μg/ml (1.42–8.45), and type 2 diabetes 5.43 (2.12–13.23). Levels in GCK-MODY were higher than in other groups (P < 10−4 vs. each group). HNF1A-MODY subjects showed no difference in unadjusted 1,5AG levels from type 2 diabetes, type 1 diabetes, and LADA. Adjusting for A1C revealed a difference between HNF1A-MODY and type 2 diabetes (P = 0.001). The discriminative accuracy of unadjusted 1,5AG levels was 0.79 for GCK-MODY versus type 2 diabetes and 0.86 for GCK-MODY versus HNF1A-MODY but was only 0.60 for HNF1A-MODY versus type 2 diabetes.

CONCLUSIONS

In our dataset, serum 1,5AG performed well in discriminating GCK-MODY from other diabetes subtypes, particularly HNF1A-MODY. Measurement of 1,5AG levels could inform decisions regarding MODY diagnostic testing.

Summary: The Sample avAILability system—SAIL—is a web based application for searching, browsing and annotating biological sample collections or biobank entries. By providing individual-level information on the availability of specific data types (phenotypes, genetic or genomic data) and samples within a collection, rather than the actual measurement data, resource integration can be facilitated. A flexible data structure enables the collection owners to provide descriptive information on their samples using existing or custom vocabularies. Users can query for the available samples by various parameters combining them via logical expressions. The system can be scaled to hold data from millions of samples with thousands of variables.

Availability: SAIL is available under Aferro-GPL open source license: https://github.com/sail.

Contact: gostev@ebi.ac.uk, support@simbioms.org

Supplementary information: Supplementary data are available at Bioinformatics online and from http://www.simbioms.org.

The association between variation in the fat mass and obesity-associated (FTO) gene and adulthood body mass index (BMI; weight (kg)/height (m)2) is well-replicated. More thorough analyses utilizing phenotypic data over the life course may deepen our understanding of the development of BMI and thus help in the prevention of obesity. The authors used a structural equation modeling approach to explore the network of variables associated with BMI from the prenatal period to age 31 years (1965–1997) in 4,435 subjects from the Northern Finland Birth Cohort 1966. The use of structural equation modeling permitted the easy inclusion of variables with missing values in the analyses without separate imputation steps, as well as differentiation between direct and indirect effects. There was an association between the FTO single nucleotide polymorphism rs9939609 and BMI at age 31 years that persisted after controlling for several relevant factors during the life course. The total effect of the FTO variant on adult BMI was mostly composed of the direct effect, but a notable part was also arising indirectly via its effects on earlier BMI development. In addition to well-established genetic determinants, many life-course factors such as physical activity, in spite of not showing mediation or interaction, had a strong independent effect on BMI.

Background

Vitamin D is crucial for maintaining musculoskeletal health. Recently, vitamin D insufficiency has been linked to a number of extraskeletal disorders, including diabetes, cancer, and cardiovascular disease. Determinants of circulating 25-hydroxyvitamin D (25-OH D) include sun exposure and dietary intake, but its high heritability suggests that genetic determinants may also play a role.

Methods

We performed a genome-wide association study of 25-OH D among ∼30,000 individuals of European descent from 15 cohorts. Five cohorts were designated as discovery cohorts (n=16,125), five as in silico replication cohorts (n=9,366), and five as de novo replication cohorts (n=8,378). Association results were combined using z-score-weighted meta-analysis. Vitamin D insufficiency was defined as 25-OH D <75 nmol/L or <50 nmol/L.

Findings

Variants at three loci reached genome-wide significance in the discovery cohorts, and were confirmed in the replication cohorts: 4p12 (overall P=1.9 × 10-109 for rs2282679, in GC); 11q12 (P=2.1 × 10-27 for rs12785878, near DHCR7); 11p15 (P=3.3 × 10-20 for rs10741657, near CYP2R1). Variants at an additional locus (20q13, CYP24A1) were genome-wide significant in the pooled sample (P=6.0 × 10-10 for rs6013897). A genotype score was constructed using the three confirmed variants. Those in the top quartile of genotype scores had 2- to 2.5-fold elevated odds of vitamin D insufficiency (P≤1 × 10-26).

Interpretation

Variants near genes involved in cholesterol synthesis (DHCR7), hydroxylation (CYP2R1, CYP24A1), and vitamin D transport (GC) influence vitamin D status. Genetic variation at these loci identifies individuals of European descent who have substantially elevated risk of vitamin D insufficiency.

Recent multi-dimensional approaches to the study of complex disease have revealed powerful insights into how genetic and epigenetic factors may underlie their aetiopathogenesis. We examined genotype-epigenotype interactions in the context of Type 2 Diabetes (T2D), focussing on known regions of genomic susceptibility. We assayed DNA methylation in 60 females, stratified according to disease susceptibility haplotype using previously identified association loci. CpG methylation was assessed using methylated DNA immunoprecipitation on a targeted array (MeDIP-chip) and absolute methylation values were estimated using a Bayesian algorithm (BATMAN). Absolute methylation levels were quantified across LD blocks, and we identified increased DNA methylation on the FTO obesity susceptibility haplotype, tagged by the rs8050136 risk allele A (p = 9.40×10−4, permutation p = 1.0×10−3). Further analysis across the 46 kb LD block using sliding windows localised the most significant difference to be within a 7.7 kb region (p = 1.13×10−7). Sequence level analysis, followed by pyrosequencing validation, revealed that the methylation difference was driven by the co-ordinated phase of CpG-creating SNPs across the risk haplotype. This 7.7 kb region of haplotype-specific methylation (HSM), encapsulates a Highly Conserved Non-Coding Element (HCNE) that has previously been validated as a long-range enhancer, supported by the histone H3K4me1 enhancer signature. This study demonstrates that integration of Genome-Wide Association (GWA) SNP and epigenomic DNA methylation data can identify potential novel genotype-epigenotype interactions within disease-associated loci, thus providing a novel route to aid unravelling common complex diseases.

We performed a genome-wide association study (GWAS) in 1705 Parkinson's disease (PD) UK patients and 5175 UK controls, the largest sample size so far for a PD GWAS. Replication was attempted in an additional cohort of 1039 French PD cases and 1984 controls for the 27 regions showing the strongest evidence of association (P< 10−4). We replicated published associations in the 4q22/SNCA and 17q21/MAPT chromosome regions (P< 10−10) and found evidence for an additional independent association in 4q22/SNCA. A detailed analysis of the haplotype structure at 17q21 showed that there are three separate risk groups within this region. We found weak but consistent evidence of association for common variants located in three previously published associated regions (4p15/BST1, 4p16/GAK and 1q32/PARK16). We found no support for the previously reported SNP association in 12q12/LRRK2. We also found an association of the two SNPs in 4q22/SNCA with the age of onset of the disease.

It has been recently hypothesized that many of the signals detected in genome-wide association studies (GWAS) to T2D and other diseases, despite being observed to common variants, might in fact result from causal mutations that are rare. One prediction of this hypothesis is that the allelic associations should be population-specific, as the causal mutations arose after the migrations that established different populations around the world. We selected 19 common variants found to be reproducibly associated to T2D risk in European populations and studied them in a large multiethnic case-control study (6,142 cases and 7,403 controls) among men and women from 5 racial/ethnic groups (European Americans, African Americans, Latinos, Japanese Americans, and Native Hawaiians). In analysis pooled across ethnic groups, the allelic associations were in the same direction as the original report for all 19 variants, and 14 of the 19 were significantly associated with risk. In summing the number of risk alleles for each individual, the per-allele associations were highly statistically significant (P<10−4) and similar in all populations (odds ratios 1.09–1.12) except in Japanese Americans the estimated effect per allele was larger than in the other populations (1.20; Phet = 3.8×10−4). We did not observe ethnic differences in the distribution of risk that would explain the increased prevalence of type 2 diabetes in these groups as compared to European Americans. The consistency of allelic associations in diverse racial/ethnic groups is not predicted under the hypothesis of Goldstein regarding “synthetic associations” of rare mutations in T2D.

Author Summary

Single rare causal alleles and/or collections of multiple rare alleles have been suggested to create “synthetic associations” with common variants in genome-wide association studies (GWAS). This model predicts that associations with common variants will not be consistent across populations. In this study, we examined 19 T2D variants for association with T2D risk in 6,142 cases and 7,403 controls from five racial/ethnic populations in the Multiethnic Cohort (European Americans, African Americans, Latinos, Japanese Americans, and Native Hawaiians). In racial/ethnic pooled analysis, all 19 variants were associated with T2D risk in the same direction as previous reports in Europeans, and the sum total of risk variants was significantly associated with T2D risk in each racial/ethnic group. The consistent associations across populations do not support the Goldstein hypothesis that rare causal alleles underlie GWAS signals. We also did not find evidence that these markers underlie racial/ethnic disparities in T2D prevalence. Large-scale GWAS and sequencing studies in these populations are necessary in order to both improve the current set of markers at these risk loci and identify new risk variants for T2D that may be difficult, or impossible, to detect in European populations.

Chronic kidney disease (CKD) is a worldwide public health problem that is associated with substantial morbidity and mortality. To search for sequence variants that associate with CKD, we conducted a genome-wide association study (GWAS) that included a total of 3,203 Icelandic cases and 38,782 controls. We observed an association between CKD and a variant with 80% population frequency, rs4293393-T, positioned next to the UMOD gene (GeneID: 7369) on chromosome 16p12 (OR = 1.25, P = 4.1×10−10). This gene encodes uromodulin (Tamm-Horsfall protein), the most abundant protein in mammalian urine. The variant also associates significantly with serum creatinine concentration (SCr) in Icelandic subjects (N = 24,635, P = 1.3×10−23) but not in a smaller set of healthy Dutch controls (N = 1,819, P = 0.39). Our findings validate the association between the UMOD variant and both CKD and SCr recently discovered in a large GWAS. In the Icelandic dataset, we demonstrate that the effect on SCr increases substantially with both age (P = 3.0×10−17) and number of comorbid diseases (P = 0.008). The association with CKD is also stronger in the older age groups. These results suggest that the UMOD variant may influence the adaptation of the kidney to age-related risk factors of kidney disease such as hypertension and diabetes. The variant also associates with serum urea (P = 1.0×10−6), uric acid (P = 0.0064), and suggestively with gout. In contrast to CKD, the UMOD variant confers protection against kidney stones when studied in 3,617 Icelandic and Dutch kidney stone cases and 43,201 controls (OR = 0.88, P = 5.7×10−5).

Author Summary

Chronic kidney disease (CKD) is a common condition that is associated with substantial morbidity and mortality and has been recognized as a major public health problem worldwide. Common causes of CKD include hypertension, diabetes, and inflammatory disorders. Previous studies have shown a significant genetic contribution to kidney disease and a recent genome-wide association study yielded a variant in the UMOD gene that affects the risk of CKD. Here, we replicate the association between UMOD and CKD in an independent analysis. We also demonstrate for the first time an interaction between the UMOD variant and age that suggests that this variant may adversely affect the aging kidney and its adaptation to age-related risk factors of kidney disease, such as hypertension and diabetes. Furthermore, we show that the UMOD variant that affects risk of CKD also provides protection against kidney stone disease.

Background

Genome-wide association studies have found type 2 diabetes-associated variants in the HNF1B gene to exhibit reciprocal associations with prostate cancer risk. We aimed to identify whether these variants may have an effect on cancer risk in general versus a specific effect on prostate cancer only.

Methodology/Principal Findings

In a collaborative analysis, we collected data from GWAS of cancer phenotypes for the frequently reported variants of HNF1B, rs4430796 and rs7501939, which are in linkage disequilibrium (r2 = 0.76, HapMap CEU). Overall, the analysis included 16 datasets on rs4430796 with 19,640 cancer cases and 21,929 controls; and 21 datasets on rs7501939 with 26,923 cases and 49,085 controls. Malignancies other than prostate cancer included colorectal, breast, lung and pancreatic cancers, and melanoma. Meta-analysis showed large between-dataset heterogeneity that was driven by different effects in prostate cancer and other cancers. The per-T2D-risk-allele odds ratios (95% confidence intervals) for rs4430796 were 0.79 (0.76, 0.83)] per G allele for prostate cancer (p<10−15 for both); and 1.03 (0.99, 1.07) for all other cancers. Similarly for rs7501939 the per-T2D-risk-allele odds ratios (95% confidence intervals) were 0.80 (0.77, 0.83) per T allele for prostate cancer (p<10−15 for both); and 1.00 (0.97, 1.04) for all other cancers. No malignancy other than prostate cancer had a nominally statistically significant association.

Conclusions/Significance

The examined HNF1B variants have a highly specific effect on prostate cancer risk with no apparent association with any of the other studied cancer types.

Previous studies have successfully identified genetic variants in several genes associated with human iris (eye) color; however, they all used simplified categorical trait information. Here, we quantified continuous eye color variation into hue and saturation values using high-resolution digital full-eye photographs and conducted a genome-wide association study on 5,951 Dutch Europeans from the Rotterdam Study. Three new regions, 1q42.3, 17q25.3, and 21q22.13, were highlighted meeting the criterion for genome-wide statistically significant association. The latter two loci were replicated in 2,261 individuals from the UK and in 1,282 from Australia. The LYST gene at 1q42.3 and the DSCR9 gene at 21q22.13 serve as promising functional candidates. A model for predicting quantitative eye colors explained over 50% of trait variance in the Rotterdam Study. Over all our data exemplify that fine phenotyping is a useful strategy for finding genes involved in human complex traits.

Author Summary

We measured human eye color to hue and saturation values from high-resolution, digital, full-eye photographs of several thousand Dutch Europeans. This quantitative approach, which is extremely cost-effective, portable, and time efficient, revealed that human eye color varies along more dimensions than the one represented by the blue-green-brown categories studied previously. Our work represents the first genome-wide study of quantitative human eye color. We clearly identified 3 new loci, LYST, 17q25.3, TTC3/DSCR9, in contributing to the natural and subtle eye color variation along multiple dimensions, providing new leads towards a more detailed understanding of the genetic basis of human eye color. Our quantitative prediction model explained over 50% of eye color variance, representing the highest accuracy achieved so far in genomic prediction of human complex and quantitative traits, with relevance for future forensic applications.

Background

The TNF/LTA locus has been a long-standing T2D candidate gene. Several studies have examined association of TNF/LTA SNPs with T2D but the majority have been small-scale and produced no convincing evidence of association. The purpose of this study is to examine T2D association of tag SNPs in the TNF/LTA region capturing the majority of common variation in a large-scale sample set of UK/Irish origin.

Methods

This study comprised a case-control (1520 cases and 2570 control samples) and a family-based component (423 parent-offspring trios). Eleven tag SNPs (rs928815, rs909253, rs746868, rs1041981 (T60N), rs1800750, rs1800629 (G-308A), rs361525 (G-238A), rs3093662, rs3093664, rs3093665, and rs3093668) were selected across the TNF/LTA locus and genotyped using a fluorescence-based competitive allele specific assay. Quality control of the obtained genotypes was performed prior to single- and multi-point association analyses under the additive model.

Results

We did not find any consistent SNP associations with T2D in the case-control or family-based datasets.

Conclusions

The present study, designed to analyse a set of tag SNPs specifically selected to capture the majority of common variation in the TNF/LTA gene region, found no robust evidence for association with T2D. To investigate the presence of smaller effects of TNF/LTA gene variation with T2D, a large-scale meta-analysis will be required.

Lipoprotein-associated phospholipase A2 (Lp-PLA2) is an emerging risk factor and therapeutic target for cardiovascular disease. The activity and mass of this enzyme are heritable traits, but major genetic determinants have not been explored in a systematic, genome-wide fashion. We carried out a genome-wide association study of Lp-PLA2 activity and mass in 6,668 Caucasian subjects from the population-based Framingham Heart Study. Clinical data and genotypes from the Affymetrix 550K SNP array were obtained from the open-access Framingham SHARe project. Each polymorphism that passed quality control was tested for associations with Lp-PLA2 activity and mass using linear mixed models implemented in the R statistical package, accounting for familial correlations, and controlling for age, sex, smoking, lipid-lowering-medication use, and cohort. For Lp-PLA2 activity, polymorphisms at four independent loci reached genome-wide significance, including the APOE/APOC1 region on chromosome 19 (p = 6×10−24); CELSR2/PSRC1 on chromosome 1 (p = 3×10−15); SCARB1 on chromosome 12 (p = 1×10−8) and ZNF259/BUD13 in the APOA5/APOA1 gene region on chromosome 11 (p = 4×10−8). All of these remained significant after accounting for associations with LDL cholesterol, HDL cholesterol, or triglycerides. For Lp-PLA2 mass, 12 SNPs achieved genome-wide significance, all clustering in a region on chromosome 6p12.3 near the PLA2G7 gene. Our analyses demonstrate that genetic polymorphisms may contribute to inter-individual variation in Lp-PLA2 activity and mass.

Author Summary

Blood levels of lipoprotein-associated phospholipase A2 (Lp-PLA2) show a strong association with atherosclerosis in humans. This enzyme is made by certain cells of the immune system, associates with lipoproteins (HDL and LDL), and is thought to be involved in inflammation. Studies have shown that Lp-PLA2 is a good predictor of cardiovascular disease, independent of HDL and LDL cholesterol levels. This has led to the development of drugs aimed at inhibiting Lp-PLA2 as a way to treat or prevent cardiovascular disease. The activity and mass of Lp-PLA2 are heritable traits, but major genetic determinants have not been explored in a systematic fashion. We examined genetic variants across the human genome to identify genes influencing Lp-PLA2 activity and mass. We studied 6,668 Caucasian subjects from the population-based Framingham Heart Study. Clinical data and genetic data on 550,000 genetic variants were available for association analysis. There was no overlap in the most significantly associated SNPs for activity and mass. We identified four distinct gene regions showing highly significant associations with Lp-PLA2 activity, all of which are known to include genes involved in cholesterol metabolism. The only locus associated with Lp-PLA2 mass was a region harboring PLA2G7, the gene that encodes lipoprotein-associated phospholipase A2.

Meta-analyses of population-based genome-wide association studies (GWAS) in adults have recently led to the detection of new genetic loci for obesity. Here we aimed to discover additional obesity loci in extremely obese children and adolescents. We also investigated if these results generalize by estimating the effects of these obesity loci in adults and in population-based samples including both children and adults. We jointly analysed two GWAS of 2,258 individuals and followed-up the best, according to lowest p-values, 44 single nucleotide polymorphisms (SNP) from 21 genomic regions in 3,141 individuals. After this DISCOVERY step, we explored if the findings derived from the extremely obese children and adolescents (10 SNPs from 5 genomic regions) generalized to (i) the population level and (ii) to adults by genotyping another 31,182 individuals (GENERALIZATION step). Apart from previously identified FTO, MC4R, and TMEM18, we detected two new loci for obesity: one in SDCCAG8 (serologically defined colon cancer antigen 8 gene; p = 1.85×10−8 in the DISCOVERY step) and one between TNKS (tankyrase, TRF1-interacting ankyrin-related ADP-ribose polymerase gene) and MSRA (methionine sulfoxide reductase A gene; p = 4.84×10−7), the latter finding being limited to children and adolescents as demonstrated in the GENERALIZATION step. The odds ratios for early-onset obesity were estimated at ∼1.10 per risk allele for both loci. Interestingly, the TNKS/MSRA locus has recently been found to be associated with adult waist circumference. In summary, we have completed a meta-analysis of two GWAS which both focus on extremely obese children and adolescents and replicated our findings in a large followed-up data set. We observed that genetic variants in or near FTO, MC4R, TMEM18, SDCCAG8, and TNKS/MSRA were robustly associated with early-onset obesity. We conclude that the currently known major common variants related to obesity overlap to a substantial degree between children and adults.

Author Summary

Genome-wide association studies (GWAS) have successfully contributed to the detection of genetic variants involved in body-weight regulation. We jointly analysed two GWAS for early-onset extreme obesity in 2,258 individuals of European origin and followed-up the findings in 3,141 individuals. Evidence for association of markers in two new genetic loci was shown (SDCCAG8 on chromosome 1q43–q44 and between TNKS/MSRA on chromosome 8p23.1). We also re-identified variants in or near FTO, MC4R, and TMEM18 to be associated with extreme obesity. In addition, we assessed the effect of the markers in 31,182 obese, lean, normal weight, and unselected individuals from population-based samples and showed that the variants near FTO, MC4R, TMEM18, and SDCCAG8 were consistently associated with obesity. For variants of TNKS/MSRA, the obesity association was limited to children and adolescents. In summary, we detected two new obesity loci and confirmed that the currently known major common variants related to obesity overlap to a substantial degree between children and adults.

Tooth development is a highly heritable process which relates to other growth and developmental processes, and which interacts with the development of the entire craniofacial complex. Abnormalities of tooth development are common, with tooth agenesis being the most common developmental anomaly in humans. We performed a genome-wide association study of time to first tooth eruption and number of teeth at one year in 4,564 individuals from the 1966 Northern Finland Birth Cohort (NFBC1966) and 1,518 individuals from the Avon Longitudinal Study of Parents and Children (ALSPAC). We identified 5 loci at P<5×10−8, and 5 with suggestive association (P<5×10−6). The loci included several genes with links to tooth and other organ development (KCNJ2, EDA, HOXB2, RAD51L1, IGF2BP1, HMGA2, MSRB3). Genes at four of the identified loci are implicated in the development of cancer. A variant within the HOXB gene cluster associated with occlusion defects requiring orthodontic treatment by age 31 years.

Author Summary

Genome-wide association studies have been used to identify genetic variants conferring susceptibility to diseases, intermediate phenotypes, and physiological traits such as height, hair color, and age at menarche. Here we analyze the NFBC1966 and ALSPAC birth cohorts to investigate the genetic determinants of a key developmental process: primary tooth development. The prospective nature of our studies allows us to exploit accurate measurements of age at first tooth eruption and number of teeth at one year, and also provides the opportunity to assess whether genetic variants affecting these traits are associated with dental problems later in the life course. Of the genes that we find to be associated with primary tooth development, several have established roles in tooth development and growth, and almost half have proposed links with the development of cancer. We find that one of the variants is also associated with occlusion defects requiring orthodontic treatment later in life. Our findings should provide a strong foundation for the study of the genetic architecture of tooth development, which as well as its relevance to medicine and dentistry, may have implications in evolutionary biology since teeth represent important markers of evolution.

It has been suggested that genetic susceptibility plays an important role in the pathogenesis of diabetic nephropathy. A large-scale genotyping analysis of gene-based single nucleotide polymorphisms (SNPs) in Japanese patients with type 2 diabetes identified the gene encoding acetyl-coenzyme A carboxylase beta (ACACB) as a candidate for a susceptibility to diabetic nephropathy; the landmark SNP was found in the intron 18 of ACACB (rs2268388: intron 18 +4139 C > T, p = 1.4×10−6, odds ratio = 1.61, 95% confidence interval [CI]: 1.33–1.96). The association of this SNP with diabetic nephropathy was examined in 9 independent studies (4 from Japan including the original study, one Singaporean, one Korean, and two European) with type 2 diabetes. One case-control study involving European patients with type 1 diabetes was included. The frequency of the T allele for SNP rs2268388 was consistently higher among patients with type 2 diabetes and proteinuria. A meta-analysis revealed that rs2268388 was significantly associated with proteinuria in Japanese patients with type 2 diabetes (p = 5.35×10−8, odds ratio = 1.61, 95% Cl: 1.35–1.91). Rs2268388 was also associated with type 2 diabetes–associated end-stage renal disease (ESRD) in European Americans (p = 6×10−4, odds ratio = 1.61, 95% Cl: 1.22–2.13). Significant association was not detected between this SNP and nephropathy in those with type 1 diabetes. A subsequent in vitro functional analysis revealed that a 29-bp DNA fragment, including rs2268388, had significant enhancer activity in cultured human renal proximal tubular epithelial cells. Fragments corresponding to the disease susceptibility allele (T) had higher enhancer activity than those of the major allele. These results suggest that ACACB is a strong candidate for conferring susceptibility for proteinuria in patients with type 2 diabetes.

Author Summary

Although cumulative epidemiological findings have suggested that genetic susceptibility plays an important role in the pathogenesis of diabetic nephropathy, no gene conferring susceptibility to diabetic nephropathy has been definitively identified. In a large-scale association study of 1,312 Japanese subjects with type 2 diabetes using SNPs from a Japanese SNP database, we show that the T-allele of ACACB rs2268388 is associated with diabetic nephropathy. We also show that the association is consistently observed in patients with type 2 diabetes and proteinuria across different ethnic groups, including populations of European descent. Because a DNA fragment corresponding to the disease susceptibility allele is shown to have higher enhancer activity, we hypothesize that the increase in the expression and/or activity of the encoded acetyl-coenzyme A carboxylase beta contributes to the development and progression of diabetic nephropathy. Our present analysis provides novel insight into the pathogenesis of diabetic nephropathy. This finding is important because diabetic nephropathy is a leading cause of end-stage renal disease and affects life expectancy in subjects with type 2 diabetes.

Background

Readily accessible samples such as peripheral blood or cell lines are increasingly being used in large cohorts to characterise gene expression differences between a patient group and healthy controls. However, cell and RNA isolation procedures and the variety of cell types that make up whole blood can affect gene expression measurements. We therefore systematically investigated global gene expression profiles in peripheral blood from six individuals collected during two visits by comparing five of the following cell and RNA isolation methods: whole blood (PAXgene), peripheral blood mononuclear cells (PBMCs), lymphoblastoid cell lines (LCLs), CD19 and CD20 specific B-cell subsets.

Results

Gene expression measurements were clearly discriminated by isolation method although the reproducibility was high for all methods (range ρ = 0.90-1.00). The PAXgene samples showed a decrease in the number of expressed genes (P < 1*10-16) with higher variability (P < 1*10-16) compared to the other methods. Differentially expressed probes between PAXgene and PBMCs were correlated with the number of monocytes, lymphocytes, neutrophils or erythrocytes. The correlations (ρ = 0.83; ρ = 0.79) of the expression levels of detected probes between LCLs and B-cell subsets were much lower compared to the two B-cell isolation methods (ρ = 0.98). Gene ontology analysis of detected genes showed that genes involved in inflammatory responses are enriched in B-cells CD19 and CD20 whereas genes involved in alcohol metabolic process and the cell cycle were enriched in LCLs.

Conclusion

Gene expression profiles in blood-based samples are strongly dependent on the predominant constituent cell type(s) and RNA isolation method. It is crucial to understand the differences and variability of gene expression measurements between cell and RNA isolation procedures, and their relevance to disease processes, before application in large clinical studies.