A catalog of associations between rare coding variants and COVID-19 outcomes.

Severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) causes coronavirus disease-19 (COVID-19), a respiratory illness that can result in hospitalization or death. We investigated associations between rare genetic variants and seven COVID-19 outcomes in 543,213 individuals, including 8,248 with COVID-19. After accounting for multiple testing, we did not identify any clear associations with rare variants either exome-wide or when specifically focusing on (i) 14 interferon pathway genes in which rare deleterious variants have been reported in severe COVID-19 patients; (ii) 167 genes located in COVID-19 GWAS risk loci; or (iii) 32 additional genes of immunologic relevance and/or therapeutic potential. Our analyses indicate there are no significant associations with rare protein-coding variants with detectable effect sizes at our current sample sizes. Analyses will be updated as additional data become available, with results publicly browsable at https://rgc-covid19.regeneron.com.

Genome-wide analysis in 756,646 individuals provides first genetic evidence that ACE2 expression influences COVID-19 risk and yields genetic risk scores predictive of severe disease.

SARS-CoV-2 enters host cells by binding angiotensin-converting enzyme 2 (ACE2). Through a genome-wide association study, we show that a rare variant (MAF = 0.3%, odds ratio 0.60, P=4.5×10-13) that down-regulates ACE2 expression reduces risk of COVID-19 disease, providing human genetics support for the hypothesis that ACE2 levels influence COVID-19 risk. Further, we show that common genetic variants define a risk score that predicts severe disease among COVID-19 cases.

Characterization of bipolar disorder patient-specific induced pluripotent stem cells from a family reveals neurodevelopmental and mRNA expression abnormalities.

Bipolar disorder (BD) is a common neuropsychiatric disorder characterized by chronic recurrent episodes of depression and mania. Despite evidence for high heritability of BD, little is known about its underlying pathophysiology. To develop new tools for investigating the molecular and cellular basis of BD, we applied a family-based paradigm to derive and characterize a set of 12 induced pluripotent stem cell (iPSC) lines from a quartet consisting of two BD-affected brothers and their two unaffected parents. Initially, no significant phenotypic differences were observed between iPSCs derived from the different family members. However, upon directed neural differentiation, we observed that CXCR4 (CXC chemokine receptor-4) expressing central nervous system (CNS) neural progenitor cells (NPCs) from both BD patients compared with their unaffected parents exhibited multiple phenotypic differences at the level of neurogenesis and expression of genes critical for neuroplasticity, including WNT pathway components and ion channel subunits. Treatment of the CXCR4(+) NPCs with a pharmacological inhibitor of glycogen synthase kinase 3, a known regulator of WNT signaling, was found to rescue a progenitor proliferation deficit in the BD patient NPCs. Taken together, these studies provide new cellular tools for dissecting the pathophysiology of BD and evidence for dysregulation of key pathways involved in neurodevelopment and neuroplasticity. Future generation of additional iPSCs following a family-based paradigm for modeling complex neuropsychiatric disorders in conjunction with in-depth phenotyping holds promise for providing insights into the pathophysiological substrates of BD and is likely to inform the development of targeted therapeutics for its treatment and ideally prevention.

Copy number variation in schizophrenia in Sweden.

Schizophrenia (SCZ) is a highly heritable neuropsychiatric disorder of complex genetic etiology. Previous genome-wide surveys have revealed a greater burden of large, rare copy number variations (CNVs) in SCZ cases and identified multiple rare recurrent CNVs that increase risk of SCZ although with incomplete penetrance and pleiotropic effects. Identification of additional recurrent CNVs and biological pathways enriched for SCZ CNVs requires greater sample sizes. We conducted a genome-wide survey for CNVs associated with SCZ using a Swedish national sample (4719 cases and 5917 controls). High-confidence CNV calls were generated using genotyping array intensity data, and their effect on risk of SCZ was measured. Our data confirm increased burden of large, rare CNVs in SCZ cases as well as significant associations for recurrent 16p11.2 duplications, 22q11.2 deletions and 3q29 deletions. We report a novel association for 17q12 duplications (odds ratio=4.16, P=0.018), previously associated with autism and mental retardation but not SCZ. Intriguingly, gene set association analyses implicate biological pathways previously associated with SCZ through common variation and exome sequencing (calcium channel signaling and binding partners of the fragile X mental retardation protein). We found significantly increased burden of the largest CNVs (>500 kb) in genes present in the postsynaptic density, in genomic regions implicated via SCZ genome-wide association studies and in gene products localized to mitochondria and cytoplasm. Our findings suggest that multiple lines of genomic inquiry--genome-wide screens for CNVs, common variation and exonic variation--are converging on similar sets of pathways and/or genes.

Evidence that duplications of 22q11.2 protect against schizophrenia.

