Evaluation of Bone Growth around Autotransplanted Teeth Using Cone-Beam Computed Tomographic Images.

INTRODUCTION: The aim of this observational study was to evaluate bone growth at the recipient site of autotransplanted teeth using cone-beam computed tomography (CBCT) images and to compare it to the control tooth. METHODS: Bone growth was measured on CBCT images obtained from the anterior region of the maxilla, which received a single transplant of immature premolars due to early loss of the central incisor. The results were compared to the bone growth of contralateral central incisor that served as control. The sample consisted of 13 autotransplanted teeth and 13 control teeth. It was determined the thickness of the buccal bone of the marginal ridge, the thickness of the palatal bone of the marginal ridge, the height of the maxillary bone, and the thickness of the maxillary bone at its mid-height. The measurements were analyzed statistically, adopting a significance level of 5%. RESULTS: The mean interval between the dental autotransplantation surgery and the CBCT scan was 16.2 years. The mean thickness of the buccal bone of the marginal ridge, thickness of the palatal bone of the marginal ridge, height of the maxillary bone, and thickness of the maxillary bone at its mid-height were, respectively, 0.85, 0.94, 19.6, and 8.85 mm for autotransplanted teeth and 0.79, 1.09, 18.88, and 8.5 mm for control teeth. There were no significant differences in the bone measurements between autotransplanted and control teeth (P > .05). CONCLUSIONS: The growth of maxillary bone in the region that received the autotransplanted tooth was proportional to that of maxillary bone that supports the contralateral tooth. Thus, autogenous transplantation is efficient in maintaining the bone architecture at the recipient site.

A new HaCV-EBHSV recombinant lagovirus circulating in European brown hares (Lepus europaeus) from Catalonia, Spain.

In 2020/2021, several European brown hare syndrome virus (EBHSV) outbreaks were recorded in European hares (Lepus europaeus) from Catalonia, Spain. Recombination analysis combined with phylogenetic reconstruction and estimation of genetic distances of the complete coding sequences revealed that 5 strains were recombinants. The recombination breakpoint is located within the non-structural protein 2C-like RNA helicase (nucleotide position ~ 1889). For the genomic fragment upstream of the breakpoint, a non-pathogenic EBHSV-related strain (hare calicivirus, HaCV; GII.2) was the most closely related sequence; for the rest of the genome, the most similar strains were the European brown hare syndrome virus (EBHSV) strains recovered from the same 2020/2021 outbreaks, suggesting a recent origin. While the functional impact of the atypical recombination breakpoint remains undetermined, the novel recombinant strain was detected in different European brown hare populations from Catalonia, located 20-100 km apart, and seems to have caused a fatal disease both in juvenile and adult animals, confirming its viability and ability to spread and establish infection. This is the first report of a recombination event involving HaCV and EBHSV and, despite the recombination with a non-pathogenic strain, it appears to be associated with mortality in European brown hares, which warrants close monitoring.

An Ancestral Major Histocompatibility Complex Organization in Cartilaginous Fish: Reconstructing MHC Origin and Evolution.

Cartilaginous fish (sharks, rays, and chimeras) comprise the oldest living jawed vertebrates with a mammalian-like adaptive immune system based on immunoglobulins (Ig), T-cell receptors (TCRs), and the major histocompatibility complex (MHC). Here, we show that the cartilaginous fish "adaptive MHC" is highly regimented and compact, containing (i) a classical MHC class Ia (MHC-Ia) region containing antigen processing (antigen peptide transporters and immunoproteasome) and presenting (MHC-Ia) genes, (ii) an MHC class II (MHC-II) region (with alpha and beta genes) with linkage to beta-2-microglobulin (ß2m) and bromodomain-containing 2, (iii) nonclassical MHC class Ib (MHC-Ib) regions with 450 million-year-old lineages, and (iv) a complement C4 associated with the MHC-Ia region. No MHC-Ib genes were found outside of the elasmobranch MHC. Our data suggest that both MHC-I and MHC-II genes arose after the second round of whole-genome duplication (2R) on a human chromosome (huchr) 6 precursor. Further analysis of MHC paralogous regions across early branching taxa from all jawed vertebrate lineages revealed that Ig/TCR genes likely arose on a precursor of the huchr9/12/14 MHC paralog. The ß2m gene is linked to the Ig/TCR genes in some vertebrates suggesting that it was present at 1R, perhaps as the donor of C1 domain to the primordial MHC gene. In sum, extant cartilaginous fish exhibit a conserved and prototypical MHC genomic organization with features found in various vertebrates, reflecting the ancestral arrangement for the jawed vertebrates.

