A loss-of-function mutation in RORB disrupts saltatorial locomotion in rabbits.

Saltatorial locomotion is a type of hopping gait that in mammals can be found in rabbits, hares, kangaroos, and some species of rodents. The molecular mechanisms that control and fine-tune the formation of this type of gait are unknown. Here, we take advantage of one strain of domesticated rabbits, the sauteur d'Alfort, that exhibits an abnormal locomotion behavior defined by the loss of the typical jumping that characterizes wild-type rabbits. Strikingly, individuals from this strain frequently adopt a bipedal gait using their front legs. Using a combination of experimental crosses and whole genome sequencing, we show that a single locus containing the RAR related orphan receptor B gene (RORB) explains the atypical gait of these rabbits. We found that a splice-site mutation in an evolutionary conserved site of RORB results in several aberrant transcript isoforms incorporating intronic sequence. This mutation leads to a drastic reduction of RORB-positive neurons in the spinal cord, as well as defects in differentiation of populations of spinal cord interneurons. Our results show that RORB function is required for the performance of saltatorial locomotion in rabbits.

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.

Functional cross-species conservation of guanylate-binding proteins in innate immunity.

Guanylate binding proteins (GBPs) represent an evolutionary ancient protein family widely distributed among eukaryotes. They are interferon (IFN)-inducible guanosine triphosphatases that belong to the dynamin superfamily. GBPs are known to have a major role in the cell-autonomous innate immune response against bacterial, parasitic and viral infections and are also involved in inflammasome activation. Evolutionary studies depicted that GBPs present a pattern of gain and loss of genes in each family with several genes pseudogenized and some genes more divergent, indicative for the birth-and-death evolution process. Most species harbor large GBP gene clusters encoding multiple paralogs. Previous functional studies mainly focused on mouse and human GBPs, but more data are becoming available, broadening the understanding of this multifunctional protein family. In this review, we will provide new insights and give a broad overview about GBP evolution, conservation and their roles in all studied species, including plants, invertebrates and vertebrates, revealing how far the described features of GBPs can be transferred to other species.

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.

Evidence for Extensive Duplication and Subfunctionalization of FCRL6 in Armadillo (Dasypus novemcinctus).

The control of infections by the vertebrate adaptive immune system requires careful modulation to optimize defense and minimize harm to the host. The Fc receptor-like (FCRL) genes encode immunoregulatory molecules homologous to the receptors for the Fc portion of immunoglobulin (FCR). To date, nine different genes (FCRL1-6, FCRLA, FCRLB and FCRLS) have been identified in mammalian organisms. FCRL6 is located at a separate chromosomal position from the FCRL1-5 locus, has conserved synteny in mammals and is situated between the SLAMF8 and DUSP23 genes. Here, we show that this three gene block underwent repeated duplication in Dasypus novemcinctus (nine-banded armadillo) resulting in six FCRL6 copies, of which five appear functional. Among 21 mammalian genomes analyzed, this expansion was unique to D. novemcinctus. Ig-like domains that derive from the five clustered FCRL6 functional gene copies show high structural conservation and sequence identity. However, the presence of multiple non-synonymous amino acid changes that would diversify individual receptor function has led to the hypothesis that FCRL6 endured subfunctionalization during evolution in D. novemcinctus. Interestingly, D. novemcinctus is noteworthy for its natural resistance to the Mycobacterium leprae pathogen that causes leprosy. Because FCRL6 is chiefly expressed by cytotoxic T and NK cells, which are important in cellular defense responses against M. leprae, we speculate that FCRL6 subfunctionalization could be relevant for the adaptation of D. novemcinctus to leprosy. These findings highlight the species-specific diversification of FCRL family members and the genetic complexity underlying evolving multigene families critical for modulating adaptive immune protection.

TLR7 and TLR8 evolution in lagomorphs: different patterns in the different lineages.

Toll-like receptors (TLRs) are one of the most ancient and widely studied innate immune receptors responsible for host defense against invading pathogens. Among the known TLRs, TLR7 and TLR8 sense and recognize single-stranded (ss) RNAs with a dynamic evolutionary history. While TLR8 was lost in birds and duplicated in turtles and crocodiles, TLR7 is duplicated in some birds, but in other tetrapods, there is only one copy. In mammals, with the exception of lagomorphs, TLR7 and TLR8 are highly conserved. Here, we aim to study the evolution of TLR7 and TLR8 in mammals, with a special focus in the order Lagomorpha. By searching public sequence databases, conducting evolutionary analysis, and evaluating gene expression, we were able to confirm that TLR8 is absent in hares but widely expressed in the European rabbit. In contrast, TLR7 is absent in the European rabbit and quite divergent in hares. Our results suggest that, in lagomorphs, more in particular in leporids, TLR7 and TLR8 genes have evolved faster than in any other mammalian group. The long history of interaction with viruses and their location in highly dynamic telomeric regions might explain the pattern observed.

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.

