Parallel detection of single nucleotide variants and copy number variants with exome analysis: Validation in a cohort of 700 undiagnosed patients

Biallelic structural variants in three patients with ERCC8-related Cockayne syndrome and a potential pitfall of copy number variation analysis

Cockayne syndrome (CS) is a rare autosomal recessive disorder caused by mutations in ERCC8 or ERCC6. Most pathogenic variants in ERCC8 are single nucleotide substitutions. Structural variants (SVs) have been reported in patients with ERCC8-related CS. However, comprehensive molecular detection, including SVs of ERCC8, in CS patients remains problematic. Herein, we present three Japanese patients with ERCC8-related CS in whom causative SVs were identified using whole-exome-based copy number variation (CNV) detection tools. One patient showed compound heterozygosity for a 259-kb deletion and a deletion of exon 4 which has previously been reported as an Asia-specific variant. The other two patients were homozygous for the same exon 4 deletion. The exon 4 deletion was detected only by the ExomeDepth software. Intrigued by the discrepancy in the detection capability of various tools for the SVs, we evaluated the analytic performance of four whole-exome-based CNV detection tools using an exome data set from 337 healthy individuals. A total of 1,278,141 exons were predicted as being affected by the 4 CNV tools. Interestingly 95.1% of these affected exons were detected by one tool alone. Thus, we expect that the use of multiple tools may improve the detection rate of SVs from aligned exome data.

Reanalysis of Chromosomal Microarray Data Using a Smaller Copy Number Variant Call Threshold Identifies Four Cases with Heterozygous Multiexon Deletions of ARID1B, EHMT1, and FOXP1 Genes

Introduction

Chromosomal microarray (CMA) is a highly accurate and established method for detecting copy number variations (CNVs) in clinical genetic testing. CNVs are important etiological factors for disorders such as intellectual disability, developmental delay, and multiple congenital anomalies. Recently developed analytical methods have facilitated the identification of smaller CNVs. Therefore, reanalyzing CMA data using a smaller CNV calling threshold may yield useful information. However, this method was left to the discretion of each institution.

Methods

We reanalyzed the CMA data of 131 patients using a smaller CNV call threshold: 50 kb 50 probes for gain and 25 kb 25 probes for loss. We interpreted the reanalyzed CNVs based on the most recently available information. In the reanalysis, we filtered the data using the Clinical Genome Resource dosage sensitivity gene list as an index to quickly and efficiently check morbid genes.

Results

The number of copy number loss was approximately 20 times greater, and copy number gain was approximately three times greater compared to those in the previous analysis. We detected new likely pathogenic CNVs in four participants: a 236.5 kb loss within ARID1B, a 50.6 kb loss including EHMT1, a 46.5 kb loss including EHMT1, and an 89.1 kb loss within the FOXP1 gene.

Conclusion

The method employed in this study is simple and effective for CMA data reanalysis using a smaller CNV call threshold. Thus, this method is efficient for both ongoing and repeated analyses. This study may stimulate further discussion of reanalysis methodology in clinical laboratories.

Café-au-lait Spots and Cleft Palate: Not a Chance Association

The recognition of syndromic forms of cleft palate is important for condition-specific management. Here, we report a patient with cleft palate, congenital heart disease, intellectual disability, and café-au-lait spots who had a deletion of chromosome 15q14. The identification of the precise breakpoints using a Nanopore-based long-read sequencer showed that the deletion spanned MEIS2 and SPRED1 loci. Cleft palate and café-au-lait spots can be ascribed to MEIS2 and SPRED1, respectively. Patients with cleft palate and café-au-lait spots should be encouraged to undergo a detailed genomic evaluation, including screening for a 15q14 deletion, to enable appropriate anticipatory medico-surgical management and genetic counseling.

Genome Analysis in Sick Neonates and Infants: High-yield Phenotypes and Contribution of Small Copy Number Variations

OBJECTIVE

To delineate the diagnostic efficacy of medical exome, whole exome, and whole genome sequencing according to primary symptoms, the contribution of small copy number variations, and the impact of molecular diagnosis on clinical management.

STUDY DESIGN

This was a prospective study of 17 tertiary care centers in Japan, conducted between April 2019 and March 2021. Critically ill neonates and infants less than 6 months of age were recruited in neonatal intensive care units and in outpatient clinics. The patients underwent medical exome, whole exome, or whole genome sequencing as the first tier of testing. Patients with negative results after medical exome or whole exome sequencing subsequently underwent whole genome sequencing. The impact of molecular diagnosis on clinical management was evaluated through contacting primary care physicians.

RESULTS

