Genetic diversity, evolution and selection in the major histocompatibility complex DRB and DQB loci in the family Equidae.

BACKGROUND: The mammalian Major Histocompatibility Complex (MHC) is a genetic region containing highly polymorphic genes with immunological functions. MHC class I and class II genes encode antigen-presenting molecules expressed on the cell surface. The MHC class II sub-region contains genes expressed in antigen presenting cells. The antigen binding site is encoded by the second exon of genes encoding antigen presenting molecules. The exon 2 sequences of these MHC genes have evolved under the selective pressure of pathogens. Interspecific differences can be observed in the class II sub-region. The family Equidae includes a variety of domesticated, and free-ranging species inhabiting a range of habitats exposed to different pathogens and represents a model for studying this important part of the immunogenome. While equine MHC class II DRA and DQA loci have received attention, the genetic diversity and effects of selection on DRB and DQB loci have been largely overlooked. This study aimed to provide the first in-depth analysis of the MHC class II DRB and DQB loci in the Equidae family. RESULTS: Three DRB and two DQB genes were identified in the genomes of all equids. The genes DRB2, DRB3 and DQB3 showed high sequence conservation, while polymorphisms were more frequent at DRB1 and DQB1 across all species analyzed. DQB2 was not found in the genome of the Asiatic asses Equus hemionus kulan and E. h. onager. The bioinformatic analysis of non-zero-coverage-bases of DRB and DQB genes in 14 equine individual genomes revealed differences among individual genes. Evidence for recombination was found for DRB1, DRB2, DQB1 and DQB2 genes. Trans-species allele sharing was identified in all genes except DRB1. Site-specific selection analysis predicted genes evolving under positive selection both at DRB and DQB loci. No selected amino acid sites were identified in DQB3. CONCLUSIONS: The organization of the MHC class II sub-region of equids is similar across all species of the family. Genomic sequences, along with phylogenetic trees suggesting effects of selection as well as trans-species polymorphism support the contention that pathogen-driven positive selection has shaped the MHC class II DRB/DQB sub-regions in the Equidae.

Comparative genomics of the Natural Killer Complex in carnivores.

Background: The 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. Methods: This 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. Results: The 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. Conclusions: Knowledge 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.

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.

Major histocompatibility complex and other allergy-related candidate genes associated with insect bite hypersensitivity in Icelandic horses.

Insect bite hypersensitivity (IBH) is an allergic dermatitis of horses caused by bites of insects. IBH is a multifactorial disease with contribution of genetic and environmental factors. Candidate gene association analysis of IBH was performed in a group of 89 Icelandic horses all born in Iceland and imported to Europe. Horses were classified in IBH-affected and non-affected based on clinical signs and history of recurrent dermatitis, and on the results of an in vitro sulfidoleukotriene (sLT)-release assay with Culicoides nubeculosus and Simulium vittatum extract. Different genetic markers were tested for association with IBH by the Fisher's exact test. The effect of the major histocompatibility complex (MHC) gene region was studied by genotyping five microsatellites spanning the MHC region (COR112, COR113, COR114, UM011 and UMN-JH34-2), and exon 2 polymorphisms of the class II Eqca-DRA gene. Associations with Eqca-DRA and COR113 were identified (p < 0.05). In addition, a panel of 20 single nucleotide polymorphisms (SNPs) in 17 candidate allergy-related genes was tested. During the initial screen, no marker from the panel was significantly (p < 0.05) associated with IBH. Five SNPs associated with IBH at p < 0.10 were therefore used for analysis of combined genotypes. Out of them, SNPs located in the genes coding for the CD14 receptor (CD14), interleukin 23 receptor (IL23R), thymic stromal lymphopoietin (TSLP) and transforming growth factor beta 3 (TGFB3) molecules were associated with IBH as parts of complex genotypes. These results are supported by similar associations and by expression data from different horse populations and from human studies.

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.

Newly identified variability of the antigen binding site coding sequences of the equine major histocompatibility complex class I and class II genes.