Physiological Responses and Technical-Tactical Performance of Youth Basketball Players: A Brief Comparison between 3x3 and 5x5 Basketball.

This study aims to examine youth players' physiological responses and technical-tactical performance when playing simulated 3x3 and 5x5 basketball games. Fifteen well-trained male basketball players (16.6 ± 0.2 years old) participated in scrimmage basketball games under two different conditions: 3x3 (half-court) and 5x5 (full-court). The players' heart rate, muscle oxygen saturation and total hemoglobin data were collected and computed to describe physiological responses, while video analysis was used to characterize their technical-tactical performance. A Bayesian one-way analysis of variance (ANOVA) was used to quantify the predictive influence of both game conditions on the physiological and the technical-tactical variables. The results indicated that different game conditions influenced the players' physiological responses slightly, as only hemoglobin sample entropy increased between the 3x3 and 5x5 game scenarios. Conversely, statistical differences in most of the technical-tactical variables were moderate and decisive in favour of the game condition model. Overall, this study emphasizes that playing 3x3 and 5x5 basketball games lead to relatively negligible differences in the players' physiological response but pronounced variations in their technical-tactical performance. Therefore, important implications may be drawn to the applied field as the specificity of technical-tactical adaptations when playing 3x3 or 5x5 formats should be considered by basketball coaches to better design the training sessions for players that fall within our sample age category.

How manipulation of playing area dimensions in ball possession games constrains physical effort and technical actions in under-11, under-15 and under-23 soccer players.

The aim of this study was to investigate the effects of playing area manipulation (20 × 15 m, 25 × 20 m and 30 × 25 m) on external workloads (total distance covered, distance covered while walking, running and sprinting, number of sprints, maximum sprint speed), internal load perceptions (rating of perceived exertion) and technical actions of passing (number of passes with dominant and non-dominant foot, and maximum passing speed) during 4v4 ball possession small-sided and conditioned games in under-11, under-15 and under-23 soccer players. Results showed higher values in the large playing area for under-11 in the distance covered in different speed zones, sprint number and RPE (all p <.001) for under-15 in sprints number (p <.01) and maximum sprint speed (p =.02), and for under-23 in both RPE and sprint numbers (p <.01). Although no significant differences were found on technical actions, it was still possible to notice some effects through pairwise comparison. High-intensity running was promoted on larger playing areas, where under-11 s were also able to perform more technical actions of passing. Opposite, under-23s were able to perform more passing on smaller playing areas, where under-11 s perceived the exercise more intense. The impact of different playing areas was reduced for the under-15.

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.

Arousal changes and delirium in acute medically-ill male older patients with and without dementia: a prospective study during hospitalization.

Objectives: Previous research has characterized the prevalence, natural course and outcomes of delirium superimposed in dementia but much less is known about the relation between preexisting dementia and the emergence of altered arousal (such as drowsiness, obtundation, stupor or agitation) during acute medical illness. This study aimed to determine the natural course of delirium and abnormal arousal states in acute medically-ill older patients with and without prior dementia during hospital stay.Methods: Observational prospective study in an acute male geriatric ward. Patients aged ≥ 65 years old were assessed by a psychiatrist within the first 72h of admission and in every other day until discharge to determine the level of arousal and the presence of delirium. Prior cognitive impairment, sociodemographic data, chronic comorbidities, psychotropic prescription and functional status were assessed at baseline.Results: 43.5% of participants in the final sample (n= 269) had dementia. Prior dementia was associated with higher rates of moderate/severe hypoarousal (29.9% vs. 4.6%; p<0.001) and delirium (20.5% vs. 7.2%; p<0.001) at admission. RASS ≤ -3 at admission predicted a 4-fold increased intra-hospital mortality risk and RASS ≠ 0 had a sensitivity of 82.8% and a specificity of85.9% for delirium.Conclusions: Moderate/severe hypoarousal is associated with adverse outcomes and should be assessed as part of delirium spectrum, particularly in subjects with prior dementia.

The evolution of S100A7: an unusual gene expansion in Myotis bats.

