Multiple Genes Influence BMI on Chromosome 7q31–34: The NHLBI Family Heart Study

The National Heart, Lung, and Blood Institute Family Heart Study (FHS) genome‐wide linkage scan identified a region of chromosome 7q31–34 with a lod score of 4.9 for BMI at D7S1804 (131.9 Mb). We report the results of linkage and association to BMI in this region for two independent FHS samples. The first sample includes 225 FHS pedigrees with evidence of linkage to 7q31–34, using 1,132 single‐nucleotide polymorphisms (SNPs) and 7 microsatellites. The second represents a case–control sample (318 cases; BMI >25 and 325 controls; BMI <25) derived from unrelated FHS participants who were not part of the genome scan. The latter set was genotyped for 606 SNPs, including 37 SNPs with prior evidence for association in the linked families. Although variance components linkage analysis using only SNPs generated a peak lod score that coincided with the original linkage scan at 131.9 Mb, a conditional linkage analysis showed evidence of a second quantitative trait locus (QTL) near 143 cM influencing BMI. Three SNPs (rs161339, rs12673281, and rs1993068) located near the three genes pleiotrophin (PTN), diacylglycerol (DAG) kinase iota (DGKι), and cholinergic receptor, muscarinic 2 (CHRM2) demonstrated significant association in both linked families (P = 0.0005, 0.002, and 0.03, respectively) and the case–control sample (P = 0.01, 0.0003, and 0.03, respectively), regardless of the genetic model tested. These findings suggest that several genes may be associated with BMI in the 7q31–34 region.     

Single-nucleotide polymorphisms and haplotype LD analysis of the 29-kb IGF2 region on chromosome 11p15.5 in the Korean population

OBJECTIVE: We investigated sequence variations of the 29-kb insulin-like growth factor 2 (IGF2) region in human chromosome region 11p15.5 in the Korean population. This region consists of IGF2, insulin-like growth factor 2 antisense (IGF2AS), and the insulin gene, all important candidate genes for various diseases, including cancer, obesity, diabetes, and coronary disease. While single nucleotide polymorphisms (SNPs) have been identified for this region and used in association studies, ethnic differences in genetic variation at this site have not been addressed. To date, SNPs for the entire 29-kb region in the Korean population have not been reported. METHODS: We surveyed a population of 108 Koreans for SNPs in the 29-kb IGF2 region. RESULTS: We identified 62 SNPs, consisting of 6 SNPs in the promoter region, 17 in the untranslated region, 19 in introns, and 20 in the intergenic region. We also analyzed linkage disequilibrium (LD) patterns and haplotypes using 36 high-frequency (> 5%)SNPs and found a well-defined LD block spanning about 13 kb that includes 8 kb of the IGF2AS gene, with two hot-spot regions flanking the LD block. CONCLUSION: These SNPs may be useful as genetic markers in disease association studies in the Korean population.     

Strong signature of natural selection within an FHIT intron implicated in prostate cancer risk

