Transcriptomic profiling of different developmental stages reveals parasitic strategies of Wohlfahrtia magnifica, a myiasis-causing flesh fly.

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.

Whole-genome SNP genotyping unveils ancestral and recent introgression in wild and domestic goats.

After the domestication of goats around 10,000 years before the present (BP), humans transported goats far beyond the range of their wild ancestor, the bezoar goat. This brought domestic goats into contact with many wild goat species such as ibex and markhor, enabling introgression between domestic and wild goats. To investigate this, while shedding light on the taxonomic status of wild and domestic goats, we analysed genome-wide SNP data of 613 specimens from 14 taxonomic units, including Capra hircus, C. pyrenaica, C. ibex (from Switzerland, Austria, Germany and Slovenia), C. aegagrus aegagrus, C. a. cretica, C. h. dorcas, C. caucasica caucasica, C. c. severtzovi, C. c. cylindricornis, C. falconeri, C. sibirica sibirica, C. s. alaiana and C. nubiana, as well as Oreamnos americanus (mountain goat) as an outgroup. To trace gene flow between domestic and wild goats, we integrated genotype data of local goat breeds from the Alps as well as from countries such as Spain, Greece, Türkiye, Egypt, Sudan, Iran, Russia (Caucasus and Altai) and Pakistan. Our phylogenetic analyses displayed a clear separation between bezoar-type and ibex-type clades with wild goats from the Greek islands of Crete and Youra clustered within domestic goats, confirming their feral origin. Our analyses also revealed gene flow between the lineages of Caucasian tur and domestic goats that most likely occurred before or during early domestication. Within the clade of domestic goats, analyses inferred gene flow between African and Iberian goats. The detected events of introgression were consistent with previous reports and offered interesting insights into the historical relationships among domestic and wild goats.

Participatory investigation of goat farmers' breeding practices, trait preference, and selection criteria in Burkina Faso.

The community-based breeding program (CBBP) is an innovative approach recommended for genetic improvement and sustainable use of animal genetic resources in extensive farming systems. Successful implementation of this approach requires an understanding of the characteristics of production systems, breeding objectives, and farmers' trait preference. This study aimed to identify the selection criteria of goat farmers in rural areas of Burkina Faso and their potential implications in establishing CBBP. Following focus group discussions, a well-structured questionnaire was designed and administered to 372 randomly selected goat farmers in two different agro-ecological zones. A list of traits obtained during focus group discussions was provided to farmers individually, and they were asked to rank the ones they preferentially use to select breeding animals. Statistical tests were conducted to compare data between the two agro-ecological zones. The results showed that the average goat flock per household was higher (P < 0.05) in the Sudanian (15.68 ± 13.76), compared to the Sudano-Sahelian area (12.93 ± 13.3). Adult females were the dominant age-sex group in both areas. Reasons for culling, keeping breeding bucks, and castration practice were significantly different (P < 0.05) among agro-ecological zones. The most important common criterion for selection in the two zones was body size, coat color, and growth rate for the bucks and does, while fertility (0.06) parameters including twining ability (0.18), kidding frequency (0.11), and mothering ability (0.15) were furthermore considered for breeding does selection. These findings provide valuable insights for developing CBBPs tailored to goat production in the study areas.

A chromosome-scale high-contiguity genome assembly of the cheetah (Acinonyx jubatus).

The cheetah (Acinonyx jubatus, SCHREBER 1775) is a large felid and is considered the fastest land animal. Historically, it inhabited open grassland across Africa, the Arabian Peninsula, and southwestern Asia; however, only small and fragmented populations remain today. Here, we present a de novo genome assembly of the cheetah based on PacBio continuous long reads and Hi-C proximity ligation data. The final assembly (VMU_Ajub_asm_v1.0) has a total length of 2.38 Gb, of which 99.7% are anchored into the expected 19 chromosome-scale scaffolds. The contig and scaffold N50 values of 96.8 Mb and 144.4 Mb, respectively, a BUSCO completeness of 95.4% and a k-mer completeness of 98.4%, emphasize the high quality of the assembly. Furthermore, annotation of the assembly identified 23,622 genes and a repeat content of 40.4%. This new highly contiguous and chromosome-scale assembly will greatly benefit conservation and evolutionary genomic analyses and will be a valuable resource, e.g., to gain a detailed understanding of the function and diversity of immune response genes in felids.

The potential and shortcomings of mitochondrial DNA analysis for cheetah conservation management.

