Determinants of telomere length across human tissues

Telomere shortening is a hallmark of aging. Telomere length (TL) in blood cells has been studied extensively as a biomarker of human aging and disease; however, little is known regarding variability in TL in nonblood, disease-relevant tissue types. Here, we characterize variability in TLs from 6391 tissue samples, representing >20 tissue types and 952 individuals from the Genotype-Tissue Expression (GTEx) project. We describe differences across tissue types, positive correlation among tissue types, and associations with age and ancestry. We show that genetic variation affects TL in multiple tissue types and that TL may mediate the effect of age on gene expression. Our results provide the foundational knowledge regarding TL in healthy tissues that is needed to interpret epidemiological studies of TL and human health.

Effect of core strength on the measure of power in the extremities

The purpose of this study was to (a) develop a functional field test to assess the role of the core musculature and its impact on sport performance in an athletic population and (b) develop a functional field test to determine how well the core can transfer forces from the lower to the upper extremities. Twenty-five DI collegiate football players performed medicine ball throws (forward, reverse, right, and left) in static and dynamic positions. The results of the medicine ball throws were compared with several athletic performance measurements: 1 repetition maximum (1RM) squat, squat kg/bw, 1RM bench press, bench kg/bw, countermovement vertical jump (CMJ), 40-yd dash (40 yd), and proagility (PrA). Push press power (PWR) was used to measure the transfer of forces through the body. Several correlations were found in both the static and dynamic medicine ball throws when compared with the performance measures. Static reverse correlated with CMJ (r = 0.44), 40 yd (r = 0.5), and PrA (r = 0.46). Static left correlated with bench kg/bw (0.42), CMJ (0.44), 40 yd (0.62), and PrA (0.59). Static right also correlated with bench kg/bw (0.41), 40 yd (0.44), and PrA (0.65). Dynamic forward (DyFw) correlated with the 1RM squat (r = 0.45) and 1RM bench (0.41). Dynamic left and Dynamic right correlated with CMJ, r = 0.48 and r = 0.40, respectively. Push press power correlated with bench kg/bw (0.50), CMJ (0.48), and PrA (0.48). A stepwise regression for PWR prediction identified 1RM squat as the best predictor. The results indicate that core strength does have a significant effect on an athlete's ability to create and transfer forces to the extremities. Currently, plank exercises are considered an adequate method of training the core for athletes to improve core strength and stability. This is a problem because it puts the athletes in a nonfunctional static position that is very rarely replicated in the demands of sport-related activities. The core is the center of most kinetic chains in the body and should be trained accordingly.

Association of Arsenic Exposure with Whole Blood DNA Methylation: An Epigenome-Wide Study of Bangladeshi Adults

BACKGROUND

Arsenic exposure affects [Formula: see text] people worldwide, including [Formula: see text] in Bangladesh. Arsenic exposure increases the risk of cancer and other chronic diseases, and one potential mechanism of arsenic toxicity is epigenetic dysregulation.

OBJECTIVE

We assessed associations between arsenic exposure and genome-wide DNA methylation measured at baseline among 396 Bangladeshi adults participating in the Health Effects of Arsenic Longitudinal Study (HEALS) who were exposed by drinking naturally contaminated well water.

METHODS

Methylation in whole blood DNA was measured at [Formula: see text] using the Illumina InfiniumMethylationEPIC (EPIC) array. To assess associations between arsenic exposure and CpG methylation, we used linear regression models adjusted for covariates and surrogate variables (SVs) (capturing unknown technical and biologic factors). We attempted replication and conducted a meta-analysis using an independent dataset of [Formula: see text] from 400 Bangladeshi individuals with arsenical skin lesions.

RESULTS

We identified 34 CpGs associated with [Formula: see text] creatinine-adjusted urinary arsenic [[Formula: see text]]. Sixteen of these CpGs annotated to the [Formula: see text] array, and 10 associations were replicated ([Formula: see text]). The top two CpGs annotated upstream of the ABR gene (cg01912040, cg10003262 ). All urinary arsenic-associated CpGs were also associated with arsenic concentration measured in drinking water ([Formula: see text]). Meta-analysis ([Formula: see text] samples) identified 221 urinary arsenic-associated CpGs ([Formula: see text]). The arsenic-associated CpGs from the meta-analysis were enriched in non-CpG islands and shores ([Formula: see text]) and depleted in promoter regions ([Formula: see text]). Among the arsenic-associated CpGs ([Formula: see text]), we observed significant enrichment of genes annotating to the reactive oxygen species pathway, inflammatory response, and tumor necrosis factor [Formula: see text] ([Formula: see text]) signaling via nuclear factor kappa-B ([Formula: see text]) hallmarks ([Formula: see text]).

