Learning epistatic polygenic phenotypes with Boolean interactions.

Detecting epistatic drivers of human phenotypes is a considerable challenge. Traditional approaches use regression to sequentially test multiplicative interaction terms involving pairs of genetic variants. For higher-order interactions and genome-wide large-scale data, this strategy is computationally intractable. Moreover, multiplicative terms used in regression modeling may not capture the form of biological interactions. Building on the Predictability, Computability, Stability (PCS) framework, we introduce the epiTree pipeline to extract higher-order interactions from genomic data using tree-based models. The epiTree pipeline first selects a set of variants derived from tissue-specific estimates of gene expression. Next, it uses iterative random forests (iRF) to search training data for candidate Boolean interactions (pairwise and higher-order). We derive significance tests for interactions, based on a stabilized likelihood ratio test, by simulating Boolean tree-structured null (no epistasis) and alternative (epistasis) distributions on hold-out test data. Finally, our pipeline computes PCS epistasis p-values that probabilisticly quantify improvement in prediction accuracy via bootstrap sampling on the test set. We validate the epiTree pipeline in two case studies using data from the UK Biobank: predicting red hair and multiple sclerosis (MS). In the case of predicting red hair, epiTree recovers known epistatic interactions surrounding MC1R and novel interactions, representing non-linearities not captured by logistic regression models. In the case of predicting MS, a more complex phenotype than red hair, epiTree rankings prioritize novel interactions surrounding HLA-DRB1, a variant previously associated with MS in several populations. Taken together, these results highlight the potential for epiTree rankings to help reduce the design space for follow up experiments.

Assessing the potential of polygenic scores to strengthen medical risk prediction models of COVID-19.

As findings on the epidemiological and genetic risk factors for coronavirus disease-19 (COVID-19) continue to accrue, their joint power and significance for prospective clinical applications remains virtually unexplored. Severity of symptoms in individuals affected by COVID-19 spans a broad spectrum, reflective of heterogeneous host susceptibilities across the population. Here, we assessed the utility of epidemiological risk factors to predict disease severity prospectively, and interrogated genetic information (polygenic scores) to evaluate whether they can provide further insights into symptom heterogeneity. A standard model was trained to predict severe COVID-19 based on principal component analysis and logistic regression based on information from eight known medical risk factors for COVID-19 measured before 2018. In UK Biobank participants of European ancestry, the model achieved a relatively high performance (area under the receiver operating characteristic curve ~90%). Polygenic scores for COVID-19 computed from summary statistics of the Covid19 Host Genetics Initiative displayed significant associations with COVID-19 in the UK Biobank (p-values as low as 3.96e-9, all with R2 under 1%), but were unable to robustly improve predictive performance of the non-genetic factors. However, error analysis of the non-genetic models suggested that affected individuals misclassified by the medical risk factors (predicted low risk but actual high risk) display a small but consistent increase in polygenic scores. Overall, the results indicate that simple models based on health-related epidemiological factors measured years before COVID-19 onset can achieve high predictive power. Associations between COVID-19 and genetic factors were statistically robust, but currently they have limited predictive power for translational settings. Despite that, the outcomes also suggest that severely affected cases with a medical history profile of low risk might be partly explained by polygenic factors, prompting development of boosted COVID-19 polygenic models based on new data and tools to aid risk-prediction.

Oligogenic Architecture of Rare Noncoding Variants Distinguishes 4 Congenital Heart Disease Phenotypes.

