High-throughput genomic technologies yield about 20,000 variants in the protein-coding exome of each individual. A commonly used approach to select candidate disease-causing variants is to test whether the associated gene has been previously reported to be disease-causing. In the absence of known disease-causing genes, it can be challenging to associate candidate genes with specific genetic diseases. To facilitate the discovery of novel gene-disease associations, we determined the putative biologically closest known genes and their associated diseases for 13,005 human genes not currently reported to be disease-associated. We used these data to construct the closest disease-causing genes (CDG) server, which can be used to infer the closest genes with an associated disease for a user-defined list of genes or diseases. We demonstrate the utility of the CDG server in five immunodeficiency patient exomes across different diseases and modes of inheritance, where CDG dramatically reduced the number of candidate genes to be evaluated. This resource will be a considerable asset for ascertaining the potential relevance of genetic variants found in patient exomes to specific diseases of interest. The CDG database and online server are freely available to non-commercial users at: http://lab.rockefeller.edu/casanova/CDG.

Globally, autosomal recessive IFNAR1 deficiency is a rare inborn error of immunity underlying susceptibility to live attenuated vaccine and wild-type viruses. We report seven children from five unrelated kindreds of western Polynesian ancestry who suffered from severe viral diseases. All the patients are homozygous for the same nonsense IFNAR1 variant (p.Glu386*). This allele encodes a truncated protein that is absent from the cell surface and is loss-of-function. The fibroblasts of the patients do not respond to type I IFNs (IFN-α2, IFN-ω, or IFN-β). Remarkably, this IFNAR1 variant has a minor allele frequency >1% in Samoa and is also observed in the Cook, Society, Marquesas, and Austral islands, as well as Fiji, whereas it is extremely rare or absent in the other populations tested, including those of the Pacific region. Inherited IFNAR1 deficiency should be considered in individuals of Polynesian ancestry with severe viral illnesses.

Cellular senescence is a mechanism used by mitotic cells to prevent uncontrolled cell division. As senescent cells persist in tissues, they cause local inflammation and are harmful to surrounding cells, contributing to aging. Generally, neurodegenerative diseases, such as Parkinson's, are disorders of aging. The contribution of cellular senescence to neurodegeneration is still unclear. SATB1 is a DNA binding protein associated with Parkinson's disease. We report that SATB1 prevents cellular senescence in post-mitotic dopaminergic neurons. Loss of SATB1 causes activation of a cellular senescence transcriptional program in dopamine neurons both in human stem cell-derived dopaminergic neurons and in mice. We observed phenotypes that are central to cellular senescence in SATB1 knockout dopamine neurons in vitro and in vivo. Moreover, we found that SATB1 directly represses expression of the pro-senescence factor p21 in dopaminergic neurons. Our data implicate senescence of dopamine neurons as a contributing factor in the pathology of Parkinson's disease.

Interindividual clinical variability in the course of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection is vast. We report that at least 101 of 987 patients with life-threatening coronavirus disease 2019 (COVID-19) pneumonia had neutralizing immunoglobulin G (IgG) autoantibodies (auto-Abs) against interferon-ω (IFN-ω) (13 patients), against the 13 types of IFN-α (36), or against both (52) at the onset of critical disease; a few also had auto-Abs against the other three type I IFNs. The auto-Abs neutralize the ability of the corresponding type I IFNs to block SARS-CoV-2 infection in vitro. These auto-Abs were not found in 663 individuals with asymptomatic or mild SARS-CoV-2 infection and were present in only 4 of 1227 healthy individuals. Patients with auto-Abs were aged 25 to 87 years and 95 of the 101 were men. A B cell autoimmune phenocopy of inborn errors of type I IFN immunity accounts for life-threatening COVID-19 pneumonia in at least 2.6% of women and 12.5% of men.

