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

SwissGenVar: A Platform for Clinical-Grade Interpretation of Genetic Variants to Foster Personalized Healthcare in Switzerland

Large-scale next-generation sequencing (NGS) germline testing is technically feasible today, but variant interpretation represents a major bottleneck in analysis workflows. This includes extensive variant prioritization, annotation, and time-consuming evidence curation. The scale of the interpretation problem is massive, and variants of uncertain significance (VUSs) are a challenge to personalized medicine. This challenge is further compounded by the complexity and heterogeneity of the standards used to describe genetic variants and the associated phenotypes when searching for relevant information to support clinical decision making. To address this, all five Swiss academic institutions for Medical Genetics joined forces with the Swiss Institute of Bioinformatics (SIB) to create SwissGenVar as a user-friendly nationwide repository and sharing platform for genetic variant data generated during routine diagnostic procedures and research sequencing projects. Its aim is to provide a protected environment for expert evidence sharing about individual variants to harmonize and upscale their significance interpretation at the clinical grade according to international standards. To corroborate the clinical assessment, the variant-related data will be combined with consented high-quality clinical information. Broader visibility will be achieved by interfacing with international databases, thus supporting global initiatives in personalized healthcare.

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.

Genetic analysis of nephrogenic diabetes insipidus patients: A study on the Iranian population

INTRODUCTION

Nephrogenic diabetes insipidus (NDI) is a rare genetic disease that causes water imbalance. The kidneys play a crucial role in regulating body fluids by controlling water balance through urine excretion. This highlights their essential function in managing the body's water levels, but individuals with NDI may have excess urine production (polyuria), that leads to excessive thirst (polydipsia). Untreated affected individuals may exhibit poor feeding and failure to thrive. This disease is caused by mutations in the AVPR2 and the AQP2 genes which have the X-linked and autosomal recessive/dominant inheritance, respectively. Both of these genes are expressed in the kidney.

METHODS

Twelve Iranian patients from 10 consanguineous families were studied in this project. DNA was extracted from the whole blood samples of the patients and their parents. All coding exons and exon-intron boundaries of the AVPR2 and AQP2 genes were sequenced in the affected individuals, and the identified variants were investigated in the parents. All variants were analyzed according to the ACMG (American College of Medical Genetics and Genomics) guidelines.

RESULTS

In this study, 6 different mutations were identified in the patients, including 5 in the AQP2 gene (c.439G>A, c.538G>A, c.140C>T, c.450T>A, and the novel c.668T>C) and 1 in the AVPR2 gene (c.337C>T) in the present study.

DISCUSSION

As expected, all the detected mutations in this study were missense. According to the ACMG guideline, the identified mutations were categorized as pathogenic or likely pathogenic. Unlike previous studies which showed more than 90% of mutations were in the AVPR2 gene, and only less than 10% of the mutations were in the AQP2 gene, it was found that more than 90% of our identified mutations located in the AQP2 gene, and only one mutation was observed in the AVPR2 gene, which seems it may be a result of the high rate of consanguineous marriages in the Iranian population. We observed genotype-phenotype correlation in some of our affected individuals, and some of the mutations were observed in unrelated families from same ethnicity which could be suggestive of a founder mutation.

Genetics and clinical phenotypes in common variable immunodeficiency

Common variable immunodeficiency (CVID) is one of the most common symptomatic groups of inborn errors of immunity. In addition to infections resulting from insufficient levels of immune globulins and antibodies, many patients develop inflammatory or autoimmune conditions, which are associated with increased mortality. This aspect of CVID has been the focus of many studies, and dissecting the clinical phenotypes of CVID, has had the goal of providing biomarkers to identify these subjects, potentially at the time of diagnosis. With the application of whole exome (WES) and whole genome analyses, an increasing number of monogenic causes of CVID have been elucidated. From the standpoint of the practicing physician, an important question is whether the clinical phenotype, particularly the occurrence of autoinflammation of autoimmunity, might suggest the likelihood of identifying a causative mutation, and if possible the gene most likely to underlie CVID. We addressed this question in a patient group of 405 subjects diagnosed with CVID from one medical center.

Exploiting public databases of genomic variation to quantify evolutionary constraint on the branch point sequence in 30 plant and animal species

The branch point sequence is a degenerate intronic heptamer required for the assembly of the spliceosome during pre-mRNA splicing. Disruption of this motif may promote alternative splicing and eventually cause phenotype variation. Despite its functional relevance, the branch point sequence is not included in most genome annotations. Here, we predict branch point sequences in 30 plant and animal species and attempt to quantify their evolutionary constraints using public variant databases. We find an implausible variant distribution in the databases from 16 of 30 examined species. Comparative analysis of variants from whole-genome sequencing shows that variants submitted from exome sequencing or false positive variants are widespread in public databases and cause these irregularities. We then investigate evolutionary constraint with largely unbiased public variant databases in 14 species and find that the fourth and sixth position of the branch point sequence are more constrained than coding nucleotides. Our findings show that public variant databases should be scrutinized for possible biases before they qualify to analyze evolutionary constraint.

Accurate phenotypic classification and exome sequencing allow identification of novel genes and variants associated with adult-onset hearing loss

