Evidence for a survival-driven traveling wave in a keystone boreal predator population.

Cyclical population dynamics are a common phenomenon in populations worldwide, yet the spatial organization of these cycles remains poorly understood. In this study, we investigated the spatial form and timing of a population collapse from 2018 to 2022 in Canada lynx (Lynx canadensis) across the northwest boreal forest. We analyzed survival, reproduction, and dispersal data from 143 individual global positioning system (GPS) collared lynx from populations across five study sites spanning interior Alaska to determine whether lynx displayed characteristics of a population wave following a concurrent wave in snowshoe hare (Lepus americanus) abundance. Reproductive rates declined across the study sites; however, site-level reproduction declined first in our easternmost study sites, supporting the idea of a population wave. Despite a clear increase in percent of dispersing lynx, there was no evidence of directional bias in dispersal following a hare population wave. Analysis did show increasingly poor survival for lynx dispersing to the east compared to combined resident and westward dispersal. This pattern is consistent with a survival-mediated population wave in lynx as the driver of the theorized population wave. The combination of these factors supports the idea of a hierarchical response to snowshoe hare population declines with a drop in lynx reproduction followed by increased dispersal, and finally reduced survival. All of this evidence is consistent with the expected characteristics of a population undergoing a traveling wave and supports the hypothesis that lynx presence may facilitate and mirror the underlying wave patterns in snowshoe hare.

HAT: de novo variant calling for highly accurate short-read and long-read sequencing data.

MOTIVATION: de novo variants (DNVs) are variants that are present in offspring but not in their parents. DNVs are both important for examining mutation rates as well as in the identification of disease-related variation. While efforts have been made to call DNVs, calling of DNVs is still challenging from parent-child sequenced trio data. We developed Hare And Tortoise (HAT) as an automated DNV detection workflow for highly accurate short-read and long-read sequencing data. Reliable detection of DNVs is important for human genomics and HAT addresses this need. RESULTS: HAT is a computational workflow that begins with aligned read data (i.e. CRAM or BAM) from a parent-child sequenced trio and outputs DNVs. HAT detects high-quality DNVs from Illumina short-read whole-exome sequencing, Illumina short-read whole-genome sequencing, and highly accurate PacBio HiFi long-read whole-genome sequencing data. The quality of these DNVs is high based on a series of quality metrics including number of DNVs per individual, percent of DNVs at CpG sites, and percent of DNVs phased to the paternal chromosome of origin. AVAILABILITY AND IMPLEMENTATION: https://github.com/TNTurnerLab/HAT.

Effect of Myxoma Virus Species Jump on Iberian Hare Populations.

The myxoma virus species jump from European rabbits (Oryctolagus cuniculus) to Iberian hares (Lepus granatensis) has raised concerns. We assess the decline suffered by Iberian hare populations on the Iberian Peninsula and discuss the association between the effect of myxomatosis and the average abundance index, which we estimated by using hunting bags.

Patterns of evolution in MHC class II DQA and DQB exon 2 genes of Alpine mountain hares, Lepus timidus varronis, and sympatric and parapatric brown hares, L. europaeus, from Switzerland.

In natural populations, hybridization is known to occur between a wide range of species. However, its evolutionary significance is less clear. Genes involved in fighting pathogens are considered excellent candidates for studying adaptive introgression, although both introgression and balancing selection can generate similar patterns of diversity and differentiation. Here, we compared DQA and DQB MHC class II and microsatellite allelic diversity of sympatric and parapatric mountain (Lepus timidus) and brown hare (L. europaeus) populations from Switzerland. We detected higher genetic diversity in brown hares compared to mountain hares at both MHC and microsatellite loci. We consider the observed patterns of microsatellite diversity both for L. europaeus and L. timidus as result of stochastic demographic processes while the pattern of MHC polymorphism of the studied hare populations can be explained by pathogen-driven selection. Rare bidirectional gene flow between both hare species seems to occur specifically for MHC alleles. However, the high number of shared alleles showing similar high frequency in both species suggests that reciprocally exchanged MHC alleles are being maintained via balancing selection. Adaptation to similar pathogen communities can also lead to parallel selection of MHC alleles. Positive selection, recombination and mutations have played different roles in shaping the patterns of MHC allelic diversity in and differentiation between both species. Results for the latter evolutionary forces do not show a better matching between the sympatric populations compared to the parapatric ones, suggesting a minor role of introgression for the observed evolutionary patterns of the studied hare species.

Vegetation biomass and topography are associated with seasonal habitat selection and fall translocation behavior in Arctic hares.

