Publications

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.

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.

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.

BACKGROUND: Arsenic exposure has been related to numerous adverse cardiovascular outcomes. The aim of this study was to investigate the cross-sectional and prospective association between arsenic exposure with echocardiographic measures of left ventricular (LV) geometry and functioning.

METHODS: A total of 1337 young adult participants free of diabetes mellitus and cardiovascular disease were recruited from the SHFS (Strong Heart Family Study). The sum of inorganic and methylated arsenic concentrations in urine (ΣAs) at baseline was used as a biomarker of arsenic exposure. LV geometry and functioning were assessed using transthoracic echocardiography at baseline and follow-up.

RESULTS: Mean follow-up was 5.6 years, and median (interquartile range) of ΣAs was 4.2 (2.8-6.9) µg/g creatinine. Increased arsenic exposure was associated with prevalent LV hypertrophy, with an odds ratio (95% CI) per a 2-fold increase in ΣAs of 1.47 (1.05-2.08) in all participants and of 1.58 (1.04-2.41) among prehypertensive or hypertensive individuals. Measures of LV geometry, including LV mass index, left atrial systolic diameter, interventricular septum, and LV posterior wall thickness, were positively and significantly related to arsenic exposure. Among measures of LV functioning, stroke volume, and ejection fraction were associated with arsenic exposure.

CONCLUSIONS: Arsenic exposure was related to an increase in LV wall thickness and LV hypertrophy in young American Indians with a low burden of cardiovascular risk factors. The relationship was stronger in participants with prehypertension or hypertension, suggesting that potential cardiotoxic effects of arsenic might be more pronounced in individuals already undergoing cardiovascular adaptive mechanisms following elevated systemic blood pressure.

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.