Pamela A. Burger

Phylogeography, genetic structure and population divergence time of cheetahs in Africa and Asia: evidence for long-term geographic isolates

The cheetah (Acinonyx jubatus) has been described as a species with low levels of genetic variation. This has been suggested to be the consequence of a demographic bottleneck 10 000–12 000 years ago (ya) and also led to the assumption that only small genetic differences exist between the described subspecies. However, analysing mitochondrial DNA and microsatellites in cheetah samples from most of the historic range of the species we found relatively deep phylogeographic breaks between some of the investigated populations, and most of the methods assessed divergence time estimates predating the postulated bottleneck. Mitochondrial DNA monophyly and overall levels of genetic differentiation support the distinctiveness of Northern‐East African cheetahs (Acinonyx jubatus soemmeringii). Moreover, combining archaeozoological and contemporary samples, we show that Asiatic cheetahs (Acinonyx jubatus venaticus) are unambiguously separated from African subspecies. Divergence time estimates from mitochondrial and nuclear data place the split between Asiatic and Southern African cheetahs (Acinonyx jubatus jubatus) at 32 000–67 000 ya using an average mammalian microsatellite mutation rate and at 4700–44 000 ya employing human microsatellite mutation rates. Cheetahs are vulnerable to extinction globally and critically endangered in their Asiatic range, where the last 70–110 individuals survive only in Iran. We demonstrate that these extant Iranian cheetahs are an autochthonous monophyletic population and the last representatives of the Asiatic subspecies A. j. venaticus. We advocate that conservation strategies should consider the uncovered independent evolutionary histories of Asiatic and African cheetahs, as well as among some African subspecies. This would facilitate the dual conservation priorities of maintaining locally adapted ecotypes and genetic diversity.

Permanent Genetic Resources added to Molecular Ecology Resources Database 1 June 2010 - 31 July 2010.

This article documents the addition of 205 microsatellite marker loci to the Molecular Ecology Resources Database. Loci were developed for the following species: Bagassa guianensis, Bulweria bulwerii, Camelus bactrianus, Chaenogobius annularis, Creontiades dilutus, Diachasmimorpha tryoni, Dioscorea alata, Euhrychiopsis lecontei, Gmelina arborea, Haliotis discus hannai, Hirtella physophora, Melanaphis sacchari, Munida isos, Thaumastocoris peregrinus and Tuberolachnus salignus. These loci were cross‐tested on the following species: Halobaena caerulea, Procellaria aequinoctialis, Oceanodroma monteiroi, Camelus ferus, Creontiades pacificus, Dioscorea rotundata, Dioscorea praehensilis, Dioscorea abyssinica, Dioscorea nummularia, Dioscorea transversa, Dioscorea esculenta, Dioscorea pentaphylla, Dioscorea trifida, Hirtella bicornis, Hirtella glandulosa, Licania alba, Licania canescens, Licania membranaceae, Couepia guianensis and 7 undescribed Thaumastocoris species.

Analysis of the mitochondrial genome of cheetahs (Acinonyx jubatus) with neurodegenerative disease

Abstract The complete mitochondrial genome of Acinonyx jubatus was sequenced and mitochondrial DNA (mtDNA) regions were screened for polymorphisms as candidates for the cause of a neurodegenerative demyelinating disease affecting captive cheetahs. The mtDNA reference sequences were established on the basis of the complete sequences of two diseased and two nondiseased animals as well as partial sequences of 26 further individuals. The A. jubatus mitochondrial genome is 17,047-bp long and shows a high sequence similarity (91%) to the domestic cat. Based on single nucleotide polymorphisms (SNPs) in the control region (CR) and pedigree information, the 18 myelopathic and 12 non-myelopathic cheetahs included in this study were classified into haplotypes I, II and III. In view of the phenotypic comparability of the neurodegenerative disease observed in cheetahs and human mtDNA-associated diseases, specific coding regions including the tRNAs leucine UUR, lysine, serine UCN, and partial complex I and V sequences were screened. We identified a heteroplasmic and a homoplasmic SNP at codon 507 in the subunit 5 (MTND5) of complex I. The heteroplasmic haplotype I-specific valine to methionine substitution represents a nonconservative amino acid change and was found in 11 myelopathic and eight non-myelopathic cheetahs with levels ranging from 29% to 79%. The homoplasmic conservative amino acid substitution valine to alanine was identified in two myelopathic animals of haplotype II. In addition, a synonymous SNP in the codon 76 of the MTND4L gene was found in the single haplotype III animal. The amino acid exchanges in the MTND5 gene were not associated with the occurrence of neurodegenerative disease in captive cheetahs.

Ancient and modern DNA reveal dynamics of domestication and cross-continental dispersal of the dromedary

