Are There Rearrangement Hotspots in the Human Genome?

In a landmark paper, Nadeau and Taylor [18] formulated the random breakage model (RBM) of chromosome evolution that postulates that there are no rearrangement hotspots in the human genome. In the next two decades, numerous studies with progressively increasing levels of resolution made RBM the de facto theory of chromosome evolution. Despite the fact that RBM had prophetic prediction power, it was recently refuted by Pevzner and Tesler [4], who introduced the fragile breakage model (FBM), postulating that the human genome is a mosaic of solid regions (with low propensity for rearrangements) and fragile regions (rearrangement hotspots). However, the rebuttal of RBM caused a controversy and led to a split among researchers studying genome evolution. In particular, it remains unclear whether some complex rearrangements (e.g., transpositions) can create an appearance of rearrangement hotspots. We contribute to the ongoing debate by analyzing multi-break rearrangements that break a genome into multiple fragments and further glue them together in a new order. In particular, we demonstrate that (1) even if transpositions were a dominant force in mammalian evolution, the arguments in favor of FBM still stand, and (2) the "gene deletion" argument against FBM is flawed.     

Estimation of the true evolutionary distance under the fragile breakage model

Background

The ability to estimate the evolutionary distance between extant genomes plays a crucial role in many phylogenomic studies. Often such estimation is based on the parsimony assumption, implying that the distance between two genomes can be estimated as the rearrangement distance equal the minimal number of genome rearrangements required to transform one genome into the other. However, in reality the parsimony assumption may not always hold, emphasizing the need for estimation that does not rely on the rearrangement distance. The distance that accounts for the actual (rather than minimal) number of rearrangements between two genomes is often referred to as the true evolutionary distance. While there exists a method for the true evolutionary distance estimation, it however assumes that genomes can be broken by rearrangements equally likely at any position in the course of evolution. This assumption, known as the random breakage model, has recently been refuted in favor of the more rigorous fragile breakage model postulating that only certain “fragile” genomic regions are prone to rearrangements.

Results

We propose a new method for estimating the true evolutionary distance between two genomes under the fragile breakage model. We evaluate the proposed method on simulated genomes, which show its high accuracy. We further apply the proposed method for estimation of evolutionary distances within a set of five yeast genomes and a set of two fish genomes.

Conclusions

The true evolutionary distances between the five yeast genomes estimated with the proposed method reveals that some pairs of yeast genomes violate the parsimony assumption. The proposed method further demonstrates that the rearrangement distance between the two fish genomes underestimates their evolutionary distance by about 20%. These results demonstrate how drastically the two distances can differ and justify the use of true evolutionary distance in phylogenomic studies.

     

Chromosome breakage and recombination at fragile sites.

Chromosomal fragile sites are points on chromosomes that usually appear as nonstaining chromosome or chromatid gaps. It has frequently been suggested that fragile sites may be involved in chromosome breakage and recombination events. We and others have previously shown that fragile sites predispose to intrachromosomal recombination as measured by sister-chromatid exchanges. These findings suggested that fragile site expression often, if not always, is accompanied by DNA strand breakage. In the present report, fragile sites are shown to predispose to deletions and interchromosomal recombination. By use of somatic cell hybrids containing either human chromosome 3 or the fragile X chromosome, deletions and translocations were induced by FUdR or aphidicolin with breakpoints at the fragile sites Xq27 or 3p14.2 (FRA3B) or at points so close to the fragile sites as to be cytogenetically indistinguishable. Southern blot analysis of DNA from a panel of chromosome 3 deletion and translocation hybrids was then utilized to detect loss or retention of markers flanking FRA3B and to corroborate the cytogenetic evidence that the breakpoints were at this fragile site. One cell line with a reciprocal translocation between human chromosome 3 (with breakpoint at 3p14.2) and a hamster chromosome showed cytogenetically that the fragile site was expressed on both derivative chromosomes, supporting the hypothesis that the fragile site represents a repeated sequence. The approach described provides a means of generating specific rearrangements in somatic cell hybrids with a breakpoint at a fragile site.(ABSTRACT TRUNCATED AT 250 WORDS)     