BACKGROUND: The S100A7 gene, also called psoriasin, was first described as an upregulated protein in psoriatic skin. For the past years, the importance of this protein as a key effector of innate immunity has been clearly established, not only due to its importance protecting against bacteria skin insult in humans, but also because of its important role in amplifying inflammatory processes. Given the importance of S100A7 in host defense, S100A7 genes have been mostly studied in humans. Here we provide a detailed analysis of the evolution of the gene family encoding for the S100A7 protein in mammals. RESULTS: Examination of several mammalian genomes revealed an unexpected variation in the copy number of S100A7. Among the most representative mammalian groups, we report that multiple events of duplication, gene loss and high mutation rates are shaping the evolution of this gene family. An unexpected result comes from Myotis species (order Chiroptera), where we found an outstanding S100A7 gene radiation, resulting in more than 10 copies in M. lucifugus and 5 copies in M. brandtii. These findings suggest a unique adaptive road in these species and are suggestive of special role of this protein in their immune system. CONCLUSIONS: We found different evolutionary histories among different mammalian groups. Overall, our results suggest that this gene family is evolving under the birth-and-death model of evolution. To our knowledge, this work represents the first detailed analysis of phylogenetic relationships of S100A7 within mammals and therefore will pave the way to further clarify their unique function in the immune system.

Analysis of substitution rates showed that TLR5 is evolving at different rates among mammalian groups.

BACKGROUND: Toll-like receptors (TLRs) are the most widely studied innate immunity receptors responsible for recognition of invading pathogens. Among the TLR family, TLR5 is the only that senses and recognizes flagellin, the major protein of bacterial flagella. TLR5 has been reported to be under overall purifying selection in mammals, with a small proportion of codons under positive selection. However, the variation of substitution rates among major mammalian groups has been neglected. Here, we studied the evolution of TLR5 in mammals, comparing the substitution rates among groups. RESULTS: In this study we analysed the TLR5 substitution rates in Euungulata, Carnivora, Chiroptera, Primata, Rodentia and Lagomorpha, groups. For that, Tajima's relative rate test, Bayesian inference of evolutionary rates and genetic distances were estimated with CODEML's branch model and RELAX. The combined results showed that in the Lagomorpha, Rodentia, Carnivora and Chiroptera lineages TLR5 is evolving at a higher substitution rate. The RELAX analysis further suggested a significant relaxation of selective pressures for the Lagomorpha (K = 0.22, p < 0.01), Rodentia (K = 0.58, p < 0.01) and Chiroptera (K = 0.65, p < 0.01) lineages and for the Carnivora ancestral branches (K = 0.13, p < 0.01). CONCLUSIONS: Our results show that the TLR5 substitution rate is not uniform among mammals. In fact, among the different mammal groups studied, the Lagomorpha, Rodentia, Carnivora and Chiroptera are evolving faster. This evolutionary pattern could be explained by 1) the acquisition of new functions of TLR5 in the groups with higher substitution rate, i.e. TLR5 neofunctionalization, 2) by the beginning of a TLR5 pseudogenization in these groups due to some redundancy between the TLRs genes, or 3) an arms race between TLR5 and species-specific parasites.

Evolution of CCL16 in Glires (Rodentia and Lagomorpha) shows an unusual random pseudogenization pattern.

BACKGROUND: The C-C motif chemokine ligand 16 (CCL16) is a potent pro-inflammatory chemokine and a chemoattractant for monocytes and lymphocytes. In normal plasma, it is present at high concentrations and elicits its effects on cells by interacting with cell surface chemokine receptors. In the European rabbit and in rodents such as mouse, rat and guinea pig, CCL16 was identified as a pseudogene, while in the thirteen-lined ground squirrel it appears to be potentially functional. To gain insight into the evolution of this gene in the superorder Glires (rodents and lagomorphs), we amplified the CCL16 gene from eleven Leporidae and seven Ochotonidae species. RESULTS: We compared our sequences with CCL16 sequences of twelve rodent species retrieved from public databases. The data show that for all leporid species studied CCL16 is a pseudogene. This is primarily due to mutations at the canonical Cys Cys motif, creating either premature stop codons, or disrupting amino acid replacements. In the Mexican cottontail, CCL16 is pseudogenized due to a frameshift deletion. Additionally, in the exon 1 (signal peptide), there are frameshift deletions present in all leporids studied. In contrast, in Ochotona species, CCL16 is potentially functional, except for an allele in Hoffmann's pika. In rodents, CCL16 is functional in a number of species, but patterns of pseudogenization similar to those observed in lagomorphs also exist. CONCLUSIONS: Our results suggest that while functional in the Glires ancestor, CCL16 underwent pseudogenization in some species. This process occurred stochastically or in specific lineages at different moments in the evolution of Glires. These observations suggest that the CCL16 had different evolutionary constrains in the Glires group that could be associated with the CCL16 biological function.

