Inversion Breakpoints Research Articles

Characterization of the chromosomal inversion associated with the Koa mutation in the mouse revealed the cause of skeletal abnormalities

BackgroundKoala (Koa) is a dominant mutation in mice causing bushy muzzle and pinna, and is associated with a chromosomal inversion on the distal half of chromosome 15. To identify the gene responsible for the Koa phenotypes, we investigated phenotypes of Koa homozygous mice and determined the breakpoints of the inversion with a genetic method using recombination between two different chromosomal inversions.ResultsSkeletal preparation of Koa homozygotes showed marked deformity of the ribs and a wider skull with extended zygomatic arches, in addition to a general reduction in the lengths of long bones. They also had open eyelids at birth caused by a defect in the extension of eyelid anlagen during the embryonic stages. The proximal and distal breakpoints of the Koa inversion were determined to be 0.8-Mb distal to the Trsps1 gene and to 0.1-Mb distal to the Hoxc4 gene, respectively, as previously reported. The phenotypes of mice with the recombinant inverted chromosomes revealed the localization of the gene responsible the Koa phenotype in the vicinity of the proximal recombinant breakpoint. Expression of the Trsps1 gene in this region was significantly reduced in the Koa homozygous and heterozygous embryos.ConclusionWhile no gene was disrupted by the chromosomal inversion, an association between the Koa phenotype and the proximal recombinant breakpoint, phenotypic similarities with Trps1-deficient mice or human patients with TRSP1 mutations, and the reduced expression of the Trsps1 gene in Koa mice, indicated that the phenotypes of the Koa mice are caused by the altered expression of the Trps1 gene.

X-Linked Congenital Adrenal Hypoplasia with Hypogonadotropic Hypogonadism Caused by an Inversion Disrupting a Conserved Noncoding Element Upstream of theNR0B1(DAX1) Gene

X-linked congenital adrenal hypoplasia with hypogonadotropic hypogonadism (AHCH) is known to be caused by coding mutations in the nuclear receptor subfamily 0, group B, member 1 (NR0B1) gene, encoding the transcriptional repressor dosage-sensitive sex-reversal adrenal hypoplasia critical region on the X chromosome protein 1 (DAX1). Four males in a family were affected by AHCH. Our aim was to locate the genetic cause of their disease, knowing that they had no mutation in the obvious candidate gene, NR0B1. Linkage analysis of the X chromosome and mutational screening of conserved noncoding regions upstream of NR0B1 were performed. To functionally characterize the genetic defect, studies of transcription and expression of DAX1 and steroidogenic factor 1 (SF-1) were done. A 60 Mb inversion on the X chromosome with one of the inversion breakpoints located in a conserved noncoding region 4 kb upstream of NR0B1 was detected. The inversion causes relocation of a putative SF-1 binding site implicated in murine gonadal development. A reporter construct lacking this enhancer element upstream of NR0B1 was unresponsive to SF-1 transcriptional activation. Immunohistochemistry suggested that the inversion leads to SF-1 silencing in the patients' testes both in childhood and in adult life. We report a noncoding mutation causing AHCH, an inversion resulting in a phenotype similar to what is caused by intragenic NR0B1 null mutations. The inversion seems to disrupt and/or relocate regulatory sites crucial in DAX1 expression.

Disruption of CNTNAP2 and additional structural genome changes in a boy with speech delay and autism spectrum disorder

Patients with autism spectrum disorder (ASD) frequently harbour chromosome rearrangements and segmental aneuploidies, which allow us to identify candidate genes. In a boy with mild facial dysmorphisms, speech delay and ASD, we reconstructed by karyotyping, FISH and SNP array-based segmental aneuploidy profiling a highly complex chromosomal rearrangement involving at least three breaks in chromosome 1 and seven breaks in chromosome 7. Chromosome banding revealed an inversion of region 7q32.1-7q35 on the derivative chromosome 7. FISH with region-specific BACs mapped both inversion breakpoints and revealed additional breaks and structural changes in the CNTNAP2 gene. Two gene segments were transposed and inserted into the 1q31.2 region, while the CNTNAP2 segment between the two transposed parts as well as intron 13 to the 5-UTR were retained on the der(7). SNP array analysis revealed an additional de novo deletion encompassing the distal part of intron1 and exon 2 of CNTNAP2, which contains FOXP2 binding sites. Second, we found another de novo deletion on chromosome 1q41, containing 15 annotated genes, including KCTD3 and USH2A. Disruptions of the CNTNAP2 gene have been associated with ASD and with Gilles de la Tourette syndrome (GTS). Comparison of disruptions of CNTNAP2 in patients with GTS and ASD suggests that large proximal disruptions result in either GTS or ASD, while relatively small distal disruptions may be phenotypically neutral. For full-blown ASD to develop, a proximal disruption of CNTNAP2 may have to occur concomitantly with additional genome mutations such as hemizygous deletions of the KCTD3 and USH2A genes.

