Non-invasive prenatal testing for autosomal recessive disorders: A new promising approach

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.

Loss-of-function variant in the LRR domain of SLITRK2 implicated in a neurodevelopmental disorder

Background: Neurodevelopmental disorders are characterized by different combinations of intellectual disability (ID), communication and social skills deficits, and delays in achieving motor or language milestones. SLITRK2 is a postsynaptic cell-adhesion molecule that promotes neurite outgrowth and excitatory synapse development. Methods and Results: In the present study, we investigated a single patient segregating Neurodevelopmental disorder. SLITRK2 associated significant neuropsychological issues inherited in a rare X-linked fashion have recently been reported. Whole-exome sequencing and data analysis revealed a novel nonsense variant [c.789T>A; p.(Cys263*); NM_032539.5; NP_115928.1] in exon 5 of the SLITRK2 gene (MIM# 300561). Three-dimensional protein modeling revealed substantial changes in the mutated SLITRK2 protein, which might lead to nonsense-medicated decay. Conclusion: This study confirms the role of SLITRK2 in neuronal development and highlights the importance of including the SLITRK2 gene in the screening of individuals presenting neurodevelopmental disorders.

Novel compound heterozygous MYO15A splicing variants in autosomal recessive non-syndromic hearing loss

BACKGROUND

Hereditary hearing loss is a highly heterogeneous disorder. This study aimed to identify the genetic cause of a Chinese family with autosomal recessive non-syndromic sensorineural hearing loss (ARNSHL).

METHODS

Clinical information and peripheral blood samples were collected from the proband and its parents. Two-step high-throughput next-generation sequencing on the Ion Torrent platform was applied to detect variants as follows. First, long-range PCR was performed to amplify all the regions of the GJB2, GJB3, SLC26A4, and MT-RNR1 genes, followed by next-generation sequencing. If no candidate pathogenetic variants were found, the targeted exon sequencing with AmpliSeq technology was employed to examine another 64 deafness-associated genes. Sanger sequencing was used to identify variants and the lineage co-segregation. The splicing of the MYO15A gene was assessed by in silico bioinformatics prediction and minigene assays.

RESULTS

Two candidate MYO15A gene (OMIM, #602,666) heterozygous splicing variants, NG_011634.2 (NM_016239.3): c.6177 + 1G > T and c.9690 + 1G > A, were identified in the proband, and these two variants were both annotated as pathogenic according to the American College of Medical Genetics and Genomics (ACMG) guidelines. Further bioinformatic analysis predicted that the c.6177 + 1G > T variant might cause exon skipping and that the c.9690 + 1G > A variant might activate a cryptic splicing donor site in the downstream intronic region. An in vitro minigene assay confirmed the above predictions.

CONCLUSIONS

We identified a compound heterozygous splicing variant in the MYO15A gene in a Han Chinese family with ARNSHL. Our results broaden the spectrum of MYO15A variants, potentially benefiting the early diagnosis, prevention, and treatment of the disease.

Molecular insight into CREBBP and TANGO2 variants causing intellectual disability

BACKGROUND

Intellectual disability (ID) can be associated with different syndromes such as Rubinstein-Taybi syndrome (RSTS) and can also be related to conditions such as metabolic encephalomyopathic crises, recurrent,with rhabdomyolysis, cardiac arrhythmias and neurodegeneration. Rare congenital RSTS1 (OMIM 180849) is characterized by mental and growth retardation, significant and duplicated distal phalanges of thumbs and halluces, facial dysmorphisms, and an elevated risk of malignancies. Microdeletions and point mutations in the CREB-binding protein (CREBBP) gene, located at 16p13.3, have been reported to cause RSTS. By contrast, TANGO2-related metabolic encephalopathy and arrhythmia (TRMEA) is a rare metabolic condition that causes repeated metabolic crises, hypoglycemia, lactic acidosis, rhabdomyolysis, arrhythmias and encephalopathy with cognitive decline. Clinicians need more clinical and genetic evidence to detect and comprehend the phenotypic spectrum of this disorder.

METHODS

Exome sequencing was used to identify the disease-causing variants in two affected families A and B from District Kohat and District Karak, Khyber Pakhtunkhwa. Affected individuals from both families presented symptoms of ID, developmental delay and behavioral abnormalities. The validation and co-segregation analysis of the filtered variant was carried out using Sanger sequencing.

RESULTS

In the present study, two families (A and B) exhibiting various forms of IDs were enrolled. In Family A, exome sequencing revealed a novel missense variant (NM 004380.3: c.4571A>G; NP_004371.2: p.Lys1524Arg) in the CREBBP gene, whereas, in Family B, a splice site variant (NM 152906.7: c.605 + 1G>A) in the TANGO2 gene was identified. Sanger sequencing of both variants confirmed their segregation with ID in both families. The in silico tools verified the aberrant changes in the CREBBP protein structure. Wild-type and mutant CREBBP protein structures were superimposed and conformational changes were observed likely altering the protein function.

CONCLUSIONS

RSTS and TRMEA are exceedingly rare disorders for which specific clinical characteristics have been clearly established, but more investigations are underway and required. Multicenter studies are needed to increase our understanding of the clinical phenotypes, mainly showing the genotype-phenotype associations.

Identification of novel homozygous variants in FOXE3 and AP4M1 underlying congenital syndromic anophthalmia and microphthalmia

BACKGROUND

Anophthalmia and microphthalmia are severe developmental ocular disorders that affect the size of the ocular globe and can be unilateral or bilateral. The disease is found in syndromic as well as non-syndromic forms. It is genetically caused by chromosomal aberrations, copy number variations and single gene mutations, along with non-genetic factors such as viral infections, deficiency of vitamin A and an exposure to alcohol or drugs during pregnancy. To date, more than 30 genes having different modes of inheritance patterns are identified as causing anophthalmia and microphthalmia.

METHODS

In the present study, a clinical and genetic analysis was performed of six patients with anophthalmia and microphthalmia and/or additional phenotypes of intellectual disability, developmental delay and cerebral palsy from a large consanguineous Pakistani family. Whole exome sequencing followed by data analysis for variants prioritization and validation through Sanger sequencing was performed to identify the disease causing variant(s). American College of Medical Genetics and Genomics (ACMG) guidelines were applied to classify clinical interpretation of the prioritized variants.

RESULTS

Clinical investigations revealed that the affected individuals are afflicted with anophthalmia. Three of the patients showed additional phenotype of intellectual disability, developmental delays and other neurological symptoms. Whole exome sequencing of the DNA samples of the affected members in the family identified a novel homozygous stop gain mutation (NM_012186: c.106G>T: p.Glu36*) in Forkhead Box E3 (FOXE3) gene shared by all affected individuals. Moreover, patients segregating additional phenotypes of spastic paraplegia, intellectual disability, hearing loss and microcephaly showed an additional homozygous sequence variant (NM_004722: c.953G>A: p.Arg318Gln) in AP4M1. Sanger sequencing validated the correct segregation of the identified variants in the affected family. ACMG guidelines predicted the variants to be pathogenic.

CONCLUSIONS

We have investigated first case of syndromic anophthalmia caused by variants in the FOXE3 and AP4M1. The present findings are helpful for understanding pathological role of the mutations of the genes in syndromic forms of anophthalmia. Furthermore, the study signifies searching for the identification of second variant in families with patients exhibiting variable phenotypes. In addition, the findings will help clinical geneticists, genetic counselors and the affected family with respect to prenatal testing, family planning and genetic counseling.

A variant in sperm-specific glycolytic enzyme enolase 4 (ENO4) causes human male infertility

BACKGROUND

Although defects in sperm morphology and physiology lead to male infertility, in many instances, the exact disruption of molecular pathways in a given patient is often unknown. The glycolytic pathway is an essential process to supply energy in sperm cell motility. Enolase 4 (ENO4) is crucial for the glycolytic process, which provides the energy for sperm cells in motility. ENO4 is located in the sperm principal piece and is essential for the motility and organization of the sperm flagellum. In the present study, we characterized a family with asthenozoospermia and abnormal sperm morphology as a result of a variant in the enolase 4 (ENO4) gene.

METHODS

Computer-assisted semen analysis, papanicolaou smear staining and scanning electron microscopy were used to examine sperm motility and morphology for semen analysis in patients. For genetic analysis, whole-exome sequencing followed by Sanger sequencing was performed.

RESULTS

Two brothers in a consanguineous family were being clinically investigated for sperm motility and morphology issues. Genetic analysis by whole-exome sequencing revealed a homozygous variant [c.293A>G, p.(Lys98Arg)] in the ENO4 gene that segregated with infertility in the family, shared by affected but not controls.

CONCLUSIONS

In view of the association of asthenozoospermia and abnormal sperm morphology in Eno4 knockout mice, we consider this to be the first report describing the involvement of ENO4 gene in human male infertility. We also explore the possible involvement of another variant in explaining other phenotypic features in this family.

Clinical features and genetic analysis of a case series of skeletal ciliopathies in a prenatal setting

BACKGROUND

Short-rib polydactyly syndrome (SRPS) refers to a group of lethal skeletal dysplasias that can be difficult to differentiate between subtypes or from other non-lethal skeletal dysplasias such as Ellis-van Creveld syndrome and Jeune syndrome in a prenatal setting. We report the ultrasound and genetic findings of four unrelated fetuses with skeletal dysplasias.

METHODS

Systemic prenatal ultrasound examination was performed in the second or third trimester. Genetic tests including GTG-banding, single nucleotide polymorphism (SNP) array and exome sequencing were performed with amniocytes or aborted fetal tissues.

RESULTS

The major and common ultrasound anomalies for the four unrelated fetuses included short long bones of the limbs and narrow thorax. No chromosomal abnormalities and pathogenic copy number variations were detected. Exome sequencing revealed three novel variants in the DYNC2H1 gene, namely NM_001080463.2:c.6809G > A p.(Arg2270Gln), NM_001080463.2:3133C > T p.(Gln1045Ter), and NM_001080463.2:c.337C > T p.(Arg113Trp); one novel variant in the IFT172 gene, NM_015662.3:4540-5 T > A; and one novel variant in the WDR19 gene, NM_025132.4:c.2596G > C p.(Gly866Arg). The genotypes of DYNC2H1, IFT172 and WDR19 and the phenotypes of the fetuses give hints for the diagnosis of short-rib thoracic dysplasia (SRTD) with or without polydactyly 3, 10, and 5, respectively.

CONCLUSION

Our findings expand the mutation spectrum of DYNC2H1, IFT172 and WDR19 associated with skeletal ciliopathies, and provide useful information for prenatal diagnosis and genetic counseling on rare skeletal disorders.

Truncated DNM1 variant underlines developmental delay and epileptic encephalopathy

Background

Developmental and epileptic encephalopathies (DEEs) signify a group of heterogeneous neurodevelopmental disorder associated with early-onset seizures accompanied by developmental delay, hypotonia, mild to severe intellectual disability, and developmental regression. Variants in the DNM1 gene have been associated with autosomal dominant DEE type 31A and autosomal recessive DEE type 31B.

Methods

In the current study, a consanguineous Pakistani family consisting of a proband (IV-2) was clinically evaluated and genetically analyzed manifesting in severe neurodevelopmental phenotypes. WES followed by Sanger sequencing was performed to identify the disease-causing variant. Furthermore, 3D protein modeling and dynamic simulation of wild-type and mutant proteins along with reverse transcriptase (RT)-based mRNA expression were checked using standard methods.

Results

Data analysis of WES revealed a novel homozygous non-sense variant (c.1402G>T; p. Glu468*) in exon 11 of the DNM1 gene that was predicted as pathogenic class I. Variants in the DNM1 gene have been associated with DEE types 31A and B. Different bioinformatics prediction tools and American College of Medical Genetics guidelines were used to verify the identified variant. Sanger sequencing was used to validate the disease-causing variant. Our approach validated the pathogenesis of the variant as a cause of heterogeneous neurodevelopmental disorders. In addition, 3D protein modeling showed that the mutant protein would lose most of the amino acids and might not perform the proper function if the surveillance non-sense-mediated decay mechanism was skipped. Molecular dynamics analysis showed varied trajectories of wild-type and mutant DNM1 proteins in terms of root mean square deviation, root mean square fluctuation and radius of gyration. Similarly, RT-qPCR revealed a substantial reduction of the DNM1 gene in the index patient.

Conclusion

Our finding further confirms the association of homozygous, loss-of-function variants in DNM1 associated with DEE type 31B. The study expands the genotypic and phenotypic spectrum of pathogenic DNM1 variants related to DNM1-associated pathogenesis.

Overlapping neurological phenotypes in two extended consanguineous families with novel variants in the CNTNAP1 and ADGRG1 genes

BACKGROUND

Population diversity is important and rare disease isolates can frequently reveal novel homozygous or biallelic mutations that lead to expanded clinical heterogeneity, with diverse clinical presentations.

METHODS

The present study describes two consanguineous families with a total of seven affected individuals suffering from a clinically similar severe syndromic neurological disorder, with abnormal development and central nervous system (CNS) and peripheral nervous system (PNS) abnormalities. Whole exome sequencing (WES) and Sanger sequencing followed by 3D protein modeling was performed to identify the disease-causing gene. RNA was extracted from the fresh blood of both families affected and healthy individuals.

RESULTS

The families were clinically assessed in the field in different regions of Khyber Pakhtunkhwa. Magnetic resonance imagining was obtained in the probands and blood was collected for DNA extraction and WES was performed. Sanger sequencing confirmed a homozygous, likely pathogenic mutation (GRCh38: chr17:42684199G>C; (NM_003632.3): c.333G>C);(NP_003623.1): p.(Trp111Cys) in the CNTNAP1 gene in family A, previously associated with Congenital Hypo myelinating Neuropathy 3 (CHN3; OMIM # 618186) and a novel nonsense variant in family B, (GRCh38: chr16: 57654086C>T; NC_000016.10 (NM_001370440.1): c.721C>T); (NP_001357369.1): p.(Gln241Ter) in the ADGRG1 gene previously associated with bilateral frontoparietal polymicrogyria (OMIM # 606854); both families have extended CNS and PNS clinical manifestations. In addition, 3D protein modeling was performed for the missense variant, p.(Trp111Cys), identified in the CNTNAP1, suggesting extensive secondary structure changes that might lead to improper function or downstream signaling. No RNA expression was observed in both families affected and healthy individuals hence showing that these genes are not expressed in blood.

CONCLUSIONS

In the present study, two novel biallelic variants in the CNTNAP1 and ADGRG1 genes in two different consanguineous families with a clinical overlap in the phenotype were identified. Thus, the clinical and mutation spectrum is expanded to provide further evidence that CNTNAP1 and ADGRG1 are very important for widespread neurological development.

Two novel homozygous variants of ATP6V0A2 and ALDH18A1 lead to autosomal recessive cutis laxa type 2 and 3 in two Pakistani families

BACKGROUND

Autosomal recessive cutis laxa type 2A (ARCL2A; OMIM: 219200) is characterized by neurovegetative, developmental and progeroid elastic skin anomalies. It is caused by biallelic variation in ATPase, H+ transporting V0 subunit A2 (ATP6V0A2; OMIM: 611716) located on chromosome 12q24.31. Autosomal recessive cutis laxa type 3A (ARCL3A; OMIM: 219150) is another subclinical type characterized by short stature, ophthalmological abnormalities and a progeria-like appearance. The ARCL3A is caused by loss of function alterations in the aldehyde dehydrogenase 18 family member A1 (ALDH18A1; OMIM: 138250) gene located at chromosome 10q24.1.

METHODS

Whole-exome sequencing (WES), and Sanger sequencing were performed for molecular diagnosis. 3D protein modeling was performed to investigate the deleterious effect of the variant on protein structure.

RESULTS

In this study, clinical and molecular diagnosis were performed for two families, ED-01 and DWF-41, which displayed hallmark features of ARCL2A and ARCL3A, respectively. Three affected individuals in the ED-01 family (IV-4, IV-5 and V-3) displayed sagging loose skin, down-slanting palpebral fissures, excessive wrinkles on the abdomen, hands and feet, and prominent veins on the trunk. Meanwhile the affected individuals in the DWF-41 family (V-2 and V-3) had progeroid skin, short stature, dysmorphology, low muscle tone, epilepsy, lordosis, scoliosis, delayed puberty and internal genitalia. WES in the index patient (ED-01: IV-4) identified a novel homozygous deletion (NM_012463.3: c.1977_1980del; p.[Val660LeufsTer23]) in exon 16 of the ATP6V0A2 while in DWF-41 a novel homozygous missense variant (NM_001323413.1:c.1867G>A; p.[Asp623Asn]) in exon 15 of the ALDH18A1 was identified. Sanger validation in all available family members confirmed the autosomal recessive modes of inheritances in each family. Three dimensional in-silico protein modeling suggested deleterious impact of the identified variants. Furthermore, these variants were assigned class 1 or "pathogenic" as per guidelines of American College of Medical Genetics 2015. Screening of ethnically matched healthy controls (n = 200 chromosomes), excluded the presence of these variations in general population.

CONCLUSIONS

To the best of our knowledge, this is the first report of ATP6V0A2 and ALDH18A1 variations in the Pakhtun ethnicity of Pakistani population. The study confirms that WES can be used as a first-line diagnostic test in patients with cutis laxa, and provides basis for population screening and premarital testing to reduce the diseases burden in future generations.

Variants in FREM1 and trisomy 18 identified in a neonatal progeria patient

BACKGROUND

Neonatal progeroid disorders are rare disorders with clinical features including low body mass index, proptosis, aged and dysmorphic facial features at the time of birth, prominent veins, sparse scalp hairs, and severe growth retardation. Very few cases have been identified with an unknown genetic cause. Here, we report clinical and genetic findings of a proband with hallmark features of neonatal progeria.

METHODS

Microarray comparative genomic hybridization, whole exome sequencing (WES) and Sanger sequencing were performed using standard methods.

RESULTS

Array combined genome hybridization data revealed trisomy 18 in the proband (II-1), and WES data identified novel compound heterozygous variants (c.247 C > T; p.H83Y and c.14769868InsA) in the FREM1 gene.

CONCLUSION

We report a novel complex case of neonatal progeria with atrial septal defects, trisomy 18 without typical features of Edward syndrome. The phenotype of the patient was more consistent with neonatal progeria, thus we speculate it to be caused by the FREM1 variants.

Genetic and molecular dynamics analysis of two variants of the COL4A5 gene causing Alport syndrome

BACKGROUND

Alport syndrome (AS; OMIM#308,940) is a hereditary kidney disease that progresses over time and is distinguished by hearing loss and ocular irregularities. The syndrome has three subtypes, namely X-linked (XL; OMIM#301,050), autosomal recessive (AR; OMIM#203,780), and autosomal dominant (AD; OMIM#104,200), which are categorized based on their respective modes of inheritance. XLAS is attributed to a pathogenic variant in the COL4A5 (OMIM*303,630) gene, which encodes the α5(IV) chain of type IV collagen (Col-IV). In contrast, ADAS and ARAS are the result of variants in the COL4A3 (OMIM*120,070) and COL4A4 (OMIM*120,131) genes, which encode the α3(IV) and α4(IV) chains of Col-IV, respectively. Typically, the diagnosis of AS necessitates hereditary or pathological assessments. The determination of splicing variants as pathogenic or non-pathogenic based on gene sequencing outcomes is challenging.

METHODS

In this study, we conducted exome sequencing and Sanger sequencing on two unrelated Chinese patients with AS. We identified a deletion variant c.4414delG in the COL4A5 gene and a splicing variant c.4298-20T > A in the same gene. In order to ascertain the impact of c.4298-20T > A on the synthesis of COL4A5 mRNA, we performed experiments involving minigene splicing. Additionally, we predicted the ability of these two variants to affect triple helix formation of α345(IV) using molecular dynamics methods.

RESULTS

The c.4414delG deletion variant caused a change in the genetic code of the COL4A5 gene. Specifically, it caused a shift in codon 1472 from encoding aspartate to encoding methionine. This shift resulted in a change of 75 amino acids in the protein sequence, ultimately leading to an early stop codon. This premature stop codon caused the production of a truncated α5(IV) chain with a predicted protein effect of p.D1472Mfs. The mRNA of the COL4A5 gene experienced intron 46 retention due to the splicing variant c.4298-20T > A, leading to the inclusion of six additional amino acids between amino acids 1432 and 1433 of the α5(IV) chain. This variant is predicted to have a protein effect of p.(P1432_G1433insDYFVEI). The impact of two variants, c.4414delG and c.4298-20T > A, on the aggregation region for α3(IV), α4(IV), and α5(IV) trimerisation were studied using molecular dynamics simulations. Results showed that the deletion variant c.4414delG had a significantly stronger disruption on NC1, compared to the splicing variant c.4298-20T > A. This difference in impact is consistent with the varying clinical phenotypes observed in the two patients. Based on the American College of Medical Genetics and Genomics (ACMG) classification criteria and guidelines for genetic variants, the deletion variant c.4414delG was rated as pathogenic while the splicing variant c.4298-20T > A was rated as likely-pathogenic.

CONCLUSION

Our study has identified two novel pathogenic loci, the deletion variant c.4414delG and the splicing variant c.4298-20T > A, associated with XLAS. This finding expands the genetic spectrum of XLAS. We suggest that molecular dynamics can effectively model the effect of genetic variation on α345(IV) trimerization, which may offer valuable insights into the mechanisms of XLAS pathogenesis.

Whole exome sequencing identified five novel variants in CNTN2, CARS2, ARSA, and CLCN4 leading to epilepsy in consanguineous families

Introduction: Epilepsy is a group of neurological disorders characterized by recurring seizures and fits. The Epilepsy genes can be classified into four distinct groups, based on involvement of these genes in different pathways leading to Epilepsy as a phenotype. Genetically the disease has been associated with various pathways, leading to pure epilepsy-related disorders caused by CNTN2 variations, or involving physical or systemic issues along with epilepsy caused by CARS2 and ARSA, or developed by genes that are putatively involved in epilepsy lead by CLCN4 variations. Methods: In this study, five families of Pakistani origin (EP-01, EP-02, EP-04, EP-09, and EP-11) were included for molecular diagnosis. Results: Clinical presentations of these patients included neurological symptoms such as delayed development, seizures, regression, myoclonic epilepsy, progressive spastic tetraparesis, vision and hearing impairment, speech problems, muscle fibrillation, tremors, and cognitive decline. Whole exome sequencing in index patients and Sanger sequencing in all available individuals in each family identified four novel homozygous variants in genes CARS2: c.655G>A p.Ala219Thr (EP-01), ARSA: c.338T>C: p.Leu113Pro (EP-02), c.938G>T p.Arg313Leu (EP-11), CNTN2: c.1699G>T p.Glu567Ter (EP-04), and one novel hemizygous variant in gene CLCN4: c.2167C>T p.Arg723Trp (EP-09). Conclusion: To the best of our knowledge these variants were novel and had not been reported in familial epilepsy. These variants were absent in 200 ethnically matched healthy control chromosomes. Three dimensional protein analyses revealed drastic changes in the normal functions of the variant proteins. Furthermore, these variants were designated as "pathogenic" as per guidelines of American College of Medical Genetics 2015. Due to overlapping phenotypes, among the patients, clinical subtyping was not possible. However, whole exome sequencing successfully pinpointed the molecular diagnosis which could be helpful for better management of these patients. Therefore, we recommend that exome sequencing be performed as a first-line molecular diagnostic test in familial cases.

Biallelic Variants in Seven Different Genes Associated with Clinically Suspected Bardet-Biedl Syndrome

Bardet-Biedl syndrome (BBS) is a rare clinically and genetically heterogeneous autosomal recessive multi-systemic disorder with 22 known genes. The primary clinical and diagnostic features include six different hallmarks, such as rod-cone dystrophy, learning difficulties, renal abnormalities, male hypogonadism, post-axial polydactyly, and obesity. Here, we report nine consanguineous families and a non-consanguineous family with several affected individuals presenting typical clinical features of BBS. In the present study, 10 BBS Pakistani families were subjected to whole exome sequencing (WES), which revealed novel/recurrent gene variants, including a homozygous nonsense mutation (c.94C>T; p.Gln32Ter) in the IFT27 (NM_006860.5) gene in family A, a homozygous nonsense mutation (c.160A>T; p.Lys54Ter) in the BBIP1 (NM_001195306.1) gene in family B, a homozygous nonsense variant (c.720C>A; p.Cys240Ter) in the WDPCP (NM_015910.7) in family C, a homozygous nonsense variant (c.505A>T; p.Lys169Ter) in the LZTFL1 (NM_020347.4) in family D, pathogenic homozygous 1 bp deletion (c.775delA; p.Thr259Leufs*21) in the MKKS/BBS5 (NM_170784.3) gene in family E, a pathogenic homozygous missense variant (c.1339G>A; p.Ala447Thr) in BBS1 (NM_024649.4) in families F and G, a pathogenic homozygous donor splice site variant (c.951+1G>A; p?) in BBS1 (NM_024649.4) in family H, a pathogenic bi-allelic nonsense variant in MKKS (NM_170784.3) (c.119C>G; p.Ser40*) in family I, and homozygous pathogenic frameshift variants (c.196delA; p.Arg66Glufs*12) in BBS5 (NM_152384.3) in family J. Our findings extend the mutation and phenotypic spectrum of four different types of ciliopathies causing BBS and also support the importance of these genes in the development of multi-systemic human genetic disorders.

Mutational Spectrum of the ABCA12 Gene and Genotype–Phenotype Correlation in a Cohort of 64 Patients with Autosomal Recessive Congenital Ichthyosis

Autosomal recessive congenital ichthyosis (ARCI) is a non-syndromic congenital disorder of cornification characterized by abnormal scaling of the skin. The three major phenotypes are lamellar ichthyosis, congenital ichthyosiform erythroderma, and harlequin ichthyosis. ARCI is caused by biallelic mutations in ABCA12, ALOX12B, ALOXE3, CERS3, CYP4F22, NIPAL4, PNPLA1, SDR9C7, SULT2B1, and TGM1. The most severe form of ARCI, harlequin ichthyosis, is caused by mutations in ABCA12. Mutations in this gene can also lead to congenital ichthyosiform erythroderma or lamellar ichthyosis. We present a large cohort of 64 patients affected with ARCI carrying biallelic mutations in ABCA12. Our study comprises 34 novel mutations in ABCA12, expanding the mutational spectrum of ABCA12-associated ARCI up to 217 mutations. Within these we found the possible mutational hotspots c.4541G>A, p.(Arg1514His) and c.4139A>G, p.(Asn1380Ser). A correlation of the phenotype with the effect of the genetic mutation on protein function is demonstrated. Loss-of-function mutations on both alleles generally result in harlequin ichthyosis, whereas biallelic missense mutations mainly lead to CIE or LI.

Biallelic mutations in FLG, TGM1, and STS genes segregated with different types of ichthyoses in eight families of Pakistani origin

BACKGROUND

Congenital ichthyosis is a diverse group of keratinization disorders associated with generalized scaling of skin of varying severity. The non-syndromic forms of congenital ichthyosis are further grouped into common ichthyosis (ichthyosis vulgaris and X-linked ichthyosis), autosomal recessive congenital ichthyosis, and keratopathic ichthyosis.

OBJECTIVE

To identify sequence variants involved in different forms of hereditary ichthyoses.

METHODS

We studied eight families with different types of ichthyosis including four families with autosomal recessive congenital ichthyosis and four families with common ichthyosis. Whole exome sequencing and PCR based genotyping was carried out to find out the molecular basis of disease.

RESULTS

In one family, a novel duplication sequence variant NM_002016.2:c.2767dupT; NP_002007.1:p.Ser923PhefsTer2 was identified in FLG gene; in four families a previously reported nonsense sequence variant NM_000359.3:c.232C>T; NP_002007.1:p.Arg78Ter was identified in TGM1 gene, while, in three families of X-linked recessive ichthyosis, the whole STS gene (NM_001320752.2; NP_001307681.2) regions were deleted.

STUDY LIMITATION

Gene expression studies have not been performed that would have strengthened the findings of computational analysis.

CONCLUSION

This study highlights the significance of the c.232C>T variant in the TGM1 gene as a possible founder mutation, complete STS gene deletion as reported previously in Pakistani population, while novel sequence variant in the FLG gene expands the spectrum of variations in this gene. These findings may be used for genetic counseling of the studied families.

Spondyloocular Syndrome: A Report of an Additional Family and Phenotypic Spectrum Delineation

Spondyloocular syndrome (SOS, OMIM # 605822) is a rare genetic disorder characterized by osseous and ocular manifestations, including generalized osteoporosis, multiple long bones fractures, platyspondyly, dense cataracts and retinal detachment, and dysmorphic facial features, with or without short stature, cardiopathy, hearing impairment, and intellectual disability. Biallelic mutations in the XYLT2 gene (OMIM * 608125), encoding the xylosyltransferase II, were shown to be responsible for this disease. To date, 22 cases with SOS have been described, with varying clinical presentations and a yet-to-be-established genotypic-phenotypic correlation. Two patients from a consanguineous Lebanese family that presented with SOS were included in this study. Whole exome sequencing revealed a novel homozygous nonsense mutation in XYLT2 (p.Tyr414*) in these patients. We review all previously reported cases with SOS, describe the second nonsense mutation in XYLT2, and contribute to a better delineation of the phenotypic spectrum of the disease.

Homozygous Missense Variant in the Solute Carrier Organic Anion Transporter 2A1 (SLCO2A1) Gene Underlies Isolated Nail Clubbing

BACKGROUND

Inherited isolated nail clubbing is a very rare Mendelian condition in humans, characterized by enlargement of the terminal segments of fingers and toes with thickened nails. Mutations in two genes have been reported to cause isolated nail clubbing in humans, which are the SLCO2A1 gene and the HPGD gene.

OBJECTIVES

An extended Pakistani family having two affected siblings born of unaffected consanguineous union was included in the study. Predominant isolated congenital nail clubbing (ICNC) without any other systemic abnormalities was observed, which we aimed to characterize at clinico-genetic level.

METHODS

Whole exome coupled with Sanger sequencing were employed to uncover the sequence variant as a cause of the disease. Furthermore, protein modeling was carried out to reveal the predicted possible effect of the mutation at the protein level.

RESULTS

Whole exome sequencing data analysis revealed a novel biallelic sequence variant (c.155T>A; p.Phe52Tyr) in the SLCO2A1 gene. Further, Sanger sequencing analysis validated and confirmed the segregation of the novel variant in the entire family. Subsequently, protein modeling of the wild-type and mutated SLCO2A1 revealed broad-scale change, which might compromise the proteins' secondary structure and function.

CONCLUSION

The present study adds another mutation to the SLCO2A1-related pathophysiology. The involvement of SLCO2A1 in the pathogenesis of ICNC may open exciting perceptions of this gene in nail development/morphogenesis.

Novel Variants in MPV17, PRX, GJB1, and SACS Cause Charcot-Marie-Tooth and Spastic Ataxia of Charlevoix-Saguenay Type Diseases

Charcot-Marie-Tooth disease (CMT) and autosomal recessive spastic ataxia of Charlevoix-Saguenay type (ARSACS) are large heterogeneous groups of sensory, neurological genetic disorders characterized by sensory neuropathies, muscular atrophies, abnormal sensory conduction velocities, and ataxia. CMT2EE (OMIM: 618400) is caused by mutations in MPV17 (OMIM: 137960), CMT4F (OMIM: 614895) is caused by PRX (OMIM: 605725), CMTX1 (OMIM: 302800) is caused by mutations in GJB1 (OMIM: 304040), and ARSACS (OMIM: 270550) is caused by mutations in SACS (OMIM: 604490). In this study, we enrolled four families: DG-01, BD-06, MR-01, and ICP-RD11, with 16 affected individuals, for clinical and molecular diagnoses. One patient from each family was analyzed for whole exome sequencing and Sanger sequencing was done for the rest of the family members. Affected individuals of families BD-06 and MR-01 show complete CMT phenotypes and family ICP-RD11 shows ARSACS type. Family DG-01 shows complete phenotypes for both CMT and ARSACS types. The affected individuals have walking difficulties, ataxia, distal limb weakness, axonal sensorimotor neuropathies, delayed motor development, pes cavus, and speech articulations with minor variations. The WES analysis in an indexed patient of family DG-01 identified two novel variants: c.83G>T (p.Gly28Val) in MPV17 and c.4934G>C (p.Arg1645Pro) in SACS. In family ICP-RD11, a recurrent mutation that causes ARSACS, c.262C>T (p.Arg88Ter) in SACS, was identified. Another novel variant, c.231C>A (p.Arg77Ter) in PRX, which causes CMT4F, was identified in family BD-06. In family MR-01, a hemizygous missense variant c.61G>C (p.Gly21Arg) in GJB1 was identified in the indexed patient. To the best of our knowledge, there are very few reports on MPV17, SACS, PRX, and GJB1 causing CMT and ARSACS phenotypes in the Pakistani population. Our study cohort suggests that whole exome sequencing can be a useful tool in diagnosing complex multigenic and phenotypically overlapping genetic disorders such as Charcot-Marie-Tooth disease (CMT) and spastic ataxia of Charlevoix-Saguenay type.

Case report: A novel de novo loss of function variant in the DNA-binding domain of TBX2 causes severe osteochondrodysplasia

Background: T-box family members are transcription factors characterized by highly conserved residues corresponding to the DNA-binding domain known as the T-box. TBX2 has been implicated in several developmental processes, such as coordinating cell fate, patterning, and morphogenesis of a wide range of tissues and organs, including lungs, limbs, heart, kidneys, craniofacial structures, and mammary glands. Methods: In the present study, we have clinically and genetically characterized a proband showing a severe form of chondrodysplasia with developmental delay. Whole-exome sequencing (WES), Sanger sequencing, and 3D protein modeling were performed in the present investigation. Results: Whole-exome sequencing revealed a novel nonsense variant (c.529A>T; p.Lys177*; NM_005994.4) in TBX2. 3D-TBX2 protein modeling revealed a substantial reduction of the mutated protein, which might lead to a loss of function (LOF) or nonsense-mediated decay (NMD). Conclusion: This study has not only expanded the mutation spectrum in the gene TBX2 but also facilitated the diagnosis and genetic counseling of related features in affected families.

Report of Hermansky-Pudlak Syndrome in Two Families with Novel Variants in HPS3 and HPS4 Genes

Background: Hermansky-Pudlak syndrome (HSP) was first reported in 1959 as oculocutaneous albinism with bleeding abnormalities, and now consists of 11 distinct heterogenic genetic disorders that are caused by mutations in four protein complexes: AP-3, BLOC1, BLOC2, and BLOC3. Most of the patients show albinism and a bleeding diathesis; additional features may present depending on the nature of a defective protein complex. The subtypes 3 and 4 have been known for mutations in HSP3 and HSP4 genes, respectively. Methods: In this study, two Pakhtun consanguineous families, ALB-09 and ALB-10, were enrolled for clinical and molecular diagnoses. Whole-exome sequencing (WES) of the index patient in each family followed by Sanger sequencing of all available samples was performed using 3Billion. Inc South Korea rare disease diagnostics services. Results: The affected individuals of families ALB-09 and ALB-10 showed typical phenotypes of HPS such as oculocutaneous albinism, poor vision, nystagmus, nystagmus-induced involuntary head nodding, bleeding diathesis, and enterocolitis; however, immune system weakness was not recorded. WES analyses of one index patient revealed a novel nonsense variant (NM_032383.4: HSP3; c.2766T > G) in family ALB-09 and a five bp deletion (NM_001349900.2: HSP4; c.1180_1184delGTTCC) variant in family ALB-10. Sanger sequencing confirmed homozygous segregation of the disease alleles in all affected individuals of the respective family. Conclusions: The substitution c.2766T > G creates a premature protein termination at codon 922 in HPS3, replacing tyrosine amino acid with a stop codon (p.Tyr922Ter), while the deletion mutation c.1180_1184delGTTCC leads to a reading frameshift and a premature termination codon adding 23 abnormal amino acids to HSP4 protein (p:Val394Pro395fsTer23). To the best of our knowledge, the two novel variants identified in HPS3 and HPS4 genes causing Hermansky-Pudlak syndrome are the first report from the Pakhtun Pakistani population. Our work expands the pathogenic spectrum of HPS3 and HPS4 genes, provides successful molecular diagnostics, and helps the families in genetic counselling and reducing the disease burden in their future generations.

Genetic variations in the DYNC2H1 gene causing SRTD3 (short-rib thoracic dysplasia 3 with or without polydactyly)

Background and aims: Short-rib thoracic dysplasia 3 with or without polydactyly (SRTD3) represents a type of severe fetal skeletal dysplasia (SD) characterized by shortened limbs, narrow thorax with or without polydactyly, which is caused by the homozygous or compound heterozygous mutations in the DYNC2H1 gene. SRTD3 is a recessive disorder, identification of the responsible genetic variation would be beneficial to an accurate prenatal diagnosis and well-grounded counseling for the affected families. Material and methods: Two families having experienced recurrent fetal SDs were recruited and submitted to a multiplatform genetic investigation. Whole-exome sequencing (WES) was performed with samples collected from the probands. Sanger sequencing and fluorescent quantitative PCR (qPCR) were conducted as validation assays for suspected variations. Results: WES identified two compound heterozygous variations in the DYNC2H1(NM_001080463.2) gene, namely c.2386C>T (p.Arg796Trp) and c.7289T>C (p.Ile2430Thr) for one; and exon (64-83)del and c.8190G>T (p.Leu2730Phe) for the other, respectively. One variant in them, exon (64-83)del, was novelly identified. Conclusion: The study detected two compound heterozygous variation in DYNC2H1 including one novel deletion: exon (64-83) del. Our findings clarified the cause of fetal skeletal dysplasia in the subject families, provided guidance for their future pregnancies, and highlighted the value of WES in diagnosis of skeletal dysplasia with unclear prenatal indications.

A novel homozygous truncating variant in PPFIBP1 further delineates PPFIBP1-associated neurodevelopmental disorder

Neurodevelopmental disorders (NDDs) are classified as a group of disorders affecting function and development of the brain and having wide clinical variability. Herein, we describe two affected individuals segregating a recessive NDD. The affected individuals exhibited phenotypes such as global developmental delay (GDD), intellectual disability (ID), microcephaly and speech delay. Whole-exome sequencing (WES) followed by bidirectional Sanger sequencing techniques identified a homozygous nonsense variant (c.466C > T; p.Gln156*) in the PPFIBP1 gene (NM_003622.4) that segregated with the disease phenotype. Further, to elucidate the effect of the variant on protein structure, 3D protein modelling was performed for the mutant and normal protein that suggested substantial reduction of the mutant protein. Our data support the evidence that PPFIBP1 has a pivotal role in neurodevelopment in humans, and loss-of-function variants cause clinically variable neurodevelopmental phenotypes.

Congenital Defects in a Patient Carrying a Novel Homozygous AEBP1 Variant: Further Expansion of the Phenotypic Spectrum of Ehlers-Danlos Syndrome Classical-like Type 2?

In 2018, a new clinical subtype, caused by biallelic variants in the AEBP1 gene, encoding the ACLP protein, was added to the current nosological classification of the Ehlers-Danlos Syndromes (EDS). This new phenotype, provisionally termed EDS classical-like type 2 (clEDS2), has not yet been fully characterized, as only nine cases have been reported to date. Here we describe a patient, homozygous for a novel AEBP1 pathogenic variant (NM_001129.5 c.2123_2124delTG (p.Val708AlafsTer5)), whose phenotype is reminiscent of classical EDS but also includes previously unreported multiple congenital malformations. Furthermore, we briefly summarize the current principal clinical manifestations of clEDS2 and the molecular evidence surrounding the role of AEBP1 in the context of extracellular matrix homeostasis and connective tissue development. Although a different coexisting etiology for the multiple congenital malformations of our patient cannot be formally excluded, the emerging role of ACLP in TGF-β and WNT pathways may explain their occurrence and the phenotypical variability of clEDS2.

Clinical genetics of spondylocostal dysostosis: A mini review

Spondylocostal dysostosis is a genetic defect associated with severe rib and vertebrae malformations. In recent years, extensive clinical and molecular diagnosis advancements enabled us to identify disease-causing variants in different genes for such severe conditions. The identification of novel candidate genes enabled us to understand the developmental biology and molecular and cellular mechanisms involved in the etiology of these rare diseases. Here, we discuss the clinical and molecular targets associated with spondylocostal dysostosis, including clinical evaluation, genes, and pathways involved. This review might help us understand the basics of such a severe disorder, which might help in proper clinical characterization and help in future therapeutic strategies.