Search

Demographic history reconstruction is based on the estimation of effective population size (Ne), which is inferred and interpreted in various fields of evolutionary and conservation biology. Interest in Ne estimation is growing, as the key evolutionary forces and their are linked to Ne, and genetic data become increasingly accessible. However, what is effective population size, and how can we obtain an estimate of effective population size? In this review, we describe the history of the term Ne and explore existing methods for obtaining historical and contemporary estimates of changes in effective population size. We provide a detailed overview of methods based on sequential Markovian coalescence (SMC), generalized phylogenetic coalescence (G-PhoCS), identity by descent (IBD) and identity by state (IBS) similarity, as well as methods using allele frequency spectrum (AFS). For each method, we briefly summarize the underlying theory and note its advantages and disadvantages. In the final section of the review, we present examples of the use of these methods for various non-model species with conservation status.

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.

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.

DNA hybridization-capture techniques allow researchers to focus their sequencing efforts on preselected genomic regions. This feature is especially useful when analys- ing ancient DNA (aDNA) extracts, which are often dominated by exogenous environ- mental sources. Here, we assessed, for the first time, the performance of hyRAD as an inexpensive and design-free alternative to commercial capture protocols to obtain authentic aDNA data from osseous remains. HyRAD relies on double enzymatic re- striction of fresh DNA extracts to produce RNA probes that cover only a fraction ofthe genome and can serve as baits for capturing homologous fragments from aDNA li- braries. We found that this approach could retrieve sequence data from horse remains coming from a range of preservation environments, including beyond radiocarbon range, yielding up to 146.5-fold on-target enrichment for aDNA extracts showing ex- tremely low endogenous content (20%¿30%), while the fraction of endogenous reads mapping on- and off-target was relatively insensi- tive to the original endogenous DNA content. Procedures based on two instead of a single round of capture increased on-target coverage up to 3.6-fold. Additionally, we used methylation-sensitive restriction enzymes to produce probes targeting hypo- methylated regions, which improved data quality by reducing post-mortem DNA dam- age and mapping within multicopy regions. Finally, we developed a fully automated hyRAD protocol utilizing inexpensive robotic platforms to facilitate capture process- ing. Overall, our work establishes hyRAD as a cost-effective strategy to recover a set of shared orthologous variants across multiple ancient samples., This project received funding from: the University Paul Sabatier IDEX Chaire d’Excellence (OURASI); the CNRS Programme de Recherche Conjoint (PRC); the CNRS International Research Project (IRP AMADEUS); the European Union’s Horizon 2020 research and innovation programme under the Marie Skłodowska-Curie grant agreement No. 797449; the Russian Foundation for Basic Research, project No. 19-59-15001 “Horses and their importance in the life of the ancient population of Altai and adjacent territories: interdisciplinary research and reconstruction”; and the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation programme (grant agreement 681605).

The brown bear (Ursus arctos) is an iconic carnivoran species of the Northern Hemisphere. Its population history has been studied extensively using mitochondrial markers, which demonstrated signatures of multiple waves of migration, arguably connected with glaciation periods. Among Eurasian brown bears, Siberian populations remain understudied. We have sequenced complete mitochondrial genomes of four ancient (~4.5–40 kya) bears from South Siberia and 19 modern bears from South Siberia and the Russian Far East. Reconstruction of phylogenetic relationships between haplotypes and evaluation of modern population structure have demonstrated that all the studied samples belong to the most widespread Eurasian clade 3. One of the ancient haplotypes takes a basal position relative to the whole of clade 3; the second is basal to the haplogroup 3a (the most common subclade), and two others belong to clades 3a1 and 3b. Modern Siberian bears retain at least some of this diversity; apart from the most common haplogroup 3a, we demonstrate the presence of clade 3b, which was previously found mainly in mainland Eurasia and Northern Japan. Our findings highlight the importance of South Siberia as a refugium for northern Eurasian brown bears and further corroborate the hypothesis of several waves of migration in the Pleistocene.

