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.     

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.     

Abundance and chromosomal distribution of six Drosophila buzzatii transposons: BuT1, BuT2, BuT3, BuT4, BuT5, and BuT6

The abundance and chromosomal distribution of six class-II transposable elements (TEs) of Drosophila buzzatii have been analyzed by Southern blotting and in situ hybridization. These six transposons had been previously found at the breakpoints of inversions 2j and 2q ⁷ of D. buzzatii. These two polymorphic inversions were generated by an ectopic recombination event between two copies of Galileo, a Foldback element. The four breakpoints became hotspots for TE insertions after the generation of the inversion and the transposons analyzed in this work are considered to be secondary invaders of these regions. Insertions of the six transposons are present in the euchromatin but show an increased density in the pericentromeric euchromatin–heterochromatin transition region and the dot chromosome. They are also more abundant in the inverted segments of chromosome 2 rearrangements. We further observed that the accumulation of TE insertions varies between elements and is correlated between dot, proximal regions, and inverted segments. These observations fully agree with previous data in Drosophila melanogaster and support recombination rate as the chief force explaining the chromosomal distribution of TEs.Electronic Supplementary Material Supplementary material is available in the online version of this article at and is accessible for authorized users.Sequence data from this article have been deposited in the EMBL/GenBank Data Libraries under accession number DQ402469.     

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.     

Evaluation of the genomic extent of effects of fixed inversion differences on intraspecific variation and interspecific gene flow in Drosophila pseudoobscura and D. persimilis

There is increasing evidence that chromosomal inversions may facilitate the formation or persistence of new species by allowing genetic factors conferring species-specific adaptations or reproductive isolation to be inherited together and by reducing or eliminating introgression. However, the genomic domain of influence of the inverted regions on introgression has not been carefully studied. Here, we present a detailed study on the consequences that distance from inversion breakpoints has had on the inferred level of gene flow and divergence between Drosophila pseudoobscura and D. persimilis. We identified the locations of the inversion breakpoints distinguishing D. pseudoobscura and D. persimilis in chromosomes 2, XR, and XL. Population genetic data were collected at specific distances from the inversion breakpoints of the second chromosome and at two loci inside the XR and XL inverted regions. For loci outside the inverted regions, we found that distance from the nearest inversion breakpoint had a significant effect on several measures of divergence and gene flow between D. pseudoobscura and D. persimilis. The data fitted a logarithmic relationship, showing that the suppression of crossovers in inversion heterozygotes also extends to loci located outside the inversion but close to it (within 1-2 Mb). Further, we detected a significant reduction in nucleotide variation inside the inverted second chromosome region of D. persimilis and near one breakpoint, consistent with a scenario in which this inversion arose and was fixed in this species by natural selection.     

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.     

Evolution versus constitution: differences in chromosomal inversion

We compared the chromosomal breakpoints of evolutionary conserved and constitutional inversions. Multicolor banding and human-specific bacterial artificial chromosomes were applied to map the breakpoints of constitutional pericentric inversions on human chromosomes 2 and 9. For the first time, we present a high-resolution analysis of the breakpoint regions, which are characterized by gene destitution, co-localization with fragile sites, multitude repeats as well as pseudogenes and, remarkably, a large sequence homology to the opposite breakpoint. In contrast, evolutionary inversion breakpoints lack such extensive cross-hybridizing regions and are often associated with fragile sites of the genome and low-copy repeats. These molecular characteristics gave evidence for different types of inversion formation and indicate that evolutionary inversions cannot originate from constitutional inversions like those of chromosomes 2 and 9. Finally, the constitutional inversion breakpoints were investigated on three different great ape species and on four test persons each bearing the same cytogenetically determined inversion on chromosomes 2 and 9, respectively. Our data indicate the existence of different molecular breakpoints for the two variant chromosomes.     

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.     

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).     

Are chromosomal inversions induced by transposable elements? A paradigm from the malaria mosquito Anopheles gambiae

Chromosomal rearrangements abound in nature and can be studied in detail in organisms with polytene chromosomes. In Drosophila and in Anopheline mosquitoes most speciation processes seem to be associated with the establishment of chromosomal rearrangements, particularly of paracentric inversions. It is not known what triggers inversions in natural populations. In the laboratory inversions are commonly generated by X-rays, mutagens or after the activity of certain transposable elements (TEs). The Anopheles gambiae complex is comprised of six sibling species, each one characterized by the presence of fixed paracentric inversions on their chromosomes. Two of these, An. gambiae s.s. and An. arabiensis, are the most important vectors of human malaria and are structured into sub-populations, each carrying a characteristic set of polymorphic chromosomal inversions. We have cloned the breakpoints of the naturally occurring polymorphic inversion In(2R)d' of An. arabiensis. Analysis of the surrounding sequences demonstrated that adjacent to the distal breakpoint lies a transposable element that we called Odysseus. Characteristics of Odysseus' terminal region and its cytological distribution in different strains as well as within the same strain indicate that Odysseus is an actively transposing element. The presence of Odysseus at the junction of the naturally occurring inversion In(2R)d' suggests that the inversion may be the result of the TEs activity. Cytological evidence from Drosophila melanogaster has also implicated the hobo transposable element in the generation of certain Hawaiian endemic inversions. This picture supports the hypothesis of the important role of TEs in generating natural inversions.     

