DOMINO: Using Machine Learning to Predict Genes Associated with Dominant Disorders

A de novo ARIH2 gene mutation was detected in a patient with autism spectrum disorders and intellectual disability

E3 ubiquitin protein ligase encoded by ARIH2 gene catalyses the ubiquitination of target proteins and plays a crucial role in posttranslational modifications across various cellular processes. As prior documented, mutations in genes involved in the ubiquitination process are often associated with autism spectrum disorder (ASD) and/or intellectual disability (ID). In the current study, a de novo heterozygous mutation was identified in the splicing intronic region adjacent to the last exon of the ARIH2 gene using whole exome sequencing (WES). We hypothesize that this mutation, found in an ASD/ID patient, disrupts the protein Ariadne domain which is involved in the autoinhibition of ARIH2 enzyme. Predictive analyses elucidated the implications of the novel mutation in the splicing process and confirmed its autosomal dominant inheritance model. Nevertheless, we cannot exclude the possibility that other genetic factors, undetectable by WES, such as mutations in non-coding regions and polygenic risk in inter-allelic complementation, may contribute to the patient's phenotype. This work aims to suggest potential relationship between the detected mutation in ARIH2 gene and both ASD and ID, even though functional studies combined with new sequencing approaches will be necessary to validate this hypothesis.

Artificial Intelligence in Point-of-Care Biosensing: Challenges and Opportunities

The integration of artificial intelligence (AI) into point-of-care (POC) biosensing has the potential to revolutionize diagnostic methodologies by offering rapid, accurate, and accessible health assessment directly at the patient level. This review paper explores the transformative impact of AI technologies on POC biosensing, emphasizing recent computational advancements, ongoing challenges, and future prospects in the field. We provide an overview of core biosensing technologies and their use at the POC, highlighting ongoing issues and challenges that may be solved with AI. We follow with an overview of AI methodologies that can be applied to biosensing, including machine learning algorithms, neural networks, and data processing frameworks that facilitate real-time analytical decision-making. We explore the applications of AI at each stage of the biosensor development process, highlighting the diverse opportunities beyond simple data analysis procedures. We include a thorough analysis of outstanding challenges in the field of AI-assisted biosensing, focusing on the technical and ethical challenges regarding the widespread adoption of these technologies, such as data security, algorithmic bias, and regulatory compliance. Through this review, we aim to emphasize the role of AI in advancing POC biosensing and inform researchers, clinicians, and policymakers about the potential of these technologies in reshaping global healthcare landscapes.

Genetic variants for head size share genes and pathways with cancer

The size of the human head is highly heritable, but genetic drivers of its variation within the general population remain unmapped. We perform a genome-wide association study on head size (N = 80,890) and identify 67 genetic loci, of which 50 are novel. Neuroimaging studies show that 17 variants affect specific brain areas, but most have widespread effects. Gene set enrichment is observed for various cancers and the p53, Wnt, and ErbB signaling pathways. Genes harboring lead variants are enriched for macrocephaly syndrome genes (37-fold) and high-fidelity cancer genes (9-fold), which is not seen for human height variants. Head size variants are also near genes preferentially expressed in intermediate progenitor cells, neural cells linked to evolutionary brain expansion. Our results indicate that genes regulating early brain and cranial growth incline to neoplasia later in life, irrespective of height. This warrants investigation of clinical implications of the link between head size and cancer.

Herpes simplex encephalitis due to a mutation in an E3 ubiquitin ligase

Encephalitis is a rare and potentially fatal manifestation of herpes simplex type 1 infection. Following genome-wide genetic analyses, we identified a previously uncharacterized and very rare heterozygous variant in the E3 ubiquitin ligase WWP2, in a 14-month-old girl with herpes simplex encephalitis. The p.R841H variant (NM_007014.4:c.2522G > A) impaired TLR3 mediated signaling in inducible pluripotent stem cells-derived neural precursor cells and neurons; cells bearing this mutation were also more susceptible to HSV-1 infection compared to control cells. The p.R841H variant increased TRIF ubiquitination in vitro. Antiviral immunity was rescued following the correction of p.R841H by CRISPR-Cas9 technology. Moreover, the introduction of p.R841H in wild type cells reduced such immunity, suggesting that this mutation is linked to the observed phenotypes.

Biallelic loss-of-function variants of SLC12A9 cause lysosome dysfunction and a syndromic neurodevelopmental disorder

PURPOSE

Pathogenic variants of FIG4 generate enlarged lysosomes and neurological and developmental disorders. To identify additional genes regulating lysosomal volume, we carried out a genome-wide activation screen to detect suppression of enlarged lysosomes in FIG4-/- cells.

METHODS

The CRISPR-a gene activation screen utilized sgRNAs from the promoters of protein-coding genes. Fluorescence-activated cell sorting separated cells with correction of the enlarged lysosomes from uncorrected cells. Patient variants of SLC12A9 were identified by exome or genome sequencing and studied by segregation analysis and clinical characterization.

RESULTS

Overexpression of SLC12A9, a solute co-transporter, corrected lysosomal swelling in FIG4-/- cells. SLC12A9 (NP_064631.2) colocalized with LAMP2 at the lysosome membrane. Biallelic variants of SLC12A9 were identified in 3 unrelated probands with neurodevelopmental disorders. Common features included intellectual disability, skeletal and brain structural abnormalities, congenital heart defects, and hypopigmented hair. Patient 1 was homozygous for nonsense variant p.(Arg615∗), patient 2 was compound heterozygous for p.(Ser109Lysfs∗20) and a large deletion, and proband 3 was compound heterozygous for p.(Glu290Glyfs∗36) and p.(Asn552Lys). Fibroblasts from proband 1 contained enlarged lysosomes that were corrected by wild-type SLC12A9 cDNA. Patient variant p.(Asn552Lys) failed to correct the lysosomal defect.

CONCLUSION

Impaired function of SLC12A9 results in enlarged lysosomes and a recessive disorder with a recognizable neurodevelopmental phenotype.

Scrutinizing the human TEX genes in the context of human male infertility

BACKGROUND

Infertility affects around 15% of all couples worldwide and is increasingly linked to variants in genes specifically expressed in the testis. Well-established causes of male infertility include pathogenic variants in the genes TEX11, TEX14, and TEX15, while few studies have recently reported variants in TEX13B, TEX13C, FAM9A (TEX39A), and FAM9B (TEX39B).

OBJECTIVES

We aimed at screening for novel potential candidate genes among the human TEX ("testis expressed") genes as well as verifying previously described disease associations in this set of genes.

MATERIALS AND METHODS

To this end, we screened the exome sequencing data of 1305 men, including 1056 crypto- and azoospermic individuals, and determined cell-specific expression by analyzing testis-specific single-cell RNA sequencing data for genes with identified variants. To investigate the overarching role in male fertility, we generated testis-specific knockdown (KD) models of all 10 orthologous TEX genes in Drosophila melanogaster.

RESULTS

We detected rare potential disease-causing variants in TEX10, TEX13A, TEX13B, TEX13C, TEX13D, ZFAND3 (TEX27), TEX33, FAM9A (TEX39A), and FAM9B (TEX39B), in 28 infertile men, of which 15 men carried variants in TEX10, TEX27, and TEX33. The KD of TEX2, TEX9, TEX10, TEX13, ZFAND3 (TEX27), TEX28, TEX30, NFX1 (TEX42), TEX261, and UTP4 (TEX292) in Drosophila resulted in normal fertility.

DISCUSSION

Based on our findings, the autosomal dominant predicted genes TEX10 and ZFAND3 (TEX27) and the autosomal recessive predicted gene TEX33, which all three are conceivably required for germ cell maturation, were identified as novel potential candidate genes for human non-obstructive azoospermia. We additionally identified hemizygous loss-of-function (LoF) variants in TEX13B, TEX13C, and FAM9A (TEX39A) as unlikely monogenic culprits of male infertility as LoF variants were also found in control men.

CONCLUSION

Our findings concerning the X-linked genes TEX13B, TEX13C, and FAM9A (TEX39A) contradict previous reports and will decrease false-positive reports in genetic diagnostics of azoospermic men.

UNC5C: Novel Gene Associated with Psychiatric Disorders Impacts Dysregulation of Axon Guidance Pathways

