Alexey I. Makunin

Genome-wide Evidence Reveals that African and Eurasian Golden Jackals Are Distinct Species

The golden jackal of Africa (Canis aureus) has long been considered a conspecific of jackals distributed throughout Eurasia, with the nearest source populations in the Middle East. However, two recent reports found that mitochondrial haplotypes of some African golden jackals aligned more closely to gray wolves (Canis lupus), which is surprising given the absence of gray wolves in Africa and the phenotypic divergence between the two species. Moreover, these results imply the existence of a previously unrecognized phylogenetically distinct species despite a long history of taxonomic work on African canids. To test the distinct-species hypothesis and understand the evolutionary history that would account for this puzzling result, we analyzed extensive genomic data including mitochondrial genome sequences, sequences from 20 autosomal loci (17 introns and 3 exon segments), microsatellite loci, X- and Y-linked zinc-finger protein gene (ZFX and ZFY) sequences, and whole-genome nuclear sequences in African and Eurasian golden jackals and gray wolves. Our results provide consistent and robust evidence that populations of golden jackals from Africa and Eurasia represent distinct monophyletic lineages separated for more than one million years, sufficient to merit formal recognition as different species: C. anthus (African golden wolf) and C. aureus (Eurasian golden jackal). Using morphologic data, we demonstrate a striking morphologic similarity between East African and Eurasian golden jackals, suggesting parallelism, which may have misled taxonomists and likely reflects uniquely intense interspecific competition in the East African carnivore guild. Our study shows how ecology can confound taxonomy if interspecific competition constrains size diversification.

Genomic Legacy of the African Cheetah, Acinonyx jubatus

BackgroundPatterns of genetic and genomic variance are informative in inferring population history for human, model species and endangered populations.ResultsHere the genome sequence of wild-born African cheetahs reveals extreme genomic depletion in SNV incidence, SNV density, SNVs of coding genes, MHC class I and II genes, and mitochondrial DNA SNVs. Cheetah genomes are on average 95 % homozygous compared to the genomes of the outbred domestic cat (24.08 % homozygous), Virunga Mountain Gorilla (78.12 %), inbred Abyssinian cat (62.63 %), Tasmanian devil, domestic dog and other mammalian species. Demographic estimators impute two ancestral population bottlenecks: one >100,000 years ago coincident with cheetah migrations out of the Americas and into Eurasia and Africa, and a second 11,084–12,589 years ago in Africa coincident with late Pleistocene large mammal extinctions. MHC class I gene loss and dramatic reduction in functional diversity of MHC genes would explain why cheetahs ablate skin graft rejection among unrelated individuals. Significant excess of non-synonymous mutations in AKAP4 (p<0.02), a gene mediating spermatozoon development, indicates cheetah fixation of five function-damaging amino acid variants distinct from AKAP4 homologues of other Felidae or mammals; AKAP4 dysfunction may cause the cheetah’s extremely high (>80 %) pleiomorphic sperm.ConclusionsThe study provides an unprecedented genomic perspective for the rare cheetah, with potential relevance to the species’ natural history, physiological adaptations and unique reproductive disposition.

Genes on B chromosomes of vertebrates

BackgroundThere is a growing body of evidence that B chromosomes, once regarded as totally heterochromatic and genetically inert, harbor multiple segmental duplications containing clusters of ribosomal RNA genes, processed pseudogenes and protein-coding genes. Application of novel molecular approaches further supports complex composition and possible phenotypic effects of B chromosomes.ResultsHere we review recent findings of gene-carrying genomic segments on B chromosomes from different vertebrate groups. We demonstrate that the genetic content of B chromosomes is highly heterogeneous and some B chromosomes contain multiple large duplications derived from various chromosomes of the standard karyotype. Although B chromosomes seem to be mostly homologous to each other within a species, their genetic content differs between species. There are indications that some genomic regions are more likely to be located on B chromosomes.ConclusionsThe discovery of multiple autosomal genes on B chromosomes opens a new discussion about their possible effects ranging from sex determination to fitness and adaptation, their complex interactions with host genome and role in evolution.

Annotated features of domestic cat – Felis catus genome