We demonstrate that stress differentially regulates glutamate homeostasis in the dorsal and ventral hippocampus and identify a role for the astroglial xCT in ventral dentate gyrus (vDG) in stress and antidepressant responses. We provide an RNA-seq roadmap for the stress-sensitive vDG. The transcription factor REST binds to xCT promoter in co-occupancy with the epigenetic marker H3K27ac to regulate expression of xCT, which is also reduced in a genetic mouse model of inherent susceptibility to depressive-like behavior. Pharmacologically, modulating histone acetylation with acetyl-L-carnitine (LAC) or acetyl-N-cysteine (NAC) rapidly increases xCT and activates a network with mGlu2 receptors to prime an enhanced glutamate homeostasis that promotes both pro-resilient and antidepressant-like responses. Pharmacological xCT blockage counteracts NAC prophylactic effects. GFAP+-Cre-dependent overexpression of xCT in vDG mimics pharmacological actions in promoting resilience. This work establishes a mechanism by which vDG protection leads to stress resilience and antidepressant responses via epigenetic programming of an xCT-mGlu2 network.

Autosomal dominant (AD) NFKB1 deficiency is thought to be the most common genetic etiology of common variable immunodeficiency (CVID). However, the causal link between NFKB1 variants and CVID has not been demonstrated experimentally and genetically, and there has been insufficient biochemical characterization and enrichment analysis. We show that the cotransfection of NFKB1-deficient HEK293T cells (lacking both p105 and its cleaved form p50) with a κB reporter, NFKB1/p105, and a homodimerization-defective RELA/p65 mutant results in p50:p65 heterodimer-dependent and p65:p65 homodimer-independent transcriptional activation. We found that 59 of the 90 variants in patients with CVID or related conditions were loss of function or hypomorphic. By contrast, 258 of 260 variants in the general population or patients with unrelated conditions were neutral. None of the deleterious variants displayed negative dominance. The enrichment in deleterious NFKB1 variants of patients with CVID was selective and highly significant (P = 2.78 × 10-15). NFKB1 variants disrupting NFKB1/p50 transcriptional activity thus underlie AD CVID by haploinsufficiency, whereas neutral variants in this assay should not be considered causal.

Autosomal inborn errors of type I IFN immunity and autoantibodies against these cytokines underlie at least 10% of critical COVID-19 pneumonia cases. We report very rare, biochemically deleterious X-linked TLR7 variants in 16 unrelated male individuals aged 7 to 71 years (mean: 36.7 years) from a cohort of 1,202 male patients aged 0.5 to 99 years (mean: 52.9 years) with unexplained critical COVID-19 pneumonia. None of the 331 asymptomatically or mildly infected male individuals aged 1.3 to 102 years (mean: 38.7 years) tested carry such TLR7 variants (p = 3.5 × 10-5). The phenotypes of five hemizygous relatives of index cases infected with SARS-CoV-2 include asymptomatic or mild infection (n=2, 5 and 38 years), or moderate (n=1, 5 years), severe (n=1, 27 years), or critical (n=1, 29 years) pneumonia. Two boys (aged 7 and 12 years) from a cohort of 262 male patients with severe COVID-19 pneumonia (mean: 51.0 years) are hemizygous for a deleterious TLR7 variant. The cumulative allele frequency for deleterious TLR7 variants in the male general population is < 6.5x10-4 We also show that blood B cell lines and myeloid cell subsets from the patients do not respond to TLR7 stimulation, a phenotype rescued by wild-type TLR7 The patients' blood plasmacytoid dendritic cells (pDCs) produce low levels of type I IFNs in response to SARS-CoV-2. Overall, X-linked recessive TLR7 deficiency is a highly penetrant genetic etiology of critical COVID-19 pneumonia, in about 1.8% of male patients below the age of 60 years. Human TLR7 and pDCs are essential for protective type I IFN immunity against SARS-CoV-2 in the respiratory tract.