BACKGROUND: Congenital heart disease (CHD) is highly heritable, but the power to identify inherited risk has been limited to analyses of common variants in small cohorts. METHODS: We performed reimputation of 4 CHD cohorts (n=55 342) to the TOPMed reference panel (freeze 5), permitting meta-analysis of 14 784 017 variants including 6 035 962 rare variants of high imputation quality as validated by whole genome sequencing. RESULTS: Meta-analysis identified 16 novel loci, including 12 rare variants, which displayed moderate or large effect sizes (median odds ratio, 3.02) for 4 separate CHD categories. Analyses of chromatin structure link 13 of the genome-wide significant loci to key genes in cardiac development; rs373447426 (minor allele frequency, 0.003 [odds ratio, 3.37 for Conotruncal heart disease]; P=1.49×10-8) is predicted to disrupt chromatin structure for 2 nearby genes BDH1 and DLG1 involved in Conotruncal development. A lead variant rs189203952 (minor allele frequency, 0.01 [odds ratio, 2.4 for left ventricular outflow tract obstruction]; P=1.46×10-8) is predicted to disrupt the binding sites of 4 transcription factors known to participate in cardiac development in the promoter of SPAG9. A tissue-specific model of chromatin conformation suggests that common variant rs78256848 (minor allele frequency, 0.11 [odds ratio, 1.4 for Conotruncal heart disease]; P=2.6×10-8) physically interacts with NCAM1 (PFDR=1.86×10-27), a neural adhesion molecule acting in cardiac development. Importantly, while each individual malformation displayed substantial heritability (observed h2 ranging from 0.26 for complex malformations to 0.37 for left ventricular outflow tract obstructive disease) the risk for different CHD malformations appeared to be separate, without genetic correlation measured by linkage disequilibrium score regression or regional colocalization. CONCLUSIONS: We describe a set of rare noncoding variants conferring significant risk for individual heart malformations which are linked to genes governing cardiac development. These results illustrate that the oligogenic basis of CHD and significant heritability may be linked to rare variants outside protein-coding regions conferring substantial risk for individual categories of cardiac malformation.

Genetically-informed prediction of short-term Parkinson's disease progression.

Parkinson's disease (PD) treatments modify disease symptoms but have not been shown to slow progression, characterized by gradual and varied motor and non-motor changes overtime. Variation in PD progression hampers clinical research, resulting in long and expensive clinical trials prone to failure. Development of models for short-term PD progression prediction could be useful for shortening the time required to detect disease-modifying drug effects in clinical studies. PD progressors were defined by an increase in MDS-UPDRS scores at 12-, 24-, and 36-months post-baseline. Using only baseline features, PD progression was separately predicted across all timepoints and MDS-UPDRS subparts in independent, optimized, XGBoost models. These predictions plus baseline features were combined into a meta-predictor for 12-month MDS UPDRS Total progression. Data from the Parkinson's Progression Markers Initiative (PPMI) were used for training with independent testing on the Parkinson's Disease Biomarkers Program (PDBP) cohort. 12-month PD total progression was predicted with an F-measure 0.77, ROC AUC of 0.77, and PR AUC of 0.76 when tested on a hold-out PPMI set. When tested on PDBP we achieve a F-measure 0.75, ROC AUC of 0.74, and PR AUC of 0.73. Exclusion of genetic predictors led to the greatest loss in predictive accuracy; ROC AUC of 0.66, PR AUC of 0.66-0.68 for both PPMI and PDBP testing. Short-term PD progression can be predicted with a combination of survey-based, neuroimaging, physician examination, and genetic predictors. Dissection of the interplay between genetic risk, motor symptoms, non-motor symptoms, and longer-term expected rates of progression enable generalizable predictions.

Pan-ancestry exome-wide association analyses of COVID-19 outcomes in 586,157 individuals.

Severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) causes coronavirus disease 2019 (COVID-19), a respiratory illness that can result in hospitalization or death. We used exome sequence data to investigate associations between rare genetic variants and seven COVID-19 outcomes in 586,157 individuals, including 20,952 with COVID-19. After accounting for multiple testing, we did not identify any clear associations with rare variants either exome wide or when specifically focusing on (1) 13 interferon pathway genes in which rare deleterious variants have been reported in individuals with severe COVID-19, (2) 281 genes located in susceptibility loci identified by the COVID-19 Host Genetics Initiative, or (3) 32 additional genes of immunologic relevance and/or therapeutic potential. Our analyses indicate there are no significant associations with rare protein-coding variants with detectable effect sizes at our current sample sizes. Analyses will be updated as additional data become available, and results are publicly available through the Regeneron Genetics Center COVID-19 Results Browser.

Genetic control of variability in subcortical and intracranial volumes.

