A genome-wide association study of alloimmunization in the TOPMed OMG-SCD cohort identifies a locus on chromosome 12.

BACKGROUND: Red cell alloimmunization after exposure to donor red cells is a very common complication of transfusion for patients with sickle cell disease (SCD), resulting frequently in accelerated donor red blood cell destruction. Patients show substantial differences in their predisposition to alloimmunization, and genetic variability is one proposed component. Although several genetic association studies have been conducted for alloimmunization, the results have been inconsistent, and the genetic determinants of alloimmunization remain largely unknown. STUDY DESIGN AND METHODS: We performed a genome-wide association study (GWAS) in 236 African American (AA) SCD patients from the Outcome Modifying Genes in Sickle Cell Disease (OMG-SCD) cohort, which is part of Trans-Omics for Precision Medicine (TOPMed), with whole-genome sequencing data available. We also performed sensitivity analyses adjusting for different sets of covariates and applied different sample grouping strategies based on the number of alloantibodies patients developed. RESULTS: We identified one genome-wide significant locus on chr12 (p = 3.1e-9) with no evidence of genomic inflation (lambda = 1.003). Further leveraging QTL evidence from GTEx whole blood and/or Jackson Heart Study PBMC RNA-Seq data, we identified a number of potential genes, such as ARHGAP9, STAT6, and ATP23, that may be driving the association signal. We also discovered some suggestive loci using different analysis strategies. DISCUSSION: We call for the community to collect additional alloantibody information within SCD cohorts to further the understanding of the genetic basis of alloimmunization in order to improve transfusion outcomes.

Single-nucleus multi-omics of Parkinson's disease reveals a glutamatergic neuronal subtype susceptible to gene dysregulation via alteration of transcriptional networks.

The genetic architecture of Parkinson's disease (PD) is complex and multiple brain cell subtypes are involved in the neuropathological progression of the disease. Here we aimed to advance our understanding of PD genetic complexity at a cell subtype precision level. Using parallel single-nucleus (sn)RNA-seq and snATAC-seq analyses we simultaneously profiled the transcriptomic and chromatin accessibility landscapes in temporal cortex tissues from 12 PD compared to 12 control subjects at a granular single cell resolution. An integrative bioinformatic pipeline was developed and applied for the analyses of these snMulti-omics datasets. The results identified a subpopulation of cortical glutamatergic excitatory neurons with remarkably altered gene expression in PD, including differentially-expressed genes within PD risk loci identified in genome-wide association studies (GWAS). This was the only neuronal subtype showing significant and robust overexpression of SNCA. Further characterization of this neuronal-subpopulation showed upregulation of specific pathways related to axon guidance, neurite outgrowth and post-synaptic structure, and downregulated pathways involved in presynaptic organization and calcium response. Additionally, we characterized the roles of three molecular mechanisms in governing PD-associated cell subtype-specific dysregulation of gene expression: (1) changes in cis-regulatory element accessibility to transcriptional machinery; (2) changes in the abundance of master transcriptional regulators, including YY1, SP3, and KLF16; (3) candidate regulatory variants in high linkage disequilibrium with PD-GWAS genomic variants impacting transcription factor binding affinities. To our knowledge, this study is the first and the most comprehensive interrogation of the multi-omics landscape of PD at a cell-subtype resolution. Our findings provide new insights into a precise glutamatergic neuronal cell subtype, causal genes, and non-coding regulatory variants underlying the neuropathological progression of PD, paving the way for the development of cell- and gene-targeted therapeutics to halt disease progression as well as genetic biomarkers for early preclinical diagnosis.

A study protocol for risk stratification in children with concussion (RSiCC): Theoretical framework, design, and methods.