Principles of genome evolution in the Drosophila melanogaster species group

That closely related species often differ by chromosomal inversions was discovered by Sturtevant and Plunkett in 1926. Our knowledge of how these inversions originate is still very limited, although a prevailing view is that they are facilitated by ectopic recombination events between inverted repetitive sequences. The availability of genome sequences of related species now allows us to study in detail the mechanisms that generate interspecific inversions. We have analyzed the breakpoint regions of the 29 inversions that differentiate the chromosomes of Drosophila melanogaster and two closely related species, D. simulans and D. yakuba, and reconstructed the molecular events that underlie their origin. Experimental and computational analysis revealed that the breakpoint regions of 59% of the inversions (17/29) are associated with inverted duplications of genes or other nonrepetitive sequences. In only two cases do we find evidence for inverted repetitive sequences in inversion breakpoints. We propose that the presence of inverted duplications associated with inversion breakpoint regions is the result of staggered breaks, either isochromatid or chromatid, and that this, rather than ectopic exchange between inverted repetitive sequences, is the prevalent mechanism for the generation of inversions in the melanogaster species group. Outgroup analysis also revealed evidence for widespread breakpoint recycling. Lastly, we have found that expression domains in D. melanogaster may be disrupted in D. yakuba, bringing into question their potential adaptive significance.     

Principles of Genome Evolution in the Drosophila melanogaster Species Group

That closely related species often differ by chromosomal inversions was discovered by Sturtevant and Plunkett in 1926. Our knowledge of how these inversions originate is still very limited, although a prevailing view is that they are facilitated by ectopic recombination events between inverted repetitive sequences. The availability of genome sequences of related species now allows us to study in detail the mechanisms that generate interspecific inversions. We have analyzed the breakpoint regions of the 29 inversions that differentiate the chromosomes of Drosophila melanogaster and two closely related species, D. simulans and D. yakuba, and reconstructed the molecular events that underlie their origin. Experimental and computational analysis revealed that the breakpoint regions of 59% of the inversions (17/29) are associated with inverted duplications of genes or other nonrepetitive sequences. In only two cases do we find evidence for inverted repetitive sequences in inversion breakpoints. We propose that the presence of inverted duplications associated with inversion breakpoint regions is the result of staggered breaks, either isochromatid or chromatid, and that this, rather than ectopic exchange between inverted repetitive sequences, is the prevalent mechanism for the generation of inversions in the melanogaster species group. Outgroup analysis also revealed evidence for widespread breakpoint recycling. Lastly, we have found that expression domains in D. melanogaster may be disrupted in D. yakuba, bringing into question their potential adaptive significance.     

Variation of Spontaneous Occurrence Rates of Chromosomal Aberrations in the Second Chromosomes of DROSOPHILA MELANOGASTER

After accumulating mutations by the aid of marked inversions, spontaneous occurrence rates of chromosome aberrations were estimated for 1148 chromosome lines that originated from five stem line second chromosomes of Drosophila melanogaster. In chromosome lines originating from three stem chromosomes (CH, PQ, and RT), mutations were accumulated for 7550, 7252, and 7256 chromosome generations, respectively, but no structural change was detected. For the chromosome lines that originated from the other two stem chromosomes, the situation was different: Twenty aberrations (19 paracentric inversions and 1 translocation between the second and the third chromosomes) during 45990 chromosome generations took place in the 500 chromosome lines derived from stem line chromosome (AW), and 92 aberrations (83 paracentric inversions, 6 pericentric inversions, 2 translocations between the second and the third chromosomes and 1 transposition) arose during 45006 chromosome generations in the 500 chromosome lines derived from stem line chromosome (JH). For the AW group the occurrence rate becomes 0.00043 per chromosome per generation for all aberrations and 0.00041 for inversions. For the JH group the corresponding rates are 0.00204 and 0.00198, respectively.-A non-random distribution of the breakpoint on the salivary gland chromosome was observed and the breakpoints were concentrated in the regions 26, 29, 33, and 34.-The cytoplasms and the chromosomes (other than the second chromosomes) were made approximately uniform throughout the experiments. Thus, this remarkable variability in the occurrence rate is most probably due to the differences in one or more chromosomal elements on the original five stem chromosomes. The mutable chromosomes (AW and JH) appear to carry a kind of mutator factor such as hi (Ives 1950).     

