Evolutionary Breakpoints Research Articles

Comparative analysis of chromosome-level genomes provides insights into chromosomal evolution in Chiroptera.

Chiroptera (bats) presents a fascinating model due to its remarkable variation in chromosome numbers, which range from 14 to 62. This astonishing diversity makes bats an excellent subject for studying chromosome evolution. The black-bearded tomb bat (Taphozous melanopogon) occupies a pivotal phylogenetic position within Chiroptera, emphasizing its crucial role in the systematic examination of bat chromosome evolution. In this study, we present the first chromosome-level genome of T. melanopogon within the family Emballonuridae. Together with previously published genomes, we construct a strongly supported phylogenetic tree of bats, which supports that Emballonuridae forms a basal group within Yangochiroptera. Furthermore, we reconstruct ancestral karyotypes at key nodes along the bat phylogeny and conduct a synteny analysis among the genomes of 12 bat species. Our findings identified evolutionary breakpoint regions (EBRs) that are of particular interest. Notably, some bat genomes exhibit an enrichment of genes related to host defense against microbial pathogens within EBRs. Remarkably, one species possesses multiple copies of some β-defensin genes, while six other species have experienced the loss of some β-defensin genes due to EBRs. Furthermore, some olfactory receptor genes are located in EBRs of 12 species, 4 of which have a significant enrichment in sensory perception of smell. Together, our comparative genomic analysis underscores the potential link between chromosome rearrangements and the adaptation of bats to defend against microbial pathogens.

A Bird's-Eye View of Chromosomic Evolution in the Class Aves.

Birds (Aves) are the most speciose of terrestrial vertebrates, displaying Class-specific characteristics yet incredible external phenotypic diversity. Critical to agriculture and as model organisms, birds have adapted to many habitats. The only extant examples of dinosaurs, birds emerged ~150 mya and >10% are currently threatened with extinction. This review is a comprehensive overview of avian genome ("chromosomic") organization research based mostly on chromosome painting and BAC-based studies. We discuss traditional and contemporary tools for reliably generating chromosome-level assemblies and analyzing multiple species at a higher resolution and wider phylogenetic distance than previously possible. These results permit more detailed investigations into inter- and intrachromosomal rearrangements, providing unique insights into evolution and speciation mechanisms. The 'signature' avian karyotype likely arose ~250 mya and remained largely unchanged in most groups including extinct dinosaurs. Exceptions include Psittaciformes, Falconiformes, Caprimulgiformes, Cuculiformes, Suliformes, occasional Passeriformes, Ciconiiformes, and Pelecaniformes. The reasons for this remarkable conservation may be the greater diploid chromosome number generating variation (the driver of natural selection) through a greater possible combination of gametes and/or an increase in recombination rate. A deeper understanding of avian genomic structure permits the exploration of fundamental biological questions pertaining to the role of evolutionary breakpoint regions and homologous synteny blocks.

Insights into the adaptive evolution of chromosome and essential traits through chromosome-level genome assembly of Gekko japonicus

Gekko japonicus possesses flexible climbing and detoxification abilities under insectivorous habits. Still, the evolutionary mechanisms behind these traits remain unclarified. This study presents a chromosome-level G.japonicus genome, revealing that its evolutionary breakpoint regions were enriched with specific repetitive elements and defense response genes. Gene families unique to G.japonicus and positively selected genes are mainly enriched in immune, sensory, and nervous pathways. Expansion of bitter taste receptor type 2 primarily in insectivorous species could be associated with toxin clearance. Detox cytochrome P450 in G. japonicus has undergone more birth and death processes than biosynthesis-type P450 genes. Proline, cysteine, glycine, and serine in corneous beta proteins of G.japonicus might influence flexibility and setae adhesiveness. Certain thermosensitive transient receptor potential channels under relaxed purifying selection or positive selection in G. japonicus might enhance adaptation to climate change. This genome assembly offers insights into the adaptive evolution of gekkotans.

Patterns of chromosome evolution in ruminants

