Exome-wide benchmark of difficult-to-sequence regions using short-read next-generation DNA sequencing.

Next-generation DNA sequencing (NGS) in short-read mode has recently been used for genetic testing in various clinical settings. NGS data accuracy is crucial in clinical settings, and several reports regarding quality control of NGS data, primarily focusing on establishing NGS sequence read accuracy, have been published thus far. Variant calling is another critical source of NGS errors that remains unexplored at the single-nucleotide level despite its established significance. In this study, we used a machine-learning-based method to establish an exome-wide benchmark of difficult-to-sequence regions at the nucleotide-residue resolution using 10 genome sequence features based on real-world NGS data accumulated in The Genome Aggregation Database (gnomAD) of the human reference genome sequence (GRCh38/hg38). The newly acquired metric, designated the 'UNMET score,' along with additional lines of structural information from the human genome, allowed us to assess the sequencing challenges within the exonic region of interest using conventional short-read NGS. Thus, the UNMET score could provide a basis for addressing potential sequential errors in protein-coding exons of the human reference genome sequence GRCh38/hg38 in clinical sequencing.

Nanopore long-read sequencing analysis reveals ZIC1 dysregulation caused by a de novo 3q inversion with a breakpoint located 7 kb downstream of ZIC1.

Zic family member 1 (ZIC1), a gene located on chromosome 3q24, encodes a transcription factor with zinc finger domains that is essential for the normal development of the cerebellum. Heterozygous loss-of-function of ZIC1 causes Dandy-Walker malformation, while heterozygous gain-of-function leads to a multiple congenital anomaly syndrome characterized by craniosynostosis, brain abnormalities, facial features, and learning disability. In this study, we present the results of genetic analysis of a male patient with clinically suspected Gomez-Lopez-Hernandez syndrome. The patient displayed multiple congenital abnormalities, including bicoronal craniosynostosis, characteristic facial features, cerebellar malformation with rhombencephalosynapsis, and temporal alopecia, and a de novo inversion of chromosome 3q. Breakpoint analysis using a Nanopore long-read sequencer revealed a breakpoint in the distal centromere of 3q24 located 7 kb downstream of the 3' untranslated region of ZIC1. On the basis of the clinical similarities, we concluded that the abnormalities in this patient were caused by the transcriptional dysregulation of ZIC1. We hypothesize the underlying molecular mechanisms of transcriptional dysregulation of ZIC1 such as the abnormalities in topologically associated domains encompassing ZIC1. This study highlights the usefulness of long-read sequencing in the analysis of de novo balanced chromosomal abnormalities.

A Japanese patient with Teebi hypertelorism syndrome and a novel CDH11 EC1 domain variant.

The gene CDH11 encodes cadherin-11, a Type II cadherin superfamily member that contains five extracellular cadherin (EC) domains. Cadherin-11 undergoes trans-dimerization via the EC1 domain to generate cadherin complexes. Compound heterozygous and homozygous loss-of-function CDH11 variants are observed in Elsahy-Waters syndrome (EWS), which shows characteristic craniofacial features, vertebral abnormalities, cutaneous syndactyly in 2-3 digits, genitourinary anomalies, and intellectual disability. Heterozygous CDH11 variants can cause Teebi hypertelorism syndrome (THS), which features widely spaced eyes and hypospadias. We report a THS patient with a novel CDH11 variant involving the EC1 domain. The patient was a 10-month-old male with normal developmental milestones, but had widely spaced eyes, strabismus, hypospadias, shawl scrotum, broad thumbs (right bifid thumb in x-ray), polysyndactyly of the left fourth finger, and cutaneous syndactyly of left third/fourth fingers. Exome sequencing identified a de novo heterozygous CDH11 variant (NM_001797.4:c.229C > T [p.Leu77Phe] NC_000016.9:g.64998856G > A). Clinical features were consistent with previously reported THS patients, but polysyndactyly, broad thumb, and cutaneous syndactyly overlapped phenotypic features of EWS. THS and EWS may represent a spectrum of CDH11-related disorders. Residue Leu77 in this novel CDH11 variant lines a large hydrophobic pocket where side chains of the partner cadherin-11 insert to trans-dimerize, suggesting that the cadherin-11 structure might be altered in this variant.

