Regional genomic instability predisposes to complex dystrophin gene rearrangements

Mutations in the dystrophin gene (DMD) cause Duchenne and Becker muscular dystrophies and the majority of cases are due to DMD gene rearrangements. Despite the high incidence of these aberrations, little is known about their causative molecular mechanism(s). We examined 792 DMD/BMD clinical samples by oligonucleotide array-CGH and report on the junction sequence analysis of 15 unique deletion cases and three complex intragenic rearrangements to elucidate potential underlying mechanism(s). Furthermore, we present three cases with intergenic rearrangements involving DMD and neighboring loci. The cases with intragenic rearrangements include an inversion with flanking deleted sequences; a duplicated segment inserted in direct orientation into a deleted region; and a splicing mutation adjacent to a deletion. Bioinformatic analysis demonstrated that 7 of 12 breakpoints combined among 3 complex cases aligned with repetitive sequences, as compared to 4 of 30 breakpoints for the 15 deletion cases. Moreover, the inversion/deletion case may involve a stem-loop structure that has contributed to the initiation of this rearrangement. For the duplication/deletion and splicing mutation/deletion cases, the presence of the first mutation, either a duplication or point mutation, may have elicited the deletion events in an attempt to correct preexisting mutations. While NHEJ is one potential mechanism for these complex rearrangements, the highly complex junction sequence of the inversion/deletion case suggests the involvement of a replication-based mechanism. Our results support the notion that regional genomic instability, aided by the presence of repetitive elements, a stem-loop structure, and possibly preexisting mutations, may elicit complex rearrangements of the DMD gene.     

Molecular cytogenetic characterization of a constitutional complex intrachromosomal 4q rearrangement in a patient with multiple congenital anomalies

Constitutional Complex Chromosomal Rearrangements (CCRs) are very rare. While the vast majority of CCRs involve more than one chromosome, only seven cases describe CCRs with four or more breakpoints within a single chromosome. Here, we present a patient with multiple congenital anomalies and mental retardation. Array Comparative Genomic Hybridisation (array CGH), FISH and Multicolour Banding FISH revealed a de novo complex rearrangement with two deletions, a duplication and an inversion of 4q. This CCR involving at least seven breakpoints is one of the most complex rearrangements of a single chromosome reported thus far. Potential mechanisms generating such complex rearrangements are discussed.     

Rearrangements at a hobo element inserted into the first intron of the singed gene in the unstable sn 49 system of Drosophila melanogaster

The cytological structure of the X chromosome and the DNA organisation of the singed locus were examined in five singed bristle mutants of Drosophila melanogaster. These mutants are all derived from the unstable mutant singed-49, isolated from a wild population in the Russian Far East in 1975. Rearrangements were found at a site within the first intron of the singed gene, where a hobo element is inserted in these mutants. One rearrangement, which is associated with a strong bristle phenotype, has an inversion between 2D and the location of singed at 7D, which separates the singed promoter from the singed coding region. Two phenotypically wild-type derivatives have smaller rearrangements within the first intron which do not appear to interfere with singed expression. Two derivatives with bristle phenotypes have more complex rearrangements, and one of them shows a dominant or antimorphic phenotype. DNA blotting and in situ hybridisation experiments show that, in addition to these rearrangements at a hobo element inserted at singed, other hobo elements in these strains have been mobilised. This system is therefore similar to others in which functional hobo elements continue to transpose, resulting in elevated rates of mutation and chromosome rearrangement.     

Reassessment of the evolution of wheat chromosomes 4A, 5A,and 7B

Key message

Comparison of genome sequences of wild emmer wheat and Aegilops tauschii suggests a novel scenario of the evolution of rearranged wheat chromosomes 4A, 5A, and 7B.

Abstract

