The aim of logistic regression is to estimate genetic effects on disease risk, while survival analysis aims to determine effects on age of onset. In practice, genetic variants may affect both types of outcomes. A cure survival model analyzes logistic and survival effects simultaneously. The aim of this simulation study is to assess the performance of logistic regression and traditional survival analysis under a cure model and to investigate the benefits of cure survival analysis. We simulated data under a cure model and varied the percentage of subjects at risk for disease (cure fraction), the logistic and survival effect sizes, and the contribution of genetic background risk factors. We then computed the error rates and estimation bias of logistic, Cox proportional hazards (PH), and cure PH analysis, respectively. The power of logistic and Cox PH analysis is sensitive to the cure fraction and background heritability. Our results show that traditional Cox PH analysis may erroneously detect age of onset effects if no such effects are present in the data. In the presence of genetic background risk even the cure model results in biased estimates of both the odds ratio and the hazard ratio. Cure survival analysis takes cure fractions into account and can be used to simultaneously estimate the effect of genetic variants on disease risk and age of onset. Since genome-wide cure survival analysis is not computationally feasible, we recommend this analysis for genetic variants that are significant in a traditional survival analysis.

Bipolar disorder (BD) patients show aberrant white matter microstructure compared to healthy controls but little is known about the relation with clinical characteristics. We therefore investigated the relation of white matter microstructure with the main pharmacological treatments as well its relation with IQ. Patients with BD (N = 257) and controls (N = 167) underwent diffusion tensor imaging (DTI) and comprehensive clinically assessments including IQ estimates. DTI images were analyzed using tract-based spatial statistics. Fractional anisotropy (FA) and Mean Diffusivity (MD) were determined. Patients had significantly lower FA and higher MD values throughout the white matter skeleton compared to controls. Within the BD patients, lithium use was associated with higher FA and lower MD. Antipsychotic medication use in the BD patients was not associated with FA but, in contrast to lithium, was associated with higher MD. IQ was significantly positively correlated with FA and negatively with MD in patients as well as in controls. In this large DTI study we found evidence for marked differences in FA and MD particularly in (but not restricted to) corpus callosum, between BD patients and controls. This effect was most pronounced in lithium-free patients, implicating that lithium affects white matter microstructure and attenuates differences associated with bipolar disorder. Effects of antipsychotic medication intake were absent in FA and only subtle in MD relative to those of lithium. The abnormal white matter microstructure was associated with IQ but not specifically for either group.

Thousands of common single nucleotide polymorphisms (SNPs) are weakly associated with schizophrenia. It is likely that subsets of disease-associated SNPs are associated with distinct heritable disease-associated phenotypes. Therefore, we examined the shared genetic susceptibility modulating schizophrenia and brain volume.

Adolescence is a transitional period between childhood and adulthood and is characterized by emotional instability. Underlying this behavior may be an imbalance between the limbic subcortical areas and the prefrontal cortex. Here, we investigated differences in these regions during adolescence and young adulthood. Fifty subjects aged 10 to 24 viewed and rated neutral, negative, and positive pictures (IAPS: International Affective Picture System), while being scanned with functional MRI. Only those trials in which there was a match between the subject's response and the IAPS rating were included in the analyses. Task performance (matching accuracy, reaction times) did not differ across age. Activity in the amygdala and hippocampus decreased with age when processing emotional salient stimuli versus neutral stimuli. In contrast, activation in the ventrolateral prefrontal cortex increased with age. Importantly, we show for the first time that these age-related changes are paralleled by an increase in functional coupling of the amygdala and hippocampus with the orbitofrontal cortex and ventrolateral prefrontal cortex. These findings are in line with the general notion that brain development from childhood to adulthood is characterized by a gradual increase in frontal control over subcortical regions. Understanding these developmental changes is important as these may underlie typical adolescent behavior.

