Chromosomes of Damaliscus (Artiodactyla, Bovidae): simple and complex centric fusion rearrangements

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.

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.

Joint identification of sex and sex-linked scaffolds in non-model organisms using low depth sequencing data

Being able to assign sex to individuals and identify autosomal and sex-linked scaffolds are essential in most population genomic analyses. Non-model organisms often have genome assemblies at scaffold-level and lack characterization of sex-linked scaffolds. Previous methods to identify sex and sex-linked scaffolds have relied on synteny between the non-model organism and a closely related species or prior knowledge about the sex of the samples to identify sex-linked scaffolds. In the latter case, the difference in depth of coverage between the autosomes and the sex chromosomes are used. Here, we present "sex assignment through coverage" (SATC), a method to assign sex to samples and identify sex-linked scaffolds from next generation sequencing (NGS) data. The method works for species with a homogametic/heterogametic sex determination system and only requires a scaffold-level reference assembly and sampling of both sexes with whole genome sequencing (WGS) data. We use the sequencing depth distribution across scaffolds to jointly identify: (i) male and female individuals, and (ii) sex-linked scaffolds. This is achieved through projecting the scaffold depths into a low-dimensional space using principal component analysis (PCA) and subsequent Gaussian mixture clustering. We demonstrate the applicability of our method using data from five mammal species and a bird species complex. The method is freely available at https://github.com/popgenDK/SATC as R code and a graphical user interface (GUI).

Genetic diversity in natural range remnants of the critically endangered hirola antelope

The hirola antelope (Beatragus hunteri) is considered to be the most endangered antelope in the world. In the ex situ translocated population at Tsavo East National Park, calf mortality and the critically low population numbers might suggest low genetic diversity and inbreeding depression. Consequently, a genetic study of the wild population is pivotal to gain an understanding of diversity and differentiation within its range before designing future translocation plans to increase the genetic diversity of the ex situ population. For that purpose, we assessed 55 individuals collected across five localities in eastern Kenya, covering its entire natural range. We used the complete mitochondrial DNA control region and microsatellite genotyping to estimate genetic diversity and differentiation across its range. Nuclear genetic diversity was moderate in comparison to other endangered African antelopes, with no signals of inbreeding. However, the mitochondrial data showed low nucleotide diversity, few haplotypes and low haplotypic differentiation. Overall, the inferred low degree of genetic differentiation and population structure suggests a single population of hirola across the natural range. An overall stable population size was inferred over the recent history of the species, although signals of a recent genetic bottleneck were found. Our results show hope for ongoing conservation management programmes and that there is a future for the hirola in Kenya.

Multi-locus phylogeny of the tribe Tragelaphini (Mammalia, Bovidae) and species delimitation in bushbuck: Evidence for chromosomal speciation mediated by interspecific hybridization

The bushbuck is the most widespread bovid species in Africa. Previous mitochondrial studies have revealed a polyphyletic pattern suggesting the possible existence of two distinct species. To assess this issue, we have sequenced 16 nuclear genes and one mitochondrial fragment (cytochrome b gene + control region) for most species of the tribe Tragelaphini, including seven bushbuck individuals belonging to the two divergent mtDNA haplogroups, Scriptus and Sylvaticus. Our phylogenetic analyses show that the Scriptus lineage is a sister-group of Sylvaticus in the nuclear tree, whereas it is related to Tragelaphus angasii in the mitochondrial tree. This mito-nuclear discordance indicates that the mitochondrial genome of Scriptus was acquired by introgression after one or several past events of hybridization between bushbuck and an extinct species closely related to T. angasii. The division into two bushbuck species is supported by the analyses of nuclear markers and by the karyotype here described for T. scriptus (2n = 57 M/58F), which is strikingly distinct from the one previously found for T. sylvaticus (2n = 33 M/34F). Molecular dating estimates suggest that the two species separated during the Early Pleistocene after an event of interspecific hybridization, which may have mediated massive chromosomal rearrangements in the common ancestor of T. scriptus.

