Prioritization of oligogenic variant combinations in whole exomes.

MOTIVATION: Whole exome sequencing (WES) has emerged as a powerful tool for genetic research, enabling the collection of a tremendous amount of data about human genetic variation. However, properly identifying which variants are causative of a genetic disease remains an important challenge, often due to the number of variants that need to be screened. Expanding the screening to combinations of variants in two or more genes, as would be required under the oligogenic inheritance model, simply blows this problem out of proportion. RESULTS: We present here the High-throughput oligogenic prioritizer (Hop), a novel prioritization method that uses direct oligogenic information at the variant, gene and gene pair level to detect digenic variant combinations in WES data. This method leverages information from a knowledge graph, together with specialized pathogenicity predictions in order to effectively rank variant combinations based on how likely they are to explain the patient's phenotype. The performance of Hop is evaluated in cross-validation on 36 120 synthetic exomes for training and 14 280 additional synthetic exomes for independent testing. Whereas the known pathogenic variant combinations are found in the top 20 in approximately 60% of the cross-validation exomes, 71% are found in the same ranking range when considering the independent set. These results provide a significant improvement over alternative approaches that depend simply on a monogenic assessment of pathogenicity, including early attempts for digenic ranking using monogenic pathogenicity scores. AVAILABILITY AND IMPLEMENTATION: Hop is available at https://github.com/oligogenic/HOP.

A milder form of NSRP1-associated neurodevelopmental disorder, caused by a missense variant in the nuclear localization signal.

Nuclear Speckle Splicing Regulator Protein 1 (NSRP1) is a splice factor found in nuclear speckles, which are small membrane-free organelles implicated in epigenetic regulation, chromatin organization, DNA repair, and RNA modification. Bi-allelic loss-of-function variants in NSRP1 have recently been identified in patients suffering from a severe neurodevelopmental disorder, presenting with neurodevelopmental delay, epilepsy, microcephaly, hypotonia, and spastic cerebral palsy. Described patients acquired neither independent walking nor speech and often showed anomalies on cerebral MRI. Here we describe the case of a 14-year-old girl with motor and language delay as well as intellectual disability, who presents an ataxic gait but walks without assistance and speaks in short sentences. Whole-genome sequencing revealed the compound heterozygous NSRP1 variants c.114 + 2T > G and c.1595T > A (p.Val532Glu). Functional validation using HEK293T cells transfected with either wild-type or mutated GFP-tagged Nsrp1 suggests that the Val532Glu variant interferes with the function of the nuclear localization signal, and leads to mislocalization of NSRP1 in the cytosol, thus confirming the pathogenicity of the observed variant. This case helps to expand the phenotypic and genetic spectrum associated with pathogenic NSRP1 variants and indicates that this diagnosis should also be suspected in patients with milder phenotypes.

Comprehensive evaluation of the implementation of episignatures for diagnosis of neurodevelopmental disorders (NDDs).

Episignatures are popular tools for the diagnosis of rare neurodevelopmental disorders. They are commonly based on a set of differentially methylated CpGs used in combination with a support vector machine model. DNA methylation (DNAm) data often include missing values due to changes in data generation technology and batch effects. While many normalization methods exist for DNAm data, their impact on episignature performance have never been assessed. In addition, technologies to quantify DNAm evolve quickly and this may lead to poor transposition of existing episignatures generated on deprecated array versions to new ones. Indeed, probe removal between array versions, technologies or during preprocessing leads to missing values. Thus, the effect of missing data on episignature performance must also be carefully evaluated and addressed through imputation or an innovative approach to episignatures design. In this paper, we used data from patients suffering from Kabuki and Sotos syndrome to evaluate the influence of normalization methods, classification models and missing data on the prediction performances of two existing episignatures. We compare how six popular normalization methods for methylarray data affect episignature classification performances in Kabuki and Sotos syndromes and provide best practice suggestions when building new episignatures. In this setting, we show that Illumina, Noob or Funnorm normalization methods achieved higher classification performances on the testing sets compared to Quantile, Raw and Swan normalization methods. We further show that penalized logistic regression and support vector machines perform best in the classification of Kabuki and Sotos syndrome patients. Then, we describe a new paradigm to build episignatures based on the detection of differentially methylated regions (DMRs) and evaluate their performance compared to classical differentially methylated cytosines (DMCs)-based episignatures in the presence of missing data. We show that the performance of classical DMC-based episignatures suffers from the presence of missing data more than the DMR-based approach. We present a comprehensive evaluation of how the normalization of DNA methylation data affects episignature performance, using three popular classification models. We further evaluate how missing data affect those models' predictions. Finally, we propose a novel methodology to develop episignatures based on differentially methylated regions identification and show how this method slightly outperforms classical episignatures in the presence of missing data.

