FecalSeq enrichment with RAD Sequencing from non-invasive environmental samples holds promise for genetic monitoring of an imperiled lagomorph.

Despite advances in genomic sequencing and bioinformatics, conservation genomics is still often hindered by a reliance on non-invasive samples. The presence of exogenous DNA and the low quantity and poor quality of DNA in non-invasive samples have been a roadblock to sequencing, thereby limiting the potential for genomic monitoring of endangered species. Recent molecular advances, such as host DNA enrichment, hold promise for facilitating sequencing from non-invasive samples. We used the FecalSeq method to enrich DNA extracted from wild-collected fecal pellets of the imperiled New England cottontail and identified SNPs from 3RAD Sequencing. We obtained SNPs from rabbit pellets, including pellets that were collected in poor environmental conditions and samples that performed poorly with microsatellites. Measures of sequencing success improved with greater amounts of starting DNA and 32% of samples generated SNP genotypes that passed quality control filtering. Genotyping error rates were high, however, and the approach was unable to consistently distinguish unique individuals or matching genotypes, while it was suitable for recovering the expected population structure. Pairing FecalSeq enrichment with RADseq is a promising low-cost method for monitoring wild populations using non-invasive samples in an environmental context, but it may be better suited for informing conservation through population genomics.

Putative climate adaptation in American pikas (Ochotona princeps) is associated with copy number variation across environmental gradients.

Improved understanding of the genetic basis of adaptation to climate change is necessary for maintaining global biodiversity moving forward. Studies to date have largely focused on sequence variation, yet there is growing evidence that suggests that changes in genome structure may be an even more significant source of adaptive potential. The American pika (Ochotona princeps) is an alpine specialist that shows some evidence of adaptation to climate along elevational gradients, but previous work has been limited to single nucleotide polymorphism based analyses within a fraction of the species range. Here, we investigated the role of copy number variation underlying patterns of local adaptation in the American pika using genome-wide data previously collected across the entire species range. We identified 37-193 putative copy number variants (CNVs) associated with environmental variation (temperature, precipitation, solar radiation) within each of the six major American pika lineages, with patterns of divergence largely following elevational and latitudinal gradients. Genes associated (n = 158) with independent annotations across lineages, variables, and/or CNVs had functions related to mitochondrial structure/function, immune response, hypoxia, olfaction, and DNA repair. Some of these genes have been previously linked to putative high elevation and/or climate adaptation in other species, suggesting they may serve as important targets in future studies.

Phylogenomics of American pika (Ochotona princeps) lineage diversification.

Quaternary climate oscillations have profoundly influenced current species distributions. For many montane species, these fluctuations were a prominent driver in species range shifts, often resulting in intraspecific diversification, as has been the case for American pikas (Ochotona princeps). Range shifts and population declines in this thermally-sensitive lagomorph have been linked to historical and contemporary environmental changes across its western North American range, with previous research reconstructing five mitochondrial DNA lineages. Here, we paired genome-wide data (25,244 SNPs) with range-wide sampling to re-examine the number and distribution of intra-specific lineages, and investigate patterns of within- and among-lineage divergence and diversity. Our results provide genomic evidence of O. princeps monophyly, reconstructing six distinct lineages that underwent multiple rounds of divergence (0.809-2.81 mya), including a new Central Rocky Mountain lineage. We further found evidence for population differentiation across multiple spatial scales, and reconstructed levels of standing variation comparable to those found in other small mammals. Overall, our findings demonstrate the influence of past glacial cycles on O. princeps lineage diversification, suggest that current subspecific taxonomy may need to be revisited, and provide an important framework for investigations of American pika adaptive potential in the face of anthropogenic climate change.

Occurrence of Anaplasma spp. in wild lagomorphs from Southern Spain: Molecular detection of new Anaplasma bovis lineages.

