Machine learning-enhanced noninvasive prenatal testing of monogenic disorders.

OBJECTIVE: Single-nucleotide variants (SNVs) are of great significance in prenatal diagnosis as they are the leading cause of inherited single-gene disorders (SGDs). Identifying SNVs in a non-invasive prenatal screening (NIPS) scenario is particularly challenging for maternally inherited SNVs. We present an improved method to predict inherited SNVs from maternal or paternal origin in a genome-wide manner. METHODS: We performed SNV-NIPS based on the combination of fragments of cell free DNA (cfDNA) features, Bayesian inference and a machine-learning (ML) prediction refinement step using random forest (RF) classifiers trained on millions of non-pathogenic variants. We next evaluate the real-world performance of our refined method in a clinical setting by testing our models on 16 families with singleton pregnancies and varying fetal fraction (FF) levels, and validate the results over millions of inherited variants in each fetus. RESULTS: The average area under the ROC curve (AUC) values are 0.996 over all families for paternally inherited variants, 0.81 for the challenging maternally inherited variants, 0.86 for homozygous biallelic variants and 0.95 for compound heterozygous variants. Discriminative AUCs were achieved even in families with a low FF. We further investigate the performance of our method in correctly predicting SNVs in coding regions of clinically relevant genes and demonstrate significantly improved AUCs in these regions. Finally, we focus on the pathogenic variants in our cohort and show that our method correctly predicts if the fetus is unaffected or affected in all (10/10, 100%) of the families containing a pathogenic SNV. CONCLUSIONS: Overall, we demonstrate our ability to perform genome-wide NIPS for maternal and homozygous biallelic variants and showcase the utility of our method in a clinical setting.

PNPT1, MYO15A, PTPRQ, and SLC12A2-associated genetic and phenotypic heterogeneity among hearing impaired assortative mating families in Southern India.

The study was conducted between 2018 and 2020. From a cohort of 113 hearing impaired (HI), five non-DFNB12 probands identified with heterozygous CDH23 variants were subjected to exome analysis. This resolved the etiology of hearing loss (HL) in four South Indian assortative mating families. Six variants, including three novel ones, were identified in four genes: PNPT1 p.(Ala46Gly) and p.(Asn540Ser), MYO15A p.(Leu1485Pro) and p.(Tyr1891Ter), PTPRQ p.(Gln1336Ter), and SLC12A2 p.(Pro988Ser). Compound heterozygous PNPT1 variants were associated with DFNB70 causing prelingual profound sensorineural hearing loss (SNHL), vestibular dysfunction, and unilateral progressive vision loss in one family. In the second family, MYO15A variants in the myosin motor domain, including a novel variant, causing DFNB3, were found to be associated with prelingual profound SNHL. A novel PTPRQ variant was associated with postlingual progressive sensorineural/mixed HL and vestibular dysfunction in the third family with DFNB84A. In the fourth family, the SLC12A2 novel variant was found to segregate with severe-to-profound HL causing DFNA78, across three generations. Our results suggest a high level of allelic, genotypic, and phenotypic heterogeneity of HL in these families. This study is the first to report the association of PNPT1, PTPRQ, and SLC12A2 variants with HL in the Indian population.

Molecular epidemiology of non-syndromic autosomal recessive congenital ichthyosis in a Middle-Eastern population.

Autosomal recessive congenital ichthyosis (ARCI) is a rare and heterogeneous skin cornification disorder presenting with generalized scaling and varying degrees of erythema. Clinical manifestations range from lamellar ichthyosis (LI), congenital ichthyosiform erythroderma (CIE) through the most severe form of ARCI, Harlequin ichthyosis (HI). We used homozygosity mapping, whole-exome and direct sequencing to delineate the relative distribution of pathogenic variants as well as identify genotype-phenotype correlations in a cohort of 62 Middle Eastern families with ARCI of various ethnic backgrounds. Pathogenic variants were identified in most ARCI-associated genes including TGM1 (21%), CYP4F22 (18%), ALOX12B (14%), ABCA12 (10%), ALOXE3 (6%), NIPAL4 (5%), PNPLA1 (3%), LIPN (2%) and SDR9C7 (2%). In 19% of cases, no mutation was identified. Our cohort revealed a higher prevalence of CYP4F22 and ABCA12 pathogenic variants and a lower prevalence of TGM1 and NIPAL4 variants, as compared to data obtained in other regions of the world. Most variants (89%) in ALOX12B were associated with CIE and were the most common cause of ARCI among patients of Muslim origin (26%). Palmoplantar keratoderma associated with fissures was exclusively a result of pathogenic variants in TGM1. To our knowledge, this is the largest cohort study of ARCI in the Middle-Eastern population reported to date. Our data demonstrate the importance of population-tailored mutation screening strategies and shed light upon specific genotype-phenotype correlations.

