Multi-ancestry Genome-Wide Association Study of Early Childhood Caries

Early childhood caries (ECC) is the most common non-communicable childhood disease. It is an important health problem with known environmental and social/behavioral influences that lacks evidence for specific associated genetic risk loci. To address this knowledge gap, we conducted a genome-wide association study of ECC in a multi-ancestry population of U.S. preschool-age children (n=6,103) participating in a community-based epidemiologic study of early childhood oral health. Calibrated examiners used ICDAS criteria to measure ECC with the primary trait using the dmfs index with decay classified as macroscopic enamel loss (ICDAS ≥3). We estimated heritability, concordance rates, and conducted genome-wide association analyses to estimate overall genetic effects; the effects stratified by sex, household water fluoride, and dietary sugar; and leveraged the combined gene/gene-environment effects using the 2-degree-of-freedom (2df) joint test. The common genetic variants explained 24% of the phenotypic variance (heritability) of the primary ECC trait and the concordance rate was higher with a higher degree of relatedness. We identified 21 novel non-overlapping genome-wide significant loci for ECC. Two loci, namely RP11-856F16 . 2 (rs74606067) and SLC41A3 (rs71327750) showed evidence of association with dental caries in external cohorts, namely the GLIDE consortium adult cohort (n=∼487,000) and the GLIDE pediatric cohort (n=19,000), respectively. The gene-based tests identified TAAR6 as a genome-wide significant gene. Implicated genes have relevant biological functions including roles in tooth development and taste. These novel associations expand the genomics knowledge base for this common childhood disease and underscore the importance of accounting for sex and pertinent environmental exposures in genetic investigations of oral health.

Large trans-ethnic meta-analysis identifies AKR1C4 as a novel gene associated with age at menarche

STUDY QUESTION

Does the expansion of genome-wide association studies (GWAS) to a broader range of ancestries improve the ability to identify and generalise variants associated with age at menarche (AAM) in European populations to a wider range of world populations?

SUMMARY ANSWER

By including women with diverse and predominantly non-European ancestry in a large-scale meta-analysis of AAM with half of the women being of African ancestry, we identified a new locus associated with AAM in African-ancestry participants, and generalised loci from GWAS of European ancestry individuals.

WHAT IS KNOWN ALREADY

AAM is a highly polygenic puberty trait associated with various diseases later in life. Both AAM and diseases associated with puberty timing vary by race or ethnicity. The majority of GWAS of AAM have been performed in European ancestry women.

STUDY DESIGN, SIZE, DURATION

We analysed a total of 38 546 women who did not have predominantly European ancestry backgrounds: 25 149 women from seven studies from the ReproGen Consortium and 13 397 women from the UK Biobank. In addition, we used an independent sample of 5148 African-ancestry women from the Southern Community Cohort Study (SCCS) for replication.

PARTICIPANTS/MATERIALS, SETTING, METHODS

Each AAM GWAS was performed by study and ancestry or ethnic group using linear regression models adjusted for birth year and study-specific covariates. ReproGen and UK Biobank results were meta-analysed using an inverse variance-weighted average method. A trans-ethnic meta-analysis was also carried out to assess heterogeneity due to different ancestry.

MAIN RESULTS AND THE ROLE OF CHANCE

We observed consistent direction and effect sizes between our meta-analysis and the largest GWAS conducted in European or Asian ancestry women. We validated four AAM loci (1p31, 6q16, 6q22 and 9q31) with common genetic variants at P < 5 × 10-7. We detected one new association (10p15) at P < 5 × 10-8 with a low-frequency genetic variant lying in AKR1C4, which was replicated in an independent sample. This gene belongs to a family of enzymes that regulate the metabolism of steroid hormones and have been implicated in the pathophysiology of uterine diseases. The genetic variant in the new locus is more frequent in African-ancestry participants, and has a very low frequency in Asian or European-ancestry individuals.

LARGE SCALE DATA

N/A.

LIMITATIONS, REASONS FOR CAUTION

Extreme AAM (<9 years or >18 years) were excluded from analysis. Women may not fully recall their AAM as most of the studies were conducted many years later. Further studies in women with diverse and predominantly non-European ancestry are needed to confirm and extend these findings, but the availability of such replication samples is limited.

WIDER IMPLICATIONS OF THE FINDINGS

Expanding association studies to a broader range of ancestries or ethnicities may improve the identification of new genetic variants associated with complex diseases or traits and the generalisation of variants from European-ancestry studies to a wider range of world populations.

STUDY FUNDING/COMPETING INTEREST(S)

Funding was provided by CHARGE Consortium grant R01HL105756-07: Gene Discovery For CVD and Aging Phenotypes and by the NIH grant U24AG051129 awarded by the National Institute on Aging (NIA). The authors have no conflict of interest to declare.

Pharmacogenomics study of thiazide diuretics and QT interval in multi-ethnic populations: the cohorts for heart and aging research in genomic epidemiology

Thiazide diuretics, commonly used antihypertensives, may cause QT interval (QT) prolongation, a risk factor for highly fatal and difficult to predict ventricular arrhythmias. We examined whether common single-nucleotide polymorphisms (SNPs) modified the association between thiazide use and QT or its component parts (QRS interval, JT interval) by performing ancestry-specific, trans-ethnic and cross-phenotype genome-wide analyses of European (66%), African American (15%) and Hispanic (19%) populations (N=78 199), leveraging longitudinal data, incorporating corrected standard errors to account for underestimation of interaction estimate variances and evaluating evidence for pathway enrichment. Although no loci achieved genome-wide significance (P<5 × 10-8), we found suggestive evidence (P<5 × 10-6) for SNPs modifying the thiazide-QT association at 22 loci, including ion transport loci (for example, NELL1, KCNQ3). The biologic plausibility of our suggestive results and simulations demonstrating modest power to detect interaction effects at genome-wide significant levels indicate that larger studies and innovative statistical methods are warranted in future efforts evaluating thiazide-SNP interactions.

BMI loci and longitudinal BMI from adolescence to young adulthood in an ethnically diverse cohort

Objective

The association of obesity susceptibility variants with change in body mass index (BMI) across the life course is not well understood.

Subjects

