Relevance of sequence and structure elements for deletion events in the dystrophin gene major hot-spot

Non Random Distribution of DMD Deletion Breakpoints and Implication of Double Strand Breaks Repair and Replication Error Repair Mechanisms

BACKGROUND

Dystrophinopathies are mostly caused by copy number variations, especially deletions, in the dystrophin gene (DMD). Despite the large size of the gene, deletions do not occur randomly but mainly in two hot spots, the main one involving exons 45 to 55. The underlying mechanisms are complex and implicate two main mechanisms: Non-homologous end joining (NHEJ) and micro-homology mediated replication-dependent recombination (MMRDR).

OBJECTIVE

Our goals were to assess the distribution of intronic breakpoints (BPs) in the genomic sequence of the main hot spot of deletions within DMD gene and to search for specific sequences at or near to BPs that might promote BP occurrence or be associated with DNA break repair.

METHODS

Using comparative genomic hybridization microarray, 57 deletions within the intron 44 to 55 region were mapped. Moreover, 21 junction fragments were sequenced to search for specific sequences.

RESULTS

Non-randomly distributed BPs were found in introns 44, 47, 48, 49 and 53 and 50% of BPs clustered within genomic regions of less than 700bp. Repeated elements (REs), known to promote gene rearrangement via several mechanisms, were present in the vicinity of 90% of clustered BPs and less frequently (72%) close to scattered BPs, illustrating the important role of such elements in the occurrence of DMD deletions. Palindromic and TTTAAA sequences, which also promote DNA instability, were identified at fragment junctions in 20% and 5% of cases, respectively. Micro-homologies (76%) and insertions or deletions of small sequences were frequently found at BP junctions.

CONCLUSIONS

Our results illustrate, in a large series of patients, the important role of RE and other genomic features in DNA breaks, and the involvement of different mechanisms in DMD gene deletions: Mainly replication error repair mechanisms, but also NHEJ and potentially aberrant firing of replication origins. A combination of these mechanisms may also be possible.

Precise mapping of 17 deletion breakpoints within the central hotspot deletion region (introns 50 and 51) of the DMD gene

Exon deletions in the human DMD gene, which encodes the dystrophin protein, are the molecular defect in 50-70% of cases of Duchenne/Becker muscular dystrophies. Deletions are preferentially clustered in the 5' (exons 2-20) and the central (exons 45-53) region of DMD, likely because local DNA structure predisposes to specific breakage or recombination events. Notably, innovative therapeutic strategies may rescue dystrophin function by homology-based specific targeting of sequences within the central DMD hot spot deletion region. To further study molecular mechanisms that generate such frequent genome variations and to identify residual intronic sequences, we sequenced 17 deletion breakpoints within introns 50 and 51 of DMD and analyzed the surrounding genomic architecture. There was no breakpoint clustering within the introns nor extensive homology between sequences adjacent to each junction. However, at or near the breakpoint, we found microhomology, short tandem repeats, interspersed repeat elements and short sequence stretches that predispose to DNA deletion or bending. Identification of such structural elements contributes to elucidate general mechanisms generating deletion within the DMD gene. Moreover, precise mapping of deletion breakpoints and localization of repeated elements are of interest, because residual intronic sequences may be targeted by therapeutic strategies based on genome editing correction.

Genome Editing Gene Therapy for Duchenne Muscular Dystrophy

Duchenne muscular dystrophy (DMD) is a severe genetic disorder caused by loss of function of the dystrophin gene on the X chromosome. Gene augmentation of dystrophin is challenging due to the large size of the dystrophin cDNA. Emerging genome editing technologies, such as TALEN and CRISPR-Cas9 systems, open a new erain the restoration of functional dystrophin and are a hallmark of bona fide gene therapy. In this review, we summarize current genome editing approaches, properties of target cell types for ex vivo gene therapy, and perspectives of in vivo gene therapy including genome editing in human zygotes. Although technical challenges, such as efficacy, accuracy, and delivery of the genome editing components, remain to be further improved, yet genome editing technologies offer a new avenue for the gene therapy of DMD.

Preimplantational genetic diagnosis and mutation detection in a family with duplication mutation of DMD gene

