Functional significance of single nucleotide polymorphisms in the lactase gene in diverse US patients and evidence for a novel lactase persistence allele at -13909 in those of European ancestry.

OBJECTIVES: Recent data from mainly homogeneous European and African populations implicate a 140-bp region 5' to the transcriptional start site of LCT (the lactase gene) as a regulatory site for lactase persistence and nonpersistence. Because there are no studies of US nonhomogeneous populations, we performed genotype/phenotype analysis of the -13910 and -22018 LCT single nucleotide polymorphisms (SNPs) in New England children, mostly of European ancestry. METHODS: Duodenal biopsies were processed for disaccharidase activities, RNA quantification by reverse transcription polymerase chain reaction (RT-PCR), allelic expression ratios by PCR, and genotyping and SNP analysis. Results were compared with clinical information. RESULTS: Lactase activity and mRNA levels, and sucrase-to-lactase ratios of enzyme activity and mRNA, showed robust correlations with genotype. None of the other LCT SNPs showed as strong a correlation with enzyme or mRNA levels as did -13910. Data were consistent, with the -13910 being the causal sequence variant instead of -22018. Four individuals heterozygous for -13910T/C had allelic expression patterns similar to individuals with -13910C/C genotypes; of these, 2 showed equal LCT expression from the 2 alleles and a novel variant (-13909C>A) associated with lactase persistence. CONCLUSIONS: The identification of -13910C/C genotype is likely to predict lactase nonpersistence, consistent with prior published studies. A -13910T/T genotype will frequently, but not perfectly, predict lactase persistence in this mixed European-ancestry population; a -13910T/C genotype will not predict the phenotype. A long, rare haplotype in 2 individuals with -13910T/C genotype but equal allele-specific expression contains a novel lactase persistence allele present at -13909.

Genome-wide association analysis identifies variants associated with nonalcoholic fatty liver disease that have distinct effects on metabolic traits.

Nonalcoholic fatty liver disease (NAFLD) clusters in families, but the only known common genetic variants influencing risk are near PNPLA3. We sought to identify additional genetic variants influencing NAFLD using genome-wide association (GWA) analysis of computed tomography (CT) measured hepatic steatosis, a non-invasive measure of NAFLD, in large population based samples. Using variance components methods, we show that CT hepatic steatosis is heritable (∼26%-27%) in family-based Amish, Family Heart, and Framingham Heart Studies (n = 880 to 3,070). By carrying out a fixed-effects meta-analysis of genome-wide association (GWA) results between CT hepatic steatosis and ∼2.4 million imputed or genotyped SNPs in 7,176 individuals from the Old Order Amish, Age, Gene/Environment Susceptibility-Reykjavik study (AGES), Family Heart, and Framingham Heart Studies, we identify variants associated at genome-wide significant levels (p<5×10(-8)) in or near PNPLA3, NCAN, and PPP1R3B. We genotype these and 42 other top CT hepatic steatosis-associated SNPs in 592 subjects with biopsy-proven NAFLD from the NASH Clinical Research Network (NASH CRN). In comparisons with 1,405 healthy controls from the Myocardial Genetics Consortium (MIGen), we observe significant associations with histologic NAFLD at variants in or near NCAN, GCKR, LYPLAL1, and PNPLA3, but not PPP1R3B. Variants at these five loci exhibit distinct patterns of association with serum lipids, as well as glycemic and anthropometric traits. We identify common genetic variants influencing CT-assessed steatosis and risk of NAFLD. Hepatic steatosis associated variants are not uniformly associated with NASH/fibrosis or result in abnormalities in serum lipids or glycemic and anthropometric traits, suggesting genetic heterogeneity in the pathways influencing these traits.

The efficacy of detecting variants with small effects on the Affymetrix 6.0 platform using pooled DNA.