Clinical and molecular diagnosis of genodermatoses: Review and perspectives.

Genodermatoses are a complex and heterogeneous group of genetic skin disorders characterized by variable expression and clinical and genetic heterogeneity, rendering their diagnosis challenging. DNA-based techniques, like whole-exome sequencing, can establish a diagnosis in 50% of cases. RNA-sequencing is emerging as an attractive tool that can obtain information regarding gene expression while integrating functional genomic data with regard to the interpretation of variants. This increases the diagnostic rate by an additional 10-15%. In the present review, we detail the clinical steps involved in the diagnosis of genodermatoses, as well as the current DNA-based technologies available to clinicians. Herein, the intention is to facilitate a better understanding of the possibilities and limitations of these diagnostic technologies. In addition, this review could guide dermatologists through new emerging techniques, such as RNA-sequencing and its applications to familiarizing them with future techniques. Currently, this multi-omics approach is likely the best strategy designed to promote the diagnosis of patients with genodermatoses and discover new skin disease genes that could result in novel targeted therapies.

Implementation of fetal clinical exome sequencing: Comparing prospective and retrospective cohorts.

PURPOSE: We compared the diagnostic yield of fetal clinical exome sequencing (fCES) in prospective and retrospective cohorts of pregnancies presenting with anomalies detected using ultrasound. We evaluated factors that led to a higher diagnostic efficiency, such as phenotypic category, clinical characterization, and variant analysis strategy. METHODS: fCES was performed for 303 fetuses (183 ongoing and 120 ended pregnancies, in which chromosomal abnormalities had been excluded) using a trio/duo-based approach and a multistep variant analysis strategy. RESULTS: fCES identified the underlying genetic cause in 13% (24/183) of prospective and 29% (35/120) of retrospective cases. In both cohorts, recessive heterozygous compound genotypes were not rare, and trio and simplex variant analysis strategies were complementary to achieve the highest possible diagnostic rate. Limited prenatal phenotypic information led to interpretation challenges. In 2 prospective cases, in-depth analysis allowed expansion of the spectrum of prenatal presentations for genetic syndromes associated with the SLC17A5 and CHAMP1 genes. CONCLUSION: fCES is diagnostically efficient in fetuses presenting with cerebral, skeletal, urinary, or multiple anomalies. The comparison between the 2 cohorts highlights the importance of providing detailed phenotypic information for better interpretation and prenatal reporting of genetic variants.

Comparative study of keratinocyte primary culture methods from paediatric skin biopsies for RNA-sequencing.