Duchenne muscular dystrophy (DMD) is an X-linked recessive neuromuscular disease caused by mutation in the DMD gene. A 38-year-old woman was referred to our hospital with her son who was diagnosed with DMD. Multiplex PCR failed to detect DMD mutations in the affected child. The female carrier underwent preimplantation genetic diagnosis by linkage analysis and gender determination. Eight embryos were biopsied after in vitro fertilization. Two healthy embryos determined both as females (E1 and E3) were transferred. Although the paternal allele was absent in E3, it was considered to be a result of allele dropout for the STR-49 marker. Surprisingly, a female and a male baby were delivered at 38 gestational weeks, suggesting that E3 was a male embryo with the allele dropout occurring at the SRY gene. Exon 2 duplication was detected in the DMD patient and the carrier mother using next-generation sequencing and multiple ligation-dependent probe amplification. Next, we verified the duplication of exon 2 by real-time PCR, using a special primer at 3' of intron 1, very close to exon 2. Finally, we confirmed that both newborns inherited the normal allele, using quantitative real-time PCR and linkage analysis.

Detection of truncated dystrophin lacking the C-terminal domain in a Chinese pedigree by next-generation sequencing

Dystrophin (DMD) gene is the largest gene containing 79 exons involving various mutation types and regions, and targeted next-generation sequencing (NGS) was employed in detecting DMD gene mutation in the present study. A literature-annotated disease nonsense mutation (c.10141C>T, NM_004006.1) in exon 70 that has been reported as Duchenne Muscular Dystrophy (DMD)-causing mutation was found in our two patients, the proband and his cousin. In the present study two main methods were used, the next-generation sequencing and the classic Sanger sequencing. The exon capture followed by HiSeq2000 sequencing was specifically used in this study. Combined applications of the next-generation sequencing platform and bioinformatics are proved to be effective methods for DMD diagnosis.

Mechanisms of genomic instabilities underlying two common fragile-site-associated loci, PARK2 and DMD, in germ cell and cancer cell lines

Common fragile sites (CFSs) are specific chromosome regions that exhibit an increased frequency of breaks when cells are exposed to a DNA-replication inhibitor such as aphidicolin. PARK2 and DMD, the causative genes for autosomal-recessive juvenile Parkinsonism and Duchenne and Becker muscular dystrophy, respectively, are two very large genes that are located within aphidicolin-induced CFSs. Gross rearrangements within these two genes are frequently observed as the causative mutations for these diseases, and similar alterations within the large fragile sites that surround these genes are frequently observed in cancer cells. To elucidate the molecular mechanisms underlying this fragility, we performed a custom-designed high-density comparative genomic hybridization analysis to determine the junction sequences of approximately 500 breakpoints in germ cell lines and cancer cell lines involving PARK2 or DMD. The sequence signatures where these breakpoints occur share some similar features both in germ cell lines and in cancer cell lines. Detailed analyses of these structures revealed that microhomologies are predominantly involved in rearrangement processes. Furthermore, breakpoint-clustering regions coincide with the latest-replicating region and with large nuclear-lamina-associated domains and are flanked by the highest-flexibility peaks and R/G band boundaries, suggesting that factors affecting replication timing collectively contribute to the vulnerability for rearrangement in both germ cell and somatic cell lines.

Characterization of deletion breakpoints in patients with dystrophinopathy carrying a deletion of exons 45-55 of the Duchenne muscular dystrophy (DMD) gene

Deletion of exons 45-55 (del45-55) in the Duchenne muscular dystrophy gene (DMD) has gained particular interest in the field of molecular therapy, because it causes a milder phenotype than DMD, and therefore, may represent a good candidate for the goal of a multiple exon-skipping strategy. We have precisely characterized deletion breakpoints in three patients with del45-55 in DMD. Two of them were young adult males of the X-linked dilated cardiomyopathy phenotype, and the third patient revealed the mild Becker muscular dystrophy phenotype of late onset. The deletion breakpoints differed among patients. The deletion started at nt 226 604, 231 518, 117 284 in intron 44, and ended at nt 64 994, 59 314, 71 806 in intron 55, respectively. Deletion junctions showed no significant homology between the sequences adjacent to the distal and proximal end joints in these patients. Deletion breakpoints were not primarily associated with any particular sequence element, or with a matrix attachment region. However, there were several palindromic sequences and short tandem repeats at or near the breakpoints. These sequences, with a marked propensity to form secondary DNA structure intermediates, may predispose local DNA to breakage and intragenic recombination in these patients.

