Publications

The objective of this study was to investigate whether leukocyte telomere length (LTL) predicts the risk for cancer mortality among American Indians participating in the Strong Heart Study (1989-1991). Participants (aged 45-74 years) were followed annually until December 2015 to collect information on morbidity/mortality. LTL was measured by qPCR using genomic DNA isolated from peripheral blood. The association between LTL and risk for cancer mortality was examined using a multivariable Cox proportional hazard model, adjusting for age, gender, education, study site, smoking, alcohol use, physical activity, systolic blood pressure, fasting blood glucose, obesity, and low- and high-density lipoprotein. Of 1945 participants (mean age 56.10 ± 8.17 at baseline, 57% women) followed for an average 20.5 years, 220 died of cancer. Results showed that longer LTL at baseline significantly predicts an increased risk of cancer death among females (HR 1.57, 95% CI 1.08-2.30), but not males (HR 0.74, 95% CI 0.49-1.12) (p for interaction 0.009). Specifically, compared with the women with the longest LTL (fourth quartile), those in the third, second, and first quartiles showed 53%, 41%, and 44% reduced risk for cancer death, respectively. The findings highlight the importance of sex-specific analysis in future telomere research.

BACKGROUND: Diet plays a key role in development of diabetes, and there has been recent interest in better understanding the association of dairy food intake with diabetes.

OBJECTIVE: This study examined the associations of full-fat and low-fat dairy food intake with incident diabetes among American Indians-a population with a high burden of diabetes.

METHODS: The study included participants from the Strong Heart Family Study (SHFS), a family-based study of cardiovascular disease in American Indians, free of diabetes at baseline (2001-2003) (n = 1623). Participants were 14-86-y-old at baseline and 60.8% were female. Dairy food intake was assessed using a Block food frequency questionnaire. Incident diabetes was defined using American Diabetes Association criteria. Parametric survival models with a Weibull distribution were used to evaluate the associations of full-fat and low-fat dairy food intake with incident diabetes. Serving sizes were defined as 250 mL for milk and 42.5 g for cheese.

RESULTS: We identified 277 cases of diabetes during a mean follow-up of 11 y. Reported intake of dairy foods was low [median full-fat dairy food intake: 0.11 serving/1000 kcal; median low-fat dairy food intake: 0.03 serving/1000 kcal]. Participants who reported the highest full-fat dairy food intake had a lower risk of diabetes compared to those who reported the lowest full-fat food dairy intake [HR (95% CI): 0.79 (0.59, 1.06); P-trend = 0.03, comparing extreme tertiles, after adjustment for age, sex, site, physical activity, education, smoking, diet quality, and low-fat dairy food intake]. Low-fat dairy food intake was not associated with diabetes.

CONCLUSIONS: American Indians who participated in the SHFS reported low dairy food intake. Participants who reported higher full-fat dairy food intake had a lower risk of diabetes than participants who reported lower intake. These findings may be of interest to populations with low dairy food intake.

Although hundreds of genome-wide association studies-implicated loci have been reported for adult obesity-related traits, less is known about the genetics specific for early-onset obesity and with only a few studies conducted in non-European populations to date. Searching for additional genetic variants associated with childhood obesity, we performed a trans-ancestral meta-analysis of 30 studies consisting of up to 13 005 cases (≥95th percentile of body mass index (BMI) achieved 2-18 years old) and 15 599 controls (consistently <50th percentile of BMI) of European, African, North/South American and East Asian ancestry. Suggestive loci were taken forward for replication in a sample of 1888 cases and 4689 controls from seven cohorts of European and North/South American ancestry. In addition to observing 18 previously implicated BMI or obesity loci, for both early and late onset, we uncovered one completely novel locus in this trans-ancestral analysis (nearest gene, METTL15). The variant was nominally associated with only the European subgroup analysis but had a consistent direction of effect in other ethnicities. We then utilized trans-ancestral Bayesian analysis to narrow down the location of the probable causal variant at each genome-wide significant signal. Of all the fine-mapped loci, we were able to narrow down the causative variant at four known loci to fewer than 10 single nucleotide polymorphisms (SNPs) (FAIM2, GNPDA2, MC4R and SEC16B loci). In conclusion, an ethnically diverse setting has enabled us to both identify an additional pediatric obesity locus and further fine-map existing loci.