GPC4 truncating variant associated with Keipert syndrome and lacrimal punctal agenesis.

Lacrimal punctal agenesis is an extremely rare condition with an unclear genetic basis. Here, we report a 3-year-old male patient harboring a hemizygous variant in glypican 4 (GPC4), which causes Keipert syndrome, who presented with complete lacrimal punctal agenesis, distinctive craniofacial features, mild developmental delay, mild intellectual disability, and autism. The craniofacial features included a prominent forehead, epicanthus, depressed and broad nasal bridge, hypoplastic columella, midface hypoplasia, tented upper lip, and low-set ears. Proband exome sequencing identified a hemizygous variant in GPC4: NM_001448.3:c.1051C > T (p.Arg351*). The GPC4 variant was inherited from his heterozygous mother; X-inactivation followed a skewed pattern in his mother. This patient demonstrated clinical features consistent with Keipert syndrome including craniofacial features, brachydactyly, broad distal phalanx, broad first toe, and mild developmental delay; however, agenesis of the lacrimal puncta has not been reported previously in Keipert syndrome. Our findings suggest that GPC4, which encodes a heparan-sulfate proteoglycan, may play an important role in lacrimal morphogenesis. Our observations also suggest that Keipert syndrome should be considered in patients with lacrimal punctal agenesis.

Noonan syndrome-like phenotype associated with an ERF frameshift variant.

Noonan syndrome is a so-called "RASopathy," that is characterized by short stature, distinctive facial features, congenital heart defects, and developmental delay. Of individuals with a clinical diagnosis of Noonan syndrome, 80%-90% have pathogenic variants in the known genes implicated in the disorder, but the molecular mechanism is unknown in the remaining cases. Heterozygous pathogenic variants of ETS2 repressor factor (ERF), which functions as a repressor in the RAS/MAPK signaling pathway, cause syndromic craniosynostosis. Here, we report an ERF frameshift variant cosegregating with a Noonan syndrome-like phenotype in a family. The proband was a 3-year-old female who presented with dysmorphic facial features, including proptosis, hypertelorism, slightly down slanted palpebral fissures, low-set posteriorly rotated ears, depressed nasal bridge, short stature, and developmental delay. Exome sequencing of the proband identified a heterozygous ERF variant [NM_006494.4: c.185del p.(Glu62Glyfs*15)]. Her mother and sister showed a similar phenotype and had the same heterozygous ERF variant. A large proportion of the previously reported patients with syndromic craniosynostosis and pathogenic ERF variants also showed characteristic features that overlap with those of Noonan syndrome. The present finding supports an association between heterozygous ERF variants and a Noonan syndrome-like phenotype.

Two siblings with acute necrotizing encephalopathy associated with variants of LARS1.

Acute necrotizing encephalopathy (ANE) is a rapidly progressive encephalopathy of unknown etiology. The underlying mechanisms are highly heterogeneous, often including genetic backgrounds. Variants of LARS1, encoding the leucyl-tRNA synthetase 1, are responsible for infantile liver failure syndrome 1. We describe two siblings with ANE caused by compound heterozygous variants of LARS1. Patient 1 was a 17-month-old girl. She presented with generalized seizure and liver dysfunction due to influenza type A infection. Brain magnetic resonance imaging on day 4 of onset showed diffuse high-intensity signals consistent with ANE. She died on day 10. Patient 2, a younger male sibling of patient 1, had mild to moderate developmental delay and growth failure at the age of 18 months. He showed a markedly elevated level of transaminases triggered by infection with human herpesvirus 6. On day 4 of onset, he had generalized seizures. Brain computed tomography showed a diffuse symmetrical hypodensity consistent with ANE. He died on day 7. Whole exome sequencing identified the compound heterozygous variants in LARS1 (NM_020117.11) as c.83_88delinsAATGGGATA, p.(Arg28_Phe30delinsLysTryAspIle) and c.1283C>T, p.(Pro428Leu) in both siblings. The severe neurologic phenotype, found in our patients, reflects the complicated pathogenesis of LARS1-related disorder.