Significance The dromedary is one of the largest domesticates, sustainably used in arid and hostile environments. It provides food and transport to millions of people in marginal agricultural areas. We show how important long-distance and back-and-forth movements in ancient caravan routes shaped the species’ genetic diversity. Using a global sample set and ancient mitochondrial DNA analyses, we describe the population structure in modern dromedaries and their wild extinct ancestors. Phylogenetic analyses of ancient and modern dromedaries suggest a history of restocking from wild animals from the southeast coast of the Arabian Peninsula. Dromedaries now extend the list of species for which classic models of domestication from a single center and from wild conspecific individuals in isolation are rejected. Dromedaries have been fundamental to the development of human societies in arid landscapes and for long-distance trade across hostile hot terrains for 3,000 y. Today they continue to be an important livestock resource in marginal agro-ecological zones. However, the history of dromedary domestication and the influence of ancient trading networks on their genetic structure have remained elusive. We combined ancient DNA sequences of wild and early-domesticated dromedary samples from arid regions with nuclear microsatellite and mitochondrial genotype information from 1,083 extant animals collected across the species’ range. We observe little phylogeographic signal in the modern population, indicative of extensive gene flow and virtually affecting all regions except East Africa, where dromedary populations have remained relatively isolated. In agreement with archaeological findings, we identify wild dromedaries from the southeast Arabian Peninsula among the founders of the domestic dromedary gene pool. Approximate Bayesian computations further support the “restocking from the wild” hypothesis, with an initial domestication followed by introgression from individuals from wild, now-extinct populations. Compared with other livestock, which show a long history of gene flow with their wild ancestors, we find a high initial diversity relative to the native distribution of the wild ancestor on the Arabian Peninsula and to the brief coexistence of early-domesticated and wild individuals. This study also demonstrates the potential to retrieve ancient DNA sequences from osseous remains excavated in hot and dry desert environments.

The de novo genome assembly and annotation of a female domestic dromedary of North African origin

The single‐humped dromedary (Camelus dromedarius) is the most numerous and widespread of domestic camel species and is a significant source of meat, milk, wool, transportation and sport for millions of people. Dromedaries are particularly well adapted to hot, desert conditions and harbour a variety of biological and physiological characteristics with evolutionary, economic and medical importance. To understand the genetic basis of these traits, an extensive resource of genomic variation is required. In this study, we assembled at 65× coverage, a 2.06 Gb draft genome of a female dromedary whose ancestry can be traced to an isolated population from the Canary Islands. We annotated 21 167 protein‐coding genes and estimated ~33.7% of the genome to be repetitive. A comparison with the recently published draft genome of an Arabian dromedary resulted in 1.91 Gb of aligned sequence with a divergence of 0.095%. An evaluation of our genome with the reference revealed that our assembly contains more error‐free bases (91.2%) and fewer scaffolding errors. We identified ~1.4 million single‐nucleotide polymorphisms with a mean density of 0.71 × 10−3 per base. An analysis of demographic history indicated that changes in effective population size corresponded with recent glacial epochs. Our de novo assembly provides a useful resource of genomic variation for future studies of the camels adaptations to arid environments and economically important traits. Furthermore, these results suggest that draft genome assemblies constructed with only two differently sized sequencing libraries can be comparable to those sequenced using additional library sizes, highlighting that additional resources might be better placed in technologies alternative to short‐read sequencing to physically anchor scaffolds to genome maps.

Genetic diversity and population structure of Mongolian domestic Bactrian camels (Camelus bactrianus)

The tradition of animal husbandry in the context of a nomadic lifestyle has been of great significance in the Mongolian society. Both Bactrian camels and horses have been invaluable for the survival and development of human activities in the harsh arid environment of the Mongolian steppe. As camels offer unique and sustainable opportunities for livestock production in marginal agro-ecological zones, we investigated the current genetic diversity of three local Mongolian camel breeds and compared their levels of variation with common native Mongolian camels distributed throughout the country. Based on mitochondrial and nuclear markers, we found levels of genetic diversity in Mongolian populations similar to that reported for Chinese Bactrian camels and for dromedaries. Little differentiation was detected between single breeds, except for a small group originating from the northwestern Mongolian Altai. We found neither high inbreeding levels in the different breeds nor evidence for a population decline. Although the Mongolian camel census size has severely declined over the past 20 years, our analyses suggest that there still exists a stable population with adequate genetic variation for continued sustainable utilization.

High mitochondrial differentiation levels between wild and domestic Bactrian camels: a basis for rapid detection of maternal hybridization

Hybridization between wild species and their domestic congeners often threatens the gene pool of the wild species. The last wild Bactrian camel (Camelus ferus) populations in Mongolia and China are examples of populations facing such a hybridization threat. To address this key issue in the conservation of wild camels, we analysed wild, hybrid and domestic Bactrian camels (Camelus bactrianus) originating from Mongolia, China and Austria. Through screening of an 804-base-pair mitochondrial fragment, we identified eight mitochondrial haplotypes and found high sequence divergence (1.9%) between C. ferus and C. bactrianus. On the basis of a mitochondrial DNA sequence fixed difference, we developed a diagnostic PCR restriction fragment length polymorphism (PCR-RFLP) assay to differentiate between wild and domestic camel samples. We applied the assay to 81 individuals and confirmed the origin of all samples including five hybrids with known maternal ancestry. The PCR-RFLP system was effective for both traditional (blood, skin) and non-invasive samples (faeces, hair), as well as for museum specimens. Our results demonstrate high levels of mitochondrial differentiation between wild and domestic Bactrian camels and that maternal hybridization can be detected by a rapid and reliable PCR-RFLP system.