In ancestry stratified models of 5,962 European American (EA), 2,080 African American (AA), and 1,582 Hispanic American (HA) individuals from the National Longitudinal Study of Adolescent to Adult Health (Add Health), we examined associations between 34 obesity SNPs with per year change in BMI, measured by the slope from a growth-curve analysis of two or more BMI measurements between adolescence and young adulthood. For SNPs nominally associated with BMI change (p<0.05), we interrogated age differences within data collection Wave and time differences between age categories that overlapped between Waves.

Results

We found SNPs in/near FTO, MC4R, MTCH2, TFAP2B, SEC16B, and TMEM18 were significantly associated (p<0.0015 ≈ 0.05/34) with BMI change in EA and the ancestry-combined meta-analysis. Rs9939609 in FTO met genome-wide significance at p<5e-08 in the EA and ancestry combined analysis, respectively [Beta(se)=0.025(0.004);Beta(se)=0.021(0.003)]. No SNPs were significant after Bonferroni correction in AA or HA, although 5 SNPs in AA and 4 SNPs in HA were nominally significant (p<0.05). In EA and the ancestry-combined meta-analysis, rs3817334 near MTCH2 showed larger effects in younger respondents, while rs987237 near TFAP2B, showed larger effects in older respondents across all Waves. Differences in effect estimates across time for MTCH2 and TFAP2B are suggestive of either era or cohort effects.

Conclusion

The observed association between variants in/near FTO, MC4R, MTCH2, TFAP2B, SEC16B, and TMEM18 with change in BMI from adolescence to young adulthood suggest that the genetic effect of BMI loci varies over time in a complex manner, highlighting the importance of investigating loci influencing obesity risk across the life course.

The interaction between physical activity and obesity gene variants in association with BMI: Does the obesogenic environment matter?

Little is known about how obesity susceptibility single nucleotide polymorphisms (SNPs) interact with moderate to vigorous physical activity (MVPA) in relation to BMI during adolescence, once obesogenic neighborhood factors are accounted for. In race stratified models, including European (EA; N=4977), African (AA; N=1726), and Hispanic Americans (HA; N=1270) from the National Longitudinal Study of Adolescent to Adult Health (1996; ages 12-21), we assessed the evidence for a SNPxMVPA interaction with BMI-for-age Z score, once accounting for obesogenic neighborhood factors including physical activity amenities, transportation and recreation infrastructure, poverty and crime. Eight SNPxMVPA interactions with suggestive significance (p<0.10; three in each EA, and AA, two in HA) were observed showing attenuation on BMI-for-age Z score in adolescents with ≥5 versus <5 bouts/week MVPA, except for rs10146997 (near NRXN3). Findings were robust to the inclusion of neighborhood-level variables as covariates. These findings suggest that any attenuation from MVPA on a genetic susceptibility to obesity during adolescence is likely not operating through obesogenic neighborhood factors.

Chronic Periodontitis Genome-wide Association Study in the Hispanic Community Health Study / Study of Latinos

Chronic periodontitis (CP) has a genetic component, particularly its severe forms. Evidence from genome-wide association studies (GWASs) has highlighted several potential novel loci. Here, the authors report the first GWAS of CP among a large community-based sample of Hispanics/Latinos. The authors interrogated a quantitative trait of CP (mean interproximal clinical attachment level determined by full-mouth periodontal examinations) among 10,935 adult participants (mean age: 45 y, range: 18 to 76 y) from the Hispanic Community Health Study / Study of Latinos. Genotyping was done with a custom Illumina Omni2.5M array, and imputation to approximately 20 million single-nucleotide polymorphisms was based on the 1000 Genomes Project phase 1 reference panel. Analyses were based on linear mixed models adjusting for sex, age, study design features, ancestry, and kinship and employed a conventional P < 5 × 10^-8 statistical significance threshold. The authors identified a genome-wide significant association signal in the 1q42.2 locus ( TSNAX-DISC1 noncoding RNA, lead single-nucleotide polymorphism: rs149133391, minor allele [C] frequency = 0.01, P = 7.9 × 10^-9) and 4 more loci with suggestive evidence of association ( P < 5 × 10^-6): 1q22 (rs13373934), 5p15.33 (rs186066047), 6p22.3 (rs10456847), and 11p15.1 (rs75715012). We tested these loci for replication in independent samples of European-American ( n = 4,402) and African-American ( n = 908) participants of the Atherosclerosis Risk in Communities study. There was no replication among the European Americans; however, the TSNAX-DISC1 locus replicated in the African-American sample (rs149133391, minor allele frequency = 0.02, P = 9.1 × 10^-3), while the 1q22 locus was directionally concordant and nominally significant (rs13373934, P = 4.0 × 10^-2). This discovery GWAS of interproximal clinical attachment level-a measure of lifetime periodontal tissue destruction-was conducted in a large, community-based sample of Hispanic/Latinos. It identified a genome-wide significant locus that was independently replicated in an African-American population. Identifying this genetic marker offers direction for interrogation in subsequent genomic and experimental studies of CP.

Influence of SNP*SNP interaction on BMI in European American adolescents: findings from the National Longitudinal Study of Adolescent Health

BACKGROUND

Adolescent obesity is predictive of future weight gain, obesity and adult onset severe obesity (body mass index [BMI] ≥40 kg m(-2) ). Despite successful efforts to identify Single Nucleotide Polymorphisms (SNPs) influencing BMI, <5% of the 40-80% heritability of the phenotype has been explained. Identification of gene-gene (G-G) interactions between known variants can help explain this hidden heritability as well as identify potential biological mechanisms affecting weight gain during this critical developmental period.

OBJECTIVE

We have recently shown distinct genetic effects on BMI across the life course, and thus it is important to examine the evidence for epistasis in adolescence.

METHODS

In adolescent participants of European descent from wave II of the National Longitudinal Study of Adolescent Health (Add Health, n = 5072, ages 12-21, 52.5% female), we tested 34 established BMI-related SNPs for G-G interaction effects on BMI z-score. We used mixed-effects regression, assuming multiplicative interaction models adjusting for age, sex and geographic region, with random effects for family and school.

RESULTS

For 28 G-G interactions that were nominally significant (P < 0.05), we attempted to replicate our results in an adolescent sample from the Childhood European American Cohort from Philadelphia. In the replication study, one interaction (PRKD1-FTO) was significant after correction for multiple testing.

