Attention-deficit/hyperactivity disorder and dopamine receptor D4 (DRD4) exon 3 variable number of tandem repeats (VNTR) 2-repeat allele.

To investigate the association of attention-deficit/hyperactivity disorder (ADHD) with the 48-base pair (bp) variable number of tandem repeats (VNTR) in exon 3 of the dopamine receptor D4 (DRD4) gene, we genotyped 240 ADHD patients and their parents from Hong Kong. The 4R allele was most common, followed by 2R. We examined association between the 2R allele (relative to 4R) and ADHD by Transmission Disequilibrium Test (TDT). The odds ratio (OR) (95% confidence interval) was 0.90 (0.64-1.3). The p-value was 0.6. Examining subgroups revealed nominally significant association of 2R with inattentive ADHD: OR = 0.33 (0.12-0.92) and p = 0.03. Because our study used TDT analysis, we meta-analyzed the association of 2R with ADHD in Asians (1329 patient alleles), revealing results similar to ours: OR = 0.97 (0.80-1.2) and p = 0.8. To examine the association of 2R with inattentive ADHD, we meta-analyzed all studies (regardless of analysis type or ethnicity, in order to increase statistical power): 702 patient alleles, 1420 control alleles, OR = 0.81 (0.57-1.1) and p = 0.2. Overall, there is no evidence of association between ADHD and the 2R allele, but the suggestive association with the inattentive type warrants further investigation.

Mendelian randomization: causal inference leveraging genetic data.

Mendelian randomization (MR) leverages genetic information to examine the causal relationship between phenotypes allowing for the presence of unmeasured confounders. MR has been widely applied to unresolved questions in epidemiology, making use of summary statistics from genome-wide association studies on an increasing number of human traits. However, an understanding of essential concepts is necessary for the appropriate application and interpretation of MR. This review aims to provide a non-technical overview of MR and demonstrate its relevance to psychiatric research. We begin with the origins of MR and the reasons for its recent expansion, followed by an overview of its statistical methodology. We then describe the limitations of MR, and how these are being addressed by recent methodological advances. We showcase the practical use of MR in psychiatry through three illustrative examples - the connection between cannabis use and psychosis, the link between intelligence and schizophrenia, and the search for modifiable risk factors for depression. The review concludes with a discussion of the prospects of MR, focusing on the integration of multi-omics data and its extension to delineating complex causal networks.

Whole-exome sequencing in a Chinese sample provides preliminary evidence for the link between rare/low-frequency immune-related variants and early-onset schizophrenia.

Rare and low-frequency variants contribute to schizophrenia (SCZ), and may influence its age-at-onset (AAO). We examined the association of rare or low-frequency deleterious coding variants in Chinese patients with SCZ. We collected DNA samples in 197 patients with SCZ spectrum disorder and 82 healthy controls (HC), and performed exome sequencing. The AAO variable was ascertained in the majority of SCZ participants for identify the early-onset (EOS, AAO<=18) and adult-onset (AOS, AAO>18) subgroups. We examined the overall association of rare/low-frequency, damaging variants in SCZ versus HC, EOS versus HC, and AOS versus HC at the gene and gene-set levels using Sequence Kernel Association Test. The quantitative rare-variant association test of AAO was conducted. Resampling was used to obtain empirical p-values and to control for family-wise error rate (FWER). In binary-trait association tests, we identified 5 potential candidate risk genes and 10 gene ontology biological processes (GOBP) terms, among which PADI2 reached FWER-adjusted significance. In quantitative rare-variant association tests, we found marginally significant correlations of AAO with alterations in 4 candidate risk genes, and 5 GOBP pathways. Together, the biological and functional profiles of these genes and gene sets supported the involvement of perturbations of neural systems in SCZ, and altered immune functions in EOS.

Computational Retinal Microvascular Biomarkers from an OCTA Image in Clinical Investigation.

