Shared loci of telomere length with brain morphology, and pleiotropy in transcriptomic and epigenomic profiles of brain

Gita A Pathak, Frank R Wendt, Daniel F Levey, Adam P Mecca, Christopher H van Dyck, Joel Gelernter, Renato Polimanti
doi: https://doi.org/10.1101/2020.08.03.20167726

ABSTRACT

Premature shortening of telomere length is observed in neuropsychiatric disorders. We tested genetic colocalization of seven and nine leukocyte telomere length (LTL) loci in two ancestry groups, European (EUR) and East-Asian (EAS), respectively (total n=60,601) with brain morphology measures for 101 region-of-interests (ROI) (n=21,821). The posterior probability (>90%) was observed for ‘fourth ventricle’, ‘gray matter’ and ‘cerebellar vermal lobules I-IV’ volumes. We found regulatory genes (p ≤2.47 x10−6) by integrating transcriptomic (EAS=4;EUR=5) and methylation data (EUR=17; EAS=4) of brain tissues using Summary-based Mendelian Randomization (SMR). The LTL SNP associations were tested for brain-based chromatin profiles using H-MAGMA (EUR=50; EAS=97; p ≤ 1.02 x10−6). Pathway enrichment of tissue-specific genes highlighted calcium ion transport (fetal brain) and G2/M cell cycle transition (adult brain). This study provides evidence that previously reported LTL associations with neuropsychiatric disorders could be related to a shared genetic relationship between LTL and brain structural and regulatory traits.

1 Introduction

Telomeres are short sequences of “TTAGGG” at the 3’ end of the chromosomes. With replicating cell cycles, telomere length continues to shorten and represents the mitotic history of the cell. The telomere length shortens (i.e., the number of TTAGGG repeat motifs decreases) over chronological age (i.e., actual age of the person) (Koliada et al., 2015). Additionally, under the effects of influences such as, inflammation, oxidative stress, elevated hormone levels under stress or pathological conditions, telomere length can shorten prematurely (Vakonaki et al., 2018). Perturbations in telomere lengthening enzymes such as telomerase reverse transcriptase (TERT) and telomerase RNA component result in cellular damage and genomic instability, further accelerating damage to physiological processes termed cellular aging (Bär and Blasco, 2016).

Psychiatric traits such as major depressive disorder (MDD), schizophrenia, anxiety, bipolar disorder (Muneer and Minhas, 2019), and neurodegenerative disorders (e.g., Alzheimer’s disease) have documented associations with shortened leukocyte telomere length (LTL)(Forero et al., 2016). Childhood anxiety and adversity contribute to lifetime diagnosis of mood disorders(McLaughlin et al., 2012) such as MDD(Kaczmarczyk et al., 2018; Wang et al., 2017), bipolar disorder(Aas et al., 2016), schizophrenia(Aas et al., 2019), and posttraumatic stress disorder(McLaughlin et al., 2017). Childhood adversities and chronic illness are correlated with shorter leukocyte telomere length in individuals with anxiety diagnosis (Kananen et al., 2010). In addition to brain-based disorders, LTL is associated with brain morphology of such structures as hippocampal, temporal, and cingulate region volumes among others (King et al., 2014). Teenagers with MDD were reported to have shorter telomeres and smaller right hippocampal volume (Henje Blom et al., 2015). A longitudinal study measuring the effects of various mental training modules (to improve cognition and behavioral traits) reported a significant association between reduced telomere length and thinning of the cortex structure highlighting neuronal plasticity variation between telomere and certain brain structures(Puhlmann et al., 2019). The association of telomere length variability with psychiatric disorders has led to the hypothesis of accelerated cellular or biological aging mediating physiological deterioration (Aas et al., 2019).

Telomerase mediates cell differentiation and is highly expressed in neural stem and progenitor cells, and declines with increasing age (Liu et al., 2018). Dysfunction in telomerase impairs neural circuitry in the hippocampus, which is one of many factors that leads to memory decline (Zhou et al., 2017). The expression of TERT is lower in astrocytes than oligodendrocytes (Szebeni et al., 2014), and telomere attrition induces cellular senescence of astrocytes leading to neurological dysregulation (Cohen and Torres, 2019). Telomere shortening disrupts the neuronal differentiation cell cycle (Ferrón et al., 2009).Neurogenesis in both adult and fetal tissue is influenced by functioning telomerase and interaction with sheltrin complex which protects the telomeric caps (i.e., prevents reduction in telomeric TTAGGG repeat number) (Lobanova et al., 2017). Even though telomere length is generally expected to be associated with age-related disorders, telomere length attrition can be influenced by maternal and intrauterine stressors during fetal development exhibiting its effects during an individual’s life course (Edlow et al., 2019).

Twin-based genetic heritability of LTL is estimated at around 30-60% (Dorajoo et al., 2019). Genome-wide association studies (GWAS) measure relative differences in frequencies of genome-wide single nucleotide polymorphisms (SNPs) for a trait of interest. A relatively high proportion of SNP associations are found in non-coding regions and may have tissue specific regulatory effects in local or distant regions (Giral et al., 2018). Tissue- and cell-type specific regulatory effects may be investigated by comparing the association of SNPs with gene expression (termed expression quantitative trait loci, eQTL) and methylation levels of CpG (cytosine-phosphate-guanosine) sites termed methylation QTL (mQTL) (T. Zhao et al., 2019). Leveraging large studies with genetic, transcriptomic, and epigenetic data allows us to identify putative tissue and cell-specific regulatory effects of trait-specific SNPs. To our knowledge, no study has investigated shared causal loci for brain morphology and brain tissue based multi-omic regulatory consequences of telomere length SNP associations.

In this study we used large scale LTL GWAS of 37,505 European and 23,096 Chinese individuals (ntotal=60,601) (Dorajoo et al., 2019) to identify a) LTL causal loci shared with brain volume measures, b) transcriptomic and epigenetic regulatory effects, and c) chromatin profiles in brain tissue (Figure 1). First we analyzed whether LTL risk loci are shared with 101 T1-MRI (magnetic resonance imaging weighted for longitudinal relaxation time) brain region-of-interest phenotypes (n=21,821) (B. Zhao et al., 2019).Then, we integrated meta-analyzed mQTL (Qi et al., 2018) (n = 1,160) and eQTL associations from different regions of the brain (n = 1,194) obtained from Genotype-Tissue Expression GTEx (GTEx Consortium et al., 2017), CommonMind Consortium (CMC) (Fromer et al., 2016), the Religious Orders Study and the Memory and Aging Project (ROSMAP) (Ng et al., 2017) and prefrontal cortex-based eQTL (n = 1,387) of PsychENCODE consortia (Gandal et al., 2018) were used to conduct gene based association study for LTL using Summary-based Mendelian Randomization (Zhu et al., 2016). Chromatin profiles from human brain tissue of two developmental stages – fetal and adult, and two cell types – astrocytes and neurons were investigated to identify genes specific to tissue and cell types using H-MAGMA (Hi-C-coupled multimarker analysis of genomic annotation) (Sey et al., 2020).

Figure 1:
Figure 1: Study Overview:

Schematic presenting the hypothesis and four major domains studied – i) shared genomic regions between leukocyte telomere length (LTL) and brain volume measures, and gene-based associations for LTL in brain tissue-based ii) transcriptomics – combined brain tissue (meta-analyzed from different regions of the brains) and prefrontal cortex, iii) methylomics – combined brain tissues (meta-analyzed), and iv) Hi-C chromatin profiles – adult and fetal brain, astrocytes and neurons.

2 Results

2.1 Genetic colocalization with brain morphology

We performed colocalization analysis using coloc (Giambartolomei et al., 2014) between LTL and 101 brain volume measures to identify causal genomic loci that with respect to telomere length and measures of brain morphology. The coloc method measures the posterior probability for five hypotheses (H0-H4), where H3 and H4 posterior probabilities (PP) represent sharing two and one genetic variant, respectively. For the EUR population, we observed four brain volumes with PP higher than 30%: left postcentral [PPH3= 1%, PPH4= 47% ; chr20(q13.33)], left caudal middle frontal [PPH3= 1%, PPH4= 51%; chr20(q13.33)], right precentral [PPH3= 73%, PPH4= 1%; chr2(p16.2)] and gray matter [PPH3= 99%, PPH4= 1% ; chr10(q24.33)]. The top variants associated with both traits in this locus [chr10(q24.33)] are rs11191849 and rs11191848, both mapping toOBFC1 (alternatively known as STN1).

