Inherited deficiency of the RNA lariat–debranching enzyme 1 (DBR1) is a rare etiology of brainstem viral encephalitis. The cellular basis of disease and the range of viral predisposition are unclear. We report inherited DBR1 deficiency in a 14-year-old boy who suffered from isolated SARS-CoV-2 brainstem encephalitis. The patient is homozygous for a previously reported hypomorphic and pathogenic DBR1 variant (I120T). Consistently, DBR1 I120T/I120T fibroblasts from affected individuals from this and another unrelated kindred have similarly low levels of DBR1 protein and high levels of RNA lariats. DBR1 I120T/I120T human pluripotent stem cell (hPSC)–derived hindbrain neurons are highly susceptible to SARS-CoV-2 infection. Exogenous WT DBR1 expression in DBR1 I120T/I120T fibroblasts and hindbrain neurons rescued the RNA lariat accumulation phenotype. Moreover, expression of exogenous RNA lariats, mimicking DBR1 deficiency, increased the susceptibility of WT hindbrain neurons to SARS-CoV-2 infection. Inborn errors of DBR1 impair hindbrain neuron–intrinsic antiviral immunity, predisposing to viral infections of the brainstem, including that by SARS-CoV-2.
Introduction
The first cases of SARS-CoV-2 encephalitis were described in March 2020, about 5 mo after the start of the COVID-19 pandemic (Moriguchi et al., 2020; Ye et al., 2020). It progressively became clear that encephalitis following SARS-CoV-2 infection comprises several different types of encephalopathy, probably due to different mechanisms of disease (Aubart et al., 2022; Cho et al., 2023; Ellul et al., 2020; Siow et al., 2021). These SARS-CoV-2–related forms of encephalitis have an estimated overall prevalence of about 2/10,000. In most cases, encephalitis occurs weeks or months after the acute respiratory viral infection and presents as post-infectious autoimmune encephalitis (Holroyd and Conway, 2023). However, encephalitis may also occur in severe cases of COVID-19 pneumonia or multisystem inflammatory syndrome in children (MIS-C) with the inflammation of multiple organs, including the brain (Olivotto et al., 2021). No particular distribution by age, sex, or ancestry has been observed in either of these types of SARS-CoV-2–related encephalitis, and overall mortality for these two forms is about 13%. Isolated encephalitis, which defines a third class of SARS-CoV-2–related encephalopathy, is extremely severe and has been reported only rarely (Cho et al., 2023; Ellul et al., 2020; Moriguchi et al., 2020; Siow et al., 2021; Ye et al., 2020). Acute SARS-CoV-2 encephalitis often involves brainstem lesions, as revealed by brain imaging data, but lesions in other parts of the brain have also been reported. Over the last 3 years, we and others have shown that critical COVID-19 pneumonia results from inborn errors of, or autoantibodies against, type I interferons (IFNs) in about 15% of cases (Bastard et al., 2021a; Casanova and Abel, 2021, 2022; Zhang et al., 2022). We also recently discovered recessive deficiencies of the OAS-RNase L pathway resulting in unchecked inflammatory responses to SARS-CoV-2 in mononuclear phagocytes, underlying MIS-C in about 1% of the patients studied (Lee et al., 2023). The pathogenesis of SARS-CoV-2 encephalitis remains much less clear. Based on our previous studies of other types of viral encephalitis, including herpes simplex virus 1 encephalitis in particular (Zhang et al., 2021), we hypothesized that inborn errors of brain-intrinsic immunity might underlie isolated SARS-CoV-2 encephalitis. We tested this hypothesis by performing whole-exome/genome sequencing for a cohort of 16 patients from the COVID Human Genetic Effort (https://www.covidhge.com) cohort who developed isolated encephalitis during acute SARS-CoV-2 infection.
Results and discussion
We analyzed the 16 exomes for candidate genotypes for the 19 known viral encephalitis–causing genes (TLR3, UNC93B1, TRIF, TRAF3, TBK1, IRF3, NEMO, GTF3A, MDA5, DOCK2, POL3A, POL3C, IFNAR1, STAT1, STAT2, TYK2, SNORA31, DBR1, and RIPK3) (Casanova and Abel, 2021; Liu et al., 2023; Zhang et al., 2021), and for another three genes implicated in inborn errors of immunity (IEI) and with functions closely related to those of the known viral encephalitis–causing genes (IFNAR2, IRF9, and JAK1) (Meyts and Casanova, 2021). We searched for rare (minor allele frequency <0.01 in the Genome Aggregation Database and in our in-house whole exome sequencing database containing about 20,000 exomes or genomes) non-synonymous or splicing (affecting essential-splicing or intronic branch-point sites) variants with a combined annotation-dependent depletion score (Kircher et al., 2014) above the mutation significance cutoff (Itan et al., 2016) for the 22 genes, all of which have a gene damage index below 13.83 (Itan et al., 2015). We considered monoallelic and/or biallelic variants in accordance with the known modes of inheritance of the genes concerned. Interestingly, this search revealed that one patient (P1) was homozygous for a missense mutation (I120T) of DBR1, which encodes the RNA lariat–debranching enzyme 1 (DBR1), a previously reported genetic etiology of brainstem viral encephalitis (BVE) (Zhang et al., 2018). More specifically, a detailed analysis of the whole-genome sequencing (WGS) data for P1 and his two parents revealed that P1 was homozygous or compound heterozygous for 18 and 1 genes, respectively (Table S1). In addition, P1 was found to be heterozygous for a known familial Mediterranean fever (FMF)–causing pathogenic MEFV gene variant, M694V (Tirosh et al., 2021), and another MEFV variant, E148Q, previously reported to confer mild genotype modifications (Fig. 1 A). The DBR1 variant of P1 was previously shown to be pathogenic, underlying HSV-1 BVE (Zhang et al., 2018). An EstiAge analysis of sequencing data for P1 and a previously described patient with HSV-1 BVE showed that they had a haplotype of about 6.61 Mb in common, estimated to have originated from a common ancestor about 21 generations ago (95% confidence interval: 7–85 generations), corresponding to a period of about 567 years (297–2,295 years) (Fig. 1 B). None of the other genes were deemed plausible candidates to explain SARS-CoV-2 encephalitis based on the biochemical nature of the variant, the reported tissue expression of the corresponding gene, or the known function of its products.
P1 is a 14-year-old boy born to nonconsanguineous parents of Syrian origin living in Sweden. P1 and one of his six siblings (S6) had a clinical history of FMF, with recurrent episodes of serositis. The familial segregation of the DBR1 I120T and MEFV M694V variants suggested incomplete penetrance of the two genotypes for the BVE and FMF clinical phenotypes, respectively (Fig. 1 A). In February 2022, P1 presented with a sore throat and fever at a time at which other members of his family had mild COVID-19. 1 wk later, he was admitted to an intensive care unit at Umeå University Hospital, Sweden, for acute dizziness, vomiting, and loss of consciousness. Magnetic resonance imaging (MRI) showed brainstem encephalitis, with lesions in the pons, mesencephalon, and cerebellum, and evidence of an increase in intracranial pressure (Fig. 1 C). At admission, PCR on a nasopharyngeal sample was positive for SARS-CoV-2, but no IgG against SARS-CoV-2 was detected in the blood or cerebral spinal fluid (CSF). SARS-CoV-2 was therefore considered to be the virus responsible for encephalitis (Table S2). Extensive searches for other pathogens in the blood and CSF, and screening for neuronal autoantibodies in the blood yielded negative results (Table S3). P1 was treated with an external ventricular drain, corticosteroids, remdesivir, and IL1 blockade, and has no sequelae 15 mo after the encephalitis episode. Immunological studies 6 mo after the recovery of the patient showed leukocyte subset counts and proliferation to be normal. The activation of T and B lymphocytes, and the functions of neutrophils, phagocytes, and the complement system were also normal (Table S3). P1 tested positive for SARS-CoV-2 IgG in serological tests performed at this time point and presented a robust T cell response to the SARS-CoV-2 spike 1 (S1), spike (S), membrane (M), and nucleocapsid (N) peptides, further confirming the history of SARS-CoV-2 infection in this patient.
DBR1 has only one known cellular function, as the only intronic RNA lariat debranching enzyme (Chapman and Boeke, 1991). We previously reported autosomal recessive DBR1 deficiency in five children from three kindreds with brainstem encephalitis due to infection with HSV-1, influenza B virus, or norovirus, with complete clinical penetrance for at least one type of viral encephalitis (Zhang et al., 2018). The I120T mutation was previously found in the homozygous state in two patients with HSV-1 BVE. The I120T variant has been shown to be biochemically deleterious, as it results in abnormally low levels of the corresponding protein and RNA lariat–debranching activity (Zhang et al., 2018). P1 and three of his siblings, aged 6 (S6), 10 (S4), and 20 (S2) years, were found to be homozygous for the I120T variant. None of these three siblings of P1 developed SARS-CoV-2 encephalitis or another severe infection, despite living in the same household as P1 and probably being exposed to the same infectious agents as P1, including viruses (Fig. S1 A). The 6-year-old sibling (S6) was the only sibling with a positive T cell response to SARS-CoV-2 S1 and SMN peptides at the latest follow-up, and this child also had a positive result in the anti-SARS-CoV-2 spike protein serological test (Table S4). The parents were heterozygous and the other three siblings tested (S1, S3, and S5) were heterozygous or WT at the I120 position (Fig. 1 A and Fig. S1 B). These findings are suggestive of autosomal recessive DBR1 deficiency underlying SARS-CoV-2 brainstem encephalitis with incomplete clinical penetrance, consistent with previous reports of IEIs underlying sporadic severe viral diseases (Casanova and Abel, 2021, 2022).
For confirmation of the DBR1 deficiency in P1 at the cellular level, we first studied SV40-transformed fibroblasts (SV40-fibroblasts) as a surrogate cellular model for tissue-resident cells, as in our previous studies (Zhang et al., 2018). SV40-fibroblasts from P1 and one of his siblings (S2) homozygous for DBR1 I120T contain low levels of DBR1 protein and high levels of DKK1 and ID1 RNA lariats, like SV40-fibroblasts from one of the previously reported I120T homozygotes (Zhang et al., 2018). This was not the case for the siblings of P1 heterozygous for I120T or WT at this position (Fig. 2, A–C and Fig. S1 B). Exogenous WT DBR1 expression in SV40-fibroblasts from P1 rescued DBR1 protein levels and decreased DKK1 and ID1 RNA lariat levels to values similar to those in healthy control cells (Fig. 2, D–F). A previous study reported abnormally high levels of IFN-γ in the blood of FMF patients even during attack-free periods (Köklü et al., 2005). However, the carriers of the FMF genotype (M694V) from this family had basal circulating IFN-α, -β, and -γ levels similar to those in individuals not carrying this variant (Fig. S1 C). The FMF genotype of this family does not, therefore, appear to affect the penetrance or expressivity of DBR1 deficiency through changes in blood type I or II IFN levels in the context of viral infection. Autosomal recessive DBR1 deficiency may, therefore, underlie brainstem SARS-CoV-2 encephalitis in P1.
We previously showed that inherited DBR1 deficiency results in high levels of RNA lariat accumulation, particularly during viral infection, thereby impairing cell-intrinsic antiviral immunity in human fibroblasts (Zhang et al., 2018). The brain-specific mechanisms of DBR1 deficiency impairing cell-intrinsic immunity to viruses remain unclear. DBR1 expression is strongest in the brainstem in humans, consistent with the clinical presentation of the known DBR1-deficient patients, who display selective susceptibility to BVE. We hypothesized that DBR1 deficiency due to homozygosity for the I120T variant would result in uncontrolled SARS-CoV-2 infection in brainstem-resident neuronal cells due to the accumulation of RNA lariats. Hindbrain neurons differentiated from human pluripotent stem cells (hPSC) express angiotensin-converting enzyme 2 (ACE2) and should therefore be permissive for ACE2-mediated SARS-CoV-2 entry (Fig. 3 A). We first assessed the levels of RNA lariats in hindbrain neurons derived from hPSCs from a previously reported DBR1-deficient patient homozygous for the same I120T variant and for similar cells derived from a healthy control (H9). We then assessed the susceptibility to SARS-CoV-2 infection of these cells relative to previously reported patients with forebrain HSV-1 encephalitis and deficiencies of TLR3 or IFNAR1 (TLR3−/−, IFNAR1−/−) (Bastard et al., 2021b; Guo et al., 2011) and a healthy control (Fig. S2, A and B). Like P1-derived SV40-fibroblasts, DBR1 I120T/I120T hPSC-derived hindbrain neurons had higher levels of ID1 and DKK1 RNA lariats than healthy control neurons (Fig. 3, B and C), and exogenous WT DBR1 expression decreased the levels of these RNA lariats (Fig. 3 D). SARS-CoV-2 infection also resulted in an increase in DKK1 lariat levels in both healthy control and DBR1 I120T/I120T hindbrain neurons (Fig. 3 E), as previously reported for the infection of DBR1-mutated dermal fibroblasts with HSV-1 (Zhang et al., 2018). It is, therefore, plausible that SARS-CoV-2 infection is uncontrolled in DBR1-deficient hindbrain neurons due to the accumulation of intronic RNA lariats.
Finally, we studied the susceptibility to SARS-CoV-2 of DBR1-deficient hindbrain neurons relative to cells from TLR3- or IFNAR1-deficient patients with forebrain herpes simplex encephalitis (Bastard et al., 2021b; Guo et al., 2011) and a healthy control. Interestingly, deficiencies of TLR3 or IFNAR1 did not render hindbrain neurons more susceptible to SARS-CoV-2 replication (Fig. 4, A–C; and Fig. S3 A), despite the demonstrated role of such deficiencies in defects of the control of HSV-1 infection in hPSC-derived cortical neurons and of the control of HSV-1, SARS-CoV-2, and other viral infections in human SV40-fibroblasts (Guo et al., 2011; Lim et al., 2019; Zhang et al., 2020; Bastard et al., 2021b). By contrast, DBR1 I120T/I120T hindbrain neurons displayed markedly higher rates of SARS-CoV-2 replication from 24 to 96 h after infection than hindbrain neurons derived from a healthy control or from TLR3−/− or IFNAR1−/− patients (Fig. 4, A–C; and Fig. S3 B). The TLR3−/− and IFNAR1−/− hindbrain neurons are normally resistant to SARS-CoV-2 infection. Moreover, despite normal cellular responses to IFN-β (Fig. 4 D and Fig. S3 C), prior treatment with IFN-β only partially restricted the replication of the virus in DBR1 I120T/I120T hindbrain neurons at a low multiplicity of infection (MOI 0.1), with no detectable restriction at a high MOI (MOI 10) (Fig. 4, A–C; and Fig. S3 A). IFN-β treatment also had no detectable effect in control neurons (Fig. 4, A–C; and Fig. S3 A). Importantly, exogenous DKK1 RNA lariat expression rendered DBR1 WT healthy control hindbrain neurons susceptible to SARS-CoV-2 infection, mimicking DBR1-deficient neurons (Fig. 4, E and F). These findings suggest that type I IFNs are not necessary to control SARS-CoV-2 in hindbrain neurons and that DBR1 is a critical viral restriction factor in these cells.
This study thus identifies inherited DBR1 deficiency as a genetic cause of SARS-CoV-2 brainstem encephalitis. Mechanistically, DBR1 deficiency disrupts intrinsic immunity to SARS-CoV-2 in the human hindbrain, resulting in uncontrolled viral replication and brainstem encephalitis. SARS-CoV-2 is known to be able to invade the central nervous system via the olfactory bulb (Yang et al., 2020; Ziegler et al., 2020), but SARS-CoV-2 encephalitis has generally been considered to be more of an inflammatory condition than a viral condition due to the lack of positive results for SARS-CoV-2 in PCR on the CSF, the observed responsiveness to immunomodulatory or immunosuppressive treatments, such as immunoglobulins or corticosteroids, and the delayed onset relative to infection in some patients (Cho et al., 2023; Ellul et al., 2020; Siow et al., 2021). Our in vitro findings suggest that brain inflammation may stem from uncontrolled SARS-CoV-2 replication, at least in this patient with inherited DBR1 deficiency who developed encephalitis during acute SARS-CoV-2 infection, but we cannot exclude a role for other cellular mechanisms in vivo. It is possible that similar viral replication–related disease mechanisms occurred in other patients with isolated SARS-CoV-2 encephalitis or encephalitis together with a severe infection of the lung or other organs. It is therefore advisable to administer an antiviral treatment, particularly during the early stage of SARS-CoV-2 encephalitis. Future studies should search for mutations of the genes encoding DBR1 and related molecules in other patients with SARS-CoV-2 encephalitis. Our findings also confirm that DBR1 is a gatekeeper of the human brainstem against various viruses, including not only HSV-1, influenza B virus, and norovirus, but also SARS-CoV-2. Patients with DBR1 deficiency should be vaccinated not only against SARS-CoV-2 but probably also against a broader range of viruses; live-attenuated vaccines may be contraindicated and should be avoided. Despite the recessive mode of inheritance of DBR1 deficiency, its penetrance for viral encephalitis appears to be incomplete. It will now be important to decipher the detailed molecular mechanisms by which the accumulation of intronic lariats impairs cell-intrinsic immunity to viruses in the brainstem. Such studies may pave the way for the development of effective preventive or therapeutic measures for patients prone to brainstem encephalitis.
Materials and methods
Human subjects
Informed consent was obtained in Sweden, in accordance with local regulations and a protocol for research on human subjects approved by the Swedish Ethical Review Authority (Dnr 2021-06541-01). Experiments were conducted in the United States and France, in accordance with local regulations and with the approval of the institutional review board of the Rockefeller University and the Institut National de la Santé et de la Recherche Médicale, respectively. Approval was obtained from the French Ethics Committee (Comité de Protection des Personnes), the French National Agency for Medicine and Health Product Safety, the Institut National de la Santé et de la Recherche Médicale in Paris, France (protocol no. C10-13), and the Rockefeller University Institutional Review Board in New York, USA (protocol no. JCA-0700).
WGS
WGS was performed with the Truseq DNA PCR-free protocol (IIIumina) according to the manufacturer’s instructions. Briefly, 1,100 ng genomic DNA was fragmented into fragments of about 350 bp in length with a Covaris E220. Fragmentation was controlled with a Tapestation 4200 (Agilent). The fragments were subjected to end repair and an A-tail and dual-index adaptors (IDT for Illumina; TruSeq DNA UDI) were ligated to the fragments, which were then subjected to double-sided purification to create a narrow fragment-size distribution. Library quantification was performed with the KAPA Library Quantification Kit (Roche). Sequencing was performed on S4 flow cells with a NovaSeq 6000 sequencer (Illumina) in paired-end 150-bp readout mode with the aim of obtaining 400 million read pairs (a mean coverage of about 30×). Demultiplexing was performed with bcl2fastq2 Conversion Software v2.20 (Illumina). The sequences were aligned with the reference human genome sequence (GRCh37) with BWA. Downstream processing was performed with the Genome Analysis Toolkit (GATK), SAM tools, and Picard Tools (https://broadinstitute.github.io/picard/). Variants were called with GATK Unified Genotyper. All calls with a Phred-scaled SNP quality ≤20 were filtered out.
Pan-viral serology using VirScan
The VirScan methodology employs programmable phage immunoprecipitation-sequencing (PhIP-Seq), where viral peptides are epitomized on the outer surface of the T7 bacteriophage for precise antibody detection, followed by next-generation deep sequencing as described in the previous studies (Schubert et al., 2019; Xu et al., 2015). The extensive VirScan library comprises 481,966 62-amino acid peptides, designed with a 14-amino acid overlap, spanning comprehensively across full-length vertebrate, mosquito-borne, and tick-borne viral genomes. Our approach to phage immunoprecipitation and sequencing incorporated nuanced modification of well-established PhIP-Seq protocols (Mandel-Brehm et al., 2019; O’Donovan et al., 2020; Schubert et al., 2019). Briefly, 1 μl of human sera was incubated with 500 μl of the VirScan library for 12–18 h at 4°C. All samples were run in two technical replicates. Antibody-bound phages were further subjected to two rounds of immunoprecipitation utilizing a mix of protein A and protein G magnetic beads (Thermo Fisher Scientific), followed by elution and sequencing to unveil the unknown antigen(s).
Western blot
Total cell extracts were prepared from SV40-fibroblasts from patients or healthy controls. Equal amounts of protein from each sample were separated by SDS-PAGE and blotted onto polyvinylidene difluoride membranes (Bio-Rad). The membranes were then probed with an anti-human DBR1 antibody (ProteinTech). They were then stripped and reprobed with an anti-GAPDH antibody (Santa-Cruz) to control for protein loading. Antibody binding was detected by enhanced chemiluminescence (ECL; Amersham-Pharmacia-Biotech) with an Amersham Imager 600 (GE Life Sciences).
Reverse transcription–quantitative PCR (RT-qPCR)
Total RNA was isolated from the patient’s SV-40 fibroblasts with the Quick-RNA Microprep kit (Zymo Research). We reverse-transcribed the extracted total RNA with random hexamers and the SuperScript III First-Strand Synthesis system (#18080051; Thermo Fisher Scientific). RT-qPCR was performed with Applied Biosystems TaqMan assays with Hs01113902_m1 (spanning DBR1 exons 2–3) and Hs01113907_m1 (exons 7–8) probes for DBR1 and with the β-glucuronidase (#4310888E; GUS) housekeeping gene used for normalization. Results were expressed according to the ΔΔCt method, performed in accordance with the manufacturer’s instructions.
For the quantification of ID1 and DKK1 intron lariats, we reverse-transcribed the extracted total RNA with 0.5 µM branch point-specific reverse primer for ID1 or DKK1 and a reverse primer for β-glucuronidase (GUS) for normalization. The ID1 and DKK1 mRNA transcripts were reverse-transcribed with random hexamers and quantified with primers for ID1 or DKK1 transcripts in the Fast SYBR Green System (#4385616; Thermo Fisher Scientific). The levels of ID1 and DKK1 lariats and transcripts were normalized relative to that of GUS transcripts and calculated according to the ∆CT method. The sequences of the primers used have been reported elsewhere (Zhang et al., 2018).
SIMOA digital ELISA
Pan-IFNα, IFNγ (duplex), and IFNβ (single-plex) protein concentrations were quantified in SIMOA digital ELISA assays developed as Quanterix Homebrews according to the manufacturer’s instructions. The limit of detection of these assays was 0.8 fg ml−1 for IFN-α, 20 fg ml−1 for IFN-γ, and 0.2 pg ml−1 for IFN-β, considering the dilution factor applied.
Fibroblast cell culture
Primary fibroblasts were isolated from skin punch biopsy specimens under sterile conditions and were cultured in DMEM (GIBCO BRL; Invitrogen) supplemented with 10% fetal calf serum (FCS) (GIBCO BRL; Invitrogen). Immortalized SV40-transformed fibroblast cell lines (SV40-fibroblasts) were created by the electroporation of about five million cells with 4 mg of a plasmid containing T-antigen DNA. The transfected cells were transferred to two fresh 75-cm2 flasks, each containing 12 ml DMEM (GIBCO BRL; Invitrogen) supplemented with 10% FCS (GIBCO BRL; Invitrogen). SV40-fibroblast clones appeared after about 15 days. These clones were cultured and passaged for experimental use.
Plasmids
The DBR1 (accession #Q9UK59) cDNA was inserted into the pDONOR vector. Site-directed mutagenesis was performed to obtain the mutant I120T DBR1 construct. WT DBR1 and I120T DBR1 constructs were then inserted into the pTRIP vector. To generate DKK1 lariat copGFP-split plasmid, the copGFP fragment was PCR-amplified from the PTY-copGFP plasmid and assembled with DKK1 into a PTY plasmid. DKK1 lariat copGFP constructs were then inserted into the pTRIP vector. All primers used for site-directed mutagenesis or subcloning were generated by SnapGene software (version 7). For lentiviral vector production, envelope plasmid pCMV-VSV-G, packaging plasmid PsPAX2, and transfer plasmid pTRIP were used. Lentivirus was concentrated with a Lenti-X concentrator (Takara Bio). All constructs were sequenced to ensure that no adventitious mutations were generated during the cloning process.
Patient-specific induced pluripotent stem cell (iPSC) reprogramming, culture, and characterization
Patient-specific iPSCs were obtained by reprogramming the patient’s primary fibroblasts by infection with the nonintegrating CytoTune Sendai viral vector kit (Life Technologies). All reprogrammed cells were karyotyped to ensure that they carried no chromosomal abnormalities. Patient-specific DBR1 mutations were confirmed by Sanger sequencing of genomic DNA extracted from the iPSC lines. Human iPSC cultures were maintained in Essential 8 medium (#A1517001; Life Technologies) on vitronectin-N (VTN-N, #A14700; Thermo Fisher Scientific)-coated plates. Healthy control hESC line H9 and two experimental control iPSC lines from a TLR3−/− patient and an IFNAR1−/− patient with established deficiencies of the TLR3-type-I IFN circuit (Bastard et al., 2021b; Guo et al., 2011) were used in this study.
Hindbrain neuron differentiation from hPSCs
The differentiation of hESCs or iPSCs (referred to jointly as hPSCs) into hindbrain neurons was induced by dual-SMAD inhibition (Chambers et al., 2009), with correct anterior-to-posterior patterning achieved by Wnt activation. Briefly, hPSCs were dissociated into a single-cell suspension with Accutase (#AT104; Innovative Cell Technologies) and used to seed Geltrex (#A1413202; Thermo Fisher Scientific)-coated plastic plates at a density of 250,000 cells/cm2 in E8 medium supplemented with 10 μM ROCK inhibitor (Y-27632; 10 µM #1254/10; R&D Systems). The cells were then transferred to neural induction medium for 11 days. This medium consisted of E6 medium supplemented with LDN193189 (100 nM #04-0074; Reprocell) and SB431542 (10 μM #1614/50; R&D Systems), with the addition of CHIR99021 (3 μM #4423; Tocris Bioscience) for the first 2 days. After 11 days, this medium was replaced by neural differentiation medium consisting of 1:1 DMEM/F12 and Neurobasal, 1× N2 supplement (#17502-048; Thermo Fisher Scientific), 1× B27 without vitamin A (#12587010; Thermo Fisher Scientific), and 1× penicillin/streptomycin (#15-140-122; Thermo Fisher Scientific). The cells were incubated in this neural differentiation medium for 4 days, and hindbrain neural progenitors were then either cryopreserved in Stem Cellbanker (#11924; Amsbio) or replated onto poly-ornithine/laminin/fibronectin plates at a density of 1.5 × 105 cells/cm2 in maturation medium consisting of Neurobasal, medium 1× B27 without vitamin A, 1× penicillin/streptomycin, 2 mM L-glutamine (#25030081; Thermo Fisher Scientific), dibutyryl cAMP (#D0627; Sigma-Aldrich), 10 μM DAPT (#2634; Tocris), 250 μM ascorbic acid (#A4034; Sigma-Aldrich), 10 ng/ml glial cell line–derived neurotrophic factor, and 10 ng/ml brain-derived neurotrophic factor. We added 10 μM Y-27632 at replating. The medium was replaced every 5 days until day 30, when the hindbrain neurons were used for experiments.
Virus propagation
The SARS-CoV-2 NYC isolate was obtained from the saliva of a deidentified patient on July 28, 2020. The sequence of the virus is publicly available (GenBank OM345241). The virus isolate was initially amplified in Caco-2 cells (passage 1, or P#1 stock). For the generation of P#2 and P#3 working stocks, Caco-2 cells were infected with the P#1 and P#2 viruses, respectively, at a MOI of 0.05 plaque-forming units (PFU)/cell and incubated for 6 and 7 days, respectively, at 37°C. The virus-containing supernatant was then harvested, clarified by centrifugation (3,000 × g for 10 min), and filtered through a disposable vacuum filter system with 0.22 mm pores. The P#3 stock used in this study had a titer of 3.4 × 106 PFU/ml, as determined on Vero E6 cells with a 1% methylcellulose overlay, as previously described (Mendoza et al., 2020).
Quantification of SARS-CoV-2 infection
Hindbrain neuron progenitors were used to seed 96-well plates at a density of 1.5 × 105 cells/cm2 and were differentiated into dorsal hindbrain neurons. After 4 wk in culture, the neurons were left untreated or were treated with IFN-β (1,000 IU/ml) for 18 h and were then infected with SARS-CoV-2 for 24, 48, 72, and 96 h. The cells were fixed with 10% neutral-buffered formalin and stained for SARS-CoV-2 nucleocapsid (#GTX135357; GeneTex) and neuron-specific cytoskeletal protein microtubule-associated protein 2 (MAP2, Ab11267; Abcam). Alexa Fluor 647– and Alexa Fluor 488–conjugated secondary antibodies (Invitrogen), respectively, were used for counterstaining. Plates were imaged with the ImageXpress micro XL High-Content Screening System, and data were analyzed with MetaXpress (Molecular Devices).
Statistical analysis
The data for hindbrain neuron infection with SARS-CoV-2 were obtained from three biological replicates in three independent experiments. For each biological replicate, three technical replicates were performed and averaged for downstream analyses. Statistical analysis was performed with Kruskal–Wallis tests with Dunn’s correction, and the results are indicated in the corresponding figures and legends (ns, not significant; *P < 0.05; **P < 0.01; ***P < 0.001).
Online supplemental material
Fig. S1 shows detailed VirScan, DBR1 Sanger sequencing, and levels of circulating type I and II IFNs of P1, his family members, and healthy controls enrolled in the study. Fig. S2 shows the characterization of hPSC-derived hindbrain neurons. Fig. S3 shows the evaluation of SARS-CoV-2 infection in hPSC-derived hindbrain neurons. Table S1 shows the homozygous or compound heterozygous rare nonsynonymous or essential splicing variants found in P1’s WGS data. Table S2 shows the viral PCR and viral antibody serological results of P1. Table S3 shows the leukocyte immunological functional tests for the patient. Table S4 shows the viral serological data from P1 and the family members of the patient.
Data availability
All data supporting the findings of this study are available in the published article and its online supplemental material or available from corresponding authors upon reasonable request in accordance with local regulations and ethical approvals related to studies of human subjects. The raw RNA sequencing (RNAseq) data generated from this study are deposited in the NCBI database under the NCBI-SRA project PRJNA1123312.
Acknowledgments
We thank the patients and their families for participating in our research. We thank all members of both branches of the Laboratory of Human Genetics of Infectious Diseases for discussions and technical and administrative support. The authors would like to acknowledge the support for next-generation sequencing and bioinformatics services provided by the Clinical Genomics Stockholm facility at the Science for Life Laboratory, hosted jointly by the Department of Microbiology, Tumor and Cell Biology at Karolinska Institutet and the School of Engineering Sciences in Chemistry, Biotechnology and Health at the Kungliga Tekniska Högskolan Royal Institute of Technology. We thank Ann-Charlotte Wikström from the Department of Clinical Immunology and Transfusion Medicine, Karolinska University Hospital, Stockholm and Valtteri Wirta, Science for Life Laboratory, Department of Microbiology, Tumor and Cell Biology, Karolinska Institutet, Stockholm; as well as the pediatric nurses Anita Jansson and Anna Höglund from the Department of Pediatrics, Härnösand Hospital, Region Västernorrland.
The Laboratory of Human Genetics of Infectious Diseases is supported by the Howard Hughes Medical Institute, the Rockefeller University, the St. Giles Foundation, the National Institutes of Health (NIH) (R01AI088364, R01AI163029, and R21AI160576), the National Center for Advancing Translational Sciences, the NIH Clinical and Translational Science Award program (UL1TR001866), the Yale Center for Mendelian Genomics and the Genome Sequencing Program Coordinating Center funded by the National Human Genome Research Institute (UM1HG006504 and U24HG008956), the Yale High-Performance Computing Center (S10OD018521), the Fisher Center for Alzheimer’s Research Foundation, the Meyer Foundation, the JPB Foundation, the French National Research Agency (ANR) under the “Investments for the Future” program (ANR-10-IAHU-01), the Integrative Biology of Emerging Infectious Diseases Laboratory of Excellence (ANR-10-LABX-62-IBEID), the French Foundation for Medical Research (FRM) (EQU201903007798), the ANR CNSVIRGEN (ANR-19-CE15-0009-01), the ANR SEAeHostFactors (ANR-18-CE15-0020-02), the ANR GenMISC (ANR-21-COVR-039), the ANRS-COV05, ANR GENVIR (ANR-20-CE93-003) and ANR AABIFNCOV (ANR-20-CO11-0001) projects, the ANR-Recherche Hospitalo-Universitaire en santé (RHU) program (ANR-21-RHUS-08), the European Union’s Horizon 2020 research and innovation program under grant agreement no. 824110 (EASI-genomics), the HORIZON-HLTH-2021-DISEASE-04 program under grant agreement 101057100 (UNDINE), the ANR-RHU Program ANR-21-RHUS-08 (COVIFERON), the Square Foundation, Grandir - Fonds de solidarité pour l’enfance, the Fondation du Souffle, the SCOR Corporate Foundation for Science, the Battersea & Bowery Advisory Group, the French Ministry of Higher Education, Research, and Innovation (MESRI-COVID-19), Institut National de la Santé et de la Recherche Médicale, and Paris Cité University. Y.-H. Chan was supported by an A*STAR International Fellowship. D. Lee was supported by a fellowship from the FRM for medical residents and fellows and the European Society for Immunodeficiencies 2023 Juniors Bridge Grant. Work in the Laboratory of Virology and Infectious Disease was supported by NIH grants P01AI138398-S1, 2U19AI111825, R01AI091707-10S1, and R01AI161444; a George Mason University Fast Grant; the G. Harold and Leila Y. Mathers Charitable Foundation; the Meyer Foundation; and the Bawd Foundation. J. Le Pen was supported by the Francois Wallace Monahan Postdoctoral Fellowship at The Rockefeller University and a European Molecular Biology Organization Long-Term Fellowship (ALTF 380-2018). I. Meyts is a senior clinical investigator at Fonds Wetenschappelijk Onderzoek (FWO) Vlaanderen and is supported by a KU Leuven C1 Grant (C16/18/007; severe infections), by FWO Grant G0B5120N and by the Jeffrey Modell Foundation. This work was supported in part by the Intramural Research Program of the National Institute of Allergy and Infectious Diseases, NIH (grant AI00122 and AI001270 to L.D. Notarangelo; grant AI001265-04 to H. Su, a member of the COVID Human Genetic Effort). O. Ekwall’s laboratory is supported by grants from The Swedish Research Council (2018-02752 and 2022-00781) and the Swedish state under an agreement between the Swedish government and the county councils, the ALF-agreement (ALFGBG-965795). Open Access funding provided by Rockefeller University.
Author contributions: Y.-H. Chan: Conceptualization, Formal analysis, Investigation, Methodology, Resources, Visualization, Writing—original draft, Writing—review & editing, V. Lundberg: Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Project administration, Resources, Validation, Visualization, Writing—original draft, Writing—review & editing, J. Le Pen: Funding acquisition, Investigation, Writing—review & editing, J. Yuan: Investigation, Methodology, D. Lee: Conceptualization, Investigation, Writing—review & editing, F. Pinci: Investigation, Methodology, S. Volpi: Resources, K. Nakajima: Formal analysis, Visualization, Writing—review & editing, V. Bondet: Formal analysis, Investigation, Methodology, Resources, Validation, Writing—review & editing, S. Akesson: Conceptualization, N.V. Khobrekar: Resources, A. Bodansky: Formal analysis, Investigation, Visualization, L. Du: Investigation, T. Melander: Resources, A.-A. Mariaggi: Resources, Y. Seeleuthner: Data curation, Resources, T.S. Saleh: Investigation, Resources, D. Chakravarty: Investigation, P. Marits: Investigation, Resources, Writing—review & editing, K. Dobbs: Resources, S. Vonlanthen: Resources, Writing—review & editing, V. Hennings: Investigation, K. Thorn: Investigation, D. Rinchai: Data curation, Formal analysis, Visualization, Writing—review & editing, L. Bizien: Investigation, M. Chaldebas: Data curation, Software, Writing—review & editing, A. Sobh: Conceptualization, Data curation, Investigation, Writing—review & editing, T. Ozcelik: Investigation, Resources, Writing—review & editing, S. Keles: Data curation, Validation, S.A. AlKhater: Resources, Writing—review & editing, C. Prando: Data curation, Investigation, Resources, Writing—review & editing, I. Meyts: Investigation, Resources, Writing—review & editing, M.R. Wilson: Formal analysis, Project administration, Writing—review & editing, J. Rosain: Investigation, E. Jouanguy: Resources, Writing—review & editing, M. Aubart: Conceptualization, Writing—review & editing, L. Abel: Formal analysis, Writing—review & editing, T.H. Mogensen: Conceptualization, Data curation, Funding acquisition, Writing—review & editing, Q. Pan-Hammarstrom: Supervision, D. Gao: Writing—review & editing, D. Duffy: Investigation, Writing—review & editing, A. Cobat: Formal analysis, S. Berg: Conceptualization, Methodology, Writing—review & editing, L.D. Notarangelo: Methodology, Resources, Writing—review & editing, O. Harschnitz: Conceptualization, Investigation, Methodology, Supervision, Writing—review & editing, C.M. Rice: Funding acquisition, Writing—review & editing, L. Studer: Conceptualization, Funding acquisition, Writing—review & editing, J.-L. Casanova: Conceptualization, Funding acquisition, Project administration, Resources, Supervision, Validation, Writing—original draft, Writing—review & editing, O. Ekwall: Conceptualization, Funding acquisition, Resources, Supervision, Writing—original draft, Writing—review & editing, S.-Y. Zhang: Conceptualization, Funding acquisition, Investigation, Project administration, Resources, Supervision, Validation, Visualization, Writing—original draft, Writing—review & editing.
References
COVID Human Genetic Effort collaborators: Paul Bastard, Alessandro Borghesi, Aziz Bousfiha, Oksana Boyarchuk, Petter Brodin, Jacinta Bustamante, Giorgio Casari, Rémi Chevalier, John Christodoulou, Roger Colobran, Antonio Condino-Neto, Juan Carlos Aldave Becerra, Lisa Arkin, Evangelos Andreakos, Christian W. Thorball, Sara Espinosa, Carlos Flores, Amyrath Geraldo, Rabih Halwani, Nevin Hatipoğlu, Brahim Melaiki, Jacques Fellay, Alenka Gagro, Yuval Itan, Chandima Jeewandara, Filomeen Haerynck, Davood Mansouri, Leslie Naesens, Lisa F.P. Ng, Keisuke Okamoto, Pere Soler-Palacin, Laurent Renia, Aurora Pujol Onofre, Igor B. Resnick, José Luis Franco Restrepo, Jacques G. Rivière, Anna Scherbina, Anna Šedivá, Mikko R.J. Seppänen, Helen Su, Stuart G. Tangye, Sehime G. Temel, Ahmad Abou Tayoun, Stuart Turvey, K.M. Furkan Uddin, Diederik van de Beek, Tom Le Voyer, Donald C. Vinh, Shen-Ying Zhang, and Jean-Laurent Casanova. Affiliations are listed at the end of the PDF.
Author notes
Y.-H. Chan, V. Lundberg, J. Le Pen, and J. Yuan contributed equally to this paper.
D. Lee, F. Pinci, and S. Volpi contributed equally to this paper.
O. Harschnitz, C.M. Rice, L. Studer, J.-L. Casanova, O. Ekwall, and S.-Y. Zhang contributed equally to this paper.
Disclosures: F. Pinci and O. Harschnitz reported a patent to differentiation of hindbrain neurons from human pluripotent stem cells pending. I. Meyts reported grants from CSL-Behring, FWO Vlaanderen G0B5120N, KU Leuven, and Jeffrey Modell Foundation during the conduct of the study; and “other” from Takeda and Boehringer-Ingelheim outside the submitted work. M.R Wilson reported grants from Genentech and Novartis, and personal fees from Delve Bio outside the submitted work. S. Berg reported personal fees from speaker’s bureau, SOBI outside the submitted work. No other disclosures were reported.