Sensitivity to external demands is essential for adaptation to dynamic environments, but comes at the cost of increased risk of adverse outcomes when facing poor environmental conditions. Here, we apply a novel methodology to perform genome-wide association analysis of mean and variance in ten key brain features (accumbens, amygdala, caudate, hippocampus, pallidum, putamen, thalamus, intracranial volume, cortical surface area, and cortical thickness), integrating genetic and neuroanatomical data from a large lifespan sample (n = 25,575 individuals; 8-89 years, mean age 51.9 years). We identify genetic loci associated with phenotypic variability in thalamus volume and cortical thickness. The variance-controlling loci involved genes with a documented role in brain and mental health and were not associated with the mean anatomical volumes. This proof-of-principle of the hypothesis of a genetic regulation of brain volume variability contributes to establishing the genetic basis of phenotypic variance (i.e., heritability), allows identifying different degrees of brain robustness across individuals, and opens new research avenues in the search for mechanisms controlling brain and mental health.

Cardiac Imaging of Aortic Valve Area From 34 287 UK Biobank Participants Reveals Novel Genetic Associations and Shared Genetic Comorbidity With Multiple Disease Phenotypes.

BACKGROUND: The aortic valve is an important determinant of cardiovascular physiology and anatomic location of common human diseases. METHODS: From a sample of 34 287 white British ancestry participants, we estimated functional aortic valve area by planimetry from prospectively obtained cardiac magnetic resonance imaging sequences of the aortic valve. Aortic valve area measurements were submitted to genome-wide association testing, followed by polygenic risk scoring and phenome-wide screening, to identify genetic comorbidities. RESULTS: A genome-wide association study of aortic valve area in these UK Biobank participants showed 3 significant associations, indexed by rs71190365 (chr13:50764607, DLEU1, P=1.8×10-9), rs35991305 (chr12:94191968, CRADD, P=3.4×10-8), and chr17:45013271:C:T (GOSR2, P=5.6×10-8). Replication on an independent set of 8145 unrelated European ancestry participants showed consistent effect sizes in all 3 loci, although rs35991305 did not meet nominal significance. We constructed a polygenic risk score for aortic valve area, which in a separate cohort of 311 728 individuals without imaging demonstrated that smaller aortic valve area is predictive of increased risk for aortic valve disease (odds ratio, 1.14; P=2.3×10-6). After excluding subjects with a medical diagnosis of aortic valve stenosis (remaining n=308 683 individuals), phenome-wide association of >10 000 traits showed multiple links between the polygenic score for aortic valve disease and key health-related comorbidities involving the cardiovascular system and autoimmune disease. Genetic correlation analysis supports a shared genetic etiology with between aortic valve area and birth weight along with other cardiovascular conditions. CONCLUSIONS: These results illustrate the use of automated phenotyping of cardiac imaging data from the general population to investigate the genetic etiology of aortic valve disease, perform clinical prediction, and uncover new clinical and genetic correlates of cardiac anatomy.

Twin study designs as a tool to identify new candidate genes for depression: A systematic review of DNA methylation studies.

Monozygotic (MZ) twin studies constitute a key resource for the dissection of environmental and biological risk factors for human complex disorders. Given that epigenetic differences accumulate throughout the lifespan, the assessment of MZ twin pairs discordant for depression offers a genetically informative design to explore DNA methylation while accounting for the typical confounders of the field, shared by co-twins of a pair. In this review, we systematically evaluate all twin studies published to date assessing DNA methylation in association with depressive phenotypes. However, difficulty to recruit large numbers of MZ twin pairs fails to provide enough sample size to develop genome-wide approaches. Alternatively, region and pathway analysis revealed an enrichment for nervous system related functions; likewise, evidence supports an accumulation of methylation variability in affected subjects when compared to their co-twins. Nevertheless, longitudinal studies incorporating known risk factors for depression such as childhood trauma are required for understanding the role that DNA methylation plays in the etiology of depression.

Brain scans from 21,297 individuals reveal the genetic architecture of hippocampal subfield volumes.

