Publications

BACKGROUND: Common genetic variation in the arsenic methyltransferase (AS3MT) gene region is known to be associated with arsenic metabolism efficiency (AME), measured as the percentage of dimethylarsinic acid (DMA%) in the urine. Rare, protein-altering variants in AS3MT could have even larger effects on AME, but their contribution to AME has not been investigated.

OBJECTIVES: We estimated the impact of rare, protein-coding variation in AS3MT on AME using a multi-population approach to facilitate the discovery of population-specific and shared causal rare variants.

METHODS: We generated targeted DNA sequencing data for the coding regions of AS3MT for three arsenic-exposed cohorts with existing data on arsenic species measured in urine: Health Effects of Arsenic Longitudinal Study (HEALS, n=2,434), Strong Heart Study (SHS, n=868), and New Hampshire Skin Cancer Study (NHSCS, n=666). We assessed the collective effects of rare (allele frequency <1%), protein-altering AS3MT variants on DMA%, using multiple approaches, including a test of the association between rare allele carrier status (yes/no) and DMA% using linear regression (adjusted for common variants in 10q24.32 region, age, sex, and population structure).

RESULTS: We identified 23 carriers of rare-protein-altering AS3MT variant across all cohorts (13 in HEALS and 5 in both SHS and NHSCS), including 6 carriers of predicted loss-of-function variants. DMA% was 6-10% lower in carriers compared with noncarriers in HEALS [β=-9.4 (95% CI: -13.9, -4.8)], SHS [β=-6.9 (95% CI: -13.6, -0.2)], and NHSCS [β=-8.7 (95% CI: -15.6, -2.2)]. In meta-analyses across cohorts, DMA% was 8.7% lower in carriers [β=-8.7 (95% CI: -11.9, -5.4)].

DISCUSSION: Rare, protein-altering variants in AS3MT were associated with lower mean DMA%, an indicator of reduced AME. Although a small percentage of the population (0.5-0.7%) carry these variants, they are associated with a 6-10% decrease in DMA% that is consistent across multiple ancestral and environmental backgrounds. https://doi.org/10.1289/EHP8152.

BACKGROUND: Whole-exome sequencing (WES) can expedite research on genetic variation in non-human primate (NHP) models of human diseases. However, NHP-specific reagents for exome capture are not available. This study reports the use of human-specific capture reagents in WES for olive baboons, marmosets, and vervet monkeys.

METHODS: Exome capture was carried out using the SureSelect Human All Exon V6 panel from Agilent Technologies, followed by high-throughput sequencing. Capture of protein-coding genes and detection of single nucleotide variants were evaluated.

RESULTS: Exome capture and sequencing results showed that more than 97% of old world and 93% of new world monkey protein coding genes were detected. Single nucleotide variants were detected across the genomes and missense variants were found in genes associated with human diseases.

CONCLUSIONS: A cost-effective approach based on commercial, human-specific reagents can be used to perform WES for the discovery of genetic variants in these NHP species.

Genetic effects are suspected to influence cadmium internal dose. Our objective was to assess genetic determinants of urine cadmium in American Indian adults participating in the Strong Heart Family Study (SHFS). Urine cadmium levels and genotyped short tandem repeat (STR) markers were available on 1936 SHFS participants. We investigated heritability, including gene-by-sex and smoking interactions, and STR-based quantitative trait locus (QTL) linkage, using a variance-component decomposition approach, which incorporates the genetic information contained in the pedigrees. We also used available single nucleotide polymorphisms (SNPs) from Illumina's Metabochip and custom panel to assess whether promising QTLs associated regions could be attributed to SNPs annotated to specific genes. Median urine cadmium levels were 0.44 μg/g creatinine. The heritability of urine cadmium concentrations was 28%, with no evidence of gene-by-sex or -smoking interaction. We found strong statistical evidence for a genetic locus at chromosome 16 determining urine cadmium concentrations (Logarithm of odds score [LOD] = 3.8). Among the top 20 associated SNPs in this locus, 17 were annotated to ABCC1 (p-values from 0.0002 to 0.02), and attenuated the maximum linkage peak by a ∼40%. Suggestive QTL signals (LOD>1.9) in chromosomes 2, 6, 11, 14, and 19, showed associated SNPs in the genes NDUFA10, PDE10A, PLEKHA7, BAZ1A and CHAF1A, respectively. Our findings support that urinary cadmium levels are heritable and influenced by a QTL on chromosome 16, which was explained by genetic variation in ABCC1. Studies with extended sets of genome-wide markers are needed to confirm these findings and to identify additional metabolism and toxicity pathways for cadmium.

BACKGROUND: Epidemiologic studies often use self-report as proxy for clinical history. However, whether self-report correctly identifies prevalence in minority populations with health disparities and poor health-care access is unknown. Furthermore, overlap of clinical vascular events with covert vascular brain injury (VBI), detected by imaging, is largely unexamined.

METHODS: The Strong Heart Study recruited American Indians from 3 regions, with surveillance and adjudication of stroke events from 1989 to 2013. In 2010-2013, all 817 survivors, aged 65-95 years, underwent brain imaging, neurological history interview, and cognitive testing. VBI was defined as imaged infarct or hemorrhage.

