A variant in human AIOLOS impairs adaptive immunity by interfering with IKAROS.

In the present study, we report a human-inherited, impaired, adaptive immunity disorder, which predominantly manifested as a B cell differentiation defect, caused by a heterozygous IKZF3 missense variant, resulting in a glycine-to-arginine replacement within the DNA-binding domain of the encoded AIOLOS protein. Using mice that bear the corresponding variant and recapitulate the B and T cell phenotypes, we show that the mutant AIOLOS homodimers and AIOLOS-IKAROS heterodimers did not bind the canonical AIOLOS-IKAROS DNA sequence. In addition, homodimers and heterodimers containing one mutant AIOLOS bound to genomic regions lacking both canonical motifs. However, the removal of the dimerization capacity from mutant AIOLOS restored B cell development. Hence, the adaptive immunity defect is caused by the AIOLOS variant hijacking IKAROS function. Heterodimeric interference is a new mechanism of autosomal dominance that causes inborn errors of immunity by impairing protein function via the mutation of its heterodimeric partner.

An essential role for the Zn2+ transporter ZIP7 in B cell development.

Despite the known importance of zinc for human immunity, molecular insights into its roles have remained limited. Here we report a novel autosomal recessive disease characterized by absent B cells, agammaglobulinemia and early onset infections in five unrelated families. The immunodeficiency results from hypomorphic mutations of SLC39A7, which encodes the endoplasmic reticulum-to-cytoplasm zinc transporter ZIP7. Using CRISPR-Cas9 mutagenesis we have precisely modeled ZIP7 deficiency in mice. Homozygosity for a null allele caused embryonic death, but hypomorphic alleles reproduced the block in B cell development seen in patients. B cells from mutant mice exhibited a diminished concentration of cytoplasmic free zinc, increased phosphatase activity and decreased phosphorylation of signaling molecules downstream of the pre-B cell and B cell receptors. Our findings highlight a specific role for cytosolic Zn2+ in modulating B cell receptor signal strength and positive selection.

Genome-wide detection of human intronic AG-gain variants located between splicing branchpoints and canonical splice acceptor sites.

Human genetic variants that introduce an AG into the intronic region between the branchpoint (BP) and the canonical splice acceptor site (ACC) of protein-coding genes can disrupt pre-mRNA splicing. Using our genome-wide BP database, we delineated the BP-ACC segments of all human introns and found extreme depletion of AG/YAG in the [BP+8, ACC-4] high-risk region. We developed AGAIN as a genome-wide computational approach to systematically and precisely pinpoint intronic AG-gain variants within the BP-ACC regions. AGAIN identified 350 AG-gain variants from the Human Gene Mutation Database, all of which alter splicing and cause disease. Among them, 74% created new acceptor sites, whereas 31% resulted in complete exon skipping. AGAIN also predicts the protein-level products resulting from these two consequences. We performed AGAIN on our exome/genomes database of patients with severe infectious diseases but without known genetic etiology and identified a private homozygous intronic AG-gain variant in the antimycobacterial gene SPPL2A in a patient with mycobacterial disease. AGAIN also predicts a retention of six intronic nucleotides that encode an in-frame stop codon, turning AG-gain into stop-gain. This allele was then confirmed experimentally to lead to loss of function by disrupting splicing. We further showed that AG-gain variants inside the high-risk region led to misspliced products, while those outside the region did not, by two case studies in genes STAT1 and IRF7. We finally evaluated AGAIN on our 14 paired exome-RNAseq samples and found that 82% of AG-gain variants in high-risk regions showed evidence of missplicing. AGAIN is publicly available from https://hgidsoft.rockefeller.edu/AGAIN and https://github.com/casanova-lab/AGAIN.

Genome-wide detection of human variants that disrupt intronic branchpoints.

