Publications

OBJECTIVE: Comprehensive assessment of alterations in lipid species preceding type 2 diabetes (T2D) is largely unknown. We aimed to identify plasma molecular lipids associated with risk of T2D in American Indians.

RESEARCH DESIGN AND METHODS: Using untargeted liquid chromatography-mass spectrometry, we repeatedly measured 3,907 fasting plasma samples from 1,958 participants who attended two examinations (∼5.5 years apart) and were followed up to 16 years in the Strong Heart Family Study. Mixed-effects logistic regression was used to identify lipids associated with risk of T2D, adjusting for traditional risk factors. Repeated measurement analysis was performed to examine the association between change in lipidome and change in continuous measures of T2D, adjusting for baseline lipids. Multiple testing was controlled by false discovery rate at 0.05.

RESULTS: Higher baseline level of 33 lipid species, including triacylglycerols, diacylglycerols, phosphoethanolamines, and phosphocholines, was significantly associated with increased risk of T2D (odds ratio [OR] per SD increase in log2-transformed baseline lipids 1.50-2.85) at 5-year follow-up. Of these, 21 lipids were also associated with risk of T2D at 16-year follow-up. Aberrant lipid profiles were also observed in prediabetes (OR per SD increase in log2-transformed baseline lipids 1.30-2.19 for risk lipids and 0.70-0.78 for protective lipids). Longitudinal changes in 568 lipids were significantly associated with changes in continuous measures of T2D. Multivariate analysis identified distinct lipidomic signatures differentiating high- from low-risk groups.

CONCLUSIONS: Lipid dysregulation occurs many years preceding T2D, and novel molecular lipids (both baseline level and longitudinal change over time) are significantly associated with risk of T2D beyond traditional risk factors. Our findings shed light on the mechanisms linking dyslipidemia to T2D and may yield novel therapeutic targets for early intervention tailored to American Indians.

Rationale: Permanent lung function impairment after active tuberculosis infection is relatively common. It remains unclear which spirometric pattern is most prevalent after tuberculosis.Objectives: Our objective was to elucidate the impact of active tuberculosis survival on lung health in the Strong Heart Study (SHS), a population of American Indians historically highly impacted by tuberculosis. As arsenic exposure has also been related to lung function in the SHS, we also assessed the joint effect between arsenic exposure and past active tuberculosis.Methods: The SHS is an ongoing population-based, prospective study of cardiovascular disease and its risk factors in American Indian adults. This study uses tuberculosis data and spirometry data from the Visit 2 examination (1993-1995). Prior active tuberculosis was ascertained by a review of medical records. Forced vital capacity (FVC), forced expiratory volume in 1 second (FEV1), and FEV1/FVC were measured by spirometry. An additional analysis was conducted to evaluate the potential association between active tuberculosis and arsenic exposure.Results: A history of active tuberculosis was associated with reduced percent predicted FVC and FEV1, an increased odds of airflow obstruction (odds ratio = 1.45, 95% confidence interval = 1.08-1.95), and spirometric restrictive pattern (odds ratio = 1.73, 95% confidence interval = 1.24-2.40). These associations persisted after adjustment for diabetes and other risk factors, including smoking. We also observed the presence of cough, phlegm, and exertional dyspnea after a history of active tuberculosis. In the additional analysis, increasing urinary arsenic concentrations were associated with decreasing lung function in those with a history of active tuberculosis, but a reduced odds of active tuberculosis was found with elevated arsenic.Conclusions: Our findings support existing knowledge that a history of active tuberculosis is a risk factor for long-term respiratory impairment. Arsenic exposure, although inversely associated with prior active tuberculosis, was associated with a further decrease in lung function among those with a prior active tuberculosis history. The possible interaction between arsenic and tuberculosis, as well as the reduced odds of tuberculosis associated with arsenic exposure, warrants further investigation, as many populations at risk of developing active tuberculosis are also exposed to arsenic-contaminated water.