Wild lagomorphs can act as reservoirs of several pathogens of public and animal health concern. However, the number of studies assessing the presence of Anaplasma spp. in these species is scarce. The aim of the present study was to molecularly identify Anaplasma spp. in wild rabbits (Oryctolagus cuniculus) and Iberian hares (Lepus granatensis) from Southern Spain and assess their epidemiological role in the maintenance of the bacterium. During 2017-2021, spleen samples of 394 wild rabbits and 145 Iberian hares were collected. Anaplasma DNA was detected using different PCR assays (16S rRNA and groEL) and phylogenetic analyses were carried out by Bayesian approach. The possible influence of lagomorph species, age and sex on the prevalence of Anaplasma spp. was evaluated by a multiple logistic regression model. The 9.4% of the rabbits were positive to Anaplasma bovis, but all the hares were negative. No significant differences were found in Anaplasma spp. prevalence regarding to age or sex. This is the first report of A. bovis in lagomorphs from Europe. The phylogenetic analysis of A. bovis confirms the existence of different clusters suggesting the existence of several lineages. In addition, a high divergence of nucleotide identity was observed within the lineage 4, which could result in the under-detection of some strains when using A. bovis-specific PCR, hindering its detection and characterization. Since this analysis is based on a limited number of nucleotide bases and sequences, more studies are needed for further characterize A. bovis, as well as its relationship with other Anaplasma spp.

Integrated Analysis of lncRNA-miRNA-mRNA Regulatory Network in Rapamycin-Induced Cardioprotection against Ischemia/Reperfusion Injury in Diabetic Rabbits.

The inhibition of mammalian target of rapamycin (mTOR) with rapamycin (RAPA) provides protection against myocardial ischemia/reperfusion (I/R) injury in diabetes. Since interactions between transcripts, including long non-coding RNA (lncRNA), microRNA(miRNA) and mRNA, regulate the pathophysiology of disease, we performed unbiased miRarray profiling in the heart of diabetic rabbits following I/R injury with/without RAPA treatment to identify differentially expressed (DE) miRNAs and their predicted targets of lncRNAs/mRNAs. Results showed that among the total of 806 unique miRNAs targets, 194 miRNAs were DE after I/R in diabetic rabbits. Specifically, eight miRNAs, including miR-199a-5p, miR-154-5p, miR-543-3p, miR-379-3p, miR-379-5p, miR-299-5p, miR-140-3p, and miR-497-5p, were upregulated and 10 miRNAs, including miR-1-3p, miR-1b, miR-29b-3p, miR-29c-3p, miR-30e-3p, miR-133c, miR-196c-3p, miR-322-5p, miR-499-5p, and miR-672-5p, were significantly downregulated after I/R injury. Interestingly, RAPA treatment significantly reversed these changes in miRNAs. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis indicated the participation of miRNAs in the regulation of several signaling pathways related to I/R injury, including MAPK signaling and apoptosis. Furthermore, in diabetic hearts, the expression of lncRNAs, HOTAIR, and GAS5 were induced after I/R injury, but RAPA suppressed these lncRNAs. In contrast, MALAT1 was significantly reduced following I/R injury, with the increased expression of miR-199a-5p and suppression of its target, the anti-apoptotic protein Bcl-2. RAPA recovered MALAT1 expression with its sponging effect on miR-199-5p and restoration of Bcl-2 expression. The identification of novel targets from the transcriptome analysis in RAPA-treated diabetic hearts could potentially lead to the development of new therapeutic strategies for diabetic patients with myocardial infarction.

Life on the edge-a changing genetic landscape within an iconic American pika metapopulation over the last half century.

Declines and extirpations of American pika (Ochotona princeps) populations at historically occupied sites started being documented in the literature during the early 2000s. Commensurate with global climate change, many of these losses at peripheral and lower elevation sites have been associated with changes in ambient air temperature and precipitation regimes. Here, we report on a decline in available genetic resources for an iconic American pika metapopulation, located at the southwestern edge of the species distribution in the Bodie Hills of eastern California, USA. Composed of highly fragmented habitat created by hard rock mining, the ore dumps at this site were likely colonized by pikas around the end of the 19th century from nearby natural talus outcrops. Genetic data extracted from both contemporary samples and archived natural history collections allowed us to track population and patch-level genetic diversity for Bodie pikas across three distinct sampling points during the last half- century (1948-1949, 1988-1991, 2013-2015). Reductions in within-population allelic diversity and expected heterozygosity were observed across the full time period. More extensive sampling of extant patches during the 1988-1991 and 2013-2015 periods revealed an increase in population structure and a reduction in effective population size. Furthermore, census records from the last 51 years as well as archived museum samples collected in 1947 from a nearby pika population in the Wassuk range (Nevada, USA) provide further support of the increasing isolation and genetic coalescence occurring in this region. This study highlights the importance of museum samples and long-term monitoring in contextualizing our understanding of population viability.

