Publications

2019

Best, Lyle G, Poojitha Balakrishnan, Shelley A Cole, Karin Haack, Jonathan M Kocarnik, Nathan Pankratz, Matthew Z Anderson, et al. (2019) 2019. “Genetic Analysis of HsCRP in American Indians: The Strong Heart Family Study.”. PloS One 14 (10): e0223574. https://doi.org/10.1371/journal.pone.0223574.

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.

Nigra, Anne E, Pablo Olmedo, Maria Grau-Perez, Rae O’Leary, Marcia O’Leary, Amanda M Fretts, Jason G Umans, et al. (2019) 2019. “Dietary Determinants of Inorganic Arsenic Exposure in the Strong Heart Family Study.”. Environmental Research 177: 108616. https://doi.org/10.1016/j.envres.2019.108616.

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.

Chittoor, Geetha, Karin Haack, Poojitha Balakrishnan, Christopher Bizon, Sandra Laston, Lyle G Best, Jean W MacCluer, et al. (2019) 2019. “Fine Mapping and Identification of Serum Urate Loci in American Indians: The Strong Heart Family Study.”. Scientific Reports 9 (1): 17899. https://doi.org/10.1038/s41598-019-52924-w.

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.

2018

Balakrishnan, Poojitha, Ana Navas-Acien, Karin Haack, Dhananjay Vaidya, Jason G Umans, Lyle G Best, Walter Goessler, et al. (2018) 2018. “Arsenic-Gene Interactions and Beta-Cell Function in the Strong Heart Family Study.”. Toxicology and Applied Pharmacology 348: 123-29. https://doi.org/10.1016/j.taap.2018.03.034.

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.

Suchy-Dicey, Astrid M, Clemma J Muller, Tara M Madhyastha, Dean Shibata, Shelley A Cole, Jinying Zhao, W T Longstreth, and Dedra Buchwald. (2018) 2018. “Telomere Length and Magnetic Resonance Imaging Findings of Vascular Brain Injury and Central Brain Atrophy: The Strong Heart Study.”. American Journal of Epidemiology 187 (6): 1231-39. https://doi.org/10.1093/aje/kwx368.

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.

Spratlen, Miranda J, Maria Grau-Perez, Lyle G Best, Joseph Yracheta, Mariana Lazo, Dhananjay Vaidya, Poojitha Balakrishnan, et al. (2018) 2018. “The Association of Arsenic Exposure and Arsenic Metabolism With the Metabolic Syndrome and Its Individual Components: Prospective Evidence From the Strong Heart Family Study.”. American Journal of Epidemiology 187 (8): 1598-1612. https://doi.org/10.1093/aje/kwy048.

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.

Balakrishnan, Poojitha, Dhananjay Vaidya, Saroja Voruganti, Karin Haack, Jack W Kent, Kari E North, Sandra Laston, et al. (2018) 2018. “Genetic Variants Related to Cardiometabolic Traits Are Associated to B Cell Function, Insulin Resistance, and Diabetes Among AmeriCan Indians: The Strong Heart Family Study.”. Frontiers in Genetics 9: 466. https://doi.org/10.3389/fgene.2018.00466.

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.

Bastarrachea, Raul A, Hugo A Laviada-Molina, Edna J Nava-Gonzalez, Irene Leal-Berumen, Claudia Escudero-Lourdes, Fabiola Escalante-Araiza, Vanessa-Giselle Peschard, et al. (2018) 2018. “Deep Multi-OMICs and Multi-Tissue Characterization in a Pre- and Postprandial State in Human Volunteers: The GEMM Family Study Research Design.”. Genes 9 (11). https://doi.org/10.3390/genes9110532.

Cardiovascular disease (CVD) and type 2 diabetes (T2D) are increasing worldwide. This is mainly due to an unhealthy nutrition, implying that variation in CVD risk may be due to variation in the capacity to manage a nutritional load. We examined the genomic basis of postprandial metabolism. Our main purpose was to introduce the GEMM Family Study (Genetics of Metabolic Diseases in Mexico) as a multi-center study carrying out an ongoing recruitment of healthy urban adults. Each participant received a mixed meal challenge and provided a 5-hours' time course series of blood, buffy coat specimens for DNA isolation, and adipose tissue (ADT)/skeletal muscle (SKM) biopsies at fasting and 3 h after the meal. A comprehensive profiling, including metabolomic signatures in blood and transcriptomic and proteomic profiling in SKM and ADT, was performed to describe tendencies for variation in postprandial response. Our data generation methods showed preliminary trends indicating that by characterizing the dynamic properties of biomarkers with metabolic activity and analyzing multi-OMICS data it could be possible, with this methodology and research design, to identify early trends for molecular biology systems and genes involved in the fasted and fed states.