BACKGROUND: Keratinocyte culture is a standard method used to study gene expression, cell differentiation and proliferation. Numerous protocols exist, however their application is frequently unsuitable for small specimens, such as 4-mm punch skin biopsies. AIMS: This study compared 3 different methods of keratinocyte culture from paediatric skin biopsies to evaluate which one ensures adequate cell growth for RNA extraction and sequencing. MATERIALS AND METHODS: Thirty-six skin samples were obtained from 4-mm punch skin biopsies from residual human body material from healthy children. They were cultured in vitro according to 3 different methods: enzymatic method, epidermis explant and direct explant method. Keratinocytes were characterized by immunocytochemistry using pan-cytokeratin. RNA extraction was performed with RNeasy Mini kit. Quantity and quality of the extracted RNA was assessed to meet the requirements of library preparation for sequencing. RESULTS: The direct explant method had largely shown its superiority over the two other methods, with a 100% success rate and an average of 15 days of culture. RNA extraction yielded a mean of 8545.85 ng of RNA per sample with an RQN of 10. Cover-clip immunochemistry staining with pan-cytokeratin had confirmed the absence of fibroblast contamination. DISCUSSION: Although the enzymatic method is the most frequently used for keratinocyte culture, it is not suitable small samples required in dermatology. The direct explant method guarantees a high growth rate and the extraction of high quality RNA. Variation in the amount of RNA harvested are related to inter- and intra-individual variations and to the conditions of the experiment. CONCLUSION: This study allowed to conclude that the direct explant method is the most efficient and easy method to ensure cell growth when the samples are from 4-mm punch skin biopsies. This technique avoids fibroblasts contamination and obtains a sufficient quantity and quality of RNA to sequence it.

Predicting disease-causing variant combinations.

Notwithstanding important advances in the context of single-variant pathogenicity identification, novel breakthroughs in discerning the origins of many rare diseases require methods able to identify more complex genetic models. We present here the Variant Combinations Pathogenicity Predictor (VarCoPP), a machine-learning approach that identifies pathogenic variant combinations in gene pairs (called digenic or bilocus variant combinations). We show that the results produced by this method are highly accurate and precise, an efficacy that is endorsed when validating the method on recently published independent disease-causing data. Confidence labels of 95% and 99% are identified, representing the probability of a bilocus combination being a true pathogenic result, providing geneticists with rational markers to evaluate the most relevant pathogenic combinations and limit the search space and time. Finally, the VarCoPP has been designed to act as an interpretable method that can provide explanations on why a bilocus combination is predicted as pathogenic and which biological information is important for that prediction. This work provides an important step toward the genetic understanding of rare diseases, paving the way to clinical knowledge and improved patient care.

Evaluation of neurodevelopmental symptoms in 10 cases of neonatal ichthyosis and sclerosing cholangitis syndrome.

Neonatal ichthyosis and sclerosing cholangitis (NISCH) syndrome is an extremely rare entity with only 19 patients described in the literature. We report an extended family with the disorder and investigate the association of neurodevelopmental symptoms. Patients with CLDN1 mutations, and specifically « the Moroccan¼ c.200_201delTT deletion, may be an increased risk for neurodevelopmental symptoms such as learning disabilities, mental retardation, and language delay.

Rarity of fetal cells in exocervical samples for noninvasive prenatal diagnosis.

OBJECTIVES: The possibility to isolate fetal cells from pregnant women cervical samples has been discussed for five decades but is not currently applied in clinical practice. This study aimed at offering prenatal genetic diagnosis from fetal cells obtained through noninvasive exocervical sampling and immuno-sorted based on expression of HLA-G. METHODS: We first developed and validated robust protocols for cell detection and isolation on control cell lines expressing (JEG-3) or not (JAR) the HLA-G antigen, a specific marker for extravillous trophoblasts. We then applied these protocols to noninvasive exocervical samples collected from pregnant women between 6 and 14 weeks of gestational age. Sampling was performed through insertion and rotation of a brush at the ectocervix close to the external os of the endocervical canal. Finally, we attempted to detect and quantify trophoblasts in exocervical samples from pregnant women by ddPCR targeting the male SRY locus. RESULTS: For immunohistochemistry, a strong specific signal for HLA-G was observed in the positive control cell line and for rare cells in exocervical samples, but only in non-fixative conditions. HLA-G positive cells diluted in HLA-G negative cells were isolated by flow cytometry or magnetic cell sorting. However, no HLA-G positive cells could be recovered from exocervical samples. SRY gene was detected by ddPCR in exocervical samples from male (50%) but also female (27%) pregnancies. CONCLUSIONS: Our data suggest that trophoblasts are too rarely and inconstantly present in noninvasive exocervical samples to be reliably retrieved by standard immunoisolation techniques and therefore cannot replace the current practice for prenatal screening and diagnosis.

