Karyotype relationships among selected deer species and cattle revealed by bovine FISH probes

Comparative karyotype analysis of the red brocket deer (M. americana sensu lato and M. rufa) complex: evidence of drastic chromosomal evolution and implications on speciation process

Chromosomal rearrangements are often associated with playing a role in the speciation process. However, the underlying mechanism that favors the genetic isolation associated with chromosomal changes remains elusive. In this sense, the genus Mazama is recognized by its high level of karyotype diversity among species with similar morphology. A cryptic species complex has been identified within the genus, with the red brocket deer (Mazama americana and Mazama rufa) being the most impressive example. The chromosome variation was clustered in cytotypes with diploid numbers ranging from 42 to 53 and was correlated with geographical location. We conducted an analysis of chromosome evolution of the red brocket deer complex using comparative chromosome painting and Bacterial Artificial Chromosome (BAC) clones among different cytotypes. The aim was to deepen our understanding of the karyotypic relationships within the red brocket, thereby elucidating the significant chromosome variation among closely related species. This underscores the significance of chromosome changes as a key evolutionary process shaping their genomes. The results revealed the presence of three distinct cytogenetic lineages characterized by significant karyotypic divergence, suggesting the existence of efficient post-zygotic barriers. Tandem fusions constitute the main mechanism driving karyotype evolution, following a few centric fusions, inversion X-autosomal fusions. The BAC mapping has improved our comprehension of the karyotypic relationships within the red brocket deer complex, prompting questions regarding the role of these changes in the speciation process. We propose the red brocket as a model group to investigate how chromosomal changes contribute to isolation and explore the implications of these changes in taxonomy and conservation.

Comparative Genomics Reveal Phylogenetic Relationship and Chromosomal Evolutionary Events of Eight Cervidae Species

Cervidae represents a family that is not only rich in species diversity but also exhibits a wide range of karyotypes. The controversies regarding the phylogeny and classification of Cervidae still persist. The flourishing development of the genomic era has made it possible to address these issues at the genomic level. Here, the genomes of nine species were used to explore the phylogeny and chromosomal evolutionary events of Cervidae. By conducting whole-genome comparisons, we identified single-copy orthologous genes across the nine species and constructed a phylogenetic tree based on the single-copy orthologous genes sequences, providing new insights into the phylogeny of Cervidae, particularly the phylogenetic relationship among sika deer, red deer, wapiti and Tarim red deer. Gene family analysis revealed contractions in the olfactory receptor gene family and expansions in the histone gene family across eight Cervidae species. Furthermore, synteny analysis was used to explore the chromosomal evolutionary events of Cervidae species, revealing six chromosomal fissions during the evolutionary process from Bovidae to Cervidae. Notably, specific chromosomal fusion events were found in four species of Cervus, and a unique chromosomal fusion event was identified in Muntiacus reevesi. Our study further completed the phylogenetic relationship within the Cervidae and demonstrated the feasibility of inferring species phylogeny at the whole-genome level. Additionally, our findings on gene family evolution and the chromosomal evolutionary events in eight Cervidae species lay a foundation for comprehensive research of the evolution of Cervidae.

X Chromosome-Specific Repeats in Non-Domestic Bovidae

Repetitive sequences form a substantial and still enigmatic part of the mammalian genome. We isolated repetitive DNA blocks of the X chromosomes of three species of the family Bovidae: Kobus defassa (KDEXr sequence), Bos taurus (BTAXr sequence) and Antilope cervicapra (ACEXr sequence). The copy numbers of the isolated sequences were assessed using qPCR, and their chromosomal localisations were analysed using FISH in ten bovid tribes and in outgroup species. Besides their localisation on the X chromosome, their presence was also revealed on the Y chromosome and autosomes in several species. The KDEXr sequence abundant in most Bovidae species also occurs in distant taxa (Perissodactyla and Carnivora) and seems to be evolutionarily older than BTAXr and ACEXr. The ACEXr sequence, visible only in several Antilopini species using FISH, is probably the youngest, and arised in an ancestor common to Bovidae and Cervidae. All three repetitive sequences analysed in this study are interspersed among gene-rich regions on the X chromosomes, apparently preventing the crossing-over in their close vicinity. This study demonstrates that repetitive sequences on the X chromosomes have undergone a fast evolution, and their variation among related species can be beneficial for evolutionary studies.

