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

Geographical distribution, genetic diversity, and environmental adaptations of dromedary camel breeds in Saudi Arabia

The dromedary camel (Camelus dromedarius) in Saudi Arabia exhibits significant genetic diversity, driven by adaptation to diverse ecological niches such as deserts, mountains, and coastal areas. This study explores the genetic structure of these camel populations, correlating their genetic diversity with geographical regions rather than ecological classifications. Through whole-genome sequencing of 63 camel genomes, we identified substantial differences in heterozygosity and inbreeding across different ecotypes, particularly noting higher genetic diversity in mountainous populations and lower diversity in coastal populations. The study also revealed significant enrichment of specific gene sets associated with environmental adaptation, such as the HECT domain in desert populations, which is crucial for maintaining protein integrity under extreme conditions. Principal component and admixture analyses further highlighted the genetic distinctiveness of certain breeds, particularly the Awarik (beach ecotype), which showed signs of genetic isolation.

Protective effect of epidermal growth factor on cryopreservation of dromedary camel epididymal spermatozoa: Evidence from in vitro and in silico studies

Epidermal growth factor (EGF) plays a crucial role in maintaining male reproductive capacity in mammals, however, its protective effects on cryopreserved dromedary camel epididymal spermatozoa have not been thoroughly investigated. This study aims to investigate the potential protective role of EGF on cryopreserved camel epididymal spermatozoa, supported by evidence from a molecular docking study. We assessed sperm motility, kinematics parameters, oxidative stress, ultrastructural changes, apoptosis, and molecular docking markers in camel epididymal spermatozoa following cryopreservation. Camel epididymal spermatozoa (n = 30 pairs of testes) were collected from local slaughterhouses. The epididymal spermatozoa were diluted with a freezing medium (SHOTOR extender) supplemented with different concentrations of EGF; 0 (EGF0), 50 (EGF50), 100 (EGF100), 200 (EGF200), and 400 (EGF400) ng/mL in SHOTOR extender and cryopreserved using a standard protocol. All EGF groups showed significant improvements in sperm progressive motility, viability, and sperm membrane function after equilibration at 5 °C for 24 hours. Regarding frozen-thawed samples, sperm progressive motility and some kinematic parameters (DAP, VSL, VCL and AHL) were significantly higher in the EFG400 group compared to the other groups (P < 0.01). A significant increase in the percentage of live/acrosome-intact sperm was observed, accompanied by a significant decrease in malondialdehyde levels in all EGF groups (P < 0.05). Both the EGF200 and EGF400 groups showed significantly higher sperm viability and significantly lower percentages of apoptotic and necrotic sperm compared to the other groups (P < 0.05). EGF supplementation preserved the ultrastructural integrity and cryotolerance of epididymal camel spermatozoa. The docking analysis indicated that EGF exhibited higher binding affinity with apoptosis sperm markers, including caspase-3 and bcl-2-associated X (Bax) proteins, with binding energies of -502.0 and -621.0 kcal/mol, respectively. In conclusion, the addition of EGF to SHOTOR extender was found to have beneficial effects on sperm motility, kinematics parameters, sperm viability, acrosome integrity, sperm ultrastructural features, and reduced oxidative stress and apoptosis-like changes in cryopreserved epididymal camel spermatozoa.

The Evolution and Role of Molecular Tools in Measuring Diversity and Genomic Selection in Livestock Populations (Traditional and Up-to-Date Insights): A Comprehensive Exploration

Distinctive molecular approaches and tools, particularly high-throughput SNP genotyping, have been applied to determine and discover SNPs, potential genes of interest, indicators of evolutionary selection, genetic abnormalities, molecular indicators, and loci associated with quantitative traits (QTLs) in various livestock species. These methods have also been used to obtain whole-genome sequencing (WGS) data, enabling the implementation of genomic selection. Genomic selection allows for selection decisions based on genomic-estimated breeding values (GEBV). The estimation of GEBV relies on the calculation of SNP effects using prediction equations derived from a subset of individuals in the reference population who possess both SNP genotypes and phenotypes for target traits. Compared to traditional methods, modern genomic selection methods offer advantages for sex-limited traits, low heritability traits, late-measured traits, and the potential to increase genetic gain by reducing generation intervals. The current availability of high-density genotyping and next-generation sequencing data allow for genome-wide scans for selection. This investigation provides an overview of the essential role of advanced molecular tools in studying genetic diversity and implementing genomic selection. It also highlights the significance of adaptive selection in light of new high-throughput genomic technologies and the establishment of selective comparisons between different genomes. Moreover, this investigation presents candidate genes and QTLs associated with various traits in different livestock species, such as body conformation, meat production and quality, carcass characteristics and composition, milk yield and composition, fertility, fiber production and characteristics, and disease resistance.

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.