There are only about 7,100 adolescent and adult cheetahs (Acinonyx jubatus) remaining in the wild. With the majority occurring outside protected areas, their numbers are rapidly declining. Evidence-based conservation measures are essential for the survival of this species. Genetic data is routinely used to inform conservation strategies, e.g., by establishing conservation units (CU). A commonly used marker in conservation genetics is mitochondrial DNA (mtDNA). Here, we investigated the cheetah's phylogeography using a large-scale mtDNA data set to refine subspecies distributions and better assign individuals to CUs. Our dataset mostly consisted of historic samples to cover the cheetah's whole range as the species has been extinct in most of its former distribution. While our genetic data largely agree with geography-based subspecies assignments, several geographic regions show conflicting mtDNA signals. Our analyses support previous findings that evolutionary forces such as incomplete lineage sorting or mitochondrial capture likely confound the mitochondrial phylogeography of this species, especially in East and, to some extent, in Northeast Africa. We caution that subspecies assignments solely based on mtDNA should be treated carefully and argue for an additional standardized nuclear single nucleotide polymorphism (SNP) marker set for subspecies identification and monitoring. However, the detection of the A. j. soemmeringii specific haplogroup by a newly designed Amplification-Refractory Mutation System (ARMS) can already provide support for conservation measures. Supplementary Information: The online version contains supplementary material available at 10.1007/s10592-022-01483-1.

Genomic analyses show extremely perilous conservation status of African and Asiatic cheetahs (Acinonyx jubatus).

We live in a world characterized by biodiversity loss and global environmental change. The extinction of large carnivores can have ramifying effects on ecosystems like an uncontrolled increase in wild herbivores, which in turn can have knock-on impacts on vegetation regeneration and communities. Cheetahs (Acinonyx jubatus) serve important ecosystem functions as apex predators; yet, they are quickly heading towards an uncertain future. Threatened by habitat loss, human-wildlife conflict and illegal trafficking, there are only approximately 7100 individuals remaining in nature. We present the most comprehensive genome-wide analysis of cheetah phylogeography and conservation genomics to date, assembling samples from nearly the entire current and past species' range. We show that their phylogeography is more complex than previously thought, and that East African cheetahs (A. j. raineyi) are genetically distinct from Southern African individuals (A. j. jubatus), warranting their recognition as a distinct subspecies. We found strong genetic differentiation between all classically recognized subspecies, thus refuting earlier findings that cheetahs show only little differentiation. The strongest differentiation was observed between the Asiatic and all the African subspecies. We detected high inbreeding in the Critically Endangered Iranian (A. j. venaticus) and North-western (A. j. hecki) subspecies, and show that overall cheetahs, along with snow leopards, have the lowest genome-wide heterozygosity of all the big cats. This further emphasizes the cheetah's perilous conservation status. Our results provide novel and important information on cheetah phylogeography that can support evidence-based conservation policy decisions to help protect this species. This is especially relevant in light of ongoing and proposed translocations across subspecies boundaries, and the increasing threats of illegal trafficking.

Unlocking the potential of a validated single nucleotide polymorphism array for genomic monitoring of trade in cheetahs (Acinonyx jubatus).

Cheetahs (Acinonyx jubatus) are listed as vulnerable on the International Union for Conservation of Nature Red List of Threatened Species. Threats include loss of habitat, human-wildlife conflict and illegal wildlife trade. In South Africa, the export of wild cheetah is a restricted activity under the Convention on International Trade in Endangered Species of Wild Fauna and Flora (CITES), however, limited legal trade is permitted of animals born to captive parents. To effectively monitor the legal and illegal trade in South Africa, it was thus essential to develop a validated molecular test. Here, we designed a single nucleotide polymorphism (SNP) array for cheetah from Double Digest Restriction Associated DNA sequencing data for individual identification and parentage testing. In order to validate the array, unrelated individuals and 16 family groups consisting of both parents and one to three offspring were genotyped using the Applied Biosystems™ QuantStudio™ 12K Flex Real-Time PCR System. In addition, parentage assignments were compared to microsatellite data. Cross-species amplification was tested in various felids and cheetah sub-species in order to determine the utility of the SNP array in other species. We obtained successful genotyping results for 218 SNPs in cheetah (A. j. jubatus) with an optimal DNA input concentration ranging from 10 to 30 ng/µl. The combination of SNPs had a higher resolving power for individual identification compared to microsatellites and provided high assignment accuracy in known pedigrees. Cross-species amplification in other felids was determined to be limited. However, the SNP array demonstrated a clear genetic discrimination of two cheetah subspecies tested here. We conclude that the described SNP array is suitable for accurate parentage assignment and provides an important traceability tool for forensic investigations of cheetah trade.

Population structure, inbreeding and admixture in local cattle populations managed by community-based breeding programs in Burkina Faso.