The UNC-5 family of netrin receptor genes, predominantly expressed in brain tissues, plays a pivotal role in various neuronal processes. Mutations in genes involved in axon development contribute to a wide spectrum of human diseases, including developmental, neuropsychiatric, and neurodegenerative disorders. The NTN1/DCC signaling pathway, interacting with UNC5C, plays a crucial role in central nervous system axon guidance and has been associated with psychiatric disorders during adolescence in humans. Whole-exome sequencing analysis unveiled two compound heterozygous causative mutations within the UNC5C gene in a patient diagnosed with psychiatric disorders. In silico analysis demonstrated that neither of the observed variants affected the allosteric linkage between UNC5C and NTN1. In fact, these mutations are located within crucial cytoplasmic domains, specifically ZU5 and the region required for the netrin-mediated axon repulsion of neuronal growth cones. These domains play a critical role in forming the supramodular protein structure and directly interact with microtubules, thereby ensuring the functionality of the axon repulsion process. We emphasize that these mutations disrupt the aforementioned processes, thereby associating the UNC5C gene with psychiatric disorders for the first time and expanding the number of genes related to psychiatric disorders. Further research is required to validate the correlation of the UNC5C gene with psychiatric disorders, but we suggest including it in the genetic analysis of patients with psychiatric disorders.

HCN1 pathogenic variants associated with childhood epilepsy in a cohort of Chinese patients

OBJECTIVE

HCN ion channel family has a widespread expression in neurons, and recently, increasing studies have demonstrated their roles in epilepsies.

METHODS

Clinical data of the patients were gathered in a retrospective study. Exon sequencing was used for the patients with unexplained recurrent seizures and varying levels of developmental delay.

RESULTS

In this study, eight de novo variants of HCN1 genes were uncovered in eight patients, including six missense variants, one nonsense variant and one frameshift insertion variant; five of them were reported for the first time. The onset age for eight patients ranges from one month to one year. Their main clinical manifestations are epilepsy and varying degrees of developmental delay, and the main type of seizure is focal secondary generalized tonic-clonic seizure. Importantly, in our study, one case presented with a form of migrating focal seizure that has not been reported in the literature. Seizures from five of the eight children were effectively controlled with antiepileptic drugs including valproic acid, levetiracetam and oxcarbazepine. One child developed normally and four children developed mild delay. One child was treated with topiramate, and the convulsion was partially controlled and showed moderate to severe developmental delay. The antiepileptic treatment failed for the other two children, and the two children were treated with sodium valproate, oxcarbazepine, lamotrigine, chlorbazan, levetiracetam and nitrodiazepam successively, but their convulsions were not controlled and showed moderate to severe developmental delay.

SIGNIFICANCE

Our research reported eight variants in HCN1 gene causing epilepsy; among these variants, five variants were never reported before. HCN1-related epilepsy usually starts infantile period, and focal secondary generalized tonic-clonic seizure is the most common seizure type. Importantly, we reported the case with migrating focal seizure was rarely reported. Our study expanded both genotype and phenotype for HCN1-related epilepsy.

Next Generation Sequencing and Electromyography Reveal the Involvement of the P2RX6 Gene in Myopathy

Ion channelopathies result from impaired ion channel protein function, due to mutations affecting ion transport across cell membranes. Over 40 diseases, including neuropathy, pain, migraine, epilepsy, and ataxia, are associated with ion channelopathies, impacting electrically excitable tissues and significantly affecting skeletal muscle. Gene mutations affecting transmembrane ionic flow are strongly linked to skeletal muscle disorders, particularly myopathies, disrupting muscle excitability and contraction. Electromyography (EMG) analysis performed on a patient who complained of weakness and fatigue revealed the presence of primary muscular damage, suggesting an early-stage myopathy. Whole exome sequencing (WES) did not detect potentially causative variants in known myopathy-associated genes but revealed a novel homozygous deletion of the P2RX6 gene likely disrupting protein function. The P2RX6 gene, predominantly expressed in skeletal muscle, is an ATP-gated ion channel receptor belonging to the purinergic receptors (P2RX) family. In addition, STRING pathways suggested a correlation with more proteins having a plausible role in myopathy. No previous studies have reported the implication of this gene in myopathy. Further studies are needed on patients with a defective ion channel pathway, and the use of in vitro functional assays in suppressing P2RX6 gene expression will be required to validate its functional role.

Artificial intelligence and database for NGS-based diagnosis in rare disease

Rare diseases (RDs) are rare complex genetic diseases affecting a conservative estimate of 300 million people worldwide. Recent Next-Generation Sequencing (NGS) studies are unraveling the underlying genetic heterogeneity of this group of diseases. NGS-based methods used in RDs studies have improved the diagnosis and management of RDs. Concomitantly, a suite of bioinformatics tools has been developed to sort through big data generated by NGS to understand RDs better. However, there are concerns regarding the lack of consistency among different methods, primarily linked to factors such as the lack of uniformity in input and output formats, the absence of a standardized measure for predictive accuracy, and the regularity of updates to the annotation database. Today, artificial intelligence (AI), particularly deep learning, is widely used in a variety of biological contexts, changing the healthcare system. AI has demonstrated promising capabilities in boosting variant calling precision, refining variant prediction, and enhancing the user-friendliness of electronic health record (EHR) systems in NGS-based diagnostics. This paper reviews the state of the art of AI in NGS-based genetics, and its future directions and challenges. It also compare several rare disease databases.

De novo variants in PLCG1 are associated with hearing impairment, ocular pathology, and cardiac defects

Phospholipase C isozymes (PLCs) hydrolyze phosphatidylinositol 4,5-bisphosphate into inositol 1,4,5-trisphosphate and diacylglycerol, important signaling molecules involved in many cellular processes. PLCG1 encodes the PLCγ1 isozyme that is broadly expressed. Hyperactive somatic mutations of PLCG1 are observed in multiple cancers, but only one germline variant has been reported. Here we describe three unrelated individuals with de novo heterozygous missense variants in PLCG1 (p.Asp1019Gly, p.His380Arg, and p.Asp1165Gly) who exhibit variable phenotypes including hearing loss, ocular pathology and cardiac septal defects. To model these variants in vivo, we generated the analogous variants in the Drosophila ortholog, small wing (sl). We created a null allele slT2A and assessed the expression pattern. sl is broadly expressed, including in wing discs, eye discs, and a subset of neurons and glia. Loss of sl causes wing size reductions, ectopic wing veins and supernumerary photoreceptors. We document that mutant flies exhibit a reduced lifespan and age-dependent locomotor defects. Expressing wild-type sl in slT2A mutant rescues the loss-of-function phenotypes whereas expressing the variants causes lethality. Ubiquitous overexpression of the variants also reduces viability, suggesting that the variants are toxic. Ectopic expression of an established hyperactive PLCG1 variant (p.Asp1165His) in the wing pouch causes severe wing phenotypes, resembling those observed with overexpression of the p.Asp1019Gly or p.Asp1165Gly variants, further arguing that these two are gain-of-function variants. However, the wing phenotypes associated with p.His380Arg overexpression are mild. Our data suggest that the PLCG1 de novo heterozygous missense variants are pathogenic and contribute to the features observed in the probands.

GNB1-Related Rod-Cone Dystrophy: A Case Report

Introduction

The GNB1 (guanine nucleotide-binding protein, β1) gene encodes for the ubiquitous β1 subunit of heterotrimeric G proteins, which are associated with G-protein-coupled receptors (GPCRs). GNB1 mutations cause a neurodevelopmental disorder characterized by a broad clinical spectrum. A novel variant has recently been confirmed in a case of rod-cone dystrophy.

Case Presentation

We describe the second confirmed case of a classical rod-cone dystrophy associated with a mutation located in exon 6 of GNB1 [NM_002074.5:c.217G>C, p.(Ala73Pro)] in a 56-year-old patient also presenting mild intellectual disability, attention deficit/hyperactivity disorder, and truncal obesity.

Conclusion

This paper confirms the role of GNB1 in the pathogenesis of a classic rod-cone dystrophy and highlights the importance of including this gene in the genetic analysis panel for inherited retinal diseases.

Rare de novo gain-of-function missense variants in DOT1L are associated with developmental delay and congenital anomalies

Misregulation of histone lysine methylation is associated with several human cancers and with human developmental disorders. DOT1L is an evolutionarily conserved gene encoding a lysine methyltransferase (KMT) that methylates histone 3 lysine-79 (H3K79) and was not previously associated with a Mendelian disease in OMIM. We have identified nine unrelated individuals with seven different de novo heterozygous missense variants in DOT1L through the Undiagnosed Disease Network (UDN), the SickKids Complex Care genomics project, and GeneMatcher. All probands had some degree of global developmental delay/intellectual disability, and most had one or more major congenital anomalies. To assess the pathogenicity of the DOT1L variants, functional studies were performed in Drosophila and human cells. The fruit fly DOT1L ortholog, grappa, is expressed in most cells including neurons in the central nervous system. The identified DOT1L variants behave as gain-of-function alleles in flies and lead to increased H3K79 methylation levels in flies and human cells. Our results show that human DOT1L and fly grappa are required for proper development and that de novo heterozygous variants in DOT1L are associated with a Mendelian disease.