The major histocompatibility complex (MHC) with its class I and II genes plays a crucial role in the immune response to pathogens by presenting oligopeptide antigens to various immune response effector cells. In order to counteract the vast variability of infectious agents, MHC class I and II genes usually retain high levels of SNPs mainly concentrated in the exons encoding the antigen binding sites. The aim of the study was to reveal new variability of selected MHC genes with a special focus on MHC class I physical haplotypes. Long-range NGS to was used to identify exon 2-exon 3 alleles in three genetically distinct horse breeds. A total of 116 allelic variants were found in the MHC class I genes Eqca-1, Eqca-2, Eqca-7 and Eqca-Ψ, 112 of which were novel. The MHC class II DRA locus was confirmed to comprise five exon 2 alleles, and no new sequences were observed. Additional variability in terms of 15 novel exon 2 alleles was identified in the DQA1 locus. Extensive overall variability across the entire MHC region was confirmed by an analysis of MHC-linked microsatellite loci. Both diversifying and purifying selection were detected within the MHC class I and II loci analyzed.

The deadly face of felid killer cells: The cytotoxic proteins and their genes.

Natural killer cells and cytotoxic T lymphocytes are the main cell populations of the immune system able to directly kill target cells via cytotoxic granules. Different mammalian species may differ in specific features of their pore-forming protein (perforin) and granule-bound serine proteases (granzymes). One perforin gene (PRF1) and four genes encoding granzymes A, B, H, and K (GZMA, GZMB, GZMH, GZMK) were identified in the reference genomes of felids. The objective of this work was to characterize the genes PRF1, GZMA and GZMB in a panel of 17 felid species by next-generation re-sequencing. A search of available felid genomes (17 species) retrieved the coding sequences of these genes for comparison to our data. Both sets of sequences or their combinations (23 species) were used for phylogenetic and selection analyses. Nucleotide PRF1, GZMA and GZMB sequences showed high similarities between felid species (over 95% identity). All trees derived from coding sequences expressed phylogenetic relationships corresponding to the zoological taxonomy of the Felidae, except GZMA. No effects of positive selection were detected in the genes studied, however, effects of purifying selection were observed for PRF1 and GZMA. The conservation of PRF1 is in agreement with its critical biological function. The differentiation observed between granzyme sub-families may reflect an adaptation to pathogen variation. The need to maintain important gene functions and at the same time cope with various pathogens may lead to an equilibrium between positive and negative selective pressures acting on GZMB. The within-species variability in wild felid populations merits further investigation.

Comparative Genomics of the Major Histocompatibility Complex (MHC) of Felids.

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 natural cytotoxicity receptor genes in the family Felidae.

Natural killer (NK) cells belong to the innate immune system. The germline-encoded natural killer cell receptors represent activating and inhibitory receptors regulating multiple NK cell activities. The natural cytotoxicity receptors (NCRs) are activating natural cytotoxicity triggering receptors 1, 2, and 3 (NKp46, NKp44, and NKp30), encoded by the genes NCR1, NCR2, and NCR3, respectively. NCRs may be expressed in different cell types engaged in mechanisms of innate and adaptive immunity. The family Felidae, comprising the domestic cat and a wide variety of free-ranging species represents a well-suited model for biomedical and evolutionary studies. We characterized the NCR1, NCR2, and NCR3 genes in a panel of felid species. We confirmed the presence of potentially functional genes NCR1, NCR2, and NCR3 in all species. All three genes are conserved within the family and are similar to other phylogenetically related mammalian families. The NCR1 and NCR2 phylogenetic trees based on both nucleotide and protein sequences corresponded to the current zoological taxonomy, with some exceptions suggesting effects of different selection pressures in some species. Highly conserved NCR3 sequences did not allow a robust phylogenetic analysis. Most interspecific differences both at the nucleotide and protein level were found in NCR2. Within species, the most polymorphic CDS was detected in NCR1. Selection analyses indicated the effects of purifying selection on individual amino acid sites in all three genes. In stray cats, a rather high intraspecific diversity was observed.

[Hypophosphatasia--biochemical and clinical manifestations, molecular genetic principles]. / Hypofosfatázia--biochemické a klinické prejavy, molekulovo-genetická podstata.

Hypophosphatasia is a rare hereditary metabolic disorder accompanying deficit of tissue nonspecific serum alkaline phosphatase. The incidence of overt forms is estimated about 1:100000 live births. In the prenatal manifestation the disease may cause severe damage to the foetus with intrauterine death. In children there is a defect of mineralization with rickets signs and the subsequent hypercalcaemia a hypercalciuria may lead to death. In adults the main manifestation is osteomalacia, skeletal deformities and fractures, early arthritis. In severe forms the heredity is autosomal recessive type. In mild forms the heredity may be dominant or recessive. In two case reports we present clinical course of the disease in two adult sisters, where diagnosis of hypophosphatasia was first time confirmed in Slovak population using molecular genetic methods.

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

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

