Genome-wide DNA methylation analysis of cannabis use disorder in a veteran cohort enriched for posttraumatic stress disorder.

Cannabis use has been increasing over the past decade, not only in the general US population, but particularly among military veterans. With this rise in use has come a concomitant increase in cannabis use disorder (CUD) among veterans. Here, we performed an epigenome-wide association study for lifetime CUD in an Iraq/Afghanistan era veteran cohort enriched for posttraumatic stress disorder (PTSD) comprising 2,310 total subjects (1,109 non-Hispanic black and 1,201 non-Hispanic white). We also investigated CUD interactions with current PTSD status and examined potential indirect effects of DNA methylation (DNAm) on the relationship between CUD and psychiatric diagnoses. Four CpGs were associated with lifetime CUD, even after controlling for the effects of current smoking (AHRR cg05575921, LINC00299 cg23079012, VWA7 cg22112841, and FAM70A cg08760398). Importantly, cg05575921, a CpG strongly linked to smoking, remained associated with lifetime CUD even when restricting the analysis to veterans who reported never smoking cigarettes. Moreover, CUD interacted with current PTSD to affect cg05575921 and cg23079012 such that those with both CUD and PTSD displayed significantly lower DNAm compared to the other groups. Finally, we provide preliminary evidence that AHRR cg05575921 helps explain the association between CUD and any psychiatric diagnoses, specifically mood disorders.

GWAS Meta-Analysis of Suicide Attempt: Identification of 12 Genome-Wide Significant Loci and Implication of Genetic Risks for Specific Health Factors.

OBJECTIVE: Suicidal behavior is heritable and is a major cause of death worldwide. Two large-scale genome-wide association studies (GWASs) recently discovered and cross-validated genome-wide significant (GWS) loci for suicide attempt (SA). The present study leveraged the genetic cohorts from both studies to conduct the largest GWAS meta-analysis of SA to date. Multi-ancestry and admixture-specific meta-analyses were conducted within groups of significant African, East Asian, and European ancestry admixtures. METHODS: This study comprised 22 cohorts, including 43,871 SA cases and 915,025 ancestry-matched controls. Analytical methods across multi-ancestry and individual ancestry admixtures included inverse variance-weighted fixed-effects meta-analyses, followed by gene, gene-set, tissue-set, and drug-target enrichment, as well as summary-data-based Mendelian randomization with brain expression quantitative trait loci data, phenome-wide genetic correlation, and genetic causal proportion analyses. RESULTS: Multi-ancestry and European ancestry admixture GWAS meta-analyses identified 12 risk loci at p values <5×10-8. These loci were mostly intergenic and implicated DRD2, SLC6A9, FURIN, NLGN1, SOX5, PDE4B, and CACNG2. The multi-ancestry SNP-based heritability estimate of SA was 5.7% on the liability scale (SE=0.003, p=5.7×10-80). Significant brain tissue gene expression and drug set enrichment were observed. There was shared genetic variation of SA with attention deficit hyperactivity disorder, smoking, and risk tolerance after conditioning SA on both major depressive disorder and posttraumatic stress disorder. Genetic causal proportion analyses implicated shared genetic risk for specific health factors. CONCLUSIONS: This multi-ancestry analysis of suicide attempt identified several loci contributing to risk and establishes significant shared genetic covariation with clinical phenotypes. These findings provide insight into genetic factors associated with suicide attempt across ancestry admixture populations, in veteran and civilian populations, and in attempt versus death.

Genetic background determines severity of Loxl1-mediated systemic and ocular elastosis in mice.

Pseudoexfoliation syndrome (PEX) is a systemic, age-related disorder characterized by elastosis and extracellular matrix deposits. Its most significant ocular manifestation is an aggressive form of glaucoma associated with variants in the gene encoding lysyl oxidase-like 1 (LOXL1). Depending upon the population, variants in LOXL1 can impart risk or protection for PEX, suggesting the importance of genetic context. As LOXL1 protein levels are lower and the degree of elastosis is higher in people with PEX, we studied Loxl1-deficient mice on three different genetic backgrounds: C57BL/6 (BL/6), 129S×C57BL/6 (50/50) and 129S. Early onset and high prevalence of spontaneous pelvic organ prolapse in BL/6 Loxl1-/- mice necessitated the study of mice that were <2 months old. Similar to pelvic organ prolapse, most elastosis endpoints were the most severe in BL/6 Loxl1-/- mice, including skin laxity, pulmonary tropoelastin accumulation, expansion of Schlemm's canal and dilation of intrascleral veins. Interestingly, intraocular pressure was elevated in 50/50 Loxl1-/- mice, depressed in BL/6 Loxl1-/- mice and unchanged in 129S Loxl1-/- mice compared to that of control littermates. Overall, the 129S background was protective against most elastosis phenotypes studied. Thus, repair of elastin-containing tissues is impacted by the abundance of LOXL1 and genetic context in young animals.