CONCLUSIONS

The novel and replicable associations between arsenic exposure and DNA methylation at specific CpGs observed in this work suggest that epigenetic alterations should be further investigated as potential mediators in arsenic toxicity and as biomarkers of exposure and effect in exposed populations. https://doi.org/10.1289/EHP3849.

Genome-wide association study of telomere length among South Asians identifies a second RTEL1 association signal

Background

Leucocyte telomere length (TL) is a potential biomarker of ageing and risk for age-related disease. Leucocyte TL is heritable and shows substantial differences by race/ethnicity. Recent genome-wide association studies (GWAS) report ~10 loci harbouring SNPs associated with leucocyte TL, but these studies focus primarily on populations of European ancestry.

Objective

This study aims to enhance our understanding of genetic determinants of TL across populations.

Methods

We performed a GWAS of TL using data on 5075 Bangladeshi adults. We measured TL using one of two technologies (qPCR or a Luminex-based method) and used standardised variables as TL phenotypes.

Results

Our results replicate previously reported associations in the TERC and TERT regions (P=2.2×10⁻⁸ and P=6.4×10⁻⁶, respectively). We observed a novel association signal in the RTEL1 gene (intronic SNP rs2297439; P=2.82×10⁻⁷) that is independent of previously reported TL-associated SNPs in this region. The minor allele for rs2297439 is common in South Asian populations (≥0.25) but at lower frequencies in other populations (eg, 0.07 in Northern Europeans). Among the eight other previously reported association signals, all were directionally consistent with our study, but only rs8105767 (ZNF208) was nominally significant (P=0.003). SNP-based heritability estimates were as high as 44% when analysing close relatives but much lower when analysing distant relatives only.

Conclusions

In this first GWAS of TL in a South Asian population, we replicate some, but not all, of the loci reported in prior GWAS of individuals of European ancestry, and we identify a novel second association signal at the RTEL1 locus.

A study of telomere length, arsenic exposure, and arsenic toxicity in a Bangladeshi cohort

BACKGROUND

Chronic arsenic exposure is associated with increased risk for arsenical skin lesions, cancer, and other adverse health outcomes. One potential mechanism of arsenic toxicity is telomere dysfunction. However, prior epidemiological studies of arsenic exposure, telomere length (TL), and skin lesion are small and cross-sectional. We investigated the associations between arsenic exposure and TL and between baseline TL and incident skin lesion risk among individuals participating in the Health Effects of Arsenic Longitudinal Study in Bangladesh (2000-2009).

METHODS

Quantitative PCR was used to measure the average TL of peripheral blood DNA collected at baseline. The association between baseline arsenic exposure (well water and urine) and TL was estimated in a randomly-selected subcohort (n = 1469). A nested case-control study (466 cases and 464 age- and sex-matched controls) was used to estimate the association between baseline TL and incident skin lesion risk (diagnosed < 8 years after baseline).

RESULTS

No association was observed between arsenic exposure (water or urine) and TL. Among incident skin lesion cases and matched controls, we observed higher skin lesion risk among individuals with shorter TL (P_trend = 1.5 × 10^-5) with odds ratios of 2.60, 1.59, and 1.10 for the first (shortest), second, and third TL quartiles compared to the fourth (longest).

CONCLUSIONS

Arsenic exposure was not associated with TL among Bangladeshi adults, suggesting that leukocyte TL may not reflect a primary mode of action for arsenic's toxicity. However, short TL was associated with increased skin lesion risk, and may be a biomarker of arsenic susceptibility modifying arsenic's effect on skin lesion risk.

The contribution of parent-to-offspring transmission of telomeres to the heritability of telomere length in humans