Research shows that one in five children will experience a concussion by age 16. Compared to adults, children experience longer and more severe postconcussive symptoms (PCS), with severity and duration varying considerably among children and complicating management of these patients. Persistent PCS can result in increased school absenteeism, social isolation, and psychological distress. Although early PCS diagnosis and access to evidence-based interventions are strongly linked to positive health and academic outcomes, symptom severity and duration are not fully explained by acute post-injury symptoms. Prior research has focused on the role of neuroinflammation in mediating PCS and associated fatigue; however relationship between inflammatory biomarkers and PCS severity, has not examined longitudinally. To identify which children are at high risk for persistent PCS and poor health, academic, and social outcomes, research tracking PCS trajectories and describing school-based impacts across the entire first year postinjury is critically needed. This study will 1) define novel PCS trajectory typologies in a racially/ethnically diverse population of 500 children with concussion (11-17 years, near equal distribution by sex), 2) identify associations between these typologies and patterns of inflammatory biomarkers and genetic variants, 3) develop a risk stratification model to identify children at risk for persistent PCS; and 4) gain unique insights and describe PCS impact, including fatigue, on longer-term academic and social outcomes. We will be the first to use NIH's symptom science model and patient-reported outcomes to explore the patterns of fatigue and other physical, cognitive, psychological, emotional and academic responses to concussion in children over a full year. Our model will enable clinicians and educators to identify children most at risk for poor long-term health, social, and academic outcomes after concussion. This work is critical to meeting our long-term goal of developing personalized concussion symptom-management strategies to improve outcomes and reduce disparities in the health and quality of life of children.

Risk of meningomyelocele mediated by the common 22q11.2 deletion.

Meningomyelocele is one of the most severe forms of neural tube defects (NTDs) and the most frequent structural birth defect of the central nervous system. We assembled the Spina Bifida Sequencing Consortium to identify causes. Exome and genome sequencing of 715 parent-offspring trios identified six patients with chromosomal 22q11.2 deletions, suggesting a 23-fold increased risk compared with the general population. Furthermore, analysis of a separate 22q11.2 deletion cohort suggested a 12- to 15-fold increased NTD risk of meningomyelocele. The loss of Crkl, one of several neural tube-expressed genes within the minimal deletion interval, was sufficient to replicate NTDs in mice, where both penetrance and expressivity were exacerbated by maternal folate deficiency. Thus, the common 22q11.2 deletion confers substantial meningomyelocele risk, which is partially alleviated by folate supplementation.

Validation of human telomere length multi-ancestry meta-analysis association signals identifies POP5 and KBTBD6 as human telomere length regulation genes.

Genome-wide association studies (GWAS) have become well-powered to detect loci associated with telomere length. However, no prior work has validated genes nominated by GWAS to examine their role in telomere length regulation. We conducted a multi-ancestry meta-analysis of 211,369 individuals and identified five novel association signals. Enrichment analyses of chromatin state and cell-type heritability suggested that blood/immune cells are the most relevant cell type to examine telomere length association signals. We validated specific GWAS associations by overexpressing KBTBD6 or POP5 and demonstrated that both lengthened telomeres. CRISPR/Cas9 deletion of the predicted causal regions in K562 blood cells reduced expression of these genes, demonstrating that these loci are related to transcriptional regulation of KBTBD6 and POP5. Our results demonstrate the utility of telomere length GWAS in the identification of telomere length regulation mechanisms and validate KBTBD6 and POP5 as genes affecting telomere length regulation.

Genome-wide association analyses identify 95 risk loci and provide insights into the neurobiology of post-traumatic stress disorder.

Post-traumatic stress disorder (PTSD) genetics are characterized by lower discoverability than most other psychiatric disorders. The contribution to biological understanding from previous genetic studies has thus been limited. We performed a multi-ancestry meta-analysis of genome-wide association studies across 1,222,882 individuals of European ancestry (137,136 cases) and 58,051 admixed individuals with African and Native American ancestry (13,624 cases). We identified 95 genome-wide significant loci (80 new). Convergent multi-omic approaches identified 43 potential causal genes, broadly classified as neurotransmitter and ion channel synaptic modulators (for example, GRIA1, GRM8 and CACNA1E), developmental, axon guidance and transcription factors (for example, FOXP2, EFNA5 and DCC), synaptic structure and function genes (for example, PCLO, NCAM1 and PDE4B) and endocrine or immune regulators (for example, ESR1, TRAF3 and TANK). Additional top genes influence stress, immune, fear and threat-related processes, previously hypothesized to underlie PTSD neurobiology. These findings strengthen our understanding of neurobiological systems relevant to PTSD pathophysiology, while also opening new areas for investigation.