Comparative studies of X chromosomes in Cervidae family

The family Cervidae is the second most diverse in the infraorder Pecora and is characterized by variability in the diploid chromosome numbers among species. X chromosomes in Cervidae evolved through complex chromosomal rearrangements of conserved segments within the chromosome, changes in centromere position, heterochromatic variation, and X-autosomal translocations. The family Cervidae consists of two subfamilies: Cervinae and Capreolinae. Here we build a detailed X chromosome map with 29 cattle bacterial artificial chromosomes of representatives of both subfamilies: reindeer (Rangifer tarandus), gray brocket deer (Mazama gouazoubira), Chinese water deer (Hydropotes inermis) (Capreolinae); black muntjac (Muntiacus crinifrons), tufted deer (Elaphodus cephalophus), sika deer (Cervus nippon) and red deer (Cervus elaphus) (Cervinae). To track chromosomal rearrangements during Cervidae evolution, we summarized new data, and compared them with available X chromosomal maps and chromosome level assemblies of other species. We demonstrate the types of rearrangements that may have underlined the variability of Cervidae X chromosomes. We detected two types of cervine X chromosome-acrocentric and submetacentric. The acrocentric type is found in three independent deer lineages (subfamily Cervinae and in two Capreolinae tribes-Odocoileini and Capreolini). We show that chromosomal rearrangements on the X-chromosome in Cervidae occur at a higher frequency than in the entire Ruminantia lineage: the rate of rearrangements is 2 per 10 million years.

Assessing the Taxonomic Status of the Gray Brocket Mazama simplicicornis argentina Lönnberg, 1919 (Artiodactyla: Cervidae)

Mazama simplicicornis argentina is the name that was given to describe a gray brocket collected by Lönberg in 1919 in the central Chaco region of Argentina. Subsequent authors, based on morphological similarities, considered this name to be a synonym for the species Subulo gouazoubira Fischer, 1814 from Paraguay. In the absence of genetic analyses to compare the Argentinian and Paraguayan gray brockets, we aimed to clarify the taxonomy of M. simplicicornis argentina through an integrative assessment using morphological, cytogenetical, and molecular data from its holotype and a current topotype. Qualitative skull features and cranio-morphometric results of M. simplicicornis argentina showed a great similarity with the S. gouazoubira neotype characters. The diploid chromosome number of M. simplicicornis argentina topotype corresponded with the karyotypical pattern of S. gouazoubira with 2n = 70 and FN = 70, showing a great similarity in all classic and molecular cytogenetic results and revealing the homologies between karyotypes. The phylogenetic analysis of mitochondrial genes used in this study (concatenated partial ND5 and Cytb gene) allocated the M. simplicicornis argentina specimens in the monophyletic clade of S. gouazoubira with a branch value of 100%. These results show that there is no discontinuity between the Argentinian and Paraguayan gray brockets. Therefore, the individuals originally described as M. simplicicornis argentina should be recognized as S. gouazoubira.

Revalidation of Passalites Gloger, 1841 for the Amazon brown brocket deer P.nemorivagus (Cuvier, 1817) (Mammalia, Artiodactyla, Cervidae)

Mazamanemorivaga (Cuvier, 1817) is a gray brocket deer that inhabits the Amazon region. An assessment of previous studies revealed inconsistencies in its current taxonomic classification, suggesting the need for an update in its genus classification. A taxonomic repositioning of this species is proposed through the collection of a specimen from its type locality (French Guiana) with subsequent morphological (coloring pattern, body measurements, and craniometry), cytogenetics (G Band, C Band, conventional Giemsa, Ag-NOR staining, and BAC probe mapping), and molecular phylogenetic analysis (mitochondrial genes Cyt B of 920 bp, COI I of 658 bp, D-loop 610 bp), and comparisons with other specimens of the same taxon, as well as other Neotropical deer species. The morphological and cytogenetic differences between this and other Neotropical Cervidae confirm the taxon as a unique and valid species. The phylogenetic analysis evidenced the basal position of the M.nemorivaga specimens within the Blastocerina clade. This shows early diversification and wide divergence from the other species, suggesting that the taxon should be transferred to a different genus. A taxonomic update of the genus name is proposed through the validation of Passalites Gloger, 1841, with Passalitesnemorivagus (Cuvier, 1817) as the type species. Future research should focus on evaluating the potential existence of other species within the genus Passalites, as suggested in the literature.