Bands involved in primary chromosome rearrangements in sarcomas are not constitutionally liable to breakage in sarcoma patients

Summary The localization of breakpoints in spontaneous chromosome aberrations, i.e., chromatid and chromosome gaps, breaks, and exchanges, has been studied in cultured skin fibroblasts from 34 untreated patients with musculoskeletal sarcoma and 38 controls. A total of 325 aberrations in the sarcoma group and 251 in the control group could be assigned to particular bands. The distribution was non-random (P<0.001) in both groups. Twenty-one bands in the sarcoma group and 20 in the control group appeared as hot spots, with 11 represented in both groups. Only three hot spots, all of which were present among both patients and controls, coincided with bands involved in primary sarcoma-associated chromosome rearrangements. The results indicate that the chromosome breakage pattern of non-malignant cells is similar in sarcoma patients and controls. Hence, the occurrence of primary structural rearrangements in sarcomas cannot be accounted for by any constitutional proneness to chromosome breakage at these bands.     

129 The molecular mechanism of breakage at fragile site FRA16D

Replication stress induces physical breakage at discrete loci in chromosomes, which can be visualized on a metaphase chromosome spread. These common fragile sites (CFS) are conserved across species and are hotspots for sister chromatid recombination, viral integration, rearrangements, translocations, and deletions (Glover et al 2005). Despite multiple theories, the molecular mechanisms of CFS expression and genomic instability are still not well understood. The fragile site FRA16D is of special interest because it is the second most highly expressed fragile site and is located within the WWOX tumor suppressor gene. Previous data identified a polymorphic AT repeat within a FRA16D subregion called F1 that causes chromosome fragility and replication fork stalling in a yeast model (Zhang and Freudenreich 2007). Recently, we have found that breakage increases in an AT repeat length-dependent manner. Our results suggest that the AT repeat in the context of F1 forms a secondary structure, making the region more vulnerable to breakage.     

Chromosomal breakpoint reuse in genome sequence rearrangement.

In order to apply gene-order rearrangement algorithms to the comparison of genome sequences, Pevzner and Tesler bypass gene finding and ortholog identification and use the order of homologous blocks of unannotated sequence as input. The method excludes blocks shorter than a threshold length. Here we investigate possible biases introduced by eliminating short blocks, focusing on the notion of breakpoint reuse introduced by these authors. Analytic and simulation methods show that reuse is very sensitive to the proportion of blocks excluded. As is pertinent to the comparison of mammalian genomes, this exclusion risks randomizing the comparison partially or entirely.     

The random versus fragile breakage models of chromosome evolution: a matter of resolution

Conserved synteny––the sharing of at least one orthologous gene by a pair of chromosomes from two species––can, in the strictest sense, be viewed as sequence conservation between chromosomes of two related species, irrespective of whether coding or non-coding sequence is examined. The recent sequencing of multiple vertebrate genomes indicates that certain chromosomal segments of considerable size are conserved in gene order as well as underlying non-coding sequence across all vertebrates. Some of these segments lost genes or non-coding sequence and/or underwent breakage only in teleost genomes, presumably because evolutionary pressure acting on these regions to remain intact were relaxed after an additional round of whole genome duplication. Random reporter insertions into zebrafish chromosomes combined with computational genome-wide analysis indicate that large chromosomal areas of multiple genes contain long-range regulatory elements, which act on their target genes from several gene distances away. In addition, computational breakpoint analyses suggest that recurrent evolutionary breaks are found in “fragile regions” or “hotspots”, outside of the conserved blocks of synteny. These findings cannot be accommodated by the random breakage model and suggest that this view of genome and chromosomal evolution requires substantial reassessment.     

Congenital chromosome breakage clusters within Giemsa-light bands and identifies sites of chromatin instability

