STRUCTURE AND POPULATION GENETICS OF THE BREAKPOINTS OF A POLYMORPHIC INVERSION IN DROSOPHILA SUBOBSCURA

Drosophila subobscura is a paleartic species of the obscura group with a rich chromosomal polymorphism. To further our understanding on the origin of inversions and on how they regain variation, we have identified and sequenced the two breakpoints of a polymorphic inversion of D. subobscura—inversion 3 of the O chromosome—in a population sample. The breakpoints could be identified as two rather short fragments (～300 bp and 60 bp long) with no similarity to any known transposable element family or repetitive sequence. The presence of the ～300‐bp fragment at the two breakpoints of inverted chromosomes implies its duplication, an indication of the inversion origin via staggered double‐strand breaks. Present results and previous findings support that the mode of origin of inversions is neither related to the inversion age nor species‐group specific. The breakpoint regions do not consistently exhibit the lower level of variation within and stronger genetic differentiation between arrangements than more internal regions that would be expected, even in moderately small inversions, if gene conversion were greatly restricted at inversion breakpoints. Comparison of the proximal breakpoint region in species of the obscura group shows that this breakpoint lies in a small high‐turnover fragment within a long collinear region (～300 kb).     

Nucleotide sequence comparison of a chromosome rearrangement on human chromosome 12 and the corresponding ape chromosomes

Chromosome rearrangement has been considered to be important in the evolutionary process. Here, we demonstrate the evolutionary relationship of the rearranged human chromosome 12 and the corresponding chromosome XII in apes (chimpanzee, bonobo, gorilla, orangutan, and gibbon) by examining PCR products derived from the breakpoints of inversions and by conducting shotgun sequencing of a gorilla fosmid clone containing the breakpoint and a "duplicated segment" (duplicon). We confirmed that a pair of 23-kb duplicons flank the breakpoints of inversions on the long and short arms of chimpanzee chromosome XII. Although only the 23-kb duplicon on the long arm of chimpanzee chromosome XII and its telomeric flanking sequence are found to be conserved among the hominoids (human, great apes, and gibbons), the duplicon on the short arm of chimpanzee chromosome XII is suggested to be the result of a duplication from that on the long arm. Furthermore, the shotgun sequencing of a gorilla fosmid indicated that the breakpoint on the long arm of the gorilla is located at a different position 1.9 kb from that of chimpanzee. The region is flanked by a sequence homologous to that of human chromosome 6q22. Our findings and sequence analysis suggest a close relationship between segmental duplication and chromosome rearrangement (or breakpoint of inversion) in Hominoidea. The role of the chromosome rearrangement in speciation is also discussed based on our new results.     

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.     

The foldback-like transposon Galileo is involved in the generation of two different natural chromosomal inversions of Drosophila buzzatii

Chromosomal inversions are the most common type of genome rearrangement in the genus Drosophila. Although the potential of transposable elements (TEs) for generating inversions has been repeatedly demonstrated in the laboratory, little is known on their role in the generation of natural inversions, which are those effectively contributing to the adaptation and/or evolution of species. We have cloned and sequenced the two breakpoints of the polymorphic inversion 2q7 of D. buzzatii. The sequence analysis of the breakpoint regions revealed the presence in the inverted chromosomes of large insertions, formed by complex assemblies of transposons, that are absent from the chromosomes without the inversion. Among the transposons inserted, the Foldback-like element Galileo, that was previously found responsible of the generation of the widespread inversion 2j of D. buzzatii, is present at both 2q7 breakpoints and is the most likely inducer of the inversion. A detailed study of the nucleotide and structural variation in the breakpoint regions of six chromosomal lines with the 2q7 inversion detected no nucleotide differences between them, which suggests a monophyletic and recent origin. In contrast, a remarkable degree of structural variation was observed in the same six chromosomal lines. It thus appears that the two breakpoints of the inverted chromosomes have become genetically unstable hotspots, as was previously found for the 2j inversion breakpoints. The possibility that this instability is caused by structural properties of Foldback elements is discussed.     

Testing chromosomal phylogenies and inversion breakpoint reuse in Drosophila

A combination of cytogenetic and bioinformatic procedures was used to test the chromosomal phylogeny relating Drosophila buzzatii with D. repleta. Chromosomes X and 2, harboring most of the inversions fixed between these two species, were analyzed. First, chromosomal segments conserved during the divergence of the two species were identified by comparative in situ hybridization to the D. repleta chromosomes of 180 BAC clones from a BAC-based physical map of the D. buzzatii genome. These conserved segments were precisely delimited with the aid of clones containing inversion breakpoints. Then GRIMM software was used to estimate the minimum number of rearrangements necessary to transform one genome into the other and identify all possible rearrangement scenarios. Finally, the most plausible inversion trajectory was tested by hybridizing 12 breakpoint-bearing BAC clones to the chromosomes of seven other species in the repleta group. The results show that chromosomes X and 2 of D. buzzatii and D. repleta differ by 12 paracentric inversions. Nine of them are fixed in chromosome 2 and entail two breakpoint reuses. Our results also show that the cytological relationship between D. repleta and D. mercatorum is closer than that between D. repleta and D. peninsularis, and we propose that the phylogenetic relationships in this lineage of the repleta group be reconsidered. We also estimated the rate of rearrangement between D. repleta and D. buzzatii and conclude that rates within the genus Drosophila vary substantially between lineages, even within a single species group.     