The Genomics of Speciation in Drosophila: Diversity, Divergence, and Introgression Estimated Using Low-Coverage Genome Sequencing

In nature, closely related species may hybridize while still retaining their distinctive identities. Chromosomal regions that experience reduced recombination in hybrids, such as within inversions, have been hypothesized to contribute to the maintenance of species integrity. Here, we examine genomic sequences from closely related fruit fly taxa of the Drosophila pseudoobscura subgroup to reconstruct their evolutionary histories and past patterns of genic exchange. Partial genomic assemblies were generated from two subspecies of Drosophila pseudoobscura (D. ps.) and an outgroup species, D. miranda. These new assemblies were compared to available assemblies of D. ps. pseudoobscura and D. persimilis, two species with overlapping ranges in western North America. Within inverted regions, nucleotide divergence among each pair of the three species is comparable, whereas divergence between D. ps. pseudoobscura and D. persimilis in non-inverted regions is much lower and closer to levels of intraspecific variation. Using molecular markers flanking each of the major chromosomal inversions, we identify strong crossover suppression in F1 hybrids extending over 2 megabase pairs (Mbp) beyond the inversion breakpoints. These regions of crossover suppression also exhibit the high nucleotide divergence associated with inverted regions. Finally, by comparison to a geographically isolated subspecies, D. ps. bogotana, our results suggest that autosomal gene exchange between the North American species, D. ps. pseudoobscura and D. persimilis, occurred since the split of the subspecies, likely within the last 200,000 years. We conclude that chromosomal rearrangements have been vital to the ongoing persistence of these species despite recent hybridization. Our study serves as a proof-of-principle on how whole genome sequencing can be applied to formulate and test hypotheses about species formation in lesser-known non-model systems.

Cruciform-forming inverted repeats appear to have mediated many of the microinversions that distinguish the human and chimpanzee genomes

Submicroscopic inversions have contributed significantly to the genomic divergence between humans and chimpanzees over evolutionary time. Those microinversions which are flanked by segmental duplications (SDs) are presumed to have originated via non-allelic homologous recombination between SDs arranged in inverted orientation. However, the nature of the mechanisms underlying those inversions which are not flanked by SDs remains unclear. We have investigated 35 such inversions, ranging in size from 51-nt to 22056-nt, with the goal of characterizing the DNA sequences in the breakpoint-flanking regions. Using the macaque genome as an outgroup, we determined the lineage specificity of these inversions and noted that the majority (N = 31; 89%) were associated with deletions (of length between 1-nt and 6754-nt) immediately adjacent to one or both inversion breakpoints. Overrepresentations of both direct and inverted repeats, >or= 6-nt in length and capable of non-B DNA structure formation, were noted in the vicinity of breakpoint junctions suggesting that these repeats could have contributed to double strand breakage. Inverted repeats capable of cruciform structure formation were also found to be a common feature of the inversion breakpoint-flanking regions, consistent with these inversions having originated through the resolution of Holliday junction-like cruciforms. Sequences capable of non-B DNA structure formation have previously been implicated in promoting gross deletions and translocations causing human genetic disease. We conclude that non-B DNA forming sequences may also have promoted the occurrence of mutations in an evolutionary context, giving rise to at least some of the inversion/deletions which now serve to distinguish the human and chimpanzee genomes.

Split hand/foot malformation due to chromosome 7q aberrations(SHFM1): additional support for functional haploinsufficiency as the causative mechanism

We report on three patients with split hand/foot malformation type 1 (SHFM1). We detected a deletion in two patients and an inversion in the third, all involving chromosome 7q21q22. We performed conventional chromosomal analysis, array comparative genomic hybridization and fluorescence in situ hybridization. Both deletions included the known genes associated with SHFM1 (DLX5, DLX6 and DSS1), whereas in the third patient one of the inversion break points was located just centromeric to these genes. These observations confirm that haploinsufficiency due to either a simultaneous deletion of these genes or combined downregulation of gene expression due to a disruption in the region between these genes and a control element could be the cause of the syndrome. We review previously reported studies that support this hypothetical mechanism.