Analysis of the distribution of published chromosome breaks in cells with constitutional chromosome aberrations showed a nonrandom distribution of breaks among chromosomes and chromosome regions. A significant amount of breakage occurred at Giemsa-negative bands. In addition, chromosome sites associated with a number of fragile sites and cellular oncogene sites were affected nonrandomly. The data are consistent with the hypothesis that chromosome breakage occurs in somatic or germ cells as a result of recombinational errors involving actively transcribing chromatin regions or regions of unstable DNA sequence structure placed in proximity during interphase.     

Chromosome breakage at a major fragile site associated with P-glycoprotein gene amplification in multidrug-resistant CHO cells.

Recent studies of several drug-resistant Chinese hamster cell lines suggested that a breakage-fusion-bridge mechanism is frequently involved in the amplification of drug resistance genes. These observations underscore the importance of chromosome breakage in the initiation of DNA amplification in mammalian cells. However, the mechanism of this breakage is unknown. Here, we propose that the site of chromosome breakage consistent with the initial event of P-glycoprotein (P-gp) gene amplification via the breakage-fusion-bridge cycle in three independently established multidrug-resistant CHO cells was located at 1q31. This site is a major chromosome fragile site that can be induced by methotrexate and aphidicolin treatments. Pretreatments of CHO cells with methotrexate or aphidicolin enhanced the frequencies of resistance to vinca alkaloid and amplification of the P-gp gene. These observations suggest that chromosome fragile sites play a pivotal role in DNA amplification in mammalian cells. Our data are also consistent with the hypothesis that gene amplification can be initiated by stress-induced chromosome breakage that is independent of modes of action of cytotoxic agents. Drug-resistant variants may arise by their growth advantage due to overproduction of cellular target molecules via gene amplification.     

Folic acid and chromosome breakage. II. A methionine effect similar to that in fragile X expression

Summary The absence of folic acid results in an increase in spontaneous chromosome breakage in the presence of, but not in the absence of, methionine. This is similar to the methionine effect seen in studies of the fragile X chromosome. The results support the idea that fragile sites and at least some types of spontaneous chromosome breaks may share a common mechanism.     

Inhibition of topoisomerase I prevents chromosome breakage at common fragile sites

Common fragile sites are loci that preferentially form gaps and breaks on metaphase chromosomes when DNA synthesis is perturbed, particularly after treatment with the DNA polymerase inhibitor, aphidicolin. We and others have identified several cell cycle checkpoint and DNA repair proteins that influence common fragile site stability. However, the initial events underlying fragile site breakage remain poorly understood. We demonstrate here that aphidicolin-induced gaps and breaks at fragile sites are prevented when cells are co-treated with low concentrations of the topoisomerase I inhibitor, camptothecin. This reduction in breakage is accompanied by a reduction in aphidicolin-induced RPA foci, CHK1 and RPA2 phosphorylation, and PCNA monoubiquitination, indicative of reduced levels of single stranded DNA. Furthermore, camptothecin reduces spontaneous fragile site breakage seen in cells lacking ATR, even in the absence of aphidicolin. These data from cultured human cells demonstrate that topoisomerase I activity is required for DNA common fragile site breaks and suggest that polymerase–helicase uncoupling is a key initial event in this process.     

SCJ: a breakpoint-like distance that simplifies several rearrangement problems

The breakpoint distance is one of the most straightforward genome comparison measures. Surprisingly, when it comes to defining it precisely for multichromosomal genomes with both linear and circular chromosomes, there is more than one way to go about it. Pevzner and Tesler gave a definition in a 2003 paper, Tannier et al. defined it differently in 2008, and in this paper we provide yet another alternative, calling it SCJ for single-cut-or-join, in analogy to the popular double cut and join (DCJ) measure. We show that several genome rearrangement problems, such as median and halving, become easy for SCJ, and provide linear and higher polynomial time algorithms for them. For the multichromosomal linear genome median problem, this is the first polynomial time algorithm described, since for other distances this problem is NP-hard. In addition, we show that small parsimony under SCJ is also easy, and can be solved by a variant of Fitch's algorithm. In contrast, big parsimony is NP-hard under SCJ. This new distance measure may be of value as a speedily computable, first approximation to distances based on more realistic rearrangement models.     