Disrupted circadian rhythms in the plateau pika.

The plateau pika (Ochotona curzoniae) is the most populous mammal on the 'third pole', the Qinghai-Tibet Plateau, and is presumed to have inhabited the region before the plateau rose up from sea level. Herein we discuss the disrupted circadian rhythm in the plateau pika and the gene polymorphism behind this phenotype, placing these findings in the broader context of circadian rhythms under extreme conditions.

Species-Level Taxonomic Characterization of Uncultured Core Gut Microbiota of Plateau Pika.

Rarely has the vast diversity of bacteria on Earth been profiled, particularly on inaccessible plateaus. These uncultured microbes, which are also known as "microbial dark matter," may play crucial roles in maintaining the ecosystem and are linked to human health, regarding pathogenicity and prebioticity. The plateau pika (Ochotona curzoniae) is a small burrowing steppe lagomorph that is endemic to the Qinghai-Tibetan Plateau and is a keystone species in the maintenance of ecological balance. We used a combination of full-length 16S rRNA amplicon sequencing, shotgun metagenomics, and metabolomics to elucidate the species-level community structure and the metabolic potential of the gut microbiota of the plateau pika. Using a full-length 16S rRNA metataxonomic approach, we clustered 618 (166 ± 35 per sample) operational phylogenetic units (OPUs) from 105 plateau pika samples and assigned them to 215 known species, 226 potentially new species, and 177 higher hierarchical taxa. Notably, 39 abundant OPUs (over 60% total relative abundance) are found in over 90% of the samples, thereby representing a "core microbiota." They are all classified as novel microbial lineages, from the class to the species level. Using metagenomic reads, we independently assembled and binned 109 high-quality, species-level genome bins (SGBs). Then, a precise taxonomic assignment was performed to clarify the phylogenetic consistency of the SGBs and the 16S rRNA amplicons. Thus, the majority of the core microbes possess their genomes. SGBs belonging to the genus Treponema, the families Muribaculaceae, Lachnospiraceae, and Oscillospiraceae, and the order Eubacteriales are abundant in the metagenomic samples. In addition, multiple CAZymes are detected in these SGBs, indicating their efficient utilization of plant biomass. As the most widely connected metabolite with the core microbiota, tryptophan may relate to host environmental adaptation. Our investigation allows for a greater comprehension of the composition and functional capacity of the gut microbiota of the plateau pika. IMPORTANCE The great majority of microbial species remain uncultured, severely limiting their taxonomic characterization and biological understanding. The plateau pika (Ochotona curzoniae) is a small burrowing steppe lagomorph that is endemic to the Qinghai-Tibetan Plateau and is considered to be the keystone species in the maintenance of ecological stability. We comprehensively investigated the gut microbiota of the plateau pika via a multiomics endeavor. Combining full-length 16S rRNA metataxonomics, shotgun metagenomics, and metabolomics, we elucidated the species-level taxonomic assignment of the core uncultured intestinal microbiota of the plateau pika and revealed their correlation to host nutritional metabolism and adaptation. Our findings provide insights into the microbial diversity and biological significance of alpine animals.

Ancient DNA of the Don-Hares Assumes the Existence of Two Distinct Mitochondrial Clades in Northeast Asia.

Paleoclimatic changes during the Pleistocene-Holocene transition is suggested as a main factor that led to species extinction, including the woolly mammoth (Mammuthus primigenius), Steller's sea cow (Hydrodamalis gigas) and the Don-hare (Lepus tanaiticus). These species inhabited the territory of Eurasia during the Holocene, but eventually went extinct. The Don-hare is an extinct species of the genus Lepus (Leporidae, Lagomorpha), which lived in the Late Pleistocene-Early Holocene in Eastern Europe and Northern Asia. For a long time, the Don-hare was considered a separate species, but at the same time, its species status was disputed, taking into account both morphological data and mitochondrial DNA. In this study, mitochondrial genomes of five Don-hares, whose remains were found on the territory of Northeastern Eurasia were reconstructed. Firstly, we confirm the phylogenetic proximity of the "young" specimens of Don-hare and mountain or white hare, and secondly, that samples older than 39 Kya form a completely distinct mitochondrial clade.