Retinal structural and functional changes in humans can be manifestations of different physiological or pathological conditions. Retinal imaging is the only way to directly inspect blood vessels and their pathological changes throughout the whole body non-invasively. Various quantitative analysis metrics have been used to measure the abnormalities of retinal microvasculature in the context of different retinal, cerebral and systemic disorders. Recently developed optical coherence tomography angiography (OCTA) is a non-invasive imaging tool that allows high-resolution three-dimensional mapping of the retinal microvasculature. The identification of retinal biomarkers from OCTA images could facilitate clinical investigation in various scenarios. We provide a framework for extracting computational retinal microvasculature biomarkers (CRMBs) from OCTA images through a knowledge-driven computerized automatic analytical system. Our method allows for improved identification of the foveal avascular zone (FAZ) and introduces a novel definition of vessel dispersion in the macular region. Furthermore, retinal large vessels and capillaries of the superficial and deep plexus can be differentiated, correlating with retinal pathology. The diagnostic value of OCTA CRMBs was demonstrated by a cross-sectional study with 30 healthy subjects and 43 retinal vein occlusion (RVO) patients, which identified strong correlations between OCTA CRMBs and retinal function in RVO patients. These OCTA CRMBs generated through this "all-in-one" pipeline may provide clinicians with insights about disease severity, treatment response and prognosis, aiding in the management and early detection of various disorders.

Association of neurotransmitter pathway polygenic risk with specific symptom profiles in psychosis.

A primary goal of psychiatry is to better understand the pathways that link genetic risk to psychiatric symptoms. Here, we tested association of diagnosis and endophenotypes with overall and neurotransmitter pathway-specific polygenic risk in patients with early-stage psychosis. Subjects included 205 demographically diverse cases with a psychotic disorder who underwent comprehensive psychiatric and neurological phenotyping and 115 matched controls. Following genotyping, we calculated polygenic scores (PGSs) for schizophrenia (SZ) and bipolar disorder (BP) using Psychiatric Genomics Consortium GWAS summary statistics. To test if overall genetic risk can be partitioned into affected neurotransmitter pathways, we calculated pathway PGSs (pPGSs) for SZ risk affecting each of four major neurotransmitter systems: glutamate, GABA, dopamine, and serotonin. Psychosis subjects had elevated SZ PGS versus controls; cases with SZ or BP diagnoses had stronger SZ or BP risk, respectively. There was no significant association within psychosis cases between individual symptom measures and overall PGS. However, neurotransmitter-specific pPGSs were moderately associated with specific endophenotypes; notably, glutamate was associated with SZ diagnosis and with deficits in cognitive control during task-based fMRI, while dopamine was associated with global functioning. Finally, unbiased endophenotype-driven clustering identified three diagnostically mixed case groups that separated on primary deficits of positive symptoms, negative symptoms, global functioning, and cognitive control. All clusters showed strong genome-wide risk. Cluster 2, characterized by deficits in cognitive control and negative symptoms, additionally showed specific risk concentrated in glutamatergic and GABAergic pathways. Due to the intensive characterization of our subjects, the present study was limited to a relatively small cohort. As such, results should be followed up with additional research at the population and mechanism level. Our study suggests pathway-based PGS analysis may be a powerful path forward to study genetic mechanisms driving psychiatric endophenotypes.

Temporal changes in brain morphology related to inflammation and schizophrenia: an omnigenic Mendelian randomization study.

BACKGROUND: Over the past several decades, more research focuses have been made on the inflammation/immune hypothesis of schizophrenia. Building upon synaptic plasticity hypothesis, inflammation may contribute the underlying pathophysiology of schizophrenia. Yet, pinpointing the specific inflammatory agents responsible for schizophrenia remains a complex challenge, mainly due to medication and metabolic status. Multiple lines of evidence point to a wide-spread genetic association across genome underlying the phenotypic variations of schizophrenia. METHOD: We collected the latest genome-wide association analysis (GWAS) summary data of schizophrenia, cytokines, and longitudinal change of brain. We utilized the omnigenic model which takes into account all genomic SNPs included in the GWAS of trait, instead of traditional Mendelian randomization (MR) methods. We conducted two round MR to investigate the inflammatory triggers of schizophrenia and the resulting longitudinal changes in the brain. RESULTS: We identified seven inflammation markers linked to schizophrenia onset, which all passed the Bonferroni correction for multiple comparisons (bNGF, GROA(CXCL1), IL-8, M-CSF, MCP-3 (CCL7), TNF-ß, CRP). Moreover, CRP were found to significantly influence the linear rate of brain morphology changes, predominantly in the white matter of the cerebrum and cerebellum. CONCLUSION: With an omnigenic approach, our study sheds light on the immune pathology of schizophrenia. Although these findings need confirmation from future studies employing different methodologies, our work provides substantial evidence that pervasive, low-level neuroinflammation may play a pivotal role in schizophrenia, potentially leading to notable longitudinal changes in brain morphology.

