Search

Schizophrenia is a mental disorder constituting different symptom clusters. Its high heterogeneity in both pathophysiology and clinical manifestations hampered effective prevention and treatment. It has long been recognized that one of the core features of schizophrenia is its intellectual decline. Using the proxy-phenotype method (PPM), we tried to identify core genes, the expression of which in the dorsal lateral prefrontal cortex (DLPFC) showed a genetic dependence between intelligence (IT) and schizophrenia (SCZ). The result revealed ten genes of genetic dependence in their genetic expression in DLPFC between IT and schizophrenia. Further, a clustering analysis using the expression matrix of these ten genes identified four biotypes in our patient group. Subsequent phenotypic profiling of these four biotypes indicated a significant difference in working memory capacity, the gray matter volume (GMV) of five brain regions (lLimbicA_TempPole_2, rLimbicA_TempPole_2 rLimbicB_OFC_1, rContA_IPS_1 and rContB_PFClv_1), structural network and psychopathology. An in-vitro investigation of the biological functions of these core genes indicated their potentially critical role in neuronal growth, especially in dendritic spines. Our current study employed a novel statistical approach to identify the core genes associated with IT and explore the possibility of using the expression knowledge of these core genes to reduce the heterogeneity of schizophrenia. The results pinpointed one biotype that exhibited significant deficits in working memory, GMV in limbic and prefrontal areas, and also showed psychopathology of core negative symptom and worse outcomes.

The last decades of research in schizophrenia witnessed a shift of etiological speculation from neurotransmitters to inflammation. However, identifying definite inflammatory effectors of schizophrenia remains elusive due to confounding factors such as medication and metabolic status. To tackle this issue, we carried out omnigenic-based Mendelian randomization (MR) analysis to explore the inflammatory responses of schizophrenia and the brain morphological consequences caused by these SCZ-triggering inflammation responses. Our results identified seven SCZ-triggering inflammation markers, with P values surviving the Bonferroni multiple comparisons (B_NGF, P = 1.45 × 10−8; GROA (CXCL1) P = 1.15 × 10−4; IL8, P = 3.64 × 10−7; MCSF, P = 9.30 × 10−4; MCP3 (CCL7), P = 1.3 × 10−6; TNF_β, P = 3.63 × 10−4; CRP, P = 1.71 × 10−32). Further, three of them, GROA (CXCL1), IL8 and CRP, could lead to significant linear change rate of brain morphologies, especially white matter in both cerebral and cerebellum. Our study is the first to use an omnigenic conceptual framework to capture the immune pathology of schizophrenia. Although future studies adopting a different methodology are needed to validate our results, our study provides another piece of evidence that extensive and low-grade neuroinflammation exists in schizophrenia and that some of these inflammation markers could be potential targets for the precise diagnosis and treatment of schizophrenia.

Deciphering the relationships between genes and complex traits could help us better understand the biological mechanisms leading to phenotypic variations and disease onset. Univariate gene-based analyses are widely used to characterize gene-phenotype relationships, but are subject to the influence of confounders. Furthermore, while some genes directly contribute to traits variations, others may exert their effects through other genes. How to quantify individual genes’ direct and indirect effects on complex traits remains an important yet challenging question. We presented a novel framework (BN-GWAS) to decipher the total and direct causal effects of individual genes using imputed gene expression data from GWAS and raw gene expression from GTEx. The study was partially motivated by the quest to differentiate “core” genes (genes with direct causal effect on the phenotype) from “peripheral” ones. Our proposed framework is based on a Bayesian network (BN) approach, which produces a directed graph showing the relationship between genes and the phenotype. The approach aims to uncover the overall causal structure, to examine the role of individual genes and quantify the direct and indirect effects by each gene. An important advantage and novelty of the proposed framework is that it allows gene expression and disease trait(s) to be evaluated in different samples, significantly improving the flexibility and applicability of the approach. It uses IDA and jointIDA incorporating a novel p-value-based regularization approach to quantify the causal effects (including total causal effects, direct causal effects, and medication effects) of genes. The proposed approach can be extended to decipher the joint causal network of 2 or more traits, and has high specificity and precision (a.k.a., positive predictive value), making it particularly useful for selecting genes for follow-up studies. We verified the feasibility and validity of the proposed framework by extensive simulations and applications to 52 traits across 14 tissues in the UK Biobank (UKBB). Split-half replication and stability selection analyses were performed to demonstrate the accuracy and efficiency of our proposed method to identify causally relevant genes. The identified (direct) causal genes were found to be significantly enriched for genes highlighted in the OpenTargets database, and the enrichment was stronger than achieved by conventional univariate gene-based tests. Encouragingly, many enriched pathways were supported by the literature, and some of the enriched drugs have been tested or used to treat patients in clinical practice. Our proposed framework provides a powerful way to prioritize genes with large direct or indirect causal effects and to estimate the ‘variable importance’ of such genes.

