Determining the order of genes,centromeres, and rearrangement breakpoints inNeurospora by tests of duplication coverage

Nontandem segmental duplications provide a useful alternative to conventional recombination mapping for determining gene order in a haploid organism such asNeurospora. When an insertional or terminal rearrangement is crossed by Normal sequence, a class of progeny is produced that have a precisely delimited chromosome segment duplicated. In such Duplication progeny, a recessive gene in the Normal-sequence donor chromosome may or may not be masked (“covered”) by its dominant wild-type allele in the translocation-sequence recipient chromosome. Coverage depends upon whether the gene in question is left or right of the rearrangement breakpoint. The recessive gene will be heterozygous and covered (not expressed) if its locus is within the duplicated segment, but it will be haploid and expressed if the locus is outside the segment. Not only genes but also centromeres can be mapped by means of duplications, because genes included in. the same viable duplication must reside in the same chromosome arm. - Numerous sequences in the current genetic maps ofN. crassa have been determined using duplications. Gene order in the albino region and in the centromere region of linkage group I provide examples. Over 50 insertional or terminal rearrangements are available from which nontandem duplications of defined content can be obtained at will; collectively these cover about 75% of the genome. - Intercrosses between partially overlapping chromosome rearrangements also produce Duplication progeny containing two copies of regions between the breakpoints. The 180 mapped reciprocal translocations and inversions include numerous overlapping combinations that can be used for duplication mapping.     

Small deletion variants have stable breakpoints commonly associated with alu elements

Copy number variants (CNVs) contribute significantly to human genomic variation, with over 5000 loci reported, covering more than 18% of the euchromatic human genome. Little is known, however, about the origin and stability of variants of different size and complexity. We investigated the breakpoints of 20 small, common deletions, representing a subset of those originally identified by array CGH, using Agilent microarrays, in 50 healthy French Caucasian subjects. By sequencing PCR products amplified using primers designed to span the deleted regions, we determined the exact size and genomic position of the deletions in all affected samples. For each deletion studied, all individuals carrying the deletion share identical upstream and downstream breakpoints at the sequence level, suggesting that the deletion event occurred just once and later became common in the population. This is supported by linkage disequilibrium (LD) analysis, which has revealed that most of the deletions studied are in moderate to strong LD with surrounding SNPs, and have conserved long-range haplotypes. Analysis of the sequences flanking the deletion breakpoints revealed an enrichment of microhomology at the breakpoint junctions. More significantly, we found an enrichment of Alu repeat elements, the overwhelming majority of which intersected deletion breakpoints at their poly-A tails. We found no enrichment of LINE elements or segmental duplications, in contrast to other reports. Sequence analysis revealed enrichment of a conserved motif in the sequences surrounding the deletion breakpoints, although whether this motif has any mechanistic role in the formation of some deletions has yet to be determined. Considered together with existing information on more complex inherited variant regions, and reports of de novo variants associated with autism, these data support the presence of different subgroups of CNV in the genome which may have originated through different mechanisms.     

Next-Generation Sequencing of Duplication CNVs Reveals that Most Are Tandem and Some Create Fusion Genes at Breakpoints

Interpreting the genomic and phenotypic consequences of copy-number variation (CNV) is essential to understanding the etiology of genetic disorders. Whereas deletion CNVs lead obviously to haploinsufficiency, duplications might cause disease through triplosensitivity, gene disruption, or gene fusion at breakpoints. The mutational spectrum of duplications has been studied at certain loci, and in some cases these copy-number gains are complex chromosome rearrangements involving triplications and/or inversions. However, the organization of clinically relevant duplications throughout the genome has yet to be investigated on a large scale. Here we fine-mapped 184 germline duplications (14.7 kb–25.3 Mb; median 532 kb) ascertained from individuals referred for diagnostic cytogenetics testing. We performed next-generation sequencing (NGS) and whole-genome sequencing (WGS) to sequence 130 breakpoints from 112 subjects with 119 CNVs and found that most (83%) were tandem duplications in direct orientation. The remainder were triplications embedded within duplications (8.4%), adjacent duplications (4.2%), insertional translocations (2.5%), or other complex rearrangements (1.7%). Moreover, we predicted six in-frame fusion genes at sequenced duplication breakpoints; four gene fusions were formed by tandem duplications, one by two interconnected duplications, and one by duplication inserted at another locus. These unique fusion genes could be related to clinical phenotypes and warrant further study. Although most duplications are positioned head-to-tail adjacent to the original locus, those that are inverted, triplicated, or inserted can disrupt or fuse genes in a manner that might not be predicted by conventional copy-number assays. Therefore, interpreting the genetic consequences of duplication CNVs requires breakpoint-level analysis.     