Repeat mediated gene duplication in the Drosophila pseudoobscura genome

Genetic mutations can occur on a wide variety of scales, including those that change single nucleotides, those that add or remove content to/from a genome, and those that change the organization of a genome. Gene duplications are a specific class of mutations that add content to a genome, and they can arise via a wide variety of mechanisms. I examined the mechanisms responsible for recently duplicated genes in the D. pseudoobscura genome, and I observed both retroposed and DNA duplications. Many duplicated genes lack signatures of either retroposition or DNA-based mechanisms, but other features of these ambiguously duplicated genes suggest that most were generated via retroposition. Furthermore, close examination of sequences flanking DNA duplications and those found at the breakpoints of chromosomal inversions suggests a connection between these two events. In Drosophila, duplicated genes near inversion breakpoints can arise via unequal genetic exchange during the non-allelic crossing over event giving rise to the inversion. I observed one duplicated gene in the D. pseudoobscura genome that appears to have been generated by this mechanism. Additionally, many DNA duplications in the D. pseudoobscura genome are flanked by a repetitive sequence also found at the breakpoints of chromosomal inversions. This suggests that the molecular mechanisms responsible for chromosomal rearrangements and some duplicated genes have overlapping processes.

A chromosome 10 variant with a 12 Mb inversion [inv(10)(q11.22q21.1)] identical by descent and frequent in the Swedish population

We identified a paracentric inversion of chromosome 10 [inv(10)(q11.22q21.1)] in 0.20% of Swedish individuals (15/7,439) referred for cytogenetic analysis. A retrospective analysis of 8,896 karyotypes from amniocenteses in Sweden revealed a carrier frequency of 0.079% (7/8,896) for the inversion. Cloning and detailed analysis of the inversion breakpoint regions show enrichment for interspersed repeat elements and AT-stretches. The centromeric breakpoint coincides with that of a predicted inversion from HapMap data, which suggests that this region is involved in several chromosome 10 variants. No known gene or predicted transcript are disrupted by the inversion which spans approximately 12 Mb. Carriers from four non-related Swedish families have identical inversion breakpoints and haplotype analysis confirmed that the rearrangement is identical by descent. Diagnosis was retrieved in 6 out of the 15 carriers referred for cytogenetic analysis. No consistent phenotype was found to be associated with the inversion. Our study demonstrates that the inv(10)(q11.22q21.1) is a rare and inherited chromosome variant with a broad geographical distribution in Sweden.

Identification and Characterization of a Novel Melanoma Tumor Suppressor Gene on Human Chromosome 6q21

By characterizing a complex chromosome rearrangement involving 6q and 17p in melanoma cell line UACC-930, we isolated a candidate tumor suppressor gene at 6q21, named prenyl diphosphate synthase subunit 2 (PDSS2), which was interrupted by an inversion breakpoint. The purpose of this study was to determine the tumor-suppressive potential of PDSS2 in the development of melanoma. To isolate the rearranged 6q in UACC-930 cells, a bacterial artificial chromosome clone (RP1-67A8) covering the breakpoint at 6q21 was digested with HindIII and each DNA fragment was used as the probe for the breakpoint in Southern blotting. The HindIII fragment probe covering the breakpoint was then used to screen an EcoRI-digested DNA library generated from UACC-930. To characterize the tumor-suppressive potential of PDSS2, PDSS2 was stably transfected into a highly tumorigenic melanoma cell line, UACC-903. The tumor-suppressive function of PDSS2 was shown by both in vitro and in vivo assays. The differential expression of PDSS2 in benign nevi and primary melanoma samples was also studied. Down-regulation of PDSS2 was observed in 59 of 87 (67.8%) primary melanomas, which was significantly higher than that in benign nevi (7 of 66, 10.6%; P < 0.001). In addition, an overexpression of the PDSS2 in UACC-903 cells could inhibit tumor cell growth, decrease the colony-forming ability in soft agar, and totally abrogate the tumorigenicity of UACC-903 in nude mice. Our results support the proposal that PDSS2 is a novel tumor suppressor gene that plays an important role in the development of malignant melanoma.

High-resolution mapping of heterochromatin redistribution in a Drosophila position-effect variegation model.