High-throughput genomic markers provide an opportunity to assess important indicators of genetic diversity for populations managed in livestock breeding programs. While well-structured breeding programs are common in developed countries, in developing country situations, especially in West Africa, on-farm performance and pedigree recordings are rare, and thus, genomic markers provide insights to the levels of genetic diversity, inbreeding and introgression by other breeds. In this study, we analysed key population parameters such as population structure, admixture and levels of inbreeding in three neighbouring populations of African taurine and taurine × Zebu crosses managed by community-based breeding programs in the South-West of Burkina Faso. The three populations were pure Baoulé (called Lobi locally) in sedentary production systems, Baoulé x Zebu crossbreds in sedentary systems and Zebu × Baoulé crossbreds in transhumant production systems, respectively. The total sample analysed included 631 animals and 38,207 single nucleotide polymorphisms after quality control. Results of principal component and admixture analyses confirmed the genetic background of two distinct ancestral populations (taurine and zebuine) and levels of admixture in all three breeding populations, including the presumably pure Baoulé group of animals. Inbreeding levels were moderate, compared to European dairy and beef cattle populations and higher than those of Brazilian Nellore cattle. Very few animals with inbreeding levels indicating parent-offspring or full sib mating were observed, and inbreeding levels indicating half sib mating were also rare. For the management of breeding populations, farmers were advised to exchange best young bulls. The crossbreeding levels of presumably pure Baoulé animals are of concern to the breeding program due to the high level of endangerment of pure African taurine cattle populations across West Africa. Future rounds of bull selection in the community-based breeding program will make use of genomic information about admixture levels.

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

BACKGROUND: Immune-response (IR) genes have an important role in the defense against highly variable pathogens, and therefore, diversity in these genomic regions is essential for species' survival and adaptation. Although current genome assemblies from Old World camelids are very useful for investigating genome-wide diversity, demography and population structure, they have inconsistencies and gaps that limit analyses at local genomic scales. Improved and more accurate genome assemblies and annotations are needed to study complex genomic regions like adaptive and innate IR genes. RESULTS: In this work, we improved the genome assemblies of the three Old World camel species - domestic dromedary and Bactrian camel, and the two-humped wild camel - via different computational methods. The newly annotated dromedary genome assembly CamDro3 served as reference to scaffold the NCBI RefSeq genomes of domestic Bactrian and wild camels. These upgraded assemblies were then used to assess nucleotide diversity of IR genes within and between species, and to compare the diversity found in immune genes and the rest of the genes in the genome. We detected differences in the nucleotide diversity among the three Old World camelid species and between IR gene groups, i.e., innate versus adaptive. Among the three species, domestic Bactrian camels showed the highest mean nucleotide diversity. Among the functionally different IR gene groups, the highest mean nucleotide diversity was observed in the major histocompatibility complex. CONCLUSIONS: The new camel genome assemblies were greatly improved in terms of contiguity and increased size with fewer scaffolds, which is of general value for the scientific community. This allowed us to perform in-depth studies on genetic diversity in immunity-related regions of the genome. Our results suggest that differences of diversity across classes of genes appear compatible with a combined role of population history and differential exposures to pathogens, and consequent different selective pressures.

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

Dromedaries have been fundamental to the development of human societies in arid landscapes and for long-distance trade across hostile hot terrains for 3,000 y. Today they continue to be an important livestock resource in marginal agro-ecological zones. However, the history of dromedary domestication and the influence of ancient trading networks on their genetic structure have remained elusive. We combined ancient DNA sequences of wild and early-domesticated dromedary samples from arid regions with nuclear microsatellite and mitochondrial genotype information from 1,083 extant animals collected across the species' range. We observe little phylogeographic signal in the modern population, indicative of extensive gene flow and virtually affecting all regions except East Africa, where dromedary populations have remained relatively isolated. In agreement with archaeological findings, we identify wild dromedaries from the southeast Arabian Peninsula among the founders of the domestic dromedary gene pool. Approximate Bayesian computations further support the "restocking from the wild" hypothesis, with an initial domestication followed by introgression from individuals from wild, now-extinct populations. Compared with other livestock, which show a long history of gene flow with their wild ancestors, we find a high initial diversity relative to the native distribution of the wild ancestor on the Arabian Peninsula and to the brief coexistence of early-domesticated and wild individuals. This study also demonstrates the potential to retrieve ancient DNA sequences from osseous remains excavated in hot and dry desert environments.

Range-Wide Snow Leopard Phylogeography Supports Three Subspecies.