AbstractStudying when and where gross genomic rearrangements occurred during evolution is key to understanding changes in genome structure with functional consequences that might eventually lead to speciation. Here we identified chromosome rearrangements in ruminants, a clade characterized by large chromosome differences. Using 26 genome assemblies, we reconstructed five ancestral karyotypes and classified the rearrangement events occurring in each lineage. With these reconstructions, we then identified evolutionary breakpoints regions (EBRs) and synteny fragments. Ruminant karyotype evolution is characterized by inversions, while interchromosomal rearrangements occurred preferentially in the oldest ancestor of ruminants. We found that EBRs are depleted of protein coding genes, including housekeeping genes. Similarly, EBRs are not enriched in high GC regions, suggesting that meiotic double strand breaks might not be their origin. Overall, our results characterize at fine detail the location of chromosome rearrangements in ruminant evolution and provide new insights into the formation of EBRs.

97 Dissection of Evolution of Cis-Regulatory Elements and Its Application on Genetic Control of Complex Traits in Farm Animals

Abstract Functional Annotation of Animal Genomes (FAANG) has made significant progress on the improvement of the annotation of the functional elements, particularly non-coding regions in farm animal genomes. Different epigenomes including ChIP-seq for 4 histone marks, ATAC-seq, CTCF-seq, RNA-seq, and Hi-C across different tissues in different animal species were generated. Thousands of cis-regulatory elements including promoters and enhancers across different tissues in different species were identified and annotated. A comprehensive characterization and comparison of cis-regulatory elements across tissues and species can help us understand the genetic and evolutionary basis of complex phenotypes. Yet, attempts to dissect the evolution of regulatory elements have generally been small, both regarding the number of tissues and species compared. We demonstrate that genomic variants associated with complex traits and adaptive evolution in farm animals are significantly enriched in active promoters and enhancers. Further analysis suggests that the regulators were significantly enriched in evolutionary breakpoint regions (EBRs) in a tissue-specific and regulator-specific manner. Furthermore, we reveal distinct tissue-specific regulatory elements play important roles in domestication processes. Finally, we provide biological insights on tissue-specific regulatory conservation, and by integrating genome-wide association studies on complex traits, we demonstrate that comparative epigenomic analysis enhances our interpretations of complex traits variations in farm animals and provides novel insights on the role of regulatory elements on the functional conservation across species.

Integrative comparative analysis of avian chromosome evolution by in-silico mapping of the gene ontology of homologous synteny blocks and evolutionary breakpoint regions

Avian chromosomes undergo more intra- than interchromosomal rearrangements, which either induce or are associated with genome variations among birds. Evolving from a common ancestor with a karyotype not dissimilar from modern chicken, two evolutionary elements characterize evolutionary change: homologous synteny blocks (HSBs) constitute common conserved parts at the sequence level, while evolutionary breakpoint regions (EBRs) occur between HSBs, defining the points where rearrangement occurred. Understanding the link between the structural organization and functionality of HSBs and EBRs provides insight into the mechanistic basis of chromosomal change. Previously, we identified gene ontology (GO) terms associated with both; however, here we revisit our analyses in light of newly developed bioinformatic algorithms and the chicken genome assembly galGal6. We aligned genomes available for six birds and one lizard species, identifying 630 HSBs and 19 EBRs. We demonstrate that HSBs hold vast functionality expressed by GO terms that have been largely conserved through evolution. Particularly, we found that genes within microchromosomal HSBs had specific functionalities relevant to neurons, RNA, cellular transport and embryonic development, and other associations. Our findings suggest that microchromosomes may have conserved throughout evolution due to the specificity of GO terms within their HSBs. The detected EBRs included those found in the genome of the anole lizard, meaning they were shared by all saurian descendants, with others being unique to avian lineages. Our estimate of gene richness in HSBs supported the fact that microchromosomes contain twice as many genes as macrochromosomes.

Can changes in 3D genome architecture create new regulatory landscapes that contribute to phenotypic evolution?