Whole-Genome Resequencing Analysis of the Camelus bactrianus (Bactrian Camel) Genome Identifies Mutations and Genes Affecting Milk Production Traits

Milk production is an important trait that influences the economic value of camels. However, the genetic regulatory mechanisms underlying milk production in camels have not yet been elucidated. We aimed to identify candidate molecular markers that affect camel milk production. We classified Junggar Bactrian camels (9-10-year-old) as low-yield (<1.96 kg/d) or high-yield (>2.75 kg/d) based on milk production performance. Milk fat (5.16 ± 0.51 g/100 g) and milk protein (3.59 ± 0.22 g/100 g) concentrations were significantly lower in high-yielding camels than those in low-yielding camels (6.21 ± 0.59 g/100 g, and 3.93 ± 0.27 g/100 g, respectively) (p < 0.01). There were no apparent differences in gland tissue morphology between the low- and high-production groups. Whole-genome resequencing of 12 low- and 12 high-yield camels was performed. The results of selection mapping methods, performed using two methods (FST and θπ), showed that 264 single nucleotide polymorphism sites (SNPs) overlapped between the two methods, identifying 181 genes. These genes were mainly associated with the regulation of oxytocin, estrogen, ErbB, Wnt, mTOR, PI3K-Akt, growth hormone synthesis/secretion/action, and MAPK signaling pathways. A total of 123 SNPs were selected, based on significantly associated genomic regions and important pathways for SNP genotyping, for verification in 521 additional Bactrian camels. This analysis showed that 13 SNPs were significantly associated with camel milk production yield and 18 SNPs were significantly associated with camel milk composition percentages. Most of these SNPs were located in coding regions of the genome. However, five and two important mutation sites were found in the introns of CSN2 (β-casein) and CSN3 (κ-casein), respectively. Among the candidate genes, NR4A1, ADCY8, PPARG, CSN2, and CSN3 have previously been well studied in dairy livestock. These observations provide a basis for understanding the molecular mechanisms underlying milk production in camels as well as genetic markers for breeding programs aimed at improving milk production.

A Review on Camel Milk Composition, Techno-Functional Properties and Processing Constraints

Camel milk plays a critical role in the diet of peoples belongs to the semi-arid and arid regions. Since prehistoric times, camel milk marketing was limited due to lacking the processing facilities in the camel-rearing areas, nomads practiced the self-consumption of raw and fermented camel milk. A better understanding of the techno-functional properties of camel milk is required for product improvement to address market and customer needs. Despite the superior nutraceutical and health promoting potential, limited camel dairy products are available compared to other bovines. It is a challenging impetus for the dairy industry to provide diversified camel dairy products to consumers with superior nutritional and functional qualities. The physicochemical behavior and characteristics of camel milk is different than the bovine milk, which poses processing and technological challenges. Traditionally camel milk is only processed into various fermented and non-fermented products; however, the production of commercially important dairy products (cheese, butter, yogurt, and milk powder) from camel milk still needs to be processed successfully. Therefore, the industrial processing and transformation of camel milk into various products, including fermented dairy products, pasteurized milk, milk powder, cheese, and other products, require the development of new technologies based on applied research. This review highlights camel milk's processing constraints and techno-functional properties while presenting the challenges associated with processing the milk into various dairy products. Future research directions to improve product quality have also been discussed.

Camelus knoblochi genome reveals the complex evolutionary history of Old World camels

Extant Old World camels (genus Camelus) contributed to the economic and cultural exchanges between the East and West for thousands of years.1,2 Although many remains have been unearthed,3,4,5 we know neither whether the prevalent hybridization observed between extant Camelus species2,6,7 also occurred between extinct lineages and the ancestors of extant Camelus species nor why some populations became extinct while others survived. To investigate these questions, we generated paleogenomic and stable isotope data from an extinct two-humped camel species, Camelus knoblochi. We find that in the mitochondrial phylogeny, all C. knoblochi form a paraphyletic group that nests within the diversity of modern, wild two-humped camels (Camelus ferus). In contrast, they are clearly distinguished from both wild and domesticated (Camelus bactrianus) two-humped camels on the nuclear level. Moreover, the divergence pattern of the three camel species approximates a trifurcation, because the most common topology is only slightly more frequent than the two other possible topologies. This mito-nuclear phylogenetic discordance likely arose due to interspecific gene flow between all three species, suggesting that interspecific hybridization is not exclusive to modern camels but a recurrent phenomenon throughout the evolutionary history of the genus Camelus. These results suggest that the genomic complexity of Old World camels' evolutionary history is underestimated when considering data from only modern species. Finally, we find that C. knoblochi populations began declining prior to the last glacial maximum and, by integrating palaeoecological evidence and stable isotope data, suggest that this was likely due to failure to adapt to a changing environment.