The hippocampus is a heterogeneous structure, comprising histologically distinguishable subfields. These subfields are differentially involved in memory consolidation, spatial navigation and pattern separation, complex functions often impaired in individuals with brain disorders characterized by reduced hippocampal volume, including Alzheimer's disease (AD) and schizophrenia. Given the structural and functional heterogeneity of the hippocampal formation, we sought to characterize the subfields' genetic architecture. T1-weighted brain scans (n = 21,297, 16 cohorts) were processed with the hippocampal subfields algorithm in FreeSurfer v6.0. We ran a genome-wide association analysis on each subfield, co-varying for whole hippocampal volume. We further calculated the single-nucleotide polymorphism (SNP)-based heritability of 12 subfields, as well as their genetic correlation with each other, with other structural brain features and with AD and schizophrenia. All outcome measures were corrected for age, sex and intracranial volume. We found 15 unique genome-wide significant loci across six subfields, of which eight had not been previously linked to the hippocampus. Top SNPs were mapped to genes associated with neuronal differentiation, locomotor behaviour, schizophrenia and AD. The volumes of all the subfields were estimated to be heritable (h2 from 0.14 to 0.27, all p < 1 × 10-16) and clustered together based on their genetic correlations compared with other structural brain features. There was also evidence of genetic overlap of subicular subfield volumes with schizophrenia. We conclude that hippocampal subfields have partly distinct genetic determinants associated with specific biological processes and traits. Taking into account this specificity may increase our understanding of hippocampal neurobiology and associated pathologies.

Publisher Correction: Common brain disorders are associated with heritable patterns of apparent aging of the brain.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Association between the 4p16 genomic locus and different types of congenital heart disease: results from adult survivors in the UK Biobank.

Congenital heart disease is the most common birth defect in newborns and the leading cause of death in infancy, affecting nearly 1% of live births. A locus in chromosome 4p16, adjacent to MSX1 and STX18, has been associated with atrial septal defects (ASD) in multiple European and Chinese cohorts. Here, genotyping data from the UK Biobank was used to test for associations between this locus and congenital heart disease in adult survivors of left ventricular outflow tract obstruction (n = 164) and ASD (n = 223), with a control sample of 332,788 individuals, and a meta-analysis of the new and existing ASD data was performed. The results show an association between the previously reported markers at 4p16 and risk for either ASD or left ventricular outflow tract obstruction, with effect sizes similar to the published data (OR between 1.27-1.45; all p < 0.05). Differences in allele frequencies remained constant through the studied age range (40-70 years), indicating that the variants themselves do not drive lethal genetic defects. Meta-analysis shows an OR of 1.35 (95% CI: 1.25-1.46; p < 10-4) for the association with ASD. The findings show that the genetic associations with ASD can be generalized to adult survivors of both ASD and left ventricular lesions. Although the 4p16 associations are statistically compelling, the mentioned alleles confer only a small risk for disease and their frequencies in this adult sample are the same as in children, likely limiting their clinical significance. Further epidemiological and functional studies may elicit factors triggering disease in interaction with the risk alleles.

Common brain disorders are associated with heritable patterns of apparent aging of the brain.

Common risk factors for psychiatric and other brain disorders are likely to converge on biological pathways influencing the development and maintenance of brain structure and function across life. Using structural MRI data from 45,615 individuals aged 3-96 years, we demonstrate distinct patterns of apparent brain aging in several brain disorders and reveal genetic pleiotropy between apparent brain aging in healthy individuals and common brain disorders.

Brain Heterogeneity in Schizophrenia and Its Association With Polygenic Risk.