Leukocyte telomere length (LTL) is a heritable trait with two potential sources of heritability (h²): inherited variation in non-telomeric regions (e.g., SNPs that influence telomere maintenance) and variability in the lengths of telomeres in gametes that produce offspring zygotes (i.e., "direct" inheritance). Prior studies of LTL h² have not attempted to disentangle these two sources. Here, we use a novel approach for detecting the direct inheritance of telomeres by studying the association between identity-by-descent (IBD) sharing at chromosome ends and phenotypic similarity in LTL. We measured genome-wide SNPs and LTL for a sample of 5069 Bangladeshi adults with substantial relatedness. For each of the 6318 relative pairs identified, we used SNPs near the telomeres to estimate the number of chromosome ends shared IBD, a proxy for the number of telomeres shared IBD (T_shared). We then estimated the association between T_shared and the squared pairwise difference in LTL ((ΔLTL)²) within various classes of relatives (siblings, avuncular, cousins, and distant), adjusting for overall genetic relatedness (ϕ). The association between T_shared and (ΔLTL)² was inverse among all relative pair types. In a meta-analysis including all relative pairs (ϕ > 0.05), the association between T_shared and (ΔLTL)² (P = 0.01) was stronger than the association between ϕ and (ΔLTL)² (P = 0.43). Our results provide strong evidence that telomere length (TL) in parental germ cells impacts TL in offspring cells and contributes to LTL h² despite telomere "reprogramming" during embryonic development. Applying our method to larger studies will enable robust estimation of LTL h² attributable to direct transmission of telomeres.

Treatment of liver injuries at level I and level II centers in a multi-institutional metropolitan trauma system. The Midwest Trauma Society Liver Trauma Study Group

OBJECTIVE

The development of trauma systems and trauma centers has had a major impact on the fate of the critically injured patient. However, some have suggested that care may be compromised if too many trauma centers are designated for a given area. As of 1987, the state of Missouri had designated six adult trauma centers, two Level I and four Level II, for the metropolitan Kansas City, Mo, area, serving a population of approximately 1 million people. To determine whether care was comparable between the Level I and II centers, we conducted a concurrent evaluation of the fate of patients with a sentinel injury, hepatic trauma, over a 6-year period (1987-1992) who were treated at these six trauma centers.

METHODS

All patients during the 6-year study period who suffered liver trauma and who survived long enough to be evaluated by computerized tomography or celiotomy were entered into the study. Patients with central nervous system trauma were excluded from analysis. Information concerning mechanism of injury, RTS, Injury Severity Score (ISS), presence of shock, liver injury scoring, mode of treatment, mortality, and length of stay were recorded on abstract forms for analysis. Care was evaluated by mortality, time to the operating room (OR), and intensive care unit (ICU) and hospital length of stay.

RESULTS

Over the 6-year period 300 patients with non-central nervous system liver trauma were seen. Level I centers cared for 195 patients and Level II centers cared for 105. There was no difference in mean ISS or ISS > 25 between Level I and II centers. Fifty-five (28%) patients arrived in shock at Level I centers and 24 (23%) at Level II centers. Forty-eight patients (16%) died. Thirty-two (16%) died at Level I centers, and 16 (15%) died at Level II centers. Twenty of 55 patients (36%) in shock died at Level I centers, and 11 of 24 (46%) died at Level II centers (p = 0.428). Forty-three patients (22%) had liver scaling scores of IV-VI at Level I centers, and 10 (10%) had similar scores at Level II centers (p < 0.01). With liver scores IV-VI, 22 of 43 (51%) died at Level I centers and 10 of 14 (71%) died at Level II centers (p = 0.184). There was no difference in mean time or in delays beyond 1 hour to the OR for those patients in shock between Level I and II centers. There was a longer ICU stay at Level II centers (5.0 +/- 8.3 vs. 2.8 +/- 8.4 days, p = 0.04). This difference was confined to penetrating injuries. There was no difference in hospital length of stay.

CONCLUSIONS

In a metropolitan trauma system, when Level I and II centers were compared for their ability to care for victims of hepatic trauma, there was no discernible difference in care rendered with respect to severity of injury, mortality, delays to the OR, or hospital length of stay. It was observed that more severe liver injuries were seen at Level I centers, but this did not seem to significantly affect care at Level II centers. There was a longer ICU stay observed at Level II centers owing to penetrating injuries, possibly because there were fewer penetrating injuries treated at these facilities. Although the bulk of patients were seen at Level I centers, care throughout the system was equivalent.

