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.

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.

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.

Exploring the mediation of DNA methylation across the epigenome between childhood adversity and First Episode of Psychosis-findings from the EU-GEI study.

ABTRACT: 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. Adversity 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) 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 = 1.68; p = <0.001; abuse: OR = 2.16; p < 0.001; neglect: OR = 2.27; p = <0.001). None of the CpG sites significantly mediated the adversity-psychosis association after Bonferroni correction (p < 8.1 × 10-8). However, 28, 34 and 29 differentially methylated probes associated with 21, 27, 20 genes passed a less stringent discovery threshold (p < 5 × 10-5) for composite, abuse and neglect respectively, with a lack of overlap between abuse and neglect. These included genes previously associated to psychosis in EWAS studies, such as PANK1, SPEG TBKBP1, TSNARE1 or H2R. Downstream gene ontology analyses did not reveal any biological pathways that survived false discovery rate correction. Although at a non-significant level, DNAm changes in genes previously associated with schizophrenia in EWAS studies may mediate the CA-psychosis association. These results and associated involved processes such as mitochondrial or histaminergic disfunction, immunity or neural signalling requires replication in well powered samples. The lack of overlap between mediating genes associated with abuse and neglect suggests differential biological trajectories linking CA subtypes and psychosis.

Clarifying the causes of consistent and inconsistent findings in genetics.

As research in genetics has advanced, some findings have been unexpected or shown to be inconsistent between studies or datasets. The reasons these inconsistencies arise are complex. Results from genetic studies can be affected by various factors including statistical power, linkage disequilibrium, quality control, confounding and selection bias, as well as real differences from interactions and effect modifiers, which may be informative about the mechanisms of traits and disease. Statistical artefacts can manifest as differences between results but they can also conceal underlying differences, which implies that their critical examination is important for understanding the underpinnings of traits. In this review, we examine these factors and outline how they can be identified and conceptualised with structural causal models. We explain the consequences they have on genetic estimates, such as genetic associations, polygenic scores, family- and genome-wide heritability, and describe methods to address them to aid in the estimation of true effects of genetic variation. Clarifying these factors can help researchers anticipate when results are likely to diverge and aid researchers' understanding of causal relationships between genes and complex traits.

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.

Use of multiple polygenic risk scores for distinguishing schizophrenia-spectrum disorder and affective psychosis categories in a first-episode sample; the EU-GEI study.

BACKGROUND: Schizophrenia (SZ), bipolar disorder (BD) and depression (D) run in families. This susceptibility is partly due to hundreds or thousands of common genetic variants, each conferring a fractional risk. The cumulative effects of the associated variants can be summarised as a polygenic risk score (PRS). Using data from the EUropean Network of national schizophrenia networks studying Gene-Environment Interactions (EU-GEI) first episode case-control study, we aimed to test whether PRSs for three major psychiatric disorders (SZ, BD, D) and for intelligent quotient (IQ) as a neurodevelopmental proxy, can discriminate affective psychosis (AP) from schizophrenia-spectrum disorder (SSD). METHODS: Participants (842 cases, 1284 controls) from 16 European EU-GEI sites were successfully genotyped following standard quality control procedures. The sample was stratified based on genomic ancestry and analyses were done only on the subsample representing the European population (573 cases, 1005 controls). Using PRS for SZ, BD, D, and IQ built from the latest available summary statistics, we performed simple or multinomial logistic regression models adjusted for 10 principal components for the different clinical comparisons. RESULTS: In case-control comparisons PRS-SZ, PRS-BD and PRS-D distributed differentially across psychotic subcategories. In case-case comparisons, both PRS-SZ [odds ratio (OR) = 0.7, 95% confidence interval (CI) 0.54-0.92] and PRS-D (OR = 1.31, 95% CI 1.06-1.61) differentiated AP from SSD; and within AP categories, only PRS-SZ differentiated BD from psychotic depression (OR = 2.14, 95% CI 1.23-3.74). CONCLUSIONS: Combining PRS for severe psychiatric disorders in prediction models for psychosis phenotypes can increase discriminative ability and improve our understanding of these phenotypes. Our results point towards the potential usefulness of PRSs in specific populations such as high-risk or early psychosis phases.

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.

The association between reasons for first using cannabis, later pattern of use, and risk of first-episode psychosis: the EU-GEI case-control study.

