Karyotype evolution of giraffes (Giraffa camelopardalis) revealed by cross-species chromosome painting with Chinese muntjac (Muntiacus reevesi) and human (Homo sapiens) paints

A chromosome-scale genome assembly of the okapi (Okapia johnstoni)

The okapi (Okapia johnstoni), or forest giraffe, is the only species in its genus and the only extant sister group of the giraffe within the family Giraffidae. The species is one of the remaining large vertebrates surrounded by mystery because of its elusive behavior as well as the armed conflicts in the region where it occurs, making it difficult to study. Deforestation puts the okapi under constant anthropogenic pressure, and it is currently listed as “Endangered” on the IUCN Red List. Here, we present the first annotated de novo okapi genome assembly based on PacBio continuous long reads, polished with short reads, and anchored into chromosome-scale scaffolds using Hi-C proximity ligation sequencing. The final assembly (TBG_Okapi_asm_v1) has a length of 2.39 Gbp, of which 98% are represented by 28 scaffolds > 3.9 Mbp. The contig N50 of 61 Mbp and scaffold N50 of 102 Mbp, together with a BUSCO score of 94.7%, and 23 412 annotated genes, underline the high quality of the assembly. This chromosome-scale genome assembly is a valuable resource for future conservation of the species and comparative genomic studies among the giraffids and other ruminants.

A towering genome: Experimentally validated adaptations to high blood pressure and extreme stature in the giraffe

The suite of adaptations associated with the extreme stature of the giraffe has long interested biologists and physiologists. By generating a high-quality chromosome-level giraffe genome and a comprehensive comparison with other ruminant genomes, we identified a robust catalog of giraffe-specific mutations. These are primarily related to cardiovascular, bone growth, vision, hearing, and circadian functions. Among them, the giraffe FGFRL1 gene is an outlier with seven unique amino acid substitutions not found in any other ruminant. Gene-edited mice with the giraffe-type FGFRL1 show exceptional hypertension resistance and higher bone mineral density, both of which are tightly connected with giraffe adaptations to high stature. Our results facilitate a deeper understanding of the molecular mechanism underpinning distinct giraffe traits, and may provide insights into the study of hypertension in humans.

An integrated chromosome-scale genome assembly of the Masai giraffe (Giraffa camelopardalis tippelskirchi)

BACKGROUND

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

FINDINGS

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

CONCLUSIONS

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

Giraffe genome sequence reveals clues to its unique morphology and physiology

The origins of giraffe's imposing stature and associated cardiovascular adaptations are unknown. Okapi, which lacks these unique features, is giraffe's closest relative and provides a useful comparison, to identify genetic variation underlying giraffe's long neck and cardiovascular system. The genomes of giraffe and okapi were sequenced, and through comparative analyses genes and pathways were identified that exhibit unique genetic changes and likely contribute to giraffe's unique features. Some of these genes are in the HOX, NOTCH and FGF signalling pathways, which regulate both skeletal and cardiovascular development, suggesting that giraffe's stature and cardiovascular adaptations evolved in parallel through changes in a small number of genes. Mitochondrial metabolism and volatile fatty acids transport genes are also evolutionarily diverged in giraffe and may be related to its unusual diet that includes toxic plants. Unexpectedly, substantial evolutionary changes have occurred in giraffe and okapi in double-strand break repair and centrosome functions.

Giraffe's unique anatomy and physiology include its stature and associated cardiovascular adaptation. Here, Douglas Cavener and colleagues provide de novo genome assemblies of giraffe and its closest relative okapi and provide comparative analyses to infer insights into evolution and adaptation.

The binucleate cell of okapi and giraffe placenta shows distinctive glycosylation compared with other ruminants: a lectin histochemical study

The placenta of ruminants contains characteristic binucleate cells (BNC) with a highly conserved glycan structure which evolved early in Ruminant phylogenesis. Giraffe and Okapi placentae also contain these cells and it is not known whether they have a similar glycan array. We have used lectin histochemistry to examine the glycosylation of these cells in these species and compare them with bovine BNC which have a typical ruminant glycan composition. Two placentae, mid and near term, from Giraffe (Giraffa camelopardalis) and two term placenta of Okapi (Okapia johnstoni) were embedded in resin and stained with a panel of 23 lectins and compared with near-term bovine (Bos taurus) placenta. Significant differences were found in the glycans of Giraffe and Okapi BNC compared with those from the bovine, with little or no expression of terminal αN-acetylgalactosamine bound by Dolichos biflorus and Vicia villosa agglutinins which instead bound to placental blood vessels. Higher levels of N-acetylglucosamine bound by Lycopersicon esculentum and Phytolacca americana agglutinins were also apparent. Some differences between Okapi and Giraffe were evident. Most N-linked glycans were similarly expressed in all three species as were fucosyl residues. Interplacentomal areas in Giraffe and Bovine showed differences from the placentomal cells though no intercotyledonary BNC were apparent in Okapi. In conclusion, Giraffidae BNC developed different glycan biosynthetic pathways following their split from the Bovidae with further differences evolving as Okapi and Giraffe diverged from each other, affecting both inter and placentomal BNC which may have different functions during development.