The association between telomere length and mortality in Bangladesh

Telomeres are tandem repeat sequences at the end of chromosomes that bind proteins to protect chromosome ends. Telomeres shorten with age, and shorter leukocyte telomere length (TL) has been associated with overall mortality in numerous studies. However, this association has not been tested in populations outside of Europe and the U.S. We assessed the association between TL and subsequent mortality using data on 744 mortality cases and 761 age-/sex-matched controls sampled from >27,000 participants from three longitudinal Bangladeshi cohorts: Health Effects of Arsenic Longitudinal Study (HEALS), HEALS Expansion (HEALS-E), and Bangladesh Vitamin E and Selenium Trial (BEST). We used conditional logistic regression to estimate odds ratios (ORs) for the association between a standardized TL variable and overall mortality, as well as mortality from chronic diseases, respiratory diseases, circulatory diseases, and cancer. In HEALS and BEST, we observed an association between shorter TL and increased overall mortality (P=0.03 and P=0.03), mortality from chronic disease (P=0.01 and P=0.03) and mortality from circulatory disease (P=0.03 and P=0.04). Results from pooled analyses of all cohorts were consistent with HEALS and BEST. This is the first study demonstrating an association between short TL and increased mortality in a population of non-European ancestry.

Association between genome-wide copy number variation and arsenic-induced skin lesions: a prospective study

BACKGROUND

Exposure to arsenic in drinking water is a global health problem and arsenic-induced skin lesions are hallmark of chronic arsenic toxicity. We and others have reported germline genetic variations as risk factors for such skin lesions. The role of copy number variation (CNV) in the germline DNA in this regard is unknown.

METHODS

From a large prospectively followed-up cohort, exposed to arsenic, we randomly selected 2171 subjects without arsenic-induced skin lesions at enrollment and genotyped their whole blood DNA samples on Illumina Cyto12v2.1 SNP chips to generate DNA copy number. Participants were followed up every 2 years for a total of 8 years, especially for the development of skin lesions. In Cox regression models, each CNV segment was used as a predictor, accounting for other potential covariates, for incidence of skin lesions.

RESULT

The presence of genomic deletion(s) in a number of genes (OR5J2, GOLGA6L7P, APBA2, GALNTL5, VN1R31P, PHKG1P2, SGCZ, ZNF658) and lincRNA genes (RP11-76I14.1, CTC-535 M15.2, RP11-73B2.2) were associated with higher risk [HR between 1.67 (CI 1.3-2.1) and 2.15 (CI 1.5-2.9) for different CNVs] for development of skin lesions independent of gender, age, and arsenic exposure. Some deletions had stronger effect in a specific gender (ZNF658 in males, SGCZ in females) and some had stronger effect in higher arsenic exposure (lincRNA CTD-3179P9.1) suggesting a possible gene-environment interaction.

CONCLUSION

This first genome-wide CNV study in a prospectively followed-up large cohort, exposed to arsenic, suggests that DNA deletion in several genes and lincRNA genes may predispose an individual to a higher risk of development of arsenic-induced skin lesions.

A novel pooled-sample multiplex luminex assay for high-throughput measurement of relative telomere length

OBJECTIVES

Relative telomere length (RTL) is a potential biomarker of aging and risk for chronic disease. Previously, we developed a probe-based RTL assay on Luminex platform, where probes for Telomere (T) and reference gene (R) for a given DNA sample were tested in a single well. Here, we describe a method of pooling multiple samples in one well to increase the throughput and cost-effectiveness.

METHODS

We used four different microbeads for the same T-probe and four different microbeads for the same R-probe. Each pair of probe sets were hybridized to DNA in separate plates and then pooled in a single plate for all the subsequent steps. We used DNA samples from 60 independent individuals and repeated in multiple batches to test the precision.

RESULTS