Adult-onset progressive hearing loss is a common, complex disease with a strong genetic component. Although to date over 150 genes have been identified as contributing to human hearing loss, many more remain to be discovered, as does most of the underlying genetic diversity. Many different variants have been found to underlie adult-onset hearing loss, but they tend to be rare variants with a high impact upon the gene product. It is likely that combinations of more common, lower impact variants also play a role in the prevalence of the disease. Here we present our exome study of hearing loss in a cohort of 532 older adult volunteers with extensive phenotypic data, including 99 older adults with normal hearing, an important control set. Firstly, we carried out an outlier analysis to identify genes with a high variant load in older adults with hearing loss compared to those with normal hearing. Secondly, we used audiometric threshold data to identify individual variants which appear to contribute to different threshold values. We followed up these analyses in a second cohort. Using these approaches, we identified genes and variants linked to better hearing as well as those linked to worse hearing. These analyses identified some known deafness genes, demonstrating proof of principle of our approach. However, most of the candidate genes are novel associations with hearing loss. While the results support the suggestion that genes responsible for severe deafness may also be involved in milder hearing loss, they also suggest that there are many more genes involved in hearing which remain to be identified. Our candidate gene lists may provide useful starting points for improved diagnosis and drug development.

Novel homozygous frameshift variant in the ATCAY gene in an Iranian patient with Cayman cerebellar ataxia; expanding the neuroimaging and clinical features: a case report

BACKGROUND

Pathogenic variants in the ATCAY gene are associated with a rare autosomal recessive disorder called Cayman cerebellar ataxia. Affected individuals display psychomotor retardation, cerebellar dysfunction, nystagmus, intention tremor, ataxic gait and dysarthric in some cases.

CASE PRESENTATION

Whole exome sequencing was performed for a 21-month-old girl suffering from developmental delay specifically in motor and language aspects, hypotonia, nystagmus, pes planus and strabismus. The detected variant in the patient was confirmed by family segregation analysis by Sanger sequencing in both of her parents. Previously three homozygous variants in the ATCAY gene (missense, splice site and frameshift deletion) have been reported in patients with Cayman cerebellar ataxia. Here we report the fourth homozygous variant and the second homozygous frameshift deletion in this gene to be associated with autosomal recessive Cayman cerebellar ataxia.

CONCLUSION

The identification of this novel homozygous frameshift deletion in the ATCAY gene expands our understanding of the genetic landscape underlying Cayman cerebellar ataxia. Furthermore, the occurrence of this variant in Iran, in addition to Pakistan, signifies the importance of considering genotypic and phenotypic factors beyond ethnicity when studying this disorder. These findings contribute to the ongoing efforts to unravel the molecular basis of Cayman cerebellar ataxia and improve diagnostic approaches and potential therapeutic interventions.

Accurate phenotypic classification and exome sequencing allow identification of novel genes and variants associated with adult-onset hearing loss

Adult-onset progressive hearing loss is a common, complex disease with a strong genetic component. Although to date over 150 genes have been identified as contributing to human hearing loss, many more remain to be discovered, as does most of the underlying genetic diversity. Many different variants have been found to underlie adult-onset hearing loss, but they tend to be rare variants with a high impact upon the gene product. It is likely that combinations of more common, lower impact variants also play a role in the prevalence of the disease. Here we present our exome study of hearing loss in a cohort of 532 older adult volunteers with extensive phenotypic data, including 99 older adults with normal hearing, an important control set. Firstly, we carried out an outlier analysis to identify genes with a high variant load in older adults with hearing loss compared to those with normal hearing. Secondly, we used audiometric threshold data to identify individual variants which appear to contribute to different threshold values. We followed up these analyses in a second cohort. Using these approaches, we identified genes and variants linked to better hearing as well as those linked to worse hearing. These analyses identified some known deafness genes, demonstrating proof of principle of our approach. However, most of the candidate genes are novel associations with hearing loss. While the results support the suggestion that genes responsible for severe deafness may also be involved in milder hearing loss, they also suggest that there are many more genes involved in hearing which remain to be identified. Our candidate gene lists may provide useful starting points for improved diagnosis and drug development.

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.

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.

Identifying shared genetic factors underlying epilepsy and congenital heart disease in Europeans

Epilepsy (EP) and congenital heart disease (CHD) are two apparently unrelated diseases that nevertheless display substantial mutual comorbidity. Thus, while 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 as 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. Finally, 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.

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.

Role of Moonlighting Proteins in Disease: Analyzing the Contribution of Canonical and Moonlighting Functions in Disease Progression

The term moonlighting proteins refers to those proteins that present alternative functions performed by a single polypeptide chain acquired throughout evolution (called canonical and moonlighting, respectively). Over 78% of moonlighting proteins are involved in human diseases, 48% are targeted by current drugs, and over 25% of them are involved in the virulence of pathogenic microorganisms. These facts encouraged us to study the link between the functions of moonlighting proteins and disease. We found a large number of moonlighting functions activated by pathological conditions that are highly involved in disease development and progression. The factors that activate some moonlighting functions take place only in pathological conditions, such as specific cellular translocations or changes in protein structure. Some moonlighting functions are involved in disease promotion while others are involved in curbing it. The disease-impairing moonlighting functions attempt to restore the homeostasis, or to reduce the damage linked to the imbalance caused by the disease. The disease-promoting moonlighting functions primarily involve the immune system, mesenchyme cross-talk, or excessive tissue proliferation. We often find moonlighting functions linked to the canonical function in a pathological context. Moonlighting functions are especially coordinated in inflammation and cancer. Wound healing and epithelial to mesenchymal transition are very representative. They involve multiple moonlighting proteins with a different role in each phase of the process, contributing to the current-phase phenotype or promoting a phase switch, mitigating the damage or intensifying the remodeling. All of this implies a new level of complexity in the study of pathology genesis, progression, and treatment. The specific protein function involved in a patient's progress or that is affected by a drug must be elucidated for the correct treatment of diseases.