The taxonomy of the genus Calomyscus remains controversial. According to the latest systematics the genus includes eight species with great karyotypic variation. Here, we studied karyotypes of 14 Calomyscus individuals from different regions of Iran and Turkmenistan using a new set of chromosome painting probes from a Calomyscus sp. male (2n = 46, XY; Shahr-e-Kord-Soreshjan-Cheshme Maiak Province). We showed the retention of large syntenic blocks in karyotypes of individuals with identical chromosome numbers. The only rearrangement (fusion 2/21) differentiated Calomyscus elburzensis, Calomyscus mystax mystax, and Calomyscus sp. from Isfahan Province with 2n = 44 from karyotypes of C. bailwardi, Calomyscus sp. from Shahr-e-Kord, Chahar Mahal and Bakhtiari-Aloni, and Khuzestan-Izeh Provinces with 2n = 46. The individuals from Shahdad tunnel, Kerman Province with 2n = 51–52 demonstrated non-centric fissions of chromosomes 4, 5, and 6 of the 46-chromosomal form with the formation of separate small acrocentrics. A heteromorphic pair of chromosomes in a specimen with 2n = 51 resulted from a fusion of two autosomes. C-banding and chromomycin A3-DAPI staining after G-banding showed extensive heterochromatin variation between individuals.

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.

The genus status of Urocricetus was defined recently based on morphological and molecular data. Even though the amount of evidence for a separate phylogenetic position of this genus among Cricetinae continues to increase, there is still no consensus on its relationship to other groups. Here we give the first comprehensive description of the U. kamensis karyotype (2n = 30, NFa = 50) including results of comparative cytogenetic analysis and detailed examination of its phylogenetic position by means of numerous molecular markers. The molecular data strongly indicated that Urocricetus is a distant sister group to Phodopus. Comparative cytogenetic data showed significant reorganization of the U. kamensis karyotype compared to karyotypes of all other hamsters investigated earlier. The totality of findings undoubtedly means that Urocricetus belongs to a separate divergent lineage of Cricetinae.

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.

Domestication of horses fundamentally transformed long-range mobility and warfare. However, modern domesticated breeds do not descend from the earliest domestic horse lineage associated with archaeological evidence of bridling, milking and corralling at Botai, Central Asia around 3500 bc3. Other longstanding candidate regions for horse domestication, such as Iberia and Anatolia, have also recently been challenged. Thus, the genetic, geographic and temporal origins of modern domestic horses have remained unknown. Here we pinpoint the Western Eurasian steppes, especially the lower Volga-Don region, as the homeland of modern domestic horses. Furthermore, we map the population changes accompanying domestication from 273 ancient horse genomes. This reveals that modern domestic horses ultimately replaced almost all other local populations as they expanded rapidly across Eurasia from about 2000 bc, synchronously with equestrian material culture, including Sintashta spoke-wheeled chariots. We find that equestrianism involved strong selection for critical locomotor and behavioural adaptations at the GSDMC and ZFPM1 genes. Our results reject the commonly held association between horseback riding and the massive expansion of Yamnaya steppe pastoralists into Europe around 3000 bc driving the spread of Indo-European languages. This contrasts with the scenario in Asia where Indo-Iranian languages, chariots and horses spread together, following the early second millennium bc Sintashta culture., The work by G. Boeskorov is done on state assignment of DPMGI SB RAS. This project was supported by the University Paul Sabatier IDEX Chaire d’Excellence (OURASI); Villum Funden miGENEPI research programme; the CNRS ‘Programme de Recherche Conjoint’ (PRC); the CNRS International Research Project (IRP AMADEUS); the France Génomique Appel à Grand Projet (ANR-10-INBS-09-08, BUCEPHALE project); IB10131 and IB18060, both funded by Junta de Extremadura (Spain) and European Regional Development Fund; Czech Academy of Sciences (RVO:67985912); the Zoological Institute ZIN RAS (АААА-А19-119032590102-7); and King Saud University Researchers Supporting Project (NSRSP–2020/2). The research was carried out with the financial support of the Russian Foundation for Basic Research (19-59-15001 and 20-04-00213), the Russian Science Foundation (16-18-10265, 20-78-10151, and 21-18-00457), the Government of the Russian Federation (FENU-2020-0021), the Estonian Research Council (PRG29), the Estonian Ministry of Education and Research (PRG1209), the Hungarian Scientific Research Fund (Project NF 104792), the Hungarian Academy of Sciences (Momentum Mobility Research Project of the Institute of Archaeology, Research Centre for the Humanities); and the Polish National Science Centre (2013/11/B/HS3/03822). This project has received funding from the European Union’s Horizon 2020 research and innovation programme under the Marie Skłodowska-Curie (grant agreement 797449). This project has received funding from the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation programme (grant agreements 681605, 716732 and 834616).