Genetic variability among and within domestic Old and New World camels at the α-lactalbumin gene (LALBA) reveals new alleles and polymorphisms responsible for differential expression

α-Lactalbumin (α-LA), which is encoded by the LALBA gene, is a major whey protein that binds to Ca2+ and facilitates lactose synthesis as a regulatory subunit of the synthase enzyme complex. In addition, it has been shown to play central roles in immune modulation, cell-growth regulation, and antimicrobial activity. In this study, a multitechnical approach was used to fully characterize the LALBA gene and its variants in both coding and regulatory regions for domestic camelids (dromedary, Bactrian camel, alpaca, and llama). The gene analysis revealed a conserved structure among the camelids, but a slight difference in size (2,012 bp on average) due to intronic variations. Promoters were characterized for the transcription factor binding sites (11 found in total). Intraspecies sequence comparison showed 36 SNPs in total (2 in the dromedary, none in the Bactrian camel, 22 in the alpaca, and 12 in the llama), whereas interspecies comparison showed 86 additional polymorphic sites. Eight SNPs were identified as trans-specific polymorphisms, and 2 of them (g.112A>G and g.1229A>G) were particularly interesting in the New World camels. The first creates a new binding site for transcription factor SP1. An enhancing effect of the g.112G variant on the expression was demonstrated by 3 independent pGL3 gene reporter assays. The latter is responsible for the p.78Ile>Val AA replacement and represents novel allelic variants (named LALBA A and B). A link to protein variants has been established by isoelectric focusing (IEF), and bioinformatics analysis revealed that carriers of valine (g.1229G) have a higher glycosylation rate. Genotyping methods based on restriction fragment length polymorphism (PCR-RFLP) were set up for both SNPs. Overall, adenine was more frequent (0.54 and 0.76) at both loci. Four haplotypes were found, and the AA and GA were the most common with a frequency of 0.403 and 0.365, respectively. Conversely, a putative biological gain characterizes the haplotype GG. Therefore, opportunities for rapid directional selection can be realized if this haplotype is associated with favorable milk protein properties. This study adds knowledge at the gene and protein level for α-LA (LALBA) in camelids and importantly contributes to a relatively unexplored research area in these species.

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.

Structure and Dynamics Guiding Design of Antibody Therapeutics and Vaccines

Antibodies and other new antibody-like formats have emerged as one of the most rapidly growing classes of biotherapeutic proteins. Understanding the structural features that drive antibody function and, consequently, their molecular recognition is critical for engineering antibodies. Here, we present the structural architecture of conventional IgG antibodies alongside other formats. We emphasize the importance of considering antibodies as conformational ensembles in solution instead of focusing on single-static structures because their functions and properties are strongly governed by their dynamic nature. Thus, in this review, we provide an overview of the unique structural and dynamic characteristics of antibodies with respect to their antigen recognition, biophysical properties, and effector functions. We highlight the numerous technical advances in antibody structure prediction and design, enabled by the vast number of experimentally determined high-quality structures recorded with cryo-EM, NMR, and X-ray crystallography. Lastly, we assess antibody and vaccine design strategies in the context of structure and dynamics.

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.

Exploiting morphobiometric and genomic variability of African indigenous camel populations-A review

Camels (Camelus dromedarius) in Africa are adapted to arid and the semi-arid environmental conditions, and are valuable for meat, milk and fiber production. On account of the growing demand for camels in this continent, there is a need for knowledge on their phenotypic and genetic diversity. This is fundamental to sustainable herd management and utilization including the design of appropriate breeding and conservation strategies. We reviewed studies on the phenotypic and genetic characterization, breeding objectives, systems of production, productive and reproductive performances, and pathways for the sustainable rearing and use of camels in Africa. The morphological and genetic diversity, productive and reproductive abilities of African camels suggest the existence of genetic variations that can be utilized for breeds/ecotypes' genetic improvement and conservation. Possible areas of intervention include the establishment of open nucleus and community-based breeding schemes and utilization of modern reproductive technologies for the genetic improvement of milk and meat yields, sustainable management of rangelands, capacity building of the pastoralists and agro-pastoralists, institutional supports, formation of centralized conservation centres and efficient and effective marketing systems.

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.

