Publications

BACKGROUND: American Indian (AI) communities are affected by uranium exposure from abandoned mines and naturally contaminated drinking water. Few studies have evaluated geographical differences across AI communities and the role of dietary exposures.

OBJECTIVE: We evaluated differences in urinary uranium levels by diet and geographical area among AI participants from the Northern Plains, the Southern Plains, and the Southwest enrolled in the Strong Heart Family Study (SHFS).

METHODS: We used food frequency questionnaires to determine dietary sources related to urinary uranium levels for 1,682 SHFS participants in 2001-2003. We calculated adjusted geometric mean ratios (GMRs) of urinary uranium for an interquartile range (IQR) increase in self-reported food group consumption accounting for family clustering and adjusting for sociodemographic variables and other food groups. We determined the percentage of variability in urinary uranium explained by diet.

RESULTS: Median (IQR) urinary uranium levels were 0.027 (0.012, 0.057) μg/g creatinine. Urinary uranium levels were higher in Arizona (median 0.039 μg/g) and North Dakota and South Dakota (median 0.038 μg/g) and lower in Oklahoma (median 0.019 μg/g). The adjusted percent increase (95% confidence interval) of urinary uranium levels per IQR increase in reported food intake was 20% (5%, 36%) for organ meat, 11% (1%, 23%) for cereals, and 14% (1%, 29%) for alcoholic drinks. In analyses stratified by study center, the association with organ meat was specific to North Dakota and South Dakota participants. An IQR increase in consumption of fries and chips was inversely associated with urinary uranium levels -11% (-19%, -3%). Overall, we estimated that self-reported dietary exposures explained 1.71% of variability in urine uranium levels.

IMPACT: Our paper provides a novel assessment of self-reported food intake and urinary uranium levels in a cohort of American Indian participants. We identify foods (organ meat, cereals, and alcohol) positively associated with urinary uranium levels, find that organ meat consumption is only associated with urine uranium in North Dakota and South Dakota, and estimate that diet explains relatively little variation in total urinary uranium concentrations. Our findings contribute meaningful data toward a more comprehensive estimation of uranium exposure among Native American communities and support the need for high-quality assessments of water and dust uranium exposures in SHFS communities.

BACKGROUND: Chronic kidney disease (CKD) is highly prevalent in Central America, and genetic factors may contribute to CKD risk. To understand the influences of genetic admixture on CKD susceptibility, we conducted an admixture mapping screening of CKD traits and risk factors in US Hispanic and Latino individuals from Central America country of origin.