Pinnipedia karyotype evolution was studied here using human, domestic dog, and stone marten whole-chromosome painting probes to obtain comparative chromosome maps among species of Odobenidae (Odobenus rosmarus), Phocidae (Phoca vitulina, Phoca largha, Phoca hispida, Pusa sibirica, Erignathus barbatus), and Otariidae (Eumetopias jubatus, Callorhinus ursinus, Phocarctos hookeri, and Arctocephalus forsteri). Structural and functional chromosomal features were assessed with telomere repeat and ribosomal-DNA probes and by CBG (C-bands revealed by barium hydroxide treatment followed by Giemsa staining) and CDAG (Chromomycin A3-DAPI after G-banding) methods. We demonstrated diversity of heterochromatin among pinniped karyotypes in terms of localization, size, and nucleotide composition. For the first time, an intrachromosomal rearrangement common for Otariidae and Odobenidae was revealed. We postulate that the order of evolutionarily conserved segments in the analyzed pinnipeds is the same as the order proposed for the ancestral Carnivora karyotype (2n = 38). The evolution of conserved genomes of pinnipeds has been accompanied by few fusion events (less than one rearrangement per 10 million years) and by novel intrachromosomal changes including the emergence of new centromeres and pericentric inversion/centromere repositioning. The observed interspecific diversity of pinniped karyotypes driven by constitutive heterochromatin variation likely has played an important role in karyotype evolution of pinnipeds, thereby contributing to the differences of pinnipeds&rsquo, chromosome sets.

A growing number of researchers studying horse domestication come to a conclusion that this process happened in multiple locations and involved multiple wild maternal lines. The most promising approach to address this problem involves mitochondrial haplotype comparison of wild and domestic horses from various locations coupled with studies of possible migration routes of the ancient shepherds. Here, we sequenced complete mitochondrial genomes of six horses from burials of the Ukok plateau (Russia, Altai Mountains) dated from 2.7 to 1.4 thousand years before present and a single late Pleistocene wild horse from the neighboring region (Denisova cave). Sequencing data indicates that the wild horse belongs to an extinct pre-domestication lineage. Integration of the domestic horse data with known Eurasian haplotypes of a similar age revealed two distinct groups: the first one widely distributed in Europe and presumably imported to Altai, and the second one specific for Altai Mountains and surrounding area.

Euchromatic segments of the X chromosomes of placental mammals are the most conservative elements of the karyotype, only rarely subjected to either inter- or intrachromosomal rearrangements. Here, using microdissection-derived set of region-specific probes of Terricola savii we detailed the evolutionary rearrangements found in X chromosomes in 20 vole species (Arvicolinae, Rodentia). We show that the evolution of X chromosomes in this taxon was accompanied by multiple para- and pericentric inversions and centromere shifts. The contribution of intrachromosomal rearrangements to the karyotype evolution of Arvicolinae species was approximately equivalent in both the separate autosomal conserved segments and the X chromosomes. Intrachromosmal rearrangements and structural reorganization of the X chromosomes was likely accompanied by an accumulation, distribution, and evolution of repeated sequences.

The mandarin vole, Lasiopodomys mandarinus, is one of the most intriguing species among mammals with non-XX/XY sex chromosome system. It combines polymorphism in diploid chromosome numbers, variation in the morphology of autosomes, heteromorphism of X chromosomes, and several sex chromosome systems the origin of which remains unexplained. Here we elucidate the sex determination system in Lasiopodomys mandarinus vinogradovi using extensive karyotyping, crossbreeding experiments, molecular cytogenetic methods, and single chromosome DNA sequencing. Among 205 karyotyped voles, one male and three female combinations of sex chromosomes were revealed. The chromosome segregation pattern and karyomorph-related reproductive performances suggested an aberrant sex determination with almost half of the females carrying neo-X/neo-Y combination. The comparative chromosome painting strongly supported this proposition and revealed the mandarin vole sex chromosome systems originated due to at least two de novo autosomal translocations onto the ancestral X chromosome. The polymorphism in autosome 2 was not related to sex chromosome variability and was proved to result from pericentric inversions. Sequencing of microdissection derived of sex chromosomes allowed the determination of the coordinates for syntenic regions but did not reveal any Y-specific sequences. Several possible sex determination mechanisms as well as interpopulation karyological differences are discussed.

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.