Schizophrenia has been conceptualized as a disorder of brain connectivity. Recent studies suggest that brain connectivity may be disproportionally impaired among the so-called rich club. This small core of densely interconnected hub regions has been hypothesized to form an important infrastructure for global brain communication and integration of information across different systems of the brain. Given the heritable nature of the illness, we hypothesized that connectivity disturbances, including abnormal rich club connectivity, may be related to familial vulnerability for schizophrenia. To test this hypothesis, both schizophrenia patients and unaffected siblings of patients were investigated. Rich club organization was examined in networks derived from diffusion-weighted imaging in 40 schizophrenia patients, 54 unaffected siblings of patients, and 51 healthy control subjects. Connectivity between rich club hubs was differentially reduced across groups (P = .014), such that it was highest in controls, intermediate in siblings (7.9% reduced relative to controls), and lowest in patients (19.6% reduced compared to controls). Furthermore, in patients, lower levels of rich club connectivity were found to be related to longer duration of illness and worse overall functioning. Together, these findings suggest that impaired rich club connectivity is related to familial, possibly reflecting genetic, vulnerability for schizophrenia. Our findings support a central role for abnormal rich club organization in the etiology of schizophrenia.

Genetic polymorphisms can shape the global landscape of DNA methylation, by either changing substrates for DNA methyltransferases or altering the DNA binding affinity of cis-regulatory proteins. The interactions between CpG methylation and genetic polymorphisms have been previously investigated by methylation quantitative trait loci (mQTL) and allele-specific methylation (ASM) analysis. However, it remains unclear whether these approaches can effectively and comprehensively identify all genetic variants that contribute to the inter-individual variation of DNA methylation levels. Here we used three independent approaches to systematically investigate the influence of genetic polymorphisms on variability in DNA methylation by characterizing the methylation state of 96 whole blood samples in 52 parent-child trios from 22 nuclear pedigrees. We performed targeted bisulfite sequencing with padlock probes to quantify the absolute DNA methylation levels at a set of 411,800 CpG sites in the human genome. With mid-parent offspring analysis (MPO), we identified 10,593 CpG sites that exhibited heritable methylation patterns, among which 70.1% were SNPs directly present in methylated CpG dinucleotides. We determined the mQTL analysis identified 49.9% of heritable CpG sites for which regulation occurred in a distal cis-regulatory manner, and that ASM analysis was only able to identify 5%. Finally, we identified hundreds of clusters in the human genome for which the degree of variation of CpG methylation, as opposed to whether or not CpG sites were methylated, was associated with genetic polymorphisms, supporting a recent hypothesis on the genetic influence of phenotypic plasticity. These results show that cis-regulatory SNPs identified by mQTL do not comprise the full extent of heritable CpG methylation, and that ASM appears overall unreliable. Overall, the extent of genome-methylome interactions is well beyond what is detectible with the commonly used mQTL and ASM approaches, and is likely to include effects on plasticity.

Genome-wide association studies have identified several common risk loci for schizophrenia (SCZ). In parallel, neuroimaging studies have shown consistent findings of widespread white matter disconnectivity in patients with SCZ.

Cytosine-5 methylation within CpG dinucleotides is a potentially important mechanism of epigenetic influence on human traits and disease. In addition to influences of age and gender, genetic control of DNA methylation levels has recently been described. We used whole blood genomic DNA in a twin set (23 MZ twin-pairs and 23 DZ twin-pairs, N = 92) as well as healthy controls (N = 96) to investigate heritability and relationship with age and gender of selected DNA methylation profiles using readily commercially available GoldenGate bead array technology. Despite the inability to detect meaningful methylation differences in the majority of CpG loci due to tissue type and locus selection issues, we found replicable significant associations of DNA methylation with age and gender. We identified associations of genetically heritable single nucleotide polymorphisms with large differences in DNA methylation levels near the polymorphism (cis effects) as well as associations with much smaller differences in DNA methylation levels elsewhere in the human genome (trans effects). Our results demonstrate the feasibility of array-based approaches in studies of DNA methylation and highlight the vast differences between individual loci. The identification of CpG loci of which DNA methylation levels are under genetic control or are related to age or gender will facilitate further studies into the role of DNA methylation and disease.

