Publications

Telomeres play a central role in cellular aging, and shorter telomere length has been associated with age-related disorders including diabetes. However, a causal link between telomere shortening and diabetes risk has not been established. In a well-characterized longitudinal cohort of American Indians participating in the Strong Heart Family Study, we examined whether leukocyte telomere length (LTL) at baseline predicts incident diabetes independent of known diabetes risk factors. Among 2,328 participants free of diabetes at baseline, 292 subjects developed diabetes during an average 5.5 years of follow-up. Compared with subjects in the highest quartile (longest) of LTL, those in the lowest quartile (shortest) had an almost twofold increased risk of incident diabetes (hazard ratio [HR] 1.83 [95% CI 1.26-2.66]), whereas the risk for those in the second (HR 0.87 [95% CI 0.59-1.29]) and the third (HR 0.95 [95% CI 0.65-1.38]) quartiles was statistically nonsignificant. These findings suggest a nonlinear association between LTL and incident diabetes and indicate that LTL could serve as a predictive marker for diabetes development in American Indians, who suffer from disproportionately high rates of diabetes.

Inorganic arsenic is methylated in the body by arsenic (III) methyltransferase (AS3MT). Arsenic methylation is thought to play a role in arsenic-related epigenetic phenomena, including aberrant DNA and histone methylation. However, it is unclear whether the promoter of the AS3MT gene, which codes for AS3MT, is differentially methylated as a function of arsenic exposure. In this study, we evaluated AS3MT promoter methylation according to exposure, assessed by urinary arsenic excretion in a stratified random sample of 48 participants from the Strong Heart Study who had urine arsenic measured at baseline and DNA available from 1989 to 1991 and 1998-1999. For this study, all data are from the 1989-1991 visit. We measured AS3MT promoter methylation at its 48 CpG loci by bisulphite sequencing. We compared mean  % methylation at each CpG locus by arsenic exposure group using linear regression adjusted for study centre, age and sex. A hypomethylated region in the AS3MT promoter was associated with higher arsenic exposure. In vitro, arsenic induced AS3MT promoter hypomethylation, and it increased AS3MT expression in human peripheral blood mononuclear cells. These findings may suggest that arsenic exposure influences the epigenetic regulation of a major arsenic metabolism gene.

OBJECTIVE: Gene × diet interaction plays an important role in atherosclerosis, an inflammatory disorder. Leukotrienes are the most potent inflammatory mediators, and genetic variants encoding leukotriene genes have been implicated in atherosclerosis. This study tests nutrigenetic interaction of a previously defined leukotriene haplotype on carotid artery hypertrophy and atherosclerosis in American Indians.

METHODS: This study included 3402 American Indians participating in the Strong Heart Family Study (SHFS). Carotid artery measurements, including intima-media thickness (IMT), vascular mass, and plaque, were assessed using ultrasound. Eleven tagSNPs in the leukotriene A4 hydrolase (LTA4H) gene were genotyped in all subjects. Main haplotype effect and haplotype × diet interaction were examined by generalized estimating equation, adjusting for known risk factors.

RESULTS: There was no significant main effect of haplotype or diet on any of the carotid artery measures. However, a previously defined LTA4H haplotype, called HapE, significantly interacted with dietary intake of n-3 and n-6 fatty acids on both IMT (P(HapE × n3) = 0.018, P(HapE × n6) = 0.040) and vascular mass (P(HapE × n3) = 0.012, P(HapE × n6) = 0.018), but not plaque. The direction of this nutrigenetic interaction on IMT was consistent with that reported in a recent study of Caucasian twins.

CONCLUSION: Dietary intake of polyunsaturated fatty acids significantly modifies the effect of a leukotriene haplotype on carotid artery hypertrophy but not atherosclerosis in American Indians, independent of established cardiovascular risk factors. Replication of nutrigenetic interaction in two distinct ethnic groups suggests the robustness and generalizability of our findings to diverse populations.

