Search

BackgroundThe mammalian Natural Killer Complex (NKC) harbors genes and gene families encoding a variety of C-type lectin-like proteins expressed on various immune cells. The NKC is a complex genomic region well-characterized in mice, humans and domestic animals. The major limitations of automatic annotation of the NKC in non-model animals include short-read based sequencing, methods of assembling highly homologous and repetitive sequences, orthologues missing from reference databases and weak expression. In this situation, manual annotations of complex genomic regions are necessary.MethodsThis study presents a manual annotation of the genomic structure of the NKC region in a high-quality reference genome of the domestic cat and compares it with other felid species and with representatives of other carnivore families. Reference genomes of Carnivora, irrespective of sequencing and assembly methods, were screened by BLAST to retrieve information on their killer cell lectin-like receptor (KLR) gene content. Phylogenetic analysis of in silico translated proteins of expanded subfamilies was carried out.ResultsThe overall genomic structure of the NKC in Carnivora is rather conservative in terms of its C-type lectin receptor gene content. A novel KLRH-like gene subfamily (KLRL) was identified in all Carnivora and a novel KLRJ-like gene was annotated in the Mustelidae. In all six families studied, one subfamily (KLRC) expanded and experienced pseudogenization. The KLRH gene subfamily expanded in all carnivore families except the Canidae. The KLRL gene subfamily expanded in carnivore families except the Felidae and Canidae, and in the Canidae it eroded to fragments.ConclusionsKnowledge of the genomic structure and gene content of the NKC region is a prerequisite for accurate annotations of newly sequenced genomes, especially of endangered wildlife species. Identification of expressed genes, pseudogenes and gene fragments in the context of expanded gene families would allow the assessment of functionally important variability in particular species.

Abstract The anthropogenic impact on wildlife is ever increasing. With shrinking habitats, wild populations are being pushed to co-exist in proximity to humans leading to an increased threat of infectious diseases. Therefore, understanding the immune system of a species is key to assess its resilience in a changing environment. The innate immune system (IIS) is the body’s first line of defense against pathogens. High variability in IIS genes, like toll-like receptor (TLR) genes, appears to be associated with resistance to infectious diseases. However, few studies have investigated diversity in TLR genes in vulnerable species for conservation. Large predators are threatened globally including leopards and cheetahs, both listed as 'vulnerable' by IUCN. To examine IIS diversity in these sympatric species, we used next-generation-sequencing to compare selected TLR genes in African leopards and cheetahs. Despite differences, both species show some TLR haplotype similarity. Historic cheetahs from all subspecies exhibit greater genetic diversity than modern Southern African cheetahs. The diversity in investigated TLR genes is lower in modern Southern African cheetahs than in African leopards. Compared to historic cheetah data and other subspecies, a more recent population decline might explain the observed genetic impoverishment of TLR genes in modern Southern African cheetahs. However, this may not yet impact the health of this cheetah subspecies.