BACKGROUND: Increased serum levels of C-reactive protein (CRP), an important component of the innate immune response, are associated with increased risk of cardiovascular disease (CVD). Multiple single nucleotide polymorphisms (SNP) have been identified which are associated with CRP levels, and Mendelian randomization studies have shown a positive association between SNPs increasing CRP expression and risk of colon cancer (but thus far not CVD). The effects of individual genetic variants often interact with the genetic background of a population and hence we sought to resolve the genetic determinants of serum CRP in a number of American Indian populations.

METHODS: The Strong Heart Family Study (SHFS) has serum CRP measurements from 2428 tribal members, recruited as large families from three regions of the United States. Microsatellite markers and MetaboChip defined SNP genotypes were incorporated into variance components, decomposition-based linkage and association analyses.

RESULTS: CRP levels exhibited significant heritability (h2 = 0.33 ± 0.05, p<1.3 X 10-20). A locus on chromosome (chr) 6, near marker D6S281 (approximately at 169.6 Mb, GRCh38/hg38) showed suggestive linkage (LOD = 1.9) to CRP levels. No individual SNPs were found associated with CRP levels after Bonferroni adjustment for multiple testing (threshold <7.77 x 10-7), however, we found nominal associations, many of which replicate previous findings at the CRP, HNF1A and 7 other loci. In addition, we report association of 46 SNPs located at 7 novel loci on chromosomes 2, 5, 6(2 loci), 9, 10 and 17, with an average of 15.3 Kb between SNPs and all with p-values less than 7.2 X 10-4.

CONCLUSION: In agreement with evidence from other populations, these data show CRP serum levels are under considerable genetic influence; and include loci, such as near CRP and other genes, that replicate results from other ethnic groups. These findings also suggest possible novel loci on chr 6 and other chromosomes that warrant further investigation.

BACKGROUND: Chronic exposure to inorganic arsenic (iAs) in the US occurs mainly through drinking water and diet. Although American Indian (AI) populations have elevated urinary arsenic concentrations compared to the general US population, dietary sources of arsenic exposure in AI populations are not well characterized.

METHODS: We evaluated food frequency questionnaires to determine the major dietary sources of urinary arsenic concentrations (measured as the sum of arsenite, arsenate, monomethylarsonate, and dimethylarsinate, ΣAs) for 1727 AI participants in the Strong Heart Family Study (SHFS). We compared geometric mean ratios (GMRs) of urinary ΣAs for an interquartile range (IQR) increase in reported food group consumption. Exploratory analyses were stratified by gender and study center.

RESULTS: In fully adjusted generalized estimating equation models, the percent increase (95% confidence interval) of urinary ΣAs per increase in reported food consumption corresponding to the IQR was 13% (5%, 21%) for organ meat, 8% (4%, 13%) for rice, 7% (2%, 13%) for processed meat, and 4% (1%, 7%) for non-alcoholic drinks. In analyses stratified by study center, the association with organ meat was only observed in North/South Dakota. Consumption of red meat [percent increase -7% (-11%, -3%)] and fries and chips [-6% (-10%, -2%)] was inversely associated with urinary ΣAs.

CONCLUSIONS: Organ meat, processed meat, rice, and non-alcoholic drinks contribute to ΣAs exposure in the SHFS population. Organ meat is a unique source of ΣAs exposure for North and South Dakota participants and may reflect local food consumption. Further studies should comprehensively evaluate drinking water arsenic in SHFS communities and determine the relative contribution of diet and drinking water to total arsenic exposure.

