Chromosome-length genome assemblies and cytogenomic analyses of pangolins reveal remarkable chromosome counts and plasticity.

We report the first chromosome-length genome assemblies for three species in the mammalian order Pholidota: the white-bellied, Chinese, and Sunda pangolins. Surprisingly, we observe extraordinary karyotypic plasticity within this order and, in female white-bellied pangolins, the largest number of chromosomes reported in a Laurasiatherian mammal: 2n = 114. We perform the first karyotype analysis of an African pangolin and report a Y-autosome fusion in white-bellied pangolins, resulting in 2n = 113 for males. We employ a novel strategy to confirm the fusion and identify the autosome involved by finding the pseudoautosomal region (PAR) in the female genome assembly and analyzing the 3D contact frequency between PAR sequences and the rest of the genome in male and female white-bellied pangolins. Analyses of genetic variability show that white-bellied pangolins have intermediate levels of genome-wide heterozygosity relative to Chinese and Sunda pangolins, consistent with two moderate declines of historical effective population size. Our results reveal a remarkable feature of pangolin genome biology and highlight the need for further studies of these unique and endangered mammals.

Chromosome-length genome assembly and karyotype of the endangered black-footed ferret (Mustela nigripes).

The black-footed ferret (Mustela nigripes) narrowly avoided extinction to become an oft-cited example of the benefits of intensive management, research, and collaboration to save a species through ex situ conservation breeding and reintroduction into its former range. However, the species remains at risk due to possible inbreeding, disease susceptibility, and multiple fertility challenges. Here, we report the de novo genome assembly of a male black-footed ferret generated through a combination of linked-read sequencing, optical mapping, and Hi-C proximity ligation. In addition, we report the karyotype for this species, which was used to anchor and assign chromosome numbers to the chromosome-length scaffolds. The draft assembly was ~2.5 Gb in length, with 95.6% of it anchored to 19 chromosome-length scaffolds, corresponding to the 2n = 38 chromosomes revealed by the karyotype. The assembly has contig and scaffold N50 values of 148.8 kbp and 145.4 Mbp, respectively, and is up to 96% complete based on BUSCO analyses. Annotation of the assembly, including evidence from RNA-seq data, identified 21,406 protein-coding genes and a repeat content of 37.35%. Phylogenomic analyses indicated that the black-footed ferret diverged from the European polecat/domestic ferret lineage 1.6 million yr ago. This assembly will enable research on the conservation genomics of black-footed ferrets and thereby aid in the further restoration of this endangered species.

Demographic History, Adaptation, and NRAP Convergent Evolution at Amino Acid Residue 100 in the World Northernmost Cattle from Siberia.

Native cattle breeds represent an important cultural heritage. They are a reservoir of genetic variation useful for properly responding to agriculture needs in the light of ongoing climate changes. Evolutionary processes that occur in response to extreme environmental conditions could also be better understood using adapted local populations. Herein, different evolutionary histories of the world northernmost native cattle breeds from Russia were investigated. They highlighted Kholmogory as a typical taurine cattle, whereas Yakut cattle separated from European taurines approximately 5,000 years ago and contain numerous ancestral and some novel genetic variants allowing their adaptation to harsh conditions of living above the Polar Circle. Scans for selection signatures pointed to several common gene pathways related to adaptation to harsh climates in both breeds. But genes affected by selection from these pathways were mostly different. A Yakut cattle breed-specific missense mutation in a highly conserved NRAP gene represents a unique example of a young amino acid residue convergent change shared with at least 16 species of hibernating/cold-adapted mammals from six distinct phylogenetic orders. This suggests a convergent evolution event along the mammalian phylogenetic tree and fast fixation in a single isolated cattle population exposed to a harsh climate.

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.

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 Tandemly Arranged Repetitive DNAs in Three Species of Cyprinoidei with Different Ploidy Levels.