Genome-Wide Noninvasive Prenatal Diagnosis of SNPs and Indels.

Noninvasive prenatal diagnosis (NIPD) is an emerging field, that enables testing for diseases in the fetus with no risk to the pregnancy, compared to invasive methods (e.g., amniocentesis). The procedure is based on the presence of fetal DNA within the mother's plasma cell-free DNA (cfDNA). Today, NIPD is performed for chromosomal abnormalities (e.g., Down syndrome) and some large deletions/duplications. It is also available for point mutations but is limited for one mutation or up to several genes simultaneously. Genome-wide detection of fetal point mutations was presented in a few studies, and the first software tool for this task, Hoobari, has recently become available. Here we describe the necessary steps in genome-wide noninvasive fetal genotyping, including examples using the Hoobari software. We discuss the various materials, software, computational infrastructure, and samples required for this analysis. Genome-wide analysis of point mutations in the fetus is not widely studied, albeit much space for algorithmic improvements exists. Here we suggest practical solutions for challenges along the process. Our work assists bioinformaticians in accessing NIPD data analysis and can eventually be utilized for other cfDNA-related fields.

Genome-Wide Noninvasive Prenatal Diagnosis of De Novo Mutations.

Noninvasive prenatal diagnosis (NIPD) has become a common, safe, and effective procedure for detection of inherited diseases early in pregnancy. It is based on the analysis of fetal cell-free DNA (cffDNA) derived from the placenta, circulating in the maternal plasma. De novo mutations, although rare, cause a considerable number of dominant genetic disorders. Due to the sparse representation of fetal-derived sequences in the blood, the challenge of detecting low frequency fetal de novo mutations becomes preponderant. Hence, this detection type requires deep genome-wide sequencing of cffDNA from maternal plasma and a unique analysis approach. Here we suggest and discuss a method for identifying de novo mutations based on whole genome sequencing (WGS) of cell-free DNA (cfDNA) from maternal plasma samples. Our method consists of an augmented pipeline for analysis of de novo mutation candidates. It begins with an enhanced noninvasive fetal variant calling step, followed by a candidate de novo mutation filtration, and then finally, a supervised machine learning approach is utilized for reduction of false positive rates. Overall, this study provides a basis for genome-wide de novo mutation analysis in NIPD procedures, which could be used in any procedure where rare de novo mutations should be carefully picked out of a sea of data.

Improved noninvasive fetal variant calling using standardized benchmarking approaches.

The technology of noninvasive prenatal testing (NIPT) enables risk-free detection of genetic conditions in the fetus, by analysis of cell-free DNA (cfDNA) in maternal blood. For chromosomal abnormalities, NIPT often effectively replaces invasive tests (e.g. amniocentesis), although it is considered as screening rather than diagnostics. Most recently, the NIPT has been applied to genome-wide, comprehensive genotyping of the fetus using cfDNA, i.e. identifying all its genetic variants and mutations. Previously, we suggested that NIPD should be treated as a special case of variant calling, and presented Hoobari, the first software tool for noninvasive fetal variant calling. Using a unique pipeline, we were able to comprehensively decipher the inheritance of SNPs and indels. A few caveats still exist in this pipeline. Performance was lower for indels and biparental loci (i.e. where both parents carry the same mutation), and performance was not uniform across the genome. Here we utilized standardized methods for benchmarking of variant calling pipelines and applied them to noninvasive fetal variant calling. By using the best performing pipeline and by focusing on coding regions, we showed that noninvasive fetal genotyping greatly improves performance, particularly in indels and biparental loci. These results emphasize the importance of using widely accepted concepts to describe the challenge of genome-wide NIPT of point mutations; and demonstrate a benchmarking process for the first time in this field. This study brings genome-wide and complete NIPD closer to the clinic; while potentially alleviating uncertainty and anxiety during pregnancy, and promoting informed choices among families and physicians.

Genome-wide noninvasive prenatal diagnosis of monogenic disorders: Current and future trends.

Noninvasive prenatal diagnosis (NIPD) is a risk-free alternative to invasive methods for prenatal diagnosis, e.g. amniocentesis. NIPD is based on the presence of fetal DNA within the mother's plasma cell-free DNA (cfDNA). Though currently available for various monogenic diseases through detection of point mutations, NIPD is limited to detecting one mutation or up to several genes simultaneously. Noninvasive prenatal whole exome/genome sequencing (WES/WGS) has demonstrated genome-wide detection of fetal point mutations in a few studies. However, Genome-wide NIPD of monogenic disorders currently has several challenges and limitations, mainly due to the small amounts of cfDNA and fetal-derived fragments, and the deep coverage required. Several approaches have been suggested for addressing these hurdles, based on various technologies and algorithms. The first relevant software tool, Hoobari, recently became available. Here we review the approaches proposed and the paths required to make genome-wide monogenic NIPD widely available in the clinic.