Different lines of evidence support a causal role for microglia in the pathogenesis of schizophrenia. However, how schizophrenia patient-derived microglia are affected at the phenotypic and functional level is still largely unknown. We used a recently described model to induce patient-derived microglia-like cells and used this to analyze changes in the molecular phenotype and function of myeloid cells in schizophrenia. We isolated monocytes from twenty recent-onset schizophrenia patients and twenty non-psychiatric controls. We cultured the cells towards an induced microglia-like phenotype (iMG), analyzed the phenotype of the cells by RNA sequencing and mass cytometry, and their response to LPS. Mass cytometry showed a high heterogeneity of iMG in cells derived from patients as well as controls. The prevalence of two iMG clusters was significantly higher in schizophrenia patients (adjusted p-value < 0.001). These subsets are characterized by expression of ApoE, Ccr2, CD18, CD44, and CD95, as well as IRF8, P2Y12, Cx3cr1 and HLA-DR. In addition, we found that patient-derived iMG show an enhanced response to LPS, with increased secretion of TNF-α. Further studies are needed to replicate these findings, to determine whether similar subclusters are present in schizophrenia patients in vivo, and to address how these subclusters are related to the increased response to LPS, as well as other microglial functions.

Microglia, the immune cells of the brain, are important for neurodevelopment and have been hypothesized to play a role in the pathogenesis of schizophrenia (SCZ). Although previous postmortem studies pointed toward presence of microglial activation, this view has been challenged by more recent hypothesis-driven and hypothesis-free analyses. The aim of the present study is to further understand the observed microglial changes in SCZ. We first performed a detailed meta-analysis on studies that analyzed microglial cell density, microglial morphology, and expression of microglial-specific markers. We then further explored findings from the temporal cortex by performing immunostainings and qPCRs on an additional dataset. A random effect meta-analysis showed that the density of microglial cells was unaltered in SCZ (ES: 0.144 95% CI: 0.102 to 0.390, p = .250), and clear changes in microglial morphology were also absent. The expression of several microglial specific genes, such as CX3CR1, CSF1R, IRF8, OLR1, and TMEM119 was decreased in SCZ (ES: -0.417 95% CI: -0.417 to -0.546, p < .0001), consistent with genome-wide transcriptome meta-analysis results. These results indicate a change in microglial phenotype rather than density, which was validated with the use of TMEM119/Iba1 immunostainings on temporal cortex of a separate cohort. Changes in microglial gene expression were overlapping between SCZ and other psychiatric disorders, but largely opposite from changes reported in Alzheimer's disease. This distinct microglial phenotype provides a crucial molecular hallmark for future research into the role of microglia in SCZ and other psychiatric disorders.

Suicide is a leading cause of death worldwide, and nonfatal suicide attempts, which occur far more frequently, are a major source of disability and social and economic burden. Both have substantial genetic etiology, which is partially shared and partially distinct from that of related psychiatric disorders.

Genomic loci where recurrent pathogenic copy number variants are associated with psychiatric phenotypes in the population may also be sensitive to the collective impact of multiple functional low-frequency single nucleotide variants (SNVs).

DNA methylation has recently been identified as a mediator between in utero famine exposure and a range of metabolic and psychiatric traits. However, genome-wide analyses are scarce and cross-sectional analyses are hampered by many potential confounding factors. Moreover, causal relations are hard to identify due to the lack of controlled experimental designs. In the current study, we therefore combined a comprehensive assessment of genome-wide DNA methylation differences in people exposed to the great Chinese famine in utero with an in vitro study in which we deprived fibroblasts of nutrition.

