1. Introduction
The DMD gene encodes the dystrophin protein, which binds the extracellular matrix to the cytoskeleton in the skeletal muscle fibres and is involved in sarcolemma integrity during muscle contraction [1]. Pathogenic variants in the DMD gene produce a weakening of the sarcolemma and muscle damage that results in muscle disease. The term dystrophinopathies includes Duchenne muscular dystrophy (DMD), Becker muscular dystrophy (BMD), and dilated cardiomyopathy (DCM) [2,3,4]. In 90% of these patients, the clinical presentation can be predicted by the reading frame rule [5,6]. DMD, the most severe phenotype with early childhood onset and rapid progression, is characterised by the absence of dystrophin expression due to out-of-frame variants. In contrast, BMD patients have a milder phenotype, later onset, and slower progression, showing reduced dystrophin staining due to a partially functional dystrophin protein caused by in-frame deletions or missense variants.
Dystrophinopathies follow an X-linked recessive inheritance pattern, in which males are affected by the disease and female carriers are usually asymptomatic. It is estimated that 2.5–19% of female DMD carriers may manifest muscle symptoms, ranging from mild muscle weakness and myalgia to a severe DMD phenotype [7,8,9,10,11,12]. Non-muscular manifestations such as cardiomyopathy are present in 7.3–16.7% of female carriers [12]. The most common explanation for the clinical manifestation in DMD female carriers is due to a skewed X chromosome inactivation (XCI) [13,14]; however, a lack of correlation between the XCI pattern, clinical manifestations, and dystrophin alterations in muscle biopsy has been described [15].
The most prevalent pathogenic variants in the DMD gene are exonic deletions (60%) and duplications (10%), followed by point mutations [16,17]. However, a small number of patients (2–7%) harbour deep-intronic variants or complex chromosomal rearrangements, such as inversions, insertions, or translocations, which cannot be identified by current diagnostic techniques, such as Multiplex Ligation-dependent Probe Amplification (MLPA) and panel, or even exome, sequencing. The DMD gene is the largest known gene in the human genome encoding several isoforms expressed in a variety of tissues. The full-length muscle isoform Dp417m consists of 79 exons corresponding to an open reading frame of 11.3 kb [1], but 99% of the DMD gene consists of intronic sequences. Large introns in the DMD gene harbour a substantial proportion of repetitive and transposable elements (TEs) that can result in recombination events or replication errors, thus predisposing the patient to chromosomal rearrangements. This can explain the high incidence of de novo variants in the DMD gene, which has been estimated at 33% [17,18].
Here, we present a female patient with a de novo translocation between Xp21 and 17q23.2 resulting in a DMD-like phenotype. This structural variant was characterised using an integrated multi-omics approach. This study highlights the utility of transcriptome and whole genome analysis in the detection and characterisation of structural variants, in particular, in patients where clinical diagnosis might have been misguided and standard cytogenetic testing overlooked.
2. Results
2.1. Clinical Presentation
The patient was born from healthy non-consanguineous parents. Early developmental milestones were achieved at the expected times, and she achieved independent walking at 16 months old. At the age of 7 years old, she presented with a history of frequent falls and difficulty running and standing up from the floor. On examination at that time, she had prominent calf hypertrophy and pelvic girdle weakness (Medical Research Council (MRC) grade 3+ to 4/5) and got up from the floor with a Gowers manoeuvre. Her serum creatine kinase (CK) levels were 14,050 u/L. The weakness was progressive, and at 10 years old she was commenced on steroid treatment with deflazacort, which was continued until she was approximately 25 years old. She was diagnosed with osteoporosis at 21 years old. At 18 years old, the patient was able to walk very short distances holding onto furniture indoors but used a powered wheelchair outside her home. This continued until she was 24 years old when she sustained fractures to her femur and ankle, leading to a loss of independent ambulation, which resulted in a significant weight gain.
At the most recent assessment, the patient was 29 years old and had a BMI of 39.1. She was a full-time powered wheelchair user and used a hoist for transfers. She had limited arm and leg movement with limited residual distal function preserved (MRC grade 3+ for wrists and 2 for ankle dorsiflexors). She was able to bring her right hand to her mouth with some compensation. She reported that her left arm was slightly weaker than the right one. Her spine was flexible with a slight single curvature on the back with a convexity to the right side. Her respiratory assessment showed a forced vital capacity of 0.86 L (25% predicted value), with no postural drop. Peak cough flow was also limited (208 L/min); however, an overnight pulse oximeter showed no sign of significant desaturations (Mean SPO2 98%, oxygen desaturation index for >4% of 0.7/h). Cardiac assessments have shown no evidence of impaired function, with ejection fraction >50%. Cognition was normal.
2.2. Histopathological Findings
Histological analysis showed dystrophic features including fibre size variability, internalised nuclei, and severe adipose and connective tissue replacement (Figure 1A). Immunohistochemical labelling for dystrophin was absent on many fibres, with only a few small groups showing labelling (Figure 1G,H). The same groups of fibres showed an overall reduction in labelling for α, β, γ, δ sarcoglycan, and β dystroglycan (Figure 1B–F). Upregulation of utrophin was also seen (Figure 1I). This pattern of labelling was not considered typical of a manifesting DMD carrier, where a mosaic of dystrophin-positive and -negative fibres is usually observed, as the primary deficiency of dystrophin was much more severe.
Western blotting was performed alongside immunohistochemistry (Figure 2). No C-terminal dystrophin labelling could be seen but a faint, clear band was detected for the dystrophin rod domain. This band was full-sized but severely reduced in abundance in comparison to the normal control. In addition, the abundance of the sarcoglycan and dystroglycan bands was also reduced but still visible.
2.3. Routine Genetic Studies
Initial standard single gene screening of the DMD (MLPA and sequencing), SGCA, SGCB, SCGG, SGCD, CAPN3, ANO5, and FRKP genes was negative, and an SNP array did not detect any alteration. Later, a neuromuscular disease in silico gene panel (via whole genome sequencing) of 201 genes was negative. Structural variants (SV) were not analysed. Trio whole exome sequencing (WES) was performed in a research setting, but no candidate variant was identified amongst 429 neuromuscular-related genes (as per Töpf A. et al. 2020 [19]). In addition, the X-chromosome inactivation pattern was considered to be within the normal range (73–27%).
2.4. Transcriptomics and Whole Genome Sequencing
RNA sequencing (RNAseq) from muscle biopsy and trio whole genome sequencing(WGS) were performed in parallel as part of the Solve-RD project (
Following this finding, copy number variants (CNVs) and SVs were evaluated in the trio-WGS data. No CNVs were identified within the DMD gene, but a putative translocation event between chromosome 17 and intron 16 of the DMD gene was detected in the proband. This translocation had the highest quality score and was found exclusively in this patient amongst the WGS Solve-RD cohort (n = 2303). Furthermore, this translocation was not present in the gnomAD SVs database (v4.1.0, last accession 14 June 2024).
Raw trio-WGS data were visualised on an Integrative Genomics Viewer (IGV) and revealed in the proband the presence of chimeric reads aligning to both chromosome 17 and chromosome X (Figure 4). Allele specific PCR primers were designed to validate the exact sequence of the breakpoints and to confirm the occurrence of this de novo balanced reciprocal translocation between the DMD gene (chrX) and the BCAS3 gene (chr17) (Figure 4). The breakpoints of the translocation t(X;17)(p21.1;q23.2) were chr17:58,861,454-chrX:32,578,628 (Figure 5A), and chrX:32,578,622-chr17:58,861,455 (Figure 5B), leading to two derivative chromosomes: one encoding exons 1 to 6 of the BCAS3 gene and exons 1 to 16 of the DMD gene and the second one encoding exons 7 to 25 of the BCAS3 gene and exons 17 to 79 of the DMD gene. A 3 bp microhomology at breakpoint 1 and an 11 bp insertion at breakpoint 2 with no homology in the nearby breakpoint sequence (300 bp) were detected. The breakpoint located in intron 16 of the DMD gene lies within the LINE element (L1PA10).
Despite this translocation disrupting two protein-coding genes, DMD on chromosome X and BCAS3 on chromosome 17, only DMD expression was affected, with normal BCAS3 expression and splicing in RNAseq analysis.
2.5. Karyotype
To further confirm the translocation, a karyotype was performed on patient’s peripheral blood (Figure 6). G-band analysis shows a female karyotype with 46 chromosomes, including an apparently balanced reciprocal translocation between the short arm of a chromosome X (break at band p21) and the long arm of a chromosome 17 (break at band q23).
3. Discussion
During the last decades, the increasing availability of next-generation sequencing (NGS), mainly exome and gene panel sequencing, is believed to have revolutionised the diagnostic journey of patients with rare diseases. Nonetheless, it is estimated that 50% of patients with muscle diseases remain undiagnosed after NGS testing [19,22,23]. This may be attributed to the difficulties in reclassifying variants of uncertain significance (VUS) or detecting pathogenic changes not covered by exome sequencing, such as non-coding or structural variants. In particular, RNA sequencing has been useful in the detection and interpretation of intronic cryptic variants and splicing defects [24,25,26]. In this study, we report a female patient with a DMD-like phenotype caused by a de novo reciprocal translocation t(X;17)(p21.1;q23.2) detected through a transcriptomics and genomics integrated approach.
The proband was a female patient with childhood-onset progressive proximal weakness and elevated creatine kinase and reduced staining of the dystrophin–glycoprotein complex (DGC) on muscle biopsy (Figure 1), in addition to a reduction in abundance of DGC proteins on Western blot (Figure 2). Given the clinical presentation and histopathological findings, MLPA and sequencing of the DMD gene were indubitably the first line of analysis, but all proved negative. By that time, novel NGS techniques became available in the diagnostic setting and trio-WES was consequently performed but failed to identify a candidate variant, leaving the patient still undiagnosed. Years later, RNAseq and trio-WGS were performed simultaneously in a research setting. RNAseq allowed the rapid identification of DMD as a gene with outlier expression (Figure 3; consequently, the SVs and CNVs analysis from the trio-WGS data were focused on the DMD gene revealing a de novo balanced translocation t(X;17)(p21.1;q23.2) (Figure 4 and Figure 5). This translocation completely disrupts the Dp472 muscle isoform of the DMD gene and the BCAS3 gene, which is an important cytoskeletal protein during embryogenesis, angiogenesis, and tumorigenesis, and has been recently identified to cause the autosomal, recessive Hengel–Maroofian–Schols syndrome [27].
Previous reports have proposed that the underlying mechanisms of chromosomal rearrangements in the DMD gene may involve a double-stranded break (DSB) followed by non-homologous end joining (NHEJ) [28], microhomology-mediated end joining (MMEJ) (also known as alternative end joining) [29], or a replication-based mechanism [30]. In the patient described here, the presence of a small microhomology (3 bp) in one of the translocation breakpoints (Figure 5) and the insertion of 11 bp not originating from the nearby sequence suggests that NHEJ, which acts on microhomologies of less than 4 bp [31], is likely to be the mechanism that caused the translocation. We hypothesise that the fact that the DMD breakpoint lies within a LINE element has occurred by chance, as transposable elements and repetitive sequences have been detected in chromosomal rearrangement breakpoints in a similar proportion to that found in the human genome [32,33].
The presence of clinical manifestations in female DMD carriers can be explained by (1) a DMD pathogenic variant in both alleles [34], (2) one DMD variant in the only X chromosome due to uniparental disomy or Turner syndrome (45,X) [35,36], or (3) a skewed XCI resulting in the inactivation of the wild-type allele and expression of the mutant allele, harbouring either a single nucleotide or structural variant [37]. To date, 28 balanced translocations disrupting the DMD gene have been reported in female dystrophinopathy patients. All the reported patients showed a DMD-like phenotype with or without cardiomyopathy and cognitive impairment and usually, but not always, present skewed XCI.
In the past, genomic rearrangements such as the one presented here were detected by karyotyping; however, small chromosomal rearrangements could not be detected, and breakpoint characterisation was not performed. The underdiagnosed prevalence of a DMD-like phenotype in female patients often leads to a misdiagnosis and/or diagnostic odyssey of female patients with muscular dystrophy and alterations in the DGC staining, as in the case described here. A secondary reduction in dystrophin expression has been described in sarcoglycanopathies, particularly in patients with pathogenic variants in β-, γ-, and δ-sarcoglycan [38]. Nevertheless, in female patients with high CK levels, a positive Gower’s sign, cardiomyopathy, and no candidate variants on exome sequencing, the presence of chromosomal rearrangements, such as deletions, insertions, inversions, or translocations within the DMD gene, should, be investigated by standard karyotyping or, otherwise, WGS.
As shown here, RNAseq can frequently help in the prioritisation of WGS findings through the detection of aberrantly spliced or expressed genes, or genes with monoallelic expression. However, RNAseq alone would not have sufficed to reach the genetic diagnosis of the patient. A combination of RNAseq and WGS or, alternatively, targeted-long-read sequencing of DMD should be the strategy in those females with a DMD-like phenotype without a clear genetic diagnosis after routine genetic testing. Overall, this case illustrates the diagnostic odyssey in female patients with complex variants truncating the DMD gene that are not detected by exome or panel sequencing.
4. Materials and Methods
Patient. The patient was clinically assessed at the NHS England Highly Specialised Service for Rare Neuromuscular Disorders, UK. Informed consent for medical research was obtained and biological samples were submitted to the Newcastle Medical Research Council (MRC) Centre Biobank for Neuromuscular Diseases for which ethical approval was granted by the National Research Ethics Service (NRES) Committee North East–Newcastle & North Tyneside 1 (reference 19/NE/0028).
Muscle biopsy. At the age of 9 years, a muscle biopsy from the quadriceps was taken.
Immunohistochemical labelling. Labelling was performed on 6 μm sections of frozen skeletal muscle. The sections were equilibrated to room temperature after being removed from storage at −80 °C, before being permeabilized in 0.1% Triton X-100 in PBS (TBS). Samples were then incubated at 4 °C overnight in primary antibodies; dystrophin (N- and C- terminal domains), sarcoglycans (α, β, γ, and δ), dystroglycans (α and β), and utrophin. This was followed by washing in TBS. X-Cell Plus Universal probe was applied, followed by X-Cell Plus Polymer HRP, and developed with Liquid Stable DAB according to the manufacturer’s instructions. Nuclei were counterstained with Carazzi’s haematoxylin.
Western blotting. BioRad Protean II equipment was used to cast two 16 cm gels, 1.5 cm thick. The resolving gel was poured in two phases to give a lower phase of 7% and upper phase of 5.5% polyacrylamide. Once the resolving gel was set, a 1cm deep 3% stacking gel was poured to create the sample lanes. Frozen muscle samples were homogenised by pipette with treatment buffer. The samples were placed in boiling water for 2 min and centrifuged at 1300× g for 5 min. A total of 46 μL of supernatant was loaded and gels were run at 21 mA overnight at 10 °C. The gels were blotted for 7 h onto 0.45 mm nitrocellulose in transfer buffer at −13.5 °C. Immunolabelling was performed by blocking the nitrocellulose in 5% milk powder, incubating in a multiplex cocktail of primary antibodies diluted in Tris buffered saline with 0.05% Tween 20 (TBST), washing in TBST for 15 min, repeat milk blocking for 15 min, washing in TBST for 15 min, incubation in peroxidase-conjugated anti-mouse secondary antibody for 1 h, and visualisation with Supersignal West Pico in accordance with manufacturer’s instructions.
Genetic Studies
DNA. DNA from the proband and relatives was extracted from peripheral blood. MLPA of the DMD gene was performed using the P034 and P035 Salsa Kit (MRC-Holland, Amsterdam, Netherlands). Trio whole exome sequencing (WES) was performed at the Broad Institute of MIT and Harvard’s Genomics Platform (Cambridge, MA, USA), as described previously [39]. Single nucleotide variants (SNVs) and copy number variants (CNVs) were analysed on the seqr platform (seqr.broadinstitute.org). Trio whole genome sequencing (WGS) was performed as part of the Solve-RD project [40] at BGI Europe (Copenhagen, Denmark). Reads were aligned to the GRCh37 reference genome with BWA-MEM v0.7.8. The data were uploaded and analysed on the Genome-Phenome Analysis Platform (GPAP; platform.rd-connect.eu). Variants identified in trio-WES and trio-WGS were interpreted according to the American College of Medical Genetics and Genomics (ACMG) guidelines [41] and filtered by population allele frequency (<0.1%). ClinCNV and Manta v1.6.0 were used to detect CNVs and structural variants, respectively, in the trio-WGS [42,43]. Breakpoints were manually inspected using Integrative Genomics Viewer v2.16.2 (IGV). The RepeatMasker track from the UCSC Genome Browser was used to detect transposable elements or repetitive sequences nearby translocation breakpoints. Two primer pairs were used to validate the breakpoints through PCR and Sanger Sequencing: 5′ AGACCCTTTCTTTCCTGCGT 3′/5′ ACGGATGCTGGGCTCAATAT 3′ and 5′ TCAACAGCACATGTGATTTCAGTC 3′/5′ TGGGCAGCTGTAGTGAACAA 3′.
SNP Array. The patient’s DNA was hybridised in a Infinium CytoSNP-850K v1.3 BeadChip array (Illumina, San Diego, CA, USA), following the manufacturer’s instructions. BlueFuse Multi v4.5 software (Illumina, San Diego, CA, USA) was used to analyse the data with an average backbone resolution of approximately 50kb and an average targeted gene resolution of approximately 10kb. Classification of the copy number variant followed the ACMG guidelines [44].
RNA sequencing. Total RNA was extracted from proband’s muscle biopsy. A Qubit RNA BR Assay kit (Thermo Fisher Scientific, Waltham, MA, USA) was used to assess the quantity and an Agilent Fragment Analyzer DNF-471 RNA Kit (15 nt) (Agilent Technologies, Santa Clara, CA, USA) was used to evaluate integrity. RNA-Seq libraries were generated from total RNA using the TruSeq Stranded mRNA Library Prep Kit (Illumina, San Diego, CA, USA). In brief, mRNA was enriched with oligo-dT magnetic beads from 500 ng of total RNA. The resulting blunt-ended, double-stranded cDNA was 3′ adenylated, and Illumina-platform-compatible adaptors with unique dual indexes and unique molecular identifiers (Integrated DNA Technologies, CA, USA) were ligated. The ligation product underwent enrichment with 15 PCR cycles. The final library was validated using an Agilent Bioanalyzer DNA 7500 assay (Agilent Technologies, Santa Clara, CA, USA). Libraries were sequenced on a NovaSeq 6000 (Illumina, San Diego, CA, USA) in paired-end mode with a read length of 2 × 151 bp, following the manufacturer’s protocol for dual indexing. Image analysis, base calling, and quality scoring were processed using the manufacturer’s software Real Time Analysis (RTA 3.4.4), followed by the generation of FASTQ sequence files. RNA-seq reads were trimmed with trim_galore/0.6.7 [45] and mapped against the hg19 human reference with STAR/2.7.8a using the - -twopassMode = BASIC parameter [46]. Mapping quality metrics were calculated with Qualimap [47], featureCounts [48], STAR log files, and custom scripts. RNAseq data were then analysed using the Detection of RNA Outliers Pipeline (DROP, v.1.3.3) [20]. This analysis was performed within a cohort of 116 muscle samples sequenced and processed using the same pipeline within the Solve-RD project. DROP was also used to verify that the RNA and DNA samples belonged to the same individual.
X-chromosome inactivation assay. DNA from the proband extracted from peripheral blood was subjected to an initial PCR to confirm that the patient was heterozygous for the CAG repeat in exon 1 of the androgen receptor (AR) gene using primers 5’-TCCAGAATCTGTTCCAGAGCGT (forward) and 5′-(FAM) GGCTGTGAAGGTTGCTGTTCCTCAT (reverse). Female DNA was then spiked with male DNA (2:1 ratio) to act as a control for restriction digestion. One-half of the spiked DNA was digested overnight at 37 °C with the methylation-sensitive enzymes HhaI and HpaII. The digested spiked DNA and the undigested spiked DNA were then used as templates in the AR (fluorescent) PCR. A PCR product should only be seen for an inactive (methylated) X chromosome, which cannot be digested by the enzymes, HhaI and HpaII, and can, therefore, act as a template for the PCR. The presence of the male allele in the digested sample suggests that digestion has been compromised. Following capillary electrophoresis, Genemarker (
Conceptualization, A.S.-C. and A.T.; Formal analysis, A.S.-C., V.A.Y., G.D., S.L., A.E.-C., Y.P., E.H. and A.T.; Funding acquisition, J.G., V.S. and A.T.; Investigation, R.M.-L., E.H. and C.M.-B.; Methodology, A.S.-C., V.A.Y., G.D., S.L., A.E.-C., Y.P. and A.T.; Writing—original draft, A.S.-C., Y.P., R.M.-L., E.H. and A.T.; Writing—review and editing, A.S.-C., V.A.Y., G.D., S.L., A.E.-C., J.G., Y.P., R.M.-L., E.H., C.M.-B., V.S. and A.T. All authors have read and agreed to the published version of the manuscript.
This study was conducted in accordance with the Declaration of Helsinki and approved by the National Research Ethics Service (NRES) Committee North East–Newcastle & North Tyneside (reference 19/NE/0028).
Informed consent was obtained from the patient involved in the study.
Further details and data are available upon request from the corresponding author.
We would like to thank the patient and their relatives for their participation in the study. We thank NHS England Highly Specialised Service for Rare Neuromuscular Disorders, UK and Stephen Morris and the cytogenetics facility at Leeds Teaching Hospitals NHS Trust, UK.
The authors declare no conflicts of interest.
Footnotes
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
References
1. Koenig, M.; Monaco, A.P.; Kunkel, L.M. The Complete Sequence of Dystrophin Predicts a Rod-Shaped Cytoskeletal Protein. Cell; 1988; 53, pp. 219-228. [DOI: https://dx.doi.org/10.1016/0092-8674(88)90383-2] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/3282674]
2. Hoffman, E.P.; Brown, R.H.; Kunkel, L.M. Dystrophin: The Protein Product of the Duchenne Muscular Dystrophy Locus. Cell; 1987; 51, pp. 919-928. [DOI: https://dx.doi.org/10.1016/0092-8674(87)90579-4] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/3319190]
3. Muntoni, F.; Torelli, S.; Ferlini, A. Dystrophin and Mutations: One Gene, Several Proteins, Multiple Phenotypes. Lancet Neurol.; 2003; 2, pp. 731-740. [DOI: https://dx.doi.org/10.1016/S1474-4422(03)00585-4]
4. Duan, D.; Goemans, N.; Takeda, S.; Mercuri, E.; Aartsma-Rus, A. Duchenne Muscular Dystrophy. Nat. Rev. Dis. Primers; 2021; 7, 13. [DOI: https://dx.doi.org/10.1038/s41572-021-00248-3] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/33602943]
5. Monaco, A.P.; Bertelson, C.J.; Liechti-Gallati, S.; Moser, H.; Kunkel, L.M. An Explanation for the Phenotypic Differences between Patients Bearing Partial Deletions of the DMD Locus. Genomics; 1988; 2, pp. 90-95. [DOI: https://dx.doi.org/10.1016/0888-7543(88)90113-9]
6. Aartsma-Rus, A.; Van Deutekom, J.C.T.; Fokkema, I.F.; Van Ommen, G.J.B.; Den Dunnen, J.T. Entries in the Leiden Duchenne Muscular Dystrophy Mutation Database: An Overview of Mutation Types and Paradoxical Cases That Confirm the Reading-Frame Rule. Muscle Nerve; 2006; 34, pp. 135-144. [DOI: https://dx.doi.org/10.1002/mus.20586]
7. Norman, A.; Harper, P. A Survey of Manifesting Carriers of Duchenne and Becker Muscular Dystrophy in Wales. Clin. Genet.; 1989; 36, pp. 31-37. [DOI: https://dx.doi.org/10.1111/j.1399-0004.1989.tb03363.x]
8. Politano, L.; Nigro, V.; Nigro, G.; Petretta, V.R.; Passamano, L.; Papparella, S.; Di Somma, S.; Comi, L.I. Development of Cardiomyopathy in Female Carriers of Duchenne and Becker Muscular Dystrophies. JAMA; 1996; 275, pp. 1335-1338. [DOI: https://dx.doi.org/10.1001/jama.1996.03530410049032]
9. Hoogerwaard, E.M.; Bakker, E.; Ippel, P.F.; Oosterwijk, J.C.; Majoor-Krakauer, D.F.; Leschot, N.J.; Van Essen, A.J.; Brunner, H.G.; van der Wouw, P.A.; Wilde, A.A. et al. Signs and Symptoms of Duchenne Muscular Dystrophy and Becker Muscular Dystrophy among Carriers in The Netherlands: A Cohort Study. Lancet; 1999; 353, pp. 2116-2119. [DOI: https://dx.doi.org/10.1016/S0140-6736(98)10028-4]
10. Grain, L.; Cortina-Borja, M.; Forfar, C.; Hilton-Jones, D.; Hopkin, J.; Burch, M. Cardiac Abnormalities and Skeletal Muscle Weakness in Carriers of Duchenne and Becker Muscular Dystrophies and Controls. Neuromuscul. Disord.; 2001; 11, pp. 186-191. [DOI: https://dx.doi.org/10.1016/S0960-8966(00)00185-1]
11. Mccaffrey, T.; Guglieri, M.; Murphy, A.P.; Bushby, K.; Johnson, A.; Bourke, J.P. Cardiac Involvement in Female Carriers of Duchenne or Becker Muscular Dystrophy. Muscle Nerve; 2017; 55, pp. 810-818. [DOI: https://dx.doi.org/10.1002/mus.25445] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/27761893]
12. Ishizaki, M.; Kobayashi, M.; Adachi, K.; Matsumura, T.; Kimura, E. Female Dystrophinopathy: Review of Current Literature. Neuromuscul. Disord.; 2018; 28, pp. 572-581. [DOI: https://dx.doi.org/10.1016/j.nmd.2018.04.005] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/29801751]
13. Soltanzadeh, P.; Friez, M.J.; Dunn, D.; von Niederhausern, A.; Gurvich, O.L.; Swoboda, K.J.; Sampson, J.B.; Pestronk, A.; Connolly, A.M.; Florence, J.M. et al. Clinical and Genetic Characterization of Manifesting Carriers of DMD Mutations. Neuromuscul. Disord.; 2010; 20, pp. 499-504. [DOI: https://dx.doi.org/10.1016/j.nmd.2010.05.010]
14. Juan-Mateu, J.; Rodríguez, M.J.; Nascimento, A.; Jiménez-Mallebrera, C.; González-Quereda, L.; Rivas, E.; Paradas, C.; Madruga, M.; Sánchez-Ayaso, P.; Jou, C. et al. Prognostic Value of X-Chromosome Inactivation in Symptomatic Female Carriers of Dystrophinopathy. Orphanet. J. Rare Dis.; 2012; 7, 82. [DOI: https://dx.doi.org/10.1186/1750-1172-7-82] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/23092449]
15. Brioschi, S.; Gualandi, F.; Scotton, C.; Armaroli, A.; Bovolenta, M.; Falzarano, M.S.; Sabatelli, P.; Selvatici, R.; D’Amico, A.; Pane, M. et al. Genetic Characterization in Symptomatic Female DMD Carriers: Lack of Relationship between X-Inactivation, Transcriptional DMD Allele Balancing and Phenotype. BMC Med. Genet.; 2012; 13, 73. [DOI: https://dx.doi.org/10.1186/1471-2350-13-73]
16. Flanigan, K.M.; Dunn, D.M.; von Niederhausern, A.; Soltanzadeh, P.; Gappmaier, E.; Howard, M.T.; Sampson, J.B.; Mendell, J.R.; Wall, C.; King, W.M. et al. Mutational Spectrum of DMD Mutations in Dystrophinopathy Patients: Application of Modern Diagnostic Techniques to a Large Cohort. Hum. Mutat.; 2009; 30, pp. 1657-1666. [DOI: https://dx.doi.org/10.1002/humu.21114] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/19937601]
17. Tuffery-Giraud, S.; Béroud, C.; Leturcq, F.; Yaou, R.B.; Hamroun, D.; Michel-Calemard, L.; Moizard, M.P.; Bernard, R.; Cossée, M.; Boisseau, P. et al. Genotype-Phenotype Analysis in 2,405 Patients with a Dystrophinopathy Using the UMD-DMD Database: A Model of Nationwide Knowledgebase. Hum. Mutat.; 2009; 30, pp. 934-945. [DOI: https://dx.doi.org/10.1002/humu.20976]
18. Bladen, C.L.; Salgado, D.; Monges, S.; Foncuberta, M.E.; Kekou, K.; Kosma, K.; Dawkins, H.; Lamont, L.; Roy, A.J.; Chamova, T. et al. The TREAT-NMD DMD Global Database: Analysis of More than 7,000 Duchenne Muscular Dystrophy Mutations. Hum. Mutat.; 2015; 36, pp. 395-402. [DOI: https://dx.doi.org/10.1002/humu.22758]
19. Töpf, A.; Johnson, K.; Bates, A.; Phillips, L.; Chao, K.R.; England, E.M.; Laricchia, K.M.; Mullen, T.; Valkanas, E.; Xu, L. et al. Sequential Targeted Exome Sequencing of 1001 Patients Affected by Unexplained Limb-Girdle Weakness. Genet. Med.; 2020; 22, pp. 1478-1488. [DOI: https://dx.doi.org/10.1038/s41436-020-0840-3]
20. Yépez, V.A.; Mertes, C.; Müller, M.F.; Klaproth-Andrade, D.; Wachutka, L.; Frésard, L.; Gusic, M.; Scheller, I.F.; Goldberg, P.F.; Prokisch, H. et al. Detection of Aberrant Gene Expression Events in RNA Sequencing Data. Nat. Protoc.; 2021; 16, pp. 1276-1296. [DOI: https://dx.doi.org/10.1038/s41596-020-00462-5]
21. Brechtmann, F.; Mertes, C.; Matusevičiūtė, A.; Yépez, V.A.; Avsec, Ž.; Herzog, M.; Bader, D.M.; Prokisch, H.; Gagneur, J. OUTRIDER: A Statistical Method for Detecting Aberrantly Expressed Genes in RNA Sequencing Data. Am. J. Hum. Genet.; 2018; 103, pp. 907-917. [DOI: https://dx.doi.org/10.1016/j.ajhg.2018.10.025]
22. Nallamilli, B.R.R.; Chaubey, A.; Valencia, C.A.; Stansberry, L.; Behlmann, A.M.; Ma, Z.; Mathur, A.; Shenoy, S.; Ganapathy, V.; Jagannathan, L. et al. A Single NGS-Based Assay Covering the Entire Genomic Sequence of the DMD Gene Facilitates Diagnostic and Newborn Screening Confirmatory Testing. Hum. Mutat.; 2021; 42, pp. 626-638. [DOI: https://dx.doi.org/10.1002/humu.24191]
23. Gonzalez-Quereda, L.; Rodriguez, M.J.; Diaz-Manera, J.; Alonso-Perez, J.; Gallardo, E.; Nascimento, A.; Ortez, C.; Natera-De Benito, D.; Olive, M.; Gonzalez-Mera, L. et al. Targeted Next-Generation Sequencing in a Large Cohort of Genetically Undiagnosed Patients with Neuromuscular Disorders in Spain. Genes; 2020; 11, 539. [DOI: https://dx.doi.org/10.3390/genes11050539] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/32403337]
24. Cummings, B.B.; Marshall, J.L.; Tukiainen, T.; Lek, M.; Donkervoort, S.; Foley, A.R.; Bolduc, V.; Waddell, L.B.; Sandaradura, S.A.; O’Grady, G.L. et al. Improving Genetic Diagnosis in Mendelian Disease with Transcriptome Sequencing. Sci. Transl. Med.; 2017; 9, eaal5209. [DOI: https://dx.doi.org/10.1126/scitranslmed.aal5209] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/28424332]
25. Gonorazky, H.D.; Naumenko, S.; Ramani, A.K.; Nelakuditi, V.; Mashouri, P.; Wang, P.; Kao, D.; Ohri, K.; Viththiyapaskaran, S.; Tarnopolsky, M.A. et al. Expanding the Boundaries of RNA Sequencing as a Diagnostic Tool for Rare Mendelian Disease. Am. J. Hum. Genet.; 2019; 104, pp. 466-483. [DOI: https://dx.doi.org/10.1016/j.ajhg.2019.01.012]
26. Yépez, V.A.; Gusic, M.; Kopajtich, R.; Mertes, C.; Smith, N.H.; Alston, C.L.; Ban, R.; Beblo, S.; Berutti, R.; Blessing, H. et al. Clinical Implementation of RNA Sequencing for Mendelian Disease Diagnostics. Genome Med.; 2022; 14, 38. [DOI: https://dx.doi.org/10.1186/s13073-022-01019-9] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/35379322]
27. Hengel, H.; Hannan, S.B.; Dyack, S.; MacKay, S.B.; Schatz, U.; Fleger, M.; Kurringer, A.; Balousha, G.; Ghanim, Z.; Alkuraya, F.S. et al. Bi-Allelic Loss-of-Function Variants in BCAS3 Cause a Syndromic Neurodevelopmental Disorder. Am. J. Hum. Genet.; 2021; 108, pp. 1069-1082. [DOI: https://dx.doi.org/10.1016/j.ajhg.2021.04.024] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/34022130]
28. Luce, L.; Abelleyro, M.M.; Carcione, M.; Mazzanti, C.; Rossetti, L.; Radic, P.; Szijan, I.; Menazzi, S.; Francipane, L.; Nevado, J. et al. Analysis of Complex Structural Variants in the DMD Gene in One Family. Neuromuscul. Disord.; 2021; 31, pp. 253-263. [DOI: https://dx.doi.org/10.1016/j.nmd.2020.11.015]
29. Ankala, A.; Kohn, J.N.; Hegde, A.; Meka, A.; Ephrem, C.L.H.; Askree, S.H.; Bhide, S.; Hegde, M.R. Aberrant Firing of Replication Origins Potentially Explains Intragenic Nonrecurrent Rearrangements within Genes, Including the Human DMD Gene. Genome Res.; 2012; 22, pp. 25-34. [DOI: https://dx.doi.org/10.1101/gr.123463.111] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/22090376]
30. Baskin, B.; Gibson, W.T.; Ray, P.N. Duchenne Muscular Dystrophy Caused by a Complex Rearrangement between Intron 43 of the DMD Gene and Chromosome 4. Neuromuscul. Disord.; 2011; 21, pp. 178-182. [DOI: https://dx.doi.org/10.1016/j.nmd.2010.11.008]
31. Chang, H.H.Y.; Pannunzio, N.R.; Adachi, N.; Lieber, M.R. Non-Homologous DNA End Joining and Alternative Pathways to Double-Strand Break Repair. Nat. Rev. Mol. Cell Biol.; 2017; 18, pp. 495-506. [DOI: https://dx.doi.org/10.1038/nrm.2017.48] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/28512351]
32. Schluth-Bolard, C.; Diguet, F.; Chatron, N.; Rollat-Farnier, P.-A.; Bardel, C.; Afenjar, A.; Amblard, F.; Amiel, J.; Blesson, S.; Callier, P. et al. Whole Genome Paired-End Sequencing Elucidates Functional and Phenotypic Consequences of Balanced Chromosomal Rearrangement in Patients with Developmental Disorders. J. Med. Genet.; 2019; 56, pp. 526-535. [DOI: https://dx.doi.org/10.1136/jmedgenet-2018-105778] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/30923172]
33. Nilsson, D.; Pettersson, M.; Gustavsson, P.; Förster, A.; Hofmeister, W.; Wincent, J.; Zachariadis, V.; Anderlid, B.-M.; Nordgren, A.; Mäkitie, O. et al. Whole-Genome Sequencing of Cytogenetically Balanced Chromosome Translocations Identifies Potentially Pathological Gene Disruptions and Highlights the Importance of Microhomology in the Mechanism of Formation. Hum. Mutat.; 2017; 38, pp. 180-192. [DOI: https://dx.doi.org/10.1002/humu.23146] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/27862604]
34. Fujii, K.; Minami, N.; Hayashi, Y.; Nishino, I.; Nonaka, I.; Tanabe, Y.; Takanashi, J.; Kohno, Y. Homozygous Female Becker Muscular Dystrophy. Am. J. Med. Genet. A; 2009; 149A, pp. 1052-1055. [DOI: https://dx.doi.org/10.1002/ajmg.a.32808] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/19396825]
35. Quan, F.; Janas, J.; Toth-Fejel, S.; Johnson, D.B.; Wolford, J.K.; Popovich, B.W. Uniparental Disomy of the Entire X Chromosome in a Female with Duchenne Muscular Dystrophy. Am. J. Hum. Genet.; 1997; 60, pp. 160-165.
36. Baiget, M.; Tizzano, E.; Volpini, V.; del Rio, E.; Pérez-Vidal, T.; Gallano, P. DMD Carrier Detection in a Female with Mosaic Turner’s Syndrome. J. Med. Genet.; 1991; 28, pp. 209-210. [DOI: https://dx.doi.org/10.1136/jmg.28.3.209]
37. Boyd, Y.; Buckle, V.; Holt, S.; Munro, E.; Hunter, D.; Craig, I. Muscular Dystrophy in Girls with X;Autosome Translocations. J. Med. Genet.; 1986; 23, pp. 484-490. [DOI: https://dx.doi.org/10.1136/jmg.23.6.484] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/3806636]
38. Vainzof, M.; Passos-Bueno, M.R.; Canovas, M.; Moreira, E.S.; Pavanello, R.C.; Marie, S.K.; Anderson, L.V.; Bonnemann, C.G.; McNally, E.M.; Nigro, V. et al. The Sarcoglycan Complex in the Six Autosomal Recessive Limb-Girdle Muscular Dystrophies. Hum. Mol. Genet.; 1996; 5, pp. 1963-1969. [DOI: https://dx.doi.org/10.1093/hmg/5.12.1963]
39. Harris, E.; Topf, A.; Barresi, R.; Hudson, J.; Powell, H.; Tellez, J.; Hicks, D.; Porter, A.; Bertoli, M.; Evangelista, T. et al. Exome Sequences versus Sequential Gene Testing in the UK Highly Specialised Service for Limb Girdle Muscular Dystrophy. Orphanet. J. Rare Dis.; 2017; 12, 151. [DOI: https://dx.doi.org/10.1186/s13023-017-0699-9]
40. Zurek, B.; Ellwanger, K.; Vissers, L.E.L.M.; Schüle, R.; Synofzik, M.; Töpf, A.; de Voer, R.M.; Laurie, S.; Matalonga, L.; Gilissen, C. et al. Solve-RD: Systematic Pan-European Data Sharing and Collaborative Analysis to Solve Rare Diseases. Eur. J. Hum. Genet.; 2021; 29, pp. 1325-1331. [DOI: https://dx.doi.org/10.1038/s41431-021-00859-0] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/34075208]
41. Richards, S.; Aziz, N.; Bale, S.; Bick, D.; Das, S.; Gastier-Foster, J.; Grody, W.W.; Hegde, M.; Lyon, E.; Spector, E. 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, pp. 405-424. [DOI: https://dx.doi.org/10.1038/gim.2015.30] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/25741868]
42. Chen, X.; Schulz-Trieglaff, O.; Shaw, R.; Barnes, B.; Schlesinger, F.; Källberg, M.; Cox, A.J.; Kruglyak, S.; Saunders, C.T. Manta: Rapid Detection of Structural Variants and Indels for Germline and Cancer Sequencing Applications. Bioinformatics; 2016; 32, pp. 1220-1222. [DOI: https://dx.doi.org/10.1093/bioinformatics/btv710] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/26647377]
43. Demidov, G.; Sturm, M.; Ossowski, S. ClinCNV: Multi-Sample Germline CNV Detection in NGS Data. BioRxiv; 2022; 2022-06. [DOI: https://dx.doi.org/10.1101/2022.06.10.495642]
44. Riggs, E.R.; Andersen, E.F.; Cherry, A.M.; Kantarci, S.; Kearney, H.; Patel, A.; Raca, G.; Ritter, D.I.; South, S.T.; Thorland, E.C. et al. Technical Standards for the Interpretation and Reporting of Constitutional Copy-Number Variants: A Joint Consensus Recommendation of the American College of Medical Genetics and Genomics (ACMG) and the Clinical Genome Resource (ClinGen). Genet. Med.; 2020; 22, pp. 245-257. [DOI: https://dx.doi.org/10.1038/s41436-019-0686-8] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/31690835]
45. Felix Krueger Trim Galore. Available online: https://github.com/FelixKrueger/TrimGalore (accessed on 2 February 2023).
46. Dobin, A.; Davis, C.A.; Schlesinger, F.; Drenkow, J.; Zaleski, C.; Jha, S.; Batut, P.; Chaisson, M.; Gingeras, T.R. STAR: Ultrafast Universal RNA-Seq Aligner. Bioinformatics; 2013; 29, pp. 15-21. [DOI: https://dx.doi.org/10.1093/bioinformatics/bts635]
47. Okonechnikov, K.; Conesa, A.; García-Alcalde, F. Qualimap 2: Advanced Multi-Sample Quality Control for High-Throughput Sequencing Data. Bioinformatics; 2016; 32, pp. 292-294. [DOI: https://dx.doi.org/10.1093/bioinformatics/btv566] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/26428292]
48. Liao, Y.; Smyth, G.K.; Shi, W. featureCounts: An Efficient General Purpose Program for Assigning Sequence Reads to Genomic Features. Bioinformatics; 2014; 30, pp. 923-930. [DOI: https://dx.doi.org/10.1093/bioinformatics/btt656]
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer
© 2024 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/). Notwithstanding the ProQuest Terms and Conditions, you may use this content in accordance with the terms of the License.
Abstract
Duchenne and Becker muscular dystrophies, caused by pathogenic variants in DMD, are the most common inherited neuromuscular conditions in childhood. These diseases follow an X-linked recessive inheritance pattern, and mainly males are affected. The most prevalent pathogenic variants in the DMD gene are copy number variants (CNVs), and most patients achieve their genetic diagnosis through Multiplex Ligation-dependent Probe Amplification (MLPA) or exome sequencing. Here, we investigated a female patient presenting with muscular dystrophy who remained genetically undiagnosed after MLPA and exome sequencing. RNA sequencing (RNAseq) from the patient’s muscle biopsy identified an 85% reduction in DMD expression compared to 116 muscle samples included in the cohort. A de novo balanced translocation between chromosome 17 and the X chromosome (t(X;17)(p21.1;q23.2)) disrupting the DMD and BCAS3 genes was identified through trio whole genome sequencing (WGS). The combined analysis of RNAseq and WGS played a crucial role in the detection and characterisation of the disease-causing variant in this patient, who had been undiagnosed for over two decades. This case illustrates the diagnostic odyssey of female DMD patients with complex structural variants that are not detected by current panel or exome sequencing analysis.
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer
Details






1 John Walton Muscular Dystrophy Research Centre, Translational and Clinical Research Institute, Newcastle University and Newcastle Hospitals NHS Foundation Trust, Newcastle upon Tyne NE1 3BZ, UK;
2 School of Computation, Information and Technology, Technical University of Munich, 85748 Garching, Germany
3 Universitätsklinikum Tübingen—Institut für Medizinische Genetik und angewandte Genomik, 72076 Tübingen, Germany
4 Centro Nacional de Análisis Genómico (CNAG), 08028 Barcelona, Spain
5 Centro Nacional de Análisis Genómico (CNAG), 08028 Barcelona, Spain; Universitat de Barcelona (UB), 08007 Barcelona, Spain
6 School of Computation, Information and Technology, Technical University of Munich, 85748 Garching, Germany; Institute of Human Genetics, School of Medicine, Technical University of Munich, 81675 Munich, Germany; Computational Health Center, Helmholtz Center Munich, 85764 Neuherberg, Germany
7 Muscle Immunoanalysis Unit, Newcastle upon Tyne Hospitals NHS Foundation Trust, Newcastle upon Tyne NE7 7DN, UK
8 John Walton Muscular Dystrophy Research Centre, Translational and Clinical Research Institute, Newcastle University and Newcastle Hospitals NHS Foundation Trust, Newcastle upon Tyne NE1 3BZ, UK;