BACKGROUND: While cannabis use is a well-established risk factor for psychosis, little is known about any association between reasons for first using cannabis (RFUC) and later patterns of use and risk of psychosis. METHODS: We used data from 11 sites of the multicentre European Gene-Environment Interaction (EU-GEI) case-control study. 558 first-episode psychosis patients (FEPp) and 567 population controls who had used cannabis and reported their RFUC.We ran logistic regressions to examine whether RFUC were associated with first-episode psychosis (FEP) case-control status. Path analysis then examined the relationship between RFUC, subsequent patterns of cannabis use, and case-control status. RESULTS: Controls (86.1%) and FEPp (75.63%) were most likely to report 'because of friends' as their most common RFUC. However, 20.1% of FEPp compared to 5.8% of controls reported: 'to feel better' as their RFUC (χ2 = 50.97; p < 0.001). RFUC 'to feel better' was associated with being a FEPp (OR 1.74; 95% CI 1.03-2.95) while RFUC 'with friends' was associated with being a control (OR 0.56; 95% CI 0.37-0.83). The path model indicated an association between RFUC 'to feel better' with heavy cannabis use and with FEPp-control status. CONCLUSIONS: Both FEPp and controls usually started using cannabis with their friends, but more patients than controls had begun to use 'to feel better'. People who reported their reason for first using cannabis to 'feel better' were more likely to progress to heavy use and develop a psychotic disorder than those reporting 'because of friends'.

Estimating indirect parental genetic effects on offspring phenotypes using virtual parental genotypes derived from sibling and half sibling pairs.

Indirect parental genetic effects may be defined as the influence of parental genotypes on offspring phenotypes over and above that which results from the transmission of genes from parents to their children. However, given the relative paucity of large-scale family-based cohorts around the world, it is difficult to demonstrate parental genetic effects on human traits, particularly at individual loci. In this manuscript, we illustrate how parental genetic effects on offspring phenotypes, including late onset conditions, can be estimated at individual loci in principle using large-scale genome-wide association study (GWAS) data, even in the absence of parental genotypes. Our strategy involves creating "virtual" mothers and fathers by estimating the genotypic dosages of parental genotypes using physically genotyped data from relative pairs. We then utilize the expected dosages of the parents, and the actual genotypes of the offspring relative pairs, to perform conditional genetic association analyses to obtain asymptotically unbiased estimates of maternal, paternal and offspring genetic effects. We apply our approach to 19066 sibling pairs from the UK Biobank and show that a polygenic score consisting of imputed parental educational attainment SNP dosages is strongly related to offspring educational attainment even after correcting for offspring genotype at the same loci. We develop a freely available web application that quantifies the power of our approach using closed form asymptotic solutions. We implement our methods in a user-friendly software package IMPISH (IMputing Parental genotypes In Siblings and Half Siblings) which allows users to quickly and efficiently impute parental genotypes across the genome in large genome-wide datasets, and then use these estimated dosages in downstream linear mixed model association analyses. We conclude that imputing parental genotypes from relative pairs may provide a useful adjunct to existing large-scale genetic studies of parents and their offspring.

Potential role of regulatory DNA variants in modifying the risk of severe cutaneous reactions induced by aromatic anti-seizure medications.

OBJECTIVE: Stevens-Johnson syndrome (SJS) and toxic epidermal necrolysis (TEN) are severe cutaneous adverse drug reactions. Antiseizure medications (ASMs) with aromatic ring structure, including carbamazepine, are among the most common culprits. Screening for human leukocyte antigen (HLA) allele HLA-B*15:02 is recommended prior to initiating treatment with carbamazepine in Asians, but this allele has low positive predictive value. METHODS: We performed whole genome sequencing and analyzed 6 199 696 common variants among 113 aromatic ASM-induced SJS/TEN cases and 84 tolerant controls of Han Chinese ethnicity. RESULTS: In the primary analysis, nine variants reached genome-wide significance (p < 5e-08), one in the carbamazepine subanalysis (85 cases vs. 77 controls) and a further eight identified in HLA-B*15:02-negative subanalysis (35 cases and 53 controls). Interaction analysis between each novel variant from the primary analysis found that five increased risk irrespective of HLA-B*15:02 status or zygosity. HLA-B*15:02-positive individuals were found to have reduced risk if they also carried a chromosome 12 variant, chr12.9426934 (heterozygotes: relative risk = .71, p = .001; homozygotes: relative risk = .23, p < .001). All significant variants lie within intronic or intergenic regions with poorly understood functional consequence. In silico functional analysis of suggestive variants (p < 5e-6) identified through the primary and subanalyses (stratified by HLA-B*15:02 status and drug exposure) suggests that genetic variation within regulatory DNA may contribute to risk indirectly by disrupting the regulation of pathology-related genes. The genes implicated were specific either to the primary analysis (CD9), HLA-B*15:02 carriers (DOCK10), noncarriers (ABCA1), carbamazepine exposure (HLA-E), or phenytoin exposure (CD24). SIGNIFICANCE: We identified variants that could explain why some carriers of HLA-B*15:02 tolerate treatment, and why some noncarriers develop ASM-induced SJS/TEN. Additionally, this analysis suggests that the mixing of HLA-B*15:02 carrier status in previous studies might have masked variants contributing to susceptibility, and that inheritance of risk for ASM-induced SJS/TEN is complex, likely involving multiple risk variants.