Alternated selection mechanisms maintain adaptive diversity in different demographic scenarios of a large carnivore.

BACKGROUND: Different population trajectories are expected to impact the signature of neutral and adaptive processes at multiple levels, challenging the assessment of the relative roles of different microevolutionary forces. Here, we integrate adaptive and neutral variability patterns to disentangle how adaptive diversity is driven under different demographic scenarios within the Iberian wolf (Canis lupus) range. We studied the persistent, the expanding and a small, isolated group within the Iberian wolf population, using 3 MHC class II genes (DRB1, DQA1, and DQB1), which diversity was compared with 39 microsatellite loci. RESULTS: Both the persistent and the expanding groups show evidence of balancing selection, revealed by a significant departure from neutrality at MHC loci, significant higher observed and expected heterozygosity and lower differentiation at MHC than at neutral loci, and signs of positive selection. However, despite exhibiting a significantly higher genetic diversity than the isolated group, the persistent group did not show significant excess of MHC heterozygotes. The expanding group, while showing a similar level of genetic diversity than the persistent group, displays by contrast a significant excess of MHC heterozygotes, which is compatible with the heterozygote advantage mechanism. Results are not clear regarding the role of drift and selection in the isolated group due to the small size of this population. Although diversity indices of MHC loci correspond to neutral expectations in the isolated group, accelerated MHC divergence, revealed by a higher differentiation at MHC than neutral loci, may indicate diversifying selection. CONCLUSION: Different selective pressures were observed in the three different demographic scenarios, which are possibly driven by different selection mechanisms to maintain adaptive diversity.

ICTV Virus Taxonomy Profile: Caliciviridae.

The family Caliciviridae includes viruses with single-stranded, positive-sense RNA genomes of 7.4-8.3 kb. The most clinically important representatives are human noroviruses, which are a leading cause of acute gastroenteritis in humans. Virions are non-enveloped with icosahedral symmetry. Members of seven genera infect mammals (Lagovirus, Norovirus, Nebovirus, Recovirus, Sapovirus, Valovirus and Vesivirus), members of two genera infect birds (Bavovirus and Nacovirus), and members of two genera infect fish (Minovirus and Salovirus). This is a summary of the International Committee on Taxonomy of Viruses (ICTV) Report on the family Caliciviridae, which is available at ictv.global/report/caliciviridae.

The antiviral activity of rodent and lagomorph SERINC3 and SERINC5 is counteracted by known viral antagonists.

A first step towards the development of a human immunodeficiency virus (HIV) animal model has been the identification and surmounting of species-specific barriers encountered by HIV along its replication cycle in cells from small animals. Serine incorporator proteins 3 (SERINC3) and 5 (SERINC5) were recently identified as restriction factors that reduce HIV-1 infectivity. Here, we compared the antiviral activity of SERINC3 and SERINC5 among mice, rats and rabbits, and their susceptibility to viral counteraction to their human counterparts. In the absence of viral antagonists, rodent and lagomorph SERINC3 and SERINC5 displayed anti-HIV activity in a similar range to human controls. Vesicular stomatitis virus G protein (VSV-G) pseudotyped virions were considerably less sensitive to restriction by all SERINC3/5 orthologs. Interestingly, HIV-1 Nef, murine leukemia virus (MLV) GlycoGag and equine infectious anemia virus (EIAV) S2 counteracted the antiviral activity of all SERINC3/5 orthologs with similar efficiency. Our results demonstrate that the antiviral activity of SERINC3/5 proteins is conserved in rodents and rabbits, and can be overcome by all three previously reported viral antagonists.

The evolution of S100A7 in primates: a model of concerted and birth-and-death evolution.