Background: Intellectual disability (ID) is an important medical and social problem that can be caused by different genetic and environmental factors. One such factor could be rDNA amplification and changes in rRNA expression and maturation. Aim of the study: The aim of the present study was to investigate rRNA levels in patients with heteromorphism of the p-arms of acrocentric chromosomes bearing nucleolus organizer regions compared to a healthy control group. Material and methods: Frequencies of p-arms enlargements in patients with ID and in healthy people were analyzed by G-banding screening. rRNA gene copy numbers on affected acrocentric chromosomes in peripheral blood lymphocytes were evaluated in ID patients and healthy bearers using FISH, and in immortalized lymphocytes of one patient – using FISH and real time PCR. Simultaneously, levels of 18S, 28S and 5,8S rRNA in both groups by means of qRT-PCR were investigated. Results: No difference in acrocentric chromosome heteromorphism frequency in patients versus the healthy group were found. However, we found an amplification of rDNA, a significant elevation in 28S and 5.8S rRNA expression and changes in the 28S/18S rRNA ratio in ID patients compared to healthy controls. At the same time, FISH appeared to be not reliable enough for copy number evaluation, but RT-PCR showed rDNA copy changes in heteromorphic cells compared to normal. Conclusion: Our findings indicate a loss of the correct regulation of rDNA activity and processing after amplification. This could disturb the ribosomal apparatus and thus lead to intellectual disability via at least two mechanisms. Keywords: Intellectual disability, Acrocentric chromosomes, Nucleolus organizer region – NOR, Ribosomal DNA – rDNA, Ribosomal RNA – rRNA

The taxonomy of the genus Calomyscus remains controversial. According to the latest systematics the genus includes eight species with great karyotypic variation. Here, we studied karyotypes of 14 Calomyscus individuals from different regions of Iran and Turkmenistan using a new set of chromosome painting probes from a Calomyscus sp. male (2n = 46, XY, Shahr-e-Kord-Soreshjan-Cheshme Maiak Province). We showed the retention of large syntenic blocks in karyotypes of individuals with identical chromosome numbers. The only rearrangement (fusion 2/21) differentiated Calomyscus elburzensis, Calomyscus mystax mystax, and Calomyscus sp. from Isfahan Province with 2n = 44 from karyotypes of C. bailwardi, Calomyscus sp. from Shahr-e-Kord, Chahar Mahal and Bakhtiari-Aloni, and Khuzestan-Izeh Provinces with 2n = 46. The individuals from Shahdad tunnel, Kerman Province with 2n = 51–52 demonstrated non-centric fissions of chromosomes 4, 5, and 6 of the 46-chromosomal form with the formation of separate small acrocentrics. A heteromorphic pair of chromosomes in a specimen with 2n = 51 resulted from a fusion of two autosomes. C-banding and chromomycin A3-DAPI staining after G-banding showed extensive heterochromatin variation between individuals.

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&rsquo, 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.

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.