Importance: Between-individual variability in brain structure is determined by gene-environment interactions, possibly reflecting differential sensitivity to environmental and genetic perturbations. Magnetic resonance imaging (MRI) studies have revealed thinner cortices and smaller subcortical volumes in patients with schizophrenia. However, group-level comparisons may mask considerable within-group heterogeneity, which has largely remained unnoticed in the literature. Objectives: To compare brain structural variability between individuals with schizophrenia and healthy controls and to test whether respective variability reflects the polygenic risk score (PRS) for schizophrenia in an independent sample of healthy controls. Design, Setting, and Participants: This case-control and polygenic risk analysis compared MRI-derived cortical thickness and subcortical volumes between healthy controls and patients with schizophrenia across 16 cohorts and tested for associations between PRS and MRI features in a control cohort from the UK Biobank. Data were collected from October 27, 2004, through April 12, 2018, and analyzed from December 3, 2017, through August 1, 2018. Main Outcomes and Measures: Mean and dispersion parameters were estimated using double generalized linear models. Vertex-wise analysis was used to assess cortical thickness, and regions-of-interest analyses were used to assess total cortical volume, total surface area, and white matter, subcortical, and hippocampal subfield volumes. Follow-up analyses included within-sample analysis, test of robustness of the PRS threshold, population covariates, outlier removal, and control for image quality. Results: A comparison of 1151 patients with schizophrenia (mean [SD] age, 33.8 [10.6] years; 68.6% male [n = 790] and 31.4% female [n = 361]) with 2010 healthy controls (mean [SD] age, 32.6 [10.4] years; 56.0% male [n = 1126] and 44.0% female [n = 884]) revealed higher heterogeneity in schizophrenia for cortical thickness and area (t = 3.34), cortical (t = 3.24) and ventricle (t range, 3.15-5.78) volumes, and hippocampal subfields (t range, 2.32-3.55). In the UK Biobank sample of 12 490 participants (mean [SD] age, 55.9 [7.5] years; 48.2% male [n = 6025] and 51.8% female [n = 6465]), higher PRS was associated with thinner frontal and temporal cortices and smaller left CA2/3 (t = -3.00) but was not significantly associated with dispersion. Conclusions and Relevance: This study suggests that schizophrenia is associated with substantial brain structural heterogeneity beyond the mean differences. These findings may reflect higher sensitivity to environmental and genetic perturbations in patients, supporting the heterogeneous nature of schizophrenia. A higher PRS was associated with thinner frontotemporal cortices and smaller hippocampal subfield volume, but not heterogeneity. This finding suggests that brain variability in schizophrenia results from interactions between environmental and genetic factors that are not captured by the PRS. Factors contributing to heterogeneity in frontotemporal cortices and hippocampus are key to furthering our understanding of how genetic and environmental factors shape brain biology in schizophrenia.

Cerebellar Gray Matter Volume Is Associated With Cognitive Function and Psychopathology in Adolescence.

BACKGROUND: Accumulating evidence supports cerebellar involvement in mental disorders, such as schizophrenia, bipolar disorder, depression, anxiety disorders, and attention-deficit/hyperactivity disorder. However, little is known about the cerebellum in developmental stages of these disorders. In particular, whether cerebellar morphology is associated with early expression of specific symptom domains remains unclear. METHODS: We used machine learning to test whether cerebellar morphometric features could robustly predict general cognitive function and psychiatric symptoms in a large and well-characterized developmental community sample centered on adolescence (Philadelphia Neurodevelopmental Cohort, n = 1401, age 8-23 years). RESULTS: Cerebellar morphology was associated with both general cognitive function and general psychopathology (mean correlations between predicted and observed values: r = .20 and r = .13; p < .001). Analyses of specific symptom domains revealed significant associations with rates of norm-violating behavior (r = .17; p < .001) as well as psychosis (r = .12; p < .001) and anxiety (r = .09; p = .012) symptoms. In contrast, we observed no associations with attention deficits or depressive, manic, or obsessive-compulsive symptoms. Crucially, across 52 brain-wide anatomical features, cerebellar features emerged as the most important for prediction of general psychopathology, psychotic symptoms, and norm-violating behavior. Moreover, the association between cerebellar volume and psychotic symptoms and, to a lesser extent, norm-violating behavior remained significant when adjusting for several potentially confounding factors. CONCLUSIONS: The robust associations with psychiatric symptoms in the age range when these typically emerge highlight the cerebellum as a key brain structure in the development of severe mental disorders.

Oxytocin pathway gene networks in the human brain.