Pre-messenger RNA splicing is initiated with the recognition of a single-nucleotide intronic branchpoint (BP) within a BP motif by spliceosome elements. Forty-eight rare variants in 43 human genes have been reported to alter splicing and cause disease by disrupting BP. However, until now, no computational approach was available to efficiently detect such variants in massively parallel sequencing data. We established a comprehensive human genome-wide BP database by integrating existing BP data and generating new BP data from RNA sequencing of lariat debranching enzyme DBR1-mutated patients and from machine-learning predictions. We characterized multiple features of BP in major and minor introns and found that BP and BP-2 (two nucleotides upstream of BP) positions exhibit a lower rate of variation in human populations and higher evolutionary conservation than the intronic background, while being comparable to the exonic background. We developed BPHunter as a genome-wide computational approach to systematically and efficiently detect intronic variants that may disrupt BP recognition. BPHunter retrospectively identified 40 of the 48 known pathogenic BP variants, in which we summarized a strategy for prioritizing BP variant candidates. The remaining eight variants all create AG-dinucleotides between the BP and acceptor site, which is the likely reason for missplicing. We demonstrated the practical utility of BPHunter prospectively by using it to identify a novel germline heterozygous BP variant of STAT2 in a patient with critical COVID-19 pneumonia and a novel somatic intronic 59-nucleotide deletion of ITPKB in a lymphoma patient, both of which were validated experimentally. BPHunter is publicly available from https://hgidsoft.rockefeller.edu/BPHunter and https://github.com/casanova-lab/BPHunter.

Loss of C2orf69 defines a fatal autoinflammatory syndrome in humans and zebrafish that evokes a glycogen-storage-associated mitochondriopathy.

Human C2orf69 is an evolutionarily conserved gene whose function is unknown. Here, we report eight unrelated families from which 20 children presented with a fatal syndrome consisting of severe autoinflammation and progredient leukoencephalopathy with recurrent seizures; 12 of these subjects, whose DNA was available, segregated homozygous loss-of-function C2orf69 variants. C2ORF69 bears homology to esterase enzymes, and orthologs can be found in most eukaryotic genomes, including that of unicellular phytoplankton. We found that endogenous C2ORF69 (1) is loosely bound to mitochondria, (2) affects mitochondrial membrane potential and oxidative respiration in cultured neurons, and (3) controls the levels of the glycogen branching enzyme 1 (GBE1) consistent with a glycogen-storage-associated mitochondriopathy. We show that CRISPR-Cas9-mediated inactivation of zebrafish C2orf69 results in lethality by 8 months of age due to spontaneous epileptic seizures, which is preceded by persistent brain inflammation. Collectively, our results delineate an autoinflammatory Mendelian disorder of C2orf69 deficiency that disrupts the development/homeostasis of the immune and central nervous systems.

Immunodeficiency, autoinflammation and amylopectinosis in humans with inherited HOIL-1 and LUBAC deficiency.

We report the clinical description and molecular dissection of a new fatal human inherited disorder characterized by chronic autoinflammation, invasive bacterial infections and muscular amylopectinosis. Patients from two kindreds carried biallelic loss-of-expression and loss-of-function mutations in HOIL1 (RBCK1), a component of the linear ubiquitination chain assembly complex (LUBAC). These mutations resulted in impairment of LUBAC stability. NF-κB activation in response to interleukin 1ß (IL-1ß) was compromised in the patients' fibroblasts. By contrast, the patients' mononuclear leukocytes, particularly monocytes, were hyper-responsive to IL-1ß. The consequences of human HOIL-1 and LUBAC deficiencies for IL-1ß responses thus differed between cell types, consistent with the unique association of autoinflammation and immunodeficiency in these patients. These data suggest that LUBAC regulates NF-κB-dependent IL-1ß responses differently in different cell types.

Negative selection on human genes underlying inborn errors depends on disease outcome and both the mode and mechanism of inheritance.