For the EAS population, we found six brain volume measures: gray matter [PPH3= 10%, PPH4= 32% ; chr11 (q22.3)]; right cerebellum white matter [PPH3=1%, PPH4 = 42% ; chr10 (q24 -q25.1)]; right basal forebrain [PPH3= 3%, PPH4 =39% ; chr20 (q13.33)]; cerebellar vermal lobules I-IV [PPH3= 42%, PPH4 =50%; chr7 (q31.33)]; right supramarginal - [PPH3= 14%, PPH4 =65% ; chr10 (q24 -q25.1)], and fourth ventricle [PPH3= 4%, PPH4 = 96% ; chr20 (q13.33)](Figure 2). The strongest PP was seen on chr20 (q13.33) for sharing risk variants between LTL and fourth ventricle which is located dorsal to the brain stem and anterior to the cerebellum, and drains the cerebrospinal fluid (CSF) into the spinal cord (Roesch and Tadi, 2020). The top variant in this locus rs7361098 is in the intronic region of the RTEL1 (regulator of telomere elongation helicase 1; hg19 and 38) and upstream to the RTEL1-TNFRSF6B. On chr7 (q31.33), the cerebellar vermal lobules I-IV had high aggregated probability (92%) of shared genomic risk locus between LTL and volume of cerebellar vermal lobule. (Supplementary file 2). The top variants in this locus were rs34204767 inGPR37 and rs6968500 in C7orf77 (upstream). The vermal lobules are located in the cerebellum midline and affect behavioral and cognitive phenotypes (Yucel et al., 2013).

Figure 2:
Figure 2: Genetic colocalization

A) The ideogram highlights genomic regions wherein genetic variants (SNPs) are shared between leukocyte telomere length (LTL) and respective brain volume measures whose PP >30%. Posterior probability >90% of sharing a causal variants for each LTL-brain volume measure pair is bolded. B) Regional association plot for LTL and fourth ventricle on chr20 (q13.33) and the top significant variants for both traits in the locus is shown followed by genes in the region (hg19 build). C) Regional association plot for LTL and gray matter on chr10 (q24.33) and the top significant variants for both traits in the locus is shown followed by genes in the region (hg19 build). D) Regional association plot for both traits on chr7 (q31.33) and the top two most significant variants in the locus are shown followed by genes in the region (hg19 build).

2.2 Brain tissue based transcriptomic and epigenetic QTLs

We integrated expression meta-analyzed cis-QTL data from different tissues of the brain (brainmeta) with SNPs associated with LTL considering allele frequency and LD structure of the EUR and EAS population separately. For brainmeta eQTL, we found one bonferroni significant gene for the EAS population (DHRS1; p=1.55 x 10−6and pHEIDI = 1.25 x 10−21). None of the genes were bonferroni significant for the EUR population. Significant p-value HEIDI test indicates a linkage model i.e. more than one variant that are probably in linkage, are associated with eQTL and telomere length, whereas non-significant HEIDI p-value indicate pleiotropy under single causal variant model (Figure 3A & 3B).

Figure 3:
Figure 3: Gene-based associations integrating brain tissue eQTL and mQTL profiles

with genetic variants associated with telomere length in EUR and EAS populations. A,B) Bokeh plots showing LTL-associated genes for gene expression (eQTLs) of combined brain tissues performed using SMR. C,D) Bokeh plots showing LTL-associated genes for expression in prefrontal cortex E,F) Bokeh plots showing LTL-associated genes for methylation in combined brain tissues. The left side panels show association in EUR population and EAS associations are shown on the right. The x-axis shows genomic coordinate for the gene/CpG sites and the y-axis shows the -log10 of p-value of SMR test, so the order of increasing significance moves vertically. The dotted line shows the Bonferroni significance and annotated with respective value on top of the line. The size of the dots are scaled to the -log10(p-value) of HEIDI test which tests for pleiotropy under single causal variant model defined as p-valueHEIDI >0.05 or -log10 (p-value) = 1.3 highlighted in pink text label.. Note: Bokeh plots are intersection between manhattan and bubble plots.

Integrating eQTL data from prefrontal cortex (eQTLpfc), we found five genes - LRRC34 (p= 1.15×10−11;pHEIDI = 9.64 x10−4), RP11-362K14.5 (p= 4.38 x10−10;pHEIDI = 0.13), MYNN (p= 1.11 x10−9;pHEIDI = 0.17), RP11-362K14.6 (p= 4.67 x10−7pHEIDI = 0.17), and RP11-541N10.3 (p= 1.18 x10−6; pHEIDI = 0.10) in the EUR population. All the significant genes, except LRRC34 had a non-significant HEIDI p-value, suggesting pleiotropic variants in the risk region sharing association with both eQTLpfc and LTL. We found three genes associated in the EAS population; GPR37 (p=9.75 x10−7; pHEIDI =0.22), ATM (p= 1.38 x10−6; pHEIDI =0.06) and DHRS1 (p= 2.47 x10−6; pHEIDI = 1.11 x10−5). Genes GPR37 and ATM have non-significant HEIDI test statistics, suggesting pleiotropy(single-variant) in the region (Figure 3C & 3D) (Supplementary file 3).

For the brain tissue mQTL data integrated with LTL association statistics, we found 17 significant CpG sites (p ≤4.54 x10-7) in EUR ancestry subjects. Six of the 17 CpG sites were associated under the pleiotropic model - cg00624418 (pHEIDI =0.82), cg00832555 (pHEIDI =0.06), cg01920524 (pHEIDI =0.27), cg21291985 (pHEIDI =0.22), cg08095637 (pHEIDI=0.14)and cg01861555 (pHEIDI =0.40). In the EAS population, we found four CpG sites; cg00624418 (p = 1.77 x10−17; pHEIDI =0.01), cg03935379 (p= 4.60 x10−10; pHEIDI =0.09), cg07380026 (p= 4.37×10−9; pHEIDI =0.01) and cg06294279 (p= 2.65×10−7; pHEIDI =0.35). Of these associations, cg03935379 (TERT) and cg06294279 (GPR37) are associated under the pleiotropic model (Figure 3E & 3F) (Supplementary file 3). Since methylation studies are mostly performed in blood, and methylation levels vary between tissues, we interpreted the CpG sites in blood and different brain tissues using BECon (Edgar et al., 2017) (Supplementary file 4 and Figure 4). Using the EWAS Atlas (Li et al., 2019), we investigated the eight CpG sites that showed pleiotropy, for their disease-based tissue variability of methylation levels (Supplementary file4).

Figure 4:
Figure 4: CpG site attributes and blood-brain tissue correlation.

To aid in interpretation of the identified CpG sites, we used BECon for visualizing the inter-individual and tissue-based variability. (Top) The individual variability across three brain and blood tissue methylation levels is shown in each panel for respective CpG site from BECon data. (Bottom) The panel shows genomic location and annotation of the CpG sites and correlation value of blood methylation levels with each of the three brain tissues methylation. Tabular information is listed in Supplementary file 4.

2.3 Brain chromatin association for tissue and cell-types

We performed association analysis for LTL, by combining Hi-C chromatin SNP-gene associations for developing brain – fetal brain, dorsolateral prefrontal cortex – adult brain, and cell types – astrocytes and neurons using H-MAGMA (Sey et al., 2020). For fetal brain, we observed 21 genes (EUR) and 50 genes (EAS); adult brain (EUR-23; EAS-53); astrocytes (EUR-23; EAS – 53); neurons positive [cells sorted as NeuN+] (EUR – 21; EAS-54) and neurons negative [cells sorted as NeuN-] (EUR-25; EAS-55) (Figure 5).Several genes overlapped across tissue and cell types (unique number of genes - EUR= 50, significant genes lower than p 1.02 x10−06; EAS=97; significant genes lower than p 9.71 x10−07), and their specific distribution is shown in the Venn diagram (Figure 5) (Supplementary file5). Combining genes across both populations and considering those unique to each tissue and cell type, we found unique genes for fetal brain (eight genes), adult brain (nine genes), astrocytes (nine genes), neurons ; NeuN+ (11 genes) and NeuN-(15 genes) (Table 1). These genes were further investigated for enriched gene ontology pathways. We observed processes; calcium ion transport (pFDR =0.012; fetal brain) and G2/M cell cycle transition (pFDR = 0.036; adult brain) at pFDR<0.05 (Figure 5 top-right panel) (Supplementary file5 Table S6 and S7). To further the interpretation, we analyzed these genes using the whole-brain specific protein-protein interaction of all unique genes to tissue or cell types. The network was filtered to retain at least two-degrees, meaning that gene (node) should be connected to at least two other genes (nodes), to avoid selecting orphan/singular branches. This approach prioritizes genes that interact the most in a given network (Figure 5 bottom-right panel). This network revealed 12 intermediate genes that connected our query genes.

View this table:
Table 1: List of genes unique or non-overlapping among the tissue and cell types.

The genes highlighted in bold drive enrichment results discussed in the main text.

Figure 5:
Figure 5: (A)Genes associated with LTL for tissue (fetal and adult brain) and cells (astrocytes and neurons).