Inherited human ITK deficiency impairs IFN-γ immunity and underlies tuberculosis

Inborn errors of IFN-γ immunity can underlie tuberculosis (TB). We report three patients from two kindreds without EBV viremia or disease but with severe TB and inherited complete ITK deficiency, a condition associated with severe EBV disease that renders immunological studies challenging. They have CD4+ αβ T lymphocytopenia with a concomitant expansion of CD4-CD8- double-negative (DN) αβ and Vδ2- γδ T lymphocytes, both displaying a unique CD38+CD45RA+T-bet+EOMES- phenotype. Itk-deficient mice recapitulated an expansion of the γδ T and DN αβ T lymphocyte populations in the thymus and spleen, respectively. Moreover, the patients' T lymphocytes secrete small amounts of IFN-γ in response to TCR crosslinking, mitogens, or forced synapse formation with autologous B lymphocytes. Finally, the patients' total lymphocytes secrete small amounts of IFN-γ, and CD4+, CD8+, DN αβ T, Vδ2+ γδ T, and MAIT cells display impaired IFN-γ production in response to BCG. Inherited ITK deficiency undermines the development and function of various IFN-γ-producing T cell subsets, thereby underlying TB.

Exome-wide association study to identify rare variants influencing COVID-19 outcomes: Results from the Host Genetics Initiative

Host genetics is a key determinant of COVID-19 outcomes. Previously, the COVID-19 Host Genetics Initiative genome-wide association study used common variants to identify multiple loci associated with COVID-19 outcomes. However, variants with the largest impact on COVID-19 outcomes are expected to be rare in the population. Hence, studying rare variants may provide additional insights into disease susceptibility and pathogenesis, thereby informing therapeutics development. Here, we combined whole-exome and whole-genome sequencing from 21 cohorts across 12 countries and performed rare variant exome-wide burden analyses for COVID-19 outcomes. In an analysis of 5,085 severe disease cases and 571,737 controls, we observed that carrying a rare deleterious variant in the SARS-CoV-2 sensor toll-like receptor TLR7 (on chromosome X) was associated with a 5.3-fold increase in severe disease (95% CI: 2.75-10.05, p = 5.41x10-7). This association was consistent across sexes. These results further support TLR7 as a genetic determinant of severe disease and suggest that larger studies on rare variants influencing COVID-19 outcomes could provide additional insights.

Computational approaches for predicting variant impact: An overview from resources, principles to applications

One objective of human genetics is to unveil the variants that contribute to human diseases. With the rapid development and wide use of next-generation sequencing (NGS), massive genomic sequence data have been created, making personal genetic information available. Conventional experimental evidence is critical in establishing the relationship between sequence variants and phenotype but with low efficiency. Due to the lack of comprehensive databases and resources which present clinical and experimental evidence on genotype-phenotype relationship, as well as accumulating variants found from NGS, different computational tools that can predict the impact of the variants on phenotype have been greatly developed to bridge the gap. In this review, we present a brief introduction and discussion about the computational approaches for variant impact prediction. Following an innovative manner, we mainly focus on approaches for non-synonymous variants (nsSNVs) impact prediction and categorize them into six classes. Their underlying rationale and constraints, together with the concerns and remedies raised from comparative studies are discussed. We also present how the predictive approaches employed in different research. Although diverse constraints exist, the computational predictive approaches are indispensable in exploring genotype-phenotype relationship.

DPP9 deficiency: An inflammasomopathy that can be rescued by lowering NLRP1/IL-1 signaling

Dipeptidyl peptidase 9 (DPP9) is a direct inhibitor of NLRP1, but how it affects inflammasome regulation in vivo is not yet established. Here, we report three families with immune-associated defects, poor growth, pancytopenia, and skin pigmentation abnormalities that segregate with biallelic DPP9 rare variants. Using patient-derived primary cells and biochemical assays, these variants were shown to behave as hypomorphic or knockout alleles that failed to repress NLRP1. The removal of a single copy of Nlrp1a/b/c, Asc, Gsdmd, or Il-1r, but not Il-18, was sufficient to rescue the lethality of Dpp9 mutant neonates in mice. Similarly, dpp9 deficiency was partially rescued by the inactivation of asc, an obligate downstream adapter of the NLRP1 inflammasome, in zebrafish. These experiments suggest that the deleterious consequences of DPP9 deficiency were mostly driven by the aberrant activation of the canonical NLRP1 inflammasome and IL-1β signaling. Collectively, our results delineate a Mendelian disorder of DPP9 deficiency driven by increased NLRP1 activity as demonstrated in patient cells and in two animal models of the disease.