Molecular phylogeny and chromosomal evolution of Alcelaphini (Antilopinae)

Robertsonian (Rb) translocations, in particular centric fusions, are thought to play a primary role in evolution and speciation of the Bovidae family. However, Rb fusions are often polymorphic within species, being suggested as phylogenetically uninformative characters. This work studies chromosome variation in 72 captive individuals of 6 species of Alcelaphini (Antilopinae): The hartebeest (genus Alcelaphus), hirola (Beatragus), black and blue wildebeests (Connochaetes), and the topi and bontebok (Damaliscus). We infer the phylogenic relationships among Alcelaphini species and determine patterns of chromosomal evolution using G-banded karyotypes and complete mitochondrial genome sequences. The molecular phylogeny showed an early divergence of Connochaetes, followed by the split of Alcelaphus plus Beatragus + Damaliscus as sister taxa. Mitochondrial and chromosomal phylogenies only differed in the position of the critically endangered Beatragus, likely due to homoplasic chromosome characters. Patterns of chromosome evolution, reconstructed using a probabilistic approach, suggest that chromosome changes leading to speciation in Alcelaphini do not exclusively involve consecutive reduction of diploid number through centric fusion but also the losses and reversions of Rb translocations in Beatragus and Damaliscus lineages. Our results provide evidence that complex scenarios of chromosomal rearrangements can be detected in relatively recent-diverged bovids, as in this group of antelopes.

Chromosome homologies between tsessebe (Damaliscus lunatus) and Chinese muntjac (Muntiacus reevesi) facilitate tracing the evolutionary history of Damaliscus (Bovidae, Antilopinae, Alcelaphini)

Genome-wide homologies between the tsessebe (Damaliscus lunatus, 2n = 36) and Chinese muntjac (Muntiacus reevesi, 2n = 46) have been established by cross-species painting with Chinese muntjac chromosome paints. Twenty-two autosomal painting probes detected 35 orthologous segments in the tsessebe. Hybridization results confirmed that: (i) D. lunatus carries the (9;14) reciprocal translocation that has been proposed to be a derived chromosomal landmark shared by all species of the Antilopinae; (ii) the karyotype of D. lunatus can be derived almost exclusively from the bovid ancestral karyotype through 12 Robertsonian translocations involving 24 ancestral acrocentric autosomes; (iii) in addition to the Rb fusions, pericentric heterochromatic amplification has shaped the morphology of several of the D. lunatus chromosomes. Integrated analysis of these and published cytogenetic data on pecorans has allowed us to accurately discern the karyotype history of Damaliscus (D. lunatus; D. pygargus, 2n = 38; D. hunteri, 2n = 44). The phylogenomic relationships of 3 species reflected by specific chromosomal rearrangements were consistent with published phylogenies based on morphology, suggesting that chromosomal rearrangements have played an important role in speciation within the Alcelaphini, and that karyotype characters are valuable phylogenetic markers in this group.

Species identification and chromosome variation of captive two-toed sloths

Two-toed sloth species, Linnaeus's and Hoffmman's, are frequent residents of zoo collections in North America. However, species identification has always been problematic because of their large overlap in external morphology, which represents an obstacle to the captive breeding program. We describe here a PCR-based technique that allows species identification of two-toed sloths without requiring sequencing, by using a mitochondrial marker (COI gene) and restriction enzyme assay. We also report intra- and inter-specific patterns of chromosome variation in captive two-toed sloths. Molecularly, we identified 22 samples of Linnaeus's and Hoffmman's two-toed sloths corresponding to 14 and 8 individuals, respectively. One animal was identified as a hybrid using the nuclear gene Enam having alleles derived from both species. The chromosome number in Hoffman's two-toed sloths showed low variation ranging only between 50 and 51. In contrast, Linnaeus's two-toed sloths appeared to vary widely, with diploid numbers ranging from 53 to 67, suggesting distinct geographic groups. The species identification method presented here represents a low-cost easy-to-use tool that will help to improve management of the captive population of two-toed sloths.