Animal genomes are compartmentalized into insulated regulatory units named topology-associated domains (TADs). TADs insulate gene promoters from enhancers that occupy neighboring TADs. Chromosomal rearrangements that disrupt TAD structure can generate new regulatory interactions between enhancers and promoters that were once separated into different TADs, which might lead to new gene expression patterns. On the one hand, TAD rearrangements are known to cause deleterious phenotypes, but, on the other hand, rearrangements can also create novel expression patterns that may be selected during evolution because they generate advantageous phenotypes. Here, we review recent studies that explore the effects of chromosomal rearrangements and genetic perturbations on TAD structure and gene regulation in the context of development and evolution. We discuss the possible contribution of evolutionary breakpoints (EBRs) that affect TAD structure to the evolution of gene regulation and the phenotype.

Comparative Studies of Karyotypes in the Cervidae Family

The family Cervidae is the second most diverse family in the infraorder Pecora and is characterized by a striking variability in the diploid chromosome numbers among species, ranging from 6 to 70. Chromosomal rearrangements in Cervidae have been studied in detail by chromosome painting. There are many comparative cytogenetic data for both subfamilies (Cervinae and Capreolinae) based on homologies with chromosomes of cattle and Chinese muntjac. Previously it was found that interchromosomal rearrangements are the major type of rearrangements occurring in the Cervidae family. Here, we build a detailed chromosome map of a female reindeer (Rangifer tarandus, 2n = 70, Capreolinae) and a female black muntjac (Muntiacus crinifrons, 2n = 8, Cervinae) with dromedary homologies to find out what other types of rearrangements may have underlined the variability of Cervidae karyotypes. To track chromosomal rearrangements and the distribution of nucleolus organizer regions not only during Cervidae but also Pecora evolution, we summarized new data and compared them with chromosomal maps of other already studied species. We discuss changes in the pecoran ancestral karyotype in the light of new painting data. We show that intrachromosomal rearrangements in autosomes of Cervidae are more frequent than previously thought: at least 13 inversions in evolutionary breakpoint regions were detected.

Evolution of the ancestral mammalian karyotype and syntenic regions

Decrypting the rearrangements that drive mammalian chromosome evolution is critical to understanding the molecular bases of speciation, adaptation, and disease susceptibility. Using 8 scaffolded and 26 chromosome-scale genome assemblies representing 23/26 mammal orders, we computationally reconstructed ancestral karyotypes and syntenic relationships at 16 nodes along the mammalian phylogeny. Three different reference genomes (human, sloth, and cattle) representing phylogenetically distinct mammalian superorders were used to assess reference bias in the reconstructed ancestral karyotypes and to expand the number of clades with reconstructed genomes. The mammalian ancestor likely had 19 pairs of autosomes, with nine of the smallest chromosomes shared with the common ancestor of all amniotes (three still conserved in extant mammals), demonstrating a striking conservation of synteny for ∼320 My of vertebrate evolution. The numbers and types of chromosome rearrangements were classified for transitions between the ancestral mammalian karyotype, descendent ancestors, and extant species. For example, 94 inversions, 16 fissions, and 14 fusions that occurred over 53 My differentiated the therian from the descendent eutherian ancestor. The highest breakpoint rate was observed between the mammalian and therian ancestors (3.9 breakpoints/My). Reconstructed mammalian ancestor chromosomes were found to have distinct evolutionary histories reflected in their rates and types of rearrangements. The distributions of genes, repetitive elements, topologically associating domains, and actively transcribed regions in multispecies homologous synteny blocks and evolutionary breakpoint regions indicate that purifying selection acted over millions of years of vertebrate evolution to maintain syntenic relationships of developmentally important genes and regulatory landscapes of gene-dense chromosomes.

Genomic Organization of Microsatellites and LINE-1-like Retrotransposons: Evolutionary Implications for Ctenomys minutus (Rodentia: Ctenomyidae) Cytotypes.

