Is mammalian chromosomal evolution driven by regions of genome fragility?

Background

A fundamental question in comparative genomics concerns the identification of mechanisms that underpin chromosomal change. In an attempt to shed light on the dynamics of mammalian genome evolution, we analyzed the distribution of syntenic blocks, evolutionary breakpoint regions, and evolutionary breakpoints taken from public databases available for seven eutherian species (mouse, rat, cattle, dog, pig, cat, and horse) and the chicken, and examined these for correspondence with human fragile sites and tandem repeats.

Results

Our results confirm previous investigations that showed the presence of chromosomal regions in the human genome that have been repeatedly used as illustrated by a high breakpoint accumulation in certain chromosomes and chromosomal bands. We show, however, that there is a striking correspondence between fragile site location, the positions of evolutionary breakpoints, and the distribution of tandem repeats throughout the human genome, which similarly reflect a non-uniform pattern of occurrence.

Conclusion

These observations provide further evidence that certain chromosomal regions in the human genome have been repeatedly used in the evolutionary process. As a consequence, the genome is a composite of fragile regions prone to reorganization that have been conserved in different lineages, and genomic tracts that do not exhibit the same levels of evolutionary plasticity.     

Is mammalian chromosomal evolution driven by regions of genome fragility?

Background  

A fundamental question in comparative genomics concerns the identification of mechanisms that underpin chromosomal change. In an attempt to shed light on the dynamics of mammalian genome evolution, we analyzed the distribution of syntenic blocks, evolutionary breakpoint regions, and evolutionary breakpoints taken from public databases available for seven eutherian species (mouse, rat, cattle, dog, pig, cat, and horse) and the chicken, and examined these for correspondence with human fragile sites and tandem repeats.     

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.     

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.     

A high-resolution cat radiation hybrid and integrated FISH mapping resource for phylogenomic studies across Felidae

We describe the construction of a high-resolution radiation hybrid (RH) map of the domestic cat genome, which includes 2662 markers, translating to an estimated average intermarker distance of 939 kilobases (kb). Targeted marker selection utilized the recent feline 1.9x genome assembly, concentrating on regions of low marker density on feline autosomes and the X chromosome, in addition to regions flanking interspecies chromosomal breakpoints. Average gap (breakpoint) size between cat-human ordered conserved segments is less than 900 kb. The map was used for a fine-scale comparison of conserved syntenic blocks with the human and canine genomes. Corroborative fluorescence in situ hybridization (FISH) data were generated using 129 domestic cat BAC clones as probes, providing independent confirmation of the long-range correctness of the map. Cross-species hybridization of BAC probes on divergent felids from the genera Profelis (serval) and Panthera (snow leopard) provides further evidence for karyotypic conservation within felids, and demonstrates the utility of such probes for future studies of chromosome evolution within the cat family and in related carnivores. The integrated map constitutes a comprehensive framework for identifying genes controlling feline phenotypes of interest, and to aid in assembly of a higher coverage feline genome sequence.     

Intrachromosomal rearrangements in avian genome evolution: evidence for regions prone to breakpoints

It is generally believed that the organization of avian genomes remains highly conserved in evolution as chromosome number is constant and comparative chromosome painting demonstrated there to be very few interchromosomal rearrangements. The recent sequencing of the zebra finch (Taeniopygia guttata) genome allowed an assessment of the number of intrachromosomal rearrangements between it and the chicken (Gallus gallus) genome, revealing a surprisingly high number of intrachromosomal rearrangements. With the publication of the turkey (Meleagris gallopavo) genome it has become possible to describe intrachromosomal rearrangements between these three important avian species, gain insight into the direction of evolutionary change and assess whether breakpoint regions are reused in birds. To this end, we aligned entire chromosomes between chicken, turkey and zebra finch, identifying syntenic blocks of at least 250 kb. Potential optimal pathways of rearrangements between each of the three genomes were determined, as was a potential Galliform ancestral organization. From this, our data suggest that around one-third of chromosomal breakpoint regions may recur during avian evolution, with 10% of breakpoints apparently recurring in different lineages. This agrees with our previous hypothesis that mechanisms of genome evolution are driven by hotspots of non-allelic homologous recombination.     