Genome-level evolution of resistance genes in Arabidopsis thaliana

Pathogen resistance genes represent some of the most abundant and diverse gene families found within plant genomes. However, evolutionary mechanisms generating resistance gene diversity at the genome level are not well understood. We used the complete Arabidopsis thaliana genome sequence to show that most duplication of individual NBS-LRR sequences occurs at close physical proximity to the parent sequence and generates clusters of closely related NBS-LRR sequences. Deploying the statistical strength of phylogeographic approaches and using chromosomal location as a proxy for spatial location, we show that apparent duplication of NBS-LRR genes to ectopic chromosomal locations is largely the consequence of segmental chromosome duplication and rearrangement, rather than the independent duplication of individual sequences. Although accounting for a smaller fraction of NBS-LRR gene duplications, segmental chromosome duplication and rearrangement events have a large impact on the evolution of this multigene family. Intergenic exchange is dramatically lower between NBS-LRR sequences located in different chromosome regions as compared to exchange between sequences within the same chromosome region. Consequently, once translocated to new chromosome locations, NBS-LRR gene copies have a greater likelihood of escaping intergenic exchange and adopting new functions than do gene copies located within the same chromosomal region. We propose an evolutionary model that relates processes of genome evolution to mechanisms of evolution for the large, diverse, NBS-LRR gene family.     

Recent duplication of the common carp (Cyprinus carpio L.) genome as revealed by analyses of microsatellite loci

Genome duplications may have played a role in the early stages of vertebrate evolution, near the time of divergence of the lamprey lineage. Additional genome duplication, specifically in ray-finned fish, may have occurred before the divergence of the teleosts. The common carp (Cyprinus carpio) has been considered tetraploid because of its chromosome number (2n = 100) and its high DNA content. We studied variation using 59 microsatellite primer pairs to better understand the ploidy level of the common carp. Based on the number of PCR amplicons per individual, about 60% of these primer pairs are estimated to amplify duplicates. Segregation patterns in families suggested a partially duplicated genome structure and disomic inheritance. This could suggest that the common carp is tetraploid and that polyploidy occurred by hybridization (allotetraploidy). From sequences of microsatellite flanking regions, we estimated the difference per base between pairs of alleles and between pairs of paralogs. The distribution of differences between paralogs had two distinct modes suggesting one whole-genome duplication and a more recent wave of segmental duplications. The genome duplication was estimated to have occurred about 12 MYA, with the segmental duplications occurring between 2.3 and 6.8 MYA. At 12 MYA, this would be one of the most recent genome duplications among vertebrates. Phylogenetic analysis of several cyprinid species suggests an evolutionary model for this tetraploidization, with a role for polyploidization in speciation and diversification.     

Genomic rearrangements resulting in PLP1 deletion occur by nonhomologous end joining and cause different dysmyelinating phenotypes in males and females

In the majority of patients with Pelizaeus-Merzbacher disease, duplication of the proteolipid protein gene PLP1 is responsible, whereas deletion of PLP1 is infrequent. Genomic mechanisms for these submicroscopic chromosomal rearrangements remain unknown. We identified three families with PLP1 deletions (including one family described elsewhere) that arose by three distinct processes. In one family, PLP1 deletion resulted from a maternal balanced submicroscopic insertional translocation of the entire PLP1 gene to the telomere of chromosome 19. PLP1 on the 19qtel is probably inactive by virtue of a position effect, because a healthy male sibling carries the same der(19) chromosome along with a normal X chromosome. Genomic mapping of the deleted segments revealed that the deletions are smaller than most of the PLP1 duplications and involve only two other genes. We hypothesize that the deletion is infrequent, because only the smaller deletions can avoid causing either infertility or lethality. Analyses of the DNA sequence flanking the deletion breakpoints revealed Alu-Alu recombination in the family with translocation. In the other two families, no homologous sequence flanking the breakpoints was found, but the distal breakpoints were embedded in novel low-copy repeats, suggesting the potential involvement of genome architecture in stimulating these rearrangements. In one family, junction sequences revealed a complex recombination event. Our data suggest that PLP1 deletions are likely caused by nonhomologous end joining.     