Genetic variants underlying life-threatening diseases, being unlikely to be transmitted to the next generation, are gradually and selectively eliminated from the population through negative selection. We study the determinants of this evolutionary process in human genes underlying monogenic diseases by comparing various negative selection scores and an integrative approach, CoNeS, at 366 loci underlying inborn errors of immunity (IEI). We find that genes underlying autosomal dominant (AD) or X-linked IEI have stronger negative selection scores than those underlying autosomal recessive (AR) IEI, whose scores are not different from those of genes not known to be disease causing. Nevertheless, genes underlying AR IEI that are lethal before reproductive maturity with complete penetrance have stronger negative selection scores than other genes underlying AR IEI. We also show that genes underlying AD IEI by loss of function have stronger negative selection scores than genes underlying AD IEI by gain of function, while genes underlying AD IEI by haploinsufficiency are under stronger negative selection than other genes underlying AD IEI. These results are replicated in 1,140 genes underlying inborn errors of neurodevelopment. Finally, we propose a supervised classifier, SCoNeS, which predicts better than state-of-the-art approaches whether a gene is more likely to underlie an AD or AR disease. The clinical outcomes of monogenic inborn errors, together with their mode and mechanisms of inheritance, determine the levels of negative selection at their corresponding loci. Integrating scores of negative selection may facilitate the prioritization of candidate genes and variants in patients suspected to carry an inborn error.

Multibatch Cytometry Data Integration for Optimal Immunophenotyping.