PHACES-like syndrome with TMEM260 compound heterozygous variants.

PHACES syndrome is a multiple congenital disorder with unknown etiology that is characterized by Posterior fossa anomalies, Hemangioma, Arterial lesions, Cardiac abnormalities/coarctation of the aorta, Eye anomalies, and Sternal cleft. Compound heterozygous or homozygous TMEM260 variants cause structural heart defects and renal anomalies syndrome (SHDRA). We describe a 10-year-old male patient with a PHACES-like syndrome and TMEM260 compound heterozygous variants who demonstrated overlapping phenotypes between the two syndromes. He presented with truncus arteriosus, supraumbilical raphe, ophthalmological abnormality, vertebral abnormality, borderline intellectual disability, and hearing loss. He had normal serum creatinine. In proband exome sequencing, compound heterozygous TMEM260 variants (NM_017799.4 c.1617delG p.(Trp539Cysfs*9)/c.1858C > T p.(Gln620*)) were identified. Twelve patients have been reported with TMEM260-related SHDRA: 10 had truncus arteriosus and 6 had renal failure. One previously reported patient had facial port wine nevus and another patient had supraumbilical raphe, which are the cardinal signs for PHACES syndrome. TMEM260-related SHDRA could share overlapping clinical features with PHACES syndrome. This report expands the phenotypic spectrum of a TMEM260-related disorder.

Complex congenital cardiovascular anomaly in a patient with AGO1-associated disorder.

Pathogenic AGO1 variants have been associated with neurodevelopmental disorders, including autism spectrum disorder, developmental delay, intellectual disability, and dysmorphic facial appearance. In mammalian models, defects in microRNA (miRNA) biogenesis are associated with congenital heart disease and dilated cardiomyopathy. We describe the case of a patient with partial anomalous pulmonary venous return, hypoplastic left lung, bilateral pulmonary sequestration, and dilated myocardiopathy. We identified a de novo pathogenic variant of AGO1, which encodes an Argonaute protein forming a gene-silencing complex with microRNAs. The patient was diagnosed with dilated cardiomyopathy with no apparent cause at 3 years of age. She was started on enalapril and carvedilol, and her heart failure was well controlled. We expanded the AGO1-associated phenotype to include complex congenital cardiovascular anomaly and dilated cardiomyopathy in humans.

Divergent variant patterns among 19 patients with Rubinstein-Taybi syndrome uncovered by comprehensive genetic analysis including whole genome sequencing.

Rubinstein-Taybi syndrome (RSTS) is characterized by dysmorphic facial features, broad thumbs, and intellectual disability. CREB-binding protein (CREBBP) or E1A-binding protein P300 (EP300) are causative genes. To elucidate the underlying genetic and genomic architecture related to the RSTS phenotype, we performed comprehensive genetic analysis targeting CREBBP and/or EP300 in 22 clinically diagnosed patients. During the 11-year study period, we used several analysis methods including high-resolution melting, array-based comparative genomic hybridization, panel-based exome sequencing, whole exome sequencing, and whole genome sequencing (WGS). We identified the causative variants in 19 patients (86.3%), but they were variable and complex, so we must combine multiple analysis methods. Notably, we found genetic alterations in the non-coding regions of two patients (10.5%, 2/19): scattered deletions including a partial 5'-untranslated region of CREBBP in one patient (all coding exons were intact), and a deep 229-bp intronic deletion in another patient, resulting in a splicing error. Furthermore, we identified rare clinical findings: two patients with an EP300 variant showed abnormal development of the neural tube, and one patient with a CREBBP variant had anorectal atresia with a cloaca. Our findings expand the allelic heterogeneity of RSTS, underscore the utility of comprehensive genetic analysis, and suggest that WGS may be a practical diagnostic strategy.