Abstract Background Wohlfahrtia magnifica is an obligatory parasite that causes myiasis in several warm-blooded vertebrates. Adult females deposit the first-stage larvae directly onto wounds or natural body orifices (e.g., genitalia) of the host, from where they quickly colonize the host tissue and feed on it for development. The infestation of W. magnifica can lead to health issues, welfare concerns, and substantial economic losses. To date, little is known about the molecular mechanisms of the W. magnifica-causing myiasis. Results In this study, we collected parasitic-stage larvae of W. magnifica from wounds of naturally infested Bactrian camels, as well as pupae and adult flies reared in vitro from the wound-collected larvae, for investigating the gene expression profiles of the different developmental stages of W. magnifica, with a particular focus on examining gene families closely related to the parasitism of the wound-collected larvae. As key proteins related to the parasite-host interaction, 2049 excretory/secretory (ES) proteins were identified in W. magnifica through the integration of multiple bioinformatics approaches. Functional analysis indicates that these ES proteins are primarily involved in cuticle development, peptidase activity, immune response, and metabolic processes. The global investigation of gene expression at different developmental stages using pairwise comparisons and weighted correlation network analysis (WGCNA) showed that the upregulated genes during second-stage larvae were related to cuticle development, peptidase activity, and RNA transcription and translation; during third-stage larvae to peptidase inhibitor activity and nutrient reservoir activity; during pupae to cell and tissue morphogenesis and cell and tissue development; and during adult flies to signal perception, many of them involved in light perception, and adult behavior, e.g., feeding, mating, and locomotion. Specifically, the expression level analysis of the likely parasitism-related genes in parasitic wound-collected larvae revealed a significant upregulation of 88 peptidase genes (including 47 serine peptidase genes), 110 cuticle protein genes, and 21 heat shock protein (hsp) genes. Interestingly, the expression of 2 antimicrobial peptide (AMP) genes, including 1 defensin and 1 diptericin, was also upregulated in the parasitic larvae. Conclusions We identified ES proteins in W. magnifica and investigated their functional distribution. In addition, gene expression profiles at different developmental stages of W. magnifica were examined. Specifically, we focused on gene families closely related to parasitism of wound-collected larvae. These findings shed light on the molecular mechanisms underlying the life cycle of the myiasis-causing fly, especially during the parasitic larval stages, and provide guidance for the development of control measures against W. magnifica.

IntroductionThe two-humped Bactrian camel (Camelus bactrianus) is a large, even-toed ungulate native to the steppes of Central Asia. Domestic Bactrian camels are economically important in Mongolia and other Central Asian countries. These animals are used for transport, milk and meat production, and camel racing which is a great culture of nomads. Eimeriosis, also known as coccidiosis, is considered as an economically important parasitic diseases in Bactrian camels. There is still considerable lack of data concerning the spectrum of monoxenous Eimeria species, their epizootiology as well as their precise life cycles in Bactrian camels. This study was performed to determine the prevalence of Eimeria species in camelids from southern part of Mongolia.MethodsA total of 536 fresh camel fecal samples (n = 536) collected from herds located in five different Aimags (provinces) of Mongolia were examined. Eimeria spp. oocysts were isolated using the sugar flotation technique, and after sporulation, oocysts were identified by morphometric evaluation.ResultsWe identified the most common Eimeria species infecting Mongolian Bactrian camels: Eimeria cameli (22.3%), Eimeria rajasthani (37.3%) and Eimeria dromedarii (27.7%). Interestingly, mixed infections were detected in 24.8% (n = 133) of the samples, while 39.0% (n = 209) were negative for coccidian stages. To investigate the immunogenetic response of the Mongolian Bactrian camels to Eimeria spp. infection, we screened the genetic diversity in a functional important immune response gene of the major histocompatibility complex (MHC). We detected two polymorphic sites in the MHC class II DRA exon 2, which translated into one non-synonymous and one synonymous amino acid (aa) change.DiscussionThe resulting aa alleles were not significantly associated with any of the three detected Eimeria species infections, nor could we show heterozygote advantage in non-infected Mongolian Bactrian camels. Further investigations on molecular epidemiology, in vitro culture, pathogenicity and host–parasite interactions will be necessary to better understand the impact of eimeriosis in Bactrian camels.

Common names allow species diversity to be acknowledged by experts and non-specialists alike; they are descriptors with both scientific and cultural implications. However, a lack of clarity when using a common name could risk altering perceptions of threatened species. This is the case for the Critically Endangered wild camel Camelus ferus, which, despite extensive evidence of its species status, is frequently referred to in English as wild Bactrian camel. However, the wild camel (Mongolian: хавтгай, khavtgai; Chinese: 野骆驼, ye luo tuo) is not a wild version of the domestic Bactrian camel Camelus bactrianus but a separate species near extinction, with an estimated population of c. 950. Failure to clearly separate Bactrian and wild camels in name risks masking the plight of the few remaining wild camels with the visible abundance of the domesticated species. Here we advocate the use of the accurate English common name wild camel for C. ferus ideally alongside its Indigenous names to correctly represent its cultural and conservation importance.