Comparative Analyses of Turkish Variome and Widely Used Genomic Variation Databases for the Evaluation of Rare Sequence Variants in Turkish Individuals: Idiopathic Hypogonadotropic Hypogonadism as a Disease Model

Objective

With the increasing use of whole-exome sequencing, one of the challenges in identifying the causal allele for a Mendelian disease is the lack of availability of population-specific human genetic variation reference databases. The people of Turkey were not represented in GnomAD or other publicly available large databases until recently, when the first comprehensive genomic variation database, Turkish Variome (TRV), was published. The aim of this study was to evaluate whether TRV or other publicly available large genomic variation databases can reliably be used for rare disease variant evaluation in Turkish individuals.

Methods

Sixty non-disease-causing, non-synonymous variants (minor allele frequencies >1%) were identified in 58 genes that are known to be associated with idiopathic hypogonadotropic hypogonadism from a large Turkish patient cohort. The allelic frequencies of these variants were then compared with those in various public genomic variation databases, including TRV.

Results

Our cohort variants showed the highest correlations with those in the TRV, Iranome, and The Greater Middle East Variome, in decreasing order.

Conclusion

These results suggest that the TRV is the appropriate database to use for rare genomic variant evaluations in the Turkish population. Our data also suggest that variomes from geographic neighborhoods may serve as substitute references for populations devoid of their own genomic variation databases.

Aggregated Genomic Data as Cohort-Specific Allelic Frequencies can Boost Variants and Genes Prioritization in Non-Solved Cases of Inherited Retinal Dystrophies

The introduction of NGS in genetic diagnosis has increased the repertoire of variants and genes involved and the amount of genomic information produced. We built an allelic-frequency (AF) database for a heterogeneous cohort of genetic diseases to explore the aggregated genomic information and boost diagnosis in inherited retinal dystrophies (IRD). We retrospectively selected 5683 index-cases with clinical exome sequencing tests available, 1766 with IRD and the rest with diverse genetic diseases. We calculated a subcohort’s IRD-specific AF and compared it with suitable pseudocontrols. For non-solved IRD cases, we prioritized variants with a significant increment of frequencies, with eight variants that may help to explain the phenotype, and 10/11 of uncertain significance that were reclassified as probably pathogenic according to ACMG. Moreover, we developed a method to highlight genes with more frequent pathogenic variants in IRD cases than in pseudocontrols weighted by the increment of benign variants in the same comparison. We identified 18 genes for further studies that provided new insights in five cases. This resource can also help one to calculate the carrier frequency in IRD genes. A cohort-specific AF database assists with variants and genes prioritization and operates as an engine that provides a new hypothesis in non-solved cases, augmenting the diagnosis rate.

Multifocal organoids reveal clonal associations between synchronous intestinal tumors with pervasive heterogeneous drug responses

Multifocal colorectal cancer (CRC) comprises both clonally independent primary tumors caused by inherited predisposition and clonally related tumors mainly due to intraluminal spreading along an intact basement membrane. The distinction between these multifocal CRCs is essential because therapeutic strategies vary according to the clonal association of multiple tumor masses. Here, we report one unique case of synchronous intestinal cancer (SIC) with tumors occurring along the entire bowel tract, including the small intestine. We established six patient-derived organoids (PDOs), and patient-derived cell lines (PDCs) from each site of the SIC, which were subjected to extensive genomic, transcriptomic, and epigenomic sequencing. We also estimated the drug responses of each multifocal SIC to 25 clinically relevant therapeutic compounds to validate how the clinically actionable alternations between SICs were associated with drug sensitivity. Our data demonstrated distinct clonal associations across different organs, which were consistently supported by multi-omics analysis, as well as the accordant responses to various therapeutic compounds. Our results indicated the imminent drawback of a single tumor-based diagnosis of multifocal CRC and suggested the necessity of an in-depth molecular analysis of all tumor regions to avoid unexpected resistance to the currently available targeted therapies.

Premature termination codon: a tunable protein translation approach

Cellular protein-protein interactions are largely dependent on the activities of signaling proteins. Here, we present a technique to tune gene expression at translation level based on G418-inducible readthrough premature termination codon (PTC-on). To demonstrate how this PTC-on can control the expression level of a cellular signaling protein to regulate signal transduction, we settled a p53 PTC-on system in p53-null H1299 cells. After treating with G418, the cells expressed full-length p53 protein in a dose-dependent manner. We further demonstrated to use this PTC-on approach to dissect p53-dependent and p53-independent apoptosis in response to the DNA double strand breaks in H1299 cells. In principle, the PTC-on can be used as a general approach for exploring the functions of any other signaling proteins.

Investigating the characteristics of genes and variants associated with self-reported hearing difficulty in older adults in the UK Biobank

BACKGROUND

Age-related hearing loss is a common, heterogeneous disease with a strong genetic component. More than 100 loci have been reported to be involved in human hearing impairment to date, but most of the genes underlying human adult-onset hearing loss remain unknown. Most genetic studies have focussed on very rare variants (such as family studies and patient cohort screens) or very common variants (genome-wide association studies). However, the contribution of variants present in the human population at intermediate frequencies is hard to quantify using these methods, and as a result, the landscape of variation associated with adult-onset hearing loss remains largely unknown.

RESULTS

