Chromosome-level assembly of wild Bactrian camel genome reveals organization of immune gene loci

Transcriptomics-Based Study of Immune Genes Associated with Subclinical Mastitis in Bactrian Camels

The significant increase in demand for camel milk has led to a rapid increase in the number of Bactrian camels. However, the widespread occurrence of mastitis significantly impacts the development of the Bactrian camel milk industry and poses a public health risk. Despite this, there is a lack of research on the transcriptional response, immune response pathways, and changes in core genes of Bactrian camels with subclinical mastitis. This study aimed to reveal the changes in immune-related response pathways and gene transcription levels in Bactrian camels with subclinical mastitis by analyzing the blood transcriptional response after the occurrence of subclinical mastitis in natural conditions. This study focused on 7-year-old Bactrian camels and collected 2 mL of blood from the camels that tested positive with a 4-peak California Mastitis Test (CMT) and those that tested negative with a 3-peak CMT. RNA sequencing (RNA-Seq) technology was used to analyze gene expression in the blood samples. Gene expression was verified using quantitative reverse transcription polymerase chain reaction (RT-qPCR). Overall, 1722 differentially expressed genes were sequenced in the blood samples of CMT-positive and CMT-negative Bactrian camels, including 1061 upregulated and 661 downregulated genes. After conducting gene ontology functional enrichment, 453 differentially expressed genes were identified. We also discovered pathways such as immune response, the G-protein-coupled receptor signaling pathway, and internal signal transmission. Kyoto Encyclopedia of Genes and Genomes (KEGG) signaling pathway enrichment detected 668 differentially expressed genes annotated onto 309 metabolic pathways, with significantly enriched immune pathways including cytokine-cytokine receptor interaction, complex and coalescence cascades, natural killer cell-mediated cytotoxicity, and T helper type 17 cell differentiation, among others. Through a STRING protein interaction database and cytoscape analysis, it was found that core differentially expressed genes related to immunity included IL10, CCL5, IL1B, OSM, TNFRSF1B, IL7, and CCR3, among others. The RT-qPCR results for six randomly selected core differentially expressed genes showed that the RT-qPCR expression pattern was consistent with the RNA Seq results. The immune-related genes in Bactrian camels affected by subclinical mastitis are primarily concentrated in the immune response and the cytokine-cytokine receptor interaction pathway. Given the importance of these pathways and the connections among related genes, the immune genes within these pathways may play a crucial role in the pathogenesis of subclinical mastitis in Bactrian camels. This study provides a valuable reference for investigating the immune regulatory mechanisms of subclinical mastitis in Bactrian camels.

Single-cell 5' RNA sequencing of camelid peripheral B cells provides insights into cellular basis of heavy-chain antibody production

Camelids produce both conventional tetrameric antibodies (Abs) and dimeric heavy-chain antibodies (HCAbs). Although B cells that generate these two types of Abs exhibit distinct B cell receptors (BCRs), whether these two B cell populations differ in their phenotypes and developmental processes remains unclear. Here, we performed single-cell 5' RNA profiling of peripheral blood mononuclear cell samples from Bactrian camels before and after immunization. We characterized the functional subtypes and differentiation trajectories of circulating B cells in camels, and reconstructed single-cell BCR sequences. We found that in contrast to humans, the proportion of T-bet+ B cells was high among camelid peripheral B cells. Several marker genes of human B cell subtypes, including CD27 and IGHD, were expressed at low levels in the corresponding camel B cell subtypes. Camelid B cells expressing variable genes of HACbs (VHH) were widely present in various functional subtypes and showed highly overlapping differentiation trajectories with B cells expressing variable genes of conventional Abs (VH). After immunization, the transcriptional changes in VHH+ and VH+ B cells were largely consistent. Through structure modeling, we identified a variety of scaffold types among the reconstructed VHH sequences. Our study provides insights into the cellular context of HCAb production in camels and lays the foundation for developing single-B cell-based camelid single-domain Ab screening.