Idiopathic Parkinson's Disease (PD) is characterized by the loss of dopaminergic neurons in the substantia nigra pars compacta, which is associated with neuroinflammation and reactive gliosis. The underlying cause of PD and the concurrent neuroinflammation are not well understood. In this study, we utilized human and murine neuronal lines, stem cell-derived dopaminergic neurons, and mice to demonstrate that three previously identified genetic risk factors for PD, namely SATB1, MIR22HG, and GBA, are components of a single gene regulatory pathway. Our findings indicate that dysregulation of this pathway leads to the upregulation of glucocerebrosides (GluCer), which triggers a cellular senescence-like phenotype in dopaminergic neurons. Specifically, we discovered that downregulation of the transcriptional repressor SATB1 results in the derepression of the microRNA miR-22-3p, leading to decreased GBA expression and subsequent accumulation of GluCer. Furthermore, our results demonstrate that an increase in GluCer alone is sufficient to impair lysosomal and mitochondrial function, thereby inducing cellular senescence dependent on S100A9 and stress factors. Dysregulation of the SATB1-MIR22-GBA pathway, observed in both PD patients and normal aging, leads to lysosomal and mitochondrial dysfunction due to the GluCer accumulation, ultimately resulting in a cellular senescence-like phenotype in dopaminergic neurons. Therefore, our study highlights a novel pathway involving three genetic risk factors for PD and provides a potential mechanism for the senescence-induced neuroinflammation and reactive gliosis observed in both PD and normal aging.

Most humans are exposed to Tropheryma whipplei (Tw). Whipple's disease (WD) strikes only a small minority of individuals infected with Tw (<0.01%), whereas asymptomatic chronic carriage is more common (<25%). We studied a multiplex kindred, containing four WD patients and five healthy Tw chronic carriers. We hypothesized that WD displays autosomal dominant (AD) inheritance, with age-dependent incomplete penetrance. We identified a single very rare non-synonymous mutation in the four patients: the private R98W variant of IRF4, a transcription factor involved in immunity. The five Tw carriers were younger, and also heterozygous for R98W. We found that R98W was loss-of-function, modified the transcriptome of heterozygous leukocytes following Tw stimulation, and was not dominant-negative. We also found that only six of the other 153 known non-synonymous IRF4 variants were loss-of-function. Finally, we found that IRF4 had evolved under purifying selection. AD IRF4 deficiency can underlie WD by haploinsufficiency, with age-dependent incomplete penetrance.

SARS-CoV-2 is responsible for the ongoing world-wide pandemic which has already taken more than two million lives. Effective treatments are urgently needed. The enzymatic activity of the HECT-E3 ligase family members has been implicated in the cell egression phase of deadly RNA viruses such as Ebola through direct interaction of its VP40 Protein. Here we report that HECT-E3 ligase family members such as NEDD4 and WWP1 interact with and ubiquitylate the SARS-CoV-2 Spike protein. Furthermore, we find that HECT family members are overexpressed in primary samples derived from COVID-19 infected patients and COVID-19 mouse models. Importantly, rare germline activating variants in the NEDD4 and WWP1 genes are associated with severe COVID-19 cases. Critically, I3C, a natural NEDD4 and WWP1 inhibitor from Brassicaceae, displays potent antiviral effects and inhibits viral egression. In conclusion, we identify the HECT family members of E3 ligases as likely novel biomarkers for COVID-19, as well as new potential targets of therapeutic strategy easily testable in clinical trials in view of the established well-tolerated nature of the Brassicaceae natural compounds.

Clinical outcome upon infection with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) ranges from silent infection to lethal coronavirus disease 2019 (COVID-19). We have found an enrichment in rare variants predicted to be loss-of-function (LOF) at the 13 human loci known to govern Toll-like receptor 3 (TLR3)- and interferon regulatory factor 7 (IRF7)-dependent type I interferon (IFN) immunity to influenza virus in 659 patients with life-threatening COVID-19 pneumonia relative to 534 subjects with asymptomatic or benign infection. By testing these and other rare variants at these 13 loci, we experimentally defined LOF variants underlying autosomal-recessive or autosomal-dominant deficiencies in 23 patients (3.5%) 17 to 77 years of age. We show that human fibroblasts with mutations affecting this circuit are vulnerable to SARS-CoV-2. Inborn errors of TLR3- and IRF7-dependent type I IFN immunity can underlie life-threatening COVID-19 pneumonia in patients with no prior severe infection.