Simple SummaryIn animals, several species contain substantial chromosomal and genomic variation among their populations, but as to what could have driven such diversification is still a puzzle for most cases. Here, we used molecular cytogenetic analysis to expose the main genomic elements involved in the population variation observed in the Neotropical underground rodents of the genus Ctenomys (Rodentia: Ctenomyidae), which harbor the most significant chromosomal variation among mammals (2n = 10 to 2n = 70). These data provide evidence for a correlation between repetitive genomic content and localization of evolutionary breakpoint regions (EBRs) and highlight their direct impact in promoting chromosomal rearrangements.The Neotropical underground rodents of the genus Ctenomys (Rodentia: Ctenomyidae) comprise about 65 species, which harbor the most significant chromosomal variation among mammals (2n = 10 to 2n = 70). Among them, C. minutus stands out with 45 different cytotypes already identified, among which, seven parental ones, named A to G, are parapatrically distributed in the coastal plains of Southern Brazil. Looking for possible causes that led to such extensive karyotype diversification, we performed chromosomal mapping of different repetitive DNAs, including microsatellites and long interspersed element-1 (LINE-1) retrotransposons in the seven parental cytotypes. Although microsatellites were found mainly in the centromeric and telomeric regions of the chromosomes, different patterns occur for each cytotype, thus revealing specific features. Likewise, the LINE-1-like retrotransposons also showed a differential distribution for each cytotype, which may be linked to stochastic loss of LINE-1 in some populations. Here, microsatellite motifs (A)30, (C)30, (CA)15, (CAC)10, (CAG)10, (CGG)10, (GA)15, and (GAG)10 could be mapped to fusion of chromosomes 20/17, fission and inversion in the short arm of chromosome 2, fusion of chromosomes 23/19, and different combinations of centric and tandem fusions of chromosomes 22/24/16. These data provide evidence for a correlation between repetitive genomic content and localization of evolutionary breakpoints and highlight their direct impact in promoting chromosomal rearrangements.

Enriched tandem repeats in chromosomal fusion points of Rineloricaria latirostris (Boulenger, 1900) (Siluriformes: Loricariidae).

Cytogenetic data showed the enrichment of repetitive DNAs in chromosomal rearrangement points between closely related species in armored catfishes. Still, few studies integrated cytogenetic and genomic data aiming to identify their prone-to-break DNA sites. Here, we aimed to obtain the repetitive fraction in Rineloricaria latirostris to recognize the microsatellite and homopolymers flanking the regions previously described as chromosomal fusion points. The results indicated that repetitive DNAs in R. latirostris are predominantly DNA transposons, and considering the microsatellite and homopolymers, A/T-rich expansions were the most abundant. The in situ localization demonstrated the A/T-rich repetitive sequences were scattered on the chromosomes, while A/G-rich microsatellite units were accumulated in some regions. The DNA transposon hAT, the 5S rDNA, and 45S rDNA (previously identified in Robertsonian fusion points in R. latirostris) were clusterized with some microsatellites, especially (CA)n, (GA)n, and poly-A, which were also enriched in regions of chromosomal fusions. Our findings demonstrated that repetitive sequences such as rDNAs, hAT transposons, and microsatellite units flank probable evolutionary breakpoint regions in R. latirostris. However, due to the sequence unit homologies in different chromosomal sites, these repeat DNAs only may facilitate chromosome fusion events in R. latirostris rather than working as a double-strand breakpoint site.

3D chromatin remodelling in the germ line modulates genome evolutionary plasticity

Chromosome folding has profound impacts on gene regulation, whose evolutionary consequences are far from being understood. Here we explore the relationship between 3D chromatin remodelling in mouse germ cells and evolutionary changes in genome structure. Using a comprehensive integrative computational analysis, we (i) reconstruct seven ancestral rodent genomes analysing whole-genome sequences of 14 species representatives of the major phylogroups, (ii) detect lineage-specific chromosome rearrangements and (iii) identify the dynamics of the structural and epigenetic properties of evolutionary breakpoint regions (EBRs) throughout mouse spermatogenesis. Our results show that EBRs are devoid of programmed meiotic DNA double-strand breaks (DSBs) and meiotic cohesins in primary spermatocytes, but are associated in post-meiotic cells with sites of DNA damage and functional long-range interaction regions that recapitulate ancestral chromosomal configurations. Overall, we propose a model that integrates evolutionary genome reshuffling with DNA damage response mechanisms and the dynamic spatial genome organisation of germ cells.