Improved nonparametric penalized maximum likelihood estimation for arbitrarily censored survival data.

Nonparametric maximum likelihood estimation encompasses a group of classic methods to estimate distribution-associated functions from potentially censored and truncated data, with extensive applications in survival analysis. These methods, including the Kaplan-Meier estimator and Turnbull's method, often result in overfitting, especially when the sample size is small. We propose an improvement to these methods by applying kernel smoothing to their raw estimates, based on a BIC-type loss function that balances the trade-off between optimizing model fit and controlling model complexity. In the context of a longitudinal study with repeated observations, we detail our proposed smoothing procedure and optimization algorithm. With extensive simulation studies over multiple realistic scenarios, we demonstrate that our smoothing-based procedure provides better overall accuracy in both survival function estimation and individual-level time-to-event prediction (imputation) by reducing overfitting. Our smoothing procedure decreases the bias (discrepancy between the estimated and true simulated survival function) using interval-censored data by up to 48% compared to the raw un-smoothed estimate, with similar improvements of up to 34% and 23% in within-sample and out-of-sample prediction, respectively. Our smoothing algorithm also demonstrates significant overall improvement across all three metrics when compared to a popular semiparametric B-splines estimation method. Finally, we apply our method to real data on censored breast cancer diagnosis, which similarly shows improvement when compared to empirical survival estimates from uncensored data. We provide an R package, SISE, for implementing our penalized likelihood method.

Integrative analysis of metabolomic, genomic, and imaging-based phenotypes identify very-low-density lipoprotein as a potential risk factor for lumbar Modic changes.

PURPOSE: Modic changes (MC) on magnetic resonance imaging (MRI) have been associated with the development and severity of low back pain (LBP). The etiology of MC remains elusive, but it has been suggested that altered metabolism may be a risk factor. As such, this study aimed to identify metabolomic biomarkers for MC phenotypes of the lumbar spine via a combined metabolomic-genomic approach. METHODS: A population cohort of 3,584 southern Chinese underwent lumbar spine MRI. Blood samples were genotyped with single-nucleotide polymorphisms (SNP) arrays (n = 2,482) and serum metabolomics profiling using magnetic resonance spectroscopy (n = 757), covering 130 metabolites representing three molecular windows, were assessed. Genome-wide association studies (GWAS) were performed on each metabolite, to construct polygenic scores for predicting metabolite levels in subjects who had GWAS but not metabolomic data. Associations between predicted metabolite levels and MC phenotypes were assessed using linear/logistic regression and least absolute shrinkage and selection operator (LASSO). Two-sample Mendelian randomization analysis tested for causal relationships between metabolic biomarkers and MC. RESULTS: 20.4% had MC (10.6% type 1, 67.2% type 2, 22.2% mixed types). Significant MC metabolomic biomarkers were mean diameter of very-low-density lipoprotein (VLDL)/low-density lipoprotein (LDL) particles and cholesterol esters/phospholipids in large LDL. Mendelian randomization indicated that decreased VLDL mean diameter may lead to MC. CONCLUSIONS: This large-scale study is the first to address metabolomics in subject with/without lumbar MC. Causality studies implicate VLDL related to MC, noting a metabolic etiology. Our study substantiates the field of "spino-metabolomics" and illustrates the power of integrating metabolomics-genomics-imaging phenotypes to discover biomarkers for spinal disorders, paving the way for more personalized spine care for patients.

Intermediate confounding in trio relationships: The importance of complete data in effect size estimation.

We present an important characteristic of trio models which may lead to bias and loss of power when one parent is unmodeled in trio analyses. Motivated by recent interest in estimating parental effects on postnatal and later-life phenotypes, we consider a causal model where each parent has both an effect on their child's phenotype which is mediated through the genotype transmitted to the child and a direct effect on the phenotype through the parentally provided environment. We derive the power and bias of models in which one parent's genotype is not modeled, showing that while the effect of the child's genotype is biased in the direction of the unmodeled parent's effect as expected, the estimated effect of the observed parent's genotype is also biased in the opposite direction. While this phenomenon may not be intuitive under the assumption of random mating, it can be explained by intermediate confounding of the child's genotype-phenotype effect. These observations have implications for the accurate estimation of maternal and paternal effects in trio data sets with missing genotype data.