A de novo variant of BICRA results in Coffin-Siris syndrome 12

BACKGROUND

BICRA, a transcript regulator, was identified as the genetic factor of Coffin-Siris syndrome 12 (CSS12) recently, which was characterized by diverse neurodevelopmental delays. Up to now, limited studies of BICRA in neurodevelopmental delay have been reported.

METHODS

Clinical data such as EEGs, MRIs, routine blood, and physical examination were collected. Trio whole exome sequencing (WES) of the family was performed, and all variants with a minor allele frequency (<0.01) in exon and canonical splicing sites were selected for further pathogenic evaluation. Candidate variants were validated by Sanger sequencing. The BICRA-related literature was reviewed and the clinical characteristics were summarized.

RESULTS

We reported a CSS12 proband with a narrow and slightly clinical phenotype who only exhibited language developmental delay, hypotonia, and slight gastrointestinal features. WES revealed a de novo variant in exon 6 of BICRA [NM_015711.3: c.1666C>T, p.Gln556*]. This variant resulted in an early translation termination at 556th of BICRA, not collected in the public population database (gnomAD), and classified as pathogenic according to the ACMG guideline.

CONCLUSION

Our results expanded the pathogenic genetic and clinical spectrum of BICRA-related diseases.

De novo CLPTM1 variants with reduced GABAA R current response in patients with epilepsy

OBJECTIVE

To investigate the clinical features and potential pathogenesis mechanism of de novo CLPTM1 variants associated with epilepsy.

METHODS

Identify de novo genetic variants associated with epilepsy by reanalyzing trio-based whole-exome sequencing data. We analyzed the clinical characteristics of patients with these variants and performed functional in vitro studies in cells expressing mutant complementary DNA for these variants using whole-cell voltage-clamp current recordings and outside-out patch-clamp recordings from transiently transfected human embryonic kidney (HEK) cells.

RESULTS

Two de novo missense variants related to epilepsy were identified in the CLPTM1 gene. Functional studies indicated that CLPTM1-p.R454H and CLPTM1-p.R568Q variants reduced the γ-aminobutyric acid A receptor (GABAA R) current response amplitude recorded under voltage clamp compared to the wild-type receptors. These variants also reduced the charge transfer and altered the time course of desensitization and deactivation following rapid removal of GABA. The surface expression of the GABAA R γ2 subunit from the CLPTM1-p.R568Q group was significantly reduced compared to CLPTM1-WT.

SIGNIFICANCE

This is the first report of functionally relevant variants within the CLPTM1 gene. Patch-clamp recordings showed that these de novo CLPTM1 variants reduce GABAA R currents and charge transfer, which should promote excitation and hypersynchronous activity. This study may provide insights into the molecular mechanisms of the CLPTM1 variants underlying the patients' phenotypes, as well as for exploring potential therapeutic targets for epilepsy.

Whole-exome sequencing reveals candidate high-risk susceptibility genes for endometriosis

BACKGROUND

Endometriosis is a common, chronic disease among fertile-aged women. Disease course may be highly invasive, requiring extensive surgery. The etiology of endometriosis remains elusive, though a high level of heritability is well established. Several low-penetrance predisposing loci have been identified, but high-risk susceptibility remains undetermined. Endometriosis is known to increase the risk of epithelial ovarian cancers, especially of endometrioid and clear cell types. Here, we have analyzed a Finnish family where four women have been diagnosed with surgically verified, severely symptomatic endometriosis and two of the patients also with high-grade serous carcinoma.

RESULTS

Whole-exome sequencing revealed three rare candidate predisposing variants segregating with endometriosis. The variants were c.1238C>T, p.(Pro413Leu) in FGFR4, c.5065C>T, p.(Arg1689Trp) in NALCN, and c.2086G>A, p.(Val696Met) in NAV2. The only variant predicted deleterious by in silico tools was the one in FGFR4. Further screening of the variants in 92 Finnish endometriosis and in 19 endometriosis-ovarian cancer patients did not reveal additional carriers. Histopathology, positive p53 immunostaining, and genetic analysis supported the high-grade serous subtype of the two tumors in the family.

CONCLUSIONS

Here, we provide FGFR4, NALCN, and NAV2 as novel high-risk candidate genes for familial endometriosis. Our results also support the association of endometriosis with high-grade serous carcinoma. Further studies are required to validate the findings and to reveal the exact pathogenesis mechanisms of endometriosis. Elucidating the genetic background of endometriosis defines the etiology of the disease and provides opportunities for expedited diagnostics and personalized treatments.

XGDAG: explainable gene-disease associations via graph neural networks

MOTIVATION

Disease gene prioritization consists in identifying genes that are likely to be involved in the mechanisms of a given disease, providing a ranking of such genes. Recently, the research community has used computational methods to uncover unknown gene-disease associations; these methods range from combinatorial to machine learning-based approaches. In particular, during the last years, approaches based on deep learning have provided superior results compared to more traditional ones. Yet, the problem with these is their inherent black-box structure, which prevents interpretability.

RESULTS

We propose a new methodology for disease gene discovery, which leverages graph-structured data using graph neural networks (GNNs) along with an explainability phase for determining the ranking of candidate genes and understanding the model's output. Our approach is based on a positive-unlabeled learning strategy, which outperforms existing gene discovery methods by exploiting GNNs in a non-black-box fashion. Our methodology is effective even in scenarios where a large number of associated genes need to be retrieved, in which gene prioritization methods often tend to lose their reliability.

AVAILABILITY AND IMPLEMENTATION

The source code of XGDAG is available on GitHub at: https://github.com/GiDeCarlo/XGDAG. The data underlying this article are available at: https://www.disgenet.org/, https://thebiogrid.org/, https://doi.org/10.1371/journal.pcbi.1004120.s003, and https://doi.org/10.1371/journal.pcbi.1004120.s004.

A Dominant-Negative Mutant of ANXA7 Impairs Calcium Signaling and Enhances the Proliferation of Prostate Cancer Cells by Downregulating the IP3 Receptor and the PI3K/mTOR Pathway

Annexin A7/ANXA7 is a calcium-dependent membrane fusion protein with tumor suppressor gene (TSG) properties, which is located on chromosome 10q21 and is thought to function in the regulation of calcium homeostasis and tumorigenesis. However, whether the molecular mechanisms for tumor suppression are also involved in the calcium- and phospholipid-binding properties of ANXA7 remain to be elucidated. We hypothesized that the 4 C-terminal endonexin-fold repeats in ANXA7 (GX(X)GT), which are contained within each of the 4 annexin repeats with 70 amino acids, are responsible for both calcium- and GTP-dependent membrane fusion and the tumor suppressor function. Here, we identified a dominant-negative triple mutant (DNTM/DN-ANXA7J) that dramatically suppressed the ability of ANXA7 to fuse with artificial membranes while also inhibiting tumor cell proliferation and sensitizing cells to cell death. We also found that the [DNTM]ANA7 mutation altered the membrane fusion rate and the ability to bind calcium and phospholipids. In addition, in prostate cancer cells, our data revealed that variations in phosphatidylserine exposure, membrane permeabilization, and cellular apoptosis were associated with differential IP3 receptor expression and PI3K/AKT/mTOR modulation. In conclusion, we discovered a triple mutant of ANXA7, associated with calcium and phospholipid binding, which leads to the loss of several essential functions of ANXA7 pertinent to tumor protection and highlights the importance of the calcium signaling and membrane fusion functions of ANXA7 for preventing tumorigenesis.

Loss of Neuron Navigator 2 Impairs Brain and Cerebellar Development