First detection and molecular characterization of rabbit hemorrhagic disease virus (RHDV) in Algeria.

Since the first detection of rabbit hemorrhagic disease (RHD), the rabbit hemorrhagic disease virus (RHDV) has been responsible for high morbidity and mortality worldwide, both in domestic and in wild rabbits. Despite the apparent control of RHD in rabbitries through vaccination, several studies highlighted the rapid evolution of RHDV by recombination, which may facilitate the emergence of new pathogenic strains. The aim of this study was to confirm the presence and characterize RHDV in Algeria. For this, rabbit samples were collected in the north of Algeria, between 2018 and 2021, from small farms where the virus was suspected after the sudden death of a high number of rabbits, and from healthy hunted wild rabbits. The domestic rabbits revealed clinical signs and lesions that were suggestive of RHD. RT-PCR showed that 79.31% of the domestic rabbit samples were positive for RHDV, while in 20.69%, including the hunted rabbits, the virus was not detected. Phylogenetic analysis of the Algerian strains allowed the confirmation and identification as GI.2 (RHDV2), and showed a close relation to GI.3P-GI.2 recombinant strains, suggesting a potential introduction from other countries, with an older strain potentially originated from neighboring Tunisia, while more recent isolates grouped with strains from North America. Our study reports for the first time the presence of GI.2 (RHDV2) in Algeria with multiple routes of introduction. Consequently, we propose that RHDV control in Algeria should be based on epidemiological surveys in association with an adequate prophylactic program.

Corrigendum: Uncovering a 500 million year old history and evidence of pseudogenization for TLR15.

[This corrects the article DOI: 10.3389/fimmu.2022.1020601.].

The potential role of scavenging flies as mechanical vectors of Lagovirus europaeus/GI.2.

The European rabbit (Oryctolagus cuniculus) populations of the Iberian Peninsula have been severely affected by the emergence of the rabbit haemorrhagic disease virus (RHDV) Lagovirus europaeus/GI.2 (RHDV2/b). Bushflies and blowflies (Muscidae and Calliphoridae families, respectively) are important RHDV vectors in Oceania, but their epidemiological role is unknown in the native range of the European rabbit. In this study, scavenging flies were collected between June 2018 and February 2019 in baited traps at one site in southern Portugal, alongside a longitudinal capture-mark-recapture study of a wild European rabbit population, aiming to provide evidence of mechanical transmission of GI.2 by flies. Fly abundance, particularly from Calliphoridae and Muscidae families, peaked in October 2018 and in February 2019. By employing molecular tools, we were able to detect the presence of GI.2 in flies belonging to the families Calliphoridae, Muscidae, Fanniidae and Drosophilidae. The positive samples were detected during an RHD outbreak and absent in samples collected when no evidence of viral circulation in the local rabbit population was found. We were able to sequence a short viral genomic fragment, confirming its identity as RHDV GI.2. The results suggest that scavenging flies may act as mechanical vectors of GI.2 in the native range of the southwestern Iberian subspecies O. cuniculus algirus. Future studies should better assess their potential in the epidemiology of RHD and as a tool for monitoring viral circulation in the field.

Extensive MHC class IIß diversity across multiple loci in the small-spotted catshark (Scyliorhinus canicula).

The major histocompatibility complex (MHC) is a multigene family responsible for pathogen detection, and initiation of adaptive immune responses. Duplication, natural selection, recombination, and their resulting high functional genetic diversity spread across several duplicated loci are the main hallmarks of the MHC. Although these features were described in several jawed vertebrate lineages, a detailed MHC IIß characterization at the population level is still lacking for chondrichthyans (chimaeras, rays and sharks), i.e. the most basal lineage to possess an MHC-based adaptive immune system. We used the small-spotted catshark (Scyliorhinus canicula, Carcharhiniformes) as a case-study species to characterize MHC IIß diversity using complementary molecular tools, including publicly available genome and transcriptome datasets, and a newly developed high-throughput Illumina sequencing protocol. We identified three MHC IIß loci within the same genomic region, all of which are expressed in different tissues. Genetic screening of the exon 2 in 41 individuals of S. canicula from a single population revealed high levels of sequence diversity, evidence for positive selection, and footprints of recombination. Moreover, the results also suggest the presence of copy number variation in MHC IIß genes. Thus, the small-spotted catshark exhibits characteristics of functional MHC IIß genes typically observed in other jawed vertebrates.