Single domain antibody: Development and application in biotechnology and biopharma

Heavy-chain antibodies (HCAbs) are a unique type of antibodies devoid of light chains, and comprised of two heavy chains-only that recognize their cognate antigen by virtue of a single variable domain also referred to as VHH, single domain antibody (sdAb), or nanobody (Nb). These functional HCAbs, serendipitous discovered about three decades ago, are exclusively found in camelids, comprising dromedaries, camels, llamas, and vicugnas. Nanobodies have become an essential tool in biomedical research and medicine, both in diagnostics and therapeutics due to their beneficial properties: small size, high stability, strong antigen-binding affinity, low immunogenicity, low production cost, and straightforward engineering into more potent affinity reagents. The occurrence of HCAbs in camelids remains intriguing. It is believed to be an evolutionary adaptation, equipping camelids with a robust adaptive immune defense suitable to respond to the pressure from a pathogenic invasion necessitating a more profound antigen recognition and neutralization. This evolutionary innovation led to a simplified HCAb structure, possibly supported by genetic mutations and drift, allowing adaptive mutation and diversification in the heavy chain variable gene and constant gene regions. Beyond understanding their origins, the application of nanobodies has significantly advanced over the past 30 years. Alongside expanding laboratory research, there has been a rapid increase in patent application for nanobodies. The introduction of commercial nanobody drugs such as Cablivi, Nanozora, Envafolimab, and Carvykti has boosted confidence among in their potential. This review explores the evolutionary history of HCAbs, their ontogeny, and applications in biotechnology and pharmaceuticals, focusing on approved and ongoing medical research pipelines.

Mitochondrial DNA of the Arabian Camel Camelus dromedarius

The Camelidae family, ranging from southwest Asia to north Africa, South America, and Australia, includes key domesticated species adapted to diverse environments. Among these, the Arabian camel (Camelus dromedarius) is vital to the cultural and economic landscape of the Arabian Peninsula. This review explores the mitochondrial DNA of the dromedary camel, focusing on the D-loop region to understand its genetic diversity, maternal inheritance, and evolutionary history. We aim to investigate the unique characteristics of Arabian camel mtDNA, analyze the D-loop for genetic diversity and maternal lineage patterns, and explore the implications of mitochondrial genomic studies for camel domestication and adaptation. Key findings on mtDNA structure and variation highlight significant genetic differences and adaptive traits. The D-loop, essential for mtDNA replication and transcription, reveals extensive polymorphisms and haplotypes, providing insights into dromedary camel domestication and breeding history. Comparative analyses with other camelid species reveal unique genetic signatures in the Arabian camel, reflecting its evolutionary and adaptive pathways. Finally, this review integrates recent advancements in mitochondrial genomics, demonstrating camel genetic diversity and potential applications in conservation and breeding programs. Through comprehensive mitochondrial genome analysis, we aim to enhance the understanding of Camelidae genetics and contribute to the preservation and improvement of these vital animals.

The immunoglobulin A isotype of the Arabian camel (Camelus dromedarius) preserves the dualistic structure of unconventional single-domain and canonical heavy chains

Introduction

The evolution of adaptive immunity in Camelidae resulted in the concurrent expression of classic heterotetrameric and unconventional homodimeric heavy chain-only IgG antibodies. Heavy chain-only IgG bears a single variable domain and lacks the constant heavy (CH) γ1 domain required for pairing with the light chain. It has not been reported whether this distinctive feature of IgG is also observed in the IgA isotype.

Methods

Gene-specific primers were used to generate an IgA heavy chain cDNA library derived from RNA extracted from the dromedary's third eyelid where isolated lymphoid follicles and plasma cells abound at inductive and effector sites, respectively.

Results