Past research suggested that wheat chromosome 4A was subjected to a reciprocal translocation T(4AL;5AL)1 that occurred in the diploid progenitor of the wheat A subgenome and to three major rearrangements that occurred in polyploid wheat: pericentric inversion Inv(4AS;4AL)1, paracentric inversion Inv(4AL;4AL)1, and reciprocal translocation T(4AL;7BS)1. Gene collinearity along the pseudomolecules of tetraploid wild emmer wheat (Triticum turgidum ssp. dicoccoides, subgenomes AABB) and diploid Aegilops tauschii (genomes DD) was employed to confirm these rearrangements and to analyze the breakpoints. The exchange of distal regions of chromosome arms 4AS and 4AL due to pericentric inversion Inv(4AS;4AL)1 was detected, and breakpoints were validated with an optical Bionano genome map. Both breakpoints contained satellite DNA. The breakpoints of reciprocal translocation T(4AL;7BS)1 were also found. However, the breakpoints that generated paracentric inversion Inv(4AL;4AL)1 appeared to be collocated with the 4AL breakpoints that had produced Inv(4AS;4AL)1 and T(4AL;7BS)1. Inv(4AS;4AL)1, Inv(4AL;4AL)1, and T(4AL;7BS)1 either originated sequentially, and Inv(4AL;4AL)1 was produced by recurrent chromosome breaks at the same breakpoints that generated Inv(4AS;4AL)1 and T(4AL;7BS)1, or Inv(4AS;4AL)1, Inv(4AL;4AL)1, and T(4AL;7BS)1 originated simultaneously. We prefer the latter hypothesis since it makes fewer assumptions about the sequence of events that produced these chromosome rearrangements.

     

Molecular cytogenetic characterization of the first reported case of an inv dup (4p)(p15.1-pter) with a concomitant 4q35.1-qter deletion and normal parents

Inverted duplications associated with terminal deletions are complex anomalies described in an increasing of chromosome ends. We report on the cytogenetic characterization of the first de novo inv dup del(4) with partial 4p duplication and 4q deletion in a girl with clinical signs consistent with “recombinant 4 syndrome”. This abnormality was suspected by banding, but high-resolution molecular cytogenetic investigations allowed us to define the breakpoints of the rearrangement. The terminal duplicated region extending from 4p15.1 to the telomere was estimated to be 29.27 Mb, while the size of the terminal deletion was 3.114 Mb in the 4q35.1 region. Until now, 10 patients with duplicated 4p14-p15 and deleted 4q35 chromosome 4 have been described. In all cases the abnormal chromosome 4 was derived from a pericentric inversion inherited from one of the parents. In conclusion, we have identified the first case of inv dup del(4) with normal parents suggesting that, often, terminal duplications or terminal deletions mask complex rearrangements.     

Breakpoint junctions of chromosome 9p deletions in two human glioma cell lines.

Interstitial deletions of the short arm of chromosome 9 are associated with glioma, acute lymphoblastic leukemia, melanoma, mesothelioma, lung cancer, and bladder cancer. The distal breakpoints of the deletions (in relation to the centromere) in 14 glioma and leukemia cell lines have been mapped within the 400 kb IFN gene cluster located at band 9p21. To obtain information about the mechanism of these deletions, we have isolated and analyzed the nucleotide sequences at the breakpoint junctions in two glioma-derived cell lines. The A1235 cell line has a complex rearrangement of chromosome 9, including a deletion and an inversion that results in two breakpoint junctions. Both breakpoints of the distal inversion junction occurred within AT-rich regions. In the A172 cell line, a tandem heptamer repeat was found on either side of the deletion breakpoint junction. The distal breakpoint occurred 5' of IFNA2; the 256 bp sequenced from the proximal side of the breakpoint revealed 95% homology to long interspersed nuclear elements. One- and two-base-pair overlaps were observed at these junctions. The possible role of sequence overlaps, and repetitive sequences, in the rearrangement is discussed.     

Decoding NF1 Intragenic Copy-Number Variations