Insulin resistance (IR) is a metabolic dysfunction often co-morbid with major depressive disorder (MDD). The paths to development of MDD remain largely unspecified, highlighting a need for identification of risk factors. Here, we tested whether specific subscales of childhood trauma as well as family history of type-2 diabetes (Fam-Hx-Dm2) are risk factors for development of metabolic dysfunction and severity of depressive symptoms.

Major depressive disorder (MDD) is a primary psychiatric illness worldwide; there is a dearth of new mechanistic models for the development of better therapeutic strategies. Although we continue to discover individual biological factors, a major challenge is the identification of integrated, multidimensional traits underlying the complex heterogeneity of depression and treatment outcomes. Here, we set out to ascertain the emergence of the novel mitochondrial mediator of epigenetic function acetyl-L-carnitine (LAC) in relation to previously described individual predictors of antidepressant responses to the insulin-sensitizing agent pioglitazone. Herein, we report that i) subjects with MDD and shorter leukocyte telomere length (LTL) show decreased levels of LAC, increased BMI, and a history of specific types of childhood trauma; and that ii) these multidimensional factors spanning mitochondrial metabolism, cellular aging, metabolic function, and childhood trauma provide more detailed signatures to predict longitudinal changes in depression severity in response to pioglitazone than individual factors. The findings of multidimensional signatures involved in the pathophysiology of depression and their role in predicting treatment outcomes provide a starting point for the development of a mechanistic framework linking biological networks and environmental factors to clinical outcomes in pursuit of personalized medicine strategies to effectively treat MDD.

Population stratification is a confounder of genetic association studies. In analyses of rare variants, corrections based on principal components (PCs) and linear mixed models (LMMs) yield conflicting conclusions. Studies evaluating these approaches generally focused on limited types of structure and large sample sizes. We investigated the properties of several correction methods through a large simulation study using real exome data, and several within- and between-continent stratification scenarios. We considered different sample sizes, with situations including as few as 50 cases, to account for the analysis of rare disorders. Large samples showed that accounting for stratification was more difficult with a continental than with a worldwide structure. When considering a sample of 50 cases, an inflation of type-I-errors was observed with PCs for small numbers of controls (≤ 100), and with LMMs for large numbers of controls (≥ 1000). We also tested a novel local permutation method (LocPerm), which maintained a correct type-I-error in all situations. Powers were equivalent for all approaches pointing out that the key issue is to properly control type-I-errors. Finally, we found that power of analyses including small numbers of cases can be increased, by adding a large panel of external controls, provided an appropriate stratification correction was used.

Autosomal recessive IRF7 and IRF9 deficiencies impair type I and III IFN immunity and underlie severe influenza pneumonitis. We report three unrelated children with influenza A virus (IAV) infection manifesting as acute respiratory distress syndrome (IAV-ARDS), heterozygous for rare TLR3 variants (P554S in two patients and P680L in the third) causing autosomal dominant (AD) TLR3 deficiency. AD TLR3 deficiency can underlie herpes simplex virus-1 (HSV-1) encephalitis (HSE) by impairing cortical neuron-intrinsic type I IFN immunity to HSV-1. TLR3-mutated leukocytes produce normal levels of IFNs in response to IAV. In contrast, TLR3-mutated fibroblasts produce lower levels of IFN-β and -λ, and display enhanced viral susceptibility, upon IAV infection. Moreover, the patients' iPSC-derived pulmonary epithelial cells (PECs) are susceptible to IAV. Treatment with IFN-α2b or IFN-λ1 rescues this phenotype. AD TLR3 deficiency may thus underlie IAV-ARDS by impairing TLR3-dependent, type I and/or III IFN-mediated, PEC-intrinsic immunity. Its clinical penetrance is incomplete for both IAV-ARDS and HSE, consistent with their typically sporadic nature.