INSIDER: Interpretable sparse matrix decomposition for RNA expression data analysis.

RNA sequencing (RNA-Seq) is widely used to capture transcriptome dynamics across tissues, biological entities, and conditions. Currently, few or no methods can handle multiple biological variables (e.g., tissues/ phenotypes) and their interactions simultaneously, while also achieving dimension reduction (DR). We propose INSIDER, a general and flexible statistical framework based on matrix factorization, which is freely available at https://github.com/kai0511/insider. INSIDER decomposes variation from different biological variables and their interactions into a shared low-rank latent space. Particularly, it introduces the elastic net penalty to induce sparsity while considering the grouping effects of genes. It can achieve DR of high-dimensional data (of > = 3 dimensions), as opposed to conventional methods (e.g., PCA/NMF) which generally only handle 2D data (e.g., sample × expression). Besides, it enables computing 'adjusted' expression profiles for specific biological variables while controlling variation from other variables. INSIDER is computationally efficient and accommodates missing data. INSIDER also performed similarly or outperformed a close competing method, SDA, as shown in simulations and can handle complex missing data in RNA-Seq data. Moreover, unlike SDA, it can be used when the data cannot be structured into a tensor. Lastly, we demonstrate its usefulness via real data analysis, including clustering donors for disease subtyping, revealing neuro-development trajectory using the BrainSpan data, and uncovering biological processes contributing to variables of interest (e.g., disease status and tissue) and their interactions.

Improved estimation of functional enrichment in SNP heritability using feasible generalized least squares.

Functional enrichment results typically implicate tissue or cell-type-specific biological pathways in disease pathogenesis and as therapeutic targets. We propose generalized linkage disequilibrium score regression (g-LDSC) that requires only genome-wide association studies (GWASs) summary-level data to estimate functional enrichment. The method adopts the same assumptions and regression model formulation as stratified linkage disequilibrium score regression (s-LDSC). Although s-LDSC only partially uses LD information, our method uses the whole LD matrix, which accounts for possible correlated error structure via a feasible generalized least-squares estimation. We demonstrate through simulation studies under various scenarios that g-LDSC provides more precise estimates of functional enrichment than s-LDSC, regardless of model misspecification. In an application to GWAS summary statistics of 15 traits from the UK Biobank, estimates of functional enrichment using g-LDSC were lower and more realistic than those obtained from s-LDSC. In addition, g-LDSC detected more significantly enriched functional annotations among 24 functional annotations for the 15 traits than s-LDSC (118 vs. 51).

Bidirectional two-sample Mendelian randomization study of differential white blood cell counts and schizophrenia.