RESULTS: Adjudicated stroke was prevalent in 4% of participants and separately collected, self-reported stroke in 8%. Imaging-defined VBI was detected in 51% and not associated with any stroke event in 47%. Compared with adjudication, self-report had 76% sensitivity and 95% specificity. Participants with adjudicated or self-reported stroke had the poorest performance on cognitive testing; those with imaging-only (covert) VBI had intermediate performance.

CONCLUSION: In this community-based cohort, self-report for prior stroke had good performance metrics. A majority of participants with VBI did not have overt, clinically recognized events but did have neurological or cognitive symptoms. Data collection methodology for studies in a resource-limited setting must balance practical limitations in costs, accuracy, feasibility, and research goals.

BACKGROUND: Cancer mortality among American Indian (AI) people varies widely, but factors associated with cancer mortality are infrequently assessed.

METHODS: Cancer deaths were identified from death certificate data for 3516 participants of the Strong Heart Study, a population-based cohort study of AI adults ages 45-74 years in Arizona, Oklahoma, and North and South Dakota. Cancer mortality was calculated by age, sex and region. Cox proportional hazards model was used to assess independent associations between baseline factors in 1989 and cancer death by 2010.

RESULTS: After a median follow-up of 15.3 years, the cancer death rate per 1000 person-years was 6.33 (95 % CI 5.67-7.04). Cancer mortality was highest among men in North/South Dakota (8.18; 95 % CI 6.46-10.23) and lowest among women in Arizona (4.57; 95 % CI 2.87-6.92). Factors independently associated with increased cancer mortality included age, current or former smoking, waist circumference, albuminuria, urinary cadmium, and prior cancer history. Factors associated with decreased cancer mortality included Oklahoma compared to Dakota residence, higher body mass index and total cholesterol. Sex was not associated with cancer mortality. Lung cancer was the leading cause of cancer mortality overall (1.56/1000 person-years), but no lung cancer deaths occurred among Arizona participants. Mortality from unspecified cancer was relatively high (0.48/100 person-years; 95 % CI 0.32-0.71).

CONCLUSIONS: Regional variation in AI cancer mortality persisted despite adjustment for individual risk factors. Mortality from unspecified cancer was high. Better understanding of regional differences in cancer mortality, and better classification of cancer deaths, will help healthcare programs address cancer in AI communities.

IMPORTANCE: American Indian communities experience a high burden of coronary heart disease (CHD). Strategies are needed to identify individuals at risk and implement preventive interventions.

OBJECTIVE: To investigate the association of blood DNA methylation (DNAm) with incident CHD using a large number of methylation sites (cytosine-phosphate-guanine [CpG]) in a single model.

DESIGN, SETTING, AND PARTICIPANTS: This prospective study, including a discovery cohort (the Strong Heart Study [SHS]) and 4 additional cohorts (the Women's Health Initiative [WHI], the Framingham Heart Study [FHS], the Atherosclerosis Risk in Communities Study ([ARIC]-Black, and ARIC-White), evaluated 12 American Indian communities in 4 US states; African American women, Hispanic women, and White women throughout the US; White men and White women from Massachusetts; and Black men and women and White men and women from 4 US communities. A total of 2321 men and women (mean [SD] follow-up, 19.1 [9.2] years) were included in the SHS, 1874 women (mean [SD] follow-up, 15.8 [5.9] years) in the WHI, 2128 men and women (mean [SD] follow-up, 7.7 [1.8] years) in the FHS, 2114 men and women (mean [SD] follow-up, 20.9 [7.2] years) in the ARIC-Black, and 931 men and women (mean [SD] follow-up, 20.9 [7.2] years) in the ARIC-White. Data were collected from May 1989 to December 2018 and analyzed from February 2019 to May 2021.

EXPOSURE: Blood DNA methylation.

MAIN OUTCOME AND MEASURE: Using a high-dimensional time-to-event elastic-net model for the association of 407 224 CpG sites with incident CHD in the SHS (749 events), this study selected the differentially methylated CpG positions (DMPs) selected in the SHS and evaluated them in the WHI (531 events), FHS (143 events), ARIC-Black (350 events), and ARIC-White (121 events) cohorts.

RESULTS: The median (IQR) age of participants in SHS was 55 (49-62) years, and 1359 participants (58.6%) were women. Elastic-net models selected 505 DMPs associated with incident CHD in the SHS beyond established risk factors, center, blood cell counts, and genetic principal components. Among those DMPs, 33 were commonly selected in 3 or 4 of the other cohorts and the pooled hazard ratios from the standard Cox models were significant at P < .05 for 10 of the DMPs. For example, the hazard ratio (95% CI) for CHD comparing the 90th and 10th percentiles of differentially methylated CpGs was 0.86 (0.78-0.95) for cg16604233 (tagged to COL11A2) and 1.23 (1.08-1.39) for cg09926486 (tagged to FRMD5). Some of the DMPs were consistent in the direction of the association; others showed associations in opposite directions across cohorts. Untargeted independent elastic-net models of CHD showed distinct DMPs, genes, and network of genes in the 5 cohorts.

CONCLUSIONS AND RELEVANCE: In this multi-cohort study, blood-based DNAm findings supported an association between a complex blood epigenomic signature and CHD that was largely different across populations.

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.