Analysis of Micro-Rearrangements in 25 Eukaryotic Species Pairs by SyntenyMapper

High-quality mapping of genomic regions and genes between two organisms is an indispensable prerequisite for evolutionary analyses and comparative genomics. Existing approaches to this problem focus on either delineating orthologs or finding extended sequence regions of common evolutionary origin (syntenic blocks). We propose SyntenyMapper, a novel tool for refining predefined syntenic regions. SyntenyMapper creates a set of blocks with conserved gene order between two genomes and finds all minor rearrangements that occurred since the evolutionary split of the two species considered. We also present TrackMapper, a SyntenyMapper-based tool that allows users to directly compare genome features, such as histone modifications, between two organisms, and identify genes with highly conserved features. We demonstrate SyntenyMapper''s advantages by conducting a large-scale analysis of micro-rearrangements within syntenic regions of 25 eukaryotic species. Unsurprisingly, the number and length of syntenic regions is correlated with evolutionary distance, while the number of micro-rearrangements depends only on the size of the harboring region. On the other hand, the size of rearranged regions remains relatively constant regardless of the evolutionary distance between the organisms, implying a length constraint in the rearrangement process. SyntenyMapper is a useful software tool for both large-scale and gene-centric genome comparisons.     

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.     

Assessing the role of tandem repeats in shaping the genomic architecture of great apes

Background

Ancestral reconstructions of mammalian genomes have revealed that evolutionary breakpoint regions are clustered in regions that are more prone to break and reorganize. What is still unclear to evolutionary biologists is whether these regions are physically unstable due solely to sequence composition and/or genome organization, or do they represent genomic areas where the selection against breakpoints is minimal.

Methodology and Principal Findings

Here we present a comprehensive study of the distribution of tandem repeats in great apes. We analyzed the distribution of tandem repeats in relation to the localization of evolutionary breakpoint regions in the human, chimpanzee, orangutan and macaque genomes. We observed an accumulation of tandem repeats in the genomic regions implicated in chromosomal reorganizations. In the case of the human genome our analyses revealed that evolutionary breakpoint regions contained more base pairs implicated in tandem repeats compared to synteny blocks, being the AAAT motif the most frequently involved in evolutionary regions. We found that those AAAT repeats located in evolutionary regions were preferentially associated with Alu elements.

Significance

Our observations provide evidence for the role of tandem repeats in shaping mammalian genome architecture. We hypothesize that an accumulation of specific tandem repeats in evolutionary regions can promote genome instability by altering the state of the chromatin conformation or by promoting the insertion of transposable elements.     

progressiveMauve: Multiple Genome Alignment with Gene Gain,Loss and Rearrangement

Background

Multiple genome alignment remains a challenging problem. Effects of recombination including rearrangement, segmental duplication, gain, and loss can create a mosaic pattern of homology even among closely related organisms.

Methodology/Principal Findings

We describe a new method to align two or more genomes that have undergone rearrangements due to recombination and substantial amounts of segmental gain and loss (flux). We demonstrate that the new method can accurately align regions conserved in some, but not all, of the genomes, an important case not handled by our previous work. The method uses a novel alignment objective score called a sum-of-pairs breakpoint score, which facilitates accurate detection of rearrangement breakpoints when genomes have unequal gene content. We also apply a probabilistic alignment filtering method to remove erroneous alignments of unrelated sequences, which are commonly observed in other genome alignment methods. We describe new metrics for quantifying genome alignment accuracy which measure the quality of rearrangement breakpoint predictions and indel predictions. The new genome alignment algorithm demonstrates high accuracy in situations where genomes have undergone biologically feasible amounts of genome rearrangement, segmental gain and loss. We apply the new algorithm to a set of 23 genomes from the genera Escherichia, Shigella, and Salmonella. Analysis of whole-genome multiple alignments allows us to extend the previously defined concepts of core- and pan-genomes to include not only annotated genes, but also non-coding regions with potential regulatory roles. The 23 enterobacteria have an estimated core-genome of 2.46Mbp conserved among all taxa and a pan-genome of 15.2Mbp. We document substantial population-level variability among these organisms driven by segmental gain and loss. Interestingly, much variability lies in intergenic regions, suggesting that the Enterobacteriacae may exhibit regulatory divergence.