Genome-wide genotyping of a cohort using pools rather than individual samples has long been proposed as a cost-saving alternative for performing genome-wide association (GWA) studies. However, successful disease gene mapping using pooled genotyping has thus far been limited to detecting common variants with large effect sizes, which tend not to exist for many complex common diseases or traits. Therefore, for DNA pooling to be a viable strategy for conducting GWA studies, it is important to determine whether commonly used genome-wide SNP array platforms such as the Affymetrix 6.0 array can reliably detect common variants of small effect sizes using pooled DNA. Taking obesity and age at menarche as examples of human complex traits, we assessed the feasibility of genome-wide genotyping of pooled DNA as a single-stage design for phenotype association. By individually genotyping the top associations identified by pooling, we obtained a 14- to 16-fold enrichment of SNPs nominally associated with the phenotype, but we likely missed the top true associations. In addition, we assessed whether genotyping pooled DNA can serve as an inexpensive screen as the second stage of a multi-stage design with a large number of samples by comparing the most cost-effective 3-stage designs with 80% power to detect common variants with genotypic relative risk of 1.1, with and without pooling. Given the current state of the specific technology we employed and the associated genotyping costs, we showed through simulation that a design involving pooling would be 1.07 times more expensive than a design without pooling. Thus, while a significant amount of information exists within the data from pooled DNA, our analysis does not support genotyping pooled DNA as a means to efficiently identify common variants contributing small effects to phenotypes of interest. While our conclusions were based on the specific technology and study design we employed, the approach presented here will be useful for evaluating the utility of other or future genome-wide genotyping platforms in pooled DNA studies.

Association of linear growth impairment in pediatric Crohn's disease and a known height locus: a pilot study.

The etiology of growth impairment in Crohn's disease (CD) has been inadequately explained by nutritional, hormonal, and/or disease-related factors, suggesting that genetics may be an additional contributor. The aim of this cross-sectional study was to investigate genetic variants associated with linear growth in pediatric-onset CD. We genotyped 951 subjects (317 CD patient-parent trios) for 64 polymorphisms within 14 CD-susceptibility and 23 stature-associated loci. Patient height-for-age Z-score < -1.64 was used to dichotomize probands into growth-impaired and nongrowth-impaired groups. The transmission disequilibrium test (TDT) was used to study association to growth impairment. There was a significant association between growth impairment in CD (height-for-age Z-score < -1.64) and a stature-related polymorphism in the dymeclin gene DYM (rs8099594) (OR = 3.2, CI [1.57-6.51], p = 0.0007). In addition, there was nominal over-transmission of two CD-susceptibility alleles, 10q21.1 intergenic region (rs10761659) and ATG16L1 (rs10210302), in growth-impaired CD children (OR = 2.36, CI [1.26-4.41] p = 0.0056 and OR = 2.45, CI [1.22-4.95] p = 0.0094, respectively). Our data indicate that genetic influences due to stature-associated and possibly CD risk alleles may predispose CD patients to alterations in linear growth. This is the first report of a link between a stature-associated locus and growth impairment in CD.

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

UNLABELLED: 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 x 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 x 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.

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.

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.

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).

Use of weighted reference panels based on empirical estimates of ancestry for capturing untyped variation.

Many association methods use a subset of genotyped single nucleotide polymorphisms (SNPs) to capture or infer genotypes at other untyped SNPs. We and others previously showed that tag SNPs selected to capture common variation using data from The International HapMap Consortium (Nature 437:1299-1320, 2005), The International HapMap Consortium (Nature 449:851-861, 2007) could also capture variation in populations of similar ancestry to HapMap reference populations (de Bakker et al. in Nat Genet 38:1298-1303, 2006; González-Neira et al. in Genome Res 16:323-330, 2006; Montpetit et al. in PLoS Genet 2:282-290, 2006; Mueller et al. in Am J Hum Genet 76:387-398, 2005). To capture variation in admixed populations or populations less similar to HapMap panels, a "cosmopolitan approach," in which all samples from HapMap are used as a single reference panel, was proposed. Here we refine this suggestion and show that use of a "weighted reference panel," constructed based on empirical estimates of ancestry in the target population (relative to available reference panels), is more efficient than the cosmopolitan approach. Weighted reference panels capture, on average, only slightly fewer common variants (minor allele frequency > 5%) than the cosmopolitan approach (mean r (2) = 0.977 vs. 0.989, 94.5% variation captured vs. 96.8% at r (2) > 0.8), across the five populations of the Multiethnic Cohort, but entail approximately 25% fewer tag SNPs per panel (average 538 vs. 718). These results extend a recent study in two Indian populations (Pemberton et al. in Ann Hum Genet 72:535-546, 2008). Weighted reference panels are potentially useful for both the selection of tag SNPs in diverse populations and perhaps in the design of reference panels for imputation of untyped genotypes in genome-wide association studies in admixed populations.

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.

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 x 10(-7) to 8 x 10(-22)). Together, these 12 loci account for approximately 2% of the population variation in height. Individuals with < or =8 height-increasing alleles and > or =16 height-increasing alleles differ in height by approximately 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.