Here we present a study based on exome sequencing and self-reported hearing difficulty in the UK Biobank, a large-scale biomedical database. We have carried out variant load analyses using different minor allele frequency and impact filters, and compared the resulting gene lists to a manually curated list of nearly 700 genes known to be involved in hearing in humans and/or mice. An allele frequency cutoff of 0.1, combined with a high predicted variant impact, was found to be the most effective filter setting for our analysis. We also found that separating the participants by sex produced markedly different gene lists. The gene lists obtained were investigated using gene ontology annotation, functional prioritisation and expression analysis, and this identified good candidates for further study.

CONCLUSIONS

Our results suggest that relatively common as well as rare variants with a high predicted impact contribute to age-related hearing impairment and that the genetic contributions to adult hearing difficulty may differ between the sexes. Our manually curated list of deafness genes is a useful resource for candidate gene prioritisation in hearing loss.

A clinical, molecular genetics and pathological study of a FTDP-17 family with a heterozygous splicing variant c.823-10G>T at the intron 9/exon 10 of the MAPT gene

We report the first clinical-radiological-genetic-molecular-pathological study of a kindred with c.823-10G>T MAPT intronic variant (rs63749974) associated with frontotemporal dementia and parkinsonism linked to chromosome 17 (FTDP-17). We describe the clinical spectrum within this family and emphasize the association between MAPT gene variants and motor neuron disease. This report of a second family with FTDP-17 associated with c.823-10G>T MAPT variant strongly supports pathogenicity of the variant and confirms it is a 4-repeat (4R) tauopathy. This intronic point mutation, probably strengthens the polypyrimidine tract and alters the splicing of exon 10 (10 nucleotides into intron 9) close to the 3' splice site.

Inherited human c-Rel deficiency disrupts myeloid and lymphoid immunity to multiple infectious agents

We studied a child with severe viral, bacterial, fungal, and parasitic diseases, who was homozygous for a loss-of-function mutation of REL, encoding c-Rel, which is selectively expressed in lymphoid and myeloid cells. The patient had low frequencies of NK, effector memory cells reexpressing CD45RA (Temra) CD8+ T cells, memory CD4+ T cells, including Th1 and Th1*, Tregs, and memory B cells, whereas the counts and proportions of other leukocyte subsets were normal. Functional deficits of myeloid cells included the abolition of IL-12 and IL-23 production by conventional DC1s (cDC1s) and monocytes, but not cDC2s. c-Rel was also required for induction of CD86 expression on, and thus antigen-presenting cell function of, cDCs. Functional deficits of lymphoid cells included reduced IL-2 production by naive T cells, correlating with low proliferation and survival rates and poor production of Th1, Th2, and Th17 cytokines by memory CD4+ T cells. In naive CD4+ T cells, c-Rel is dispensable for early IL2 induction but contributes to later phases of IL2 expression. The patient's naive B cells displayed impaired MYC and BCL2L1 induction, compromising B cell survival and proliferation and preventing their differentiation into Ig-secreting plasmablasts. Inherited c-Rel deficiency disrupts the development and function of multiple myeloid and lymphoid cells, compromising innate and adaptive immunity to multiple infectious agents.

Inherited PD-1 deficiency underlies tuberculosis and autoimmunity in a child

The pathophysiology of adverse events following programmed cell death protein 1 (PD-1) blockade, including tuberculosis (TB) and autoimmunity, remains poorly characterized. We studied a patient with inherited PD-1 deficiency and TB who died of pulmonary autoimmunity. The patient's leukocytes did not express PD-1 or respond to PD-1-mediated suppression. The patient's lymphocytes produced only small amounts of interferon (IFN)-γ upon mycobacterial stimuli, similarly to patients with inborn errors of IFN-γ production who are vulnerable to TB. This phenotype resulted from a combined depletion of Vδ2⁺ γδ T, mucosal-associated invariant T and CD56^bright natural killer lymphocytes and dysfunction of other T lymphocyte subsets. Moreover, the patient displayed hepatosplenomegaly and an expansion of total, activated and RORγT⁺ CD4^-CD8^- double-negative αβ T cells, similar to patients with STAT3 gain-of-function mutations who display lymphoproliferative autoimmunity. This phenotype resulted from excessive amounts of STAT3-activating cytokines interleukin (IL)-6 and IL-23 produced by activated T lymphocytes and monocytes, and the STAT3-dependent expression of RORγT by activated T lymphocytes. Our work highlights the indispensable role of human PD-1 in governing both antimycobacterial immunity and self-tolerance, while identifying potentially actionable molecular targets for the diagnostic and therapeutic management of TB and autoimmunity in patients on PD-1 blockade.

Recurrent Pregnancy Loss and Concealed Long-QT Syndrome

Background

Recurrent pregnancy loss affects 1% to 2% of couples attempting childbirth. A large fraction of all cases remains idiopathic, which warrants research into monogenic causes of this distressing disorder.

Methods and Results