A number of large, rare copy number variants (CNVs) are deleterious for neurodevelopmental disorders, but large, rare, protective CNVs have not been reported for such phenotypes. Here we show in a CNV analysis of 47 005 individuals, the largest CNV analysis of schizophrenia to date, that large duplications (1.5-3.0 Mb) at 22q11.2--the reciprocal of the well-known, risk-inducing deletion of this locus--are substantially less common in schizophrenia cases than in the general population (0.014% vs 0.085%, OR=0.17, P=0.00086). 22q11.2 duplications represent the first putative protective mutation for schizophrenia.

The one and the many: effects of the cell adhesion molecule pathway on neuropsychological function in psychosis.

BACKGROUND: Genetic studies of single gene variants have been criticized as providing a simplistic characterization of the genetic basis of illness risk that ignores the effects of other variants within the same biological pathways. Of candidate biological pathways for schizophrenia (SZ), the cell adhesion molecule (CAM) pathway has repeatedly been linked to both psychosis and neurocognitive dysfunction. Here we tested, using risk allele scores derived from the Schizophrenia Psychiatric Genome-Wide Association Study Consortium (PGC-SCZ), whether alleles within the CAM pathway were correlated with poorer neuropsychological function in patients. METHOD: In total, 424 patients with psychosis were assessed in areas of cognitive ability typically found to be impaired in SZ: intelligence quotient, memory, working memory and attentional control. CAM pathway genes were identified using the Kyoto Encyclopedia of Genes and Genomes (KEGG) database. Alleles within these genes identified as significantly associated with SZ risk in the PGC-SCZ were then used to calculate a CAM pathway-based polygenic risk allele score for each patient and these scores were tested for association with cognitive ability. RESULTS: Increased CAM pathway polygenic risk scores were significantly associated with poorer performance on measures of memory and attention, explaining 1-3% of variation on these measures. Notably, the most strongly associated single nucleotide polymorphism (SNP) in the CAM pathway (rs9272105 within HLA-DQA1) explained a similar amount of variance in attentional control, but not memory, as the polygenic risk analysis. CONCLUSIONS: These data support a role for the CAM pathway in cognitive function, both at the level of individual SNPs and the wider pathway. In so doing these data highlight the value of pathway-based polygenic risk score studies as well as single gene studies for understanding SZ-associated deficits in cognition.

Detecting large copy number variants using exome genotyping arrays in a large Swedish schizophrenia sample.

Although copy number variants (CNVs) are important in genomic medicine, CNVs have not been systematically assessed for many complex traits. Several large rare CNVs increase risk for schizophrenia (SCZ) and autism and often demonstrate pleiotropic effects; however, their frequencies in the general population and other complex traits are unknown. Genotyping large numbers of samples is essential for progress. Large cohorts from many different diseases are being genotyped using exome-focused arrays designed to detect uncommon or rare protein-altering sequence variation. Although these arrays were not designed for CNV detection, the hybridization intensity data generated in each experiment could, in principle, be used for gene-focused CNV analysis. Our goal was to evaluate the extent to which CNVs can be detected using data from one particular exome array (the Illumina Human Exome Bead Chip). We genotyped 9100 Swedish subjects (3962 cases with SCZ and 5138 controls) using both standard genome-wide association study (GWAS) and exome arrays. In comparison with CNVs detected using GWAS arrays, we observed high sensitivity and specificity for detecting genic CNVs 400 kb including known pathogenic CNVs along with replicating the literature finding that cases with SCZ had greater enrichment for genic CNVs. Our data confirm the association of SCZ with 16p11.2 duplications and 22q11.2 deletions, and suggest a novel association with deletions at 11q12.2. Our results suggest the utility of exome-focused arrays in surveying large genic CNVs in very large samples; and thereby open the door for new opportunities such as conducting well-powered CNV assessment and comparisons between different diseases. The use of a single platform also minimizes potential confounding factors that could impact accurate detection.

Meta-analysis of association between obsessive-compulsive disorder and the 3' region of neuronal glutamate transporter gene SLC1A1.

The neuronal glutamate transporter gene SLC1A1 is a candidate gene for obsessive-compulsive disorder (OCD) based on linkage studies and convergent evidence implicating glutamate in OCD etiology. The 3' end of SLC1A1 is the only genomic region with consistently demonstrated OCD association, especially when analyzing male-only probands. However, specific allele associations have not been consistently replicated, and recent OCD genome-wide association and meta-analysis studies have not incorporated all previously associated SLC1A1 SNPs. To clarify the nature of association between SLC1A1 and OCD, pooled analysis was performed on all available relevant raw study data, comprising a final sample of 815 trios, 306 cases and 634 controls. This revealed weak association between OCD and one of nine tested SLC1A1 polymorphisms (rs301443; uncorrected P = 0.046; non-significant corrected P). Secondary analyses of male-affecteds only (N = 358 trios and 133 cases) demonstrated modest association between OCD and a different SNP (rs12682807; uncorrected P = 0.012; non-significant corrected P). Findings of this meta-analysis are consistent with the trend of previous candidate gene studies in psychiatry and do not clarify the putative role of SLC1A1 in OCD pathophysiology. Nonetheless, it may be important to further examine the potential associations demonstrated in this amalgamated sample, especially since the SNPs with modest associations were not included in the more highly powered recent GWAS or in a past meta-analysis including five SLC1A1 polymorphisms. This study underscores the need for much larger sample sizes in future genetic association studies and suggests that next-generation sequencing may be beneficial in examining the potential role of rare variants in OCD.

