Is mammalian chromosomal evolution driven by regions of genome fragility?

Addressing chromosome evolution in the whole-genome sequence era

The evolution of karyotypes has been the subject of intensive study since the middle of the 20th century. This was motivated by the observation that the karyotypes of related species showed remarkable conservation. The recent emergence of whole-genome sequencing projects gives the opportunity to complement the cytogenetic approaches by addressing the conservation of karyotypes using chromosome sequence comparison. In this short review we present a description of recent advances in computational biology methods dedicated to the study of chromosome evolution and more specifically ancestral karyotype reconstruction in an attempt to provide an integrated overview of both cytogenetic and computational approaches.

Primate chromosome evolution: ancestral karyotypes, marker order and neocentromeres

In 1992 the Japanese macaque was the first species for which the homology of the entire karyotype was established by cross-species chromosome painting. Today, there are chromosome painting data on more than 50 species of primates. Although chromosome painting is a rapid and economical method for tracking translocations, it has limited utility for revealing intrachromosomal rearrangements. Fortunately, the use of BAC-FISH in the last few years has allowed remarkable progress in determining marker order along primate chromosomes and there are now marker order data on an array of primate species for a good number of chromosomes. These data reveal inversions, but also show that centromeres of many orthologous chromosomes are embedded in different genomic contexts. Even if the mechanisms of neocentromere formation and progression are just beginning to be understood, it is clear that these phenomena had a significant impact on shaping the primate genome and are fundamental to our understanding of genome evolution. In this report we complete and integrate the dataset of BAC-FISH marker order for human syntenies 1, 2, 4, 5, 8, 12, 17, 18, 19, 21, 22 and the X. These results allowed us to develop hypotheses about the content, marker order and centromere position in ancestral karyotypes at five major branching points on the primate evolutionary tree: ancestral primate, ancestral anthropoid, ancestral platyrrhine, ancestral catarrhine and ancestral hominoid. Current models suggest that between-species structural rearrangements are often intimately related to speciation. Comparative primate cytogenetics has become an important tool for elucidating the phylogeny and the taxonomy of primates. It has become increasingly apparent that molecular cytogenetic data in the future can be fruitfully combined with whole-genome assemblies to advance our understanding of primate genome evolution as well as the mechanisms and processes that have led to the origin of the human genome.

Chromosomal expression and localization of aphidicolin-induced fragile sites in the standard karyotype of river buffalo (Bubalus bubalis)

The present study reports on the chromosomal expression and localization of aphidicolin-induced fragile sites in the standard karyotype of river buffalo (Bubalus bubalis, 2n = 50) with the aim of establishing a 'fragile site map' of the species. Totally, 400 aphidicolin-induced breakages were analyzed from eight young and clinically healthy animals, four males and four females; these breakages were localized in 106 RBG-negative chromosome bands or at the band-interband regions. The number of breakages per chromosome did not vary statistically 'among' the animals investigated but the differences among individual chromosomes were highly significant thus indicating that the chromosomal distribution of the breakages is not random and appears only partially related to chromosome length. Fragile sites were statistically determined as those chromosomal bands showing three or more breakages. In the river buffalo karyotype, 51 fragile sites were detected and localized on the standardized ideogram of the species. The most fragile bands were as follows: 9q213 with 24 breakages out of 400; 19q21 with 16, 17q21 and inacXq24 with 15, 15q23 with 13 and 13q23 with 12 breaks, respectively. Previous gene mapping analysis in this species has revealed that the closest loci to these fragile sites contain genes such as RASA1 and CAST (9q214), NPR3 and C9 (19q19), PLP and BTK (Xq24-q25), OarCP09 (15q24), and EDNRB (13q22) whose mutations are responsible for severe phenotypic malformations and immunodeficiency in humans as well as in mice and meat quality in pigs. Further cytogenetic and molecular studies are needed to fully exploit the biological significance of the fragile sites in karyotype evolution of domestic animals and their relationships with productive and reproductive efficiency of livestock.

Genomic instability within centromeres of interspecific marsupial hybrids