Generation and characterization of transgenic mice with the full-length human DMD gene

We report the generation of mice with an intact and functional copy of the 2.3-megabase human dystrophin gene (hDMD), the largest functional stretch of human DNA thus far integrated into a mouse chromosome. Yeast spheroplasts containing an artificial chromosome with the full-length hDMD gene were fused with mouse embryonic stem cells and were subsequently injected into mouse blastocysts to produce transgenic hDMD mice. Human-specific PCR, Southern blotting, and fluorescent in situ hybridization techniques demonstrated the intactness and stable chromosomal integration of the hDMD gene on mouse chromosome 5. Expression of the transgene was confirmed by RT-PCR and Western blotting. The tissue-specific expression pattern of the different DMD transcripts was maintained. However, the human Dp427p and Dp427m transcripts were expressed at 2-fold higher levels and human Dp427c and Dp260 transcripts were expressed at 2- and 4-fold lower levels than their endogenous counterparts. Ultimate functional proof of the hDMD transgene was obtained by crossing of hDMD mice with dystrophin-deficient mdx mice and dystrophin and utrophin-deficient mdx x Utrn-/- mice. The hDMD transgene rescued the lethal dystrophic phenotype of the mdx x Utrn-/- mice. All signs of muscular dystrophy disappeared in the rescued mice, as demonstrated by histological staining of muscle sections and gene expression profiling experiments. Currently, hDMD mice are extensively used for preclinical testing of sequence-specific therapeutics for the treatment of Duchenne muscular dystrophy. In addition, the hDMD mouse can be used to study the influence of the genomic context on deletion and recombination frequencies, genome stability, and gene expression regulation.

Copy number variation in the genome; the human DMD gene as an example

Recent developments have yielded new technologies that have greatly simplified the detection of deletions and duplications, i.e., copy number variants (CNVs). These technologies can be used to screen for CNVs in and around specific genomic regions, as well as genome-wide. Several genome-wide studies have demonstrated that CNV in the human genome is widespread and may include millions of nucleotides. One of the questions that emerge is which sequences, structures and/or processes are involved in their generation. Using as an example the human DMD gene, mutations in which cause Duchenne and Becker muscular dystrophy, we review the current data, determine the deletion and duplication profile across the gene and summarize the information that has been collected regarding their origin. In addition we discuss the methods most frequently used for their detection, in particular MAPH and MLPA.

A region in the dystrophin gene major hot spot harbors a cluster of deletion breakpoints and generates double-strand breaks in yeast

Deletions within the dystrophin gene (DMD) account for >70% of mutations leading to Duchenne and Becker muscular dystrophies (DMD and BMD). Deletion breakpoints were reported to be scattered within regions that also represent meiotic recombination hot spots. Recent studies indicates that deletion junctions arise from nonhomologous end joining (NHEJ), a major pathway for repairing DNA double-strand breaks (DSBs) in mammals. Here we show that a region in intron 47 (i.e., a major deletion hot spot in the DMD gene) generates DSBs during meiosis in yeast and harbors a cluster of previously sequenced deletion breaks. Mapping of breakpoints in 26 BMD/DMD patients indicated that the frequency of breakpoint occurrence around this region is 3-fold higher than expected by chance. These findings suggest that DSBs mediate deletion formation in intron 47 and possibly account for the high frequency of meiotic recombination in the region. Statistical analysis indicated the presence of at least one other breakpoint cluster in intron 47. Taken together, these results suggest that the primary events in deletion formation occur within discrete regions and that the scattered breakpoint distribution reflects both a variable degree of DSB end processing and the availability of a small (compared to the huge regions involved) deletion junction sample.

Gross genomic rearrangements involving deletions in the CFTR gene: characterization of six new events from a large cohort of hitherto unidentified cystic fibrosis chromosomes and meta-analysis of the underlying mechanisms