The genes are shown as circle data points categorized by tissue and cell types (x-axis) and arranged by genomic position (ascending position; y-axis). The dark green circles show significant genes for the EUR reference panel and EAS is shown as light green. The venn diagram on the right shows the distribution of distinct and overlapping genes for each population. Details are provided in Supplementary file 5. The Venn diagram on top right panel shows the distribution of combined significant genes in tissue and cell types. (B) Gene enrichment analysis of distinct genes. The genes that were unique to each tissue/cell type were analyzed for gene ontology enrichment. The processes are shown as dendograms with their pFDR value on the blue nodes. (C) PPI network for all tissue/cell-type distinct genes. The whole-brain specific PPI (protein-protein interaction) network for neuron genes is shown on the right, with query genes circled. Details of gene enrichment are reported in Supplementary file5 Tables S6 and S7.

2.4 Gene-drug interaction and enrichment of gene functions

We analyzed all significant genes identified from analyses for gene-drug interaction. There were 14 genes on whom the drug interaction was available (Table S1 in Supplementary file 6). To connect the influence of these drugs on gene functions weperformed drug set enrichment analysis. Out of 80 drugs, 10 drugs had associated information with gene function sets. There were 24 nominally significant (p<0.05) gene functions processes and one process i.e. transmission across chemical synapses was FDR significant (p= 0.00028; pFDR=0.019) (Fig S1 and Table S2 in Supplementary file 6).

3 Discussion

In this study, we investigated the role of genetic variants associated with leukocyte telomere length (LTL) in the context of brain-based expression, methylation, and chromatin profiles of two tissue types and three cell types. Analyzing different brain morphological traits, we discovered that the volumes of global gray matter, fourth ventricle and cerebellar vermal lobules I-IV, have genomic causal loci in common with LTL. The shared locus with global gray matter reported in our results corroborates previous findings of LTL association with cortical gray matter (King et al., 2014). A study of 15,892 individuals including subjects diagnosed with six psychiatric disorders showed gray matter loss across segmental brain regions in comparison to controls (Goodkind et al., 2015). The fourth ventricle bordered by the brainstem and the cerebellum, is one of the sites for CSF production and is the last ventricle connecting to the subarachnoid space (Roesch and Tadi, 2020). Pathologies involving the fourth ventricle are associated with abnormal CSF volume and flow (Whedon and Glassey, 2009). Different neuropsychiatric disorders such as schizophrenia have demonstrated enlarged fourth ventricle – which reflects a deficit in bordering tissue when compared to brain samples from healthy subjects (Juuhl-Langseth et al., 2012). Schizophrenia also has association with reduced telomere length (Russo et al., 2018). Furthermore, an extreme pathology (i.e. congenital deformity of the fourth ventricle) causes Dandy-Walker malformation which exhibits some symptoms with similarities to those seen in schizophrenia, obsessive compulsive disorder, and type-I and type-II bipolar disorder (Lingeswaran et al., 2009). Furthermore, the cerebellar vermal lobules I-IV are the top four anterior lobes connecting both hemispheres of the cerebellum, and are involved in motor, cognitive, and behavioral processing in autism spectrum disorders(D’Mello et al., 2015).Taken together, our results further corroborate the relationship between CNS morphology and telomere length.

The brain transcriptomic integrated association studies highlighted genes that are involved with both LTL and neuronal effects. GPR37 (G-protein coupled receptor 37) is expressed in oligodendrocytes in different tissues of the brain contributing to myelination regulation (Smith et al., 2017). GPR37 is associated with telomere length and serves as stress response neuroprotective gene. In major depression disorder, GPR37 expression levels in brain tissue were lower compared to healthy participants (Mamdani et al., 2015). Furthermore, CpG site - cg06294279 located in the GPR37 intragenic region was also associated under pleiotropic genomic risk (i.e. sharing association with both methylation level of the site and telomere length). The top variants in the colocalization analysis between LTL and cerebellar vermal lobules overlap with GPR37. The BECon databank shows that Broadman areas 7, 10, and 20 share little to low negative correlation with blood-based methylation in healthy individuals. To understand the role of the CpG site in diseases, we mined the EWAS Datahub (Li et al., 2019); CpG site - cg06294279 showed high variability in different brain tissues between cases and controls, for schizophrenia and Alzheimer’s disease (Supplementaryfile4).LRRC34 (leucine rich repeat containing 34) gene expression is high in stem cells which lowers upon cell differentiation(Lührig et al., 2014). Individuals with Bipolar disorder who were on lithium showed changes in SNP-based transcriptome effect of LRRC34 (Coutts et al., 2019). Furthermore, we found six CpG sites in the LRRC34 gene that displayed mQTL associations for SNPs associated with LTL, out of which, site-cg08095637 had single-variant pleiotropic association.

Among genes identified from epigenomic integrated associations, PARP1 (Poly(ADP-Ribose) Polymerase 1) is a regulator of DNA repair and cell death taking part in neurogenesis of neuronal and embryonic stem cells (Hong et al., 2019). Lowering the expression levels of PARP1 is linked with abnormal ventricle volume in the brain and exhibits psychiatric symptoms including anxiety, depression, and behavioral impairment (Hong et al., 2019; Sriram et al., 2016). ATM (Ataxia Telangiectasia Mutated) is also involved in cellular DNA repair and differentiating pluripotent stem cells (Corti et al., 2019). ATM regulates GABAergic based neuronal plasticity and brain development and its dysregulation is associated with onset of different psychopathologies (Pizzamiglio et al., 2016). The mQTL analysis of brain tissue found CpG sites in addition to LRRC34, in the TERT gene. Methylation based biological age acceleration is associated with TERT and genetic overlap with schizophrenia, bipolar disorder, and dementia (Lu et al., 2018). TERT is also involved in paternally inherited epigenetic imprint observed in advanced paternal age which is a risk for several neuropsychiatric disorder (Hehar et al., 2017).Methylation levels of ZBTB46 are associated with cognitive function (Starnawska et al., 2017) and schizophrenia in a RNA-sequencing gene expression study(Maycox et al., 2009; Zhang et al., 2020).

By integrating chromatin profiles for brain at two stages of development and different cell types, we found several overlapping genes across brain tissue and cell specific genes. For gene specific to Hi-C associations from fetal brain, we observed calcium ion transport process; this enrichment was driven by JPH4, and CUL5. JPH4 belongs to the junctophilin family, is primarily expressed in neurons and regulates calcium signaling in normal brain development (Landstrom et al., 2014). CUL5 (cullin5) encodes the vasopressin-activated calcium-mobilizing receptor (VACM-1) and is responsible for chromosomal stability and DNA repair mechanisms (Tapia-Laliena et al., 2019). Neural stem cells are precursors of neurons and glial cells in the brain and periodic responses from calcium signaling drive cell-specific neurogenesis (Toth et al., 2016). In the adultbrain, we observed enrichment of the G2/M cell cycle transition process. While neurons are generally not believed to enter M-phase, cell cycle dysregulation arising from pathological conditions can lead to undifferentiated neurons back in to the cell cycle and increase vulnerability to cell death (Walton et al., 2019). Reactivation of the cell cycle in adult neurons occurs in the circumstance of degeneration or injury to the central nervous system(Frade and Ovejero-Benito, 2015). Two genes driving cell cycle ontology were EPHX1 and PSME1. Both EPHX1 and PSME1 are high expressed in the brain, EPHX1 regulates endocannabinoid 2-arachidonoylglycerol in depression, stress, Alzheimer’s disease, and methamphetamine-dependence (Nithipatikom et al.,2014. Václavíková et al., 2015). Lower transcript levels of PSME1 (proteasome activator subunit 1) have been observed in individuals with major depression disorder compared to healthy/control population (Redei et al., 2014).

We also observed several brain-based mQTLs which showed pleiotropic (single variant) association with LTL SNPs and several genes that were significant when integrating chromatin profiles. These results are in line with the genomic link between TERT and accelerated biological aging derived from methylation expression (Lu et al., 2018). The drug-gene interaction of significant genes showed enrichment for different transporter systems of chemical synapses, aquaporins and mitochondrial tRNA aminoacylation and transmission in the postsynaptic cell (Fig S1: Supplementary file 6). Together, these findings indicate the set of LTL-associated genes identified herein are associated with brain-related traits due to their regulatory effects on brain tissues and cells. Therefore, we believe telomere length genomics is an important avenue for future studies of mental health.

Most significant genes from one population are nominally significant in the other. However, there were few regions that are significant in one population, but not the other. These differences are due to recombination of the SNPs and allele frequency differences in the locus. Five loci on chromosomes 3, 4, 5, 10 and 20 overlap in both populations, however, the SNPs in these regions show differences in linkage disequilibrium highlighting the importance of studying the population based genetic association. These differences in SNP-peaks are visualized in genome-wide and regional plots (Supplementary file 1).