Polyploid species represent a challenge for both cytogenetic and genomic studies due to their high chromosome numbers and the morphological similarity between their paralogous chromosomes. This paper describes the use of low-coverage high-throughput sequencing to identify the 14 most abundant tandemly arranged repetitive elements in the paleotetraploid genome of the crucian carp (Carassius carassius, 2n = 100). These repetitive elements were then used for molecular cytogenetic studies of a closely related functionally triploid form of the Prussian carp (Carassius gibelio, 3n = 150 + Bs) and a relatively distant diploid species, the tench (Tinca tinca, 2n = 48). According to their distribution on the chromosomes of the 3 aforementioned species, the repetitive elements here identified can be divided into 5 groups: (1) those specific to a single genomic locus in both Carassius species, despite the recent carp-specific genome duplication; (2) those located in a single genomic locus of T. tinca, but amplified in one or both Carassius species; (3) those massively amplified in the B chromosomes of C. gibelio; (4) those located in a single locus in C. gibelio, but amplified in many blocks in C. carassius; and (5) those located in multiple pericentromeric loci in both Carassius species. Our data indicate that some of the repetitive elements are highly conserved in cyprinoid species and may serve as good cytogenetic and genomic markers for discriminating paralogous chromosomes, while others are evolutionarily recent, and their amplification may be related to the last whole-genome duplication event.

Draft de novo Genome Assembly of the Elusive Jaguarundi, Puma yagouaroundi.

The Puma lineage within the family Felidae consists of 3 species that last shared a common ancestor around 4.9 million years ago. Whole-genome sequences of 2 species from the lineage were previously reported: the cheetah (Acinonyx jubatus) and the mountain lion (Puma concolor). The present report describes a whole-genome assembly of the remaining species, the jaguarundi (Puma yagouaroundi). We sequenced the genome of a male jaguarundi with 10X Genomics linked reads and assembled the whole-genome sequence. The assembled genome contains a series of scaffolds that reach the length of chromosome arms and is similar in scaffold contiguity to the genome assemblies of cheetah and puma, with a contig N50 = 100.2 kbp and a scaffold N50 = 49.27 Mbp. We assessed the assembled sequence of the jaguarundi genome using BUSCO, aligned reads of the sequenced individual and another published female jaguarundi to the assembled genome, annotated protein-coding genes, repeats, genomic variants and their effects with respect to the protein-coding genes, and analyzed differences of the 2 jaguarundis from the reference mitochondrial genome. The jaguarundi genome assembly and its annotation were compared in quality, variants, and features to the previously reported genome assemblies of puma and cheetah. Computational analyzes used in the study were implemented in transparent and reproducible way to allow their further reuse and modification.

Low-pass single-chromosome sequencing of human small supernumerary marker chromosomes (sSMCs) and Apodemus B chromosomes.

Supernumerary chromosomes sporadically arise in many eukaryotic species as a result of genomic rearrangements. If present in a substantial part of species population, those are called B chromosomes, or Bs. This is the case for 70 mammalian species, most of which are rodents. In humans, the most common types of extra chromosomes, sSMCs (small supernumerary marker chromosomes), are diagnosed in approximately 1 of 2000 postnatal cases. Due to low frequency in population, human sSMCs are not considered B chromosomes. Genetic content of both B-chromosomes and sSMCs in most cases remains understudied. Here, we apply microdissection of single chromosomes with subsequent low-pass sequencing on Ion Torrent PGM and Illumina MiSeq to identify unique and repetitive DNA sequences present in a single human sSMC and several B chromosomes in mice Apodemus flavicollis and Apodemus peninsulae. The pipeline for sequencing data analysis was made available in Galaxy interface as an addition to previously published command-line version. Human sSMC was attributed to the proximal part of chromosome 15 long arm, and breakpoints leading to its formation were located into satellite DNA arrays. Genetic content of Apodemus B chromosomes was species-specific, and minor alterations were observed in both species. Common features of Bs in these Apodemus species were satellite DNA and ERV enrichment, as well as the presence of the vaccinia-related kinase gene Vrk1. Understanding of the non-essential genome elements content provides important insights into genome evolution in general.

First report on B chromosome content in a reptilian species: the case of Anolis carolinensis.

Supernumerary elements of the genome are often called B chromosomes. They usually consist of various autosomal sequences and, because of low selective pressure, are mostly pseudogenized and contain many repeats. There are numerous reports on B chromosomes in mammals, fish, invertebrates, plants, and fungi, but only a few of them have been studied using sequencing techniques. However, reptilian supernumerary chromosomes have been detected only cytogenetically and never sequenced or analyzed at the molecular level. One model squamate species with available genome sequence is Anolis carolinensis. The scope of the present article is to describe the genetic content of A. carolinensis supernumerary chromosomes. In this article, we confirm the presence of B chromosomes in this species by reverse painting and synaptonemal complex analysis. We applied low-pass high-throughput sequencing to analyze flow-sorted B chromosomes. Anole B chromosomes exhibit similar traits to other supernumerary chromosomes from different taxons: they contain two genes related to cell division control (INCENP and SPIRE2), are enriched in specific repeats, and show a high degree of pseudogenization. Therefore, the present study confirms that reptilian B chromosomes resemble supernumerary chromosomes of other taxons.

