Comprehensive structural variation genome map of individuals carrying complex chromosomal rearrangements

A national long-read sequencing study on chromosomal rearrangements uncovers hidden complexities

Clinical genetic laboratories often require a comprehensive analysis of chromosomal rearrangements/structural variants (SVs), from large events like translocations and inversions to supernumerary ring/marker chromosomes and small deletions or duplications. Understanding the complexity of these events and their clinical consequences requires pinpointing breakpoint junctions and resolving the derivative chromosome structure. This task often surpasses the capabilities of short-read sequencing technologies. In contrast, long-read sequencing techniques present a compelling alternative for clinical diagnostics. Here, Genomic Medicine Sweden-Rare Diseases has explored the utility of HiFi Revio long-read genome sequencing (lrGS) for digital karyotyping of SVs nationwide. The 16 samples from 13 families were collected from all Swedish healthcare regions. Prior investigations had identified 16 SVs, ranging from simple to complex rearrangements, including inversions, translocations, and copy number variants. We have established a national pipeline and a shared variant database for variant calling and filtering. Using lrGS, 14 of the 16 known SVs are detected. Of these, 13 are mapped at nucleotide resolution, and one complex rearrangement is only visible by read depth. Two Chromosome 21 rearrangements, one mosaic, remain undetected. Average read lengths are 8.3-18.8 kb with coverage exceeding 20× for all samples. De novo assembly results in a limited number of phased contigs per individual (N50 6-86 Mb), enabling direct characterization of the chromosomal rearrangements. In a national pilot study, we demonstrate the utility of HiFi Revio lrGS for analyzing chromosomal rearrangements. Based on our results, we propose a 5-year plan to expand lrGS use for rare disease diagnostics in Sweden.

Characterization of cryptic complex chromosome rearrangements in balanced chromosomal rearrangement carriers and their PGT-SR clinical outcome assessments

Several reports have presented that balanced chromosomal rearrangements (BCRs) carriers with normal phenotypes may be carriers of complex rearrangements. However, the incidence and PGT clinical outcomes of cryptic complex chromosome rearrangements (CCCRs) in individuals with BCRs is remain unknown. We recruited a cohort of 1,264 individuals with BCR carriers from 2016 to 2021 at the Reproductive and Genetic Hospital of CITIC Xiangya. Peripheral blood was collected for karyotyping and genomic DNA extraction and the PGT-SR clinical outcomes of CCCRs carriers were analyzed and compared with those of BCR carriers. Our findings revealed that 3.6% (45/1,264) of BCR carriers had CCCRs, involving 3-25 breakpoints on 1-3 chromosomes. Furthermore, when mate-pair sequencing was employed, 63.3% (19/30) of CCCR carriers were found to have chromosome rearrangements that were different from those identified by the MicroSeq technique. And the transferable embryo rate of CCCR carriers with 3 chromosomes was significantly lower than that of CCCR carriers with only 1-2 chromosomes. In this research, we revealed that some of the BCR carriers were actually CCCR carriers, and the prognosis of PGT in CCCR carriers with one or two chromosomes is better than that of CCCR carriers with three chromosomes.

A combination of long- and short-read genomics reveals frequent p-arm breakpoints within chromosome 21 complex genomic rearrangements

Purpose

Although chromosome 21 is the smallest human chromosome, it is highly relevant in the pathogenicity of both cancer and congenital diseases, including Alzheimer disease and trisomy 21 (Down syndrome). In addition, cases with rare structural variants (SVs) of chromosome 21 have been reported. These events vary in size and include large chromosomal events, such as ring chromosomes and small partial aneuploidies. The p-arm of the acrocentric chromosome 21 was devoid of reference genomic sequence in GRCh37 and GRCh38, which hampered our ability to solve genomic rearrangements and find the mechanism of formation of disease-causing SVs. We hypothesize that conserved satellite structures and segmental duplications located on the p-arm play an important role in the formation of complex SVs involving chromosome 21.

Methods

Three cases with complex chromosome 21 rearrangements were studied with a combination of short-read and long-read genome sequencing, as well as optical genome mapping. The data were aligned to the T2T-CHM13 assembly.

Results

We were able to resolve all 3 complex chromosome 21 rearrangements in which 15, 8, and 26 breakpoints were identified, respectively. By comparing the identified SV breakpoints, we were able to pinpoint a region between 21p13 and 21p12 that appears to be frequently involved in chromosome 21 rearrangements. Importantly, we observed acrocentric satellite DNA at several breakpoint junctions suggesting an important role for those elements in the formation of complex SVs.

Conclusion

Taken together, our results provide further insights into the architecture and underlying mechanisms of complex rearrangements on acrocentric chromosomes.

Long-read sequencing and optical mapping generates near T2T assemblies that resolves a centromeric translocation

Long-read genome sequencing (lrGS) is a promising method in genetic diagnostics. Here we investigate the potential of lrGS to detect a disease-associated chromosomal translocation between 17p13 and the 19 centromere. We constructed two sets of phased and non-phased de novo assemblies; (i) based on lrGS only and (ii) hybrid assemblies combining lrGS with optical mapping using lrGS reads with a median coverage of 34X. Variant calling detected both structural variants (SVs) and small variants and the accuracy of the small variant calling was compared with those called with short-read genome sequencing (srGS). The de novo and hybrid assemblies had high quality and contiguity with N50 of 62.85 Mb, enabling a near telomere to telomere assembly with less than a 100 contigs per haplotype. Notably, we successfully identified the centromeric breakpoint of the translocation. A concordance of 92% was observed when comparing small variant calling between srGS and lrGS. In summary, our findings underscore the remarkable potential of lrGS as a comprehensive and accurate solution for the analysis of SVs and small variants. Thus, lrGS could replace a large battery of genetic tests that were used for the diagnosis of a single symptomatic translocation carrier, highlighting the potential of lrGS in the realm of digital karyotyping.