Causal relationships between blood lipids and depression phenotypes: a Mendelian randomisation analysis.

BACKGROUND: The etiology of depression remains poorly understood. Changes in blood lipid levels were reported to be associated with depression and suicide, however study findings were mixed. METHODS: We performed a two-sample Mendelian randomisation (MR) analysis to investigate the causal relationship between blood lipids and depression phenotypes, based on large-scale GWAS summary statistics (N = 188 577/480 359 for lipid/depression traits respectively). Five depression-related phenotypes were included, namely major depression (MD; from PGC), depressive symptoms (DS; from SSGAC), longest duration and number of episodes of low mood, and history of deliberate self-harm (DSH)/suicide (from UK Biobank). MR was conducted with inverse-variance weighted (MR-IVW), Egger and Generalised Summary-data-based MR (GSMR) methods. RESULTS: There was consistent evidence that triglyceride (TG) is causally associated with DS (MR-IVW ß for one-s.d. increase in TG = 0.0346, 95% CI 0.0114-0.0578), supported by MR-IVW and GSMR and multiple r2 clumping thresholds. We also observed relatively consistent associations of TG with DSH/suicide (MR-Egger OR = 2.514, CI 1.579-4.003). There was moderate evidence for positive associations of TG with MD and the number of episodes of low mood. For HDL-c, we observed moderate evidence for causal associations with DS and MD. LDL-c and TC did not show robust causal relationships with depression phenotypes, except for weak evidence that LDL-c is inversely related to DSH/suicide. We did not detect significant associations when depression phenotypes were treated as exposures. CONCLUSIONS: This study provides evidence to a causal relationship between TG, and to a lesser extent, altered cholesterol levels with depression phenotypes. Further studies on its mechanistic basis and the effects of lipid-lowering therapies are warranted.

Jumping to conclusions, general intelligence, and psychosis liability: findings from the multi-centre EU-GEI case-control study.

BACKGROUND: The 'jumping to conclusions' (JTC) bias is associated with both psychosis and general cognition but their relationship is unclear. In this study, we set out to clarify the relationship between the JTC bias, IQ, psychosis and polygenic liability to schizophrenia and IQ. METHODS: A total of 817 first episode psychosis patients and 1294 population-based controls completed assessments of general intelligence (IQ), and JTC, and provided blood or saliva samples from which we extracted DNA and computed polygenic risk scores for IQ and schizophrenia. RESULTS: The estimated proportion of the total effect of case/control differences on JTC mediated by IQ was 79%. Schizophrenia polygenic risk score was non-significantly associated with a higher number of beads drawn (B = 0.47, 95% CI -0.21 to 1.16, p = 0.17); whereas IQ PRS (B = 0.51, 95% CI 0.25-0.76, p < 0.001) significantly predicted the number of beads drawn, and was thus associated with reduced JTC bias. The JTC was more strongly associated with the higher level of psychotic-like experiences (PLEs) in controls, including after controlling for IQ (B = -1.7, 95% CI -2.8 to -0.5, p = 0.006), but did not relate to delusions in patients. CONCLUSIONS: Our findings suggest that the JTC reasoning bias in psychosis might not be a specific cognitive deficit but rather a manifestation or consequence, of general cognitive impairment. Whereas, in the general population, the JTC bias is related to PLEs, independent of IQ. The work has the potential to inform interventions targeting cognitive biases in early psychosis.

Daily use of high-potency cannabis is associated with more positive symptoms in first-episode psychosis patients: the EU-GEI case-control study.