Large Inverted Duplications in the Human Genome Form via a Fold-Back Mechanism

Inverted duplications are a common type of copy number variation (CNV) in germline and somatic genomes. Large duplications that include many genes can lead to both neurodevelopmental phenotypes in children and gene amplifications in tumors. There are several models for inverted duplication formation, most of which include a dicentric chromosome intermediate followed by breakage-fusion-bridge (BFB) cycles, but the mechanisms that give rise to the inverted dicentric chromosome in most inverted duplications remain unknown. Here we have combined high-resolution array CGH, custom sequence capture, next-generation sequencing, and long-range PCR to analyze the breakpoints of 50 nonrecurrent inverted duplications in patients with intellectual disability, autism, and congenital anomalies. For half of the rearrangements in our study, we sequenced at least one breakpoint junction. Sequence analysis of breakpoint junctions reveals a normal-copy disomic spacer between inverted and non-inverted copies of the duplication. Further, short inverted sequences are present at the boundary of the disomic spacer and the inverted duplication. These data support a mechanism of inverted duplication formation whereby a chromosome with a double-strand break intrastrand pairs with itself to form a “fold-back” intermediate that, after DNA replication, produces a dicentric inverted chromosome with a disomic spacer corresponding to the site of the fold-back loop. This process can lead to inverted duplications adjacent to terminal deletions, inverted duplications juxtaposed to translocations, and inverted duplication ring chromosomes.     

Patterns of segmental duplication in the human genome

We analyzed the completed human genome for recent segmental duplications (size > or = 1 kb and sequence similarity > or = 90%). We found that approximately 4% of the genome is covered by duplications and that the extent of segmental duplication varies from 1% to 14% among the 24 chromosomes. Intrachromosomal duplication is more frequent than interchromosomal duplication in 15 chromosomes. The duplication frequencies in pericentromeric and subtelomeric regions are greater than the genome average by approximately threefold and fourfold. We examined factors that may affect the frequency of duplication in a region. Within individual chromosomes, the duplication frequency shows little correlation with local gene density, repeat density, recombination rate, and GC content, except chromosomes 7 and Y. For the entire genome, the duplication frequency is correlated with each of the above factors. Based on known genes and Ensembl genes, the proportion of duplications containing complete genes is 3.4% and 10.7%, respectively. The proportion of duplications containing genes is higher in intrachromosomal than in interchromosomal duplications, and duplications containing genes have a higher sequence similarity and tend to be longer than duplications containing no genes. Our simulation suggests that many duplications containing genes have been selectively maintained in the genome.     

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.     

The Effects of Translocations on Recombination Frequency in Caenorhabditis Elegans

In the nematode Caenorhabditis elegans, recombination suppression in translocation heterozygotes is severe and extensive. We have examined the meiotic properties of two translocations involving chromosome I, szT1(I;X) and hT1(I;V). No recombination was observed in either of these translocation heterozygotes along the left (let-362-unc-13) 17 map units of chromosome I. Using half-translocations as free duplications, we mapped the breakpoints of szT1 and hT1. The boundaries of crossover suppression coincided with the physical breakpoints. We propose that DNA sequences at the right end of chromosome I facilitate pairing and recombination. We use the data from translocations of other chromosomes to map the location of pairing sites on four other chromosomes. hT1 and szT1 differed markedly in their effect on recombination adjacent to the crossover suppressed region. hT1 had no effect on recombination in the adjacent interval. In contrast, the 0.8 map unit interval immediately adjacent to the szT1(I;X) breakpoint on chromosome I increased to 2.5 map units in translocation heterozygotes. This increase occurs in a chromosomal interval which can be expanded by treatment with radiation. These results are consistent with the suggestion that the szT1(I) breakpoint is in a region of DNA in which meiotic recombination is suppressed relative to the genomic average. We propose that DNA sequences disrupted by the szT1 translocation are responsible for determining the frequency of meiotic recombination in the vicinity of the breakpoint.     