CONTEXT: Osteoporosis is a major public health burden with significant economic costs. However, the correlates of bone health in Hispanic children are understudied.

OBJECTIVE: We aimed to identify genetic variants associated with bone mineral density (BMD) and bone mineral content (BMC) at multiple skeletal sites in Hispanic children.

METHODS: We conducted a cross-sectional genome-wide linkage analysis, genome-wide and exome-wide association analysis of BMD and BMC. The Viva La Familia Study is a family-based cohort with a total of 1030 Hispanic children (4-19 years old at baseline) conducted in Houston, TX. BMD and BMC were measured by Dual-energy X-ray absorptiometry.

RESULTS: Significant heritability were observed for BMC and BMD at multiple skeletal sites ranging between 44 and 68% (P < 2.8 × 10-9). Significant evidence for linkage was found for BMD of pelvis and left leg on chromosome 7p14, lumbar spine on 20q13 and left rib on 6p21, and BMC of pelvis on chromosome 20q12 and total body on 14q22-23 (logarithm of odds score > 3). We found genome-wide significant association between BMC of right arm and rs762920 at PVALB (P = 4.6 × 10-8), and between pelvis BMD and rs7000615 at PTK2B (P = 7.4 × 10-8). Exome-wide association analysis revealed novel association of variants at MEGF10 and ABRAXAS2 with left arm and lumber spine BMC, respectively (P < 9 × 10-7).

CONCLUSIONS: We identified novel loci associated with BMC and BMD in Hispanic children, with strongest evidence for PTK2B. These findings provide better understanding of bone genetics and shed light on biological mechanisms underlying BMD and BMC variation.

BACKGROUND: Chronic exposure to arsenic (As), a human toxicant and carcinogen, remains a global public health problem. Health risks persist after As exposure has ended, suggesting epigenetic dysregulation as a mechanistic link between exposure and health outcomes.

OBJECTIVES: We investigated the association between total urinary As and locus-specific DNA methylation in the Strong Heart Study, a cohort of American Indian adults with low-to-moderate As exposure [total urinary As, mean (±SD) μg/g creatinine: 11.7 (10.6)].

METHODS: DNA methylation was measured in 2,325 participants using the Illumina MethylationEPIC array. We implemented linear models to test differentially methylated positions (DMPs) and the DMRcate method to identify regions (DMRs) and conducted gene ontology enrichment analysis. Models were adjusted for estimated cell type proportions, age, sex, body mass index, smoking, education, estimated glomerular filtration rate, and study center. Arsenic was measured in urine as the sum of inorganic and methylated species.

RESULTS: In adjusted models, methylation at 20 CpGs was associated with urinary As after false discovery rate (FDR) correction (FDR< 0.05). After Bonferroni correction, 5 CpGs remained associated with total urinary As (pBonferroni<0.05), located in SLC7A11, ANKS3, LINGO3, CSNK1D, ADAMTSL4. We identified one DMR on chromosome 11 (chr11:2,322,050-2,323,247), annotated to C11orf2; TSPAN32 genes.

DISCUSSION: This is one of the first epigenome-wide association studies to investigate As exposure and locus-specific DNA methylation using the Illumina MethylationEPIC array and the largest epigenome-wide study of As exposure. The top DMP was located in SLC7A11A, a gene involved in cystine/glutamate transport and the biosynthesis of glutathione, an antioxidant that may protect against As-induced oxidative stress. Additional DMPs were located in genes associated with tumor development and glucose metabolism. Further research is needed, including research in more diverse populations, to investigate whether As-related DNA methylation signatures are associated with gene expression or may serve as biomarkers of disease development. https://doi.org/10.1289/EHP6263.

INTRODUCTION: We examined whether abuse in childhood and/or adolescence was associated with shorter telomere length in a pooled analysis of 3,232 participants from five diverse cohorts. We also assessed whether religion or spirituality (R/S) could buffer deleterious effects of abuse.