Of the 85 patients, 41 (48%) had positive results. Based on the primary symptoms, patients with metabolic phenotypes had the highest diagnostic yield (67%, 4/6 patients), followed by renal (60%, 3/5 patients), and neurologic phenotypes (58%, 14/24 patients). Among them, 4 patients had pathogenic small copy number variations identified using whole genome sequencing. In the 41 patients with a molecular diagnosis, 20 (49%) had changes in clinical management.

CONCLUSIONS

Genome analysis for critically ill neonates and infants had a high diagnostic yield for metabolic, renal, and neurologic phenotypes. Small copy number variations detected using whole genome sequencing contributed to the overall molecular diagnosis in 5% of all the patients. The resulting molecular diagnoses had a significant impact on clinical management.

Methods to Improve Molecular Diagnosis in Genomic Cold Cases in Pediatric Neurology

During the last decade, genetic testing has emerged as an important etiological diagnostic tool for Mendelian diseases, including pediatric neurological conditions. A genetic diagnosis has a considerable impact on disease management and treatment; however, many cases remain undiagnosed after applying standard diagnostic sequencing techniques. This review discusses various methods to improve the molecular diagnostic rates in these genomic cold cases. We discuss extended analysis methods to consider, non-Mendelian inheritance models, mosaicism, dual/multiple diagnoses, periodic re-analysis, artificial intelligence tools, and deep phenotyping, in addition to integrating various omics methods to improve variant prioritization. Last, novel genomic technologies, including long-read sequencing, artificial long-read sequencing, and optical genome mapping are discussed. In conclusion, a more comprehensive molecular analysis and a timely re-analysis of unsolved cases are imperative to improve diagnostic rates. In addition, our current understanding of the human genome is still limited due to restrictions in technologies. Novel technologies are now available that improve upon some of these limitations and can capture all human genomic variation more accurately. Last, we recommend a more routine implementation of high molecular weight DNA extraction methods that is coherent with the ability to use and/or optimally benefit from these novel genomic methods.

A patient with compound heterozygosity of SMPD4: Another example of utility of exome-based copy number analysis in autosomal recessive disorders

For the efficient diagnosis of rare and undiagnosed diseases, the parallel detection of copy number variants (CNVs) and single nucleotide variants using exome analysis is required. Recently, our group reported the usefulness of a program called EXCAVATOR2, which screens for CNVs from aligned exome data in bam format. This method is expected to contribute to the identification of structural variants and to improve the diagnosis rate, especially for the diagnosis of autosomal recessive disease, when a conventional exome analysis identifies a pathogenic variant in one allele but not the other. Here we report a 2-year-old Japanese boy with an undiagnosed disease. He had severe neonatal asphyxia, severe intellectual disability, intractable seizures, cerebellar and brainstem hypoplasia and dysmorphic features including a prominent supraorbital ridge, thin upper lip, and prominent antihelix. An exome analysis reinforced with a copy number analysis using the EXCAVATOR2 method revealed that the patient had a hemizygous variant in chr2(GRCh37):g.130925108G>A, NM_017951.4 c.832C>T, p.(Arg278*) in SMPD4 that was derived from his father and a deletion of SMPD4 derived from his mother. The presence of the deletion spanning SMPD4 was confirmed by short-read and long-read whole-genome sequencing. The successful diagnosis of this reported patient demonstrates the diagnostic utility of EXCAVATOR2 and overcomes the weakness of exome analysis for the detection of autosomal recessive diseases in nonconsanguineous families, significantly impacting genetic counseling for family planning.

A Recurrent De Novo Terminal Duplication of 14q32 in Korean Siblings Associated with Developmental Delay and Intellectual Disability, Growth Retardation, Facial Dysmorphism, and Cerebral Infarction: A Case Report and Literature Review

The terminal 14q32 duplication has been reported often in association with other cytogenetic abnormalities, and individuals with this specific duplication showed varying degrees of developmental delay/intellectual disability (DD/ID) and growth retardation (GR), and distinct facial dysmorphisms. Herein, based on the limited cases of terminal duplication of 14q32 known to date, we present new affected siblings presenting with DD/ID, GR, and facial dysmorphism, as well as cerebral infarction caused by recurrent de novo der(14)t(14;14)(p11.2;q32.1) leading to terminal duplication of 14q32. We used coverage analysis generated via duo exome sequencing, performed chromosomal microarray (CMA) as a confirmatory test, and compared our findings with those reported previously. Coverage analysis generated via duo exome sequencing revealed a 17.2 Mb heterozygous duplication at chromosome 14q32.11-q32.33 with a Z ratio ranging between 0.5 and 1 in the proband and her elder brother. As a complementary method, CMA established a terminal duplication described as the arr[hg19]14q32.11q32.33(90,043,558_107,258,824)x3 in the proband and her elder brother; however, the parents and other siblings showed normal karyotyping and no abnormal gain or loss of CMA results. Five candidate genes, BCL11B, CCNK, YY1, DYNC1H1, and PACS2, were associated with the clinical phenotypes in our cases. Although the parents had normal chromosomes, two affected cases carrying terminal duplication of 14q32 can be explained by gonadal mosaicism. Further studies are needed to establish the association between cerebrovascular events and terminal duplication of chromosome 14q32, including investigation into the cytogenetics of patients with precise clinical descriptions.