Epigenetic Aging Associations With Psychoneurological Symptoms and Social Functioning in Adults With Sickle Cell Disease.

Objective: Sickle cell disease (SCD), the most common inherited blood disorder in the United States, is associated with severe psychoneurological symptoms. While epigenetic age acceleration has been linked to psychoneurological symptom burden in other diseases, this connection is unexplored in SCD. This study aimed to assess the association between epigenetic age acceleration and psychoneurological symptom burden in SCD. Methods: In this cross-sectional study, emotional impact, pain impact, sleep impact, social functioning, and cognitive function were assessed in 87 adults living with SCD. DNA methylation data were generated from blood specimens and used to calculate epigenetic age using five clocks (Horvath, Hannum, PhenoAge, GrimAge, & DunedinPACE). Associations between epigenetic age acceleration and symptoms were assessed. Results: The sample (N = 87) had a mean (SD) chronologic age was 30.6 (8.1) years. Epigenetic age acceleration was associated with several symptom outcomes. GrimAge age acceleration (ß = -0.49, p = .03) and increased DunedinPACE (ß = -2.23, p = .004) were associated with worse emotional impact scores. PhenoAge (ß = -0.32, p = .04) and the GrimAge (ß = -0.48, p = .05) age acceleration were associated with worse pain impact scores. Increased DunedinPACE (ß = -2.07 p = .04) were associated with worse sleep impact scores. Increased DunedinPACE (ß = -2.87, p = .005) was associated with worse social functioning scores. We did not find associations between epigenetic age acceleration and cognitive function in this sample. Conclusion: Epigenetic age acceleration was associated with worse symptom experiences, suggesting the potential for epigenetic age acceleration as a biomarker to aid in risk stratification or targets for intervention to mitigate symptom burden in SCD.

Deployment-related toxic exposures are associated with worsening mental and physical health after military service: Results from a self-report screening of veterans deployed after 9/11.

Exposure to toxins-such as heavy metals and air pollution-can result in poor health and wellbeing. Recent scientific and media attention has highlighted negative health outcomes associated with toxic exposures for U.S. military personnel deployed overseas. Despite established health risks, less empirical work has examined whether deployment-related toxic exposures are associated with declines in mental and physical health after leaving military service, particularly among the most recent cohort of veterans deployed after September 11, 2001. Using data from 659 U.S. veterans in the VISN 6 MIRECC Post-Deployment Mental Health Study, we tested whether self-reported toxic exposures were associated with poorer mental and physical health. At baseline, veterans who reported more toxic exposures also reported more mental health, ß = 0.14, 95% CI [0.04, 0.23], p = 0.004, and physical health symptoms, ß = 0.21, 95% CI [0.11, 0.30], p < 0.001. Over the next ten years, veterans reporting more toxic exposures also had greater increases in mental health symptoms, ß = 0.23, 95% CI [0.15, 0.31], p < 0.001, physical health symptoms, ß = 0.22, 95% CI [0.14, 0.30], p < 0.001, and chronic disease diagnoses, ß = 0.15, 95% CI [0.07, 0.23], p < 0.001. These associations accounted for demographic and military covariates, including combat exposure. Our findings suggest that toxic exposures are associated with worsening mental and physical health after military service, and this recent cohort of veterans will have increased need for mental health and medical care as they age into midlife and older age.

Demographic characteristics and epigenetic biological aging among post-9/11 veterans: Associations of DunedinPACE with sex, race, and age.

