<figure><img src="image_1.png" /></figure>

HAL
open science

Structural variant calling and clinical interpretation in
6224 unsolved rare disease exomes

German Demidov, Steven Laurie, Annalaura Torella, Giulio Piluso, Marcello Scala, Manuela Morleo, Vincenzo Nigro, Holm Graessner, Siddharth Banka, Alfons Macaya, et al.

<figure><img src="image_2.png" /></figure>

► To cite this version:

German Demidov, Steven Laurie, Annalaura Torella, Giulio Piluso, Marcello Scala, et al. Structural variant calling and clinical interpretation in 6224 unsolved rare disease exomes. European Journal of Human Genetics, 2024, 32 (8), pp.998 - 1004. 10.1038/s41431-024-01637-4. hal-05163100

HAL Id: hal-05163100

https://hal.science/hal-05163100v1

Submitted on 15 Jul 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_3.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
ARTICLE **OPEN**

<figure><img src="image_4.png" /></figure>

Check for updates

# Structural variant calling and clinical interpretation in 6224 unsolved rare disease exomes

German Demidov ID<fnref n="1" />,*, Steven Laurie ID², Annalaura Torella³,⁴, Giulio Piluso ID³, Marcello Scala ID<fnref n="5" />,⁶, Manuela Morleo³,⁴, Vincenzo Nigro³,⁴, Holm Graessner ID<fnref n="1" />,⁷, Siddharth Banka ID⁸,⁹, Solve-RD consortium*, Katja Lohmann ID¹⁰ and Stephan Ossowski ID<fnref n="1" />

© The Author(s) 2024, corrected publication 2025

Structural variants (SVs), including large deletions, duplications, inversions, translocations, and more complex events have the potential to disrupt gene function resulting in rare disease. Nevertheless, current pipelines and clinical decision support systems for exome sequencing (ES) tend to focus on small alterations such as single nucleotide variants (SNVs) and insertions-deletions shorter than 50 base pairs (indels). Additionally, detection and interpretation of large copy-number variants (CNVs) are frequently performed. However, detection of other types of SVs in ES data is hampered by the difficulty of identifying breakpoints in off-target (intergenic or intronic) regions, which makes robust identification of SVs challenging. In this paper, we demonstrate the utility of SV calling in ES resulting in a diagnostic yield of 0.4% (23 out of 5825 probands) for a large cohort of unsolved patients collected by the Solve-RD consortium. Remarkably, 8 out of 23 pathogenic SV were not found by comprehensive read-depth-based CNV analysis, resulting in a 0.13% increased diagnostic value.

European Journal of Human Genetics (2024) 32:998–1004; https://doi.org/10.1038/s41431-024-01637-4

## INTRODUCTION

Next-generation sequencing (NGS) is a method widely used for the detection of single-nucleotide variants (SNVs) and short indels [1] in the clinical diagnosis of rare genetic diseases. Detection of larger or complex variants with short read NGS is challenging as reads with a length of 100–150 bp cannot span across distantly located breakpoints. Nonetheless, methods for the detection of structural variants (SVs), including copy-number variants (CNVs) as well as copy-number neutral variants (such as inversions and balanced translocations) have been developed for genome sequencing (GS). Three main signals are used to detect SVs from genome data: (1) paired-end (PE) orientation and abnormal insert size (distance between read one and read two in a pair), (2) the presence of split and soft-clipped (SC) reads at the breakpoints of SVs, and (3) abnormal read depths (RD) in CNVs [2].

Exome sequencing (ES) is a cost-effective alternative to GS. ES protocols include exon-enrichment by probe hybridization, followed by high-depth sequencing of the enriched exonic regions (typically 100x coverage for clinical exome sequencing). Thus, ES allows reliable detection of small coding and near-splice-site variants within these targeted regions, which are causative for the largest fraction of genetic diseases based on current knowledge. However, depending on the enrichment strategy, ES almost completely lacks coverage of deep-intronic and intergenic variants, making the detection of variants, including SV

breakpoints, in these genomic regions essentially impossible. Since it is much more likely for breakpoints of SVs to occur in the >98% non-exonic regions of the human genome, SVs are exceedingly hard to detect with ES data. Hence, usage of PE and SC signals is restricted to breakpoint-detection in regions covered by ES reads [3] and thus, SV detection in ES data is mainly limited to detection of CNVs (deletions and duplications) using the normalized RD signature. It is worth noting that the RD signal in ES suffers from many biases such as GC-content correlated coverage bias and does not allow a robust detection of short coding CNVs affecting only one to several exons, while short CNVs displaying PE or SC signals in addition to changes in RD are more reliably detectable. Therefore, a frequently used approach for SV detection in patients with negative ES results is to perform additional NGS analyses, such as short or long read GS, which increases both cost and time to diagnosis.