BackgroundPosition-effect variegation (PEV) is the stochastic transcriptional silencing of a gene positioned adjacent to heterochromatin. white-mottled X-chromosomal inversions in Drosophila are classic PEV models that show variegation of the eye color gene white due to its relocation next to pericentric heterochromatin. It has been suggested that in these models the spreading of heterochromatin across the rearrangement breakpoint causes the silencing of white. However, the extent of this spreading and the precise pattern of heterochromatin redistribution have remained unclear. To obtain insight into the mechanism of PEV, we constructed high-resolution binding maps of Heterochromatin Protein 1 (HP1) on white-mottled chromosomes.ResultsWe find that HP1 invades euchromatin across the inversion breakpoints over ~175 kb and ~30 kb, causing de novo association of HP1 with 20 genes. However, HP1 binding levels in these regions show substantial local variation, and white is the most strongly bound gene. Remarkably, white is also the only gene that is detectably repressed by heterochromatin. Furthermore, we find that HP1 binding to the invaded region is particularly sensitive to the dosage of the histone methyltransferase Su(var)3-9, indicating that the de novo formed heterochromatin is less stable than naturally occurring constitutive heterochromatin.ConclusionOur molecular maps demonstrate that heterochromatin can invade a normally euchromatic region, yet the strength of HP1 binding and effects on gene expression are highly dependent on local context. Our data suggest that the white gene has an unusual intrinsic affinity for heterochromatin, which may cause this gene to be more sensitive to PEV than most other genes.

Efficient sampling of parsimonious inversion histories with application to genome rearrangement in Yersinia.

Inversions are among the most common mutations acting on the order and orientation of genes in a genome, and polynomial-time algorithms exist to obtain a minimal length series of inversions that transform one genome arrangement to another. However, the minimum length series of inversions (the optimal sorting path) is often not unique as many such optimal sorting paths exist. If we assume that all optimal sorting paths are equally likely, then statistical inference on genome arrangement history must account for all such sorting paths and not just a single estimate. No deterministic polynomial algorithm is known to count the number of optimal sorting paths nor sample from the uniform distribution of optimal sorting paths.Here, we propose a stochastic method that uniformly samples the set of all optimal sorting paths. Our method uses a novel formulation of parallel Markov chain Monte Carlo. In practice, our method can quickly estimate the total number of optimal sorting paths. We introduce a variant of our approach in which short inversions are modeled to be more likely, and we show how the method can be used to estimate the distribution of inversion lengths and breakpoint usage in pathogenic Yersinia pestis.The proposed method has been implemented in a program called “MC4Inversion.” We draw comparison of MC4Inversion to the sampler implemented in BADGER and a previously described importance sampling (IS) technique. We find that on high-divergence data sets, MC4Inversion finds more optimal sorting paths per second than BADGER and the IS technique and simultaneously avoids bias inherent in the IS technique.

Chromosomal Inversions between Human and Chimpanzee Lineages Caused by Retrotransposons

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

The mouse hairy ears mutation exhibits an extended growth (anagen) phase in hair follicles and altered Hoxc gene expression in the ears

The mouse In(15)2Rl (hairy ears, Eh) mutation is a paracentric inversion of the distal half of chromosome 15 (Chr 15). Heterozygous Eh/+ mice display misshaped and hairy ears that have more and longer hair than the ears of their wild-type littermates. We mapped, cloned and sequenced both inversion breakpoints. No protein-coding transcript was disrupted by either breakpoint. The proximal breakpoint is located between syntrophin basic 1 (Sntb1) and hyaluronan synthase 2 (Has2), and the distal breakpoint maps between homeobox C4 (Hoxc4) and single-strand selective monofunctional uracil DNA glycosylase (Smug1), near the middle and the telomere ends of Chr 15, respectively. The inversion spans ~47 megabases. Our genetic analysis suggests that the hairy-ear phenotype is caused by the proximal breakpoint of the inversion-bearing Chr 15. Quantitative RNA analysis by real-time polymerase chain reaction for the genes flanking the breakpoint indicated no changes in expression levels except for some homeobox C (Hoxc) genes whose expression was elevated in developing and mature skin of the ears but not of other body regions. The increased hair length on the ears of Eh/+ mice was due to an extension of the anagen stage in the hair cycle, as determined by histological analysis. Our data indicate that the Eh phenotype arises from mis-expression of Hoxc genes.