The evolutionary history of Drosophila buzzatii. XXXV. Inversion polymorphism and nucleotide variability in different regions of the second chromosome

Inversions are portions of a chromosome where the gene order is reversed relative to a standard reference orientation. Because of reduced levels of recombination in heterokaryotypes, inversions have a potentially important effect on patterns of nucleotide variability in those genomic regions close to, or included in, the inverted fragments. Here we report sequence variation at three anonymous regions (STSs) located at different positions in relation to second-chromosome inversion breakpoints in 29 isochromosomal lines derived from an Argentinean population of Drosophila buzzatii. In agreement with previous findings in Drosophila, gene flux (crossing over and/or gene conversion) between arrangements seems to appreciably increase as we approach the middle sections of inversion 2j, and patterns of nucleotide variability within, as well as genetic differentiation between chromosome arrangements, are comparable to those observed at the molecular marker outside the inverted fragments. On the other hand, nucleotide diversity near the proximal breakpoint of inversion 2j is reduced when contrasted with that found at the other regions, particularly in the case of derived inverted chromosomes. Using the data from the breakpoint, we estimate that the inversion polymorphism is approximately 1.63 N generations old, where N is the effective population size. An excess of low-frequency segregating polymorphisms is detected; mostly in the ancestral 2st arrangement and probably indicating a population expansion that predates the coalescent time of inversion 2j. Heterogeneity in mutation rates between the markers linked to the inversions may be sufficient to explain the different levels of nucleotide diversity observed. When considered in the context of other studies on patterns of variation relative to physical distance to inversion breakpoints, our data appear to be consistent with the conclusion that inversions are unlikely to be "long-lived" balanced polymorphisms.     

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.     

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.     

The evolutionary history of human chromosome 7

We report on a comparative molecular cytogenetic and in silico study on evolutionary changes in human chromosome 7 homologs in all major primate lineages. The ancestral mammalian homologs comprise two chromosomes (7a and 7b/16p) and are conserved in carnivores. The subchromosomal organization of the ancestral primate segment 7a shared by a lemur and higher Old World monkeys is the result of a paracentric inversion. The ancestral higher primate chromosome form was then derived by a fission of 7b/16p, followed by a centric fusion of 7a/7b as observed in the orangutan. In hominoids two further inversions with four distinct breakpoints were described in detail: the pericentric inversion in the human/African ape ancestor and the paracentric inversion in the common ancestor of human and chimpanzee. FISH analysis employing BAC probes confined the 7p22.1 breakpoint of the pericentric inversion to 6.8 Mb on the human reference sequence map and the 7q22.1 breakpoint to 97.1 Mb. For the paracentric inversion the breakpoints were found in 7q11.23 between 76.1 and 76.3 Mb and in 7q22.1 at 101.9 Mb. All four breakpoints were flanked by large segmental duplications. Hybridization patterns of breakpoint-flanking BACs and the distribution of duplicons suggest their presence before the origin of both inversions. We propose a scenario by which segmental duplications may have been the cause rather than the result of these chromosome rearrangements.     

Philadelphia chromosomal breakpoints are clustered within a limited region, bcr, on chromosome 22