The precision was good to excellent with Intraclass correlation coefficient (ICC) of 0.908 (95% CI 0.856-0.942). More than 67% of the variation in the RTL could be explained by sample-to-sample variation; less than 0.1% variation was due to batch-to-batch variation and 0.3% variation was explained by bead-to-bead variation. We increased the throughput of RTL Luminex assay from 60 to 240 samples per run. The new assay was validated against the original Luminex assay without pooling (r = 0.79, P = 1.44 × 10-15 ). In an independent set of samples (n = 550), the new assay showed a negative correlation of RTL with age (r = -0.41), a result providing external validation for the method.

CONCLUSION

We describe a novel high throughput pooled-sample multiplex Luminex assay for RTL with good to excellent precision suitable for large-scale studies.

Rare, Protein-Altering Variants in AS3MT and Arsenic Metabolism Efficiency: A Multi-Population Association Study

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.

Assessing the impact of arsenic metabolism efficiency on DNA methylation using Mendelian randomization

BACKGROUND

Arsenic exposure affects >100 million people globally and increases risk for chronic diseases. One possible toxicity mechanism is epigenetic modification. Previous epigenome-wide association studies (EWAS) have identified associations between arsenic exposure and CpG-specific DNA methylation. To provide additional evidence that observed associations represent causal relationships, we examine the association between genetic determinants of arsenic metabolism efficiency (percent dimethylarsinic acid, DMA%, in urine) and DNA methylation among individuals from the Health Effects of Arsenic Longitudinal Study (n = 379) and Bangladesh Vitamin E and Selenium Trial (n = 393).

METHODS

We used multivariate linear models to assess the association of methylation at 221 arsenic-associated CpGs with DMA% and measures of genetically predicted DMA% derived from three SNPs (rs9527, rs11191527, and rs61735836). We also conducted two-sample Mendelian randomization analyses to estimate the association between arsenic metabolism efficiency and CpG methylation.

RESULTS

Among the associations between DMA% and methylation at each of 221 CpGs, 64% were directionally consistent with associations observed between arsenic exposure and the 221 CpGs from a prior EWAS. Similarly, among the associations between genetically predicted DMA% and each CpG, 62% were directionally consistent with the prior EWAS results. Two-sample Mendelian randomization analyses produced similar conclusions.

CONCLUSION

Our findings support the hypothesis that arsenic exposure effects DNA methylation at specific CpGs in whole blood. Our novel approach for assessing the impact of arsenic exposure on DNA methylation requires larger samples in order to draw more robust conclusions for specific CpG sites.

The ventral hippocampus and nucleus accumbens as neural substrates for cocaine contextual memory reconsolidation

Drug craving triggered by cues that were once associated with drug intoxication is a major contributor to continued drug-seeking behaviors. Addictive drugs engage molecular pathways of associative learning and memory. Reactivated memories are vulnerable to disruption by interference with the process of reconsolidation, hence targeting reconsolidation could be a strategy to reduce cue-induced drug craving and relapse. Here we examined the circuitry of cocaine contextual memory reconsolidation and explored neuroplasticity following memory reactivation. Mice underwent chemogenetic inhibition of either nucleus accumbens (NA) neurons or the glutamatergic projection neurons from the ventral hippocampus (vHPC) to NA using inhibitory designer receptors exclusively activated by designer drugs (iDREADD). Mice underwent cocaine conditioned place preference followed by reactivation of the cocaine contextual memory. Clozapine-N-oxide (CNO) was administered after memory reactivation to inhibit either NA neurons or the accumbens-projecting vHPC neurons during the reconsolidation period. When retested 3 days later, a significant reduction in the previously established preference for the cocaine context was found in both conditions. FosTRAP2-Ai14 mice were used to identify neurons activated by cocaine memory recall and to evaluate plasticity in NA medium spiny neurons (MSNs) and vHPC pyramidal neurons upon recall of cocaine memories. Results indicate a significant increase in dendritic spine density in NA MSNs activated by cocaine memory recall, particularly of the thin spine type. Sholl analysis indicated longer dendritic length and more branching of NA MSNs after cocaine memory recall than without memory reactivation. vHPC neurons showed increased spine density, with the most robust change in stubby spines. These results implicate a circuit involving glutamatergic projections from the vHPC onto NA neurons which is necessary for the reconsolidation of cocaine memories. Interruption of cocaine memory reconsolidation reduced drug-seeking behavior.