We investigated a nonconsanguineous Estonian family who had experienced 5 live births, intersected by 3 early pregnancy losses, and 6 fetal deaths, 3 of which occurred during the second trimester. No fetal malformations were described at the autopsies performed in 3 of 6 cases of fetal death. Parental and fetal chromosomal abnormalities (including submicroscopic) and maternal risk factors were excluded. Material for genetic testing was available from 4 miscarried cases (gestational weeks 11, 14, 17, and 18). Exome sequencing in 3 pregnancy losses and the mother identified no rare variants explicitly shared by the miscarried conceptuses. However, the mother and 2 pregnancy losses carried a heterozygous nonsynonymous variant, resulting in p.Val173Asp (rs199472695) in the ion channel gene KCNQ1. It is expressed not only in heart, where mutations cause type 1 long‐QT syndrome, but also in other tissues, including uterus. The p.Val173Asp variant has been previously identified in a patient with type 1 long‐QT syndrome, but not reported in the Genome Aggregation Database. With heterologous expression in CHO cells, our in vitro electrophysiologic studies indicated that the mutant slowly activating voltage‐gated K+ channel (I_Ks) is dysfunctional. It showed reduced total activating and deactivating currents (P<0.01), with dramatically positive shift of voltage dependence of activation by ≈10 mV (P<0.05).

Conclusions

The current study uncovered concealed maternal type 1 long‐QT syndrome as a potential novel cause behind recurrent fetal loss.

A computational approach for detecting physiological homogeneity in the midst of genetic heterogeneity

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.

Detection of homozygous and hemizygous complete or partial exon deletions by whole-exome sequencing

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

Interstitial Lung Disease in Common Variable Immunodeficiency

Interstitial lung disease (ILD) is a common complication in patients with common variable immunodeficiency (CVID) and often associated with other features, such as bronchiectasis and autoimmunity. As the ILD term encompasses different acute and chronic pulmonary conditions, the diagnosis is commonly made based on imaging features; histopathology is less frequently available. From a cohort of 637 patients with CVID followed at our center over 4 decades, we reviewed the data for 46 subjects (30 females, 16 males) who had lung biopsies with proven ILD. They had a median age at CVID diagnosis of 26 years old, with a median IgG level at diagnosis of 285.0 mg/dL with average isotype switched memory B cells of 0.5%. Lung biopsy pathology revealed granulomas in 25 patients (54.4%), lymphoid interstitial pneumonia in 13 patients (28.3%), lymphoid hyperplasia not otherwise specified in 7 patients (15.2%), cryptogenic organizing pneumonia in 7 patients (15.2%), follicular bronchitis in 4 patients (8.7%), and predominance of pulmonary fibrosis in 4 patients (8.7%). Autoimmune manifestations were common and were present in 28 (60.9%) patients. Nine patients (19.6%) died, with a median age at death of 49-years-old. Lung transplant was done in 3 of these patients (6.5%) who are no longer alive. These analyses reveal the high burden of this complication, with almost one-fifth of the group deceased in this period. Further understanding of the causes of the development and progression of ILD in CVID patients is required to define the best management for this patient population.

A call to eradicate non-inclusive terms from the life sciences

Since the Black Lives Matter movement rose to mainstream prominence, the academic enterprise has started recognizing the systematic racism present in science. However, there have been relatively few efforts to make sure that the language used to communicate science is inclusive. Here, I quantify the number of research articles published between 2000 and 2020 that contained non-inclusive terms with racial connotations, such as "blacklist" and "whitelist", or "master" and "slave". This reveals that non-inclusive language is being increasingly used in the life sciences literature, and I urge the global academic community to expunge these archaic terms to make science inclusive for everyone.

Spectrum-frequency and genotype-phenotype analysis of rhodopsin variants

Mutations in RHO are the most common cause of autosomal dominant retinitis pigmentosa. However, the pathogenicity of many RHO variants is questionable. This study was designed to investigate the genotype-phenotype correlation for RHO variants. These RHO variants were collected from the in-house exome sequencing data of 7092 probands suffering from different types of eye conditions. The variants were classified using bioinformatics tools, family segregation, and clinical phenotypes. The RHO variants were assessed using multiple online tools and a genotype-phenotype analysis based on the data collected from of ours, gnomAD, and published literature. Totally, 52 heterozygous variants of RHO were detected in the 7092 probands. Of these 52, 17 were potentially pathogenic, were present in 35 families, and comprised 15 missense variants, one inframe deletion and one nonsense variant. All the 15 missense variants were predicted to be damaging by five different online tools. The analysis of the clinical data of the patients from the 35 families revealed certain common features, of an early damage to both the rods and the cones, relatively preserved visual acuity in adulthood, and mid-peripheral tapetoretinal degeneration with pigmentation or RPE atrophy. Our data, the data from gnomAD, and the systematic review of the 246 previously reported variants suggest that approximately two-thirds of the rare missense variants and most of the truncated variants involving upstream of K296 are likely benign. This study provides a brief summary of the characteristics of the pathogenic RHO variants. It emphasizes that the systematic evaluation of these variants at the individual-gene level is crucial in the current era of clinical genetic testing even for a well-known gene such as RHO.

Current genetic landscape in common variable immune deficiency

