Johannah L. Butler

Identification of ten loci associated with height highlights new biological pathways in human growth

Height is a classic polygenic trait, reflecting the combined influence of multiple as-yet-undiscovered genetic factors. We carried out a meta-analysis of genome-wide association study data of height from 15,821 individuals at 2.2 million SNPs, and followed up the strongest findings in >10,000 subjects. Ten newly identified and two previously reported loci were strongly associated with variation in height (P values from 4 × 10−7 to 8 × 10−22). Together, these 12 loci account for ∼2% of the population variation in height. Individuals with ≤8 height-increasing alleles and ≥16 height-increasing alleles differ in height by ∼3.5 cm. The newly identified loci, along with several additional loci with strongly suggestive associations, encompass both strong biological candidates and unexpected genes, and highlight several pathways (let-7 targets, chromatin remodeling proteins and Hedgehog signaling) as important regulators of human stature. These results expand the picture of the biological regulation of human height and of the genetic architecture of this classical complex trait.

Discerning the ancestry of European Americans in genetic association studies

European Americans are often treated as a homogeneous group, but in fact form a structured population due to historical immigration of diverse source populations. Discerning the ancestry of European Americans genotyped in association studies is important in order to prevent false-positive or false-negative associations due to population stratification and to identify genetic variants whose contribution to disease risk differs across European ancestries. Here, we investigate empirical patterns of population structure in European Americans, analyzing 4,198 samples from four genome-wide association studies to show that components roughly corresponding to northwest European, southeast European, and Ashkenazi Jewish ancestry are the main sources of European American population structure. Building on this insight, we constructed a panel of 300 validated markers that are highly informative for distinguishing these ancestries. We demonstrate that this panel of markers can be used to correct for stratification in association studies that do not generate dense genotype data.

PNPLA3 variants specifically confer increased risk for histologic nonalcoholic fatty liver disease but not metabolic disease

Single nucleotide polymorphisms (SNPs) near 7 loci have been associated with liver function tests or with liver steatosis by magnetic resonance spectroscopy. In this study we aim to test whether these SNPs influence the risk of histologically‐confirmed nonalcoholic fatty liver disease (NAFLD). We tested the association of histologic NAFLD with SNPs at 7 loci in 592 cases of European ancestry from the Nonalcoholic Steatohepatitis Clinical Research Network and 1405 ancestry‐matched controls. The G allele of rs738409 in PNPLA3 was associated with increased odds of histologic NAFLD (odds ratio [OR] = 3.26, 95% confidence intervals [CI] = 2.11‐7.21; P = 3.6 × 10−43). In a case only analysis of G allele of rs738409 in PNPLA3 was associated with a decreased risk of zone 3 centered steatosis (OR = 0.46, 95% CI = 0.36‐0.58; P = 5.15 × 10−11). We did not observe any association of this variant with body mass index, triglyceride levels, high‐ and low‐density lipoprotein levels, or diabetes (P > 0.05). None of the variants at the other 6 loci were associated with NAFLD. Conclusion: Genetic variation at PNPLA3 confers a markedly increased risk of increasingly severe histological features of NAFLD, without a strong effect on metabolic syndrome component traits. (HEPATOLOGY 2010)

Transferability of tag SNPs in genetic association studies in multiple populations

A general question for linkage disequilibrium–based association studies is how power to detect an association is compromised when tag SNPs are chosen from data in one population sample and then deployed in another sample. Specifically, it is important to know how well tags picked from the HapMap DNA samples capture the variation in other samples. To address this, we collected dense data uniformly across the four HapMap population samples and eleven other population samples. We picked tag SNPs using genotype data we collected in the HapMap samples and then evaluated the effective coverage of these tags in comparison to the entire set of common variants observed in the other samples. We simulated case-control association studies in the non-HapMap samples under a disease model of modest risk, and we observed little loss in power. These results demonstrate that the HapMap DNA samples can be used to select tags for genome-wide association studies in many samples around the world.

The Association of a SNP Upstream of INSIG2 with Body Mass Index is Reproduced in Several but Not All Cohorts