Oxytocin is a neuropeptide involved in animal and human reproductive and social behavior. Three oxytocin signaling genes have been frequently implicated in human social behavior: OXT (structural gene for oxytocin), OXTR (oxytocin receptor), and CD38 (oxytocin secretion). Here, we characterized the distribution of OXT, OXTR, and CD38 mRNA across the human brain by creating voxel-by-voxel volumetric expression maps, and identified putative gene pathway interactions by comparing gene expression patterns across 20,737 genes. Expression of the three selected oxytocin pathway genes was enriched in subcortical and olfactory regions and there was high co-expression with several dopaminergic and muscarinic acetylcholine genes, reflecting an anatomical basis for critical gene pathway interactions. fMRI meta-analysis revealed that the oxytocin pathway gene maps correspond with the processing of anticipatory, appetitive, and aversive cognitive states. The oxytocin signaling system may interact with dopaminergic and muscarinic acetylcholine signaling to modulate cognitive state processes involved in complex human behaviors.

Data-Driven Clustering Reveals a Link Between Symptoms and Functional Brain Connectivity in Depression.

BACKGROUND: Depression is a complex disorder with large interindividual variability in symptom profiles that often occur alongside symptoms of other psychiatric domains, such as anxiety. A dimensional and symptom-based approach may help refine the characterization of depressive and anxiety disorders and thus aid in establishing robust biomarkers. We use resting-state functional magnetic resonance imaging to assess the brain functional connectivity correlates of a symptom-based clustering of individuals. METHODS: We assessed symptoms using the Beck Depression and Beck Anxiety Inventories in individuals with or without a history of depression (N = 1084) and high-dimensional data clustering to form subgroups based on symptom profiles. We compared dynamic and static functional connectivity between subgroups in a subset of the total sample (n = 252). RESULTS: We identified five subgroups with distinct symptom profiles, which cut across diagnostic boundaries with different total severity, symptom patterns, and centrality. For instance, inability to relax, fear of the worst, and feelings of guilt were among the most severe symptoms in subgroups 1, 2, and 3, respectively. The distribution of individuals was 32%, 25%, 22%, 10%, and 11% in subgroups 1 to 5, respectively. These subgroups showed evidence of differential static brain-connectivity patterns, in particular comprising a frontotemporal network. In contrast, we found no significant associations with clinical sum scores, dynamic functional connectivity, or global connectivity. CONCLUSIONS: Adding to the pursuit of individual-based treatment, subtyping based on a dimensional conceptualization and unique constellations of anxiety and depression symptoms is supported by distinct patterns of static functional connectivity in the brain.

Loss of function, missense, and intronic variants in NOTCH1 confer different risks for left ventricular outflow tract obstructive heart defects in two European cohorts.

Loss of function variants in NOTCH1 cause left ventricular outflow tract obstructive defects (LVOTO). However, the risk conferred by rare and noncoding variants in NOTCH1 for LVOTO remains largely uncharacterized. In a cohort of 49 families affected by hypoplastic left heart syndrome, a severe form of LVOTO, we discovered predicted loss of function NOTCH1 variants in 6% of individuals. Rare or low-frequency missense variants were found in 16% of families. To make a quantitative estimate of the genetic risk posed by variants in NOTCH1 for LVOTO, we studied associations of 400 coding and noncoding variants in NOTCH1 in 1,085 cases and 332,788 controls from the UK Biobank. Two rare intronic variants in strong linkage disequilibrium displayed significant association with risk for LVOTO amongst European-ancestry individuals. This result was replicated in an independent analysis of 210 cases and 68,762 controls of non-European and mixed ancestry. In conclusion, carrying rare predicted loss of function variants in NOTCH1 confer significant risk for LVOTO. In addition, the two intronic variants seem to be associated with an increased risk for these defects. Our approach demonstrates the utility of population-based data sets in quantifying the specific risk of individual variants for disease-related phenotypes.

Epigenetic outlier profiles in depression: A genome-wide DNA methylation analysis of monozygotic twins.