Measures of epigenetic aging derived from DNA methylation (DNAm) have enabled the assessment of biological aging in new populations and cohorts. In the present study, we used an epigenetic measure of aging, DunedinPACE, to examine rates of aging across demographic groups in a sample of 2,309 United States military veterans from the VISN 6 MIRECC's Post-Deployment Mental Health Study. As assessed by DunedinPACE, female veterans were aging faster than male veterans (ß = 0.39, 95 % CI [0.29, 0.48], p < .001), non-Hispanic Black veterans were aging faster than non-Hispanic White veterans (ß = 0.58, 95 % CI [0.50, 0.66], p < .001), and older veterans were biologically aging faster than younger veterans (ß = 0.21, 95 % CI [0.18, 0.25], p < .001). In secondary analyses, these differences in rates of aging were not explained by a variety of biopsychosocial covariates. In addition, the percentage of European genetic admixture in non-Hispanic Black veterans was not associated with DunedinPACE. Our findings suggest that female and non-Hispanic Black veterans are at greater risk of accelerated aging among post-9/11 veterans. Interventions that slow aging might provide relatively greater benefit among veterans comprising these at-risk groups.

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.

Nontargeted Plasma Proteomic Analysis of Renal Disease and Pulmonary Hypertension in Patients with Sickle Cell Disease.

Sickle cell disease (SCD) is characterized by red blood cell sickling, vaso-occlusion, hemolytic anemia, damage to multiple organ systems, and, as a result, shortened life expectancy. Sickle cell disease nephropathy (SCDN) and pulmonary hypertension (pHTN) are common and frequently co-occurring complications of SCD; both are associated with markedly accelerated mortality. To identify candidate circulating biomarkers of SCDN and pHTN, we used mass spectrometry to quantify the relative abundance of >1000 proteins in plasma samples from 189 adults with SCD from the Outcome Modifying Genes in SCD (OMG-SCD) cohort (ProteomeXchange identifier PXD048716). Forty-four proteins were differentially abundant in SCDN, most significantly cystatin-C and collagen α-1(XVIII) chain (COIA1), and 55 proteins were dysregulated in patients with SCDN and pHTN, most significantly insulin-like growth factor-binding protein 6 (IBP6). Network analysis identified a module of 133 coregulated proteins significantly associated with SCDN, that was enriched for extracellular matrix proteins, insulin-like growth factor binding proteins, cell adhesion proteins, EGF-like calcium binding proteins, and several cadherin family members. Collectively, these data provide a comprehensive understanding of plasma protein changes in SCDN and pHTN which validate numerous studies of chronic kidney disease and suggest shared profiles of protein disruption in kidney dysfunction and pHTN among SCD patients.

Posttraumatic stress disorder, trauma, and accelerated biological aging among post-9/11 veterans.

People who experience trauma and develop posttraumatic stress disorder (PTSD) are at increased risk for poor health. One mechanism that could explain this risk is accelerated biological aging, which is associated with the accumulation of chronic diseases, disability, and premature mortality. Using data from 2309 post-9/11 United States military veterans who participated in the VISN 6 MIRECC's Post-Deployment Mental Health Study, we tested whether PTSD and trauma exposure were associated with accelerated rate of biological aging, assessed using a validated DNA methylation (DNAm) measure of epigenetic aging-DunedinPACE. Veterans with current PTSD were aging faster than those who did not have current PTSD, ß = 0.18, 95% CI [0.11, 0.27], p < .001. This effect represented an additional 0.4 months of biological aging each year. Veterans were also aging faster if they reported more PTSD symptoms, ß = 0.13, 95% CI [0.09, 0.16], p < 0.001, or higher levels of trauma exposure, ß = 0.09, 95% CI [0.05, 0.13], p < 0.001. Notably, veterans with past PTSD were aging more slowly than those with current PTSD, ß = -0.21, 95% CI [-0.35, -0.07], p = .003. All reported results accounted for age, gender, self-reported race/ethnicity, and education, and remained when controlling for smoking. Our findings suggest that an accelerated rate of biological aging could help explain how PTSD contributes to poor health and highlights the potential benefits of providing efficacious treatment to populations at increased risk of trauma and PTSD.