Further delineation of SET-related intellectual disability syndrome.

A loss-of-function mutation of SET causes nonsyndromic intellectual disability, often associated with mild facial dysmorphic features, including plagiocephaly, facial asymmetry, broad and high forehead, a wide mouth, and a prominent mandible. We report a male individual with a 2.0 Mb deletion within 9q34.11, involving SET and SPTAN1, but not STXBP1. Among the genes with a high probability of being loss-of-function intolerant in the deletion interval, only SPTAN1 and SET had haploinsufficiency score (%HI) <10, indicating a high likelihood of haploinsufficiency. Pathogenic variants in SPTAN1 are responsible for early-onset epileptic encephalopathy by exerting a dominant-negative effect. However, whether haploinsufficiency of SPTAN1 alone also causes the severe phenotype remained unknown. SET is a regulator of cell differentiation in early human development and a component of the inhibitor of histone acetyltransferases complex. Therefore, combining the previously reported patients, our patient delineated the phenotypic spectrum of SET-related nonsyndromic intellectual disability with mild facial dysmorphism.

Update of the genotype and phenotype of KMT2D and KDM6A by genetic screening of 100 patients with clinically suspected Kabuki syndrome.

Kabuki syndrome is characterized by a variable degree of intellectual disability, characteristic facial features, and complications in various organs. Many variants have been identified in two causative genes, that is, lysine methyltransferase 2D (KMT2D) and lysine demethylase 6A (KDM6A). In this study, we present the results of genetic screening of 100 patients with a suspected diagnosis of Kabuki syndrome in our center from July 2010 to June 2018. We identified 76 variants (43 novel) in KMT2D and 4 variants (3 novel) in KDM6A as pathogenic or likely pathogenic. Rare variants included a deep splicing variant (c.14000-8C>G) confirmed by RNA sequencing and an 18% mosaicism level for a KMT2D mutation. We also characterized a case with a blended phenotype consisting of Kabuki syndrome, osteogenesis imperfecta, and 16p13.11 microdeletion. We summarized the clinical phenotypes of 44 patients including a patient who developed cervical cancer of unknown origin at 16 years of age. This study presents important details of patients with Kabuki syndrome including rare clinical cases and expands our genetic understanding of this syndrome, which will help clinicians and researchers better manage and understand patients with Kabuki syndrome they may encounter.

An efficient genetic test flow for multiple congenital anomalies and intellectual disability.

BACKGROUND: Genetic testing has enabled the diagnosis of multiple congenital anomalies and/or intellectual disabilities. However, because of the phenotypic variability in these disorders, selection of an appropriate genetic test can be difficult and complex. For clinical examination, particularly in clinical facilities, a simple and standardized system is needed. METHODS: We compared microarray comparative genomic hybridization and clinical exome sequencing with regard to diagnostic yield, cost, and time required to reach a definitive diagnosis. After first performing G-banding for 200 patients with multiple congenital anomalies and/or intellectual disability, as a subsequent genetic test, microarray and clinical exome sequencing were compared with regard to diagnostic yield, cost, and time required. RESULTS: There was no obvious difference in the diagnostic rate between the two methods; however, clinical exome sequencing was superior in terms of cost and time. In addition, clinical exome sequencing could sufficiently identify copy number variants, and even smaller copy number variants could be identified. CONCLUSIONS: Clinical exome sequencing should be implemented earlier as a genetic test for undiagnosed patients with multiple congenital anomalies and/or intellectual disabilities. Our results can be used to establish inspection methods in clinical facilities.

