Chromoanagenesis Event Underlies a de novo Pericentric and Multiple Paracentric Inversions in a Single Chromosome Causing Coffin-Siris Syndrome

Genomic Balancing Act: deciphering DNA rearrangements in the complex chromosomal aberration involving 5p15.2, 2q31.1, and 18q21.32

Despite extensive research into the genetic underpinnings of neurodevelopmental disorders (NDD), many clinical cases remain unresolved. We studied a female proband with a NDD, mildly dysmorphic facial features, and brain stem hypoplasia on neuroimaging. Comprehensive genomic analyses revealed a terminal 5p loss and a terminal 18q gain in the proband while a diploid copy number for chromosomes 5 and 18 in both parents. Genomic investigations in the proband identified an unbalanced translocation t(5;18) with additional genetic material from chromosome 2 (2q31.3) inserted at the breakpoint, pointing to a complex chromosomal rearrangement (CCR) involving 5p15.2, 2q31.3, and 18q21.32. Breakpoint junction analyses enabled by long-read genome sequencing unveiled the presence of four distinct junctions in the father, who is a carrier of a balanced CCR. The proband inherited from the father both the abnormal chromosome 5 resulting in segmental aneusomies of chr5 (loss) and chr18 (gain) and a der(2) homologue. Evidences suggest a chromoplexy mechanism for this CCR derivation, involving double-strand breaks (DSBs) repaired by non-homologous end joining (NHEJ) or alternative end joining (alt-EJ). The complexity of the CCR and the segregation of homologues elucidate the genetic model for this family. This study demonstrates the importance of combining multiple genomic technologies to uncover genetic causes of complex neurodevelopmental syndromes and to better understand genetic disease mechanisms.

Investigation of a pathogenic inversion in UNC13D and comprehensive analysis of chromosomal inversions across diverse datasets

Inversions are known contributors to the pathogenesis of genetic diseases. Identifying inversions poses significant challenges, making it one of the most demanding structural variants (SVs) to detect and interpret. Recent advancements in sequencing technologies and the development of publicly available SV datasets have substantially enhanced our capability to explore inversions. However, a cross-comparison in those datasets remains unexplored. In this study, we reported a proband with familial hemophagocytic lymphohistiocytosis type-3 carrying a splicing variant (c.1389+1G>A) in trans with an inversion present in 0.006345% of individuals in gnomAD (v4.0) that disrupts UNC13D. Based on this result, we investigate the features of potentially pathogenic inversions in gnomAD which revealed 98.9% of them are rare and disrupt 5% of protein-coding genes associated with a phenotype in OMIM. We then conducted a comparative analysis of additional public datasets, including DGV, 1KGP, and two recent studies from the Human Genome Structural Variation Consortium which revealed common and dataset-specific inversion characteristics suggesting methodology detection biases. Next, we investigated the genetic features of inversions disrupting the protein-coding genes. Notably, we found that the majority of protein-coding genes in OMIM disrupted by inversions are associated with autosomal recessive phenotypes supporting the hypothesis that inversions in trans with other variants are potential hidden causes of monogenic diseases. This effort aims to fill the gap in our understanding of the molecular characteristics of inversions with low frequency in the population and highlight the importance of identifying them in rare disease studies.

Leveraging the T2T assembly to resolve rare and pathogenic inversions in reference genome gaps

