<figure><img src="image_1.png" /></figure>

HAL
open science

# Climate Anomalies and Spillover of Bat-Borne Viral Diseases in the Asia-Pacific Region and the Arabian Peninsula

Alice Latinne, Serge Morand

► To cite this version:

Alice Latinne, Serge Morand. Climate Anomalies and Spillover of Bat-Borne Viral Diseases in the Asia-Pacific Region and the Arabian Peninsula. Viruses, 2022, 14 (5), pp.1100. 10.3390/v14051100. hal-04787436

HAL Id: hal-04787436

https://hal.science/hal-04787436v1

Submitted on 26 Feb 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
Article

# Climate Anomalies and Spillover of Bat-Borne Viral Diseases in the Asia-Pacific Region and the Arabian Peninsula

Alice Latinne <fnref n="1" />,<fnref n="2" />,<fnref n="3" />,<fnref n="4" />,* and Serge Morand <fnref n="3" />,<fnref n="4" />,<fnref n="5" />

<figure><img src="image_3.png" /></figure>

check for
updates

**Citation:** Latinne, A.; Morand, S.
Climate Anomalies and Spillover of
Bat-Borne Viral Diseases in the
Asia-Pacific Region and the Arabian
Peninsula. *Viruses* **2022**, *14*, 1100.
https://doi.org/10.3390/v14051100

Academic Editor: Peng Zhou

Received: 5 April 2022

Accepted: 16 May 2022

Published: 20 May 2022

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_4.png" /></figure>

**Copyright:** © 2022 by the authors.
Licensee MDPI, Basel, Switzerland.
This article is an open access article
distributed under the terms and
conditions of the Creative Commons
Attribution (CC BY) license (https://
creativecommons.org/licenses/by/
4.0/).

**Abstract:** Climate variability and anomalies are known drivers of the emergence and outbreaks of infectious diseases. In this study, we investigated the potential association between climate factors and anomalies, including El Niño Southern Oscillation (ENSO) and land surface temperature anomalies, as well as the emergence and spillover events of bat-borne viral diseases in humans and livestock in the Asia-Pacific region and the Arabian Peninsula. Our findings from time series analyses, logistic regression models, and structural equation modelling revealed that the spillover patterns of the Nipah virus in Bangladesh and the Hendra virus in Australia were differently impacted by climate variability and with different time lags. We also used event coincidence analysis to show that the emergence events of most bat-borne viral diseases in the Asia-Pacific region and the Arabian Peninsula were statistically associated with ENSO climate anomalies. Spillover patterns of the Nipah virus in Bangladesh and the Hendra virus in Australia were also significantly associated with these events, although the pattern and co-influence of other climate factors differed. Our results suggest that climate factors and anomalies may create opportunities for virus spillover from bats to livestock and humans. Ongoing climate change and the future intensification of El Niño events will therefore potentially increase the emergence and spillover of bat-borne viral diseases in the Asia-Pacific region and the Arabian Peninsula.

**Keywords:** bat-borne virus; spillover; SARS-CoV-2; Nipah virus; Hendra virus; climate change; El Niño Southern Oscillation; event coincidence analysis; temporal analysis; structural equation modelling

## 1. Introduction

Most zoonotic and vector-borne diseases are climate-sensitive, particularly to temperature or precipitations [1], and several are also sensitive to climate variability and anomalies [2–4]. The El Niño Southern Oscillation (ENSO), with its alternating warming (El Niño), cooling (La Niña), and neutral phases, is one of the most important climate phenomena due to its ability to modify the global atmospheric circulation and the temperature and precipitation patterns across the globe [5]. The ENSO has significant cascade effects on ecosystems [6–8] and agriculture productivity [9]. ENSO-related climate variability is also a known driver of the emergence and outbreaks of infectious diseases [3,10]. Outbreaks of numerous infectious diseases, such as cholera [11], Rift Valley fever [12], visceral leishmaniasis [13], dengue [14,15], Zika virus [16], and malaria [17], among others, have been linked to the ENSO. The emergences of some viral diseases of bat origin, such as the Hendra virus (HeV) in Australia and the Nipah virus (NiV) in Malaysia, have also been associated with El Niño events [18,19]. Several mechanisms by which the ENSO affects and facilitates the transmission of zoonotic diseases have been suggested, including modifications of seasonal cycles, population dynamics, and distribution ranges of vectors and hosts of

<footnote n="1">¹ Wildlife Conservation Society, Viet Nam Country Program, Ha Noi 100000, Vietnam</footnote>

<footnote n="2">² Wildlife Conservation Society, Global Conservation Program, Bronx, NY 10460, USA</footnote>

<footnote n="3">³ MIVEGEC, CNRS—IRD—Montpellier Université, 911 Avenue Agropolis, BP 6450, 34394 Montpellier, France; 
serge.morand@umontpellier.fr</footnote>

<footnote n="4">⁴ Faculty of Veterinary Technology, University of Kasetsart, Bangkok 10900, Thailand</footnote>

<footnote n="5">⁵ Faculty of Tropical Medicine, University of Mahidol, Bangkok 10400, Thailand</footnote>

<footnote>* Correspondence: alice.latinne@gmail.com</footnote>
zoonotic pathogens, as well as alterations of the replication and transmission patterns of these pathogens [1,3,4].

In this study, we reviewed the emergence events and recurring spillover events of bat-borne viral diseases in humans and livestock in the Asia-Pacific region and the Arabian Peninsula, i.e., two regions highly affected by El Niño/La Niña events [20]. Furthermore, we investigated the potential association between climate anomalies, El Niño/La Niña events, and these emergence and spillover events. First, we tested the potential association between the spillover events of HeV in Australia and NiV in Bangladesh, i.e., two bat-borne viruses characterized by a high number of recurring spillover events in these two regions, and climate factors (temperature, rainfall) and anomalies (ENSO and land surface temperature anomalies) using time-series analyses, logistic regression models, and structural equation modelling. Second, we assessed potential simultaneities between the emergence events of bat-borne viruses in human and livestock populations and El Niño/La Niña events using event coincidence analysis (ECA) [21]. We then discussed the potential ecological mechanisms that may explain the emergence and spillover events of bat-borne viruses in relation to climate factors and events of El Niño/La Niña.

## 2. Materials and Methods

### 2.1. Data