Gross genomic rearrangements involving deletions in the CFTR gene have recently been found to account for approximately 20% of unidentified cystic fibrosis (CF) chromosomes in both French and Italian patients. Using QMPSF and walking quantitative DHPLC, six novel mutations (three simple deletions, two complex deletions with short insertions of 3-6 bp, and a complex deletion with a 182 bp inverted downstream sequence) were characterized by screening 274 unidentified CF chromosomes from 10 different countries. These lesions increase the total number of fully characterized large CFTR genomic rearrangements involving deletions to 21. Systematic analysis of the 42 associated breakpoints indicated that all 21 events were caused by nonhomologous recombination. Whole gene complexity analysis revealed a significant correlation between regions of low sequence complexity and the locations of the deletion breakpoints. Known recombination-promoting motifs were noted in the vicinity of the breakpoints. A total of 11 simple deletions were potentially explicable in terms of the classical model of replication slippage. However, the complex deletions appear to have arisen via multiple mechanisms; three of the five complex deletions with short insertions and both examples of large inverted insertions (299 and 182 bp, respectively) can be explained by either a model of serial replication slippage in cis (SRScis) or SRS in trans (SRStrans). Finally, the nature and distribution of large genomic rearrangements in the CFTR gene were compared and contrasted with those of two other genes, DMD and MSH2, with a view to gaining a broader understanding of DNA sequence context in mediating the diverse underlying mutational mechanisms.

Splicing analysis disclosed a determinant single nucleotide for exon skipping caused by a novel intraexonic four-nucleotide deletion in the dystrophin gene

BACKGROUND

Mutations in exonic splicing enhancer sequences are known to cause splicing errors. Although exonic splicing enhancers have been identified as a stretch of purine-rich sequences, it has been difficult to precisely pinpoint the determinant nucleotides in these sequences. This article reports that a 4-bp deletion in exon 38 of the dystrophin gene induced complete exon 38 skipping in vivo. Moreover, the third nucleotide of the deletion was shown to be determinant for the exonic splicing enhancer activity in in vivo splicing analysis of hybrid minigenes encoding mutant exons.

METHOD

Genomic DNA analysis of a 2-year-old boy with a raised level of serum creatine kinase yielded a 4-bp deletion 11 bp upstream of the 3' end of exon 38 of the dystrophin gene (c. 5434-5437del TTCA), disrupting a predicted SC35-binding site.

RESULT

Interestingly, his dystrophin mRNA was shown to completely lack exon 38 (exon 38- transcript). As the exon 38- transcript coded for a truncated dystrophin protein, this exon skipping was determined to be a modifying factor of his phenotype. In an in vivo splicing assay, a hybrid minigene encoding exon 38 with the 4-bp deletion was shown to induce complete exon 38 skipping, confirming the deleted region as a splicing enhancer sequence. Site-directed mutagenesis of the deleted sequence showed that the complete exon 38 skipping was caused by mutation of the third nucleotide position of the deletion (C5436), whereas mutations at the other three nucleotide positions induced partial exon skipping.

CONCLUSION

Our results underline the potential of understanding the regulation of exonic splicing enhancer sequences and exon skipping therapy for treatment of Duchenne's muscular dystrophy.

Double-strand DNA break repair with replication slippage on two strands: a novel mechanism of deletion formation

We have characterized an unusual family with two different androgen receptor (AR) gene deletions, in which we propose a novel mechanism of deletion formation has occurred. Affected individuals have the X-linked disorder androgen insensitivity syndrome, and we previously showed that different family members have deletions of different exons of the AR gene. We have now fully sequenced the deletions from affected individuals, and confirmed the presence of different deletions in different affected family members. Most affected and heterozygote individuals have a 4,430-bp deletion of exon 5 that occurred between repeated GTGGCAT motifs in introns 4 and 5. One affected hemizygous individual has a 4,033-bp deletion of exons 6 and 7 that occurred between repeated CCTC motifs in introns 5 and 7. The intron 5 breakpoint junctions of the two deletions are only 11 bp apart. Surprisingly, the maternal grandmother of the original index case was found to be mosaic for both deletional events, as well as having the normal AR gene. Karyotyping ruled out 47,XXX trisomy, indicating triple mosaicism for the two different deleted AR alleles and a normal AR allele. This triple mosaicism must have occurred early in embryonic development, as both deletions were passed on to different children. Based on these findings, we propose a novel mechanism of deletion formation. We suggest that during AR gene replication, a double strand DNA break occurred in intron 5, and that a variant of replication slippage occurred on both newly synthesized strands between the repeat motifs of microhomology, leading to the formation of the two different AR gene deletions.