The Population Diversity of Candidate Genes for Resistance/Susceptibility to Coronavirus Infection in Domestic Cats: An Inter-Breed Comparison.

Feline coronavirus (FCoV) is a complex pathogen causing feline infectious peritonitis (FIP). Host genetics represents a factor contributing to the pathogenesis of the disease. Differential susceptibility of various breeds to FIP was reported with controversial results. The objective of this study was to compare the genetic diversity of different breeds on a panel of candidate genes potentially affecting FCoV infection. One hundred thirteen cats of six breeds were genotyped on a panel of sixteen candidate genes. SNP allelic/haplotype frequencies were calculated; pairwise FST and molecular variance analyses were performed. Principal coordinate (PCoA) and STRUCTURE analyses were used to infer population structure. Interbreed differences in allele frequencies were observed. PCoA analysis performed for all genes of the panel indicated no population substructure. In contrast to the full marker set, PCoA of SNP markers associated with FCoV shedding (NCR1 and SLX4IP) showed three clusters containing only alleles associated with susceptibility to FCoV shedding, homozygotes and heterozygotes for the susceptibility alleles, and all three genotypes, respectively. Each cluster contained cats of multiple breeds. Three clusters of haplotypes were identified by PCoA, two clusters by STRUCTURE. Haplotypes of a single gene (SNX5) differed significantly between the PCoA clusters.

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

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

Crimean-Congo Hemorrhagic Fever Virus Past Infections Are Associated with Two Innate Immune Response Candidate Genes in Dromedaries.

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.

Complex variation in the KLRA (LY49) immunity-related genomic region in horses.

Natural killer (NK) cells play important roles in innate and adaptive immunity, as well as in the reproduction of placental mammals. Ly49 (KLRA) molecules represent a lectin-like type of NK cell receptor encoded within a complex genomic region, the NK cell complex. In rodents and horses, an expansion of the genes encoding Ly49 receptors leading to the formation of a gene family was observed. High sequence similarities and frequent high polymorphism of multiple family members represent an obstacle both for their individual identification and for annotation in the reference genomes of their respective species. Here, we focused on resolving complex variation of the KLRA gene family observed in domestic and Przewalski's horses. The KLRA (LY49) genomic region contains six genes (KLRA2-KLRA7) and one putative pseudogene, KLRA1. Two types of polymorphism were observed in the horses analyzed. Copy number variation between haplotypes was documented for the gene KLRA7 by polymerase chain reaction. As expected, the major source of variation of all KLRA genes, including KLRA7, is because of single nucleotide polymorphisms, many of them being nonsynonymous substitutions. Extensive allelic variability of the expanded KLRA (LY49) genes was observed. For four out of the six functional KLRA, high numbers of novel allelic amino acid sequence variants were identified in the genes studied, suggesting that this variation might be of functional importance, especially in the context of high polymorphism of their presumed ligands encoded by major histocompatibility complex class I genes. In fact, polymorphic amino acid sites were mostly found in the ligand-binding C-type lectin-like domain of the putative receptor molecule.

Candidate Gene Markers Associated with Fecal Shedding of the Feline Enteric Coronavirus (FECV).

The Feline coronavirus (FCoV) can cause a fatal disease, the Feline Infectious Peritonitis. Persistent shedders represent the most important source of infection. The role of the host in FCoV fecal shedding is unknown. The objective of this study was to develop gene markers and to test their associations with FCoV shedding patterns. Fecal samples were taken from 57 cats of 12 breeds on the day 0 and after 2, 4 and 12 months. Variation from persistent and/or high-intensity shedding to no shedding was observed. Thirteen immunity-related genes were selected as functional and positional/functional candidates. Positional candidates were selected in a candidate region detected by a GWAS analysis. Tens to hundreds of single nucleotide polymorphisms (SNPs) per gene were identified using next generation sequencing. Associations with different phenotypes were assessed by chi-square and Fisher's exact tests. SNPs of one functional and one positional candidate (NCR1 and SLX4IP, respectively) and haplotypes of four genes (SNX5, NCR2, SLX4IP, NCR1) were associated with FCoV shedding at pcorected < 0.01. Highly significant associations were observed for extreme phenotypes (persistent/high-intensity shedders and non-shedders) suggesting that there are two major phenotypes associated with different genotypes, highly susceptible cats permanently shedding high amounts of viral particles and resistant non-shedders.

Spectrum of DHCR7 mutations in Slovak patients with Smith-Lemli-Opitz syndrome and detection of common mutations by PCR-based assays.

