Sex-Dependent Shared and Nonshared Genetic Architecture Across Mood and Psychotic Disorders

BACKGROUND: Sex differences in incidence and/or presentation of schizophrenia (SCZ), major depressive disorder (MDD), and bipolar disorder (BIP) are pervasive. Previous evidence for shared genetic risk and sex differences in brain abnormalities across disorders suggest possible shared sex-dependent genetic risk. METHODS: We conducted the largest to date genome-wide genotype-by-sex (GxS) interaction of risk for these disorders using 85,735 cases (33,403 SCZ, 19,924 BIP, and 32,408 MDD) and 109,946 controls from the PGC (Psychiatric Genomics Consortium) and iPSYCH. RESULTS: Across disorders, genome-wide significant single nucleotide polymorphism-by-sex interaction was detected for a locus encompassing NKAIN2 (rs117780815, p = 3.2 x 10(-8)), which interacts with sodium/potassium-transporting ATPase (adenosine triphosphatase) enzymes, implicating neuronal excitability. Three additional loci showed evidence (p < 1 x 10(-8)) for cross-disorder GxS interaction (rs7302529, p = 1.6 x 10(-7); rs73033497, p = 8.8 x 10(-7); rs7914279, p = 6.4 x 10(-7)), implicating various functions. Gene-based analyses identified GxS interaction across disorders (p = 8.97 x 10(-7)) with transcriptional inhibitor SLTM. Most significant in SCZ was a MOCOS gene locus (rs11665282, p = 1.5 x 10(-7)), implicating vascular endothelial cells. Secondary analysis of the PGC-SCZ dataset detected an interaction (rs13265509, p = 1.1 x 10(-7)) in a locus containing IDO2, a kynurenine pathway enzyme with immunoregulatory functions implicated in SCZ, BIP, and MDD. Pathway enrichment analysis detected significant GxS interaction of genes regulating vascular endothelial growth factor receptor signaling in MDD (false discovery rate-corrected p < .05). CONCLUSIONS: In the largest genome-wide GXS analysis of mood and psychotic disorders to date, there was substantial genetic overlap between the sexes. However, significant sex-dependent effects were enriched for genes related to neuronal development and immune and vascular functions across and within SCZ, BIP, and MDD at the variant, gene, and pathway levels.Peer reviewe

Genetic correlations of psychiatric traits with body composition and glycemic traits are sex- and age-dependent

Body composition is often altered in psychiatric disorders. Using genome-wide common genetic variation data, we calculate sex-specific genetic correlations amongst body fat %, fat mass, fat-free mass, physical activity, glycemic traits and 17 psychiatric traits (up to N = 217,568). Two patterns emerge: (1) anorexia nervosa, schizophrenia, obsessive-compulsive disorder, and education years are negatively genetically correlated with body fat % and fat-free mass, whereas (2) attention-deficit/hyperactivity disorder (ADHD), alcohol dependence, insomnia, and heavy smoking are positively correlated. Anorexia nervosa shows a stronger genetic correlation with body fat % in females, whereas education years is more strongly correlated with fat mass in males. Education years and ADHD show genetic overlap with childhood obesity. Mendelian randomization identifies schizophrenia, anorexia nervosa, and higher education as causal for decreased fat mass, with higher body fat % possibly being a causal risk factor for ADHD and heavy smoking. These results suggest new possibilities for targeted preventive strategies

A recessive genetic model and runs of homozygosity in major depressive disorder