Outcome of publicly funded nationwide first-tier noninvasive prenatal screening.

PURPOSE: Noninvasive prenatal screening (NIPS) using cell-free DNA has transformed prenatal care. Belgium was the first country to implement and fully reimburse NIPS as a first-tier screening test offered to all pregnant women. A consortium consisting of all Belgian genetic centers report the outcome of two years genome-wide NIPS implementation. METHODS: The performance for the common trisomies and for secondary findings was evaluated based on 153,575 genome-wide NIP tests. Furthermore, the evolution of the number of invasive tests and the incidence of Down syndrome live births was registered. RESULTS: Trisomies 21, 18, and 13 were detected in respectively 0.32%, 0.07%, and 0.06% of cases, with overall positive predictive values (PPVs) of 92.4%, 84.6%, and 43.9%. Rare autosomal trisomies and fetal segmental imbalances were detected in respectively 0.23% and 0.07% of cases with PPVs of 4.1% and 47%. The number of invasive obstetric procedures decreased by 52%. The number of trisomy 21 live births dropped to 0.04%. CONCLUSION: Expanding the scope of NIPS beyond trisomy 21 fetal screening allows the implementation of personalized genomic medicine for the obstetric population. This genome-wide NIPS approach has been embedded successfully in prenatal genetic care in Belgium and might serve as a framework for other countries offering NIPS.

ORVAL: a novel platform for the prediction and exploration of disease-causing oligogenic variant combinations.

A tremendous amount of DNA sequencing data is being produced around the world with the ambition to capture in more detail the mechanisms underlying human diseases. While numerous bioinformatics tools exist that allow the discovery of causal variants in Mendelian diseases, little to no support is provided to do the same for variant combinations, an essential task for the discovery of the causes of oligogenic diseases. ORVAL (the Oligogenic Resource for Variant AnaLysis), which is presented here, provides an answer to this problem by focusing on generating networks of candidate pathogenic variant combinations in gene pairs, as opposed to isolated variants in unique genes. This online platform integrates innovative machine learning methods for combinatorial variant pathogenicity prediction with visualization techniques, offering several interactive and exploratory tools, such as pathogenic gene and protein interaction networks, a ranking of pathogenic gene pairs, as well as visual mappings of the cellular location and pathway information. ORVAL is the first web-based exploration platform dedicated to identifying networks of candidate pathogenic variant combinations with the sole ambition to help in uncovering oligogenic causes for patients that cannot rely on the classical disease analysis tools. ORVAL is available at https://orval.ibsquare.be.

Novel homozygous variant of carbonic anhydrase 8 gene expanding the phenotype of cerebellar ataxia, mental retardation, and disequilibrium syndrome subtype 3.

We report the case of an 11-year-old Syrian girl born to consanguineous parents, who presents an ataxic gait from early childhood. On clinical examination, she presented a severe static - kinetic cerebellar syndrome, walking without support is possible for short distances only. Strikingly, three consecutive MRIs did not show any sign of cerebellar abnormalities, but a brain positron emission tomography (PET) using [18F]-fluorodeoxyglucose (FDG) demonstrated a clear decrease in glucose metabolism in the cerebellum as well as the anterior and medial temporal lobe bilaterally. A clinical exome analysis identified a novel homozygous c.251A > G (p.Asn84Ser) likely pathogenic variant in the carbonic anhydrase 8 (CA8) gene. CA8 mutations cause cerebellar ataxia, mental retardation, and disequilibrium syndrome subtype 3 (CAMRQ3), a rare genetically autosomal recessive disorder, only described in four families, so far with the frequent observation of quadrupedal gait. The proband differed with other reported CA8 mutations by the absence of clear cerebellar signs on brain MRI and the presence of focal seizures. This report expands the clinical spectrum associated with mutations in CA8 and illustrates the possible discrepancy between (mild) neuro-radiological images (MRI) and (severe) clinical phenotype in young individuals. In contrast, the observation of clear cerebellar abnormal metabolic findings suggests that the FDG-PET scan may be used as an early marker for hereditary ataxia.