A combined banding method that allows the reliable identification of chromosomes as well as differentiation of AT- and GC-rich heterochromatin.

Сonstitutive heterochromatin areas are revealed by differential staining as C-positive chromosomal regions. These C-positive bands may greatly vary by location, size, and nucleotide composition. CBG-banding is the most commonly used method to detect structural heterochromatin in animals. The difficulty in identification of individual chromosomes represents an unresolved problem of this method as the body of the chromosome is stained uniformly and does not have banding pattern beyond C-bands. Here, we present the method that we called CDAG for sequential heterochromatin staining after differential GTG-banding. The method uses G-banding followed by heat denaturation in the presence of formamide with consecutive fluorochrome staining. The new technique is valid for the concurrent revealing of heterochromatin position due to differential banding of chromosomes and heterochromatin composition (AT-/GC-rich) in animal karyotyping.

Chromosome Translocations as a Driver of Diversification in Mole Voles Ellobius (Rodentia, Mammalia).

The involvement of chromosome changes in the initial steps of speciation is controversial. Here we examine diversification trends within the mole voles Ellobius, a group of subterranean rodents. The first description of their chromosome variability was published almost 40 years ago. Studying the G-band structure of chromosomes in numerous individuals revealed subsequent homologous, step-by-step, Robertsonian translocations, which changed diploid numbers from 54 to 30. Here we used a molecular cytogenetic strategy which demonstrates that chromosomal translocations are not always homologous; consequently, karyotypes with the same diploid number can carry different combinations of metacentrics. We further showed that at least three chromosomal forms with 2n = 34 and distinct metacentrics inhabit the Pamir-Alay mountains. Each of these forms independently hybridized with E. tancrei, 2n = 54, forming separate hybrid zones. The chromosomal variations correlate slightly with geographic barriers. Additionally, we confirmed that the emergence of partial or monobrachial homology appeared to be a strong barrier for hybridization in nature, in contradistinction to experiments which we reported earlier. We discuss the possibility of whole arm reciprocal translocations for mole voles. Our findings suggest that chromosomal translocations lead to diversification and speciation.

De novo assembling and primary analysis of genome and transcriptome of gray whale Eschrichtius robustus.

BACKGROUND: Gray whale, Eschrichtius robustus (E. robustus), is a single member of the family Eschrichtiidae, which is considered to be the most primitive in the class Cetacea. Gray whale is often described as a "living fossil". It is adapted to extreme marine conditions and has a high life expectancy (77 years). The assembly of a gray whale genome and transcriptome will allow to carry out further studies of whale evolution, longevity, and resistance to extreme environment. RESULTS: In this work, we report the first de novo assembly and primary analysis of the E. robustus genome and transcriptome based on kidney and liver samples. The presented draft genome assembly is complete by 55% in terms of a total genome length, but only by 24% in terms of the BUSCO complete gene groups, although 10,895 genes were identified. Transcriptome annotation and comparison with other whale species revealed robust expression of DNA repair and hypoxia-response genes, which is expected for whales. CONCLUSIONS: This preliminary study of the gray whale genome and transcriptome provides new data to better understand the whale evolution and the mechanisms of their adaptation to the hypoxic conditions.

Evolutionary plasticity of acipenseriform genomes.

Acipenseriformes is an order of ray-finned fishes, comprising 27 extant species of sturgeons and paddlefishes inhabiting waters of the Northern Hemisphere. The order has a basal position within Actinopteri (ray-finned fish minus polypterids) and is characterized by many specific morphological and genomic features, including high diploid chromosome numbers, various levels of ploidy between species, unclear sex determination, and propensity to interspecific hybridization. Recent advances in molecular genetics, genomics, and comparative cytogenetics produced novel data on different aspects of acipenseriform biology, including improved phylogenetic reconstructions and better understanding of genome structure. Here, we discuss the cytogenetic and genomic traits of acipenseriforms and their connection with polyploidization and tolerance to interspecific hybridization.

Genomic Organization and Physical Mapping of Tandemly Arranged Repetitive DNAs in Sterlet (Acipenser ruthenus).