Using whole-exome sequencing to examine the genetic causes of immune deficiency in 235 common variable immunodeficiency (CVID) patients seen in the United States (Mount Sinai, New York), 128 patients from Sweden, and 208 from Iran revealed 68 known disease-causing genes underlying this heterogeneous immune defect. The patients at the time of study ranged from 4 to 90 years of age. Overall, 31%, 36%, and 54% of the patients in the US, Swedish, or Iranian cohorts had mutations. The multiplicity of genes identified in the 571 subjects reflects the complex requirements of B-cell antigen signaling, activation, survival, migration, maturation, and maintenance of antibody-secreting memory B-cell populations to the plasma cell stage. For the US and Swedish cohorts, CVID subjects with noninfectious complications, lymphoid infiltrations, inflamatory conditions, or autoimmunity were somewhat more likely to have an identifiable gene, but in both cohorts, numerous subjects with these medical conditions had no potential gene that could be assigned. Specific clinical patterns of illnesses were also not linked to any given gene defect as there was considerable overlap in clinical presentations. These observations led to a new perspective on the complexity of the immunologic phenotype found in CVID syndrome.

Exome sequencing identification of susceptibility genes in Chinese patients with keratoconus

PURPOSE

Keratoconus (KC) is a corneal ectasia disease with complex genetic heterogeneity. The present study aimed to identify susceptibility genes in Chinese patients with KC.

METHODS

Exome sequencing (ES) was performed in 28 Chinese KC patients to search for susceptibility genes of the disease. The candidate variants were filtered out by multi-step bioinformatics analysis and validated by Sanger sequencing. Another 100 individuals with KC were also recruited to verify those variants by Sanger sequencing.

RESULTS

By filtering out nonsynonymous variants located in exon, selecting variants which were presented in two or more samples and applying public databases to remove common variants, along with the inclusion of missense SNVs located in differential expressed genes and protein damaging variants (stop gain/stop loss SNVs and InDels), we have identified 6 SNVs (4 missense SNVs: c.1168 T > C in TRANK1, c.341A>T in ERMP1, c.4346 T > C in SDK2, c.1730A>C in COL6A1; 2 stop gain SNVs: c.1138 C > T in CNBD1, c.241 C > T in KRT82) and 2 InDels (c.193_195del in NSUN5, c.1690_1698del in COL9A3) as candidate variants for KC. The verifying results showed that c.341A>T in ERMP1 and c.193_195del in NSUN5 was found in one and two samples, respectively.

CONCLUSIONS

Our study suggested that a total of six SNVs in six genes and two InDels in two genes might be considered as candidate variants in Chinese patients with KC.

Children's rare disease cohorts: an integrative research and clinical genomics initiative

While genomic data is frequently collected under distinct research protocols and disparate clinical and research regimes, there is a benefit in streamlining sequencing strategies to create harmonized databases, particularly in the area of pediatric rare disease. Research hospitals seeking to implement unified genomics workflows for research and clinical practice face numerous challenges, as they need to address the unique requirements and goals of the distinct environments and many stakeholders, including clinicians, researchers and sequencing providers. Here, we present outcomes of the first phase of the Children’s Rare Disease Cohorts initiative (CRDC) that was completed at Boston Children’s Hospital (BCH). We have developed a broadly sharable database of 2441 exomes from 15 pediatric rare disease cohorts, with major contributions from early onset epilepsy and early onset inflammatory bowel disease. All sequencing data is integrated and combined with phenotypic and research data in a genomics learning system (GLS). Phenotypes were both manually annotated and pulled automatically from patient medical records. Deployment of a genomically-ordered relational database allowed us to provide a modular and robust platform for centralized storage and analysis of research and clinical data, currently totaling 8516 exomes and 112 genomes. The GLS integrates analytical systems, including machine learning algorithms for automated variant classification and prioritization, as well as phenotype extraction via natural language processing (NLP) of clinical notes. This GLS is extensible to additional analytic systems and growing research and clinical collections of genomic and other types of data.

Non-canonical RNA-DNA differences and other human genomic features are enriched within very short tandem repeats

Very short tandem repeats bear substantial genetic, evolutional, and pathological significance in genome analyses. Here, we compiled a census of tandem mono-nucleotide/di-nucleotide/tri-nucleotide repeats (MNRs/DNRs/TNRs) in GRCh38, which we term "polytracts" in general. Of the human genome, 144.4 million nucleotides (4.7%) are occupied by polytracts, and 0.47 million single nucleotides are identified as polytract hinges, i.e., break-points of tandem polytracts. Preliminary exploration of the census suggested polytract hinge sites and boundaries of AAC polytracts may bear a higher mapping error rate than other polytract regions. Further, we revealed landscapes of polytract enrichment with respect to nearly a hundred genomic features. We found MNRs, DNRs, and TNRs displayed noticeable difference in terms of locational enrichment for miscellaneous genomic features, especially RNA editing events. Non-canonical and C-to-U RNA-editing events are enriched inside and/or adjacent to MNRs, while all categories of RNA-editing events are under-represented in DNRs. A-to-I RNA-editing events are generally under-represented in polytracts. The selective enrichment of non-canonical RNA-editing events within MNR adjacency provides a negative evidence against their authenticity. To enable similar locational enrichment analyses in relation to polytracts, we developed a software Polytrap which can handle 11 reference genomes. Additionally, we compiled polytracts of four model organisms into a Track Hub which can be integrated into USCS Genome Browser as an official track for convenient visualization of polytracts.

The MUC6/AP2A2 Locus and Its Relevance to Alzheimer's Disease: A Review