BACKGROUND: C-reactive protein (CRP) is a biomarker of inflammation. Genome-wide association studies (GWAS) have identified single-nucleotide polymorphisms (SNPs) associated with CRP concentrations and inflammation-related traits such as cardiovascular disease, type 2 diabetes mellitus, and obesity. We aimed to replicate previous CRP-SNP associations, assess whether these associations generalize to additional race/ethnicity groups, and evaluate inflammation-related SNPs for a potentially pleiotropic association with CRP.

METHODS AND RESULTS: We selected and analyzed 16 CRP-associated and 250 inflammation-related GWAS SNPs among 40 473 African American, American Indian, Asian/Pacific Islander, European American, and Hispanic participants from 7 studies collaborating in the Population Architecture using Genomics and Epidemiology (PAGE) study. Fixed-effect meta-analyses combined study-specific race/ethnicity-stratified linear regression estimates to evaluate the association between each SNP and high-sensitivity CRP. Overall, 18 SNPs in 8 loci were significantly associated with CRP (Bonferroni-corrected P<3.1×10(-3) for replication, P<2.0×10(-4) for pleiotropy): Seven of these were specific to European Americans, while 9 additionally generalized to African Americans (1), Hispanics (5), or both (3); 1 SNP was seen only in African Americans and Hispanics. Two SNPs in the CELSR2/PSRC1/SORT1 locus showed a potentially novel association with CRP: rs599839 (P=2.0×10(-6)) and rs646776 (P=3.1×10(-5)).

CONCLUSIONS: We replicated 16 SNP-CRP associations, 10 of which generalized to African Americans and/or Hispanics. We also identified potentially novel pleiotropic associations with CRP for two SNPs previously associated with coronary artery disease and/or low-density lipoprotein-cholesterol. These findings demonstrate the benefit of evaluating genotype-phenotype associations in multiple race/ethnicity groups and looking for pleiotropic relationships among SNPs previously associated with related phenotypes.

BACKGROUND: Cigarette smoking negatively affects kidney function. Genetic variants in the nicotinic acetylcholine receptor (nAChR) genes have been associated with nicotine dependence, and are likely to influence renal function and related traits. Whereas each single variant may only exert a small effect, the joint contribution of multiple variants to the risk of disease could be substantial.

METHODS: Using a gene-family approach, we investigated the joint association of 61 tagging SNPs in seven genes encoding the nAChRs with kidney function in 3620 American Indians participating in the Strong Heart Family Study, independent of known risk factors. Kidney function was evaluated by estimated glomerular filtration rate, urinary albumin/creatinine ratio, albuminuria and chronic kidney disease. The joint impact of smoking-related variants was assessed using the weighted truncated product method.

RESULTS: Multiple SNPs showed marginal individual effect on renal function variability, and only a few survive multiple comparison correction. In contrast, a gene-family analysis considering the joint impact of all 61 SNPs reveals significant associations of the nAChR gene family with kidney function variables including estimated glomerular filtration rate, urinary albumin/creatinine ratio, and albuminuria (all Ps ≤ 0.0001) after adjusting for established risk factors including cigarette smoking.

CONCLUSION: Genetic variants in nAChR genes jointly contribute to renal function or kidney damage in American Indians. The effects of these genetic variants on kidney function or damage are independent of traditional risk factors including cigarette smoking per se.