Camel Proteins and Enzymes: A Growing Resource for Functional Evolution and Environmental Adaptation

In less agroecological parts of the Asian, Arabian, and African deserts, Camelus dromedarius play an important role in human survival. For many years, camels have been employed as a source of food, a tool of transportation, and a means of defense. They are becoming increasingly important as viable livestock animals in many desert climates. With the help of camel genetics, genomics and proteomics known so far, this review article will summarize camel enzymes and proteins, which allow them to thrive under varied harsh environmental situations. An in-depth study of the dromedary genome revealed the existence of protein-coding and fast-developing genes that govern a variety of metabolic responses including lipid and protein metabolism, glucoamylase, flavin-containing monooxygenase and guanidinoacetate methyltransferase are other metabolic enzymes found in the small intestine, liver, pancreas, and spleen. In addition, we will discuss the handling of common medications by camel liver cytochrome p 450, which are different from human enzymes. Moreover, camels developed several paths to get optimum levels of trace elements like copper, zinc, selenium, etc., which have key importance in their body for normal regulation of metabolic events. Insulin tolerance, carbohydrate and energy metabolism, xenobiotics metabolizing enzymes, vimentin functions, behavior during the rutting season, resistance to starvation and changes in blood composition and resistance to water loss were among the attractive aspects of camel enzymes and proteins peculiarities in the camels. Resolving the enigma of the method of adaptation and the molecular processes linked with camel life is still a developing repository full of mysteries that need additional exploration.

Comparative analysis of transposable elements provides insights into genome evolution in the genus Camelus

BACKGROUND

Transposable elements (TEs) are common features in eukaryotic genomes that are known to affect genome evolution critically and to play roles in gene regulation. Vertebrate genomes are dominated by TEs, which can reach copy numbers in the hundreds of thousands. To date, details regarding the presence and characteristics of TEs in camelid genomes have not been made available.

RESULTS

