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.

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.

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.

Distinguishing early and late brain aging from the Alzheimer's disease spectrum: consistent morphological patterns across independent samples.

Alzheimer's disease (AD) is a debilitating age-related neurodegenerative disorder. Accurate identification of individuals at risk is complicated as AD shares cognitive and brain features with aging. We applied linked independent component analysis (LICA) on three complementary measures of gray matter structure: cortical thickness, area and gray matter density of 137 AD, 78 mild (MCI) and 38 subjective cognitive impairment patients, and 355 healthy adults aged 18-78 years to identify dissociable multivariate morphological patterns sensitive to age and diagnosis. Using the lasso classifier, we performed group classification and prediction of cognition and age at different age ranges to assess the sensitivity and diagnostic accuracy of the LICA patterns in relation to AD, as well as early and late healthy aging. Three components showed high sensitivity to the diagnosis and cognitive status of AD, with different relationships with age: one reflected an anterior-posterior gradient in thickness and gray matter density and was uniquely related to diagnosis, whereas the other two, reflecting widespread cortical thickness and medial temporal lobe volume, respectively, also correlated significantly with age. Repeating the LICA decomposition and between-subject analysis on ADNI data, including 186 AD, 395 MCI and 220 age-matched healthy controls, revealed largely consistent brain patterns and clinical associations across samples. Classification results showed that multivariate LICA-derived brain characteristics could be used to predict AD and age with high accuracy (area under ROC curve up to 0.93 for classification of AD from controls). Comparison between classifiers based on feature ranking and feature selection suggests both common and unique feature sets implicated in AD and aging, and provides evidence of distinct age-related differences in early compared to late aging.

Altered amygdalar resting-state connectivity in depression is explained by both genes and environment.

Recent findings indicate that alterations of the amygdalar resting-state fMRI connectivity play an important role in the etiology of depression. While both depression and resting-state brain activity are shaped by genes and environment, the relative contribution of genetic and environmental factors mediating the relationship between amygdalar resting-state connectivity and depression remain largely unexplored. Likewise, novel neuroimaging research indicates that different mathematical representations of resting-state fMRI activity patterns are able to embed distinct information relevant to brain health and disease. The present study analyzed the influence of genes and environment on amygdalar resting-state fMRI connectivity, in relation to depression risk. High-resolution resting-state fMRI scans were analyzed to estimate functional connectivity patterns in a sample of 48 twins (24 monozygotic pairs) informative for depressive psychopathology (6 concordant, 8 discordant and 10 healthy control pairs). A graph-theoretical framework was employed to construct brain networks using two methods: (i) the conventional approach of filtered BOLD fMRI time-series and (ii) analytic components of this fMRI activity. Results using both methods indicate that depression risk is increased by environmental factors altering amygdalar connectivity. When analyzing the analytic components of the BOLD fMRI time-series, genetic factors altering the amygdala neural activity at rest show an important contribution to depression risk. Overall, these findings show that both genes and environment modify different patterns the amygdala resting-state connectivity to increase depression risk. The genetic relationship between amygdalar connectivity and depression may be better elicited by examining analytic components of the brain resting-state BOLD fMRI signals.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.