NOVOPlasty: de novo assembly of organelle genomes from whole genome data.

The evolution in next-generation sequencing (NGS) technology has led to the development of many different assembly algorithms, but few of them focus on assembling the organelle genomes. These genomes are used in phylogenetic studies, food identification and are the most deposited eukaryotic genomes in GenBank. Producing organelle genome assembly from whole genome sequencing (WGS) data would be the most accurate and least laborious approach, but a tool specifically designed for this task is lacking. We developed a seed-and-extend algorithm that assembles organelle genomes from whole genome sequencing (WGS) data, starting from a related or distant single seed sequence. The algorithm has been tested on several new (Gonioctena intermedia and Avicennia marina) and public (Arabidopsis thaliana and Oryza sativa) whole genome Illumina data sets where it outperforms known assemblers in assembly accuracy and coverage. In our benchmark, NOVOPlasty assembled all tested circular genomes in less than 30 min with a maximum memory requirement of 16 GB and an accuracy over 99.99%. In conclusion, NOVOPlasty is the sole de novo assembler that provides a fast and straightforward extraction of the extranuclear genomes from WGS data in one circular high quality contig. The software is open source and can be downloaded at https://github.com/ndierckx/NOVOPlasty.

Understanding mutational effects in digenic diseases.

To further our understanding of the complexity and genetic heterogeneity of rare diseases, it has become essential to shed light on how combinations of variants in different genes are responsible for a disease phenotype. With the appearance of a resource on digenic diseases, it has become possible to evaluate how digenic combinations differ in terms of the phenotypes they produce. All instances in this resource were assigned to two classes of digenic effects, annotated as true digenic and composite classes. Whereas in the true digenic class variants in both genes are required for developing the disease, in the composite class, a variant in one gene is sufficient to produce the phenotype, but an additional variant in a second gene impacts the disease phenotype or alters the age of onset. We show that a combination of variant, gene and higher-level features can differentiate between these two classes with high accuracy. Moreover, we show via the analysis of three digenic disorders that a digenic effect decision profile, extracted from the predictive model, motivates why an instance was assigned to either of the two classes. Together, our results show that digenic disease data generates novel insights, providing a glimpse into the oligogenic realm.

BCL11A frameshift mutation associated with dyspraxia and hypotonia affecting the fine, gross, oral, and speech motor systems.

We report the case of a 7-year-old male of Western European origin presenting with moderate intellectual disability, severe childhood apraxia of speech in the presence of oral and manual dyspraxia, and hypotonia across motor systems including the oral and speech motor systems. Exome sequencing revealed a de novo frameshift protein truncating mutation in the fourth exon of BCL11A, a gene recently demonstrated as being involved in cognition and language development. Making parallels with a previously described patient with a 200 kb 2p15p16.1 deletion encompassing the entire BCL11A gene and displaying a similar phenotype, we characterize in depth how BCL11A is involved in clinical aspects of language development and oral praxis.

DIDA: A curated and annotated digenic diseases database.

DIDA (DIgenic diseases DAtabase) is a novel database that provides for the first time detailed information on genes and associated genetic variants involved in digenic diseases, the simplest form of oligogenic inheritance. The database is accessible via http://dida.ibsquare.be and currently includes 213 digenic combinations involved in 44 different digenic diseases. These combinations are composed of 364 distinct variants, which are distributed over 136 distinct genes. The web interface provides browsing and search functionalities, as well as documentation and help pages, general database statistics and references to the original publications from which the data have been collected. The possibility to submit novel digenic data to DIDA is also provided. Creating this new repository was essential as current databases do not allow one to retrieve detailed records regarding digenic combinations. Genes, variants, diseases and digenic combinations in DIDA are annotated with manually curated information and information mined from other online resources. Next to providing a unique resource for the development of new analysis methods, DIDA gives clinical and molecular geneticists a tool to find the most comprehensive information on the digenic nature of their diseases of interest.

FGFR1 mutations cause Hartsfield syndrome, the unique association of holoprosencephaly and ectrodactyly.