Chromosome-level assembly of the Rangifer tarandus genome and validation of cervid and bovid evolution insights

BACKGROUND

Genome assembly into chromosomes facilitates several analyses including cytogenetics, genomics and phylogenetics. Despite rapid development in bioinformatics, however, assembly beyond scaffolds remains challenging, especially in species without closely related well-assembled and available reference genomes. So far, four draft genomes of Rangifer tarandus (caribou or reindeer, a circumpolar distributed cervid species) have been published, but none with chromosome-level assembly. This emblematic northern species is of high interest in ecological studies and conservation since most populations are declining.

RESULTS

We have designed specific probes based on Oligopaint FISH technology to upgrade the latest published reindeer and caribou chromosome-level genomes. Using this oligonucleotide-based method, we found six mis-assembled scaffolds and physically mapped 68 of the largest scaffolds representing 78% of the most recent R. tarandus genome assembly. Combining physical mapping and comparative genomics, it was possible to document chromosomal evolution among Cervidae and closely related bovids.

CONCLUSIONS

Our results provide validation for the current chromosome-level genome assembly as well as resources to use chromosome banding in studies of Rangifer tarandus.

Resurrection of the genus Subulo Smith, 1827 for the gray brocket deer, with designation of a neotype

The gray brocket deer, Mazama gouazoubiraG. Fischer, 1814, occurs in South America and presents an extensive degree of morphological and genetic variability. Previous phylogenetic research showed that the genus Mazama is polyphyletic and imposed the designation of a different genus-group name for M. gouazoubira. We aimed to review and clarify the taxonomy of M. gouazoubira through the proposal of updating the nomenclature for this taxon and by the characterization of specimens collected close to the original type locality (topotypes). The topotypes were characterized by morphological (general characterization and morphometry), cytogenetic (conventional staining, Ag-NOR, G- and C-banding, and fluorescence in situ hybridization), and phylogenetic (mitogenomes) approaches. We revealed chromosome homologies between cattle and M. gouazoubira using an entire set of cattle whole chromosome painting probes and propose an updated G-band idiogram for the species. The morphometric analysis did not discriminate the individuals of M. gouazoubira, including the topotypes, from other small brocket deer species. However, the phylogenetic analysis, based on a Bayesian inference tree of the mitogenomes, confirmed the polyphyly of the genus Mazama and supported the need to change the gray brocket deer genus-group name. Based on our revision, we validated the genus SubuloSmith, 1827, and fixed a type species for the genus. In the absence of the holotype, we denominated a neotype described by the collection of a male topotype in Paraguay. The nomenclature rearrangement presented here is a starting point that will assist in the taxonomic resolution of Neotropical deer.

Revalidation of Mazama rufa (Illiger 1815) (Artiodactyla: Cervidae) as a Distinct Species out of the Complex Mazama americana (Erxleben 1777)