Cerebellar hypoplasia and dysplasia encompass a group of clinically and genetically heterogeneous disorders frequently associated with neurodevelopmental impairment. The Neuron Navigator 2 (NAV2) gene (MIM: 607,026) encodes a member of the Neuron Navigator protein family, widely expressed within the central nervous system (CNS), and particularly abundant in the developing cerebellum. Evidence across different species supports a pivotal function of NAV2 in cytoskeletal dynamics and neurite outgrowth. Specifically, deficiency of Nav2 in mice leads to cerebellar hypoplasia with abnormal foliation due to impaired axonal outgrowth. However, little is known about the involvement of the NAV2 gene in human disease phenotypes. In this study, we identified a female affected with neurodevelopmental impairment and a complex brain and cardiac malformations in which clinical exome sequencing led to the identification of NAV2 biallelic truncating variants. Through protein expression analysis and cell migration assay in patient-derived fibroblasts, we provide evidence linking NAV2 deficiency to cellular migration deficits. In model organisms, the overall CNS histopathology of the Nav2 hypomorphic mouse revealed developmental anomalies including cerebellar hypoplasia and dysplasia, corpus callosum hypo-dysgenesis, and agenesis of the olfactory bulbs. Lastly, we show that the NAV2 ortholog in Drosophila, sickie (sick) is widely expressed in the fly brain, and sick mutants are mostly lethal with surviving escapers showing neurobehavioral phenotypes. In summary, our results unveil a novel human neurodevelopmental disorder due to genetic loss of NAV2, highlighting a critical conserved role of the NAV2 gene in brain and cerebellar development across species.

The first genetic landscape of inherited retinal dystrophies in Portuguese patients identifies recurrent homozygous mutations as a frequent cause of pathogenesis

Inherited retinal diseases (IRDs) are a group of ocular conditions characterized by an elevated genetic and clinical heterogeneity. They are transmitted almost invariantly as monogenic traits. However, with more than 280 disease genes identified so far, association of clinical phenotypes with genotypes can be very challenging, and molecular diagnosis is essential for genetic counseling and correct management of the disease. In addition, the prevalence and the assortment of IRD mutations are often population-specific. In this work, we examined 230 families from Portugal, with individuals suffering from a variety of IRD diagnostic classes (270 subjects in total). Overall, we identified 157 unique mutations (34 previously unreported) in 57 distinct genes, with a diagnostic rate of 76%. The IRD mutational landscape was, to some extent, different from those reported in other European populations, including Spanish cohorts. For instance, the EYS gene appeared to be the most frequently mutated, with a prevalence of 10% among all IRD cases. This was, in part, due to the presence of a recurrent and seemingly founder mutation involving the deletion of exons 13 and 14 of this gene. Moreover, our analysis highlighted that as many as 51% of our cases had mutations in a homozygous state. To our knowledge, this is the first study assessing a cross-sectional genotype-phenotype landscape of IRDs in Portugal. Our data reveal a rather unique distribution of mutations, possibly shaped by a small number of rare ancestral events that have now become prevalent alleles in patients.

Dendritic Spine in Autism Genetics: Whole-Exome Sequencing Identifying De Novo Variant of CTTNBP2 in a Quad Family Affected by Autism Spectrum Disorder

Autism spectrum disorder (ASD) affects around 1% of children with no effective blood test or cure. Recent studies have suggested that these are neurological disorders with a strong genetic basis and that they are associated with the abnormal formation of dendritic spines. Chromosome microarray (CMA) together with high-throughput sequencing technology has been used as a powerful tool to identify new candidate genes for ASD. In the present study, CMA was first used to scan for genome-wide copy number variants in a proband, and no clinically significant copy number variants were found. Whole-exome sequencing (WES) was used further for genetic testing of the whole quad family affected by ASD, including the proband, his non-autistic sister, and his parents. Sanger sequencing and MassARRAY-based validation were used to identify and confirm variants associated with ASD. WES yielded a 151-fold coverage depth for each sample. A total of 98.65% of the targeted whole-exome region was covered at >20-fold depth. A de novo variant in CTTNBP2, p.M115T, was identified. The CTTNBP2 gene belongs to a family of ankyrin repeat domain-containing proteins associated with dendritic spine formation. Although CTTNBP2 has been associated with ASD, limited studies have been developed to identify clinically relevant de novo mutations of CTTNBP2 in children with ASD; family-based WES successfully identified a clinically relevant mutation in the CTTNBP2 gene in a quad family affected by ASD. Considering the neuron-specific expression of CTTNBP2 and its role in dendritic spine formation, our results suggest a correlation between the CTTNBP2 mutation and ASD, providing genetic evidence for ASD spine pathology. Although the present study is currently insufficient to support the assertion that the de novo mutation M115T in CTTNBP2 directly causes the autism phenotype, our study provides support for the assertion that this mutation is a candidate clinically relevant variant in autism.

Predicting Genetic Disorder and Types of Disorder Using Chain Classifier Approach

Genetic disorders are the result of mutation in the deoxyribonucleic acid (DNA) sequence which can be developed or inherited from parents. Such mutations may lead to fatal diseases such as Alzheimer's, cancer, Hemochromatosis, etc. Recently, the use of artificial intelligence-based methods has shown superb success in the prediction and prognosis of different diseases. The potential of such methods can be utilized to predict genetic disorders at an early stage using the genome data for timely treatment. This study focuses on the multi-label multi-class problem and makes two major contributions to genetic disorder prediction. A novel feature engineering approach is proposed where the class probabilities from an extra tree (ET) and random forest (RF) are joined to make a feature set for model training. Secondly, the study utilizes the classifier chain approach where multiple classifiers are joined in a chain and the predictions from all the preceding classifiers are used by the conceding classifiers to make the final prediction. Because of the multi-label multi-class data, macro accuracy, Hamming loss, and α-evaluation score are used to evaluate the performance. Results suggest that extreme gradient boosting (XGB) produces the best scores with a 92% α-evaluation score and a 84% macro accuracy score. The performance of XGB is much better than state-of-the-art approaches, in terms of both performance and computational complexity.

Analysis of copy number variation in men with non-obstructive azoospermia

BACKGROUND

Recent findings demonstrate that single nucleotide variants can cause non-obstructive azoospermia (NOA). In contrast, copy number variants (CNVs) were only analysed in few studies in infertile men. Some have reported a higher prevalence of CNVs in infertile versus fertile men.

OBJECTIVES

This study aimed to elucidate if CNVs are associated with NOA.

MATERIALS AND METHODS

We performed array-based comparative genomic hybridisation (aCGH) in 37 men with meiotic arrest, 194 men with Sertoli cell-only phenotype, and 21 control men. We filtered our data for deletions affecting genes and prioritised the affected genes according to the literature search. Prevalence of CNVs was compared between all groups. Exome data of 2,030 men were screened to detect further genetic variants in prioritised genes. Modelling was performed for the protein encoded by the novel candidate gene TEKT5 and we stained for TEKT5 in human testicular tissue.

RESULTS

We determined the cause of infertility in two individuals with homozygous deletions of SYCE1 and in one individual with a heterozygous deletion of SYCE1 combined with a likely pathogenic missense variant on the second allele. We detected heterozygous deletions affecting MLH3, EIF2B2, SLX4, CLPP and TEKT5, in one subject each. CNVs were not detected more frequently in infertile men compared with controls.

DISCUSSION

While SYCE1 and MLH3 encode known meiosis-specific proteins, much less is known about the proteins encoded by the other identified candidate genes, warranting further analyses. We were able to identify the cause of infertility in one out of the 231 infertile men by aCGH and in two men by using exome sequencing data.

CONCLUSION

As aCGH and exome sequencing are both expensive methods, combining both in a clinical routine is not an effective strategy. Instead, using CNV calling from exome data has recently become more precise, potentially making aCGH dispensable.

Mendelian gene identification through mouse embryo viability screening

BACKGROUND

The diagnostic rate of Mendelian disorders in sequencing studies continues to increase, along with the pace of novel disease gene discovery. However, variant interpretation in novel genes not currently associated with disease is particularly challenging and strategies combining gene functional evidence with approaches that evaluate the phenotypic similarities between patients and model organisms have proven successful. A full spectrum of intolerance to loss-of-function variation has been previously described, providing evidence that gene essentiality should not be considered as a simple and fixed binary property.

METHODS

Here we further dissected this spectrum by assessing the embryonic stage at which homozygous loss-of-function results in lethality in mice from the International Mouse Phenotyping Consortium, classifying the set of lethal genes into one of three windows of lethality: early, mid, or late gestation lethal. We studied the correlation between these windows of lethality and various gene features including expression across development, paralogy and constraint metrics together with human disease phenotypes. We explored a gene similarity approach for novel gene discovery and investigated unsolved cases from the 100,000 Genomes Project.

RESULTS

We found that genes in the early gestation lethal category have distinct characteristics and are enriched for genes linked with recessive forms of inherited metabolic disease. We identified several genes sharing multiple features with known biallelic forms of inborn errors of the metabolism and found signs of enrichment of biallelic predicted pathogenic variants among early gestation lethal genes in patients recruited under this disease category. We highlight two novel gene candidates with phenotypic overlap between the patients and the mouse knockouts.

CONCLUSIONS