Genomic rearrangements can cause both Mendelian and complex disorders. Currently, several major mechanisms causing genomic rearrangements, such as non-allelic homologous recombination (NAHR), non-homologous end joining (NHEJ), fork stalling and template switching (FoSTeS), and microhomology-mediated break-induced replication (MMBIR), have been proposed. However, to what extent these mechanisms contribute to gene-specific pathogenic copy-number variations (CNVs) remains understudied. Furthermore, few studies have resolved these pathogenic alterations at the nucleotide-level. Accordingly, our aim was to explore which mechanisms contribute to a large, unique set of locus-specific non-recurrent genomic rearrangements causing the genetic neurocutaneous disorder neurofibromatosis type 1 (NF1). Through breakpoint-spanning PCR as well as array comparative genomic hybridization, we have identified the breakpoints in 85 unrelated individuals carrying an NF1 intragenic CNV. Furthermore, we characterized the likely rearrangement mechanisms of these 85 CNVs, along with those of two additional previously published NF1 intragenic CNVs. Unlike the most typical recurrent rearrangements mediated by flanking low-copy repeats (LCRs), NF1 intragenic rearrangements vary in size, location, and rearrangement mechanisms. We propose the DNA-replication-based mechanisms comprising both FoSTeS and/or MMBIR and serial replication stalling to be the predominant mechanisms leading to NF1 intragenic CNVs. In addition to the loop within a 197-bp palindrome located in intron 40, four Alu elements located in introns 1, 2, 3, and 50 were also identified as intragenic-rearrangement hotspots within NF1.     

RET rearrangements in radiation-induced papillary thyroid carcinomas: high prevalence of topoisomerase I sites at breakpoints and microhomology-mediated end joining in ELE1 and RET chimeric genes

Children exposed to radioactive iodine after the Chernobyl reactor accident frequently developed papillary thyroid carcinomas (PTC). The predominant molecular lesions in these tumors are rearrangements of the RET receptor tyrosine kinase gene. Various types of RET rearrangements have been described. More than 90% of PTC with RET rearrangement exhibit a PTC1 or PTC3 type of rearrangement with an inversion of the H4 or ELE1 gene, respectively, on chromosome 10. To obtain closer insight into the mechanisms underlying PTC3 inversions, we analyzed the genomic breakpoints of 22 reciprocal and 4 nonreciprocal ELE1 and RET rearrangements in 26 post-Chernobyl tumor samples. In contrast to previous assumptions, an accumulation of breakpoints at the two Alu elements in the ELE1 sequence was not observed. Instead, breakpoints are distributed in the affected introns of both genes without significant clustering. When compared to the corresponding wildtype sequences, the majority of breakpoints (92%) do not contain larger deletions or insertions. Most remarkably, at least one topoisomerase I site was found exactly at or in close vicinity to all breakpoints, indicating a potential role for this enzyme in the formation of DNA strand breaks and/or ELE1 and RET inversions. The presence of short regions of sequence homology (microhomologies) and short direct and inverted repeats at the majority of breakpoints furthermore indicates a nonhomologous DNA end-joining mechanism in the formation of chimeric ELE1/Ret and Ret/ELE1 genes.     

Array CGH detection of a cryptic deletion in a complex chromosome rearrangement

Balanced complex chromosome rearrangements (CCR) are extremely rare in humans. They are usually ascertained either by abnormal phenotype or reproductive failure in carriers. These abnormalities are attributed to disruption of genes at the breakpoints, position effect or cryptic imbalances in the genome. However, little is known about possible imbalances at the junction points. We report here a patient with a CCR involving three chromosomes (2;10;11) and eight breakpoints. The patient presented with behavioural problems as the sole phenotypic abnormality. The rearrangement, which is apparently balanced in G-banding and multicolour FISH, was shown by genomic array analysis to include a deletion of 0.15–1.5 Mb associated with one of the breakpoints. To explain the formation of this rearrangement through the smallest possible number of breakage-and-reunion events, one has to assume that the breaks have not occurred simultaneously, but in a temporal order within the span of a single cell division. We demonstrate that array comparative genomic hybridisation (CGH) is a useful complementary tool to cytogenetic analysis for detecting and mapping cryptic imbalances associated with chromosome rearrangement.     

Comparative genomics, marker density and statistical analysis of chromosome rearrangements