Loss-of-Function Variants in SERPINA12 Underlie Autosomal Recessive Palmoplantar Keratoderma.

Inherited palmoplantar keratodermas refer to a large and heterogeneous group of conditions resulting from abnormal epidermal differentiation and featuring thickening of the skin of the palms and soles. Here, we aimed at delineating the genetic basis of an autosomal recessive form of palmoplantar keratodermas manifesting with erythematous hyperkeratotic plaques over the palms and soles, extending to non-palmoplantar areas. Whole-exome sequencing in affected individuals revealed homozygous nonsense variants in the SERPINA12 gene. SERPINA12 encodes the visceral adipose tissue-derived serpin A12, a serine protease inhibitor. The pathogenic variants were found to result in reduced visceral adipose tissue-derived serpin A12 expression in patients' skin biopsies in comparison to healthy controls. In addition, SERPINA12 downregulation in three-dimensional skin equivalents was associated with marked epidermal acanthosis and hyperkeratosis, replicating the human phenotype. Moreover, decreased SERPINA12 expression resulted in reduced visceral adipose tissue-derived serpin A12-mediated inhibition of kallikrein 7 activity as well as decreased levels of desmoglein-1 and corneodesmosin, two known kallikrein 7 substrates, which are required for normal epidermal differentiation. The present data, taken collectively, demarcate a unique type of autosomal recessive palmoplantar keratodermas, attribute to visceral adipose tissue-derived serpin A12 a role in skin biology, and emphasize the importance of mechanisms regulating proteolytic activity for normal epidermal differentiation.

[LINEAR MYCOSIS FUNGOIDES IN A BLASCHKOID DISTRIBUTION].

INTRODUCTION: Mycosis fungoides (MF) is the most common type of primary cutaneous T cell lymphoma. Many clinicopathological variants of MF have been described in the literature, though only a few presented in a segmental pattern. There are several unique patterns of distribution of skin diseases, one of which is the Blaschko Lines. Congenital skin diseases develop in a Blaschkoid pattern due to mosaicism. In contrast, according to Happle, the development of acquired skin diseases in a similar pattern is explained by superimposed segmental manifestation - a process which involves mosaicism overlapping a preexisting congenital mutation. The theories by which previous case reports explained the segmental appearance of MF did not cover the molecular basis for their development. We report a case of a patient who presented with MF in a unique segmental distribution consistent with the Blaschko lines. The patient was found to have an acquired mosaic mutation in GNAS gene exclusively in the involved skin which represents a superimposed segmental manifestation according to Happle's theory. This case demonstrates the hidden potential of these rare cases which allows a better understanding of the pathogenesis by which acquired diseases develop. This is a basis for further research that could help identify new therapeutic targets for MF and other diseases that share its genetic etiology.

A novel intronic mutation of PDE6B is a major cause of autosomal recessive retinitis pigmentosa among Caucasus Jews.

Purpose: To identify the genetic basis for retinitis pigmentosa (RP) in a cohort of Jewish patients from Caucasia. Methods: Patients underwent a detailed ophthalmic evaluation, including funduscopic examination, visual field testing, optical coherence tomography (OCT), and electrophysiological tests, electroretinography (ERG) and visual evoked potentials (VEP). Genetic analysis was performed with a combination of whole exome sequencing (WES) and Sanger sequencing. Bioinformatic analysis of the WES results was performed via a customized pipeline. Pathogenicity of the identified intronic variant was evaluated in silico using the web tool Human Splicing Finder, and in vitro, using a minigene-based splicing assay. Linkage disequilibrium (LD) analysis was used to demonstrate a founder effect, and the decay of LD over generations around the mutation in Caucasus Jewish chromosomes was modeled to estimate the age of the most recent common ancestor. Results: In eight patients with RP from six unrelated families, all of Caucasus Jewish ancestry, we identified a novel homozygous intronic variant, located at position -9 of PDE6B intron 15. The c.1921-9C>G variant was predicted to generate a novel acceptor splice site, nine bases upstream of the original splice site of intron 15. In vitro splicing assay demonstrated that this novel acceptor splice site is used instead of the wild-type site, leading to an 8-bp insertion into exon 16, which is predicted to cause a frameshift. The presence of a common ancestral haplotype in mutation-bearing chromosomes was compatible with a founder effect. Conclusions: The PDE6B c.1921-9C>G intronic mutation is a founder mutation that accounts for at least 40% (6/15 families) of autosomal recessive RP among Caucasus Jews. This result is highly important for molecular diagnosis, carrier screening, and genetic counseling in this population.

Bayesian-based noninvasive prenatal diagnosis of single-gene disorders.