METHODS: We analyzed 1023 participants of HCHS/SOL (Hispanic Community Health Study/Study of Latinos) who reported 4 grandparents originating from the same Central America country. Ancestry admixture findings were validated on 8191 African Americans from WHI (Women's Health Initiative), 3141 American Indians from SHS (Strong Heart Study), and over 1.1 million European individuals from a multistudy meta-analysis.

RESULTS: We identified 3 novel genomic regions for albuminuria (chromosome 14q24.2), CKD (chromosome 6q25.3), and type 2 diabetes (chromosome 3q22.2). The 14q24.2 locus driven by a Native American ancestry had a protective effect on albuminuria and consisted of 2 nearby regions spanning the RGS6 gene. Variants at this locus were validated in American Indians. The 6q25.3 African ancestry-derived locus, encompassing the ARID1B gene, was associated with increased risk for CKD and replicated in African Americans through admixture mapping. The European ancestry type 2 diabetes locus at 3q22.2, encompassing the EPHB1 and KY genes, was validated in European individuals through variant association.

CONCLUSIONS: US Hispanic/Latino populations are culturally and genetically diverse. This study focusing on Central America grandparent country of origin provides new loci discovery and insights into the ancestry-of-origin influences on CKD and risk factors in US Hispanic and Latino individuals.

The last couple of decades have highlighted the importance of studying hybridization, particularly among primate species, as it allows us to better understand our own evolutionary trajectory. Here, we report on genetic ancestry estimates using dense, full genome data from 881 olive (Papio anubus), yellow (Papio cynocephalus), or olive-yellow crossed captive baboons from the Southwest National Primate Research Center. We calculated global and local ancestry information, imputed low coverage genomes (n = 830) to improve marker quality, and updated the genetic resources of baboons available to assist future studies. We found evidence of historical admixture in some putatively purebred animals and identified errors within the Southwest National Primate Research Center pedigree. We also compared the outputs between two different phasing and imputation pipelines along with two different global ancestry estimation software. There was good agreement between the global ancestry estimation software, with R2 > 0.88, while evidence of phase switch errors increased depending on what phasing and imputation pipeline was used. We also generated updated genetic maps and created a concise set of ancestry informative markers (n = 1,747) to accurately obtain global ancestry estimates.

BACKGROUND: A growing body of research indicates that associations of ceramides and sphingomyelins with mortality depend on the chain length of the fatty acid acylated to the backbone sphingoid base. We examined associations of 8 ceramide and sphingomyelin species with mortality among an American Indian population.

METHODS AND RESULTS: The analysis comprised 2688 participants from the SHFS (Strong Heart Family Study). Plasma ceramide and sphingomyelin species carrying long-chain (ie, 16:0) and very-long-chain (ie, 20:0, 22:0, 24:0) saturated fatty acids were measured by sequential liquid chromatography and mass spectroscopy using samples from 2001 to 2003. Participants were followed for 18.8 years (2001-2020). Associations of ceramides and sphingomyelins with mortality were assessed using Cox models. The mean age of participants was 40.8 years. There were 574 deaths during a median 17.4-year follow-up. Ceramides and sphingomyelins carrying fatty acid 16:0 were positively associated with mortality. Ceramides and sphingomyelins carrying longer fatty acids were inversely associated with mortality. Per SD difference in each ceramide and sphingomyelin species, hazard ratios for death were: 1.68 (95% CI, 1.44-1.96) for ceramide-16 (Cer-16), 0.82 (95% CI, 0.71-0.95) for Cer-20, 0.60 (95% CI, 0.51-0.70) for Cer-22, 0.67 (95% CI, 0.56-0.79) for Cer-24, 1.80 (95% CI-1.57, 2.05) for sphingomyelin-16 (SM-16), 0.54 (95% CI, 0.47-0.62) for SM-20, 0.50 (95% CI, 0.44-0.57) for SM-22, and 0.59 (95% CI, 0.52-0.67) for SM-24.

CONCLUSIONS: The direction/magnitude of associations of ceramides and sphingomyelins with mortality differs according to the length of the fatty acid acylated to the backbone sphingoid base.

REGISTRATION: URL: https://www.clinicatrials.gov; Unique identifier: NCT00005134.

The marmoset is a fundamental nonhuman primate model for the study of aging, neurobiology, and many other topics. Genetic management of captive marmoset colonies is complicated by frequent chimerism in the blood and other tissues, a lack of tools to enable cost-effective, genome-wide interrogation of variation, and historic mergers and migrations of animals between colonies. We implemented genotype-by-sequencing (GBS) of hair follicle derived DNA (a minimally chimeric DNA source) of 82 marmosets housed at the Southwest National Primate Research Center (SNPRC). Our primary goals were the genetic characterization of our marmoset population for pedigree verification and colony management and to inform the scientific community of the functional genetic makeup of this valuable resource. We used the GBS data to reconstruct the genetic legacy of recent mergers between colonies, to identify genetically related animals whose relationships were previously unknown due to incomplete pedigree information, and to show that animals in the SNPRC colony appear to exhibit low levels of inbreeding. Of the >99,000 single-nucleotide variants (SNVs) that we characterized, >9800 are located within gene regions known to harbor pathogenic variants of clinical significance in humans. Overall, we show the combination of low-resolution (sparse) genotyping using hair follicle DNA is a powerful strategy for the genetic management of captive marmoset colonies and for identifying potential SNVs for the development of biomedical research models.

BACKGROUND: Although COVID-19 infection has been associated with a number of clinical and environmental risk factors, host genetic variation has also been associated with the incidence and morbidity of infection. The CRP gene codes for a critical component of the innate immune system and CRP variants have been reported associated with infectious disease and vaccination outcomes. We investigated possible associations between COVID-19 outcome and a limited number of candidate gene variants including rs1205.

METHODOLOGY/PRINCIPAL FINDINGS: The Strong Heart and Strong Heart Family studies have accumulated detailed genetic, cardiovascular risk and event data in geographically dispersed American Indian communities since 1988. Genotypic data and 91 COVID-19 adjudicated deaths or hospitalizations from 2/1/20 through 3/1/23 were identified among 3,780 participants in two subsets. Among 21 candidate variants including genes in the interferon response pathway, APOE, TMPRSS2, TLR3, the HLA complex and the ABO blood group, only rs1205, a 3' untranslated region variant in the CRP gene, showed nominally significant association in T-dominant model analyses (odds ratio 1.859, 95%CI 1.001-3.453, p = 0.049) after adjustment for age, sex, center, body mass index, and a history of cardiovascular disease. Within the younger subset, association with the rs1205 T-Dom genotype was stronger, both in the same adjusted logistic model and in the SOLAR analysis also adjusting for other genetic relatedness.

CONCLUSION: A T-dominant genotype of rs1205 in the CRP gene is associated with COVID-19 death or hospitalization, even after adjustment for relevant clinical factors and potential participant relatedness. Additional study of other populations and genetic variants of this gene are warranted.

BACKGROUND: Identifying lipidomic markers of diet quality is needed to inform the development of biomarkers of diet, and to understand the mechanisms driving the diet- coronary heart disease (CHD) association.

OBJECTIVES: This study aimed to identify lipidomic markers of diet quality and examine whether these lipids are associated with incident CHD.

METHODS: Using liquid chromatography-mass spectrometry, we measured 1542 lipid species from 1694 American Indian adults (aged 18-75 years, 62% female) in the Strong Heart Family Study. Participants were followed up for development of CHD through 2020. Information on the past year diet was collected using the Block Food Frequency Questionnaire, and diet quality was assessed using the Alternative Healthy Eating Index-2010 (AHEI). Mixed-effects linear regression was used to identify individual lipids cross-sectionally associated with AHEI. In prospective analysis, Cox frailty model was used to estimate the hazard ratio (HR) of each AHEI-related lipid for incident CHD. All models were adjusted for age, sex, center, education, body mass index, smoking, alcohol drinking, level of physical activity, energy intake, diabetes, hypertension, and use of lipid-lowering drugs. Multiple testing was controlled at a false discovery rate of <0.05.

RESULTS: Among 1542 lipid species measured, 71 lipid species (23 known), including acylcarnitine, cholesterol esters, glycerophospholipids, sphingomyelins and triacylglycerols, were associated with AHEI. Most of the identified lipids were associated with consumption of ω-3 (n-3) fatty acids. In total, 147 participants developed CHD during a mean follow-up of 17.8 years. Among the diet-related lipids, 10 lipids [5 known: cholesterol ester (CE)(22:5)B, phosphatidylcholine (PC)(p-14:0/22:1)/PC(o-14:0/22:1), PC(p-38:3)/PC(o-38:4)B, phosphatidylethanolamine (PE)(p-18:0/20:4)/PE(o-18:0/20:4), and sphingomyelin (d36:2)A] were associated with incident CHD. On average, each standard deviation increase in the baseline level of these 5 lipids was associated with 17%-23% increased risk of CHD (from HR: 1.17; 95% CI: 1, 1.36; to HR: 1.23; 95% CI: 1.05, 1.43).

CONCLUSIONS: In this study, lipidomic markers of diet quality in American Indian adults are found. Some diet-related lipids are associated with risk of CHD beyond established risk factors.

BACKGROUND: Although many studies on the association between dyslipidemia and cardiovascular disease (CVD) exist in older adults, data on the association among adolescents and young adults living with disproportionate burden of cardiometabolic disorders are scarce.

METHODS AND RESULTS: The SHFS (Strong Heart Family Study) is a multicenter, family-based, prospective cohort study of CVD in an American Indian populations, including 12 communities in central Arizona, southwestern Oklahoma, and the Dakotas. We evaluated SHFS participants, who were 15 to 39 years old at the baseline examination in 2001 to 2003 (n=1440). Lipids were measured after a 12-hour fast. We used carotid ultrasounds to detect plaque at baseline and follow-up in 2006 to 2009 (median follow-up=5.5 years). We identified incident CVD events through 2020 with a median follow-up of 18.5 years. We used shared frailty proportional hazards models to assess the association between dyslipidemia and subclinical or clinical CVD, while controlling for covariates. Baseline dyslipidemia prevalence was 55.2%, 73.6%, and 78.0% for participants 15 to 19, 20 to 29, and 30 to 39 years old, respectively. Approximately 2.8% had low-density lipoprotein cholesterol ≥160 mg/dL, which is higher than the recommended threshold for lifestyle or medical interventions in young adults of 20 to 39 years old. During follow-up, 9.9% had incident plaque (109/1104 plaque-free participants with baseline and follow-up ultrasounds), 11.0% had plaque progression (128/1165 with both baseline and follow-up ultrasounds), and 9% had incident CVD (127/1416 CVD-free participants at baseline). Plaque incidence and progression were higher in participants with total cholesterol ≥200 mg/dL, low-density lipoprotein cholesterol ≥160 mg/dL, or non-high-density lipoprotein cholesterol ≥130 mg/dL, while controlling for covariates. CVD risk was independently associated with low-density lipoprotein cholesterol ≥160 mg/dL.

CONCLUSIONS: Dyslipidemia is a modifiable risk factor that is associated with both subclinical and clinical CVD, even among the younger American Indian population who have unexpectedly high rates of significant CVD events. Therefore, this population is likely to benefit from a variety of evidence-based interventions including screening, educational, lifestyle, and guideline-directed medical therapy at an early age.

The generalized estimating equations method (GEE) is commonly applied to analyze data obtained from family studies. GEE is well known for its robustness on misspecification of correlation structure. However, the unbalanced distribution of family sizes and complicated genetic relatedness structure within each family may challenge GEE performance. We focused our research on binary outcomes. To evaluate the performance of GEE, we conducted a series of simulations, on data generated adopting the kinship matrix (correlation structure within each family) from the Strong Heart Family Study (SHFS). We performed a fivefold cross-validation to further evaluate the GEE predictive power on data from the SHFS. A Bayesian modeling approach, with direct integration of the kinship matrix, was also included to contrast with GEE. Our simulation studies revealed that GEE performs well on a binary outcome from families having a relatively simple kinship structure. However, data with a binary outcome generated from families with complex kinship structures, especially with a large genetic variance, can challenge the performance of GEE.

BACKGROUND: Cardiovascular disease (CVD) is a leading cause of morbidity and mortality in American Indian people. In 2022, the American Heart Association developed the Life's Essential 8 goals to promote cardiovascular health (CVH) for Americans, composed of diet, physical activity, nicotine exposure, sleep, body mass index, blood lipids, blood pressure, and blood glucose. We examined whether achievement of Life's Essential 8 goals was associated with incident CVD among SHFS (Strong Heart Family Study) participants.

METHODS AND RESULTS: A total of 2139 SHFS participants without CVD at baseline were included in analyses. We created a composite CVH score based on achievement of Life's Essential 8 goals, excluding sleep. Scores of 0 to 49 represented low CVH, 50 to 69 represented moderate CVH, and 70 to 100 represented high CVH. Incident CVD was defined as incident myocardial infarction, coronary heart disease, congestive heart failure, or stroke. Cox proportional hazard models were used to examine the relationship of CVH and incident CVD. The incidence rate of CVD at the 20-year follow-up was 7.43 per 1000 person-years. Compared with participants with low CVH, participants with moderate and high CVH had a lower risk of incident CVD; the hazard ratios and 95% CIs for incident CVD for moderate and high CVH were 0.52 (95% CI, 0.40-0.68) and 0.25 (95% CI, 0.14-0.44), respectively, after adjustment for age, sex, education, and study site.

CONCLUSIONS: Better CVH was associated with lower CVD risk which highlights the need for comprehensive public health interventions targeting CVH promotion to reduce CVD risk in American Indian communities.