Estimates of the number of chromosomal breakpoints that have arisen (e.g., by translocation and inversion) in the evolutionary past between two species and their common ancestor can be made by comparing map positions of marker loci. Statistical methods for doing so are based on a random-breakage model of chromosomal rearrangement. The model treats all modes of chromosome rearrangement alike, and it assumes that chromosome boundaries and breakpoints are distributed randomly along a single genomic interval. Here we use simulation and numerical analysis to test the validity of these model assumptions. Mean estimates of numbers of breakpoints are close to those expected under the random-breakage model when marker density is high relative to the amount of chromosomal rearrangement and when rearrangements occur by translocation alone. But when marker density is low relative to the number of chromosomes, and when rearrangements occur by both translocation and inversion, the number of breakpoints is underestimated. The underestimate arises because rearranged segments may contain markers, yet the rearranged segments may, nevertheless, be undetected. Variances of the estimate of numbers of breakpoints decrease rapidly as markers are added to the comparative maps, but are less influenced by the number or type of chromosomal rearrangement separating the species. Variances obtained with simulated genomes comprised of chromosomes of equal length are substantially lower than those obtained when chromosome size is unconstrained. Statistical power for detecting heterogeneity in the rate of chromosomal rearrangement is also investigated. Results are interpreted with respect to the amount of marker information required to make accurate inferences about chromosomal evolution.     

Distribution of induced chromosome rearrangement breakpoints along the chromosome length and the problem of intercalary heterochromatin in Drosophila melanogaster

Historically, the term “intercalary heterochromatin” was based on the finding that induced chromosome rearrangements occur at a higher frequency in the corresponding regions. The available molecular genetic data and, in particular, the results of the Drosophila Genome Project made it possible to decide between two possible explanations of the preferential location of chromosome rearrangement breakpoints in intercalary heterochromatin regions. Namely, a higher frequency of radiation-induced rearrangements in these regions correlates with the DNA content and probably lacks an association with the features of chromatin organization.     

Locating rearrangement events in a phylogeny based on highly fragmented assemblies

Background

The inference of genome rearrangement operations requires complete genome assemblies as input data, since a rearrangement can involve an arbitrarily large proportion of one or more chromosomes. Most genome sequence projects, especially those on non-model organisms for which no physical map exists, produce very fragmented assembles, so that a rearranged fragment may be impossible to identify because its two endpoints are on different scaffolds. However, breakpoints are easily identified, as long as they do not coincide with scaffold ends. For the phylogenetic context, in comparing a fragmented assembly with a number of complete assemblies, certain combinatorial constraints on breakpoints can be derived. We ask to what extent we can use breakpoint data between a fragmented genome and a number of complete genomes to recover all the arrangements in a phylogeny.

Results

We simulate genomic evolution via chromosomal inversion, fragmenting one of the genomes into a large number of scaffolds to represent the incompleteness of assembly. We identify all the breakpoints between this genome and the remainder. We devise an algorithm which takes these breakpoints into account in trying to determine on which branch of the phylogeny a rearrangement event occurred. We present an analysis of the dependence of recovery rates on scaffold size and rearrangement rate, and show that the true tree, the one on which the rearrangement simulation was performed, tends to be most parsimonious in estimating the number of true events inferred.

Conclusions

It is somewhat surprising that the breakpoints identified just between the fragmented genome and each of the others suffice to recover most of the rearrangements produced by the simulations. This holds even in parts of the phylogeny disjoint from the lineage of the fragmented genome.

     

Determination of the Genome Sequence of Porphyromonas gingivalis Strain ATCC 33277 and Genomic Comparison with Strain W83 Revealed Extensive Genome Rearrangements in P. gingivalis