Spratlen, Miranda J, Maria Grau-Perez, Jason G Umans, Joseph Yracheta, Lyle G Best, Kevin Francesconi, Walter Goessler, et al. (2018) 2018. “Arsenic, One Carbon Metabolism and Diabetes-Related Outcomes in the Strong Heart Family Study.”. Environment International 121 (Pt 1): 728-40. https://doi.org/10.1016/j.envint.2018.09.048.

BACKGROUND: Inorganic arsenic exposure and inter-individual differences in its metabolism have been associated with cardiometabolic risk. A more efficient arsenic metabolism profile (lower MMA%, higher DMA%) has been associated with reduced risk for arsenic-related health outcomes; however, this profile has also been associated with increased risk for diabetes-related outcomes. The mechanism behind these contrasting associations is equivocal; we hypothesized one carbon metabolism (OCM) may play a role.

METHODS: We evaluated the association between OCM-related variables (nutrient intake and genetic variants) and both arsenic metabolism biomarkers (iAs%, MMA% and DMA%) and diabetes-related outcomes (metabolic syndrome, diabetes, HOMA2-IR and waist circumference) in 935 participants free of prevalent diabetes and metabolic syndrome from the Strong Heart Family Study, a family-based prospective cohort comprised of American Indian tribal members aged 14+ years.

RESULTS: Of the 935 participants free of both diabetes and metabolic syndrome at baseline, 279 (29.8%) developed metabolic syndrome over a median of 5.3 years of follow-up and of the 1458 participants free of diabetes at baseline, 167 (11.3%) developed diabetes over follow-up. OCM nutrients were not associated with arsenic metabolism, however, higher vitamin B6 was associated with diabetes-related outcomes (higher HOMA2-IR and increased risk for diabetes and metabolic syndrome). A polymorphism in an OCM-related gene, methionine synthase (MTR), was associated with both higher MMA% (β = 2.57, 95% CI: 0.22, 4.92) and lower HOMA2-IR (GMR = 0.79, 95% CI = 0.66, 0.93 per 5 years of follow-up). Adjustment for OCM variables did not affect previously reported associations between arsenic metabolism and diabetes-related outcomes; however, the association between the MTR variant and diabetes-related outcomes were attenuated after adjustment for arsenic metabolism.

CONCLUSIONS: Our findings suggest MMA% may be a partial mediator in the association between OCM and diabetes-related outcomes. Additional mediation analyses with longer follow-up period are needed to confirm this finding. Further research is needed to determine whether excess B vitamin intake is associated with increased risk for diabetes-related outcomes.

2017

Peng, Hao, Mihriye Mete, Sameer Desale, Amanda M Fretts, Shelley A Cole, Lyle G Best, Jue Lin, et al. (2017) 2017. “Leukocyte Telomere Length and Ideal Cardiovascular Health in American Indians: The Strong Heart Family Study.”. European Journal of Epidemiology 32 (1): 67-75. https://doi.org/10.1007/s10654-016-0199-6.

Telomere length, a marker of biological aging, has been associated with cardiovascular disease (CVD) and its risk factors. Ideal cardiovascular health (CVH), defined by the American Heart Association (AHA), has also been associated with a reduced risk of CVD, but the relationship between telomere length and ideal CVH is unclear. We measured leukocyte telomere length (LTL) by qPCR in 2568 American Indians in the Strong Heart Family Study (SHFS). All participants were free of overt CVD at enrollment (2001-2003). CVH indices included four behavioral factors (smoking, physical activity, diet, BMI) and three health factors (blood pressure, cholesterol, fasting glucose). Each index was categorized as poor, intermediate, or ideal according to the AHA's guideline. CVH was further categorized into below average (0-1), average (2-3) and above average (≥4) based on the total number of ideal indices. Results showed that, 29, 50 and 21 % of study participants had below average, average, and above average CVH, respectively. Participants with above average CVH had significantly longer LTL than those with below average CVH (β = 0.034, P = 0.042) after adjusting for age, sex, education level, marital status, processed meat consumption, alcohol consumption, and study site. Compared to the U.S. general population, American Indians achieved lower rates for five out of the seven ideal CVH metrics, including smoking, BMI, physical activity, diet, and blood pressure. Achieving four or more ideal CVH metrics was significantly associated with longer LTL. This finding suggests that achieving an ideal CVH may prevent or delay CVD, probably through promoting healthy aging.