Transmission dynamics of lyssavirus in : mechanistic modelling study based on longitudinal seroprevalence data.

We investigated the transmission dynamics of lyssavirus in and , using serological, virological, demographic and ecological data collected between 2015 and 2022 from two maternity colonies in northern Italian churches. Despite no lyssavirus detection in 556 bats sampled over 11 events by reverse transcription-polymerase chain reaction (RT-PCR), 36.3% of 837 bats sampled over 27 events showed neutralizing antibodies to , with a significant increase in summers. By fitting sets of mechanistic models to seroprevalence data, we investigated factors that influenced lyssavirus transmission within and between years. Five models were selected as a group of final models: in one model, a proportion of exposed bats (median model estimate: 5.8%) became infectious and died while the other exposed bats recovered with immunity without becoming infectious; in the other four models, all exposed bats became infectious and recovered with immunity. The final models supported that the two colonies experienced seasonal outbreaks driven by: (i) immunity loss particularly during hibernation, (ii) density-dependent transmission, and (iii) a high transmission rate after synchronous birthing. These findings highlight the importance of understanding ecological factors, including colony size and synchronous birthing timing, and potential infection heterogeneities to enable more robust assessments of lyssavirus spillover risk.