Comprehensive evaluation of ESR2 variation and ovarian cancer risk.

Studies indicate that estrogen receptor beta, encoded by the ESR2 gene on chromosome 14q, may play a role in ovarian carcinogenesis. Using the genetic structure data generated by the Breast and Prostate Cohort Consortium (BPC3), we have comprehensively characterized the role of haplotype diversity in ESR2 and risk of ovarian cancer. Five haplotypes with a frequency of > or =5% were observed in White subjects and five haplotype tagging SNPs (htSNP) were selected to capture the locus diversity with a minimum R(h)(2) of 0.81. The htSNPs were genotyped in 574 White controls, 417 White invasive ovarian cancer cases, and 123 White borderline ovarian cancer cases from case-control studies carried out in Los Angeles County from 1994 through 2004. No statistically significant association was observed between the five htSNPs and related haplotypes and risk of ovarian cancer overall. Haplotype D was associated with a nonstatistically significant increased risk of invasive ovarian cancer overall (odds ratio, 1.38; 95% confidence interval, 0.93-2.02; P = 0.11) relative to the most common haplotype and a statistically significant increased risk of invasive clear cell ovarian cancer (odds ratio, 3.88; 95% confidence interval, 1.28-11.73; P = 0.016). Haplotype D was also reported by the BPC3 to be associated with increased risk of breast cancer. This haplotype warrants further investigation to rule out any effect with invasive ovarian cancer risk.

Common genetic variation in eight genes of the GH/IGF1 axis does not contribute to adult height variation.

Stature (adult height) is one of the most heritable human traits, yet few genes, if any, have been convincingly associated with adult height variation in the general population. Here, we selected 150 tag SNPs from eight candidate genes in the growth hormone (GH)/insulin-like growth factor-1 (IGF1) axis (GHR, GHRH, GHRHR, IGF1, IGFALS, IGFBP3, JAK2, STAT5B), and genotyped them in approximately 2,200 individuals ascertained for short or tall stature. Nominally significant tag SNPs were then tested in three additional replication cohorts, including a family-based panel to rule out spurious associations owing to population stratification. Across the four height cohorts (N = 6,075 individuals), we did not observe any consistent associations between stature and common variants (> or =5% minor allele frequency) in these eight genes, including a common deletion of the growth hormone receptor gene exon 3. Tests of epistatic interactions between these genes did not yield any results beyond those expected by chance. Although we have not tested all genes in the GH/IGF1 axis, our results indicate that common variation in these GH/IGF1 axis genes is not a major determinant of stature, and suggest that if common variation contributes to adult height variation in the general population, the variants are in other, possibly unanticipated genes.

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.

Determination of sequence variation and haplotype structure for the gonadotropin-releasing hormone (GnRH) and GnRH receptor genes: investigation of role in pubertal timing.

Because GnRH and its receptor (GnRHR) are pivotal regulators of the reproductive endocrine axis and mutations in GNRHR lead to hypogonadotropic hypogonadism, we investigated whether genetic variation in GNRHR or GNRH1 affects pubertal timing in the general population. To screen for missense mutations in these genes that might affect pubertal timing, we resequenced the coding regions of these genes in 48 probands with late but otherwise normal pubertal development. No missense variants were found in either gene, except for a previously identified single nucleotide polymorphism (SNP) in GNRH1 that was not associated with late pubertal development. To search for common variants that might affect pubertal timing, we took a haplotype-based association approach. To identify common haplotypes in these genes, we genotyped 41 SNPs in DNA from commercially available European-derived multigenerational pedigrees and participants in a multiethnic cohort (MEC). Two blocks of strong linkage disequilibrium were identified that spanned GNRHR and one was identified spanning GNRH1; within each block, more than 80% of chromosomes carried one of a few common haplotypes. A set of haplotype-tagging SNPs that mark these common haplotypes in all five ethnic groups within the MEC were defined and used to perform association studies among 125 trios (probands with late pubertal development and their parents) and 506 women from the MEC who had early (menarche < 11 yr of age, n = 216) or late (menarche > or = 15 yr of age, n = 290) pubertal development. Three SNPs in GNRHR showed modest association with late pubertal development in the trios; among the 506 women, a different SNP was associated with late menarche, and one rare haplotype was associated with early age of menarche. All of the observed associations were relatively modest and only nominally statistically significant; replication is needed to determine their validity. We conclude that genetic variation in GNRH1 and GNRHR is not likely to be a substantial modulator of pubertal timing in the general population.