A 2.1 Mb deletion adjacent but distal to a 14q21q23 paracentric inversion in a family with spherocytosis and severe learning difficulties

A familial q21.1q23.2-inversion on chromosome 14 that co-segregated with spherocytosis and learning difficulties or mild mental retardation was extensively investigated by bacterial artificial chromosome fluorescence in situ hybridization and array-comparative genomic hybridization. As expected, a deletion of the beta-spectrin gene SPTB, a known cause of spherocytosis, was found. More unexpectedly, this deletion was approximately 1.6 Mb distal to the 14q23.2-inversion breakpoint. The deletion spanned approximately 2.1 Mb and contained 15 annotated genes in addition to SPTB, among them PLEKHG3, a guanide nucleotide exchange factor for Rho GTPases. This gene is highly expressed in the brain and our best candidate for causing the mild mental retardation. The case illustrates that inversions can be associated with microdeletions close to but not including one of the inversion breakpoints.

Reconstructing Ancestral Autosomal Arrangements in theAnopheles gambiaeComplex

Members of the Anopheles gambiae complex have remarkably distinct ecological adaptations, behaviors, and degrees of vectorial capacity. Inferring phylogenetic relationships in the complex is crucial for identifying the genomic changes associated with the origin and loss of epidemiologically important traits. However, the high level of sequence similarity, genetic introgression, and shared molecular ancestral polymorphisms makes reconstruction of the A. gambiae complex phylogeny difficult. Phylogenetic relationships among the members of species complexes can be inferred from the distribution of fixed chromosomal inversions if outgroup arrangements are known. The aim of this work is to test a possibility of determining ancestral autosomal arrangements in the A. gambiae complex using outgroup chromosomes and a combination of bioinformatic and cytogenetic approaches. The minimum number of inversions between members of the A. gambiae complex and the outgroup species A. funestus and A. stephensi was calculated using the Multiple Genome Rearrangements (MGR) and Sorting Permutation by Reversals and block-INterchanGes (SPRING) programs. The physical mapping of A. merus chromosomes identified molecular coordinates of the proximal 2Ro+ inversion breakpoint in A. gambiae. DNA probes from 2La+ and 2Ro+ inversion breakpoints of the A. gambiae were mapped to the A. stephensi chromosomes. Assuming monophyletic origin of the inversions, this study concludes that physical mapping of ingroup and outgroup species can be used for identifying inversion breakpoints and ancestral autosomal arrangements within species complexes. Molecular characterization of the breakpoints in both ingroup and outgroup species will provide a solid basis for reconstructing the inversion history in the A. gambiae complex.

A position effect on TRPS1 is associated with Ambras syndrome in humans and the Koala phenotype in mice

Ambras syndrome (AS) is a rare form of congenital hypertrichosis with excessive hair on the shoulders, face and ears. Cytogenetic studies have previously implicated an association with rearrangements of chromosome 8. Here we define an 11.5 Mb candidate interval for AS on chromosome 8q based on cytogenetic breakpoints in three patients. TRPS1, a gene within this interval, was deleted in a patient with an 8q23 chromosomal rearrangement, while its expression was significantly downregulated in another patient with an inversion breakpoint 7.3 Mb downstream of TRPS1. Here, we describe the first potential long-range position effect on the expression of TRPS1. To gain insight into the mechanisms by which Trps1 affects the hair follicle, we performed a detailed analysis of the hair abnormalities in Koa mice, a mouse model of hypertrichosis. We found that the proximal breakpoint of the Koa inversion is located 791 kb upstream of Trps1. Quantitative real-time polymerase chain reaction, in situ hybridization and immunofluorescence analysis revealed that Trps1 expression levels are reduced in Koa mutant mice at the sites of pathology for the phenotype. We determined that the Koa inversion creates a new Sp1 binding site and translocates additional Sp1 binding sites within a highly conserved stretch spanning the proximal breakpoint, providing a potential mechanism for the position effect. Collectively, these results describe a position effect that downregulates TRPS1 expression as the probable cause of hypertrichosis in AS in humans and the Koa phenotype in mice.

An old bilbo-like non-LTR retroelement insertion provides insight into the relationship of species of the virilis group