Uncovering a 500 million year old history and evidence of pseudogenization for TLR15.

Introduction: Toll like receptors (TLRs) are at the front line of pathogen recognition and host immune response. Many TLR genes have been described to date with some being found across metazoans while others are restricted to specific lineages. A cryptic member of the TLR gene family, TLR15, has a unique phylogenetic distribution. Initially described in extant species of birds and reptiles, an ortholog has been reported for cartilaginous fish. Methods: Here, we significantly expanded the evolutionary analysis of TLR15 gene evolution, taking advantage of large genomic and transcriptomic resources available from different lineages of vertebrates. Additionally, we objectively search for TLR15 in lobe-finned and ray-finned fish, as well as in cartilaginous fish and jawless vertebrates. Results and discussion: We confirm the presence of TLR15 in early branching jawed vertebrates - the cartilaginous fish, as well as in basal Sarcopterygii - in lungfish. However, within cartilaginous fish, the gene is present in Holocephalans (all three families) but not in Elasmobranchs (its sister-lineage). Holocephalans have long TLR15 protein sequences that disrupt the typical TLR structure, and some species display a pseudogene sequence due to the presence of frameshift mutations and early stop codons. Additionally, TLR15 has low expression levels in holocephalans when compared with other TLR genes. In turn, lungfish also have long TLR15 protein sequences but the protein structure is not compromised. Finally, TLR15 presents several sites under negative selection. Overall, these results suggest that TLR15 is an ancient TLR gene and is experiencing ongoing pseudogenization in early-branching vertebrates.

Identification of the Fc-alpha/mu receptor in Xenopus provides insight into the emergence of the poly-Ig receptor (pIgR) and mucosal Ig transport.

The polyimmunoglobulin receptor (pIgR) transcytoses J chain-containing antibodies through mucosal epithelia. In mammals, two cis-duplicates of PIGR, FCMR, and FCAMR, flank the PIGR gene. A PIGR duplication is first found in amphibians, previously annotated as PIGR2 (herein xlFCAMR), and is expressed by APCs. We demonstrate that xlFcamR is the equivalent of mammalian FcamR. It has been assumed that pIgR is the oldest member of this family, yet our data could not distinguish whether PIGR or FCAMR emerged first; however, FCMR was the last family member to emerge. Interestingly, bony fish "pIgR" is not an orthologue of tetrapod pIgR, and possibly acquired its function via convergent evolution. PIGR/FCAMR/FCMR are members of a larger superfamily, including TREM, CD300, and NKp44, which we name the "double-disulfide Ig superfamily" (ddIgSF). Domains related to each ddIgSF family were identified in cartilaginous fish (sharks, chimeras) and encoded in a single gene cluster syntenic to the human pIgR locus. Thus, the ddIgSF families date back to the earliest antibody-based adaptive immunity, but apparently not before. Finally, our data strongly suggest that the J chain arose in evolution only for Ig multimerization. This study provides a framework for further studies of pIgR and the ddIgSF in vertebrates.

A Highly Complex, MHC-Linked, 350 Million-Year-Old Shark Nonclassical Class I Lineage.

Cartilaginous fish, or Chondrichthyes, are the oldest extant vertebrates to possess the MHC and the Ig superfamily-based Ag receptors, the defining genes of the gnathostome adaptive immune system. In this work, we have identified a novel MHC lineage, UEA, a complex multigene nonclassical class I family found in sharks (division Selachii) but not detected in chimaeras (subclass Holocephali) or rays (division Batoidea). This new lineage is distantly related to the previously reported nonclassical class I lineage UCA, which appears to be present only in dogfish sharks (order Squaliformes). UEA lacks conservation of the nine invariant residues in the peptide (ligand)-binding regions (PBR) that bind to the N and C termini of bound peptide in most vertebrate classical class I proteins, which are replaced by relatively hydrophobic residues compared with the classical UAA. In fact, UEA and UCA proteins have the most hydrophobic-predicted PBR of all identified chondrichthyan class I molecules. UEA genes detected in the whale shark and bamboo shark genome projects are MHC linked. Consistent with UEA comprising a very large gene family, we detected weak expression in different tissues of the nurse shark via Northern blotting and RNA sequencing. UEA genes fall into three sublineages with unique characteristics in the PBR. UEA shares structural and genetic features with certain nonclassical class I genes in other vertebrates, such as the highly complex XNC nonclassical class I genes in Xenopus, and we anticipate that each shark gene, or at least each sublineage, will have a unique function, perhaps in bacterial defense.