CONCLUSIONS

Our results are suggestive of epistatic effects on BMI during adolescence and point to potentially interactive effects between genes in biological pathways important in obesity.

Lipoprotein lipase variants interact with polyunsaturated fatty acids for obesity traits in women: replication in two populations

BACKGROUND AND AIMS

Lipoprotein lipase (LPL) is a candidate gene for obesity based on its role in triglyceride hydrolysis and the partitioning of fatty acids towards storage or oxidation. Whether dietary fatty acids modify LPL associated obesity risk is unknown.

METHODS AND RESULTS

We examined five single nucleotide polymorphisms (SNPs) (rs320, rs2083637, rs17411031, rs13702, rs2197089) for potential interaction with dietary fatty acids for obesity traits in 1171 participants (333 men and 838 women, aged 45-75 y) of the Boston Puerto Rican Health Study (BPRHS). In women, SNP rs320 interacted with dietary polyunsaturated fatty acids (PUFA) for body mass index (BMI) (P = 0.002) and waist circumference (WC) (P = 0.001) respectively. Higher intake of PUFA was associated with lower BMI and WC in homozygotes of the major allele (TT) (P = 0.01 and 0.005) but not in minor allele carriers (TG and GG). These interactions were replicated in an independent population, African American women of the Atherosclerosis Risk in Communities (ARIC) study (n = 1334).

CONCLUSION

Dietary PUFA modulated the association of LPL rs320 with obesity traits in two independent populations. These interactions may be relevant to the dietary management of obesity, particularly in women.

Chronic periodontitis genome-wide association studies: gene-centric and gene set enrichment analyses

Recent genome-wide association studies (GWAS) of chronic periodontitis (CP) offer rich data sources for the investigation of candidate genes, functional elements, and pathways. We used GWAS data of CP (n = 4,504) and periodontal pathogen colonization (n = 1,020) from a cohort of adult Americans of European descent participating in the Atherosclerosis Risk in Communities study and employed a MAGENTA approach (i.e., meta-analysis gene set enrichment of variant associations) to obtain gene-centric and gene set association results corrected for gene size, number of single-nucleotide polymorphisms, and local linkage disequilibrium characteristics based on the human genome build 18 (National Center for Biotechnology Information build 36). We used the Gene Ontology, Ingenuity, KEGG, Panther, Reactome, and Biocarta databases for gene set enrichment analyses. Six genes showed evidence of statistically significant association: 4 with severe CP (NIN, p = 1.6 × 10(-7); ABHD12B, p = 3.6 × 10(-7); WHAMM, p = 1.7 × 10(-6); AP3B2, p = 2.2 × 10(-6)) and 2 with high periodontal pathogen colonization (red complex-KCNK1, p = 3.4 × 10(-7); Porphyromonas gingivalis-DAB2IP, p = 1.0 × 10(-6)). Top-ranked genes for moderate CP were HGD (p = 1.4 × 10(-5)), ZNF675 (p = 1.5 × 10(-5)), TNFRSF10C (p = 2.0 × 10(-5)), and EMR1 (p = 2.0 × 10(-5)). Loci containing NIN, EMR1, KCNK1, and DAB2IP had showed suggestive evidence of association in the earlier single-nucleotide polymorphism-based analyses, whereas WHAMM and AP2B2 emerged as novel candidates. The top gene sets included severe CP ("endoplasmic reticulum membrane," "cytochrome P450," "microsome," and "oxidation reduction") and moderate CP ("regulation of gene expression," "zinc ion binding," "BMP signaling pathway," and "ruffle"). Gene-centric analyses offer a promising avenue for efficient interrogation of large-scale GWAS data. These results highlight genes in previously identified loci and new candidate genes and pathways possibly associated with CP, which will need to be validated via replication and mechanistic studies.

Genetic risk score and adiposity interact to influence triglyceride levels in a cohort of Filipino women

Background/Objectives:

Individually, genetic variants only moderately influence cardiometabolic (CM) traits, such as lipid and inflammatory markers. In this study we generated genetic risk scores from a combination of previously reported variants influencing CM traits, and used these scores to explore how adiposity levels could mediate genetic contributions to CM traits.

Subjects/Methods:

Participants included 1649 women from the 2005 Cebu Longitudinal Health and Nutrition Survey. Three genetic risk scores were constructed for C-reactive protein (CRP), high-density lipoprotein cholesterol (HDL-C) and triglycerides (TGs). We used linear regression models to assess the association between each genetic risk score and its related trait. We also tested for interactions between each score and measures of adiposity.

Results:

Each genetic risk score explained a greater proportion of variance in trait levels than any individual genetic variant. We found an interaction between the TG genetic risk score (2.29–14.34 risk alleles) and waist circumference (WC) (P_interaction=1.66 × 10⁻²). Based on model predictions, for individuals with a higher TG genetic risk score (75th percentile=12), having an elevated WC (⩾80 cm) increased TG levels from 1.32 to 1.71 mmol l⁻¹. However, for individuals with a lower score (25th percentile=7), having an elevated WC did not significantly change TG levels.

Conclusions:

The TG genetic risk score interacted with adiposity to synergistically influence TG levels. For individuals with a genetic predisposition to elevated TG levels, our results suggest that reducing adiposity could possibly prevent further increases in TG levels and thereby lessen the likelihood of adverse health outcomes such as cardiovascular disease.

Moderate to vigorous physical activity interactions with genetic variants and body mass index in a large US ethnically diverse cohort

BACKGROUND

Little is known about the interaction between genetic and behavioural factors during lifecycle risk periods for obesity and how associations vary across race/ethnicity.

OBJECTIVE

The objective of this study was to examine joint associations of adiposity-related single-nucleotide polymorphisms (SNPs) and moderate to vigorous physical activity (MVPA) with body mass index (BMI) in a diverse adolescent cohort.

METHODS

Using data from the National Longitudinal Study of Adolescent Health (n = 8113: Wave II 1996; ages 12-21, Wave III; ages 18-27), we assessed interactions of 41 well-established SNPs and MVPA with BMI-for-age Z-scores in European Americans (EA; n = 5077), African-Americans (AA; n = 1736) and Hispanic Americans (HA; n = 1300).

RESULTS