Comparative chromosome painting of pronghorn (Antilocapra americana) and saola (Pseudoryx nghetinhensis) karyotypes with human and dromedary camel probes

BACKGROUND

Pronghorn (Antilocapridae, 2n = 58) and saola (Bovidae, 2n = 50) are members of Pecora, a highly diversified group of even-toed hoofed mammals. Karyotypes of these species were not involved in chromosome painting studies despite their intriguing phylogenetic positions in Pecora.

RESULTS

To trace the chromosome evolution during very fast radiation of main families from the common Pecoran ancestor, high-resolution comparative chromosome maps of pronghorn and saola with human (HSA) and dromedary camel (CDR) painting probes were established. The human and dromedary camel painting probes revealed 50 and 64 conserved segments respectively in the pronghorn genome, while 51 and 63 conserved segments respectively in the saola genome. Integrative analysis with published comparative maps showed that inversions in chromosomes homologous to CDR19/35/19 (HSA 10/20/10), CDR12/34/12 (HSA12/22/12/22), CDR10/33/10 (HSA 11) are present in representatives of all five living Pecoran families. The pronghorn karyotype could have formed from a putative 2n = 58 Pecoran ancestral karyotype by one fission and one fusion and that the saola karyotype differs from the presumed 2n = 60 bovid ancestral karyotype (2n = 60) by five fusions.

CONCLUSION

The establishment of high-resolution comparative maps for pronghorn and saola has shed some new insights into the putative ancestral karyotype, chromosomal evolution and phylogenic relationships in Pecora. No cytogenetic signature rearrangements were found that could unite the Antilocapridae with Giraffidae or with any other Pecoran families. Our data on the saola support a separate position of Pseudorigyna subtribe rather than its affinity to either Bovina or Bubalina, but the saola phylogenetic position within Bovidae remains unresolved.

The telomeric sync model of speciation: species-wide telomere erosion triggers cycles of transposon-mediated genomic rearrangements, which underlie the saltatory appearance of nonadaptive characters

Charles Darwin knew that the fossil record is not overwhelmingly supportive of genetic and phenotypic gradualism; therefore, he developed the core of his theory on the basis of breeding experiments. Here, I present evidence for the existence of a cell biological mechanism that strongly points to the almost forgotten European concept of saltatory evolution of nonadaptive characters, which is in perfect agreement with the gaps in the fossil record. The standard model of chromosomal evolution has always been handicapped by a paradox, namely, how speciation can occur by spontaneous chromosomal rearrangements that are known to decrease the fertility of heterozygotes in a population. However, the hallmark of almost all closely related species is a differing chromosome complement and therefore chromosomal rearrangements seem to be crucial for speciation. Telomeres, the caps of eukaryotic chromosomes, erode in somatic tissues during life, but have been thought to remain stable in the germline of a species. Recently, a large human study spanning three healthy generations clearly found a cumulative telomere effect, which is indicative of transgenerational telomere erosion in the human species. The telomeric sync model of speciation presented here is based on telomere erosion between generations, which leads to identical fusions of chromosomes and triggers a transposon-mediated genomic repatterning in the germline of many individuals of a species. The phenotypic outcome of the telomere-triggered transposon activity is the saltatory appearance of nonadaptive characters simultaneously in many individuals. Transgenerational telomere erosion is therefore the material basis of aging at the species level.

Molecular cytogenetic insights to the phylogenetic affinities of the giraffe (Giraffa camelopardalis) and pronghorn (Antilocapra americana)

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.

Comparative molecular cytogenetics in Cetartiodactyla

Cetartiodactyla comprises Artiodactyla (even-toed ungulates) and Cetacea (whales, dolphins and porpoises). Artiodactyla is a large taxon represented by about 200 living species ranked in 10 families. Cetacea are classified into 13 families with almost 80 species. Many publications concerning karyotypic relationships in Cetartiodactyla have been published in previous decades. Formerly, the karyotypes of closely related species were compared by chromosome banding. Introduction of molecular cytogenetic methods facilitated comparative mapping between species with highly rearranged karyotypes and distantly related species. Such information is a prerequisite for the understanding of karyotypic phylogeny and the reconstruction of the karyotypes of common ancestors. This study summarizes the data on chromosome evolution in Cetartiodactyla, mainly derived from molecular cytogenetic studies. Traditionally, phylogenetic relationships of most groups have been estimated using morphological data. However, the results of some molecular studies of mammalian phylogeny are discordant with traditional conceptions of phylogeny. Cetartiodactyls provide several examples of incongruence between traditional morphological and molecular data. Such cases of conflict include the relationships of the major clades of artiodactyls, the relationships among the extant families of the suborder Ruminantia or the phylogeny of the family Bovidae. The most unexpected aspect of the molecular phylogeny was the recognition that Cetacea is a deeply nested member of Artiodactyla. The largest living order of terrestrial hoofed mammals is the even-toed hoofed mammals, or Artiodactyla. The artiodactyls are composed of over 190 living species including pigs, peccaries, hippos, camels, llamas, deer, pronghorns, giraffes, sheep, goats, cattle and antelopes. Cetacea is an order of wholly aquatic mammals, which include whales, dolphins and porpoises. Cetartiodactyla has become the generally accepted name for the clade containing both of these orders.