Majority of the cDNA clones revealed hallmarks of heavy chain-only antibodies, i.e. camelid-specific amino acid substitutions in framework region 1 and 2, broad length distribution of complementarity determining region 3, and the absence of the CHα1 domain. In a few clones, however, the cDNA of the canonical IgA heavy chain was amplified which included the CHα1 domain, analogous to CHγ1 domain in IgG1 subclass. Moreover, we noticed a short, proline-rich hinge, and, at the N-terminal end of the CHα3 domain, a unique, camelid-specific pentapeptide of undetermined function, designated as the inter-α region. Immunoblots using rabbit anti-camel IgA antibodies raised against CHα2 and CHα3 domains as well as the inter-α region revealed the expression of a ~52 kDa and a ~60 kDa IgA species, corresponding to unconventional and canonical IgA heavy chain, respectively, in the third eyelid, trachea, small and large intestine. In contrast, the leporine anti-CHα1 antibody detected canonical, but not unconventional IgA heavy chain, in all the examined tissues, milk, and serum, in addition to another hitherto unexplored species of ~45 kDa in milk and serum. Immunohistology using anti-CHα domain antibodies confirmed the expression of both variants of IgA heavy chains in plasma cells in the third eyelid's lacrimal gland, conjunctiva, tracheal and intestinal mucosa.

Conclusion

We found that in the dromedary, the IgA isotype has expanded the immunoglobulin repertoire by co-expressing unconventional and canonical IgA heavy chains, comparable to the IgG class, thus underscoring the crucial role of heavy chain-only antibodies not only in circulation but also at the mucosal frontiers.

OmicsSuite: a customized and pipelined suite for analysis and visualization of multi-omics big data

