Search

Since the beginning of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2)/coronavirus disease 2019 (COVID-19) pandemic, global sequencing efforts have led in the field of inborn errors of immunity, and inspired particularly by previous research on life-threatening influenza, they have revealed that known and novel inborn errors affecting type I interferon immunity underlie critical COVID-19 in up to 5% of cases. In addition, neutralizing autoantibodies against type I interferons have been identified in up to 20% of patients with critical COVID-19 who are older than 80 years and 20% of fatal cases, with a higher prevalence in men and individuals older than 70 years. Also, inborn errors impairing regulation of type I interferon responses and RNA degradation have been found as causes of multisystem inflammatory syndrome in children, a life-threatening hyperinflammatory condition complicating otherwise mild initial SARS-CoV-2 infection in children and young adults. Better understanding of these immunologic mechanisms can aid in designing treatments for severe COVID-19, multisystem inflammatory syndrome in children, long COVID, and neuro-COVID.

Epilepsy (EP) and congenital heart disease (CHD) are two apparently unrelated diseases that nevertheless display substantial mutual comorbidity. Thus, whilst congenital heart defects are associated with an elevated risk of developing epilepsy, the incidence of epilepsy in CHD patients correlates with CHD severity. Although genetic determinants have been postulated to underlie the comorbidity of EP and CHD, the precise genetic etiology is unknown. We performed variant and gene association analyses on EP and CHD patients separately, using whole exomes of genetically identified Europeans from the UK Biobank and Mount Sinai BioMe Biobank. We prioritized biologically plausible candidate genes and investigated the enriched pathways and other identified comorbidities by biological proximity calculation, pathway analyses, and gene-level phenome-wide association studies. Our variant- and gene-level results point to the Voltage-Gated Calcium Channels (VGCC) pathway being a unifying framework for EP and CHD comorbidity. Additionally, pathway-level analyses indicated that the functions of disease-associated genes partially overlap between the two disease entities. Lastly, phenome-wide association analyses of prioritized candidate genes revealed that cerebral blood flow and ulcerative colitis constitute the two main traits associated with both EP and CHD.

Summary Identifying whether a given genetic mutation results in a gene product with increased (gain-of-function; GOF) or diminished (loss-of-function; LOF) activity is an important step toward understanding disease mechanisms because they may result in markedly different clinical phenotypes. Here, we generated an extensive database of documented germline GOF and LOF pathogenic variants by employing natural language processing (NLP) on the available abstracts in the Human Gene Mutation Database. We then investigated various gene- and protein-level features of GOF and LOF variants and applied machine learning and statistical analyses to identify discriminative features. We found that GOF variants were enriched in essential genes, for autosomal-dominant inheritance, and in protein binding and interaction domains, whereas LOF variants were enriched in singleton genes, for protein-truncating variants, and in protein core regions. We developed a user-friendly web-based interface that enables the extraction of selected subsets from the GOF/LOF database by a broad set of annotated features and downloading of up-to-date versions. These results improve our understanding of how variants affect gene/protein function and may ultimately guide future treatment options.

Since the arrival of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in December 2019, its characterization as a novel human pathogen, and the resulting coronavirus disease 2019 (COVID-19) pandemic, over 6.5 million people have died worldwide—a stark and sobering reminder of the fundamental and nonredundant roles of the innate and adaptive immune systems in host defense against emerging pathogens. Inborn errors of immunity (IEI) are caused by germline variants, typically in single genes. IEI are characterized by defects in development and/or function of cells involved in immunity and host defense, rendering individuals highly susceptible to severe, recurrent, and sometimes fatal infections, as well as immune dysregulatory conditions such as autoinflammation, autoimmunity, and allergy. The study of IEI has revealed key insights into the molecular and cellular requirements for immune-mediated protection against infectious diseases. Indeed, this has been exemplified by assessing the impact of SARS-CoV-2 infection in individuals with previously diagnosed IEI, as well as analyzing rare cases of severe COVID-19 in otherwise healthy individuals. This approach has defined fundamental aspects of mechanisms of disease pathogenesis, immunopathology in the context of infection with a novel pathogen, and therapeutic options to mitigate severe disease. This review summarizes these findings and illustrates how the study of these rare experiments of nature can inform key features of human immunology, which can then be leveraged to improve therapies for treating emerging and established infectious diseases.