Despite the aforementioned issues with detection of SVs in ES data, for approximately 2-4% of SVs (depending on the size of the exome kit's target region) we would expect the breakpoints to occur close enough to a targeted region, to be detected. In this paper, we evaluate the increase in diagnostic yield by PE- and SC-based SV calling in a large dataset of more than six thousand individuals with rare diseases (Solve-RD data freezes 1 and 2 [4]) who remained undiagnosed after standard ES analysis. Despite the large data heterogeneity, we show a valuable, albeit small

<footnote n="1">¹Institute of Medical Genetics and Applied Genomics, University of Tübingen, Tübingen, Germany. ²Centro Nacional de Análisis Genómico (CNAG), C/Baldiri Reixac 4, 08028 Barcelona, Spain. ³Department of Precision Medicine, University of Campania 'Luigi Vanvitelli', Naples, Italy. ⁴Telethon Institute of Genetics and Medicine, Pozzuoli, Italy.</footnote>

<footnote n="5">⁵Department of Neurosciences, Rehabilitation, Ophthalmology, Genetics, Maternal and Child Health, Università Degli Studi di Genova, Genoa, Italy. ⁶Medical Genetics Unit, IRCCS Istituto Giannina Gaslini, Genoa, Italy. ⁷Centre for Rare Diseases, University of Tübingen, Tübingen, Germany. ⁸Division of Evolution & Genomic Sciences, School of Biological Sciences, Faculty of Biology, Medicine and Health, University of Manchester, Manchester, UK. ⁹Manchester Centre for Genomic Medicine, St. Mary's Hospital, Manchester University NHS Foundation Trust, Health Innovation Manchester, Manchester, UK. ¹⁰Institute of Neurogenetics, University of Lübeck and University Hospital Schleswig-Holstein, Lübeck, Germany. *A list of authors and their affiliations appears at the end of the paper. email: German.Demidov@med.uni-tuebingen.de</footnote>
increase in the diagnostic yield following SV calling, quality
filtering, functional annotation, and clinical interpretation. We
demonstrate the benefit of SV calling as part of ES reanalysis for
patients with no previously identified molecular cause, under-
lining the value of applying SV calling algorithms in ES (re-)analysis
projects.

MATERIALS AND METHODS

The recruitment of individuals investigated by the Solve-RD consortium has been described in detail elsewhere [4]. The part of the Solve-RD cohort used in this study includes 9351 exome samples from 9314 individuals (Solve-RD data freezes 1 and 2), including 6224 affected individuals from 5825 families and 3090 unaffected relatives, usually parents. Samples were collected from multiple centers across Europe within four European Reference Networks (ERNs). Our cohort includes (1) 1,892 (30.4%) affected individuals from 1821 families from ERN ITHACA ("Intellectual disability, TeleHealth And Congenital Anomalies", including 26 affected individuals (0.4%) from 21 families from Undiagnosed Rare Diseases Program (Spain)), (2) 2,343 (37.6%) patients from 2200 families from ERN RND (Rare Neurological Diseases), (3) 1632 (26.2%) patients from 1499 families from EURO-NMD (ERN for NeuroMuscular Diseases) and (4) 357 (5.7%) patients from 340 families from ERN-GENTURIS (GENetic TUMour RISK Syndromes) (see Laurie et al., 2024, under review, for more details). Patients were submitted by 44 clinical research groups from these 4 ERNs, who used 28 different exome enrichment kits (various vendors and kit versions). Human phenotype ontology (HPO) terms, family relationship information, and raw ES data for all patients were submitted to the RD-Connect GPAP [5] by the corresponding clinicians and submitters.

Reads were aligned to the GRCh37 hs37d5 human reference genome using BWA-MEM 0.7.8. Manta SV caller [6] was used to detect SVs using default parameters and exome flag on. The population allele frequency (AF) of a detected breakpoint was estimated across the complete dataset using the number of breakpoints detected within ±20 base pairs to the focal breakpoint. This vicinity was defined empirically based on the confidence intervals provided by Manta for each detected breakpoint. For approximately 60% of the SVs the confidence intervals are narrower than 40 base pairs, with a median of 17, thus, larger windows would lead to overestimation, while smaller windows would lead to an underestimation of the allele frequency of SVs. Filtering based on breakpoint frequency in the cohort allowed removal of frequently observed events (polymorphisms) and artifacts. We kept only SVs with a breakpoint frequency of less than or equal to 20 out of 9351 exome datasets from affected and unaffected individuals. This AF threshold of ~0.21% was selected empirically, considering the largest family sizes identified via the relatedness analysis between all samples (largest family has 18 members) and previously reported highest frequencies of variants involved in recessive disorders. Moreover, only variants flagged as "PASSED" by Manta were kept for further interpretation.

Samples submitted by ITHACA, RND, NMD, and GENTURIS were filtered using lists of known disease genes provided by clinical experts from each ERN (3081, 1820, 611, 230 genes, respectively, see Laurie et al. 2024). Only SVs affecting at least one exon or located within 5 bp of an exon of these genes, were reported to the corresponding submitters for clinical interpretation. Intersection of SVs with short HIGH impact variants [7] and known pathogenic missense variants [8] previously identified in probands, for the identification of compound heterozygous events yielded no additional candidate SVs in disease gene lists. Finally, samples with candidate SVs on more than 5 chromosomes were filtered out, since high numbers of candidate SVs likely suggest low quality of the DNA or sequencing data rather than the presence of several possibly disease-causing rare SVs on different chromosomes. Annotation of SVs reported to clinicians included the following features: chromosome and position of the breakpoints, allele frequency within the Solve-RD exome cohort, affected status of the sequenced individual, HPO terms of the index case, ORDO code, genes potentially affected by the SV, other samples showing the same SV.

Further evaluation of technical quality and potential causality was
performed by clinical researchers from the 42 submitting groups. All SV
calls affecting genes reported as autosomal dominant (AD) and X-linked
(XL, dominant [XLD] or recessive [XLR]) in OMIM [9] with AF less than 4 per
cohort (thus, likely affecting members of only one family or a few
independent patients) were evaluated. Investigation of potential biallelic
variants in recessive disease genes, i.e., a combination of an SV with a small

heterozygous variant on the other allele, was performed by submitters if
sufficient expert time was available, which yielded one solved case.
Technical validation of calls was undertaken at corresponding facilities.
Clear-cut variants were considered as "validation not required". We define
"clear-cut" as events matching all the following criteria: (1) SVs supported
by multiple lines of evidence such as normalized coverage depth (CNV
calling), split reads, paired-end distance and orientation, and B-allele
frequency, (2) phenotypic description (HPO terms) matching to an affected
gene, (3) segregation confirmed in all available family members, and (4)
visual inspection of the area did not indicate an abundance of "random"
split reads or unusual paired end distances, which may indicate some sort
of sequencing quality failure. Hence, such SVs required multiple sources of
sequencing signal, such as PE and RD in samples with no obvious quality
issues, a match of one affected gene with the disease phenotype and the
inheritance patterns, as well as no alternative variants explaining the
phenotype.

Visual inspection and quality assessment were undertaken using the IGV
genome browser. We generated screenshots of the left and right
breakpoints and the complete SV as described in [10]. Screenshots
accompanying the filtered variants were returned to the data submitters for
inspection.

In addition to Manta, we tested InDelible [11], an SV caller developed specifically for ES data. The suggested pipeline was run according to the authors' recommendations and further filtered according to the recommended "strict" routine (https://github.com/HurlesGroupSanger/indelible), with the exception of trio-specific filters since most of our samples were singletons. In addition to the author recommended filters, all samples with more than 80,000 detected SVs were filtered out resulting in removal of 7.5% of cases. Results from InDelible were not submitted for expert evaluation, as no additional candidates were obtained, while several good candidates identified by Manta were missed.

Short variant analysis and phenotypic data collection was achieved using the RD-Connect Genome-Phenome Analysis Platform (GPAP) as described in [12]. Parallel RD-based CNV analysis using ClinCNV, ExomeDepth and Conifer was performed as described in [10].

RESULTS

We detected, quality-filtered, and annotated SVs in 9351 ES
datasets collected within the Solve-RD project. SV callsets from
unaffected relatives were used for segregation analysis,
population-AF based filtering and identification of systematic
errors, which helped to dramatically reduce the number of
potentially causal SV candidates for in-depth inspection. Following
automated quality- and annotation-based filtering of the raw SV
callsets generated by Manta, 1404 SV in 868 samples remained in
ITHACA, 798 SV in 487 samples in RND, 1519 SV in 606 samples in
NMD and 15 SV in 15 samples in GENTURIS callsets. The
distribution of SVs per sample was not uniform, and some
samples contained a high number of candidates for clinical
interpretation.

Upon expert-inspection by experts from the corresponding European Research Networks participating in Solve-RD, 23 distinct SVs were considered to be causal in 32 out of 6,224 (0.51%) affected individuals, pertaining to 23 unrelated families (Table 1). Eight of these 23 SVs (0.13% of the affected individuals) were not reported by read depth (RD) based CNV detection using ClinCNV [13], ExomeDepth [14] and Conifer [15] performed in parallel [10] (Table 2). To evaluate the added diagnostic value in comparison with conventional read-depth CNV analysis for ES, we present the identified SVs in three categories.

First, a detected SV can be classified as a simple deletion or
duplication (CNV). 15 distinct SVs in 19 patients were simple CNVs,
however, 5 of them were not detected by the CNV detection
approaches due to their small size (typically affecting a single exon
or only part of an exon affected by deletions with overall length
ranging from 66 to 3077 base pairs). The longest CNV, not
detected by RD methods, was a 3077 base-pair long deletion,
which only affects 36 base pairs of exon 23 of *SHANK3*
(NM_001372044.1), an exon with a total length of 2.2 kb. One
simple deletion affecting the recessive gene *B4GALNT1* was not
<table><caption>**Table 1.** The total number of structural variant calls, the number of evaluated calls and the diagnostic value increase per ERN. Total number of affected individuals denotes all the affected family members including index cases.</caption><thead><tr><th></th><th>ERN RND</th><th>ERN ITHACA</th><th>ERN NMD</th><th>ERN GENTURIS</th></tr></thead><tbody><tr><td>Number of affected individuals</td><td>2.343</td><td>1.892</td><td>1.632</td><td>357</td></tr><tr><td>Number of index patients</td><td>2.2</td><td>1.821</td><td>1.499</td><td>340</td></tr><tr><td>Known disease genes in gene list</td><td>1.82</td><td>3.081</td><td>611</td><td>230</td></tr><tr><td>Number of candidate variants, after filtering</td><td>798</td><td>1.404</td><td>1.519</td><td>15</td></tr><tr><td>Number of samples with SVs, after filtering</td><td>487</td><td>868</td><td>606</td><td>15</td></tr><tr><td>Number of solved index patients/all affected patients</td><td>7 (0.32%) / 11</td><td>9 (0.49%) / 9</td><td>6 (0.4%) / 9</td><td>1 (0.29%) / 3</td></tr><tr><td>Percentage of causal SVs among investigated SVs</td><td>1.37%</td><td>0.64%</td><td>0.59%</td><td>20%</td></tr></tbody></table>

<table><caption>**Table 2.** Variants detected via paired-end or soft-clipped signal based SV analysis (Manta) in exomes, considered to be causative for the corresponding rare diseases.</caption><thead><tr><th colspan="9">SIMPLE CNVS</th></tr><tr><th>ERN</th><th>Chr</th><th>Start</th><th>End</th><th>Manta Type</th><th>IGV evaluation</th><th>Gene(s)</th><th>CNV analysis</th><th>Length</th></tr></thead><tbody><tr><td>GENTURIS</td><td>16</td><td>68845794</td><td>68999240</td><td>BND_0/1</td><td>DEL</td><td>CDH1</td><td>YES</td><td>153446</td></tr><tr><td>NMD</td><td>10</td><td>69917991</td><td>69950535</td><td>BND_0/1</td><td>TANDEM_DUP (within gene)</td><td>MYPN</td><td>YES</td><td>32544</td></tr><tr><td>ITHACA</td><td>6</td><td>31630177</td><td>31657902</td><td>BND_0/1</td><td>DEL</td><td>CSNK2B</td><td>YES</td><td>27725</td></tr><tr><td>NMD</td><td>2</td><td>179448847</td><td>179459005</td><td>BND_0/1</td><td>DEL</td><td>TTN</td><td>YES</td><td>10158</td></tr><tr><td>RND</td><td>12</td><td>58014702</td><td>58024265</td><td>BND_0/1</td><td>DEL</td><td>B4GALNT1</td><td>YES</td><td>9563</td></tr><tr><td>RND</td><td>16</td><td>89610020</td><td>89619318</td><td>BND_1/1</td><td>DEL</td><td>SPG7</td><td>YES</td><td>9298</td></tr><tr><td>ITHACA</td><td>2</td><td>162269403</td><td>162274096</td><td>BND_0/1</td><td>DEL</td><td>TBR1</td><td>YES</td><td>4693</td></tr><tr><td>ITHACA</td><td>22</td><td>51160830</td><td>51163907</td><td>DEL_0/1</td><td>DEL</td><td>SHANK3</td><td>NO</td><td>3077</td></tr><tr><td>RND</td><td>11</td><td>66472730</td><td>66475186</td><td>BND_0/1</td><td>DEL</td><td>SPTBN2</td><td>YES</td><td>2456</td></tr><tr><td>RND</td><td>3</td><td>38938611</td><td>38939384</td><td>BND_0/1</td><td>DEL</td><td>SCN11A</td><td>NO</td><td>773</td></tr><tr><td>NMD</td><td>17</td><td>48247360</td><td>48247704</td><td>DEL_1/1</td><td>DEL</td><td>SGCA</td><td>YES</td><td>344</td></tr><tr><td>ITHACA</td><td>6</td><td>1.58E+08</td><td>1.58E+08</td><td>BND_0/1</td><td>DEL</td><td>ARID1B</td><td>NO</td><td>245</td></tr><tr><td>NMD</td><td>6</td><td>152485331</td><td>152485416</td><td>DUP_0/1</td><td>TANDEM_DUP (within gene)</td><td>SYNE1</td><td>NO</td><td>85</td></tr><tr><td>ITHACA</td><td>X</td><td>73749063</td><td>73749133</td><td>DUP_0/1</td><td>TANDEM_DUP (within gene)</td><td>SLC16A2</td><td>NO</td><td>70</td></tr><tr><td>ITHACA</td><td>14</td><td>36987112</td><td>36987178</td><td>DEL_0/1</td><td>DEL</td><td>NKX2-1</td><td>NO</td><td>66</td></tr><tr><th colspan="9">PARTS OF COMPLEX EVENTS</th></tr><tr><td>RND</td><td>14</td><td>50878079</td><td>51132277</td><td>BND_0/1</td><td>INVERSION</td><td>ATL1</td><td>YES</td><td>254198</td></tr><tr><td>RND</td><td>X</td><td>153630035</td><td>153668349</td><td>BND_0/1</td><td>COMPLEX_DUP</td><td>ATP6AP1;GDI1;RPL10;TAZ</td><td>YES</td><td>38314</td></tr><tr><td>ITHACA</td><td>16</td><td>2214362</td><td>2229919</td><td>BND_0/1</td><td>DUP (partial)</td><td>TRAF7</td><td>YES</td><td>15557</td></tr><tr><td>RND</td><td>3</td><td>11070530</td><td>11075542</td><td>BND_0/1</td><td>TANDEM_DUP (within gene) plus DEL</td><td>SLC6A1</td><td>YES</td><td>5012</td></tr><tr><td>NMD</td><td>2</td><td>179553482</td><td>179557174</td><td>DUP_0/1</td><td>DEL+DUP</td><td>TTN</td><td>YES</td><td>3692</td></tr><tr><td>ITHACA</td><td>X</td><td>152958716</td><td>152959469</td><td>Mix</td><td>Retroduplication</td><td>SLC6A8</td><td>YES</td><td>753</td></tr><tr><th colspan="9">COPY-NUMBER NEUTRAL</th></tr><tr><td>ITHACA</td><td>9</td><td>130887682</td><td>140727114</td><td>BND_0/1</td><td>INVERSION</td><td>EHMT1</td><td>NO</td><td>9839432</td></tr><tr><td>NMD</td><td>X</td><td>30960717</td><td>31140070</td><td>BND_1/1</td><td>INVERSION</td><td>DMD</td><td>NO</td><td>179353</td></tr></tbody></table>

reported in RD-based CNV analysis due to an excessive number of CNV candidates in the sample, but was prioritized in our SV analysis, which produced a much smaller number of calls for expert evaluation. Further inspection of this 9563 bp deletion, removing six exons, revealed that it was detected by RD-based CNV callers. Moreover, a missense variant, known to be pathogenic, was detected in the other allele, thus forming a compound heterozygote explaining the phenotype [Laurie at al 2024, under review]. It is worth mentioning that even for SVs detected in parallel by RD-based CNV callers the accurate information on breakpoint coordinates can facilitate the functional interpretation of the variant, as in e.g. [16].

The second category comprises variants for which a detected SV can be a part of a complex SV (a combination of deletions, duplications, and inversions), part of which is detectable via standard RD-based CNV analysis. It is rarely possible to detect all breakpoints in this situation using ES data. However, detected SVs in combination with CNVs detected by RD-based methods can indicate and potentially better resolve the complex nature of an

event. Six variants which were considered causal in 10 patients were complex events, such as a deletion followed by a duplication, or an inversion followed by a duplication.

Finally, a detected SV can be copy-number neutral. We were able to find causal inversions, but no translocations. Two pathogenic inversions were found in two patients. The first one presented an inversion of almost 10 Mb in length affecting *EHMT1* (MIM * 607001) in intron 25 of 26 (g:9:130,887,682-140,727,115, hg19) (ERN-ITHACA; Fig. 1C) This subject meets the criteria for Kleefstra syndrome type 1 (KS1) [17], a well-described syndromic neurodevelopmental condition characterized by psychomotor delay, cognitive impairment, behavioral disorders, facial dysmorphism, abnormal skull shape, abnormalities of hands, and congenital heart defects [17, 18]. Indeed, this patient presented with severe intellectual disability including absence of speech, hand stereotypies, aggressivity, and hypotonia. Additional clinical manifestations consistent with a possible diagnosis of KS1 included dysmorphic features (e.g., sparse eyebrows, short nose, protruding tongue, absence of lateral incisors) and hand
<figure><img src="image_5.png" /><figcaption>**Fig. 1** Visualisation of read alignments for various types of pathogenic SVs in the tool IGV. A Simple duplication of 70 base pairs missed by RD-based CNV analysis and short variant calling. B A complex SV: paired-end distance indicates the presence of a tandem duplication. RD-based CNV calling (top 2 tracks) indicates the presence of a deletion next to the detected duplication, together forming a complex SV. C A 9.8Mb-long genomic inversion affecting the penultimate intron in *EHMT1* was detected via structural variant calling. D 179Kb long inversion involving the last exon of the DMD gene. Screenshots produced using Integrative Genomics Viewer [21]. Red arrows indicate the breakpoints identified via SV analysis.</figcaption></figure>

abnormalities (syndactyly and drumstick fingers). Considering the suggestive phenotype, genetic testing for KS1 had previously been performed through array comparative genomic hybridization (aCGH) and Sanger sequencing of the *EHMT1* gene, according to the suggestions from the diagnostic guidelines [17]. However, these tests did not detect any possibly pathogenic variant in the candidate gene, leading to several years of a delay in diagnosing the patient despite the strongly suggestive clinical and phenotypic observations. Following the identification of the 10 Mb inversion identified here, karyotyping and FISH were performed to furtherly validate the molecular diagnosis. Karyotyping at 400 band resolution did not identify the inversion that was correctly observed by FISH (probe utilized RP11-31M4 (9q34.3) Empire Genomics) (data not shown). The second patient had an inversion of approximately 180 kb in length affecting *DMD* (ERN-NMD; Fig. 1D). Both inversions occurred in patients with strikingly fitting phenotypes which had been extensively examined for variants in the corresponding genes as candidate genes for several years before being enrolled into the Solve-RD project. The genes are directly affected by these inversions, leading to loss-of-function phenotypes.

Thus, the added diagnostic value gained solely from paired-end mapping-based SV detection in the reanalysis of unsolved individuals, was 8 out of 5825 index patients (0.14%) for pipelines that include a comprehensive read-depth based CNV detection

method. This value increases to 0.4% (23/5,825) for pipelines detecting only short variants, hence lacking CNV calls.

In addition to Manta, we initially tested InDelible [11], an SV caller developed specifically for ES data. Despite the lower allele frequency filtering threshold, recommended by the authors (0.04% in comparison with 0.2% we used in our analysis), 684 variants passed the filtering (Manta: 1976, Table 1). A preliminary examination of the calls revealed no additional candidates not reported by Manta, but InDelible missed all the clear rare candidates detectable with paired-end information. Thus, InDelible results were not submitted for expert evaluation and were excluded from further use in this study.

## DISCUSSION

Here we demonstrate the utility of PE- and SC-based SV calling in ES data for undiagnosed individuals with a rare disease using a large cohort of almost 10,000 individuals from the Solve-RD project. Despite the modest overall increase in the diagnostic yield, each successfully diagnosed patient represents a family whose diagnostic odyssey is coming to an end. Especially for patients initially sequenced many years ago, alternative ways of investigating structural variants, such as GS or long-read sequencing, may be unavailable due to financial or logistical reasons e.g. lack of DNA samples. Thus, re-analysis approaches
using existing ES data to solve unsolved patients can benefit from PE-based SV calling. In a global effort to solve previously unsolved patients, SV calling in ES data could help hundreds of people to find their diagnosis after many years of being undiagnosed.

Even though the majority of our findings (15 out of 23 unique causal SVs) were also discovered in parallel via read-depth based CNV analysis, SV identification using PE and SC signals should not be considered redundant in such cases. In some patients (6 out of 23 variants) the detected SVs occur next to CNVs detected by RD callers, but knowledge of the SV's presence is crucial for the interpretation of the rearrangement, uncovering its underlying complexity, and may indicate the necessity for additional molecular analyses such as long read DNA sequencing or RNA sequencing. Furthermore, SV detection in ES, even when it identifies the same variant also found by RD methods, increases the confidence that the variant is a true positive. PE-based SV calling facilitates the validation of CNVs since it provides information about the exact breakpoint coordinates and enables locus-specific PCR amplification and Sanger sequencing. Moreover, we were able to identify eight SVs that were not identified by RD methods [10], five of which were simple deletions and duplications, two were inversions, and one was a larger deletion forming a compound heterozygote with a SNV in an autosomal recessive disease gene (which was not reported based on RD-calling due to an excessive number of CNV calls in recessive genes).

Notably, both inversions were identified in patients with strikingly matching phenotypes, who had undergo multiple rounds of genetic testing of the respective genes before the start of Solve-RDIn this regard, the reanalysis of the genetic data focused on the investigation of structural variants (SV) through the Solve-RD platform was fundamental for the identification of a balanced rearrangement causing the disruption of *EHMT1*, thus providing the pathophysiological link to the genetic condition displayed by the proband of this family. This achievement ended the diagnostic odyssey of this family, finally offering a genetic diagnosis and the chance of targeted clinical management.

In comparison with widely used RD based CNV callers, we cannot estimate the recall rate for detecting variants with PE-based SV calling from ES data, since breakpoints need to affect a targeted region of the exome. For this reason, the chance to successfully identify an SV is independent of its length, as its breakpoints must be covered with a significant number of short reads exhibiting abnormal orientation or insert size. Nonetheless, among our more than 6000 index patients, five CNVs were detected that were missed by CNV analysis since they were too short. Furthermore, RD-based CNV-calling cannot detect copy-number-neutral events, such as two inversions in our cohort. 16 out of 105 unique CNVs identified as being pathogenic or likely-pathogenic [10] were also detected in parallel by SV analysis (24 out of 115 solved, candidate or partially explaining the phenotype patients), allowing to better define their breakpoints. Interestingly, this fraction is higher than expected when considering the genomic region covered by ES (< 2%), likely because we focus our analysis on SVs and CNVs overlapping known disease genes. Moreover, the chance of detected CNVs to be true positives increases, if they are supported by paired-end and split-read signatures, allowing genetic analysts to evaluate the clinical relevance of the variant without doubts regarding the technical quality. Finally, CNVs detected in ES data do not provide accurate breakpoint information as SV calling can, which further simplifies the clinical interpretation of variants. Hence, we conclude that SV calling in ES often provides valuable additional information even for variants already detected by CNV analysis.

A previously published study by [11] in which SVs were analyzed in a cohort of 13,438 probands from the Deciphering Developmental Disorders (DDD) study came to comparable conclusions regarding the fraction of detectable causal SVs. Gardner et al. used a

combination of the tool XHMM [19] for RD-based CNV calling and InDelible for split read (SR)-based SV calling, while we combined three RD-based tools with the SV caller Manta. Gardner et al. report 30 unique pathogenic SVs identified by InDelible (0.22%) in addition to 128 CNVs detected by XHMM (0.95%). Interestingly, the fraction of causal CNVs found via the RD-method XHMM is substantially lower compared to the combination of three RD-Methods in Solve-RD (0.95% vs. 1.6%) [10], but on the other hand the additional diagnostic value of SVs is slightly higher (0.22% vs. 0.14% in our study). Of note, applying InDelible to our cohort (data not shown) did not reveal additional SVs except for 3 mobile element insertions (MEI), which we had already detected using specialized tools [20], but it did miss several causal SVs identified by Manta. We conclude that despite using different patient cohorts and different approaches for the detection of SVs and CNVs, the DDD and Solve-RD initiatives produced comparable results regarding the fraction of causal SVs detected.

The categorization of the identified variants into three groups also allows us to specify the technical and analytical challenges encountered while evaluating them. For the first group (simple deletions and duplications) the main challenge was to evaluate the technical quality of the calls which were missed by RD methods. In these cases, we had to rely only on paired-end and split-read information, since the coverage did not change visibly beyond the standard level of noise. The evaluation of the SVs that we defined as "complex" necessarily involved joint visual analysis of RD and PE signals. Since we could often only detect one pair of breakpoints in "complex" SVs, only the presence of abnormal coverage next to the identified SV indicated the complex nature of the event. In the case of inversions, in most cases Manta did not report the type of variant as INV (inversion) but BND (breakpoint) in ES data, and hence all breakpoints had to be visually explored in order to identify inversions. This weakness can only be overcome by genome sequencing, which will typically result in clear identification of an inversion.

As discussed, the PE-based SV calling approach using ES has some weaknesses. We did not perform an evaluation of SV calling sensitivity in ES since the sensitivity is low, as expected by the nature of targeted sequencing. Usage of Manta or analogous tools for SV detection cannot replace RD-based CNV detection. SV calling should only be considered as an addition, which may result in a slight increase indiagnostic yield but does not guarantee robust detection even of long rearrangements if their breakpoints are intergenic or deep-intronic. Furthermore, improvements and automation of the clinical interpretation of SVs occurring in recessive genes are necessary to reduce the number of reported calls for expert evaluation, as all tools report many false positives. One way to achieve this is through automated phenotypic matching procedures, which we plan to implement during the analysis of the final Solve-RD data freeze and will report the results to the medical genetics community.

In summary, PE and SC based SV calling is a valuable addition to RD-based CNV calling, providing a diagnosis to a small but important fraction of rare disease patients, who would otherwise remain undiagnosed.

## DATA AVAILABILITY

All raw and processed data files are available at the EGA (Datasets EGAD00001009767, EGAD00001009768, EGAD00001009769, and EGAD00001009770, under Solve-RD study EGAS00001003851). The family (FAM) and participant (P) identifiers used in this manuscript are pseudonymized and known only to the researchers involved In Solve-RD.

## REFERENCES

1. Deciphering Developmental Disorders Study. Large-scale discovery of novel genetic causes of developmental disorders. Nature. 2015;519:223-8.
2. Tattini L, D'Aurizio R, Magi A. Detection of genomic structural variants from next-generation sequencing data. Front Bioeng Biotechnol. 2015;3:92.

3. Sadedin SP, Ellis JA, Masters SL, Oshlack A. Zimmer: a system for improving accuracy and consistency of CNV calling from exome data. Gigascience. 2018;7:giy112.

4. Zurek B, Ellwanger K, Vissers LELM, Schüle R, Synofzik M, Töpf A, et al. Solve-RD: systematic pan-European data sharing and collaborative analysis to solve rare diseases. Eur J Hum Genet. 2021;29:1325–31.

5. Laurie S, Piscia D, Matalonga L, Corvó A, Fernández-Callejo M, Garcia-Linares C, et al. The RD-connect genome-phenome analysis platform: accelerating diagnosis, research, and gene discovery for rare diseases. Hum Mutat. 2022;43:717–33.

6. Chen X, Schulz-Trie glaff O, Shaw R, Barnes B, Schlesinger F, Källberg M, et al. Manta: rapid detection of structural variants and indels for germline and cancer sequencing applications. Bioinformatics. 2016;32:1220–2.

7. McLaren W, Gil L, Hunt SE, Riat HS, Ritchie GRS, Thormann A, et al. The ensemble effect predictor. Genome Biol. 2016;17:122.

8. Landrum MJ, Chitipiralla S, Brown GR, Chen C, Gu B, Hart J, et al. ClinVar: improvements to accessing data. Nucleic Acids Res. 2020;48:D835–44.

9. Amberger JS, Bocchini CA, Schiettecatte F, Scott AF, Hamosh A. OMIM.org: Online Mendelian Inheritance in Man (OMIM®), an online catalog of human genes and genetic disorders. Nucleic Acids Res. 2015;43:D789–98.

10. Demidov G, Yaldiz B, Garcia-Pelaez J, de Boer E, Van de Vondel L, Paramonov I, et al. Comprehensive reanalysis for CNVs in ES data from unsolved rare disease cases results in new diagnoses. medRxiv [Internet]. 2023; Available from: https://doi.org/10.1101/2023.10.22.32396993

11. Gardner EJ, Sifrim A, Lindsay SJ, Prigmore E, Rajan D, Danecek P, et al. Detecting cryptic clinically relevant structural variation in exome-sequencing data increases diagnostic yield for developmental disorders. Am J Hum Genet. 2021;108:2186–94.

12. Matalonga L, Hernández-Ferrer C, Piscia D, Solve-RD SNV-indel working group, Schüle R, Synofzik M, et al. Solving patients with rare diseases through programmatic reanalysis of genome-phenome data. Eur J Hum Genet. 2021;29:1337–47.

13. Demidov G, Sturm M, Ossowski S ClinCNV: multi-sample germline CNV detection in NGS data. bioRxiv [Internet]. 2022; Available from: https://doi.org/10.1101/2022.06.10.495642

14. Plagnol V, Curtis J, Epstein M, Mok KY, Stebbings E, Grigoriadou S, et al. A robust model for read count data in exome sequencing experiments and implications for copy number variant calling. Bioinformatics. 2012;28:2747–54.

15. Krumm N, Sudmant PH, Ko A, O’Roak BJ, Malig M, Coe BP, et al. Copy number variation detection and genotyping from exome sequence data. Genome Res. 2012;22:1525–32.

16. São José C, Garcia-Pelaez J, Ferreira M, Arrieta O, André A, Martins N, et al. Combined loss of CDH1 and downstream regulatory sequences drive early-onset diffuse gastric cancer and increase penetrance of hereditary diffuse gastric cancer. Gastric Cancer. 2023;26:653–66.

17. Kleefstra T, van Zelst-Stams WA, Nillesen WM, Cormier-Daire V, Houge G, Foulds N, et al. Further clinical and molecular delineation of the 9q subtelomeric deletion syndrome supports a major contribution of EHMT1 haploinsufficiency to the core phenotype. J Med Genet. 2009;46:598–606.

18. Willemse MH, Vulto-van Silfhout AT, Nillesen WM, Wissink-Lindhout WM, van Bokhoven H, Philip N, et al. Update on Kleefstra Syndrome. Mol Syndromol. 2012;2:202–12.

19. Fromer M, Purcell SM. Using XHMM software to detect copy number variation in whole-exome sequencing data. Curr Protoc Hum Genet. 2014;81:7.23–7.23.21.

20. Wijngaard R, Demidov G, O’Gorman L, Corominas-Galbany J, Yaldiz B, Steyaert W, et al. Mobile element insertions in rare diseases: a comparative benchmark and reanalysis of 60,000 exome samples. Eur J Hum Genet. 2024;32:200–208.

21. Robinson JT, Thorvaldsdóttir H, Winckler W, Guttman M, Lander ES, Getz G, et al. Integrative genomics viewer. Nat Biotechnol. 2011;29:24–6.

## ACKNOWLEDGEMENTS

We want to gratefully acknowledge the interpretation, evaluation and validation of the candidate structural variants, performed by the following experts: Miquel Raspall-Chaure (Pediatric Neurology Research Group, Vall d'Hebron Research Institute, Universitat Autònoma de Barcelona, Barcelona, Spain), Swati Naik (Birmingham Women's and Children's Hospital NHS Trust, Birmingham, UK), Lukáš Ryba, (Department of Biology and Medical Genetics, 2nd Faculty of Medicine, Charles University in Prague and Motol University Hospital, Prague, Czech Republic), Patrick Callier (FHU TRANSLAD, CHU Dijon, France and GAD team, INSERM UMR1231, Université de Bourgogne-Franche Comté, Dijon, France), Valeria Capra (Medical

Genetics Unit, IRCCS Istituto Giannina Gaslini, Genoa, Italy), Flavia Pivitera, Maria Antonietta Mencarelli, Ilaria Longo, Francesca Ariani, Kristina Zguro (Medical Genetics, Med Biotech Hub and Competence Center, Department of Medical Biotechnologies, University of Siena, Italy), David Gómez-Andrés (Pediatric Neurology Research Group, Vall d'Hebron Research Institute, Barcelona, Spain). We acknowledge bioinformatic analysis undertaken by Burcu Yaldiz (Department of Human Genetics, Radboud Institute for Molecular Life Sciences, Radboud University Medical Center), Leon Schuetz (Institute of Medical Genetics and Applied Genomics, University of Tübingen, Tübingen, Germany).

## AUTHOR CONTRIBUTIONS

Solve-RD consortium: contribution of data, expertise and infrastructure. GD, SL performed the bioinformatic analysis. AT, JP, MS, MM, VN, SB, KL evaluated the variants, performed their validation and clinical reexamination if necessary. GD, HG, SL, KL, SO developed the study design. HG, SL, KL and SO supervised the project. GD wrote the manuscript, which was reviewed and approved by all authors. All other Solve-RD authors: data acquisition, the interpretation, evaluation and validation of the candidate structural variants.

## FUNDING

The Solve-RD project has received funding from the European Union's Horizon 2020 research and innovation programme under grant agreement No 779257. This work was supported in part by Telethon Undiagnosed Diseases Program (TUDP, GSP15001). Open Access funding enabled and organized by Projekt DEAL.

## COMPETING INTERESTS

The authors declare no competing interests.

## ETHICS STATEMENT

The Ethics committee of the Eberhard Karl University of Tubingen gave ethical approval for this work. The responsibility of checking the data is suitable for submission to the RD-Connect GPAP and Solve-RD, including informed consent, lies within the data submitter as required by their Code of Conduct and Data Sharing Policy, respectively. In some cases, individuals had to be re-consented prior to data submission. This study adheres to the principles set out in the Declaration of Helsinki.

## ADDITIONAL INFORMATION

**Supplementary information** The online version contains supplementary material available at https://doi.org/10.1038/s41431-024-01637-4.

**Correspondence** and requests for materials should be addressed to German Demidov.

**Reprints and permission information** is available at http://www.nature.com/reprints

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_6.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2024, corrected publication 2025
SOLVE-RD CONSORTIUM

German Demidov<sup>1,3,4</sup>, Steven Laurie<sup>2</sup>, Annalaura Torella<sup>3,4</sup>, Giulio Piluso<sup>1,3</sup>, Marcello Scala<sup>5,6</sup>, Manuela Morleo<sup>3,4</sup>, Vincenzo Nigro<sup>3,4</sup>, Holm Graessner<sup>1,7</sup>, Siddharth Banka<sup>8,9</sup>, Alfons Macaya<sup>11,12</sup>, Belén Pérez-Dueñas<sup>11</sup>, Adam Jackson<sup>8,9</sup>, Giovanni Stevanin<sup>13,14,15,16,17</sup>, Jean-Madeleine de Sainte Agathe<sup>18</sup>, Markéta Havlovicová<sup>19</sup>, Rita Horvath<sup>20</sup>, Michele Pinelli<sup>4</sup>, Nienke J. H. van Os<sup>21</sup>, Bart P. C. van de Warrenburg<sup>21</sup>, Anne-Sophie Denommé-Pichon<sup>22</sup>, Marco Savarese<sup>23</sup>, Mridul Johari<sup>23</sup>, Bruno Dallapiccola<sup>24</sup>, Marco Tartaglia<sup>24</sup>, Martje G. Pauly<sup>10</sup>, Anna Katharina Sommer<sup>25</sup>, Tobias B. Haack<sup>1</sup>, Ana Töpf<sup>26</sup>, Lacombe Didier<sup>27</sup>, Chiara Fallerini<sup>28,29</sup>, Alessandra Renieri<sup>28,29,30</sup>, Patrick F. Chinnery<sup>20,31</sup>, Daniel Natera-de Benito<sup>32</sup>, Andres Nascimento<sup>32</sup>, Aurélien Trimouille<sup>33</sup>, Francina Munell<sup>11</sup>, Anna Marcé-Grau<sup>11</sup>, Ben Yaou Rabah<sup>34,35,36</sup>, Gisèle Bonne<sup>34</sup>, Liedewei Van de Vondel<sup>37,38,39</sup>, Katja Lohmann<sup>10</sup> and Stephan Ossowski<sup>1</sup>

<footnote n="11"><sup>11</sup>Pediatric Neurology Research Group, Vall d'Hebron Research Institute, Universitat Autònoma de Barcelona, Barcelona, Spain. <sup>12</sup>Institute of Neuroscience, Universitat Autònoma de Barcelona, Barcelona, Spain. <sup>13</sup>Institut Nacional de la Santé et de la Recherche Médicale (INSERM) U1127, Paris, France. <sup>14</sup>Centre National de la Recherche Scientifique, Unité Mixte de Recherche (UMR) 7225, Paris, France. <sup>15</sup>Unité Mixte de Recherche en Santé 1127, Université Pierre et Marie Curie (Paris 06), Sorbonne Universités, Paris, France. <sup>16</sup>Institut du Cerveau - ICM, Paris, France. <sup>17</sup>Ecole Pratique des Hautes Etudes, Paris Sciences et Lettres Research University, Paris, France. <sup>18</sup>Department of Genetics, Assistance Publique-Hôpitaux de Paris—Sorbonne Université, Pitié-Salpétrière University Hospital, 83 Boulevard de l'Hôpital, Paris, France. <sup>19</sup>Department of Biology and Medical Genetics, Charles University Prague-2nd Faculty of Medicine and University Hospital Motol, Prague, Czech Republic. <sup>20</sup>Department of Clinical Neurosciences, University of Cambridge, Cambridge, UK. <sup>21</sup>Department of Neurology, Radboud University Medical Center, Nijmegen, The Netherlands. <sup>22</sup>Inserm - University of Burgundy-Franche Comté, UMR1231 GAD, Dijon, France.</footnote>

<footnote n="23"><sup>23</sup>Folkhälsan Research Centre and Medicum, University of Helsinki, Helsinki, Finland. <sup>24</sup>Molecular Genetics and Functional Genomics, Ospedale Pediatrico Bambino Gesù, IRCCS, Rome, Italy. <sup>25</sup>Institute of Human Genetics, Medical Faculty, University of Bonn, Bonn, Germany. <sup>26</sup>John Walton Muscular Dystrophy Research Centre, Translational and Clinical Research Institute, Newcastle University and Newcastle Hospitals NHS Foundation Trust, Newcastle upon Tyne, UK. <sup>27</sup>Univ. Bordeaux, MRGM INSERM U1211, CHU de Bordeaux, Service de Génétique Médicale, F-33000 Bordeaux, France. <sup>28</sup>Med Biotech Hub and Competence Center, Department of Medical Biotechnologies, University of Siena, Siena, Italy.</footnote>

<footnote n="29"><sup>29</sup>Medical Genetics, University of Siena, Siena, Italy. <sup>30</sup>Genetica Medica, Azienda Ospedaliero-Universitaria Senese, Siena, Italy. <sup>31</sup>Medical Research Council Mitochondrial Biology Unit, University of Cambridge, Cambridge, UK. <sup>32</sup>Neuromuscular Disorders Unit, Department of Pediatric Neurology. Hospital Sant Joan de Déu, Barcelona, Spain. <sup>33</sup>Laboratoire de Génétique Moléculaire, Service de Génétique Médicale, CHU Bordeaux—Hôpital Pellegrin, Place Amélie Raba Léon, 33076 Bordeaux Cede, France. <sup>34</sup>Sorbonne Université, Inserm, Institut de Myologie, Centre de Recherche en Myologie, F-75013 Paris, France. <sup>35</sup>AP-HP, Centre de Référence de Pathologie Neuromusculaire Nord, Est, Ile-de-France, Institut de Myologie, G.H. Pitié-Salpétrière, F-75013 Paris, France. <sup>36</sup>Institut de Myologie, Equipe Bases de données, G.H. Pitié-Salpétrière, F-75013 Paris, France. <sup>37</sup>Peripheral Neuropathy Research Group, University of Antwerp, Antwerp, Belgium. <sup>38</sup>Laboratory of Neuromuscular Pathology, Institute Born-Bunge, University of Antwerp, Antwerpen, Belgium.</footnote>

<footnote n="39"><sup>39</sup>Translational Neurosciences, Faculty of Medicine and Health Sciences, University of Antwerp, Antwerpen, Belgium. A full list of members and their affiliations appears in the Supplementary Information.</footnote>