The human gene damage index as a gene-level approach to prioritizing exome variants

FDPSM: Feature-Driven Prediction Modeling of Pathogenic Synonymous Mutations

Synonymous mutations, once considered to be biologically neutral, are now recognized to affect protein expression and function by altering the RNA splicing, stability, or translation efficiency. These effects can contribute to disease, making the prediction of the pathogenicity a crucial task. Computational methods have been developed to analyze the sequence features and biological functions of synonymous mutations, but existing methods face limitations, including scarcity of labeled data, reliance on other prediction tools, and insufficient representation of feature interrelationships. Here, we present FDPSM, a novel prediction method specifically designed to predict pathogenic synonymous mutations. FDPSM was trained on a robust data set of 4251 positive and negative training samples to enhance predictive accuracy. The method leveraged a comprehensive set of features, including genomic context, conservation, splicing effects, functional effects, and epigenomics, without relying on prediction scores from other mutation pathogenicity tools. Recognizing that original features alone may not fully capture the distinctions between pathogenic and benign synonymous mutations, we enhanced the feature set by extracting effective information from the interactions and distribution of these features. The experimental results showed that FDPSM significantly outperformed existing methods in predicting the pathogenicity of synonymous mutations, offering a more accurate and reliable tool for this important task. FDPSM is available at https://github.com/xialab-ahu/FDPSM.

Deciphering the digenic architecture of congenital heart disease using trio exome sequencing data

Congenital heart disease (CHD) is the most common congenital anomaly and a leading cause of infant morbidity and mortality. Despite extensive exploration of the monogenic causes of CHD over the last decades, ∼55% of cases still lack a molecular diagnosis. Investigating digenic interactions, the simplest form of oligogenic interactions, using high-throughput sequencing data can elucidate additional genetic factors contributing to the disease. Here, we conducted a comprehensive analysis of digenic interactions in CHD by utilizing a large CHD trio exome sequencing cohort, comprising 3,910 CHD and 3,644 control trios. We extracted pairs of presumably deleterious rare variants observed in CHD-affected and unaffected children but not in a single parent. Burden testing of gene pairs derived from these variant pairs revealed 29 nominally significant gene pairs. These gene pairs showed a significant enrichment for known CHD genes (p < 1.0 × 10-4) and exhibited a shorter average biological distance to known CHD genes than expected by chance (p = 3.0 × 10-4). Utilizing three complementary biological relatedness approaches including network analyses, biological distance calculations, and candidate gene prioritization methods, we prioritized 10 final gene pairs that are likely to underlie CHD. Analysis of bulk RNA-sequencing data showed that these genes are highly expressed in the developing embryonic heart (p < 1 × 10-4). In conclusion, our findings suggest the potential role of digenic interactions in CHD pathogenesis and provide insights into unresolved molecular diagnoses. We suggest that the application of the digenic approach to additional disease cohorts will significantly enhance genetic discovery rates.

A graph neural network approach for accurate prediction of pathogenicity in multi-type variants

Accurate prediction of pathogenic variants in human disease-associated genes would have a profound effect on clinical decision-making; however, it remains a significant challenge due to the overwhelming number of these variants. We propose graph neural network for multimodal annotation-based pathogenicity prediction (GNN-MAP), a novel deep learning framework that effectively integrates multimodal annotations and similarity relationships among variants to predict the pathogenicity of multi-type variants. Trained on the ClinVar dataset, GNN-MAP exhibits superior predictive performance in internal validation and orthogonal test datasets, accurately predicting variant pathogenicity. Notably, GNN-MAP enables accurate prediction of the pathogenicity of rare variants and highly imbalanced datasets. Furthermore, it achieves high performance in the pathogenicity prediction of inherited retinal disease-specific variants, highlighting its effectiveness in disease-specific variant prediction. These findings suggest that the robust capability of GNN-MAP to predict pathogenicity across multiple variant types and datasets holds significant potential for applications in research and clinical settings.

Genetic Landscape and Mitochondrial Metabolic Dysregulation in Patients Suffering From Severe Long COVID

Long COVID represents a significant global health challenge with an unclear etiology. Alongside accumulating evidence of mitochondrial dysfunction in patients with acute SARS-CoV-2 infection, a symptomatic overlap exists between long COVID and mitochondrial disorders. However, the genetic underpinnings of mitochondrial dysfunction in long COVID have not been previously explored. We employed whole genome sequencing to analyze 13 patients with severe long COVID to identify genetic defects related to mitochondrial function. We performed extracellular bioenergetics flux analysis on peripheral blood mononuclear cells and proteomics to evaluate cellular bioenergetics and compared the results to those of healthy controls. Our investigation identified 10 variants classified as pathogenic or likely pathogenic and 83 variants of unknown significance affecting a wide range of mitochondria-associated biological functions. Bioenergetics flux analysis in peripheral blood mononuclear cells revealed an altered ATP production rate in four long COVID patients compared to healthy controls. This study presents initial evidence of a potential underlying genetic predisposition to mitochondrial dysfunction in long COVID while demonstrating altered cellular energy capacity in a subset of these patients. These findings open avenues for further research into the role of mitochondrial dysfunction and pathology in patients suffering from long COVID and may pave the way for targeted therapeutic strategies aimed at mitigating mitochondrial dysfunction.

Integrating rare pathogenic variant prioritization with gene-based association analysis to identify novel genes and relevant multimodal traits for Alzheimer's disease

INTRODUCTION

Increasing evidence has highlighted rare variants in Alzheimer's disease (AD). However, insufficient sample sizes, especially in underrepresented ethnic groups, hinder their investigation. Additionally, their impact on endophenotypes remains largely unexplored.

METHODS

We prioritized rare likely-deleterious variants based on whole-genome sequencing data from a Chinese AD cohort (n = 988). Gene-based optimal sequence kernel association tests were conducted between AD cases and normal controls to identify AD-related genes. Network clustering, endophenotype association, and cellular experiments were conducted to evaluate their functional consequences.

RESULTS

We identified 11 novel AD candidate genes, which captured AD-related pathways and enhanced AD risk prediction performance. Key genes (RABEP1, VIPR1, RPL3L, and CABIN1) were linked to cognitive decline and brain atrophy. Experiments showed RABEP1 p.R845W inducing endocytosis dysregulation and exacerbating toxic amyloid β accumulation, underscoring its therapeutic potential.

DISCUSSION

Our findings highlighted the contributions of rare variants to AD and provided novel insights into AD therapeutics.

HIGHLIGHTS

Identified 11 novel AD candidate genes in a Chinese AD cohort. Correlated candidate genes with AD-related cognitive and brain imaging traits. Indicated RABEP1 p.R845W as a critical AD contributor in the endocytic pathway.

GoFCards: an integrated database and analytic platform for gain of function variants in humans