BACKGROUND: Schizophrenia and white blood cell counts (WBC) are both complex and polygenic traits. Previous evidence suggests that increased WBC are associated with higher all-cause mortality, and other studies have found elevated WBC in first-episode psychosis and chronic schizophrenia. However, these observational findings may be confounded by antipsychotic exposures and their effects on WBC. Mendelian randomization (MR) is a useful method for examining the directions of genetically-predicted relationships between schizophrenia and WBC. METHODS: We performed a two-sample MR using summary statistics from genome-wide association studies (GWAS) conducted by the Psychiatric Genomics Consortium Schizophrenia Workgroup (N = 130,644) and the Blood Cell Consortium (N = 563,946). The MR methods included inverse variance weighted (IVW), MR Egger, weighted median, MR-PRESSO, contamination mixture, and a novel approach called mixture model reciprocal causal inference (MRCI). False discovery rate was employed to correct for multiple testing. RESULTS: Multiple MR methods supported bidirectional genetically-predicted relationships between lymphocyte count and schizophrenia: IVW (b = 0.026; FDR p-value = 0.008), MR Egger (b = 0.026; FDR p-value = 0.008), weighted median (b = 0.013; FDR p-value = 0.049), and MR-PRESSO (b = 0.014; FDR p-value = 0.010) in the forward direction, and IVW (OR = 1.100; FDR p-value = 0.021), MR Egger (OR = 1.231; FDR p-value < 0.001), weighted median (OR = 1.136; FDR p-value = 0.006) and MRCI (OR = 1.260; FDR p-value = 0.026) in the reverse direction. MR Egger (OR = 1.171; FDR p-value < 0.001) and MRCI (OR = 1.154; FDR p-value = 0.026) both suggested genetically-predicted eosinophil count is associated with schizophrenia, but MR Egger (b = 0.060; FDR p-value = 0.010) and contamination mixture (b = -0.013; FDR p-value = 0.045) gave ambiguous results on whether genetically predicted liability to schizophrenia would be associated with eosinophil count. MR Egger (b = 0.044; FDR p-value = 0.010) and MR-PRESSO (b = 0.009; FDR p-value = 0.045) supported genetically predicted liability to schizophrenia is associated with elevated monocyte count, and the opposite direction was also indicated by MR Egger (OR = 1.231; FDR p-value = 0.045). Lastly, unidirectional genetic liability from schizophrenia to neutrophil count were proposed by MR-PRESSO (b = 0.011; FDR p-value = 0.028) and contamination mixture (b = 0.011; FDR p-value = 0.045) method. CONCLUSION: This MR study utilised multiple MR methods to obtain results suggesting bidirectional genetic genetically-predicted relationships for elevated lymphocyte counts and schizophrenia risk. In addition, moderate evidence also showed bidirectional genetically-predicted relationships between schizophrenia and monocyte counts, and unidirectional effect from genetic liability for eosinophil count to schizophrenia and from genetic liability for schizophrenia to neutrophil count. The influence of schizophrenia to eosinophil count is less certain. Our findings support the role of WBC in schizophrenia and concur with the hypothesis of neuroinflammation in schizophrenia.

A Bayesian fine-mapping model using a continuous global-local shrinkage prior with applications in prostate cancer analysis.

The aim of fine mapping is to identify genetic variants causally contributing to complex traits or diseases. Existing fine-mapping methods employ Bayesian discrete mixture priors and depend on a pre-specified maximum number of causal variants, which may lead to sub-optimal solutions. In this work, we propose a Bayesian fine-mapping method called h2-D2, utilizing a continuous global-local shrinkage prior. We also present an approach to define credible sets of causal variants in continuous prior settings. Simulation studies demonstrate that h2-D2 outperforms current state-of-the-art fine-mapping methods such as SuSiE and FINEMAP in accurately identifying causal variants and estimating their effect sizes. We further applied h2-D2 to prostate cancer analysis and discovered some previously unknown causal variants. In addition, we inferred 369 target genes associated with the detected causal variants and several pathways that were significantly over-represented by these genes, shedding light on their potential roles in prostate cancer development and progression.

SumVg: Total Heritability Explained by All Variants in Genome-Wide Association Studies Based on Summary Statistics with Standard Error Estimates.

Genome-wide association studies (GWAS) are commonly employed to study the genetic basis of complex traits/diseases, and a key question is how much heritability could be explained by all single nucleotide polymorphisms (SNPs) in GWAS. One widely used approach that relies on summary statistics only is linkage disequilibrium score regression (LDSC); however, this approach requires certain assumptions about the effects of SNPs (e.g., all SNPs contribute to heritability and each SNP contributes equal variance). More flexible modeling methods may be useful. We previously developed an approach recovering the "true" effect sizes from a set of observed z-statistics with an empirical Bayes approach, using only summary statistics. However, methods for standard error (SE) estimation are not available yet, limiting the interpretation of our results and the applicability of the approach. In this study, we developed several resampling-based approaches to estimate the SE of SNP-based heritability, including two jackknife and three parametric bootstrap methods. The resampling procedures are performed at the SNP level as it is most common to estimate heritability from GWAS summary statistics alone. Simulations showed that the delete-d-jackknife and parametric bootstrap approaches provide good estimates of the SE. In particular, the parametric bootstrap approaches yield the lowest root-mean-squared-error (RMSE) of the true SE. We also explored various methods for constructing confidence intervals (CIs). In addition, we applied our method to estimate the SNP-based heritability of 12 immune-related traits (levels of cytokines and growth factors) to shed light on their genetic architecture. We also implemented the methods to compute the sum of heritability explained and the corresponding SE in an R package SumVg. In conclusion, SumVg may provide a useful alternative tool for calculating SNP heritability and estimating SE/CI, which does not rely on distributional assumptions of SNP effects.