Habitat selection theory suggests that environmental features selected at coarse scales reveal fundamental factors affecting animal fitness. When these factors vary across seasons, they may lead to large-scale movements, including long-distance seasonal migrations. We analyzed the seasonal habitat selection of 25 satellite-tracked Arctic hares from a population on Ellesmere Island (Nunavut, Canada) that relocated over 100 km in the fall. Since no other lagomorph is known to perform such extensive movements, this population offered an ideal setting to test animal movement and habitat selection theory. On summer grounds hares selected low elevation areas, while on winter grounds they selected high vegetation biomass, high elevation, and steep slopes. During fall relocation, they alternated between stopover and traveling behavioral states (ratio 2:1). Stopover locations were characterized by higher vegetation heterogeneity and lower rugosity than traveling locations, while vegetation biomass and elevation interacted to explain stopover locations in a more complex way. The selected combination of environmental features thus varied across seasons and behavioral states, in a way broadly consistent with predictions based on the changing food and safety needs of hares. Although causality was not demonstrated, our results improve our understanding of long-distance movements and habitat selection in Arctic hares, as well as herbivore ecology in the polar desert. Results also provide strong support to animal movement and habitat selection theory, by showing how some important hypotheses hold when tested in a species phylogenetically distinct from most animal models used in this research field.

Typing of Leishmania isolates from vectors and leporids of the Madrid (Spain) outbreak.

In 2009, a large outbreak of leishmaniasis, associated with environmental changes, was declared near Madrid (Spain), in which Phlebotomus perniciosus was the vector, whereas the main reservoirs were hares and rabbits. Analysis of isolates from humans, vectors and leporids from the focus identified the Leishmania infantum ITS-Lombardi genotype. However, multilocus enzyme electrophoresis (MLEE), the reference technique for Leishmania typing, and sequencing of the hsp70 gene, a commonly used marker, were not performed. In the present study, 19 isolates from P. perniciosus (n = 11), hares (n = 5) and rabbits (n = 3) from the outbreak area, all characterized as ITS-Lombardi in previous studies, were analysed by MLEE and hsp70 sequencing. The hsp70 results confirmed that all the analysed strains are L. infantum. However, by MLEE, 4 different zymodemes of L. infantum were identified based on variable mobilities of the NP1 enzyme: MON-34 (NP1100, n = 11), MON-80 (NP1130, n = 6), MON-24 (NP1140, n = 1) and MON-331 (NP1150, n = 1). The relative frequency of these zymodemes does not correspond to their usual occurrence in Spain. Moreover, MON-34 and MON-80 were found in P. perniciosus, hares and rabbits for the first time. These findings continue to provide insights into the outbreak and call for further studies with a higher number of strains.

Estrogen regulated transcription of the non-estrogen-regulated hamster uteroglobin gene in MCF-7 cells.

To study the estrogen regulated transcription of the uteroglobin (UG) gene, the founding member of the secretoglobin family widely expressed in many different mammalian species, we re-created functional estrogen response elements (EREs) in the UG gene promoter from a species where UG expression is not regulated by estrogens: the hamster Mesocricetus auratus (Ma), to ascertain if the lack of functional EREs is the real cause of its estrogen insensitivity. Functional EREs in the hamster promoter, including the consensus ERE (cERE), failed to respond to an appropriate estrogen stimulus compared with its estrogen regulated ortholog from the brown hare Lepus capensis (Lc). As the nucleotide sequence is the only difference between genetic constructs from these two species, we suspected that the UG promoter from the hamster probably contains cis-acting genetic elements that negatively impairs the estrogen-regulated transcription mediated by the functional ERE. Accordingly, we prepared chimeric DNA constructs which eventually allowed to identify a region located 29 base pairs (bp) downstream of the ERE as responsible for the lack of estrogen-responsiveness of the Ma-UG gene in the breast cancer cell line MCF-7. This region contains the sequence ACACCCC which has been identified as the core sequence of the Sp/ Krüppel-like factor (KLF) family of transcription factors. This finding is relevant, not only due to the observation on a novel mechanism that control estrogen-induced transcription, but also because it may encourage further investigation for better defining specific genes with an ERE that do not respond to estrogen signaling in MCF-7 cells, a cell line widely employed as an in vitro model in breast cancer research.

Selective sweep analysis of the adaptability of the Yarkand hare (Lepus yarkandensis) to hot arid environments using SLAF-seq.