Bipolar disorder is a heritable mental illness with complex etiology. We performed a genome-wide association study of 41,917 bipolar disorder cases and 371,549 controls of European ancestry, which identified 64 associated genomic loci. Bipolar disorder risk alleles were enriched in genes in synaptic signaling pathways and brain-expressed genes, particularly those with high specificity of expression in neurons of the prefrontal cortex and hippocampus. Significant signal enrichment was found in genes encoding targets of antipsychotics, calcium channel blockers, antiepileptics and anesthetics. Integrating expression quantitative trait locus data implicated 15 genes robustly linked to bipolar disorder via gene expression, encoding druggable targets such as HTR6, MCHR1, DCLK3 and FURIN. Analyses of bipolar disorder subtypes indicated high but imperfect genetic correlation between bipolar disorder type I and II and identified additional associated loci. Together, these results advance our understanding of the biological etiology of bipolar disorder, identify novel therapeutic leads and prioritize genes for functional follow-up studies.

Left-right asymmetry is an important organizing feature of the healthy brain that may be altered in schizophrenia, but most studies have used relatively small samples and heterogeneous approaches, resulting in equivocal findings. We carried out the largest case-control study of structural brain asymmetries in schizophrenia, with MRI data from 5,080 affected individuals and 6,015 controls across 46 datasets, using a single image analysis protocol. Asymmetry indexes were calculated for global and regional cortical thickness, surface area, and subcortical volume measures. Differences of asymmetry were calculated between affected individuals and controls per dataset, and effect sizes were meta-analyzed across datasets. Small average case-control differences were observed for thickness asymmetries of the rostral anterior cingulate and the middle temporal gyrus, both driven by thinner left-hemispheric cortices in schizophrenia. Analyses of these asymmetries with respect to the use of antipsychotic medication and other clinical variables did not show any significant associations. Assessment of age- and sex-specific effects revealed a stronger average leftward asymmetry of pallidum volume between older cases and controls. Case-control differences in a multivariate context were assessed in a subset of the data (N = 2,029), which revealed that 7% of the variance across all structural asymmetries was explained by case-control status. Subtle case-control differences of brain macrostructural asymmetry may reflect differences at the molecular, cytoarchitectonic, or circuit levels that have functional relevance for the disorder. Reduced left middle temporal cortical thickness is consistent with altered left-hemisphere language network organization in schizophrenia.

There is currently no quantifiable method to predict long-term clinical outcomes in patients presenting with a first episode of psychosis. A major barrier to developing useful markers for this is biological heterogeneity, where many different pathological mechanisms may underly the same set of symptoms in different individuals. Normative modelling has been used to quantify this heterogeneity in established psychotic disorders by identifying regions of the cortex which are thinner than expected based on a normative healthy population range. These brain atypicalities are measured at the individual level and therefore potentially useful in a clinical setting. However, it is still unclear whether alterations in individual brain structure can be detected at the time of the first psychotic episode, and whether they are associated with subsequent clinical outcomes. We applied normative modelling of cortical thickness to a sample of first-episode psychosis patients, with the aim of quantifying heterogeneity and to use any pattern of cortical atypicality to predict symptoms and response to antipsychotic medication at timepoints from baseline up to 95 weeks (median follow-ups = 4). T1-weighted brain magnetic resonance images from the GAP and OPTiMiSE samples were processed with Freesurfer V6.0.0 yielding 148 cortical thickness features. An existing normative model of cortical thickness (n = 37,126) was adapted to integrate data from each clinical site and account for effects of gender and site. Our test sample consisted of control participants (n = 149, mean age = 26, SD = 6.7) and patient data (n = 295, mean age = 26, SD = 6.7), this sample was used for estimating deviations from the normative model and subsequent statistical analysis. For each individual, the 148 cortical thickness features were mapped to centiles of the normative distribution and converted to z-scores reflecting the distance from the population mean. Individual cortical thickness metrics of +/- 2.6 standard deviations from the mean were considered extreme deviations from the norm. We found that no more than 6.4% of psychosis patients had extreme deviations in a single brain region (regional overlap) demonstrating a high degree of heterogeneity. Mann-Whitney U tests were run on z-scores for each region and significantly lower z-scores were observed in FEP patients in the frontal, temporal, parietal and occipital lobes. Finally, linear mixed-effects modelling showed that negative deviations in cortical thickness in parietal and temporal regions at baseline are related to more severe negative symptoms over the medium-term. This study shows that even at the early stage of symptom onset normative modelling provides a framework to identify individualised cortical markers which can be used for early personalised intervention and stratification.