Cannabis use as a potential mediator between childhood adversity and first-episode psychosis: results from the EU-GEI case-control study.

BACKGROUND: Childhood adversity and cannabis use are considered independent risk factors for psychosis, but whether different patterns of cannabis use may be acting as mediator between adversity and psychotic disorders has not yet been explored. The aim of this study is to examine whether cannabis use mediates the relationship between childhood adversity and psychosis. METHODS: Data were utilised on 881 first-episode psychosis patients and 1231 controls from the European network of national schizophrenia networks studying Gene-Environment Interactions (EU-GEI) study. Detailed history of cannabis use was collected with the Cannabis Experience Questionnaire. The Childhood Experience of Care and Abuse Questionnaire was used to assess exposure to household discord, sexual, physical or emotional abuse and bullying in two periods: early (0-11 years), and late (12-17 years). A path decomposition method was used to analyse whether the association between childhood adversity and psychosis was mediated by (1) lifetime cannabis use, (2) cannabis potency and (3) frequency of use. RESULTS: The association between household discord and psychosis was partially mediated by lifetime use of cannabis (indirect effect coef. 0.078, s.e. 0.022, 17%), its potency (indirect effect coef. 0.059, s.e. 0.018, 14%) and by frequency (indirect effect coef. 0.117, s.e. 0.038, 29%). Similar findings were obtained when analyses were restricted to early exposure to household discord. CONCLUSIONS: Harmful patterns of cannabis use mediated the association between specific childhood adversities, like household discord, with later psychosis. Children exposed to particularly challenging environments in their household could benefit from psychosocial interventions aimed at preventing cannabis misuse.

Dynamic structure-function coupling across three major psychiatric disorders.

BACKGROUND: Convergent evidence has suggested atypical relationships between brain structure and function in major psychiatric disorders, yet how the abnormal patterns coincide and/or differ across different disorders remains largely unknown. Here, we aim to investigate the common and/or unique dynamic structure-function coupling patterns across major depressive disorder (MDD), bipolar disorder (BD), and schizophrenia (SZ). METHODS: We quantified the dynamic structure-function coupling in 452 patients with psychiatric disorders (MDD/BD/SZ = 166/168/118) and 205 unaffected controls at three distinct brain network levels, such as global, meso-, and local levels. We also correlated dynamic structure-function coupling with the topological features of functional networks to examine how the structure-function relationship facilitates brain information communication over time. RESULTS: The dynamic structure-function coupling is preserved for the three disorders at the global network level. Similar abnormalities in the rich-club organization are found in two distinct functional configuration states at the meso-level and are associated with the disease severity of MDD, BD, and SZ. At the local level, shared and unique alterations are observed in the brain regions involving the visual, cognitive control, and default mode networks. In addition, the relationships between structure-function coupling and the topological features of functional networks are altered in a manner indicative of state specificity. CONCLUSIONS: These findings suggest both transdiagnostic and illness-specific alterations in the dynamic structure-function relationship of large-scale brain networks across MDD, BD, and SZ, providing new insights and potential biomarkers into the neurodevelopmental basis underlying the behavioral and cognitive deficits observed in these disorders.

Early-onset schizophrenia is associated with immune-related rare variants in a Chinese sample.