A novel gene (FAM20B encoding glycosaminoglycan xylosylkinase) for neonatal short limb dysplasia resembling Desbuquois dysplasia.

Desbuquois dysplasia (DBQD) is an autosomal recessive heterogeneous disorder characterized by joint laxity and skeletal changes, including a distinctive monkey-wrench appearance of the femora, advanced carpal ossification, and abnormal patterning of the preaxial digits. Two genes for DBQD (CANT1 encoding calcium-activated nucleotidase-1 and XYLT1 encoding xylosyltransferase-1) have been reported. We propose a novel gene for neonatal short limb dysplasia resembling DBQD, based on the phenotype and genotype of two affected siblings. The affected boy and girl died in early infancy and shortly after birth, respectively. The clinical hallmarks included mid-face hypoplasia, thoracic hypoplasia with respiratory failure, very short stature (approximately -7 SD of birth length) with mesomelic shortening of the limbs, and multiple dislocations of the large joints. Radiological examinations showed prominent lesser trochanter, flared metaphyses of the long bones, and joint dislocations. The affected boy had preaxial digital hypoplasia, and the affected girl showed overlapping and syndactyly of the preaxial digits. Molecular analyses of the girl showed compound heterozygous variants in FAM20B (NM_014864: c.174_178delTACCT p.T59Afs*19/c.1038delG p.N347Mfs*4). FAM20B encodes glycosaminoglycan xylosylkinase, which acts downstream of xylosyltransferase-1. Given the fact that FAM20B deficiency causes skeletal phenotypes in mice and zebrafish, these variants are highly probable to be pathogenic.

Familial total anomalous pulmonary venous return with 15q11.2 (BP1-BP2) microdeletion.

A 15q11.2 microdeletion (BP1-BP2) is associated with congenital heart diseases (CHDs), developmental delay, and epilepsy. This deletion co-occurs with CHD in 20-30% patients, but a familial case of CHD and a 15q11.2 deletion has not been identified. Here we report the first familial (three siblings) case of total anomalous pulmonary venous return associated with 15q11.2 deletion. Array comparative genomic hybridization identified a ~395 kb deletion at 15q11.2 in patient 1. This deletion was confirmed by fluorescence in situ hybridization in patients 1 and 3 and their asymptomatic father. No deleterious mutation was identified by proband-only exome sequencing of patient 1. One healthy sibling and their mother did not carry the deletion. This deletion is often inherited from asymptomatic parents with an estimated low penetrance of 10.4%. Conversely, we observed high penetrance of this deletion, but secondary copy-number variants or pathogenic variants were not detected in this family.

Delineation of the KIAA2022 mutation phenotype: two patients with X-linked intellectual disability and distinctive features.

Next-generation sequencing has enabled the screening for a causative mutation in X-linked intellectual disability (XLID). We identified KIAA2022 mutations in two unrelated male patients by targeted sequencing. We selected 13 Japanese male patients with severe intellectual disability (ID), including four sibling patients and nine sporadic patients. Two of thirteen had a KIAA2022 mutation. Patient 1 was a 3-year-old boy. He had severe ID with autistic behavior and hypotonia. Patient 2 was a 5-year-old boy. He also had severe ID with autistic behavior, hypotonia, central hypothyroidism, and steroid-dependent nephrotic syndrome. Both patients revealed consistent distinctive features, including upswept hair, narrow forehead, downslanting eyebrows, wide palpebral fissures, long nose, hypoplastic alae nasi, open mouth, and large ears. De novo KIAA2022 mutations (p.Q705X in Patient 1, p.R322X in Patient 2) were detected by targeted sequencing and confirmed by Sanger sequencing. KIAA2022 mutations and alterations have been reported in only four families with nonsyndromic ID and epilepsy. KIAA2022 is highly expressed in the fetal and adult brain and plays a crucial role in neuronal development. These additional patients support the evidence that KIAA2022 is a causative gene for XLID.