While our study has provided key evidence of shared genomic causal loci among telomere length, brain morphology, and regulatory traits, it has limitations. The results showing genetic co-localizing cannot be interpreted as mediating the direction of causality. Furthermore, colocalization of two traits under a single variant is a conservative model of testing, and there is biological plausibility of more than one variant exhibiting small effects. There are several SNPs in a locus exhibiting different effect directions for LTL and brain volume (Fig S2: Supplementary file 6). Therefore, it is difficult to assess the phenotypic consequences of each SNP because variants may alter physiological states in tissue or cell-type specific manner. We assessed regulatory changes using two populations and therefore the results may not apply to other major populations and will require separate investigation. Further studies using animal models would narrow targets with specific cellular changes. Our next steps will involve investigating the role of telomere length with respect to neuropsychiatric brain-disorders to deepen the understanding of cellular aging in neurological outcomes.

4 Conclusion

We identified shared genomic causal loci among those affecting telomere length, brain morphology and regulatory traits. The identified genes show involvement in the central nervous system, and provide evidence that shared genetics could explain, in part, the association of LTL with several brain-based outcomes reported previously. These results highlight that telomere length may serve as a marker for certain brain-based diseases such as psychiatric disorders. These findings additionally corroborate that in addition to methylation based biological aging, telomere-based cellular aging has the potential to provide further resolution on biological mechanisms of brain disorders.

5 Methods

All data used in this study were made publicly available by the cited consortia. This study was exempted from requiring institutional review board approval due to its use of de-identified aggregated data.

5.1 Cross-ancestry comparison of LTL-associations

We used the largest available LTL GWAS for EUR and EAS population. There are seven genome loci in EUR population, and nine independent loci in EAS population associated with leukocyte telomere length (LTL) (Dorajoo et al., 2019). Genomic risk regions for LTL were identified using FUMA (default values; LD -0.6 and 250kb locus) (Watanabe et al., 2017). A locus contains several SNPs independent significant SNPs by merging LD blocks within 250kb. Loci on chromosomes 2 and 19 that are observed in EURs are not present in EAS population. Similarly, there are loci specific to EAS population are on chromosomes 1, 7 and 14. Regions which are positionally similar in both populations are on chromosomes 3, 4, 5, 10 and 20. Length of these genomic regions varies by population due to LD structure and different number of significant variants within these regions. Therefore, in this study we investigated pleiotropic associations of SNPs with LTL and quantitative expression of brain tissues in both populations (Supplementary file1).

5.2 Brain morphology colocalization

To evaluate which loci are shared between LTL and brain volume measures, we used the largest available GWAS of 101 phenotypes of brain morphology(B. Zhao et al., 2019). In this study, the brain volume measures or region-of-interest (ROI) were quantified using magnetic resonance imaging and mapped to catalogued ROI to identify (i.e. map them to) structural measures (B. Zhao et al., 2019). The GWAS of 101 brain morphology/region-of-interest measures includes meta-analysis of five cohorts – the Human Connectome Project (HCP) study, the Pediatric Imaging, Neurocognition, and Genetics (PING) study, the Philadelphia Neurodevelopmental Cohort (PNC) study, and the Alzheimer’s Disease Neuroimaging Initiative (ADNI) study (n=21,821). The GWAS association statistics for each of the 101 brain morphology measures (B. Zhao et al., 2019) were evaluated with respect to GWAS of LTL (Dorajoo et al., 2019) to identify which causal risk loci were shared across both traits using coloc (Giambartolomei et al., 2014). We performed pairwise colocalization analysis for each of the risk locus between LTL and each brain volume phenotypes using the appropriate 1000 Genomes Phase 3 reference panels (i.e., EUR and EAS). The coloc method in R-v3.6, evaluates posterior probability (PP) for four alternative hypothesis – A) H1: association with trait 1, not with trait 2, B) H2: association with trait 2, not with trait 1, C) H3: association with trait 1 and trait 2, two independent SNPs, and D) H4: association with trait 1 and trait 2, one shared SNP (Giambartolomei et al., 2014). We used the coloc.abf function using p-values, respective ancestry’s allele frequencies and sample size with default priors. Here we report results for sharing single causal variant with both traits (LTL and brain morphology) having H4 PP higher than 30% in the main text, detailed results can be found in supplementary file 2.

5.3 Transcriptomic and epigenetic QTLs in brain tissues

The association of LTL was integrated with eQTL and mQTL information derived from Genotype-Tissue Expression (GTEx) (GTEx Consortium et al., 2017), CommonMind Consortium (CMC) (Fromer et al., 2016), Religious Orders Study and the Memory and Aging Project (ROSMAP) and PsychENCODE consortia (Gandal et al., 2018). Data were included from two transcriptomic studies: a) multiple brain tissue i.e. meta-analysis of gene expression from different regions of the brain (brainmeta) and b) prefrontal cortex (brainpfc), and methylation association study using meta-analyzed methylation expression from multiple regions of the brain using Summary-based Mendelian Randomization (SMR) (Zhu et al., 2016). Bonferroni correction (0.05/number of associations) was applied to identify significant gene associations for each of the six QTL analysis. We performed heterogeneity in dependent instruments (HEIDI) association analysis to test if a single variant is associated with phenotype and QTL (single variant model of pleiotropy). Significant HEIDI associations (p< 0.05) depart from this assumption of single variant association while non-significant associations (p > 0.05) do not reject the null, indicating a linkage model (i.e. different SNPs in linkage) (Zhu et al., 2016). The analyses were performed in each ancestry group (i.e., EUR and EAS) separately. Two sets of files were created using LD and SNP allele frequencies of each population from the 1000 Genome Phase 3 reference panel. This resulted in different number of SNPs considered for each population. Misaligned allele direction between QTL and LTL-SNP dataset were removed. The CpG sites which showed pleiotropic evidence i.e. single causal variant between LTL and brain mQTL were interpreted using BECon for methylation level correlation between brain tissues and blood (Edgar et al., 2017) and extracted from EWAS DataHub for different brain disorders (Li et al., 2019).

5.4 Brain chromatin profiles using H-MAGMA

The gene-based association study for brain-derived chromatin profiles for fetal and adult brain regions and cells (i.e. astrocytes and neurons) was performed using H-MAGMA (Sey et al., 2020). This approach aggregates genetic variants to nearest genes derived from Hi-C (high-resolution chromatin conformation) data of fetal developing cortex, adult dorsolateral prefrontal cortex, and cell-types - iPSC derived astrocytes and cells sorted as NeuN+ and NeuN-from anterior cingulate cortex(Rajarajan et al., 2018). The H-MAGMA approach aids in providing functional/regulatory effects of non-coding SNPs, and is a complementary to transcriptome-wide association studies and the coloc approach used in this study (Sey et al., 2020). Due to telomeres/TERT’s role in the central nervous system and strong PP seen with brain volume measures, we tested all five tissue profiles with LTL GWAS association statistics. Each tissue/cell type association was performed for EUR and EAS population reference panels. Bonferroni correction (0.05/number of associations) was applied to identify gene associations for each the five tissue/cell type analyses. The genes that were unique to each tissue or cell type were analyzed for pathway enrichment using ShinyGO (Ge et al., 2019). Additionally, all unique genes were investigated for their role in whole-brain specific protein-protein interaction network using NetworkAnalyst3.0 (Zhou et al., 2019) retaining genes with at least two genes. The visualizations were created in ggplot, Gviz and PhenoGram (Wolfe et al., 2013). All significant genes from aforementioned analyses were investigated for gene-drug interactions using DGidb(Griffith et al., 2013). Then analyzed for gene function enrichment using DSEA(Napolitano et al., 2016).

Data Availability