The human S100A7 resides in the epidermal differentiation complex (EDC) and has been described as a key effector of innate immunity. In humans, there are five S100A7 genes located in tandem-S100A7A, S100A7P1, S100AL2, S100A7, and S100AP2. The presence of several retroelements in the S100A7A/S100A7P1 and S100A7/S100A7P2 clusters suggests that these genes were originated from a duplication around ~ 35 million years ago, during or after the divergence of Platyrrhini and Catarrhini primates. To test this hypothesis, and taking advantage of the high number of genomic sequences available in the public databases, we retrieved S100A7 gene sequences of 12 primates belonging to the Cercopithecoidea and Hominoidea (Catarrhini species). Our results support the duplication theory, with at least one gene of each cluster being identified in both Cercopithecoidea and Hominoidea species. Moreover, given the presence of an ongoing gene conversion event between S100A7 and S100A7A, a high rate of mutation in S100A7L2 and the presence of pseudogenes, we proposed a model of concerted and birth-and-death evolution to explain the evolution of S100A7 gene family. Indeed, our results suggest that S100A7L2 most likely suffered a neofunctionalization in the Catarrhini group. Being S100A7 a major protein in innate defense, we believe that our findings could open new doors in the study of this gene family in immune system.

Strong selection of the TLR2 coding region among the Lagomorpha suggests an evolutionary history that differs from other mammals.

Toll-like receptors (TLRs) are one of the first lines of defense against pathogens and are crucial for triggering an appropriate immune response. Among TLRs, TLR2 is functional in all vertebrates and has high ability in detecting bacterial and viral pathogen ligands. The mammals' phylogenetic tree of TLR2 showed longer branches for the Lagomorpha clade, raising the hypothesis that lagomorphs experienced an acceleration of the mutation rate. This hypothesis was confirmed by (i) Tajima's test of neutrality that revealed different evolutionary rates between lagomorphs and the remaining mammals with lagomorphs presenting higher nucleotide diversity; (ii) genetic distances were similar among lagomorphs and between lagomorphs and other mammals; and (iii) branch models reinforced the existence of an acceleration of the mutation rate in lagomorphs. These results suggest that the lagomorph TLR2 has been strongly involved in pathogen recognition, which probably caused a host-pathogen arms race that led to the observed acceleration of the mutation rate.

Host-Specific Glycans Are Correlated with Susceptibility to Infection by Lagoviruses, but Not with Their Virulence.

Rabbit hemorrhagic disease virus (RHDV) and European brown hare syndrome virus (EBHSV) are two lagoviruses from the family Caliciviridae that cause fatal diseases in two leporid genera, Oryctolagus and Lepus, respectively. In the last few years, several examples of host jumps of lagoviruses among leporids were recorded. In addition, a new pathogenic genotype of RHDV emerged, and many nonpathogenic strains of lagoviruses have been described. The molecular mechanisms behind host shifts and the emergence of virulence are unknown. Since RHDV uses glycans of the histo-blood group antigen type as attachment factors to initiate infection, we studied if glycan specificities of the new pathogenic RHDV genotype, nonpathogenic lagoviruses, and EBHSV potentially play a role in determining the host range and virulence of lagoviruses. We observed binding to A, B, or H antigens of the histo-blood group family for all strains known to primarily infect European rabbits (Oryctolagus cuniculus), which have recently been classified as GI strains. However, we could not explain the emergence of virulence, since similar glycan specificities were found in several pathogenic and nonpathogenic strains. In contrast, EBHSV, recently classified as GII.1, bound to terminal ß-linked N-acetylglucosamine residues of O-glycans. Expression of these attachment factors in the upper respiratory and digestive tracts in three lagomorph species (Oryctolagus cuniculus, Lepuseuropaeus, and Sylvilagus floridanus) showed species-specific patterns regarding susceptibility to infection by these viruses, indicating that species-specific glycan expression is likely a major contributor to lagovirus host specificity and range.IMPORTANCE Lagoviruses constitute a genus of the family Caliciviridae comprising highly pathogenic viruses, RHDV and EBHSV, that infect rabbits and hares, respectively. Recently, nonpathogenic strains were discovered and new pathogenic strains have emerged. In addition, host jumps between lagomorphs have been observed. The mechanisms responsible for the emergence of pathogenicity and host species range are unknown. Previous studies showed that RHDV strains attach to glycans expressed in the upper respiratory and digestive tracts of rabbits, the likely portals of virus entry. Here, we studied the glycan-binding properties of novel pathogenic and nonpathogenic strains looking for a link between glycan binding and virulence or between glycan specificity and host range. We found that glycan binding did not correlate with virulence. However, expression of glycan motifs in the upper respiratory and digestive tracts of lagomorphs revealed species-specific patterns associated with the host ranges of the virus strains, suggesting that glycan diversity contributes to lagovirus host ranges.