Studies conducted in psychotic disorders have shown that DNA-methylation (DNAm) is sensitive to the impact of Childhood Adversity (CA). However, whether it mediates the association between CA and psychosis is yet to be explored. Epigenome-wide association studies (EWAS) using the Illumina Infinium-Methylation EPIC array in peripheral blood tissue from 366 First-episode of psychosis and 517 healthy controls was performed. Polyvictimization scores were created for abuse, neglect, and composite adversity with the Childhood Trauma Questionnaire (CTQ). Regressions examining (I) CTQ scores with psychosis; (II) with DNAm EWAS level and (III), then between DNAm and caseness, adjusted for a variety of confounders were conducted. Divide-Aggregate Composite-null Test for the composite null hypothesis of no mediation effect was conducted. Enrichment analyses were conducted with missMethyl package and the KEGG database. Our results show that CA was associated with psychosis (Composite: OR = 3.09; p =

Despite many methodological studies on the global heritability estimation of complex traits, very few methods have been designed to dissect local heritability accurately. A precise local heritability would provide a high-resolution map for precision medicine. We report an effective heritability estimator (EHE) for local heritability estimation using p values from genome-wide association studies (GWAS). With a directly converting marginal heritability to the effective one, EHE shows higher accuracy and is the only method correctly calculating standard errors of local heritability estimation among the four compared methods. Importantly, EHE can be directly used to estimate the conditional heritability of nearby regions or genes in which the redundant heritability can be cleaned. The conditional heritability estimation procedure can be supervised by tissue-specific expression profiles to prioritize and quantify functionally more important genes of complex phenotypes. Applying EHE to 42 complex traits from the UK biobank, we found that heritability density in lncRNA genes tends to be higher than protein-coding genes, and genes pleiotropic to more phenotypes are not enriched for more heritability among candidate susceptibility genes. EHE provides an efficient way to finely and precisely dissect the genetic architecture of complex phenotypes.

Hirschsprung disease (HSCR) is characterized by the absence of enteric neurons caused by the defects of enteric neural crest cells (ENCCs), leading to intestinal obstruction. Here, using iPSC-based models of HSCR and single-cell transcriptomic analysis, we identified a core gene set of 118 genes commonly dysregulated in all HSCR-ENCCs, with HDAC1 found to be a master regulator of these genes. Furthermore, upregulation of RNA splicing mediators and enhanced alternative splicing events were associated with severe form of HSCR. In particular, the higher inclusion rate of exon 9 in PTBP1 and the perturbed expression of a PTBP1-target, PKM, were significantly enriched in these HSCR-ENCCs, and associated with the defective oxidative-phosphorylation (OXPHOS) and impaired neurogenesis. Hedgehog-induced OXPHOS significantly enhanced the survival and differentiation capacity of HSCR-ENCCs. In sum, we have defined the core factors underpinning HSCR disease and demonstrated the implications of OXPHOS in ENCC development and HSCR pathogenesis.

Background Jumping to conclusions (JTC) is a well-established reasoning and data gathering bias found in patients with psychosis even at illness onset (First Episode Psychosis, FEP). Preliminary work in this field focused primarily on the association with delusions, although jumping to conclusions has also been found in non-deluded schizophrenia patients after remission, and in individual with at risk mental state. Moreover, psychotic patients tend to show impairments in social cognition, struggling in identifying, processing and interpreting social clues. Deficits in facial emotion recognition (FER) – a key component of the construct – represent a well-replicated finding in schizophrenia. Furthermore, deficits in global facial affect recognition have been found in FEP with the same severity as at further stages, especially for anger recognition. The present study aims to measure JTC and FER bias in a sample of FEP recruited across 5 European countries, compared with healthy controls. Methods Data on JTC (Beads task 60:40), FER (Degraded Facial Recognition task – DFAR) and socio-demographics have been analysed in a sample of 643 FEP and 1019 population controls recruited as part as the EU-GEI study across UK, Netherlands, France, Spain, and Italy. IQ scores were used to exclude cases and controls with current IQ

Background Current diagnostic models of psychosis have been questioned since Kraepelin’s original dichotomy of dementia praecox and manic depression. Indeed, increasing evidence has suggested that a dimensional approach might be a valid alternative platform for research. However, while an increasing number of studies have investigated how environmental risk factors for affective and non-affective psychosis map onto symptom dimensions, only a few have examined these dimensions in relation to genetic variants as summarised by Polygenic Risk Score (PRS). Furthermore, no studies have examined the putative effect of PRS for Schizophrenia (SZ), Bipolar Disorder (BP), and Major Depressive Disorder (MDD) on previously identified general and specific symptom dimensions. At the same time, only one study has investigated how symptoms vary according to epidemiological factors such as living in urban neighbourhoods. The objectives of this study were to: 1) test whether a bi-factor model statistically fits the conceptualization of psychosis as composed of general and specific dimensions; 2) examine the extent to which SZ, BP, and MDD PRSs explain the phenotypic variance due to general and specific dimensions; 3) test the hypothesis that the general psychosis dimension would be more severe in highly urban environments. Methods We used clinical and epidemiological data from the EUropean network of national schizophrenia networks studying Gene-Environment Interactions (EUGEI) study, including 2322 First Episode Psychosis (FEP) patients recruited in 17 sites across 6 countries. Genetic variants were collectively analyzed for 800 individuals. The following analysis steps were performed: 1) Psychopathology items were analysed using multidimensional item response modelling in MPlus to estimate unidimensional, multidimensional, and bi-factor models of psychosis. Model fit statistics included Log-Likelihood, and Akaike and Bayesian Information Criteria to compare these models. 2) SZ, BP, and MDD PRSs for general and specific dimensions were built using PRSice. Summary statistics from large case-control mega-analyses from the Psychiatric Genomics Consortium were used as base data sets and general and specific dimension scores were used as discovery data sets. Individuals’ number of risk alleles in the discovery sample was weighted by the log odds ratio from the base samples, accounting for population stratification, and summed into the three PRSs. 3) Multilevel regression analysis was used in STATA 14 to examine the variance in general dimension due to the population density levels across the sites. Results A bi-factor solution, composed of one general and five specific symptom dimensions, showed the best model fit statistics. Higher SZ PRS score was associated with higher scores on positive dimensions (β= 0.27, t=2.11, p

Background: Robust evidence has demonstrated that cannabis use increases the risk to develop psychotic disorders. However, a limited number of studies have investigated if and how cannabis use influences psychopathology profiles at first episode psychosis (FEP). Based on the evidence that dopamine dysfunction contributes to explain positive symptoms in psychosis, and that the main cannabis’ psychoactive component, Δ9-Tetrahydrocannabinol (THC), modulates the dopamine system, we hypothesise that: 1) positive symptoms at FEP are more common among psychotic patients who used cannabis compared with never users; 2) this association is a dose-response relationship. Methods: We analyzed a sample of 1130 FEP patients as part of the EUGEI study, recruited across six countries. The MRC Socio-demographic Schedule was used to collect sociodemographic information. Psychopathology was assessed with the OPerational CRITeria (OPCRIT), and symptom items were analyzed using Mplus to estimate a multidimensional model of psychosis. The Cannabis Experience Questionnaire modified version (CEQmv) was administered to collect information on cannabis, and different patterns of use were computed based on frequency of consumption and type of cannabis, as a proxy of exposure to THC. Results: The lifetime rate of cannabis use was 63%. Fifty-five percent of cannabis users consumed mostly high-potency cannabis, and 46% showed a daily frequency. Mixed-effects linear regression revealed that frequency of cannabis use was associated with the positive symptom dimension score. Daily users of high-potency cannabis presented with the strongest association (Β=0.19, 95%CI=0.02–0.38), even after gender, age, ethnicity, other drug use, and study site were controlled for. Discussion: Our results show that patients with a history of daily use of high potency cannabis express more positive symptoms at psychosis onset, even after taking into account other substance use and relevant sociodemographic factors.