DNA methylation is involved in regulatory processes throughout the animal kingdom. Here, the authors map DNA methylation in 535 vertebrates and 45 invertebrates, establishing a reference dataset for cross-species analysis and exploring epigenetic variation across vertebrate evolution.

Comparative transcriptomic analyses between dromedary camels and rats provide insight into water homeostasis gene responses in animals under water stress conditions.

BackgroundThe mammalian Leukocyte Receptor Complex (LRC) chromosomal region may contain gene families for the killer cell immunoglobulin-like receptor (KIR) and/or leukocyte immunoglobulin-like receptor (LILR) collections as well as various framing genes. This complex region is well described in humans, mice, and some domestic animals. Although single KIR genes are known in some Carnivora, their complements of LILR genes remain largely unknown due to obstacles in the assembly of regions of high homology in short-read based genomes.MethodsAs part of the analysis of felid immunogenomes, this study focuses on the search for LRC genes in reference genomes and the annotation of LILR genes in Felidae. Chromosome-level genomes based on single-molecule long-read sequencing were preferentially sought and compared to representatives of the Carnivora.ResultsSeven putatively functional LILR genes were found across the Felidae and in the Californian sea lion, four to five genes in Canidae, and four to nine genes in Mustelidae. They form two lineages, as seen in the Bovidae. The ratio of functional genes for activating LILRs to inhibitory LILRs is slightly in favor of inhibitory genes in the Felidae and the Canidae; the reverse is seen in the Californian sea lion. This ratio is even in all of the Mustelidae except the Eurasian otter, which has a predominance of activating LILRs. Various numbers of LILR pseudogenes were identified.ConclusionsThe structure of the LRC is rather conservative in felids and the other Carnivora studied. The LILR sub-region is conserved within the Felidae and has slight differences in the Canidae, but it has taken various evolutionary paths in the Mustelidae. Overall, the process of pseudogenization of LILR genes seems to be more frequent for activating receptors. Phylogenetic analysis found no direct orthologues across the Carnivora which corroborate the rapid evolution of LILRs seen in mammals.

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.

Myiasis caused by Wohlfahrtia magnifica is a widespread parasitic infestation in mammals. The infested host suffers from damage as the developing larvae feed on its tissues. For the control of myiasis infestation, genetic methods have been shown to be effective and promising as an alternative to insecticides. Combining genome, isoform sequencing (Iso-Seq), and RNA sequencing (RNA-seq) data, we isolated and characterized two sex-determination genes, W. magnifica transformer (Wmtra) and W. magnifica transformer2 (Wmtra2), whose orthologs in a number of insect pests have been utilized to develop genetic control approaches. Wmtra transcripts are sex-specifically spliced; only the female transcript encodes a full-length functional protein, while the male transcript encodes a truncated and non-functional polypeptide due to the presence of the male-specific exon containing multiple in-frame stop codons. The existence of five predicted TRA/TRA2 binding sites in the male-specific exon and the surrounding intron of Wmtra, as well as the presence of an RNA-recognition motif in WmTRA2 may suggest the auto-regulation of Wmtra by its own protein interacting with WmTRA2. This results in the skipping of the male-specific exon and translation of the full-length functional protein only in females. Our comparative study in dipteran species showed that both the WmTRA and WmTRA2 proteins exhibit a high degree of similarity to their orthologs in the myiasis-causing blow flies. Additionally, transcriptome profiling performed between adult females and adult males reported 657 upregulated and 365 downregulated genes. Functional analysis showed that among upregulated genes those related to meiosis and mitosis Gene Ontology (GO) terms were enriched, while, among downregulated genes, those related to muscle cell development and aerobic metabolic processes were enriched. Among the female-biased gene set, we detected five candidate genes, vasa (vas), nanos (nanos), bicoid (bcd), Bicaudal C (BicC), and innexin5 (inx5). The promoters of these genes may be able to upregulate Cas9 expression in the germline in Cas9-based homing gene drive systems as established in some flies and mosquitoes. The isolation and characterization of these genes is an important step toward the development of genetic control programs against W. magnifica infestation.