Integrating Optical Genome Mapping and Whole Genome Sequencing in Somatic Structural Variant Detection

Structural variants drive tumorigenesis by disrupting normal gene function through insertions, inversions, translocations, and copy number changes, including deletions and duplications. Detecting structural variants is crucial for revealing their roles in tumor development, clinical outcomes, and personalized therapy. Presently, most studies rely on short-read data from next-generation sequencing that aligns back to a reference genome to determine if and, if so, where a structural variant occurs. However, structural variant discovery by short-read sequencing is challenging, primarily because of the difficulty in mapping regions of repetitive sequences. Optical genome mapping (OGM) is a recent technology used for imaging and assembling long DNA strands to detect structural variations. To capture the structural variant landscape more thoroughly in the human genome, we developed an integrated pipeline that combines Bionano OGM and Illumina whole-genome sequencing and applied it to samples from 29 pediatric B-ALL patients. The addition of OGM allowed us to identify 511 deletions, 506 insertions, 93 duplications/gains, and 145 translocations that were otherwise missed in the short-read data. Moreover, we identified several novel gene fusions, the expression of which was confirmed by RNA sequencing. Our results highlight the benefit of integrating OGM and short-read detection methods to obtain a comprehensive analysis of genetic variation that can aid in clinical diagnosis, provide new therapeutic targets, and improve personalized medicine in cancers driven by structural variation.

Prenatal Diagnosis of a de novo 2q14.3-q22.1 Deletion with Complex Chromosomal Rearrangement

Introduction

Chromosomal aberrations due to complex chromosomal rearrangements (CCRs) can cause abnormal phenotypes if accompanied by microdeletions or microduplications near the breakpoint, or gene breaks.

Case Presentation

We report a prenatal diagnostic case of 2q14.3-q22.1 deletion with ultrasound suggestive of absent nasal bone accompanied by CCRs involving 6 chromosomes. Cytogenetic analysis revealed a karyotype of 46,XY,der(1)t(1;2)(p13.3;p11.2),der(2)t(1;2)inv(2)(q12q14.2)del(2)(q14.3q22.1),t(12;16)(q21.2;q12.1),t(13;21)(q32;q22.1). Chromosomal microarray analysis identified a 14.90 Mb deletion on 2q14.3q22.1. The copy number variant was de novo, as determined by karyotype analysis of the parents' peripheral blood G-banding.

Conclusion

The region contains haploinsufficient genes that can cause different phenotypes, mainly associated with neurodevelopmental and autism spectrum disorders. However, the genotype-phenotype correlation is limited in prenatal evaluation. Therefore, the combined use of multiple diagnostic techniques has an important role in the assessment of CCRs and genetic counseling.

Analysis of complex chromosomal rearrangement involving chromosome 6 via the integration of optical genomic mapping and molecular cytogenetic methodologies

Complex chromosomal rearrangements (CCRs) can result in spontaneous abortions, infertility, and malformations in newborns. In this study, we explored a familial CCR involving chromosome 6 by combining optical genomic mapping (OGM) and molecular cytogenetic methodologies. Within this family, the father and the paternal grandfather were both asymptomatic carriers of an identical balanced CCR, while the two offspring with an unbalanced paternal-origin CCR and two microdeletions presented with clinical manifestation. The first affected child, a 5-year-old boy, exhibited neurodevelopmental delay, while the second, a fetus, presented with hydrops fetalis. SNP-genotype analysis revealed a recombination event during gamete formation in the father that may have contributed to the deletion in his offspring. Meanwhile, the couple's haplotypes will facilitate the selection of normal gametes in the setting of assisted reproduction. Our study demonstrated the potential of OGM in identifying CCRs and its ability to work with current methodologies to refine precise breakpoints and construct accurate haplotypes for couples with a CCR.

Identification of complex and cryptic chromosomal rearrangements by optical genome mapping

BACKGROUND

Optical genome mapping (OGM) has developed into a highly promising method for detecting structural variants (SVs) in human genomes. Complex chromosomal rearrangements (CCRs) and cryptic translocations are rare events that are considered difficult to detect by routine cytogenetic methods. In this study, OGM was applied to delineate the precise chromosomal rearrangements in three cases with uncertain or unconfirmed CCRs detected by conventional karyotyping and one case with a cryptic translocation suggested by fetal chromosomal microarray analysis (CMA).

RESULTS

In the three cases with CCRs, OGM not only confirmed or revised the original karyotyping results but also refined the precise chromosomal structures. In the case with a suspected translocation not detected by karyotyping, OGM efficiently identified the cryptic translocation and defined the genomic breakpoints with relatively high accuracy.

CONCLUSIONS

Our study confirmed OGM as a robust alternative approach to karyotyping for the detection of chromosomal structural rearrangements, including CCRs and cryptic translocations.

Optical Genome Mapping for Comprehensive Assessment of Chromosomal Aberrations and Discovery of New Fusion Genes in Pediatric B-Acute Lymphoblastic Leukemia

PURPOSE

To assess the potential added value of Optical Genomic Mapping (OGM) for identifying chromosomal aberrations.

METHODS

We utilized Optical Genomic Mapping (OGM) to determine chromosomal aberrations in 46 children with B-cell Acute lymphoblastic leukemia ALL (B-ALL) and compared the results of OGM with conventional technologies. Partial detection results were verified by WGS and PCR.