The red brocket deer Mazama americana Erxleben, 1777 is considered a polyphyletic complex of cryptic species with wide chromosomal divergence. Evidence indicates that the observed chromosomal divergences result in reproductive isolation. The description of a neotype for M. americana allowed its genetic characterization and represented a comparative basis to resolve the taxonomic uncertainties of the group. Thus, we designated a neotype for the synonym Mazama rufa Illiger, 1815 and tested its recognition as a distinct species from the M. americana complex with the analysis of morphological, cytogenetic and molecular data. We also evaluated its distribution by sampling fecal DNA in the wild. Morphological data from craniometry and body biometry indicated an overlap of quantitative measurements between M. rufa and the entire M. americana complex. The phylogenetic hypothesis obtained through mtDNA confirmed the reciprocal monophyly relationship between M. americana and M. rufa, and both were identified as distinct molecular operational taxonomic units by the General Mixed Yule Coalescent species delimitation analysis. Finally, classic cytogenetic data and fluorescence in situ hybridization with whole chromosome painting probes showed M. rufa with a karyotype of 2n = 52, FN = 56. Comparative analysis indicate that at least fifteen rearrangements separate M. rufa and M. americana (sensu stricto) karyotypes, which confirmed their substantial chromosomal divergence. This divergence should represent an important reproductive barrier and allow its characterization as a distinct and valid species. Genetic analysis of fecal samples demonstrated a wide distribution of M. rufa in the South American continent through the Atlantic Forest, Cerrado and south region of Amazon. Thus, we conclude for the revalidation of M. rufa as a distinct species under the concept of biological isolation, with its karyotype as the main diagnostic character. The present work serves as a basis for the taxonomic review of the M. americana complex, which should be mainly based on cytogenetic characterization and directed towards a better sampling of the Amazon region, the evaluation of available names in the species synonymy and a multi-locus phylogenetic analysis.

CNVs with adaptive potential in Rangifer tarandus: genome architecture and new annotated assembly

Rangifer tarandus has experienced recent drastic population size reductions throughout its circumpolar distribution and preserving the species implies genetic diversity conservation. To facilitate genomic studies of the species populations, we improved the genome assembly by combining long read and linked read and obtained a new highly accurate and contiguous genome assembly made of 13,994 scaffolds (L90 = 131 scaffolds). Using de novo transcriptome assembly of RNA-sequencing reads and similarity with annotated human gene sequences, 17,394 robust gene models were identified. As copy number variations (CNVs) likely play a role in adaptation, we additionally investigated these variations among 20 genomes representing three caribou ecotypes (migratory, boreal and mountain). A total of 1,698 large CNVs (length > 1 kb) showing a genome distribution including hotspots were identified. 43 large CNVs were particularly distinctive of the migratory and sedentary ecotypes and included genes annotated for functions likely related to the expected adaptations. This work includes the first publicly available annotation of the caribou genome and the first assembly allowing genome architecture analyses, including the likely adaptive CNVs reported here.

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 Summary

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

Abstract

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.

Chromosomal Aberrations in Cattle

Chromosomal aberrations and their mechanisms have been studied for many years in livestock. In cattle, chromosomal abnormalities are often associated with serious reproduction-related problems, such as infertility of carriers and early mortality of embryos. In the present work, we review the mechanisms and consequences of the most important bovine chromosomal aberrations: Robertsonian translocations and reciprocal translocations. We also discuss the application of bovine cell cultures in genotoxicity studies.

Analysis of human chromosomes by imaging flow cytometry

Chromosomal analysis is traditionally performed by karyotyping on metaphase spreads, or by fluorescent in situ hybridization (FISH) on interphase cells or metaphase spreads. Flow cytometry was introduced as a new method to analyze chromosomes number (ploidy) and structure (telomere length) in the 1970s with data interpretation largely based on fluorescence intensity. This technology has had little uptake for human cytogenetic applications primarily due to analytical challenges. The introduction of imaging flow cytometry, with the addition of digital images to standard multi-parametric flow cytometry quantitative tools, has added a new dimension. The ability to visualize the chromosomes and FISH signals overcomes the inherent difficulties when the data is restricted to fluorescence intensity. This field is now moving forward with methods being developed to assess chromosome number and structure in whole cells (normal and malignant) in suspension. A recent advance has been the inclusion of immunophenotyping such that antigen expression can be used to identify specific cells of interest for specific chromosomes and their abnormalities. This capability has been illustrated in blood cancers, such as chronic lymphocytic leukemia and plasma cell myeloma. The high sensitivity and specificity achievable highlights the potential imaging flow cytometry has for cytogenomic applications (i.e., diagnosis and disease monitoring). This review introduces and describes the development, current status, and applications of imaging flow cytometry for chromosomal analysis of human chromosomes.