Bat-borne viral pathogens were defined as viruses with bats (Chiroptera) as their natural animal hosts, i.e., the long-term ecological niche of a viral population [22], or viruses whose closest viral relatives have bats as their natural hosts. The zoonotic sources of these viruses in human populations was either bats or another host species involved as an intermediate host in their emergence. Emergence events were defined as the first detection of a bat-borne viral pathogen in human or in livestock populations or the first detection of a bat-borne viral pathogen in a region significantly distant from any other regions where it was previously observed (e.g., Nipah virus emergence events in Malaysia, India, and the Philippines). Recurring spillover events correspond to subsequent pathogen detections after its first emergence in the same region. Data on the emergence of bat-borne viral diseases in human and livestock populations in the study area during the period 1990–2020 were gathered from original sources (Table 1) and from several databases (Emerging Infectious Diseases Repository (EIDR), https://eidr.ecohealthalliance.org/ (accessed on 17 May 2022); World Animal Health Information System (OIE-WAHIS), https://wahis.oie.int/#/home (accessed on 17 May 2022); PROMED). Data on recurring spillover events of the Nipah virus in Bangladesh and India (Table 2) were obtained from three studies [23–25], while data on recurring spillover events of the Hendra virus in Australia were obtained from the Queensland Government database (https://www.business.qld.gov.au/industries/service-industries-professionals/service-industries/veterinary-surgeons/guidelines-hendra/incident-summary (accessed on 17 May 2022)) (Table 3).

<table><caption>**Table 1.** Emergence of bat-borne viruses in the Asia-Pacific region and the Arabian Peninsula in relation to El Niño Southern Oscillation (ENSO)-driven climate anomalies (data on ENSO were retrieved from NOAA).</caption><thead><tr><td>Emergence</td><td>Viral Family</td><td>Natural Reservoir</td><td>Intermediate host</td><td>Date, Location</td><td>ENSO Phase</td><td>References</td></tr></thead><tbody><tr><td>Hendra virus</td><td>Paramyxoviridae</td><td>Pteropodid bats</td><td>Horse</td><td>Aug 1994, Australia</td><td>Warm Phase/El Niño</td><td>Giles et al., 2018 [18]</td></tr><tr><td>Australian bat lyssavirus</td><td>Rhabdoviridae</td><td>Pteropodid bats</td><td>None</td><td>Oct 1996, Australia</td><td>Neutral Phase</td><td>Field et al., 1999 [26]</td></tr><tr><td>Menangle virus</td><td>Paramyxoviridae</td><td>Pteropodid bats</td><td>Pig</td><td>Jun 1997, Australia</td><td>Warm Phase/El Niño</td><td>Chant et al., 1998 [27]</td></tr><tr><td>Nipah virus</td><td>Paramyxoviridae</td><td>Pteropodid bats</td><td>Pig</td><td>Sep 1998, Malaysia</td><td>Cool Phase/La Niña</td><td>Ang et al., 2018 [25]</td></tr><tr><td>Nipah virus</td><td>Paramyxoviridae</td><td>Pteropodid bats</td><td>None</td><td>Jan 2001, India</td><td>Cool Phase/La Niña</td><td>Ang et al., 2018 [25]</td></tr></tbody></table>
<table><caption>**Table 1.** *Cont.*</caption><thead><tr><th>Emergence</th><th>Viral Family</th><th>Natural Reservoir</th><th>Intermediate host</th><th>Date, Location</th><th>ENSO Phase</th><th>References</th></tr></thead><tbody><tr><td>SARS-CoV-1</td><td>Coronaviridae</td><td>Rhinolophid bats</td><td>Small carnivores</td><td>Nov 2002, China</td><td>Warm Phase/El Niño</td><td>Ge et al., 2013 [28]</td></tr><tr><td>Melaka virus</td><td>Reoviridae</td><td>Pteropodid bats</td><td>None</td><td>Mar 2006, Malaysia</td><td>Cool Phase/La Niña</td><td>Chua et al., 2008 [29]</td></tr><tr><td>Kampar virus</td><td>Reoviridae</td><td>Pteropodid bats</td><td>None</td><td>Aug 2006, Malaysia</td><td>Neutral Phase</td><td>Chua et al., 2008 [29]</td></tr><tr><td>MERS-CoV</td><td>Coronaviridae</td><td>Vespertilionid bats</td><td>Camel</td><td>Apr 2012, Middle East</td><td>Cool Phase/La Niña</td><td>Zaki et al., 2012 [30]</td></tr><tr><td>Nipah virus</td><td>Paramyxoviridae</td><td>Pteropodid bats</td><td>Horse</td><td>Mar 2014, The Philippines</td><td>Neutral phase</td><td>Ching et al., 2015 [31]</td></tr><tr><td>SADS-CoV</td><td>Coronaviridae</td><td>Rhinolophid bats</td><td>Pig (no human cases)</td><td>Oct 2016, China</td><td>Cool Phase/La Nina</td><td>Gong et al., 2017 [32]</td></tr><tr><td>SARS-CoV-2</td><td>Coronaviridae</td><td>Rhinolophid bats</td><td>?</td><td>December 2019, China</td><td>Warm Phase/El Niño</td><td>Zhu et al., 2020 [33]</td></tr></tbody></table>

<table><caption>**Table 2.** Recurring spillover events of the Nipah virus after its first emergence in South Asia (India and Bangladesh) in relation to El Niño Southern Oscillation (ENSO)-driven climate anomalies (data on ENSO were retrieved from NOAA). Data on outbreaks were retrieved from Rahman and Chakraborty (2012) [23], Ang et al. (2018) [25], and Rahman et al. (2021) [24].</caption><thead><tr><th>Country</th><th>Date (Month/Year)</th><th>ENSO Phase</th></tr></thead><tbody><tr><td>Bangladesh</td><td>April 2001</td><td>Neutral Phase</td></tr><tr><td>Bangladesh</td><td>January 2003</td><td>Warm Phase/El Niño</td></tr><tr><td>Bangladesh</td><td>January 2004</td><td>Neutral Phase</td></tr><tr><td>Bangladesh</td><td>April 2004</td><td>Neutral Phase</td></tr><tr><td>Bangladesh</td><td>January 2005</td><td>Warm Phase/El Niño</td></tr><tr><td>Bangladesh</td><td>January 2007</td><td>Warm Phase/El Niño</td></tr><tr><td>Bangladesh</td><td>March 2007</td><td>Neutral Phase</td></tr><tr><td>Bangladesh</td><td>April 2007</td><td>Neutral Phase</td></tr><tr><td>India</td><td>April 2007</td><td>Neutral Phase</td></tr><tr><td>Bangladesh</td><td>February 2008</td><td>Cool Phase/La Niña</td></tr><tr><td>Bangladesh</td><td>April 2008</td><td>Cool Phase/La Niña</td></tr><tr><td>Bangladesh</td><td>January 2009</td><td>Cool Phase/La Niña</td></tr><tr><td>Bangladesh</td><td>February 2010</td><td>Warm Phase/El Niño</td></tr><tr><td>Bangladesh</td><td>January 2011</td><td>Cool Phase/La Niña</td></tr><tr><td>Bangladesh</td><td>January 2012</td><td>Cool Phase/La Niña</td></tr><tr><td>Bangladesh</td><td>January 2013</td><td>Neutral Phase</td></tr><tr><td>Bangladesh</td><td>January 2014</td><td>Warm Phase/El Niño</td></tr><tr><td>Bangladesh</td><td>January 2015</td><td>Warm Phase/El Niño</td></tr><tr><td>Bangladesh</td><td>February 2015</td><td>Neutral Phase</td></tr><tr><td>Bangladesh</td><td>March 2015</td><td>Warm Phase/El Niño</td></tr><tr><td>Bangladesh</td><td>February 2017</td><td>Neutral Phase</td></tr><tr><td>Bangladesh</td><td>February 2018</td><td>Cool Phase/La Niña</td></tr><tr><td>Bangladesh</td><td>April 2018</td><td>Cool Phase/La Niña</td></tr><tr><td>India</td><td>May 2018</td><td>Neutral Phase</td></tr><tr><td>India</td><td>June 2019</td><td>Neutral Phase</td></tr></tbody></table>
<table><caption>**Table 3.** Recurring spillover events of the Hendra virus after its first emergence in Australia in relation to El Niño Southern Oscillation (ENSO)-driven climate anomalies (data on ENSO were retrieved from NOAA). Data on outbreaks were retrieved from the Queensland Government database (https://www.business.qld.gov.au/industries/service-industries-professionals/service-industries/veterinary-surgeons/guidelines-hendra/incident-summary, accessed on 17 May 2022).</caption><thead><tr><th>Country</th><th>Date (Month/Year)</th><th>ENSO Phase</th></tr></thead><tbody><tr><td>Australia</td><td>September 1994</td><td>Warm Phase/El Niño</td></tr><tr><td>Australia</td><td>January 1999</td><td>Cool Phase/La Niña</td></tr><tr><td>Australia</td><td>October 2004</td><td>Warm Phase/El Niño</td></tr><tr><td>Australia</td><td>December 2004</td><td>Warm Phase/El Niño</td></tr><tr><td>Australia</td><td>June 2006</td><td>Neutral Phase</td></tr><tr><td>Australia</td><td>October 2006</td><td>Warm Phase/El Niño</td></tr><tr><td>Australia</td><td>June 2007</td><td>Neutral Phase</td></tr><tr><td>Australia</td><td>July 2007</td><td>Cool Phase/La Niña</td></tr><tr><td>Australia</td><td>June 2008</td><td>Cool Phase/La Niña</td></tr><tr><td>Australia</td><td>July 2008,</td><td>Neutral Phase</td></tr><tr><td>Australia</td><td>July 2009</td><td>Neutral Phase</td></tr><tr><td>Australia</td><td>September 2009</td><td>Warm Phase/El Niño</td></tr><tr><td>Australia</td><td>May 2010</td><td>Neutral Phase</td></tr><tr><td>Australia</td><td>June 2011 (4 events)</td><td>Cool Phase/La Niña</td></tr><tr><td>Australia</td><td>July 2011 (8 events)</td><td>Neutral Phase</td></tr><tr><td>Australia</td><td>August 2011 (5 events)</td><td>Cool Phase/La Niña</td></tr><tr><td>Australia</td><td>October 2011</td><td>Cool Phase/La Niña</td></tr><tr><td>Australia</td><td>January 2012</td><td>Cool Phase/La Niña</td></tr><tr><td>Australia</td><td>May 2012 (2 events)</td><td>Neutral Phase</td></tr><tr><td>Australia</td><td>June 2012</td><td>Neutral Phase</td></tr><tr><td>Australia</td><td>July 2012 (2 events)</td><td>Neutral Phase</td></tr><tr><td>Australia</td><td>September 2012</td><td>Neutral Phase</td></tr><tr><td>Australia</td><td>October 2012</td><td>Neutral Phase</td></tr><tr><td>Australia</td><td>January 2013</td><td>Neutral Phase</td></tr><tr><td>Australia</td><td>February 2013</td><td>Neutral Phase</td></tr><tr><td>Australia</td><td>June 2013 (2 events)</td><td>Neutral Phase</td></tr><tr><td>Australia</td><td>July 2013 (4 events)</td><td>Neutral Phase</td></tr><tr><td>Australia</td><td>March 2014</td><td>Neutral Phase</td></tr><tr><td>Australia</td><td>June 2014 (2 events)</td><td>Neutral Phase</td></tr><tr><td>Australia</td><td>July 2014</td><td>Neutral Phase</td></tr><tr><td>Australia</td><td>June 2015</td><td>Warm Phase/El Niño</td></tr><tr><td>Australia</td><td>July 2015</td><td>Warm Phase/El Niño</td></tr><tr><td>Australia</td><td>September 2015</td><td>Warm Phase/El Niño</td></tr><tr><td>Australia</td><td>December 2016</td><td>Cool Phase/La Niña</td></tr><tr><td>Australia</td><td>May 2017</td><td>Neutral Phase</td></tr><tr><td>Australia</td><td>July 2017</td><td>Neutral Phase</td></tr><tr><td>Australia</td><td>August 2017 (2 events)</td><td>Neutral Phase</td></tr><tr><td>Australia</td><td>September 2018</td><td>Neutral Phase</td></tr><tr><td>Australia</td><td>June 2019</td><td>Warm Phase/El Niño</td></tr><tr><td>Australia</td><td>June 2020</td><td>Neutral Phase</td></tr></tbody></table>

Data on ENSO values were retrieved from the National Oceanic and Atmospheric Administration (NOAA, https://www.noaa.gov, accessed on 17 May 2022). The 'NINO 3.4' index is the most commonly used index to define El Niño and La Niña events and to study climate-rainfall or climate-disease connections [3]. The NINO 3.4 index is based on a 5-month running mean of the sea surface temperature (SST) in the region bounded by 5° N to 5° S, from 170° W to 120° W. El Niño (warm phase) and La Niña (cool phase) are defined when anomalies in the NINO 3.4 index exceeds +0.4 °C or -0.4 °C, respectively (NOAA, https://www.noaa.gov, accessed on 17 May 2022). The R package rsi0 [34] was used to import the NINO 3.4 index values for the period 1990-2020 and the corresponding defined El Niño (warm phase) and La Niña (cool phase) phases from the NOAA website.
Data on the average monthly temperature and rainfall data in Australia and Bangladesh were also gathered from the World Bank database (https://climateknowledgeportal.worldbank.org, accessed on 17 May 2022) and data on the global land surface temperature anomalies were obtained from the NOAA (https://www.ncei.noaa.gov/access/monitoring/global-temperature-anomalies, accessed on 17 May 2022) to investigate the impact of other climate factors (rainfall, temperature) and climate variability (global land surface temperature anomalies), in addition to the ENSO anomalies, on the occurrence of recurring spillover events of HeV and NiV viruses.

## 2.2. Statistical Analyses

First, we investigated the potential association between the spillover events of HeV in Australia and NiV in Bangladesh and several climate factors. For this, we used: (1) time-series analyses to investigate the temporal association between HeV/NiV spillover events, and temperature, rainfall, ENSO, and land surface temperature anomalies and to estimate their time lag values; (2) logistic regression models to determine the significant factors (temperature, rainfall, ENSO, land surface temperature anomalies) explaining the spillover events using the time lag values computed from the results of the time-series analyses; and (3) structural equation modelling to test a causal chain of correlation that may explain the spillover events using the results of the logistic regression analyses. Second, we assessed potential simultaneities between the emergence events of all bat-borne viruses in human and livestock populations and El Niño/La Niña events using event coincidence analysis to test the hypothesis that HeV/NiV spillover events were statistically preceded by an event of El Niño/La Niña.

### 2.2.1. Time-Series Analysis

Time-series analyses were used to study the temporal patterns of ENSO anomalies (using the NINO 3.4 index), as well as the average monthly temperature and rainfall in Australia and Bangladesh using the ncf function implemented in R [35]. The time series included 330 months in total from January 1993, i.e., one year before the first spillover event recorded in our dataset, to June 2020. The residual autocorrelation function (ACF) was examined to determine the general form of the model to be fitted. A wavelet analysis was used to decompose a time series to reveal periodic signals at each time point in the series. The wavelet analysis coefficients show the correlation magnitudes of ENSO anomalies (NINO 3.4 index), temperature, or rainfall for each year and period length of the time series (i.e., 1993 to 2020), displayed using a power spectrum over the full time series using the biwavelet and WaveletComp packages [36,37] implemented in R [38]. The ccf function was then used to compute the cross-correlation or cross-covariance between univariate series, i.e., ENSO (NINO 3.4 index); the average monthly temperature; the average monthly rainfall; land surface temperature anomalies; and either HeV or NiV recurring spillover events in Australia and in Bangladesh, respectively.

### 2.2.2. Logistic Regression with Time Lag Analysis

Logistic regression modelling with a logit function and lag was used to test the significant effects of monthly rainfall; monthly temperature; anomalies in land surface temperature; and ENSO anomalies (NINO 3.4 index) on the recurring spillover events of HeV and NiV in Australia and Bangladesh, respectively. The most significant lag values computed by time-series cross-correlation analysis, as described above, and the glm function implemented in R with the family binomial [38] were used.

The initial general linear model with logit function was of the form:

Recurrent spillover of HeV/NiV ~ lag(NINO 3.4 index, lag value)

+ lag(average temperature Australia/Bangladesh, lag value)

+ lag(average rainfall Australia/Bangladesh, lag value)

+ lag(global land surface temperature anomalies, lag value)
Initial models included variables with significant lag values obtained from cross-correlation time-series analysis. Final models were selected using backward selection and AIC criterion using the stepAIC function of the MASS package [39] implemented in R.

### 2.2.3. Structural Equation Modelling

Structural equation modelling (SEM) was used to investigate the temporal relationships between recurring outbreaks of HeV and NiV, respectively, in Australia and Bangladesh, in relation to monthly rainfall, monthly temperature, anomalies in land surface temperature, and the NINO 3.4 index values. SEM combines measurement models (e.g., reliability) with structural models (e.g., regression), thus testing a chain of causality (path analysis) between outbreaks of HeV and NiV in Australia and Bangladesh and these climatic factors. SEM was performed using the 'piecewiseSEM' package [40]. The following structural equation model was tested for the period from January 1993 to June 2020 of the dataset:

$$F (\text{outbreaks of HeV or NiV}) = f_1 (\text{lag NINO 3.4 index}) + f_2 (\text{lag temperature Australia or Bangladesh}) + f_3 (\text{lag rainfall Australia or Bangladesh}) + f_4 (\text{lag global land surface temperature anomalies}) + b_1$$

$$G (\text{lag temperature Australia or Bangladesh}) = g_1 (\text{lag NINO 3.4 index}) + b_2$$

$$H (\text{lag rainfall Australia or Bangladesh}) = h_1 (\text{lag NINO 3.4 index}) + b_3$$

$$I (\text{lag anomalies of the land surface temperature}) = i (\text{lag NINO 3.4 index}) + b_4$$

with lag values computed by time-series cross-correlation analysis.

### 2.2.4. Event Coincidence Analysis

Finally, event coincidence analysis (ECA) was used to test if events of a given type are causally influenced by the timing of events of second type [41] and to investigate the statistical interdependence between emergence and spillover events of bat-borne viruses and El Niño/La Niña events. ECA was implemented in the CoinCalc R package [42] to test whether the observed coincidence rates are significantly different from two independent random events [41]. ECA defines the precursor coincidence rate (pcr) and the trigger coincidence rate (tcr). The pcr describes the fraction of first-type events, i.e., emergence/spillover events, preceded by at least one second-type event, i.e., El Niño/La Niña events. The tcr describes the fraction of second-type events, i.e., El Niño/La Niña events, followed by at least one first-type event, i.e., emergence/spillover events (see [41]). CoinCalc computed the probability of the precursor and trigger coincidence rates occurring by chance, with the null hypothesis that the observed precursor and trigger coincidence rates can be explained by two independent series of randomly distributed events [42]. The p-value of the corresponding analytical significance test corresponds to the probability that the two types of events are randomly distributed and independent of each other (following two independent Poisson processes) and sufficiently rare.

We tested the hypothesis that the pcr describing the emergence of bat-borne viruses or recurrent spillover events of NiV and HeV were statistically preceded by an El Niño/La Niña event, while the tcr did not depart from a random association. For the analysis of the recurrent spillover events of NiV in Bangladesh and HeV in Australia, the time lag values estimated in months using cross-correlation among the time series of outbreak events and NINO 3.4 index values were used. The time lag values were also moved around their estimates to explore the stability and constancy of the association given by ECA.

## 3. Results

A total of ten bat-borne viruses, belonging to the Coronaviridae (n = 4), Paramyxoviridae (n = 3), Reoviridae (n = 2), and Rhabdoviridae (n = 1) families, emerged in the Asia-Pacific region and the Arabian Peninsula in the period 1990–2020 (Table 1; Figure 1). Nine of these viruses emerged in humans, while the swine acute diarrhea syndrome coronavirus (SADS-CoV) emerged in swine populations but was never detected in humans. Natural bat hosts of the Coronaviridae viruses were vespertilionid and rhinolophid bats,
while their intermediate hosts included several mammal species. Pteropodid bats were the hosts of the emerging Paramyxoviridae, Reoviridae, and Rhabdoviridae, and livestock (horses and pigs) were involved as intermediate hosts in the emergence of HeV and NiV. Most of these viruses emerged in a single geographic area; only NiV emerged in several distant locations over a 16-year period. It first emerged in Malaysia in 1998, before then emerging in Bangladesh and India in 2001 and in the Philippines in 2014 (Table 1; Figure 1). After their first emergence, NiV and HeV then regularly spilled over in Bangladesh and Australia, respectively (Tables 2 and 3).

<figure><img src="image_5.png" /><figcaption>**Figure 1.** (A) Map showing the locations of emergence of bat-borne viruses in the Asia-Pacific region and the Arabian Peninsula (see Table 1) and the bat reservoir of each virus. Virus names are colored</figcaption></figure>
according to the ENSO phase at the time of their emergence: neutral phase (black), cool-phase La Niña (blue), or warm-phase El Niño (red). (B) Variations of the NINO 3.4 index characterizing the El Niño Southern Oscillation (ENSO) retrieved from the National Oceanic and Atmospheric Administration (NOAA, https://www.noaa.gov, accessed on 17 May 2022) from 1990 to 2020. Red and blue threshold lines indicate warming El Niño or cooling La Niña climate anomalies, respectively. Arrows indicate the emergence time of new bat-borne viruses in the Asia-Pacific region and the Arabian Peninsula (see Table 1). Virus names are colored according to the ENSO phase at the time of their emergence: neutral phase (black), cool-phase La Niña (blue), or warm-phase El Niño (red).

Five emergence events, including NiV in Malaysia and India, Melaka virus, Middle East respiratory syndrome coronavirus (MERS-CoV), and SARS-CoV, occurred during a cool phase (La Niña event), while four of them, i.e., HeV, Menangle virus, and severe acute respiratory syndrome coronavirus 1 and 2 (SARS-CoV-1 and -2), occurred during a warm phase (El Niño event) (Figure 1; Table 1). The remaining three emergence events of bat-borne viruses, Kampar virus, Australian bat lyssavirus, and NiV in the Philippines occurred during a neutral phase. However, the spillover of the Australian bat lyssavirus to a human might not be considered as a natural emergence since the infection was acquired from a pet fruit bat living in captivity [26]. It should also be noted that the emergence of NiV in the Philippines in March 2014 followed the major volcanic activity of Mayon volcano that started in May 2013 [43].

### 3.1. *Time-Series Analyses for HeV and NiV*

There were strong and significant seasonal patterns of 12 months for temperature and rainfall, both in Bangladesh and Australia, as shown by the ACF and wavelet analysis (Figure 2A–D). Significant patterns over 24 and 36 months were observed for NINO 3.4 index (Figure 2E). Moreover, an increasing trend in the global land surface temperature anomalies was observed from 1990 to 2020 (Figure S1).

Cross-correlation analysis among pairs of temporal series of spillover events of NiV and HeV revealed several significant correlations with climate variables (Figure 3; Table 4). Significant correlations were observed for recurring spillover events of NiV with monthly rainfall (lag of 1 month), monthly temperature (lag of 1 month), and land surface temperature anomalies (lag of 10 months) (Figure 3A,B,D; Table 4). No significant correlation was observed for recurring spillover events of NiV with NINO 3.4 index values, although the best correlation was observed for no lag (Figure 3C; Table 4). Significant correlations were observed for recurring spillover events of HeV with NINO 3.4 index values (lag of 7 months), monthly rainfall (lag of 1 month), monthly temperature (no lag), and land surface temperature anomalies (lag of 3 months) (Figure 3E–H; Table 4). Cross-correlation analysis among pairs of temporal series of monthly rainfall, monthly temperature, and land surface temperature anomalies revealed few significant correlations with NINO 3.4 index values (Table 4). There was a significant correlation between land surface temperature anomalies and NINO 3.4 index values (lag of 3 months) and a significant correlation between monthly rainfall in Australia and NINO 3.4 index values (no lag) (Table 4). Non-significant correlations were observed for the monthly temperature in Australia (lag of 7 months), as well as for the monthly temperature and monthly rainfall in Bangladesh (with lags of 10 and 11 months, respectively) (Table 4).

### 3.2. *Logistic Regression Analyses*

The above results were used to build two initial logistic regression models based on the lag values obtained by the time-series cross-correlation analysis (see Figure 3 and Table 4). The selected model of the recurring spillover events of HeV in Australia show the significant effects of temperature (with no lag), the global land surface temperature anomalies (with a lag of 3 months), and NINO 3.4 index values (with a lag of 7 months). Rainfall was retained as a variable in the best explanatory model for HeV but had no significant effect (Table 5). The selected model of the recurrent spillover events of NiV in Bangladesh show the only significant effect of rainfall (with a lag of one month), but no
effects of global land surface temperature anomalies, the NINO 3.4 index, and the mean
temperature (Table 5).

<figure><img src="image_6.png" /><figcaption>**Figure 2.** Time series and residual autocorrelation function (ACF) with significant auto-correlation values in dashed lines (left column) and wavelet power spectrum (right column) from January 1993 to</figcaption></figure>
June 2020 (330 months) of (A) monthly temperature in Bangladesh, (B) monthly rainfall in Bangladesh, (C) monthly temperature in Australia, (D) monthly rainfall in Australia, and (E) NINO 3.4 index values decomposed in smooth trend and seasonal effect. Wavelet power values increased from blue to red, and black contour lines indicate a 5% significance level.

<figure><img src="image_7.png" /><figcaption>**Figure 3.** Temporal correlation from January 1993 to June 2020 (330 months) with significant auto-correlation values in dashed lines between spillover events of the Nipah virus in Bangladesh and the Hendra virus in Australia, as well as the monthly temperature (A,E), the monthly rainfall (B,F), the NINO 3.4 index (C,G), and anomalies of the land surface temperature (D,H).</figcaption></figure>
<table><caption>**Table 4.** Results of temporal cross-association between spillover events of the Hendra virus (HeV), the Nipah virus (NiV), NINO 3.4 index, the average monthly temperature in Australia and Bangladesh, the average monthly rainfall in Australia and Bangladesh, and land surface temperature anomalies (monthly lag values were obtained from time-series cross-correlation analysis). Significant correlations are highlighted in bold.</caption><thead><tr><th>First Time-Series</th><th>Second Time-Series</th><th>Lag</th><th>Correlation (p Value)</th></tr></thead><tbody><tr><td rowspan="3">Spillover events of HeV (Australia)</td><td>NINO 3.4 index</td><td>7 months</td><td><strong>0.13 (0.018)</strong></td></tr><tr><td>Rainfall (Australia)</td><td>1 month</td><td><strong>0.18 (0.002)</strong></td></tr><tr><td>Temperature (Australia)</td><td>0 month</td><td><strong>0.24 (&lt; 0.001)</strong></td></tr><tr><td rowspan="6">Spillover events of NiV (Bangladesh)</td><td>Land surface temperature anomalies</td><td>3 months</td><td><strong>0.14 (0.012)</strong></td></tr><tr><td>NINO 3.4 index</td><td>0 month</td><td>0.05 (0.35)</td></tr><tr><td>Rainfall (Bangladesh)</td><td>1 month</td><td><strong>0.22 (0.008)</strong></td></tr><tr><td>Temperature (Bangladesh)</td><td>1 month</td><td><strong>0.32 (0.008)</strong></td></tr><tr><td>Land surface temperature anomalies</td><td>10 months</td><td><strong>0.10 (0.047)</strong></td></tr><tr><td>Rainfall (Australia)</td><td>0 month</td><td><strong>0.16 (0.007)</strong></td></tr><tr><td rowspan="5">NINO 3.4 index</td><td>Temperature (Australia)</td><td>7 months</td><td>0.08 (0.35)</td></tr><tr><td>Rainfall (Bangladesh)</td><td>10 months</td><td>0.03 (0.58)</td></tr><tr><td>Temperature (Bangladesh)</td><td>11 months</td><td>0.06 (0.27)</td></tr><tr><td>Land surface temperature anomalies</td><td>3 months</td><td><strong>0.40 (&lt;0.0001)</strong></td></tr></tbody></table>

<table><caption>**Table 5.** Results of the logistic regression modelling with lags to explore the temporal association between spillover events of the Hendra virus (HeV) and the Nipah virus (NiV) from January 1993 to June 2020. The initial models included the following variables: the NINO 3.4 index, the average monthly temperature in Australia/Bangladesh, the average monthly rainfall in Australia/Bangladesh, and the land surface temperature anomalies with lag values obtained from time-series cross-correlation analysis (Table 4). The best explanatory models were selected using a backward procedure with AIC criterion. Significant *p*-values are highlighted in bold.</caption><thead><tr><th>Response Variable</th><th>Predictor Variable</th><th>Estimate (Std Err)</th><th>Odds Ratio (2.5–97.5 %)</th><th>p</th><th>R2 (Global)</th></tr></thead><tbody><tr><td rowspan="4">HeV spillover events (Australia)</td><td>NINO 3.4 index (lag = 7 months)</td><td>−0.72 (0.23)</td><td>0.49 (0.30–0.75)</td><td><strong>0.002</strong></td><td></td></tr><tr><td>Rainfall (lag = 1 month)</td><td>−0.01 (0.01)</td><td>0.99 (0.96–1.00)</td><td>0.16</td><td></td></tr><tr><td>Temperature (lag = 0 month)</td><td>−0.11 (0.04)</td><td>0.90 (0.82–0.98)</td><td><strong>0.013</strong></td><td></td></tr><tr><td>Land surface temperature anomalies (lag = 3 months)</td><td>3.43 (1.03)</td><td>30.96 (4.32–254.30)</td><td><strong>0.001</strong></td><td>0.21</td></tr><tr><td>NiV spillover events (Bangladesh)</td><td>Rainfall (lag = 1 month)</td><td>−0.03 (0.01)</td><td>0.97 (0.95–0.99)</td><td><strong>0.008</strong></td><td>0.30</td></tr></tbody></table>

### 3.3. *Structural Equation Modelling*

SEM confirmed the above results for HeV. Significant correlations were observed between the series of recurring HeV spillovers and the mean monthly temperatures in Australia (with no lag, *p* = 0.003), the anomalies in the land surface temperature (with a lag of 3 months, *p* < 0.001), and NINO 3.4 index values (with a lag of 7 months, *p* < 0.001) (Figure 4A; Table 6). The mean monthly rainfall (with a lag of 1 month) and anomalies of the land surface temperature (with a lag of 3 months) were also significantly correlated with NINO 3.4 index values (*p* = 0.003 and *p* < 0.001, respectively) (Figure 4A; Table 6), with respective lags taking into account their estimated values given in Table 4.
A. Hendra virus in Australia

<figure><img src="image_8.png" /></figure>

B. Nipah virus in Bangladesh

<figure><img src="image_9.png" /><figcaption>**Figure 4.** Results of structural equation modelling of (A) spillover events of the Hendra virus in Australia and (B) spillover events of the Nipah virus in Bangladesh on temporal trends from January 1993 to June 2020 (330 months) based on results obtained from temporal series analyses (Figure 2) and logistic regression analyses with lags (Table 4). Significant partial correlations are presented in continuous lines and non-significant partial correlations are presented in dashed lines with values of standardized estimates.</figcaption></figure>
<table><caption>**Table 6.** Results of the structural equation modelling (SEM) to explore the temporal associations between spillover events of the Hendra virus (HeV) or the Nipah virus (NiV), the monthly temperature, the monthly rainfall, the NINO 3.4 index, and land surface temperature (LST) anomalies from January 1991 to June 2020. Lag values were obtained from the time-series analyses and the logistic regression analyses (Table 4). Significant *p*-values are highlighted in bold.</caption><thead><tr><th>Model</th><th>Response Variable</th><th>Predictor Variable</th><th>Estimate (Std Err), df</th><th>Standardized Estimate</th><th>p</th><th>R2 (Individual)</th></tr></thead><tbody><tr><td rowspan="5">HeV</td><td rowspan="5">Spillover events of HeV</td><td>'NINO 3.4' (lag = 7 months)</td><td>-0.077 (0.022), 313</td><td>-0.204</td><td><strong>&lt;0.001</strong></td><td></td></tr><tr><td>Rainfall (lag = 1 month)</td><td>-0.001 (0.001), 313</td><td>-0.067</td><td>0.28</td><td></td></tr><tr><td>Temperature (lag = 0 month)</td><td>-0.078 (0.022), 313</td><td>-0.190</td><td><strong>0.003</strong></td><td></td></tr><tr><td>LST anom (lag = 3 months)</td><td>0.3123 (0.095), 313</td><td>0.201</td><td><strong>&lt;0.001</strong></td><td>0.11</td></tr><tr><td>Rainfall (lag = 1 month)</td><td>'NINO 3.4' (lag = 1 month)</td><td>-5.223 (1.7670), 316</td><td>-0.1638</td><td><strong>0.003</strong></td><td>0.027</td></tr><tr><td rowspan="3">NiV</td><td rowspan="3">Spillover events of NiV</td><td>Temperature</td><td>'NINO 3.4' (lag = 7 months)</td><td>0.271 (0.313), 316</td><td>0.049</td><td>0.389</td><td>0.002</td></tr><tr><td>LST anomalies (lag = 3 months)</td><td>'NINO 3.4' (lag = 6 months)</td><td>0.101 (0.012), 316</td><td>0.434</td><td><strong>&lt;0.001</strong></td><td>0.19</td></tr><tr><td>Rainfall (lag = 1 month)</td><td>'NINO 3.4' (lag = 0 month)</td><td>-0.014 (0.015), 312</td><td>-0.887</td><td>0.38</td><td></td></tr><tr><td rowspan="5">NiV</td><td rowspan="5">Spillover events of NiV</td><td>Temperature (lag = 1 month)</td><td>-0.020 (0.005), 312</td><td>-3.689</td><td><strong>0.003</strong></td><td></td></tr><tr><td>LST anom (lag = 10 month)</td><td>0.098 (0.067), 312</td><td>1.443</td><td>0.15</td><td>0.12</td></tr><tr><td>Rainfall (lag = 1 month)</td><td>NINO 3.4' (lag = 11 months)</td><td>6.730 (11.178), 315</td><td>0.602</td><td>0.55</td><td>0.001</td></tr><tr><td>Temperature (lag = 1 month)</td><td>NINO 3.4' (lag = 12 months)</td><td>0.254 (0.231), 315</td><td>1.100</td><td>0.272</td><td>0.004</td></tr><tr><td>LST anomalies (lag = 10 months)</td><td>'NINO 3.4' (lag = 13 months)</td><td>0.092 (0.012), 315</td><td>7.774</td><td><strong>&lt;0.001</strong></td><td>0.16</td></tr></tbody></table>

We included the same variables in SEM for NiV and HeV using the lag with the best correlation, although some were non-significant (Table 4). SEM shows that the series of recurring NiV spillovers was only correlated with the mean monthly temperature in Bangladesh (with a lag of one month, *p* < 0.001) (Table 6; Figure 4B), while the monthly rainfall had no significant effect, contrary to what was suggested by the GLM (Table 5). The land surface temperature anomalies were evidently significantly correlated with NINO 3.4 index values (*p* < 0.001), with lag taking into account values given in Table 4.

### 3.4. Event Coincidence Analysis

Using the locations and dates of the emergence, as well as the series of spillover events of bat-borne viruses (Tables 1–3) and the corresponding ENSO phases (warm, neutral, cool) (Figure 1B; Table 1), we tested the hypothesis that the outbreaks of bat-borne viral diseases were directly preceded (no lag) by an ENSO-driven El Niño/La Niña climate event. The results of the ECA show a random association between an emergence event of bat-borne viral disease (Table 1) following an El Niño/La Niña event (n = 12) given
by the non-significant value of the precursor coincidence rate (0.67, *p* = 0.066) and the non-significant value of the trigger coincidence rate (0.05, *p* = 0.28) (Figure 5A). However, a non-random association was observed with a significant value of precursor coincidence rate (0.80, *p* = 0.014) when the Australian bat lyssavirus (1996) and the Nipah virus in the Philippines (2014) were removed, as other factors may have impacted these emergence events, as explained above.

<figure><img src="image_10.png" /><figcaption>**Figure 5.** Event coincidence analyses of the association between (A) an emergence event of bat-borne virus in the Asia-Pacific region and the Arabian Peninsula, (B) a spillover event of the Nipah virus in Bangladesh, and (C) a spillover event of the Hendra virus in Australia with an ENSO event (El Niño or La Niña phases) with values of precursor coincidence rate (pcr) and its associated probability and lag value between paired events (lags values used for Nipah virus and Hendra virus analyses correspond to the results of cross-temporal series correlations, shown in Figure 3 and Table 4).</figcaption></figure>

Even if there were no significant associations between NINO 3.4 index values and the NiV outbreak events using time-series analyses or SEM (see above), we explored a possible association using event coincidence analysis. We observed a random statistical relationship between an outbreak event of NiV in Bangladesh (Table 2) and events of El Niño/La Niña (n = 22) with no lag (best correlation observed in our cross-correlation analysis, as shown in Table 4). However, a non-random statistical and highly significant relationship between an outbreak event of NiV in Bangladesh and events of El Niño/La Niña using a lag of 3 months was observed, with a significant precursor coincidence rate (0.73, *p* = 0.003) and the non-significant value of a trigger coincidence rate (0.10, *p* = 0.09), suggesting a global lag effect of an ENSO event, whatever its phase (El Niño or La Niña) (Figure 5B). There was a non-random statistical relationship observed between an outbreak
event of HeV in Australia (n = 40) following an event of El Niño/La Niña (a significant precursor coincidence rate = 0.63, *p* = 0.006; a trigger coincidence rate = 0.14, *p* = 0.15) (Figure 5C)), with a lag of 7 months estimated by cross-correlation time-series analysis (Figure 3G; Table 4).

**4. Discussion**

Numerous studies have investigated the origins and drivers of emergence of bat-borne viruses [44–48], with several pinpointing the importance of climate factors and their variability [49,50]. Abnormal rainfall, temperature, and vegetation development associated with ENSO climatic anomalies are known to create appropriate ecological conditions for pathogen emergence, transmission, and propagation [3]. Here, our main objective was to further investigate the potential statistical correlation between emergence and spillover events of several bat-borne viruses and climate factors, such as rainfall, temperature, global surface temperature anomalies, and ENSO events, using a single analytical framework. We used diverse methodologies, such as time-series analyses, logistic regression, SEM, and ECA, to better depict the complex relationships between these climate factors, variability, and anomalies, as well as the emergence and spillover events of bat-borne viruses in the Asia-Pacific region and the Arabian Peninsula.

While it has already been suggested that the emergences of HeV in Australia in 1994 and NiV in Malaysia in 1998 were associated with El Niño events [18,19,51], the long-term surveillance data of NiV in Bangladesh and HeV in Australia provided a better assessment of the influence of climate variability on the recurring spillovers of these bat-borne viruses. Our findings revealed that the spillover patterns of these two closely related paramyxoviruses, both belonging to the genus henipavirus and with closely related flying fox hosts (*Pteropus* spp.), are differently impacted by climate variability and with different time lags, according to our time-series cross-correlation analysis.

NiV outbreaks occurred almost annually in Bangladesh since 2001, following a seasonal pattern with most outbreaks occurring during the winter months [52]. Even if some livestock and domestic animals were found infected by NiV in Bangladesh [53], contact with sick livestock or domestic animals is not considered an important risk factor of spillover infections in Bangladesh [54]. The seasonal timing and spatial distribution of outbreaks coincide with patterns of raw date palm sap production and consumption [52], suggesting that human behavior and the consumption of date palm sap contaminated by *Pteropus* bats play an important role in these spillovers [54]. However, the fact that the number of NiV spillover events varies greatly from year to year suggests that additional factors influencing bat ecology and movement must be at play [55–57]. This was confirmed by a serological survey of *Pteropus medius* bats in Bangladesh which indicated that NiV viral shedding by bats can happen at any time of year and that viral dynamics are cyclical, but not annual or seasonal [58].

So far, a single climatic factor, winter temperature, was shown to be linked to NiV spillover, with colder winter temperatures being associated with more spillovers [56]. Our study confirms the influence of temperature, but our logistic regression models also suggest a correlation between NiV spillover and monthly rainfall, with a short time lag of one month, with lower rainfall being associated with more spillover events. Our cross-correlation analysis, logistic regression models, and SEM did not detect any significant direct correlation between the NINO 3.4 index and NiV spillover events. However, our ECA findings suggest that NiV spillovers in Bangladesh were significantly associated with ENSO events, either El Niño or La Niña phases, with a time lag of three months. This suggests that climate anomalies related to both warm and cool ENSO events may be linked to increased risk of NiV spillover from bats in Bangladesh, which may be explained by the fact that El Niño and La Niña events are characterized by similar rainfall and temperature anomalies in large regions of Bangladesh [59]. ENSO events are associated with the incidence of other diseases in Bangladesh, such as dengue and cholera [60,61]. Climate change will lead to warmer winter temperatures in Bangladesh over the next few decades [62], which
may reduce the number of NiV spillover events in the country, as the negative correlation between temperature and NiV outbreaks was clearly demonstrated in this study and others [56]. However, droughts are expected to significantly increase in some regions of Bangladesh under global warming [62,63]. This could negatively impact bat food resources, induce increased bat movement, and potentially lead to more NiV spillover events, as shown by the negative correlation between rainfall and NiV spillover events observed in this study.

Our study revealed a stronger correlation between climate variability and the spillover pattern of HeV in eastern Australia than NiV in Bangladesh. HeV prevalence in flying foxes in Australia has shown multi-year inter-epidemic periods, suggesting that viral dynamics are not annual, but the ecological drivers and the climate influence behind this pattern remain unclear [64–66]. Numerous factors including food shortage, low concentration of nectar-based resources, extreme temperatures, dry conditions, phenology of eucalypt forests, physiological stress, flying fox foraging behavior, and use of wintering roosts in urban and agricultural areas were all suggested to be associated with increased HeV shedding in Australian flying foxes [18,66,67]. Our logistic regression models and SEM show that seasonal climate factors (monthly temperature), but also multi-annual climate variability (ENSO 3.4 index) and long-trend climate anomalies (land surface temperature anomalies), significantly influence the complex pattern of HeV spillover events in Australia. Our ECA also confirmed the hypothesis that HeV outbreaks were preceded by an ENSO-driven El Niño/La Niña climate event with a time lag of seven months. Interestingly, our results show no time lag between the mean monthly temperature recorded in Australia and the spillovers of HeV, suggesting a direct influence of climate seasonality. McMichael et al. [68] hypothesized that this correlation between lower winter temperature and increased HeV shedding in flying foxes could be mediated by the physiological cost of thermoregulation. The temporal lags observed between the ENSO 3.4 index (7 months) or the anomalies of the land surface temperature (3 months) suggest an indirect effect of the climate variability through ecological cascades that may affect food availability, bat migration patterns, and physiological stresses. A previous study [18] showed that the significant impact of ENSO on the flowering phenology of eucalypt, and consequently on bat foraging activities, was characterized by a time lag (3–8 months) similar to the one we observed between ENSO and HeV spillover events in this study (7 months). Lower eucalypt flowering and bat foraging activities induced by an El Niño event may lead to increased HeV prevalence a few months later [18].

Stress induced by climate variability can have a profound effect on disease dynamics in wild animal populations, mostly in relation to immune changes [69] or behavioral changes, such as climate-driven temporary migrations [18]. Bats undergo seasonal physiological changes, including immunological functions, which affect viral shedding [69]. Flying fox immunocompetence is challenged during food shortages driven by climatic anomalies, and HeV (sero) prevalence in Australian pteropid bats increased when their body condition decreased [67,70]. Immunological stress caused by physiological and behavioral changes during the breeding season has been suggested as a contributing factor in HeV shedding in some studies [67], but has been found to have no effect in others [70].

Beyond NiV and HeV, our findings suggest that the emergence of most viral diseases of bat origin was likely driven by ENSO climatic anomalies, as 9 out of 12 bat-borne viruses emerged in the Asia-Pacific region and the Arabian Peninsula after an ENSO event over the last three decades (Table 1; Figure 1). Removing the emergence of the Australian bat lyssavirus and NiV in the Philippines, given that other factors may have impacted these emergence events, gave a high prior probability for the emergence of bat virus after an event of El Niño/La Niña in our ECA (Figure 5). The recent emergence of SARS-CoV-2, responsible of the coronavirus disease COVID-19 in China in late 2019, also followed an important El Niño event, which had particularly affected China [71].

The emergence of Nipah in the Philippines in March 2014 was not linked to a warm or cool ENSO phase, but did occur following major volcanic activity at Mt Mayon that started
one year prior to the emergence [43]. A study conducted after the more recent eruption of Mt. Mayon in 2018 showed large vegetation and environmental impacts of this eruption, which have affected the whole archipelago up to the northern part of Borneo [72]. Studies have also stressed the likely impacts of volcanic activities on disease outbreaks [73].

Retrospectively, our findings question the absence of emergence reports during major El Niño-La Niña events before 1994. This may be related to an important limitation of our study and the fact that several past spillover events of bat-borne viruses likely remained undetected. The successful detection of spillover events requires an efficient surveillance system adapted to wildlife, or bats in this case, and many countries are still lacking such a wildlife and human health surveillance system [74]. It is also important to note that climate is not the only factor influencing the emergence of bat-borne viruses. Several additional key drivers that promote cross-species transmission and emergence of zoonotic pathogens have been identified in Asia and include rapidly urbanizing populations, widespread wildlife trade and wildlife consumption, intensive livestock production, deforestation, habitat fragmentation, land-use change, and biodiversity loss [75,76]. Therefore, the absence of bat-borne virus emergence report before 1994 may also be explained by increased contact between bat and human/livestock populations over the last three decades due to land-use change, deforestation, intensification of farming practice, and the expansion of the distribution range of certain bat species linked to climate changes [49,75,76].

Climate modelling strongly suggests an intensification of extreme El Niño events in the future [77,78], which will potentially increase the occurrence and outbreaks of infectious diseases and the emergence of bat-borne viral diseases. Climate change will also continue to shift the global distribution of bats and drive changes in bat richness which will increase the risk of bat-borne coronavirus emergence in the near future [49]. This study and its findings also stress the necessity of improving our knowledge of bat ecology. Close monitoring of bat populations will improve our understanding of viral spillover mechanisms, as demonstrated for HeV in eastern Australia and NiV in Bangladesh, and will contribute to better prediction and prevention strategies.

**Supplementary Materials:** The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/v14051100/s1. Figure S1: Global land surface temperature anomalies 1993–2020. Table S1: Emergence and spillover data used in the study. R code S1: R code used to conduct the analyses described in the study. R code S2: R code used to obtain Figure 1B.

**Author Contributions:** Both authors contributed equally. All authors have read and agreed to the published version of the manuscript.

**Funding:** This work was part of the FutureHealthSEA project funded by the French ANR (ANR-17-CE35-0003-01).

**Institutional Review Board Statement:** Not applicable.

**Informed Consent Statement:** Not applicable.

**Data Availability Statement:** Data and codes used in the study are available as Supplementary Materials (Table S1, R code S1, R code S2).

**Conflicts of Interest:** The authors declare no conflict of interest.

## References

1. McIntyre, K.M.; Setzkorn, C.; Hepworth, P.J.; Morand, S.; Morse, A.P.; Baylis, M. Systematic Assessment of the Climate Sensitivity of Important Human and Domestic Animals Pathogens in Europe. *Sci. Rep.* **2017**, *7*, 7134. [CrossRef] [PubMed]

2. Anyamba, A.; Chretien, J.-P.; Small, J.; Tucker, C.J.; Linthicum, K.J. Developing global climate anomalies suggest potential disease risks for 2006–2007. *Int. J. Health Geogr.* **2006**, *5*, 60. [CrossRef] [PubMed]

3. Anyamba, A.; Chretien, J.-P.; Britch, S.C.; Soebiyanto, R.P.; Small, J.L.; Jepsen, R.; Forshey, B.M.; Sanchez, J.L.; Smith, R.D.; Harris, R.; et al. Global Disease Outbreaks Associated with the 2015–2016 El Niño Event. *Sci. Rep.* **2019**, *9*, 1930. [CrossRef] [PubMed]

4. Morand, S.; Owers, K.A.; Waret-Szkuta, A.; McIntyre, K.M.; Baylis, M. Climate variability and outbreaks of infectious diseases in Europe. *Sci. Rep.* **2013**, *3*, 1774. [CrossRef] [PubMed]
5. Hanley, D.E.; Bourassa, M.A.; O'Brien, J.J.; Smith, S.R.; Spade, E.R. A Quantitative Evaluation of ENSO Indices. *J. Clim.* **2003**, *16*, 1249–1258. [CrossRef]

6. Kogan, F.N. Satellite-Observed Sensitivity of World Land Ecosystems to El Niño/La Niña. *Remote Sens. Environ.* **2000**, *74*, 445–462. [CrossRef]

7. Rifai, S.W.; Girardin, C.A.J.; Berenguer, E.; del Aguila-Pasquel, J.; Dahlsjö, C.A.L.; Doughty, C.E.; Jeffery, K.J.; Moore, S.; Oliveras, I.; Riutta, T.; et al. ENSO Drives interannual variation of forest woody growth across the tropics. *Philos. Trans. R. Soc. B Biol. Sci.* **2018**, *373*, 20170410. [CrossRef]

8. McPhaden, M.J.; Santoso, A.; Cai, W. *El Niño Southern Oscillation in a Changing Climate*; John Wiley & Sons: Hoboken, NJ, USA, 2020.

9. Iizumi, T.; Luo, J.-J.; Challinor, A.J.; Sakurai, G.; Yokozawa, M.; Sakuma, H.; Brown, M.; Yamagata, T. Impacts of El Niño Southern Oscillation on the global yields of major crops. *Nat. Commun.* **2014**, *5*, 3712. [CrossRef]

10. Kovats, R.S.; Bouma, M.J.; Hajat, S.; Worrall, E.; Haines, A. El Niño and health. *Lancet* **2003**, *362*, 1481–1489. [CrossRef]

11. Rodó, X.; Pascual, M.; Fuchs, G.; Faruque, A. ENSO and cholera: A nonstationary link related to climate change? *Proc. Natl. Acad. Sci. USA* **2002**, *99*, 12901–12906. [CrossRef]

12. Nicholls, N. A method for predicting murray valley encephalitis in southeast australia using the southern oscillation. *Aust. J. Exp. Biol. Med Sci.* **1986**, *64*, 587–594. [CrossRef]

13. Chaves, L.F.; Calzada, J.E.; Valderrama, A.; Saldana, A. Cutaneous leishmaniasis and sand fly fluctuations are associated with El Niño in Panamá. *PLoS Negl. Trop. Dis.* **2014**, *8*, e3210. [CrossRef] [PubMed]

14. Vincenti-Gonzalez, M.F.; Tami, A.; Lizarazo, E.; Grillet, M.E. ENSO-driven climate variability promotes periodic major outbreaks of dengue in Venezuela. *Sci. Rep.* **2018**, *8*, 5727. [CrossRef] [PubMed]

15. Pramanik, M.; Singh, P.; Kumar, G.; Ojha, V.P.; Dhiman, R.C. El Niño Southern Oscillation as an early warning tool for dengue outbreak in India. *BMC Public Health* **2020**, *20*, 1498. [CrossRef] [PubMed]

16. Caminade, C.; Turner, J.; Metelmann, S.; Hesson, J.C.; Blagrove, M.S.C.; Solomon, T.; Morse, A.P.; Baylis, M. Global risk model for vector-borne transmission of Zika virus reveals the role of El Niño 2015. *Proc. Natl. Acad. Sci. USA* **2017**, *114*, 119–124. [CrossRef] [PubMed]

17. Dhiman, R.C.; Sarkar, S. El Niño Southern Oscillation as an early warning tool for malaria outbreaks in India. *Malar. J.* **2017**, *16*, 122. [CrossRef]

18. Giles, J.R.; Eby, P.; Parry, H.; Peel, A.J.; Plowright, R.K.; Westcott, D.A.; McCallum, H. Environmental drivers of spatiotemporal foraging intensity in fruit bats and implications for Hendra virus ecology. *Sci. Rep.* **2018**, *8*, 9555. [CrossRef]

19. Daszak, P.; Zambrana-Torrelio, C.; Bogich, T.L.; Fernandez, M.; Epstein, J.H.; Murray, K.A.; Hamilton, H. Interdisciplinary approaches to understanding disease emergence: The past, present, and future drivers of Nipah virus emergence. *Proc. Natl. Acad. Sci. USA* **2012**, *110* (Suppl. 1), 3681–3688. [CrossRef]

20. Rojas, O.; Li, Y.; Cumani, R. Understanding the Drought Impact of El Niño on the Global Agricultural Areas: An Assessment Using FAO's Agricultural Stress Index (ASI); Food and Agriculture Organization of the United Nations (FAO): Rome, Italy, 2014.

21. Donges, J.F.; Schleussner, C.-F.; Siegmund, J.F.; Donner, R.V. Event coincidence analysis for quantifying statistical interrelationships between event time series. *Eur. Phys. J. Spéc. Top.* **2016**, *225*, 471–487. [CrossRef]

22. Corman, V.M.; Muth, D.; Niemeyer, D.; Drosten, C. Chapter Eight—Hosts and Sources of Endemic Human Coronaviruses. In *Advances in Virus Research 100*; Kielian, M., Mettenleiter, T.C., Roossinck, M.J., Eds.; Academic Press: Cambridge, MA, USA, 2018; pp. 163–188.

23. Rahman, M.; Chakraborty, A. Nipah virus outbreaks in Bangladesh: A deadly infectious disease. *WHO South-East Asia J. Public Health* **2012**, *1*, 208–212. [CrossRef]

24. Rahman, M.Z.; Islam, M.M.; Hossain, M.E.; Rahman, M.M.; Islam, A.; Siddika, A.; Sultana, S.; Klena, J.; Flora, M.; Daszak, P.; et al. Genetic diversity of Nipah virus in Bangladesh. *Int. J. Infect. Dis.* **2021**, *102*, 144–151. [CrossRef] [PubMed]

25. Ang, B.S.P.; Lim, T.C.C.; Wang, L.; Kraft, C.S. Nipah Virus Infection. *J. Clin. Microbiol.* **2018**, *56*, e01875-17. [CrossRef] [PubMed]

26. Field, H.; McCall, B.; Barrett, J. Australian Bat Lyssavirus Infection in a Captive Juvenile Black Flying Fox. *Emerg. Infect. Dis.* **1999**, *5*, 438–440. [CrossRef] [PubMed]

27. Chant, K.; Chan, R.; Smith, M.; Dwyer, D.E.; Kirkland, P. Probable human infection with a newly described virus in the family Paramyxoviridae. *Emerg. Infect. Dis.* **1998**, *4*, 273–275. [CrossRef]

28. Ge, X.-Y.; Li, J.-L.; Yang, X.-L.; Chmura, A.A.; Zhu, G.; Epstein, J.H.; Mazet, J.K.; Hu, B.; Zhang, W.; Peng, C.; et al. Isolation and characterization of a bat SARS-like coronavirus that uses the ACE2 receptor. *Nature* **2013**, *503*, 535–538. [CrossRef]

29. Chua, K.B.; Voon, K.; Crameri, G.; Tan, H.S.; Rosli, J.; McEachern, J.A.; Suluraju, S.; Yu, M.; Wang, L.-F. Identification and Characterization of a New Orthoreovirus from Patients with Acute Respiratory Infections. *PLoS ONE* **2008**, *3*, e3803. [CrossRef]

30. Zaki, A.M.; Van Boheemen, S.; Bestebroer, T.M.; Osterhaus, A.D.M.E.; Fouchier, R.A.M. Isolation of a Novel Coronavirus from a Man with Pneumonia in Saudi Arabia. *N. Engl. J. Med.* **2012**, *367*, 1814–1820. [CrossRef]

31. Ching, P.K.G.; de los Reyes, V.C.; Sucaldito, M.N.; Tayag, E.; Columna-Vingno, A.B.; Malbas, F.F.; Bolo, G.C.; Sejvar, J.J.; Eagles, D.; Playford, G.; et al. Outbreak of Henipavirus Infection, Philippines, 2014. *Emerg. Infect. Dis.* **2015**, *21*, 328–331. [CrossRef]

32. Gong, L.; Li, J.; Zhou, Q.; Xu, Z.; Chen, L.; Zhang, Y.; Xue, C.; Wen, Z.; Cao, Y. A New Bat-HKU2-like Coronavirus in Swine, China, 2017. *Emerg. Infect. Dis.* **2017**, *3*, 1607–1609. [CrossRef]
33. Zhu, N.; Zhang, D.; Wang, W.; Li, X.; Yang, B.; Song, J.; Zhao, X.; Huang, B.; Shi, W.; Lu, R.; et al. A Novel Coronavirus from Patients with Pneumonia in China, 2019. *N. Engl. J. Med.* **2020**, *382*, 727–733. [CrossRef]

34. Albers, S.; Campitelli, E. Rsoi: Import Various Northern and Southern Hemisphere Climate Indices. R Package Version 05. Available online: https://cran.r-project.org/web/packages/rsoi (accessed on 17 May 2022).

35. Bjornstad, O.N.; Cai, J. ncf: Spatial Covariance Functions. R Package Version 12-8. Available online: https://cran.r-project.org/web/packages/ncf (accessed on 17 May 2022).

36. Rösch, A.; Schmidbauer, H. WaveletComp 1.1: A Guided Tour through the R Package. Available online: http://www.hs-stat.com/projects/WaveletComp/WaveletComp_guided_tour.pdf (accessed on 17 May 2022).

37. Gouhier, T.; Grinsted, A.; Simko, V. biwavelet: Conduct Univariate and Bivariate Wavelet Analyses. R Package Version 02017. Available online: https://cran.r-project.org/web/packages/biwavelet/ (accessed on 17 May 2022).

38. R Core Team. *R: A language and Environment for Statistical Computing*; R Foundation for Statistical Computing: Vienna, Austria, 2020.

39. Venables, W.; Ripley, B.D. *Statistics Complements to Modern Applied Statistics with S*, 4th ed.; Springer: Berlin/Heidelberg, Germany, 2002.

40. Lefcheck, J.S. piecewiseSEM: Piecewise structural equation modelling in r for ecology, evolution, and systematics. *Methods Ecol. Evol.* **2016**, *7*, 573–579. [CrossRef]

41. Siegmund, J.F.; Sanders, T.G.; Heinrich, I.; Van der Maaten, E.; Simard, S.; Helle, G.; Donner, R. Meteorological Drivers of Extremes in Daily Stem Radius Variations of Beech, Oak, and Pine in Northeastern Germany: An Event Coincidence Analysis. *Front. Plant Sci.* **2016**, *7*, 733. [CrossRef] [PubMed]

42. Siegmund, J.F.; Siegmund, N.; Donner, R.V. CoinCalc—A new R package for quantifying simultaneities of event series. *Comput. Geosci.* **2017**, *98*, 64–72. [CrossRef]

43. Global Volcanism Program. Report on Mayon (Philippines). In *Bulletin of the Global Volcanism Network*; Wunderman, R., Ed.; Smithsonian Institution: Washington, DC, USA, 2013; Volume 38. [CrossRef]

44. Brierley, L.; Vonhof, M.J.; Olival, K.J.; Daszak, P.; Jones, K.E. Quantifying Global Drivers of Zoonotic Bat Viruses: A Process-Based Perspective. *Am. Nat.* **2016**, *187*, E53–E64. [CrossRef] [PubMed]

45. Olival, K.J.; Hosseini, P.R.; Zambrana-Torrelio, C.; Ross, N.; Bogich, T.L.; Daszak, P. Host and viral traits predict zoonotic spillover from mammals. *Nature* **2017**, *546*, 646–650. [CrossRef]

46. Letko, M.; Seifert, S.N.; Olival, K.J.; Plowright, R.K.; Munster, V.J. Bat-borne virus diversity, spillover and emergence. *Nat. Rev. Genet.* **2020**, *18*, 461–471. [CrossRef]

47. Latinne, A.; Hu, B.; Olival, K.J.; Zhu, G.; Zhang, L.; Li, H.; Chmura, A.A.; Field, H.E.; Zambrana-Torrelio, C.; Epstein, J.H.; et al. Origin and cross-species transmission of bat coronaviruses in China. *Nat. Commun.* **2020**, *11*, 4235. [CrossRef]

48. Wells, K.; Morand, S.; Wardeh, M.; Baylis, M. Distinct spread of DNA and RNA viruses among mammals amid prominent role of domestic species. *Glob. Ecol. Biogeogr.* **2020**, *29*, 470–481. [CrossRef]

49. Beyer, R.M.; Manica, A.; Mora, C. Shifts in global bat diversity suggest a possible role of climate change in the emergence of SARS-CoV-1 and SARS-CoV-2. *Sci. Total Environ.* **2021**, *767*, 145413. [CrossRef]

50. Buceta, J.; Johnson, K. Modeling the Ebola zoonotic dynamics: Interplay between enviroclimatic factors and bat ecology. *PLoS ONE* **2017**, *12*, e0179559. [CrossRef]

51. McFarlane, R.; Becker, N.; Field, H. Investigation of the Climatic and Environmental Context of Hendra Virus Spillover Events 1994–2010. *PLoS ONE* **2011**, *6*, e28374. [CrossRef] [PubMed]

52. Gurley, E.S.; Hegde, S.T.; Hossain, K.; Sazzad, H.M.S.; Hossain, M.J.; Rahman, M.; Sharker, M.A.Y.; Salje, H.; Islam, M.S.; Epstein, J.H.; et al. Convergence of Humans, Bats, Trees, and Culture in Nipah Virus Transmission, Bangladesh. *Emerg. Infect. Dis.* **2017**, *23*, 1446–1453. [CrossRef] [PubMed]

53. Chowdhury, S.; Khan, S.U.; Crameri, G.; Epstein, J.H.; Broder, C.C.; Islam, A.; Peel, A.J.; Barr, J.; Daszak, P.; Wang, L.-F.; et al. Serological Evidence of Henipavirus Exposure in Cattle, Goats and Pigs in Bangladesh. *PLOS Neglected Trop. Dis.* **2014**, *8*, e3302. [CrossRef] [PubMed]

54. Hegde, S.T.; Sazzad, H.M.S.; Hossain, M.J.; Alam, M.-U.; Kenah, E.; Daszak, P.; Rollin, P.; Rahman, M.; Luby, S.; Gurley, E.S. Investigating Rare Risk Factors for Nipah Virus in Bangladesh: 2001–2012. *EcoHealth* **2016**, *13*, 720–728. [CrossRef]

55. Cortes, M.C.; Cauchemez, S.; Lefrancq, N.; Luby, S.P.; Hossain, M.J.; Sazzad, H.; Rahman, M.; Daszak, P.; Salje, H.; Gurley, E.S. Characterization of the Spatial and Temporal Distribution of Nipah Virus Spillover Events in Bangladesh, 2007–2013. *J. Infect. Dis.* **2018**, *217*, 1390–1394. [CrossRef]

56. McKee, C.D.; Islam, A.; Luby, S.P.; Salje, H.; Hudson, P.J.; Plowright, R.K.; Gurley, E. The Ecology of Nipah Virus in Bangladesh: A Nexus of Land-Use Change and Opportunistic Feeding Behavior in Bats. *Viruses* **2021**, *13*, 169. [CrossRef]

57. Olival, K.J.; Latinne, A.; Islam, A.; Epstein, J.H.; Hersch, R.; Engstrand, R.C.; Gurley, E.S.; Amato, G.; Luby, S.P.; Daszak, P. Population genetics of fruit bat reservoir informs the dynamics, distribution and diversity of Nipah virus. *Mol. Ecol.* **2020**, *29*, 970–985. [CrossRef]

58. Epstein, J.H.; Anthony, S.J.; Islam, A.; Kilpatrick, A.M.; Khan, S.A.; Balkey, M.D.; Ross, N.; Smith, I.; Zambrana-Torrelio, C.; Tao, Y.; et al. Nipah virus dynamics in bats and implications for spillover to humans. *Proc. Natl. Acad. Sci. USA* **2020**, *117*, 29190–29201. [CrossRef]

59. Wahiduzzaman, M.; Luo, J.-J. A statistical analysis on the contribution of El Niño–Southern Oscillation to the rainfall and temperature over Bangladesh. *Meteorol. Atmos. Phys.* **2021**, *133*, 55–68. [CrossRef]
60. Cash, B.A.; Rodó, X.; Kinter, J.L.; Yunus, M. Disentangling the Impact of ENSO and Indian Ocean Variability on the Regional Climate of Bangladesh: Implications for Cholera Risk. *J. Clim.* **2010**, *23*, 2817–2831. [CrossRef]

61. Sharmin, S.; Glass, K.; Viennet, E.; Harley, D. Interaction of Mean Temperature and Daily Fluctuation Influences Dengue Incidence in Dhaka, Bangladesh. *PLOS Neglected Trop. Dis.* **2015**, *9*, e0003901. [CrossRef] [PubMed]

62. Chowdhury, M.R.; Ndiaye, O. Climate change and variability impacts on the forests of Bangladesh—A diagnostic discussion based on CMIP5 GCMs and ENSO. *Int. J. Clim.* **2017**, *37*, 4768–4782. [CrossRef]

63. Islam, A.R.M.T.; Salam, R.; Yeasmin, N.; Kamruzzaman, M.; Shahid, S.; Fattah, M.A.; Uddin, A.S.; Shahariar, M.H.; Mondol, A.H.; Jhajharia, D.; et al. Spatiotemporal distribution of drought and its possible associations with ENSO indices in Bangladesh. *Arab. J. Geosci.* **2021**, *14*, 2681. [CrossRef]

64. Plowright, R.K.; Peel, A.J.; Streicker, D.G.; Gilbert, A.T.; McCallum, H.; Wood, J.; Baker, M.; Restif, O. Transmission or Within-Host Dynamics Driving Pulses of Zoonotic Viruses in Reservoir-Host Populations. *PLOS Neglected Trop. Dis.* **2016**, *10*, e0004796. [CrossRef] [PubMed]

65. Páez, D.J.; Giles, J.; Mccallum, H.; Field, H.; Jordan, D.; Peel, A.J.; Plowright, R.K. Conditions affecting the timing and magnitude of Hendra virus shedding across pteropodid bat populations in Australia. *Epidemiol. Infect.* **2017**, *145*, 3143–3153. [CrossRef]

66. Becker, D.; Eby, P.; Madden, W.; Peel, A.; Plowright, R. Ecological conditions experienced by bat reservoir hosts predict the intensity of Hendra virus excretion over space and time. *bioRxiv* **2021**. [CrossRef]

67. Plowright, R.K.; Field, H.E.; Smith, C.; Divljan, A.; Palmer, C.; Tabor, G.; Daszak, P.; Foley, J.E. Reproduction and nutritional stress are risk factors for Hendra virus infection in little red flying foxes (*Pteropus scapulatus*). *Proc. R. Soc. B Boil. Sci.* **2008**, *275*, 861–869. [CrossRef]

68. McMichael, L.; Edson, D.; Smith, C.; Mayer, D.; Smith, I.; Kopp, S.; Meers, J.; Field, H. Physiological stress and Hendra virus in flying-foxes (*Pteropus* spp.), Australia. *PLoS ONE* **2017**, *12*, e0182171. [CrossRef]

69. Banerjee, A.; Baker, M.L.; Kulcsar, K.; Misra, V.; Plowright, R.; Mossman, K. Novel Insights into Immune Systems of Bats. *Front. Immunol.* **2020**, *11*, 26. [CrossRef]

70. Edson, D.; Peel, A.J.; Huth, L.; Mayer, D.G.; Vidgen, M.E.; McMichael, L.; Broos, A.; Melville, D.; Kristoffersen, J.; de Jong, C.; et al. Time of year, age class and body condition predict Hendra virus infection in Australian black flying foxes (*Pteropus alecto*). *Epidemiol. Infect.* **2019**, *147*, e240. [CrossRef]

71. Liqiang, H. El Nino makes trouble for broad areas of China, may hurt rice farmers. *China Daily*. 2019. Available online: http://www.chinadaily.com.cn/a/201903/01/WS5c7888eba3106c65c34ec1b7.html (accessed on 17 May 2022).

72. Tiwari, A.; Singh, S.; Soni, V.; Kumar, R. Environmental Impact of Recent Volcanic Eruption from Mt. Mayon Over 1 South-East Asia. *J. Geogr. Nat. Disast.* **2021**, *10*, 547.

73. Batumbo Boloweti, D.; Giraudoux, P.; Deniel, C.; Garnier, E.; Mauny, F.; Kasereka, C.M.; Kizungu, R.; Muyembe, J.J.; Bompangue, D.; Bornette, G. Volcanic activity controls cholera outbreaks in the East African Rift. *PLoS Negl. Trop. Dis.* **2020**, *14*, e0008406. [CrossRef] [PubMed]

74. Miranda, A.V.; Wiyono, L.; Rocha, I.C.N.; Cedeño, T.D.D.; Lucero-Prisno, D.E.I. Strengthening Virology Research in the Association of Southeast Asian Nations: Preparing for Future Pandemics. *Am. J. Trop. Med. Hyg.* **2021**, *105*, 1141–1143. [CrossRef] [PubMed]

75. Coker, R.J.; Hunter, B.M.; Rudge, J.W.; Liverani, M.; Hanvoravongchai, P. Emerging infectious diseases in southeast Asia: Regional challenges to control. *Lancet* **2011**, *377*, 599–609. [CrossRef]

76. Keesing, F.; Ostfeld, R.S. Impacts of biodiversity and biodiversity loss on zoonotic diseases. *Proc. Natl. Acad. Sci. USA* **2021**, *118*, e2023540118. [CrossRef]

77. Cai, W.; Ng, B.; Wang, G.; Santoso, A.; Wu, L.; Yang, K. Increased ENSO sea surface temperature variability under four IPCC emission scenarios. *Nat. Clim. Chang.* **2022**, *12*, 228–231. [CrossRef]

78. Hu, K.; Huang, G.; Huang, P.; Kosaka, Y.; Xie, S.-P. Intensification of El Niño-induced atmospheric anomalies under greenhouse warming. *Nat. Geosci.* **2021**, *14*, 377–382. [CrossRef]