BACKGROUND: Daily use of high-potency cannabis has been reported to carry a high risk for developing a psychotic disorder. However, the evidence is mixed on whether any pattern of cannabis use is associated with a particular symptomatology in first-episode psychosis (FEP) patients. METHOD: We analysed data from 901 FEP patients and 1235 controls recruited across six countries, as part of the European Network of National Schizophrenia Networks Studying Gene-Environment Interactions (EU-GEI) study. We used item response modelling to estimate two bifactor models, which included general and specific dimensions of psychotic symptoms in patients and psychotic experiences in controls. The associations between these dimensions and cannabis use were evaluated using linear mixed-effects models analyses. RESULTS: In patients, there was a linear relationship between the positive symptom dimension and the extent of lifetime exposure to cannabis, with daily users of high-potency cannabis having the highest score (B = 0.35; 95% CI 0.14-0.56). Moreover, negative symptoms were more common among patients who never used cannabis compared with those with any pattern of use (B = -0.22; 95% CI -0.37 to -0.07). In controls, psychotic experiences were associated with current use of cannabis but not with the extent of lifetime use. Neither patients nor controls presented differences in depressive dimension related to cannabis use. CONCLUSIONS: Our findings provide the first large-scale evidence that FEP patients with a history of daily use of high-potency cannabis present with more positive and less negative symptoms, compared with those who never used cannabis or used low-potency types.

Modeling Parent-Specific Genetic Nurture in Families with Missing Parental Genotypes: Application to Birthweight and BMI.

Disaggregation and estimation of genetic effects from offspring and parents has long been of interest to statistical geneticists. Recently, technical and methodological advances have made the genome-wide and loci-specific estimation of direct offspring and parental genetic nurture effects more possible. However, unbiased estimation using these methods requires datasets where both parents and at least one child have been genotyped, which are relatively scarce. Our group has recently developed a method and accompanying software (IMPISH; Hwang et al. in PLoS Genet 16:e1009154, 2020) which is able to impute missing parental genotypes from observed data on sibships and estimate their effects on an offspring phenotype conditional on the effects of genetic transmission. However, this method is unable to disentangle maternal and paternal effects, which may differ in magnitude and direction. Here, we introduce an extension to the original IMPISH routine which takes advantage of all available nuclear families to impute parent-specific missing genotypes and obtain asymptotically unbiased estimates of genetic effects on offspring phenotypes. We apply this this method to data from related individuals in the UK Biobank, showing concordance with previous estimates of maternal genetic effects on offspring birthweight. We also conduct the first GWAS jointly estimating offspring-, maternal-, and paternal-specific genetic effects on body-mass index.

WITER: a powerful method for estimation of cancer-driver genes using a weighted iterative regression modelling background mutation counts.

Genomic identification of driver mutations and genes in cancer cells are critical for precision medicine. Due to difficulty in modelling distribution of background mutation counts, existing statistical methods are often underpowered to discriminate cancer-driver genes from passenger genes. Here we propose a novel statistical approach, weighted iterative zero-truncated negative-binomial regression (WITER, http://grass.cgs.hku.hk/limx/witer or KGGSeq,http://grass.cgs.hku.hk/limx/kggseq/), to detect cancer-driver genes showing an excess of somatic mutations. By fitting the distribution of background mutation counts properly, this approach works well even in small or moderate samples. Compared to alternative methods, it detected more significant and cancer-consensus genes in most tested cancers. Applying this approach, we estimated 229 driver genes in 26 different types of cancers. In silico validation confirmed 78% of predicted genes as likely known drivers and many other genes as very likely new drivers for corresponding cancers. The technical advances of WITER enable the detection of driver genes in TCGA datasets as small as 30 subjects and rescue of more genes missed by alternative tools in moderate or small samples.

De novo mutations in Caudal Type Homeo Box transcription Factor 2 (CDX2) in patients with persistent cloaca.

The cloaca is an embryonic cavity that is divided into the urogenital sinus and rectum upon differentiation of the cloacal epithelium triggered by tissue-specific transcription factors including CDX2. Defective differentiation leads to persistent cloaca in humans (PC), a phenotype recapitulated in Cdx2 mutant mice. PC is linked to hypo/hyper-vitaminosis A. Although no gene has ever been identified, there is a strong evidence for a genetic contribution to PC. We applied whole-exome sequencing and copy-number-variants analyses to 21 PC patients and their unaffected parents. The damaging p.Cys132* and p.Arg237His de novo CDX2 variants were identified in two patients. These variants altered the expression of CYP26A1, a direct CDX2 target encoding the major retinoic acid (RA)-degrading enzyme. Other RA genes, including the RA-receptor alpha, were also mutated. Genes governing the development of cloaca-derived structures were recurrently mutated and over-represented in the basement-membrane components set (q-value < 1.65 × 10-6). Joint analysis of the patients' profile highlighted the extracellular matrix-receptor interaction pathway (MsigDBID: M7098, FDR: q-value < 7.16 × 10-9). This is the first evidence that PC is genetic, with genes involved in the RA metabolism at the lead. Given the CDX2 de novo variants and the role of RA, our observations could potentiate preventive measures. For the first time, a gene recapitulating PC in mouse models is found mutated in humans.