Analysis of shark NCR3 family genes reveals primordial features of vertebrate NKp30.

Natural killer (NK) cells play major roles in innate immunity against viruses and cancer. Natural killer receptors (NKR) expressed by NK cells recognize foreign- or self-ligands on infected and transformed cells as well as healthy cells. NKR genes are the most rapidly evolving loci in vertebrates, and it is generally difficult to detect orthologues in different taxa. The unique exception is NKp30, an activating NKR in mammals that binds to the self-ligand B7H6. The NKp30-encoding gene, NCR3, has been found in most vertebrates including sharks, the oldest vertebrates with human-type adaptive immunity. NCR3 has a special, non-rearranging VJ-type immunoglobulin superfamily (IgSF) domain that predates the emergence of the rearranging antigen receptors. Herein we show that NCR3 loci are linked to the shark major histocompatibility complex (MHC), proving NCR3's primordial association with the MHC. We identified eight subtypes of differentially expressed highly divergent shark NCR3 family genes. Using in situ hybridization, we detected one subtype, NS344823, to be expressed by predominantly single cells outside of splenic B cell zones. The expression by non-B cells was also confirmed by PCR in peripheral blood lymphocytes. Surprisingly, high expression of NS344823 was detected in the thymic cortex, demonstrating NS344823 expression in developing T cells. Finally, we show for the first time that shark T cells are found as single cells or in small clusters in the splenic red pulp, also unassociated with the large B cell follicles we previously identified.

Cartilaginous fish class II genes reveal unprecedented old allelic lineages and confirm the late evolutionary emergence of DM.

Cartilaginous fish (chimaeras, rays and sharks) are the most basal extant jawed vertebrates with an adaptive immune system based on the Major Histocompatibility Complex (MHC). Despite being a key taxon in the evolution of vertebrate adaptive immunity, no comprehensive characterization of MHC class II genes has been undertaken for the group. We performed extensive bioinformatic searches on a taxonomically diverse dataset of transcriptomes and genomes of cartilaginous fish targeting MHC class II sequences. Class IIα and IIß sequences were retrieved from all taxa analyzed and showed typical features of classical class II genes. Phylogenetic trees of the immunoglobulin superfamily domain showed two divergent and remarkably ancient lineages of class II genes in Selachians (sharks), originating >350 million years ago. Close linkage of lineage-specific pairs of IIα and IIß genes was found, confirming previous results, with genes from distinct lineages segregating as alleles. Nonclassical class II DM sequences were not retrieved from these data and classical class II sequences lacked the conserved residues shown to interact with DM molecules, supporting claims that the DM system arose only in the lobe-finned fish lineage leading to tetrapods. Based on our search methods, other divergent class II genes are unlikely in cartilaginous fish.

Cartilaginous fishes offer unique insights into the evolution of the nuclear receptor gene repertoire in gnathostomes.

Nuclear receptors (NRs) are key transcription factors that originated in the common ancestor of metazoans. The vast majority of NRs are triggered by binding to either endogenous (e.g. retinoic acid) or exogenous (e.g. xenobiotics) ligands, and their evolution and expansion is tightly linked to the function of endocrine systems. Importantly, they represent classic targets of physiological exploitation by endocrine disrupting chemicals. The NR gene repertoire in different lineages has been shaped by gene loss, duplication and mutation, denoting a dynamic evolutionary route. As the earliest diverging class of gnathostomes (jawed vertebrates), cartilaginous fishes offer an exceptional opportunity to address the early diversification of NR gene families and the evolution of the endocrine system in jawed vertebrates. Here we provide an exhaustive analysis into the NR gene composition in five elasmobranch (sharks and rays) and two holocephalan (chimaeras) species. For this purpose, we generated also a low coverage draft genome assembly of the chimaera small-eyed rabbitfish, Hydrolagus affinis. We show that cartilaginous fish retain an archetypal NR gene repertoire, similar to that of mammals and coincident with the two rounds of whole genome duplication that occurred in the gnathostome ancestor. Furthermore, novel gene members of the non-canonical NR0B receptors were found in the genomes of this lineage. Our findings provide an essential view into the early diversification of NRs in gnathostomes, paving the way for functional studies.

An Ancient, MHC-Linked, Nonclassical Class I Lineage in Cartilaginous Fish.