Genomic Consequences of and Demographic Response to Pervasive Hybridization Over Time in Climate-Sensitive Pikas.

Rare and geographically restricted species may be vulnerable to genetic effects from inbreeding depression in small populations or from genetic swamping through hybridization with common species, but a third possibility is that selective gene flow can restore fitness (genetic rescue). Climate-sensitive pikas (Ochotona spp.) of the Qinghai-Tibetan Plateau (QHTP) and its vicinity have been reduced to residual populations through the movement of climatic zones during the Pleistocene and recent anthropogenic disturbance, whereas the plateau pika (O. curzoniae) remains common. Population-level whole-genome sequencing (n = 142) of six closely related species in the subgenus Ochotona revealed several phases of ancient introgression, lineage replacement, and bidirectional introgression. The strength of gene flow was the greatest from the dominant O. curzoniae to ecologically distinct species in areas peripheral to the QHTP. Genetic analyses were consistent with environmental reconstructions of past population movements. Recurrent periods of introgression throughout the Pleistocene revealed an increase in genetic variation at first but subsequent loss of genetic variation in later phases. Enhanced dispersion of introgressed genomic regions apparently contributed to demographic recovery in three peripheral species that underwent range shifts following climate oscillations on the QHTP, although it failed to drive recovery of northeastern O. dauurica and geographically isolated O. sikimaria. Our findings highlight differences in timescale and environmental background to determine the consequence of hybridization and the unique role of the QHTP in conserving key evolutionary processes of sky island species.

Molecular phenotyping uncovers differences in basic housekeeping functions among closely related species of hares (Lepus spp., Lagomorpha: Leporidae).

Speciation is a fundamental evolutionary process, which results in genetic differentiation of populations and manifests as discrete morphological, physiological and behavioural differences. Each species has travelled its own evolutionary trajectory, influenced by random drift and driven by various types of natural selection, making the association of genetic differences between the species with the phenotypic differences extremely complex to dissect. In the present study, we have used an in vitro model to analyse in depth the genetic and gene regulation differences between fibroblasts of two closely related mammals, the arctic/subarctic mountain hare (Lepus timidus Linnaeus) and the temperate steppe-climate adapted brown hare (Lepus europaeus Pallas). We discovered the existence of a species-specific expression pattern of 1623 genes, manifesting in differences in cell growth, cell cycle control, respiration, and metabolism. Interspecific differences in the housekeeping functions of fibroblast cells suggest that speciation acts on fundamental cellular processes, even in these two interfertile species. Our results help to understand the molecular constituents of a species difference on a cellular level, which could contribute to the maintenance of the species boundary.

A homotetrameric hemoglobin expressed in alveolar epithelial cells increases blood oxygenation in high-altitude plateau pika (Ochotona curzoniae).

The plateau pika (Ochotona curzoniae) is native to the Qinghai-Tibet Plateau. In this study, the gene that encodes a heme-binding protein in the pulmonary surfactant (PS) of the pika is identified. The protein is a homotetrameric hemoglobin (δ4) encoded by HBD (δ). HBD is expressed in alveolar epithelial type II (ATII) and type I (ATI) cells, upregulated by hypoxia. δ4 is secreted into alveolar cavities through osmiophilic multilamellar bodies. HBD expression is downregulated by RNAi, which significantly increases hypoxia-inducible factor 1α expression in lung tissue and red blood cells and hemoglobin and blood lactate concentrations but significantly decreases arterial partial pressure of oxygen (PaO2). Our results indicate that plateau pikas physiologically show hypoxemia when HBD expression is downregulated. Therefore, specific HBD expression in the lungs helps plateau pikas to obtain oxygen under hypoxia by maintaining higher PaO2. These findings provide insights into the adaptive mechanisms of plateau pikas to withstand high-altitude environments.

Comparative genomics reveals putative evidence for high-elevation adaptation in the American pika (Ochotona princeps).