Heterogeneous duplications in patients with Pelizaeus-Merzbacher disease suggest a mechanism of coupled homologous and nonhomologous recombination

We describe genomic structures of 59 X-chromosome segmental duplications that include the proteolipid protein 1 gene (PLP1) in patients with Pelizaeus-Merzbacher disease. We provide the first report of 13 junction sequences, which gives insight into underlying mechanisms. Although proximal breakpoints were highly variable, distal breakpoints tended to cluster around low-copy repeats (LCRs) (50% of distal breakpoints), and each duplication event appeared to be unique (100 kb to 4.6 Mb in size). Sequence analysis of the junctions revealed no large homologous regions between proximal and distal breakpoints. Most junctions had microhomology of 1-6 bases, and one had a 2-base insertion. Boundaries between single-copy and duplicated DNA were identical to the reference genomic sequence in all patients investigated. Taken together, these data suggest that the tandem duplications are formed by a coupled homologous and nonhomologous recombination mechanism. We suggest repair of a double-stranded break (DSB) by one-sided homologous strand invasion of a sister chromatid, followed by DNA synthesis and nonhomologous end joining with the other end of the break. This is in contrast to other genomic disorders that have recurrent rearrangements formed by nonallelic homologous recombination between LCRs. Interspersed repetitive elements (Alu elements, long interspersed nuclear elements, and long terminal repeats) were found at 18 of the 26 breakpoint sequences studied. No specific motif that may predispose to DSBs was revealed, but single or alternating tracts of purines and pyrimidines that may cause secondary structures were common. Analysis of the 2-Mb region susceptible to duplications identified proximal-specific repeats and distal LCRs in addition to the previously reported ones, suggesting that the unique genomic architecture may have a role in nonrecurrent rearrangements by promoting instability.

Molecular characterization of the spectrum of genomic deletions in the mismatch repair genes MSH2, MLH1, MSH6, and PMS2 responsible for hereditary nonpolyposis colorectal cancer (HNPCC)

A systematic search by Southern blot analysis in a cohort of 439 hereditary nonpolyposis colorectal cancer (HNPCC) families for genomic rearrangements in the main mismatch repair (MMR) genes, namely, MSH2, MLH1, MSH6, and PMS2, identified 48 genomic rearrangements causative of this inherited predisposition to colorectal cancer in 68 unrelated kindreds. Twenty-nine of the 48 rearrangements were found in MSH2, 13 in MLH1, 2 in MSH6, and 4 in PMS2. The vast majority were deletions, although one previously described large inversion, an intronic insertion, and a more complex rearrangement also were found. Twenty-four deletion breakpoints have been identified and sequenced in order to determine the underlying recombination mechanisms. Most fall within repetitive sequences, mainly Alu repeats, in agreement with the differential distribution of deletions between the MSH2 and MLH1 genes: the higher number and density of Alu repeats in MSH2 corresponded with a higher incidence of genomic rearrangement at this disease locus when compared with other MMR genes. Long interspersed nuclear element (LINE) repeats, relatively abundant in, for example, MLH1, did not seem to contribute to the genesis of the deletions, presumably because of their older evolutionary age and divergence among individual repeat units when compared with short interspersed nuclear element (SINE) repeats, including Alu repeats. Moreover, Southern blot analysis of the introns and the genomic regions flanking the MMR genes allowed us to detect 6 novel genomic rearrangements that left the coding region of the disease-causing gene intact. These rearrangements comprised 4 deletions upstream of the coding region of MSH2 (3 cases) and MSH6 (1 case), a 2-kb insertion in intron 7 of PMS2, and a small (459-bp) deletion in intron 13 of MLH1. The characterization of these genomic rearrangements underlines the importance of genomic deletions in the etiology of HNPCC and will facilitate the development of PCR-based tests for their detection in diagnostic laboratories.

Molecular structure of double-minute chromosomes bearing amplified copies of the epidermal growth factor receptor gene in gliomas