Sperm chromosome segregation of rob(4;16) and rob(4;16)inv(4) in the brown brocket deer (Mazama gouazoubira)

The genus Mazama stands out among the Neotropical deer due to their wide intra and interspecific karyotypic diversification, which is associated with an accentuated chromosomal fragility. There are reports of heterozygous Robertsonian translocation (RT) carriers in a free-range population of Mazama gouazoubira (brown brocket deer), as well as in captive animals of this and other species of the genus. To analyze possible negative impacts of heterozygous chromosome rearrangements on reproductive fitness of the carriers, we performed an analysis of sperm meiotic segregation in four brown brocket bucks, carriers of a rob(4;16), and compared the results with those of a normal buck. We established a reliable FISH and sperm-FISH protocol for the brown brocket deer using bovine (Bos taurus; diploid number, 2n = 60) whole chromosome painting (WCP) and BAC probes. Using BAC probes, we revealed the presence of a paracentric inversion (PAI) of the fused chromosome 4 in two of the four analyzed RT carriers. The mean frequency of normal/balanced sperm in the translocation carriers was significantly lower than in the normal buck (94.78% vs 98.40%). The mean value of total unbalanced spermatozoa was almost doubled in the RT/PAI carriers (6.68%) when compared to RT carriers (3.76%), but the difference was not statistically significant. This study demonstrated the efficiency of FISH with bovine WCP and BAC probes in the characterization of chromosome rearrangements and gametic segregation patterns in brown brocket deer. Our results indicate a low to moderate increase in the rates of unbalanced meiotic segregation products in brown brocket bucks heterozygous for RT and RT/PAIs.

Analysis of muntjac deer genome and chromatin architecture reveals rapid karyotype evolution

Closely related muntjac deer show striking karyotype differences. Here we describe chromosome-scale genome assemblies for Chinese and Indian muntjacs, Muntiacus reevesi (2n = 46) and Muntiacus muntjak vaginalis (2n = 6/7), and analyze their evolution and architecture. The genomes show extensive collinearity with each other and with other deer and cattle. We identified numerous fusion events unique to and shared by muntjacs relative to the cervid ancestor, confirming many cytogenetic observations with genome sequence. One of these M. muntjak fusions reversed an earlier fission in the cervid lineage. Comparative Hi-C analysis showed that the chromosome fusions on the M. muntjak lineage altered long-range, three-dimensional chromosome organization relative to M. reevesi in interphase nuclei including A/B compartment structure. This reshaping of multi-megabase contacts occurred without notable change in local chromatin compaction, even near fusion sites. A few genes involved in chromosome maintenance show evidence for rapid evolution, possibly associated with the dramatic changes in karyotype.

Genetic architecture of individual variation in recombination rate on the X chromosome in cattle

Meiotic recombination is an essential biological process that ensures proper chromosome segregation and creates genetic diversity. Individual variation in global recombination rates has been shown to be heritable in several species, and variants significantly associated with this trait have been identified. Recombination on the sex chromosome has often been ignored in these studies although this trait may be particularly interesting as it may correspond to a biological process distinct from that on autosomes. For instance, recombination in males is restricted to the pseudo-autosomal region (PAR). We herein used a large cattle pedigree with more than 100,000 genotyped animals to improve the genetic map of the X chromosome and to study the genetic architecture of individual variation in recombination rate on the sex chromosome (XRR). The length of the genetic map was 46.4 and 121.2 cM in males and females, respectively, but the recombination rate in the PAR was six times higher in males. The heritability of CO counts on the X chromosome was comparable to that of autosomes in males (0.011) but larger than that of autosomes in females (0.024). XRR was highly correlated (0.76) with global recombination rate (GRR) in females, suggesting that both traits might be governed by shared variants. In agreement, a set of eleven previously identified variants associated with GRR had correlated effects on female XRR (0.86). In males, XRR and GRR appeared to be distinct traits, although more accurate CO counts on the PAR would be valuable to confirm these results.

Sequence Analysis and FISH Mapping of Four Satellite DNA Families among Cervidae