The gram-negative anaerobic bacterium Porphyromonas gingivalis is a major causative agent of chronic periodontitis. Porphyromonas gingivalis strains have been classified into virulent and less-virulent strains by mouse subcutaneous soft tissue abscess model analysis. Here, we present the whole genome sequence of P. gingivalis ATCC 33277, which is classified as a less-virulent strain. We identified 2090 protein-coding sequences (CDSs), 4 RNA operons, and 53 tRNA genes in the ATCC 33277 genome. By genomic comparison with the virulent strain W83, we identified 461 ATCC 33277-specific and 415 W83-specific CDSs. Extensive genomic rearrangements were observed between the two strains: 175 regions in which genomic rearrangements have occurred were identified. Thirty-five of those genomic rearrangements were inversion or translocation and 140 were simple insertion, deletion, or replacement. Both strains contained large numbers of mobile elements, such as insertion sequences, miniature inverted-repeat transposable elements (MITEs), and conjugative transposons, which are frequently associated with genomic rearrangements. These findings indicate that the mobile genetic elements have been deeply involved in the extensive genome rearrangement of P. gingivalis and the occurrence of many of the strain-specific CDSs. We also describe here a very unique feature of MITE400, which we renamed MITEPgRS (MITE of P. gingivalis with Repeating Sequences).Key words: Porphyromonas gingivalis, whole genome sequence, genome rearrangement, conjugative transposon, MITE     

Chromosomal Inversions between Human and Chimpanzee Lineages Caused by Retrotransposons

The long interspersed element-1 (LINE-1 or L1) and Alu elements are the most abundant mobile elements comprising 21% and 11% of the human genome, respectively. Since the divergence of human and chimpanzee lineages, these elements have vigorously created chromosomal rearrangements causing genomic difference between humans and chimpanzees by either increasing or decreasing the size of genome. Here, we report an exotic mechanism, retrotransposon recombination-mediated inversion (RRMI), that usually does not alter the amount of genomic material present. Through the comparison of the human and chimpanzee draft genome sequences, we identified 252 inversions whose respective inversion junctions can clearly be characterized. Our results suggest that L1 and Alu elements cause chromosomal inversions by either forming a secondary structure or providing a fragile site for double-strand breaks. The detailed analysis of the inversion breakpoints showed that L1 and Alu elements are responsible for at least 44% of the 252 inversion loci between human and chimpanzee lineages, including 49 RRMI loci. Among them, three RRMI loci inverted exonic regions in known genes, which implicates this mechanism in generating the genomic and phenotypic differences between human and chimpanzee lineages. This study is the first comprehensive analysis of mobile element bases inversion breakpoints between human and chimpanzee lineages, and highlights their role in primate genome evolution.     

Construction and analysis of a two-chromosome double reciprocal translocation that functions as a ‘balancer’ inNeurospora crassa

T(IIL; VL;IIR; VR) BLNC-1 is a compound chromosome rearrangement inNeurospora crassa that combines two reciprocal translocations:T(IIL; VL) AR30 which interchanges the left end of linkage group II with the left end of linkage group V, andT(IIR;VR) ALS154 which interchanges the right end of linkage group II with the right end of linkage group V.BLNC-1 acts as a crossover suppressor for most of both linkage groups II and V since single crossovers between the rearrangement breakpoints result in progeny with lethal unbalanced duplications and deficiencies. The integrity ofBLNC-1 following meiosis was tested in crosses of markedBLNC-1 by marked Normal sequence, with markers located at critical points on linkage groups II and V. Although recombination between distal markers in the four arms was reduced markedly, double crossovers in the long intervening regions occurred with a frequency of 21%. Of these double crossovers, most were coincidental crossovers, one in each of the long intervening regions, resulting in the resolution of the complex into its component rearrangements (16%), while a minority of the double crossovers (5%) were crossovers involving only one of the two component linkage groups, and resulted in the insertion of a segment between the breakpoints. - TheBLNC-1 balancer can be used for: (1) mapping new loci to linkage groups II and V, especially for identifying markers mapping near the tips of the linkage groups; (2) for isolating genetically intact chromosomes from natural populations or for quantitative genetic studies; and (3) for studying recombinational hot-spots which can be detected as escapes from crossover suppression. -Based on experience withBLNC-1, future two-chromosome balancers should be designed with two breakpoints near, but not at, the opposite ends of the chromosome to be balanced, and the other two breakpoints close to, but spanning, the centromere of a second chromosome. Such a construction when combined with appropriately placed selective markers should prevent breakdown of the complex, and should resemble an inversion in eliminating crossover products. Contribution no. 85-218-J from the Department of Plant Pathology, Kansas Agricultural Experiment Station, Kansas State University, Manhattan.     