Eucaryotic genome evolution through the spontaneous duplication of large chromosomal segments

There is growing evidence that duplications have played a major role in eucaryotic genome evolution. Sequencing data revealed the presence of large duplicated regions in the genomes of many eucaryotic organisms, and comparative studies have suggested that duplication of large DNA segments has been a continuing process during evolution. However, little experimental data have been produced regarding this issue. Using a gene dosage assay for growth recovery in Saccharomyces cerevisiae, we demonstrate that a majority of the revertant strains (58%) resulted from the spontaneous duplication of large DNA segments, either intra- or interchromosomally, ranging from 41 to 655 kb in size. These events result in the concomitant duplication of dozens of genes and in some cases in the formation of chimeric open reading frames at the junction of the duplicated blocks. The types of sequences at the breakpoints as well as their superposition with the replication map suggest that spontaneous large segmental duplications result from replication accidents. Aneuploidization events or suppressor mutations that do not involve large-scale rearrangements accounted for the rest of the reversion events (in 26 and 16% of the strains, respectively).     

Plant noncoding RNA gene discovery by "single-genome comparative genomics"

Plant genomes have undergone multiple rounds of duplications that contributed massively to the growth of gene families. The structure of resulting families has been studied in depth for protein-coding genes. However, little is known about the impact of duplications on noncoding RNA (ncRNA) genes. Here we perform a systematic analysis of duplicated regions in the rice genome in search of such ncRNA repeats. We observe that, just like their protein counterparts, most ncRNA genes have undergone multiple duplications that left visible sequence conservation footprints. The extent of ncRNA gene duplication in plants is such that these sequence footprints can be exploited for the discovery of novel ncRNA gene families on a large scale. We developed an SVM model that is able to retrieve likely ncRNA candidates among the 100,000+ repeat families in the rice genome, with a reasonably low false-positive discovery rate. Among the nearly 4000 ncRNA families predicted by this means, only 90 correspond to putative snoRNA or miRNA families. About half of the remaining families are classified as structured RNAs. New candidate ncRNAs are particularly enriched in UTR and intronic regions. Interestingly, 89% of the putative ncRNA families do not produce a detectable signal when their sequences are compared to another grass genome such as maize. Our results show that a large fraction of rice ncRNA genes are present in multiple copies and are species-specific or of recent origin. Intragenome comparison is a unique and potent source for the computational annotation of this major class of ncRNA.     

Resolving the breakpoints of the 17q21.31 microdeletion syndrome with next-generation sequencing

Recurrent deletions have been associated with numerous diseases and genomic disorders. Few, however, have been resolved at the molecular level because their breakpoints often occur in highly copy-number-polymorphic duplicated sequences. We present an approach that uses a combination of somatic cell hybrids, array comparative genomic hybridization, and the specificity of next-generation sequencing to determine breakpoints that occur within segmental duplications. Applying our technique to the 17q21.31 microdeletion syndrome, we used genome sequencing to determine copy-number-variant breakpoints in three deletion-bearing individuals with molecular resolution. For two cases, we observed breakpoints consistent with nonallelic homologous recombination involving only H2 chromosomal haplotypes, as expected. Molecular resolution revealed that the breakpoints occurred at different locations within a 145 kbp segment of >99% identity and disrupt KANSL1 (previously known as KANSL1). In the remaining case, we found that unequal crossover occurred interchromosomally between the H1 and H2 haplotypes and that this event was mediated by a homologous sequence that was once again missing from the human reference. Interestingly, the breakpoints mapped preferentially to gaps in the current reference genome assembly, which we resolved in this study. Our method provides a strategy for the identification of breakpoints within complex regions of the genome harboring high-identity and copy-number-polymorphic segmental duplication. The approach should become particularly useful as high-quality alternate reference sequences become available and genome sequencing of individuals'' DNA becomes more routine.     