Our study aimed to identify possible single nucleotide polymorphisms (SNPs) via a genome-wide association study (GWAS) approach and a candidate gene platform that were associated with lumbar developmental spinal stenosis (DSS). Southern Chinese population-based study volunteers were assessed (age range: 18-55 years). DSS was defined as the anteroposterior bony spinal canal diameter on T1-weighted axial MRI of L1 to S1. Genotyping was performed using the Illumina HumanOmniZhongHua-8 BeadChip. Using the canal diameter as the quantitative trait, genomic statistical analyses was performed. A total of 469 subjects were recruited. The mean axial AP measurements noted were: L1: 21.8 mm, L2: 21.9 mm, L3: 22.4 mm, L4: 20.2 mm, L5: 19.6 mm, and S1: 17.3 mm. Q-Q plots of genome-wide associations found significant differences in L4 and L5 measurements. More significant SNPs were found on chromosomes 8, 11, and 18. Low-density lipoprotein receptor-related protein 5 on chromosome 11 was found to be an important functional gene in canal bony development via candidate gene approach. We found two clusters in the findings with one including the upper levels (L1-L4) and the other the lower levels (L5 and S1). This is the first GWAS addressing DSS. The presence of multiple SNPs suggests a multi-factorial origin of DSS. Further analyses noted region-specific genetic predisposition, delineating distinct upper to lower lumbar regions of DSS. With better understanding of the DSS phenotype and genetic markers, the at-risk population can be identified early, preventative measures can be initiated, lifestyle/activity modification can be implemented, and more novel and precision-based therapeutics can be developed. © 2017 Orthopaedic Research Society. Published by Wiley Periodicals, Inc. J Orthop Res 36:1262-1268, 2018.

Supplemental Figures and Tables. Supplemental Figures S1â S14 and Tables S1â S14. (PDF 3992 kb)

Supplementary information. (DOC 653 kb)

Supplementary Methods and Material. (DOCX 51Â kb)