Of 97 assessed, we found nominally significant SNP-MVPA interactions on BMI-for-age Z-score in EA at GNPDA2 and FTO and in HA at LZTR2/SEC16B. In EA, the estimated effect of the FTO risk allele on BMI-for-age Z-score was lower (β = -0.13; 95% confidence interval [CI]: 0.08, 0.18) in individuals with ≥5 vs. <5 (β = 0.24; CI: 0.16, 0.32) bouts of MVPA per week (P for interaction 0.02). Race/ethnicity-pooled meta-analysis showed nominally significant interactions for SNPs at TFAP2B, POC5 and LYPLAL1.

CONCLUSIONS

High MVPA may attenuate underlying genetic risk for obesity during adolescence, a high-risk period for adult obesity.

Screen time behaviours may interact with obesity genes, independent of physical activity, to influence adolescent BMI in an ethnically diverse cohort

BACKGROUND

There has been little investigation of gene-by-environment interactions related to sedentary behaviour, a risk factor for obesity defined as leisure screen time (ST; i.e. television, video and computer games).

OBJECTIVE

To test the hypothesis that limiting ST use attenuates the genetic predisposition to increased body mass index (BMI), independent of physical activity.

DESIGN

Using 7642 wave II participants of the National Longitudinal Study of Adolescent Health, (Add Health; mean = 16.4 years, 52.6% female), we assessed the interaction of ST (h week(-1) ) and 41 established obesity single nucleotide polymorphisms (SNPs) with age- and sex-specific BMI Z-scores in 4788 European-American (EA), 1612 African-American (AA) and 1242 Hispanic American (HA) adolescents.

RESULTS

Nominally significant SNP*ST interaction were found for FLJ35779 in EA, GNPDA2 in AA and none in HA (EA: beta [SE] = 0.016[0.007]), AA: beta [SE] = 0.016[0.011]) per 7 h week(-1) ST and one risk allele in relation to BMI Z-score.

CONCLUSIONS

While for two established BMI loci, we find evidence that high levels of ST exacerbate the influence of obesity susceptibility variants on body mass; overall, we do not find strong evidence for interactions between the majority of established obesity loci. However, future studies with larger sample sizes, or that may build on our current study and the growing published literature, are clearly warranted.

Epidemiology and genetic determinants of progressive deterioration of glycaemia in American Indians: the Strong Heart Family Study

AIMS/HYPOTHESIS

Type 2 diabetes is a chronic, heterogeneous disease and a major risk factor for cardiovascular diseases. The underlying mechanisms leading to progression to type 2 diabetes are not fully understood and genetic tools may help to identify important pathways of glycaemic deterioration.

METHODS

Using prospective data on American Indians from the Strong Heart Family Study, we identified 373 individuals defined as progressors (diabetes incident cases), 566 individuals with transitory impaired fasting glucose (IFG) and 1,011 controls (normal fasting glycaemia at all visits). We estimated the heritability (h(2)) of the traits and the evidence for association with 16 known variants identified in type 2 diabetes genome-wide association studies.

RESULTS

We noted high h(2) for diabetes progression (h(2) = 0.65 ± 0.16, p = 2.7 × 10(-6)) but little contribution of genetic factors to transitory IFG (h(2) = 0.09 ± 0.10, p = 0.19) for models adjusted for multiple risk factors. At least three variants (in WFS1, TSPAN8 and THADA) were nominally associated with diabetes progression in age- and sex-adjusted analyses with estimates showing the same direction of effects as reported in the discovery European ancestry studies.

CONCLUSIONS/INTERPRETATION

Our findings do not exclude these loci for diabetes susceptibility in American Indians and suggest phenotypic heterogeneity of the IFG trait, which may have implications for genetic studies when diagnosis is based on a single time-point measure.

Generalization of adiposity genetic loci to US Hispanic women

BACKGROUND:

Obesity is a public health concern. Yet the identification of adiposity-related genetic variants among United States (US) Hispanics, which is the largest US minority group, remains largely unknown.

OBJECTIVE:

To interrogate an a priori list of 47 (32 overall body mass and 15 central adiposity) index single-nucleotide polymorphisms (SNPs) previously studied in individuals of European descent among 3494 US Hispanic women in the Women's Health Initiative SNP Health Association Resource (WHI SHARe).

DESIGN:

Cross-sectional analysis of measured body mass index (BMI), waist circumference (WC) and waist-to-hip ratio (WHR) were inverse normally transformed after adjusting for age, smoking, center and global ancestry. WC and WHR models were also adjusted for BMI. Genotyping was performed using the Affymetrix 6.0 array. In the absence of an a priori selected SNP, a proxy was selected (r²⩾0.8 in CEU).

RESULTS:

Six BMI loci (TMEM18, NUDT3/HMGA1, FAIM2, FTO, MC4R and KCTD15) and two WC/WHR loci (VEGFA and ITPR2-SSPN) were nominally significant (P<0.05) at the index or proxy SNP in the corresponding BMI and WC/WHR models. To account for distinct linkage disequilibrium patterns in Hispanics and further assess generalization of genetic effects at each locus, we interrogated the evidence for association at the 47 surrounding loci within 1 Mb region of the index or proxy SNP. Three additional BMI loci (FANCL, TFAP2B and ETV5) and five WC/WHR loci (DNM3-PIGC, GRB14, ADAMTS9, LY86 and MSRA) displayed Bonferroni-corrected significant associations with BMI and WC/WHR. Conditional analyses of each index SNP (or its proxy) and the most significant SNP within the 1 Mb region supported the possible presence of index-independent signals at each of these eight loci as well as at KCTD15.

CONCLUSION:

This study provides evidence for the generalization of nine BMI and seven central adiposity loci in Hispanic women. This study expands the current knowledge of common adiposity-related genetic loci to Hispanic women.

Common genetic variation near the connexin-43 gene is associated with resting heart rate in African Americans: a genome-wide association study of 13,372 participants

BACKGROUND

Genome-wide association studies have identified several genetic loci associated with variation in resting heart rate in European and Asian populations. No study has evaluated genetic variants associated with heart rate in African Americans.

OBJECTIVE

To identify novel genetic variants associated with resting heart rate in African Americans.

METHODS