Background: Rare variants are likely to contribute to schizophrenia (SCZ), given the large discrepancy between the heritability estimated from twin and GWAS studies. Furthermore, the nature of the rare-variant contribution to SCZ may vary with the "age-at-onset" (AAO), since early-onset has been suggested as being indicative of neurodevelopment deviance. Objective: To examine the association of rare deleterious coding variants in early- and adult-onset SCZ in a Chinese sample. Method: Exome sequencing was performed on DNA from 197 patients with SCZ spectrum disorder and 82 healthy controls (HC) of Chinese ancestry recruited in Hong Kong. We also gathered AAO information in the majority of SCZ samples. Patients were classified into early-onset (EOS, AAO<18) and adult-onset (AOS, AAO>18). We collapsed the rare variants to improve statistical power and examined the overall association of rare variants in SCZ versus HC, EOS versus HC, and AOS versus HC at the gene and gene-set levels by Sequence Kernel Association Test. The quantitative rare-variant association test of AAO was also conducted. We focused on variants which were predicted to have a medium or high impact on the protein-encoding process as defined by Ensembl. We applied a 100000-time permutation test to obtain empirical p-values, with significance threshold set at p < 1e -3 to control family-wise error rates. Moreover, we compared the burden of targeted rare variants in significant risk genes and gene sets in cases and controls. Results: Based on several binary-trait association tests (i.e., SCZ vs HC, EOS vs HC and AOS vs HC), we identified 7 candidate risk genes and 20 gene ontology biological processes (GOBP) terms, which exhibited higher burdens in SCZ than in controls. Based on quantitative rare-variant association tests, we found that alterations in 5 candidate risk genes and 7 GOBP pathways were significantly correlated with AAO. Based on biological and functional profiles of the candidate risk genes and gene sets, our findings suggested that, in addition to the involvement of perturbations in neural systems in SCZ in general, altered immune responses may be specifically implicated in EOS. Conclusion: Disrupted immune responses may exacerbate abnormal perturbations during neurodevelopment and trigger the early onset of SCZ. We provided evidence of rare variants increasing SCZ risk in the Chinese population.

Pathway-Specific Polygenic Scores Improve Cross-Ancestry Prediction of Psychosis and Clinical Outcomes.

Psychotic disorders are debilitating conditions with disproportionately high public health burden. Genetic studies indicate high heritability, but current polygenic scores (PGS) account for only a fraction of variance in psychosis risk. PGS often show poor portability across ancestries, performing significantly worse in non-European populations. Pathway-specific PGS (pPGS), which restrict PGS to genomic locations within distinct biological units, could lead to increased mechanistic understanding of pathways that lead to risk and improve cross-ancestry prediction by reducing noise in genetic predictors. This study examined the predictive power of genome-wide PGS and nine pathway-specific pPGS in a unique Chinese-ancestry sample of deeply-phenotyped psychosis patients and non-psychiatric controls. We found strong evidence for the involvement of schizophrenia-associated risk variants within "nervous system development" (p=2.5e-4) and "regulation of neuron differentiation" pathways (p=3.0e-4) in predicting risk for psychosis. We also found the "ion channel complex" pPGS, with weights derived from GWAS of bipolar disorder, to be strongly associated with the number of inpatient psychiatry admissions a patient experiences (p=1.5e-3) and account for a majority of the signal in the overall bipolar PGS. Importantly, although the schizophrenia genome-wide PGS alone explained only 3.7% of the variance in liability to psychosis in this Chinese ancestry sample, the addition of the schizophrenia-weighted pPGS for "nervous system development" and "regulation of neuron differentiation" increased the variance explained to 6.9%, which is on-par with the predictive power of PGS in European ancestry samples. Thus, not only can pPGS provide greater insight into mechanisms underlying genetic risk for disease and clinical outcomes, but may also improve cross-ancestry risk prediction accuracy.

Inferring CTCF-binding patterns and anchored loops across human tissues and cell types.

CCCTC-binding factor (CTCF) is a transcription regulator with a complex role in gene regulation. The recognition and effects of CTCF on DNA sequences, chromosome barriers, and enhancer blocking are not well understood. Existing computational tools struggle to assess the regulatory potential of CTCF-binding sites and their impact on chromatin loop formation. Here we have developed a deep-learning model, DeepAnchor, to accurately characterize CTCF binding using high-resolution genomic/epigenomic features. This has revealed distinct chromatin and sequence patterns for CTCF-mediated insulation and looping. An optimized implementation of a previous loop model based on DeepAnchor score excels in predicting CTCF-anchored loops. We have established a compendium of CTCF-anchored loops across 52 human tissue/cell types, and this suggests that genomic disruption of these loops could be a general mechanism of disease pathogenesis. These computational models and resources can help investigate how CTCF-mediated cis-regulatory elements shape context-specific gene regulation in cell development and disease progression.

Dissecting the high-resolution genetic architecture of complex phenotypes by accurately estimating gene-based conditional heritability.