We conducted a genome-wide comparative analysis of camelid TEs, focusing on the identification of TEs and elucidation of transposition histories in four species: Camelus dromedarius, C. bactrianus, C. ferus, and Vicugna pacos. Our TE library was created using both de novo structure-based and homology-based searching strategies ( https://github.com/kacst-bioinfo-lab/TE_ideintification_pipeline ). Annotation results indicated a similar proportion of each genomes comprising TEs (35-36%). Class I LTR retrotransposons comprised 16-20% of genomes, and mostly consisted of the endogenous retroviruses (ERVs) groups ERVL, ERVL-MaLR, ERV_classI, and ERV_classII. Non-LTR elements comprised about 12% of genomes and consisted of SINEs (MIRs) and the LINE superfamilies LINE1, LINE2, L3/CR1, and RTE clades. Least represented were the Class II DNA transposons (2%), consisting of hAT-Charlie, TcMar-Tigger, and Helitron elements and comprising about 1-2% of each genome.

CONCLUSIONS

The findings of the present study revealed that the distribution of transposable elements across camelid genomes is approximately similar. This investigation presents a characterization of TE content in four camelid to contribute to developing a better understanding of camelid genome architecture and evolution.

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.

Multiomic analysis of the Arabian camel (Camelus dromedarius) kidney reveals a role for cholesterol in water conservation

The Arabian camel (Camelus dromedarius) is the most important livestock animal in arid and semi-arid regions and provides basic necessities to millions of people. In the current context of climate change, there is renewed interest in the mechanisms that enable camelids to survive in arid conditions. Recent investigations described genomic signatures revealing evolutionary adaptations to desert environments. We now present a comprehensive catalogue of the transcriptomes and proteomes of the dromedary kidney and describe how gene expression is modulated as a consequence of chronic dehydration and acute rehydration. Our analyses suggested an enrichment of the cholesterol biosynthetic process and an overrepresentation of categories related to ion transport. Thus, we further validated differentially expressed genes with known roles in water conservation which are affected by changes in cholesterol levels. Our datasets suggest that suppression of cholesterol biosynthesis may facilitate water retention in the kidney by indirectly facilitating the AQP2-mediated water reabsorption.

Intra-horn Penile Intromission in the Alpaca Vicugna pacos and Consequences to Genital Morphology

Copulatory behavior and genital morphology interact to deliver sperm more effectively during mating, but the nature of this interaction has not been explored in depth in most vertebrates. Alpacas have unusually long copulations lasting 15-20 min, and a unique copulatory behavior, where the penis intromits all the way past the cervix, into the uterine horns. Here we describe the morphology of male and female genitalia and report unique morphological characteristics that may be associated with this unusual insemination mode. Vaginal shape is highly variable, and seemingly not associated with age or parity. The cranial vagina varies between bulbous and cylindrical, while the caudal vagina is typically narrower and always cylindrical. The cervix consists of a series of two to three spirals or rings, and it is often found in a relaxed state that may prevent damage caused by the cartilaginous penis tip as it pushes through the cervix to reach the uterine horns. The uterus and uterine horns have a complex shape with multiple constrictions. The cartilaginous penis tip has a sharp urethral process that may help to push against these constrictions. The diameter of the vaginal lumen is much greater than the diameter of the penis suggesting that there is little direct interaction between them, and that female vaginal shape is not under strong copulatory selection. In effect, the entire female reproductive tract of the female is interacting with the penis during copulation.

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.

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.

Genome-Wide Identification and Analysis of Variants in Domestic and Wild Bactrian Camels Using Whole-Genome Sequencing Data

The population size of Bactrian camels is smaller than dromedary, and they are distributed in cold and mountain regions and are also at the risk of extinction in some countries such as Iran. To identify and investigate the genome-wide variations, whole-genome sequencing of two Iranian Bactrian camels were performed with 37.4- and 42.6-fold coverage for the first time. Along with Iranian Bactrian camels, sequencing data from two Mongolian domestic and two wild Bactrian camels deposited in the NCBI were reanalyzed. The analysis eventuated to the identification of 4,908,998, 4,485,725, and 4,706,654 SNPs for Iranian, Mongolian domestic, and wild Bactrian camels, respectively. Also, INDEL variations ranged from 358,311 to 533,188 in all six camels. Results of variants annotation in all samples revealed that more than 88 percent of SNPs and INDELs were located in the intergenic and intronic regions. We found that 800,530 SNPs were common among all studied camels, containing 4,046 missense variants that affected 2,428 genes. Investigation of common genes among all camels containing the missense SNPs showed that there are 98 zinc finger and 4 fertility-related genes (ZP1, ZP2, ZP4, and ZPBP) in this set.

Genome-wide diversity and global migration patterns in dromedaries follow ancient caravan routes

Dromedaries have been essential for the prosperity of civilizations in arid environments and the dispersal of humans, goods and cultures along ancient, cross-continental trading routes. With increasing desertification their importance as livestock species is rising rapidly, but little is known about their genome-wide diversity and demographic history. As previous studies using few nuclear markers found weak phylogeographic structure, here we detected fine-scale population differentiation in dromedaries across Asia and Africa by adopting a genome-wide approach. Global patterns of effective migration rates revealed pathways of dispersal after domestication, following historic caravan routes like the Silk and Incense Roads. Our results show that a Pleistocene bottleneck and Medieval expansions during the rise of the Ottoman empire have shaped genome-wide diversity in modern dromedaries. By understanding subtle population structure we recognize the value of small, locally adapted populations and appeal for securing genomic diversity for a sustainable utilization of this key desert species.

The major histocompatibility complex of Old World camelids: Class I and class I-related genes

The genomic structure of the Major Histocompatibility Complex (MHC) region and variation in selected MHC class I related genes in Old World camels, Camelus bactrianus and Camelus dromedaries were studied. The overall genomic organization of the camel MHC region follows a general pattern observed in other mammalian species and individual MHC loci appear to be well conserved. Selected MHC class I genes B-67 and BL3-7 exhibited unexpectedly low variability, even when compared to other camel MHC class I related genes MR1 and MICA. Interspecific SNP and allele sharing are relatively common, and frequencies of heterozygotes are usually low. Such a low variation in a genomic region generally considered as one of the most polymorphic in vertebrate genomes is unusual. Evolutionary relationships between MHC class I related genes and their counterparts from other species seem to be rather complex. Often, they do not follow the general evolutionary history of the species concerned. Close evolutionary relationships of individual MHC class I loci between camels, humans and dogs were observed. Based on the results of this study and on our data on MHC class II genes, the extent and the pattern of polymorphism of the MHC region of Old World camelids differed from most mammalian groups studied so far. Camels thus seem to be an important model for our understanding of the role of genetic diversity in immune functions, especially in the context of unique features of their immunoglobulin and T-cell receptor genes.

Effect of Research Impact on Emerging Camel Husbandry, Welfare and Social-Related Awareness

The lack of applied scientific research on camels, despite them being recognized as production animals, compels the reorganization of emerging camel breeding systems with the aim of achieving successful camel welfare management strategies all over the world. Relevant and properly-framed research widely impacts dissemination of scientific contents and drives public willingness to enhance ethically acceptable conditions for domestic animals. Consumer perception of this livestock industry will improve and high-quality products will be obtained. This paper draws on bibliometric indicators as promoting factors for camel-related research advances, tracing historical scientific publications indexed in ScienceDirect directory from 1880-2019. Camel as a species did not affect Journal Citation Reports (JCR) impact (p > 0.05) despite the journal, author number, corresponding author origin, discipline and publication year affecting it (p < 0.001). Countries with traditionally well-established camel farming are also responsible for the papers with the highest academic impact. However, camel research advances may have only locally and partially influenced welfare related laws, so intentional harming acts and basic needs neglect may persist in these species. A sustainable camel industry requires those involved in camel research to influence business stakeholders and animal welfare advocacies by highlighting the benefits of camel wellbeing promotion, co-innovation partnership establishment and urgent enhancement of policy reform.

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

Camelids are characterized by their unique adaptive immune system that exhibits the generation of homodimeric heavy-chain immunoglobulins, somatic hypermutation of T-cell receptors, and low genetic diversity of major histocompatibility complex (MHC) genes. However, short-read assemblies are typically highly fragmented in these gene loci owing to their repetitive and polymorphic nature. Here, we constructed a chromosome-level assembly of wild Bactrian camel genome based on high-coverage long-read sequencing and chromatin interaction mapping. The assembly with a contig N50 of 5.37 Mb and a scaffold N50 of 76.03 Mb, represents the most contiguous camelid genome to date. The genomic organization of immunoglobulin heavy-chain locus was similar between the wild Bactrian camel and alpaca, and genes encoding for conventional and heavy-chain antibodies were intermixed. The organizations of two immunoglobulin light-chain loci and four T cell receptor loci were also fully deciphered using the new assembly. Additionally, the complete classical MHC region was resolved into a single contig. The high-quality assembly presented here provides an essential reference for future investigations examining the camelid immune system.

Whole-genome sequencing of 128 camels across Asia reveals origin and migration of domestic Bactrian camels

The domestic Bactrian camels were treated as one of the principal means of locomotion between the eastern and western cultures in history. However, whether they originated from East Asia or Central Asia remains elusive. To address this question, we perform whole-genome sequencing of 128 camels across Asia. The extant wild and domestic Bactrian camels show remarkable genetic divergence, as they were split from dromedaries. The wild Bactrian camels also contribute little to the ancestry of domestic ones, although they share close habitat in East Asia. Interestingly, among the domestic Bactrian camels, those from Iran exhibit the largest genetic distance and the earliest split from all others in the phylogeny, despite evident admixture between domestic Bactrian camels and dromedaries living in Central Asia. Taken together, our study support the Central Asian origin of domestic Bactrian camels, which were then immigrated eastward to Mongolia where native wild Bactrian camels inhabit.

Ming, Yuan et al. performed whole-genome sequencing on 128 wild and domesticated Bactrian camels across Asia. They show that wild and domestic Bactrian camels are genetically diverged from dromedaries, and that wild camels contributed little to domestic camel ancestry despite sharing a habitat in East Asia.

Zoonotic parasites of dromedary camels: so important, so ignored

With a global population of about 35 million in 47 countries, dromedary camels play a crucial role in the economy of many marginal, desert areas of the world where they survive under harsh conditions. Nonetheless, there is scarce knowledge regarding camels' parasite fauna which can reduce their milk and meat productions. In addition, only scattered information is available about zoonotic parasites transmitted to humans via contamination (e.g. Cryptosporidium spp., Giardia duodenalis, Balantidium coli, Blastocystis spp. and Enterocytozoon bieneusi), as foodborne infections (e.g. Toxoplasma gondii, Trichinella spp. and Linguatula serrata) or by arthropod vectors (Trypanosoma spp.). Herein, we draw attention of the scientific community and health policy-making organizations to the role camels play in the epidemiology of parasitic zoonotic diseases also in the view of an increase in their farming in desert areas worldwide.

The Major Histocompatibility Complex of Old World Camels-A Synopsis

This study brings new information on major histocompatibility complex (MHC) class III sub-region genes in Old World camels and integrates current knowledge of the MHC region into a comprehensive overview for Old World camels. Out of the MHC class III genes characterized, TNFA and the LY6 gene family showed high levels of conservation, characteristic for MHC class III loci in general. For comparison, an MHC class II gene TAP1, not coding for antigen presenting molecules but functionally related to MHC antigen presenting functions was studied. TAP1 had many SNPs, even higher than the MHC class I and II genes encoding antigen presenting molecules. Based on this knowledge and using new camel genomic resources, we constructed an improved genomic map of the entire MHC region of Old World camels. The MHC class III sub-region shows a standard organization similar to that of pig or cattle. The overall genomic structure of the camel MHC is more similar to pig MHC than to cattle MHC. This conclusion is supported by differences in the organization of the MHC class II sub-region, absence of functional DY genes, different organization of MIC genes in the MHC class I sub-region, and generally closer evolutionary relationships of camel and porcine MHC gene sequences analyzed so far.

From Desert to Medicine: A Review of Camel Genomics and Therapeutic Products

Camels have an important role in the lives of human beings, especially in arid regions, due to their multipurpose role and unique ability to adapt to harsh conditions. In spite of its enormous economic, cultural, and biological importance, the camel genome has not been widely studied. The size of camel genome is roughly 2.38 GB, containing over 20,000 genes. The unusual genetic makeup of the camel is the main reason behind its ability to survive under extreme environmental conditions. The camel genome harbors several unique variations which are being investigated for the treatment of several disorders. Various natural products from camels have also been tested and prescribed as adjunct therapy to control the progression of ailments. Interestingly, the camel employs unique immunological and molecular mechanisms against pathogenic agents and pathological conditions. Here, we broadly review camel classification, distribution and breed as well as recent progress in the determination of the camel genome, its size, genetic distribution, response to various physiological conditions, immunogenetics and the medicinal potential of camel gene products.

Cdrom Archive: A Gateway to Study Camel Phenotypes

Camels are livestock that exhibit unique morphological, biochemical, and behavioral traits, which arose by natural and artificial selection. Investigating the molecular basis of camel traits has been limited by: (1) the absence of a comprehensive record of morphological trait variation (e.g., diseases) and the associated mode of inheritance, (2) the lack of extended pedigrees of specific trait(s), and (3) the long reproductive cycle of the camel, which makes the cost of establishing and maintaining a breeding colony (i.e., monitoring crosses) prohibitively high. Overcoming these challenges requires (1) detailed documentation of phenotypes/genetic diseases and their likely mode of inheritance (and collection of related DNA samples), (2) conducting association studies to identify phenotypes/genetic diseases causing genetic variants (instead of classical linkage analysis, which requires extended pedigrees), and (3) validating likely causative variants by screening a large number of camel samples from different populations. We attempt to address these issues by establishing a systematic way of collecting camel DNA samples, and associated phenotypic information, which we call the "Cdrom Archive." Here, we outline the process of building this archive to introduce it to other camel researchers (as an example). Additionally, we discuss the use of this archive to study the phenotypic traits of Arabian Peninsula camel breeds (the "Mezayen" camels). Using the Cdrom Archive, we report variable phenotypic traits related to the coat (color, length, and texture), ear and tail lengths, along with other morphological measurements.

Quality and quantity of dromedary camel DNA sampled from whole-blood, saliva, and tail-hair

Camels are livestock with unique adaptations to hot-arid regions. To effectively study camel traits, a biobank of camel DNA specimens with associated biological information is needed. We examined whole-blood, saliva (buccal swabs), and tail-hair follicle samples to determine which is the best source for establishing a DNA biobank. We inspected five amounts of each of whole-blood, buccal swabs, and tail-hair follicles in nine camels, both qualitatively via gel electrophoresis and quantitatively using a NanoDrop spectrophotometer. We also tested the effects of long term-storage on the quality and quantity of DNA, and measured the rate of degradation, by analyzing three buccal swab samples and 30 tail-hair follicles over a period of nine months. Good quality DNA, in the form of visible large size DNA bands, was extracted from all three sources, for all five amounts. The five volumes of whole-blood samples (20–100μl) provided ~0.4–3.6 μg, the five quantities of buccal swabs (1–5) produced ~0.1–12 μg, while the five amounts of tail-hair follicles (10–50) resulted in ~0.7–25 μg. No differences in the rate of degradation of buccal swab and tail-hair follicle DNA were detected, but there was clearly greater deterioration in the quality of DNA extracted from buccal swabs when compared to tail-hair follicles. We recommend using tail-hair samples for camel DNA biobanking, because it resulted in both an adequate quality and quantity of DNA, along with its ease of collection, transportation, and storage. Compared to its success in studies of other domesticated animals, we anticipate that using ~50 tail-hair follicles will provide sufficient DNA for sequencing or SNP genotyping.

Genetic characterization of the oxytocin-neurophysin I gene (OXT) and its regulatory regions analysis in domestic Old and New World camelids

Oxytocin is a neurohypophysial peptide linked to a wide range of biological functions, including milk ejection, temperament and reproduction. Aims of the present study were a) the characterization of the OXT (Oxytocin-neurophysin I) gene and its regulatory regions in Old and New world camelids; b) the investigation of the genetic diversity and the discovery of markers potentially affecting the gene regulation. On average, the gene extends over 814 bp, ranging between 825 bp in dromedary, 811 bp in Bactrian and 810 bp in llama and alpaca. Such difference in size is due to a duplication event of 21 bp in dromedary. The main regulatory elements, including the composite hormone response elements (CHREs), were identified in the promoter, whereas the presence of mature microRNAs binding sequences in the 3'UTR improves the knowledge on the factors putatively involved in the OXT gene regulation, although their specific biological effect needs to be still elucidated. The sequencing of genomic DNA allowed the identification of 17 intraspecific polymorphisms and 69 nucleotide differences among the four species. One of these (MF464535:g.622C>G) is responsible, in alpaca, for the loss of a consensus sequence for the transcription factor SP1. Furthermore, the same SNP falls within a CpG island and it creates a new methylation site, thus opening future possibilities of investigation to verify the influence of the novel allelic variant in the OXT gene regulation. A PCR-RFLP method was setup for the genotyping and the frequency of the allele C was 0.93 in a population of 71 alpacas. The obtained data clarify the structure of OXT gene in domestic camelids and add knowledge to the genetic variability of a genomic region, which has received little investigation so far. These findings open the opportunity for new investigations, including association studies with productive and reproductive traits.

