Estimating the population mutation rate from a de novo assembled Bactrian camel genome and cross-species comparison with dromedary ESTs

Analysis of Genetic Diversity and Population Structure of Tarim and Junggar Bactrian Camels Based on Simplified GBS Genome Sequencing

In view of the severe reduction in Bactrian camel germplasm resources, scientific evaluation, protection, and utilization is particularly important. Therefore, it is necessary to investigate the genetic diversity and genetic structure of this species, and identify the genes that have played important roles in its evolution. In this study, 21,971 SNPs were identified in 118 domestic Bactrian camels from the Tarim (n = 60) and Junggar (n = 58) populations using simplified GBS genome sequencing. The results show that Tarim and Junggar Bactrian camels have high nucleotide diversity. A phylogenetic tree constructed using structural analysis, principal component analysis (PCA), and the adjacency method (NJ) showed that Tarim and Junggar Bactrian camels were clustered together. The selection signals revealed that the Tarim and Junggar Bactrian camels shared 108 genes under positive selection, including WNT1, WNT10B, CD14, SEC61A2, DPAGT1, FOXO6, etc. These selected genes were widely involved in the immune system, embryonic development, lipid metabolism, and other processes. From a genomic analysis perspective, the genetic relationship between TLM and ZGE camels is close, with an average Fst of 0.048 and a relatively low average differentiation coefficient between the two populations. In addition, shared selected genes in the long-term depression pathway were significantly enriched in Tarim and Junggar. These findings will offer support and assistance for the exploration of genetic resource preservation, economically significant traits, and the mechanisms underlying biological characteristics, molecular breeding, and disease.

A Study of the Genetic Structure of Hybrid Camels in Kazakhstan

Camel farming is gaining scientific interest due to its unique agricultural characteristics. Camels are versatile for milk and meat production, wool, racing, transport, and tourism. To use their full potential, it is essential to improve our understanding of the genetic structure of these animals. One-humped and two-humped camels have received detailed genetic descriptions, while there is no such information for their hybrids, which outperform their parent species in several agricultural characteristics. Thus, in this study, for the first time, the whole genome sequencing data (WGS) of five hybrid camels bred in the Almaty region of Kazakhstan are presented in comparison with the WGS data of one-humped, two-humped, and wild camels. A total of 43,552,164 single-nucleotide polymorphisms were found across the studied groups. Further comparison of these SNPs showed the following number of private SNPs among the populations: hybrid camels (3,271,083), wild camels (2,515,591), Bactrians (1,244,694), and dromedaries (531,224). The genetic structure of the studied animals was described, and a phylogenetic tree was built to assess their genetic distance. It was found that the studied hybrids are genetically closer to dromedaries since they were on the close branch of the phylogenetic tree.

Aspects of Molecular Genetics in Dromedary Camel

Dromedary camels are unique in their morphological and physiological characteristics and are capable of providing milk and meat even under extreme environmental conditions. Like other species, the dromedary camel has also benefitted from the development of the molecular genetics to increase the knowledge about different aspect in camel genetics (genetic variation, molecular marker, parentage control, gene of interest, whole genome, dating…etc.). In this paper we review the different molecular genetic technics used in this particular species and future prospects. Dromedary genetic studies started in the end of the 1980s with phenotypic evaluation and the attempts to highlight the protein and biochemical diversity. In the 2000s, with the development of molecular markers such as microsatellites, genetic diversity of different types in several countries were estimated and microsatellites were also used for parentage control. In terms of genetic characterization, microsatellites revealed a defined global structure, differentiating East African and South Arabian dromedaries from North African, North Arabian, and South Asian individuals, respectively. Also, mitochondrialDNA sequence analysis of ancient DNA proved to be crucial in resolving domestication processes in dromedaries. Ancient and modern DNA revealed dynamics of domestication and cross-continental dispersion of the dromedary. Nuclear SNPs, single nucleotide polymorphisms changes that occur approximately each 1000 bps in the mammalian genome were also applied in some studies in dromedary. These markers are a very useful alternative to microsatellites and have been employed in some studies on genetic diversity and relevant phenotypic traits in livestock. Finally, thanks to the use of Next Generation Sequencing (NGS) the whole-genome assemblies of the dromedary (Camelus dromedarius) and a work to establish the organization of the dromedary genome at chromosome level were recently published.