INVERSION LENGTH AND BREAKPOINT DISTRIBUTION IN THE DROSOPHILA BUZZATII SPECIES COMPLEX: IS INVERSION LENGTH A SELECTED TRAIT?

Length and position of breakpoints are characteristics of inversions that can be precisely determined on the polytene chromosomes of Drosophila species, and they provide crucial information about the processes that govern the origin and evolution of inversions. Eighty-six paracentric inversions described in the Drosophila buzzatii species complex and 18 inversions induced by introgressive hybridization in D. buzzatii were analyzed. In contrast to previous studies, inversion length and breakpoint distribution have been considered simultaneously. We conclude that: (1) inversion length is a selected trait; rare inversions are predominantly small while evolutionarily successful inversions, polymorphic and fixed, are predominantly intermediate in length; a nearly continuous variation in length, from small to medium sized, is found between less and more successful inversions; (2) there exists a significant negative correlation between length and number of polymorphic inversions per species which explains 39% of the inversion length variance; (3) natural selection on inversion length seems the main factor determining the relative position of breakpoints along the chromosomes; (4) the distribution of breakpoints according to their band location is non-random, with chromosomal segments that accumulate up to eight breakpoints.     

Contrasting Histories of Three Gene Regions Associated with In(3l)payne of Drosophila Melanogaster

In the present report, we studied nucleotide variation in three gene regions of Drosophila melanogaster, spanning >5 kb and showing different degrees of association with the cosmopolitan inversion In(3-L)Payne. The analysis of sequence variation in the regions surrounding the breakpoints and the heat shock 83 (Hsp83) gene locus, located close to the distal breakpoint, revealed the absence of shared polymorphisms and the presence of a number of fixed differences between arrangements, indicating absence of genetic exchange. In contrast, for the esterase-6 gene region, located in the center of the inversion, we observed the presence of shared polymorphisms between arrangements suggesting genetic exchange. In the regions close to the breakpoints, the common St arrangement is 10 times more polymorphic than inverted chromosomes. We propose that the lack of recombination between arrangements in these regions coupled with genetic hitchhiking is the best explanation for the low heterozygosity observed in inverted lines. Using the data for the breakpoints, we estimate that this inversion polymorphism is around 0.36 million yr old. Although it is widely accepted that inversions are examples of balanced polymorphisms, none of the current neutrality tests including our Monte Carlo simulations showed significant departure from neutral expectations.     

Molecular population genetics of the rp49 gene region in different chromosomal inversions of Drosophila subobscura.

Nucleotide variation at the ribosomal protein 49 (rp49) gene region has been studied in 75 lines of Drosophila subobscura belonging to four chromosomal arrangements (Ost, O3+4, O3+4+8, and O3+4+23). The location of the rp49 gene region within the inversion loop differs among heterokaryotypes: it is very close to one of the breakpoints in heterozygotes involving Ost chromosomes, while it is in a more central position in all other heterokaryotypes. The distribution of nucleotide polymorphism in the different arrangements is consistent with a monophyletic origin of the inversions. The data also provide evidence that gene conversion and possibly double crossover are involved in shuffling nucleotide variation among gene arrangements. The analyses reveal that the level of genetic exchange is higher when the region is located in a more central position of the inverted fragment than when it is close to the breakpoints. The pairwise difference distributions as well as the negative values of Tajima''s and Fu and Li''s statistics further support the hypothesis that nucleotide variation within chromosomal arrangements still reflects expansion after the origin of the inversions. Under the expansion model, we have estimated the time of origin of the studied inversions.     

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.     

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.     

Investigation of the origins of human autosomal inversions

A significant proportion of both pericentric and paracentric inversions have recurrent breakpoints and so could either have arisen through multiple independent events or be identical by descent (IBD) with a single common ancestor. Of two common variant inversions previously studied, the inv(2)(p11q13) was genuinely recurrent while the inv(10)(p11.2q21.2) was IBD in all cases tested. Excluding these two variants we have ascertained 257 autosomal inversion probands at the Wessex Regional Genetics Laboratory. There were 104 apparently recurrent inversions, representing 35 different breakpoint combinations and we speculated that at least some of these had arisen on more than one occasion. However, haplotype analysis identified no recurrent cases among eight inversions tested, including the variant inv(5)(p13q13). The cases not IBD were shown to have different breakpoints at the molecular cytogenetic level. No crossing over was detected within any of the inversions and the founder haplotypes extended for variable distances beyond the inversion breakpoints. Defining breakpoint intervals by FISH mapping identified no obvious predisposing elements in the DNA sequence. In summary the vast majority of human inversions arise as unique events. Even apparently recurrent inversions, with the exception of the inv(2)(p12q13), are likely to be either derived from a common ancestor or to have subtly different breakpoints. Presumably the lack of selection against most inversions allows them to accumulate and disperse amongst different populations over time. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.