While studies have reported genetic loci affecting serum urate (SU) concentrations, few studies have been conducted in minority populations. Our objective for this study was to identify genetic loci regulating SU in a multigenerational family-based cohort of American Indians, the Strong Heart Family Study (SHFS). We genotyped 162,718 single nucleotide polymorphisms (SNPs) in 2000 SHFS participants using an Illumina MetaboChip array. A genome-wide association analysis of SU was conducted using measured genotype analysis approach accounting for kinships in SOLAR, and meta-analysis in METAL. Our results showed strong association of SU with rs4481233, rs9998811, rs7696092 and rs13145758 (minor allele frequency (MAF) = 25-44%; P < 3 × 10-14) of solute carrier family 2, member 9 (SLC2A9) and rs41481455, rs2231142 and rs1481012 (MAF = 29%; p < 3 × 10-9) of ATP-binding cassette protein, subfamily G, member 2 (ABCG2). Carriers of G alleles of rs9998811, rs4148155 and rs1481012 and A alleles of rs4481233, rs7696092 and rs13145758 and rs2231142 had lower SU concentrations as compared to non-carriers. Genetic analysis of SU conditional on significant SLC2A9 and ABCG2 SNPs revealed new loci, nucleobindin 1 (NUCB1) and neuronal PAS domain protein 4 (NPAS4) (p <6× 10-6). To identify American Indian-specific SNPs, we conducted targeted sequencing of key regions of SLC2A9. A total of 233 SNPs were identified of which 89 were strongly associated with SU (p < 7.1 × 10-10) and 117 were American Indian specific. Analysis of key SNPs in cohorts of Mexican-mestizos, European, Indian and East Asian ancestries showed replication of common SNPs, including our lead SNPs. Our results demonstrate the association of SU with uric acid transporters in a minority population of American Indians and potential novel associations of SU with neuronal-related genes which warrant further investigation.

We explored arsenic-gene interactions influencing pancreatic beta-cell activity in the Strong Heart Family Study (SHFS). We considered 42 variants selected for associations with either beta-cell function (31 variants) or arsenic metabolism (11 variants) in the SHFS. Beta-cell function was calculated as homeostatic model - beta corrected for insulin resistance (cHOMA-B) by regressing homeostatic model - insulin resistance (HOMA-IR) on HOMA-B and adding mean HOMA-B. Arsenic exposure was dichotomized at the median of the sum of creatinine-corrected inorganic and organic arsenic species measured by high performance liquid chromatography-inductively coupled plasma mass spectrometry (HPLC-ICPMS). Additive GxE models for cHOMA-B were adjusted for age and ancestry, and accounted for family relationships. Models were stratified by center (Arizona, Oklahoma, North Dakota and South Dakota) and meta-analyzed. The two interactions between higher vs. lower arsenic and SNPs for cHOMA-B that were nominally significant at P < 0.05 were with rs10738708 (SNP overall effect -3.91, P = 0.56; interaction effect with arsenic -31.14, P = 0.02) and rs4607517 (SNP overall effect +16.61, P = 0.03; interaction effect with arsenic +27.02, P = 0.03). The corresponding genes GCK and TUSC1 suggest oxidative stress and apoptosis as possible mechanisms for arsenic impacts on beta-cell function. No interactions were Bonferroni-significant (1.16 × 10-3). Our findings are suggestive of oligogenic moderation of arsenic impacts on pancreatic β-cell endocrine function, but were not Bonferroni-significant.