Previously, a candidate gene linkage approach on brother pairs affected with prostate cancer identified a locus of prostate cancer susceptibility at D3S1234 within the fragile histidine triad gene (FHIT), a tumor suppressor that induces apoptosis. Subsequent association tests on 16 SNPs spanning approximately 381 kb surrounding D3S1234 in Americans of European descent revealed significant evidence of association for a single SNP within intron 5 of FHIT. In the current study, re-sequencing and genotyping within a 28.5 kb region surrounding this SNP further delineated the association with prostate cancer risk to a 15 kb region. Multiple SNPs in sequences under evolutionary constraint within intron 5 of FHIT defined several related haplotypes with an increased risk of prostate cancer in European-Americans. Strong associations were detected for a risk haplotype defined by SNPs 138543, 142413, and 152494 in all cases (Pearson's chi(2) = 12.34, df 1, P = 0.00045) and for the homozygous risk haplotype defined by SNPs 144716, 142413, and 148444 in cases that shared 2 alleles identical by descent with their affected brothers (Pearson's chi(2) = 11.50, df 1, P = 0.00070). In addition to highly conserved sequences encompassing SNPs 148444 and 152413, population studies revealed strong signatures of natural selection for a 1 kb window covering the SNP 144716 in two human populations, the European American (pi = 0.0072, Tajima's D = 3.31, 14 SNPs) and the Japanese (pi = 0.0049, Fay & Wu's H = 8.05, 14 SNPs), as well as in chimpanzees (Fay & Wu's H = 8.62, 12 SNPs). These results strongly support the involvement of the FHIT intronic region in an increased risk of prostate cancer.     

Investigation of a genome wide association signal for obesity: synthetic association and haplotype analyses at the melanocortin 4 receptor gene locus

Background

Independent genome-wide association studies (GWAS) showed an obesogenic effect of two single nucleotide polymorphisms (SNP; rs12970134 and rs17782313) more than 150 kb downstream of the melanocortin 4 receptor gene (MC4R). It is unclear if the SNPs directly influence MC4R function or expression, or if the SNPs are on a haplotype that predisposes to obesity or includes functionally relevant genetic variation (synthetic association). As both exist, functionally relevant mutations and polymorphisms in the MC4R coding region and a robust association downstream of the gene, MC4R is an ideal model to explore synthetic association.

Methodology/Principal Findings

We analyzed a genomic region (364.9 kb) encompassing the MC4R in GWAS data of 424 obesity trios (extremely obese child/adolescent and both parents). SNP rs12970134 showed the lowest p-value (p = 0.004; relative risk for the obesity effect allele: 1.37); conditional analyses on this SNP revealed that 7 of 78 analyzed SNPs provided independent signals (p≤0.05). These 8 SNPs were used to derive two-marker haplotypes. The three best (according to p-value) haplotype combinations were chosen for confirmation in 363 independent obesity trios. The confirmed obesity effect haplotype includes SNPs 3′ and 5′ of the MC4R. Including MC4R coding variants in a joint model had almost no impact on the effect size estimators expected under synthetic association.

Conclusions/Significance

A haplotype reaching from a region 5′ of the MC4R to a region at least 150 kb from the 3′ end of the gene showed a stronger association to obesity than single SNPs. Synthetic association analyses revealed that MC4R coding variants had almost no impact on the association signal. Carriers of the haplotype should be enriched for relevant mutations outside the MC4R coding region and could thus be used for re-sequencing approaches. Our data also underscore the problems underlying the identification of relevant mutations depicted by GWAS derived SNPs.     

Power and sample size calculations for genetic case/control studies using gene-centric SNP maps: application to human chromosomes 6, 21, and 22 in three populations

Power and sample size calculations are critical parts of any research design for genetic association. We present a method that utilizes haplotype frequency information and average marker-marker linkage disequilibrium on SNPs typed in and around all genes on a chromosome. The test statistic used is the classic likelihood ratio test applied to haplotypes in case/control populations. Haplotype frequencies are computed through specification of genetic model parameters. Power is determined by computation of the test's non-centrality parameter. Power per gene is computed as a weighted average of the power assuming each haplotype is associated with the trait. We apply our method to genotype data from dense SNP maps across three entire chromosomes (6, 21, and 22) for three different human populations (African-American, Caucasian, Chinese), three different models of disease (additive, dominant, and multiplicative) and two trait allele frequencies (rare, common). We perform a regression analysis using these factors, average marker-marker disequilibrium, and the haplotype diversity across the gene region to determine which factors most significantly affect average power for a gene in our data. Also, as a 'proof of principle' calculation, we perform power and sample size calculations for all genes within 100 kb of the PSORS1 locus (chromosome 6) for a previously published association study of psoriasis. Results of our regression analysis indicate that four highly significant factors that determine average power to detect association are: disease model, average marker-marker disequilibrium, haplotype diversity, and the trait allele frequency. These findings may have important implications for the design of well-powered candidate gene association studies. Our power and sample size calculations for the PSORS1 gene appear consistent with published findings, namely that there is substantial power (>0.99) for most genes within 100 kb of the PSORS1 locus at the 0.01 significance level.     

Genome-Wide Association Studies Using Haplotypes and Individual SNPs in Simmental Cattle

Recent advances in high-throughput genotyping technologies have provided the opportunity to map genes using associations between complex traits and markers. Genome-wide association studies (GWAS) based on either a single marker or haplotype have identified genetic variants and underlying genetic mechanisms of quantitative traits. Prompted by the achievements of studies examining economic traits in cattle and to verify the consistency of these two methods using real data, the current study was conducted to construct the haplotype structure in the bovine genome and to detect relevant genes genuinely affecting a carcass trait and a meat quality trait. Using the Illumina BovineHD BeadChip, 942 young bulls with genotyping data were introduced as a reference population to identify the genes in the beef cattle genome significantly associated with foreshank weight and triglyceride levels. In total, 92,553 haplotype blocks were detected in the genome. The regions of high linkage disequilibrium extended up to approximately 200 kb, and the size of haplotype blocks ranged from 22 bp to 199,266 bp. Additionally, the individual SNP analysis and the haplotype-based analysis detected similar regions and common SNPs for these two representative traits. A total of 12 and 7 SNPs in the bovine genome were significantly associated with foreshank weight and triglyceride levels, respectively. By comparison, 4 and 5 haplotype blocks containing the majority of significant SNPs were strongly associated with foreshank weight and triglyceride levels, respectively. In addition, 36 SNPs with high linkage disequilibrium were detected in the GNAQ gene, a potential hotspot that may play a crucial role for regulating carcass trait components.     

A family-based association study identified CYP17 as a candidate gene for obesity susceptibility in Caucasians

The cytochrome P450c17α gene (CYP17) encodes a key biosynthesis enzyme of estrogen, which is critical in regulating adipogenesis and adipocyte development in humans. We therefore hypothesized that CYP17 is a candidate gene for predicting obesity. In order to test this hypothesis, we performed a family-based association test to investigate the relationship between the CYP17 gene and obesity phenotypes in a large sample comprising 1873 subjects from 405 Caucasian nuclear families of European origin recruited by the Osteoporosis Research Center of Creighton University, USA. Both single SNPs and haplotypes were tested for associations with obesity-related phenotypes, including body mass index (BMI) and fat mass. We identified three SNPs to be significantly associated with BMI, including rs3740397, rs6163, and rs619824. We further characterized the linkage disequilibrium structure for CYP17 and found that the whole CYP17 gene was located in a single-linkage disequilibrium block. This block was observed to be significantly associated with BMI. A major haplotype in this block was significantly associated with both BMI and fat mass. In conclusion, we suggest that the CYP17 gene has an effect on obesity in the Caucasian population. Further independent studies will be needed to confirm our findings.     

Haplotypes in the promoter region of the ADIPOQ gene are associated with increased diabetes risk in a German Caucasian population.

Adiponectin, which is encoded by the ADIPOQ gene, has been shown to modulate insulin sensitivity and glucose homeostasis. Plasma adiponectin levels are decreased in type 2 diabetes and obesity. Genetic variations within the ADIPOQ gene are associated with decreased adiponectin hormone levels. To analyze specific single-nucleotide polymorphisms (SNPs) and their association with T2D, 365 German subjects with T2D and 323 control subjects were screened. Three common SNPs - +45T>G in exon 2, and 2 promoter variants SNPs -11391G>A and -11377C>G - were analyzed. We found that the variant allele of SNP -11391G>A was significantly more frequent in the diabetic patient group than in the control group (p=0.003). Carrying the haplotype of SNP -11391A and SNP -11377C was associated with a 1.50-fold (p=0.03) increase in diabetes risk. The combination of the A-C haplotype and the G-C haplotype was associated with significantly elevated diabetes risk (OR=2.82 (95% CI: 1.35-5.91), p=0.006) after correction for BMI and age. Our observations suggest that diploid combinations of haplotype in the adiponectin gene promoter region contribute to the genetic risk of T2D in individuals from a German Caucasian population.     

A common haplotype in NAPEPLD is associated with severe obesity in a Norwegian population-based cohort (the HUNT study)

Obesity has a strong genetic etiology involving numerous identified metabolic pathways and others not yet examined. We investigated the association between severe obesity and genetic variation in selected candidate genes, including three drug-related genes: cannabinoid receptor 1 (CNR1), N-acyl phosphatidylethanolamine phospholipase D (NAPEPLD), and gastric lipase (LIPF); and three genes related to inflammation: nicotinamide phosphoribosyltransferase, six-transmembrane epithelial antigen of the prostate 4 (STEAP4) and interleukin 18 (IL-18). Subjects were 1,632 individuals with severe obesity (BMI ≥ 35 kg/m2) and 3,379 controls (BMI 20-24.9 kg/m2) that took part in a Norwegian population based cohort study. Tagging single-nucleotide polymorphisms (SNPs) of the coding region of these genes were analyzed. SNP-haplotypes for each gene were constructed in order to analyze allelic, genotypic, and haplotypic association to obesity. A single SNPs rs17605251 in NAPEPLD was nominally associated with BMI ≥ 35 kg/m2 (P = 0.035). A common haplotype in NAPEPLD was associated with BMI ≥ 35 kg/m2 after correction for multiple testing. The allele frequency was 56.8% in cases and 60.3% in controls, giving an odds ratio (OR) of 0.87 (95% confidence interval (CI) 0.79, 0.95; P = 0.0016). Homozygosity for this haplotype was protective against obesity (OR 0.79 (CI 0.70-0.91); P = 0.00059). The SNP rs7913071 in LIPF was associated with obesity, but the association lost statistical significance after correction for multiple testing. The CNR1, IL-18, STEAP4, and nicotinamide phosphoribosyltransferase genes were not associated with obesity. In conclusion a common haplotype in NAPEPLD, an enzyme involved in endocannabinoid synthesis, was protective against obesity.     

Further evidence for the role of ENPP1 in obesity: association with morbid obesity in Finns

The aim of this study was to investigate a series of single-nucleotide polymorphisms (SNPs) in the genes MC2R, MC3R, MC4R, MC5R, POMC, and ENPP1 for association with obesity. Twenty-five SNPs (2-7 SNPs/gene) were genotyped in 246 Finns with extreme obesity (BMI > or = 40 kg/m2) and in 481 lean subjects (BMI 20-25 kg/m2). Of the obese subjects, 23% had concomitant type 2 diabetes. SNPs and SNP haplotypes were tested for association with obesity and type 2 diabetes. Allele frequencies differed between obese and lean subjects for two SNPs in the ENPP1 gene, rs1800949 (P = 0.006) and rs943003 (P = 0.0009). These SNPs are part of a haplotype (rs1800949 C-rs943003 A), which was observed more frequently in lean subjects compared to obese subjects (P = 0.0007). Weaker associations were detected between the SNPs rs1541276 in the MC5R, rs1926065 in the MC3R genes and obesity (P = 0.04 and P = 0.03, respectively), and between SNPs rs2236700 in the MC5R, rs2118404 in the POMC, rs943003 in the ENPP1 genes and type 2 diabetes (P = 0.03, P = 0.02 and P = 0.02, respectively); these associations did not, however, remain significant after correction for multiple testing. In conclusion, a previously unexplored ENPP1 haplotype composed of SNPs rs1800949 and rs943003 showed suggestive evidence for association with adult-onset morbid obesity in Finns. In this study, we did not find association between the frequently studied ENPP1 K121Q variant, nor SNPs in the MCR or POMC genes and obesity or type 2 diabetes.     

Single-nucleotide polymorphisms in the TNF gene are associated with obesity-related phenotypes in vervet monkeys

Tumor necrosis factor (TNF) promoter single-nucleotide polymorphisms (SNPs) have been extensively characterized in humans, with numerous reports of associations with obesity-related phenotypes as well an array of infectious, immune-mediated, and inflammatory disease phenotypes. Controlling for the multitude of environmental risk factors in human studies has been a major confounder of efforts to elucidate the role and relative contribution of TNF promoter SNPs. As part of an ongoing initiative to further genetically and phenotypically characterize the St Kitts-origin vervet monkey (Chlorocebus aethiops ssp.) as an animal model of human obesity, we have conducted association analyses between TNF SNPs and previously defined obesity-related phenotypes in 265 pedigreed vervets. We report eight SNPs (-809G, -756A, -352C, -322A, +1285T, +2133T, +2362A, +2405), all contained within the same haplotype block and comprising a single haplotype, to be significantly associated with BMI, waist circumference, total plasma cholesterol (P < 0.05), and high-density lipoprotein-cholesterol (HDL-C) (P < 0.01). This study provides additional validation of the St Kitts-origin vervet model of obesity by demonstrating genetic associations analogous to that shown in humans.     

Polymorphisms in the trace amine receptor 4 (TRAR4) gene on chromosome 6q23.2 are associated with susceptibility to schizophrenia

Several linkage studies across multiple population groups provide convergent support for a susceptibility locus for schizophrenia--and, more recently, for bipolar disorder--on chromosome 6q13-q26. We genotyped 192 European-ancestry and African American (AA) pedigrees with schizophrenia from samples that previously showed linkage evidence to 6q13-q26, focusing on the MOXD1-STX7-TRARs gene cluster at 6q23.2, which contains a number of prime candidate genes for schizophrenia. Thirty-one screening single-nucleotide polymorphisms (SNPs) were selected, providing a minimum coverage of at least 1 SNP/20 kb. The association observed with rs4305745 (P=.0014) within the TRAR4 (trace amine receptor 4) gene remained significant after correction for multiple testing. Evidence for association was proportionally stronger in the smaller AA sample. We performed database searches and sequenced genomic DNA in a 30-proband subsample to obtain a high-density map of 23 SNPs spanning 21.6 kb of this gene. Single-SNP analyses and also haplotype analyses revealed that rs4305745 and/or two other polymorphisms in perfect linkage disequilibrium (LD) with rs4305745 appear to be the most likely variants underlying the association of the TRAR4 region with schizophrenia. Comparative genomic analyses further revealed that rs4305745 and/or the associated polymorphisms in complete LD with rs4305745 could potentially affect gene expression. Moreover, RT-PCR studies of various human tissues, including brain, confirm that TRAR4 is preferentially expressed in those brain regions that have been implicated in the pathophysiology of schizophrenia. These data provide strong preliminary evidence that TRAR4 is a candidate gene for schizophrenia; replication is currently being attempted in additional clinical samples.     

Multiple SNPs in intron 41 of thyroglobulin gene are associated with autoimmune thyroid disease in the Japanese population

BACKGROUND: The etiology of the autoimmune thyroid diseases (AITDs), Graves' disease (GD) and Hashimoto's thyroiditis (HT), is largely unknown. However, genetic susceptibility is believed to play a major role. Two whole genome scans from Japan and from the US identified a locus on chromosome 8q24 that showed evidence for linkage with AITD and HT. Recent studies have demonstrated an association between thyroglobulin (Tg) polymorphisms and AITD in Caucasians, suggesting that Tg is a susceptibility gene on 8q24. OBJECTIVES: The objective of the study was to refine Tg association with AITD, by analyzing a panel of 25 SNPs across an extended 260 kb region of the Tg. METHODS: We studied 458 Japanese AITD patients (287 GD and 171 HT patients) and 221 matched Japanese control subjects in association studies. Case-control association studies were performed using 25 Tg single nucleotide polymorphisms (SNPs) chosen from a database of the Single Nucleotide Polymorphism Database (dbSNP). Haplotype analysis was undertaken using the computer program SNPAlyze version 7.0. PRINCIPAL FINDINGS AND CONCLUSIONS: In total, 5 SNPs revealed association with GD (P<0.05), with the strongest SNP associations at rs2256366 (P?=?0.002) and rs2687836 (P?=?0.0077), both located in intron 41 of the Tg gene. Because of the strong LD between these two strongest associated variants, we performed the haplotype analysis, and identified a major protective haplotype for GD (P?=?0.001). These results suggested that the Tg gene is involved in susceptibility for GD and AITD in the Japanese.     

Genetic Origins of Lactase Persistence and the Spread of Pastoralism in Africa

In humans, the ability to digest lactose, the sugar in milk, declines after weaning because of decreasing levels of the enzyme lactase-phlorizin hydrolase, encoded by LCT. However, some individuals maintain high enzyme amounts and are able to digest lactose into adulthood (i.e., they have the lactase-persistence [LP] trait). It is thought that selection has played a major role in maintaining this genetically determined phenotypic trait in different human populations that practice pastoralism. To identify variants associated with the LP trait and to study its evolutionary history in Africa, we sequenced MCM6 introns 9 and 13 and ∼2 kb of the LCT promoter region in 819 individuals from 63 African populations and in 154 non-Africans from nine populations. We also genotyped four microsatellites in an ∼198 kb region in a subset of 252 individuals to reconstruct the origin and spread of LP-associated variants in Africa. Additionally, we examined the association between LP and genetic variability at candidate regulatory regions in 513 individuals from eastern Africa. Our analyses confirmed the association between the LP trait and three common variants in intron 13 (C-14010, G-13907, and G-13915). Furthermore, we identified two additional LP-associated SNPs in intron 13 and the promoter region (G-12962 and T-956, respectively). Using neutrality tests based on the allele frequency spectrum and long-range linkage disequilibrium, we detected strong signatures of recent positive selection in eastern African populations and the Fulani from central Africa. In addition, haplotype analysis supported an eastern African origin of the C-14010 LP-associated mutation in southern Africa.     

Association of the ABCA1 gene polymorphisms with type 2 DM in a Japanese population

To examine the association of the ATP-binding cassette transporter 1 (ABCA1) gene with type 2 diabetes (DM), we studied genetic polymorphisms of the ABCA1 gene including its linkage disequilibrium (LD) and haplotype analyses using a Japanese population. A sample set (DM:72, IGT:75, and NGT:227) was genotyped with 34 SNPs distributed from the promoter region to the last exon of the ABCA1 gene. LD between SNPs was assessed in pairwise manner. Among 13 LD blocks constructed, an LD block at the 5'-region showed a significant difference in the haplotype distribution between the study groups (NGT vs. IGT + DM: overall p = 0.0180; NGT vs. DM: 0.0001). Fisher's exact probability test (NGT vs. DM) showed a significant association of the haplotype 2 of the LD block (p = 0.0001), with an odds ratio (OR) of 2.53 (95%CI:1.62-4.12). Diplotype analysis also showed a significant association of the diplotypes with the haplotype 2 (OR:2.59, 95%CI:1.48-4.54, p = 0.0013).     

Haplotype association analysis of AGT variants with hypertension-related traits: the HyperGEN study

OBJECTIVE: Function of the renin-angiotensin system is important to human hypertension, but its genetic etiology remains elusive. We set out to examine a hypothesis that multiple genetic variants in the system act together in blood pressure regulation, via intermediate phenotypes such as blood pressure reactivity. METHODS: A sample of 531 hypertensive cases and 417 controls was selected from the HyperGEN study. Hypertension-related traits including blood pressure responses to challenges to math test, handgrip and postural change (mathBP, gripBP, and postBP), and body mass index (BMI) were analyzed for association with 10 single nucleotide polymorphisms (SNPs) in the angiotensinogen (AGT) gene. Single-marker and haplotype analyses were performed to examine the effects of both individual and multiple variants. Multiple-trait profiling was used to assess interaction of latent intermediate factors with susceptible haplotypes. RESULTS: In Blacks, two SNPs in exon 5 and 3'UTR showed significant association with gripBP, and two promoter SNPs were strongly associated with postBP. In Whites, only borderline association was found for 2 promoter SNPs with mathBP. Haplotype analyses in Blacks confirmed association with gripBP, and detected significant association of a haplotype to BMI (p=0.029). With the interactions modeled, haplotype associations found in Blacks remain significant, while significant associations to BMI (p=0.009) and gripSBP emerged in Whites. CONCLUSION: Genetic variants in regulatory regions of AGT showed strong association with blood pressure reactivity. Interaction of promoter and genic SNPs in AGT revealed collective action of multiple variants on blood pressure reactivity and BMI both in Blacks and in Whites, possibly following different pathways.