The Yarkand hare (Lepus yarkandensis) inhabits arid desert areas and is endemic to China. It has evolved various adaptations to survive in hot arid environments, including stress responses, the ability to maintain water homeostasis and heat tolerance. Here, we performed a selective sweep analysis to identify the candidate genes for adaptation to hot arid environments in the Yarkand hare. A total of 397 237 single-nucleotide polymorphisms were obtained from 80 Yarkand hares, which inhabit hot arid environments, and 36 Tolai hares (Lepus tolai), which inhabit environments with a mild climate, via specific-locus amplified fragment sequencing. We identified several candidate genes that were associated with the heat stress response (HSPE1), oxidative stress response (SLC23A and GLRX2), immune response (IL1R1 and IRG1), central nervous system development (FGF13, THOC2, FMR1 and MECP2) and regulation of water homeostasis (CDK1) according to fixation index values and θπ ratios in the selective sweep analysis, and six of these genes (GLRX2, IRG1, FGF13, FMR1, MECP2 and CDK1) are newly discovered genes. To the best of our knowledge, this is the first study to identify candidate genes for adaptation to hot arid environments in the Yarkand hare. The results of this study enhance our understanding of the adaptation of the Yarkand hare to hot arid environments and will aid future studies aiming to functionally verify these candidate genes.

Toxoplasma gondii infection in European brown hares (Lepus europaeus) in Brazil.

Toxoplasmosis is a worldwide zoonosis that affects warm-blooded animals, including humans. Wild animals can act as intermediate hosts of this pathogen; thus, this study aims to detect Toxoplasma gondii infection in invasive European brown hares in Brazil. For this, 72 wild European brown hares were captured from July 2020 to June 2022 in three Brazilian states: São Paulo, Paraná, and Rio Grande do Sul. The diagnostic of Toxoplasma gondii infection was performed by bioassay in mouse, histopathology in Hematoxylin-Eosin-stained tissue sections (brain, liver, lungs, kidneys, and small intestine), serology by IFAT, and molecular techniques by conventional PCR and qPCR. The combined prevalence of the different diagnostic methods was 51.4% (37/72, CI= 40.1 - 62.6 %), and there was no statistical difference between sexes, age range, or geographical region of the hosts. Mouse bioassay was the technique that detected more positive hares. To our knowledge, this is the first confirmation of Toxoplasma gondii infection in invasive European brown hares in Brazil. These animals act as reservoirs and potential infection source for carnivores and other wild and domestic animals, including humans, thus contributing to perpetuate the disease cycle in São Paulo, Paraná, and Rio Grande do Sul States. Research such as the present study is necessary to raise awareness about the role of animals in the disease cycle.

The Risk, Need, and Responsivity Relevance of Working Alliance in a Sexual Offense Treatment Program: Its Intersection With Psychopathy, Diversity, and Treatment Change.

The current study examined the self-reported working alliance of men attending a high intensity sexual offense treatment program and its associations with psychopathy, sexual violence risk, treatment change, and recidivism, in a Canadian sample of 317 incarcerated men followed up an average of approximately 10 years post release. Working Alliance Inventory (WAI; Horvath & Greenberg, 1989) self-reported total, Task, Bond, and Goal scores were positively correlated with treatment related changes in risk, and inversely associated with Psychopathy Checklist-Revised (PCL-R; Hare, 1991; Wang & Hare, 2003) scores. The Affective facet of the PCL-R, representing the callous-unemotional features of the syndrome, uniquely predicted lower Bond and Goal scores controlling for the other facets. Cox regression survival analyses demonstrated that sexual violence risk predicted increased sexual recidivism while change predicted decreased sexual recidivism controlling for PCL-R total score; however, WAI scores (particularly the Goal component) were also unexpectedly associated with increased sexual recidivism. For violent recidivism, psychopathy, risk, and change incremented the prediction of general violence, while the WAI was not significantly associated with this outcome. A set of parallel analyses, stratified by Indigenous ethnocultural heritage, demonstrated some continuity, but also potential areas of difference, in substantive findings. Risk, need, responsivity implications of the working alliance for the treatment of high psychopathy correctional clientele, and how this may intersect with Indigenous heritage, are discussed.

Potentially toxic metals in seawater, sediment and seaweeds: bioaccumulation, ecological and human health risk assessment.

This study assesses the bioaccumulation, ecological, and health risks associated with potentially toxic metals (PTMs), including Pb, Hg, Cd, As, and Cr in Hare Island, Thoothukudi. The results revealed that the concentration of PTMs in sediment, seawater, and S. wightii ranged from 0.095 to 2.81 mg kg-1, 0.017 to 1.515 mg L-1, and 0.076 to 5.713 mg kg-1, respectively. The highest concentrations of PTMs were found in the S. wightii compared to seawater and sediment. The high bioaccumulation of Hg and As in S. wightii suggests that it can be used as a bioindicator for these elements in this region. The ecological risk indices, which include individual, complex, biological, and ecological pollution indices, suggest that Hare Island had moderate contamination with Hg and Cd. However, there are no human health risks associated with PTMs. This study examines the current ecological and health risks associated with PTMs and emphasizes the importance of regular monitoring.