Multidirectional chromosome painting between the Hirola antelope (Damaliscus hunteri, Alcelaphini, Bovidae), sheep and human

Chromosome specific painting probes of human, sheep and the Hirola antelope ( Damaliscus hunteri ) derived by flow sorting of chromosomes were used in multi directional chromosome painting experiments to better define the karyological relationship within Bovidae species (specifically, Caprini and Alcelaphini tribes) and humans. Although not all chromosomes of Damaliscus hunteri could be resolved into single peaks by flow-sorting we managed to present a complete homology map for chromosomes between the three species. When comparing the karyotype of Damaliscus hunteri with human all of the main known motives in mammalian chromosome evolution are present (i.e. associations of human homologous chromosomes 3-21, 4-8, 7-16, 14-15, 16-19 and two forms of 12-22) which were also confirmed with the sheep paint probes. Further, we observed those patterns that have been described as common derived traits for artiodactyls (i.e. associations of human homologous chromosomes 5/19 and a complex alternating pattern of hybridizations with human chromosome 14 and 15 probes). As known from classical karyotyping some of the Damaliscus chromosomes are biarmed and were supposedly involved in Robertsonian translocations frequently found in karyotype evolution of bovids. We refined these rearrangements with the molecular probes and also delineated a chromosome painting pattern that should be the result of a paracentric inversion in the Damaliscus hunteri karyotype. This study demonstrates that multidirectional chromosome painting will be a valuable tool for the investigation of the dynamics of chromosome evolution in exotic bovid species.

Centromeric heterochromatin in the cattle rob(1;29) translocation: alpha-satellite I sequences, in-situ MspI digestion patterns, chromomycin staining and C-bands

The centromeric regions and alpha-satellite I sequence were studied on chromosomes 1, 29 and the rob(1;29) translocation in a Portuguese breed of cattle, Barrosa, carrying the translocation. Rob(1;29) centromeric regions showed heterochromatic bands with propidium iodide but, unlike the acrocentric autosomes, no strong centromeric bands were revealed with chromomycin A3. An alpha-satellite I sequence was not found at the centromeres of the X, Y and rob(1;29) chromosomes in the breed, although it was present at the centromeres of all acrocentric chromosomes including 1 and 29. Restriction enzyme banding with MspI revealed polymorphisms between different rob(1;29) chromosomes in both centromeric and intercalary regions. The data show that the centromeric region of the rob(1;29) chromosome has lost the alpha-satellite I sequences, while retaining other heterochromatin, and suggest that this important and widespread translocation has occurred multiple times.

A karyotypic analysis of nilgai, Boselaphus tragocamelus (Artiodactyla: Bovidae)

A combination of chromosomal banding and fluorescence in situ hybridization (FISH) was used to characterize the karyotype of Boselaphus tragocamelus (nilgai) relative to the domestic cattle standard karyotype. G-, Q- and C-band karyotypes of nilgai are presented, and the chromosomal complement of nilgai is determined to be 2n=46 (female FN=60, male FN=59; NAA=56), consistent with previous reports for the species. Comparisons with cattle identified extensive monobrachial homologies with some noteworthy exceptions. Chromosome 25 is centrically fused to 24, and chromosome 16 is acrocentric. Both appear to have additional pericentromeric material not seen in the equivalent cattle acrocentrics. This pericentromeric chromatin may be the result of de novo additions or translocation of pericentromeric material from chromosome 6, which is shown to be centrically fused to 13 but is only about two-thirds the length of cattle 6. Comparisons with cattle demonstrated that nilgai chromosome 17 has undergone a paracentric inversion and that chromosome 20 has two blocks of interstitial constitutive heterochromatin. The identities of both chromosomes were confirmed by chromosomal FISH. Furthermore, chromosomal banding and FISH were used to determine that autosome 14 has been fused to the ancestral X and Y of nilgai to form compound neo-X and -Y chromosomes. Additional FISH analyses were conducted to confirm other proposed chromosome homologies and to identify nucleolar organizing regions within the nilgai complement.