Despite extensive research on global heritability estimation for complex traits, few methods accurately dissect local heritability. A precise local heritability estimate is crucial for high-resolution mapping in genetics. Here, we report the effective heritability estimator (EHE) that can use p values from genome-wide association studies (GWASs) for local heritability estimation by directly converting marginal heritability estimates of SNPs to a non-redundant heritability estimate of a gene or a small genomic region. EHE provides higher accuracy and precision for local heritability estimation among seven compared methods. Importantly, EHE can be applied to estimate the conditional heritability of nearby genes, where redundant heritability among the genes can also be removed further. The conditional estimation can be guided by tissue-specific expression profiles (or other functional scores) to prioritize and quantify more functionally important genes of complex phenotypes. Applying EHE to 42 complex phenotypes from the UK Biobank, we revealed the existence of two types of distinct genetic architectures for various complex phenotypes and found that highly pleiotropic genes are not enriched for more heritability compared to other candidate susceptibility genes. EHE provides an accurate and robust way to dissect the genetic architecture of complex phenotypes.

Neural variability in three major psychiatric disorders.

Across the major psychiatric disorders (MPDs), a shared disruption in brain physiology is suspected. Here we investigate the neural variability at rest, a well-established behavior-relevant marker of brain function, and probe its basis in gene expression and neurotransmitter receptor profiles across the MPDs. We recruited 219 healthy controls and 279 patients with schizophrenia, major depressive disorder, or bipolar disorders (manic or depressive state). The standard deviation of blood oxygenation level-dependent signal (SDBOLD) obtained from resting-state fMRI was used to characterize neural variability. Transdiagnostic disruptions in SDBOLD patterns and their relationships with clinical symptoms and cognitive functions were tested by partial least-squares correlation. Moving beyond the clinical sample, spatial correlations between the observed patterns of SDBOLD disruption and postmortem gene expressions, Neurosynth meta-analytic cognitive functions, and neurotransmitter receptor profiles were estimated. Two transdiagnostic patterns of disrupted SDBOLD were discovered. Pattern 1 is exhibited in all diagnostic groups and is most pronounced in schizophrenia, characterized by higher SDBOLD in the language/auditory networks but lower SDBOLD in the default mode/sensorimotor networks. In comparison, pattern 2 is only exhibited in unipolar and bipolar depression, characterized by higher SDBOLD in the default mode/salience networks but lower SDBOLD in the sensorimotor network. The expression of pattern 1 related to the severity of clinical symptoms and cognitive deficits across MPDs. The two disrupted patterns had distinct spatial correlations with gene expressions (e.g., neuronal projections/cellular processes), meta-analytic cognitive functions (e.g., language/memory), and neurotransmitter receptor expression profiles (e.g., D2/serotonin/opioid receptors). In conclusion, neural variability is a potential transdiagnostic biomarker of MPDs with a substantial amount of its spatial distribution explained by gene expressions and neurotransmitter receptor profiles. The pathophysiology of MPDs can be traced through the measures of neural variability at rest, with varying clinical-cognitive profiles arising from differential spatial patterns of aberrant variability.

The genetic basis of onset age in schizophrenia: evidence and models.

Schizophrenia is a heritable neurocognitive disorder affecting about 1% of the population, and usually has an onset age at around 21-25 in males and 25-30 in females. Recent advances in genetics have helped to identify many common and rare variants for the liability to schizophrenia. Earlier evidence appeared to suggest that younger onset age is associated with higher genetic liability to schizophrenia. Clinical longitudinal research also found that early and very-early onset schizophrenia are associated with poor clinical, neurocognitive, and functional profiles. A recent study reported a heritability of 0.33 for schizophrenia onset age, but the genetic basis of this trait in schizophrenia remains elusive. In the pre-Genome-Wide Association Study (GWAS) era, genetic loci found to be associated with onset age were seldom replicated. In the post-Genome-Wide Association Study era, new conceptual frameworks are needed to clarify the role of onset age in genetic research in schizophrenia, and to identify its genetic basis. In this review, we first discussed the potential of onset age as a characterizing/subtyping feature for psychosis, and as an important phenotypic dimension of schizophrenia. Second, we reviewed the methods, samples, findings and limitations of previous genetic research on onset age in schizophrenia. Third, we discussed a potential conceptual framework for studying the genetic basis of onset age, as well as the concepts of susceptibility, modifier, and "mixed" genes. Fourth, we discussed the limitations of this review. Lastly, we discussed the potential clinical implications for genetic research of onset age of schizophrenia, and how future research can unveil the potential mechanisms for this trait.