Gain-of-function (GOF) variants, which introduce new or amplify protein functions, are essential for understanding disease mechanisms. Despite advances in genomics and functional research, identifying and analyzing pathogenic GOF variants remains challenging owing to fragmented data and database limitations, underscoring the difficulty in accessing critical genetic information. To address this challenge, we manually reviewed the literature, pinpointing 3089 single-nucleotide variants and 72 insertions and deletions in 579 genes associated with 1299 diseases from 2069 studies, and integrated these with the 3.5 million predicted GOF variants. Our approach is complemented by a proprietary scoring system that prioritizes GOF variants on the basis of the evidence supporting their GOF effects and provides predictive scores for variants that lack existing documentation. We then developed a database named GoFCards for general geneticists and clinicians to easily obtain GOF variants in humans (http://www.genemed.tech/gofcards). This database also contains data from >150 sources and offers comprehensive variant-level and gene-level annotations, with the aim of providing users with convenient access to detailed and relevant genetic information. Furthermore, GoFCards empowers users with limited bioinformatic skills to analyze and annotate genetic data, and prioritize GOF variants. GoFCards offers an efficient platform for interpreting GOF variants and thereby advancing genetic research.

Rare Genetic Variants of NLRP12 in Admixed Latino-American Children With SARS-CoV-2-Related Multisystem Inflammatory Syndrome

Multisystem inflammatory syndrome in children (MIS-C) is a rare, potentially fatal complication of SARS-CoV-2 infection. Genetic defects in inflammation-related pathways have been linked to MIS-C, but additional research is needed, especially in diverse ethnic groups. The present study aimed to identify genetic variants underlying MIS-C in Brazilian patients. Whole exome sequencing was performed, focusing on genes involved in the host immune response to SARS-CoV-2. Functional assays assessed the impact of selected variants on nuclear factor-κB signaling. Nine rare, potentially deleterious variants were found in 8 of 21 patients, located in the IL17RC, IFNA10, or NLRP12 gene. Unlike the wild type NLRP12 protein, which inhibits nuclear factor-κB activation in HEK 293T cells, the mutant NLRP12 proteins have significantly reduced inhibitory properties. In conclusion, our results indicate that rare autosomal variants in immune-related genes may underlie MIS-C, highlighting the potential role of NLRP12 in its predisposition. These findings provide new insights for the appropriate management of MIS-C.

ACK1 and BRK non-receptor tyrosine kinase deficiencies are associated with familial systemic lupus and involved in efferocytosis

Systemic lupus erythematosus (SLE) is an autoimmune disease, the pathophysiology and genetic basis of which are incompletely understood. Using a forward genetic screen in multiplex families with SLE, we identified an association between SLE and compound heterozygous deleterious variants in the non-receptor tyrosine kinases (NRTKs) ACK1 and BRK. Experimental blockade of ACK1 or BRK increased circulating autoantibodies in vivo in mice and exacerbated glomerular IgG deposits in an SLE mouse model. Mechanistically, NRTKs regulate activation, migration, and proliferation of immune cells. We found that the patients' ACK1 and BRK variants impair efferocytosis, the MERTK-mediated anti-inflammatory response to apoptotic cells, in human induced pluripotent stem cell (hiPSC)-derived macrophages, which may contribute to SLE pathogenesis. Overall, our data suggest that ACK1 and BRK deficiencies are associated with human SLE and impair efferocytosis in macrophages.

Genetic polymorphism in untranslated regions of PRKCZ influences mRNA structure, stability and binding sites

BACKGROUND

Variations in untranslated regions (UTR) alter regulatory pathways impacting phenotype, disease onset, and course of disease. Protein kinase C Zeta (PRKCZ), a serine-threonine kinase, is implicated in cardiovascular, neurological and oncological disorders. Due to limited research on PRKCZ, this study aimed to investigate the impact of UTR genetic variants' on binding sites for transcription factors and miRNA. RNA secondary structure, eQTLs, and variation tolerance analysis were also part of the study.

METHODS

The data related to PRKCZ gene variants was downloaded from the Ensembl genome browser, COSMIC and gnomAD. The RegulomeDB database was used to assess the functional impact of 5' UTR and 3'UTR variants. The analysis of the transcription binding sites (TFBS) was done through the Alibaba tool, and the Kyoto Encyclopaedia of Genes and Genomes (KEGG) was employed to identify pathways associated with PRKCZ. To predict the effect of variants on microRNA binding sites, PolymiRTS was utilized for 3' UTR variants, and the SNPinfo tool was used for 5' UTR variants.

RESULTS

The results obtained indicated that a total of 24 variants present in the 3' UTR and 25 variants present in the 5' UTR were most detrimental. TFBS analysis revealed that 5' UTR variants added YY1, repressor, and Oct1, whereas 3' UTR variants added AP-2alpha, AhR, Da, GR, and USF binding sites. The study predicted TFs that influenced PRKCZ expression. RNA secondary structure analysis showed that eight 5' UTR and six 3' UTR altered the RNA structure by either removal or addition of the stem-loop. The microRNA binding site analysis highlighted that seven 3' UTR and one 5' UTR variant altered the conserved site and also created new binding sites. eQTLs analysis showed that one variant was associated with PRKCZ expression in the lung and thyroid. The variation tolerance analysis revealed that PRKCZ was an intolerant gene.

CONCLUSION

This study laid the groundwork for future studies aimed at targeting PRKCZ as a therapeutic target.

Inherited human RelB deficiency impairs innate and adaptive immunity to infection

We report two unrelated adults with homozygous (P1) or compound heterozygous (P2) private loss-of-function variants of V-Rel Reticuloendotheliosis Viral Oncogene Homolog B (RELB). The resulting deficiency of functional RelB impairs the induction of NFKB2 mRNA and NF-κB2 (p100/p52) protein by lymphotoxin in the fibroblasts of the patients. These defects are rescued by transduction with wild-type RELB complementary DNA (cDNA). By contrast, the response of RelB-deficient fibroblasts to Tumor Necrosis Factor (TNF) or IL-1β via the canonical NF-κB pathway remains intact. P1 and P2 have low proportions of naïve CD4+ and CD8+ T cells and of memory B cells. Moreover, their naïve B cells cannot differentiate into immunoglobulin G (IgG)- or immunoglobulin A (IgA)-secreting cells in response to CD40L/IL-21, and the development of IL-17A/F-producing T cells is strongly impaired in vitro. Finally, the patients produce neutralizing autoantibodies against type I interferons (IFNs), even after hematopoietic stem cell transplantation, attesting to a persistent dysfunction of thymic epithelial cells in T cell selection and central tolerance to some autoantigens. Thus, inherited human RelB deficiency disrupts the alternative NF-κB pathway, underlying a T- and B cell immunodeficiency, which, together with neutralizing autoantibodies against type I IFNs, confers a predisposition to viral, bacterial, and fungal infections.

SARS-CoV-2 brainstem encephalitis in human inherited DBR1 deficiency

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.

Human TMEFF1 is a restriction factor for herpes simplex virus in the brain

Most cases of herpes simplex virus 1 (HSV-1) encephalitis (HSE) remain unexplained1,2. Here, we report on two unrelated people who had HSE as children and are homozygous for rare deleterious variants of TMEFF1, which encodes a cell membrane protein that is preferentially expressed by brain cortical neurons. TMEFF1 interacts with the cell-surface HSV-1 receptor NECTIN-1, impairing HSV-1 glycoprotein D- and NECTIN-1-mediated fusion of the virus and the cell membrane, blocking viral entry. Genetic TMEFF1 deficiency allows HSV-1 to rapidly enter cortical neurons that are either patient specific or derived from CRISPR-Cas9-engineered human pluripotent stem cells, thereby enhancing HSV-1 translocation to the nucleus and subsequent replication. This cellular phenotype can be rescued by pretreatment with type I interferon (IFN) or the expression of exogenous wild-type TMEFF1. Moreover, ectopic expression of full-length TMEFF1 or its amino-terminal extracellular domain, but not its carboxy-terminal intracellular domain, impairs HSV-1 entry into NECTIN-1-expressing cells other than neurons, increasing their resistance to HSV-1 infection. Human TMEFF1 is therefore a host restriction factor for HSV-1 entry into cortical neurons. Its constitutively high abundance in cortical neurons protects these cells from HSV-1 infection, whereas inherited TMEFF1 deficiency renders them susceptible to this virus and can therefore underlie HSE.

Genetic constraint at single amino acid resolution in protein domains improves missense variant prioritisation and gene discovery

BACKGROUND

One of the major hurdles in clinical genetics is interpreting the clinical consequences associated with germline missense variants in humans. Recent significant advances have leveraged natural variation observed in large-scale human populations to uncover genes or genomic regions that show a depletion of natural variation, indicative of selection pressure. We refer to this as "genetic constraint". Although existing genetic constraint metrics have been demonstrated to be successful in prioritising genes or genomic regions associated with diseases, their spatial resolution is limited in distinguishing pathogenic variants from benign variants within genes.

METHODS

We aim to identify missense variants that are significantly depleted in the general human population. Given the size of currently available human populations with exome or genome sequencing data, it is not possible to directly detect depletion of individual missense variants, since the average expected number of observations of a variant at most positions is less than one. We instead focus on protein domains, grouping homologous variants with similar functional impacts to examine the depletion of natural variations within these comparable sets. To accomplish this, we develop the Homologous Missense Constraint (HMC) score. We utilise the Genome Aggregation Database (gnomAD) 125 K exome sequencing data and evaluate genetic constraint at quasi amino-acid resolution by combining signals across protein homologues.

RESULTS

We identify one million possible missense variants under strong negative selection within protein domains. Though our approach annotates only protein domains, it nonetheless allows us to assess 22% of the exome confidently. It precisely distinguishes pathogenic variants from benign variants for both early-onset and adult-onset disorders. It outperforms existing constraint metrics and pathogenicity meta-predictors in prioritising de novo mutations from probands with developmental disorders (DD). It is also methodologically independent of these, adding power to predict variant pathogenicity when used in combination. We demonstrate utility for gene discovery by identifying seven genes newly significantly associated with DD that could act through an altered-function mechanism.

CONCLUSIONS

Grouping variants of comparable functional impacts is effective in evaluating their genetic constraint. HMC is a novel and accurate predictor of missense consequence for improved variant interpretation.

ACK1 and BRK non-receptor tyrosine kinase deficiencies are associated with familial systemic lupus and involved in efferocytosis

Systemic Lupus Erythematosus (SLE) is an autoimmune disease, the pathophysiology and genetic basis of which are incompletely understood. Using a forward genetic screen in multiplex families with systemic lupus erythematosus (SLE) we identified an association between SLE and compound heterozygous deleterious variants in the non-receptor tyrosine kinases (NRTKs) ACK1 and BRK. Experimental blockade of ACK1 or BRK increased circulating autoantibodies in vivo in mice and exacerbated glomerular IgG deposits in an SLE mouse model. Mechanistically, non-receptor tyrosine kinases (NRTKs) regulate activation, migration, and proliferation of immune cells. We found that the patients' ACK1 and BRK variants impair efferocytosis, the MERTK-mediated anti-inflammatory response to apoptotic cells, in human induced Pluripotent Stem Cell (hiPSC)-derived macrophages, which may contribute to SLE pathogenesis. Overall, our data suggest that ACK1 and BRK deficiencies are associated with human SLE and impair efferocytosis in macrophages.

The landscape of rare genetic variation associated with inflammatory bowel disease and Parkinson's disease comorbidity

BACKGROUND

Inflammatory bowel disease (IBD) and Parkinson's disease (PD) are chronic disorders that have been suggested to share common pathophysiological processes. LRRK2 has been implicated as playing a role in both diseases. Exploring the genetic basis of the IBD-PD comorbidity through studying high-impact rare genetic variants can facilitate the identification of the novel shared genetic factors underlying this comorbidity.

METHODS

We analyzed whole exomes from the BioMe BioBank and UK Biobank, and whole genomes from a cohort of 67 European patients diagnosed with both IBD and PD to examine the effects of LRRK2 missense variants on IBD, PD and their co-occurrence (IBD-PD). We performed optimized sequence kernel association test (SKAT-O) and network-based heterogeneity clustering (NHC) analyses using high-impact rare variants in the IBD-PD cohort to identify novel candidate genes, which we further prioritized by biological relatedness approaches. We conducted phenome-wide association studies (PheWAS) employing BioMe BioBank and UK Biobank whole exomes to estimate the genetic relevance of the 14 prioritized genes to IBD-PD.

RESULTS

The analysis of LRRK2 missense variants revealed significant associations of the G2019S and N2081D variants with IBD-PD in addition to several other variants as potential contributors to increased or decreased IBD-PD risk. SKAT-O identified two significant genes, LRRK2 and IL10RA, and NHC identified 6 significant gene clusters that are biologically relevant to IBD-PD. We observed prominent overlaps between the enriched pathways in the known IBD, PD, and candidate IBD-PD gene sets. Additionally, we detected significantly enriched pathways unique to the IBD-PD, including MAPK signaling, LPS/IL-1 mediated inhibition of RXR function, and NAD signaling. Fourteen final candidate IBD-PD genes were prioritized by biological relatedness methods. The biological importance scores estimated by protein-protein interaction networks and pathway and ontology enrichment analyses indicated the involvement of genes related to immunity, inflammation, and autophagy in IBD-PD. Additionally, PheWAS provided support for the associations of candidate genes with IBD and PD.

CONCLUSIONS

Our study confirms and uncovers new LRRK2 associations in IBD-PD. The identification of novel inflammation and autophagy-related genes supports and expands previous findings related to IBD-PD pathogenesis, and underscores the significance of therapeutic interventions for reducing systemic inflammation.

Helper T cell immunity in humans with inherited CD4 deficiency

CD4+ T cells are vital for host defense and immune regulation. However, the fundamental role of CD4 itself remains enigmatic. We report seven patients aged 5-61 years from five families of four ancestries with autosomal recessive CD4 deficiency and a range of infections, including recalcitrant warts and Whipple's disease. All patients are homozygous for rare deleterious CD4 variants impacting expression of the canonical CD4 isoform. A shorter expressed isoform that interacts with LCK, but not HLA class II, is affected by only one variant. All patients lack CD4+ T cells and have increased numbers of TCRαβ+CD4-CD8- T cells, which phenotypically and transcriptionally resemble conventional Th cells. Finally, patient CD4-CD8- αβ T cells exhibit intact responses to HLA class II-restricted antigens and promote B cell differentiation in vitro. Thus, compensatory development of Th cells enables patients with inherited CD4 deficiency to acquire effective cellular and humoral immunity against an unexpectedly large range of pathogens. Nevertheless, CD4 is indispensable for protective immunity against at least human papillomaviruses and Trophyrema whipplei.

Single nucleotide polymorphisms conferring susceptibility to leukemia and oral mucositis: a multi-center pilot study of patients prior to conditioning therapy for hematopoietic cell transplant

PURPOSE

Leukemias have been associated with oral manifestations, reflecting susceptibility to cancer therapy-induced oral mucositis. We sought to identify SNPs associated with both leukemia and oral mucositis (OM).

METHODS

Whole exome sequencing was performed on leukemia and non-cancer blood disorder (ncBD) patients' saliva samples (N = 50) prior to conditioning therapy. WHO OM grading scores were determined: moderate to severe (OM2-4) vs. none to mild (OM0-1). Reads were processed using Trim Galorev0.6.7, Bowtie2v2.4.1, Samtoolsv1.10, Genome Analysis Toolkit (GATK)v4.2.6.1, and DeepVariantv1.4.0. We utilized the following pipelines: P1 analysis with PLINK2v3.7, SNP2GENEv1.4.1 and MAGMAv1.07b, and P2 [leukemia (N = 42) vs. ncBDs (N = 8)] and P3 [leukemia + OM2-4 (N = 18) vs. leukemia + OM0-1 (N = 24)] with Z-tests of genotypes and protein-protein interaction determination. GeneCardsSuitev5.14 was used to identify phenotypes (P1 and P2, leukemia; P3, oral mucositis) and average disease-causing likelihood and DGIdb for drug interactions. P1 and P2 genes were analyzed with CytoScape plugin BiNGOv3.0.3 to retrieve overrepresented Gene Ontology (GO) terms and Ensembl's VEP for SNP outcomes.

RESULTS

In P1, 457 candidate SNPs (28 genes) were identified and 21,604 SNPs (1016 genes) by MAGMAv1.07b. Eighteen genes were associated with "leukemia" per VarElectv5.14 analysis and predicted to be deleterious. In P2 and P3, 353 and 174 SNPs were significant, respectively. STRINGv12.0 returned 77 and 32 genes (C.L. = 0.7) for P2 and P3, respectively. VarElectv5.14 determined 60 genes from P2 associated with "leukemia" and 11 with "oral mucositis" from P3. Overrepresented GO terms included "cellular process," "signaling," "hemopoiesis," and "regulation of immune response."

CONCLUSIONS

We identified candidate SNPs possibly conferring susceptibility to develop leukemia and oral mucositis.

The immunopathological landscape of human pre-TCRα deficiency: From rare to common variants

We describe humans with rare biallelic loss-of-function PTCRA variants impairing pre-α T cell receptor (pre-TCRα) expression. Low circulating naive αβ T cell counts at birth persisted over time, with normal memory αβ and high γδ T cell counts. Their TCRα repertoire was biased, which suggests that noncanonical thymic differentiation pathways can rescue αβ T cell development. Only a minority of these individuals were sick, with infection, lymphoproliferation, and/or autoimmunity. We also report that 1 in 4000 individuals from the Middle East and South Asia are homozygous for a common hypomorphic PTCRA variant. They had normal circulating naive αβ T cell counts but high γδ T cell counts. Although residual pre-TCRα expression drove the differentiation of more αβ T cells, autoimmune conditions were more frequent in these patients compared with the general population.

Meta-analysis of 46,000 germline de novo mutations linked to human inherited disease

BACKGROUND

De novo mutations (DNMs) are variants that occur anew in the offspring of noncarrier parents. They are not inherited from either parent but rather result from endogenous mutational processes involving errors of DNA repair/replication. These spontaneous errors play a significant role in the causation of genetic disorders, and their importance in the context of molecular diagnostic medicine has become steadily more apparent as more DNMs have been reported in the literature. In this study, we examined 46,489 disease-associated DNMs annotated by the Human Gene Mutation Database (HGMD) to ascertain their distribution across gene and disease categories.

RESULTS

Most disease-associated DNMs reported to date are found to be associated with developmental and psychiatric disorders, a reflection of the focus of sequencing efforts over the last decade. Of the 13,277 human genes in which DNMs have so far been found, the top-10 genes with the highest proportions of DNM relative to gene size were H3-3 A, DDX3X, CSNK2B, PURA, ZC4H2, STXBP1, SCN1A, SATB2, H3-3B and TUBA1A. The distribution of CADD and REVEL scores for both disease-associated DNMs and those mutations not reported to be de novo revealed a trend towards higher deleteriousness for DNMs, consistent with the likely lower selection pressure impacting them. This contrasts with the non-DNMs, which are presumed to have been subject to continuous negative selection over multiple generations.

CONCLUSION

This meta-analysis provides important information on the occurrence and distribution of disease-associated DNMs in association with heritable disease and should make a significant contribution to our understanding of this major type of mutation.

MAGPIE: accurate pathogenic prediction for multiple variant types using machine learning approach

Identifying pathogenic variants from the vast majority of nucleotide variation remains a challenge. We present a method named Multimodal Annotation Generated Pathogenic Impact Evaluator (MAGPIE) that predicts the pathogenicity of multi-type variants. MAGPIE uses the ClinVar dataset for training and demonstrates superior performance in both the independent test set and multiple orthogonal validation datasets, accurately predicting variant pathogenicity. Notably, MAGPIE performs best in predicting the pathogenicity of rare variants and highly imbalanced datasets. Overall, results underline the robustness of MAGPIE as a valuable tool for predicting pathogenicity in various types of human genome variations. MAGPIE is available at https://github.com/shenlab-genomics/magpie .

VarCards2: an integrated genetic and clinical database for ACMG-AMP variant-interpretation guidelines in the human whole genome

VarCards, an online database, combines comprehensive variant- and gene-level annotation data to streamline genetic counselling for coding variants. Recognising the increasing clinical relevance of non-coding variations, there has been an accelerated development of bioinformatics tools dedicated to interpreting non-coding variations, including single-nucleotide variants and copy number variations. Regrettably, most tools remain as either locally installed databases or command-line tools dispersed across diverse online platforms. Such a landscape poses inconveniences and challenges for genetic counsellors seeking to utilise these resources without advanced bioinformatics expertise. Consequently, we developed VarCards2, which incorporates nearly nine billion artificially generated single-nucleotide variants (including those from mitochondrial DNA) and compiles vital annotation information for genetic counselling based on ACMG-AMP variant-interpretation guidelines. These annotations include (I) functional effects; (II) minor allele frequencies; (III) comprehensive function and pathogenicity predictions covering all potential variants, such as non-synonymous substitutions, non-canonical splicing variants, and non-coding variations and (IV) gene-level information. Furthermore, VarCards2 incorporates 368 820 266 documented short insertions and deletions and 2 773 555 documented copy number variations, complemented by their corresponding annotation and prediction tools. In conclusion, VarCards2, by integrating over 150 variant- and gene-level annotation sources, significantly enhances the efficiency of genetic counselling and can be freely accessed at http://www.genemed.tech/varcards2/.

Isolated Chronic Mucocutaneous Candidiasis due to a Novel Duplication Variant of IL17RC

PURPOSE

Inborn errors of the IL-17A/F-responsive pathway lead to chronic mucocutaneous candidiasis (CMC) as a predominant clinical phenotype, without other significant clinical manifestations apart from mucocutaneous staphylococcal diseases. Among inborn errors affecting IL-17-dependent immunity, autosomal recessive (AR) IL-17RC deficiency is a rare disease with only three kindreds described to date. The lack of an in vitro functional evaluation system of IL17RC variants renders its diagnosis difficult. We sought to characterize a 7-year-old Japanese girl with CMC carrying a novel homozygous duplication variant of IL17RC and establish a simple in vitro system to evaluate the impact of this variant.

METHODS

Flow cytometry, qPCR, RNA-sequencing, and immunoblotting were conducted, and an IL17RC-knockout cell line was established for functional evaluation.

RESULTS

The patient presented with oral and mucocutaneous candidiasis without staphylococcal diseases since the age of 3 months. Genetic analysis showed that the novel duplication variant (Chr3: 9,971,476-9,971,606 dup (+131bp)) involving exon 13 of IL17RC results in a premature stop codon (p.D457Afs*16 or p.D457Afs*17). Our functional evaluation system revealed this duplication to be loss-of-function and enabled discrimination between loss-of-function and neutral IL17RC variants. The lack of response to IL-17A by the patient's SV40-immortalized fibroblasts was restored by introducing WT-IL17RC, suggesting that the genotype identified is responsible for her clinical phenotype.

CONCLUSIONS

The clinical and cellular phenotype of the current case of AR IL-17RC deficiency supports a previous report on this rare disorder. Our newly established evaluation system will be useful for the diagnosis of AR IL-17RC deficiency, providing accurate validation of unknown IL17RC variants.

Genome-wide prediction of pathogenic gain- and loss-of-function variants from ensemble learning of a diverse feature set

Gain-of-function (GOF) variants give rise to increased/novel protein functions whereas loss-of-function (LOF) variants lead to diminished protein function. Experimental approaches for identifying GOF and LOF are generally slow and costly, whilst available computational methods have not been optimized to discriminate between GOF and LOF variants. We have developed LoGoFunc, a machine learning method for predicting pathogenic GOF, pathogenic LOF, and neutral genetic variants, trained on a broad range of gene-, protein-, and variant-level features describing diverse biological characteristics. LoGoFunc outperforms other tools trained solely to predict pathogenicity for identifying pathogenic GOF and LOF variants and is available at https://itanlab.shinyapps.io/goflof/ .

Structural and non-coding variants increase the diagnostic yield of clinical whole genome sequencing for rare diseases

BACKGROUND

Whole genome sequencing is increasingly being used for the diagnosis of patients with rare diseases. However, the diagnostic yields of many studies, particularly those conducted in a healthcare setting, are often disappointingly low, at 25-30%. This is in part because although entire genomes are sequenced, analysis is often confined to in silico gene panels or coding regions of the genome.

METHODS

We undertook WGS on a cohort of 122 unrelated rare disease patients and their relatives (300 genomes) who had been pre-screened by gene panels or arrays. Patients were recruited from a broad spectrum of clinical specialties. We applied a bioinformatics pipeline that would allow comprehensive analysis of all variant types. We combined established bioinformatics tools for phenotypic and genomic analysis with our novel algorithms (SVRare, ALTSPLICE and GREEN-DB) to detect and annotate structural, splice site and non-coding variants.

RESULTS

Our diagnostic yield was 43/122 cases (35%), although 47/122 cases (39%) were considered solved when considering novel candidate genes with supporting functional data into account. Structural, splice site and deep intronic variants contributed to 20/47 (43%) of our solved cases. Five genes that are novel, or were novel at the time of discovery, were identified, whilst a further three genes are putative novel disease genes with evidence of causality. We identified variants of uncertain significance in a further fourteen candidate genes. The phenotypic spectrum associated with RMND1 was expanded to include polymicrogyria. Two patients with secondary findings in FBN1 and KCNQ1 were confirmed to have previously unidentified Marfan and long QT syndromes, respectively, and were referred for further clinical interventions. Clinical diagnoses were changed in six patients and treatment adjustments made for eight individuals, which for five patients was considered life-saving.

CONCLUSIONS

Genome sequencing is increasingly being considered as a first-line genetic test in routine clinical settings and can make a substantial contribution to rapidly identifying a causal aetiology for many patients, shortening their diagnostic odyssey. We have demonstrated that structural, splice site and intronic variants make a significant contribution to diagnostic yield and that comprehensive analysis of the entire genome is essential to maximise the value of clinical genome sequencing.

Genome Sequencing of Consanguineous Family Implicates Ubiquitin-Specific Protease 53 (USP53) Variant in Psychosis/Schizophrenia: Wild-Type Expression in Murine Hippocampal CA 1-3 and Granular Dentate with AMPA Synapse Interactions

Psychosis is a severe mental disorder characterized by abnormal thoughts and perceptions (e.g., hallucinations) occurring quintessentially in schizophrenia and in several other neuropsychiatric disorders. Schizophrenia is widely considered as a neurodevelopmental disorder that onsets during teenage/early adulthood. A multiplex consanguineous Pakistani family was afflicted with severe psychosis and apparent autosomal recessive transmission. The first-cousin parents and five children were healthy, whereas two teenage daughters were severely affected. Structured interviews confirmed the diagnosis of DSM-V schizophrenia. Probands and father underwent next-generation sequencing. All available relatives were subjected to confirmatory Sanger sequencing. Homozygosity mapping and directed a priori filtering identified only one rare variant [MAF < 5(10)-5] at a residue conserved across vertebrates. The variant was a non-catalytic deubiquitinase, USP53 (p.Cys228Arg), predicted in silico as damaging. Genome sequencing did not identify any other potentially pathogenic single nucleotide variant or structural variant. Since the literature on USP53 lacked relevance to mental illness or CNS expression, studies were conducted which revealed USP53 localization in regions of the hippocampus (CA 1-3) and granular dentate. The staining pattern was like that seen with GRIA2/GluA2 and GRIP2 antibodies. All three proteins coimmunoprecipitated. These findings support the glutamate hypothesis of schizophrenia as part of the AMPA-R interactome. If confirmed, USP53 appears to be one of the few Mendelian variants potentially causal to a common-appearing mental disorder that is a rare genetic form of schizophrenia.

Exome sequencing findings in children with annular pancreas

BACKGROUND

Annular pancreas (AP) is a congenital defect of unknown cause in which the pancreas encircles the duodenum. Theories include abnormal migration and rotation of the ventral bud, persistence of ectopic pancreatic tissue, and inappropriate fusion of the ventral and dorsal buds before rotation. The few reported familial cases suggest a genetic contribution.

METHODS

We conducted exome sequencing in 115 affected infants from the California birth defects registry.

RESULTS

Seven cases had a single heterozygous missense variant in IQGAP1, five of them with CADD scores >20; seven other infants had a single heterozygous missense variant in NRCAM, five of them with CADD scores >20. We also looked at genes previously associated with AP and found two rare heterozygous missense variants, one each in PDX1 and FOXF1.

CONCLUSION

IQGAP1 and NRCAM are crucial in cell polarization and migration. Mutations result in decreased motility which could possibly cause the ventral bud to not migrate normally. To our knowledge, this is the first study reporting a possible association for IQGAP1 and NRCAM with AP. Our findings of rare genetic variants involved in cell migration in 15% of our population raise the possibility that AP may be related to abnormal cell migration.

A knowledge graph approach to predict and interpret disease-causing gene interactions

BACKGROUND

Understanding the impact of gene interactions on disease phenotypes is increasingly recognised as a crucial aspect of genetic disease research. This trend is reflected by the growing amount of clinical research on oligogenic diseases, where disease manifestations are influenced by combinations of variants on a few specific genes. Although statistical machine-learning methods have been developed to identify relevant genetic variant or gene combinations associated with oligogenic diseases, they rely on abstract features and black-box models, posing challenges to interpretability for medical experts and impeding their ability to comprehend and validate predictions. In this work, we present a novel, interpretable predictive approach based on a knowledge graph that not only provides accurate predictions of disease-causing gene interactions but also offers explanations for these results.

RESULTS

We introduce BOCK, a knowledge graph constructed to explore disease-causing genetic interactions, integrating curated information on oligogenic diseases from clinical cases with relevant biomedical networks and ontologies. Using this graph, we developed a novel predictive framework based on heterogenous paths connecting gene pairs. This method trains an interpretable decision set model that not only accurately predicts pathogenic gene interactions, but also unveils the patterns associated with these diseases. A unique aspect of our approach is its ability to offer, along with each positive prediction, explanations in the form of subgraphs, revealing the specific entities and relationships that led to each pathogenic prediction.

CONCLUSION

Our method, built with interpretability in mind, leverages heterogenous path information in knowledge graphs to predict pathogenic gene interactions and generate meaningful explanations. This not only broadens our understanding of the molecular mechanisms underlying oligogenic diseases, but also presents a novel application of knowledge graphs in creating more transparent and insightful predictors for genetic research.

Isolated chronic mucocutaneous candidiasis due to a novel duplication variant of IL17RC

Purpose

Inborn errors of the IL-17A/F-responsive pathway lead to chronic mucocutaneous candidiasis (CMC) as a predominant clinical phenotype, without other significant clinical manifestations apart from mucocutaneous staphylococcal diseases. Amongst inborn errors affecting IL-17-dependent immunity, autosomal recessive (AR) IL-17RC deficiency is a rare disease with only three kindreds described to date. The lack of an in vitro functional evaluation system of IL17RC variants renders its diagnosis difficult. We sought to characterize a seven-year-old Japanese girl with CMC carrying a novel homozygous duplication variant of IL17RC and establish a simple in vitro system to evaluate the impact of this variant.

Methods

Flow cytometry, qPCR, RNA-sequencing, and immunoblotting were conducted, and an IL17RC-knockout cell line was established for functional evaluation.

Results

The patient presented with oral and mucocutaneous candidiasis without staphylococcal diseases since the age of three months. Genetic analysis showed that the novel duplication variant (Chr3: 9,971,476-9,971,606 dup (+ 131bp)) involving exon 13 of IL17RC results in a premature stop codon (p.D457Afs*16 or p.D457Afs*17). Our functional evaluation system revealed this duplication to be loss-of-function and enabled discrimination between loss-of-function and neutral IL17RC variants. The lack of response to IL-17A by the patient's SV40-immortalized fibroblasts was restored by introducing WT-IL17RC, suggesting that the genotype identified is responsible for her clinical phenotype.

Conclusions

The clinical and cellular phenotype of the current case of AR IL-17RC deficiency supports a previous report on this rare disorder. Our newly established evaluation system will be useful for diagnosis of AR IL-17RC deficiency, providing accurate validation of unknown IL17RC variants.

Deep structured learning for variant prioritization in Mendelian diseases

Effective computer-aided or automated variant evaluations for monogenic diseases will expedite clinical diagnostic and research efforts of known and novel disease-causing genes. Here we introduce MAVERICK: a Mendelian Approach to Variant Effect pRedICtion built in Keras. MAVERICK is an ensemble of transformer-based neural networks that can classify a wide range of protein-altering single nucleotide variants (SNVs) and indels and assesses whether a variant would be pathogenic in the context of dominant or recessive inheritance. We demonstrate that MAVERICK outperforms all other major programs that assess pathogenicity in a Mendelian context. In a cohort of 644 previously solved patients with Mendelian diseases, MAVERICK ranks the causative pathogenic variant within the top five variants in over 95% of cases. Seventy-six percent of cases were solved by the top-ranked variant. MAVERICK ranks the causative pathogenic variant in hitherto novel disease genes within the first five candidate variants in 70% of cases. MAVERICK has already facilitated the identification of a novel disease gene causing a degenerative motor neuron disease. These results represent a significant step towards automated identification of causal variants in patients with Mendelian diseases.

Germline loss-of-function PAM variants are enriched in subjects with pituitary hypersecretion

Introduction

Pituitary adenomas (PAs) are common, usually benign tumors of the anterior pituitary gland which, for the most part, have no known genetic cause. PAs are associated with major clinical effects due to hormonal dysregulation and tumoral impingement on vital brain structures. PAM encodes a multifunctional protein responsible for the essential C-terminal amidation of secreted peptides.

Methods

Following the identification of a loss-of-function variant (p.Arg703Gln) in the peptidylglycine a-amidating monooxygenase (PAM) gene in a family with pituitary gigantism, we investigated 299 individuals with sporadic PAs and 17 familial isolated PA kindreds for PAM variants. Genetic screening was performed by germline and tumor sequencing and germline copy number variation (CNV) analysis.

Results

In germline DNA, we detected seven heterozygous, likely pathogenic missense, truncating, and regulatory SNVs. These SNVs were found in sporadic subjects with growth hormone excess (p.Gly552Arg and p.Phe759Ser), pediatric Cushing disease (c.-133T>C and p.His778fs), or different types of PAs (c.-361G>A, p.Ser539Trp, and p.Asp563Gly). The SNVs were functionally tested in vitro for protein expression and trafficking by Western blotting, splicing by minigene assays, and amidation activity in cell lysates and serum samples. These analyses confirmed a deleterious effect on protein expression and/or function. By interrogating 200,000 exomes from the UK Biobank, we confirmed a significant association of the PAM gene and rare PAM SNVs with diagnoses linked to pituitary gland hyperfunction.

Conclusion

The identification of PAM as a candidate gene associated with pituitary hypersecretion opens the possibility of developing novel therapeutics based on altering PAM function.

Encephalitis and poor neuronal death-mediated control of herpes simplex virus in human inherited RIPK3 deficiency

Inborn errors of TLR3-dependent type I IFN immunity in cortical neurons underlie forebrain herpes simplex virus-1 (HSV-1) encephalitis (HSE) due to uncontrolled viral growth and subsequent cell death. We report an otherwise healthy patient with HSE who was compound heterozygous for nonsense (R422*) and frameshift (P493fs9*) RIPK3 variants. Receptor-interacting protein kinase 3 (RIPK3) is a ubiquitous cytoplasmic kinase regulating cell death outcomes, including apoptosis and necroptosis. In vitro, the R422* and P493fs9* RIPK3 proteins impaired cellular apoptosis and necroptosis upon TLR3, TLR4, or TNFR1 stimulation and ZBP1/DAI-mediated necroptotic cell death after HSV-1 infection. The patient's fibroblasts displayed no detectable RIPK3 expression. After TNFR1 or TLR3 stimulation, the patient's cells did not undergo apoptosis or necroptosis. After HSV-1 infection, the cells supported excessive viral growth despite normal induction of antiviral IFN-β and IFN-stimulated genes (ISGs). This phenotype was, nevertheless, rescued by application of exogenous type I IFN. The patient's human pluripotent stem cell (hPSC)-derived cortical neurons displayed impaired cell death and enhanced viral growth after HSV-1 infection, as did isogenic RIPK3-knockout hPSC-derived cortical neurons. Inherited RIPK3 deficiency therefore confers a predisposition to HSE by impairing the cell death-dependent control of HSV-1 in cortical neurons but not their production of or response to type I IFNs.

Identifying high-impact variants and genes in exomes of Ashkenazi Jewish inflammatory bowel disease patients

Inflammatory bowel disease (IBD) is a group of chronic digestive tract inflammatory conditions whose genetic etiology is still poorly understood. The incidence of IBD is particularly high among Ashkenazi Jews. Here, we identify 8 novel and plausible IBD-causing genes from the exomes of 4453 genetically identified Ashkenazi Jewish IBD cases (1734) and controls (2719). Various biological pathway analyses are performed, along with bulk and single-cell RNA sequencing, to demonstrate the likely physiological relatedness of the novel genes to IBD. Importantly, we demonstrate that the rare and high impact genetic architecture of Ashkenazi Jewish adult IBD displays significant overlap with very early onset-IBD genetics. Moreover, by performing biobank phenome-wide analyses, we find that IBD genes have pleiotropic effects that involve other immune responses. Finally, we show that polygenic risk score analyses based on genome-wide high impact variants have high power to predict IBD susceptibility.

Allele-dependent interaction of LRRK2 and NOD2 in leprosy

Leprosy, caused by Mycobacterium leprae, rarely affects children younger than 5 years. Here, we studied a multiplex leprosy family that included monozygotic twins aged 22 months suffering from paucibacillary leprosy. Whole genome sequencing identified three amino acid mutations previously associated with Crohn’s disease and Parkinson’s disease as candidate variants for early onset leprosy: LRRK2 N551K, R1398H and NOD2 R702W. In genome-edited macrophages, we demonstrated that cells expressing the LRRK2 mutations displayed reduced apoptosis activity following mycobacterial challenge independently of NOD2. However, employing co-immunoprecipitation and confocal microscopy we showed that LRRK2 and NOD2 proteins interacted in RAW cells and monocyte-derived macrophages, and that this interaction was substantially reduced for the NOD2 R702W mutation. Moreover, we observed a joint effect of LRRK2 and NOD2 variants on Bacillus Calmette-Guérin (BCG)-induced respiratory burst, NF-κB activation and cytokine/chemokine secretion with a strong impact for the genotypes found in the twins consistent with a role of the identified mutations in the development of early onset leprosy.

Human IL-23 is essential for IFN-γ-dependent immunity to mycobacteria

Patients with autosomal recessive (AR) IL-12p40 or IL-12Rβ1 deficiency display Mendelian susceptibility to mycobacterial disease (MSMD) due to impaired IFN-γ production and, less commonly, chronic mucocutaneous candidiasis (CMC) due to impaired IL-17A/F production. We report six patients from four kindreds with AR IL-23R deficiency. These patients are homozygous for one of four different loss-of-function IL23R variants. All six patients have a history of MSMD, but only two suffered from CMC. We show that IL-23 induces IL-17A only in MAIT cells, possibly contributing to the incomplete penetrance of CMC in patients unresponsive to IL-23. By contrast, IL-23 is required for both baseline and Mycobacterium-inducible IFN-γ immunity in both Vδ2+ γδ T and MAIT cells, probably contributing to the higher penetrance of MSMD in these patients. Human IL-23 appears to contribute to IL-17A/F-dependent immunity to Candida in a single lymphocyte subset but is required for IFN-γ-dependent immunity to Mycobacterium in at least two lymphocyte subsets.

Inborn errors of OAS-RNase L in SARS-CoV-2-related multisystem inflammatory syndrome in children

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

Human IRF1 governs macrophagic IFN-γ immunity to mycobacteria

Inborn errors of human IFN-γ-dependent macrophagic immunity underlie mycobacterial diseases, whereas inborn errors of IFN-α/β-dependent intrinsic immunity underlie viral diseases. Both types of IFNs induce the transcription factor IRF1. We describe unrelated children with inherited complete IRF1 deficiency and early-onset, multiple, life-threatening diseases caused by weakly virulent mycobacteria and related intramacrophagic pathogens. These children have no history of severe viral disease, despite exposure to many viruses, including SARS-CoV-2, which is life-threatening in individuals with impaired IFN-α/β immunity. In leukocytes or fibroblasts stimulated in vitro, IRF1-dependent responses to IFN-γ are, both quantitatively and qualitatively, much stronger than those to IFN-α/β. Moreover, IRF1-deficient mononuclear phagocytes do not control mycobacteria and related pathogens normally when stimulated with IFN-γ. By contrast, IFN-α/β-dependent intrinsic immunity to nine viruses, including SARS-CoV-2, is almost normal in IRF1-deficient fibroblasts. Human IRF1 is essential for IFN-γ-dependent macrophagic immunity to mycobacteria, but largely redundant for IFN-α/β-dependent antiviral immunity.

Severe high-molecular-weight kininogen deficiency: clinical characteristics, deficiency-causing KNG1 variants, and estimated prevalence

BACKGROUND

Severe high-molecular-weight kininogen (HK) deficiency is a poorly studied autosomal recessive contact system defect caused by pathogenic, biallelic KNG1 variants.

AIM

We performed the first comprehensive analysis of diagnostic, clinical, genetic, and epidemiological aspects of HK deficiency.

METHODS

We collected clinical information and blood samples from a newly detected HK-deficient individual and from published cases identified by a systematic literature review. Activity and antigen levels of coagulation factors were determined. Genetic analyses of KNG1 and KLKB1 were performed by Sanger sequencing. The frequency of HK deficiency was estimated considering truncating KNG1 variants from GnomAD.

RESULTS

We identified 48 cases of severe HK deficiency (41 families), of these 47 have been previously published (n = 19 from gray literature). We genotyped 3 cases and critically appraised 10 studies with genetic data. Ten HK deficiency-causing variants (one new) were identified. All of them were truncating mutations, whereas the only known HK amino acid substitution with a relevant phenotype instead causes hereditary angioedema. Conservative estimates suggest an overall prevalence of severe HK deficiency of approximately one case per 8 million population, slightly higher in Africans. Individuals with HK deficiency appeared asymptomatic and had decreased levels of prekallikrein and factor XI, which could lead to misdiagnosis.

CONCLUSION

HK deficiency is a rare condition with only few known pathogenic variants. It has an apparently good prognosis but is prone to misdiagnosis. Our understanding of its clinical implications is still limited, and an international prekallikrein and HK deficiency registry is being established to fill this knowledge gap.

Fulminant Viral Hepatitis in Two Siblings with Inherited IL-10RB Deficiency

Fulminant viral hepatitis (FVH) caused by hepatitis A virus (HAV) is a life-threatening disease that typically strikes otherwise healthy individuals. The only known genetic etiology of FVH is inherited IL-18BP deficiency, which unleashes IL-18-dependent lymphocyte cytotoxicity and IFN-γ production. We studied two siblings who died from a combination of early-onset inflammatory bowel disease (EOIBD) and FVH due to HAV. The sibling tested was homozygous for the W100G variant of IL10RB previously described in an unrelated patient with EOIBD. We show here that the out-of-frame IL10RB variants seen in other EOIBD patients disrupt cellular responses to IL-10, IL-22, IL-26, and IFN-λs in overexpression conditions and in homozygous cells. By contrast, the impact of in-frame disease-causing variants varies between cases. When overexpressed, the W100G variant impairs cellular responses to IL-10, but not to IL-22, IL-26, or IFN-λ1, whereas cells homozygous for W100G do not respond to IL-10, IL-22, IL-26, or IFN-λ1. As IL-10 is a potent antagonist of IFN-γ in phagocytes, these findings suggest that the molecular basis of FVH in patients with IL-18BP or IL-10RB deficiency may involve excessive IFN-γ activity during HAV infections of the liver. Inherited IL-10RB deficiency, and possibly inherited IL-10 and IL-10RA deficiencies, confer a predisposition to FVH, and patients with these deficiencies should be vaccinated against HAV and other liver-tropic viruses.

Germline loss-of-function PAM variants are enriched in subjects with pituitary hypersecretion

Pituitary adenomas (PAs) are common, usually benign tumors of the anterior pituitary gland which, for the most part, have no known genetic cause. PAs are associated with major clinical effects due to hormonal dysregulation and tumoral impingement on vital brain structures. Following the identification of a loss-of-function variant (p.Arg703Gln) in the PAM gene in a family with pituitary gigantism, we investigated 299 individuals with sporadic PAs and 17 familial isolated pituitary adenomas kindreds for PAM variants. PAM encodes a multifunctional protein responsible for the essential C-terminal amidation of secreted peptides. Genetic screening was performed by germline and tumor sequencing and germline copy number variation (CNV) analysis. No germline CNVs or somatic single nucleotide variants (SNVs) were identified. We detected seven likely pathogenic heterozygous missense, truncating, and regulatory SNVs. These SNVs were found in sporadic subjects with GH excess (p.Gly552Arg and p.Phe759Ser), pediatric Cushing disease (c.-133T>C and p.His778fs), or with different types of PAs (c.-361G>A, p.Ser539Trp, and p.Asp563Gly). The SNVs were functionally tested in vitro for protein expression and trafficking by Western blotting, for splicing by minigene assays, and for amidation activity in cell lysates and serum samples. These analyses confirmed a deleterious effect on protein expression and/or function. By interrogating 200,000 exomes from the UK Biobank, we confirmed a significant association of the PAM gene and rare PAM SNVs to diagnoses linked to pituitary gland hyperfunction. Identification of PAM as a candidate gene associated with pituitary hypersecretion opens the possibility of developing novel therapeutics based on altering PAM function.

Investigation of UTR Variants by Computational Approaches Reveal Their Functional Significance in PRKCI Gene Regulation

Single nucleotide polymorphisms (SNPs) are associated with many diseases including neurological disorders, heart diseases, diabetes, and different types of cancers. In the context of cancer, the variations within non-coding regions, including UTRs, have gained utmost importance. In gene expression, translational regulation is as important as transcriptional regulation for the normal functioning of cells; modification in normal functions can be associated with the pathophysiology of many diseases. UTR-localized SNPs in the PRKCI gene were evaluated using the PolymiRTS, miRNASNP, and MicroSNIper for association with miRNAs. Furthermore, the SNPs were subjected to analysis using GTEx, RNAfold, and PROMO. The genetic intolerance to functional variation was checked through GeneCards. Out of 713 SNPs, a total of thirty-one UTR SNPs (three in 3' UTR region and twenty-nine in 5' UTR region) were marked as ≤2b by RegulomeDB. The associations of 23 SNPs with miRNAs were found. Two SNPs, rs140672226 and rs2650220, were significantly linked with expression in the stomach and esophagus mucosa. The 3' UTR SNPs rs1447651774 and rs115170199 and the 5' UTR region variants rs778557075, rs968409340, and 750297755 were predicted to destabilize the mRNA structure with substantial change in free energy (∆G). Seventeen variants were predicted to have linkage disequilibrium with various diseases. The SNP rs542458816 in 5' UTR was predicted to put maximum influence on transcription factor binding sites. Gene damage index(GDI) and loss of function (o:e) ratio values for PRKCI suggested that the gene is not tolerant to loss of function variants. Our results highlight the effects of 3' and 5' UTR SNP on miRNA, transcription and translation of PRKCI. These analyses suggest that these SNPs can have substantial functional importance in the PRKCI gene. Future experimental validation could provide further basis for the diagnosis and therapeutics of various diseases.

Adaptive Evolution of the OAS Gene Family Provides New Insights into the Antiviral Ability of Laurasiatherian Mammals

Many mammals risk damage from virus invasion due to frequent environmental changes. The oligoadenylate synthesis (OAS) gene family, which is an important component of the immune system, provides an essential response to the antiviral activities of interferons by regulating immune signal pathways. However, little is known about the evolutionary characteristics of OASs in Laurasiatherian mammals. Here, we examined the evolution of the OAS genes in 64 mammals to explore the accompanying molecular mechanisms of the antiviral ability of Laurasiatherian mammals living in different environments. We found that OAS2 and OAS3 were found to be pseudogenes in Odontoceti species. This may be related to the fact that they live in water. Some Antilopinae, Caprinae, and Cervidae species lacked the OASL gene, which may be related to their habitats being at higher altitudes. The OASs had a high number of positive selection sites in Cetartiodactyla, which drove the expression of strong antiviral ability. The OAS gene family evolved in Laurasiatherian mammals at different rates and was highly correlated with the species' antiviral ability. The gene evolution rate in Cetartiodactyla was significantly higher than that in the other orders. Compared to other species of the Carnivora family, the higher selection pressure on the OAS gene and the absence of positive selection sites in Canidae may be responsible for its weak resistance to rabies virus. The OAS gene family was relatively conserved during evolution. Conserved genes are able to provide better maintenance of gene function. The rate of gene evolution and the number of positively selected sites combine to influence the resistance of a species to viruses. The positive selection sites demonstrate the adaptive evolution of the OAS gene family to the environment. Adaptive evolution combined with conserved gene function improves resistance to viruses. Our findings offer insights into the molecular and functional evolution of the antiviral ability of Laurasian mammals.

Partial loss-of-function mutations in GINS4 lead to NK cell deficiency with neutropenia

Human NK cell deficiency (NKD) is a primary immunodeficiency in which the main clinically relevant immunological defect involves missing or dysfunctional NK cells. Here, we describe a familial NKD case in which 2 siblings had a substantive NKD and neutropenia in the absence of other immune system abnormalities. Exome sequencing identified compound heterozygous variants in Go-Ichi-Ni-San (GINS) complex subunit 4 (GINS4, also known as SLD5), an essential component of the human replicative helicase, which we demonstrate to have a damaging impact upon the expression and assembly of the GINS complex. Cells derived from affected individuals and a GINS4-knockdown cell line demonstrate delayed cell cycle progression, without signs of improper DNA synthesis or increased replication stress. By modeling partial GINS4 depletion in differentiating NK cells in vitro, we demonstrate the causal relationship between the genotype and the NK cell phenotype, as well as a cell-intrinsic defect in NK cell development. Thus, biallelic partial loss-of-function mutations in GINS4 define a potentially novel disease-causing gene underlying NKD with neutropenia. Together with the previously described mutations in other helicase genes causing NKD, and with the mild defects observed in other human cells, these variants underscore the importance of this pathway in NK cell biology.

Impaired IL-23-dependent induction of IFN-γ underlies mycobacterial disease in patients with inherited TYK2 deficiency

Human cells homozygous for rare loss-of-expression (LOE) TYK2 alleles have impaired, but not abolished, cellular responses to IFN-α/β (underlying viral diseases in the patients) and to IL-12 and IL-23 (underlying mycobacterial diseases). Cells homozygous for the common P1104A TYK2 allele have selectively impaired responses to IL-23 (underlying isolated mycobacterial disease). We report three new forms of TYK2 deficiency in six patients from five families homozygous for rare TYK2 alleles (R864C, G996R, G634E, or G1010D) or compound heterozygous for P1104A and a rare allele (A928V). All these missense alleles encode detectable proteins. The R864C and G1010D alleles are hypomorphic and loss-of-function (LOF), respectively, across signaling pathways. By contrast, hypomorphic G996R, G634E, and A928V mutations selectively impair responses to IL-23, like P1104A. Impairment of the IL-23-dependent induction of IFN-γ is the only mechanism of mycobacterial disease common to patients with complete TYK2 deficiency with or without TYK2 expression, partial TYK2 deficiency across signaling pathways, or rare or common partial TYK2 deficiency specific for IL-23 signaling.

Performance evaluation of differential splicing analysis methods and splicing analytics platform construction

A proportion of previously defined benign variants or variants of uncertain significance in humans, which are challenging to identify, may induce an abnormal splicing process. An increasing number of methods have been developed to predict splicing variants, but their performance has not been completely evaluated using independent benchmarks. Here, we manually sourced ∼50 000 positive/negative splicing variants from > 8000 studies and selected the independent splicing variants to evaluate the performance of prediction methods. These methods showed different performances in recognizing splicing variants in donor and acceptor regions, reminiscent of different weight coefficient applications to predict novel splicing variants. Of these methods, 66.67% exhibited higher specificities than sensitivities, suggesting that more moderate cut-off values are necessary to distinguish splicing variants. Moreover, the high correlation and consistent prediction ratio validated the feasibility of integration of the splicing prediction method in identifying splicing variants. We developed a splicing analytics platform called SPCards, which curates splicing variants from publications and predicts splicing scores of variants in genomes. SPCards also offers variant-level and gene-level annotation information, including allele frequency, non-synonymous prediction and comprehensive functional information. SPCards is suitable for high-throughput genetic identification of splicing variants, particularly those located in non-canonical splicing regions.

Evidence of a genetically driven metabolomic signature in actively inflamed Crohn's disease

Crohn's disease (CD) is characterised by chronic inflammation. We aimed to identify a relationship between plasma inflammatory metabolomic signature and genomic data in CD using blood plasma metabolic profiles. Proton NMR spectroscopy were achieved for 228 paediatric CD patients. Regression (OPLS) modelling and machine learning (ML) approaches were independently applied to establish the metabolic inflammatory signature, which was correlated against gene-level pathogenicity scores generated for all patients and functional enrichment was analysed. OPLS modelling of metabolomic spectra from unfasted patients revealed distinctive shifts in plasma metabolites corresponding to regions of the spectrum assigned to N-acetyl glycoprotein, glycerol and phenylalanine that were highly correlated (R² = 0.62) with C-reactive protein levels. The same metabolomic signature was independently identified using ML to predict patient inflammation status. Correlation of the individual peaks comprising this metabolomic signature of inflammation with pathogenic burden across 15,854 unselected genes identified significant enrichment for genes functioning within 'intrinsic component of membrane' (p = 0.003) and 'inflammatory bowel disease (IBD)' (p = 0.003). The seven genes contributing IBD enrichment are critical regulators of pro-inflammatory signaling. Overall, a metabolomic signature of inflammation can be detected from blood plasma in CD. This signal is correlated with pathogenic mutation in pro-inflammatory immune response genes.

An accurate prediction model of digenic interaction for estimating pathogenic gene pairs of human diseases

•

Numerous pathogenic interactions are yet to be revealed efficiently due to dimension burden.

•

Existing methods are underpowered and inaccurate to estimate pathogenic interactions.

•

We developed an accurate bioinformatic method to predict digenic interaction effects based on gene-level features.

•

The new method created a valuable resource of genome-wide pathogenic digenic interactions.

•

We found that known causal genes in Mendelian and Oligogenic diseases may be enriched with interactive effects for the first time.

Increasing evidence shows that genetic interaction across the entire genome may explain a non-trivial fraction of genetic diseases. Digenic interaction is the simplest manifestation of genetic interaction among genes. However, systematic exploration of digenic interactive effects on the whole genome is often discouraged by the high dimension burden. Thus, numerous digenic interactions are yet to be identified for many diseases. Here, we propose a Digenic Interaction Effect Predictor (DIEP), an accurate machine-learning approach to identify the genome-wide pathogenic coding gene pairs with digenic interaction effects. This approach achieved high accuracy and sensitivity in independent testing datasets, outperforming another gene-level digenic predictor (DiGePred). DIEP was also able to discriminate digenic interaction effect from bi-locus effects dual molecular diagnosis (pseudo-digenic). Using DIEP, we provided a valuable resource of genome-wide digenic interactions and demonstrated the enrichment of the digenic interaction effect in Mendelian and Oligogenic diseases. Therefore, DIEP will play a useful role in facilitating the genomic mapping of interactive causal genes for human diseases.

Whole genome sequences discriminate hereditary hemorrhagic telangiectasia phenotypes by non-HHT deleterious DNA variation

The abnormal vascular structures of hereditary hemorrhagic telangiectasia (HHT) often cause severe anemia due to recurrent hemorrhage, but HHT causal genes do not predict the severity of hematological complications. We tested for chance inheritance and clinical associations of rare deleterious variants in which loss-of-function causes bleeding or hemolytic disorders in the general population. In double-blinded analyses, all 104 patients with HHT from a single reference center recruited to the 100 000 Genomes Project were categorized on new MALO (more/as-expected/less/opposite) sub-phenotype severity scales, and whole genome sequencing data were tested for high impact variants in 75 HHT-independent genes encoding coagulation factors, or platelet, hemoglobin, erythrocyte enzyme, and erythrocyte membrane constituents. Rare variants (all gnomAD allele frequencies <0.003) were identified in 56 (75%) of these 75 HHT-unrelated genes. Deleteriousness assignments by Combined Annotation Dependent Depletion (CADD) scores >15 were supported by gene-level mutation significance cutoff scores. CADD >15 variants were identified in 38/104 (36.5%) patients with HHT, found for 1 in 10 patients within platelet genes; 1 in 8 within coagulation genes; and 1 in 4 within erythrocyte hemolytic genes. In blinded analyses, patients with greater hemorrhagic severity that had been attributed solely to HHT vessels had more CADD-deleterious variants in platelet (Spearman ρ = 0.25; P = .008) and coagulation (Spearman ρ = 0.21; P = .024) genes. However, the HHT cohort had 60% fewer deleterious variants in platelet and coagulation genes than expected (Mann-Whitney test P = .021). In conclusion, patients with HHT commonly have rare variants in genes of relevance to their phenotype, offering new therapeutic targets and opportunities for informed, personalized medicine strategies.

Genomic characterization of lymphomas in patients with inborn errors of immunity

Inborn errors of immunity-associated lymphomas are characterized by distinct clinical features and genetic signatures.

Both germline and somatic alterations contribute to lymphomagenesis in patients with inborn errors of immunity.

Patients with inborn errors of immunity (IEI) have a higher risk of developing cancer, especially lymphoma. However, the molecular basis for IEI-related lymphoma is complex and remains elusive. Here, we perform an in-depth analysis of lymphoma genomes derived from 23 IEI patients. We identified and validated disease-causing or -associated germline mutations in 14 of 23 patients involving ATM, BACH2, BLM, CD70, G6PD, NBN, PIK3CD, PTEN, and TNFRSF13B. Furthermore, we profiled somatic mutations in the lymphoma genome and identified 8 genes that were mutated at a significantly higher level in IEI-associated diffuse large B-cell lymphomas (DLBCLs) than in non-IEI DLBCLs, such as BRCA2, NCOR1, KLF2, FAS, CCND3, and BRWD3. The latter, BRWD3, is furthermore preferentially mutated in tumors of a subgroup of activated phosphoinositide 3-kinase δ syndrome patients. We also identified 5 genomic mutational signatures, including 2 DNA repair deficiency-related signatures, in IEI-associated lymphomas and a strikingly high number of inter- and intrachromosomal structural variants in the tumor genome of a Bloom syndrome patient. In summary, our comprehensive genomic characterization of lymphomas derived from patients with rare genetic disorders expands our understanding of lymphomagenesis and provides new insights for targeted therapy.

Large-scale discovery of novel neurodevelopmental disorder-related genes through a unified analysis of single-nucleotide and copy number variants

Background

Previous large-scale studies of de novo variants identified a number of genes associated with neurodevelopmental disorders (NDDs); however, it was also predicted that many NDD-associated genes await discovery. Such genes can be discovered by integrating copy number variants (CNVs), which have not been fully considered in previous studies, and increasing the sample size.

Methods

We first constructed a model estimating the rates of de novo CNVs per gene from several factors such as gene length and number of exons. Second, we compiled a comprehensive list of de novo single-nucleotide variants (SNVs) in 41,165 individuals and de novo CNVs in 3675 individuals with NDDs by aggregating our own and publicly available datasets, including denovo-db and the Deciphering Developmental Disorders study data. Third, summing up the de novo CNV rates that we estimated and SNV rates previously established, gene-based enrichment of de novo deleterious SNVs and CNVs were assessed in the 41,165 cases. Significantly enriched genes were further prioritized according to their similarity to known NDD genes using a deep learning model that considers functional characteristics (e.g., gene ontology and expression patterns).

Results

We identified a total of 380 genes achieving statistical significance (5% false discovery rate), including 31 genes affected by de novo CNVs. Of the 380 genes, 52 have not previously been reported as NDD genes, and the data of de novo CNVs contributed to the significance of three genes (GLTSCR1, MARK2, and UBR3). Among the 52 genes, we reasonably excluded 18 genes [a number almost identical to the theoretically expected false positives (i.e., 380 × 0.05 = 19)] given their constraints against deleterious variants and extracted 34 “plausible” candidate genes. Their validity as NDD genes was consistently supported by their similarity in function and gene expression patterns to known NDD genes. Quantifying the overall similarity using deep learning, we identified 11 high-confidence (> 90% true-positive probabilities) candidate genes: HDAC2, SUPT16H, HECTD4, CHD5, XPO1, GSK3B, NLGN2, ADGRB1, CTR9, BRD3, and MARK2.

Conclusions

We identified dozens of new candidates for NDD genes. Both the methods and the resources developed here will contribute to the further identification of novel NDD-associated genes.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13073-022-01042-w.

Exome sequencing identifies variants in infants with sacral agenesis

BACKGROUND

Sacral agenesis (SA) consists of partial or complete absence of the caudal end of the spine and often presents with additional birth defects. Several studies have examined gene variants for syndromic forms of SA, but only one has examined exomes of children with non-syndromic SA.

METHODS

Using buccal cell specimens from families of children with non-syndromic SA, exomes of 28 child-parent trios (eight with and 20 without a maternal diagnosis of pregestational diabetes) and two child-father duos (neither with diagnosis of maternal pregestational diabetes) were exome sequenced.

RESULTS

Three children had heterozygous missense variants in ID1 (Inhibitor of DNA Binding 1), with CADD scores >20 (top 1% of deleterious variants in the genome); two children inherited the variant from their fathers and one from the child's mother. Rare missense variants were also detected in PDZD2 (PDZ Domain Containing 2; N = 1) and SPTBN5 (Spectrin Beta, Non-erythrocytic 5; N = 2), two genes previously suggested to be associated with SA etiology. Examination of variants with autosomal recessive and X-linked recessive inheritance identified five and two missense variants, respectively. Compound heterozygous variants were identified in several genes. In addition, 12 de novo variants were identified, all in different genes in different children.

CONCLUSIONS

To our knowledge, this is the first study reporting a possible association between ID1 and non-syndromic SA. Although maternal pregestational diabetes has been strongly associated with SA, the missense variants in ID1 identified in two of three children were paternally inherited. These findings add to the knowledge of gene variants associated with non-syndromic SA and provide data for future studies.