Cis-acting regulation of brain-specific ANK3 gene expression by a genetic variant associated with bipolar disorder.

Several genome-wide association studies for bipolar disorder (BD) have found a strong association of the Ankyrin 3 (ANK3) gene. This association spans numerous linked single-nucleotide polymorphisms (SNPs) in an ~250-kb genomic region overlapping ANK3. The associated region encompasses predicted regulatory elements as well as two of the six validated alternative first exons, which encode distinct protein domains at the N-terminus of the protein also known as Ankyrin-G. Using RNA ligase-mediated rapid amplification of cDNA ends to identify novel transcripts in conjunction with a highly sensitive, exon-specific multiplexed mRNA expression assay, we detected differential regulation of distinct ANK3 transcription start sites and coupling of specific 5' ends with 3' mRNA splicing events in postmortem human brain and human stem cell-derived neural progenitors and neurons. Furthermore, allelic variation at the BD-associated SNP rs1938526 correlated with a significant difference in cerebellar expression of a brain-specific ANK3 transcript. These findings suggest a brain-specific cis-regulatory transcriptional effect of ANK3 that may be relevant to BD pathophysiology.

Genome-wide association study in a Swedish population yields support for greater CNV and MHC involvement in schizophrenia compared with bipolar disorder.

Schizophrenia (SCZ) and bipolar disorder (BD) are highly heritable psychiatric disorders with overlapping susceptibility loci and symptomatology. We conducted a genome-wide association study (GWAS) of these disorders in a large Swedish sample. We report a new and independent case-control analysis of 1507 SCZ cases, 836 BD cases and 2093 controls. No single-nucleotide polymorphisms (SNPs) achieved significance in these new samples; however, combining new and previously reported SCZ samples (2111 SCZ and 2535 controls) revealed a genome-wide significant association in the major histocompatibility complex (MHC) region (rs886424, P=4.54 × 10(-8)). Imputation using multiple reference panels and meta-analysis with the Psychiatric Genomics Consortium SCZ results underscored the broad, significant association in the MHC region in the full SCZ sample. We evaluated the role of copy number variants (CNVs) in these subjects. As in prior reports, deletions were enriched in SCZ, but not BD cases compared with controls. Singleton deletions were more frequent in both case groups compared with controls (SCZ: P=0.003, BD: P=0.013), whereas the largest CNVs (>500 kb) were significantly enriched only in SCZ cases (P=0.0035). Two CNVs with previously reported SCZ associations were also overrepresented in this SCZ sample: 16p11.2 duplications (P=0.0035) and 22q11 deletions (P=0.03). These results reinforce prior reports of significant MHC and CNV associations in SCZ, but not BD.

Molecular pathways involved in neuronal cell adhesion and membrane scaffolding contribute to schizophrenia and bipolar disorder susceptibility.

Susceptibility to schizophrenia and bipolar disorder may involve a substantial, shared contribution from thousands of common genetic variants, each of small effect. Identifying whether risk variants map to specific molecular pathways is potentially biologically informative. We report a molecular pathway analysis using the single-nucleotide polymorphism (SNP) ratio test, which compares the ratio of nominally significant (P<0.05) to nonsignificant SNPs in a given pathway to identify the 'enrichment' for association signals. We applied this approach to the discovery (the International Schizophrenia Consortium (n=6909)) and validation (Genetic Association Information Network (n=2729)) of schizophrenia genome-wide association study (GWAS) data sets. We investigated each of the 212 experimentally validated pathways described in the Kyoto Encyclopaedia of Genes and Genomes in the discovery sample. Nominally significant pathways were tested in the validation sample, and five pathways were found to be significant (P=0.03-0.001); only the cell adhesion molecule (CAM) pathway withstood conservative correction for multiple testing. Interestingly, this pathway was also significantly associated with bipolar disorder (Wellcome Trust Case Control Consortium (n=4847)) (P=0.01). At a gene level, CAM genes associated in all three samples (NRXN1 and CNTNAP2), which were previously implicated in specific language disorder, autism and schizophrenia. The CAM pathway functions in neuronal cell adhesion, which is critical for synaptic formation and normal cell signaling. Similar pathways have also emerged from a pathway analysis of autism, suggesting that mechanisms involved in neuronal cell adhesion may contribute broadly to neurodevelopmental psychiatric phenotypes.