Cloning and characterization of chromosome breakpoints of Plasmodium falciparum: breakage and new telomere formation occurs frequently and randomly in subtelomeric genes.

We analysed the genetic stability of two subtelomeric genes of the human malaria parasite Plasmodium falciparum. A PCR based assay, using a telomere and a target-gene specific primer was used to detect potential chromosome rearrangements. We show that chromosome breakage and the formation of new telomeres occur frequently in the two genes coding for histidine rich proteins (HRP I and HRP II) in laboratory isolates, but remains undetectable in clinical parasite isolates. This finding suggests an essential role of these genes in vivo and that chromosome breakage is rather an accidental process than a programmed chromosome fragmentation. Cloning and sequencing of 8 chromosome breakpoints of the HRP II gene from one parasite isolate shows that the breakage occurs within a broad region in which new telomere formation appear to take place at random sites. Furthermore, this analysis revealed no obvious sequence similarities of sites of telomere addition. Finally, we show that an irregular pattern of heterogeneous telomere repeats is added at each broken end and that each healed chromosome contains a distinct pattern of repeats. We discuss a model for telomere formation in P. falciparum.     

Erratum to “Random models of Menzerath–Altmann law in genomes” (BioSystems 107(3) (2012) 167–173)

Here we improve the mathematical arguments of Baixeries et al (BioSystems 107(3) (2012) 167–173). The corrections do not alter the conclusion that the random breakage model yields an insufficient fit to the scaling of mean chromosome length as a function of chromosome number in real genomes.     

Susceptibility of heterochromatin to aphidicolin-induced chromosomal breakage

The distribution of aphidicolin-induced chromosomal lesions was analyzed to determine the relative breakage susceptibility of euchromatin and heterochromatin in the cactus mouse, Peromyscus eremicus. The observed breakage was tested against expected distributions corresponding to the karyotypic proportions of autosomal euchromatin, autosomal heterochromatin, X-chromosomal euchromatin, and X-chromosomal heterochromatin. The distribution of induced breakage was independent of sex but dependent on the individual. In all individuals, there was a highly significant (P0.0001) deficiency in the number of breaks observed as compared to expected in autosomal heterochromatin. Sparse observations in the X chromosome and the absence of breaks in the Y chromosome precluded valid statistical tests of the sex-chromosomal distribution of induced breakage. These data indicate that the autosomal heterochromatin of Peromyscus is resistant to aphidicolin-induced chromosomal breakage and argue against a simple relationship between late replication and a general mechanism for chromosomal fragility.     

Topoisomerase II– and Condensin-Dependent Breakage of MEC1ATR-Sensitive Fragile Sites Occurs Independently of Spindle Tension,Anaphase, or Cytokinesis

Fragile sites are loci of recurrent chromosome breakage in the genome. They are found in organisms ranging from bacteria to humans and are implicated in genome instability, evolution, and cancer. In budding yeast, inactivation of Mec1, a homolog of mammalian ATR, leads to chromosome breakage at fragile sites referred to as replication slow zones (RSZs). RSZs are proposed to be homologous to mammalian common fragile sites (CFSs) whose stability is regulated by ATR. Perturbation during S phase, leading to elevated levels of stalled replication forks, is necessary but not sufficient for chromosome breakage at RSZs or CFSs. To address the nature of additional event(s) required for the break formation, we examined involvement of the currently known or implicated mechanisms of endogenous chromosome breakage, including errors in replication fork restart, premature mitotic chromosome condensation, spindle tension, anaphase, and cytokinesis. Results revealed that chromosome breakage at RSZs is independent of the RAD52 epistasis group genes and of TOP3, SGS1, SRS2, MMS4, or MUS81, indicating that homologous recombination and other recombination-related processes associated with replication fork restart are unlikely to be involved. We also found spindle force, anaphase, or cytokinesis to be dispensable. RSZ breakage, however, required genes encoding condensin subunits (YCG1, YSC4) and topoisomerase II (TOP2). We propose that chromosome break formation at RSZs following Mec1 inactivation, a model for mammalian fragile site breakage, is mediated by internal chromosomal stress generated during mitotic chromosome condensation.     