Recent discoveries highlight the importance of stochastic epigenetic changes, as indexed by epigenetic outlier DNA methylation signatures, as a valuable tool to understand aberrant cell function and subsequent human pathology. There is evidence of such changes in different complex disorders as diverse as cancer, obesity and, to a lesser extent, depression. The current study was aimed at identifying outlying DNA methylation signatures of depressive psychopathology. Here, genome-wide DNA methylation levels were measured (by means of Illumina Infinium HumanMethylation450 Beadchip) in peripheral blood of thirty-four monozygotic twins informative for depressive psychopathology (lifetime DSM-IV diagnoses). This dataset was explored to identify outlying epigenetic signatures of depression, operationalized as extreme hyper- or hypo-methylation in affected co-twins from discordant pairs that is not observed across the rest of the study sample. After adjusting for blood cell count, there were thirteen CpG sites across which depressed co-twins from the discordant pairs exhibited outlying DNA methylation signatures. None of them exhibited a methylation outlier profile in the concordant and healthy pairs, and some of these loci spanned genes previously associated with neuropsychiatric phenotypes, such as GHSR and KCNQ1. This exploratory study provides preliminary proof-of-concept validation that epigenetic outlier profiles derived from genome-wide DNA methylation data may be related to depression risk.

White matter aberrations and age-related trajectories in patients with schizophrenia and bipolar disorder revealed by diffusion tensor imaging.

Supported by histological and genetic evidence implicating myelin, neuroinflammation and oligodendrocyte dysfunction in schizophrenia spectrum disorders (SZ), diffusion tensor imaging (DTI) studies have consistently shown white matter (WM) abnormalities when compared to healthy controls (HC). The diagnostic specificity remains unclear, with bipolar disorders (BD) frequently conceptualized as a less severe clinical manifestation along a psychotic spectrum. Further, the age-related dynamics and possible sex differences of WM abnormalities in SZ and BD are currently understudied. Using tract-based spatial statistics (TBSS) we compared DTI-based microstructural indices between SZ (n = 128), BD (n = 61), and HC (n = 293). We tested for age-by-group and sex-by-group interactions, computed effect sizes within different age-bins and within genders. TBSS revealed global reductions in fractional anisotropy (FA) and increases in radial (RD) diffusivity in SZ compared to HC, with strongest effects in the body and splenium of the corpus callosum, and lower FA in SZ compared to BD in right inferior longitudinal fasciculus and right inferior fronto-occipital fasciculus, and no significant differences between BD and HC. The results were not strongly dependent on age or sex. Despite lack of significant group-by-age interactions, a sliding-window approach supported widespread WM involvement in SZ with most profound differences in FA from the late 20 s.

Increased methylation at an unexplored glucocorticoid responsive element within exon 1D of NR3C1 gene is related to anxious-depressive disorders and decreased hippocampal connectivity.

Among the major psychiatric disorders, anxious-depressive disorders stand out as one of the more prevalent and more frequently associated with hypothalamic-pituitary-adrenal (HPA) axis abnormalities. Methylation at the exon 1F of the glucocorticoid receptor gene NR3C1 has been associated with both early stress exposure and risk for developing a psychiatric disorder; however, other NR3C1 promoter regions have been underexplored. Exon 1D emerges as a suggestive new target in stress-related disorders epigenetically sensitive to early adversity. After assessment of 48 monozygotic twin pairs (n=96 subjects) informative for lifetime history of anxious-depressive disorders, they were classified as concordant, discordant or healthy in function of whether both, one or neither twin in each pair had a lifetime diagnosis of anxious-depressive disorders. DNA for epigenetic analysis was extracted from peripheral blood. Exon 1F and exon 1D CpG-specific methylation was analysed by means of pyrosequencing technology. Functional magnetic resonance imaging was available for 54 subjects (n=27 twin pairs). Exon 1D CpG-specific methylation within a glucocorticoid responsive element (GRE) was correlated with familial burden of anxious-depressive disorders (r=0.35, z=2.26, p=0.02). Right hippocampal connectivity was significantly associated with CpG-specific GRE methylation (ß=-2.33, t=-2.85, p=0.01). Exon 1F was uniformly hypomethylated across all subgroups of the present sample. GRE hypermethylation at exon 1D of the NR3C1 gene in monozygotic twins concordant for anxious-depressive disorders suggests this region plays a role in increasing vulnerability to psychosocial stress, partly mediated by altered hippocampal connectivity.