Multivariate investigation of aging in mouse models expressing the Alzheimer's protective APOE2 allele: integrating cognitive metrics, brain imaging, and blood transcriptomics.

APOE allelic variation is critical in brain aging and Alzheimer's disease (AD). The APOE2 allele associated with cognitive resilience and neuroprotection against AD remains understudied. We employed a multipronged approach to characterize the transition from middle to old age in mice with APOE2 allele, using behavioral assessments, image-derived morphometry and diffusion metrics, structural connectomics, and blood transcriptomics. We used sparse multiple canonical correlation analyses (SMCCA) for integrative modeling, and graph neural network predictions. Our results revealed brain sub-networks associated with biological traits, cognitive markers, and gene expression. The cingulate cortex emerged as a critical region, demonstrating age-associated atrophy and diffusion changes, with higher fractional anisotropy in males and middle-aged subjects. Somatosensory and olfactory regions were consistently highlighted, indicating age-related atrophy and sex differences. The hippocampus exhibited significant volumetric changes with age, with differences between males and females in CA3 and CA1 regions. SMCCA underscored changes in the cingulate cortex, somatosensory cortex, olfactory regions, and hippocampus in relation to cognition and blood-based gene expression. Our integrative modeling in aging APOE2 carriers revealed a central role for changes in gene pathways involved in localization and the negative regulation of cellular processes. Our results support an important role of the immune system and response to stress. This integrative approach offers novel insights into the complex interplay among brain connectivity, aging, and sex. Our study provides a foundation for understanding the impact of APOE2 allele on brain aging, the potential for detecting associated changes in blood markers, and revealing novel therapeutic intervention targets.

Characterizing epigenetic aging in an adult sickle cell disease cohort.

ABSTRACT: Sickle cell disease (SCD) affects ∼100 000 predominantly African American individuals in the United States, causing significant cellular damage, increased disease complications, and premature death. However, the contribution of epigenetic factors to SCD pathophysiology remains relatively unexplored. DNA methylation (DNAm), a primary epigenetic mechanism for regulating gene expression in response to the environment, is an important driver of normal cellular aging. Several DNAm epigenetic clocks have been developed to serve as a proxy for cellular aging. We calculated the epigenetic ages of 89 adults with SCD (mean age, 30.64 years; 60.64% female) using 5 published epigenetic clocks: Horvath, Hannum, PhenoAge, GrimAge, and DunedinPACE. We hypothesized that in chronic disease, such as SCD, individuals would demonstrate epigenetic age acceleration, but the results differed depending on the clock used. Recently developed clocks more consistently demonstrated acceleration (GrimAge, DunedinPACE). Additional demographic and clinical phenotypes were analyzed to explore their association with epigenetic age estimates. Chronological age was significantly correlated with epigenetic age in all clocks (Horvath, r = 0.88; Hannum, r = 0.89; PhenoAge, r = 0.85; GrimAge, r = 0.88; DunedinPACE, r = 0.34). The SCD genotype was associated with 2 clocks (PhenoAge, P = .02; DunedinPACE, P < .001). Genetic ancestry, biological sex, ß-globin haplotypes, BCL11A rs11886868, and SCD severity were not associated. These findings, among the first to interrogate epigenetic aging in adults with SCD, demonstrate epigenetic age acceleration with recently developed epigenetic clocks but not older-generation clocks. Further development of epigenetic clocks may improve their predictive ability and utility for chronic diseases such as SCD.

Gene-metabolite annotation with shortest reactional distance enhances metabolite genome-wide association studies results.