High-elevation environments have lower atmospheric oxygen content, reduced temperatures, and higher levels of UV radiation than found at lower elevations. As such, species living at high elevations must overcome these challenges to survive, grow, and reproduce. American pikas (Ochotona princeps) are alpine lagomorphs that are habitat specialists typically found at elevations >2,000 m. Previous research has shown putative evidence for high-elevation adaptation; however, investigations to date have been limited to a fraction of the genome. Here, we took a comparative genomics approach to identify putative regions under selection using a chromosomal reference genome assembly for the American pika relative to 8 other mammalian species targeted based on phylogenetic relatedness and (dis)similarity in ecology. We first identified orthologous gene groups across species and then extracted groups containing only American pika genes as well as unclustered pika genes to inform functional enrichment analyses; among these, we found 141 enriched terms with many related to hypoxia, metabolism, mitochondrial function/development, and DNA repair. We identified 15 significantly expanded gene families within the American pika across all orthologous gene groups that displayed functionally enriched terms associated with hypoxia adaptation. We further detected 196 positively selected genes, 41 of which have been associated with putative adaptation to hypoxia, cold tolerance, and response to UV following a literature review. In particular, OXNAD1, NRDC, and those genes critical in DNA repair represent important targets for future research to examine their functional implications in the American pika, especially as they may relate to adaptation to rapidly changing environments.

Genomics and morphometrics reveal the adaptive evolution of pikas.

Pikas (Lagomorpha: Ochotonidae) are small mouse-like lagomorphs. To investigate their adaptation to different ecological environments during their dispersal from the Qinghai-Xizang (Tibet) Plateau (QTP), we collected 226 pikas and measured 20 morphological characteristics and recorded habitat information. We also sequenced the genome of 81 specimens, representing 27 putative pika species. The genome-wide tree based on 4 090 coding genes identified five subgenera, i.e., Alienauroa, Conothoa, Lagotona, Ochotona, and Pika, consistent with morphometric data. Morphologically, Alienauroa and Ochotona had similar traits, including smaller size and earlier divergence time compared to other pikas. Consistently, the habitats of Alienauroa and Ochotona differed from those of the remaining subgenera. Phylogenetic signal analysis detected 83 genes significantly related to morphological characteristics, including several visual and hearing-related genes. Analysis of shared amino acid substitutions and positively selected genes (PSGs) in Alienauroa and Ochotona identified two genes, i.e., mitochondrial function-related TSFM (p.Q155E) and low-light visual sensitivity-related PROM1 (p.H419Y). Functional experiments demonstrated that TSFM-155E significantly enhanced mitochondrial function compared to TSFM-155Q in other pikas, and PROM1-419Y decreased the modeling of dynamic intracellular chloride efflux upon calcium uptake. Alienauroa and Ochotona individuals mostly inhabit different environments (e.g., subtropical forests) than other pikas, suggesting that a shift from the larger ancestral type and changes in sensory acuity and energy enhancement may have been required in their new environments. This study increases our understanding of the evolutionary history of pikas.

A highland-adaptation mutation of the Epas1 protein increases its stability and disrupts the circadian clock in the plateau pika.

The Qinghai-Tibet Plateau (QTP) harbors hundreds of species well adapted to its extreme conditions, including its low-oxygen (hypoxic) atmosphere. Here, we show that the plateau pika-a keystone mammal of the QTP-lacks robust circadian rhythms. The major form of the plateau pika Epas1 protein includes a 24-residue insert caused by a point mutation at the 5' juncture site of Intron14 and is more stable than other mammalian orthologs. Biochemical studies reveal that an Epas1-Bmal1 complex with lower trans-activation activity occupies the E1/E2 motifs at the promoter of the core-clock gene Per2, thus explaining how an Epas1 mutation-selected in the hypoxic conditions of the QTP-disrupts the molecular clockwork. Importantly, experiments with hypoxic chambers show that mice expressing the plateau pika Epas1 ortholog in their suprachiasmatic nucleus have dysregulated central clocks, and pika Epas1 knockin mice reared in hypoxic conditions exhibit dramatically reduced heart damage compared with wild-type animals.

Seasonal variations in the composition and functional profiles of gut microbiota reflect dietary changes in plateau pikas.