Ten cohort studies participating in the Candidate-gene Association Resource and Continental Origins and Genetic Epidemiology Network consortia performed genome-wide genotyping of single nucleotide polymorphisms (SNPs) and imputed 2,954,965 SNPs using HapMap YRI and CEU panels in 13,372 participants of African ancestry. Each study measured the RR interval (ms) from 10-second resting 12-lead electrocardiograms and estimated RR-SNP associations using covariate-adjusted linear regression. Random-effects meta-analysis was used to combine cohort-specific measures of association and identify genome-wide significant loci (P≤2.5×10(-8)).

RESULTS

Fourteen SNPs on chromosome 6q22 exceeded the genome-wide significance threshold. The most significant association was for rs9320841 (+13 ms per minor allele; P = 4.98×10(-15)). This SNP was approximately 350 kb downstream of GJA1, a locus previously identified as harboring SNPs associated with heart rate in Europeans. Adjustment for rs9320841 also attenuated the association between the remaining 13 SNPs in this region and heart rate. In addition, SNPs in MYH6, which have been identified in European genome-wide association study, were associated with similar changes in the resting heart rate as this population of African Americans.

CONCLUSIONS

An intergenic region downstream of GJA1 (the gene encoding connexin 43, the major protein of the human myocardial gap junction) and an intragenic region within MYH6 are associated with variation in resting heart rate in African Americans as well as in populations of European and Asian origin.

Estimation of genetic effects on BMI during adolescence in an ethnically diverse cohort: The National Longitudinal Study of Adolescent Health

OBJECTIVE

The contribution of genetic variants to body mass index (BMI) during adolescence across multiethnic samples is largely unknown. We selected genetic loci associated with BMI or obesity in European-descent samples and examined them in a multiethnic adolescent sample.

DESIGN AND SAMPLE

In 5103 European American (EA), 1748 African American (AfA), 1304 Hispanic American (HA) and 439 Asian American (AsA) participants of the National Longitudinal Study of Adolescent Health (Add Health; ages 12-21 years, 47.5% male), we assessed the association between 41 established obesity-related single-nucleotide polymorphisms (SNPs) with BMI using additive genetic models, stratified by race/ethnicity, and in a pooled meta-analysis sample. We also compared the magnitude of effect for BMI-SNP associations in EA and AfA adolescents to comparable effect estimates from 11 861 EA and AfA adults in the Atherosclerosis Risk in Communities study (ages 45-64 years, 43.2% male).

RESULTS

Thirty-five of 41 BMI-SNP associations were directionally consistent with published studies in European populations, 18 achieved nominal significance (P<0.05; effect sizes from 0.19 to 0.71 kg m(-2) increase in BMI per effect allele), while 4 (FTO, TMEM18, TFAP2B, MC4R) remained significant after Bonferroni correction (P<0.0015). Of 41 BMI-SNP associations in AfA, HA and AsA adolescents, nine, three and five, respectively, were directionally consistent and nominally significant. In the pooled meta-analysis, 36 of 41 effect estimates were directionally consistent and 21 of 36 were nominally significant. In EA adolescents, BMI effect estimates were larger (P<0.05) for variants near TMEM18, PTER and MC4R and smaller for variants near MTIF3 and NRXN3 compared with EA adults.

CONCLUSION

Our findings suggest that obesity susceptibility loci may have a comparatively stronger role during adolescence than during adulthood, with variation across race/ethnic subpopulation.

The combined influence of genetic factors and sedentary activity on body mass changes from adolescence to young adulthood: the National Longitudinal Adolescent Health Study

BACKGROUND

an increase in sedentary activities is likely a major contributor to the rise in obesity over the last three decades. Little research has examined interactions between genetic variants and sedentary activity on obesity phenotypes. High levels of sedentary activity during adolescence may interact with genetic factors to influence body mass changes between adolescence and young adulthood, a high risk period for weight gain.

METHODS

in the National Longitudinal Study of Adolescent Health, siblings and twin pairs (16.5 ± 1.7 years) were followed into young adulthood (22.4 ± 1.8 years). Self-reported screen time (TV, video, and computer use in h/week) and body mass index (kg/m(2) ), calculated from measured height and weight at adolescence and at young adulthood, were available for 3795 participants. We employed a variance component approach to estimate the interaction between genotype and screen time for body mass changes. Additive genotype-by-screen time interactions were assessed using likelihood-ratio tests. Models were adjusted for race, age, sex, and age-by-sex interaction.

RESULTS

the genetic variation in body mass changes was significantly larger in individuals with low ( δ(G) = 27.59 ± 1.58) compared with high (δ(G) = 18.76 ± 2.59) levels of screen time (p < 0.003) during adolescence.

CONCLUSIONS

Our findings demonstrate that sedentary activities during adolescence may interact with genetic factors to influence body mass changes between adolescence and young adulthood. Accounting for obesity-related behaviours may improve current understanding of the genetic variation in body mass changes.

Sex-specific interaction between APOE genotype and carbohydrate intake affects plasma HDL-C levels: the Strong Heart Family Study

Low plasma levels of high-density lipoprotein cholesterol (HDL-C) are identified as a risk factor for cardiovascular disease (CVD). Sexual dimorphism, however, is widely reported in both HDL-C and CVD, with the underlying explanations of these sexual differences not fully understood. HDL-C is a complex trait influenced by both genes and dietary factors. Here we examine evidence for a sex-specific effect of APOE and the macronutrient carbohydrate on HDL-C, triglycerides (TG) and apoprotein A-1 (ApoA-1) in a sample of 326 male and 423 female participants of the Strong Heart Family Study (SHFS). Using general estimating equations in SAS to account for kinship correlations, stratifying by sex, and adjusting for age, body mass index (BMI) and SHS center, we examine the relationship between APOE genotype and carbohydrate intake on circulating levels of HDL-C, TG, and ApoA-1 through a series of carbohydrate-by-sex interactions and stratified analyses. APOE-by-carbohydrate intake shows significant sex-specific effects. All males had similar decreases in HDL-C levels associated with increased carbohydrate intake. However, only those females with APOE-4 alleles showed significantly lower HDL-C levels as their percent of carbohydrate intake increased, while no association was noted between carbohydrate intake and HDL-C in those females without an APOE-4 allele. These findings demonstrate the importance of understanding sex differences in gene-by-nutrient interaction when examining the complex architecture of HDL-C variation.