The human genetic dissection of clinical phenotypes is complicated by genetic heterogeneity. Gene burden approaches that detect genetic signals in case-control studies are underpowered in genetically heterogeneous cohorts. We therefore developed a genome-wide computational method, network-based heterogeneity clustering (NHC), to detect physiological homogeneity in the midst of genetic heterogeneity. Simulation studies showed our method to be capable of systematically converging genes in biological proximity on the background biological interaction network, and capturing gene clusters harboring presumably deleterious variants, in an efficient and unbiased manner. We applied NHC to whole-exome sequencing data from a cohort of 122 individuals with herpes simplex encephalitis (HSE), including 13 individuals with previously published monogenic inborn errors of TLR3-dependent IFN-α/β immunity. The top gene cluster identified by our approach successfully detected and prioritized all causal variants of five TLR3 pathway genes in the 13 previously reported individuals. This approach also suggested candidate variants of three reported genes and four candidate genes from the same pathway in another ten previously unstudied individuals. TLR3 responsiveness was impaired in dermal fibroblasts from four of the five individuals tested, suggesting that the variants detected were causal for HSE. NHC is, therefore, an effective and unbiased approach for unraveling genetic heterogeneity by detecting physiological homogeneity.

Host genetic susceptibility is a key risk factor for severe illness associated with COVID-19. Despite numerous studies of COVID-19 host genetics, our knowledge of COVID-19-associated variants is still limited, and there is no resource comprising all the published variants and categorizing them based on their confidence level. Also, there are currently no computational tools available to predict novel COVID-19 severity variants. Therefore, we collated 820 host genetic variants reported to affect COVID-19 susceptibility by means of a systematic literature search and confidence evaluation, and obtained 196 high-confidence variants. We then developed the first machine learning classifier of severe COVID-19 variants to perform a genome-wide prediction of COVID-19 severity for 82,468,698 missense variants in the human genome. We further evaluated the classifier’s predictions using feature importance analyses to investigate the biological properties of COVID-19 susceptibility variants, which identified conservation scores as the most impactful predictive features. The results of enrichment analyses revealed that genes carrying high-confidence COVID-19 susceptibility variants shared pathways, networks, diseases and biological functions, with the immune system and infectious disease being the most significant categories. Additionally, we investigated the pleiotropic effects of COVID-19-associated variants using phenome-wide association studies (PheWAS) in ∼40,000 BioMe BioBank genotyped individuals, revealing pre-existing conditions that could serve to increase the risk of severe COVID-19 such as chronic liver disease and thromboembolism. Lastly, we generated a web-based interface for exploring, downloading and submitting genetic variants associated with COVID-19 susceptibility for use in both research and clinical settings (https://itanlab.shinyapps.io/COVID19webpage/). Taken together, our work provides the most comprehensive COVID-19 host genetics knowledgebase to date for the known and predicted genetic determinants of severe COVID-19, a resource that should further contribute to our understanding of the biology underlying COVID-19 susceptibility and facilitate the identification of individuals at high risk for severe COVID-19.

Background & Aims Recent literature has implicated a key role for mast cells in murine models of colonic inflammation, but their role in human ulcerative colitis (UC) is not well established. A major advance has been the identification of mrgprb2 (human orthologue, MRGPX2) as mediating IgE-independent mast cell activation. We sought to define mechanisms of mast cell activation and MRGPRX2 in human UC. Methods Colon tissues were collected from patients with UC for bulk RNA sequencing and lamina propria cells were isolated for MRGPRX2 activation studies and single-cell RNA sequencing. Genetic association of all protein-altering G-protein coupled receptor single-nucleotide polymorphism was performed in an Ashkenazi Jewish UC case-control cohort. Variants of MRGPRX2 were transfected into Chinese hamster ovary (CHO) and human mast cell (HMC) 1.1 cells to detect genotype-dependent effects on β-arrestin recruitment, IP-1 accumulation, and phosphorylated extracellular signal–regulated kinase. Results Mast cell-specific mediators and adrenomedullin (proteolytic precursor of PAMP-12, an MRGPRX2 agonist) are up-regulated in inflamed compared to uninflamed UC. MRGPRX2 stimulation induces carboxypeptidase secretion from inflamed UC. Of all protein-altering GPCR alleles, a unique variant of MRGPRX2, Asn62Ser, was most associated with and was bioinformatically predicted to alter arrestin recruitment. We validated that the UC protective serine allele enhances β-arrestin recruitment, decreases IP-1, and increases phosphorylated extracellular signal–regulated kinase with MRGPRX2 agonists. Single-cell RNA sequencing defines that adrenomedullin is expressed by activated fibroblasts and epithelial cells and that interferon gamma is a key upstream regulator of mast cell gene expression. Conclusion Inflamed UC regions are distinguished by MRGPRX2-mediated activation of mast cells, with decreased activation observed with a UC-protective genetic variant. These results define cell modules of UC activation and a new therapeutic target.

Bileaflet mitral valve prolapse (Bi-MVP) is associated with increased risk for cardiac arrest. We describe a patient who presented after a cardiac arrest with Bi-MVP and variants in Lamin A/C (LMNA) and the sodium channel alpha-subunit 5a (SCN5A). Genetic variants may be the culprit for arrhythmogenesis in Bi-MVP patients. (Level of Difficulty: Intermediate.), Central Illustration

Gain-of-function (GOF) variants give rise to increased or novel protein functions whereas loss-of-function (LOF) variants lead to diminished protein function. GOF and LOF variants can result in markedly varying phenotypes, even when occurring in the same gene. However, experimental approaches for identifying GOF and LOF are generally slow and costly, whilst currently available computational methods have not been optimized to discriminate between GOF and LOF variants. We have developed LoGoFunc, an ensemble machine learning method for predicting pathogenic GOF, pathogenic LOF, and neutral genetic variants. LoGoFunc was trained on a broad range of gene-, protein-, and variant-level features describing diverse biological characteristics, as well as network features summarizing the protein-protein interactome and structural features calculated from AlphaFold2 protein models. We analyzed GOF, LOF, and neutral variants in terms of local protein structure and function, splicing disruption, and phenotypic associations, thereby revealing previously unreported relationships between various biological phenomena and variant functional outcomes. For example, GOF and LOF variants exhibit contrasting enrichments in protein structural and functional regions, whilst LOF variants are more likely to disrupt canonical splicing as indicated by splicing-related features employed by the model. Further, by performing phenome-wide association studies (PheWAS), we identified strong associations between relevant phenotypes and high-confidence predicted GOF and LOF variants. LoGoFunc outperforms other tools trained solely to predict pathogenicity or general variant impact for the identification of pathogenic GOF and LOF variants.

Attention deficit/hyperactivity disorder (ADHD) is a common and heritable phenotype frequently accompanied by insomnia, anxiety, and depression. Here, using a reverse phenotyping approach, we report heterozygous coding variations in the core circadian clock gene cryptochrome 1 in 15 unrelated multigenerational families with combined ADHD and insomnia. The variants led to functional alterations in the circadian molecular rhythms, providing a mechanistic link to the behavioral symptoms. One variant, CRY1Δ11 c.1657+3A>C, is present in approximately 1% of Europeans, therefore standing out as a diagnostic and therapeutic marker. We showed by exome sequencing in an independent cohort of patients with combined ADHD and insomnia that 8 of 62 patients and 0 of 369 controls carried CRY1Δ11. Also, we identified a variant, CRY1Δ6 c.825+1G>A, that shows reduced affinity for BMAL1/CLOCK and causes an arrhythmic phenotype. Genotype-phenotype correlation analysis revealed that this variant segregated with ADHD and delayed sleep phase disorder (DSPD) in the affected family. Finally, we found in a phenome-wide association study involving 9438 unrelated adult Europeans that CRY1Δ11 was associated with major depressive disorder, insomnia, and anxiety. These results defined a distinctive group of circadian psychiatric phenotypes that we propose to designate as "circiatric" disorders., NIH Clinical and Translational Science Award (CTSA) Program; NIH; French National Research Agency (ANR) under the “Investments for the future” Program; Integrative Biology of Emerging Infectious Diseases Laboratoire d’Excellence; IEIHSEER Grant; SEAe-Host Factors Grant; PNEUMOID Project Grant; INCA/Cancéropole Ile-de-France; Turkish Academy of Sciences (TÜBA); National Center for Advancing Translational Sciences (NCAST); Rockefeller University, INSERM; HHMI, University of Paris; St. Giles Foundation; Charles Bronfman Institute for Personalized Medicine at the Icahn School of Medicine at Mount Sinai

Over the last decade next generation sequencing (NGS) has been extensively used to identify new pathogenic mutations and genes causing rare genetic diseases. The efficient analyses of NGS data is not trivial and requires a technically and biologically rigorous pipeline that addresses data quality control, accurate variant filtration to minimize false positives and false negatives, and prioritization of the remaining genes based on disease genomics and physiological knowledge. This review provides a pipeline including all these steps, describes popular software for each step of the analysis, and proposes a general framework for the identification of causal mutations and genes in individual patients of rare genetic diseases.

The construction of population-based variomes has contributed substantially to our understanding of the genetic basis of human inherited disease. Here, we investigated the genetic structure of Turkey from 3,362 unrelated subjects whose whole exomes (n = 2,589) or whole genomes (n = 773) were sequenced to generate a Turkish (TR) Variome that should serve to facilitate disease gene discovery in Turkey. Consistent with the history of present-day Turkey as a crossroads between Europe and Asia, we found extensive admixture between Balkan, Caucasus, Middle Eastern, and European populations with a closer genetic relationship of the TR population to Europeans than hitherto appreciated. We determined that 50\% of TR individuals had high inbreeding coefficients (>= 0.0156) with runs of homozygosity longer than 4 Mb being found exclusively in the TR population when compared to 1000 Genomes Project populations. We also found that 28\% of exome and 49\% of genome variants in the very rare range (allele frequency < 0.005) are unique to the modern TR population. We annotated these variants based on their functional consequences to establish a TR Variome containing alleles of potential medical relevance, a repository of homozygous loss-of-function variants and a TR reference panel for genotype imputation using high-quality haplotypes, to facilitate genome-wide association studies. In addition to providing information on the genetic structure of the modern TR population, these data provide an invaluable resource for future studies to identify variants that are associated with specific phenotypes as well as establishing the phenotypic consequences of mutations in specific genes.

Background: Acute myocarditis (AM) is a well-known cause of sudden death and heart failure, often caused by prevalent viruses. We previously showed that some pediatric AM correlates with putatively damaging variants in genes related to cardiomyocyte structure and function. We sought to evaluate whether deleterious cardiomyopathic variants were enriched among fatal pediatric AM cases in New York City compared with ancestry-matched controls. Methods: Twenty-four children (aged 3 weeks to 20 years) with death due to AM were identified through autopsy records; histologies were reviewed to confirm that all cases met Dallas criteria for AM and targeted panel sequencing of 57 cardiomyopathic genes was performed. Controls without cardiovascular disease were identified from a pediatric database and matched by genetic ancestry to cases using principal components from exome sequencing. Rates of putative deleterious variations (DV) were compared between cases and controls. Where available, AM tissues underwent viral analysis by polymerase chain reaction. Results: DV were identified in 4 of 24 AM cases (16.7%), compared with 2 of 96 age and ancestry-matched controls (2.1%, P =0.014). Viral causes were proven for 6 of 8 AM cases (75%), including the one DV+ case where tissue was available for testing. DV+ cases were more likely to be female, have no evidence of chronic inflammation, and associate with sudden cardiac death than DV− cases. Conclusions: Deleterious variants in genes related to cardiomyocyte integrity are more common in children with fatal AM than controls, likely conferring susceptibility. Additionally, genetically mediated AM may progress more rapidly and be more severe.

BACKGROUND AND AIMS: Polygenic risk scores (PRS) may soon be used to predict inflammatory bowel disease (IBD) risk in prevention efforts. We leveraged exome-sequence and single nucleotide polymorphism (SNP) array data from 29,358 individuals in the multiethnic, randomly ascertained health system-based BioMe biobank to define effects of common and rare IBD variants on disease prediction and pathophysiology. METHODS: PRS were calculated from European, African American, and Ashkenazi Jewish (AJ) reference case-control studies, and a meta-GWAS run using all three association datasets. PRS were then combined using regression to assess which combination of scores best predicted IBD status in European, AJ, Hispanic, and African American cohorts in BioMe. Additionally, rare variants were assessed in genes associated with very early-onset IBD (VEO-IBD), by estimating genetic penetrance in each BioMe population. RESULTS: Combining risk scores based on association data from distinct ancestral populations improved IBD prediction for every population in BioMe and significantly improved prediction among European ancestry UK Biobank individuals. Lower predictive power for non-Europeans was observed, reflecting in part substantially lower African IBD case-control reference sizes. We replicated associations for two VEO-IBD genes, ADAM17 and LRBA, with high dominant model penetrance in BioMe. Autosomal recessive LRBA risk alleles are associated with severe, early-onset autoimmunity; we show that heterozygous carriage of an African-predominant LRBA protein-altering allele is associated with significantly decreased LRBA and CTLA-4 expression with T-cell activation. CONCLUSIONS: Greater genetic diversity in African populations improves prediction across populations, and generalizes some VEO-IBD genes. Increasing African American IBD case-collections should be prioritized to reduce health disparities and enhance pathophysiological insight.

Obesity is a major risk factor for cardiometabolic diseases. Nevertheless, a substantial proportion of individuals with obesity do not suffer cardiometabolic comorbidities. The mechanisms that uncouple adiposity from its cardiometabolic complications are not fully understood. Here, we identify 62 loci of which the same allele is significantly associated with both higher adiposity and lower cardiometabolic risk. Functional analyses show that the 62 loci are enriched for genes expressed in adipose tissue, and for regulatory variants that influence nearby genes that affect adipocyte differentiation. Genes prioritized in each locus support a key role of fat distribution (FAM13A, IRS1 and PPARG) and adipocyte function (ALDH2, CCDC92, DNAH10, ESR1, FAM13A, MTOR, PIK3R1 and VEGFB). Several additional mechanisms are involved as well, such as insulin–glucose signalling (ADCY5, ARAP1, CREBBP, FAM13A, MTOR, PEPD, RAC1 and SH2B3), energy expenditure and fatty acid oxidation (IGF2BP2), browning of white adipose tissue (CSK, VEGFA, VEGFB and SLC22A3) and inflammation (SH2B3, DAGLB and ADCY9). Some of these genes may represent therapeutic targets to reduce cardiometabolic risk linked to excess adiposity.

Variability in SARS-CoV-2 susceptibility and COVID-19 disease severity between individuals is partly due to genetic factors. Here, we identify 4 genomic loci with suggestive associations for SARS-CoV-2 susceptibility and 19 for COVID-19 disease severity. Four of these 23 loci likely have an ethnicity-specific component. Genome-wide association study (GWAS) signals in 11 loci colocalize with expression quantitative trait loci (eQTLs) associated with the expression of 20 genes in 62 tissues/cell types (range: 1:43 tissues/gene), including lung, brain, heart, muscle, and skin as well as the digestive system and immune system. We perform genetic fine mapping to compute 99% credible SNP sets, which identify 10 GWAS loci that have eight or fewer SNPs in the credible set, including three loci with one single likely causal SNP. Our study suggests that the diverse symptoms and disease severity of COVID-19 observed between individuals is associated with variants across the genome, affecting gene expression levels in a wide variety of tissue types., Graphical abstract, D’Antonio et al. characterize associations between GWAS signals for COVID-19 disease and eQTLs in 69 human tissues to identify causal variants and their underlying molecular mechanisms. They show that diverse symptoms and disease severity of COVID-19 are associated with variants affecting gene expression in a wide variety of tissues.

Gain-of-function (GOF) and loss-of-function (LOF) mutations in the same gene may result in markedly different clinical phenotypes and hence require different therapeutic treatments. Identifying the functional consequences of mutations is an important step toward understanding disease mechanisms. While there are numerous computational tools (e.g., CADD, SIFT, PolyPhen-2) to predict the pathogenicity of a variant, there are currently no methods to predict whether a given genetic mutation results in a gene product with increased (gain-of-function; GOF) or diminished (loss-of-function; LOF) activity. Here, we investigated various gene- and protein-level features of GOF and LOF mutations. We generated the first extensive database of all currently known germline GOF and LOF pathogenic mutations by employing natural language processing (NLP) on the available abstracts in the Human Gene Mutation Database. Machine learning and statistical analyses of gene- and protein-level features associated with GOF and LOF mutations indicated significant differences. For example, GOF mutations were enriched in essential genes, autosomal dominant inheritance, protein binding and interaction domains, whereas LOF mutations were enriched in singleton genes, protein-truncating variants, and protein core regions. These results are consistent with the notion that mutations underlying recessive and dominant disorders have significantly different structural and functional properties. These findings could ultimately improve our understanding of how mutations affect gene/protein function thereby guiding future treatment options.

Human herpes simplex virus 1 (HSV-1) encephalitis can be caused by inborn errors of the TLR3 pathway, resulting in impairment of CNS cell-intrinsic antiviral immunity. Deficiencies of the TLR3 pathway impair cell-intrinsic immunity to vesicular stomatitis virus (VSV) and HSV-1 in fibroblasts, and to HSV-1 in cortical but not trigeminal neurons. The underlying molecular mechanism is thought to involve impaired IFN-α/β induction by the TLR3 recognition of dsRNA viral intermediates or by-products. However, we show here that human TLR3 controls constitutive levels of IFNB mRNA and secreted bioactive IFN-β protein, and thereby also controls constitutive mRNA levels for IFN-stimulated genes (ISGs) in fibroblasts. Tlr3-/- mouse embryonic fibroblasts also have lower basal ISG levels. Moreover, human TLR3 controls basal levels of IFN-β secretion and ISG mRNA in induced pluripotent stem cell-derived cortical neurons. Consistently, TLR3-deficient human fibroblasts and cortical neurons are vulnerable not only to both VSV and HSV-1, but also to several other families of viruses. The mechanism by which TLR3 restricts viral growth in human fibroblasts and cortical neurons in vitro and, by inference, by which the human CNS prevents infection by HSV-1 in vivo, is therefore based on the control of early viral infection by basal IFN-β immunity. ispartof: JOURNAL OF CLINICAL INVESTIGATION vol:131 issue:1 ispartof: location:United States status: published

Inflammatory bowel disease (IBD) is a group of chronic diseases, affecting different parts of the gastrointestinal tract, that mainly comprises Crohn’s Disease (CD) and Ulcerative Colitis (UC). Most IBD genomic research to date has involved genome-wide association studies (GWAS) of common genetic variants, mostly in Europeans, resulting in the identification of over 200 risk loci. The incidence of IBD in Ashkenazi Jews (AJ) is particularly high compared to other population groups and rare protein-coding variants are significantly enriched in AJ. These variants are expected to have a larger phenotypic effect and are hypothesized to complement the missing heritability that cannot be fully addressed by GWAS in IBD. Therefore, we genetically identified 4,974 AJs IBD cases and controls from whole exome sequencing (WES) data from the NIDDK IBD Genetics Consortium (IBDGC). We selected credible rare variants with high predicted impact, aggregated them into genes, and performed gene burden and pathway enrichment analyses to identify 7 novel plausible IBD-causing genes:NCF1, CES1, ICAM1, INPP5D, ABCB1, IL33andTLR4. We further perform bulk and single-cell RNA sequencing, demonstrating the likely relatedness of the novel genes to IBD. Importantly, we demonstrate that the rare and high impact genetic architecture of AJ adult IBD displays a significant overlap with very early onset IBD (VEOIBD) genetics. At the variant level, we performed Phenome-wide association studies (PheWAS) in the UK Biobank to replicate risk sites in IBD and reveal shared risk sites with other diseases. Finally, we showed that a polygenic risk score (PRS) has high power to differentiate AJ IBD cases from controls when using rare and high impact variants.

X-linked dominant incontinentia pigmenti (IP) and X-linked recessive anhidrotic ectodermal dysplasia with immunodeficiency (EDA-ID) are caused by loss-of-function and hypomorphic IKBKG (also known as NEMO) mutations, respectively. We describe a European mother with mild IP and a Japanese mother without IP, whose 3 boys with EDA-ID died from ID. We identify the same private variant in an intron of IKBKG, IVS4+866 C>T, which was inherited from and occurred de novo in the European mother and Japanese mother, respectively. This mutation creates a new splicing donor site, giving rise to a 44-nucleotide pseudoexon (PE) generating a frameshift. Its leakiness accounts for NF-κB activation being impaired but not abolished in the boys’ cells. However, aberrant splicing rates differ between cell types, with WT NEMO mRNA and protein levels ranging from barely detectable in leukocytes to residual amounts in induced pluripotent stem cell–derived (iPSC-derived) macrophages, and higher levels in fibroblasts and iPSC-derived neuronal precursor cells. Finally, SRSF6 binds to the PE, facilitating its inclusion. Moreover, SRSF6 knockdown or CLK inhibition restores WT NEMO expression and function in mutant cells. A recurrent deep intronic splicing mutation in IKBKG underlies a purely quantitative NEMO defect in males that is most severe in leukocytes and can be rescued by the inhibition of SRSF6 or CLK.

Human inborn errors of IFN-γ immunity underlie mycobacterial diseases. We describe patients with Mycobacterium bovis (BCG) disease who are homozygous for loss-of-function mutations of SPPL2A. This gene encodes a transmembrane protease that degrades the N-terminal fragment (NTF) of CD74 (HLA invariant chain) in antigen-presenting cells. The CD74 NTF therefore accumulates in the HLA class II+ myeloid and lymphoid cells of SPPL2a-deficient patients. This toxic fragment selectively depletes IL-12- and IL-23-producing CD1c+ conventional dendritic cells (cDC2s) and their circulating progenitors. Moreover, SPPL2a-deficient memory TH1* cells selectively fail to produce IFN-γ when stimulated with mycobacterial antigens in vitro. Finally, Sppl2a–/– mice lack cDC2s, have CD4+ T cells that produce small amounts of IFN-γ after BCG infection, and are highly susceptible to infection with BCG or Mycobacterium tuberculosis. These findings suggest that inherited SPPL2a deficiency in humans underlies mycobacterial disease by decreasing the numbers of cDC2s and impairing IFN-γ production by mycobacterium-specific memory TH1* cells. Primary immunodeficiency can predispose patients to mycobacterial disease. Casanova and colleagues identify novel human mutations in the enzyme SPPL2A that result in selective accumulation of CD74 in a dendritic cell subset and lead to their death and the failure to mount effective TH1 responses.

Background Germline heterozygous mutations in human signal transducer and activator of transcription 1 (STAT1) can cause loss of function (LOF), as in patients with Mendelian susceptibility to mycobacterial diseases, or gain of function (GOF), as in patients with chronic mucocutaneous candidiasis. LOF and GOF mutations are equally rare and can affect the same domains of STAT1, especially the coiled-coil domain (CCD) and DNA-binding domain (DBD). Moreover, 6% of patients with chronic mucocutaneous candidiasis with a GOF STAT1 mutation have mycobacterial disease, obscuring the functional significance of the identified STAT1 mutations. Current computational approaches, such as combined annotation-dependent depletion, do not distinguish LOF and GOF variants. Objective We estimated variations in the CCD/DBD of STAT1. Methods We mutagenized 342 individual wild-type amino acids in the CCD/DBD (45.6% of full-length STAT1) to alanine and tested the mutants for STAT1 transcriptional activity. Results Of these 342 mutants, 201 were neutral, 30 were LOF, and 111 were GOF mutations in a luciferase assay. This assay system correctly estimated all previously reported LOF mutations (100%) and slightly fewer GOF mutations (78.1%) in the CCD/DBD of STAT1. We found that GOF alanine mutants occurred at the interface of the antiparallel STAT1 dimer, suggesting that they destabilize this dimer. This assay also precisely predicted the effect of 2 hypomorphic and dominant negative mutations, E157K and G250E, in the CCD of STAT1 that we found in 2 unrelated patients with Mendelian susceptibility to mycobacterial diseases. Conclusion The systematic alanine-scanning assay is a useful tool to estimate the GOF or LOF status and the effect of heterozygous missense mutations in STAT1 identified in patients with severe infectious diseases, including mycobacterial and fungal diseases.

Background Congenital heart disease (CHD) affects ~ 1% of live births and is the most common birth defect. Although the genetic contribution to the CHD has been long suspected, it has only been well established recently. De novo variants are estimated to contribute to approximately 8% of sporadic CHD. Methods CHD is genetically heterogeneous, making pathway enrichment analysis an effective approach to explore and statistically validate CHD-associated genes. In this study, we performed novel gene and pathway enrichment analyses of high-impact de novo variants in the recently published whole-exome sequencing (WES) data generated from a cohort of CHD 2645 parent-offspring trios to identify new CHD-causing candidate genes and mutations. We performed rigorous variant- and gene-level filtrations to identify potentially damaging variants, followed by enrichment analyses and gene prioritization. Results Our analyses revealed 23 novel genes that are likely to cause CHD, including HSP90AA1, ROCK2, IQGAP1, and CHD4, and sharing biological functions, pathways, molecular interactions, and properties with known CHD-causing genes. Conclusions Ultimately, these findings suggest novel genes that are likely to be contributing to CHD pathogenesis.

SARS-CoV-2 infection displays immense inter-individual clinical variability, ranging from silent infection to lethal disease. The role of human genetics in determining clinical response to the virus remains unclear. Studies of outliers—individuals remaining uninfected despite viral exposure and healthy young patients with life-threatening disease—present a unique opportunity to reveal human genetic determinants of infection and disease.

Summary Next-generation sequencing (NGS) generates large amounts of genomic data and reveals about 20 000 genetic coding variants per individual studied. Several mutation damage prediction scores are available to prioritize variants, but there is currently no application to help investigators to determine the relevance of the candidate genes and variants quickly and visually from population genetics data and deleteriousness scores. Here, we present PopViz, a user-friendly, rapid, interactive, mobile-compatible webserver providing a gene-centric visualization of the variants of any human gene, with (i) population-specific minor allele frequencies from the gnomAD population genetic database; (ii) mutation damage prediction scores from CADD, EIGEN and LINSIGHT and (iii) amino-acid positions and protein domains. This application will be particularly useful in investigations of NGS data for new disease-causing genes and variants, by reinforcing or rejecting the plausibility of the candidate genes, and by selecting and prioritizing, the candidate variants for experimental testing. Availability and implementation PopViz webserver is freely accessible from http://shiva.rockefeller.edu/PopViz/. Supplementary information Supplementary data are available at Bioinformatics online.

International audience; Herpes simplex virus-1 (HSV-1) encephalitis (HSE) is typically sporadic. Inborn errors of TLR3- and DBR1-mediated central nervous system cell-intrinsic immunity can account for forebrain and brainstem HSE, respectively. We report five unrelated patients with forebrain HSE, each heterozygous for one of four rare variants of SNORA31, encoding a small nucleolar RNA of the H/ACA class that are predicted to direct the isomerization of uridine residues to pseudouridine in small nuclear RNA and ribosomal RNA. We show that CRISPR/Cas9-introduced bi- and monoallelic SNORA31 deletions render human pluripotent stem cell (hPSC)-derived cortical neurons susceptible to HSV-1. Accordingly, SNORA31-mutated patient hPSC-derived cortical neurons are susceptible to HSV-1, like those from TLR3- or STAT1-deficient patients. Exogenous interferon (IFN)-β renders SNORA31- and TLR3- but not STAT1-mutated neurons resistant to HSV-1. Finally, transcriptome analysis of SNORA31-mutated neurons revealed normal responses to TLR3 and IFN-α/β stimulation but abnormal responses to HSV-1. Human SNORA31 thus controls central nervous system neuron-intrinsic immunity to HSV-1 by a distinctive mechanism.

Genetic etiologies of chronic mucocutaneous candidiasis (CMC) disrupt human IL-17A/F-dependent immunity at mucosal surfaces, whereas those of connective tissue disorders (CTDs) often impair the TGF-β-dependent homeostasis of connective tissues. The signaling pathways involved are incompletely understood. We report a three-generation family with an autosomal dominant (AD) combination of CMC and a previously undescribed form of CTD that clinically overlaps with Ehlers-Danlos syndrome (EDS). The patients are heterozygous for a private splice-site variant of MAPK8, the gene encoding c-Jun N-terminal kinase 1 (JNK1), a component of the MAPK signaling pathway. This variant is loss-of-expression and loss-of-function in the patients' fibroblasts, which display AD JNK1 deficiency by haploinsufficiency. These cells have impaired, but not abolished, responses to IL-17A and IL-17F. Moreover, the development of the patients' TH17 cells was impaired ex vivo and in vitro, probably due to the involvement of JNK1 in the TGF-β-responsive pathway and further accounting for the patients' CMC. Consistently, the patients' fibroblasts displayed impaired JNK1- and c-Jun/ATF-2-dependent induction of key extracellular matrix (ECM) components and regulators, but not of EDS-causing gene products, in response to TGF-β. Furthermore, they displayed a transcriptional pattern in response to TGF-β different from that of fibroblasts from patients with Loeys-Dietz syndrome caused by mutations of TGFBR2 or SMAD3, further accounting for the patients' complex and unusual CTD phenotype. This experiment of nature indicates that the integrity of the human JNK1-dependent MAPK signaling pathway is essential for IL-17A- and IL-17F-dependent mucocutaneous immunity to Candida and for the TGF-β-dependent homeostasis of connective tissues.

Network biology has the capability to integrate, represent, interpret, and model complex biological systems by collectively accommodating biological omics data, biological interactions and associations, graph theory, statistical measures, and visualizations. Biological networks have recently been shown to be very useful for studies that decipher biological mechanisms and disease etiologies and for studies that predict therapeutic responses, at both the molecular and system levels. In this review, we briefly summarize the general framework of biological network studies, including data resources, network construction methods, statistical measures, network topological properties, and visualization tools. We also introduce several recent biological network applications and methods for the studies of rare diseases.

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/.