Several lines of evidence suggest that, within a lineage, particular genomic regions are subject to instability that can lead to specific types of chromosome rearrangements important in species incompatibility. Within family Macropodidae (kangaroos, wallabies, bettongs, and potoroos), which exhibit recent and extensive karyotypic evolution, rearrangements involve chiefly the centromere. We propose that centromeres are the primary target for destabilization in cases of genomic instability, such as interspecific hybridization, and participate in the formation of novel chromosome rearrangements. Here we use standard cytological staining, cross-species chromosome painting, DNA probe analyses, and scanning electron microscopy to examine four interspecific macropodid hybrids (Macropus rufogriseus x Macropus agilis). The parental complements share the same centric fusions relative to the presumed macropodid ancestral karyotype, but can be differentiated on the basis of heterochromatic content, M. rufogriseus having larger centromeres with large C-banding positive regions. All hybrids exhibited the same pattern of chromosomal instability and remodeling specifically within the centromeres derived from the maternal (M. rufogriseus) complement. This instability included amplification of a satellite repeat and a transposable element, changes in chromatin structure, and de novo whole-arm rearrangements. We discuss possible reasons and mechanisms for the centromeric instability and remodeling observed in all four macropodid hybrids.

Genome-wide mapping and characterization of hypomethylated sites in human tissues and breast cancer cell lines

We have developed a method for mapping unmethylated sites in the human genome based on the resistance of TspRI-digested ends to ExoIII nuclease degradation. Digestion with TspRI and methylation-sensitive restriction endonuclease HpaII, followed by ExoIII and single-strand DNA nuclease allowed removal of DNA fragments containing unmethylated HpaII sites. We then used array comparative genomic hybridization (CGH) to map the sequences depleted by these procedures in human genomes derived from five human tissues, a primary breast tumor, and two breast tumor cell lines. Analysis of methylation patterns of the normal tissue genomes indicates that the hypomethylated sites are enriched in the 5' end of widely expressed genes, including promoter, first exon, and first intron. In contrast, genomes of the MCF-7 and MDA-MB-231 cell lines show extensive hypomethylation in the intragenic and intergenic regions whereas the primary tumor exhibits a pattern between those of the normal tissue and the cell lines. A striking characteristic of tumor cell lines is the presence of megabase-sized hypomethylated zones. These hypomethylated zones are associated with large genes, fragile sites, evolutionary breakpoints, chromosomal rearrangement breakpoints, tumor suppressor genes, and with regions containing tissue-specific gene clusters or with gene-poor regions containing novel tissue-specific genes. Correlation with microarray analysis shows that genes with a hypomethylated sequence 2 kb up- or downstream of the transcription start site are highly expressed, whereas genes with extensive intragenic and 3' untranslated region (UTR) hypomethylation are silenced. The method described herein can be used for large-scale screening of changes in the methylation pattern in the genome of interest.

The random versus fragile breakage models of chromosome evolution: a matter of resolution

Conserved synteny--the sharing of at least one orthologous gene by a pair of chromosomes from two species--can, in the strictest sense, be viewed as sequence conservation between chromosomes of two related species, irrespective of whether coding or non-coding sequence is examined. The recent sequencing of multiple vertebrate genomes indicates that certain chromosomal segments of considerable size are conserved in gene order as well as underlying non-coding sequence across all vertebrates. Some of these segments lost genes or non-coding sequence and/or underwent breakage only in teleost genomes, presumably because evolutionary pressure acting on these regions to remain intact were relaxed after an additional round of whole genome duplication. Random reporter insertions into zebrafish chromosomes combined with computational genome-wide analysis indicate that large chromosomal areas of multiple genes contain long-range regulatory elements, which act on their target genes from several gene distances away. In addition, computational breakpoint analyses suggest that recurrent evolutionary breaks are found in "fragile regions" or "hotspots", outside of the conserved blocks of synteny. These findings cannot be accommodated by the random breakage model and suggest that this view of genome and chromosomal evolution requires substantial reassessment.

Are there rearrangement hotspots in the human genome?