Meiosis-driven genome variation in plants

Meiosis includes two successive divisions of the nucleus with one round of DNA replication and leads to the formation of gametes with half of the chromosomes of the mother cell during sexual reproduction. It provides a cytological basis for gametogenesis and nheritance in eukaryotes. Meiotic cell division is a complex and dynamic process that involves a number of molecular and cellular events, such as DNA and chromosome replication, chromosome pairing, synapsis and recombination, chromosome segregation, and cytokinesis. Meiosis maintains genome stability and integrity over sexual life cycles. On the other hand, meiosis generates genome variations in several ways. Variant meiotic recombination resulting from specific genome structures induces deletions, duplications, and other rearrangements within the genic and non-genic genomic regions and has been considered a major driving force for gene and genome evolution in nature. Meiotic abnormalities in chromosome segregation lead to chromosomally imbalanced gametes and aneuploidy. Meiotic restitution due to failure of the first or second meiotic division gives rise to unreduced gametes, which triggers polyploidization and genome expansion. This paper reviews research regarding meiosis-driven genome variation, including deletion and duplication of genomic regions, aneuploidy, and polyploidization, and discusses the effect of related meiotic events on genome variation and evolution in plants. Knowledge of various meiosis-driven genome variations provides insight into genome evolution and genetic variability in plants and facilitates plant genome research.     

Analysis of homologous gene clusters in Caenorhabditis elegans reveals striking regional cluster domains

An algorithm for detecting local clusters of homologous genes was applied to the genome of Caenorhabditis elegans. Clusters of two or more homologous genes are abundant, totaling 1391 clusters containing 4607 genes, over one-fifth of all genes in C. elegans. Cluster genes are distributed unevenly in the genome, with the large majority located on autosomal chromosome arms, regions characterized by higher genetic recombination and more repeat sequences than autosomal centers and the X chromosome. Cluster genes are transcribed at much lower levels than average and very few have gross phenotypes as assayed by RNAi-mediated reduction of function. The molecular identity of cluster genes is unusual, with a preponderance of nematode-specific gene families that encode putative secreted and transmembrane proteins, and enrichment for genes implicated in xenobiotic detoxification and innate immunity. Gene clustering in Drosophila melanogaster is also substantial and the molecular identity of clustered genes follows a similar pattern. I hypothesize that autosomal chromosome arms in C. elegans undergo frequent local gene duplication and that these duplications support gene diversification and rapid evolution in response to environmental challenges. Although specific gene clusters have been documented in C. elegans, their abundance, genomic distribution, and unusual molecular identities were previously unrecognized.     

Localization and characterization of X chromosome inversion breakpoints separating Drosophila mojavensis and Drosophila arizonae

Ectopic exchange between transposable elements or other repetitive sequences along a chromosome can produce chromosomal inversions. As a result, genome sequence studies typically find sequence similarity between corresponding inversion breakpoint regions. Here, we identify and investigate the breakpoint regions of the X chromosome inversion distinguishing Drosophila mojavensis and Drosophila arizonae. We localize one inversion breakpoint to 13.7 kb and localize the other to a 1-Mb interval. Using this localization and assuming microsynteny between Drosophila melanogaster and D. arizonae, we pinpoint likely positions of the inversion breakpoints to windows of less than 3000 bp. These breakpoints define the size of the inversion to approximately 11 Mb. However, in contrast to many other studies, we fail to find significant sequence similarity between the 2 breakpoint regions. The localization of these inversion breakpoints will facilitate future genetic and molecular evolutionary studies in this species group, an emerging model system for ecological genetics.     

Drosophila duplication hotspots are associated with late-replicating regions of the genome