Fernando Alvira-Iraizoz and Benjamin T. Gillard et al. comprehensively catalogue the kidney transcriptome and proteome of dromedary camels in response to chronic dehydration stress and subsequent acute rehydration. Their results suggest that the suppression of the cholesterol biosynthesis may facilitate water retention in the dromedary camel kidney, providing further insight into how these animals survive in arid conditions.

In Burkina Faso, goats are the second most numerous ruminant livestock population, with almost exclusively indigenous breeds being reared in extensive production systems in various agroecological zones. This study was carried out to understand the morphological variation of local goat breeds in the Sudano-Sahelian and Sudanian agroecological zones. A total of 511 adult female animals belonging to two presumed populations (Mossi breed in Sudano-Sahelian zone and Djallonké breed in Sudanian zone) were sampled and body weight as well as a range of linear body measurements, following FAO guidelines, were recorded. The least squares means of body measurements of indicated that Sudano-Sahelian goats have significantly (p < 0.001) larger body measurements than Sudanian goats. Furthermore, relative high variability of the two populations in morphometric traits was observed. Principal Component Analysis (PCA) suggested structure between Mossi breed on one side and Djallonké on the other side, but no strict separation was observed, suggesting that gene flow is occurring among the different populations. A dispersion map with four clusters was built based on the first two factors. The least square means of body measurements ranked the four groups from small to large body size, namely Djallonké, Mossi × Djallonké, Mossi, and Sahelian × Mossi. Gene flow from Sahelian goat into other populations of the country, based on migration of the Fulani ethnic group from the Sahel into areas with Mossi and Djallonké breeds, could explain this configuration and confirms the continuous erosion of genetic identity of these two local breeds. The sustainable use of these adapted local goat genetic resources calls for the promotion of sustainable genetic improvement using participatory breeding approaches.

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

This review summarizes the current knowledge on the major histocompatibility complex (MHC) of the family Felidae. This family comprises an important domestic species, the cat, as well as a variety of free-living felids, including several endangered species. As such, the Felidae have the potential to be an informative model for studying different aspects of the biological functions of MHC genes, such as their role in disease mechanisms and adaptation to different environments, as well as the importance of genetic diversity for conservation issues in free-ranging or captive populations. Despite this potential, the current knowledge on the MHC in the family as a whole is fragmentary and based mostly on studies of the domestic cat and selected species of big cats. The overall structure of the domestic cat MHC is similar to other mammalian MHCs following the general scheme “centromere-MHC class I-MHC class III-MHC class II” with some differences in the gene contents. An unambiguously defined orthologue of the non-classical class I HLA-E gene has not been identified so far and the class II DQ and DP genes are missing or pseudogenized, respectively. A comparison with available genomes of other felids showed a generally high level of structural and sequence conservation of the MHC region. Very little and fragmentary information on in vitro and/or in vivo biological functions of felid MHC genes is available. So far, no association studies have indicated effects of MHC genetic diversity on a particular disease. No information is available on the role of MHC class I molecules in interactions with Natural Killer (NK) cell receptors or on the putative evolutionary interactions (co-evolution) of the underlying genes. A comparison of complex genomic regions encoding NK cell receptors (the Leukocyte Receptor Complex, LRC and the Natural Killer Cell Complex, NKC) in the available felid genomes showed a higher variability in the NKC compared to the LRC and the MHC regions. Studies of the genetic diversity of domestic cat populations and/or specific breeds have focused mainly on DRB genes. Not surprisingly, higher levels of MHC diversity were observed in stray cats compared to pure breeds, as evaluated by DRB sequencing as well as by MHC-linked microsatellite typing. Immunogenetic analysis in wild felids has only been performed on MHC class I and II loci in tigers, Namibian leopards and cheetahs. This information is important as part of current conservation tasks to assess the adaptive potential of endangered wild species at the human-wildlife interface, which will be essential for preserving biodiversity in a functional ecosystem.