In order to assess the importance of environmental and genetic risk on transition from health to psychotic disorder, a prospective study of individuals at average (n = 462) and high genetic risk (n = 810) was conducted.

Complex diseases are often highly heritable. However, for many complex traits only a small proportion of the heritability can be explained by observed genetic variants in traditional genome-wide association (GWA) studies. Moreover, for some of those traits few significant SNPs have been identified. Single SNP association methods test for association at a single SNP, ignoring the effect of other SNPs. We show using a simple multi-locus odds model of complex disease that moderate to large effect sizes of causal variants may be estimated as relatively small effect sizes in single SNP association testing. This underestimation effect is most severe for diseases influenced by numerous risk variants. We relate the underestimation effect to the concept of non-collapsibility found in the statistics literature. As described, continuous phenotypes generated with linear genetic models are not affected by this underestimation effect. Since many GWA studies apply single SNP analysis to dichotomous phenotypes, previously reported results potentially underestimate true effect sizes, thereby impeding identification of true effect SNPs. Therefore, when a multi-locus model of disease risk is assumed, a multi SNP analysis may be more appropriate.

Posttraumatic stress disorder (PTSD) is an anxiety disorder that is often diagnosed with comorbid depressive disorder. Therefore, neuroimaging studies investigating PTSD typically include both patients with and without comorbid depression. Differences in activity of the anterior cingulate cortex (ACC) and insula have been shown to differentiate PTSD patients with and without major depressive disorder (MDD). Whether or not comorbid MDD affects resting state functional connectivity of PTSD patients has not been investigated to our knowledge. Here, resting state functional connectivity of PTSD patients with (PTSD+MDD; n=27) and without (PTSD-MDD; n=23) comorbid MDD was investigated. The subgenual ACC and insula were investigated as seed regions. Connectivity between the subgenual ACC and perigenual parts of the ACC was increased in PTSD+MDD versus PTSD-MDD, which may reflect the presence of depressive specific symptoms such as rumination. Functional connectivity of the subgenual ACC with the thalamus was reduced, potentially related to more severe deficits in executive functioning in the PTSD+MDD group versus the PTSD-MDD group. In addition, the PTSD+MDD group showed reduced functional connectivity of the insula with the hippocampus compared to the PTSD-MDD group. However, this cluster was no longer significantly different when PTSD patients that were using medication were excluded from analyses. Thus, resting state functional connectivity of the subgenual ACC can distinguish PTSD+MDD from PTSD-MDD, and this may therefore be used as a neurobiological marker for comorbid MDD in the presence of PTSD. As PTSD+MDD are more treatment resistant, these findings can also guide treatment development, for example by targeting the subgenual ACC network with treatment.

The notion that healthy brain function emerges from coordinated neural activity constrained by the brain's network of anatomical connections--i.e., the connectome--suggests that alterations in the connectome's wiring pattern may underlie brain disorders. Corroborating this hypothesis, studies in schizophrenia are indicative of altered connectome architecture including reduced communication efficiency, disruptions of central brain hubs, and affected "rich club" organization. Whether similar deficits are present in bipolar disorder is currently unknown. This study examines structural connectome topology in 216 bipolar I disorder patients as compared to 144 healthy controls, focusing in particular on central regions (i.e., brain hubs) and connections (i.e., rich club connections, interhemispheric connections) of the brain's network. We find that bipolar I disorder patients exhibit reduced global efficiency (-4.4%, P =0.002) and that this deficit relates (r = 0.56, P < 0.001) to reduced connectivity strength of interhemispheric connections (-13.0%, P = 0.001). Bipolar disorder patients were found not to show predominant alterations in the strength of brain hub connections in general, or of connections spanning brain hubs (i.e., "rich club" connections) in particular (all P > 0.1). These findings highlight a role for aberrant brain network architecture in bipolar I disorder with reduced global efficiency in association with disruptions in interhemispheric connectivity, while the central "rich club" system appears not to be particularly affected.