TCF7L2 single nucleotide polymorphisms, cardiovascular disease and all-cause mortality: the Atherosclerosis Risk in Communities (ARIC) study

AIMS/HYPOTHESIS

We hypothesised that TCF7L2 single nucleotide polymorphisms (SNPs) are associated with cardiovascular disease (CVD) and that the associations differ in diabetic and non-diabetic persons.

METHODS

Our analysis included black and white participants from the Atherosclerosis Risk in Communities study who were free of prevalent CVD at baseline and had been genotyped for rs7903146, rs12255372, rs7901695, rs11196205 and rs7895340 (n=13,369). Cox proportional hazard regression was used to estimate the associations between polymorphisms and incident events; logistic and linear regression were used for associations with baseline risk factor levels.

RESULTS

TCF7L2 SNPs were not significantly associated with incident coronary heart disease, ischaemic stroke, CVD, prevalent peripheral artery disease (PAD) or all-cause mortality in the full cohort or when stratified by race.

CONCLUSIONS/INTERPRETATION

In the whole cohort, TCF7L2 SNPs were not associated with incident CVD, all-cause mortality or prevalent PAD. This result suggests that the increased health risk associated with rs7903146 genotype is specific to diabetes.

Cyclooxygenase polymorphisms and risk of cardiovascular events: the Atherosclerosis Risk in Communities (ARIC) study

Cyclooxygenase-derived prostaglandins modulate cardiovascular disease risk. We genotyped 2212 Atherosclerosis Risk in Communities study participants (1,023 incident coronary heart disease (CHD) cases; 270 incident ischemic stroke cases; 919 non-cases) with available DNA for polymorphisms in PTGS1 and PTGS2. Using a case-cohort design, associations between genotype and CHD or stroke risk were evaluated using proportional hazards regression. In Caucasians, the reduced function PTGS1 -1006A variant allele was significantly more common among stroke cases compared to non-cases (18.2 versus 10.6%, P=0.027). In African Americans, the reduced function PTGS2 -765C variant allele was significantly more common in stroke cases (61.4 versus 49.4%, P=0.032). No significant relationships with CHD risk were observed. However, aspirin utilization appeared to modify the relationship between the PTGS2 G-765C polymorphism and CHD risk (interaction P=0.072). These findings suggest that genetic variation in PTGS1 and PTGS2 may be important risk factors for the development of cardiovascular disease events. Confirmation in independent populations is necessary.

Genotype-by-sex interaction in the aetiology of type 2 diabetes mellitus: support for sex-specific quantitative trait loci in Hypertension Genetic Epidemiology Network participants

AIMS/HYPOTHESIS

While there are sex-related differences in both the prevalence of type 2 diabetes mellitus and disease risk factors, there is only limited research on sex-specific influences on type 2 diabetes aetiology within the same study population. Thus, we assessed genotype-by-sex interaction using a liability threshold model in an attempt to localise sex-specific type 2 diabetes quantitative trait loci (QTLs).

SUBJECTS, MATERIALS AND METHODS

Hypertensive siblings and their offspring and/or parents in the Hypertension Genetic Epidemiology Network of the Family Blood Pressure Program were recruited from five field centres. The diabetic phenotype was adjusted for race, study centre, age and non-linear age effects. In total, 567 diabetic individuals were identified in 385 families. Variance component linkage analyses in the combined sample and stratified by sex and race were performed (SOLAR program) using race-specific marker allele frequencies derived from a random sample of participants at each centre.

RESULTS

We observed a QTL-specific genotype-by-sex interaction (p=0.009) on chromosome 17 at 31 cM, with females displaying a robust adjusted logarithm of odds (LOD) of 3.0 compared with 0.2 in males and 1.3 in the combined sample. Three additional regions demonstrating suggestive evidence for linkage were detected: chromosomes 2 and 5 in the female sample and chromosome 22 (adjusted LOD=1.9) in the combined sample.

CONCLUSIONS/INTERPRETATION

These findings suggest that multiple genes may regulate susceptibility to type 2 diabetes, demonstrating the importance of considering the interaction of genes and environment in the aetiology of common complex traits.

QTL-specific genotype-by-smoking interaction and burden of calcified coronary atherosclerosis: the NHLBI Family Heart Study

BACKGROUND

Calcified coronary plaque (CCP) is a complex trait influenced by both genes and environment, and plausibly an interaction between the two. Because the familial aggregation of CCP has been demonstrated and smoking is a significant, independent predictor of CCP, we assessed the evidence for genotype-by-smoking interaction and conducted linkage analysis of quantitative Agatston CCP scores in participants of the NHLBI Family Heart Study (FHS).

METHODS

During standardized clinical exams smoking habits were ascertained and CCP was quantified with cardiac computed tomography (CT). Among 4387 relationship pairs from 2128 Caucasian examinees variance component analysis was implemented in SOLAR to examine: (1) additive genotype-by-smoking status interaction using a variance component approach; (2) linkage analysis in the full sample and among smoking subsets defined by individual smoking exposure; (3) QTL-specific genotype-by-smoking interaction in the regions that appeared to differentiate between smoking strata.

RESULTS

The prevalence of CCP (and median Agatston score) was 75% (184.6) in men and 48% (51.0) in women. We detected four genome-wide significant logarithm of odds (LOD) scores in samples stratified by individual smoking exposure: chromosome 4 at 122cM (nearest marker D4S2297; robust adjusted LOD=3.1; q=0.053), chromosome 6 at 99cM (nearest marker D6S1056; robust adjusted LOD=3.3; q=0.053), chromosome 11 at 19cM (nearest marker D11S199; robust adjusted LOD=4.0; q=0.02) and chromosome 13 at 77cM (nearest marker D13S892; robust adjusted LOD=3.1; q=0.053). Additive and QTL-specific genotype-by-smoking interaction was detected on chromosomes 4, 6, 11 and 13; all P<0.05. Three of the four QTLs identified in this report have been previously linked to atherosclerosis and harbor interesting candidate genes.

CONCLUSIONS

These findings demonstrate the importance of considering complex interactions in the search for genes that influence the pathogenesis of CCP.

Efficient semiparametric estimation of haplotype-disease associations in case-cohort and nested case-control studies