Acipenseriformes represent a phylogenetically basal clade of ray-finned fish characterized by unusual genomic traits, including paleopolyploid states of extant genomes with high chromosome numbers and slow rates of molecular evolution. Despite a high interest in this fish group, only a limited number of studies have been accomplished on the isolation and characterization of repetitive DNA, karyotype standardization is not yet complete, and sex chromosomes are still to be identified. Here, we applied next-generation sequencing and cluster analysis to characterize major fractions of sterlet (Acipenser ruthenus) repetitive DNA. Using FISH, we mapped 16 tandemly arranged sequences on sterlet chromosomes and found them to be unevenly distributed in the genome with a tendency to cluster in particular regions. Some of the satellite DNAs might be used as specific markers to identify individual chromosomes and their paralogs, resulting in the unequivocal identification of at least 18 chromosome pairs. Our results provide an insight into the characteristic genomic distribution of the most common sterlet repetitive sequences. Biased accumulation of repetitive DNAs in particular chromosomes makes them especially interesting for further search for cryptic sex chromosomes. Future studies of these sequences in other acipenserid species will provide new perspectives regarding the evolution of repetitive DNA within the genomes of this fish order.

Genome-wide comparative chromosome maps of Arvicola amphibius, Dicrostonyx torquatus, and Myodes rutilus.

The subfamily Arvicolinae consists of a great number of species with highly diversified karyotypes. In spite of the wide use of arvicolines in biological and medicine studies, the data on their karyotype structures are limited. Here, we made a set of painting probes from flow-sorted chromosomes of a male Palearctic collared lemming (Dicrostonyx torquatus, DTO). Together with the sets of painting probes made previously from the field vole (Microtus agrestis, MAG) and golden hamster (Mesocricetus auratus, MAU), we carried out a reciprocal chromosome painting between these three species. The three sets of probes were further hybridized onto the chromosomes of the Eurasian water vole (Arvicola amphibius) and northern red-backed vole (Myodes rutilus). We defined the diploid chromosome number in D. torquatus karyotype as 2n = 45 + Bs and showed that the system of sex chromosomes is X1X2Y1. The probes developed here provide a genomic tool-kit, which will help to investigate the evolutionary biology of the Arvicolinae rodents. Our results show that the syntenic association MAG1/17 is present not only in Arvicolinae but also in some species of Cricetinae; and thus, should not be considered as a cytogenetic signature for Arvicolinae. Although cytogenetic signature markers for the genera have not yet been found, our data provides insight into the likely ancestral karyotype of Arvicolinae. We conclude that the karyotypes of modern voles could have evolved from a common ancestral arvicoline karyotype (AAK) with 2n = 56 mainly by centric fusions and fissions.

Contrasting origin of B chromosomes in two cervids (Siberian roe deer and grey brocket deer) unravelled by chromosome-specific DNA sequencing.

BACKGROUND: B chromosomes are dispensable and variable karyotypic elements found in some species of animals, plants and fungi. They often originate from duplications and translocations of host genomic regions or result from hybridization. In most species, little is known about their DNA content. Here we perform high-throughput sequencing and analysis of B chromosomes of roe deer and brocket deer, the only representatives of Cetartiodactyla known to have B chromosomes. RESULTS: In this study we developed an approach to identify genomic regions present on chromosomes by high-throughput sequencing of DNA generated from flow-sorted chromosomes using degenerate-oligonucleotide-primed PCR. Application of this method on small cattle autosomes revealed a previously described KIT gene region translocation associated with colour sidedness. Implementing this approach to B chromosomes from two cervid species, Siberian roe deer (Capreolus pygargus) and grey brocket deer (Mazama gouazoubira), revealed dramatically different genetic content: roe deer B chromosomes consisted of two duplicated genomic regions (a total of 1.42-1.98 Mbp) involving three genes, while grey brocket deer B chromosomes contained 26 duplicated regions (a total of 8.28-9.31 Mbp) with 34 complete and 21 partial genes, including KIT and RET protooncogenes, previously found on supernumerary chromosomes in canids. Sequence variation analysis of roe deer B chromosomes revealed a high frequency of mutations and increased heterozygosity due to either amplification within B chromosomes or divergence between different Bs. In contrast, grey brocket deer B chromosomes were found to be more homogeneous and resembled autosomes in patterns of sequence variation. Similar tendencies were observed in repetitive DNA composition. CONCLUSIONS: Our data demonstrate independent origins of B chromosomes in the grey brocket and roe deer. We hypothesize that the B chromosomes of these two cervid species represent different stages of B chromosome sequences evolution: probably nascent and similar to autosomal copies in brocket deer, highly derived in roe deer. Based on the presence of the same orthologous protooncogenes in canids and brocket deer Bs we argue that genomic regions involved in B chromosome formation are not random. In addition, our approach is also applicable to the characterization of other evolutionary and clinical rearrangements.