Evolutionary breakpoint regions and chromosomal remodeling in Harttia (Siluriformes: Loricariidae) species diversification.

The Neotropical armored catfish genus Harttia presents a wide variation of chromosomal rearrangements among its representatives. Studies indicate that translocation and Robertsonian rearrangements have triggered the karyotype evolution in the genus, including differentiation of sex chromosome systems. However, few studies used powerful tools, such as comparative whole chromosome painting, to clarify this highly diversified scenario. Here, we isolated probes from the X1 (a 5S rDNA carrier) and the X2 (a 45S rDNA carrier) chromosomes of Harttia punctata, which displays an X1X1X2X2/X1X2Y multiple sex chromosome system. Those probes were applied in other Harttia species to evidence homeologous chromosome blocks. The resulting data reinforce that translocation events played a role in the origin of the X1X2Y sex chromosome system in H. punctata. The repositioning of homologous chromosomal blocks carrying rDNA sites among ten Harttia species has also been demonstrated. Anchored to phylogenetic data it was possible to evidence some events of the karyotype diversification of the studied species and to prove an independent origin for the two types of multiple sex chromosomes, XX/XY1Y2 and X1X1X2X2/X1X2Y, that occur in Harttia species. The results point to evolutionary breakpoint regions in the genomes within or adjacent to rDNA sites that were widely reused in Harttia chromosome remodeling.

Recurrent erosion of COA1/MITRAC15 exemplifies conditional gene dispensability in oxidative phosphorylation

Skeletal muscle fibers rely upon either oxidative phosphorylation or the glycolytic pathway with much less reliance on oxidative phosphorylation to achieve muscular contractions that power mechanical movements. Species with energy-intensive adaptive traits that require sudden bursts of energy have a greater dependency on glycolytic fibers. Glycolytic fibers have decreased reliance on OXPHOS and lower mitochondrial content compared to oxidative fibers. Hence, we hypothesized that gene loss might have occurred within the OXPHOS pathway in lineages that largely depend on glycolytic fibers. The protein encoded by the COA1/MITRAC15 gene with conserved orthologs found in budding yeast to humans promotes mitochondrial translation. We show that gene disrupting mutations have accumulated within the COA1 gene in the cheetah, several species of galliform birds, and rodents. The genomic region containing COA1 is a well-established evolutionary breakpoint region in mammals. Careful inspection of genome assemblies of closely related species of rodents and marsupials suggests two independent COA1 gene loss events co-occurring with chromosomal rearrangements. Besides recurrent gene loss events, we document changes in COA1 exon structure in primates and felids. The detailed evolutionary history presented in this study reveals the intricate link between skeletal muscle fiber composition and the occasional dispensability of the chaperone-like role of the COA1 gene.

Anchoring the CerEla1.0 Genome Assembly to Red Deer (Cervus elaphus) and Cattle (Bos taurus) Chromosomes and Specification of Evolutionary Chromosome Rearrangements in Cervidae.