Estimating the effects of haplotypes on the age of onset of a disease is an important step toward the discovery of genes that influence complex human diseases. A haplotype is a specific sequence of nucleotides on the same chromosome of an individual and can only be measured indirectly through the genotype. We consider cohort studies which collect genotype data on a subset of cohort members through case-cohort or nested case-control sampling. We formulate the effects of haplotypes and possibly time-varying environmental variables on the age of onset through a broad class of semiparametric regression models. We construct appropriate nonparametric likelihoods, which involve both finite- and infinite-dimensional parameters. The corresponding nonparametric maximum likelihood estimators are shown to be consistent, asymptotically normal, and asymptotically efficient. Consistent variance-covariance estimators are provided, and efficient and reliable numerical algorithms are developed. Simulation studies demonstrate that the asymptotic approximations are accurate in practical settings and that case-cohort and nested case-control designs are highly cost-effective. An application to a major cardiovascular study is provided.

Linkage analysis of LDL cholesterol in American Indian populations: the Strong Heart Family Study

Previous studies have demonstrated that low density lipoprotein cholesterol (LDL-C) concentration is influenced by both genes and environment. Although rare genetic variants associated with Mendelian causes of increased LDL-C are known, only one common genetic variant has been identified, the apolipoprotein E gene (APOE). In an attempt to localize quantitative trait loci (QTLs) influencing LDL-C, we conducted a genome-wide linkage scan of LDL-C in participants of the Strong Heart Family Study (SHFS). Nine hundred eighty men and women, age 18 years or older, in 32 extended families at three centers (in Arizona, Oklahoma, and North and South Dakota) were phenotyped for LDL-C concentration and other risk factors. Using a variance component approach and the program SOLAR, and after accounting for the effects of covariates, we detected a QTL influencing LDL-C on chromosome 19, nearest marker D19S888 at 19q13.41 [logarithm of odds (LOD) = 4.3] in the sample from the Dakotas. This region on chromosome 19 includes many possible candidate genes, including the APOE/C1/C4/C2 gene cluster. In follow-up association analyses, no significant evidence for an association was detected with the APOE*2 and APOE*4 alleles (P = 0.76 and P = 0.53, respectively). Suggestive evidence of linkage to LDL-C was detected on chromosomes 3q, 4q, 7p, 9q, 10p, 14q, and 17q. These linkage signals overlap positive findings for lipid-related traits and harbor plausible candidate genes for LDL-C.

Sex-specific findings from a genome-wide linkage analysis of human fatness in non-Hispanic whites and African Americans: the HyperGEN study

OBJECTIVE

To conduct a full genome search for genes potentially influencing two related phenotypes: body mass index (BMI, kg/m2) and percent body fat (PBF) from bioelectric impedance in men and women.

DESIGN

A total of 3383 participants, 1348 men and 2035 women; recruitment was initiated with hypertensive sibpairs and expanded to first-degree relatives in a multicenter study of hypertension genetics.

MEASUREMENTS

Genotypes for 387 highly polymorphic markers spaced to provide a 10 cM map (CHLC-8) were generated by the NHLBI Mammalian Genotyping Service (Marshfield, WI, USA). Quantitative trait loci for obesity phenotypes, BMI and PBF, were examined with a variance components method using SOLAR, adjusting for hypertensive status, ethnicity, center, age, age2, sex, and age2 x sex. As we detected a significant genotype-by-sex interaction in initial models and because of the importance of sex effects in the expression of these phenotypes, models thereafter were stratified by sex. No genotype-by-ethnicity interactions were found.

RESULTS

A QTL influencing PBF in women was detected on chromosome12q (12q24.3-12q24.32, maximum empirical LOD score=3.8); a QTL influencing this phenotype in men was found on chromosome 15q (15q25.3, maximum empirical LOD score=3.0). These QTLs were detected in African-American and white women (12q) and men (15q). QTLs influencing both BMI and PBF were found over a broad region on chromosome 3 in men. QTLs on chromosomes 3 and 12 were found in the combined sample of men and women, but with weaker significance.

CONCLUSION

The locations with highest LOD scores have been previously reported for obesity phenotypes, indicating that at least two genomic regions influence obesity-related traits. Furthermore, our results indicate the importance of considering context-dependent effects in the search for obesity QTLs.

Genotype-by-Sex Interaction in the Regulation of High-Density Lipoprotein: The Framingham Heart Study

Low levels of high-density lipoprotein (HDL) are widely documented as a risk factor for cardiovascular disease (CVD). Furthermore, there is marked sexual dimorphism in both HDL levels and the prevalence of CVD. However, the extent to which genetic factors contribute to such dimorphism has been largely unexplored. We examined the evidence for genotype-by-sex effects on HDL in a longitudinal sample of 1562 participants from 330 families in the Framingham Heart Study at three times points corresponding approximately to 1971-1974, 1980-1983, and 1988-1991. Using a variance component method, we conducted a genome scan of HDL at each time point in males and females, separately and combined, and tested for genotype-by-sex interaction at a quantitative trait locus (QTL) at each time point. Consistent findings were noted only for females on chromosome 2 near marker D2S1328, with adjusted LOD scores of 2.6, 2.2, and 2.1 across the three time points, respectively. In males suggestive linkage was detected on chromosome 16 near marker D16S3396 at the second time point and on chromosome 18 near marker D18S851 at the third time point (adjusted LOD = 2.2 and 2.4, respectively). Although the heritability of HDL is similar in males and females, sex appears to exert a substantial effect on the QTL-specific variance of HDL. When genotype-by-sex interactions exist and are not modeled, the power to detect linkage is reduced; thus our results may explain in part the paucity of significant linkage findings for HDL.

Evidence for joint action of genes on diabetes status and CVD risk factors in American Indians: the strong heart family study

OBJECTIVES

Previous research among American Indians of the strong heart family study (SHFS) has demonstrated significant heritabilities for CVD risk factors and implicated diabetes as an important predictor of several of the phenotypes. Moreover, we recently demonstrated that genetic effects on CVD risk factors differed in diabetic and nondiabetic individuals. In this paper, we investigated whether a significant genetic influence on diabetes status could be identified, and whether there is evidence for joint action of genes on diabetes status and related CVD risk factors.

METHODS AND RESULTS