Information on the developmental period at which embryonic lethality occurs in the knockout mouse may be used for novel disease gene discovery that helps to prioritise variants in unsolved rare disease cases.

Genomic testing for copy number and single nucleotide variants in spermatogenic failure

PURPOSE

To identify clinically significant genomic copy number (CNV) and single nucleotide variants (SNV) in males with unexplained spermatogenic failure (SPGF).

MATERIALS AND METHODS

Peripheral blood DNA from 97/102 study participants diagnosed with oligozoospermia, severe oligozoospermia, or non-obstructive azoospermia (NOA) was analyzed for CNVs via array comparative genomic hybridization (aCGH) and SNVs using whole-exome sequencing (WES).

RESULTS

Of the 2544 CNVs identified in individuals with SPGF, > 90% were small, ranging from 0.6 to 75 kb. Thirty, clinically relevant genomic aberrations, were detected in 28 patients (~ 29%). These included likely diagnostic CNVs in 3/41 NOA patients (~ 7%): 1 hemizygous, intragenic TEX11 deletion, 1 hemizygous DDX53 full gene deletion, and 1 homozygous, intragenic STK11 deletion. High-level mosaicism for X chromosome disomy (~ 10% 46,XY and ~ 90% 47,XXY) was also identified in 3 of 41 NOA patients who previously tested normal with conventional karyotyping. The remaining 24 CNVs detected were heterozygous, autosomal recessive carrier variants. Follow-up WES analysis confirmed 8 of 27 (30%) CNVs (X chromosome disomy excluded). WES analysis additionally identified 13 significant SNVs and/or indels in 9 patients (~ 9%) including X-linked AR, KAL1, and NR0B1 variants.

CONCLUSION

Using a combined genome-wide aCGH/WES approach, we identified pathogenic and likely pathogenic SNVs and CNVs in 15 patients (15%) with unexplained SPGF. This value equals the detection rate of conventional testing for aneuploidies and is considerably higher than the prevalence of Y chromosome microdeletions. Our results underscore the importance of comprehensive genomic analysis in emerging diagnostic testing of complex conditions like male infertility.

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.

Biallelic variants in ZNF526 cause a severe neurodevelopmental disorder with microcephaly, bilateral cataract, epilepsy and simplified gyration

BACKGROUND

Next-generation sequencing, combined with international pooling of cases, has impressively enhanced the discovery of genes responsible for Mendelian neurodevelopmental disorders, particularly in individuals affected by clinically undiagnosed diseases. To date, biallelic missense variants in ZNF526 gene, encoding a Krüppel-type zinc-finger protein, have been reported in three families with non-syndromic intellectual disability.

METHODS

Here, we describe five individuals from four unrelated families with an undiagnosed neurodevelopmental disorder in which we performed exome sequencing, on a combination of trio-based (4 subjects) or single probands (1 subject).

RESULTS

We identified five patients from four unrelated families with homozygous ZNF526 variants by whole exome sequencing. Four had variants resulting in truncation of ZNF526; they were affected by severe prenatal and postnatal microcephaly (ranging from -4 SD to -8 SD), profound psychomotor delay, hypertonic-dystonic movements, epilepsy and simplified gyral pattern on MRI. All of them also displayed bilateral progressive cataracts. A fifth patient had a homozygous missense variant and a slightly less severe disorder, with postnatal microcephaly (-2 SD), progressive bilateral cataracts, severe intellectual disability and unremarkable brain MRI.Mutant znf526 zebrafish larvae had notable malformations of the eye and central nervous system, resembling findings seen in the human holoprosencephaly spectrum.

CONCLUSION

Our findings support the role of ZNF526 biallelic variants in a complex neurodevelopmental disorder, primarily affecting brain and eyes, resulting in severe microcephaly, simplified gyral pattern, epileptic encephalopathy and bilateral cataracts.

Exome sequencing in children with clinically suspected maturity-onset diabetes of the young

OBJECTIVE

Commercial gene panels identify pathogenic variants in as low as 27% of patients suspected to have MODY, suggesting the role of yet unidentified pathogenic variants. We sought to identify novel gene variants associated with MODY.

RESEARCH DESIGN AND METHODS

We recruited 10 children with a clinical suspicion of MODY but non-diagnostic commercial MODY gene panels. We performed exome sequencing (ES) in them and their parents.

RESULTS

Mean age at diabetes diagnosis was 10 (± 3.8) years. Six were females; 4 were non-Hispanic white, 5 Hispanic, and 1 Asian. Our variant prioritization analysis identified a pathogenic, de novo variant in INS (c.94G > A, p.Gly32Ser), confirmed by Sanger sequencing, in a proband who was previously diagnosed with "autoantibody-negative type 1 diabetes (T1D)" at 3 y/o. This rare variant, absent in the general population (gnomAD database), has been reported previously in neonatal diabetes. We also identified a frameshift deletion (c.2650delC, p.Gln884AsnfsTer57) in RFX6 in a child with a previous diagnosis of "autoantibody-negative T1D" at 12 y/o. The variant was inherited from the mother, who was diagnosed with "thin type 2 diabetes" at 25 y/o. Heterozygous protein-truncating variants in RFX6 gene have been recently reported in individuals with MODY.

CONCLUSIONS

We diagnosed two patients with MODY using ES in children initially classified as "T1D". One has a likely pathogenic novel gene variant not previously associated with MODY. We demonstrate the clinical utility of ES in patients with clinical suspicion of MODY.

A domain damage index to prioritizing the pathogenicity of missense variants

Prioritizing causal variants is one major challenge for the clinical application of sequencing data. Prompted by the observation that 74.3% of missense pathogenic variants locate in protein domains, we developed an approach named domain damage index (DDI). DDI identifies protein domains depleted of rare missense variations in the general population, which can be further used as a metric to prioritize variants. DDI is significantly correlated with phylogenetic conservation, variant-level metrics, and reported pathogenicity. DDI achieved great performance for distinguishing pathogenic variants from benign ones in three benchmark datasets. The combination of DDI with the other two best approaches improved the performance of each individual method considerably, suggesting DDI provides a powerful and complementary way of variant prioritization.

De novo variants in CACNA1E found in patients with intellectual disability, developmental regression and social cognition deficit but no seizures

BACKGROUND

De novo variants in the voltage-gated calcium channel subunit α1 E gene (CACNA1E) have been described as causative of epileptic encephalopathy with contractures, macrocephaly and dyskinesias.

METHODS

Following the observation of an index patient with developmental delay and autism spectrum disorder (ASD) without seizures who had a de novo deleterious CACNA1E variant, we screened GeneMatcher for other individuals with CACNA1E variants and neurodevelopmental phenotypes without epilepsy. The spectrum of pathogenic CACNA1E variants was compared to the mutational landscape of variants in the gnomAD control population database.

RESULTS

We identified seven unrelated individuals with intellectual disability, developmental regression and ASD-like behavioral profile, and notably without epilepsy, who had de novo heterozygous putatively pathogenic variants in CACNA1E. Age of onset of clinical manifestation, presence or absence of regression and degree of severity were variable, and no clear-cut genotype-phenotype association could be recognized. The analysis of disease-associated variants and their comparison to benign variants from the control population allowed for the identification of regions in the CACNA1E protein that seem to be intolerant to substitutions and thus more likely to harbor pathogenic variants. As in a few reported cases with CACNA1E variants and epilepsy, one patient showed a positive clinical behavioral response to topiramate, a specific calcium channel modulator.

LIMITATIONS

The significance of our study is limited by the absence of functional experiments of the effect of identified variants, the small sample size and the lack of systematic ASD assessment in all participants. Moreover, topiramate was given to one patient only and for a short period of time.

CONCLUSIONS

Our results indicate that CACNA1E variants may result in neurodevelopmental disorders without epilepsy and expand the mutational and phenotypic spectrum of this gene. CACNA1E deserves to be included in gene panels for non-specific developmental disorders, including ASD, and not limited to patients with seizures, to improve diagnostic recognition and explore the possible efficacy of topiramate.

Bi-allelic PAGR1 variants are associated with microcephaly and a severe neurodevelopmental disorder: Genetic evidence from two families

Exome and genome sequencing were used to identify the genetic etiology of a severe neurodevelopmental disorder in two unrelated Ashkenazi Jewish families with three affected individuals. The clinical findings included a prenatal presentation of microcephaly, polyhydramnios and clenched hands while postnatal findings included microcephaly, severe developmental delay, dysmorphism, neurologic deficits, and death in infancy. A shared rare homozygous, missense variant (c.274A > G; p.Ser92Gly, NM_024516.4) was identified in PAGR1, a gene currently not associated with a Mendelian disease. PAGR1 encodes a component of the histone methyltransferase MLL2/MLL3 complex and may function in the DNA damage response pathway. Complete knockout of the murine Pagr1a is embryonic-lethal. Given the available evidence, PAGR1 is a strong candidate gene for a novel autosomal recessive severe syndromic neurodevelopmental disorder.