Cartilaginous fishes, or chondrichthyans, are the oldest jawed vertebrates that have an adaptive immune system based on the MHC and Ig superfamily-based AgR. In this basal group of jawed vertebrates, we identified a third nonclassical MHC class I lineage (UDA), which is present in all species analyzed within the two major cartilaginous subclasses, Holocephali (chimaeras) and Elasmobranchii (sharks, skates, and rays). The deduced amino acid sequences of UDA have eight out of nine typically invariant residues that bind to the N and C termini of bound peptide found in most vertebrae classical class I (UAA); additionally, the other predicted 28 peptide-binding residues are perfectly conserved in all elasmobranch UDA sequences. UDA is distinct from UAA in its differential tissue distribution and its lower expression levels and is mono- or oligomorphic unlike the highly polymorphic UAA UDA has a low copy number in elasmobranchs but is multicopy in the holocephalan spotted ratfish (Hydrolagus colliei). Using a nurse shark (Ginglymostoma cirratum) family, we found that UDA is MHC linked but separable by recombination from the tightly linked cluster of UAA, TAP, and LMP genes, the so-called class I region found in most nonmammalian vertebrates. UDA has predicted structural features that are similar to certain nonclassical class I genes in other vertebrates, and, unlike polymorpic classical class I, we anticipate that it may bind to a conserved set of specialized peptides.

Genetic Diversity of IGHM and IGHE in the Leporids Revealed Different Patterns of Diversity in the Two European Rabbit Subspecies (O. cuniculus algirus and O. c. cuniculus).

The European rabbit (Oryctolagus cuniculus) has been an important model for immunological studies but the study of its immunoglobulins (Ig) has been restricted to its unique IgA and IgG. Here, we studied the genetic diversity of IgM and IgE in several species of leporids and performed population genetics studies on European rabbit wild populations and domestic breeds. The leporids sequencing showed that these Ig are well conserved (98% sequence similarity among leporids), For IgM the Cµ1 and Cµ4 were the most diverse and most conserved domains, respectively, while for IgE the Cε1 was the most diverse domain and Cε2 and Cε3 the most conserved domains. The differences in the pattern of most conserved and most diverse domain between the Ig isotypes are most likely related to each isotype function. The genetic population data showed contrasting results for IgM and IgE. For both Ig, as expected, a greater diversity was observed in the original species range, the Iberian Peninsula. However, unexpectedly the genetic diversity found for IgE in the domestic animals is higher than that for the French wild populations. These results will increase knowledge of the genetic diversity of leporids and wild and domestic rabbit populations and are important tools for the management of wild populations and rabbitries.

Trichomonas Infection in a Community of Free-Ranging Domestic and Wild Columbiformes and Bonelli's Eagle (Aquila fasciata).

Trichomonas gallinae is a pathogen of conservation relevance, whose main maintenance hosts are Columbiformes, but spillover to avian predators has been described. The goal of this study was to characterize the epidemiology of Trichomonas spp. in a community of free-ranging domestic and wild Columbiformes and an endangered predator, Bonelli's eagle Aquila fasciata. We surveyed 253 live-captured Rock doves, 16 nestling Bonelli's eagles and 41 hunted Columbiformes. Oro-esophageal swabs were incubated in culture media and Trichomonas spp. isolated from Bonelli's eagle (6.3%, CI95 1.1-28.3), Turtle dove Streptopelia turtur (56.3%, CI95 39.3-71.8), Wood pigeon Columba palumbus (83.3%, CI95 43.7-97.0) and Rock dove Columba livia (68.4%, CI95 62.4-73.8). Infected Rock doves showed significantly poorer body condition than uninfected ones (p = 0.022). From a subset of 32 isolates, 18S and ITS1/5.8S/ITS2 rRNA genes were sequenced and Maximum-Likelihood trees inferred. Four ribotypes of Trichomonas spp. were identified. In this study area Trichomonas spp. seem to persist in a multi-host system involving several species of Columbiformes. Conservation actions aimed at increasing the availability of trophic resources for Bonelli's eagles through Rock dove restocking should consider the risk of pathogen transmission and of introduction of alien strains.

Sequencing of VDJ genes in Lepus americanus confirms a correlation between VHn expression and the leporid species continent of origin.