An 8.22 Mb Assembly and Annotation of the Alpaca (Vicugna pacos) Y Chromosome

The unique evolutionary dynamics and complex structure make the Y chromosome the most diverse and least understood region in the mammalian genome, despite its undisputable role in sex determination, development, and male fertility. Here we present the first contig-level annotated draft assembly for the alpaca (Vicugna pacos) Y chromosome based on hybrid assembly of short- and long-read sequence data of flow-sorted Y. The latter was also used for cDNA selection providing Y-enriched testis transcriptome for annotation. The final assembly of 8.22 Mb comprised 4.5 Mb of male specific Y (MSY) and 3.7 Mb of the pseudoautosomal region. In MSY, we annotated 15 X-degenerate genes and two novel transcripts, but no transposed sequences. Two MSY genes, HSFY and RBMY, are multicopy. The pseudoautosomal boundary is located between SHROOM2 and HSFY. Comparative analysis shows that the small and cytogenetically distinct alpaca Y shares most of MSY sequences with the larger dromedary and Bactrian camel Y chromosomes. Most of alpaca X-degenerate genes are also shared with other mammalian MSYs, though WWC3Y is Y-specific only in alpaca/camels and the horse. The partial alpaca Y assembly is a starting point for further expansion and will have applications in the study of camelid populations and male biology.

Genomic signatures of domestication in Old World camels

Domestication begins with the selection of animals showing less fear of humans. In most domesticates, selection signals for tameness have been superimposed by intensive breeding for economical or other desirable traits. Old World camels, conversely, have maintained high genetic variation and lack secondary bottlenecks associated with breed development. By re-sequencing multiple genomes from dromedaries, Bactrian camels, and their endangered wild relatives, here we show that positive selection for candidate genes underlying traits collectively referred to as 'domestication syndrome' is consistent with neural crest deficiencies and altered thyroid hormone-based signaling. Comparing our results with other domestic species, we postulate that the core set of domestication genes is considerably smaller than the pan-domestication set - and overlapping genes are likely a result of chance and redundancy. These results, along with the extensive genomic resources provided, are an important contribution to understanding the evolutionary history of camels and the genomic features of their domestication.

Old World camels in a modern world - a balancing act between conservation and genetic improvement

Old World camels have served humans in cross‐continental caravans, transporting people and goods, connecting different cultures and providing milk, meat, wool and draught since their domestication around 3000–6000 years ago. In a world of modern transport and fast connectivity, these beasts of burden seem to be out‐dated. However, a growing demand for sustainable milk and meat production, especially in countries affected by climate change and increasing desertification, brings dromedaries (Camelus dromedarius) and Bactrian camels (Camelus bactrianus) back onstage and into the focus of animal breeders and scientists. In this review on the molecular genetics of these economically important species we give an overview about the evolutionary history, domestication and dispersal of Old World camels, whereas highlighting the need for conservation of wild two‐humped camels (Camelus ferus) as an evolutionarily unique and highly endangered species. We provide cutting‐edge information on the current molecular resources and on‐going sequencing projects. We cannot emphasise enough the importance of balancing the need for improving camel production traits with maintaining the genetic diversity in two domestic species with specific physiological adaptation to a desert environment.

Comparative genomics suggests loss of keratin K24 in three evolutionary lineages of mammals

Keratins are the main cytoskeletal proteins of epithelial cells and changes in the expression of keratins have contributed to the evolutionary adaptation of epithelia to different environments. Keratin K24 was proposed to be a differentiation marker of epidermal keratinocytes but the significance of K24 expression in the epidermis versus other tissues has remained elusive. Here, we show by RT-PCR, western blot, and immunofluorescence analyses that K24 is highly expressed in the epithelium of the cornea whereas its expression levels are significantly lower in other stratified epithelia including in the epidermis. To investigate the evolutionary history of K24, we screened the genome sequences of vertebrates for orthologs of the human KRT24 gene. The results of this comparative genomics study suggested that KRT24 originated in a common ancestor of amniotes and that it was lost independently in three clades of mammals, i.e. camels, cetaceans, and a subclade of pinnipeds comprising eared seals and the walrus. Together, the results of this study identify K24 as component of the cytoskeleton in the human corneal epithelium and reveal previously unknown differences of keratin gene content among mammalian species.