The confluence of two scientific disciplines may lead to nomenclature conflicts that require new terms while respecting historical definitions. This is the situation with the current state of cytology and genomics, which offer examples of distinct nomenclature and vocabularies that require reconciliation. In this article, we propose the new terms C-scaffold (for chromosome-scale assemblies of sequenced DNA fragments, commonly named scaffolds) and scaffotype (the resulting collection of C-scaffolds that represent an organism's genome). This nomenclature avoids conflict with the historical definitions of the terms chromosome (a microscopic body made of DNA and protein) and karyotype (the collection of images of all chromosomes of an organism or species). As large-scale sequencing projects progress, adoption of this nomenclature will assist end users to properly classify genome assemblies, thus facilitating genomic analysis.

BackgroundThe Masai giraffe (Giraffa camelopardalis tippelskirchi) is the largest-bodied giraffe and the world's tallest terrestrial animal. With its extreme size and height, the giraffe's unique anatomical and physiological adaptations have long been of interest to diverse research fields. Giraffes are also critical to ecosystems of sub-Saharan Africa, with their long neck serving as a conduit to food sources not shared by other herbivores. Although the genome of a Masai giraffe has been sequenced, the assembly was highly fragmented and suboptimal for genome analysis. Herein we report an improved giraffe genome assembly to facilitate evolutionary analysis of the giraffe and other ruminant genomes.FindingsUsing SOAPdenovo2 and 170 Gbp of Illumina paired-end and mate-pair reads, we generated a 2.6-Gbp male Masai giraffe genome assembly, with a scaffold N50 of 3 Mbp. The incorporation of 114.6 Gbp of Chicago library sequencing data resulted in a HiRise SOAPdenovo + Chicago assembly with an N50 of 48 Mbp and containing 95% of expected genes according to BUSCO analysis. Using the Reference-Assisted Chromosome Assembly tool, we were able to order and orient scaffolds into 42 predicted chromosome fragments (PCFs). Using fluorescence in situ hybridization, we placed 153 cattle bacterial artificial chromosomes onto giraffe metaphase spreads to assess and assign the PCFs on 14 giraffe autosomes and the X chromosome resulting in the final assembly with an N50 of 177.94 Mbp. In this assembly, 21,621 protein-coding genes were identified using both de novo and homology-based predictions.ConclusionsWe have produced the first chromosome-scale genome assembly for a Giraffidae species. This assembly provides a valuable resource for the study of artiodactyl evolution and for understanding the molecular basis of the unique adaptive traits of giraffes. In addition, the assembly will provide a powerful resource to assist conservation efforts of Masai giraffe, whose population size has declined by 52% in recent years.

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.

The sterlet (Acipenser ruthenus Linnaeus, 1758) is a relatively small sturgeon widely distributed in Eurasian rivers from the Danube to the Yenisei. During the twentieth century, all wild sterlet populations have declined due to anthropogenic factors including: overfishing, poaching, construction of dams, and pollution. Despite the necessity of characterization both wild and captive stocks, few studies of population genetics have been performed thus far. Here we studied the genetic diversity and geographic structure of sterlet populations across the eastern range – Ob-Irtysh and Yenisei basins – by sequencing a 628-bp fragment of mitochondrial DNA control region. We identified 98 new haplotypes, delineated 12 haplogroups and estimated the time of basal haplogroup divergence within the species as over 8 million years ago. Our data suggest that Ob-Irtysh and Yenisei populations are isolated from each other and much lower genetic diversity is present in the Yenisei population than in the Ob-Irtysh population. Our data imply that sterlet populations in Siberian rivers underwent bottleneck or fragmentation, followed by subsequent population expansion. The data obtained here are important for sterlet population monitoring and restocking management.

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.

Sussemionus is an extinct subgenus of Equus first characterized and delineated in 2010. The almost complete mitochondrial genome is available only for a single specimen of Sussemionus – a 40,000 years old E. ovodovi from Proskuryakova cave (Khakassia, Russia). Our studies of ancient horses from Denisova cave (Altai, Russia) revealed mitochondrial DNA of this species in a 32,000 years old sample. Using alignments to multiple mitochondrial genomes of non-caballine equids, we recovered 100% complete mitochondrial genome of E. ovodovi for the first time. Phylogenetic analysis demonstrates close relationship between this individual and the one previously described in Khakassia.

Remarkably stable genomic chromosome elements (evolutionary conserved segments or syntenies) are the basis of large-scale chromosome architecture in vertebrate species. However, these syntenic elements harbour evolutionary important changes through intrachromosomal rearrangements such as inversions and centromere repositioning. Here, using FISH with a set of 20 region-specific probes on a wide array of 28 species, we analyzed evolution of three conserved syntenic regions of the Arvicolinae ancestral karyotype. Inside these syntenies we uncovered multiple, previously cryptic intrachromosomal rearrangements. Although in each of the three conserved blocks we found inversions and centromere repositions, the blocks experienced different types of rearrangements. In two syntenies centromere repositioning predominated, while in the third region, paracentric inversions were more frequent, whereas pericentric inversions were not detected. We found that some of the intrachromosomal rearrangements, mainly paracentric inversions, were synapomorphic for whole arvicoline genera or tribes: genera Alexandromys and Microtus, tribes Ellobini and Myodini. We hypothesize that intrachromosomal rearrangements within conserved syntenic blocks are a major evolutionary force modulating genome architecture in species-rich and rapidly-evolving rodent taxa. Inversions and centromere repositioning may impact speciation and provide a potential link between genome evolution, speciation, and biogeography.

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

There are differences in number and localization of nucleolus organizer regions (NORs) in genomes. In mammalian genomes, NORs are located on autosomes, which are often situated on short arms of acrocentric chromosomes and more rarely in telomeric, pericentromeric, or interstitial regions. In this work, we report the unique case of active NORs located on gonоsomes of a eutherian mammal, the Javan mouse-deer (Tragulus javanicus). We have investigated the position of NORs by FISH experiments with ribosomal DNA (rDNA) sequences (18S, 5.8S, and 28S) and show the presence of a single NOR site on the X and Y chromosomes. The NOR is localized interstitially on the p-arm of the X chromosome in close proximity with prominent C-positive heterochromatin blocks and in the pericentromeric area of mostly heterochromatic Y. The NOR sites are active on both the X and Y chromosomes in the studied individual and surrounded by GC enriched heterochromatin. We hypothesize that the surrounding heterochromatin might have played a role in the transfer of NORs from autosomes to sex chromosomes during the karyotype evolution of the Javan mouse-deer.

Naked mole rat (NMR) is the long-lived and tumor-resistant rodent. NMRs possess multiple adaptations that may contribute to longevity and cancer-resistance. However, whether NMRs have more efficient DNA repair have not been directly tested. Here we compared base excision repair (BER) and nucleotide excision repair (NER) systems in extracts from NMR and mouse fibroblasts after UVC irradiation. Transcript levels of the key repair enzymes demonstrated in most cases higher inducibility in the mouse vs the NMR cells. Ratios of repair enzymes activities in the extracts somewhat varied depending on post-irradiation time. NMR cell extracts were 2-3-fold more efficient at removing the bulky lesions, 1.5-3-fold more efficient at removing uracil, and about 1.4-fold more efficient at cleaving the AP-site than the mouse cells, while DNA polymerase activities being as a whole higher in the mouse demonstrate different patterns of product distribution. The level of poly(ADP-ribose) synthesis was 1.4-1.8-fold higher in the NMR cells. Furthermore, NMR cell extracts displayed higher binding of PARP1 to DNA probes containing apurinic/apyrimidinic site or photo-reactive DNA lesions. Cumulatively, our results suggest that the NMR has more efficient excision repair systems than the mouse, which may contribute to longevity and cancer resistance of this species.

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.

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. Electronic supplementary material The online version of this article (doi: 10.1186/s12862-017-1103-z) contains supplementary material, which is available to authorized users.

Several whole genome duplication (WGD) events followed by rediploidization took place in the evolutionary history of vertebrates. Acipenserids represent a convenient model group for investigation of the consequences of WGD as their representatives underwent additional WGD events in different lineages resulting in ploidy level variation between species, and these processes are still ongoing. Earlier, we obtained a set of sterlet (Acipenser ruthenus) chromosome-specific libraries by microdissection and revealed that they painted two or four pairs of whole sterlet chromosomes, as well as additional chromosomal regions, depending on rediploidization status and chromosomal rearrangements after genome duplication. In this study, we employed next generation sequencing to estimate the content of libraries derived from different paralogous chromosomes of sterlet. For this purpose, we aligned the obtained reads to the spotted gar (Lepisosteus oculatus) reference genome to reveal syntenic regions between these two species having diverged 360 Mya. We also showed that the approach is effective for synteny prediction at various evolutionary distances and allows one to clearly distinguish paralogous chromosomes in polyploid genomes. We postulated that after the acipenserid-specific WGD sterlet karyotype underwent multiple interchromosomal rearrangements, but different chromosomes were involved in this process unequally.

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.

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.

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.

Five families are traditionally recognized within higher ruminants (Pecora): Bovidae, Moschidae, Cervidae, Giraffidae and Antilocapridae. The phylogenetic relationships of Antilocapridae and Giraffidae within Pecora are, however, uncertain. While numerous fusions (mostly Robertsonian) have accumulated in the giraffe's karyotype (Giraffa camelopardalis, Giraffidae, 2n = 30), that of the pronghorn (Antilocapra americana, Antilocapridae, 2n = 58) is very similar to the hypothesised pecoran ancestral state (2n = 58). We examined the chromosomal rearrangements of two species, the giraffe and pronghorn, using a combination of fluorescence in situ hybridization painting probes and BAC clones derived from cattle (Bos taurus, Bovidae). Our data place Moschus (Moschidae) closer to Bovidae than Cervidae. Although the alternative (i.e., Moschidae + Cervidae as sister groups) could not be discounted in recent sequence-based analyses, cytogenetics bolsters conclusions that the former is more likely. Additionally, DNA sequences were isolated from the centromeric regions of both species and compared. Analysis of cenDNA show that unlike the pronghorn, the centromeres of the giraffe are probably organized in a more complex fashion comprising different repetitive sequences specific to single chromosomal pairs or groups of chromosomes. The distribution of nucleolar organiser region (NOR) sites, often an effective phylogenetic marker, were also examined in the two species. In the giraffe, the position of NORs seems to be autapomorphic since similar localizations have not been found in other species within Pecora.

Chromosome painting and G-banding analyses were used to delimit homologous chromosomal segments among 4 taxa of the hamsters genus Allocricetulus Argyropulo 1932 (Cricetidae, Murоidea, Rodentia) - A. curtatus (2n = 20), A. eversmanni eversmanni, A. eversmanni pseudocurtatus, and the hybrid A. eversmanni beljaevi × A. eversmanni pseudocurtatus (all 2n = 26). Comparative maps between the 4 karyotypes were established based on chromosome painting of chromosome-specific probes from the Syrian hamster (Mesocricetus auratus, 2n = 44). A putative ancestral karyotype for the genus Allocricetulus (AAK) was proposed and contains 12-13 ancestral autosomal elements. Integrated maps demonstrate extended conservation of syntenies within this rodent genus and show the predominant role of Robertsonian rearrangements in the karyotype evolution of the genus Allocricetulus. At the cytogenetic level, we clearly demonstrate karyological differences between karyotypes of species (A. curtatus vs. A. eversmanni) and subspecies A. e. eversmanni and A. e. beljaevi versus A. e. pseudocurtatus, but the karyotypes of A. e. eversmanni and A. e. beljaevi are identical at this level of resolution.

Zoo-FISH (fluorescence in situ hybridization) (Scherthan et al. 1994), also known as cross-species chromosome painting or comparative chromosome painting, essentially describes the use of whole chromosome- or chromosomal arm- or region-specific painting probes to delimit homologous segments (chromosome or chromosomal segments with evolutionarily conserved synteny) in other species by means of FISH. Like most conventional cytogenetic methodologies, the technique of chromosome painting was first introduced to human cytogenetics in 1988 (Pinkel et al. 1988). Wienberg and colleagues (1990) were among the pioneers who intro-duced chromosomal painting to the field of comparative cytogenetics of primates. They established the first genome-wide chromosome maps between human and a Japanese macaque. In the beginning, painting probes were derived from libraries of human chromosome-specific DNA clones. Limited by the availability of painting probes and technical difficulties in comparing distantly related species, early cross-species chromosome painting experiments primarily used apes and Old World monkeys (Jauch et al. 1992). Scherthan and colleagues (1994) were among the first to demonstrate the feasibility of comparing species as distantly related as those from different orders. Meanwhile, the invention of degenerate oligonucleotide-primed PCR (DOP-PCR, Telenius et al. 1992a), coupled with chromosomal sorting by flow cytometry, made it possible to generate painting probes for any given vertebrate species (Telenius et al 1992b; Rabbitts et al. 1995; Yang et al. 1995) and to carry out multidirectional cross-species chromosome painting. The whole set of human chromosome-specific painting probes derived from DOP-PCR were made commercially available around 1992, and the whole set of chromosome painting probes for mouse became available in late 1995. Human chromosome-specific

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. Electronic supplementary material The online version of this article (doi:10.1186/s12864-016-2933-6) contains supplementary material, which is available to authorized users.

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.