Increased serum uric acid (SUA) or hyperuricemia, a risk factor for gout, renal and cardiovascular diseases, is caused by either increased production or decreased excretion of uric acid or a mix of both. The solute carrier protein 2 family, member 9 (SLC2A9) gene encodes a transporter that mediates urate flux across the renal proximal tubule. Genome-wide association studies have consistently shown the association of single-nucleotide polymorphisms in this gene with SUA in majority populations. American Indian participants of the Strong Heart Family Study, belonging to multigenerational families, have high prevalence of hyperuricemia. We conducted measured genotype analyses, based on variance components decomposition method and accounting for family relationships, to assess whether the association between SUA and SLC2A9 gene polymorphisms generalized to American Indians (n=3604) of this study. Seven polymorphisms were selected for genotyping based on their association with SUA levels in other populations. A strong association was found between SLC2A9 gene polymorphisms and SUA in all centers combined (P-values: 1.3 × 10(-31)-5.1 × 10(-23)) and also when stratified by recruitment center; P-values: 1.2 × 10(-14)-1.0 × 10(-5). These polymorphisms were also associated with the estimated glomerular filtration rate and serum creatinine but not albumin-creatinine ratio. In summary, the association of polymorphisms in the uric acid transporter gene with SUA levels extends to a new population of American Indians.

Increased immunoglobulin G (IgG) response to dietary antigens can be associated with gastrointestinal dysfunction and autoimmunity. The underlying processes contributing to these adverse reactions remain largely unknown, and it is likely that genetic factors play a role. Here, we estimate heritability and attempt to localize genetic factors influencing IgG antibody levels against food-derived antigens using an integrative genomics approach. IgG antibody levels were determined by ELISA in >1,300 Mexican Americans for the following food antigens: wheat gliadin; bovine casein; and two forms of bovine serum albumin (BSA-a and BSA-b). Pedigree-based variance components methods were used to estimate additive genetic heritability (h(2) ), perform genome-wide association analyses, and identify transcriptional signatures (based on 19,858 transcripts from peripheral blood lymphocytes). Heritability estimates were significant for all traits (0.15-0.53), and shared environment (based on shared residency among study participants) was significant for casein (0.09) and BSA-a (0.33). Genome-wide significant evidence of association was obtained only for antibody to gliadin (P = 8.57 × 10(-8) ), mapping to the human leukocyte antigen II region, with HLA-DRA and BTNL2 as the best candidate genes. Lack of association of known celiac disease risk alleles HLA-DQ2.5 and -DQ8 with antigliadin antibodies in the studied population suggests a separate genetic etiology. Significant transcriptional signatures were found for all IgG levels except BSA-b. These results demonstrate that individual genetic differences contribute to food antigen antibody measures in this population. Further investigations may elucidate the underlying immunological processes involved.

Cigarette smoke is a strong risk factor for obesity and cardiovascular disease. The effect of genetic variants involved in nicotine metabolism on obesity or body composition has not been well studied. Though many genetic variants have previously been associated with adiposity or body fat distribution, a single variant usually confers a minimal individual risk. The goal of this study is to evaluate the joint association of multiple variants involved in cigarette smoke or nicotine dependence with obesity-related phenotypes in American Indians. To achieve this goal, we genotyped 61 tagSNPs in seven genes encoding nicotine acetylcholine receptors (nAChRs) in 3,665 American Indians participating in the Strong Heart Family Study. Single SNP association with obesity-related traits was tested using family-based association, adjusting for traditional risk factors including smoking. Joint association of all SNPs in the seven nAChRs genes were examined by gene-family analysis based on weighted truncated product method (TPM). Multiple testing was controlled by false discovery rate (FDR). Results demonstrate that multiple SNPs showed weak individual association with one or more measures of obesity, but none survived correction for multiple testing. However, gene-family analysis revealed significant associations with waist circumference (p = 0.0001) and waist-to-hip ratio (p = 0.0001), but not body mass index (p = 0.20) and percent body fat (p = 0.29), indicating that genetic variants are jointly associated with abdominal, but not general, obesity among American Indians. The observed combined genetic effect is independent of cigarette smoking per se. In conclusion, multiple variants in the nAChR gene family are jointly associated with abdominal obesity in American Indians, independent of general obesity and cigarette smoking per se.

BACKGROUND: In rodents, hypothalamic brain-derived neurotrophic factor (BDNF) expression appears to be regulated by melanocortin-4 receptor (MC4R) activity. The impact of MC4R genetic variation on circulating BDNF in humans is unknown.