High-dimensional cytometry is a powerful technique for deciphering the immunopathological factors common to multiple individuals. However, rational comparisons of multiple batches of experiments performed on different occasions or at different sites are challenging because of batch effects. In this study, we describe the integration of multibatch cytometry datasets (iMUBAC), a flexible, scalable, and robust computational framework for unsupervised cell-type identification across multiple batches of high-dimensional cytometry datasets, even without technical replicates. After overlaying cells from multiple healthy controls across batches, iMUBAC learns batch-specific cell-type classification boundaries and identifies aberrant immunophenotypes in patient samples from multiple batches in a unified manner. We illustrate unbiased and streamlined immunophenotyping using both public and in-house mass cytometry and spectral flow cytometry datasets. The method is available as the R package iMUBAC (https://github.com/casanova-lab/iMUBAC).

X-Linked TLR7 Deficiency Underlies Critical COVID-19 Pneumonia in a Male Patient with Ataxia-Telangiectasia.

BACKGROUND: Coronavirus disease 2019 (COVID-19) exhibits a wide spectrum of clinical manifestations, ranging from asymptomatic to critical conditions. Understanding the mechanism underlying life-threatening COVID-19 is instrumental for disease prevention and treatment in individuals with a high risk. OBJECTIVES: We aimed to identify the genetic cause for critical COVID-19 pneumonia in a patient with a preexisting inborn error of immunity (IEI). METHODS: Serum levels of specific antibodies against the virus and autoantibodies against type I interferons (IFNs) were measured. Whole exome sequencing was performed, and the impacts of candidate gene variants were investigated. We also evaluated 247 ataxia-telangiectasia (A-T) patients in the Iranian IEI registry. RESULTS: We report a 7-year-old Iranian boy with a preexisting hyper IgM syndrome who developed critical COVID-19 pneumonia. IgM only specific COVID-19 immune response was detected but no autoantibodies against type I IFN were observed. A homozygous deleterious mutation in the ATM gene was identified, which together with his antibody deficiency, radiosensitivity, and neurological signs, established a diagnosis of A-T. Among the 247 A-T patients evaluated, 36 had SARS-CoV-2 infection, but all had mild symptoms or were asymptomatic except the index patient. A hemizygous deleterious mutation in the TLR7 gene was subsequently identified in the patient. CONCLUSIONS: We report a unique IEI patient with combined ATM and TLR7 deficiencies. The two genetic defects underlie A-T and critical COVID-19 in this patient, respectively.

SCID and Other Inborn Errors of Immunity with Low TRECs - the Brazilian Experience.

Severe combined immunodeficiency, SCID, is a pediatric emergency that represents the most critical group of inborn errors of immunity (IEI). Affected infants present with early onset life-threatening infections due to absent or non-functional T cells. Without early diagnosis and curative treatment, most die in early infancy. As most affected infants appear healthy at birth, newborn screening (NBS) is essential to identify and treat patients before the onset of symptoms. Here, we report 47 Brazilian patients investigated between 2009 and 2020 for SCID due to either a positive family history and/or clinical impression and low TRECs. Based on clinical presentation, laboratory finding, and genetic information, 24 patients were diagnosed as typical SCID, 14 as leaky SCID, and 6 as Omenn syndrome; 2 patients had non-SCID IEI, and 1 remained undefined. Disease onset median age was 2 months, but at the time of diagnosis and treatment, median ages were 6.5 and 11.5 months, respectively, revealing considerable delay which affected negatively treatment success. While overall survival was 51.1%, only 66.7% (30/45) lived long enough to undergo hematopoietic stem-cell transplantation, which was successful in 70% of cases. Forty-three of 47 (91.5%) patients underwent genetic testing, with a 65.1% success rate. Even though our patients did not come from the NBS programs, the diagnosis of SCID improved in Brazil during the pilot programs, likely due to improved medical education. However, we estimate that at least 80% of SCID cases are still missed. NBS-SCID started to be universally implemented in the city of São Paulo in May 2021, and it is our hope that other cities will follow, leading to early diagnosis and higher survival of SCID patients in Brazil.

An ACT1 mutation selectively abolishes interleukin-17 responses in humans with chronic mucocutaneous candidiasis.

Patients with inborn errors of interleukin-17F (IL-17F) or IL-17RA display chronic mucocutaneous candidiasis (CMC). We report a biallelic missense mutation (T536I) in the adaptor molecule ACT1 in two siblings with CMC. The mutation, located in the SEFIR domain, abolished the homotypic interaction of ACT1 with IL-17 receptors, with no effect on homodimerization. The patients' fibroblasts failed to respond to IL-17A and IL-17F, and their T cells to IL-17E. By contrast, healthy individuals homozygous for the common variant D10N, located in the ACT1 tumor necrosis factor receptor-associated factor-interacting domain and previously associated with psoriasis, had impaired, but not abolished, responses to IL-17 cytokines. SEFIR-independent interactions of ACT1 with other proteins, such as CD40, heat shock protein 70 (HSP70) and HSP90, were not affected by the T536I mutation. Overall, human IL-17A and IL-17F depend on ACT1 to mediate protective mucocutaneous immunity. Moreover, other ACT1-dependent IL-17 cytokines seem to be largely redundant in host defense.

Blacklisting variants common in private cohorts but not in public databases optimizes human exome analysis.

Computational analyses of human patient exomes aim to filter out as many nonpathogenic genetic variants (NPVs) as possible, without removing the true disease-causing mutations. This involves comparing the patient's exome with public databases to remove reported variants inconsistent with disease prevalence, mode of inheritance, or clinical penetrance. However, variants frequent in a given exome cohort, but absent or rare in public databases, have also been reported and treated as NPVs, without rigorous exploration. We report the generation of a blacklist of variants frequent within an in-house cohort of 3,104 exomes. This blacklist did not remove known pathogenic mutations from the exomes of 129 patients and decreased the number of NPVs remaining in the 3,104 individual exomes by a median of 62%. We validated this approach by testing three other independent cohorts of 400, 902, and 3,869 exomes. The blacklist generated from any given cohort removed a substantial proportion of NPVs (11-65%). We analyzed the blacklisted variants computationally and experimentally. Most of the blacklisted variants corresponded to false signals generated by incomplete reference genome assembly, location in low-complexity regions, bioinformatic misprocessing, or limitations inherent to cohort-specific private alleles (e.g., due to sequencing kits, and genetic ancestries). Finally, we provide our precalculated blacklists, together with ReFiNE, a program for generating customized blacklists from any medium-sized or large in-house cohort of exome (or other next-generation sequencing) data via a user-friendly public web server. This work demonstrates the power of extracting variant blacklists from private databases as a specific in-house but broadly applicable tool for optimizing exome analysis.

A deep intronic splice mutation of STAT3 underlies hyper IgE syndrome by negative dominance.

Heterozygous in-frame mutations in coding regions of human STAT3 underlie the only known autosomal dominant form of hyper IgE syndrome (AD HIES). About 5% of familial cases remain unexplained. The mutant proteins are loss-of-function and dominant-negative when tested following overproduction in recipient cells. However, the production of mutant proteins has not been detected and quantified in the cells of heterozygous patients. We report a deep intronic heterozygous STAT3 mutation, c.1282-89C>T, in 7 relatives with AD HIES. This mutation creates a new exon in the STAT3 complementary DNA, which, when overexpressed, generates a mutant STAT3 protein (D427ins17) that is loss-of-function and dominant-negative in terms of tyrosine phosphorylation, DNA binding, and transcriptional activity. In immortalized B cells from these patients, the D427ins17 protein was 2 kDa larger and 4-fold less abundant than wild-type STAT3, on mass spectrometry. The patients' primary B and T lymphocytes responded poorly to STAT3-dependent cytokines. These findings are reminiscent of the impaired responses of leukocytes from other patients with AD HIES due to typical STAT3 coding mutations, providing further evidence for the dominance of the mutant intronic allele. These findings highlight the importance of sequencing STAT3 introns in patients with HIES without candidate variants in coding regions and essential splice sites. They also show that AD HIES-causing STAT3 mutant alleles can be dominant-negative even if the encoded protein is produced in significantly smaller amounts than wild-type STAT3.

Human intracellular ISG15 prevents interferon-α/ß over-amplification and auto-inflammation.

Intracellular ISG15 is an interferon (IFN)-α/ß-inducible ubiquitin-like modifier which can covalently bind other proteins in a process called ISGylation; it is an effector of IFN-α/ß-dependent antiviral immunity in mice. We previously published a study describing humans with inherited ISG15 deficiency but without unusually severe viral diseases. We showed that these patients were prone to mycobacterial disease and that human ISG15 was non-redundant as an extracellular IFN-γ-inducing molecule. We show here that ISG15-deficient patients also display unanticipated cellular, immunological and clinical signs of enhanced IFN-α/ß immunity, reminiscent of the Mendelian autoinflammatory interferonopathies Aicardi-Goutières syndrome and spondyloenchondrodysplasia. We further show that an absence of intracellular ISG15 in the patients' cells prevents the accumulation of USP18, a potent negative regulator of IFN-α/ß signalling, resulting in the enhancement and amplification of IFN-α/ß responses. Human ISG15, therefore, is not only redundant for antiviral immunity, but is a key negative regulator of IFN-α/ß immunity. In humans, intracellular ISG15 is IFN-α/ß-inducible not to serve as a substrate for ISGylation-dependent antiviral immunity, but to ensure USP18-dependent regulation of IFN-α/ß and prevention of IFN-α/ß-dependent autoinflammation.

SeqTailor: a user-friendly webserver for the extraction of DNA or protein sequences from next-generation sequencing data.

Human whole-genome-sequencing reveals about 4 000 000 genomic variants per individual. These data are mostly stored as VCF-format files. Although many variant analysis methods accept VCF as input, many other tools require DNA or protein sequences, particularly for splicing prediction, sequence alignment, phylogenetic analysis, and structure prediction. However, there is no existing webserver capable of extracting DNA/protein sequences for genomic variants from VCF files in a user-friendly and efficient manner. We developed the SeqTailor webserver to bridge this gap, by enabling rapid extraction of (i) DNA sequences around genomic variants, with customizable window sizes and options to annotate the splice sites closest to the variants and to consider the neighboring variants within the window; and (ii) protein sequences encoded by the DNA sequences around genomic variants, with built-in SnpEff annotator and customizable window sizes. SeqTailor supports 11 species, including: human (GRCh37/GRCh38), chimpanzee, mouse, rat, cow, chicken, lizard, zebrafish, fruitfly, Arabidopsis and rice. Standalone programs are provided for command-line-based needs. SeqTailor streamlines the sequence extraction process, and accelerates the analysis of genomic variants with software requiring DNA/protein sequences. It will facilitate the study of genomic variation, by increasing the feasibility of sequence-based analysis and prediction. The SeqTailor webserver is freely available at http://shiva.rockefeller.edu/SeqTailor/.

Incomplete penetrance for isolated congenital asplenia in humans with mutations in translated and untranslated RPSA exons.

Isolated congenital asplenia (ICA) is the only known human developmental defect exclusively affecting a lymphoid organ. In 2013, we showed that private deleterious mutations in the protein-coding region of RPSA, encoding ribosomal protein SA, caused ICA by haploinsufficiency with complete penetrance. We reported seven heterozygous protein-coding mutations in 8 of the 23 kindreds studied, including 6 of the 8 multiplex kindreds. We have since enrolled 33 new kindreds, 5 of which are multiplex. We describe here 11 new heterozygous ICA-causing RPSA protein-coding mutations, and the first two mutations in the 5'-UTR of this gene, which disrupt mRNA splicing. Overall, 40 of the 73 ICA patients (55%) and 23 of the 56 kindreds (41%) carry mutations located in translated or untranslated exons of RPSA. Eleven of the 43 kindreds affected by sporadic disease (26%) carry RPSA mutations, whereas 12 of the 13 multiplex kindreds (92%) carry RPSA mutations. We also report that 6 of 18 (33%) protein-coding mutations and the two (100%) 5'-UTR mutations display incomplete penetrance. Three mutations were identified in two independent kindreds, due to a hotspot or a founder effect. Finally, RPSA ICA-causing mutations were demonstrated to be de novo in 7 of the 23 probands. Mutations in RPSA exons can affect the translated or untranslated regions and can underlie ICA with complete or incomplete penetrance.

The genetic basis of pneumococcal and staphylococcal infections: inborn errors of human TLR and IL-1R immunity.

Many bacteria can cause pyogenic lesions in humans. Most of these bacteria are harmless in most individuals, but they, nevertheless, cause significant morbidity and mortality worldwide. The inherited and acquired immunodeficiencies underlying these pyogenic infections differ between bacteria. This short review focuses on two emblematic pyogenic bacteria: pneumococcus (Streptococcus pneumoniae) and Staphylococcus, both of which are Gram-positive encapsulated bacteria. We will discuss the contribution of human genetic studies to the identification of germline mutations of the TLR and IL-1R pathways.

From Dysgammaglobulinemia to Autosomal-Dominant Activation-Induced Cytidine Deaminase Deficiency: Unraveling an Inherited Immunodeficiency after 50 Years.

The genetic investigation of a family presenting with a dominant form of hyper IgM syndrome published in 1963 and 1975 revealed a R190X nonsense mutation in activation-induced cytidine deaminase. This report illustrates the progress made over 6 decades in the characterization of primary immunodeficiencies, from immunochemistry to whole-exome sequencing.

T-cell defects in patients with ARPC1B germline mutations account for combined immunodeficiency.

ARPC1B is a key factor for the assembly and maintenance of the ARP2/3 complex that is involved in actin branching from an existing filament. Germline biallelic mutations in ARPC1B have been recently described in 6 patients with clinical features of combined immunodeficiency (CID), whose neutrophils and platelets but not T lymphocytes were studied. We hypothesized that ARPC1B deficiency may also lead to cytoskeleton and functional defects in T cells. We have identified biallelic mutations in ARPC1B in 6 unrelated patients with early onset disease characterized by severe infections, autoimmune manifestations, and thrombocytopenia. Immunological features included T-cell lymphopenia, low numbers of naïve T cells, and hyper-immunoglobulin E. Alteration in ARPC1B protein structure led to absent/low expression by flow cytometry and confocal microscopy. This molecular defect was associated with the inability of patient-derived T cells to extend an actin-rich lamellipodia upon T-cell receptor (TCR) stimulation and to assemble an immunological synapse. ARPC1B-deficient T cells additionally displayed impaired TCR-mediated proliferation and SDF1-α-directed migration. Gene transfer of ARPC1B in patients' T cells using a lentiviral vector restored both ARPC1B expression and T-cell proliferation in vitro. In 2 of the patients, in vivo somatic reversion restored ARPC1B expression in a fraction of lymphocytes and was associated with a skewed TCR repertoire. In 1 revertant patient, memory CD8+ T cells expressing normal levels of ARPC1B displayed improved T-cell migration. Inherited ARPC1B deficiency therefore alters T-cell cytoskeletal dynamics and functions, contributing to the clinical features of CID.