The snow leopard, Panthera uncia, is an elusive high-altitude specialist that inhabits vast, inaccessible habitat across Asia. We conducted the first range-wide genetic assessment of snow leopards based on noninvasive scat surveys. Thirty-three microsatellites were genotyped and a total of 683 bp of mitochondrial DNA sequenced in 70 individuals. Snow leopards exhibited low genetic diversity at microsatellites (AN = 5.8, HO = 0.433, HE = 0.568), virtually no mtDNA variation, and underwent a bottleneck in the Holocene (∼8000 years ago) coinciding with increased temperatures, precipitation, and upward treeline shift in the Tibetan Plateau. Multiple analyses supported 3 primary genetic clusters: (1) Northern (the Altai region), (2) Central (core Himalaya and Tibetan Plateau), and (3) Western (Tian Shan, Pamir, trans-Himalaya regions). Accordingly, we recognize 3 subspecies, Panthera uncia irbis (Northern group), Panthera uncia uncia (Western group), and Panthera uncia uncioides (Central group) based upon genetic distinctness, low levels of admixture, unambiguous population assignment, and geographic separation. The patterns of variation were consistent with desert-basin "barrier effects" of the Gobi isolating the northern subspecies (Mongolia), and the trans-Himalaya dividing the central (Qinghai, Tibet, Bhutan, and Nepal) and western subspecies (India, Pakistan, Tajikistan, and Kyrgyzstan). Hierarchical Bayesian clustering analysis revealed additional subdivision into a minimum of 6 proposed management units: western Mongolia, southern Mongolia, Tian Shan, Pamir-Himalaya, Tibet-Himalaya, and Qinghai, with spatial autocorrelation suggesting potential connectivity by dispersing individuals up to ∼400 km. We provide a foundation for global conservation of snow leopard subspecies, and set the stage for in-depth landscape genetics and genomic studies.

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

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

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

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

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

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

Integrating One Health into Systems Science.

One Health Systems Science. The three subsystems of One Health (ecosystem, human and animal health) are integrated in the Systems Science concept, where objects or adaptive agents (circles) interact with a dynamic environment, and Systems Thinking can lead it intervations (Systems Design) generating a change in One Health outcomes. Real-time genomic data retrieved from the three subsystems porvide information fo Systems Thinking and Systems Design.Unlabelled Image.

Diversity of major histocompatibility complex (MHC) and natural killer cell receptor (NKR) genes and their interactions in domestic horses.

The immunogenome is the part of the genome that underlies immune mechanisms and evolves under various selective pressures. Two complex regions of the immunogenome, major histocompatibility complex (MHC) and natural killer cell receptor (NKR) genes, play an important role in the response to selective pressures of pathogens. Their importance is expressed by their genetic polymorphism at the molecular level, and their diversity associated with different types of diseases at the population level. Findings of associations between specific combinations of MHC/NKR haplotypes with different diseases in model species suggest that these gene complexes did not evolve independently. No such associations have been described in horses so far. The aim of the study was to detect associations between MHC and NKR gene/microsatellite haplotypes in three horse breed groups (Camargue, African, and Romanian) by statistical methods; chi-square test, Fisher's exact test, Pearson's goodness-of-fit test and logistic regression. Associations were detected for both MHC/NKR genes and microsatellites; the most significant associations were found between the most variable KLRA3 gene and the EQCA-1 or EQCA-2 genes. This finding supports the assumption that the KLRA3 is an important receptor for MHC I and that interactions of these molecules play important roles in the horse immunity and reproduction. Despite some limitations of the study such as low numbers of horses or lack of knowledge of the selected genes functions, the results were consistent across different statistical methods and remained significant even after overconservative Bonferroni corrections. We therefore consider them biologically plausible.

Diversity of selected toll-like receptor genes in cheetahs (Acinonyx jubatus) and African leopards (Panthera pardus pardus).

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.

Domestication of the Dromedary Revisited and Its Consequences for Legislation as to Keeping Livestock or Pet Animals.

Being in an advanced stage of domestication is a newly proposed requirement to decide which animals can be safely kept by humans. Dutch legislators were the first to apply it and other European countries may be tempted to adopt a similar approach. Unexpectedly, the Dutch assessors considered the dromedary (Camelus dromedarius) as being insufficiently domesticated and this species will therefore no longer be able to be kept as a production animal from 2024 onwards. In a recent publication on this topic, we showed that the domestication of the dromedary is actually very advanced. In this paper, we apply the same criteria that were used by the Dutch assessors to determine the degree of domestication, taking into account the most recent scientific developments in this area, even though it should be noted that these criteria have neither been peer-reviewed, nor published in an international scientific journal. For the sake of comparison, and in order to validate the procedure, we also applied these criteria to the house cat. The results confirm that the dromedary is highly domesticated, but also that the house cat (Felis silvestris catus) is at most semi-domesticated. Obviously, we agree with the decision of the Dutch legislators to place the house cat on the positive list, but our analysis demonstrates that this was decided on false grounds. Our analysis makes it clear that the requirement of being in an advanced stage of domestication is not suitable. Instead of maintaining this requirement, we recommend implementing evidence-based, peer-reviewed methods to decide which animals can be kept by humans, and to include species specific-guidelines in the legislation on how this can be achieved safely.

Comparative genomics of the Leukocyte Receptor Complex in carnivores.

Background: The 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. Methods: As 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. Results: Seven 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. Conclusions: The 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.