Chromosome painting in three-toed sloths: a cytogenetic signature and ancestral karyotype for Xenarthra

BACKGROUND

Xenarthra (sloths, armadillos and anteaters) represent one of four currently recognized Eutherian mammal supraorders. Some phylogenomic studies point to the possibility of Xenarthra being at the base of the Eutherian tree, together or not with the supraorder Afrotheria. We performed painting with human autosomes and X-chromosome specific probes on metaphases of two three-toed sloths: Bradypus torquatus and B. variegatus. These species represent the fourth of the five extant Xenarthra families to be studied with this approach.

RESULTS

Eleven human chromosomes were conserved as one block in both B. torquatus and B. variegatus: (HSA 5, 6, 9, 11, 13, 14, 15, 17, 18, 20, 21 and the X chromosome). B. torquatus, three additional human chromosomes were conserved intact (HSA 1, 3 and 4). The remaining human chromosomes were represented by two or three segments on each sloth. Seven associations between human chromosomes were detected in the karyotypes of both B. torquatus and B. variegatus: HSA 3/21, 4/8, 7/10, 7/16, 12/22, 14/15 and 17/19. The ancestral Eutherian association 16/19 was not detected in the Bradypus species.

CONCLUSIONS

Our results together with previous reports enabled us to propose a hypothetical ancestral Xenarthran karyotype with 48 chromosomes that would differ from the proposed ancestral Eutherian karyotype by the presence of the association HSA 7/10 and by the split of HSA 8 into three blocks, instead of the two found in the Eutherian ancestor. These same chromosome features point to the monophyly of Xenarthra, making this the second supraorder of placental mammals to have a chromosome signature supporting its monophyly.

Chromosome painting in Tragulidae facilitates the reconstruction of Ruminantia ancestral karyotype

Although Tragulidae, as the basal family in Ruminantia phylogenetic tree, is the key taxon for understanding the early chromosome evolution of extant ruminants, comparative molecular cytogenetic data on the tragulids are scarce. Here, we present the first genome-wide comparative map of the Java mouse deer (Tragulus javanicus, Tragulidae) revealed by chromosome painting with human and dromedary probes. Together with the published comparative maps of major representative cetartiodactyl species established with the same set of probes, our results allowed us to reconstruct a 2n = 48 Ruminantia ancestral karyotype, which is similar to the cetartiodactyl ancestral karyotype. The karyotype evolution of T. javanicus has involved multiple rearrangements, most of which appear to be apomorphic and have not found in karyotype evolution of pecoran species (i.e., Ruminantia excluding Tragulidae). The rate of chromosome evolution of the mouse deer was rather low-0.4 R/Ma, while the estimated tempo of chromosome changes on the lineages leading from Cetartiodactyla ancestor to Ruminantia and from Ruminantia to Pecora were roughly the same (about 1.2 R/Ma).

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.

Cross-species chromosome painting in Cetartiodactyla: reconstructing the karyotype evolution in key phylogenetic lineages

Recent molecular and morphological studies place Artiodactyla and Cetacea into the order Cetartiodactyla. Within the Cetartiodactyla such families as Bovidae, Cervidae, and Suidae are well studied by comparative chromosome painting, but many taxa that are crucial for understanding cetartiodactyl phylogeny remain poorly studied. Here we present the genome-wide comparative maps of five cetartiodactyl species obtained by chromosome painting with human and dromedary paint probes from four taxa: Cetacea, Hippopotamidae, Giraffidae, and Moschidae. This is the first molecular cytogenetic report on pilot whale, hippopotamus, okapi, and Siberian musk deer. Our results, when integrated with previously published comparative chromosome maps allow us to reconstruct the evolutionary pathway and rates of chromosomal rearrangements in Cetartiodactyla. We hypothesize that the putative cetartiodactyl ancestral karyotype (CAK) contained 25-26 pairs of autosomes, 2n = 52-54, and that the association of human chromosomes 8/9 could be a cytogenetic signature that unites non-camelid cetartiodactyls. There are no unambiguous cytogenetic landmarks that unite Hippopotamidae and Cetacea. If we superimpose chromosome rearrangements on the supertree generated by Price and colleagues, several homoplasy events are needed to explain cetartiodactyl karyotype evolution. Our results apparently favour a model of non-random breakpoints in chromosome evolution. Cetariodactyl karyotype evolution is characterized by alternating periods of low and fast rates in various lineages. The highest rates are found in Suina (Suidae+Tayasuidae) lineage (1.76 rearrangements per million years (R/My)) and the lowest in Cetaceans (0.07 R/My). Our study demonstrates that the combined use of human and camel paints is highly informative for revealing evolutionary karyotypic rearrangements among cetartiodactyl species.