Lysyl oxidase-like 1-antisense 1 (LOXL1-AS1) lncRNA differentially regulates gene and protein expression, signaling and morphology of human ocular cells.

Pseudoexfoliation glaucoma (PEXG) is characterized by dysregulated extracellular matrix (ECM) homeostasis that disrupts conventional outflow function and increases intraocular pressure (IOP). Prolonged IOP elevation results in optic nerve head damage and vision loss. Uniquely, PEXG is a form of open angle glaucoma that has variable penetrance, is difficult to treat and does not respond well to common IOP-lowering pharmaceuticals. Therefore, understanding modulators of disease severity will aid in targeted therapies for PEXG. Genome-wide association studies have identified polymorphisms in the long non-coding RNA lysyl oxidase-like 1-antisense 1 (LOXL1-AS1) as a risk factor for PEXG. Risk alleles, oxidative stress and mechanical stretch all alter LOXL1-AS1 expression. As a long non-coding RNA, LOXL1-AS1 binds hnRNPL and regulates global gene expression. In this study, we focus on the role of LOXL1-AS1 in the ocular cells (trabecular meshwork and Schlemm's canal) that regulate IOP. We show that selective knockdown of LOXL1-AS1 leads to cell-type-specific changes in gene expression, ECM homeostasis, signaling and morphology. These results implicate LOXL1-AS1 as a modulator of cellular homeostasis, altering cell contractility and ECM turnover, both of which are well-known contributors to PEXG. These findings support LOXL1-AS1 as a key target for modifying the disease.

Large epigenome-wide association study identifies multiple novel differentially methylated CpG sites associated with suicidal thoughts and behaviors in veterans.

Introduction: The U.S. suicide mortality rate has steadily increased during the past two decades, particularly among military veterans; however, the epigenetic basis of suicidal thoughts and behaviors (STB) remains largely unknown. Methods: To address this issue, we conducted an epigenome-wide association study of DNA methylation (DNAm) of peripheral blood samples obtained from 2,712 U.S. military veterans. Results: Three DNAm probes were significantly associated with suicide attempts, surpassing the multiple testing threshold (FDR q-value <0.05), including cg13301722 on chromosome 7, which lies between the genes SLC4A2 and CDK5; cg04724646 in PDE3A; and cg04999352 in RARRES3. cg13301722 was also found to be differentially methylated in the cerebral cortex of suicide decedents in a publicly-available dataset (p = 0.03). Trait enrichment analysis revealed that the CpG sites most strongly associated with STB in the present sample were also associated with smoking, alcohol consumption, maternal smoking, and maternal alcohol consumption, whereas pathway enrichment analysis revealed significant associations with circadian rhythm, adherens junction, insulin secretion, and RAP-1 signaling, each of which was recently associated with suicide attempts in a large, independent genome-wide association study of suicide attempts of veterans. Discussion: Taken together, the present findings suggest that SLC4A2, CDK5, PDE3A, and RARRES3 may play a role in STB. CDK5, a member of the cyclin-dependent kinase family that is highly expressed in the brain and essential for learning and memory, appears to be a particularly promising candidate worthy of future study; however, additional work is still needed to replicate these finding in independent samples.

Rare protective variants and glaucoma-relevant cell stressors modulate Angiopoietin-like 7 expression.

Rare missense and nonsense variants in the Angiopoietin-like 7 (ANGPTL7) gene confer protection from primary open-angle glaucoma (POAG), though the functional mechanism remains uncharacterized. Interestingly, a larger variant effect size strongly correlates with in silico predictions of increased protein instability (r = -0.98), suggesting that protective variants lower ANGPTL7 protein levels. Here, we show that missense and nonsense variants cause aggregation of mutant ANGPTL7 protein in the endoplasmic reticulum (ER) and decreased levels of secreted protein in human trabecular meshwork (TM) cells; a lower secreted:intracellular protein ratio strongly correlates with variant effects on intraocular pressure (r = 0.81). Importantly, accumulation of mutant protein in the ER does not increase expression of ER stress proteins in TM cells (P > 0.05 for all variants tested). Cyclic mechanical stress, a glaucoma-relevant physiologic stressor, also significantly lowers ANGPTL7 expression in primary cultures of human Schlemm's canal (SC) cells (-2.4-fold-change, P = 0.01). Collectively, these data suggest that the protective effects of ANGPTL7 variants in POAG stem from lower levels of secreted protein, which may modulate responses to physiologic and pathologic ocular cell stressors. Downregulation of ANGPTL7 expression may therefore serve as a viable preventative and therapeutic strategy for this common, blinding disease.

Genome-wide association study identifies four pan-ancestry loci for suicidal ideation in the Million Veteran Program.

Suicidal ideation (SI) often precedes and predicts suicide attempt and death, is the most common suicidal phenotype and is over-represented in veterans. The genetic architecture of SI in the absence of suicide attempt (SA) is unknown, yet believed to have distinct and overlapping risk with other suicidal behaviors. We performed the first GWAS of SI without SA in the Million Veteran Program (MVP), identifying 99,814 SI cases from electronic health records without a history of SA or suicide death (SD) and 512,567 controls without SI, SA or SD. GWAS was performed separately in the four largest ancestry groups, controlling for sex, age and genetic substructure. Ancestry-specific results were combined via meta-analysis to identify pan-ancestry loci. Four genome-wide significant (GWS) loci were identified in the pan-ancestry meta-analysis with loci on chromosomes 6 and 9 associated with suicide attempt in an independent sample. Pan-ancestry gene-based analysis identified GWS associations with DRD2, DCC, FBXL19, BCL7C, CTF1, ANNK1, and EXD3. Gene-set analysis implicated synaptic and startle response pathways (q's<0.05). European ancestry (EA) analysis identified GWS loci on chromosomes 6 and 9, as well as GWS gene associations in EXD3, DRD2, and DCC. No other ancestry-specific GWS results were identified, underscoring the need to increase representation of diverse individuals. The genetic correlation of SI and SA within MVP was high (rG = 0.87; p = 1.09e-50), as well as with post-traumatic stress disorder (PTSD; rG = 0.78; p = 1.98e-95) and major depressive disorder (MDD; rG = 0.78; p = 8.33e-83). Conditional analysis on PTSD and MDD attenuated most pan-ancestry and EA GWS signals for SI without SA to nominal significance, with the exception of EXD3 which remained GWS. Our novel findings support a polygenic and complex architecture for SI without SA which is largely shared with SA and overlaps with psychiatric conditions frequently comorbid with suicidal behaviors.

Impact of traumatic life events and polygenic risk scores for major depression and posttraumatic stress disorder on Iraq/Afghanistan Veterans.

Traumatic experiences and genetic heritability are among the most widely acknowledged risk factors leading to the development of psychopathology; including posttraumatic stress disorder (PTSD) and major depressive disorder (MDD). The purpose of this study was to investigate if polygenic risk scores (PRS) among Veterans interacted with traumatic stress to predict PTSD and MDD. 1,389 Iraq-Afghanistan military service Veterans from the Mental Illness Research Education and Clinical Center dataset were analyzed. Genome-wide association study (GWAS) statistics were utilized to generate PRS for PTSD (PRSPTSD) and PRS for MDD (PRSMDD) in order to analyze PRS-by-environment (PRSxE) with trauma exposure to predict PTSD and MDD diagnoses. Trauma exposure and PRSPTSD, were independently associated with a current PTSD diagnosis (p < 0.001 and p < 0.001, respectively). The interaction between trauma exposure and PRSMDD to predict a current diagnosis of PTSD trended towards significance (p = 0.053). Stratifying by trauma thresholds, among those within the lowest trauma load, the association of PRSMDD with PTSD was found to be nominally significant (p = 0.03). For a MDD diagnosis, there was a significant association with trauma exposure (p < 0.001); and the association with PRSMDD was found to be nominally significant (p = 0.03). No significant PRSxE effects were found with MDD. Our findings corroborate previous research highlighting trauma exposure, and genetic heritability, as risk factors for the development of PTSD and MDD in a Veteran population. Additionally, findings suggest that genetic vulnerability may be less important as trauma exposure increases, with high levels of trauma likely to result in PTSD and MDD, regardless of genetic vulnerability.

Identification of Novel, Replicable Genetic Risk Loci for Suicidal Thoughts and Behaviors Among US Military Veterans.

Importance: Suicide is a leading cause of death; however, the molecular genetic basis of suicidal thoughts and behaviors (SITB) remains unknown. Objective: To identify novel, replicable genomic risk loci for SITB. Design, Setting, and Participants: This genome-wide association study included 633 778 US military veterans with and without SITB, as identified through electronic health records. GWAS was performed separately by ancestry, controlling for sex, age, and genetic substructure. Cross-ancestry risk loci were identified through meta-analysis. Study enrollment began in 2011 and is ongoing. Data were analyzed from November 2021 to August 2022. Main Outcome and Measures: SITB. Results: A total of 633 778 US military veterans were included in the analysis (57 152 [9%] female; 121 118 [19.1%] African ancestry, 8285 [1.3%] Asian ancestry, 452 767 [71.4%] European ancestry, and 51 608 [8.1%] Hispanic ancestry), including 121 211 individuals with SITB (19.1%). Meta-analysis identified more than 200 GWS (P < 5 × 10-8) cross-ancestry risk single-nucleotide variants for SITB concentrated in 7 regions on chromosomes 2, 6, 9, 11, 14, 16, and 18. Top single-nucleotide variants were largely intronic in nature; 5 were independently replicated in ISGC, including rs6557168 in ESR1, rs12808482 in DRD2, rs77641763 in EXD3, rs10671545 in DCC, and rs36006172 in TRAF3. Associations for FBXL19 and AC018880.2 were not replicated. Gene-based analyses implicated 24 additional GWS cross-ancestry risk genes, including FURIN, TSNARE1, and the NCAM1-TTC12-ANKK1-DRD2 gene cluster. Cross-ancestry enrichment analyses revealed significant enrichment for expression in brain and pituitary tissue, synapse and ubiquitination processes, amphetamine addiction, parathyroid hormone synthesis, axon guidance, and dopaminergic pathways. Seven other unique European ancestry-specific GWS loci were identified, 2 of which (POM121L2 and METTL15/LINC02758) were replicated. Two additional GWS ancestry-specific loci were identified within the African ancestry (PET112/GATB) and Hispanic ancestry (intergenic locus on chromosome 4) subsets, both of which were replicated. No GWS loci were identified within the Asian ancestry subset; however, significant enrichment was observed for axon guidance, cyclic adenosine monophosphate signaling, focal adhesion, glutamatergic synapse, and oxytocin signaling pathways across all ancestries. Within the European ancestry subset, genetic correlations (r > 0.75) were observed between the SITB phenotype and a suicide attempt-only phenotype, depression, and posttraumatic stress disorder. Additionally, polygenic risk score analyses revealed that the Million Veteran Program polygenic risk score had nominally significant main effects in 2 independent samples of veterans of European and African ancestry. Conclusions and Relevance: The findings of this analysis may advance understanding of the molecular genetic basis of SITB and provide evidence for ESR1, DRD2, TRAF3, and DCC as cross-ancestry candidate risk genes. More work is needed to replicate these findings and to determine if and how these genes might impact clinical care.

The social construction of genomics and genetic analysis in ocular diseases in Ibadan, South-western Nigeria.

Genomics, an emerging field to improve public health practice, has potential benefits to understanding ocular diseases. This study explored the social construction of genomics in ocular diseases in the blind community in Ibadan, Nigeria, through two focus group discussions and twelve in-depth interview sessions conducted among people living with ocular disorders. The data were thematic and content-analysed. Although the participants had limited knowledge about ocular diseases, genomics, and their nexus, they maintained a positive attitude toward its potential benefits. This informed their willingness to participate in genomics testing for ocular diseases. The participants preferred saliva-based sample collection over blood-based, and expressed concern for the procedure and accrued benefits of genomics studies. Thus, public sensitisation about ocular diseases and client-centred genomics testing procedures should be engendered.

A genome-wide association study of suicide attempts in the million veterans program identifies evidence of pan-ancestry and ancestry-specific risk loci.

To identify pan-ancestry and ancestry-specific loci associated with attempting suicide among veterans, we conducted a genome-wide association study (GWAS) of suicide attempts within a large, multi-ancestry cohort of U.S. veterans enrolled in the Million Veterans Program (MVP). Cases were defined as veterans with a documented history of suicide attempts in the electronic health record (EHR; N = 14,089) and controls were defined as veterans with no documented history of suicidal thoughts or behaviors in the EHR (N = 395,064). GWAS was performed separately in each ancestry group, controlling for sex, age and genetic substructure. Pan-ancestry risk loci were identified through meta-analysis and included two genome-wide significant loci on chromosomes 20 (p = 3.64 × 10-9) and 1 (p = 3.69 × 10-8). A strong pan-ancestry signal at the Dopamine Receptor D2 locus (p = 1.77 × 10-7) was also identified and subsequently replicated in a large, independent international civilian cohort (p = 7.97 × 10-4). Additionally, ancestry-specific genome-wide significant loci were also detected in African-Americans, European-Americans, Asian-Americans, and Hispanic-Americans. Pathway analyses suggested over-representation of many biological pathways with high clinical significance, including oxytocin signaling, glutamatergic synapse, cortisol synthesis and secretion, dopaminergic synapse, and circadian rhythm. These findings confirm that the genetic architecture underlying suicide attempt risk is complex and includes both pan-ancestry and ancestry-specific risk loci. Moreover, pathway analyses suggested many commonly impacted biological pathways that could inform development of improved therapeutics for suicide prevention.

Diversity in Polygenic Risk of Primary Open-Angle Glaucoma.

Glaucoma is the leading cause of irreversible blindness worldwide. Primary open-angle glaucoma (POAG), the most common glaucoma subtype, is more prevalent and severe in individuals of African ancestry. Unfortunately, this ancestral group has been historically under-represented among genetic studies of POAG. Moreover, both genetic and polygenic risk scores (GRS, PRS) that are typically based on genetic data from European-descent populations are not transferable to individuals without a majority of European ancestry. Given the aspirations of leveraging genetic information for precision medicine, GRS and PRS demonstrate clinical potential but fall short, in part due to the lack of diversity in these studies. Prioritizing diversity in the discovery of risk variants will improve the performance and utility of GRS and PRS-derived risk estimation for disease stratification, which could bring about earlier POAG intervention and treatment for a disease that often goes undetected until significant damage has occurred.

Identification of Estrogen Signaling in a Prioritization Study of Intraocular Pressure-Associated Genes.

Elevated intraocular pressure (IOP) is the only modifiable risk factor for primary open-angle glaucoma (POAG). Herein we sought to prioritize a set of previously identified IOP-associated genes using novel and previously published datasets. We identified several genes for future study, including several involved in cytoskeletal/extracellular matrix reorganization, cell adhesion, angiogenesis, and TGF-ß signaling. Our differential correlation analysis of IOP-associated genes identified 295 pairs of 201 genes with differential correlation. Pathway analysis identified ß-estradiol as the top upstream regulator of these genes with ESR1 mediating 25 interactions. Several genes (i.e., EFEMP1, FOXC1, and SPTBN1) regulated by ß-estradiol/ESR1 were highly expressed in non-glaucomatous human trabecular meshwork (TM) or Schlemm's canal (SC) cells and specifically expressed in TM/SC cell clusters defined by single-cell RNA-sequencing. We confirmed ESR1 gene and protein expression in human TM cells and TM/SC tissue with quantitative real-time PCR and immunofluorescence, respectively. 17ß-estradiol was identified in bovine, porcine, and human aqueous humor (AH) using ELISA. In conclusion, we have identified estrogen receptor signaling as a key modulator of several IOP-associated genes. The expression of ESR1 and these IOP-associated genes in TM/SC tissue and the presence of 17ß-estradiol in AH supports a role for estrogen signaling in IOP regulation.

Gene Expression Analysis in Three Posttraumatic Stress Disorder Cohorts Implicates Inflammation and Innate Immunity Pathways and Uncovers Shared Genetic Risk With Major Depressive Disorder.

Posttraumatic stress disorder (PTSD) is a complex psychiatric disorder that can develop following exposure to traumatic events. The Psychiatric Genomics Consortium PTSD group (PGC-PTSD) has collected over 20,000 multi-ethnic PTSD cases and controls and has identified both genetic and epigenetic factors associated with PTSD risk. To further investigate biological correlates of PTSD risk, we examined three PGC-PTSD cohorts comprising 977 subjects to identify differentially expressed genes among PTSD cases and controls. Whole blood gene expression was quantified with the HumanHT-12 v4 Expression BeadChip for 726 OEF/OIF veterans from the Veterans Affairs (VA) Mental Illness Research Education and Clinical Center (MIRECC), 155 samples from the Injury and Traumatic Stress (INTRuST) Clinical Consortium, and 96 Australian Vietnam War veterans. Differential gene expression analysis was performed in each cohort separately followed by meta-analysis. In the largest cohort, we performed co-expression analysis to identify modules of genes that are associated with PTSD and MDD. We then conducted expression quantitative trait loci (eQTL) analysis and assessed the presence of eQTL interactions involving PTSD and major depressive disorder (MDD). Finally, we utilized PTSD and MDD GWAS summary statistics to identify regions that colocalize with eQTLs. Although not surpassing correction for multiple testing, the most differentially expressed genes in meta-analysis were interleukin-1 beta (IL1B), a pro-inflammatory cytokine previously associated with PTSD, and integrin-linked kinase (ILK), which is highly expressed in brain and can rescue dysregulated hippocampal neurogenesis and memory deficits. Pathway analysis revealed enrichment of toll-like receptor (TLR) and interleukin-1 receptor genes, which are integral to cellular innate immune response. Co-expression analysis identified four modules of genes associated with PTSD, two of which are also associated with MDD, demonstrating common biological pathways underlying the two conditions. Lastly, we identified four genes (UBA7, HLA-F, HSPA1B, and RERE) with high probability of a shared causal eQTL variant with PTSD and/or MDD GWAS variants, thereby providing a potential mechanism by which the GWAS variant contributes to disease risk. In summary, we provide additional evidence for genes and pathways previously reported and identified plausible novel candidates for PTSD. These data provide further insight into genetic factors and pathways involved in PTSD, as well as potential regions of pleiotropy between PTSD and MDD.

Examining Individual and Synergistic Contributions of PTSD and Genetics to Blood Pressure: A Trans-Ethnic Meta-Analysis.

Growing research suggests that posttraumatic stress disorder (PTSD) may be a risk factor for poor cardiovascular health, and yet our understanding of who might be at greatest risk of adverse cardiovascular outcomes after trauma is limited. In this study, we conducted the first examination of the individual and synergistic contributions of PTSD symptoms and blood pressure genetics to continuous blood pressure levels. We harnessed the power of the Psychiatric Genomics Consortium-PTSD Physical Health Working Group and investigated these associations across 11 studies of 72,224 trauma-exposed individuals of European (n = 70,870) and African (n = 1,354) ancestry. Genetic contributions to blood pressure were modeled via polygenic scores (PGS) for systolic blood pressure (SBP) and diastolic blood pressure (DBP) that were derived from a prior trans-ethnic blood pressure genome-wide association study (GWAS). Results of trans-ethnic meta-analyses revealed significant main effects of the PGS on blood pressure levels [SBP: ß = 2.83, standard error (SE) = 0.06, p < 1E-20; DBP: ß = 1.32, SE = 0.04, p < 1E-20]. Significant main effects of PTSD symptoms were also detected for SBP and DBP in trans-ethnic meta-analyses, though there was significant heterogeneity in these results. When including data from the largest contributing study - United Kingdom Biobank - PTSD symptoms were negatively associated with SBP levels (ß = -1.46, SE = 0.44, p = 9.8E-4) and positively associated with DBP levels (ß = 0.70, SE = 0.26, p = 8.1E-3). However, when excluding the United Kingdom Biobank cohort in trans-ethnic meta-analyses, there was a nominally significant positive association between PTSD symptoms and SBP levels (ß = 2.81, SE = 1.13, p = 0.01); no significant association was observed for DBP (ß = 0.43, SE = 0.78, p = 0.58). Blood pressure PGS did not significantly moderate the associations between PTSD symptoms and blood pressure levels in meta-analyses. Additional research is needed to better understand the extent to which PTSD is associated with high blood pressure and how genetic as well as contextual factors may play a role in influencing cardiovascular risk.

Epigenome-wide meta-analysis of PTSD across 10 military and civilian cohorts identifies methylation changes in AHRR.

Epigenetic differences may help to distinguish between PTSD cases and trauma-exposed controls. Here, we describe the results of the largest DNA methylation meta-analysis of PTSD to date. Ten cohorts, military and civilian, contribute blood-derived DNA methylation data from 1,896 PTSD cases and trauma-exposed controls. Four CpG sites within the aryl-hydrocarbon receptor repressor (AHRR) associate with PTSD after adjustment for multiple comparisons, with lower DNA methylation in PTSD cases relative to controls. Although AHRR methylation is known to associate with smoking, the AHRR association with PTSD is most pronounced in non-smokers, suggesting the result was independent of smoking status. Evaluation of metabolomics data reveals that AHRR methylation associated with kynurenine levels, which are lower among subjects with PTSD. This study supports epigenetic differences in those with PTSD and suggests a role for decreased kynurenine as a contributor to immune dysregulation in PTSD.

Genetic predictors of hippocampal subfield volume in PTSD cases and trauma-exposed controls.

Behavioural, structural, and functional neuroimaging have implicated the hippocampus as a critical brain region in posttraumatic stress disorder (PTSD) pathogenesis. Recent work in a normative, primarily European, sample identified 15 unique genetic loci contributing to structural variability in six hippocampal subfield volumes. We explored the relevance of these loci in two samples (Mental Illness Research Education and Clinical Centre [MIRECC] and Grady; n = 290) of trauma-exposed individuals enriched for PTSD and of diverse ancestry. Four of the previous loci demonstrated nominal evidence of replication in the MIRECC dataset, primarily within non-Hispanic whites (NHW). One locus replicated in the Grady cohort, which was composed exclusively of non-Hispanic blacks (NHB). Our data supported genetic interactions with diagnosis of lifetime PTSD and genetic interactions with childhood trauma in the MIRECC sample, but not the Grady sample. Given the racial, diagnostic, and trauma-exposure differences with the original genome-wide association study (GWAS) report, we conducted a full GWAS in the MIRECC and Grady datasets. Interactions between genetic variants and lifetime PTSD or childhood trauma were interrogated for single nucleotide polymorphisms (SNPs) with evidence of main effects. Genetic associations surpassed false discovery rate (FDR)-correction within hippocampal subfields in fimbria, subiculum, cornu ammonis-1 (CA1), and hippocampal amygdala transition area (HATA). One association was replicated in the Grady cohort (rs12880795 in TUNAR with left (L)-HATA volume). The most significant association in the MIRECC dataset was between rs6906714 in LINC02571 and right (R)-fimbria volume (p = 5.99×10-8, q = 0.0056). Interestingly, the effect of rs6906714 on R-fimbria volume increased with exposure to childhood trauma (gene*environment [G*E] interaction p = 0.022). These preliminary results argue for G*E interactions between genetic loci with PTSD and childhood trauma on hippocampal phenotypes. Our results underscore the need for larger neuroimaging-genetic studies in PTSD, trauma, and ancestrally diverse populations.

Las neuroimagenologia, conductual, estructural y funcional ha implicado que el hipocampo se constituye como una región cerebral critica en la patogénesis del trastorno de estrés postraumático (TEPT). Un trabajo reciente, realizado en una muestra normativa, principalmente europea, identifico 15 loci genéticos únicos que contribuyen a la variabilidad estructural de seis volúmenes de subcampos hipocampales. Exploramos la relevancia de estos loci en dos muestras enriquecidas para TEPT (del Centro de Educación y Clínica Investigación sobre Enfermedades Mentales, MIRECC por sus siglas en inglés y Grady; n=290) de individuos expuestos a trauma y de ascendencia diversa. Cuatro de los loci previos demostraron evidencia nominal de replicación en la base de datos MIRECC, principalmente en personas de raza blanca no hispánicos (NHW). Se replicó un locus en la cohorte Grady, que estaba compuesta exclusivamente por personas de raza negra no hispánicos (NHB). Nuestros datos respaldaron las interacciones genéticas con el diagnostico de TEPT a lo largo de la vida e interacciones genéticas con trauma infantil en la muestra MIRECC, pero no en la de Grady. Debido a las diferencias raciales, diagnósticas y de exposición a trauma con el reporte original del estudio de asociación del genoma completo (GWAS por sus siglas en ingles), realizamos un GWAS completo con la base de datos MIRECC y Grady. Se exploraron polimorfismos de nucleótido único (SNPs por sus siglas en ingles) en las interacciones de variantes genéticas del TEPT a lo largo de la vida y del trauma infantil, con evidencia de un efecto genético principal. Las asociaciones genéticas sobrepasaron a la corrección de tasa de falso descubrimiento (FDR) dentro de los subcampos hipocampales de la fimbria, subículo, asta de Amon-1 (CA1), y el área de transición hipocampo-amigdaliana (HATA). Una asociación se replicó en la cohorte de Grady (rs 12880795 en TUNAR con volumen izquierdo del HATA). La asociación más significativa en la base de datos de MIRECC estuvo entre rs6906714 en LINC02571 y el volumen de la fimbria derecha (p=5.99×10-8, q=0.0056). Interesantemente, el efecto rs6906714 sobre el volumen de la fimbria derecha se incrementaba con la exposición al trauma infantil (interacción gen* ambiente [G*A] p=0.022). Estos resultados preliminares orientarían a la presencia de interacciones G*A de loci genéticos con el TEPT y trauma infantil en fenotipos del hipocampo. Nuestro resultados destacan la necesidad de estudios más grandes que vinculen neuroimagenes y genética en poblaciones con TEPT, trauma, y genealogía diversa.

Differential DNA methylation patterns in human Schlemm's canal endothelial cells with glaucoma.

Purpose: Schlemm's canal (SC) endothelial cells derived from donors with or without glaucoma showed different mechanical properties and gene expression. As an important contributor to the regulation of intraocular pressure (IOP) and pathogenesis of primary open-angle glaucoma (POAG), the heritable key epigenetic changes, methylation may play an important role in the physiologic function of SC cells. This study aims to identify differentially methylated CpG sites (DMSs) in primary cultures of human SC cells with or without glaucoma. Methods: We examined the methylation pattern of seven strains of primary human cells (two glaucoma and five normal SC cell samples), which were isolated and characterized using established protocols. DNA methylation was profiled using Illumina Human Methylation 450 BeadChip. Raw data were extracted and exported using Illumina GenomeStudio software. After quantile normalization, DNA methylation data were analyzed using R package RnBeads in Bioconductor. DMSs were filtered with p ≤ 1E-5, methylation change ≥ 0.1, and false discovery rate ≤ 0.05. The closest genes and the location of each CpG site were annotated using R package FDb.InfiniumMethylation.hg19. Gene Ontology and pathway analysis was performed using WebGestalt. Selected DMSs were validated using the Zymo qMethyl kit. Results: We used five non-glaucoma and two glaucomatous SC cell samples to profile genome-wide DNA methylation using Illumina Infinium Methylation BeadChips. Principle component analysis showed the separation between the glaucoma and control samples. After quality control and differential analysis, we identified 298 highly significant DMSs (p ≤ 1E-5). Among them, 221 DMSs were within 1 kb of a nearby gene. Gene Ontology analysis demonstrated significant enrichment in positive regulation of cell migration, negative regulation of endothelial cell proliferation, and stress fiber and actin filament bundles. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis showed enrichment in cell adhesion and gap junctions. Several glaucoma-related genes were identified, including TGFBR3, THBS1, PITX2, DAXX, TBX3, TNXB, ANGPT1, and PLEKHA7. We also examined differentially methylated regions (DMRs) near these CpG sites and identified significant DMRs in TBX3, TNXB1, DAXX, and PITX2. Conclusions: This study represents the first genome-wide DNA methylation profiling in cultured human primary SC cells. The DMSs were enriched in the pathways related to outflow resistance. Several DMRs were validated in glaucoma-associated genes, further suggesting the role of DNA methylation in glaucoma development. This study could provide comprehensive understanding of DNA methylation in glaucoma and its effect on aqueous humor outflow.

Integral role for lysyl oxidase-like-1 in conventional outflow tissue function and behavior.

Lysyl oxidase-like-1 (LOXL1), a vital crosslinking enzyme in elastin fiber maintenance, is essential for the stability and strength of elastic vessels and tissues. Variants in the LOXL1 locus associate with a dramatic increase in risk of exfoliation syndrome (XFS), a systemic fibrillopathy, which often presents with ocular hypertension and exfoliation glaucoma (XFG). We examined the role of LOXL1 in conventional outflow function, the prime regulator of intraocular pressure (IOP). Using Loxl1-/- , Loxl1+/- , and Loxl1+/+ mice, we observed an inverse relationship between LOXL1 expression and IOP, which worsened with age. Elevated IOP in Loxl1-/- mice was associated with a larger globe, decreased ocular compliance, increased outflow facility, extracellular matrix (ECM) abnormalities, and dilated intrascleral veins, yet, no dilation of arteries or capillaries. Interestingly, in living Loxl1-/- mouse eyes, Schlemm's canal (SC) was less susceptible to collapse when challenged with acute elevations in IOP, suggesting elevated episcleral venous pressure (EVP). Thus, LOXL1 expression is required for normal IOP control, while ablation results in altered ECM repair/homeostasis and conventional outflow physiology. Dilation of SC and distal veins, but not arteries, is consistent with key structural and functional roles for elastin in low-pressure vessels subjected to cyclical mechanical stress.

An epigenome-wide association study of posttraumatic stress disorder in US veterans implicates several new DNA methylation loci.

BACKGROUND: Previous studies using candidate gene and genome-wide approaches have identified epigenetic changes in DNA methylation (DNAm) associated with posttraumatic stress disorder (PTSD). METHODS: In this study, we performed an EWAS of PTSD in a cohort of Veterans (n = 378 lifetime PTSD cases and 135 controls) from the Translational Research Center for TBI and Stress Disorders (TRACTS) cohort assessed using the Illumina EPIC Methylation BeadChip which assesses DNAm at more than 850,000 sites throughout the genome. Our model included covariates for ancestry, cell heterogeneity, sex, age, and a smoking score based on DNAm at 39 smoking-associated CpGs. We also examined in EPIC-based DNAm data generated from pre-frontal cortex (PFC) tissue from the National PTSD Brain Bank (n = 72). RESULTS: The analysis of blood samples yielded one genome-wide significant association with PTSD at cg19534438 in the gene G0S2 (p = 1.19 × 10-7, padj = 0.048). This association was replicated in an independent PGC-PTSD-EWAS consortium meta-analysis of military cohorts (p = 0.0024). We also observed association with the smoking-related locus cg05575921 in AHRR despite inclusion of a methylation-based smoking score covariate (p = 9.16 × 10-6), which replicates a previously observed PGC-PTSD-EWAS association (Smith et al. 2019), and yields evidence consistent with a smoking-independent effect. The top 100 EWAS loci were then examined in the PFC data. One of the blood-based PTSD loci, cg04130728 in CHST11, which was in the top 10 loci in blood, but which was not genome-wide significant, was significantly associated with PTSD in brain tissue (in blood p = 1.19 × 10-5, padj = 0.60, in brain, p = 0.00032 with the same direction of effect). Gene set enrichment analysis of the top 500 EWAS loci yielded several significant overlapping GO terms involved in pathogen response, including "Response to lipopolysaccharide" (p = 6.97 × 10-6, padj = 0.042). CONCLUSIONS: The cross replication observed in independent cohorts is evidence that DNA methylation in peripheral tissue can yield consistent and replicable PTSD associations, and our results also suggest that that some PTSD associations observed in peripheral tissue may mirror associations in the brain.