Role of chimeric transcript formation in the pathogenesis of birth defects

Chimeric transcripts are formed by chromosomal aberrations. Little is known about the role of chimeric transcripts in the pathogenesis of birth defects. We reanalyzed RNA-seq data in alignment map files from the peripheral blood of 56 patients in whom the diagnoses could not be confirmed by standard exome analysis and transcriptome analysis to screen for chimeric transcripts using a dedicated software, ChimPipe. Chimeric analysis led to a diagnosis in two of the 56 patients: (a) the first patient had a chimeric transcript spanning the causative gene ZEB2 and the GTDC1 gene in its neighboring locus. RNA-seq revealed reads spanning exon 5 of ZEB2 and exon 7 of GTDC1. Whole genome sequencing revealed a 436-kb deletion spanning intron 4 of ZEB2 and intron 7 of GTDC1 and the diagnosis of Mowat-Wilson syndrome was made. (b) The second patient had a chimeric transcript spanning the causative gene KCNK9 and the TRAPPC9 gene in its neighboring locus. RNA-seq revealed reads spanning exon 21 of TRAPPC9 and exon 1 of KCNK9. Whole genome sequencing revealed a 186-kb deletion spanning intron 20 of TRAPPC9 and intron 1 of KCNK9 in this patient. KCNK9 gene is a maternally expressed imprinted gene. The diagnosis of Birk-Barel syndrome was made. Thus, both patients had chimeric transcripts that were directly involved in the pathogenesis of the birth defects. The approach reported herein, of detecting chimeric transcripts from RNA-seq data, is unique in that the approach does not rely on any prior information on the presence of genomic deletion.

The phenotype of 15 cases with rare 8q24.13-q24.3 deletions-A new syndrome or still an enigma?

Diagnosis of rare copy number variants (CNVs) with scarce literature evidence poses a major challenge for interpretation of the clinical significance of chromosomal microarray analysis (CMA) results, especially in the prenatal setting. Bioinformatic tools can be used to assist in this issue; however, this prediction can be imprecise. Our objective was to describe the phenotype of the rare copy number losses encompassing the 8q24.13-q24.3 locus, and to find common features in terms of genomic coordinates, gene content, and clinical phenotypic characteristics. Appropriate cases were retrieved using local databases of two largest Israeli centers performing CMA analysis. In addition, literature and public databases search was performed. Local database search yielded seven new patients with del (8)(q24.13q24.3) (one of these with an additional copy number variant). Literature and public databases search yielded eight additional patients. The cases showed high phenotypic variability, ranging from asymptomatic adults and fetuses with normal ultrasound to patients with autism/developmental delay (6/11 postnatal cases, 54.5%). No clear association was noted between the specific disease-causing/high-pLI gene content of the described del (8)(q24.13q24.3) to neurodevelopmental disorders, except for a possibly relevant locus encompassing the KCNQ3 gene. We present the challenges in classification of rare variants with limited clinical information. In such cases, genotype-phenotype correlation must be assessed with extra-caution and possibly using additional methods to assist the classification, especially in the prenatal setting.

Biallelic loss of OTUD7A causes severe muscular hypotonia, intellectual disability, and seizures

The heterozygous deletion of 15q13.3 is a recurrently observed microdeletion syndrome associated with a relatively mild phenotype including learning disability and language impairment. In contrast, the homozygous deletion of 15q13.3 is extremely rare and is associated with a much severer phenotype that includes epileptic encephalopathy, profound intellectual disability, and hypotonia. Which of the genes within the deleted interval is responsible for the more severe features when biallelically deleted is currently unknown. Here, we report a patient with profound hypotonia, severe intellectual disability, and seizures who had biallelic loss-of-function variants in OTUD7A: a 15q13.3 deletion including the OTUD7A locus, and a frameshift OTUD7A variant c.1125del, p.(Glu375Aspfs*11). Unexpectedly, both aberrations occurred de novo. Our experiment using Caenorhabditis elegans showed that worms carrying a corresponding homozygous variant in the homolog OTUB-2 exhibited weakened muscle contraction suggestive of aberrant neuromuscular transmission. We concluded that the biallelic complete loss of OTUD7A in humans represents a presumably new autosomal recessive disorder characterized by profound hypotonia, severe intellectual disability, and seizures.