In the last decade, noninvasive prenatal diagnosis (NIPD) has emerged as an effective procedure for early detection of inherited diseases during pregnancy. This technique is based on using cell-free DNA (cfDNA) and fetal cfDNA (cffDNA) in maternal blood, and hence, has minimal risk for the mother and fetus compared with invasive techniques. NIPD is currently used for identifying chromosomal abnormalities (in some instances) and for single-gene disorders (SGDs) of paternal origin. However, for SGDs of maternal origin, sensitivity poses a challenge that limits the testing to one genetic disorder at a time. Here, we present a Bayesian method for the NIPD of monogenic diseases that is independent of the mode of inheritance and parental origin. Furthermore, we show that accounting for differences in the length distribution of fetal- and maternal-derived cfDNA fragments results in increased accuracy. Our model is the first to predict inherited insertions-deletions (indels). The method described can serve as a general framework for the NIPD of SGDs; this will facilitate easy integration of further improvements. One such improvement that is presented in the current study is a machine learning model that corrects errors based on patterns found in previously processed data. Overall, we show that next-generation sequencing (NGS) can be used for the NIPD of a wide range of monogenic diseases, simultaneously. We believe that our study will lead to the achievement of a comprehensive NIPD for monogenic diseases.

Variant PADI3 in Central Centrifugal Cicatricial Alopecia.

BACKGROUND: Central centrifugal cicatricial alopecia (CCCA) is the most common form of scarring alopecia among women of African ancestry. The disease is occasionally observed to affect women in families in a manner that suggests an autosomal dominant trait and usually manifests clinically after intense hair grooming. We sought to determine whether there exists a genetic basis of CCCA and, if so, what it is. METHODS: We used exome sequencing in a group of women with alopecia (discovery set), compared the results with those in a public repository, and applied other filtering criteria to identify candidate genes. We then performed direct sequencing to identify disease-associated DNA variations and RNA sequencing, protein modeling, immunofluorescence staining, immunoblotting, and an enzymatic assay to evaluate the consequences of potential etiologic mutations. We used a replication set that consisted of women with CCCA to confirm the data obtained with the discovery set. RESULTS: In the discovery set, which included 16 patients, we identified one splice site and three heterozygous missense mutations in PADI3 in 5 patients (31%). (The approximate prevalence of the disease is up to 5.6%.) PADI3 encodes peptidyl arginine deiminase, type III (PADI3), an enzyme that post-translationally modifies other proteins that are essential to hair-shaft formation. All three CCCA-associated missense mutations in PADI3 affect highly conserved residues and are predicted to be pathogenic; protein modeling suggests that they result in protein misfolding. These mutations were found to result in reduced PADI3 expression, abnormal intracellular localization of the protein, and decreased enzymatic activity - findings that support their pathogenicity. Immunofluorescence staining showed decreased expression of PADI3 in biopsy samples of scalp skin obtained from patients with CCCA. We then directly sequenced PADI3 in an additional 42 patients (replication set) and observed genetic variants in 9 of them. A post hoc analysis of the combined data sets showed that the prevalence of PADI3 mutation was higher among patients with CCCA than in a control cohort of women of African ancestry (P = 0.002 by the chi-square test; P = 0.006 by Fisher's exact test; and after adjustment for relatedness of persons, P = 0.03 and P = 0.04, respectively). CONCLUSIONS: Mutations in PADI3, which encodes a protein that is essential to proper hair-shaft formation, were associated with CCCA. (Funded by the Ram Family Foundation and others.).

Filaggrin 2 Deficiency Results in Abnormal Cell-Cell Adhesion in the Cornified Cell Layers and Causes Peeling Skin Syndrome Type A.

Peeling skin syndromes form a large and heterogeneous group of inherited disorders characterized by superficial detachment of the epidermal cornified cell layers, often associated with inflammatory features. Here we report on a consanguineous family featuring noninflammatory peeling of the skin exacerbated by exposure to heat and mechanical stress. Whole exome sequencing revealed a homozygous nonsense mutation in FLG2, encoding filaggrin 2, which cosegregated with the disease phenotype in the family. The mutation was found to result in decreased FLG2 RNA levels as well as almost total absence of filaggrin 2 in the patient epidermis. Filaggrin 2 was found to be expressed throughout the cornified cell layers and to colocalize with corneodesmosin that plays a crucial role in maintaining cell-cell adhesion in this region of the epidermis. The absence of filaggrin 2 in the patient skin was associated with markedly decreased corneodesmosin expression, which may contribute to the peeling phenotype displayed by the patients. Accordingly, using the dispase dissociation assay, we showed that FLG2 downregulation interferes with keratinocyte cell-cell adhesion. Of particular interest, this effect was aggravated by temperature elevation, consistent with the clinical phenotype. Restoration of corneodesmosin levels by ectopic expression rescued cell-cell adhesion. Taken together, the present data suggest that filaggrin 2 is essential for normal cell-cell adhesion in the cornified cell layers.