With the advancements in high-throughput sequencing technologies such as Illumina, PacBio, and 10X Genomics platforms, and gas/liquid chromatography-mass spectrometry, large volumes of biological data in multiple formats can now be obtained through multi-omics analysis. Bioinformatics is constantly evolving and seeking breakthroughs to solve multi-omics problems; however, it is challenging for most experimental biologists to analyse data using command-line interfaces, coding, and scripting. Based on experience with multi-omics, we have developed OmicsSuite, a desktop suite that comprehensively integrates statistics and multi-omics analysis and visualization. The suite has 175 sub-applications in 12 categories, including Sequence, Statistics, Algorithm, Genomics, Transcriptomics, Enrichment, Proteomics, Metabolomics, Clinical, Microorganism, Single Cell, and Table Operation. We created the user interface with Sequence View, Table View, and intelligent components based on JavaFX and the popular Shiny framework. The multi-omics analysis functions were developed based on BioJava and 300+ packages provided by the R CRAN and Bioconductor communities, and it encompasses over 3000 adjustable parameter interfaces. OmicsSuite can directly read multi-omics raw data in FastA, FastQ, Mutation Annotation Format, mzML, Matrix, and HDF5 formats, and the programs emphasize data transfer directions and pipeline analysis functions. OmicsSuite can produce pre-publication images and tables, allowing users to focus on biological aspects. OmicsSuite offers multi-omics step-by-step workflows that can be easily applied to horticultural plant breeding and molecular mechanism studies in plants. It enables researchers to freely explore the molecular information contained in multi-omics big data (Source: https://github.com/OmicsSuite/, Website: https://omicssuite.github.io, v1.3.9).

A chromosome-level genome assembly of the beet armyworm Spodoptera exigua

BACKGROUND

The beet armyworm Spodoptera exigua is a polyphagous caterpillar that causes serious damage to many species of crops and vegetables. To gain insight into how this polyphagous insect differs from less harmful oligophagous species, we generated a chromosome-level assembly and compared it to closely related species with the same or different feeding habits.

RESULTS

Based on Illumina and Pacific Biosciences data and Hi-C technology, 425.6 Mb of genome sequences were anchored and oriented into 31 linkage groups, with an N50 length of 14.8 Mb. A total of 24,649 gene models were predicted, of which 97.4% were identified in the genome assembly. Chemosensory genes are vital for locating food: of the four main families, odorant-binding proteins, chemosensory proteins and olfactory receptors showed little difference, whereas gustatory receptors are greatly expanded in S. exigua. Examination of other polyphagous insects confirmed this difference from oligophagous congeners and further identified the bitter receptor subfamily as being particularly affected.

CONCLUSION

Our high-quality genome sequence for beet armyworm identified a key expansion of the bitter gustatory receptor subfamily in this and other pests that differs crucially from more benign relatives and offers insight into the biology and possible future means of control for these economically important insects.

Characterization of heavy-chain antibody gene repertoires in Bactrian camels

Camelids are the only mammals that can produce functional heavy-chain antibodies (HCAbs). Although HCAbs were discovered over 30 years ago, the antibody gene repertoire of Bactrian camels remains largely underexplored. To characterize the diversity of variable genes of HCAbs (VHHs), germline and rearranged VHH repertoires are constructed. Phylogenetics analysis shows that all camelid VHH genes are derived from a common ancestor and the nucleotide diversity of VHHs is similar across all camelid species. While species-specific hallmark sites are identified, the non-canonical cysteines specific to VHHs are distinct in Bactrian camels and dromedaries compared with alpacas. Though low divergence at the germline repertoire between wild and domestic Bactrian camels, higher expression of VHHs is observed in some wild Bactrian camels than that of domestic ones. This study not only adds our understanding of VHH repertoire diversity across camelids, but also provides useful resources for HCAb engineering.

Best genome sequencing strategies for annotation of complex immune gene families in wildlife

BACKGROUND

The biodiversity crisis and increasing impact of wildlife disease on animal and human health provides impetus for studying immune genes in wildlife. Despite the recent boom in genomes for wildlife species, immune genes are poorly annotated in nonmodel species owing to their high level of polymorphism and complex genomic organisation. Our research over the past decade and a half on Tasmanian devils and koalas highlights the importance of genomics and accurate immune annotations to investigate disease in wildlife. Given this, we have increasingly been asked the minimum levels of genome quality required to effectively annotate immune genes in order to study immunogenetic diversity. Here we set out to answer this question by manually annotating immune genes in 5 marsupial genomes and 1 monotreme genome to determine the impact of sequencing data type, assembly quality, and automated annotation on accurate immune annotation.

RESULTS

Genome quality is directly linked to our ability to annotate complex immune gene families, with long reads and scaffolding technologies required to reassemble immune gene clusters and elucidate evolution, organisation, and true gene content of the immune repertoire. Draft-quality genomes generated from short reads with HiC or 10× Chromium linked reads were unable to achieve this. Despite mammalian BUSCOv5 scores of up to 94.1% amongst the 6 genomes, automated annotation pipelines incorrectly annotated up to 59% of manually annotated immune genes regardless of assembly quality or method of automated annotation.

CONCLUSIONS

Our results demonstrate that long reads and scaffolding technologies, alongside manual annotation, are required to accurately study the immune gene repertoire of wildlife species.

Microsatellite markers of the major histocompatibility complex genomic region of domestic camels

We identified and characterized 11 polymorphic microsatellite markers suitable for routine testing (three in the MHC class I sub-region, four in MHC class II and four in the MHC class III sub-region) of dromedaries and Bactrian camels. In total, 38 dromedaries and 33 Bactrian camels were genotyped, and interspecific differences were observed in the numbers of alleles and in allelic frequencies, as well as in the observed heterozygosity. These loci may be used as markers to study the adaptive genetic diversity of the MHC region in Old World camels.

Coevolved Canonical Loops Conformations of Single-Domain Antibodies: A Tale of Three Pockets Playing Musical Chairs

Single-domain antibodies (sdAbs) are a promising class of biotherapeutics with unique structural traits within their paratope region. The distribution of canonical conformations explored by their complementarity determining region (CDR) loops differs to some extent from conventional two-chain Fv fragments of monoclonal antibodies (mAbs). In this study, we explored in detail the canonical structures of sdAb CDR-H1 and CDR-H2 loops and compared those with mAbs from the IGHV3 and IGHV1 gene families. We surveyed the antibody structures catalogued in SAbDab and clustered the CDR canonical loops in Cartesian space. While most of the sdAb clusters were sub-populations of previously defined canonical Fv conformations of CDR-H1 and CDR-H2, our stricter clustering approach defined narrower clusters in sequence-space. Meticulous visual inspection of sub-populations allowed a clearer understanding of sequence-structure relationships. The packing densities within structural pockets contacted by CDR-H1 and CDR-H2 canonical conformations were analyzed on the premise that these pockets cannot be left vacant as they would leave exposed supportive hydrophobic residues. The fine resolution of the canonical clusters defined here revealed unique signatures within these pockets, including distinct structural complementarities between CDR-H1 and CDR-H2 canonical clusters, which could not be perceived with the previous coarser clusters. We highlight examples where a single residue change in CDR-H1 sequence is sufficient to induce a dramatic population shift in CDR-H2 conformation. This suggests that preferences in combining CDR-H1 and CDR-H2 emerged naturally during antibody evolution, leading to preferred sets of conserved amino acids at key positions in the framework as well as within the CDR loops. We outline a game of musical chairs that is necessary to maintain the integrity of the antibody structures that arose during evolution. Our study also provides refined CDR-H1 and CDR-H2 structural templates for sdAb homology modeling that could be leveraged for improved antibody design.

Nanobodies: From Serendipitous Discovery of Heavy Chain-Only Antibodies in Camelids to a Wide Range of Useful Applications

The presence of unique heavy chain-only antibodies (HCAbs) in camelids was discovered at Vrije Universiteit Brussel (VUB, Brussels, Belgium) at a time when many researchers were exploring the cloning and expression of smaller antigen-binding fragments (Fv and Fab) from hybridoma-derived antibodies. The potential importance of this discovery was anticipated, and efforts were immediately undertaken to understand the emergence and ontogeny of these HCAbs as well as to investigate the applications of the single-domain antigen-binding variable domains of HCAbs (nanobodies). Nanobodies were demonstrated to possess multiple biochemical and biophysical advantages over other antigen-binding antibody fragments and alternative scaffolds. Today, nanobodies have a significant and growing impact on research, biotechnology, and medicine.

Probiotics: their action against pathogens can be turned around

Probiotics when applied in complex evolving (micro-)ecosystems, might be selectively beneficial or detrimental to pathogens when their prophylactic efficacies are prone to ambient interactions. Here, we document a counter-intuitive phenomenon that probiotic-treated zebrafish (Danio rerio) were respectively healthy at higher but succumbed at lower level of challenge with a pathogenic Vibrio isolate. This was confirmed by prominent dissimilarities in fish survival and histology. Based upon the profiling of the zebrafish microbiome, and the probiotic and the pathogen shared gene orthogroups (genetic niche overlaps in genomes), this consequently might have modified the probiotic metabolome as well as the virulence of the pathogen. Although it did not reshuffle the architecture of the commensal microbiome of the vertebrate host, it might have altered the probiotic-pathogen inter-genus and intra-species communications. Such in-depth analyses are needed to avoid counteractive phenomena of probiotics and to optimise their efficacies to magnify human and animal well-being. Moreover, such studies will be valuable to improve the relevant guidelines published by organisations such as FAO, OIE and WHO.

Innate and Adaptive Immune Genes Associated with MERS-CoV Infection in Dromedaries

The recent SARS-CoV-2 pandemic has refocused attention to the betacoronaviruses, only eight years after the emergence of another zoonotic betacoronavirus, the Middle East respiratory syndrome coronavirus (MERS-CoV). While the wild source of SARS-CoV-2 may be disputed, for MERS-CoV, dromedaries are considered as source of zoonotic human infections. Testing 100 immune-response genes in 121 dromedaries from United Arab Emirates (UAE) for potential association with present MERS-CoV infection, we identified candidate genes with important functions in the adaptive, MHC-class I (HLA-A-24-like) and II (HLA-DPB1-like), and innate immune response (PTPN4, MAGOHB), and in cilia coating the respiratory tract (DNAH7). Some of these genes previously have been associated with viral replication in SARS-CoV-1/-2 in humans, others have an important role in the movement of bronchial cilia. These results suggest similar host genetic pathways associated with these betacoronaviruses, although further work is required to better understand the MERS-CoV disease dynamics in both dromedaries and humans.

The Genomic Organisation of the TRA/TRD Locus Validates the Peculiar Characteristics of Dromedary δ-Chain Expression

The role of γδ T cells in vertebrate immunity is still an unsolved puzzle. Species such as humans and mice display a low percentage of these T lymphocytes (i.e., “γδ low species”) with a restricted diversity of γδ T cell receptors (TR). Conversely, artiodactyl species (i.e., “γδ high species”) account for a high proportion of γδ T cells with large γ and δ chain repertoires. The genomic organisation of the TR γ (TRG) and δ (TRD) loci has been determined in sheep and cattle, noting that a wide number of germline genes that encode for γ and δ chains characterise their genomes. Taking advantage of the current improved version of the genome assembly, we have investigated the genomic structure and gene content of the dromedary TRD locus, which, as in the other mammalian species, is nested within the TR α (TRA) genes. The most remarkable finding was the identification of a very limited number of variable germline genes (TRDV) compared to sheep and cattle, which supports our previous expression analyses for which the somatic hypermutation mechanism is able to enlarge and diversify the primary repertoire of dromedary δ chains. Furthermore, the comparison between genomic and expressed sequences reveals that D genes, up to four incorporated in a transcript, greatly contribute to the increased diversity of the dromedary δ chain antigen binding-site.

A Deadly Cargo: Gene Repertoire of Cytotoxic Effector Proteins in the Camelidae

Cytotoxic T cells and natural killer cells can kill target cells based on their expression and release of perforin, granulysin, and granzymes. Genes encoding these molecules have been only poorly annotated in camelids. Based on bioinformatic analyses of genomic resources, sequences corresponding to perforin, granulysin, and granzymes were identified in genomes of camelids and related ungulate species, and annotation of the corresponding genes was performed. A phylogenetic tree was constructed to study evolutionary relationships between the species analyzed. Re-sequencing of all genes in a panel of 10 dromedaries and 10 domestic Bactrian camels allowed analyzing their individual genetic polymorphisms. The data showed that all extant Old World camelids possess functional genes for two pore-forming proteins (PRF1, GNLY) and six granzymes (GZMA, GZMB, GZMH, GZMK, GZMM, and GZMO). All these genes were represented as single copies in the genome except the GZMH gene exhibiting interspecific differences in the number of loci. High protein sequence similarities with other camelid and ungulate species were observed for GZMK and GZMM. The protein variability in dromedaries and Bactrian camels was rather low, except for GNLY and chymotrypsin-like granzymes (GZMB, GZMH).

Recent Advances in Camel Immunology

Camels are domesticated animals that are highly adapted to the extreme desert ecosystem with relatively higher resistance to a wide range of pathogens compared to many other species from the same geographical region. Recently, there has been increased interest in the field of camel immunology. As the progress in the analysis of camel immunoglobulins has previously been covered in many recent reviews, this review intends to summarize published findings related to camel cellular immunology with a focus on the phenotype and functionality of camel leukocyte subpopulations. The review also describes the impact of different physiological (age and pregnancy) and pathological (e.g. infection) conditions on camel immune cells. Despite the progress achieved in the field of camel immunology, there are gaps in our complete understanding of the camel immune system. Questions remain regarding innate recognition mechanisms, the functional characterization of antigen-presenting cells, and the characterization of camel NK and cytotoxic T cells.

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.

Applications of Nanobodies

Unique, functional, homodimeric heavy chain-only antibodies, devoid of light chains, are circulating in the blood of Camelidae. These antibodies recognize their cognate antigen via one single domain, known as VHH or Nanobody. This serendipitous discovery made three decades ago has stimulated a growing number of researchers to generate highly specific Nanobodies against a myriad of targets. The small size, strict monomeric state, robustness, and easy tailoring of these Nanobodies have inspired many groups to design innovative Nanobody-based multi-domain constructs to explore novel applications. As such, Nanobodies have been employed as an exquisite research tool in structural, cell, and developmental biology. Furthermore, Nanobodies have demonstrated their benefit for more sensitive diagnostic tests. Finally, several Nanobody-based constructs have been designed to develop new therapeutic products.

Nucleotide diversity of functionally different groups of immune response genes in Old World camels based on newly annotated and reference-guided assemblies

Background

Immune-response (IR) genes have an important role in the defense against highly variable pathogens, and therefore, diversity in these genomic regions is essential for species’ survival and adaptation. Although current genome assemblies from Old World camelids are very useful for investigating genome-wide diversity, demography and population structure, they have inconsistencies and gaps that limit analyses at local genomic scales. Improved and more accurate genome assemblies and annotations are needed to study complex genomic regions like adaptive and innate IR genes.

Results

In this work, we improved the genome assemblies of the three Old World camel species – domestic dromedary and Bactrian camel, and the two-humped wild camel – via different computational methods. The newly annotated dromedary genome assembly CamDro3 served as reference to scaffold the NCBI RefSeq genomes of domestic Bactrian and wild camels. These upgraded assemblies were then used to assess nucleotide diversity of IR genes within and between species, and to compare the diversity found in immune genes and the rest of the genes in the genome. We detected differences in the nucleotide diversity among the three Old World camelid species and between IR gene groups, i.e., innate versus adaptive. Among the three species, domestic Bactrian camels showed the highest mean nucleotide diversity. Among the functionally different IR gene groups, the highest mean nucleotide diversity was observed in the major histocompatibility complex.

Conclusions

The new camel genome assemblies were greatly improved in terms of contiguity and increased size with fewer scaffolds, which is of general value for the scientific community. This allowed us to perform in-depth studies on genetic diversity in immunity-related regions of the genome. Our results suggest that differences of diversity across classes of genes appear compatible with a combined role of population history and differential exposures to pathogens, and consequent different selective pressures.

Mitochondrial DNA variation and phylogeography of Old World camels

Objective

Old World camels are a valuable genetic resource for many countries around the world due to their adaptation to the desert environment. At present, Old World camels have encountered the challenge of unprecedented loss of genetic resources. Through our research, we would reveal the population structure and genetic variation in Old World camel populations, which provides a theoretical basis for understanding the germplasm resources and origin and evolution of different Old World camel populations.

Methods

In the present study, we assessed mtDNA control region sequences of 182 individuals from Old World camels to unravel genetic diversity, phylogeography, and demographic dynamics.

Results

Thirty-two haplotypes confirmed by 54 polymorphic sites were identified in the 156 sequences, which included 129 domestic and 27 wild Bactrian camels. Meanwhile, 14 haplotypes were defined by 47 polymorphic sites from 26 sequences in the dromedaries. The wild Bactrian camel population showed the lowest haplotype and nucleotide diversity, while the dromedaries investigated had the highest. The phylogenetic analysis suggests that there are several shared haplotypes in different Bactrian camel populations, and that there has been genetic introgression between domestic Bactrian camels and dromedaries. In addition, positive values of Tajima’s D and Fu’s Fs test demonstrated a decrease in population size and/or balancing selection in the wild Bactrian camel population. In contrast, the negative values of Tajima’s D and Fu’s Fs test in East Asian Bactrian camel populations explained the demographic expansion and/or positive selection.

Conclusion

In summary, we report novel information regarding the genetic diversity, population structure and demographic dynamics of Old World camels. The findings obtained from the present study reveal that abundant genetic diversity occurs in domestic Bactrian camel populations and dromedaries, while there are low levels of haplotype and nucleotide diversity in the wild Bactrian camel population.