Evolutionary dynamics of Anolis sex chromosomes revealed by sequencing of flow sorting-derived microchromosome-specific DNA.

Squamate reptiles show a striking diversity in modes of sex determination, including both genetic (XY or ZW) and temperature-dependent sex determination systems. The genomes of only a handful of species have been sequenced, analyzed and assembled including the genome of Anolis carolinensis. Despite a high genome coverage, only macrochromosomes of A. carolinensis were assembled whereas the content of most microchromosomes remained unclear. Most of the Anolis species have homomorphic XY sex chromosome system. However, some species have large heteromorphic XY chromosomes (e.g., A. sagrei) and even multiple sex chromosomes systems (e.g. A. pogus), that were shown to be derived from fusions of the ancestral XY with microautosomes. We applied next generation sequencing of flow sorting-derived chromosome-specific DNA pools to characterize the content and composition of microchromosomes in A. carolinensis and A. sagrei. Comparative analysis of sequenced chromosome-specific DNA pools revealed that the A. sagrei XY sex chromosomes contain regions homologous to several microautosomes of A. carolinensis. We suggest that the sex chromosomes of A. sagrei are derived by fusions of the ancestral sex chromosome with three microautosomes and subsequent loss of some genetic content on the Y chromosome.

Comparative Chromosome Map and Heterochromatin Features of the Gray Whale Karyotype (Cetacea).

Cetacean karyotypes possess exceptionally stable diploid numbers and highly conserved chromosomes. To date, only toothed whales (Odontoceti) have been analyzed by comparative chromosome painting. Here, we studied the karyotype of a representative of baleen whales, the gray whale (Eschrichtius robustus, Mysticeti), by Zoo-FISH with dromedary camel and human chromosome-specific probes. We confirmed a high degree of karyotype conservation and found an identical order of syntenic segments in both branches of cetaceans. Yet, whale chromosomes harbor variable heterochromatic regions constituting up to a third of the genome due to the presence of several types of repeats. To investigate the cause of this variability, several classes of repeated DNA sequences were mapped onto chromosomes of whale species from both Mysticeti and Odontoceti. We uncovered extensive intrapopulation variability in the size of heterochromatic blocks present in homologous chromosomes among 3 individuals of the gray whale by 2-step differential chromosome staining. We show that some of the heteromorphisms observed in the gray whale karyotype are due to distinct amplification of a complex of common cetacean repeat and heavy satellite repeat on homologous autosomes. Furthermore, we demonstrate localization of the telomeric repeat in the heterochromatin of both gray and pilot whale (Globicephala melas, Odontoceti). Heterochromatic blocks in the pilot whale represent a composite of telomeric and common repeats, while heavy satellite repeat is lacking in the toothed whale consistent with previous studies.

Chromosome Polymorphism in Microtus (Alexandromys) mujanensis (Arvicolinae, Rodentia).

The Muya Valley vole (Microtus mujanensis) has a constant diploid chromosome number of 2n = 38, but an unstable karyotype with polymorphic chromosome pairs. Here, we describe 4 karyotypic variants involving 2 polymorphic chromosome pairs, MMUJ8 and MMUJ14, in 6 animals from Buryatia using a combination of GTG-banding and chromosome painting with M. agrestis probes. We suggest that the polymorphic pairs MMUJ8 and MMUJ14 were formed through pericentric inversions that played a major role during karyotype evolution of the species. We also propose that the stable diploid number with some ongoing polymorphism in the number of chromosome arms indicates that this evolutionarily young endemic species of Russian Far East is on the way to karyotype and likely species stabilization.

A new case of an inherited reciprocal translocation in cattle: rcp(13;26) (q24;q11).

This study reports on a unique balanced reciprocal chromosome translocation detected in a phenotypically normal cattle dam and her calf. CBG-, GTG-banding and FISH using bovine whole-chromosome and telomere probes were applied. The analyses showed that the breakpoints were located near to the centromere in chromosome 26 (q11) and exceptionally close to the telomere in chromosome 13 (q24). The whole euchromatin segment of chromosome 26 was translocated onto chromosome 13. The distal end of chromosome 13 was translocated to the subcentromeric region of chromosome 26. We describe this aberration as a balanced reciprocal translocation rcp(13;26) (q24;q11). It appears that this aberration was maternally derived and may have originated de novo in the dam.