Simple SummaryThe red deer (Cervus elaphus) de novo genome assembly (CerEla1.0) has provided a great resource for genetic studies in various deer species. In this study, we used gene order comparisons between C. elaphus CerEla1.0 and B. taurus ARS-UCD1.2 genome assemblies and fluorescence in situ hybridization (FISH) with bovine BAC probes to verify the red deer-bovine chromosome relationships and anchor the CerEla1.0 C-scaffolds to karyotypes of both species. We showed the homology between bovine and deer chromosomes and determined the centromere-telomere orientation of the CerEla1.0 C-scaffolds. Using a set of BAC probes, we were able to narrow the positions of evolutionary chromosome breakpoints defining the family Cervidae. In addition, we revealed several errors in the current CerEla1.0 genome assembly. Finally, we expanded our analysis to other Cervidae and confirmed the locations of the cervid evolutionary fissions and orientation of the fused chromosomes in eight cervid species. Our results can serve as a basis for necessary improvements of the red deer genome assembly and provide support to other genetic studies in Cervidae.The family Cervidae groups a range of species with an increasing economic significance. Their karyotypes share 35 evolutionary conserved chromosomal segments with cattle (Bos taurus). Recent publication of the annotated red deer (Cervus elaphus) whole genome assembly (CerEla1.0) has provided a basis for advanced genetic studies. In this study, we compared the red deer CerEla1.0 and bovine ARS-UCD1.2 genome assembly and used fluorescence in situ hybridization with bovine BAC probes to verify the homology between bovine and deer chromosomes, determined the centromere-telomere orientation of the CerEla1.0 C-scaffolds and specified positions of the cervid evolutionary chromosome breakpoints. In addition, we revealed several incongruences between the current deer and bovine genome assemblies that were shown to be caused by errors in the CerEla1.0 assembly. Finally, we verified the centromere-to-centromere orientation of evolutionarily fused chromosomes in seven additional deer species, giving a support to previous studies on their chromosome evolution.

A map of cis-regulatory elements and 3D genome structures in zebrafish.

The zebrafish (Danio rerio) has been widely used in the study of human disease and development, and about 70% of the protein-coding genes are conserved between the two species1. However, studies in zebrafish remain constrained by the sparse annotation of functional control elements in the zebrafish genome. Here we performed RNA sequencing, assay for transposase-accessible chromatin using sequencing (ATAC-seq), chromatin immunoprecipitation with sequencing, whole-genome bisulfite sequencing, and chromosome conformation capture (Hi-C) experiments in up to eleven adult and two embryonic tissues to generate a comprehensive map of transcriptomes, cis-regulatory elements, heterochromatin, methylomes and 3D genome organization in the zebrafish Tübingen reference strain. A comparison of zebrafish, human and mouse regulatory elements enabled the identification of both evolutionarily conserved and species-specific regulatory sequences and networks. We observed enrichment of evolutionary breakpoints at topologically associating domain boundaries, which were correlated with strong histone H3 lysine 4 trimethylation (H3K4me3) and CCCTC-binding factor (CTCF) signals. We performed single-cell ATAC-seq in zebrafish brain, which delineated 25 different clusters of cell types. By combining long-read DNA sequencing and Hi-C, we assembled the sex-determining chromosome 4 de novo. Overall, our work provides an additional epigenomic anchor for the functional annotation of vertebrate genomes and the study of evolutionarily conserved elements of 3D genome organization.

Circular DNA intermediates in the generation of large human segmental duplications

BackgroundDuplications of large genomic segments provide genetic diversity in genome evolution. Despite their importance, how these duplications are generated remains uncertain, particularly for distant duplicated genomic segments.ResultsHere we provide evidence of the participation of circular DNA intermediates in the single generation of some large human segmental duplications. A specific reversion of sequence order from A-B/C-D to B-A/D-C between duplicated segments and the presence of only microhomologies and short indels at the evolutionary breakpoints suggest a circularization of the donor ancestral locus and an accidental replicative interaction with the acceptor locus.ConclusionsThis novel mechanism of random genomic mutation could explain several distant genomic duplications including some of the ones that took place during recent human evolution.

Comparative cytogenetic analyses in Ancistrus species (Siluriformes: Loricariidae)