Leporid VH genes used in the generation of their primary antibody repertoire exhibit highly divergent lineages. For the European rabbit (Oryctolagus cuniculus) four VHa lineages have been described, the a1, a2, a3 and a4. Hares (Lepus spp.) and cottontail (Sylvilagus floridanus) express one VHa lineage each, the a2L and the a5, respectively, along with a more ancient lineage, the Lepus spp. sL and S. floridanus sS. Both the European rabbit and the Lepus europaeus use a third lineage, VHn, in a low proportion of their VDJ rearrangements. The VHn genes are a conserved ancestral polymorphism that is being maintained in the leporid genome.Their usage in a low proportion of VDJ rearrangements by both European rabbit and L. europaeus but not S. floridanus has been argued to be a remnant of an ancient European leporid immunologic response to pathogens. To address this hypothesis, in this study we sequenced VDJ rearranged genes for another North American leporid, L. americanus. Our results show that L. americanus expressed these genes less frequently and in a highly modified fashion compared to the European Lepus species. Our results suggest that the American leporid species use a different VH repertoire than the European species which may be related with an immune adaptation to different environmental conditions, such as different pathogenic agents.

De novo assembly of the kidney and spleen transcriptomes of the cosmopolitan blue shark, Prionace glauca.

Cartilaginous fishes (sharks, rays and chimaeras) comprise a highly diversified group of basal vertebrates occupying a plethora of ecological aquatic niches. They represent critical components of marine ecosystems and food webs, although numerous species are threatened and almost half are poorly known. Genomic resources emerging from this basal jawed vertebrate group have offered valuable insights into the evolution of vertebrate-specific traits. Yet, the taxon remains largely understudied. Here, we generated the first high-quality de novo assembly of kidney and spleen transcriptomes of the blue shark (Prionace glauca). A total of 32,917,412 and 52,666,542 reads were obtained for spleen and kidney, respectively, using RNA-Seq Illumina technology. De novo multi-tissue assembly resulted in 97,317 unigenes with an N50 of 1975bp, in which 87,571 were assigned to a particular tissue or combination of tissues based on the sequencing read mapping. Functional annotation generated 28,564 and 19,854 open reading frames in spleen and kidney, respectively. This dataset provides a significant resource for physiological and evolutionary studies, namely into the unique osmoregulatory system of Chondrichthyes and the evolution of the immune system in vertebrates.

Proposal for a unified classification system and nomenclature of lagoviruses.

Lagoviruses belong to the Caliciviridae family. They were first recognized as highly pathogenic viruses of the European rabbit (Oryctolagus cuniculus) and European brown hare (Lepus europaeus) that emerged in the 1970-1980s, namely, rabbit haemorrhagic disease virus (RHDV) and European brown hare syndrome virus (EBHSV), according to the host species from which they had been first detected. However, the diversity of lagoviruses has recently expanded to include new related viruses with varying pathogenicity, geographic distribution and host ranges. Together with the frequent recombination observed amongst circulating viruses, there is a clear need to establish precise guidelines for classifying and naming lagovirus strains. Therefore, here we propose a new nomenclature based on phylogenetic relationships. In this new nomenclature, a single species of lagovirus would be recognized and called Lagovirus europaeus. The species would be divided into two genogroups that correspond to RHDV- and EBHSV-related viruses, respectively. Genogroups could be subdivided into genotypes, which could themselves be subdivided into phylogenetically well-supported variants. Based on available sequences, pairwise distance cutoffs have been defined, but with the accumulation of new sequences these cutoffs may need to be revised. We propose that an international working group could coordinate the nomenclature of lagoviruses and any proposals for revision.

Genetic characterization of interleukins (IL-1α, IL-1ß, IL-2, IL-4, IL-8, IL-10, IL-12A, IL-12B, IL-15 and IL-18) with relevant biological roles in lagomorphs.

ILs, as essential innate immune modulators, are involved in an array of biological processes. In the European rabbit (Oryctolagus cuniculus) IL-1α, IL-1ß, IL-2, IL-4, IL-8, IL-10, IL-12A, IL-12B, IL-15 and IL-18 have been implicated in inflammatory processes and in the immune response against rabbit hemorrhagic disease virus and myxoma virus infections. In this study we characterized these ILs in six Lagomorpha species (European rabbit, pygmy rabbit, two cottontail rabbit species, European brown hare and American pika). Overall, these ILs are conserved between lagomorphs, including in their exon/intron structure. Most differences were observed between leporids and American pika. Indeed, when comparing both, some relevant differences were observed in American pika, such as the location of the stop codon in IL-1α and IL-2, the existence of a different transcript in IL8 and the number of cysteine residues in IL-1ß. Changes at N-glycosylation motifs were also detected in IL-1, IL-10, IL-12B and IL-15. IL-1α is the protein that presents the highest evolutionary distances, which is in contrast to IL-12A where the distances between lagomorphs are the lowest. For all these ILs, sequences of human and European rabbit are more closely related than between human and mouse or European rabbit and mouse.