Genome-wide association studies (GWASs) of major depressive disorder (MDD) have yet to identify variants that surpass the threshold for genome-wide significance. A recent study reported that runs of homozygosity (ROH) are associated with schizophrenia, reflecting a novel genetic risk factor resulting from increased parental relatedness and recessive genetic effects. Here, we explore the possibility of such a recessive model in MDD. In a sample of 9,238 cases and 9,521 controls reported in a recent mega-analysis of 9 GWAS we perform an analysis of ROH and common variants under a recessive model. Since evidence for association with ROH could reflect a recessive mode of action at loci, we also conducted a genome-wide association analyses under a recessive model. The genome-wide association analysis using a recessive model found no significant associations. Our analysis of ROH suggested that there was significant heterogeneity of effect across studies in effect (P=0.001), and it was associated with genotyping platform and country of origin. The results of the ROH analysis show that differences across studies can lead to conflicting systematic genome-wide differences between cases and controls that are unaccounted for by traditional covariates. They highlight the sensitivity of the ROH method to spurious associations, and the need to carefully control for potential confounds in such analyses. We found no strong evidence for a recessive model underlying MDD

Familial aggregation and heritability of schizophrenia and co-aggregation of psychiatric illnesses in affected families

Strong familial aggregation of schizophrenia has been reported but there is uncertainty concerning the degree of genetic contribution to the phenotypic variance of the disease. This study aimed to examine the familial aggregation and heritability of schizophrenia, and the relative risks (RRs) of other psychiatric diseases, in relatives of people with schizophrenia using the Taiwan National Health Insurance Database. The study population included individuals with affected first-degree or second-degree relatives identified from all beneficiaries (n = 23 422 955) registered in 2013. Diagnoses of schizophrenia made by psychiatrists were ascertained between January 1, 1996 and December 31, 2013. Having an affected co-twin, first-degree relative, second-degree relative, or spouse was associated with an adjusted RR (95% CI) of 37.86 (30.55-46.92), 6.30 (6.09-6.53), 2.44 (1.91-3.12), and 1.88 (1.64-2.15), respectively. Compared with the general population, individuals with one affected first-degree relative had a RR (95% CI) of 6.00 (5.79-6.22) and those with 2 or more had a RR (95% CI) of 14.66 (13.00-16.53) for schizophrenia. The accountability for the phenotypic variance of schizophrenia was 47.3% for genetic factors, 15.5% for shared environmental factors, and 37.2% for non-shared environmental factors. The RR (95% CI) in individuals with a first-degree relative with schizophrenia was 3.49 (3.34-3.64) for mood disorders and 3.91 (3.35-4.57) for delusional disorders. A family history of schizophrenia is therefore associated with a higher risk of developing schizophrenia, mood disorders, and delusional disorders. Heritability and environmental factors each account for half of the phenotypic variance of schizophrenia

Phenome-wide analysis of genome-wide polygenic scores

Genome-wide polygenic scores (GPS), which aggregate the effects of thousands of DNA variants from genome-wide association studies (GWAS), have the potential to make genetic predictions for individuals. We conducted a systematic investigation of associations between GPS and many behavioral traits, the behavioral phenome. For 3152 unrelated 16-year-old individuals representative of the United Kingdom, we created 13 GPS from the largest GWAS for psychiatric disorders (for example, schizophrenia, depression and dementia) and cognitive traits (for example, intelligence, educational attainment and intracranial volume). The behavioral phenome included 50 traits from the domains of psychopathology, personality, cognitive abilities and educational achievement. We examined phenome-wide profiles of associations for the entire distribution of each GPS and for the extremes of the GPS distributions. The cognitive GPS yielded stronger predictive power than the psychiatric GPS in our UK-representative sample of adolescents. For example, education GPS explained variation in adolescents’ behavior problems (~0.6%) and in educational achievement (~2%) but psychiatric GPS were associated with neither. Despite the modest effect sizes of current GPS, quantile analyses illustrate the ability to stratify individuals by GPS and opportunities for research. For example, the highest and lowest septiles for the education GPS yielded a 0.5 s.d. difference in mean math grade and a 0.25 s.d. difference in mean behavior problems. We discuss the usefulness and limitations of GPS based on adult GWAS to predict genetic propensities earlier in development

Evidence For Genetic Heterogeneity Between Clinical Subtypes of Bipolar Disorder