Seasonal variations in gut microbiota of small mammals and how they are influenced by environmental variables are relatively poorly understood. We sampled 162 wild plateau pikas (Ochotona curzoniae) in 4 seasons over 2 and a half years and recorded the air temperature, precipitation, and nutrient content in edible vegetation at the sampling site. After conducting 16S rRNA and shotgun metagenomic sequencing, we found that the highest alpha diversity, the relative abundance of Firmicutes, and the simplest co-occurrence network occurred in winter, whereas the highest relative abundance of Proteobacteria and the most complex network structure were observed in spring. The highest relative abundance of Verrucomicrobiota and Spirochaetota was seen in summer and autumn, respectively. Air temperature, precipitation, and the contents of crude protein, crude fiber, and polysaccharide in vegetation had significant effects on the seasonal changes in gut microbiota. Diet contributed more to microbial variation than climatic factors. Metagenomic analysis revealed that the amino acid metabolism pathway and axillary activity enzymes were most abundant in summer, while abundance of carbohydrate-binding modules and carbohydrate esterases was highest in spring. These microbial variations were related to the changes in dietary nutrition, indicating that gut microbiota of plateau pika contribute to the efficient use of food resources. This study provides new evidence of how external environmental factors affect the intestinal environment of small mammals.

Lung transcriptome analysis for the identification of genes involved in the hypoxic adaptation of plateau pika (Ochotona curzoniae).

The plateau pika, a typical hypoxia-tolerant mammal lives 3000-5000 m above sea level on the Qinghai-Tibet Plateau, has acquired many physiological and morphological characteristics and strategies in its adaptation to sustained, high-altitude hypoxia. Blunted hypoxic pulmonary vasoconstriction is one such strategy, but the genes involved in this strategy have not been elucidated. Here, we investigated the genes involved and their expression profiles in the lung transcriptome of plateau pikas subjected to different hypoxic conditions (using low-pressure oxygen cabins). A slight, right ventricular hypertrophy was observed in pikas of the control group (altitude: 3200 m) vs. those exposed to 5000 m altitude conditions for one week. Our assembly identified 67,774 genes; compared with their expression in the control animals, 866 and 8364 genes were co-upregulated and co-downregulated, respectively, in pikas subjected to 5000 m altitude conditions for 1 and 4 w. We elucidated pathways that were associated with pulmonary vascular arterial pressure, including vascular smooth muscle contraction, HIF-1 signalling, calcium signalling, cGMP-PKG signalling, and PI3K-Akt signalling based on the differentially expressed genes; the top-100 pathway enrichments were found between the control group and the group exposed to 5000 m altitude conditions for 4 w. The mRNA levels of 18 candidate gene showed that more than 83% of genes were expressed and the number of transcriptome The up-regulated genes were EPAS1, Hbα, iNOS, CX40, CD31, PPM1B, HIF-1α, MYLK, Pcdh12, Surfactant protein B, the down-regulated genes were RYR2, vWF, RASA1, CLASRP, HIF-3α. Our transcriptome data are a valuable resource for future genomic studies on plateau pika.

The complete mitochondrial genome of the 'Zacatuche' Volcano rabbit (Romerolagus diazi), an endemic and endangered species from the Volcanic Belt of Central Mexico.

BACKGROUND: The 'Zacatuche', 'Teporingo', or Volcano rabbit (Romerolagus diazi) belongs to the family Leporidae, is an endemic species restricted to the Central part of the Trans-Mexican Volcanic Belt, and is considered 'endangered' by the IUCN Red List of Threatened Species. METHODS AND RESULTS: This study reports, for the first time, the complete mitochondrial genome of R. diazi and examined the phylogenetic position of R. diazi among other closely related co-familiar species using mitochondrial protein-coding genes (PCGs). The mitogenome of R. diazi was assembled from short Illumina 150 bp pair-end reads with a coverage of 189x. The AT-rich mitochondrial genome of R. diazi is 17,400 bp in length and is comprised of 13 PCGs, two ribosomal RNA genes, and 22 transfer RNA genes. The gene order observed in the mitochondrial genome of R. diazi is identical to that reported for other leporids. Phylogenetic analyses based on PCGs support the basal position of Romerolagus within the Leporidae, at least when compared to the genera Oryctolagus and Lepus. Nonetheless, additional mitochondrial genomes from species belonging to the genera Bunolagus, Sylvilagus, and Pronolagus, among others, are needed before a more robust conclusion about the derived vs basal placement of Romerolagus within the family Leporidae can be reached based on mitochondrial PCGs. CONCLUSIONS: This is the first genomic resource developed for R. diazi and it represents a tool to improve our understanding about the ecology and evolutionary biology of this iconic and endangered species.