RESULTS

OGM showed a good concordance with conventional cytogenetic techniques in identifying the reproducible and pathologically significant genomic SVs. Two new fusion genes (LMNB1::PPP2R2B and TMEM272::KDM4B) were identified by OGM and verified by WGS and RT-PCR for the first time. OGM has a greater ability to detect complex chromosomal aberrations, refine complicated karyotypes, and identify more SVs. Several novel fusion genes and single-gene alterations, associated with definite or potential pathologic significance that had not been detected by traditional methods, were also identified.

CONCLUSION

OGM addresses some of the limitations associated with conventional cytogenomic testing. This all-in-one process allows the detection of most major genomic risk markers in one test, which may have important meanings for the development of leukemia pathogenesis and targeted drugs.

Multi-omics analysis reveals multiple mechanisms causing Prader-Willi like syndrome in a family with a X;15 translocation

Prader-Willi syndrome (PWS; MIM# 176270) is a neurodevelopmental disorder caused by the loss of expression of paternally imprinted genes within the PWS region located on 15q11.2. It is usually caused by either maternal uniparental disomy of chromosome 15 (UPD15) or 15q11.2 recurrent deletion(s). Here, we report a healthy carrier of a balanced X;15 translocation and her two daughters, both with the karyotype 45,X,der(X)t(X;15)(p22;q11.2),-15. Both daughters display symptoms consistent with haploinsufficiency of the SHOX gene and PWS. We explored the architecture of the derivative chromosomes and investigated effects on gene expression in patient-derived neural cells. First, a multiplex ligation-dependent probe amplification methylation assay was used to determine the methylation status of the PWS-region revealing maternal UPD15 in daughter 2, explaining her clinical symptoms. Next, short read whole genome sequencing and 10X genomics linked read sequencing was used to pinpoint the exact breakpoints of the translocation. Finally, we performed transcriptome sequencing on neuroepithelial stem cells from the mother and from daughter 1 and observed biallelic expression of genes in the PWS region (including SNRPN) in daughter 1. In summary, our multi-omics analysis highlights two different PWS mechanisms in one family and provide an example of how structural variation can affect imprinting through long-range interactions.

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.

Analyses of rare predisposing variants of lung cancer in 6,004 whole genomes in Chinese

We present the largest whole-genome sequencing (WGS) study of non-small cell lung cancer (NSCLC) to date among 6,004 individuals of Chinese ancestry, coupled with 23,049 individuals genotyped by SNP array. We construct a high-quality haplotype reference panel for imputation and identify 20 common and low-frequency loci (minor allele frequency [MAF] ≥ 0.5%), including five loci that have never been reported before. For rare loss-of-function (LoF) variants (MAF < 0.5%), we identify BRCA2 and 18 other cancer predisposition genes that affect 5.29% of individuals with NSCLC, and 98.91% (181 of 183) of LoF variants have not been linked previously to NSCLC risk. Promoter variants of BRCA2 also have a substantial effect on NSCLC risk, and their prevalence is comparable with BRCA2 LoF variants. The associations are validated in an independent case-control study including 4,410 individuals and a prospective cohort study including 23,826 individuals. Our findings not only provide a high-quality reference panel for future array-based association studies but depict the whole picture of rare pathogenic variants for NSCLC.

Multi-Omic Investigations of a 17-19 Translocation Links MINK1 Disruption to Autism, Epilepsy and Osteoporosis

Balanced structural variants, such as reciprocal translocations, are sometimes hard to detect with sequencing, especially when the breakpoints are located in repetitive or insufficiently mapped regions of the genome. In such cases, long-range information is required to resolve the rearrangement, identify disrupted genes and, in symptomatic carriers, pinpoint the disease-causing mechanisms. Here, we report an individual with autism, epilepsy and osteoporosis and a de novo balanced reciprocal translocation: t(17;19) (p13;p11). The genomic DNA was analyzed by short-, linked- and long-read genome sequencing, as well as optical mapping. Transcriptional consequences were assessed by transcriptome sequencing of patient-specific neuroepithelial stem cells derived from induced pluripotent stem cells (iPSC). The translocation breakpoints were only detected by long-read sequencing, the first on 17p13, located between exon 1 and exon 2 of MINK1 (Misshapen-like kinase 1), and the second in the chromosome 19 centromere. Functional validation in induced neural cells showed that MINK1 expression was reduced by >50% in the patient’s cells compared to healthy control cells. Furthermore, pathway analysis revealed an enrichment of changed neural pathways in the patient’s cells. Altogether, our multi-omics experiments highlight MINK1 as a candidate monogenic disease gene and show the advantages of long-read genome sequencing in capturing centromeric translocations.

Interstitial deletion 4p15.32p16.1 and complex chromoplexy in a female proband with severe neurodevelopmental delay, growth failure and dysmorphism

Complex chromosomal rearrangements involve the restructuring of genetic material within a single chromosome or across multiple chromosomes. These events can cause serious human disease by disrupting coding DNA and gene regulatory elements via deletions, duplications, and structural rearrangements. Here we describe a 5-year-old female with severe developmental delay, dysmorphic features, multi-suture craniosynostosis, and growth failure found to have a complex series of balanced intra- and inter-chromosomal rearrangements involving chromosomes 4, 11, 13, and X. Initial clinical studies were performed by karyotype, chromosomal microarray, and FISH with research-based short-read genome sequencing coupled with sanger sequencing to precisely map her breakpoints to the base pair resolution to understand the molecular basis of her phenotype. Genome analysis revealed two pathogenic deletions at 4p16.1-p15.32 and 4q31.1, accounting for her developmental delay and dysmorphism. We identified over 60 breakpoints, many with blunt ends and limited homology, supporting a role for non-homologous end joining in restructuring and resolution of the seminal chromoplexy event. We propose that the complexity of our patient’s genomic rearrangements with a high number of breakpoints causes dysregulation of gene expression by three-dimensional chromatin interactions or topologically associating domains leading to growth failure and craniosynostosis. Our work supports an important role for genome sequencing in understanding the molecular basis of complex chromosomal rearrangements in human disease.

Linked-read whole-genome sequencing resolves common and private structural variants in multiple myeloma

Linked-read WGS can be performed without DNA purification and allows for resolution of the diverse structural variants found in MM.

Linked-read WGS can, as a standalone assay, provide comprehensive genetics in myeloma and other diseases with complex genomes.

Multiple myeloma (MM) is an incurable and aggressive plasma cell malignancy characterized by a complex karyotype with multiple structural variants (SVs) and copy-number variations (CNVs). Linked-read whole-genome sequencing (lrWGS) allows for refined detection and reconstruction of SVs by providing long-range genetic information from standard short-read sequencing. This makes lrWGS an attractive solution for capturing the full genomic complexity of MM. Here we show that high-quality lrWGS data can be generated from low numbers of cells subjected to fluorescence-activated cell sorting (FACS) without DNA purification. Using this protocol, we analyzed MM cells after FACS from 37 patients with MM using lrWGS. We found high concordance between lrWGS and fluorescence in situ hybridization (FISH) for the detection of recurrent translocations and CNVs. Outside of the regions investigated by FISH, we identified >150 additional SVs and CNVs across the cohort. Analysis of the lrWGS data allowed for resolution of the structure of diverse SVs affecting the MYC and t(11;14) loci, causing the duplication of genes and gene regulatory elements. In addition, we identified private SVs causing the dysregulation of genes recurrently involved in translocations with the IGH locus and show that these can alter the molecular classification of MM. Overall, we conclude that lrWGS allows for the detection of aberrations critical for MM prognostics and provides a feasible route for providing comprehensive genetics. Implementing lrWGS could provide more accurate clinical prognostics, facilitate genomic medicine initiatives, and greatly improve the stratification of patients included in clinical trials.

The Clinical Utility of Optical Genome Mapping for the Assessment of Genomic Aberrations in Acute Lymphoblastic Leukemia

Simple Summary

The stratification of childhood ALL is currently based on various diagnostic assays. This study investigates the feasibility of Optical Genome Mapping (OGM) to determine the genetic risk profile of ALL using fresh and frozen blood cells in an all-in-one approach. Acute lymphoblastic leukemia samples with data available from SNP-array/array-CGH, RNA-Seq, MLPA, karyotyping and FISH were compared to results obtained by OGM. We show that OGM has the potential to simplify the diagnostic workflow and to identify new structural variants helpful for classifying patients into treatment groups.

Abstract

Acute lymphoblastic leukemia (ALL) is the most prevalent type of cancer occurring in children. ALL is characterized by structural and numeric genomic aberrations that strongly correlate with prognosis and clinical outcome. Usually, a combination of cyto- and molecular genetic methods (karyotyping, array-CGH, FISH, RT-PCR, RNA-Seq) is needed to identify all aberrations relevant for risk stratification. We investigated the feasibility of optical genome mapping (OGM), a DNA-based method, to detect these aberrations in an all-in-one approach. As proof of principle, twelve pediatric ALL samples were analyzed by OGM, and results were validated by comparing OGM data to results obtained from routine diagnostics. All genomic aberrations including translocations (e.g., dic(9;12)), aneuploidies (e.g., high hyperdiploidy) and copy number variations (e.g., IKZF1, PAX5) known from other techniques were also detected by OGM. Moreover, OGM was superior to well-established techniques for resolution of the more complex structure of a translocation t(12;21) and had a higher sensitivity for detection of copy number alterations. Importantly, a new and unknown gene fusion of JAK2 and NPAT due to a translocation t(9;11) was detected. We demonstrate the feasibility of OGM to detect well-established as well as new putative prognostic markers in an all-in-one approach in ALL. We hope that these limited results will be confirmed with testing of more samples in the future.

Comparison of structural variants detected by optical mapping with long-read next-generation sequencing

MOTIVATION

Recent studies have shown the potential of using long-read whole-genome sequencing (WGS) approaches and optical mapping (OM) for the detection of clinically relevant structural variants (SVs) in cancer research. Three main long-read WGS platforms are currently in use: Pacific Biosciences (PacBio), Oxford Nanopore Technologies (ONT) and 10x Genomics. Recently, whole-genome OM technology (Bionano Genomics) has been introduced into human diagnostics. Questions remain about the accuracy of these long-read sequencing platforms, how comparable/interchangeable they are when searching for SVs and to what extent they can be replaced or supplemented by OM. Moreover, no tool can effectively compare SVs obtained by OM and WGS.

RESULTS

This study compared optical maps of the breast cancer cell line SKBR3 with AnnotSV outputs from WGS platforms. For this purpose, a software tool with comparative and filtering features was developed. The majority of SVs up to a 50 kbp distance variance threshold found by OM were confirmed by all WGS platforms, and 99% of translocations and 80% of deletions found by OM were confirmed by both PacBio and ONT, with ∼70% being confirmed by 10x Genomics in combination with PacBio and/or ONT. Interestingly, long deletions (>100 kbp) were detected only by 10x Genomics. Regarding insertions, ∼72% was confirmed by PacBio and ONT, but none by 10x Genomics. Inversions and duplications detected by OM were not detected by WGS. Moreover, the tool enabled the confirmation of SVs that overlapped in the same gene(s) and was applied to the filtering of disease-associated SVs.

AVAILABILITY

https://github.com/novosadt/om-annotsv-svc.

A New Catalog of Structural Variants in 1,301 A. thaliana Lines from Africa, Eurasia, and North America Reveals a Signature of Balancing Selection at Defense Response Genes

Genomic variation in the model plant Arabidopsis thaliana has been extensively used to understand evolutionary processes in natural populations, mainly focusing on single-nucleotide polymorphisms. Conversely, structural variation has been largely ignored in spite of its potential to dramatically affect phenotype. Here, we identify 155,440 indels and structural variants ranging in size from 1 bp to 10 kb, including presence/absence variants (PAVs), inversions, and tandem duplications in 1,301 A. thaliana natural accessions from Morocco, Madeira, Europe, Asia, and North America. We show evidence for strong purifying selection on PAVs in genes, in particular for housekeeping genes and homeobox genes, and we find that PAVs are concentrated in defense-related genes (R-genes, secondary metabolites) and F-box genes. This implies the presence of a "core" genome underlying basic cellular processes and a "flexible" genome that includes genes that may be important in spatially or temporally varying selection. Further, we find an excess of intermediate frequency PAVs in defense response genes in nearly all populations studied, consistent with a history of balancing selection on this class of genes. Finally, we find that PAVs in genes involved in the cold requirement for flowering (vernalization) and drought response are strongly associated with temperature at the sites of origin.

Whole Genome Sequencing in the Evaluation of Fetal Structural Anomalies: A Parallel Test with Chromosomal Microarray Plus Whole Exome Sequencing

Whole genome sequencing (WGS) is a powerful tool for postnatal genetic diagnosis, but relevant clinical studies in the field of prenatal diagnosis are limited. The present study aimed to prospectively evaluate the utility of WGS compared with chromosomal microarray (CMA) and whole exome sequencing (WES) in the prenatal diagnosis of fetal structural anomalies. We performed trio WGS (≈40-fold) in parallel with CMA in 111 fetuses with structural or growth anomalies, and sequentially performed WES when CMA was negative (CMA plus WES). In comparison, WGS not only detected all pathogenic genetic variants in 22 diagnosed cases identified by CMA plus WES, yielding a diagnostic rate of 19.8% (22/110), but also provided additional and clinically significant information, including a case of balanced translocations and a case of intrauterine infection, which might not be detectable by CMA or WES. WGS also required less DNA (100 ng) as input and could provide a rapid turnaround time (TAT, 18 ± 6 days) compared with that (31 ± 8 days) of the CMA plus WES. Our results showed that WGS provided more comprehensive and precise genetic information with a rapid TAT and less DNA required than CMA plus WES, which enables it as an alternative prenatal diagnosis test for fetal structural anomalies.

Complex genomic alterations and intellectual disability: an interpretative challenge

BACKGROUND

Complex chromosomal rearrangements (CCRs) are structural rearrangements involving more than three chromosomes or having more than two breaks; approximately 70% are not associated with any clinical phenotype. Here, we describe a CCR segregating in a two-generation family.

METHOD

A 4-year-old male was evaluated for developmental delay, mild intellectual disability and epicanthus. Karyotype, fluorescence in situ hybridisation (FISH) analysis and array comparative genomic hybridisation (aCGH) analysis were performed on the patient and of all family members.

RESULT

Array CGH analysis of the proband detected two non-contiguous genomic gains of chromosome 2 at bands q32.3q33.2 and bands q36.1q36.3. Both karyotype and FISH analysis revealed a recombinant chromosome 2 with a direct insertion of regions q32.3q33.2 and q36.1q36.3 into region q12. Both of these regions were also present in their original location. Karyotype and FISH analysis of the father revealed a de novo direct insertion of regions q32.3q33.2 and q36.1q36.3 into region q12. Moreover, a de novo balanced translocation involving the q arm of the same chromosome 2 and the p arm of chromosome 10 was observed in the father of the proband. The single nucleotide polymorphism (SNP) array analysis and haplotype reconstruction confirmed the paternal origin of the duplications. Karyotype, FISH analysis and array CGH analysis of other family members were all normal.

CONCLUSION

This report underlines the importance of using different methods to correctly evaluate the origin and the structure of CCRs in order to provide an appropriate management of the patients and a good estimation of the reproductive risk of the family.

Hybrid sequencing resolves two germline ultra-complex chromosomal rearrangements consisting of 137 breakpoint junctions in a single carrier

Chromoanagenesis is a genomic event responsible for the formation of complex structural chromosomal rearrangements (CCRs). Germline chromoanagenesis is rare and the majority of reported cases are associated with an affected phenotype. Here, we report a healthy female carrying two de novo CCRs involving chromosomes 4, 19, 21 and X and chromosomes 7 and 11, respectively, with a total of 137 breakpoint junctions (BPJs). We characterized the CCRs using a hybrid-sequencing approach, combining short-read sequencing, nanopore sequencing, and optical mapping. The results were validated using multiple cytogenetic methods, including fluorescence in situ hybridization, spectral karyotyping, and Sanger sequencing. We identified 137 BPJs, which to our knowledge is the highest number of reported breakpoint junctions in germline chromoanagenesis. We also performed a statistical assessment of the positioning of the breakpoints, revealing a significant enrichment of BPJ-affecting genes (96 intragenic BPJs, 26 genes, p < 0.0001), indicating that the CCRs formed during active transcription of these genes. In addition, we find that the DNA fragments are unevenly and non-randomly distributed across the derivative chromosomes indicating a multistep process of scattering and re-joining of DNA fragments. In summary, we report a new maximum number of BPJs (137) in germline chromoanagenesis. We also show that a hybrid sequencing approach is necessary for the correct characterization of complex CCRs. Through in-depth statistical assessment, it was found that the CCRs most likely was formed through an event resembling chromoplexy—a catastrophic event caused by erroneous transcription factor binding.

Integrative analysis of structural variations using short-reads and linked-reads yields highly specific and sensitive predictions

Genetic diseases are driven by aberrations of the human genome. Identification of such aberrations including structural variations (SVs) is key to our understanding. Conventional short-reads whole genome sequencing (cWGS) can identify SVs to base-pair resolution, but utilizes only short-range information and suffers from high false discovery rate (FDR). Linked-reads sequencing (10XWGS) utilizes long-range information by linkage of short-reads originating from the same large DNA molecule. This can mitigate alignment-based artefacts especially in repetitive regions and should enable better prediction of SVs. However, an unbiased evaluation of this technology is not available. In this study, we performed a comprehensive analysis of different types and sizes of SVs predicted by both the technologies and validated with an independent PCR based approach. The SVs commonly identified by both the technologies were highly specific, while validation rate dropped for uncommon events. A particularly high FDR was observed for SVs only found by 10XWGS. To improve FDR and sensitivity, statistical models for both the technologies were trained. Using our approach, we characterized SVs from the MCF7 cell line and a primary breast cancer tumor with high precision. This approach improves SV prediction and can therefore help in understanding the underlying genetics in various diseases.

Cancer and many other diseases are often driven by structural rearrangements in the patients. Their precise identification is necessary to understand evolution and cure for the disease. In this study, we have compared two sequencing technologies for the identification of structural variations i.e. Illumina’s short-reads and 10X Genomics linked-reads sequencing. Short-reads sequencing is already known to have high false discovery rate for structural variations, while, an unbiased performance evaluation of linked-reads sequencing is missing. Hence, we evaluate the performance of these two technologies using computational and PCR based methodologies. Moreover, we also present a statistical approach to increase their performance, supporting better detection of structural variations and thus further research into disease biology.

Cytogenetically visible inversions are formed by multiple molecular mechanisms

Cytogenetically detected inversions are generally assumed to be copy number and phenotypically neutral events. While nonallelic homologous recombination is thought to play a major role, recent data suggest the involvement of other molecular mechanisms in inversion formation. Using a combination of short-read whole-genome sequencing (WGS), 10X Genomics Chromium WGS, droplet digital polymerase chain reaction and array comparative genomic hybridization we investigated the genomic structure of 18 large unique cytogenetically detected chromosomal inversions and achieved nucleotide resolution of at least one chromosomal inversion junction for 13/18 (72%). Surprisingly, we observed that seemingly copy number neutral inversions can be accompanied by a copy-number gain of up to 350 kb and local genomic complexities (3/18, 17%). In the resolved inversions, the mutational signatures are consistent with nonhomologous end-joining (8/13, 62%) or microhomology-mediated break-induced replication (5/13, 38%). Our study indicates that short-read 30x coverage WGS can detect a substantial fraction of chromosomal inversions. Moreover, replication-based mechanisms are responsible for approximately 38% of those events leading to a significant proportion of inversions that are actually accompanied by additional copy-number variation potentially contributing to the overall phenotypic presentation of those patients.

A pipeline for complete characterization of complex germline rearrangements from long DNA reads

BACKGROUND

Many genetic/genomic disorders are caused by genomic rearrangements. Standard methods can often characterize these variations only partly, e.g., copy number changes or breakpoints. It is important to fully understand the order and orientation of rearranged fragments, with precise breakpoints, to know the pathogenicity of the rearrangements.

METHODS

We performed whole-genome-coverage nanopore sequencing of long DNA reads from four patients with chromosomal translocations. We identified rearrangements relative to a reference human genome, subtracted rearrangements shared by any of 33 control individuals, and determined the order and orientation of rearranged fragments, with our newly developed analysis pipeline.

RESULTS

We describe the full characterization of complex chromosomal rearrangements, by filtering out genomic rearrangements seen in controls without the same disease, reducing the number of loci per patient from a few thousand to a few dozen. Breakpoint detection was very accurate; we usually see ~ 0 ± 1 base difference from Sanger sequencing-confirmed breakpoints. For one patient with two reciprocal chromosomal translocations, we find that the translocation points have complex rearrangements of multiple DNA fragments involving 5 chromosomes, which we could order and orient by an automatic algorithm, thereby fully reconstructing the rearrangement. A rearrangement is more than the sum of its parts: some properties, such as sequence loss, can be inferred only after reconstructing the whole rearrangement. In this patient, the rearrangements were evidently caused by shattering of the chromosomes into multiple fragments, which rejoined in a different order and orientation with loss of some fragments.

CONCLUSIONS

We developed an effective analytic pipeline to find chromosomal aberration in congenital diseases by filtering benign changes, only from long read sequencing. Our algorithm for reconstruction of complex rearrangements is useful to interpret rearrangements with many breakpoints, e.g., chromothripsis. Our approach promises to fully characterize many congenital germline rearrangements, provided they do not involve poorly understood loci such as centromeric repeats.

Familial Infertility (Azoospermia and Cryptozoospermia) in Two Brothers-Carriers of t(1;7) Complex Chromosomal Rearrangement (CCR): Molecular Cytogenetic Analysis

Structural aberrations involving more than two breakpoints on two or more chromosomes are known as complex chromosomal rearrangements (CCRs). They can reduce fertility through gametogenesis arrest developed due to disrupted chromosomal pairing in the pachytene stage. We present a familial case of two infertile brothers (with azoospermia and cryptozoospermia) and their mother, carriers of an exceptional type of CCR involving chromosomes 1 and 7 and three breakpoints. The aim was to identify whether meiotic disruption was caused by CCR and/or genomic mutations. Additionally, we performed a literature survey for male CCR carriers with reproductive failures. The characterization of the CCR chromosomes and potential genomic aberrations was performed using: G-banding using trypsin and Giemsa staining (GTG banding), fluorescent in situ hybridization (FISH) (including multicolor FISH (mFISH) and bacterial artificial chromosome (BAC)-FISH), and genome-wide array comparative genomic hybridization (aCGH). The CCR description was established as: der(1)(1qter->1q42.3::1p21->1q42.3::7p14.3->7pter), der(7)(1pter->1p2 1::7p14.3->7qter). aCGH revealed three rare genes variants: ASMT, GARNL3, and SESTD1, which were ruled out due to unlikely biological functions. The aCGH analysis of three breakpoint CCR regions did not reveal copy number variations (CNVs) with biologically plausible genes. Synaptonemal complex evaluation (brother-1; spermatocytes II/oligobiopsy; the silver staining technique) showed incomplete conjugation of the chromosomes. Associations between CCR and the sex chromosomes (by FISH) were not found. A meiotic segregation pattern (brother-2; ejaculated spermatozoa; FISH) revealed 29.21% genetically normal/balanced spermatozoa. The aCGH analysis could not detect smaller intergenic CNVs of few kb or smaller (indels of single exons or few nucleotides). Since chromosomal aberrations frequently do not affect the phenotype of the carrier, in contrast to the negative influence on spermatogenesis, there is an obvious need for genomic sequencing to investigate the point mutations that may be responsible for the differences between the azoospermic and cryptozoospermic phenotypes observed in a family. Progeny from the same parents provide a unique opportunity to discover a novel genomic background of male infertility.

Catastrophic Endgames: Emerging Mechanisms of Telomere-Driven Genomic Instability

When cells progress to malignancy, they must overcome a final telomere-mediated proliferative lifespan barrier called replicative crisis. Crisis is characterized by extensive telomere fusion that drives widespread genomic instability, mitotic arrest, hyperactivation of autophagy, and cell death. Recently, it has become apparent that that the resolution of dicentric chromosomes, which arise from telomere fusions during crisis, can initiate a sequence of events that leads to chromothripsis, a form of extreme genomic catastrophe. Chromothripsis is characterized by localized genomic regions containing tens to thousands of rearrangements and it is becoming increasingly apparent that chromothripsis occurs widely across tumor types and has a clinical impact. Here we discuss how telomere dysfunction can initiate genomic complexity and the emerging mechanisms of chromothripsis.

Structural variation in the sequencing era

Identifying structural variation (SV) is essential for genome interpretation but has been historically difficult due to limitations inherent to available genome technologies. Detection methods that use ensemble algorithms and emerging sequencing technologies have enabled the discovery of thousands of SVs, uncovering information about their ubiquity, relationship to disease and possible effects on biological mechanisms. Given the variability in SV type and size, along with unique detection biases of emerging genomic platforms, multiplatform discovery is necessary to resolve the full spectrum of variation. Here, we review modern approaches for investigating SVs and proffer that, moving forwards, studies integrating biological information with detection will be necessary to comprehensively understand the impact of SV in the human genome.

Chromothripsis and DNA Repair Disorders

Chromothripsis is a mutational mechanism leading to complex and relatively clustered chromosomal rearrangements, resulting in diverse phenotypic outcomes depending on the involved genomic landscapes. It may occur both in the germ and the somatic cells, resulting in congenital and developmental disorders and cancer, respectively. Asymptomatic individuals may be carriers of chromotriptic rearrangements and experience recurrent reproductive failures when two or more chromosomes are involved. Several mechanisms are postulated to underlie chromothripsis. The most attractive hypothesis involves chromosome pulverization in micronuclei, followed by the incorrect reassembly of fragments through DNA repair to explain the clustered nature of the observed complex rearrangements. Moreover, exogenous or endogenous DNA damage induction and dicentric bridge formation may be involved. Chromosome instability is commonly observed in the cells of patients with DNA repair disorders, such as ataxia telangiectasia, Nijmegen breakage syndrome, and Bloom syndrome. In addition, germline variations of TP53 have been associated with chromothripsis in sonic hedgehog medulloblastoma and acute myeloid leukemia. In the present review, we focus on the underlying mechanisms of chromothripsis and the involvement of defective DNA repair genes, resulting in chromosome instability and chromothripsis-like rearrangements.

Refined detection and phasing of structural aberrations in pediatric acute lymphoblastic leukemia by linked-read whole-genome sequencing

Structural chromosomal rearrangements that can lead to in-frame gene-fusions are a leading source of information for diagnosis, risk stratification, and prognosis in pediatric acute lymphoblastic leukemia (ALL). Traditional methods such as karyotyping and FISH struggle to accurately identify and phase such large-scale chromosomal aberrations in ALL genomes. We therefore evaluated linked-read WGS for detecting chromosomal rearrangements in primary samples of from 12 patients diagnosed with ALL. We assessed the effect of input DNA quality on phased haplotype block size and the detectability of copy number aberrations and structural variants in the ALL genomes. We found that biobanked DNA isolated by standard column-based extraction methods was sufficient to detect chromosomal rearrangements even at low 10x sequencing coverage. Linked-read WGS enabled precise, allele-specific, digital karyotyping at a base-pair resolution for a wide range of structural variants including complex rearrangements and aneuploidy assessment. With use of haplotype information from the linked-reads, we also identified previously unknown structural variants, such as a compound heterozygous deletion of ERG in a patient with the DUX4-IGH fusion gene. We conclude that linked-read WGS allows detection of important pathogenic variants in ALL genomes at a resolution beyond that of traditional karyotyping and FISH.

Challenges in identifying large germline structural variants for clinical use by long read sequencing

Genomic structural variations, previously considered rare events, are widely recognized as a major source of inter-individual variability and hence, a major hurdle in optimum patient stratification and disease management. Herein, we focus on large complex germline structural variations and present challenges towards target treatment via the synergy of state-of-the-art approaches and information technology tools. A complex structural variation detection remains challenging, as there is no gold standard for identifying such genomic variations with long reads, especially when the chromosomal rearrangement in question is a few Mb in length. A clinical case with a large complex chromosomal rearrangement serves as a paradigm. We feel that functional validation and data interpretation are of outmost importance for information growth to be translated into knowledge growth and hence, new working practices are highlighted.

Enzyme-free optical DNA mapping of the human genome using competitive binding

Optical DNA mapping (ODM) allows visualization of long-range sequence information along single DNA molecules. The data can for example be used for detecting long range structural variations, for aiding DNA sequence assembly of complex genomes and for mapping epigenetic marks and DNA damage across the genome. ODM traditionally utilizes sequence specific marks based on nicking enzymes, combined with a DNA stain, YOYO-1, for detection of the DNA contour. Here we use a competitive binding approach, based on YOYO-1 and netropsin, which highlights the contour of the DNA molecules, while simultaneously creating a continuous sequence specific pattern, based on the AT/GC variation along the detected molecule. We demonstrate and validate competitive-binding-based ODM using bacterial artificial chromosomes (BACs) derived from the human genome and then turn to DNA extracted from white blood cells. We generalize our findings with in-silico simulations that show that we can map a vast majority of the human genome. Finally, we demonstrate the possibility of combining competitive binding with enzymatic labeling by mapping DNA damage sites induced by the cytotoxic drug etoposide to the human genome. Overall, we demonstrate that competitive-binding-based ODM has the potential to be used both as a standalone assay for studies of the human genome, as well as in combination with enzymatic approaches, some of which are already commercialized.

From cytogenetics to cytogenomics: whole-genome sequencing as a first-line test comprehensively captures the diverse spectrum of disease-causing genetic variation underlying intellectual disability

Background

Since different types of genetic variants, from single nucleotide variants (SNVs) to large chromosomal rearrangements, underlie intellectual disability, we evaluated the use of whole-genome sequencing (WGS) rather than chromosomal microarray analysis (CMA) as a first-line genetic diagnostic test.

Methods

We analyzed three cohorts with short-read WGS: (i) a retrospective cohort with validated copy number variants (CNVs) (cohort 1, n = 68), (ii) individuals referred for monogenic multi-gene panels (cohort 2, n = 156), and (iii) 100 prospective, consecutive cases referred to our center for CMA (cohort 3). Bioinformatic tools developed include FindSV, SVDB, Rhocall, Rhoviz, and vcf2cytosure.

Results

First, we validated our structural variant (SV)-calling pipeline on cohort 1, consisting of three trisomies and 79 deletions and duplications with a median size of 850 kb (min 500 bp, max 155 Mb). All variants were detected. Second, we utilized the same pipeline in cohort 2 and analyzed with monogenic WGS panels, increasing the diagnostic yield to 8%. Next, cohort 3 was analyzed by both CMA and WGS. The WGS data was processed for large (> 10 kb) SVs genome-wide and for exonic SVs and SNVs in a panel of 887 genes linked to intellectual disability as well as genes matched to patient-specific Human Phenotype Ontology (HPO) phenotypes. This yielded a total of 25 pathogenic variants (SNVs or SVs), of which 12 were detected by CMA as well. We also applied short tandem repeat (STR) expansion detection and discovered one pathologic expansion in ATXN7. Finally, a case of Prader-Willi syndrome with uniparental disomy (UPD) was validated in the WGS data.

Important positional information was obtained in all cohorts. Remarkably, 7% of the analyzed cases harbored complex structural variants, as exemplified by a ring chromosome and two duplications found to be an insertional translocation and part of a cryptic unbalanced translocation, respectively.

Conclusion

The overall diagnostic rate of 27% was more than doubled compared to clinical microarray (12%). Using WGS, we detected a wide range of SVs with high accuracy. Since the WGS data also allowed for analysis of SNVs, UPD, and STRs, it represents a powerful comprehensive genetic test in a clinical diagnostic laboratory setting.

Toward Clinical Implementation of Next-Generation Sequencing-Based Genetic Testing in Rare Diseases: Where Are We?

Next-generation sequencing (NGS) technologies have changed the landscape of genetic testing in rare diseases. However, the rapid evolution of NGS technologies has outpaced its clinical adoption. Here, we re-evaluate the critical steps in the clinical application of NGS-based genetic testing from an informatics perspective. We suggest a 'fit-for-purpose' triage of current NGS technologies. We also point out potential shortcomings in the clinical management of genetic variants and offer ideas for potential improvement. We specifically emphasize the importance of ensuring the accuracy and reproducibility of NGS-based genetic testing in the context of rare disease diagnosis. We highlight the role of artificial intelligence (AI) in enhancing understanding and prioritization of variance in the clinical setting and propose deep learning frameworks for further investigation.