We recently reported evidence of Alzheimer's disease (AD)-linked genetic variation within the mucin 6 (MUC6) gene on chromosome 11p, nearby the adaptor-related protein complex 2 subunit alpha 2 (AP2A2) gene. This locus has interesting features related to human genomics and clinical research. MUC6 gene variants have been reported to potentially influence viral-including herpesvirus-immunity and the gut microbiome. Within the MUC6 gene is a unique variable number of tandem repeat (VNTR) region. We discovered an association between MUC6 VNTR repeat expansion and AD pathologic severity, particularly tau proteinopathy. Here, we review the relevant literature. The AD-linked VNTR polymorphism may also influence AP2A2 gene expression. AP2A2 encodes a polypeptide component of the adaptor protein complex, AP-2, which is involved in clathrin-coated vesicle function and was previously implicated in AD pathogenesis. To provide background information, we describe some key knowledge gaps in AD genetics research. The "missing/hidden heritability problem" of AD is highlighted. Extensive portions of the human genome, including the MUC6 VNTR, have not been thoroughly evaluated due to limitations of existing high-throughput sequencing technology. We present and discuss additional data, along with cautionary considerations, relevant to the hypothesis that MUC6 repeat expansion influences AD pathogenesis.

Identifying disease-causing mutations in genomes of single patients by computational approaches

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.

Problems in variation interpretation guidelines and in their implementation in computational tools

Background

ACMG/AMP and AMP/ASCO/CAP have released guidelines for variation interpretation, and ESHG for diagnostic sequencing. These guidelines contain recommendations including the use of computational prediction methods. The guidelines per se and the way they are implemented cause some problems.

Methods

Logical reasoning based on domain knowledge.

Results

According to the guidelines, several methods have to be used and they have to agree. This means that the methods with the poorest performance overrule the better ones. The choice of the prediction method(s) should be made by experts  based on systematic benchmarking studies reporting all the relevant performance measures. Currently variation interpretation methods have been applied mainly to amino acid substitutions and splice site variants; however, predictors for some other types of variations are available and there will be tools for new application areas in the near future. Common problems in prediction method usage are discussed. The number of features used for method training or the number of variation types predicted by a tool are not indicators of method performance. Many published gene, protein or disease‐specific benchmark studies suffer from too small dataset rendering the results useless. In the case of binary predictors, equal number of positive and negative cases is beneficial for training, the imbalance has to be corrected for performance assessment. Predictors cannot be better than the data they are based on and used for training and testing. Minor allele frequency (MAF) can help to detect likely benign cases, but the recommended MAF threshold is apparently too high. The fact that many rare variants are disease‐causing or ‐related does not mean that rare variants in general would be harmful. How large a portion of the tested variants a tool can predict (coverage) is not a quality measure.

Conclusion

Methods used for variation interpretation have to be carefully selected. It should be possible to use only one predictor, with proven good performance or a limited number of complementary predictors with state‐of‐the‐art performance. Bear in mind that diseases and pathogenicity have a continuum and variants are not dichotomic i.e. either pathogenic or benign, either.

De novo variants in exomes of congenital heart disease patients identify risk genes and pathways

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.

Homozygous NLRP1 gain-of-function mutation in siblings with a syndromic form of recurrent respiratory papillomatosis

Juvenile-onset recurrent respiratory papillomatosis (JRRP) is a rare and debilitating childhood disease that presents with recurrent growth of papillomas in the upper airway. Two common human papillomaviruses (HPVs), HPV-6 and -11, are implicated in most cases, but it is still not understood why only a small proportion of children develop JRRP following exposure to these common viruses. We report 2 siblings with a syndromic form of JRRP associated with mild dermatologic abnormalities. Whole-exome sequencing of the patients revealed a private homozygous mutation in NLRP1, encoding Nucleotide-Binding Domain Leucine-Rich Repeat Family Pyrin Domain-Containing 1. We find the NLRP1 mutant allele to be gain of function (GOF) for inflammasome activation, as demonstrated by the induction of inflammasome complex oligomerization and IL-1β secretion in an overexpression system. Moreover, patient-derived keratinocytes secrete elevated levels of IL-1β at baseline. Finally, both patients displayed elevated levels of inflammasome-induced cytokines in the serum. Six NLRP1 GOF mutations have previously been described to underlie 3 allelic Mendelian diseases with differing phenotypes and modes of inheritance. Our results demonstrate that an autosomal recessive, syndromic form of JRRP can be associated with an NLRP1 GOF mutation.

Which genes to assess in the NGS diagnostics of intellectual disability? The case for a consensus database-driven and expert-curated approach

When deciding on which genes to assess in larger Next-Generation Sequencing (NGS) datasets for the molecular genetic diagnosis of intellectual disability (ID), geneticists today have a variety of gene-phenotype databases and expert-curated gene lists available. To quantify their respective completeness, we compare an ID gene selection auto-generated from the Human Phenotype Ontology gene-phenotype association database and expert-curated ID gene lists from three reputable sources (sysID, the DDD consortium and Genomics England) and analyse some of their differences. We give examples of what we regard as genuine gaps ("missing ID genes") for each of these and conclude that a complementary or consensus approach is needed to maximise diagnostic yield in ID patients. We propose several consensus gene lists with ID-associated genes of different confidence levels.