Complex structural variants in Mendelian disorders: identification and breakpoint resolution using short- and long-read genome sequencing
Genome Medicinevolume 10, Article number: 95 (2018)
Studies have shown that complex structural variants (cxSVs) contribute to human genomic variation and can cause Mendelian disease. We aimed to identify cxSVs relevant to Mendelian disease using short-read whole-genome sequencing (WGS), resolve the precise variant configuration and investigate possible mechanisms of cxSV formation.
We performed short-read WGS and analysis of breakpoint junctions to identify cxSVs in a cohort of 1324 undiagnosed rare disease patients. Long-read WGS and gene expression analysis were used to resolve one case.
We identified three pathogenic cxSVs: a de novo duplication-inversion-inversion-deletion affecting ARID1B, a de novo deletion-inversion-duplication affecting HNRNPU and a homozygous deletion-inversion-deletion affecting CEP78. Additionally, a de novo duplication-inversion-duplication overlapping CDKL5 was resolved by long-read WGS demonstrating the presence of both a disrupted and an intact copy of CDKL5 on the same allele, and gene expression analysis showed both parental alleles of CDKL5 were expressed. Breakpoint analysis in all the cxSVs revealed both microhomology and longer repetitive elements.
Our results corroborate that cxSVs cause Mendelian disease, and we recommend their consideration during clinical investigations. We show that resolution of breakpoints can be critical to interpret pathogenicity and present evidence of replication-based mechanisms in cxSV formation.
Structural variants (SVs) are a major source of variation in the human genome and collectively account for more differences between individuals than single nucleotide variants (SNVs) [1, 2]. SVs are categorised as canonical or complex . The canonical forms can be balanced or unbalanced and comprise inversions, insertions, translocations, deletions and duplications. More complex rearrangements are typically composed of three or more breakpoint junctions and cannot be characterised as a single canonical SV type. These are known as non-canonical or complex SVs (cxSVs) [3, 4].
Several previous studies have reported clinically relevant cxSVs in individuals with Mendelian disorders. For example, a duplication-triplication-inversion-duplication was found at the MECP2 and PLP1 loci in individuals with MECP2 duplication syndrome or Lubs syndrome (MIM: 300260) and Pelizaeus-Merzbacher disease (MIM: 312080) [5, 6], and a duplication-inversion-terminal deletion of chromosome 13 was present in foetuses with 13q deletion syndrome , among others [8,9,10]. Recently, pathogenic cxSVs associated with autism spectrum disorder and neuropsychiatric disorders have also been reported [11, 12]. Whole-genome sequencing (WGS) studies have shown that cxSVs are considerably more abundant and diverse than had been previously appreciated, representing an estimated 2% of the SVs in the human genome, and each human genome contains on average 14 cxSVs . The presence of multiple types of cxSVs has also been independently observed in several other studies [5, 12,13,14]. Extreme cases of cxSVs, such as chromothripsis, have also been identified in both cancer cells and the germline and involve hundreds of rearrangements often concerning more than one chromosome [11, 15].
Nevertheless, cxSVs are not typically considered during genomic analysis, largely due to technical challenges of identification. Complex SVs have been reported in projects such as the 1000 Genomes, but these primarily focused on the canonical types [1, 16, 17]. With the rapid expansion of high-throughput sequencing technologies including long-read WGS, genome-wide characterisation of SVs with high precision has been achieved , facilitating the study of more complex forms of SVs.
Therefore, in the present study, we sought to investigate the role of cxSVs in Mendelian disease by first identifying potentially clinically relevant cxSVs in a subset of the NIHR BioResource project using short-read WGS, second, resolving the variant configuration to base pair level resolution and, third, investigating possible mechanisms of cxSV formation by breakpoint analysis.
This cohort comprises 1324 individuals from the NIHR BioResource research study, which performs WGS of individuals with undiagnosed rare disorders. It is composed of three different subprojects: 725 were in the Inherited Retinal Disorders (IRD) project, 472 were in the Neurological and Developmental Disorders (NDD) project and 127 were in the Next Generation Children (NGC) project, which performs diagnostic trio WGS of individuals from Neonatal and Paediatric Intensive Care Units.
Short-read WGS and variant identification
We performed short-read WGS and excluded the possibility of pathogenic SNVs or indels, as part of the NIHR BioResource project as previously described . For the NDD and IRD subprojects we restricted SNVs and indel analysis to known disease-associated genes, which we assembled from sources including OMIM, RetNet and literature searches, then curated to ensure they comply with previously described criteria . The lists comprise 1423 genes (NDD) and 248 genes (IRD). For NGC participants, trio analysis focused on de novo and rare biallelic variant discovery unrestricted by a gene list.
The first stage of cxSV identification was calling and filtering canonical SVs. These initial calls comprise simple canonical SVs, and those which were potentially individual segments of cxSVs were then identified by clustering. These canonical SVs were called by Canvas , which identifies copy number gains and losses based on read depth, and Manta , which calls translocations, deletions, tandem duplications, insertions and inversions, and is based on both paired read fragment spanning and split read evidence. SVs were initially filtered to keep only those that pass standard Illumina quality filters, do not overlap previously reported CNVs in healthy cohorts  and are rare (minor allele frequency < 0.01) in the whole NIHR BioResource study (n = 9453) Schematic of the workflow can be found in Additional file 1: Figure S1.
Identification of potentially clinically relevant cxSVs
To identify potentially clinically relevant cxSVs we first identified putative cxSVs in the 1324 individuals by clustering canonical SV calls from Canvas and Manta using Bedtools cluster allowing a maximum distance between calls of 1 Kb . We then categorised the putative cxSVs into different subtypes previously described .
Next, as the scope of this study was limited to the identification of potentially clinically relevant cxSVs, we performed strict post-processing on the list of putative cxSVs. We excluded any for which visual inspection of the reads in IGV suggested the cxSV was not real but an artefact of a region of low sequencing quality and any in a gene that was not consistent with both the expected genotype and phenotype of the patient. We also filtered out retrotransposons, which are miscalled as multiple clustered intronic deletions, and dispersed duplications, which are frequently miscalled as an overlapping deletion and tandem duplication .
Sanger sequencing of the PCR product of the breakpoints was performed using standard protocols. Copy number variable segments of cxSVs and regions of homozygosity were confirmed using Illumina SNP genotyping array as previously described , and/or CytoScan® 750 K Cytogenetics Solution microarray (Affymetrix).
To resolve the configuration of the cxSV in participant 4 (P4), we performed long-read WGS with Oxford Nanopore Technologies (ONT). The sample was prepared using the 1D ligation library prep kit (SQK-LSK108), and genomic libraries were sequenced on R9 flowcell. Read sequences were extracted from base-called FAST5 files by albacore (version 2.0.2) to generate FASTQ files and then aligned against the GRCh37/hg19 human reference genome using NGMLR (version 0.2.6)  and LAST (version 912) , in order to compare results. Analysis was performed using default parameters, and for LAST, we used first last-train function to optimise alignment scoring. Variant calling was performed with Sniffles  and NanoSV , respectively.
RNA gene expression analysis of CDKL5 was as also performed on P4 and both parents. RNA was extracted from blood using the PAXgene Blood RNA Kit (QIAGEN) and retro-transcribed using the High Capacity cDNA Reverse Transcription Kit (Thermo Fisher Scientific). We performed PCR amplification and Sanger sequencing of the informative SNP rs35478150 (X:g.18638082A>C).
Breakpoint flanking sequence analysis
We analysed the sequence flanking each confirmed breakpoint (+/− 150 bps) and manually identified microhomology. The percentage of repetitive sequence was then calculated using RepeatMasker version open-4.0.7 (http://www.repeatmasker.org). In order to identify de novo SNVs and indels at the breakpoint junctions, SNV/indel calling and de novo filtering was performed with Platypus (http://github.com/andyrimmer/Platypus) for those participants for which parental WGS was available (P1 and P4).
Potentially clinically relevant cxSVs in four individuals with Mendelian disease
We identified four individuals with potentially clinically relevant cxSVs. Participant 1 (P1) presents a de novo duplication-inversion-inversion-deletion encompassing ARID1B (MIM: 135900) that causes Coffin-Siris syndrome (CSS [MIM: 135900]). This individual was a 4-month-old female who was born prematurely and presented with characteristic features of CSS as a neonate. CSS is a multiple malformation syndrome characterised by intellectual disability, severe speech impairment, coarse facial features, microcephaly, developmental delay and hypoplastic nails of the fifth digits .
A large cxSV was identified on chromosome 6, comprising a 3.3 Mb duplication, two inversions of 4.9 Kb and 3.3 Mb, and a 16.3 Mb deletion (Fig. 1a; Table 1). A total of 87 protein-coding genes were within the structural variant boundaries (Additional file 2: Table S1), of which 21 have been previously described as disease-associated in OMIM. The 16.3 Mb deletion contains 72 genes, of which only 6 have been reported as associated with autosomal dominant disease or constrained for loss-of-function (LOF) variation in ExAC  (Additional file 2: Table S1). Of these 6, only ARID1B has been previously reported as disease-associated with a LOF mechanism. Haploinsufficiency of ARID1B causes CSS, consistent with the phenotype of P1. We also looked at the 10 autosomal recessive genes within the deletion and did not find a second likely pathogenic variant in any. No disease-associated gene that was present within the duplicated region had been reported to be triplosensitive. Furthermore, the first inversion and the 3′ breakpoint of the second inversion were within CNKSR3 (MIM: 617476). However, CNKSR3 has not previously been associated with disease and is not constrained for LOF variation in ExAC; thus, the effect of this inversion on the phenotype remains unknown.
Although the LOF of ARID1B likely explains the phenotype of this individual, it is possible that other genes affected by the cxSV might contribute to the phenotype. Examination of the parental origin of the hemizygous variants in the deleted region confirmed that the cxSV occurred on the paternal chromosome, consistent with previously reported observations that ~ 80% of de novo mutations are of paternal origin .
Participant 2 (P2) has a de novo deletion-inversion-duplication encompassing HNRNPU (MIM: 602869). This individual is a 22-year-old male who presented at term with hypotonia. All his early developmental milestones were delayed, and he presented with tonic-clonic seizures at 9 months. His seizure disorder has been managed by medication but has continued episodically into adulthood. He also has significant intellectual disability, autism, and limited speech and language, and MRI showed partial agenesis of the corpus callosum and enlarged ventricles.
We identified a cxSV on chromosome 1, formed by a 1.2 Mb deletion and a 246 Kb duplication flanking an inversion of 505 Kb (Fig. 1b; Table 1). This variant encompassed eight genes (Additional file 2: Table S1), of which two were previously associated with disease: COX20 (MIM: 614698) and HNRNPU, both within the deletion boundaries. Haploinsufficiency of COX20 was not deemed likely to be pathogenic as variants in this gene have an autosomal recessive mode of inheritance and result in a mitochondrial complex IV deficiency (MIM: 220110), which is not consistent with the individual’s phenotype, and no second rare variant was identified. However, HNRNPU is a highly constrained gene for LOF variants, in which haploinsufficiency causes early infantile epileptic encephalopathy (EIEE [MIM: 617391]). Microdeletions of HNRNPU have been described in individuals with intellectual disability and other clinical features, such as seizures, corpus callosum abnormalities and microcephaly .
Participant 3 (P3), a 66-year-old male, presented with a cone-rod dystrophy and hearing loss due to a homozygous deletion-inversion-deletion overlapping CEP78 (MIM: 617110). Onset was in his fifth decade with central vision loss, photophobia and nystagmus accompanied by progressive hearing impairment, following a severe influenza-like viral infection. Two homozygous deletions in chromosome 9 of nearly 6 and 10 Kb were found flanking an inversion of 298 bp (Fig. 1c; Table 1). The second deletion intersects with the first 5 exons of CEP78. Biallelic LOF variants in this gene have been previously shown to cone-rod dystrophy and hearing loss (MIM: 617236) . Although we could not perform segregation analysis due to lack of parental DNA, it was observed to be within a copy number neutral region of homozygosity that comprised approximately Chr9:70984372-86933884.
Participant 4 (P4) presents a duplication-inversion-duplication overlapping CDKL5 (MIM: 300203) on chromosome X. This individual was a female term (41 + 1) neonate who presented with foetal bradycardia. She was diagnosed with hypoxic-ischemic encephalopathy grade 2, intrauterine hypoxia, and perinatal asphyxia, with poor cord gases. Hypothermia was induced after birth for 72 h to reduce brain injury. WGS revealed a de novo duplication-inversion-duplication, with the respective sizes of 280 Kb, 458 Kb and 283 Kb (Fig. 1d; Table 1). The inversion 3′ breakpoint is in intron 3 of 20 of CDKL5 (NM_003159). Heterozygous rare variants in X-linked CDKL5 in females cause EIEE, severe intellectual disability and Rett-like features (MIM: 300672). There are three other genes within the boundaries of this cxSV, none of them disease-associated in OMIM (Additional file 2: Table S1).
Accurate resolution of variant configuration is necessary for interpretation of pathogenicity
For each of the four cases, we validated each breakpoint in order to resolve the variant configuration to base pair level resolution. For P1, P2 and P3, all novel junctions were confirmed by Sanger sequencing, and all copy number changes were confirmed by microarrays that were performed concurrently with the WGS (Fig. 1a–c). No alternative pathogenic SNVs, indels or canonical SVs were identified, and the cxSVs were absent in our internal cohort of 9453 genomes, ClinVar or DECIPHER. The cxSVs in P1, P2 and P3 were therefore classified as pathogenic according to the ACMG guidelines .
Resolving the configuration of the cxSV in P4 was more challenging because the SV calls from short-read WGS were consistent with multiple possible configurations (Additional file 1: Figure S2). Importantly, in two of the possible configurations, there is an intact copy of CDKL5 on the non-reference allele, in addition to the disrupted copy (Additional file 1: Figure S2A–B), whereas in others there is no intact copy of CDKL5 (Additional file 1: Figure S2C–D). Therefore, resolving the configuration was essential for the interpretation of the pathogenicity of this variant. We attempted PCR amplification over the predicted new formed breakpoint junctions and could only amplify one supporting the disrupted CDKL5, due to repetitive sequence around the other breakpoints. Both duplications were confirmed by microarray.
In order to resolve the configuration, we performed long-read WGS of P4 using ONT. We obtained a median read length of 8136 bp (Additional file 1: Figure S3A), 56% of the genome was covered with a minimum coverage of 3x (Additional file 1: Figure S3B), and around 97% of the reads mapped to the human genome (GRCh37/hg19). All the breakpoints of the cxSV were covered by at least four reads. Coverage was insufficient to resolve the cxSV using long-read SV calling algorithms such as Sniffles  or NanoSV  (for which a minimum coverage of 10x is recommended). In lieu of this, we manually reviewed the split long reads across the cxSV junction breakpoints. Eight of the reads that covered the cxSV breakpoints were identified as inherited from the paternal chromosome, either by SNP phasing (Fig. 1d, J2, J3, J4 and J6) or by indirect phasing based on the assumption that breakpoint junctions occur on the same allele (Fig. 1d, J5). Therefore, ONT sequencing allowed us to identify two reads supporting the junction that was initially not possible to confirm by Sanger sequencing (J5) due to repetitive sequences. By phasing analysis, we were also able to identify three reads supporting an intact copy of CDKL5 in the allele inherited from the father (Fig. 1d, J6), confirming that the cxSV harbours an intact copy of CDKL5. Two possible configurations remain (Additional file 1: Figure S2A–B), both of which have been proposed previously [12, 34]. These are indistinguishable by short-read sequencing technology because the breakpoint junctions are identical, or even by long reads unless all junctions are crossed in the same molecule. Only one of these possible configurations is represented in Fig. 1d and Table 1 for clarity.
We performed RNA expression analysis (Sanger sequencing of one informative SNP using cDNA) and demonstrated biparental allele expression of CDKL5 in the child (Additional file 1: Figure S4). This further supports the presence of an intact copy of CDKL5 on the paternal allele and suggests that regulation of CDKL5 is probably not perturbed by the nearby cxSV. This variant was classified as VUS. The child is currently 1 year old and developmentally normal with no seizures, but remains under ongoing follow-up.
Microhomology and repetitive elements occur in conjunction at cxSV breakpoint junctions
Mutational signatures around novel breakpoint junctions of SVs can yield insights into the mechanisms by which they were formed. Therefore, we analysed the sequences of all of the novel breakpoint junctions. It has previously been reported that DNA replication-based mechanisms such as microhomology-mediated break-induced replication (MMBIR) or fork stalling and template switching (FoSTeS) are likely to be the primary mechanism responsible for the formation of cxSVs [3, 4, 35,36,37]. Our data overall support this as there is microhomology of at least 3 bp in all of the eight novel breakpoint junctions in the four individuals (Additional file 1: Figure S5 and Additional file 2: Table S2). We also observe in P2 the insertion of two sequences of 5 and 48 bp in J1 and J3 junctions, and the insertion in P4 of a 100 bp Alu sequence in J2 junction. It has been previously suggested that Alu elements could facilitate template switching and annealing via homology between replication forks .
Additional evaluation of the breakpoint sequences with RepeatMasker also identified longer repetitive elements in all of the individuals (Table 2 and Additional file 1: Figure S5). In P1, we found that sequence flanking two of the breakpoints had high similarity to SINE sequences (ERVL-MaLRs), one with LINE sequences (L2) and one with DNA/hAT-Charlie (MER3) sequences (Table 1); in P2, we noted that sequence flanking three of the breakpoints had similarity to SINE sequences (Alu and MIR); in P3, sequences surrounding all the breakpoints presented high similarity to LINEs; and in P4, one of the breakpoints had similarity to SINE/Alu sequences.
A recent study showed that a high proportion of Alu-mediated SVs contain a hybrid Alu element in the derivative chromosome . However, we do not observe these in our study. Studies have also shown that due to the error-prone nature of replication-based mechanisms of cxSV formation, de novo SNVs and indels can occur concomitantly to cxSVs . Thus, we looked for de novo SNVs or indels in 1 Kb regions around each novel breakpoint junction in those individuals for which parental WGS data was available (P1 and P4), and there were none. For P2 and P3, we considered all rare SNVs and indels in those regions and did not identify any.
In the present study, we aimed to identify cxSVs relevant to Mendelian disease using short-read WGS, to resolve the precise variant configurations and to investigate possible mechanisms of cxSV formation. We have presented three individuals with pathogenic cxSVs and one with an interesting cxSV of unknown significance. We showed that precise resolution of variant configuration can be essential for interpreting pathogenicity and presented evidence of both DNA replication based and homologous recombination mechanisms of formation.
Here, we highlight the role of cxSVs as a cause of Mendelian disease. However, cxSVs are not typically considered in analysis pipelines, in part due to the technical and analytical challenges around identification and interpretation, and when there is an associated deletion or duplication detected by microarray further analysis is rarely performed. Therefore, when not included in analyses, the full scope of genome-wide structural variation is overlooked. In this study, 0.2% (3/1324) of Mendelian disease cases were caused by cxSV. This compares to approximately 5–20% of individuals with Mendelian disorders who have a clinically relevant canonical SV [18, 39, 40]. However, 0.2% is likely to be an underestimate because short-read WGS has limited power to resolve cxSVs and because our workflow was designed to maximise specificity.
There are several technologies available for the identification of cxSVs, including short-read WGS, long-read WGS, long insert WGS (liWGS) and microarrays, each with strengths and limitations. This study demonstrates the potential of short-read WGS to identify clinically relevant cxSVs. An advantage of this is that it is a more commonly used technology than some of the other options and does not necessitate using a different technology specifically to identify cxSVs. The main limitation is that at repetitive regions mapping and variant calling algorithms have lower sensitivity. This is particularly problematic given that cxSVs are more likely to occur in repetitive regions. These limitations may now be addressed by long-read sequencing technologies such as Nanopore, either in combination with another technology as in this study or as a first line approach. These have the advantage of reads of 10–100 Kb allowing for more accurate mapping particularly over repetitive regions and facilitating phasing . Various other studies have already demonstrated the power of long-read WGS to detect SVs and cxSVs [25, 27, 41, 42]. However, the limitations are that coverage is lower and error rate is higher than short-read WGS. For these reasons, it is unlikely that, in our case, these variants would have been identified from long-read data alone without prior knowledge of the region of interest. A third possible method to identify complex rearrangements is liWGS, which has been successfully employed to detect cxSVs in other studies, and has the advantage over short-read WGS of improved mapping particularly over repetitive regions due to the large fragments, but has a lower resolution of ~ 5 Kb [11,12,13].
The segments of cxSVs that have copy number changes could in some cases be detected by microarray. These have the advantage of low cost, and that they are already the first-line test in many cases. However, they do not allow precise resolution of breakpoints and they would also miss inversions and those cxSVs where there is no CNV change. Importantly, a routinely pre-screening by microarray could potentially identify a high number of cxSVs. Actually, it has been reported that 7.6% of all rare duplications detected by microarray are part of a complex rearrangement . However, it is likely that many ‘canonical CNVs’ detected by microarray are actually misclassified cxSVs, in part due to the impossibility for detecting inversions by microarray, since inversions are involved in 84.8% of cxSVs . In this study, the copy number changes of all four cxSVs were confirmed by microarrays.
Our experience with P4, whose cxSV intersects CDKL5, demonstrates that understanding the precise configuration of a cxSV can be essential for interpreting the pathogenicity of the variant, especially if the gene of interest is disrupted by a duplication or inversion rather than a deletion. The impact of a deletion on the function of affected genes is generally assumed to be LOF. However, the consequence of a duplication can be uncertain and depends on precisely how the variant rearranges the gene, as well as gene-specific factors such as dosage sensitivity. Furthermore, duplications intersecting regulatory regions can result in a different phenotype from variants within the gene itself .
Investigating cxSVs in our cohort identified previously reported subclasses (delINVdup, delINVdel and dupINVdup in P2, P3 and P4, respectively), as well as a dupINVinvDEL in P1 . One of the limitations of our study is that certain subclasses of cxSVs such as chromothripsis and those formed by nested rather than chained breakpoints would be excluded by our filtering and clustering method. Our method was designed specifically to identify possibly pathogenic variants in Mendelian disease, and thus was optimised to maximise specificity, at the cost of sensitivity and scalability. Therefore, while a detailed, sensitive, large-scale assessment of the role of cxSVs in Mendelian disease would be valuable, this is beyond the scope of the present study. Even when a putative cxSVs is identified, sometimes it is not straightforward to define them because the distinction between canonical SVs, cxSVs and chromoanagenesis can be unclear [35, 38]. It is therefore perhaps appropriate to consider types of human genomic variation as a continuum rather than discrete classes, progressing from SNVs (that typically cause the least disruption to the genome), through indels, canonical SVs and cxSVs to the highly disruptive chromoanagenesis and aneuploidies.
The high frequency of microhomology observed at the breakpoint junctions of the cxSVs in our study and the presence of inserted sequence in three of them is consistent with the hypothesis that replication-based mechanisms such as FoSTeS/MMBIR are primarily responsible for the formation of cxSVs [3, 4, 35,36,37]. However, we also find longer repetitive elements including Alu elements in the vicinity of breakpoint junctions in all of the cxSVs, consistent with other studies [44, 45]. The exact role of repetitive elements in SV/cxSV formation is currently unclear. Repetitive elements have classically been seen as signatures of recombination-based mechanisms such as non-allelic homologous recombination, which is a well-known mechanism of formation of recurrent SVs [4, 44]. This could suggest that recombination-based and replication-based mutational mechanisms might together mediate the formation of non-recurrent cxSVs. However, it is more likely that repetitive elements facilitate replication-based SV/cxSV formation by, for example, providing the requisite microhomology islands or increasing the susceptibility of the region to the formation of secondary DNA structures that can cause replication fork collapse [46, 47].
Our work demonstrates that cxSVs contribute to rare Mendelian disorders, and provides insight into identifying and resolving both the conformation and the mechanism of formation of cxSVs by using short and long-read WGS. We demonstrate that understanding the precise configuration can be essential for interpreting the pathogenicity of cxSVs. We suggest that cxSVs should be included into research and clinical diagnosis and considered when screening SVs in the human genome. Further detailed characterisation of cxSVs in large-scale WGS studies will be essential for further unveiling the complex architecture of cxSVs and determining accurate population frequencies.
Complex structural variant
Early infantile epileptic encephalopathy
Fork stalling and template switching
Inherited retinal disorder
Long insert WGS
Microhomology-mediated break-induced replication
Neurological and developmental disorder
Oxford Nanopore Technologies
Single nucleotide variant
Sudmant PH, Rausch T, Gardner EJ, Handsaker RE, Abyzov A, Huddleston J, et al. An integrated map of structural variation in 2,504 human genomes. Nature. 2015;526(7571):75–81.
Conrad DF, Bird C, Blackburne B, Lindsay S, Mamanova L, Lee C, et al. Mutation spectrum revealed by breakpoint sequencing of human germline CNVs. Nat Genet. 2010;42(5):385–91.
Quinlan AR, Hall IM. Characterizing complex structural variation in germline and somatic genomes. Trends Genet. 2012;28(1):43–53.
Weckselblatt B, Rudd MK. Human structural variation: mechanisms of chromosome rearrangements. Trends Genet. 2015;31(10):587–99.
Carvalho CM, Ramocki MB, Pehlivan D, Franco LM, Gonzaga-Jauregui C, Fang P, et al. Inverted genomic segments and complex triplication rearrangements are mediated by inverted repeats in the human genome. Nat Genet. 2011;43(11):1074–81.
Beck CR, Carvalho CM, Banser L, Gambin T, Stubbolo D, Yuan B, et al. Complex genomic rearrangements at the PLP1 locus include triplication and quadruplication. PLoS Genet. 2015;11(3):e1005050.
Quelin C, Spaggiari E, Khung-Savatovsky S, Dupont C, Pasquier L, Loeuillet L, et al. Inversion duplication deletions involving the long arm of chromosome 13: phenotypic description of additional three fetuses and genotype-phenotype correlation. Am J Med Genet A. 2014;164A(10):2504–9.
Arno G, Agrawal SA, Eblimit A, Bellingham J, Xu M, Wang F, et al. Mutations in REEP6 cause autosomal-recessive retinitis pigmentosa. Am J Hum Genet. 2016;99(6):1305–15.
Kuroda Y, Ohashi I, Saito T, Nagai J, Ida K, Naruto T, et al. Deletion of UBE3A in brothers with Angelman syndrome at the breakpoint with an inversion at 15q11.2. Am J Med Genet A. 2014;164A(11):2873–8.
Carvalho CM, Zhang F, Liu P, Patel A, Sahoo T, Bacino CA, et al. Complex rearrangements in patients with duplications of MECP2 can occur by fork stalling and template switching. Hum Mol Genet. 2009;18(12):2188–203.
Collins RL, Brand H, Redin CE, Hanscom C, Antolik C, Stone MR, et al. Defining the diverse spectrum of inversions, complex structural variation, and chromothripsis in the morbid human genome. Genome Biol. 2017;18(1):36.
Brand H, Collins RL, Hanscom C, Rosenfeld JA, Pillalamarri V, Stone MR, et al. Paired-duplication signatures mark cryptic inversions and other complex structural variation. Am J Hum Genet. 2015;97(1):170–6.
Brand H, Pillalamarri V, Collins RL, Eggert S, O'Dushlaine C, Braaten EB, et al. Cryptic and complex chromosomal aberrations in early-onset neuropsychiatric disorders. Am J Hum Genet. 2014;95(4):454–61.
Lohmann K, Redin C, Tonnies H, Bressman SB, Subero JIM, Wiegers K, et al. Complex and dynamic chromosomal rearrangements in a family with seemingly non-Mendelian inheritance of Dopa-responsive dystonia. JAMA Neurol. 2017;74(7):806–12.
Stephens PJ, Greenman CD, Fu B, Yang F, Bignell GR, Mudie LJ, et al. Massive genomic rearrangement acquired in a single catastrophic event during cancer development. Cell. 2011;144(1):27–40.
Hehir-Kwa JY, Marschall T, Kloosterman WP, Francioli LC, Baaijens JA, Dijkstra LJ, et al. A high-quality human reference panel reveals the complexity and distribution of genomic structural variants. Nat Commun. 2016;7:12989.
Genomes Project C, Abecasis GR, Auton A, Brooks LD, MA DP, Durbin RM, et al. An integrated map of genetic variation from 1,092 human genomes. Nature. 2012;491(7422):56–65.
Carss KJ, Arno G, Erwood M, Stephens J, Sanchis-Juan A, Hull S, et al. Comprehensive rare variant analysis via whole-genome sequencing to determine the molecular pathology of inherited retinal disease. Am J Hum Genet. 2017;100(1):75–90.
Wright CF, Fitzgerald TW, Jones WD, Clayton S, McRae JF, van Kogelenberg M, et al. Genetic diagnosis of developmental disorders in the DDD study: a scalable analysis of genome-wide research data. Lancet. 2015;385(9975):1305–14.
Roller E, Ivakhno S, Lee S, Royce T, Tanner S. Canvas: versatile and scalable detection of copy number variants. Bioinformatics. 2016;32(15):2375–7.
Chen X, Schulz-Trieglaff O, Shaw R, Barnes B, Schlesinger F, Kallberg M, et al. Manta: rapid detection of structural variants and indels for germline and cancer sequencing applications. Bioinformatics. 2016;32(8):1220–2.
Zarrei M, MacDonald JR, Merico D, Scherer SW. A copy number variation map of the human genome. Nat Rev Genet. 2015;16(3):172–83.
Quinlan AR, Hall IM. BEDTools: a flexible suite of utilities for comparing genomic features. Bioinformatics. 2010;26(6):841–2.
Zhao X, Emery SB, Myers B, Kidd JM, Mills RE. Resolving complex structural genomic rearrangements using a randomized approach. Genome Biol. 2016;17(1):126.
Sedlazeck FJ, Rescheneder P, Smolka M, Fang H, Nattestad M, von Haeseler A, et al. Accurate detection of complex structural variations using single-molecule sequencing. Nat Methods. 2018;15(6):461–8.
Kielbasa SM, Wan R, Sato K, Horton P, Frith MC. Adaptive seeds tame genomic sequence comparison. Genome Res. 2011;21(3):487–93.
Cretu Stancu M, van Roosmalen MJ, Renkens I, Nieboer MM, Middelkamp S, de Ligt J, et al. Mapping and phasing of structural variation in patient genomes using nanopore sequencing. Nat Commun. 2017;8(1):1326.
Wieczorek D, Bogershausen N, Beleggia F, Steiner-Haldenstatt S, Pohl E, Li Y, et al. A comprehensive molecular study on Coffin-Siris and Nicolaides-Baraitser syndromes identifies a broad molecular and clinical spectrum converging on altered chromatin remodeling. Hum Mol Genet. 2013;22(25):5121–35.
Lek M, Karczewski KJ, Minikel EV, Samocha KE, Banks E, Fennell T, et al. Analysis of protein-coding genetic variation in 60,706 humans. Nature. 2016;536(7616):285–91.
Goldmann JM, Wong WS, Pinelli M, Farrah T, Bodian D, Stittrich AB, et al. Parent-of-origin-specific signatures of de novo mutations. Nat Genet. 2016;48(8):935–9.
Bramswig NC, Ludecke HJ, Hamdan FF, Altmuller J, Beleggia F, Elcioglu NH, et al. Heterozygous HNRNPU variants cause early onset epilepsy and severe intellectual disability. Hum Genet. 2017;136(7):821–34.
Namburi P, Ratnapriya R, Khateb S, Lazar CH, Kinarty Y, Obolensky A, et al. Bi-allelic truncating mutations in CEP78, encoding centrosomal protein 78, cause cone-rod degeneration with sensorineural hearing loss. Am J Hum Genet. 2016;99(3):777–84.
Richards S, Aziz N, Bale S, Bick D, Das S, Gastier-Foster J, et al. Standards and guidelines for the interpretation of sequence variants: a joint consensus recommendation of the American College of Medical Genetics and Genomics and the Association for Molecular Pathology. Genet Med. 2015;17(5):405–24.
Gu S, Yuan B, Campbell IM, Beck CR, Carvalho CM, Nagamani SC, et al. Alu-mediated diverse and complex pathogenic copy-number variants within human chromosome 17 at p13.3. Hum Mol Genet. 2015;24(14):4061–77.
Liu P, Erez A, Nagamani SC, Dhar SU, Kolodziejska KE, Dharmadhikari AV, et al. Chromosome catastrophes involve replication mechanisms generating complex genomic rearrangements. Cell. 2011;146(6):889–903.
Lee JA, Carvalho CM, Lupski JR. A DNA replication mechanism for generating nonrecurrent rearrangements associated with genomic disorders. Cell. 2007;131(7):1235–47.
Zhang F, Khajavi M, Connolly AM, Towne CF, Batish SD, Lupski JR. The DNA replication FoSTeS/MMBIR mechanism can generate genomic, genic and exonic complex rearrangements in humans. Nat Genet. 2009;41(7):849–53.
Carvalho CM, Pehlivan D, Ramocki MB, Fang P, Alleva B, Franco LM, et al. Replicative mechanisms for CNV formation are error prone. Nat Genet. 2013;45(11):1319–26.
Cooper GM, Coe BP, Girirajan S, Rosenfeld JA, Vu TH, Baker C, et al. A copy number variation morbidity map of developmental delay. Nat Genet. 2011;43(9):838–46.
Miller DT, Adam MP, Aradhya S, Biesecker LG, Brothman AR, Carter NP, et al. Consensus statement: chromosomal microarray is a first-tier clinical diagnostic test for individuals with developmental disabilities or congenital anomalies. Am J Hum Genet. 2010;86(5):749–64.
Merker JD, Wenger AM, Sneddon T, Grove M, Zappala Z, Fresard L, et al. Long-read genome sequencing identifies causal structural variation in a Mendelian disease. Genet Med. 2018;20(1):159–63.
Gong L, Wong CH, Cheng WC, Tjong H, Menghi F, Ngan CY, et al. Picky comprehensively detects high-resolution structural variants in nanopore long reads. Nat Methods. 2018;15(6):455–60.
Flottmann R, Kragesteen BK, Geuer S, Socha M, Allou L, Sowinska-Seidler A, et al. Noncoding copy-number variations are associated with congenital limb malformation. Genet Med. 2018;20(6):599-607.
Stankiewicz P, Lupski JR. Genome architecture, rearrangements and genomic disorders. Trends Genet. 2002;18(2):74–82.
Gu S, Posey JE, Yuan B, Carvalho CM, Luk HM, Erikson K, et al. Mechanisms for the generation of two quadruplications associated with split-hand malformation. Hum Mutat. 2016;37(2):160–4.
Carvalho CM, Lupski JR. Mechanisms underlying structural variant formation in genomic disorders. Nat Rev Genet. 2016;17(4):224–38.
Bacolla A, Jaworski A, Larson JE, Jakupciak JP, Chuzhanova N, Abeysinghe SS, et al. Breakpoints of gross deletions coincide with non-B DNA conformations. Proc Natl Acad Sci U S A. 2004;101(39):14162–7.
Collins RL, Stone MR, Brand H, Glessner JT, Talkowski ME. CNView: a visualization and annotation tool for copy number variation from whole-genome sequencing. bioRxiv. 2016. https://github.com/RCollins13/CNView.
Ordulu Z, Wong KE, Currall BB, Ivanov AR, Pereira S, Althari S, et al. Describing sequencing results of structural chromosome rearrangements with a suggested next-generation cytogenetic nomenclature. Am J Hum Genet. 2014;94(5):695–709.
We thank the participants involved in this study and their families. We thank Dr. Ernest Turro for his part in initiating our collaboration with Oxford Nanopore Technologies and for the helpful comments on this manuscript. We also thank Mateus Patricio for his contribution to the breakpoint analysis.
This work was supported by the National Institute for Health Research England (NIHR) for the NIHR BioResource project (grant number RG65966). GA and ARW are supported by the NIHR Biomedical Research Centre at Moorfields Eye Hospital NHS Trust and UCL Institute of Ophthalmology, Moorfields Eye Charity, Fight for Sight (UK), Foundation Fighting Blindness and Retinitis Pigmentosa Fighting Blindness. GA is a recipient of a Fight for Sight (UK) Early Career Investigator Award.
Availability of data and materials
Sequence data for the 1324 individuals in this work have been deposited at the European Genome-Phenome Archive. IRD (EGAD00001002656) is available at the time of publication of this article. NGC (EGAD00001004357) and NDD(EGAD00001004456) will become available in 2019. In the meantime, the data will be available upon reasonable request, and sequencing data for just the four individuals presented in this work (P1, P2, P3 and P4) have also been made available in a separate release (EGAD00001004438). The four complex structural variants have been submitted to ClinVar (SCV000807695, SCV000807696, SCV000807697 and SCV000807698).
Ethics approval and consent to participate
All participants provided written informed consent to participate in the study. The study was approved by the East of England Cambridge South national institutional review board (13/EE/0325). The research conforms with the principles of the Declaration of Helsinki.
Consent for publication
All individuals involved in this study provided consent for publication.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Figure S1. cxSV analysis workflow. Figure S2. cxSV models for P4. Figure S3. Quality control results of the long-read WGS performed on P4. Figure S4. RNA gene expression of CDKL5 from P4 and both parents. Figure S5. Proposed mechanisms of cxSV formation and breakpoint junction alignments. (PDF 1966 kb)
Table S1. List of genes present in P1, P2, P3 and P4 cxSVs). Table S2. Complex SVs breakpoint details) (XLSX 32 kb)