Chromosomal inversions (INVs) are particularly challenging to detect due to their copy-number neutral state and association with repetitive regions. Inversions represent about 1/20 of all balanced structural chromosome aberrations and can lead to disease by gene disruption or altering regulatory regions of dosage-sensitive genes in cis Short-read genome sequencing (srGS) can only resolve ∼70% of cytogenetically visible inversions referred to clinical diagnostic laboratories, likely due to breakpoints in repetitive regions. Here, we study 12 inversions by long-read genome sequencing (lrGS) (n = 9) or srGS (n = 3) and resolve nine of them. In four cases, the inversion breakpoint region was missing from at least one of the human reference genomes (GRCh37, GRCh38, T2T-CHM13) and a reference agnostic analysis was needed. One of these cases, an INV9 mappable only in de novo assembled lrGS data using T2T-CHM13 disrupts EHMT1 consistent with a Mendelian diagnosis (Kleefstra syndrome 1; MIM#610253). Next, by pairwise comparison between T2T-CHM13, GRCh37, and GRCh38, as well as the chimpanzee and bonobo, we show that hundreds of megabases of sequence are missing from at least one human reference, highlighting that primate genomes contribute to genomic diversity. Aligning population genomic data to these regions indicated that these regions are variable between individuals. Our analysis emphasizes that T2T-CHM13 is necessary to maximize the value of lrGS for optimal inversion detection in clinical diagnostics. These results highlight the importance of leveraging diverse and comprehensive reference genomes to resolve unsolved molecular cases in rare diseases.

De novo variation in ARID1B gene causes Coffin-Siris syndrome 1 in a Chinese family with excessive early-onset high myopia

Coffin-Siris syndrome (CSS) is a rare autosomal dominant inheritance disorder characterized by distinctive facial features, hypoplasia of the distal phalanx or nail of the fifth and additional digits, developmental or cognitive delay of varying degree, hypotonia, hirsutism/hypertrichosis, sparse scalp hair and varying kind of congenital anomalies. CSS can easily be misdiagnosed as other syndromes or disorders with a similar clinical picture because of their genetic and phenotypic heterogeneity. We describde the genotype-phenotype correlation of one patient from a healthy Chinese family with a novel genotype underlying CSS, who was first diagnosed in the ophthalmology department as early-onset high myopia (eoHM). Comprehensive ophthalmic tests as well as other systemic examinations were performed on participants to confirm the phenotype. The genotype was identified using whole exome sequencing, and further verified the results among other family members by Sanger sequencing. Real-time quantitative PCR (RT-qPCR) technology was used to detect the relative mRNA expression levels of candidate genes between proband and normal family members. The pathogenicity of the identified variant was determined by The American College of Medical Genetics and Genomics (ACMG) guidelines. STRING protein-protein interactions (PPIs) network analysis was used to detect the interaction of candidate gene-related proteins with high myopia gene-related proteins. The patient had excessive eoHM, cone-rod dystrophy, coarse face, excessive hair growth on the face, sparse scalp hair, developmental delay, intellectual disability, moderate hearing loss, dental hypoplasia, patent foramen ovale, chronic non-atrophic gastritis, bilateral renal cysts, cisterna magna, and emotional outbursts with aggression. The genetic assessment revealed that the patient carries a de novo heterozygous frameshift insertion variant in the ARID1B c.3981dup (p.Glu1328ArgfsTer5), which are strongly associated with the typical clinical features of CSS patients. The test results of RT-qPCR showed that mRNA expression of the ARID1B gene in the proband was approximately 30% lower than that of the normal control in the family, suggesting that the variant had an impact on the gene function at the level of mRNA expression. The variant was pathogenic as assessed by ACMG guidelines. Analysis of protein interactions in the STRING online database revealed that the ARID1A protein interacts with the high myopia gene-related proteins FGFR3, ASXL1, ERBB3, and SOX4, whereas the ARID1A protein antagonizes the ARID1B protein. Therefore, in this paper, we are the first to report a de novo heterozygous frameshift insertion variant in the ARID1B gene causing CSS with excessive eoHM. Our study extends the genotypic and phenotypic spectrums for ARID1B-CSS and supplies evidence of significant association of eoHM with variant in ARID1B gene. As CSS has high genetic and phenotypic heterogeneity, our findings highlight the importance of molecular genetic testing and an interdisciplinary clinical diagnostic workup to avoid misdiagnosis as some disorders with similar manifestations of CSS.

Mind the gap: the relevance of the genome reference to resolve rare and pathogenic inversions

Both long-read genome sequencing (lrGS) and the recently published Telomere to Telomere (T2T) reference genome provide increased coverage and resolution across repetitive regions promising heightened structural variant detection and improved mapping. Inversions (INV), intrachromosomal segments which are rotated 180° and inserted back into the same chromosome, are a class of structural variants particularly challenging to detect due to their copy-number neutral state and association with repetitive regions. Inversions represent about 1/20 of all balanced structural chromosome aberrations and can lead to disease by gene disruption or altering regulatory regions of dosage sensitive genes in cis . Here we remapped the genome data from six individuals carrying unsolved cytogenetically detected inversions. An INV6 and INV10 were resolved using GRCh38 and T2T-CHM13. Finally, an INV9 required optical genome mapping, de novo assembly of lrGS data and T2T-CHM13. This inversion disrupted intron 25 of EHMT1, confirming a diagnosis of Kleefstra syndrome 1 (MIM#610253). These three inversions, only mappable in specific references, prompted us to investigate the presence and population frequencies of differential reference regions (DRRs) between T2T-CHM13, GRCh37, GRCh38, the chimpanzee and bonobo, and hundreds of megabases of DRRs were identified. Our results emphasize the significance of the chosen reference genome and the added benefits of lrGS and optical genome mapping in solving rearrangements in challenging regions of the genome. This is particularly important for inversions and may impact clinical diagnostics.

Single variant, yet "double trouble": TSC and KBG syndrome because of a large de novo inversion

Despite the advances in high-throughput sequencing, many rare disease patients remain undiagnosed. In particular, the patients with well-defined clinical phenotypes and established clinical diagnosis, yet missing or partial genetic diagnosis, may hold a clue to more complex genetic mechanisms of a disease that could be missed by available clinical tests. Here, we report a patient with a clinical diagnosis of Tuberous sclerosis, combined with unusual secondary features, but negative clinical tests including TSC1 and TSC2 Short-read whole-genome sequencing combined with advanced bioinformatics analyses were successful in uncovering a de novo pericentric 87-Mb inversion with breakpoints in TSC2 and ANKRD11, which explains the TSC clinical diagnosis, and confirms a second underlying monogenic disorder, KBG syndrome. Our findings illustrate how complex variants, such as large inversions, may be missed by clinical tests and further highlight the importance of well-defined clinical diagnoses in uncovering complex molecular mechanisms of a disease, such as complex variants and "double trouble" effects.

Genomic Balancing Act: Deciphering DNA rearrangements in the Complex Chromosomal Aberration involving 5p15.2, 2q31.1 and 18q21.32

Despite extensive research into the genetic underpinnings of neurodevelopmental disorders (NDD), many clinical cases remain unresolved. We studied a female proband with a NDD, mildly dysmorphic facial features, and brain stem hypoplasia on neuroimaging. Comprehensive genomic analyses revealed a terminal 5p loss and terminal 18q gain in the proband while a diploid copy number for chromosomes 5 and 18 in both parents. Genomic investigations in the proband identified an unbalanced translocation t(5;18) with additional genetic material from chromosome 2 (2q31.3) inserted at the breakpoint, pointing to a complex chromosomal rearrangement (CCR) involving 5p15.2, 2q31.3, and 18q21.32. Breakpoint junction analyses enabled by long read genome sequencing unveiled the presence of four distinct junctions in the father, who is carrier of a balanced CCR. The proband inherited from the father both the abnormal chromosome 5 resulting in segmental aneusomies of chr5 (loss) and chr18 (gain) and a der(2) homologue. Evidences suggest a chromoplexy mechanism for this CCR derivation, involving double-strand breaks (DSBs) repaired by non-homologous end joining (NHEJ) or alternative end joining (alt-EJ). The complexity of the CCR and the segregation of homologues elucidate the genetic model for this family. This study demonstrates the importance of combining multiple genomic technologies to uncover genetic causes of complex neurodevelopmental syndrome and to better understand genetic disease mechanisms.

Complex genomic rearrangements: an underestimated cause of rare diseases

Complex genomic rearrangements (CGRs) are known contributors to disease but are often missed during routine genetic screening. Identifying CGRs requires (i) identifying copy number variants (CNVs) concurrently with inversions, (ii) phasing multiple breakpoint junctions incis, as well as (iii) detecting and resolving structural variants (SVs) within repeats. We demonstrate how combining cytogenetics and new sequencing methodologies is being successfully applied to gain insights into the genomic architecture of CGRs. In addition, we review CGR patterns and molecular features revealed by studying constitutional genomic disorders. These data offer invaluable lessons to individuals interested in investigating CGRs, evaluating their clinical relevance and frequency, as well as assessing their impact(s) on rare genetic diseases.

A chromoanagenesis-driven ultra-complex t(5;7;21)dn truncates neurodevelopmental genes in a disabled boy as revealed by whole-genome sequencing

Germline or constitutional chromoanagenesis-related complex chromosomal rearrangements (CCRs) are rare, apparently "all-at-once", catastrophic events that occur in a single cell cycle, exhibit an unexpected complexity, and sometimes correlate with a severe abnormal phenotype. The term chromoanagenesis encompasses three distinct phenomena, namely chromothripsis, chromoanasynthesis, and chromoplexy. Herein, we found hallmarks of chromothripsis and chromoplexy in an ultra-complex t(5; 7;21)dn involving several disordered breakpoint junctions (BPJs) accompanied by some microdeletions and the disruption of neurodevelopmental genes in a patient with a phenotype resembling autosomal dominant MRD44 (OMIM 617061). G-banded chromosomes and FISH showed that the CCR implied the translocation of the 5p15.2→pter segment onto 7q11.23; in turn, the fragment 7q11.23→qter of der(7) separated into two pieces: the segment q11.23→q32 translocated onto 5p15.2 and fused to 21q22.1→ter in the der(5) while the distal 7q32→qter segment translocated onto der(21) at q22.1. Subsequent whole-genome sequencing unveiled that CCT5, CMBL, RETREG1, MYO10, and TRIO from der(5), IMMP2L, TES, VPS37D, DUS4L, TYW1B, and FEZF1-AS1 from der(7), and TIAM1 and SOD1 from der(21), were disrupted by BPJs, whereas some other genes (predicted to be haplosufficient or inconsequential) were completely deleted. Although remarkably CCT5, TRIO, TES, MYO10, and TIAM1 (and even VPS37D) cooperate in key biological processes for normal neuronal development such as cell adhesion, migration, growth, and/or cytoskeleton formation, the disruption of TRIO most likely caused the patient's MRD44-like phenotype, including intellectual disability, microcephaly, finger anomalies, and facial dysmorphia. Our observation represents the first truncation of TRIO related to a chromoanagenesis event and therefore expands the mutational spectrum of this crucial gene. Moreover, our findings indicate that more than one mechanism is involved in modeling the architecture of ultra-complex rearrangements.

Clan genomics: From OMIM phenotypic traits to genes and biology

Clinical characterization of a patient phenotype has been the quintessential approach for elucidating a differential diagnosis and a hypothesis to explore a potential clinical diagnosis. This has resulted in a language of medicine and a semantic ontology, with both specialty- and subspecialty-specific lexicons, that can be challenging to translate and interpret. There is no 'Rosetta Stone' of clinical medicine such as the genetic code that can assist translation and interpretation of the language of genetics. Nevertheless, the information content embodied within a clinical diagnosis can guide management, therapeutic intervention, and potentially prognostic outlook of disease enabling anticipatory guidance for patients and families. Clinical genomics is now established firmly in medical practice. The granularity and informative content of a personal genome is immense. Yet, we are limited in our utility of much of that personal genome information by the lack of functional characterization of the overwhelming majority of computationally annotated genes in the haploid human reference genome sequence. Whereas DNA and the genetic code have provided a 'Rosetta Stone' to translate genetic variant information, clinical medicine, and clinical genomics provide the context to understand human biology and disease. A path forward will integrate deep phenotyping, such as available in a clinical synopsis in the Online Mendelian Inheritance in Man (OMIM) entries, with personal genome analyses.