In a landmark paper, Nadeau and Taylor [18] formulated the random breakage model (RBM) of chromosome evolution that postulates that there are no rearrangement hotspots in the human genome. In the next two decades, numerous studies with progressively increasing levels of resolution made RBM the de facto theory of chromosome evolution. Despite the fact that RBM had prophetic prediction power, it was recently refuted by Pevzner and Tesler [4], who introduced the fragile breakage model (FBM), postulating that the human genome is a mosaic of solid regions (with low propensity for rearrangements) and fragile regions (rearrangement hotspots). However, the rebuttal of RBM caused a controversy and led to a split among researchers studying genome evolution. In particular, it remains unclear whether some complex rearrangements (e.g., transpositions) can create an appearance of rearrangement hotspots. We contribute to the ongoing debate by analyzing multi-break rearrangements that break a genome into multiple fragments and further glue them together in a new order. In particular, we demonstrate that (1) even if transpositions were a dominant force in mammalian evolution, the arguments in favor of FBM still stand, and (2) the "gene deletion" argument against FBM is flawed.

Chromosomal instability in Afrotheria: fragile sites, evolutionary breakpoints and phylogenetic inference from genome sequence assemblies

BACKGROUND

Extant placental mammals are divided into four major clades (Laurasiatheria, Supraprimates, Xenarthra and Afrotheria). Given that Afrotheria is generally thought to root the eutherian tree in phylogenetic analysis of large nuclear gene data sets, the study of the organization of the genomes of afrotherian species provides new insights into the dynamics of mammalian chromosomal evolution. Here we test if there are chromosomal bands with a high tendency to break and reorganize in Afrotheria, and by analyzing the expression of aphidicolin-induced common fragile sites in three afrotherian species, whether these are coincidental with recognized evolutionary breakpoints.

RESULTS

We described 29 fragile sites in the aardvark (OAF) genome, 27 in the golden mole (CAS), and 35 in the elephant-shrew (EED) genome. We show that fragile sites are conserved among afrotherian species and these are correlated with evolutionary breakpoints when compared to the human (HSA) genome. Inddition, by computationally scanning the newly released opossum (Monodelphis domestica) and chicken sequence assemblies for use as outgroups to Placentalia, we validate the HSA 3/21/5 chromosomal synteny as a rare genomic change that defines the monophyly of this ancient African clade of mammals. On the other hand, support for HSA 1/19p, which is also thought to underpin Afrotheria, is currently ambiguous.

CONCLUSION

We provide evidence that (i) the evolutionary breakpoints that characterise human syntenies detected in the basal Afrotheria correspond at the chromosomal band level with fragile sites, (ii) that HSA 3p/21 was in the amniote ancestor (i.e., common to turtles, lepidosaurs, crocodilians, birds and mammals) and was subsequently disrupted in the lineage leading to marsupials. Its expansion to include HSA 5 in Afrotheria is unique and (iii) that its fragmentation to HSA 3p/21 + HSA 5/21 in elephant and manatee was due to a fission within HSA 21 that is probably shared by all Paenungulata.

Comparative analysis of chicken chromosome 28 provides new clues to the evolutionary fragility of gene-rich vertebrate regions

The chicken genome draft sequence has provided a valuable resource for studies of an important agricultural and experimental model species and an important data set for comparative analysis. However, some of the most gene-rich segments are missing from chicken genome draft assemblies, limiting the analysis of a substantial number of genes and preventing a closer look at regions that are especially prone to syntenic rearrangements. To facilitate the functional and evolutionary analysis of one especially gene-rich, rearrangement-prone genomic region, we analyzed sequence from BAC clones spanning chicken microchromosome GGA28; as a complement we also analyzed a gene-sparse, stable region from GGA11. In these two regions we documented the conservation and lineage-specific gain and loss of protein-coding genes and precisely mapped the locations of 31 major human-chicken syntenic breakpoints. Altogether, we identified 72 lineage-specific genes, many of which are found at or near syntenic breaks, implicating evolutionary breakpoint regions as major sites of genetic innovation and change. Twenty-two of the 31 breakpoint regions have been reused repeatedly as rearrangement breakpoints in vertebrate evolution. Compared with stable GC-matched regions, GGA28 is highly enriched in CpG islands, as are break-prone intervals identified elsewhere in the chicken genome; evolutionary breakpoints are further enriched in GC content and CpG islands, highlighting a potential role for these features in genome instability. These data support the hypothesis that chromosome rearrangements have not occurred randomly over the course of vertebrate evolution but are focused preferentially within "fragile" regions with unusual DNA sequence characteristics.