The history of Old World camelids in the light of molecular genetics

Old World camels have come into the focus as sustainable livestock species, unique in their morphological and physiological characteristics and capable of providing vital products even under extreme environmental conditions. The evolutionary history of dromedary and Bactrian camels traces back to the middle Eocene (around 40 million years ago, mya), when the ancestors of Camelus emerged on the North American continent. While the genetic status of the two domestic species has long been established, the wild two-humped camel has only recently been recognized as a separate species, Camelus ferus, based on molecular genetic data. The demographic history established from genome drafts of Old World camels shows the independent development of the three species over the last 100,000 years with severe bottlenecks occurring during the last glacial period and in the recent past. Ongoing studies involve the immune system, relevant production traits, and the global population structure and domestication of Old World camels. Based on the now available whole genome drafts, specific metabolic pathways have been described shedding new light on the camels' ability to adapt to desert environments. These new data will also be at the origin for genome-wide association studies to link economically relevant phenotypes to genotypes and to conserve the diverse genetic resources in Old World camelids.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

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 camel's 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.

Molecular analysis of some camel cytochrome P450 enzymes reveals lower evolution and drug-binding properties

Camels bear unique genotypes and phenotypes for adaptation of their harsh environment. They have unique visual systems, sniffing, water metabolism, and heat-control mechanisms that are different from other creatures. The recent announcement for the complete sequence of camel genome will allow for the discovery of many secrets of camel life. In this context, the genetic bases of camel drug-metabolizing enzymes are still unknown. Furthermore, the genomic content of camel that rendered it highly susceptible to some drugs (as monensin and salinomycin) and became easily intoxicated needs to be investigated. The objectives of this work are the annotation of camel genome and retrieval of camel for cytochrome P450 (CYP) 1A1, 2C, and 3A enzymes. This is followed by comprehensive phylogenetic, evolution, molecular modeling, and docking studies. In comparison with the human enzymes, camel CYPs showed lower evolution rate, especially CYP1A1. Furthermore, the binding of monensin, salinomycin, alfa-naphthoflavone, felodepine, and ritonavir was weaker in camel enzymes. Interestingly, rerank score indicated instable binding of monensin and salinomycin with camel CYP1A1 as well as salinomycin with camel CYP2C. The results of this work suggest that camels are more susceptible to toxicity with compounds undergoing metabolic oxidation. This conclusion was based on lower evolution rate and lower binding potency of camels compared with the human enzymes.

Camelid genomes reveal evolution and adaptation to desert environments

Bactrian camel (Camelus bactrianus), dromedary (Camelus dromedarius) and alpaca (Vicugna pacos) are economically important livestock. Although the Bactrian camel and dromedary are large, typically arid-desert-adapted mammals, alpacas are adapted to plateaus. Here we present high-quality genome sequences of these three species. Our analysis reveals the demographic history of these species since the Tortonian Stage of the Miocene and uncovers a striking correlation between large fluctuations in population size and geological time boundaries. Comparative genomic analysis reveals complex features related to desert adaptations, including fat and water metabolism, stress responses to heat, aridity, intense ultraviolet radiation and choking dust. Transcriptomic analysis of Bactrian camels further reveals unique osmoregulation, osmoprotection and compensatory mechanisms for water reservation underpinned by high blood glucose levels. We hypothesize that these physiological mechanisms represent kidney evolutionary adaptations to the desert environment. This study advances our understanding of camelid evolution and the adaptation of camels to arid-desert environments.

Assembly information services in the European Nucleotide Archive

The European Nucleotide Archive (ENA; http://www.ebi.ac.uk/ena) is a repository for the world public domain nucleotide sequence data output. ENA content covers a spectrum of data types including raw reads, assembly data and functional annotation. ENA has faced a dramatic growth in genome assembly submission rates, data volumes and complexity of datasets. This has prompted a broad reworking of assembly submission services, for which we now reach the end of a major programme of work and many enhancements have already been made available over the year to components of the submission service. In this article, we briefly review ENA content and growth over 2013, describe our rapidly developing services for genome assembly information and outline further major developments over the last year.