Telomeres are repeating regions of DNA that cap chromosomes. They shorten over the mammalian life span, especially in the presence of oxidative stress and inflammation. Telomeres may play a direct role in cell senescence, serving as markers of premature vascular aging. Leukocyte telomere length (LTL) may be associated with premature vascular brain injury and cerebral atrophy. However, reports have been inconsistent, especially among minority populations with a heavy burden of illness related to vascular aging. We examined associations between LTL and magnetic resonance imaging in 363 American Indians aged 64-93 years from the Strong Heart Study (1989-1991) and its ancillary study, Cerebrovascular Disease and Its Consequences in American Indians (2010-2013). Our results showed significant associations of LTL with ventricular enlargement and the presence of white matter hyperintensities. Secondary models indicated that renal function may mediate these associations, although small case numbers limited inference. Hypertension and diabetes showed little evidence of effect modification. Results were most extreme among participants who evinced the largest decline in LTL. Although this study was limited to cross-sectional comparisons, it represents (to our knowledge) the first consideration of associations between telomere length and brain aging in American Indians. Findings suggest a relationship between vascular aging by cell senescence and severity of brain disease.

Inorganic arsenic exposure is ubiquitous, and both exposure and interindividual differences in its metabolism have been associated with cardiometabolic risk. However, the associations of arsenic exposure and arsenic metabolism with the metabolic syndrome (MetS) and its individual components are relatively unknown. We used Poisson regression with robust variance to evaluate the associations of baseline arsenic exposure (urinary arsenic levels) and metabolism (relative percentage of arsenic species over their sum) with incident MetS and its individual components (elevated waist circumference, elevated triglycerides, reduced high-density lipoprotein cholesterol, hypertension, and elevated fasting plasma glucose) in 1,047 participants from the Strong Heart Family Study, a prospective family-based cohort study in American Indian communities (baseline visits were held in 1998-1999 and 2001-2003, follow-up visits in 2001-2003 and 2006-2009). Over the course of follow-up, 32% of participants developed MetS. An interquartile-range increase in arsenic exposure was associated with a 1.19-fold (95% confidence interval: 1.01, 1.41) greater risk of elevated fasting plasma glucose concentration but not with other individual components of the MetS or MetS overall. Arsenic metabolism, specifically lower percentage of monomethylarsonic acid and higher percentage of dimethylarsinic acid, was associated with higher risk of overall MetS and elevated waist circumference but not with any other MetS component. These findings support the hypothesis that there are contrasting and independent associations of arsenic exposure and arsenic metabolism with metabolic outcomes which may contribute to overall diabetes risk.

Background: Genetic research may inform underlying mechanisms for disparities in the burden of type 2 diabetes mellitus among American Indians. Our objective was to assess the association of genetic variants in cardiometabolic candidate genes with B cell dysfunction via HOMA-B, insulin resistance via HOMA-IR, and type 2 diabetes mellitus in the Strong Heart Family Study (SHFS). Methods and Results: We examined the association of variants, previously associated with cardiometabolic traits (∼200,000 from Illumina Cardio MetaboChip), using mixed models of HOMA-B residuals corrected for HOMA-IR (cHOMA-B), log transformed HOMA-IR, and incident diabetes, adjusted for age, sex, population stratification, and familial relatedness. Center-specific estimates were combined using fixed effect meta-analyses. We used Bonferroni correction to account for multiple testing (P < 4.13 × 10-7). We also assessed the association between variants in candidate diabetes genes with these metabolic traits. We explored the top SNPs in an independent, replication sample from Southwestern Arizona. We identified significant associations with cHOMA-B for common variants at 26 loci of which 8 were novel (PRSS7, FCRL5, PEL1, LRP12, IGLL1, ARHGEF10, PARVA, FLJ16686). The most significant variant association with cHOMA-B was observed on chromosome 5 for an intergenic variant near PARP8 (rs2961831, P = 6.39 × 10-9). In the replication study, we found a signal at rs4607517 near GCK/YKT6 (P = 0.01). Variants near candidate diabetes genes (especially GCK and KCNQ1) were also nominally associated with HOMA-IR and cHOMA-B. Conclusion: We identified variants at novel loci and confirmed those at known candidate diabetes loci associations for cHOMA-B. This study also provided evidence for association of variants at KCNQ2, CTNAA2, and KCNQ1with cHOMA-B among American Indians. Further studies are needed to account for the high heritability of diabetes among the American Indian participants of the SHFS cohort.