The worldwide dromedary milk production has increased sharply since the beginning of this century due to prolonged shelf life, improved food-safety and perceived health benefits. Scientific confirmation of health claims will expand the market of dromedary milk further. As a result, more and more dromedaries will be bred for one purpose only: the highest possible milk production. However, intensive dromedary farming systems have consequences for animal welfare and may lead to genetic changes. Tighter regulations will be implemented to restrict commercialization of raw milk. Protocols controlling welfare of dromedaries and gene databases of milk-dromedaries will prevent negative consequences of intensive farming. In countries where dromedaries have only recently been introduced as production animal, legislators have limited expertise on this species. This is exemplified by an assessment on behalf of the Dutch government, recommending prohibiting keeping this species from 2024 onwards because the dromedary was deemed to be insufficiently domesticated. Implementation of this recommendation in Dutch law would have devastating effects on existing dromedary farms and could also pave the way for adopting similar measures in other European countries. In this paper it is shown that the Dutch assessment lacks scientific rigor. Awareness of breeders and legislators for the increasing knowledge about dromedaries and their products would strengthen the position of dromedaries as one of the most adapted and sustainable animals.

The genomes of crossbred (admixed) individuals are a mosaic of ancestral haplotypes formed by recombination in each generation. The proportion of these ancestral haplotypes in certain genomic regions can be responsible for either susceptibility or tolerance against pathogens, and for performances in production traits. Using a medium-density genomic marker panel from the Illumina Bovine SNP50 BeadChip, we estimated individual admixture proportions for Baoulé x Zebu crossbred cattle in Burkina Faso, which were tested for trypanosome infection by direct ELISA from blood samples. Furthermore, we calculated local ancestry deviation from average for each SNP across 29 autosomes to identify potential regions under selection in the trypanotolerant Baoulé cattle and their crossbreds. We identified significant deviation from the local average ancestry (above 5 and 10% genome-wide thresholds) on chromosomes 8 and 19 in the positive animals, while the negative ones showed higher deviation on chromosomes 6, 19, 21, and 22. Some candidate genes on chromosome 6 (PDGFRA) and chromosome 19 (CDC6) have been found associated to trypanotolerance in West African taurines. Screening for FST outliers in trypanosome positive/negative animals we detected seven variants putatively under selection. Finally, we identified a minimum set of highly ancestry informative markers for routine admixture testing. The results of this study contribute to a better understanding of the genetic basis of trypanotolerance in Baoulé cattle and their crossbreeds. Furthermore, we provide a small informative marker set to monitor admixture in this valuable indigenous breed. As such, our results are important for conserving the genetic uniqueness and trypanotolerance of Baoulé cattle, as well as for the improvement of Baoulé and Zebu crossbreds in specific community-based breeding programs.

In this study, single-SNP GWAS analyses were conducted to find regions affecting tolerance against trypanosomosis and morphometrics traits in purebred and crossbred Baoulé cattle of Burkina Faso. The trypanosomosis status (positive and negative) and a wide set of morphological traits were recorded for purebred Baoulé and crossbred Zebu x Baoulé cattle, and genotyped with the Illumina Bovine SNP50 BeadChip. After quality control, 36,203 SNPs and 619 animals including 343 purebred Baoulé and 279 crossbreds were used for the GWAS analyses. Several important genes were found that can influence morphological parameters. Although there were no genes identified with a reported strong connection to size traits, many of them were previously identified in various growth-related studies. A re-occurring theme for the genes residing in the regions identified by the most significant SNPs was pleiotropic effect on growth of the body and the cardiovascular system. Regarding trypanosomosis tolerance, two potentially important regions were identified in purebred Baoulé on chromosomes 16 and 24, containing the CFH, CRBN, TRNT1 and, IL5RA genes, and one additional genomic region in Baoulé, x Zebu crossbreds on chromosome 5, containing MGAT4C and NTS. Almost all of these regions and genes were previously related to the trait of interest, while the CRBN gene was to our knowledge presented in the context of trypanosomiasis tolerance for the first time.