We have identified and molecularly cloned 46 kb of human DNA from chromosome 22 using a probe specific for the Philadelphia (Ph') translocation breakpoint domain of one chronic myelocytic leukemia (CML) patient. The DNAs of 19 CML patients were examined for rearrangements on chromosome 22 with probes isolated from this cloned region. In 17 patients, chromosomal breakpoints were found within a limited region of up to 5.8 kb, for which we propose the term "breakpoint cluster region" (bcr). The two patients having no rearrangements within bcr lacked the Ph' chromosome. The highly specific presence of a chromosomal breakpoint within bcr in Ph'-positive CML patients strongly suggests the involvement of bcr in this type of leukemia.     

Molecular Characterization of Hobo-Mediated Inversions in Drosophila Melanogaster

The structure of chromosomal inversions mediated by hobo transposable elements in the Uc-1 X chromosome was investigated using cytogenetic and molecular methods. Uc-1 contains a phenotypically silent hobo element inserted in an intron of the Notch locus. Cytological screening identified six independent Notch mutations resulting from chromosomal inversions with one breakpoint at cytological position 3C7, the location of Notch. In situ hybridization to salivary gland polytene chromosomes determined that both ends of each inversion contained hobo and Notch sequences. Southern blot analyses showed that both breakpoints in each inversion had hobo-Notch junction fragments indistinguishable in structure from those present in the Uc-1 X chromosome prior to the rearrangements. Polymerase chain reaction amplification of the 12 hobo-Notch junction fragments in the six inversions, followed by DNA sequence analysis, determined that each was identical to one of the two hobo-Notch junctions present in Uc-1. These results are consistent with a model in which hobo-mediated inversions result from homologous pairing and recombination between a pair of hobo elements in reverse orientation.     

A transgenic system for generation of transposon Ac/Ds-induced chromosome rearrangements in rice

The maize Activator (Ac)/Dissociation (Ds) transposable element system has been used in a variety of plants for insertional mutagenesis. Ac/Ds elements can also generate genome rearrangements via alternative transposition reactions which involve the termini of closely linked transposons. Here, we introduced a transgene containing reverse-oriented Ac/Ds termini together with an Ac transposase gene into rice (Oryza sativa ssp. japonica cv. Nipponbare). Among the transgenic progeny, we identified and characterized 25 independent genome rearrangements at three different chromosomal loci. The rearrangements include chromosomal deletions and inversions and one translocation. Most of the deletions occurred within the T-DNA region, but two cases showed the loss of 72 kilobase pairs (kb) and 79 kb of rice genomic DNA flanking the transgene. In addition to deletions, we obtained chromosomal inversions ranging in size from less than 10 kb (within the transgene DNA) to over 1 million base pairs (Mb). For 11 inversions, we cloned and sequenced both inversion breakpoints; in all 11 cases, the inversion junctions contained the typical 8 base pairs (bp) Ac/Ds target site duplications, confirming their origin as transposition products. Together, our results indicate that alternative Ac/Ds transposition can be an efficient tool for functional genomics and chromosomal manipulation in rice.     

Assaying chromosomal inversions by single-molecule haplotyping

Inversions are an important form of structural variation, but they are difficult to characterize, as their breakpoints often fall within inverted repeats. We have developed a method called 'haplotype fusion' in which an inversion breakpoint is genotyped by performing fusion PCR on single molecules of human genomic DNA. Fusing single-copy sequences bracketing an inversion breakpoint generates orientation-specific PCR products, exemplified by a genotyping assay for the int22 hemophilia A inversion on Xq28. Furthermore, we demonstrated that inversion events with breakpoints embedded within long (>100 kb) inverted repeats can be genotyped by haplotype-fusion PCR followed by bead-based single-molecule haplotyping on repeat-specific markers bracketing the inversion breakpoint. We illustrate this method by genotyping a Yp paracentric inversion sponsored by >300-kb-long inverted repeats. The generality of our methods to survey for, and genotype chromosomal inversions should help our understanding of the contribution of inversions to genomic variation, inherited diseases and cancer.     

Molecular characterization and chromosomal distribution of Galileo, Kepler and Newton, three foldback transposable elements of the Drosophila buzzatii species complex

Galileo is a foldback transposable element that has been implicated in the generation of two polymorphic chromosomal inversions in Drosophila buzzatii. Analysis of the inversion breakpoints led to the discovery of two additional elements, called Kepler and Newton, sharing sequence and structural similarities with Galileo. Here, we describe in detail the molecular structure of these three elements, on the basis of the 13 copies found at the inversion breakpoints plus 10 additional copies isolated during this work. Similarly to the foldback elements described in other organisms, these elements have long inverted terminal repeats, which in the case of Galileo possess a complex structure and display a high degree of internal variability between copies. A phylogenetic tree built with their shared sequences shows that the three elements are closely related and diverged approximately 10 million years ago. We have also analyzed the abundance and chromosomal distribution of these elements in D. buzzatii and other species of the repleta group by Southern analysis and in situ hybridization. Overall, the results suggest that these foldback elements are present in all the buzzatti complex species and may have played an important role in shaping their genomes. In addition, we show that recombination rate is the main factor determining the chromosomal distribution of these elements.     

Chromosomal distributions of breakpoints in cancer,infertility, and evolution

We extract 11 genome-wide sets of breakpoint positions from databases on reciprocal translocations, inversions and deletions in neoplasms, reciprocal translocations and inversions in families carrying rearrangements and the human-mouse comparative map, and for each set of positions construct breakpoint distributions for the 44 autosomal arms. We identify and interpret four main types of distribution: (i) a uniform distribution associated both with families carrying translocations or inversions, and with the comparative map, (ii) telomerically skewed distributions of translocations or inversions detected consequent to births with malformations, (iii) medially clustered distributions of translocation and deletion breakpoints in tumor karyotypes, and (iv) bimodal translocation breakpoint distributions for chromosome arms containing telomeric proto-oncogenes.