Metabolite genome-wide association studies (mGWAS) have advanced our understanding of the genetic control of metabolite levels. However, interpreting these associations remains challenging due to a lack of tools to annotate gene-metabolite pairs beyond the use of conservative statistical significance threshold. Here, we introduce the shortest reactional distance (SRD) metric, drawing from the comprehensive KEGG database, to enhance the biological interpretation of mGWAS results. We applied this approach to three independent mGWAS, including a case study on sickle cell disease patients. Our analysis reveals an enrichment of small SRD values in reported mGWAS pairs, with SRD values significantly correlating with mGWAS p values, even beyond the standard conservative thresholds. We demonstrate the utility of SRD annotation in identifying potential false negatives and inaccuracies within current metabolic pathway databases. Our findings highlight the SRD metric as an objective, quantitative and easy-to-compute annotation for gene-metabolite pairs, suitable to integrate statistical evidence to biological networks.

Polarized Desmosome and Hemidesmosome Shedding via Small Extracellular Vesicles is an Early Indicator of Outer Blood-Retina Barrier Dysfunction.

The retinal pigmented epithelium (RPE) constitutes the outer blood-retinal barrier, enables photoreceptor function of the eye, and is constantly exposed to oxidative stress. As such, dysfunction of the RPE underlies pathology leading to development of age-related macular degeneration (AMD), the leading cause of vision loss among the elderly in industrialized nations. A major responsibility of the RPE is to process photoreceptor outer segments, which relies on the proper functioning of its endocytic pathways and endosomal trafficking. Exosomes and other extracellular vesicles (EVs) from RPE are an essential part of these pathways and may be early indicators of cellular stress. To test the role of small EVs (sEVs) including exosomes, that may underlie the early stages of AMD, we used a polarized primary RPE cell culture model under chronic subtoxic oxidative stress. Unbiased proteomic analyses of highly purified basolateral sEVs from oxidatively stressed RPE cultures revealed changes in proteins involved in epithelial barrier integrity. There were also significant changes in proteins accumulating in the basal-side sub-RPE extracellular matrix during oxidative stress, that could be prevented with an inhibitor of sEV release. Thus, chronic subtoxic oxidative stress in primary RPE cultures induces changes in sEV content, including basal-side specific desmosome and hemidesmosome shedding via sEVs. These findings provide novel biomarkers of early cellular dysfunction and opportunity for therapeutic intervention in age-related retinal diseases (e.g., AMD).

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.

Integrative single-nucleus multi-omics analysis prioritizes candidate cis and trans regulatory networks and their target genes in Alzheimer's disease brains.

BACKGROUND: The genetic underpinnings of late-onset Alzheimer's disease (LOAD) are yet to be fully elucidated. Although numerous LOAD-associated loci have been discovered, the causal variants and their target genes remain largely unknown. Since the brain is composed of heterogenous cell subtypes, it is imperative to study the brain on a cell subtype specific level to explore the biological processes underlying LOAD. METHODS: Here, we present the largest parallel single-nucleus (sn) multi-omics study to simultaneously profile gene expression (snRNA-seq) and chromatin accessibility (snATAC-seq) to date, using nuclei from 12 normal and 12 LOAD brains. We identified cell subtype clusters based on gene expression and chromatin accessibility profiles and characterized cell subtype-specific LOAD-associated differentially expressed genes (DEGs), differentially accessible peaks (DAPs) and cis co-accessibility networks (CCANs). RESULTS: Integrative analysis defined disease-relevant CCANs in multiple cell subtypes and discovered LOAD-associated cell subtype-specific candidate cis regulatory elements (cCREs), their candidate target genes, and trans-interacting transcription factors (TFs), some of which, including ELK1, JUN, and SMAD4 in excitatory neurons, were also LOAD-DEGs. Finally, we focused on a subset of cell subtype-specific CCANs that overlap known LOAD-GWAS regions and catalogued putative functional SNPs changing the affinities of TF motifs within LOAD-cCREs linked to LOAD-DEGs, including APOE and MYO1E in a specific subtype of microglia and BIN1 in a subpopulation of oligodendrocytes. CONCLUSIONS: To our knowledge, this study represents the most comprehensive systematic interrogation to date of regulatory networks and the impact of genetic variants on gene dysregulation in LOAD at a cell subtype resolution. Our findings reveal crosstalk between epigenetic, genomic, and transcriptomic determinants of LOAD pathogenesis and define catalogues of candidate genes, cCREs, and variants involved in LOAD genetic etiology and the cell subtypes in which they act to exert their pathogenic effects. Overall, these results suggest that cell subtype-specific cis-trans interactions between regulatory elements and TFs, and the genes dysregulated by these networks contribute to the development of LOAD.