A SNP upstream of the INSIG2 gene, rs7566605, was recently found to be associated with obesity as measured by body mass index (BMI) by Herbert and colleagues. The association between increased BMI and homozygosity for the minor allele was first observed in data from a genome-wide association scan of 86,604 SNPs in 923 related individuals from the Framingham Heart Study offspring cohort. The association was reproduced in four additional cohorts, but was not seen in a fifth cohort. To further assess the general reproducibility of this association, we genotyped rs7566605 in nine large cohorts from eight populations across multiple ethnicities (total n = 16,969). We tested this variant for association with BMI in each sample under a recessive model using family-based, population-based, and case-control designs. We observed a significant (p < 0.05) association in five cohorts but saw no association in three other cohorts. There was variability in the strength of association evidence across examination cycles in longitudinal data from unrelated individuals in the Framingham Heart Study Offspring cohort. A combined analysis revealed significant independent validation of this association in both unrelated (p = 0.046) and family-based (p = 0.004) samples. The estimated risk conferred by this allele is small, and could easily be masked by small sample size, population stratification, or other confounders. These validation studies suggest that the original association is less likely to be spurious, but the failure to observe an association in every data set suggests that the effect of SNP rs7566605 on BMI may be heterogeneous across population samples.

Genome-wide association of anthropometric traits in African- and African-derived populations

Genome-wide association (GWA) studies have identified common variants that are associated with a variety of traits and diseases, but most studies have been performed in European-derived populations. Here, we describe the first genome-wide analyses of imputed genotype and copy number variants (CNVs) for anthropometric measures in African-derived populations: 1188 Nigerians from Igbo-Ora and Ibadan, Nigeria, and 743 African-Americans from Maywood, IL. To improve the reach of our study, we used imputation to estimate genotypes at approximately 2.1 million single-nucleotide polymorphisms (SNPs) and also tested CNVs for association. No SNPs or common CNVs reached a genome-wide significance level for association with height or body mass index (BMI), and the best signals from a meta-analysis of the two cohorts did not replicate in approximately 3700 African-Americans and Jamaicans. However, several loci previously confirmed in European populations showed evidence of replication in our GWA panel of African-derived populations, including variants near IHH and DLEU7 for height and MC4R for BMI. Analysis of global burden of rare CNVs suggested that lean individuals possess greater total burden of CNVs, but this finding was not supported in an independent European population. Our results suggest that there are not multiple loci with strong effects on anthropometric traits in African-derived populations and that sample sizes comparable to those needed in European GWA studies will be required to identify replicable associations. Meta-analysis of this data set with additional studies in African-ancestry populations will be helpful to improve power to detect novel associations.

Fine mapping of the association with obesity at the FTO locus in African-derived populations

Genome-wide association studies have identified many common genetic variants that are associated with polygenic traits, and have typically been performed with individuals of recent European ancestry. In these populations, many common variants are tightly correlated, with the perfect or near-perfect proxies for the functional or true variant showing equivalent evidence of association, considerably limiting the resolution of fine mapping. Populations with recent African ancestry often have less extensive and/or different patterns of linkage disequilibrium (LD), and have been proposed to be useful in fine-mapping studies. Here, we strongly replicate and fine map in populations of predominantly African ancestry the association between variation at the FTO locus and body mass index (BMI) that is well established in populations of European ancestry. We genotyped single nucleotide polymorphisms that are correlated with the signal of association in individuals of European ancestry but that have varying degrees of correlation in African-derived individuals. Most of the variants, including one previously proposed as functionally important, have no significant association with BMI, but two variants, rs3751812 and rs9941349, show strong evidence of association (P = 2.58 x 10(-6) and 3.61 x 10(-6) in a meta-analysis of 9881 individuals). Thus, we have both strongly replicated this association in African-ancestry populations and narrowed the list of potentially causal variants to those that are correlated with rs3751812 and rs9941349 in African-derived populations. This study illustrates the potential of using populations with different LD patterns to fine map associations and helps pave the way for genetically guided functional studies at the FTO locus.

Association Studies of Common Variants in 10 Hypogonadotropic Hypogonadism Genes with Age at Menarche