Chromosome breakage in human preimplantation embryos from carriers of structural chromosomal abnormalities in relation to fragile sites, maternal age, and poor sperm factors

Chromosome breakage is a fairly widespread phenomenon in preimplantation embryos affecting at least 10% of day 3 cleavage stage embryos. It may be detected during preimplantation genetic diagnosis (PGD). For carriers of structural chromosomal abnormalities, PGD involves the removal and testing of single blastomeres from cleavage stage embryos, aiming towards an unaffected pregnancy. Twenty-two such couples were referred for PGD, and biopsied blastomeres on day 3 and untransferred embryos (day 5/6) were tested using fluorescence in situ hybridisation (FISH) with appropriate probes. This study investigated whether chromosome breakage (a) was detected more frequently in cases where the breakpoint of the aberration was in the same chromosomal band as a fragile site and (b) was influenced by maternal age, sperm parameters, reproductive history, or the sex of the carrier parent. The frequency of breakage seemed to be independent of fragile sites, maternal age, reproductive history, and sex of the carrier parent. However, chromosome breakage was very significantly higher in embryos from male carriers with poor sperm parameters versus embryos from male carriers with normal sperm parameters. Consequently, embryos from certain couples were more prone to chromosome breakage, fragment loss, and hence chromosomally unbalanced embryos, independently of meiotic segregation.     

The controlling sequence for site-specific chromosome breakage in Tetrahymena

Site-specific chromosome breakage occurs in many ciliated protozoa during nuclear differentiation. We have determined the cis-acting sequence that controls this process in Tetrahymena thermophila. The Tetrahymena ribosomal RNA gene is bounded by two breakage sites. Injection of this gene into developing macronuclei leads to breakage at these sites. Deletion analysis has localized the sequences essential for breakage to a 28 bp region that includes a 15 bp sequence (Cbs) known to be present in other breakage sites. Insertions of Cbs allow breakage to occur at new sites, which is accompanied by elimination of surrounding DNAs and formation of telomeric sequences, as it is at natural sites. Thus, Cbs is the necessary and sufficient sequence signal for chromosome breakage in Tetrahymena.     

A high resolution map of mammalian X chromosome fragile regions assessed by large-scale comparative genomics

Chromosomal evolution involves multiple changes at structural and numerical levels. These changes, which are related to the variation of the gene number and their location, can be tracked by the identification of syntenic blocks (SB). First reports proposed that ~180–280 SB might be shared by mouse and human species. More recently, further studies including additional genomes have identified up to ~1,400 SB during the evolution of eutherian species. A considerable number of studies regarding the X chromosome’s structure and evolution have been undertaken because of its extraordinary biological impact on reproductive fitness and speciation. Some have identified evolutionary breakpoint regions and fragile sites at specific locations in the human X chromosome. However, mapping these regions to date has involved using low-to-moderate resolution techniques. Such scenario might be related to underestimating their total number and giving an inaccurate location. The present study included using a combination of bioinformatics methods for identifying, at base-pair level, chromosomal rearrangements occurring during X chromosome evolution in 13 mammalian species. A comparative technique using four different algorithms was used for optimizing the detection of hotspot regions in the human X chromosome. We identified a significant interspecific variation in SB size which was related to genetic information gain regarding the human X chromosome. We found that human hotspot regions were enriched by LINE-1 and Alu transposable elements, which may have led to intraspecific chromosome rearrangement events. New fragile regions located in the human X chromosome have also been postulated. We estimate that the high resolution map of X chromosome fragile sites presented here constitutes useful data concerning future studies on mammalian evolution and human disease.