Dromedaries are an important livestock, used as beasts of burden and for meat and milk production. However, they can act as an intermediate source or vector for transmitting zoonotic viruses to humans, such as the Middle East respiratory syndrome coronavirus (MERS-CoV) or Crimean–Congo hemorrhagic fever virus (CCHFV). After several outbreaks of CCHFV in the Arabian Peninsula, recent studies have demonstrated that CCHFV is endemic in dromedaries and camel ticks in the United Arab Emirates (UAE). There is no apparent disease in dromedaries after the bite of infected ticks; in contrast, fever, myalgia, lymphadenopathy, and petechial hemorrhaging are common symptoms in humans, with a case fatality ratio of up to 40%. We used the in-solution hybridization capture of 100 annotated immune genes to genotype 121 dromedaries from the UAE tested for seropositivity to CCHFV. Through univariate linear regression analysis, we identified two candidate genes belonging to the innate immune system: FCAR and CLEC2B. These genes have important functions in the host defense against viral infections and in stimulating natural killer cells, respectively. This study opens doors for future research into immune defense mechanisms in an enzootic host against an important zoonotic disease.

The adaptive immune receptors repertoire is highly plastic, with its ability to produce antigen-binding molecules and select those with high affinity for their antigen. Species have developed diverse genetic and structural strategies to create their respective repertoires required for their survival in the different environments. Camelids, until now, considered as a case of evolutionary innovation because of their only heavy-chain antibodies, represent a new mammalian model particularly useful for understanding the role of diversity in the immune system function. Here, we review the structural and functional characteristics and the current status of the genomic organization of camel immunoglobulins (IG) or antibodies, α/ß and γ/δ T cell receptors (TR), and major histocompatibility complex (MHC). In camelid humoral response, in addition to the conventional antibodies, there are IG with “only-heavy-chain” (no light chain, and two identical heavy gamma chains lacking CH1 and with a VH domain designated as VHH). The unique features of these VHH offer advantages in biotechnology and for clinical applications. The TRG and TRD rearranged variable domains of Camelus dromedarius (Arabian camel) display somatic hypermutation (SHM), increasing the intrinsic structural stability in the γ/δ heterodimer and influencing the affinity maturation to a given antigen similar to immunoglobulin genes. The SHM increases the dromedary γ/δ repertoire diversity. In Camelus genus, the general structural organization of the TRB locus is similar to that of the other artiodactyl species, with a pool of TRBV genes positioned at the 5’ end of three in tandem D-J-C clusters, followed by a single TRBV gene with an inverted transcriptional orientation located at the 3’ end. At the difference of TRG and TRD, the diversity of the TRB variable domains is not shaped by SHM and depends from the classical combinatorial and junctional diversity. The MHC locus is located on chromosome 20 in Camelus dromedarius. Cytogenetic and comparative whole genome analyses revealed the order of the three major regions “Centromere-ClassII-ClassIII-ClassI”. Unexpectedly low extent of polymorphisms and haplotypes was observed in all Old World camels despite different geographic origins.