Conclusions

The multiple genome alignments generated by our software provide a platform for comparative genomic and population genomic studies. Free, open-source software implementing the described genome alignment approach is available from http://gel.ahabs.wisc.edu/mauve.     

Comparative cytogenetics of human chromosome 3q21.3 reveals a hot spot for ectopic recombination in hominoid evolution

Fluorescence in situ hybridization mapping of fully integrated human BAC clones to primate chromosomes, combined with precise breakpoint localization by PCR analysis of flow-sorted chromosomes, was used to analyze the evolutionary rearrangements of the human 3q21.3-syntenic region in orangutan, siamang gibbon, and silvered-leaf monkey. Three independent evolutionary breakpoints were localized within a 230-kb segment contained in BACs RP11-93K22 and RP11-77P16. Approximately 200 kb of the human 3q21.3 sequence was not present on the homologous orangutan, siamang, and Old World monkey chromosomes, suggesting a genomic DNA insertion into the breakpoint region in the lineage leading to humans and African great apes. The breakpoints in the orangutan and siamang genomes were narrowed down to 12- and 20-kb DNA segments, respectively, which are enriched with endogenous retrovirus long terminal repeats and other repetitive elements. The inserted DNA segment represents part of an ancestral duplication. Paralogous sequence blocks were found at human 3q21, approximately 4 Mb proximal to the evolutionary breakpoint cluster region; at human 3p12.3, which contains an independent orangutan-specific breakpoint; and at the subtelomeric and pericentromeric regions of multiple human and orangutan chromosomes. The evolutionary breakpoint regions between human chromosome 3 and orangutan 2 as well their paralogous segments in the human genome coincide with breaks of chromosomal synteny in the mouse, rat, and/or chicken genomes. Collectively our data reveal reuse of the same short recombinogenic DNA segments in primate and vertebrate evolution, supporting a nonrandom breakage model of genome evolution.     

Precise detection of rearrangement breakpoints in mammalian chromosomes

Background  

Genomes undergo large structural changes that alter their organisation. The chromosomal regions affected by these rearrangements are called breakpoints, while those which have not been rearranged are called synteny blocks. We developed a method to precisely delimit rearrangement breakpoints on a genome by comparison with the genome of a related species. Contrary to current methods which search for synteny blocks and simply return what remains in the genome as breakpoints, we propose to go further and to investigate the breakpoints themselves in order to refine them.     

Synorth: exploring the evolution of synteny and long-range regulatory interactions in vertebrate genomes

Genomic regulatory blocks are chromosomal regions spanned by long clusters of highly conserved noncoding elements devoted to long-range regulation of developmental genes, often immobilizing other, unrelated genes into long-lasting syntenic arrangements. Synorth is a web resource for exploring and categorizing the syntenic relationships in genomic regulatory blocks across multiple genomes, tracing their evolutionary fate after teleost whole genome duplication at the level of genomic regulatory block loci, individual genes, and their phylogenetic context.     

Comparative mapping of human chromosome 10 to pig chromosomes 10 and 14