Duplications play a significant role in both extremes of the phenotypic spectrum of newly arising mutations: they can have severe deleterious effects (e.g. duplications underlie a variety of diseases) but can also be highly advantageous. The phenotypic potential of newly arisen duplications has stimulated wide interest in both the mutational and selective processes shaping these variants in the genome. Here we take advantage of the Drosophila simulans-Drosophila melanogaster genetic system to further our understanding of both processes. Regarding mutational processes, the study of two closely related species allows investigation of the potential existence of shared duplication hotspots, and the similarities and differences between the two genomes can be used to dissect its underlying causes. Regarding selection, the difference in the effective population size between the two species can be leveraged to ask questions about the strength of selection acting on different classes of duplications. In this study, we conducted a survey of duplication polymorphisms in 14 different lines of D. simulans using tiling microarrays and combined it with an analogous survey for the D. melanogaster genome. By integrating the two datasets, we identified duplication hotspots conserved between the two species. However, unlike the duplication hotspots identified in mammalian genomes, Drosophila duplication hotspots are not associated with sequences of high sequence identity capable of mediating non-allelic homologous recombination. Instead, Drosophila duplication hotspots are associated with late-replicating regions of the genome, suggesting a link between DNA replication and duplication rates. We also found evidence supporting a higher effectiveness of selection on duplications in D. simulans than in D. melanogaster. This is also true for duplications segregating at high frequency, where we find evidence in D. simulans that a sizeable fraction of these mutations is being driven to fixation by positive selection.     

Chromatin structural elements and chromosomal translocations in leukemia

Recurring chromosome abnormalities are strongly associated with certain subtypes of leukemia, lymphoma and sarcomas. More recently, their potential involvement in carcinomas, i.e. prostate cancer, has been recognized. They are among the most important factors in determining disease prognosis, and in many cases, identification of these chromosome abnormalities is crucial in selecting appropriate treatment protocols. Chromosome translocations are frequently observed in both de novo and therapy-related acute myeloid leukemia (AML) and myelodysplastic syndromes (MDS). The mechanisms that result in such chromosome translocations in leukemia and other cancers are largely unknown. Genomic breakpoints in all the common chromosome translocations in leukemia, including t(4;11), t(9;11), t(8;21), inv(16), t(15;17), t(12;21), t(1;19) and t(9;22), have been cloned. Genomic breakpoints tend to cluster in certain intronic regions of the relevant genes including MLL, AF4, AF9, AML1, ETO, CBFB, MYHI1, PML, RARA, TEL, E2A, PBX1, BCR and ABL. However, whereas the genomic breakpoints in MLL tend to cluster in the 5' portion of the 8.3 kb breakpoint cluster region (BCR) in de novo and adult patients and in the 3' portion in infant leukemia patients and t-AML patients, those in both the AML1 and ETO genes occur in the same clustered regions in both de novo and t-AML patients. These differences may reflect differences in the mechanisms involved in the formation of the translocations. Specific chromatin structural elements, such as in vivo topoisomerase II (topo II) cleavage sites, DNase I hypersensitive sites and scaffold attachment regions (SARs) have been mapped in the breakpoint regions of the relevant genes. Strong in vivo topo II cleavage sites and DNase I hypersensitive sites often co-localize with each other and also with many of the BCRs in most of these genes, whereas SARs are associated with BCRs in MLL, AF4, AF9, AML1, ETO and ABL, but not in the BCR gene. In addition, the BCRs in MLL, AML1 and ETO have the lowest free energy level for unwinding double strand DNA. Virtually all chromosome translocations in leukemia that have been analyzed to date show no consistent homologous sequences at the breakpoints, whereas a strong non-homologous end joining (NHEJ) repair signature exists at all of these chromosome translocation breakpoint junctions; this includes small deletions and duplications in each breakpoint, and micro-homologies and non-template insertions at genomic junctions of each chromosome translocation. Surprisingly, the size of these deletions and duplications in the same translocation is much larger in de novo leukemia than in therapy-related leukemia. We propose a non-homologous chromosome recombination model as one of the mechanisms that results in chromosome translocations in leukemia. The topo II cleavage sites at open chromatin regions (DNase I hypersensitive sites), SARs or the regions with low energy level are vulnerable to certain genotoxic or other agents and become the initial breakage sites, which are followed by an excision end joining repair process.