Due to production of special homodimeric heavy chain antibodies, somatic hypermutation of their T-cell receptor genes and unusually low diversity of their major histocompatibility complex genes, camels represent an important model for immunogenetic studies. Here, we analyzed genes encoding selected natural killer cell receptors with a special focus on genes encoding receptors for major histocompatibility complex (MHC) class I ligands in the two domestic camel species, Camelus dromedarius and Camelus bactrianus. Based on the dromedary genome assembly CamDro2, we characterized the genetic contents, organization, and variability of two complex genomic regions, the leukocyte receptor complex and the natural killer complex, along with the natural cytotoxicity receptor genes NCR1, NCR2, and NCR3. The genomic organization of the natural killer complex region of camels differs from cattle, the phylogenetically most closely related species. With its minimal set of KLR genes, it resembles this complex in the domestic pig. Similarly, the leukocyte receptor complex of camels is strikingly different from its cattle counterpart. With KIR pseudogenes and few LILR genes, it seems to be simpler than in the pig. The syntenies and protein sequences of the NCR1, NCR2, and NCR3 genes in the dromedary suggest that they could be human orthologues. However, only NCR1 and NCR2 have a structure of functional genes, while NCR3 appears to be a pseudogene. High sequence similarities between the two camel species as well as with the alpaca Vicugna pacos were observed. The polymorphism in all genes analyzed seems to be generally low, similar to the rest of the camel genomes. This first report on natural killer cell receptor genes in camelids adds new data to our understanding of specificities of the camel immune system and its functions, extends our genetic knowledge of the innate immune variation in dromedaries and Bactrian camels, and contributes to studies of natural killer cell receptors evolution in mammals.

Polymorphic markers on the male-specific part of the Y chromosome (MSY) provide useful information for tracking male genealogies. While maternal lineages are well studied in Old World camelids using mitochondrial DNA, the lack of a Y-chromosomal reference sequence hampers the analysis of male-driven demographics. Recently, a shotgun assembly of the horse MSY was generated based on short read next generation sequencing data. The haplotype network resulting from single copy MSY variants using the assembly as a reference revealed sufficient resolution to trace individual male lines in this species. In a similar approach we generated a 3.8 Mbp sized assembly of the MSY of Camelus bactrianus. The camel MSY assembly was used as a reference for variant calling using short read data from eight Old World camelid individuals. Based on 596 single nucleotide variants we revealed a Y-phylogenetic network with seven haplotypes. Wild and domestic Bactrian camels were clearly separated into two different haplogroups with an estimated divergence time of 26,999 ± 2,268 years. Unexpectedly, one wild camel clustered into the domestic Bactrian camels’ haplogroup. The observation of a domestic paternal lineage within the wild camel population is concerning in view of the importance to conserve the genetic integrity of these highly endangered species in their natural habitat.

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.

Myostatin, a negative regulator of skeletal muscle mass in animals, has been shown to play a role in determining muscular hypertrophy in several livestock species, and a high degree of polymorphism has been previously reported for this gene in humans and cattle. In this study, we provide a characterization of the myostatin gene in the dromedary (Camelus dromedarius) at the genomic, transcript and protein level. The gene was found to share high structural and sequence similarity with other mammals, notably Old World camelids. 3D modeling highlighted several non-conservative SNP variants compared to the bovine, as well as putative functional variants involved in the stability of the myostatin dimer. NGS data for nine dromedaries from various countries revealed 66 novel SNPs, all of them falling either upstream or downstream the coding region. The analysis also confirmed the presence of three previously described SNPs in intron 1, predicted here to alter both splicing and transcription factor binding sites (TFBS), thus possibly impacting myostatin processing and/or regulation. Several putative TFBS were identified in the myostatin upstream region, some of them belonging to the myogenic regulatory factor family. Patterns of SNP distribution across countries, as suggested by Bayesian clustering of the nine dromedaries using the 69 SNPs, pointed to weak geographic differentiation, in line with known recurrent gene flow at ancient trading centers along caravan routes. Myostatin expression was investigated in a set of 8 skeletal muscles, both at transcript and protein level, via Digital Droplet PCR and Western Blotting, respectively. No significant differences were observed at the transcript level, while, at the protein level, the only significant differences concerned the promyostatin dimer (75 kDa), in four pair-wise comparisons, all involving the tensor fasciae latae muscle. Beside the mentioned band at 75 kDa, additional bands were observed at around 40 and 25 kDa, corresponding to the promyostatin monomer and the active C-terminal myostatin dimer, respectively. Their weaker intensity suggests that the unprocessed myostatin dimers could act as important reservoirs of slowly available myostatin forms. Under this assumption, the sequential cleavage steps may contribute additional layers of control within an already complex regulatory framework.