Concentration of selected biogenic and risk elements in liver, kidneys and muscle of domestic rabbit and wild brown hare.

In the present study the concentration of selected elements in tissues of domestic rabbits and of wild brown-hares (kidneys, liver, and muscle - m. quadriceps femoris) in Slovakian habitats were determined. After mineralization the elements examined were detected using flame atomic absorption spectrophotometry/graphite furnace atomic absorption spectrophotometry. For rabbits, Fe in the liver was correlated with essential (Mn, Cu) (R2 = 0.94, p < 0.05; R2 = 0.96, p < 0.05 respectively) or toxic (Pb) elements (R2 = -0.93, p < 0.05). For hares, significant correlations were found between Cd and Cu or between Cd and Mn in the kidneys (R2 = -0.96, p < 0.05; R2 = 0.92, p < 0.05 respectively), which is the target organ for Cd. Higher concentrations of the elements were found in hare tissue, and this may be linked to pollution of their wild habitats. The xenobiotic elements as well as the essential elements were accumulated in the kidneys of the hares than rabbits. For liver, differences were less pronounced and significance was only for Fe and Cu. Muscle of hares was more contaminated than of rabbits for both biogenic and toxic elements. These results show that detectable concentrations of inorganic elements. These levels may be linked to contamination of the natural habitats of wild biota due to industry, traffic, agriculture, and urban sprawl.

Leishmania donovani Transmission Cycle Associated with Human Infection, Phlebotomus alexandri Sand Flies, and Hare Blood Meals, Israel1.

Cutaneous leishmaniasis caused by Leishmania major or L. tropica and visceral leishmaniasis caused by L. infantum have been reported in Israel. We collected Phlebotomus spp. sand flies in the Negev desert of southern Israel to identify circulating Leishmania spp. Of 22,636 trapped sand flies, 80% were P. alexandri. We sequenced Leishmania-specific internal transcribed spacer 1 fragments and K26 genes. Of 5,019 Phlebotomus female sand flies, 2.5% were Leishmania DNA-positive; 92% of infections were L. donovani. Phylogenetic analyses showed separate clustering of L. donovani and L. infantum. P. alexandri flies positive for L. donovani harbored blood meals from European hares. Leishmania DNA isolated from a patient with cutaneous leishmaniasis who lived in the survey area was identical to L. donovani from P. alexandri flies. We report circulation of L. donovani, a cause of visceral leishmaniasis, in southern Israel. Prompt diagnosis and Leishmania spp. identification are critical to prevent leishmaniasis progression.

Hepatobiliary organoids derived from leporids support the replication of hepatotropic lagoviruses.

The genus Lagovirus of the family Caliciviridae contains some of the most virulent vertebrate viruses known. Lagoviruses infect leporids, such as rabbits, hares and cottontails. Highly pathogenic viruses such as Rabbit haemorrhagic disease virus 1 (RHDV1) cause a fulminant hepatitis that typically leads to disseminated intravascular coagulation within 24-72 h of infection, killing over 95 % of susceptible animals. Research into the pathophysiological mechanisms that are responsible for this extreme phenotype has been hampered by the lack of a reliable culture system. Here, we report on a new ex vivo model for the cultivation of lagoviruses in cells derived from the European rabbit (Oryctolagus cuniculus) and European brown hare (Lepus europaeus). We show that three different lagoviruses, RHDV1, RHDV2 and RHDVa-K5, replicate in monolayer cultures derived from rabbit hepatobiliary organoids, but not in monolayer cultures derived from cat (Felis catus) or mouse (Mus musculus) organoids. Virus multiplication was demonstrated by (i) an increase in viral RNA levels, (ii) the accumulation of dsRNA viral replication intermediates and (iii) the expression of viral structural and non-structural proteins. The establishment of an organoid culture system for lagoviruses will facilitate studies with considerable implications for the conservation of endangered leporid species in Europe and North America, and the biocontrol of overabundant rabbit populations in Australia and New Zealand.

When death comes: linking predator-prey activity patterns to timing of mortality to understand predation risk.