In Drosophila, at the population and species level, fixation of a TE insertion is an unlikely fate. Of the few reported fixations at the species level most involve non-LTR retroelements. In this work we report the fixation of a non-LTR retroelement in five species ( Drosophila littoralis, Drosophila virilis, Drosophila lummei, Drosophila americana and Drosophila novamexicana) of the virilis group of Drosophila. In most species, this TE insertion is being lost through the accumulation of small deletions, but there is also evidence for the accumulation of large deletions. In the americana lineage an insertion of about 900 bp of the non-LTR retroelement is a marker for the Xc inversion. This insertion is, at most, 80 kb away from the basal Xc inversion breakpoint. The presence of a bilbo-like element in D. littoralis but not in D. kanekoi, suggests that D. littoralis is more closely related to species of the virilis phylad than to species of the montana phylad, which is in contrast with the traditional view. Nevertheless, the phylogenetic analyses here performed using a 7 gene dataset suggest that D. littoralis is indeed more closely related to species of the virilis phylad than to species of the montana phylad. The re-evaluation of the phylogenetic relationship of the species of the virilis group, under the assumption of a relaxed molecular clock, results in an estimated age of the bilbo-like element insertion of at least 7.5 Mya.

Evolutionary toggling of the MAPT 17q21.31 inversion region

Using comparative sequencing approaches, we investigated the evolutionary history of the European-enriched 17q21.31 MAPT inversion polymorphism. We present a detailed, BAC-based sequence assembly of the inverted human H2 haplotype and compare it to the sequence structure and genetic variation of the corresponding 1.5-Mb region for the noninverted H1 human haplotype and that of chimpanzee and orangutan. We found that inversion of the MAPT region is similarly polymorphic in other great ape species, and we present evidence that the inversions occurred independently in chimpanzees and humans. In humans, the inversion breakpoints correspond to core duplications with the LRRC37 gene family. Our analysis favors the H2 configuration and sequence haplotype as the likely great ape and human ancestral state, with inversion recurrences during primate evolution. We show that the H2 architecture has evolved more extensive sequence homology, perhaps explaining its tendency to undergo microdeletion associated with mental retardation in European populations.

Sequence signatures of a recent chromosomal rearrangement in Drosophila mojavensis

The X-chromosome inversion, Xe, distinguishes Drosophila mojavensis and D. arizonae. Earlier work mapped the breakpoints of this inversion to large intervals and provided hypotheses for the locations of the breakpoints within 3000-bp intergenic regions on the D. mojavensis genome sequence assembly. Here, we sequenced these regions directly in the putatively ancestral D. arizonae X-chromosome. We find that the two inversion breakpoints are near an inverted gene duplication and a common repetitive element, respectively, and these features were likely present in the non-inverted ancestral chromosome on the D. mojavensis lineage. Contrary to an earlier hypothesis, the inverted gene duplication appears to predate the inversion. We find no sequence similarity between the breakpoint regions in the D. mojavensis ancestor, excluding an ectopic-exchange model of chromosome rearrangements. We also found no evidence that staggered single-strand breaks caused the inversion. We suggest these features may have contributed to the chromosomal breakages resulting in this inversion.

Chromosomal Rearrangement Inferred From Comparisons of 12 Drosophila Genomes

The availability of 12 complete genomes of various species of genus Drosophila provides a unique opportunity to analyze genome-scale chromosomal rearrangements among a group of closely related species. This article reports on the comparison of gene order between these 12 species and on the fixed rearrangement events that disrupt gene order. Three major themes are addressed: the conservation of syntenic blocks across species, the disruption of syntenic blocks (via chromosomal inversion events) and its relationship to the phylogenetic distribution of these species, and the rate of rearrangement events over evolutionary time. Comparison of syntenic blocks across this large genomic data set confirms that genetic elements are largely (95%) localized to the same Muller element across genus Drosophila species and paracentric inversions serve as the dominant mechanism for shuffling the order of genes along a chromosome. Gene-order scrambling between species is in accordance with the estimated evolutionary distances between them and we find it to approximate a linear process over time (linear to exponential with alternate divergence time estimates). We find the distribution of synteny segment sizes to be biased by a large number of small segments with comparatively fewer large segments. Our results provide estimated chromosomal evolution rates across this set of species on the basis of whole-genome synteny analysis, which are found to be higher than those previously reported. Identification of conserved syntenic blocks across these genomes suggests a large number of conserved blocks with varying levels of embryonic expression correlation in Drosophila melanogaster. On the other hand, an analysis of the disruption of syntenic blocks between species allowed the identification of fixed inversion breakpoints and estimates of breakpoint reuse and lineage-specific breakpoint event segregation.