Estimating the Population Mutation Rate from a de novo Assembled Bactrian Camel Genome and Cross-Species Comparison with Dromedary ESTs

The Bactrian camel (Camelus bactrianus) and the dromedary (Camelus dromedarius) are among the last species that have been domesticated around 3000–6000 years ago. During domestication, strong artificial (anthropogenic) selection has shaped the livestock, creating a huge amount of phenotypes and breeds. Hence, domestic animals represent a unique resource to understand the genetic basis of phenotypic variation and adaptation. Similar to its late domestication history, the Bactrian camel is also among the last livestock animals to have its genome sequenced and deciphered. As no genomic data have been available until recently, we generated a de novo assembly by shotgun sequencing of a single male Bactrian camel. We obtained 1.6 Gb genomic sequences, which correspond to more than half of the Bactrian camel’s genome. The aim of this study was to identify heterozygous single-nucleotide polymorphisms (SNPs) and to estimate population parameters and nucleotide diversity based on an individual camel. With an average 6.6-fold coverage, we detected over 116 000 heterozygous SNPs and recorded a genome-wide nucleotide diversity similar to that of other domesticated ungulates. More than 20 000 (85%) dromedary expressed sequence tags successfully aligned to our genomic draft. Our results provide a template for future association studies targeting economically relevant traits and to identify changes underlying the process of camel domestication and environmental adaptation.

Combined Hybridization Capture and Shotgun Sequencing for Ancient DNA Analysis of Extinct Wild and Domestic Dromedary Camel

The performance of hybridization capture combined with next‐generation sequencing (NGS) has seen limited investigation with samples from hot and arid regions until now. We applied hybridization capture and shotgun sequencing to recover DNA sequences from bone specimens of ancient‐domestic dromedary (Camelus dromedarius) and its extinct ancestor, the wild dromedary from Jordan, Syria, Turkey and the Arabian Peninsula, respectively. Our results show that hybridization capture increased the percentage of mitochondrial DNA (mtDNA) recovery by an average 187‐fold and in some cases yielded virtually complete mitochondrial (mt) genomes at multifold coverage in a single capture experiment. Furthermore, we tested the effect of hybridization temperature and time by using a touchdown approach on a limited number of samples. We observed no significant difference in the number of unique dromedary mtDNA reads retrieved with the standard capture compared to the touchdown method. In total, we obtained 14 partial mitochondrial genomes from ancient‐domestic dromedaries with 17–95% length coverage and 1.27–47.1‐fold read depths for the covered regions. Using whole‐genome shotgun sequencing, we successfully recovered endogenous dromedary nuclear DNA (nuDNA) from domestic and wild dromedary specimens with 1–1.06‐fold read depths for covered regions. Our results highlight that despite recent methodological advances, obtaining ancient DNA (aDNA) from specimens recovered from hot, arid environments is still problematic. Hybridization protocols require specific optimization, and samples at the limit of DNA preservation need multiple replications of DNA extraction and hybridization capture as has been shown previously for Middle Pleistocene specimens.

The major histocompatibility complex in Old World camelids and low polymorphism of its class II genes

BackgroundThe Major Histocompatibility Complex (MHC) is a genomic region containing genes with crucial roles in immune responses. MHC class I and class II genes encode antigen-presenting molecules expressed on the cell surface. To counteract the high variability of pathogens, the MHC evolved into a region of considerable heterogeneity in its organization, number and extent of polymorphism. Studies of MHCs in different model species contribute to our understanding of mechanisms of immunity, diseases and their evolution. Camels are economically important domestic animals and interesting biomodels. Three species of Old World camels have been recognized: the dromedary (Camelus dromedarius), Bactrian camel (Camelus bactrianus) and the wild camel (Camelus ferus). Despite their importance, little is known about the MHC genomic region, its organization and diversity in camels. The objectives of this study were to identify, map and characterize the MHC region of Old World camelids, with special attention to genetic variation at selected class MHC II loci.ResultsPhysical mapping located the MHC region to the chromosome 20 in Camelus dromedarius. Cytogenetic and comparative analyses of whole genome sequences showed that the order of the three major sub-regions is “Centromere - Class II – Class III – Class I”. DRA, DRB, DQA and DQB exon 2 sequences encoding the antigen binding site of the corresponding class II antigen presenting molecules showed high degree of sequence similarity and extensive allele sharing across the three species. Unexpectedly low extent of polymorphism with low numbers of alleles and haplotypes was observed in all species, despite different geographic origins of the camels analyzed. The DRA locus was found to be polymorphic, with three alleles shared by all three species. DRA and DQA sequences retrieved from ancient DNA samples of Camelus dromedarius suggested that additional polymorphism might exist.ConclusionsThis study provided evidence that camels possess an MHC comparable to other mammalian species in terms of its genomic localization, organization and sequence similarity. We described ancient variation at the DRA locus, monomorphic in most species. The extent of molecular diversity of MHC class II genes seems to be substantially lower in Old World camels than in other mammalian species.