Centromeric and pericentromeric chromosome regions are occupied by satellite DNA. Satellite DNAs play essential roles in chromosome segregation, and, thanks to their extensive sequence variability, to some extent, they can also be used as phylogenetic markers. In this paper, we isolated and sequenced satellite DNA I-IV in 11 species of Cervidae. The obtained satellite DNA sequences and their chromosomal distribution were compared among the analysed representatives of cervid subfamilies Cervinae and Capreolinae. Only satI and satII sequences are probably present in all analysed species with high abundance. On the other hand, fluorescence in situ hybridisation (FISH) with satIII and satIV probes showed signals only in a part of the analysed species, indicating interspecies copy number variations. Several indices, including FISH patterns, the high guanine and cytosine (GC) content, and the presence of centromere protein B (CENP-B) binding motif, suggest that the satII DNA may represent the most important satellite DNA family that might be involved in the centromeric function in Cervidae. The absence or low intensity of satellite DNA FISH signals on biarmed chromosomes probably reflects the evolutionary reduction of heterochromatin following the formation of chromosome fusions. The phylogenetic trees constructed on the basis of the satellite I-IV DNA relationships generally support the present cervid taxonomy.

Comparative Chromosome Mapping of Musk Ox and the X Chromosome among Some Bovidae Species

: Bovidae, the largest family in Pecora infraorder, are characterized by a striking variability in diploid number of chromosomes between species and among individuals within a species. The bovid X chromosome is also remarkably variable, with several morphological types in the family. Here we built a detailed chromosome map of musk ox (Ovibos moschatus), a relic species originating from Pleistocene megafauna, with dromedary and human probes using chromosome painting. We trace chromosomal rearrangements during Bovidae evolution by comparing species already studied by chromosome painting. The musk ox karyotype differs from the ancestral pecoran karyotype by six fusions, one fission, and three inversions. We discuss changes in pecoran ancestral karyotype in the light of new painting data. Variations in the X chromosome structure of four bovid species nilgai bull (Boselaphus tragocamelus), saola (Pseudoryx nghetinhensis), gaur (Bos gaurus), and Kirk's Dikdik (Madoqua kirkii) were further analyzed using 26 cattle BAC-clones. We found the duplication on the X in saola. We show main rearrangements leading to the formation of four types of bovid X: Bovinae type with derived cattle subtype formed by centromere reposition and Antilopinae type with Caprini subtype formed by inversion in XSB3.

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.

Evidence for a Robertsonian fusion in Solea senegalensis (Kaup, 1858) revealed by zoo-FISH and comparative genome analysis

Background

Solea senegalensis (Kaup, 1858) is a commercially important flatfish species, belonging to the Pleuronectiformes order. The taxonomy of this group has long been controversial, and the karyotype of the order presents a high degree of variability in diploid number, derived from chromosomal rearrangements such as Robertsonian fusions. Previously it has been proposed that the large metacentric chromosome of S. senegalensis arises from this kind of chromosome rearrangement and that this is a proto-sex chromosome.

Results

In this work, the Robertsonian origin of the large metacentric chromosome of S. senegalensis has been tested by the Zoo-FISH technique applied to two species of the Soleidae family (Dicologlossa cuneata and Dagetichthys lusitanica), and by comparative genome analysis with Cynoglossus semilaevis. From the karyotypic analysis we were able to determine a chromosome complement comprising 2n = 50 (FN = 54) in D. cuneata and 2n = 42 (FN = 50) in D. lusitanica. The large metacentric painting probe gave consistent signals in four acrocentric chromosomes of the two Soleidae species; and the genome analysis proved a common origin with four chromosome pairs of C. semilaevis. As a result of the genomic analysis, up to 61 genes were annotated within the thirteen Bacterial Artificial Chromosome clones analysed.

Conclusions

These results confirm that the large metacentric chromosome of S. senegalensis originated from a Robertsonian fusion and provide new data about the chromosome evolution of S. senegalensis in particular, and of Pleuronectiformes in general.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5216-6) contains supplementary material, which is available to authorized users.