BACKGROUND: Harstfield syndrome is the rare and unique association of holoprosencephaly (HPE) and ectrodactyly, with or without cleft lip and palate, and variable additional features. All the reported cases occurred sporadically. Although several causal genes of HPE and ectrodactyly have been identified, the genetic cause of Hartsfield syndrome remains unknown. We hypothesised that a single key developmental gene may underlie the co-occurrence of HPE and ectrodactyly. METHODS: We used whole exome sequencing in four isolated cases including one case-parents trio, and direct Sanger sequencing of three additional cases, to investigate the causative variants in Hartsfield syndrome. RESULTS: We identified a novel FGFR1 mutation in six out of seven patients. Affected residues are highly conserved and are located in the extracellular binding domain of the receptor (two homozygous mutations) or the intracellular tyrosine kinase domain (four heterozygous de novo variants). Strikingly, among the six novel mutations, three are located in close proximity to the ATP's phosphates or the coordinating magnesium, with one position required for kinase activity, and three are adjacent to known mutations involved in Kallmann syndrome plus other developmental anomalies. CONCLUSIONS: Dominant or recessive FGFR1 mutations are responsible for Hartsfield syndrome, consistent with the known roles of FGFR1 in vertebrate ontogeny and conditional Fgfr1-deficient mice. Our study shows that, in humans, lack of accurate FGFR1 activation can disrupt both brain and hand/foot midline development, and that FGFR1 loss-of-function mutations are responsible for a wider spectrum of clinical anomalies than previously thought, ranging in severity from seemingly isolated hypogonadotropic hypogonadism, through Kallmann syndrome with or without additional features, to Hartsfield syndrome at its most severe end.

A novel homozygous variant in SLC25A46 gene associated with pontocerebellar hypoplasia type 1E: a case report.

Neonatal encephalopathy (NE) is a complex clinical condition with diverse etiologies. Hypoxic-ischemic encephalopathy (HIE) is a major contributor to NE cases. However, distinguishing NE subtypes, such as pontocerebellar hypoplasia type 1E (PCH1E), from HIE can be challenging due to overlapping clinical features. Here, we present a case of PCH1E in a neonate with a homozygous mutation c.72delT p. (Phe24LeufsTer20) in the SLC25A46 gene. The severity of PCH1E associated NE highlighted the significance of early recognition to guide appropriate clinical management.

POT1 tumour predisposition: a broader spectrum of associated malignancies and proposal for additional screening program.

Protection of Telomeres Protein 1 (POT1) protein is an essential subunit of the shelterin telomere binding complex, regulating telomere length. Some POT1 gene pathogenic variants (PV) lead to telomere elongation, genomic instability and higher risk of cancer. POT1 tumour predisposition syndrome (POT1-TPD) has autosomal dominant inheritance and unknown penetrance. It is associated with increased risk of cutaneous melanoma, chronic lymphocytic leukaemia, angiosarcoma and gliomas. In this work, we aim to describe a broader cancer phenotype related to POT1-TPD, in three families (two with a four generation pedigree, one with a five generation pedigree). The three index cases were referred to our oncogenetic centre for genetic counselling due to their personal history of cancer. Two underwent clinical exome sequencing of 4,867 genes associated with Mendelian genetic diseases, and another underwent gene panel sequencing including POT1, which identified three different POT1 PV: NC_000007.14(NM_015450.2):c.349C>T; NC_000007.14(NM_015450.2):c.233T>C and NC_000007.14(NM_015450.2):c.818G>A; already described in the literature. Referenced relatives, did a target genetic test (according to the POT1 PV identified in the family). In total, 37 individuals were tested (51.4% females), median age of 46 (22-81) years, with POT1 PV detected in 22. POT1-TPD was observed, but also a higher incidence of other cancers (other sarcomas, papillary thyroid cancer, early onset prostate cancer and leukaemia). These findings contribute to an increase in our knowledge about POT1 PV, and it can play a role in the definition of future POT1 PV screening criteria, POT1 carrier surveillance protocols (possibly considering screening for all types of sarcomas) and in genetic counselling.