The detection of copy number variations (CNVs) in whole-exome sequencing (WES) data is important, as CNVs may underlie a number of human genetic disorders. The recently developed HMZDelFinder algorithm can detect rare homozygous and hemizygous (HMZ) deletions in WES data more effectively than other widely used tools. Here, we present HMZDelFinder_opt, an approach that outperforms HMZDelFinder for the detection of HMZ deletions, including partial exon deletions in particular, in WES data from laboratory patient collections that were generated over time in different experimental conditions. We show that using an optimized reference control set of WES data, based on a PCA-derived Euclidean distance for coverage, strongly improves the detection of HMZ complete exon deletions both in real patients carrying validated disease-causing deletions and in simulated data. Furthermore, we develop a sliding window approach enabling HMZDelFinder_opt to identify HMZ partial deletions of exons that are undiscovered by HMZDelFinder. HMZDelFinder_opt is a timely and powerful approach for detecting HMZ deletions, particularly partial exon deletions, in WES data from inherently heterogeneous laboratory patient collections.

In eukaryotes, carnitine is best known for its ability to shuttle esterified fatty acids across mitochondrial membranes for β-oxidation. It also returns to the cytoplasm, in the form of acetyl-L-carnitine (LAC), some of the resulting acetyl groups for posttranslational protein modification and lipid biosynthesis. While dietary LAC supplementation has been clinically investigated, its effects on cellular metabolism are not well understood. To explain how exogenous LAC influences mammalian cell metabolism, we synthesized isotope-labeled forms of LAC and its analogs. In cultures of glucose-limited U87MG glioma cells, exogenous LAC contributed more robustly to intracellular acetyl-CoA pools than did β-hydroxybutyrate, the predominant circulating ketone body in mammals. The fact that most LAC-derived acetyl-CoA is cytosolic is evident from strong labeling of fatty acids in U87MG cells by exogenous 13C2-acetyl-L-carnitine. We found that the addition of d3-acetyl-L-carnitine increases the supply of acetyl-CoA for cytosolic posttranslational modifications due to its strong kinetic isotope effect on acetyl-CoA carboxylase, the first committed step in fatty acid biosynthesis. Surprisingly, whereas cytosolic carnitine acetyltransferase is believed to catalyze acetyl group transfer from LAC to coenzyme A, CRAT-/- U87MG cells were unimpaired in their ability to assimilate exogenous LAC into acetyl-CoA. We identified carnitine octanoyltransferase as the key enzyme in this process, implicating a role for peroxisomes in efficient LAC utilization. Our work has opened the door to further biochemical investigations of a new pathway for supplying acetyl-CoA to certain glucose-starved cells.

Multisystem inflammatory syndrome in children (MIS-C) is a rare and severe condition that follows benign COVID-19. We report autosomal recessive deficiencies of OAS1, OAS2, or RNASEL in five unrelated children with MIS-C. The cytosolic double-stranded RNA (dsRNA)-sensing OAS1 and OAS2 generate 2'-5'-linked oligoadenylates (2-5A) that activate the single-stranded RNA-degrading ribonuclease L (RNase L). Monocytic cell lines and primary myeloid cells with OAS1, OAS2, or RNase L deficiencies produce excessive amounts of inflammatory cytokines upon dsRNA or severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) stimulation. Exogenous 2-5A suppresses cytokine production in OAS1-deficient but not RNase L-deficient cells. Cytokine production in RNase L-deficient cells is impaired by MDA5 or RIG-I deficiency and abolished by mitochondrial antiviral-signaling protein (MAVS) deficiency. Recessive OAS-RNase L deficiencies in these patients unleash the production of SARS-CoV-2-triggered, MAVS-mediated inflammatory cytokines by mononuclear phagocytes, thereby underlying MIS-C.