Molecular analysis of a constitutional complex genome rearrangement with 11 breakpoints involving chromosomes 3, 11, 12, and 21 and a ∼0.5-Mb submicroscopic deletion in a patient with mild mental retardation

Complex chromosome rearrangements (CCRs) are extremely rare but often associated with mental retardation, congenital anomalies, or recurrent spontaneous abortions. We report a de novo apparently balanced CCR involving chromosomes 3 and 12 and a two-way translocation between chromosomes 11 and 21 in a woman with mild intellectual disability, obesity, coarse facies, and apparent synophrys without other distinctive dysmorphia or congenital anomalies. Molecular analysis of breakpoints using fluorescence in situ hybridization (FISH) with region-specific BAC clones revealed a more complex character for the CCR. The rearrangement is a result of nine breaks and involves reciprocal translocation of terminal chromosome fragments 3p24.1→pter and 12q23.1→qter, insertion of four fragments of the long arm of chromosome 12: q14.1→q21?, q21?→q22, q22→q23.1, and q23.1→q23.1 and a region 3p22.3→p24.1 into chromosome 3q26.31. In addition, we detected a ～0.5-Mb submicroscopic deletion at 3q26.31. The deletion involves the chromosome region that has been previously associated with Cornelia de Lange syndrome (CdLS) in which a novel gene NAALADL2 has been mapped recently. Other potential genes responsible for intellectual deficiency disrupted as a result of patient’s chromosomal rearrangement map at 12q14.1 (TAFA2), 12q23.1 (METAP2), and 11p14.1 (BDNF).     

Rearrangements at a hobo element inserted into the first intron of the singed gene in the unstable sn 49 system of Drosophila melanogaster

The cytological structure of the X chromosome and the DNA organisation of the singed locus were examined in five singed bristle mutants of Drosophila melanogaster. These mutants are all derived from the unstable mutant singed-49, isolated from a wild population in the Russian Far East in 1975. Rearrangements were found at a site within the first intron of the singed gene, where a hobo element is inserted in these mutants. One rearrangement, which is associated with a strong bristle phenotype, has an inversion between 2D and the location of singed at 7D, which separates the singed promoter from the singed coding region. Two phenotypically wild-type derivatives have smaller rearrangements within the first intron which do not appear to interfere with singed expression. Two derivatives with bristle phenotypes have more complex rearrangements, and one of them shows a dominant or antimorphic phenotype. DNA blotting and in situ hybridisation experiments show that, in addition to these rearrangements at a hobo element inserted at singed, other hobo elements in these strains have been mobilised. This system is therefore similar to others in which functional hobo elements continue to transpose, resulting in elevated rates of mutation and chromosome rearrangement. Received: 19 February 1997 / Accepted: 8 October 1997     

Sequence Features Contributing to Chromosomal Rearrangements in Neisseria gonorrhoeae

Through whole genome sequence alignments, breakpoints in chromosomal synteny can be identified and the sequence features associated with these determined. Alignments of the genome sequences of Neisseria gonorrhoeae strain FA1090, N. gonorrhoeae strain NCCP11945, and N. gonorrhoeae strain TCDC-NG08107 reveal chromosomal rearrangements that have occurred. Based on these alignments and dot plot pair-wise comparisons, the overall chromosomal arrangement of strain NCCP11945 and TCDC-NG08107 are very similar, with no large inversions or translocations. The insertion of the Gonococcal Genetic Island in strain NCCP11945 is the most prominent distinguishing feature differentiating these strains. When strain NCCP11945 is compared to strain FA1090, however, 14 breakpoints in chromosomal synteny are identified between these gonococcal strains. The majority of these, 11 of 14, are associated with a prophage, IS elements, or IS-like repeat enclosed elements which appear to have played a role in the rearrangements observed. Additional rearrangements of small regions of the genome are associated with pilin genes. Evidence presented here suggests that the rearrangements of blocks of sequence are mediated by activation of prophage and associated IS elements and reintegration elsewhere in the genome or by homologous recombination between IS-like elements that have generated inversions.