Mitogenome Sequencing in the Genus Camelus Reveals Evidence for Purifying Selection and Long-term Divergence between Wild and Domestic Bactrian Camels

The genus Camelus is an interesting model to study adaptive evolution in the mitochondrial genome, as the three extant Old World camel species inhabit hot and low-altitude as well as cold and high-altitude deserts. We sequenced 24 camel mitogenomes and combined them with three previously published sequences to study the role of natural selection under different environmental pressure, and to advance our understanding of the evolutionary history of the genus Camelus. We confirmed the heterogeneity of divergence across different components of the electron transport system. Lineage-specific analysis of mitochondrial protein evolution revealed a significant effect of purifying selection in the concatenated protein-coding genes in domestic Bactrian camels. The estimated dN/dS < 1 in the concatenated protein-coding genes suggested purifying selection as driving force for shaping mitogenome diversity in camels. Additional analyses of the functional divergence in amino acid changes between species-specific lineages indicated fixed substitutions in various genes, with radical effects on the physicochemical properties of the protein products. The evolutionary time estimates revealed a divergence between domestic and wild Bactrian camels around 1.1 [0.58-1.8] million years ago (mya). This has major implications for the conservation and management of the critically endangered wild species, Camelus ferus.

Evaluation of a panel of antibodies for the immunohistochemical identification of immune cells in paraffin-embedded lymphoid tissues of new- and old-world camelids

Different species of camelids play an important role in the epidemiology of various emerging infectious diseases such as Middle East respiratory syndrome. For precise investigations of the immunopathogenesis in these host species, appropriate immunohistochemical markers are highly needed in order to phenotype distinct immune cells populations in camelids. So far, specific immunohistochemical markers for camelid immune cells are rarely commercially available, and cross-reactivity studies are restricted to the use of frozen dromedary tissues. To bridge this gap, 14 commercially available primary antibodies were tested for their suitability to demonstrate immune cell populations on formalin fixed paraffin-embedded (FFPE) tissue sections of dromedaries, Bactrian camels, llamas, and alpacas in the present study. Out of these, 9 antibodies directed against CD3, CD20, CD79α, HLA-DR, Iba-1, myeloid/histiocyte antigen, CD204, CD208, and CD68 antigen exhibited distinct immunoreaction patterns to certain camelid immune cell subsets. The distribution of these antigens was comparatively evaluated in different anatomical compartments of thymus, spleen, mesenteric, and tracheobronchial lymph nodes. The presented results will provide a basis for further investigations in camelids, especially with respect to the role of the immune response in certain infectious diseases, which harbor a considerable risk to spill over to other species.