BackgroundDomestic cats enjoy an extensive veterinary medical surveillance which has described nearly 250 genetic diseases analogous to human disorders. Feline infectious agents offer powerful natural models of deadly human diseases, which include feline immunodeficiency virus, feline sarcoma virus and feline leukemia virus. A rich veterinary literature of feline disease pathogenesis and the demonstration of a highly conserved ancestral mammal genome organization make the cat genome annotation a highly informative resource that facilitates multifaceted research endeavors.FindingsHere we report a preliminary annotation of the whole genome sequence of Cinnamon, a domestic cat living in Columbia (MO, USA), bisulfite sequencing of Boris, a male cat from St. Petersburg (Russia), and light 30× sequencing of Sylvester, a European wildcat progenitor of cat domestication. The annotation includes 21,865 protein-coding genes identified by a comparative approach, 217 loci of endogenous retrovirus-like elements, repetitive elements which comprise about 55.7% of the whole genome, 99,494 new SNVs, 8,355 new indels, 743,326 evolutionary constrained elements, and 3,182 microRNA homologues. The methylation sites study shows that 10.5% of cat genome cytosines are methylated. An assisted assembly of a European wildcat, Felis silvestris silvestris, was performed; variants between F. silvestris and F. catus genomes were derived and compared to F. catus.ConclusionsThe presented genome annotation extends beyond earlier ones by closing gaps of sequence that were unavoidable with previous low-coverage shotgun genome sequencing. The assembly and its annotation offer an important resource for connecting the rich veterinary and natural history of cats to genome discovery.

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

BackgroundB 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.ResultsIn 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.ConclusionsOur data demonstrate independent origins of B chromosomes in the grey brocket and roe deer. We hypothesize that the B chromosomes of these two cervid species represent different stages of B chromosome sequences evolution: probably nascent and similar to autosomal copies in brocket deer, highly derived in roe deer. Based on the presence of the same orthologous protooncogenes in canids and brocket deer Bs we argue that genomic regions involved in B chromosome formation are not random. In addition, our approach is also applicable to the characterization of other evolutionary and clinical rearrangements.

Evolutionary plasticity of acipenseriform genomes.

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

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

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

A High-Resolution SNP Array-Based Linkage Map Anchors a New Domestic Cat Draft Genome Assembly and Provides Detailed Patterns of Recombination

High-resolution genetic and physical maps are invaluable tools for building accurate genome assemblies, and interpreting results of genome-wide association studies (GWAS). Previous genetic and physical maps anchored good quality draft assemblies of the domestic cat genome, enabling the discovery of numerous genes underlying hereditary disease and phenotypes of interest to the biomedical science and breeding communities. However, these maps lacked sufficient marker density to order thousands of shorter scaffolds in earlier assemblies, which instead relied heavily on comparative mapping with related species. A high-resolution map would aid in validating and ordering chromosome scaffolds from existing and new genome assemblies. Here, we describe a high-resolution genetic linkage map of the domestic cat genome based on genotyping 453 domestic cats from several multi-generational pedigrees on the Illumina 63K SNP array. The final maps include 58,055 SNP markers placed relative to 6637 markers with unique positions, distributed across all autosomes and the X chromosome. Our final sex-averaged maps span a total autosomal length of 4464 cM, the longest described linkage map for any mammal, confirming length estimates from a previous microsatellite-based map. The linkage map was used to order and orient the scaffolds from a substantially more contiguous domestic cat genome assembly (Felis catus v8.0), which incorporated ∼20 × coverage of Illumina fragment reads. The new genome assembly shows substantial improvements in contiguity, with a nearly fourfold increase in N50 scaffold size to 18 Mb. We use this map to report probable structural errors in previous maps and assemblies, and to describe features of the recombination landscape, including a massive (∼50 Mb) recombination desert (of virtually zero recombination) on the X chromosome that parallels a similar desert on the porcine X chromosome in both size and physical location.

X chromosome evolution in cetartiodactyla

The phenomenon of a remarkable conservation of the X chromosome in eutherian mammals has been first described by Susumu Ohno in 1964. A notable exception is the cetartiodactyl X chromosome, which varies widely in morphology and G-banding pattern between species. It is hypothesized that this sex chromosome has undergone multiple rearrangements that changed the centromere position and the order of syntenic segments over the last 80 million years of Cetartiodactyla speciation. To investigate its evolution we have selected 26 evolutionarily conserved bacterial artificial chromosome (BAC) clones from the cattle CHORI-240 library evenly distributed along the cattle X chromosome. High-resolution BAC maps of the X chromosome on a representative range of cetartiodactyl species from different branches: pig (Suidae), alpaca (Camelidae), gray whale (Cetacea), hippopotamus (Hippopotamidae), Java mouse-deer (Tragulidae), pronghorn (Antilocapridae), Siberian musk deer (Moschidae), and giraffe (Giraffidae) were obtained by fluorescent in situ hybridization. To trace the X chromosome evolution during fast radiation in specious families, we performed mapping in several cervids (moose, Siberian roe deer, fallow deer, and Pere David’s deer) and bovid (muskox, goat, sheep, sable antelope, and cattle) species. We have identified three major conserved synteny blocks and rearrangements in different cetartiodactyl lineages and found that the recently described phenomenon of the evolutionary new centromere emergence has taken place in the X chromosome evolution of Cetartiodactyla at least five times. We propose the structure of the putative ancestral cetartiodactyl X chromosome by reconstructing the order of syntenic segments and centromere position for key groups.

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

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