All the data used in this study are publicly available, hosted at the respective consortia portal. Chinese and European population telomere length GWAS (https://www.nature.com/articles/s41467-019-10443-2#data-availability); Brain Region Phenotypes GWAS (https://github.com/BIG-S2/GWAS);eQTL and mQTL data is available from SMR (https://cnsgenomics.com/software/smr/#Overview); H-MAGMA (https://github.com/thewonlab/H-MAGMA). The results generated in this study are provided in the main text and supplementary files

7 Competing Interests

Drs. Polimanti and Gelernter are paid for their editorial work for the journal ‘Complex Psychiatry’. The other authors have no conflicts of interest to declare.

8 Data Availability

All the data used in this study are publicly available, hosted at the respective consortia portal. Chinese and European population telomere length GWAS (https://www.nature.com/articles/s41467-019-10443-2#data-availability);Brain Region Phenotypes GWAS (https://github.com/BIG-S2/GWAS);eQTL and mQTL data is available from SMR (https://cnsgenomics.com/software/smr/#Overview); H-MAGMA (https://github.com/thewonlab/H-MAGMA). The results generated in this study are provided in the main text and supplementary files.

11 Supplementary files

  • Supplementaryfile1.pdf: Genome-wide and regional plots of LTL-GWAS in EUR and EAS population showing annotation of CADD scores, RegulomeDB, regulatory sites, and HiC chromatin interactions

  • Supplementaryfile2.xlsx: Details of genetic colocalization between leukocyte telomere length (LTL) of EUR and EAS populations with 101 brain morphology measures

    • Table S1: Genomic locus in each population and corresponding five posterior probability values for each region reported in the manuscript

    • Table S2: Legend of brain volume measures analyzed with LTL

    • Table S3: Details of posterior probability values for all brain regions for Region1 (chr2:54461744-54498311) of EUR population

    • Table S4: Details of posterior probability values for all brain regions for Region2 (chr3:169392323-169599567) of EUR population

    • Table S5: Details of posterior probability values for all brain regions for Region3 (chr4:164001945-164126538) of EUR population

    • Table S6: Details of posterior probability values for all brain regions for Region4 (chr5:1279790-1355859) of EUR population

    • Table S7: Details of posterior probability values for all brain regions for Region5 (chr10:105608838-105694301) of EUR population

    • Table S8: Details of posterior probability values for all brain regions for Region6 (chr19:22089575-22319421) of EUR population

    • Table S9: Details of posterior probability values for all brain regions for Region7 (chr20:62421622-62467796) of EUR population

    • Table S10: Details of posterior probability values for all brain regions for Region1 (chr1:226376883-226653478) of EAS population

    • Table S11: Details of posterior probability values for all brain regions for Region2 (chr3:169360327-169582223) of EAS population

    • Table S12: Details of posterior probability values for all brain regions for Region3 (chr4:164007820-164126538) of EAS population

    • Table S13: Details of posterior probability values for all brain regions for Region4 (chr5:1276050-1343794) of EAS population

    • Table S14: Details of posterior probability values for all brain regions for Region5 (chr7:124366188-124907090) of EAS population

    • Table S15: Details of posterior probability values for all brain regions for Region6 (chr10:105671683-105801818) of EAS population

    • Table S16: Details of posterior probability values for all brain regions for Region7 (chr11:107983015-108367453) of EAS population

    • Table S17: Details of posterior probability values for all brain regions for Region8 (chr14:24623428-24730100) of EAS population

    • Table S18: Details of posterior probability values for all brain regions for Region9 (chr20:62275844-62335293) of EAS population

  • Supplementaryfile3.pdf: Details of significant LTL associations for eQTL (gene expression) of meta-analyzed brain tissues and prefrontal cortex, and mQTL associations of meta-analyzed brain tissues in each population. Comparative results of all significant associations from one population are shown in other population as well.

  • Supplementaryfile4.pdf: Visualization of mQTL sites that showed evidence of pleiotropy (single variant model from non-significant p-value HEIDI) were mined from EWAS DataHub in EWAS Atlas to show tissue-based variability in neuropsychiatric phenotypes

  • Supplementaryfile5.xlsx: Details of chromatin-profiles associations for tissue and cell types in EUR and EAS population

    • Table S1: List of significant genes in fetal brain tissue for LTL associations in EUR and EAS population. Comparative results of all significant associations from one population are shown in other population as well.

    • Table S2: List of significant genes in adult brain tissue for LTL associations in EUR and EAS population. Comparative results of all significant associations from one population are shown in other population as well.

    • Table S3: List of significant genes in astrocyte cell types for LTL associations in EUR and EAS population. Comparative results of all significant associations from one population are shown in other population as well.

    • Table S4: List of significant genes in positive neurons cell types for LTL associations in EUR and EAS population. Comparative results of all significant associations from one population are shown in other population as well.

    • Table S5: List of significant genes in negative neuron cell type for LTL associations in EUR and EAS population. Comparative results of all significant associations from one population are shown in other population as well.

    • Table S6: List of genes that were distinct/non-overlapping in each tissue and cell-types and gene ontology enrichment results are shown next to the respective list of genes.

    • Table S7: Gene ontology enrichment results of the network from unique genes shown in Figure 6 (right panel).

    • Table S8: List of Bonferroni thresholds for each tissue and cell type in both populations and the number of significant associations.

  • Supplementary file 6. Details of gene-drug interaction, drug set enrichment and effect direction of SNPs in the loci that are colocalized

    • Table S1: Gene-Drug interaction mined from The Drug-Gene interaction database

    • Fig S1. Drug Enrichment Set Analysis (DSEA). Analyzing enriched gene functions for drugs identified from drug-gene interaction (as shown above)

    • Table S2: Tabular associations of DSEA shown in the above figure

    • Fig S2: Effect direction of SNPs in the shared causal region for LTL and brain volume measures

6 Acknowledgements

This study was supported by the National Institutes of Health (R21 DC018098, R21 DA047527, and F32 MH122058). We would like to acknowledge the participants and the investigators of all consortia’s whose data helped investigate the research question.

9 References

  1. ↵
    Aas, M., Elvsåshagen, T., Westlye, L.T., Kaufmann, T., Athanasiu, L., Djurovic, S., Melle, I., van der Meer, D., Martin-Ruiz, C., Steen, N.E., Agartz, I., Andreassen, O.A., 2019. Telomere length is associated with childhood trauma in patients with severe mental disorders. Transl. Psychiatry 9, 97. doi:10.1038/s41398-019-0432-7
    OpenUrlCrossRef
  2. ↵
    Aas, M., Henry, C., Andreassen, O.A., Bellivier, F., Melle, I., Etain, B., 2016. The role of childhood trauma in bipolar disorders. Int. J. Bipolar Disord. 4, 2. doi:10.1186/s40345-015-0042-0
    OpenUrlCrossRef
  3. ↵
    Bär, C., Blasco, M.A., 2016. Telomeres and telomerase as therapeutic targets to prevent and treat age-related diseases. [version 1; peer review: 4 approved]. F1000Res. 5. doi:10.12688/f1000research.7020.1
    OpenUrlCrossRef
  4. ↵
    Cohen, J., Torres, C., 2019. Astrocyte senescence: Evidence and significance. Aging Cell 18, e12937. doi:10.1111/acel.12937
    OpenUrlCrossRef
  5. ↵
    Corti, A., Sota, R., Dugo, M., Calogero, R.A., Terragni, B., Mantegazza, M., Franceschetti, S., Restelli, M., Gasparini, P., Lecis, D., Chrzanowska, K.H., Delia, D., 2019. DNA damage and transcriptional regulation in iPSC-derived neurons from Ataxia Telangiectasia patients. Sci. Rep. 9, 651. doi:10.1038/s41598-018-36912-0
    OpenUrlCrossRef
  6. ↵
    Coutts, F., Palmos, A.B., Duarte, R.R.R., de Jong, S., Lewis, C.M., Dima, D., Powell, T.R., 2019. The polygenic nature of telomere length and the anti-ageing properties of lithium. Neuropsychopharmacology 44, 757–765. doi:10.1038/s41386-018-0289-0
    OpenUrlCrossRef
  7. ↵
    D’Mello, A.M., Crocetti, D., Mostofsky, S.H., Stoodley, C.J., 2015. Cerebellar gray matter and lobular volumes correlate with core autism symptoms. Neuroimage Clin. 7, 631–639. doi:10.1016/j.nicl.2015.02.007
    OpenUrlCrossRefPubMed
  8. ↵
    Dorajoo, R., Chang, X., Gurung, R.L., Li, Z., Wang, L., Wang, R., Beckman, K.B., Adams-Haduch, J., M, Y., Liu, S., Meah, W.Y., Sim, K.S., Lim, S.C., Friedlander, Y., Liu, J., van Dam, R.M., Yuan, J.-M., Koh, W.-P., Khor, C.C., Heng, C.-K., 2019. Loci for human leukocyte telomere length in the Singaporean Chinese population and trans-ethnic genetic studies. Nat. Commun. 10, 2491. doi:10.1038/s41467-019-10443-2
    OpenUrlCrossRef
  9. ↵
    Edgar, R.D., Jones, M.J., Meaney, M.J., Turecki, G., Kobor, M.S., 2017. BECon: a tool for interpreting DNA methylation findings from blood in the context of brain. Transl. Psychiatry 7, e1187. doi:10.1038/tp.2017.171
    OpenUrlCrossRef
  10. ↵
    Edlow, A.G., Guedj, F., Sverdlov, D., Pennings, J.L.A., Bianchi, D.W., 2019. Significant effects of maternal diet during pregnancy on the murine fetal brain transcriptome and offspring behavior. Front. Neurosci. 13, 1335. doi:10.3389/fnins.2019.01335
    OpenUrlCrossRef
  11. ↵
    Ferrón, S.R., Marqués-Torrejón, M.A., Mira, H., Flores, I., Taylor, K., Blasco, M.A., Fariñas, I., 2009. Telomere shortening in neural stem cells disrupts neuronal differentiation and neuritogenesis. J. Neurosci. 29, 14394–14407. doi:10.1523/JNEUROSCI.3836-09.2009
    OpenUrlAbstract/FREE Full Text
  12. ↵
    Forero, D.A., González-Giraldo, Y., López-Quintero, C., Castro-Vega, L.J., Barreto, G.E., Perry, G., 2016. Meta-analysis of Telomere Length in Alzheimer’s Disease. J. Gerontol. A, Biol. Sci. Med. Sci. 71, 1069–1073. doi:10.1093/gerona/glw053
    OpenUrlCrossRefPubMed
  13. ↵
    Frade, J.M., Ovejero-Benito, M.C., 2015. Neuronal cell cycle: the neuron itself and its circumstances. Cell Cycle 14, 712–720. doi:10.1080/15384101.2015.1004937
    OpenUrlCrossRefPubMed
  14. ↵
    Fromer, M., Roussos, P., Sieberts, S.K., Johnson, J.S., Kavanagh, D.H., et al., 2016. Gene expression elucidates functional impact of polygenic risk for schizophrenia. Nat. Neurosci. 19, 1442–1453. doi:10.1038/nn.4399
    OpenUrlCrossRefPubMed
  15. ↵
    Gandal, M.J., Zhang, P., Hadjimichael, E., Walker, R.L., Chen, C., Liu, S., Won, H., van Bakel, H., Varghese, M., Wang, Y., Shieh, A.W., Haney, J., Parhami, S., Belmont, J., Kim, M.,Moran Losada, P., Khan, Z., Mleczko, J., Xia, Y., Dai, R., Wang, D., Yang, Y.T., Xu, M., Fish, K., Hof, P.R., Warrell, J., Fitzgerald, D., White, K., Jaffe, A.E., PsychENCODE Consortium, Peters, M.A., Gerstein, M., Liu, C., Iakoucheva, L.M., Pinto, D., Geschwind, D.H., 2018. Transcriptome-wide isoform-level dysregulation in ASD, schizophrenia, and bipolar disorder. Science 362. doi:10.1126/science.aat8127
    OpenUrlAbstract/FREE Full Text
  16. ↵
    Ge, S.X., Jung, D., Yao, R., 2019. ShinyGO: a graphical enrichment tool for animals and plants. Bioinformatics. doi:10.1093/bioinformatics/btz931
    OpenUrlCrossRef
  17. ↵
    Giambartolomei, C., Vukcevic, D., Schadt, E.E., Franke, L., Hingorani, A.D., Wallace, C., Plagnol, V., 2014. Bayesian test for colocalisation between pairs of genetic association studies using summary statistics. PLoS Genet. 10, e1004383. doi:10.1371/journal.pgen.1004383
    OpenUrlCrossRefPubMed
  18. ↵
    Giral, H., Landmesser, U., Kratzer, A., 2018. Into the Wild: GWAS Exploration of Non-coding RNAs. Front. Cardiovasc. Med. 5, 181. doi:10.3389/fcvm.2018.00181
    OpenUrlCrossRef
  19. ↵
    Goodkind, M., Eickhoff, S.B., Oathes, D.J., Jiang, Y., Chang, A., Jones-Hagata, L.B., Ortega, B.N., Zaiko, Y.V., Roach, E.L., Korgaonkar, M.S., Grieve, S.M., Galatzer-Levy, I., Fox, P.T., Etkin, A., 2015. Identification of a common neurobiological substrate for mental illness. JAMA Psychiatry 72, 305–315. doi:10.1001/jamapsychiatry.2014.2206
    OpenUrlCrossRefPubMed
  20. ↵
    Griffith, M., Griffith, O.L., Coffman, A.C., Weible, J.V., McMichael, J.F., Spies, N.C., Koval, J., Das, I., Callaway, M.B., Eldred, J.M., Miller, C.A., Subramanian, J., Govindan, R., Kumar, R.D., Bose, R., Ding, L., Walker, J.R., Larson, D.E., Dooling, D.J., Smith, S.M., Ley, T.J., Mardis, E.R., Wilson, R.K., 2013. DGIdb: mining the druggable genome.Nat. Methods 10, 1209–1210. doi:10.1038/nmeth.2689
    OpenUrlCrossRefPubMedWeb of Science
  21. GTEx Consortium, Laboratory, Data Analysis &Coordinating Center (LDACC)—Analysis Working Group, Statistical Methods groups—Analysis Working Group, Enhancing GTEx (eGTEx) groups, NIH Common Fund, NIH/NCI, NIH/NHGRI, NIH/NIMH, NIH/NIDA, Biospecimen Collection Source Site—NDRI, Biospecimen Collection Source Site—RPCI, Biospecimen Core Resource—VARI, Brain Bank Repository— University of Miami Brain Endowment Bank, Leidos Biomedical—Project Management, ELSI Study, Genome Browser Data Integration &Visualization—EBI, Genome Browser Data Integration &Visualization—UCSC Genomics Institute, University of California Santa Cruz, Lead analysts:, Laboratory, Data Analysis &Coordinating Center (LDACC):, NIH program management:, Biospecimen collection:, Pathology:, eQTL manuscript working group:,Battle, A., Brown, C.D., Engelhardt, B.E., Montgomery, S.B., 2017. Genetic effects on gene expression across human tissues. Nature 550, 204–213. doi:10.1038/nature24277
    OpenUrlCrossRefPubMedWeb of Science
  22. ↵
    Hehar, H., Ma, I., Mychasiuk, R., 2017. Intergenerational Transmission of Paternal Epigenetic Marks: Mechanisms Influencing Susceptibility to Post-Concussion Symptomology in a Rodent Model. Sci. Rep. 7, 7171. doi:10.1038/s41598-017-07784-7
    OpenUrlCrossRef
  23. ↵
    Henje Blom, E., Han, L.K.M., Connolly, C.G., Ho, T.C., Lin, J., LeWinn, K.Z., Simmons, A.N.,Sacchet, M.D., Mobayed, N., Luna, M.E., Paulus, M., Epel, E.S., Blackburn, E.H., Wolkowitz, O.M., Yang, T.T., 2015. Peripheral telomere length and hippocampal volume in adolescents with major depressive disorder. Transl. Psychiatry 5, e676. doi:10.1038/tp.2015.172
    OpenUrlCrossRef
  24. ↵
    Hong, S., Yi, J.H., Lee, S., Park, C.-H., Ryu, J.H., Shin, K.S., Kang, S.J., 2019. Defective neurogenesis and schizophrenia-like behavior in PARP-1-deficient mice. Cell Death Dis. 10, 943. doi:10.1038/s41419-019-2174-0
    OpenUrlCrossRef
  25. ↵
    Juuhl-Langseth, M., Rimol, L.M., Rasmussen, I.A., Thormodsen, R., Holmén, A., Emblem, K.E., Due-Tønnessen, P., Rund, B.R., Agartz, I., 2012. Comprehensive segmentation of subcortical brain volumes in early onset schizophrenia reveals limited structural abnormalities. Psychiatry Res. 203, 14–23. doi:10.1016/j.pscychresns.2011.10.005
    OpenUrlCrossRef
  26. ↵
    Kaczmarczyk, M., Wingenfeld, K., Kuehl, L.K., Otte, C., Hinkelmann, K., 2018. Childhood trauma and diagnosis of major depression: Association with memory and executive function. Psychiatry Res. 270, 880–886. doi:10.1016/j.psychres.2018.10.071
    OpenUrlCrossRef
  27. ↵
    Kananen, L., Surakka, I., Pirkola, S., Suvisaari, J., Lönnqvist, J., Peltonen, L., Ripatti, S., Hovatta, I., 2010. Childhood adversities are associated with shorter telomere length at adult age both in individuals with an anxiety disorder and controls. PLoS One 5, e10826. doi:10.1371/journal.pone.0010826
    OpenUrlCrossRefPubMed
  28. ↵
    King, K.S., Kozlitina, J., Rosenberg, R.N., Peshock, R.M., McColl, R.W., Garcia, C.K., 2014. Effect of leukocyte telomere length on total and regional brain volumes in a large population-based cohort. JAMA Neurol. 71, 1247–1254. doi:10.1001/jamaneurol.2014.1926
    OpenUrlCrossRefPubMed
  29. ↵
    Koliada, A.K., Krasnenkov, D.S., Vaiserman, A.M., 2015. Telomeric aging: mitotic clock or stress indicator? Front. Genet. 6, 82. doi:10.3389/fgene.2015.00082
    OpenUrlCrossRef
  30. ↵
    Landstrom, A.P., Beavers, D.L., Wehrens, X.H.T., 2014. The junctophilin family of proteins: from bench to bedside. Trends Mol. Med. 20, 353–362. doi:10.1016/j.molmed.2014.02.004
    OpenUrlCrossRefPubMed
  31. ↵
    Li, M., Zou, D., Li, Z., Gao, R., Sang, J., Zhang, Y., Li, R., Xia, L., Zhang, T., Niu, G., Bao, Y., Zhang, Z., 2019. EWAS Atlas: a curated knowledgebase of epigenome-wide association studies. Nucleic Acids Res. 47, D983–D988. doi:10.1093/nar/gky1027
    OpenUrlCrossRef
  32. ↵
    Lingeswaran, A., Barathi, D., Sharma, G., 2009. Dandy-Walker variant associated with bipolar affective disorder. J. Pediatr. Neurosci. 4, 131–132. doi:10.4103/1817-1745.57341
    OpenUrlCrossRefPubMed
  33. ↵
    Liu, M.-Y., Nemes, A., Zhou, Q.-G., 2018. The emerging roles for telomerase in the central nervous system. Front. Mol. Neurosci. 11, 160. doi:10.3389/fnmol.2018.00160
    OpenUrlCrossRefPubMed
  34. ↵
    Lobanova, A., She, R., Pieraut, S., Clapp, C., Maximov, A., Denchi, E.L., 2017. Different requirements of functional telomeres in neural stem cells and terminally differentiated neurons. Genes Dev. 31, 639–647. doi:10.1101/gad.295402.116
    OpenUrlAbstract/FREE Full Text
  35. ↵
    Lu, A.T., Xue, L., Salfati, E.L., Chen, B.H., Ferrucci, L., Levy, D., Joehanes, R., Murabito, J.M.,Kiel, D.P., Tsai, P.-C., Yet, I., Bell, J.T., Mangino, M., Tanaka, T., McRae, A.F., Marioni, R.E., Visscher, P.M., Wray, N.R., Deary, I.J., Levine, M.E., Quach, A., Assimes, T., Tsao, P.S., Absher, D., Stewart, J.D., Li, Y., Reiner, A.P., Hou, L., Baccarelli, A.A., Whitsel, E.A., Aviv, A., Cardona, A., Day, F.R., Wareham, N.J., Perry, J.R.B., Ong, K.K., Raj, K., Lunetta, K.L., Horvath, S., 2018. GWAS of epigenetic aging rates in blood reveals a critical role for TERT. Nat. Commun. 9, 387. doi:10.1038/s41467-017-02697-5
    OpenUrlCrossRef
  36. ↵
    Lührig, S., Siamishi, I., Tesmer-Wolf, M., Zechner, U., Engel, W., Nolte, J., 2014. Lrrc34, a novel nucleolar protein, interacts with npm1 and ncl and has an impact on pluripotent stem cells. Stem Cells Dev. 23, 2862–2874. doi:10.1089/scd.2013.0470
    OpenUrlCrossRef
  37. ↵
    Mamdani, F., Rollins, B., Morgan, L., Myers, R.M., Barchas, J.D., Schatzberg, A.F., Watson, S.J., Akil, H., Potkin, S.G., Bunney, W.E., Vawter, M.P., Sequeira, P.A., 2015. Variable telomere length across post-mortem human brain regions and specific reduction in the hippocampus of major depressive disorder. Transl. Psychiatry 5, e636. doi:10.1038/tp.2015.134
    OpenUrlCrossRef
  38. ↵
    Maycox, P.R., Kelly, F., Taylor, A., Bates, S., Reid, J., Logendra, R., Barnes, M.R., Larminie, C., Jones, N., Lennon, M., Davies, C., Hagan, J.J., Scorer, C.A., Angelinetta, C., Akbar, M.T., Hirsch, S., Mortimer, A.M., Barnes, T.R.E., de Belleroche, J., 2009. Analysis of gene expression in two large schizophrenia cohorts identifies multiple changes associated with nerve terminal function. Mol. Psychiatry 14, 1083–1094. doi:10.1038/mp.2009.18
    OpenUrlCrossRefPubMedWeb of Science
  39. ↵
    McLaughlin, K.A., Greif Green, J., Gruber, M.J., Sampson, N.A., Zaslavsky, A.M., Kessler, R.C., 2012. Childhood adversities and first onset of psychiatric disorders in a national sample of US adolescents. Arch. Gen. Psychiatry 69, 1151–1160. doi:10.1001/archgenpsychiatry.2011.2277
    OpenUrlCrossRefPubMedWeb of Science
  40. ↵
    McLaughlin, K.A., Koenen, K.C., Bromet, E.J., Karam, E.G., Liu, H., Petukhova, M., Ruscio, A.M., Sampson, N.A., Stein, D.J., Aguilar-Gaxiola, S., Alonso, J., Borges, G., Demyttenaere, K., Dinolova, R.V., Ferry, F., Florescu, S., de Girolamo, G., Gureje, O., Kawakami, N., Lee, S., Navarro-Mateu, F., Piazza, M., Pennell, B.-E., Posada-Villa, J., Ten Have, M., Viana, M.C., Kessler, R.C., 2017. Childhood adversities and post-traumatic stress disorder: evidence for stress sensitisation in the World Mental Health Surveys. Br. J. Psychiatry 211, 280–288. doi:10.1192/bjp.bp.116.197640
    OpenUrlAbstract/FREE Full Text
  41. ↵
    Muneer, A., Minhas, F.A., 2019. Telomere biology in mood disorders: an updated, comprehensive review of the literature. Clin. Psychopharmacol. Neurosci. 17, 343– 363. doi:10.9758/cpn.2019.17.3.343
    OpenUrlCrossRef
  42. ↵
    Napolitano, F., Sirci, F., Carrella, D., di Bernardo, D., 2016. Drug-set enrichment analysis: a novel tool to investigate drug mode of action. Bioinformatics 32, 235–241. doi:10.1093/bioinformatics/btv536
    OpenUrlCrossRefPubMed
  43. ↵
    Ng, B., White, C.C., Klein, H.-U., Sieberts, S.K., McCabe, C., Patrick, E., Xu, J., Yu, L., Gaiteri, C., Bennett, D.A., Mostafavi, S., De Jager, P.L., 2017. An xQTL map integrates the genetic architecture of the human brain’s transcriptome and epigenome. Nat. Neurosci. 20, 1418–1426. doi:10.1038/nn.4632
    OpenUrlCrossRefPubMed
  44. ↵
    Nithipatikom, K., Endsley, M.P., Pfeiffer, A.W., Falck, J.R., Campbell, W.B., 2014. A novel activity of microsomal epoxide hydrolase: metabolism of the endocannabinoid 2-arachidonoylglycerol. J. Lipid Res. 55, 2093–2102. doi:10.1194/jlr.M051284
    OpenUrlAbstract/FREE Full Text
  45. ↵
    Pizzamiglio, L., Focchi, E., Murru, L., Tamborini, M., Passafaro, M., Menna, E., Matteoli, M., Antonucci, F., 2016. New role of ATM in controlling gabaergic tone during development. Cereb. Cortex 26, 3879–3888. doi:10.1093/cercor/bhw125
    OpenUrlCrossRefPubMed
  46. ↵
    Puhlmann, L.M.C., Valk, S.L., Engert, V., Bernhardt, B.C., Lin, J., Epel, E.S., Vrticka, P., Singer, T., 2019. Association of Short-term Change in Leukocyte Telomere Length With Cortical Thickness and Outcomes of Mental Training Among Healthy Adults: A Randomized Clinical Trial. JAMA Netw. Open 2, e199687. doi:10.1001/jamanetworkopen.2019.9687
    OpenUrlCrossRef
  47. ↵
    Qi, T., Wu, Y., Zeng, J., Zhang, F., Xue, A., Jiang, L., Zhu, Z., Kemper, K., Yengo, L., Zheng, Z., eQTLGen Consortium, Marioni, R.E., Montgomery, G.W., Deary, I.J., Wray, N.R., Visscher, P.M., McRae, A.F., Yang, J., 2018. Identifying gene targets for brain-related traits using transcriptomic and methylomic data from blood. Nat. Commun. 9, 2282. doi:10.1038/s41467-018-04558-1
    OpenUrlCrossRef
  48. ↵
    Rajarajan, P., Borrman, T., Liao, W., Schrode, N., Flaherty, E., Casiño, C., Powell, S., Yashaswini, C., LaMarca, E.A., Kassim, B., Javidfar, B., Espeso-Gil, S., Li, A., Won, H., Geschwind, D.H., Ho, S.-M., MacDonald, M., Hoffman, G.E., Roussos, P., Zhang, B., Hahn, C.-G., Weng, Z., Brennand, K.J., Akbarian, S., 2018. Neuron-specific signatures in the chromosomal connectome associated with schizophrenia risk. Science 362. doi:10.1126/science.aat4311
    OpenUrlAbstract/FREE Full Text
  49. ↵
    Redei, E.E., Andrus, B.M., Kwasny, M.J., Seok, J., Cai, X., Ho, J., Mohr, D.C., 2014. Blood transcriptomic biomarkers in adult primary care patients with major depressivedisorder undergoing cognitive behavioral therapy. Transl. Psychiatry 4, e442. doi:10.1038/tp.2014.66
    OpenUrlCrossRef
  50. ↵
    Roesch, Z.K., Tadi, P., 2020. Neuroanatomy, Fourth Ventricle, in: StatPearls. StatPearls Publishing, Treasure Island (FL).
  51. ↵
    Russo, P., Prinzi, G., Proietti, S., Lamonaca, P., Frustaci, A., Boccia, S., Amore, R., Lorenzi, M., Onder, G., Marzetti, E., Valdiglesias, V., Guadagni, F., Valente, M.G., Cascio, G.L., Fraietta, S., Ducci, G., Bonassi, S., 2018. Shorter telomere length in schizophrenia: Evidence from a real-world population and meta-analysis of most recent literature. Schizophr. Res. 202, 37–45. doi:10.1016/j.schres.2018.07.015
    OpenUrlCrossRef
  52. ↵
    Sey, N.Y.A., Hu, B., Mah, W., Fauni, H., McAfee, J.C., Rajarajan, P., Brennand, K.J., Akbarian, S., Won, H., 2020. A computational tool (H-MAGMA) for improved prediction of brain-disorder risk genes by incorporating brain chromatin interaction profiles.Nat. Neurosci. 23, 583–593. doi:10.1038/s41593-020-0603-0
    OpenUrlCrossRef
  53. ↵
    Smith, B.M., Giddens, M.M., Neil, J., Owino, S., Nguyen, T.T., Duong, D., Li, F., Hall, R.A., 2017. Mice lacking Gpr37 exhibit decreased expression of the myelin-associated glycoprotein MAG and increased susceptibility to demyelination. Neuroscience 358, 49–57. doi:10.1016/j.neuroscience.2017.06.006
    OpenUrlCrossRef
  54. ↵
    Sriram, C.S., Jangra, A., Gurjar, S.S., Mohan, P., Bezbaruah, B.K., 2016. Edaravone abrogates LPS-induced behavioral anomalies, neuroinflammation and PARP-1. Physiol. Behav. 154, 135–144. doi:10.1016/j.physbeh.2015.10.029
    OpenUrlCrossRefPubMed
  55. ↵
    Starnawska, A., Tan, Q., McGue, M., Mors, O., Børglum, A.D., Christensen, K., Nyegaard, M., Christiansen, L., 2017. Epigenome-Wide Association Study of Cognitive Functioning in Middle-Aged Monozygotic Twins. Front. Aging Neurosci. 9, 413. doi:10.3389/fnagi.2017.00413
    OpenUrlCrossRefPubMed
  56. ↵
    Szebeni, A., Szebeni, K., DiPeri, T., Chandley, M.J., Crawford, J.D., Stockmeier, C.A., Ordway, G.A., 2014. Shortened telomere length in white matter oligodendrocytes in major depression: potential role of oxidative stress. Int. J. Neuropsychopharmacol. 17, 1579–1589. doi:10.1017/S1461145714000698
    OpenUrlCrossRefPubMed
  57. ↵
    Tapia-Laliena, M.Á., Korzeniewski, N., Peña-Llopis, S., Scholl, C., Fröhling, S., Hohenfellner, M., Duensing, A., Duensing, S., 2019. Cullin 5 is a novel candidate tumor suppressor in renal cell carcinoma involved in the maintenance of genome stability. Oncogenesis 8, 4. doi:10.1038/s41389-018-0110-2
    OpenUrlCrossRef
  58. ↵
    Toth, A.B., Shum, A.K., Prakriya, M., 2016. Regulation of neurogenesis by calcium signaling. Cell Calcium 59, 124–134. doi:10.1016/j.ceca.2016.02.011
    OpenUrlCrossRefPubMed
  59. Václavíková, R., Hughes, D.J., Soucek, P., 2015. Microsomal epoxide hydrolase 1 (EPHX1): Gene, structure, function, and role in human disease. Gene 571, 1–8.doi:10.1016/j.gene.2015.07.071
    OpenUrlCrossRefPubMed
  60. ↵
    Vakonaki, E., Tsiminikaki, K., Plaitis, S., Fragkiadaki, P., Tsoukalas, D., Katsikantami, I., Vaki, G., Tzatzarakis, M.N., Spandidos, D.A., Tsatsakis, A.M., 2018. Common mental disorders and association with telomere length. Biomed. Rep. 8, 111–116. doi:10.3892/br.2018.1040
    OpenUrlCrossRef
  61. ↵
    Walton, C.C., Zhang, W., Patiño-Parrado, I., Barrio-Alonso, E., Garrido, J.-J., Frade, J.M., 2019. Primary neurons can enter M-phase. Sci. Rep. 9, 4594. doi:10.1038/s41598-019-40462-4
    OpenUrlCrossRef
  62. Wang, X., Sundquist, K., Hedelius, A., Palmér, K., Memon, A.A., Sundquist, J., 2017.Leukocyte telomere length and depression, anxiety and stress and adjustment disorders in primary health care patients. BMC Psychiatry 17, 148. doi:10.1186/s12888-017-1308-0
    OpenUrlCrossRef
  63. ↵
    Watanabe, K., Taskesen, E., van Bochoven, A., Posthuma, D., 2017. Functional mapping and annotation of genetic associations with FUMA. Nat. Commun. 8, 1826. doi:10.1038/s41467-017-01261-5
    OpenUrlCrossRefPubMed
  64. ↵
    Whedon, J.M., Glassey, D., 2009. Cerebrospinal fluid stasis and its clinical significance. Altern Ther Health Med 15, 54–60.
    OpenUrlPubMed
  65. ↵
    Wolfe, D., Dudek, S., Ritchie, M.D., Pendergrass, S.A., 2013. Visualizing genomic information across chromosomes with PhenoGram. BioData Min 6, 18. doi:10.1186/1756-0381-6-18
    OpenUrlCrossRefPubMed
  66. ↵
    Yucel, K., Nazarov, A., Taylor, V.H., Macdonald, K., Hall, G.B., Macqueen, G.M., 2013. Cerebellar vermis volume in major depressive disorder. Brain Struct. Funct. 218, 851–858. doi:10.1007/s00429-012-0433-2
    OpenUrlCrossRef
  67. ↵
    Zhang, Y., You, X., Li, S., Long, Q., Zhu, Y., Teng, Z., Zeng, Y., 2020. Peripheral Blood Leukocyte RNA-Seq Identifies a Set of Genes Related to Abnormal Psychomotor Behavior Characteristics in Patients with Schizophrenia. Med. Sci. Monit. 26, e922426. doi:10.12659/MSM.922426
    OpenUrlCrossRef
  68. ↵
    Zhao, B., Luo, T., Li, T., Li, Y., Zhang, J., Shan, Y., Wang, X., Yang, L., Zhou, F., Zhu, Z., Alzheimer’s Disease Neuroimaging Initiative, Pediatric Imaging, Neurocognition and Genetics, Zhu, H., 2019. Genome-wide association analysis of 19,629 individuals identifies variants influencing regional brain volumes and refines their genetic co-architecture with cognitive and mental health traits. Nat. Genet. 51, 1637–1644. doi:10.1038/s41588-019-0516-6
    OpenUrlCrossRefPubMed
  69. Zhao, T., Hu, Y., Zang, T., Wang, Y., 2019. Integrate GWAS, eQTL, and mQTL Data to Identify Alzheimer’s Disease-Related Genes. Front. Genet. 10, 1021. doi:10.3389/fgene.2019.01021
    OpenUrlCrossRef
  70. ↵
    Zhou, G., Soufan, O., Ewald, J., Hancock, R.E.W., Basu, N., Xia, J., 2019. NetworkAnalyst 3.0: a visual analytics platform for comprehensive gene expression profiling and meta-analysis. Nucleic Acids Res. 47, W234–W241. doi:10.1093/nar/gkz240
    OpenUrlCrossRef
  71. ↵
    Zhou, Q.-G., Liu, M.-Y., Lee, H.-W., Ishikawa, F., Devkota, S., Shen, X.-R., Jin, Xin, Wu, H.-Y., Liu, Z., Liu, X., Jin, Xun, Zhou, H.-H., Ro, E.J., Zhang, J., Zhang, Y., Lin, Y.-H., Suh, H., Zhu, D.-Y., 2017. Hippocampal TERT Regulates Spatial Memory Formation through Modulation of Neural Development. Stem Cell Rep. 9, 543–556. doi:10.1016/j.stemcr.2017.06.014
    OpenUrlCrossRef
  72. ↵
    Zhu, Z., Zhang, F., Hu, H., Bakshi, A., Robinson, M.R., Powell, J.E., Montgomery, G.W., Goddard, M.E., Wray, N.R., Visscher, P.M., Yang, J., 2016. Integration of summary data from GWAS and eQTL studies predicts complex trait gene targets. Nat. Genet. 48, 481–487. doi:10.1038/ng.3538
    OpenUrlCrossRefPubMed
Back to top
Posted August 04, 2020.
Download PDF
Data/Code
Email
Shared loci of telomere length with brain morphology, and pleiotropy in transcriptomic and epigenomic profiles of brain
Gita A Pathak, Frank R Wendt, Daniel F Levey, Adam P Mecca, Christopher H van Dyck, Joel Gelernter, Renato Polimanti
medRxiv 2020.08.03.20167726; doi: https://doi.org/10.1101/2020.08.03.20167726
Twitter logo Facebook logo LinkedIn logo Mendeley logo
Citation Tools

Subject Area