The assumption that activity and foraging are risky for prey underlies many predator-prey theories and has led to the use of predator-prey activity overlap as a proxy of predation risk. However, the simultaneous measures of prey and predator activity along with timing of predation required to test this assumption have not been available. Here, we used accelerometry data on snowshoe hares (Lepus americanus) and Canada lynx (Lynx canadensis) to determine activity patterns of prey and predators and match these to precise timing of predation. Surprisingly we found that lynx kills of hares were as likely to occur during the day when hares were inactive as at night when hares were active. We also found that activity rates of hares were not related to the chance of predation at daily and weekly scales, whereas lynx activity rates positively affected the diel pattern of lynx predation on hares and their weekly kill rates of hares. Our findings suggest that predator-prey diel activity overlap may not always be a good proxy of predation risk, and highlight a need for examining the link between predation and spatio-temporal behaviour of predator and prey to improve our understanding of how predator-prey behavioural interactions drive predation risk.

Does coat colour influence survival? A test in a cyclic population of snowshoe hares.

Some mammal species inhabiting high-latitude biomes have evolved a seasonal moulting pattern that improves camouflage via white coats in winter and brown coats in summer. In many high-latitude and high-altitude areas, the duration and depth of snow cover has been substantially reduced in the last five decades. This reduction in depth and duration of snow cover may create a mismatch between coat colour and colour of the background environment, and potentially reduce the survival rate of species that depend on crypsis. We used long-term (1977-2020) field data and capture-mark-recapture models to test the hypothesis that whiteness of the coat influences winter apparent survival in a cyclic population of snowshoe hares (Lepus americanus) at Kluane, Yukon, Canada. Whiteness of the snowshoe hare coat in autumn declined during this study, and snowshoe hares with a greater proportion of whiteness in their coats in autumn survived better during winter. However, whiteness of the coat in spring did not affect subsequent summer survival. These results are consistent with the hypothesis that the timing of coat colour change in autumn can reduce overwinter survival. Because declines in cyclic snowshoe hare populations are strongly affected by low winter survival, the timing of coat colour change may adversely affect snowshoe hare population dynamics as climate change continues.

Adaptive cellular evolution or cellular system drift in hares.

Adaptations occur at many levels, for example, from DNA sequence of regulatory elements and cellular homeostatic systems to organismal physiology and behaviour (Mayr, 1997). Established adaptations are maintained by purifying and stabilizing selection. Students of animal diversity tend to focus on higher order traits, anatomy, physiology, organismal function and interactions. The core cellular and metabolic systems of metazoans evolved early in their history and are assumed to be rather similar between groups. The housekeeping functions and core metabolic functions of cells are generally considered relatively static, especially among closely related species. The extent to which evolution shapes core cellular metabolism and physiology in animals is largely unexplored. Ecological opportunities or strong positive selection can alter basal metabolic rate, activity levels and life-history traits (e.g., life span, age of maturity, offspring number) and potentially lead to divergence in core cellular and metabolic trait systems (Norin & Metcalfe, 2019; Speakman, 2005). Furthermore, systems under stabilizing selection can also change. Developmental systems of related species may produce the same phenotype or structure, but experience drift that can alter connections and even lead to turnover of cogs in the system (True & Haag, 2001). Are the cellular functions of animals highly constrained, subject to cellular system drift or affected by positive selection? This was tackled by a new study by Kateryna Gaertner and colleagues in a From the Cover manuscript in this issue of Molecular Ecology (Gaertner et al., 2022), using fibroblasts from the closely related but ecologically distinct brown and mountain hares.

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.

Stochastic resonance in climate reddening increases the risk of cyclic ecosystem extinction via phase-tipping.

Human activity is leading to changes in the mean and variability of climatic parameters in most locations around the world. The changing mean has received considerable attention from scientists and climate policy makers. However, recent work indicates that the changing variability, that is, the amplitude and the temporal autocorrelation of deviations from the mean, may have greater and more imminent impact on ecosystems. In this paper, we demonstrate that changes in climate variability alone could drive cyclic predator-prey ecosystems to extinction via so-called phase-tipping (P-tipping), a new type of instability that occurs only from certain phases of the predator-prey cycle. We construct a mathematical model of a variable climate and couple it to two self-oscillating paradigmatic predator-prey models. Most importantly, we combine realistic parameter values for the Canada lynx and snowshoe hare with actual climate data from the boreal forest. In this way, we demonstrate that critically important species in the boreal forest have increased likelihood of P-tipping to extinction under predicted changes in climate variability, and are most vulnerable during stages of the cycle when the predator population is near its maximum. Furthermore, our analysis reveals that stochastic resonance is the underlying mechanism for the increased likelihood of P-tipping to extinction.