Molecular identification of zoonotic Rickettsia species in Ixodidae parasitizing wild lagomorphs from Mediterranean ecosystems.

A survey study was carried out to identify tick species parasitizing wild lagomorphs in Mediterranean ecosystems in southern Spain and to determine the occurrence of Rickettsia species present in these ticks in this region. A total of 1304 European wild rabbits (Oryctolagus cuniculus) and 58 Iberian hares (Lepus granatensis) were individually examined for the presence of ticks. Ticks were found in 42.9% and 50% of the wild rabbits and hares sampled, respectively. A total of 1122 ticks were collected and five species, including Rhipicephalus pusillus, Hyalomma lusitanicum, Haemaphysalis hispanica, Ixodes ventalloi and Rhipicephalus sanguineus sensu lato (s.l.), were microscopically and molecularly identified at the 16S rRNA gene. This is the first study on Ixodidae parasitizing Iberian hares. The presence of Rickettsia DNA was assessed in 254 tick pools (according to hunting states, lagomorph species, tick species and tick development stage) using PCR assays targeting the rOmpA, rOmpB and gltA. Twenty-seven pools (10.6%) were positive to Rickettsia DNA. Five zoonotic Rickettsia species were identified, being Rickettsia massiliae the most frequent (4.7%), followed by Rickettsia sibirica subsp. mongolitimonae (2.8%), Rickettsia slovaca (2.0%), Rickettsia aeschlimannii (0.8%) and Rickettsia africae (0.4%). The results suggest that wild rabbits and Iberian hares are parasitized by a wide range of tick species and that these lagomorphs may play an important role in the sylvatic cycle of some zoonotic Rickettsia species in Mediterranean ecosystems. Our data represent the first report of R. massiliae, R. aeschlimannii, R. slovaca and R africae in ticks collected in wild lagomorphs in Europe, and the first report of not imported R. africae in this continent. Since R. slovaca and R. africae DNA was detected in tick species different to their main vectors, further studies are warranted to unravel the role of wild lagomorphs in the epidemiology of these vector-borne pathogens.

Comparative susceptibility of eastern cottontails and New Zealand white rabbits to classical rabbit haemorrhagic disease virus (RHDV) and RHDV2.

Rabbit haemorrhagic disease virus (RHDV) is associated with high morbidity and mortality in the European rabbit (Oryctolagus cuniculus). In 2010, a genetically distinct RHDV named RHDV2 emerged in Europe and spread to many other regions, including North America in 2016. Prior to this study it was unknown if eastern cottontails (ECT(s); Sylvilagus floridanus), one of the most common wild lagomorphs in the United States, were susceptible to RHDV2. In this study, 10 wild-caught ECTs and 10 New Zealand white rabbits (NZWR(s); O. cuniculus) were each inoculated orally with either RHDV (RHDVa/GI.1a; n = 5 per species) or RHDV2 (a recombinant GI.1bP-GI.2; n = 5 per species) and monitored for the development of disease. Three of the five ECTs that were infected with RHDV2 developed disease consistent with RHD and died at 4 and 6 days post-inoculation (DPI). The RHDV major capsid protein/antigen (VP60) was detected in the livers of three ECTs infected with RHDV2, but none was detected in the ECTs infected with RHDV. Additionally, RHD viral RNA was detected in the liver, spleen, intestine and blood of ECTs infected with RHDV2, but not in the ECTs infected with RHDV. RHD viral RNA was detected in urine, oral swabs and rectal swabs in at least two of five ECTs infected with RHDV2. One ECT inoculated with RHDV2 seroconverted and developed a high antibody titre by the end of the experimental period (21 DPI). ECTs inoculated with the classic RHDV did not seroconvert. In comparison, NZWRs inoculated with RHDV2 exhibited high mortality (five of five) at 2 DPI and four of five NZWRs inoculated with RHDV either died or were euthanized at 2 DPI indicating both of these viruses were highly pathogenic to this species. This experiment indicates that ECTs are susceptible to RHDV2 and can shed viral RNA, thereby suggesting this species could be involved in the epidemiology of this virus.