ABSTRACT Ancistrus is a specious genus of armored catfishes that has been extensively used for cytogenetic studies in the last 17 years. A comparison of the extensive karyotypic plasticity within this genus is presented with new cytogenetic analysis for Ancistrus cf. multispinis and Ancistrus aguaboensis. This study aims to improve our understanding of chromosomal evolution associated with changes in the diploid number (2n) and the dispersion of ribosomal DNAs (rDNAs) within Ancistrus. Ancistrus cf. multispinis and A. aguaboensis exhibit 2n of 52 and 50 chromosomes, respectively. Given that A. cf. multispinis shares a 2n = 52 also found in Pterygoplichthyini, the sister group for Ancistrini, a Robertsonian (Rb) fusion event is proposed for the 2n reduction in A. aguaboensis. 5S rDNAs pseudogenes sites have already been associated with Rb fusion in Ancistrus and our analysis suggests that the 2n reduction in A. aguaboensis was triggered by double strand breaks (DSBs) and chromosomal rearrangements at 5S rDNA sites. The presence of evolutionary breakpoint regions (EBRs) into rDNA cluster is proposed to explain part of the Rb fusion in Ancistrus. Cytogenetic data presented extends the diversity already documented in Ancistrus to further understand the role of chromosomal rearrangements in the diversification of Ancistrini.

Chromosome-level genome assembly for giant panda provides novel insights into Carnivora chromosome evolution

BackgroundChromosome evolution is an important driver of speciation and species evolution. Previous studies have detected chromosome rearrangement events among different Carnivora species using chromosome painting strategies. However, few of these studies have focused on chromosome evolution at a nucleotide resolution due to the limited availability of chromosome-level Carnivora genomes. Although the de novo genome assembly of the giant panda is available, current short read-based assemblies are limited to moderately sized scaffolds, making the study of chromosome evolution difficult.ResultsHere, we present a chromosome-level giant panda draft genome with a total size of 2.29 Gb. Based on the giant panda genome and published chromosome-level dog and cat genomes, we conduct six large-scale pairwise synteny alignments and identify evolutionary breakpoint regions. Interestingly, gene functional enrichment analysis shows that for all of the three Carnivora genomes, some genes located in evolutionary breakpoint regions are significantly enriched in pathways or terms related to sensory perception of smell. In addition, we find that the sweet receptor gene TAS1R2, which has been proven to be a pseudogene in the cat genome, is located in an evolutionary breakpoint region of the giant panda, suggesting that interchromosomal rearrangement may play a role in the cat TAS1R2 pseudogenization.ConclusionsWe show that the combined strategies employed in this study can be used to generate efficient chromosome-level genome assemblies. Moreover, our comparative genomics analyses provide novel insights into Carnivora chromosome evolution, linking chromosome evolution to functional gene evolution.

Evolution of gene regulation in ruminants differs between evolutionary breakpoint regions and homologous synteny blocks.

The role of chromosome rearrangements in driving evolution has been a long-standing question of evolutionary biology. Here we focused on ruminants as a model to assess how rearrangements may have contributed to the evolution of gene regulation. Using reconstructed ancestral karyotypes of Cetartiodactyls, Ruminants, Pecorans, and Bovids, we traced patterns of gross chromosome changes. We found that the lineage leading to the ruminant ancestor after the split from other cetartiodactyls was characterized by mostly intrachromosomal changes, whereas the lineage leading to the pecoran ancestor (including all livestock ruminants) included multiple interchromosomal changes. We observed that the liver cell putative enhancers in the ruminant evolutionary breakpoint regions are highly enriched for DNA sequences under selective constraint acting on lineage-specific transposable elements (TEs) and a set of 25 specific transcription factor (TF) binding motifs associated with recently active TEs. Coupled with gene expression data, we found that genes near ruminant breakpoint regions exhibit more divergent expression profiles among species, particularly in cattle, which is consistent with the phylogenetic origin of these breakpoint regions. This divergence was significantly greater in genes with enhancers that contain at least one of the 25 specific TF binding motifs and located near bovidae-to-cattle lineage breakpoint regions. Taken together, by combining ancestral karyotype reconstructions with analysis of cis regulatory element and gene expression evolution, our work demonstrated that lineage-specific regulatory elements colocalized with gross chromosome rearrangements may have provided valuable functional modifications that helped to shape ruminant evolution.