METHODS: Physical and sexual abuse in childhood (age <12) and adolescence (age 12-18) was assessed using the Revised Conflict Tactics Scale and questions from a 1995 Gallup survey. We measured relative leukocyte telomere lengths (RTL) using quantitative real time polymerase chain reaction. We used generalized estimating equations to assess associations of physical and sexual abuse with log-transformed RTL z-scores. Analyses were conducted in each cohort, overall, and stratified by extent of religiosity or spirituality and religious coping in adulthood. We pooled study-specific estimates using random-effects models and assessed between-study heterogeneity.

RESULTS: Compared to no abuse, severe sexual abuse was associated with lower RTL z-scores, in childhood: -15.6%, 95% CI: -25.9, -4.9; p-trend = 0.04; p-heterogeneity = 0.58 and in adolescence: -16.5%, 95% CI: -28.1, -3.0; p-trend = 0.08; p-heterogeneity = 0.68. Sexual abuse experienced in both childhood and adolescence was associated with 11.3% lower RTL z-scores after adjustment for childhood and demographic covariates (95% CI: -20.5%, -2.0%; p-trend = 0.03; p-heterogeneity = 0.62). There was no evidence of effect modification by R/S. Physical abuse was not associated with telomere length.

CONCLUSIONS: Sexual abuse in childhood or adolescence was associated with a marker of accelerated biological aging, decreased telomere length. The lack of moderation by R/S may be due to inability to capture the appropriate time period for those beliefs and practices.

BACKGROUND: Elevated adiposity is often posited by medical and public health researchers to be a risk factor associated with cardiovascular disease, diabetes, and other diseases. These health challenges are now thought to be reflected in epigenetic modifications to DNA molecules, such as DNA methylation, which can alter gene expression.

METHODS: Here we report the results of three Epigenome Wide Association Studies (EWAS) in which we assessed the differential methylation of DNA (obtained from peripheral blood) associated with three adiposity phenotypes (BMI, waist circumference, and impedance-measured percent body fat) among American Indian adult participants in the Strong Heart Study.

RESULTS: We found differential methylation at 8264 CpG sites associated with at least one of our three response variables. Of the three adiposity proxies we measured, waist circumference had the highest number of associated differentially methylated CpGs, while percent body fat was associated with the lowest. Because both waist circumference and percent body fat relate to physiology, we focused interpretations on these variables. We found a low degree of overlap between these two variables in our gene ontology enrichment and Differentially Methylated Region analyses, supporting that waist circumference and percent body fat measurements represent biologically distinct concepts.

CONCLUSIONS: We interpret these general findings to indicate that highly significant regions of the genome (DMR) and synthesis pathways (GO) in waist circumference analyses are more likely to be associated with the presence of visceral/abdominal fat than more general measures of adiposity. Our findings confirmed numerous CpG sites previously found to be differentially methylated in association with adiposity phenotypes, while we also found new differentially methylated CpG sites and regions not previously identified.

A large number of studies have shown that the baboon is one of the most commonly used non-human primate (NHP) research model for the study of immunometabolic complex traits such as type 2 diabetes (T2D), insulin resistance (IR), adipose tissue dysfunction (ATD), dyslipidemia, obesity (OB) and cardiovascular disease (CVD). This paper reports on innovative technologies and advanced research strategies for energetics and translational medicine with this NHP model. This includes the following: measuring resting energy expenditure (REE) with the mobile indirect calorimeter Breezing®; monitoring daily body temperature using subcutaneously implanted data loggers; quantifying metabolic heat with veterinary infrared thermography (IRT) imaging, and non-viral non-invasive, tissue-specific ultrasound-targeted microbubble destruction (UTMD) gene-based therapy. These methods are of broad utility; for example, they may facilitate the engineering of ectopic overexpression of brown adipose tissue (BAT) mUCP-1 via UTMD-gene therapy into baboon SKM to achieve weight loss, hypophagia and immunometabolic improvement. These methods will be valuable to basic and translational research, and human clinical trials, in the areas of metabolism, cardiovascular health, and immunometabolic and infectious diseases.