Yield of clinically reportable genetic variants in unselected cerebral palsy by whole genome sequencing

Cerebral palsy (CP) is the most common cause of childhood physical disability, with incidence between 1/500 and 1/700 births in the developed world. Despite increasing evidence for a major contribution of genetics to CP aetiology, genetic testing is currently not performed systematically. We assessed the diagnostic rate of genome sequencing (GS) in a clinically unselected cohort of 150 singleton CP patients, with CP confirmed at >4 years of age. Clinical grade GS was performed on the proband and variants were filtered, and classified according to American College of Medical Genetics and Genomics–Association for Molecular Pathology (ACMG-AMP) guidelines. Variants classified as pathogenic or likely pathogenic (P/LP) were further assessed for their contribution to CP. In total, 24.7% of individuals carried a P/LP variant(s) causing or increasing risk of CP, with 4.7% resolved by copy number variant analysis and 20% carrying single nucleotide or indel variants. A further 34.7% carried one or more rare, high impact variants of uncertain significance (VUS) in variation intolerant genes. Variants were identified in a heterogeneous group of genes, including genes associated with hereditary spastic paraplegia, clotting and thrombophilic disorders, small vessel disease, and other neurodevelopmental disorders. Approximately 1/2 of individuals were classified as likely to benefit from changed clinical management as a result of genetic findings. In addition, no significant association between genetic findings and clinical factors was detectable in this cohort, suggesting that systematic sequencing of CP will be required to avoid missed diagnoses.

CNV Detection from Exome Sequencing Data in Routine Diagnostics of Rare Genetic Disorders: Opportunities and Limitations

To assess the potential of detecting copy number variations (CNVs) directly from exome sequencing (ES) data in diagnostic settings, we developed a CNV-detection pipeline based on ExomeDepth software and applied it to ES data of 450 individuals. Initially, only CNVs affecting genes in the requested diagnostic gene panels were scored and tested against arrayCGH results. Pathogenic CNVs were detected in 18 individuals. Most detected CNVs were larger than 400 kb (11/18), but three individuals had small CNVs impacting one or a few exons only and were thus not detectable by arrayCGH. Conversely, two pathogenic CNVs were initially missed, as they impacted genes not included in the original gene panel analysed, and a third one was missed as it was in a poorly covered region. The overall combined diagnostic rate (SNVs + CNVs) in our cohort was 36%, with wide differences between clinical domains. We conclude that (1) the ES-based CNV pipeline detects efficiently large and small pathogenic CNVs, (2) the detection of CNV relies on uniformity of sequencing and good coverage, and (3) in patients who remain unsolved by the gene panel analysis, CNV analysis should be extended to all captured genes, as diagnostically relevant CNVs may occur everywhere in the genome.

Genotype-phenotype correlation in Phelan-McDermid syndrome: A comprehensive review of chromosome 22q13 deleted genes