Approximately 950 men and women, age 18 or older, in 32 extended families, were examined between 1997 and 1999. We estimated the effects of genes and environmental covariates on diabetes status using a threshold model and a maximum likelihood variance component approach. Diabetes status exhibited a residual heritability of 22% (h2=0.22). We also estimated the genetic and environmental correlations between diabetes susceptibility and eight risk factors for CVD. All eight CVD risk factors displayed significant genetic correlations with diabetes status (BMI (rhoG=0.55), fibrinogen (rhoG=0.40), HDL-C (rhoG=-0.37), ln triglycerides (rhoG=0.65), FAT (rhoG=0.38 ), PAI-1 (rhoG=0.67), SBP (rhoG=0.57), and WHR (rhoG=0.58)). Three of eight traits (HDL-C (rhoE=-0.32), ln triglycerides (rhoE=0.33), and fibrinogen (rhoE=0.20)) displayed significant environmental correlations with diabetes status.

CONCLUSIONS

These findings suggest that in the context of a high prevalence of diabetes, still unidentified diabetes genes may play an important role in influencing variation in CVD risk factors.

The utility of Bayesian model averaging for detecting known oligogenic effects

In this paper we examine the utility of the Bayesian model averaging approach proposed by Blangero et al. [Genet Epidemiol 17(Suppl. 1):S67-72, 1999] for the detection of known oligogenic effects on quantitative trait 1 (Q1) in the simulated Genetic Analysis Workshop 12 data. The results were consistent with the underlying genetic model: two major QTLs were identified on chromosomes 19 and 2, with heritabilities of 17% and 20%, respectively, located near the position of major genes 1 and 2.

Distribution of lod scores in oligogenic linkage analysis

In variance component oligogenic linkage analysis it can happen that the residual additive genetic variance bounds to zero when estimating the effect of the ith quantitative trait locus. Using quantitative trait Q1 from the Genetic Analysis Workshop 12 simulated general population data, we compare the observed lod scores from oligogenic linkage analysis with the empirical lod score distribution under a null model of no linkage. We find that zero residual additive genetic variance in the null model alters the usual distribution of the likelihood-ratio statistic.

The origins of the Irish travellers and the genetic structure of Ireland

Ireland's unique and well-documented history provides insight into the formation and origins of population subdivisions. Of particular interest, is the controversial ethnogenesis of an itinerant population of Ireland: the Travellers. The objectives of this study were: (1) to determine the genetic affinity of the Travellers to the general Irish population based on gene frequency data, subdivided by county, and (2) to explore the relationship between subpopulations of Ireland, given its turbulent history. The gene frequencies of standard genetic markers collected from populations residing in counties of Ireland and the Travellers were calculated and analysed using several multivariate methods. First, a relationship (R) matrix was used to ascertain the scaled variance covariance matrix of population similarity. Second, mean per locus heterozygosity (H) was regressed on distance of the region from the gene frequency centroid (r(ii)). The results of this study include: (1) the confirmation of Crawford's (1975, in Biosocial Interrelations in Population Adaptations, E. S. Watts et al. (eds), pp. 93-103) conclusions concerning the origins and genetic affinity of the Travellers; (2) based on several multivariate analyses, the major influence on population structure was unique historical events; and (3) Relethford and Crawford's (1995, American Journal of Physical Anthropology, 96, 25-38) hypothesis concerning the distinctiveness of the midland counties was verified by this study.

Gravidity and parity in postmenopausal American Indian women: the Strong Heart Study

The fertility of a large sample of American Indian women participating in the Strong Heart Study was examined to determine which factors are associated with variation in completed fertility among women in this population. The Strong Heart Study (SHS) is a study of cardiovascular disease (CVD) and its risk factors in American Indians living in Arizona, Oklahoma, and the Dakotas. Data were derived from a baseline examination between 1989 and 1992 of approximately 1,500 men and women, aged 45-74, from each of the 3 SHS centers. A personal interview elicited demographic information, family health history, and information on several life-style variables. A total of 1,955 ever-married, postmenopausal women were considered in these analyses. Women were considered to be postmenopausal if their menstrual cycles had stopped completely for at least 12 months, either because of natural or surgical processes. The average number of pregnancies (gravidity) for all women was 5.9, whereas the mean number of live births (parity) was 5.3. Women living in Arizona (5.6) and the Dakotas (5.8) had higher parity than those in Oklahoma (4.6). Furthermore, there was lower completed fertility in younger women: When American Indian women from all 3 centers were considered together, women born between 1910 and 1919 had a mean parity of 5.3, whereas women born between 1940 and 1949 had a mean parity of 4.0. Although previous research has suggested a relationship between parity and CVD risk factors, no linear associations between CVD risk factors and fertility were indicated in this population. We also examined the relationship of contraception, level of education, and income to fertility. While no significant relationship between contraception and the level of fertility was identified, there was a significant inverse linear relationship of both education and income with fertility. In summary, fertility rates in American Indian women are high, but appear to be decreasing in younger generations. Fertility is higher in those with less education and lower incomes.

Spatial variation of anthropometric traits in Ireland

To further elucidate the relationship between geography and genetics in Ireland, we considered variation in anthropometric traits of adult males by town using spatial autocorrelation methods. By describing and distinguishing significant patterns of anthropometric variation, we determined whether the anthropometric traits display a simple pattern of spatial variation, as predicted by the isolation by distance model, or other patterns of spatial variation. Several hypotheses were examined, including (1) whether there was spatial patterning of 20 anthropometric phenotypic distributions and 7 principal components of Irish males and (2) if there was, whether these phenotypic distributions could be explained by a simple isolation by distance model. The results of this study can be summarized by several key findings: (1) There is significant spatial patterning among towns, as detected in correlograms of 14 anthropometric traits and 2 principal component factor scores (values of Moran's I ranging from 0.7510 to -0.3616, p < or = 0.0071); (2) 4 spatial patterns were detected, including clinical patterns, long-distance differentiation, distance distinction, and regional patchiness. These results suggest several likely causes of the observed spatial patterns. First, in Ireland patterns of anthropometric variation could not be explained by a single spatial pattern (i.e., isolation by distance). Second, through an examination of the various combinations of statistical homogeneity or heterogeneity, spatial patterning or nonpatterning, and similarity or dissimilarity of spatial patterns, we conclude that several migrational events structured the genetic landscape of Ireland.