The rhesus macaque (Macaca mulatta) is the most widely studied nonhuman primate (NHP) in biomedical research. We present an updated reference genome assembly (Mmul_10, contig N50 = 46 Mbp) that increases the sequence contiguity 120-fold and annotate it using 6.5 million full-length transcripts, thus improving our understanding of gene content, isoform diversity, and repeat organization. With the improved assembly of segmental duplications, we discovered new lineage-specific genes and expanded gene families that are potentially informative in studies of evolution and disease susceptibility. Whole-genome sequencing (WGS) data from 853 rhesus macaques identified 85.7 million single-nucleotide variants (SNVs) and 10.5 million indel variants, including potentially damaging variants in genes associated with human autism and developmental delay, providing a framework for developing noninvasive NHP models of human disease.

BACKGROUND: Inorganic arsenic exposure is ubiquitous and both 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. This profile, however, has also been associated with increased risk for diabetes-related outcomes.

OBJECTIVES: The mechanism behind these conflicting associations is unclear; we hypothesized the one-carbon metabolism (OCM) pathway may play a role.

METHODS: We evaluated the influence of OCM on the relationship between arsenic metabolism and diabetes-related outcomes (HOMA2-IR, waist circumference, fasting plasma glucose) using metabolomic data from an OCM-specific and P180 metabolite panel measured in plasma, arsenic metabolism measured in urine, and HOMA2-IR and FPG measured in fasting plasma. Samples were drawn from baseline visits (2001-2003) in 59 participants from the Strong Heart Family Study, a family-based cohort study of American Indians aged ≥14 years from Arizona, Oklahoma, and North/South Dakota.

RESULTS: In unadjusted analyses, a 5% increase in DMA% was associated with higher HOMA2-IR (geometric mean ratio (GMR)= 1.13 (95% CI: 1.03, 1.25)) and waist circumference (mean difference=3.66 (0.95, 6.38). MMA% was significantly associated with lower HOMA2-IR and waist circumference. After adjustment for OCM-related metabolites (SAM, SAH, cysteine, glutamate, lysophosphatidylcholine 18.2, and three phosphatidlycholines), associations were attenuated and no longer significant.

CONCLUSIONS: These preliminary results indicate that the association of lower MMA% and higher DMA% with diabetes-related outcomes may be influenced by OCM status, either through confounding, reverse causality, or mediation.

INTRODUCTION: While several mechanisms may explain metal-related health effects, the exact cellular processes are not fully understood. We evaluated the association between leukocyte telomere length (LTL) and urine arsenic (ΣAs), cadmium (Cd) and tungsten (W) exposure in the Strong Heart Study (SHS, N = 1702) and in the Strong Heart Family Study (SHFS, N = 1793).

METHODS: Urine metal concentrations were measured using ICP-MS. Arsenic exposure was assessed as the sum of inorganic arsenic, monomethylarsonate and dimethylarsinate levels (ΣAs). LTL was measured by quantitative polymerase chain reaction.

RESULTS: In the SHS, median levels were 1.09 for LTL, and 8.8, 1.01 and 0.11 μg/g creatinine for ΣAs, Cd, and W, respectively. In the SHFS, median levels were 1.01 for LTL, and 4.3, 0.44, and 0.10 μg/g creatinine. Among SHS participants, increased urine ΣAs, Cd, and W was associated with shorter LTL. The adjusted geometric mean ratio (95% confidence interval) of LTL per an increase equal to the difference between the percentiles 90th and 10th in metal distributions was 0.85 (0.79, 0.92) for ΣAs, 0.91 (0.84, 1.00) for Cd and 0.93 (0.88, 0.98) for W. We observed no significant associations among SHFS participants. The findings also suggest that the association between arsenic and LTL might be differential depending on the exposure levels or age.

CONCLUSIONS: Additional research is needed to confirm the association between metal exposures and telomere length.