Identification of predictive markers in QTL regions that impact production traits in commercial populations of swine is dependent on construction of dense comparative maps with human and mouse genomes. Chromosomal painting in swine suggests that large genomic blocks are conserved between pig and human, while mapping of individual genes reveals that gene order can be quite divergent. High-resolution comparative maps in regions affecting traits of interest are necessary for selection of positional candidate genes to evaluate nucleotide variation causing phenotypic differences. The objective of this study was to construct an ordered comparative map of human chromosome 10 and pig chromosomes 10 and 14. As a large portion of both pig chromosomes are represented by HSA10, genes at regularly spaced intervals along this chromosome were targeted for placement in the porcine genome. A total of 29 genes from human chromosome 10 were mapped to porcine chromosomes 10 (SSC10) and 14 (SSC14) averaging about 5 Mb distance of human DNA per marker. Eighteen genes were assigned by linkage in the MARC mapping population, five genes were physically assigned with the IMpRH mapping panel and seven genes were assigned on both maps. Seventeen genes from human 10p mapped to SSC10, and 12 genes from human 10q mapped to SSC14. Comparative maps of mammalian species indicate that chromosomal segments are conserved across several species and represent syntenic blocks with distinct breakpoints. Development of comparative maps containing several species should reveal conserved syntenic blocks that will allow us to better define QTL regions in livestock.     

Chromosomal inversions associated with environmental adaptation in honeybees

Chromosomal inversions can facilitate local adaptation in the presence of gene flow by suppressing recombination between well‐adapted native haplotypes and poorly adapted migrant haplotypes. East African mountain populations of the honeybee Apis mellifera are highly divergent from neighbouring lowland populations at two extended regions in the genome, despite high similarity in the rest of the genome, suggesting that these genomic regions harbour inversions governing local adaptation. Here, we utilize a new highly contiguous assembly of the honeybee genome to characterize these regions. Using whole‐genome sequencing data from 55 highland and lowland bees, we find that the highland haplotypes at both regions are present at high frequencies in three independent highland populations but extremely rare elsewhere. The boundaries of both divergent regions are characterized by regions of high homology with each other positioned in opposite orientations and contain highly repetitive, long inverted repeats with homology to transposable elements. These regions are likely to represent inversion breakpoints that participate in nonallelic homologous recombination. Using long‐read data, we confirm that the lowland samples are contiguous across breakpoint regions. We do not find evidence for disruption of functional sequence by these breakpoints, which suggests that the inversions are likely maintained due to their allelic content conferring local adaptation in highland environments. Finally, we identify a third divergent genomic region, which contains highly divergent segregating haplotypes that also may contain inversion variants under selection. The results add to a growing body of evidence indicating the importance of chromosomal inversions in local adaptation.     

Comparative organization of wheat homoeologous group 3S and 7L using wheat-rice synteny and identification of potential markers for genes controlling xanthophyll content in wheat

EST and genomic DNA sequencing efforts for rice and wheat have provided the basis for interpreting genome organization and evolution. In this study we have used EST and genomic sequencing information and a bioinformatic approach in a two-step strategy to align portions of the wheat and rice genomes. In the first step, wheat ESTs were used to identify rice orthologs and it was shown that wheat 3S and rice 1 contain syntenic units with intrachromosomal rearrangements. Further analysis using anchored rice contiguous sequences and TBLASTX alignments in a second alignment step showed interruptions by orthologous genes that map elsewhere in the wheat genome. This indicates that gene content and order is not as conserved as large chromosomal blocks as previously predicted. Similarly, chromosome 7L contains syntenic units with rice 6 and 8 but is interrupted by combinations of intrachromosomal and interchromosomal rearrangements involving syntenic units and single gene orthologs from other rice chromosome groups. We have used the rice sequence annotations to identify genes that can be used to develop markers linked to biosynthetic pathways on 3BS controlling xanthophyll production in wheat and thus involved in determining flour colour.Electronic Supplementary Material Supplementary material is available in the online version of this article at .     

Chromosomal rearrangements breakpoints clusterization: is the clonal selection involved

Certain types of cancer are often correlate with certain chromosomal rearrangements. The chimaeric genes are formed as a result of this rearrangements, and the chimaeric proteins are the products of their expression. The breakpoints of such translocation are often clustered in the genome. Moreover, such breakpoint clusters often contain specific genomic elements like topoisomerase II consensus sites, nuclear matrix attachment regions and DNA sequences, which can make up secondary non-canonical structure. In this review we discuss whether breakpoints may be induced by chromatin structure. Furthermore, we bring up not touched upon literature question about the relation between the breakpoint clusters and the domain organization of corresponding proteins. We also consider possible mechanisms of chromosomal rearrangements.