The Smith-Lemli-Opitz syndrome (SLOS), an autosomal recessive disorder associated with multiple developmental malformations, is caused by a large spectrum of mutations in the DHCR7 gene. Mutations in the DHCR7 gene lead to a 7-dehydrocholesterol reductase deficiency, which is the final enzyme in the pathway of the cholesterol biosynthesis. Reduced cholesterol levels and elevated concentrations of its precursor 7-dehydrocholesterol in plasma and tissues are the major biochemical hallmarks of this disorder. In all patients a biochemical analysis of blood sterols using the gas chromatography/mass spectrometry was performed to confirm the clinical diagnosis of SLOS. We have also determined the mutational spectrum of DHCR7 gene in 17 Slovak patients. We identified six different mutations: nonsense mutation W151X and missense mutations V326L, L109P, G410S, R352Q, Y432C. Mutations W151X and V326L accounted for 76% of the SLOS alleles in Slovak population. The Slovak mutational spectrum is similar to that observed in other Central European countries. We also report simple polymerase chain reaction (PCR)-based assays that allow efficient and rapid mutation analysis.

Natural Killer Cell Receptor Genes in Camels: Another Mammalian Model.

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.

Natural killer cell receptor genes in the family Equidae: not only Ly49.

Natural killer (NK) cells have important functions in immunity. NK recognition in mammals can be mediated through killer cell immunoglobulin-like receptors (KIR) and/or killer cell lectin-like Ly49 receptors. Genes encoding highly variable NK cell receptors (NKR) represent rapidly evolving genomic regions. No single conservative model of NKR genes was observed in mammals. Single-copy low polymorphic NKR genes present in one mammalian species may expand into highly polymorphic multigene families in other species. In contrast to other non-rodent mammals, multiple Ly49-like genes appear to exist in the horse, while no functional KIR genes were observed in this species. In this study, Ly49 and KIR were sought and their evolution was characterized in the entire family Equidae. Genomic sequences retrieved showed the presence of at least five highly conserved polymorphic Ly49 genes in horses, asses and zebras. These findings confirmed that the expansion of Ly49 occurred in the entire family. Several KIR-like sequences were also identified in the genome of Equids. Besides a previously identified non-functional KIR-Immunoglobulin-like transcript fusion gene (KIR-ILTA) and two putative pseudogenes, a KIR3DL-like sequence was analyzed. In contrast to previous observations made in the horse, the KIR3DL sequence, genomic organization and mRNA expression suggest that all Equids might produce a functional KIR receptor protein molecule with a single non-mutated immune tyrosine-based inhibition motif (ITIM) domain. No evidence for positive selection in the KIR3DL gene was found. Phylogenetic analysis including rhinoceros and tapir genomic DNA and deduced amino acid KIR-related sequences showed differences between families and even between species within the order Perissodactyla. The results suggest that the order Perissodactyla and its family Equidae with expanded Ly49 genes and with a potentially functional KIR gene may represent an interesting model for evolutionary biology of NKR genes.

Genetic diversity and conservation in a small endangered horse population.

The Old Kladruber horses arose in the 17th century as a breed used for ceremonial purposes. Currently, grey and black coat colour varieties exist as two sub-populations with different recent breeding history. As the population underwent historical bottlenecks and intensive inbreeding, loss of genetic variation is considered as the major threat. Therefore, genetic diversity in neutral and non-neutral molecular markers was examined in the current nucleus population. Fifty microsatellites, 13 single nucleotide polymorphisms (SNPs) in immunity-related genes, three mutations in coat colour genes and one major histocompatibility (MHC-DRA) gene were studied for assessing genetic diversity after 15 years of conservation. The results were compared to values obtained in a similar study 13 years ago. The extent of genetic diversity of the current population was comparable to other breeds, despite its small size and isolation. The comparison between 1997 and 2010 did not show differences in the extent of genetic diversity and no loss of allele richness and/or heterozygosity was observed. Genetic differences identified between the black and grey sub-populations observed 13 years ago persisted. Deviations from the Hardy-Weinberg equilibrium found in 19 microsatellite loci and in five SNP loci are probably due to selective breeding. No differences between neutral and immunity-related markers were found. No changes in the frequencies of markers associated with two diseases, melanoma and insect bite hypersensitivity, were observed, due probably to the short interval of time between comparisons. It, thus, seems that, despite its small size, previous bottlenecks and inbreeding, the molecular variation of Old Kladruber horses is comparable to other horse breeds and that the current breeding policy does not compromise genetic variation of this endangered population.