We performed a genome-wide association study of 6447 bipolar disorder (BD) cases and 12 639 controls from the International Cohort Collection for Bipolar Disorder (ICCBD). Meta-analysis was performed with prior results from the Psychiatric Genomics Consortium Bipolar Group for a combined sample of 13 902 cases and 19 279 controls. We identified eight genome-wide significant, associated regions, including a novel associated region on chromosome 10 (rs10884920; P = 3.28 × 10 − 8) that includes the brain-enriched cytoskeleton protein adducin 3 (ADD3), a non-coding RNA, and a neuropeptide-specific aminopeptidase P (XPNPEP1). Our large sample size allowed us to test the heritability and genetic correlation of BD subtypes and investigate their genetic overlap with schizophrenia (SCZ) and major depressive disorder. We found a significant difference in heritability of the two most common forms of BD (BD I h2 = 0.35; BD II h2 = 0.25; P = 0.02) with a genetic correlation between BD I and BD II of 0.78,compared with a genetic correlation of 0.97 when BD cohorts containing both types were compared. In addition, we demonstrated a significantly greater load of polygenic risk alleles for SCZ and BD in patients with BD I compared with patients with BD II, and a greater load of SCZ risk alleles in the bipolar type of schizoaffective disorder (SAB) compared with both other BD subtypes. These results point to a partial difference in genetic architecture of BD subtypes, and are suggestive of a molecular correlate for the clinical division of BD into subtypes

The polygenic risk for bipolar disorder influences brain regional function relating to visual and default state processing of emotional information

Genome-wise association studies have identified a number of common single-nucleotide polymorphisms (SNPs), each of small effect, associated with risk to bipolar disorder (BD). Several risk-conferring SNPs have been individually shown to influence regional brain activation thus linking genetic risk for BD to altered brain function. The current study examined whether the polygenic risk score method, which models the cumulative load of all known risk-conferring SNPs, may be useful in the identification of brain regions whose function may be related to the polygenic architecture of BD. We calculated the individual polygenic risk score for BD (PGR-BD) in forty-one patients with the disorder, twenty-five unaffected first-degree relatives and forty-six unrelated healthy controls using the most recent Psychiatric Genomics Consortium data. Functional magnetic resonance imaging was used to define task-related brain activation patterns in response to facial affect and working memory processing. We found significant effects of the PGR-BD score on task-related activation irrespective of diagnostic group. There was a negative association between the PGR-BD score and activation in the visual association cortex during facial affect processing. In contrast, the PGR-BD score was associated with failure to deactivate the ventromedial prefrontal region of the default mode network during working memory processing. These results are consistent with the threshold-liability model of BD, and demonstrate the usefulness of the PGR-BD score in identifying brain functional alternations associated with vulnerability to BD. Additionally, our findings suggest that the polygenic architecture of BD is not regionally confined but impacts on the task-dependent recruitment of multiple brain regions

Profiling allele-specific gene expression in brains from individuals with autism spectrum disorder reveals preferential minor allele usage.

One fundamental but understudied mechanism of gene regulation in disease is allele-specific expression (ASE), the preferential expression of one allele. We leveraged RNA-sequencing data from human brain to assess ASE in autism spectrum disorder (ASD). When ASE is observed in ASD, the allele with lower population frequency (minor allele) is preferentially more highly expressed than the major allele, opposite to the canonical pattern. Importantly, genes showing ASE in ASD are enriched in those downregulated in ASD postmortem brains and in genes harboring de novo mutations in ASD. Two regions, 14q32 and 15q11, containing all known orphan C/D box small nucleolar RNAs (snoRNAs), are particularly enriched in shifts to higher minor allele expression. We demonstrate that this allele shifting enhances snoRNA-targeted splicing changes in ASD-related target genes in idiopathic ASD and 15q11-q13 duplication syndrome. Together, these results implicate allelic imbalance and dysregulation of orphan C/D box snoRNAs in ASD pathogenesis

Copy number variation in bipolar disorder.

Large (>100 kb), rare (500 kb) CNVs in BD compared with SZ, most notably for deletions >1 Mb (P=9 × 10(-4))