CONTEXT Although the timing of puberty is a highly heritable trait, little is known about the genes that regulate pubertal timing in the general population. Several genes have been identified that, when mutated, cause disorders of delayed or absent puberty such as hypogonadotropic hypogonadism (HH). OBJECTIVE Because severe variants in HH-related genes cause a severe puberty phenotype, we hypothesized that common subtle variation in these genes could contribute to the population variation in pubertal timing. DESIGN We assessed common genetic variation in 10 HH-related genes in 1801 women from the Hawaii and Los Angeles Multiethnic Cohort with either early (age<11 yr) or late (age>14 yr) menarche and in other replication samples. In addition to these common variants, we also studied the most frequently reported HH mutations to assess their role in the population variation in pubertal timing. SETTING AND PATIENTS/OTHER PARTICIPANTS: Within the general community, 1801 women from the Hawaii and Los Angeles Multiethnic Cohort participated. MAIN OUTCOME MEASURES We assessed the association of genetic variation with age at menarche. RESULTS We found no significant association between any of the variants tested and age at menarche, although we cannot rule out modest effects of these variants or of other variants at long distances from the coding region. In several self-reported racial/ethnic groups represented in our study, we observed an association between estimated genetic ancestry and age at menarche. CONCLUSIONS Our results suggest that common variants near 10 HH-related loci do not play a substantial role in the regulation of age at menarche in the general population.

Rapid assessment of genetic ancestry in populations of unknown origin by genome-wide genotyping of pooled samples.

As we move forward from the current generation of genome-wide association (GWA) studies, additional cohorts of different ancestries will be studied to increase power, fine map association signals, and generalize association results to additional populations. Knowledge of genetic ancestry as well as population substructure will become increasingly important for GWA studies in populations of unknown ancestry. Here we propose genotyping pooled DNA samples using genome-wide SNP arrays as a viable option to efficiently and inexpensively estimate admixture proportion and identify ancestry informative markers (AIMs) in populations of unknown origin. We constructed DNA pools from African American, Native Hawaiian, Latina, and Jamaican samples and genotyped them using the Affymetrix 6.0 array. Aided by individual genotype data from the African American cohort, we established quality control filters to remove poorly performing SNPs and estimated allele frequencies for the remaining SNPs in each panel. We then applied a regression-based method to estimate the proportion of admixture in each cohort using the allele frequencies estimated from pooling and populations from the International HapMap Consortium as reference panels, and identified AIMs unique to each population. In this study, we demonstrated that genotyping pooled DNA samples yields estimates of admixture proportion that are both consistent with our knowledge of population history and similar to those obtained by genotyping known AIMs. Furthermore, through validation by individual genotyping, we demonstrated that pooling is quite effective for identifying SNPs with large allele frequency differences (i.e., AIMs) and that these AIMs are able to differentiate two closely related populations (HapMap JPT and CHB).

An age-dependent diet-modified effect of the PPARγ Pro12Ala polymorphism in children

Variation in the peroxisome proliferator-activated receptor γ gene alters the risk for adiposity in adults, with evidence of interaction with diet. We investigated the age-related association between the Pro12Ala variant (rs1801282) and diet in obesity-related traits in children. The Pro12Ala variant was assayed in 2102 young children aged 1 to 6 years and in 794 periadolescent children aged 10 to 12 years of Greek origin. In both cohorts, no differences were found in obesity traits between the Ala allele carriers and Pro/Pro homozygotes. Sex-stratified analysis showed that, in periadolescent boys, Ala carriers exhibited lower measures of skinfolds (triceps: 16.9 ± 6.9 vs 19.4 ± 7.9 mm, P = .01; subscapular: 9.6 ± 4.5 vs 11.2 ± 5.4 mm, P = .02). On the other hand, young girls who were Ala carriers presented higher measures of triceps skinfold thickness (10.5 ± 3.0 vs 9.9 ± 2.8 mm, P = .04). Nominal gene-diet interactions were revealed in periadolescents for saturated fatty acid (SFA) intake and skinfolds (P for interaction = .05). In Pro/Pro homozygous young girls, SFA and total fat (TF) intake was positively associated with higher body mass index (BMI) (P = .01), waist circumference (P = .02), and skinfold thickness (triceps-SFA: P = 10⁻⁵, triceps-TF: P = 10⁻⁹, subscapular-SFA: P = 10⁻⁶, subscapular-TF: P = 10⁻⁴). For Pro/Pro homozygotes, unsaturated fat intake was inversely associated with BMI (P = .04) in young girls, and with BMI (P = .03), waist circumference (P = .03), and triceps (P = .02) in periadolescent boys. Our results suggest that adiposity in children is influenced by the Pro12Ala polymorphism in a sex-specific and age-dependent manner. We also demonstrate evidence of an age-dependent gene-diet (SFA, TF) interaction, suggesting that the type of fat intake modifies the effect of the Pro12 allele on obesity-related measures.