Phelan‐McDermid syndrome (PMS, OMIM #606232), also known as chromosome 22q13 deletion syndrome, is a rare genetic disorder characterized by intellectual disability, hypotonia, delayed or absent speech, motor impairment, autism spectrum disorder, behavioral anomalies, and minor aspecific dysmorphic features. Haploinsufficiency of SHANK3, due to intragenic deletions or point mutations, is sufficient to cause many neurobehavioral features of PMS. However, several additional genes located within larger 22q13 deletions can contribute to the great interindividual variability observed in the PMS phenotype. This review summarizes the phenotypic contributions predicted for 213 genes distributed along the largest 22q13.2‐q13.33 terminal deletion detected in our sample of 63 PMS patients by array‐CGH analysis, spanning 9.08 Mb. Genes have been grouped into four categories: (1) genes causing human diseases with an autosomal dominant mechanism, or (2) with an autosomal recessive mechanism; (3) morphogenetically relevant genes, either involved in human diseases with additive co‐dominant, polygenic, and/or multifactorial mechanisms, or implicated in animal models but not yet documented in human pathology; (4) protein coding genes either functionally nonrelevant, with unknown function, or pathogenic through mechanisms other than haploinsufficiency; piRNAs, noncoding RNAs, miRNAs, novel transcripts and pseudogenes. Our aim is to understand genotype–phenotype correlations in PMS patients and to provide clinicians with a conceptual framework to promote evidence‐based genetic work‐ups, clinical assessments, and therapeutic interventions.

Insights into the pathogenicity of missense variants in the forkhead domain of FOX proteins underlying Mendelian disorders

Forkhead box (FOX) proteins are members of a conserved family of transcription factors. Pathogenic variants in FOX genes have been shown to be responsible for several human genetic diseases. Here, we have studied the molecular and structural features of germline pathogenic variants in seven FOX proteins involved in Mendelian disorders and compared them with those of variants present in the general population (gnomAD). Our study shows that the DNA-binding domain of FOX proteins is particularly sensitive to damaging variation, although some family members show greater mutational tolerance than others. Next, we set to demonstrate that this tolerance depends on the inheritance mode of FOX-linked disorders. Accordingly, genes whose variants underlie recessive conditions are supposed to have a greater tolerance to variation. This is what we found. As expected, variants responsible for disorders with a dominant inheritance pattern show a higher degree of pathogenicity compared to those segregating in the general population. Moreover, we show that pathogenic and likely pathogenic variants tend to affect mutually exclusive sites with respect to those reported in gnomAD. The former also tend to affect sites with lower solvent exposure and a higher degree of conservation. Our results show the value of using publicly available databases and bioinformatics to gain insights into the molecular and structural bases of disease-causing genetic variation.

Developmental delay with hypotrophy associated with homozygous functionally relevant REV3L variant

REV3L encodes a catalytic subunit of DNA polymerase zeta (Pol zeta) which is essential for the tolerance of DNA damage by inducing translesion synthesis (TLS). So far, the only Mendelian disease associated with REV3L was Moebius syndrome (3 patients with dominant REV3L mutations causing monoallelic loss-of-function were reported). We describe a homozygous ultra-rare REV3L variant (T2753R) identified with whole exome sequencing in a child without Moebius syndrome but with developmental delay, hypotrophy, and dysmorphic features who was born to healthy parents (heterozygous carriers of the variant). The variant affects the amino acid adjacent to functionally important KKRY motif. By introducing an equivalent mutation (S1192R) into the REV3 gene in yeasts, we showed that, whereas it retained residual function, it caused clear dysfunction of TLS in the nucleus and instability of mitochondrial genetic information. In particular, the mutation increased UV sensitivity measured by cell survival, decreased both the spontaneous (P < 0.005) and UV-induced (P < 0.0001) mutagenesis rates of nuclear DNA and increased the UV-induced mutagenesis rates of mitochondrial DNA (P < 0.0005). We propose that our proband is the first reported case of a REV3L associated disease different from Moebius syndrome both in terms of clinical manifestations and inheritance (autosomal recessive rather than dominant). KEY MESSAGES: First description of a human recessive disorder associated with a REV3L variant. A study in yeast showed that the variant affected the enzymatic function of the protein. In particular, it caused increased UV sensitivity and abnormal mutagenesis rates.

MVP predicts the pathogenicity of missense variants by deep learning

Accurate pathogenicity prediction of missense variants is critically important in genetic studies and clinical diagnosis. Previously published prediction methods have facilitated the interpretation of missense variants but have limited performance. Here, we describe MVP (Missense Variant Pathogenicity prediction), a new prediction method that uses deep residual network to leverage large training data sets and many correlated predictors. We train the model separately in genes that are intolerant of loss of function variants and the ones that are tolerant in order to take account of potentially different genetic effect size and mode of action. We compile cancer mutation hotspots and de novo variants from developmental disorders for benchmarking. Overall, MVP achieves better performance in prioritizing pathogenic missense variants than previous methods, especially in genes tolerant of loss of function variants. Finally, using MVP, we estimate that de novo coding variants contribute to 7.8% of isolated congenital heart disease, nearly doubling previous estimates.

Accurate prediction of variant pathogenicity is essential to understanding genetic risks in disease. Here, the authors present a deep neural network method for prediction of missense variant pathogenicity, MVP, and demonstrate its utility in prioritizing de novo variants contributing to developmental disorders.

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

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

Identification of haploinsufficient genes from epigenomic data using deep forest

Haploinsufficiency, wherein a single allele is not enough to maintain normal functions, can lead to many diseases including cancers and neurodevelopmental disorders. Recently, computational methods for identifying haploinsufficiency have been developed. However, most of those computational methods suffer from study bias, experimental noise and instability, resulting in unsatisfactory identification of haploinsufficient genes. To address those challenges, we propose a deep forest model, called HaForest, to identify haploinsufficient genes. The multiscale scanning is proposed to extract local contextual representations from input features under Linear Discriminant Analysis. After that, the cascade forest structure is applied to obtain the concatenated features directly by integrating decision-tree-based forests. Meanwhile, to exploit the complex dependency structure among haploinsufficient genes, the LightGBM library is embedded into HaForest to reveal the highly expressive features. To validate the effectiveness of our method, we compared it to several computational methods and four deep learning algorithms on five epigenomic data sets. The results reveal that HaForest achieves superior performance over the other algorithms, demonstrating its unique and complementary performance in identifying haploinsufficient genes. The standalone tool is available at https://github.com/yangyn533/HaForest.

Transcriptome-directed analysis for Mendelian disease diagnosis overcomes limitations of conventional genomic testing

BACKGROUNDTranscriptome sequencing (RNA-seq) improves diagnostic rates in individuals with suspected Mendelian conditions to varying degrees, primarily by directing the prioritization of candidate DNA variants identified on exome or genome sequencing (ES/GS). Here we implemented an RNA-seq-guided method to diagnose individuals across a wide range of ages and clinical phenotypes.METHODSOne hundred fifteen undiagnosed adult and pediatric patients with diverse phenotypes and 67 family members (182 total individuals) underwent RNA-seq from whole blood and skin fibroblasts at the Baylor College of Medicine (BCM) Undiagnosed Diseases Network clinical site from 2014 to 2020. We implemented a workflow to detect outliers in gene expression and splicing for cases that remained undiagnosed despite standard genomic and transcriptomic analysis.RESULTSThe transcriptome-directed approach resulted in a diagnostic rate of 12% across the entire cohort, or 17% after excluding cases solved on ES/GS alone. Newly diagnosed conditions included Koolen-de Vries syndrome (KANSL1), Renpenning syndrome (PQBP1), TBCK-associated encephalopathy, NSD2- and CLTC-related intellectual disability, and others, all with negative conventional genomic testing, including ES and chromosomal microarray (CMA). Skin fibroblasts exhibited higher and more consistent expression of clinically relevant genes than whole blood. In solved cases with RNA-seq from both tissues, the causative defect was missed in blood in half the cases but none from fibroblasts.CONCLUSIONSFor our cohort of undiagnosed individuals with suspected Mendelian conditions, transcriptome-directed genomic analysis facilitated diagnoses, primarily through the identification of variants missed on ES and CMA.TRIAL REGISTRATIONNot applicable.FUNDINGNIH Common Fund, BCM Intellectual and Developmental Disabilities Research Center, Eunice Kennedy Shriver National Institute of Child Health & Human Development.

Intellectual disability-associated gain-of-function mutations in CERT1 that encodes the ceramide transport protein CERT

Intellectual disability (ID) is a developmental disorder that includes both intellectual and adaptive functioning deficits in conceptual, social, and practical domains. Although evidence-based interventions for patients have long been desired, their progress has been hindered due to various determinants. One of these determinants is the complexity of the origins of ID. The ceramide transport protein (CERT) encoded by CERT1 mediates inter-organelle trafficking of ceramide for the synthesis of intracellular sphingomyelin. Utilizing whole exome sequencing analysis, we identified a novel CERT variant, which substitutes a serine at position 135 (S135) for a proline in a patient with severe ID. Biochemical analysis showed that S135 is essential for hyperphosphorylation of a serine-repeat motif of CERT, which is required for down-regulation of CERT activity. Amino acid replacements of S135 abnormally activated CERT and induced an intracellular punctate distribution pattern of this protein. These results identified specific ID-associated CERT1 mutations that induced gain-of-function effects on CERT activity. These findings provide a possible molecular basis for not only new diagnostics but also a conceivable pharmaceutical intervention for ID disorders caused by gain-of-function mutations in CERT1.

BICRA, a SWI/SNF Complex Member, Is Associated with BAF-Disorder Related Phenotypes in Humans and Model Organisms

SWI/SNF-related intellectual disability disorders (SSRIDDs) are rare neurodevelopmental disorders characterized by developmental disability, coarse facial features, and fifth digit/nail hypoplasia that are caused by pathogenic variants in genes that encode for members of the SWI/SNF (or BAF) family of chromatin remodeling complexes. We have identified 12 individuals with rare variants (10 loss-of-function, 2 missense) in the BICRA (BRD4 interacting chromatin remodeling complex-associated protein) gene, also known as GLTSCR1, which encodes a subunit of the non-canonical BAF (ncBAF) complex. These individuals exhibited neurodevelopmental phenotypes that include developmental delay, intellectual disability, autism spectrum disorder, and behavioral abnormalities as well as dysmorphic features. Notably, the majority of individuals lack the fifth digit/nail hypoplasia phenotype, a hallmark of most SSRIDDs. To confirm the role of BICRA in the development of these phenotypes, we performed functional characterization of the zebrafish and Drosophila orthologs of BICRA. In zebrafish, a mutation of bicra that mimics one of the loss-of-function variants leads to craniofacial defects possibly akin to the dysmorphic facial features seen in individuals harboring putatively pathogenic BICRA variants. We further show that Bicra physically binds to other non-canonical ncBAF complex members, including the BRD9/7 ortholog, CG7154, and is the defining member of the ncBAF complex in flies. Like other SWI/SNF complex members, loss of Bicra function in flies acts as a dominant enhancer of position effect variegation but in a more context-specific manner. We conclude that haploinsufficiency of BICRA leads to a unique SSRIDD in humans whose phenotypes overlap with those previously reported.

JARID2 haploinsufficiency is associated with a clinically distinct neurodevelopmental syndrome

PURPOSE

JARID2, located on chromosome 6p22.3, is a regulator of histone methyltransferase complexes that is expressed in human neurons. So far, 13 individuals sharing clinical features including intellectual disability (ID) were reported with de novo heterozygous deletions in 6p22-p24 encompassing the full length JARID2 gene (OMIM 601594). However, all published individuals to date have a deletion of at least one other adjoining gene, making it difficult to determine if JARID2 is the critical gene responsible for the shared features. We aim to confirm JARID2 as a human disease gene and further elucidate the associated clinical phenotype.

METHODS

Chromosome microarray analysis, exome sequencing, and an online matching platform (GeneMatcher) were used to identify individuals with single-nucleotide variants or deletions involving JARID2.

RESULTS

We report 16 individuals in 15 families with a deletion or single-nucleotide variant in JARID2. Several of these variants are likely to result in haploinsufficiency due to nonsense-mediated messenger RNA (mRNA) decay. All individuals have developmental delay and/or ID and share some overlapping clinical characteristics such as facial features with those who have larger deletions involving JARID2.

CONCLUSION

We report that JARID2 haploinsufficiency leads to a clinically distinct neurodevelopmental syndrome, thus establishing gene-disease validity for the purpose of diagnostic reporting.

De Novo KAT5 Variants Cause a Syndrome with Recognizable Facial Dysmorphisms, Cerebellar Atrophy, Sleep Disturbance, and Epilepsy

KAT5 encodes an essential lysine acetyltransferase, previously called TIP60, which is involved in regulating gene expression, DNA repair, chromatin remodeling, apoptosis, and cell proliferation; but it remains unclear whether variants in this gene cause a genetic disease. Here, we study three individuals with heterozygous de novo missense variants in KAT5 that affect normally invariant residues, with one at the chromodomain (p.Arg53His) and two at or near the acetyl-CoA binding site (p.Cys369Ser and p.Ser413Ala). All three individuals have cerebral malformations, seizures, global developmental delay or intellectual disability, and severe sleep disturbance. Progressive cerebellar atrophy was also noted. Histone acetylation assays with purified variant KAT5 demonstrated that the variants decrease or abolish the ability of the resulting NuA4/TIP60 multi-subunit complexes to acetylate the histone H4 tail in chromatin. Transcriptomic analysis in affected individual fibroblasts showed deregulation of multiple genes that control development. Moreover, there was also upregulated expression of PER1 (a key gene involved in circadian control) in agreement with sleep anomalies in all of the individuals. In conclusion, dominant missense KAT5 variants cause histone acetylation deficiency with transcriptional dysregulation of multiples genes, thereby leading to a neurodevelopmental syndrome with sleep disturbance, cerebellar atrophy, and facial dysmorphisms, and suggesting a recognizable syndrome.

Community-based recruitment and exome sequencing indicates high diagnostic yield in adults with intellectual disability

BACKGROUND

Establishing a genetic diagnosis for individuals with intellectual disability (ID) benefits patients and their families as it may inform the prognosis, lead to appropriate therapy, and facilitate access to medical and supportive services. Exome sequencing has been successfully applied in a diagnostic setting, but most clinical exome referrals are pediatric patients, with many adults with ID lacking a comprehensive genetic evaluation.

METHODS

Our unique recruitment strategy involved partnering with service and education providers for individuals with ID. We performed exome sequencing and analysis, and clinical variant interpretation for each recruited family.

RESULTS

All five families enrolled in the study opted-in for the return of genetic results. In three out of five families exome sequencing analysis identified pathogenic or likely pathogenic variants in KANSL1, TUSC3, and MED13L genes. Families discussed the results and any potential medical follow-up in an appointment with a board certified clinical geneticist.

CONCLUSION

Our study suggests high yield of exome sequencing as a diagnostic tool in adult patients with ID who have not undergone comprehensive sequencing-based genetic testing. Research studies including an option of return of results through a genetic clinic could help minimize the disparity in exome diagnostic testing between pediatric and adult patients with ID.

A Synonymous Exonic Splice Silencer Variant in IRF6 as a Novel and Cryptic Cause of Non-Syndromic Cleft Lip and Palate

Missense, nonsense, splice site and regulatory region variants in interferon regulatory factor 6 (IRF6) have been shown to contribute to both syndromic and non-syndromic forms of cleft lip and/or palate (CL/P). We report the diagnostic evaluation of a complex multigeneration family of Honduran ancestry with a pedigree structure consistent with autosomal-dominant inheritance with both incomplete penetrance and variable expressivity. The proband's grandmother bore children with two partners and CL/P segregates on both sides of each lineage. Through whole-exome sequencing of five members of the family, we identified a single shared synonymous variant, located in the middle of exon 7 of IRF6 (p.Ser307Ser; g.209963979 G>A; c.921C>T). The variant was shown to segregate in the seven affected individuals and through three unaffected obligate carriers, spanning both sides of this pedigree. This variant is very rare, only being found in three (all of Latino ancestry) of 251,352 alleles in the gnomAD database. While the variant did not create a splice acceptor/donor site, in silico analysis predicted it to impact an exonic splice silencer element and the binding of major splice regulatory factors. In vitro splice assays supported this by revealing multiple abnormal splicing events, estimated to impact >60% of allelic transcripts. Sequencing of the alternate splice products demonstrated the unmasking of a cryptic splice site six nucleotides 5' of the variant, as well as variable utilization of cryptic splice sites in intron 6. The ectopic expression of different splice regulatory proteins altered the proportion of abnormal splicing events seen in the splice assay, although the alteration was dependent on the splice factor. Importantly, each alternatively spliced mRNA is predicted to result in a frame shift and prematurely truncated IRF6 protein. This is the first study to identify a synonymous variant as a likely cause of NS-CL/P and highlights the care that should be taken by laboratories when considering and interpreting variants.

Biallelic variants in ADARB1, encoding a dsRNA-specific adenosine deaminase, cause a severe developmental and epileptic encephalopathy

BACKGROUND

Adenosine-to-inosine RNA editing is a co-transcriptional/post-transcriptional modification of double-stranded RNA, catalysed by one of two active adenosine deaminases acting on RNA (ADARs), ADAR1 and ADAR2. ADARB1 encodes the enzyme ADAR2 that is highly expressed in the brain and essential to modulate the function of glutamate and serotonin receptors. Impaired ADAR2 editing causes early onset progressive epilepsy and premature death in mice. In humans, ADAR2 dysfunction has been very recently linked to a neurodevelopmental disorder with microcephaly and epilepsy in four unrelated subjects.

METHODS

We studied three children from two consanguineous families with severe developmental and epileptic encephalopathy (DEE) through detailed physical and instrumental examinations. Exome sequencing (ES) was used to identify ADARB1 mutations as the underlying genetic cause and in vitro assays with transiently transfected cells were performed to ascertain the impact on ADAR2 enzymatic activity and splicing.

RESULTS

All patients showed global developmental delay, intractable early infantile-onset seizures, microcephaly, severe-to-profound intellectual disability, axial hypotonia and progressive appendicular spasticity. ES revealed the novel missense c.1889G>A, p.(Arg630Gln) and deletion c.1245_1247+1 del, p.(Leu415PhefsTer14) variants in ADARB1 (NM_015833.4). The p.(Leu415PhefsTer14) variant leads to incorrect splicing resulting in frameshift with a premature stop codon and loss of enzyme function. In vitro RNA editing assays showed that the p.(Arg630Gln) variant resulted in a severe impairment of ADAR2 enzymatic activity.

CONCLUSION

In conclusion, these data support the pathogenic role of biallelic ADARB1 variants as the cause of a distinctive form of DEE, reinforcing the importance of RNA editing in brain function and development.

12q21 deletion syndrome: Narrowing the critical region down to 1.6 Mb including SYT1 and PPP1R12A

Deletions in the 12q21 region are rare and non-recurrent CNVs. To date, only 11 patients with deletions in this region have been reported in the literature. These patients most often presented with syndromic intellectual deficiency, ventriculomegaly or hydrocephalus, ectodermal abnormalities, growth retardation and renal and cardiac malformations, suggesting a recognizable microdeletion syndrome. We report three new patients with overlapping deletions of the 12q21 region, including the smallest deletion reported to date and the first case characterized by array CGH during pregnancy. We describe specific clinical findings and shared facial features as developmental delay, ectodermal abnormalities, ventriculomegaly or hydrocephalus, axial hypotonia or spastic diplegia, growth retardation, heart defect, hydronephrosis, ureteral reflux or horseshoe kidney, large thorax or pectus excavatum, syndactyly of 2-3 toes, pterygium coli or excess nuchal skin, large anterior fontanel, low set ears, prominent forehead, short-upturned nose with nostril hypoplasia, microretrognathia and hypertelorism. These new patients and a comprehensive review of the literature allow us to define a minimum critical region spanning 1.6 Mb in 12q21. By screening the critical region using prediction tools, we identified two candidate genes: SYT1and PPP1R12A.

Grading of Clear Cell Renal Cell Carcinomas by Using Machine Learning Based on Artificial Neural Networks and Radiomic Signatures Extracted From Multidetector Computed Tomography Images

RATIONALE AND OBJECTIVES

To evaluate the ability of artificial neural networks (ANN) fed with radiomic signatures (RSs) extracted from multidetector computed tomography images in differentiating the histopathological grades of clear cell renal cell carcinomas (ccRCCs).

MATERIALS AND METHODS

The multidetector computed tomography images of 227 ccRCCs were retrospectively analyzed. For each ccRCC, 14 conventional image features (CIFs) were extracted manually by two radiologists, and 556 texture features (TFs) were extracted by a free software application, MaZda (version 4.6). The high-dimensional dataset of these RSs was reduced using the least absolute shrinkage and selection operator. Five minimum mean squared error models (minMSEMs) for predicting the ccRCC histopathological grades were constructed from the CIFs, the TFs of the corticomedullary phase images (CMP), and the TFs of the parenchyma phase (PP) images and their combinations, respectively abbreviated as CIF-minMSEM, CMP-minMSEM, PP-minMSEM, CIF+CMP-minMSEM, and CIF+PP-minMSEM. The RSs of each model were fed 30 times consecutively into an ANN for machine learning, and the predictive accuracy of each time ML was recorded for the statistical analysis.

RESULTS

The five predictive models were constructed from 12, 19, and 10 features selected from the CIFs, the TFs of the CMP images, and that of PP images, respectively. On the basis of their accuracy across the whole cohort, the five models were ranked as follows: CIF+CMP-minMSEM (accuracy: 94.06% ± 1.14%), CIF + PP-minMSEM (accuracy: 93.32% ± 1.23%), CIF-minMSEM (accuracy: 92.26% ± 1.65%), CMP-minMSEM (accuracy: 91.76% ± 1.74%), and PP-minMSEM (accuracy: 90.89% ± 1.47%).

CONCLUSION

Machine learning based on ANN helped establish an optimal predictive model, and TFs contributed to the development of high accuracy predictive models. The CIF+CMP-minMSEM showed the greatest accuracy for differentiating low- and high-grade ccRCCs.