Amplification of the epidermal growth factor receptor gene on double minutes is recurrently observed in cells of advanced gliomas, but the structure of these extrachromosomal circular DNA molecules and the mechanisms responsible for their formation are still poorly understood. By using quantitative PCR and chromosome walking, we investigated the genetic content and the organization of the repeats in the double minutes of seven gliomas. It was established that all of the amplicons of a given tumor derive from a single founding extrachromosomal DNA molecule. In each of these gliomas, the founding molecule was generated by a simple event that circularizes a chromosome fragment overlapping the epidermal growth factor receptor gene. In all cases, the fusion of the two ends of this initial amplicon resulted from microhomology-based nonhomologous end-joining. Furthermore, the corresponding chromosomal loci were not rearranged, which strongly suggests that a postreplicative event was responsible for the formation of each of these initial amplicons.

An intragenic deletion/inversion event in the DMD gene determines a novel exon creation and results in a BMD phenotype

Duchenne and Becker Muscular Dystrophy (DMD and BMD) are caused, in the majority of cases, by deletions in the dystrophin gene ( DMD). Here we describe the unprecedented case of a BMD patient carrying a large out-of-frame intragenic deletion, together with an inversion in the DMD gene, resulting in the inclusion of a novel exon in the transcript. Multiplex PCR amplification revealed the presence of a 48-52 exon deletion, but transcript analysis identified two unexpected products, neither of them including exon 53. The shorter mRNA derived from the juxtaposition of exons 47-54 (in-frame), while the longer one resulted from the inclusion of a novel 73-bp exon between exons 47 and 54. Sequence analysis revealed that the inserted sequence derived from an inverted portion of intron 53; its inclusion is predicted to determine protein truncation. The presence of a genomic inversion involving exon 53 and flanking regions was confirmed, and inversion/deletion breakpoints were sequenced. The inverted 73-bp sequence displays splicing signals at both ends and thus it is probably recognized as a novel exon when the partially inverted hnRNA is processed. These findings highlight the importance of mRNA analysis on patients that, based on routine DNA screenings, do not follow the reading-frame rule. This is the first reported patient carrying both an intragenic deletion and inversion in the DMD locus. This case might provide further insight into both the mechanisms that determine genomic rearrangements in the DMD locus and the molecular signals that drive exon inclusion.

Dystrophin and mutations: one gene, several proteins, multiple phenotypes

A large and complex gene on the X chromosome encodes dystrophin. Many mutations have been described in this gene, most of which affect the expression of the muscle isoform, the best-known protein product of this locus. These mutations result in the Duchenne and Becker muscular dystrophies (DMD and BMD). However, there are several other tissue specific isoforms of dystrophin, some exclusively or predominantly expressed in the brain or the retina. Mutations affecting the correct expression of these tissue-specific isoforms have been associated with the CNS involvement common in DMD. Rare mutations also account for the allelic disorder X-linked dilated cardiomyopathy, in which dystrophin expression or function is affected mostly or exclusively in the heart. Genotype definition of the dystrophin gene in patients with dystrophinopathies has taught us much about functionally important domains of the protein itself and has provided insights into several regulatory mechanisms governing the gene expression profile. Here, we focus on current understanding of the genotype-phenotype relation for mutations in the dystrophin gene and their implications for gene functions.

Comparative analysis of vertebrate dystrophin loci indicate intron gigantism as a common feature

The human DMD gene is the largest known to date, spanning > 2000 kb on the X chromosome. The gene size is mainly accounted for by huge intronic regions. We sequenced 190 kb of Fugu rubripes (pufferfish) genomic DNA corresponding to the complete dystrophin gene (FrDMD) and provide the first report of gene structure and sequence comparison among dystrophin genomic sequences from different vertebrate organisms. Almost all intron positions and phases are conserved between FrDMD and its mammalian counterparts, and the predicted protein product of the Fugu gene displays 55% identity and 71% similarity to human dystrophin. In analogy to the human gene, FrDMD presents several-fold longer than average intronic regions. Analysis of intron sequences of the human and murine genes revealed that they are extremely conserved in size and that a similar fraction of total intron length is represented by repetitive elements; moreover, our data indicate that intron expansion through repeat accumulation in the two orthologs is the result of independent insertional events. The hypothesis that intron length might be functionally relevant to the DMD gene regulation is proposed and substantiated by the finding that dystrophin intron gigantism is common to the three vertebrate genes.