Genomic Structural Equation Modeling Reveals Latent Phenotypes in the Human Cortex with Distinct Genetic Architecture.

Genetic contributions to human cortical structure manifest pervasive pleiotropy. This pleiotropy may be harnessed to identify unique genetically-informed parcellations of the cortex that are neurobiologically distinct from functional, cytoarchitectural, or other cortical parcellation schemes. We investigated genetic pleiotropy by applying genomic structural equation modeling (SEM) to map the genetic architecture of cortical surface area (SA) and cortical thickness (CT) for the 34 brain regions recently reported in the ENIGMA cortical GWAS. Genomic SEM uses the empirical genetic covariance estimated from GWAS summary statistics with LD score regression (LDSC) to discover factors underlying genetic covariance, which we are denoting genetically informed brain networks (GIBNs). Genomic SEM can fit a multivariate GWAS from summary statistics for each of the GIBNs, which can subsequently be used for LD score regression (LDSC). We found the best-fitting model of cortical SA identified 6 GIBNs and CT identified 4 GIBNs. The multivariate GWASs of these GIBNs identified 74 genome-wide significant (GWS) loci (p<5×10-8), including many previously implicated in neuroimaging phenotypes, behavioral traits, and psychiatric conditions. LDSC of GIBN GWASs found that SA-derived GIBNs had a positive genetic correlation with bipolar disorder (BPD), and cannabis use disorder, indicating genetic predisposition to a larger SA in the specific GIBN is associated with greater genetic risk of these disorders. A negative genetic correlation was observed with attention deficit hyperactivity disorder (ADHD), major depressive disorder (MDD), and insomnia, indicating genetic predisposition to a larger SA in the specific GIBN is associated with lower genetic risk of these disorders. CT GIBNs displayed a negative genetic correlation with alcohol dependence. Jointly modeling the genetic architecture of complex traits and investigating multivariate genetic links across phenotypes offers a new vantage point for mapping the cortex into genetically informed networks.

Discovery of 95 PTSD loci provides insight into genetic architecture and neurobiology of trauma and stress-related disorders.

Posttraumatic stress disorder (PTSD) genetics are characterized by lower discoverability than most other psychiatric disorders. The contribution to biological understanding from previous genetic studies has thus been limited. We performed a multi-ancestry meta-analysis of genome-wide association studies across 1,222,882 individuals of European ancestry (137,136 cases) and 58,051 admixed individuals with African and Native American ancestry (13,624 cases). We identified 95 genome-wide significant loci (80 novel). Convergent multi-omic approaches identified 43 potential causal genes, broadly classified as neurotransmitter and ion channel synaptic modulators (e.g., GRIA1, GRM8, CACNA1E ), developmental, axon guidance, and transcription factors (e.g., FOXP2, EFNA5, DCC ), synaptic structure and function genes (e.g., PCLO, NCAM1, PDE4B ), and endocrine or immune regulators (e.g., ESR1, TRAF3, TANK ). Additional top genes influence stress, immune, fear, and threat-related processes, previously hypothesized to underlie PTSD neurobiology. These findings strengthen our understanding of neurobiological systems relevant to PTSD pathophysiology, while also opening new areas for investigation.