OBJECTIVE: The objective of this study is to compare BDNF concentrations of subjects with loss-of-function (LOF) and gain-of-function (GOF) MC4R variants with those of controls with common sequence MC4R.

METHODS: Circulating BDNF was measured in two cohorts with known MC4R sequence: 148 subjects of Pima Indian heritage ((mean±s.d.): age, 15.7±6.5 years; body mass index z-scores (BMI-Z), 1.63±1.03) and 69 subjects of Hispanic heritage (10.8±3.6 years; BMI-Z, 1.57±1.07). MC4R variants were characterized in vitro by cell surface expression, receptor binding and cyclic AMP response after agonist administration. BDNF single-nucleotide polymorphisms (SNPs) rs12291186, rs6265 and rs7124442 were also genotyped.

RESULTS: In the Pima cohort, no significant differences in serum BDNF was observed for 43 LOF subjects versus 65 LOF-matched controls (age, sex and BMI matched; P=0.29) or 20 GOF subjects versus 20 GOF-matched controls (P=0.40). Serum BDNF was significantly associated with genotype for BDNF rs12291186 (P=0.006) and rs6265 (P=0.009), but not rs7124442 (P=0.99); BDNF SNPs did not interact with MC4R status to predict serum BDNF. In the Hispanic cohort, plasma BDNF was not significantly different among 21 LOF subjects, 20 GOF subjects and 28 controls (P=0.79); plasma BDNF was not predicted by BDNF genotype or BDNF-x-MC4R genotype interaction.

CONCLUSIONS: Circulating BDNF concentrations were not significantly associated with MC4R functional status, suggesting that peripheral BDNF does not directly reflect hypothalamic BDNF secretion and/or that MC4R signaling is not a significant regulator of the bulk of BDNF expression in humans.

BACKGROUND: The association between human blood DNA global methylation and global hydroxymethylation has not been evaluated in population-based studies. No studies have evaluated environmental determinants of global DNA hydroxymethylation, including exposure to metals.

OBJECTIVE: We evaluated the association between global DNA methylation and global DNA hydroxymethylation in 48 Strong Heart Study participants for which selected metals had been measured in urine at baseline and DNA was available from 1989-1991 (visit 1) and 1998-1999 (visit 3).

METHODS: We measured the percentage of 5-methylcytosine (5-mC) and 5-hydroxymethylcytosine (5-hmC) in samples using capture and detection antibodies followed by colorimetric quantification. We explored the association of participant characteristics (i.e., age, adiposity, smoking, and metal exposure) with both global DNA methylation and global DNA hydroxymethylation.

RESULTS: The Spearman's correlation coefficient for 5-mC and 5-hmC levels was 0.32 (p = 0.03) at visit 1 and 0.54 (p < 0.001) at visit 3. Trends for both epigenetic modifications were consistent across potential determinants. In cross-sectional analyses, the odds ratios of methylated and hydroxymethylated DNA were 1.56 (95% CI: 0.95, 2.57) and 1.76 (95% CI: 1.07, 2.88), respectively, for the comparison of participants above and below the median percentage of dimethylarsinate. The corresponding odds ratios were 1.64 (95% CI: 1.02, 2.65) and 1.16 (95% CI: 0.70, 1.94), respectively, for the comparison of participants above and below the median cadmium level. Arsenic exposure and metabolism were consistently associated with both epigenetic markers in cross-sectional and prospective analyses. The positive correlation of 5-mC and 5-hmC levels was confirmed in an independent study population.

CONCLUSIONS: Our findings support that both epigenetic measures are related at the population level. The consistent trends in the associations between these two epigenetic modifications and the characteristics evaluated, especially arsenic exposure and metabolism, suggest the need for understanding which of the two measures is a better biomarker for environmental epigenetic effects in future large-scale epidemiologic studies.