# Longitudinal monitoring in Cambodia suggests higher circulation of alpha and betacoronaviruses in juvenile and immature bats of three species

Julien Cappelle, Neil Furey, Thavry Hoem, Tey Putita Ou, Thona Lim, Vibol

Hul, Oudam Heng, Véronique Chevalier, Philippe Dussart, Veasna Duong

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

## ► To cite this version:

Julien Cappelle, Neil Furey, Thavry Hoem, Tey Putita Ou, Thona Lim, et al.. Longitudinal monitoring in Cambodia suggests higher circulation of alpha and betacoronaviruses in juvenile and immature bats of three species. Scientific Reports, 2021, 11 (1), pp.24145. 10.1038/s41598-021-03169-z. hal-03521210

HAL Id: hal-03521210

https://hal.inrae.fr/hal-03521210

Submitted on 11 Jan 2022

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published 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
<figure><img src="image_3.png" /></figure>

OPEN

# Longitudinal monitoring in Cambodia suggests higher circulation of alpha and betacoronaviruses in juvenile and immature bats of three species

Julien Cappelle<fnref n="1" />,²✉, Neil Furey³,⁴, Thavry Hoem⁵, Tey Putita Ou⁵, Thona Lim⁶, Vibol Hul⁵,⁷, Oudam Heng⁵, Véronique Chevalier<fnref n="1" />,²,⁵, Philippe Dussart⁵,⁸ & Veasna Duong⁵

Recent studies suggest that coronaviruses circulate widely in Southeast Asian bat species and that the progenitors of the SARS-CoV-2 virus could have originated in rhinolophid bats in the region. Our objective was to assess the diversity and circulation patterns of coronavirus in several bat species in Southeast Asia. We undertook monthly live-capture sessions and sampling in Cambodia over 17 months to cover all phases of the annual reproduction cycle of bats and test specifically the association between their age and CoV infection status. We additionally examined current information on the reproductive phenology of *Rhinolophus* and other bat species presently known to occur in mainland southeast China, Vietnam, Laos and Cambodia. Results from our longitudinal monitoring (573 bats belonging to 8 species) showed an overall proportion of positive PCR tests for CoV of 4.2% (24/573) in cave-dwelling bats from Kampot and 4.75% (22/463) in flying-foxes from Kandal. Phylogenetic analysis showed that the PCR amplicon sequences of CoVs (n = 46) obtained clustered in *Alphacoronavirus* and *Betacoronavirus*. Interestingly, *Hipposideros larvatus* sensu lato harbored viruses from both genera. Our results suggest an association between positive detections of coronaviruses and juvenile and immature bats in Cambodia (OR = 3.24 [1.46–7.76], p = 0.005). Since the limited data presently available from literature review indicates that reproduction is largely synchronized among rhinolophid and hipposiderid bats in our study region, particularly in its more seasonal portions (above 16° N), this may lead to seasonal patterns in CoV circulation. Overall, our study suggests that surveillance of CoV in insectivorous bat species in Southeast Asia, including SARS-CoV-related coronaviruses in rhinolophid bats, could be targeted from June to October for species exhibiting high proportions of juveniles and immatures during these months. It also highlights the need to develop long-term longitudinal surveys of bats and improve our understanding of their ecology in the region, for both biodiversity conservation and public health reasons.

Emerging infectious diseases are a major threat to global health as illustrated by the recent emergence of the SARS-CoV-2 in China<fnref n="1" />. With more than 110 million human cases as of March 2021, the resulting disease, named COVID-19, led in a few months to the lock-down of billions of people around the world, with long-term political and economic impacts ongoing and difficult to estimate. Although the proximal origin of the virus is still unknown, phylogenetic studies suggest that progenitors of the virus could have originated in rhinolophid bats in Southeast Asia²,³. Hence, surveillance of *Rhinolophus* species, the putative reservoir of these progenitors,

<footnote n="1">¹CIRAD, UMR ASTRE (Animal, Santé, Territoires, Risques, Ecosystèmes), TA A 117/E, Campus International de Baillarguet, 34398 Montpellier CEDEX 5, France. ²ASTRE, CIRAD, INRAE, Univ Montpellier, Montpellier, France. ³Harrison Institute, Sevenoaks, Kent, England. ⁴Fauna and Flora International, Phnom Penh, Cambodia. ⁵Virology Unit, Institut Pasteur du Cambodge, Institut Pasteur International Network, Phnom Penh, Cambodia. ⁶Free the Bears, Phnom Penh, Cambodia. ⁷Unité des Virus Émergents (UVE: Aix-Marseille Univ-IRD 190-Inserm 1207), Marseille, France. ⁸Present address: Virology Unit, Institut Pasteur de Madagascar, Institut Pasteur International Network, Antananarivo, Madagascar. ☒email: julien.cappelle@cirad.fr</footnote>
could contribute to the development of prevention measures and early warning systems by identifying potentially zoonotic viruses circulating in bats and their typical circulation patterns.

The surveillance of other coronaviruses (CoVs), which are classified into four genera, *Alphacoronavirus*, *Beta-coronavirus*, *Gammacoronavirus* and *Deltacoronavirus*, may also be of interest⁴. CoVs can infect a wide variety of animals, causing diseases of varying severity in their respiratory, enteric, hepatic and neurological systems. Furthermore, viruses in the *Alphacoronavirus* and *Betacoronavirus* genera have caused disease in humans⁵. CoVs have been ubiquitously detected worldwide but only *Alphacoronavirus* and *Betacoronavirus* have been reported in bats⁴. Viruses in both genera show higher detections in faecal than oral samples, suggesting bat excretions may be more important in spill-over events⁵. In southeast Asia, both genera have been detected in a wide range of bat species in Cambodia, China, Laos, the Philippines, Taiwan and Thailand⁶⁻¹¹ and cross-species transmission among bat hosts appears to be common¹². In bats, reproduction phenology is linked to climate seasonality and has been identified as an important factor influencing pathogen circulation. When reproduction is synchronized, the resulting pulse of juveniles can lead to a sharp increase of immunologically naive individuals in a bat population, amplifying circulation of the pathogen¹³⁻¹⁵. Indeed, the number of annual birth pulses could impact the capacity of a pathogen to persist in a population as suggested by a modelling study for filoviruses¹⁶. Such seasonality related to bat reproduction cycles has been suggested for coronaviruses (CoVs) in several bat species in Africa and Europe. Juveniles in particular may play an important role in the circulation of these viruses¹⁷⁻²¹. In Southeast Asia, a large diversity of CoVs are known to circulate in bat populations, including SARS-related CoVs in horseshoe bats (*Rhinolophus* spp.)⁶,²². The impact of reproductive phenology on pathogen circulation has been examined in other pathogens such as Nipah virus but limited information is available regarding the ecology of CoVs in the region's bats²³,²⁴.

Our study was conducted before the COVID-19 pandemic and its main objectives were to assess the diversity of CoVs circulating in several bat species in two provinces of Cambodia and test the influence of reproductive phenology on levels of CoV circulation. We undertook monthly live-captures and sampling at two sites over 17 months between 2014 and 2016 to cover all phases of the annual reproduction cycle and specifically test the association between bat age and CoV infection status. We focused on cave-roosting bats in Kampot province and on tree-roosting fruit bats in Kandal province. Additionally, we summarize current knowledge on the reproductive phenology of *Rhinolophus* and other bat species in selected areas of Southeast Asia (mainland southeast China, Laos, Vietnam and Cambodia) to provide guidelines for targeting future surveillance and prevention efforts for CoVs that could potentially emerge in the region, including SARS-related CoVs.

## Material and methods

**Study sites.** The samples tested in this study originated from two sites in Cambodia where bats were previously captured and sampled to assess their diet and reproduction phenology (Fig. 1)²⁴⁻²⁶. The first site consisted of three caves (Bat Khteas, Phras Mea Kong Kea and Trai Lak) at Chhngauk hill in Kampot Province, which harbours populations of several cave-roosting bat species. Details regarding this site are provided by Thavry et al.²⁵ and Lim et al.²⁶. The second site was a roost of a single flying-fox species (*Pteropus lylei*) located in Kandal Province, for which further details are provided by Cappelle et al.²⁴.

**Bat capture and sampling.** Bat captures and sampling in Kampot Province were undertaken monthly from September 2014 to January 2016 (with the exception of October and December 2014 and January and May 2015 when the sampling team was not available). In each of the 13 sampling sessions, bats were captured using mist-nets over 3 consecutive nights (1 night per cave). Details on capture methods are provided by Thavry et al. and Lim et al.²⁵,²⁶. Bats were captured and handled in the field in accordance with guidelines approved by the American Society of Mammalogists²⁷, in addition to the requirements of the statutory study permission provided by the national authority responsible for wildlife research, the Forestry Administration of the Cambodian Ministry of Agriculture, Forest and Fisheries. Ethical approval for the study was also obtained from the Forestry Administration who participated in the field investigations and oversaw all aspects of the study. In Kandal Province, the flying-fox roost site, captures and sampling were implemented monthly from August 2015 to August 2016 with two additional sessions in May 2015 and December 2016. During each session, bats were captured using mist-nets deployed in the canopy of trees surrounding the roost over 2 consecutive nights.

All bats captured were examined to determine their sex, age and reproductive status following Anthony²⁸ and Racey²⁹. The reproductive phenology of bat species captured at the two study sites has been described previously²⁴,²⁶. The bats were grouped according to their age and sexual maturity and status as follows: females were classified as juvenile, nulliparous, parous, pregnant and/or lactating, whereas males were classified as juvenile, immature or mature.

**Coronavirus screening.** Rectal swabs were collected in the field from individual animals and placed in viral transport medium (VTM; containing tryptose phosphate Broth 2.95%, 145 mM of NaCl, 5% gelatin, 54 mM Amphotericin B, 106 U of penicillin-streptomycine per liter, 80 mg of gentamycine per liter [Sigma-Aldrich, Steinheim, Germany]) before immediate storage in liquid nitrogen for transport to the Institut Pasteur du Cambodge where they were stored at -80 °C. In the laboratory, viral RNA was extracted using Zymo Research Direct-zol RNA MiniPrep kit (Zymo Research, Cat # R2050, CA, USA) according to the manufacturer's instructions. The RNAs were screened with primary protocol by pan-coronaviridae RT-PCR targeting RdRp gene (adapted from Watanabe et al. (Watanabe-primers)³⁰ and suspected positive samples were confirmed with secondary RT-PCR adopted from Quan, et al. (Quan-primers)³¹. Reverse transcription (RT) was performed using SuperScript III First-Strand Synthesis Super-Mix (Invitrogen, San Diego, CA) according to the manufacturer's instructions.
<figure><img src="image_4.png" /><figcaption>**Figure 1.** Map of Cambodia showing the two study sites (generated with QGIS 3.18, www.qgis.org).</figcaption></figure>

The primary protocol PCR mixture (final volume: 25 µl) contained 2 µl of cDNA, 1 X PCR buffer, 0.20 mM (each) deoxynucleoside triphosphates (dNTPs), 1.5 mM of MgCl2, 0.2 µM of each primer (Appendix Table 1), and 2 U of Invitrogen Platinum TAQ DNA polymerase kit (Invitrogen, Carlsbad, CA, USA). The PCR mixture was incubated at 95 °C for 2 min, followed by 35 cycles of 94 °C for 30 s, 50 °C for 30 s and 72 °C for 30 s and finished with 72 °C for 5 min. The secondary protocol PCR mixture (final volume: 25 µl) contained 2 µl of cDNA, 1 X PCR buffer, 0.20 mM (each) deoxynucleoside triphosphates (dNTPs), 2 mM of MgCl2, 0.2 µM of each primer (Appendix Table 1), and 2 U of Invitrogen Platinum TAQ DNA polymerase kit (Invitrogen, Carlsbad, CA, USA). The PCR mixture was incubated at 95 °C for 5 min, followed by 15 cycles of 95 °C for 30 s, 65 °C for 30 s and 72 °C for 45 s, then 35 cycles of 94 °C for 30 s, 50 °C for 30 s and 72 °C for 45 s and finished with 72 °C for 7 min. Both nested PCR amplifications using round 2 forward and reverse primers (Appendix Table 1) were performed on 1 µl of the primary PCR products, using the same 1st round amplification conditions of each protocol. The amplification of sequences specific to CoVs was attested by the visualization of a 440 bp and a 434 bp fragments for modified Watanabe-primers and a 520 bp and a 328 bp for Quan-primers PCRs after the first and second PCR round, respectively. To limit the risk of contamination, RNA extraction, RT and PCR, nested-PCR and gel electrophoresis were carried out in separate rooms. In addition, negative controls (water) were included in each run of the nested-RT-PCR assay and results were validated only if these controls tested negative while the positive controls had to be positive. Positive PCR products were sequenced in both directions by direct Sanger sequencing at a commercial facility (Macrogen, Inc., Seoul, Republic of Korea).

**Sequence analysis.** Sequences were assembled and analyzed using CLC Genomics Workbench version 3.6.1 and BioEdit, version 7.0.9.1. Sequences were aligned with a representative set of CoV sequences including previously published CoVs in Cambodia from GenBank (Appendix Table 2) using MEGA version 6.06³². Phylogenetic trees based on RNA sequences were constructed using the Maximum Likelihood method with the GTR+G4+I model and bootstrap values were calculated after 1000 replicates. The best evolutionary model, phylogenetic and molecular evolutionary analyses were conducted using MEGA version 6.06³².

**Drivers of coronavirus positivity.** We employed a generalized linear mixed model (GLMM) to test the impact of several drivers on the probability of a bat being positive by PCR for CoV infection (lme4 package in R software³³). The result of the PCR test was the response variable with a binomial distribution. We used species, sex and age (classified in two categories, “mature adults” and “juveniles and immatures”) as explanatory variables with a fixed effect. To account for clustered samples collected during the same sampling session at the same
site, we included session identification code as a random effect to control for repeated measures from the same trapping session.

**Diversity and reproductive phenology of *Rhinolophus* and other bat species in mainland south-east China, Laos, Vietnam and Cambodia.** Because bats within the *Rhinolophus* genus are suspected to be the reservoir of the progenitors of both SARS-CoV-1 and SARS-CoV-2, we reviewed the literature to summarize current knowledge of their diversity and reproductive phenology in mainland southeast China (defined here as east of the Mekong River and south of the Yangtze River), Vietnam, Laos and Cambodia. Information on species occurrence was largely based on Burgin³⁴ which is the most recent and comprehensive treatment for the genus within the region, supplemented by verified records of previously undocumented rhinolophid species in Cambodia (Supplementary Materials). As Google Scholar has been found to provide superior coverage compared to major bibliographic databases³⁵, this was employed to assess the literature on reproductive phenology for our region between 2000 and 2021. The following keywords were used as search criteria: "Rhinolophus" AND "Reproductive phenology" AND "China" OR "Laos" OR "Vietnam" OR "Cambodia". Because information on bat reproductive phenology is rarely reported for our region, we also incorporated datasets from a variety of unpublished field studies. Since climate patterns have a major bearing on the reproductive phenology of insectivorous bats³⁶, we divided our region into two zones representing the two major climatic regimes present: (1) subtropics and seasonal tropics above 16° N in the north (southeast China, northern Vietnam and northern Laos) where temperatures and rainfall are noticeably seasonal (for example, in Bac Kan province [Vietnam] near the Vietnam-China border at 22° N, average monthly temperatures typically range between 12.3 and 27.5 °C and a wet summer season occurs from May to September, when ca. 75% of annual rainfall occurs³⁷), and (2) increasingly tropical below 16° N (southern Vietnam, southern Laos and Cambodia) where seasonal fluctuations in temperature progressively lessen and a wet season also occurs³⁸ (for instance, in Tay Ninh province [Vietnam] near the Vietnam-Cambodia border at 11° N, average monthly temperatures typically range between 25.2 and 28.9 °C, whereas ca. 85% of annual rainfall occurs from May to October³⁹).

The same process was repeated for genera that tested positive for coronavirus infections in our study (*Hipposideros*, *Eonycteris* and *Pteropus*), whereby their genus names were used instead of "Rhinolophus" as keywords in search criteria during their literature reviews (all other criteria remaining unchanged). As with the review for rhinolophid bats, datasets from unpublished field studies were incorporated in each instance.

**Results**

**Coronavirus screening.** A total of 573 bats belonging to 8 species were captured and sampled at the first site in Kampot province and 463 *Pteropus lylei* were captured and sampled at the second site in Kandal province (Appendix Table 3). Among the samples tested, 46 different rectal swabs tested positive for CoV, including 45 with Watanabe-primers and 21 with Quan-primers (meaning that 20 samples tested positive with both sets of primers and only 1 sample tested positive with Quan-primers only). The overall proportion of positive PCR tests for CoV was 4.2% (24/573) in the cave-dwelling bats from Kampot and 4.75% (22/463) in the flying-foxes from Kandal. Three species were dominant in the caves in Kampot: *Hipposideros larvatus* sensu lato (s.l.) with 277 individuals (11 PCR positive, 10 with modified Watanabe-primers and 1 with Quan-primers, 4.0%), *Eonycteris spelaea* with 128 individuals (13 PCR positive with modified Watanabe-primers, 10.2%) and *Taphozous melanopogon* with 110 individuals (no PCR positive). Only 22 individuals from a single *Rhinolophus* species were captured (*R. malayanus*, no PCR positive). Another 36 samples were collected from 4 other bat species (22 *Hipposideros armiger*, 8 *Lyroderma lyra*, 5 *Megaderma spasma* and 1 *Cynopterus* sp., no PCR positive).

Figure 2 shows monthly variation in coronavirus prevalence in the three most-captured species with at least one positive individual (detailed results are available in Appendix Tables 3, 4). For both *H. larvatus* s.l. and *E. spelaea*, prevalence was higher between July and November, i.e. 3–7 months after the birthing period. For both species, sexually immature individuals were predominant among the PCR-positive individuals: 76.9% (10/13) for *E. spelaea* and 72.7% (8/11) for *H. larvatus* s.l. For *P. lylei*, a peak in prevalence was observed in June 2016 with 17 positive individuals out of 45 sampled (37.8%), including 14 juveniles (82.4%), 2 other non-sexually mature individuals and only one sexually mature bat. An increase in prevalence was also observed in August and September of the preceding year.

**Drivers of coronavirus positivity.** Results of the GLMM confirmed the significantly higher positivity of juveniles and immatures individuals and did not show any significant effect of species (*Hipposideros*, *Eonycteris* and *Pteropus*, the three most-captured species with at least one positive individual) or sex on CoV positivity (Table 1, Appendix Fig. 2).

**Sequence analysis.** Phylogenetic analysis using sequences (n=45) detected with Watanabe-primers showed that CoVs detected in this study clustered in both alphacoronavirus and betacoronavirus groups (Fig. 3). CoVs detected in *H. larvatus* s.l. belonged to alphacoronavirus (n=6) in subcluster 4 along with CoV detected in *Hipposideros spp.* from Laos, Thailand and Hong Kong. Another 4 CoVs sequences from *H. larvatus* s.l. clustered in a separate branch from other SARS-CoV related viruses within group B of betacoronavirus. All 13 CoV sequences from *E. spelaea* grouped with other bat CoVs detected in Asia (group D) including Cambodia, China, Thailand and Singapore. Interestingly, *P. lylei* (n=22) harbored a unique and separate cluster with a CoV strain in *Cynopterus sphinx* from Thailand within betacoronavirus. Another phylogenetic tree using sequences generated with Quan-primers showed a similar clustering of Cambodian CoV in alpha and betacoronavirus (Appendix Fig. 1) and the additional CoV sequence from *H. larvatus* s.l. (MW507203) detected only with Quan-primers grouped with other sequences from the same bat species within alphacoronavirus.
<figure><img src="image_5.png" /><figcaption>**Figure 2.** Results of the longitudinal sampling of three bat species for coronaviruses. The first graph shows estimation of the coronavirus prevalence (with CI95) for *Pteropus lylei* at the second study site in Kandal province. The second and third graph show the same for *Hipposideros larvatus* s.l. and *Eonycteris spelaea* (respectively) at the first study site in Kampot province. The light grey areas show periods when pregnant females were captured, whereas dark grey areas show periods when lactating females were captured.</figcaption></figure>

<table><caption>**Table 1.** Results of the GLMM with the following references for the three explanatory variables: *E. spelaea* for species, female for sex and mature adults for age. OR odds ratio, CI95 95% confidence interval.</caption><thead><tr><td>Variables</td><td>Odd ratios</td><td>p value</td><td>OR CI95</td></tr></thead><tbody><tr><td>Intercept</td><td>0.014</td><td>&lt; 10<sup>-5</sup></td><td>0.00123-0.067</td></tr><tr><td>Species = H. larvatus</td><td>0.436</td><td>0.081</td><td>0.169-1.108</td></tr><tr><td>Species = P. lylei</td><td>0.117</td><td>0.053</td><td>0.007-0.929</td></tr><tr><td>Sex = male</td><td>1.566</td><td>0.209</td><td>0.778-3.18</td></tr><tr><td>Age = juveniles and immatures</td><td>3.245</td><td>0.005</td><td>1.462-7.758</td></tr></tbody></table>

**Diversity and reproductive phenology of Rhinolophus, Hipposideros, Eonycteris and Pteropus species.** Existing literature indicates that at least 27 rhinolophid and 16 hipposiderid species are currently documented in mainland southeast China, Vietnam, Laos and Cambodia, including 19 and six (respectively) in mainland southeast China where SARS-CoV and SARS-CoV-2 are thought to have emerged in humans (Supplementary Materials). A single species of *Eonycteris* (*E. spelaea*) occurs throughout this region, whereas three *Pteropus* species (*P. lylei*, *P. vampyrus*, *P. hypomelanus*) occur in its southerly portions (Cambodia and south Vietnam).

Our review of published and unpublished datasets (Supplementary Materials) reveals that basic information on reproductive phenology is lacking for most of these species and is summarized for the northern and southern portions of the region in Fig. 4. Overall, 16 studies and unpublished datasets including reproductive phenology data were obtained for rhinolophid bats (13 species), whereas 14 were obtained for hipposiderid bats (14 species)
<figure><img src="image_6.png" /><figcaption>**Figure 3.** Phylogenetic tree of (A) alphaCoV and (B) betaCoV partial RdRp gene with modified primers. The sequences detected in Cambodia are marked by red triangles for strains in the current study and black diamonds for strains previously described by Lacroix et al.⁶. The tree was built using the maximum likelihood method based on the GTR+G4+I model. The robustness of nodes was assessed with 1000 bootstrap replicates. Bootstrap values < 70 are not shown. The GenBank accession numbers of CoV detected in this study and included in the phylogenetic analysis are: MW507190-MW507207 and MW507209-MW507235.</figcaption></figure>
<figure><img src="image_7.png" /><figcaption>**Figure 3.** (continued)</figcaption></figure>

and four for *E. spelaea*. Aside from the present study, no data was available in the literature on the reproductive
phenology of the three *Pteropus* species within our study region.

Though data are presently only available for a limited number of species, these indicate that the reproductive
patterns of rhinolophid and hipposiderid bats in the northern zone (above 16° N) are characterized by restricted
seasonal monoestry (one litter within a 2-month period each year), whereby gestation occurs from March to May
and lactation largely from May to July, with species numbers for juveniles peaking in June–July. In contrast, the
pattern is less clear in the southern zone where even fewer data are presently available, but nonetheless suggests
<table><thead><tr><th colspan="16">Rhinolophus species</th></tr><tr><th>Northern Area</th><th>n (species)</th><th>Jan <sup>(2)</sup></th><th>Feb <sup>(1)</sup></th><th>Mar <sup>(3)</sup></th><th>Apr <sup>(3)</sup></th><th>May <sup>(7)</sup></th><th>Jun <sup>(5)</sup></th><th>Jul <sup>(3)</sup></th><th>Aug <sup>(4)</sup></th><th>Sep <sup>(2)</sup></th><th>Oct <sup>(3)</sup></th><th>Nov <sup>(2)</sup></th><th>Dec <sup>(3)</sup></th></tr></thead><tbody><tr><td>Gestation</td><td>10</td><td></td><td></td><td>6</td><td>7</td><td>6</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Lactation</td><td>12</td><td></td><td></td><td></td><td></td><td>8</td><td>10</td><td>5</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Volant Juveniles</td><td>9</td><td></td><td></td><td></td><td></td><td>3</td><td>6</td><td>5</td><td>1</td><td></td><td></td><td></td><td></td></tr><tr><th colspan="16">Hipposideros species</th></tr><tr><th>Northern Area</th><th>n (species)</th><th>Jan <sup>(2)</sup></th><th>Feb <sup>(1)</sup></th><th>Mar <sup>(3)</sup></th><th>Apr <sup>(3)</sup></th><th>May <sup>(6)</sup></th><th>Jun <sup>(4)</sup></th><th>Jul <sup>(3)</sup></th><th>Aug <sup>(4)</sup></th><th>Sep <sup>(2)</sup></th><th>Oct <sup>(3)</sup></th><th>Nov <sup>(2)</sup></th><th>Dec <sup>(3)</sup></th></tr><tr><td>Gestation</td><td>9</td><td></td><td></td><td>2</td><td>5</td><td>7</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Lactation</td><td>9</td><td></td><td></td><td></td><td></td><td>6</td><td>8</td><td>6</td><td>1</td><td>1</td><td></td><td></td><td></td></tr><tr><td>Volant Juveniles</td><td>3</td><td></td><td></td><td></td><td></td><td>2</td><td>2</td><td>3</td><td></td><td></td><td></td><td></td><td></td></tr><tr><th colspan="16">Eonycteris spelaea</th></tr><tr><th>Northern Area</th><th>n (species)</th><th>Jan <sup>(2)</sup></th><th>Feb <sup>(1)</sup></th><th>Mar <sup>(2)</sup></th><th>Apr <sup>(2)</sup></th><th>May <sup>(2)</sup></th><th>Jun <sup>(2)</sup></th><th>Jul <sup>(2)</sup></th><th>Aug <sup>(2)</sup></th><th>Sep <sup>(2)</sup></th><th>Oct <sup>(2)</sup></th><th>Nov <sup>(2)</sup></th><th>Dec <sup>(2)</sup></th></tr><tr><td>Gestation</td><td>1</td><td></td><td></td><td>1</td><td></td><td></td><td></td><td>1</td><td>1</td><td>1</td><td></td><td></td><td></td></tr><tr><td>Lactation</td><td>1</td><td></td><td></td><td></td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td></td></tr><tr><td>Volant Juveniles</td><td>1</td><td></td><td>1</td><td></td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td></tr><tr><th colspan="16">Southern Area</th></tr><tr><th>n (species)</th><th>Jan <sup>(3)</sup></th><th>Feb <sup>(2)</sup></th><th>Mar <sup>(5)</sup></th><th>Apr <sup>(3)</sup></th><th>May <sup>(3)</sup></th><th>Jun <sup>(4)</sup></th><th>Jul <sup>(3)</sup></th><th>Aug <sup>(3)</sup></th><th>Sep <sup>(3)</sup></th><th>Oct <sup>(3)</sup></th><th>Nov <sup>(3)</sup></th><th>Dec <sup>(3)</sup></th></tr><tr><td>Gestation</td><td>1</td><td>1</td><td>1</td><td></td><td>1</td><td>1</td><td></td><td>1</td><td>1</td><td></td><td>1</td><td>1</td><td></td></tr><tr><td>Lactation</td><td>1</td><td>1</td><td>1</td><td>1</td><td></td><td>1</td><td></td><td>1</td><td>1</td><td>1</td><td></td><td>1</td><td></td></tr><tr><td>Volant Juveniles</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td></td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td></tr></tbody></table>

**Figure 4.** Summary of available data on the reproductive phenology of *Rhinolophus*, *Hipposideros* and *Eonycteris* species in mainland southeast China, Vietnam, Laos and Cambodia. The chart shows the number of species for which a reproduction phase has been observed in a given month of the year for the northern (above 16° N) and southern (below 16° N) portions of the region. Each cell is shaded proportionally according to the number of species for which data were retrieved during the literature review, whereas superscript figures refer to the number of studies for which reproductive data were available for the genus in a given month.

that breeding is less seasonally constrained and may therefore qualify as extended seasonal monoestry (one litter per year within 2–7 months) or possibly even aseasonal monoestry (one litter at any time of the year) for some species. A rather different scenario applies to *E. spelaea*, whereby the species reproduces throughout the year. Within the northern zone, its reproductive patterns are characterized by seasonal bimodal polyoestry (two litters per year in close synchrony) with the majority of births occurring in March–April and August–September (and juveniles in May–July and September–November), whereas in the southern zone, the limited data currently available suggests that while births evidently occur throughout the year, these may peak in January and June–August.

Discussion

Our longitudinal monitoring of several bat species in Cambodia suggests that prevalence of different CoVs is higher in juvenile and immature individuals than sexually mature bats of three species (*E. spelaea*, *H. larvatus* s.l. and *P. lylei*) in Cambodia. The same finding has been observed in other regions and this may be attributable to the greater susceptibility of these immunologically naive individuals compared to adults<sup>13,17,18,21</sup>. When reproduction is highly seasonal and synchronized within a species, the pulses of juveniles entering the general population can lead to sharp seasonal increases in pathogen circulation as observed for Marburg virus in *Rousettus aegyptiacus* in Uganda<sup>14</sup>. Reproduction phenology can thus shape the circulation pattern of a pathogen and influence its capacity to persist in a population<sup>15,16</sup>.

While our results suggest a higher circulation of CoVs in juvenile and immature individuals, the limitation of our sampling to 17 months does not allow us to assess potential seasonality of CoVs in these bat populations. Furthermore, our results also suggest some level of variability in the occurrence of a seasonal peak between years, with a peak observed in juveniles in June 2016 in *P. lylei* but an increased prevalence among immatures in September 2015. Similarly, the increased CoV prevalence observed between July and November 2015 in both *H. larvatus* s.l. and *E. spelaea* were not observed in September and November of the preceding year, although this may be due to limited sampling in 2014. These inter-annual differences may be associated with variations in population-level immunity which in turn impact the level of maternal immunity transferred to juveniles<sup>40</sup>. High levels of maternal immunity might thus explain a several-month delay in the peak of CoV circulation observed in young bats, a period corresponding to the waning of maternal antibodies<sup>41</sup>. Irrespectively, multiple consecutive years of monitoring would be needed to confirm these potential variations in a general seasonal circulation
pattern of CoVs. Notwithstanding this, the association between CoV positivity and juvenile and immature individuals in our study suggests that surveillance of CoV could be targeted to the periods when juveniles and then sexually immature individuals form a higher proportion of bat populations.

The molecular characterization of CoVs detected during our study shows that several strains were circulating concomitantly in different bat species during the key period associated with an increased presence of young bats in the population. Hence, targeted sampling of coronaviruses during this period may help to survey the diversity of coronaviruses circulating in bat populations, including potentially zoonotic forms such as SARS-CoV-2 related viruses. Interestingly, bat species roosting in the same cave in Kampot Province did not share the same coronaviruses, which is consistent with previous associations between CoV clades and bat species observed in Cambodia⁶. As such, surveillance of CoVs should be conducted on all bat species present at a site, even when the different species roost together in close contact.

The primary message from our examination of the literature is that limited data are available on the reproductive phenology of rhinolophid and hipposiderid bats in southern China and Southeast Asia. More data should be collected to better understand the ecology of these potentially important species regarding the risk of future CoV emergences. Because these species and other insectivorous bats play important roles in Southeast Asian ecosystems and provide services that significantly benefit human populations in the region such as pest control and guano production⁴²,⁴³, improved knowledge of their ecology would also benefit public health and their conservation. Despite providing limited data, our literature review indicates that reproductive phenology is synchronized among *Rhinolophus* and *Hipposideros* species in the northern portion of the region (mainland southeast China, northern Laos and Vietnam above 16°N), with births peaking within a 2-month period in late spring (April-May). The same applies to many vespertilionid bat taxa that frequently share the same roost sites³⁷ and this period corresponds with the start of the summer rains in the region³⁸,³⁹ when insect biomass increases markedly⁴⁴,⁴⁵, providing optimal conditions for the survival and development of young bats. In contrast, the breeding season appears to be longer in the southern portion of the region (below 16°N in southern Laos, southern Vietnam and Cambodia), although births of some rhinolophid and hipposiderid taxa and other insectivorous cave bats still appear to peak in April-May²⁶,⁴³. A similar situation may also apply nearer to the equator, for instance in Peninsular Malaysia (1°-6°N), where although some *Rhinolophus* species evidently breed all year round, the majority mostly produce young in the first half of the year, with parturition in several cave-roosting rhinolophids also confined to April-May⁴⁶,⁴⁷. More broadly, a pattern of late dry season pregnancies and wet season production of young is common for insectivorous bats in many parts of the tropics worldwide⁴⁸⁻⁵⁰, although it is also possible that the months of parturition could vary depending on the seasonality and modality of rainfall patterns. As such, these factors should be accounted for when planning surveillance of rhinolophid, hipposiderid and other cave roosting insectivorous bats.

To conclude, our study highlights the benefits of undertaking longitudinal monitoring of bat populations, including species hosting SARS-related CoVs. Ideally, such programs should be permanent to allow the collection of long-term data and observation of inter-annual trends that may reveal years with higher prevalence of viruses in bats and a subsequent higher risk of spillover to humans⁵¹,⁵². Although more robust data are needed, the association observed between CoV positivity and juvenile and immature bats could be used to guide surveillance and maximize the probability of detecting more strains of CoV. In southern China, Vietnam, Laos and Cambodia and potentially elsewhere in Southeast Asia, the period from June to November could be prioritized for CoV surveillance in young rhinolophid and hipposiderid bats if the objective is to maximize the probability of detecting the circulation of these potentially zoonotic strains. Although crude, this first level of targeted surveillance may be particularly relevant when the financial resources available to conduct such programs are limited. Similarly, additional data on bat reproductive phenology will help to refine the strategies for the surveillance of CoV in Southeast Asia, including the surveillance of SARS-CoV-2 related viruses in rhinolophid bats.

Received: 15 April 2021; Accepted: 25 November 2021

Published online: 17 December 2021

## References

1. Wang, C., Horby, P. W., Hayden, F. G. & Gao, G. F. A novel coronavirus outbreak of global health concern. *Lancet* **395**, 470–473 (2020).

2. Andersen, K. G., Rambaut, A., Lipkin, W. I., Holmes, E. C. & Garry, R. F. The proximal origin of SARS-CoV-2. *Nat. Med.* **26**, 450–452 (2020).

3. Zhou, P. et al. A pneumonia outbreak associated with a new coronavirus of probable bat origin. *Nature* **579**, 270–273 (2020).

4. Woo, P. C. Y. et al. Discovery of seven novel mammalian and avian coronaviruses in the genus deltacoronavirus supports bat coronaviruses as the gene source of alphacoronavirus and betacoronavirus and avian coronaviruses as the gene source of gammapacoronavirus and deltacoronavirus. *J. Virol.* **86**, 3995–4008 (2012).

5. Wong, A. C. P., Li, X., Lau, S. K. P. & Woo, P. C. Y. Global epidemiology of bat coronaviruses. *Viruses* **11**, 174 (2019).

6. Lacroix, A. et al. Genetic diversity of coronaviruses in bats in Lao PDR and Cambodia. *Infect. Genet. Evol.* **48**, 10–18 (2017).

7. Tsuda, S. et al. Genomic and serological detection of bat coronavirus from bats in the Philippines. *Arch. Virol.* **157**, 2349–2355 (2012).

8. Han, Y. et al. Identification of diverse bat alphacoronaviruses and betacoronaviruses in China provides new insights into the evolution and origin of coronavirus-related diseases. *Front. Microbiol.* **10**, 20 (2019).

9. Xu, L. et al. Detection and characterization of diverse alpha- and betacoronaviruses from bats in China. *Virol. Sin.* **31**, 69–77 (2016).

10. Chen, Y.-N. et al. Detection of the severe acute respiratory syndrome-related coronavirus and alphacoronavirus in the bat population of Taiwan. *Zoonoses Public Health* **63**, 608–615 (2016).

11. Wacharapluesadee, S. et al. Group C betacoronavirus in bat guano fertilizer, Thailand. *Emerg. Infect. Dis.* **19**, 1349–1351 (2013).

12. Latinne, A. et al. Origin and cross-species transmission of bat coronaviruses in China. *Nat. Commun.* **11**, 4235 (2020).

13. Drexler, J. F. et al. Amplification of emerging viruses in a bat colony. *Emerg. Infect. Dis.* **17**, 449–456 (2011).
14. Amman, B. R. et al. Seasonal pulses of marburg virus circulation in Juvenile *Rousettus aegyptiacus* bats coincide with periods of increased risk of human infection. *PLoS Pathog.* **8**, 25 (2012).

15. Peel, A. J. et al. The effect of seasonal birth pulses on pathogen persistence in wild mammal populations. *Proc. R. Soc. Lond. B Biol. Sci.* **281**, 20132962 (2014).

16. Hayman, D. T. S. Biannual birth pulses allow filoviruses to persist in bat populations. *Proc. R. Soc. Lond. B Biol. Sci.* **282**, 20142591 (2015).

17. Gloza-Rausch, F. et al. Detection and prevalence patterns of Group I Coronaviruses in Bats, Northern Germany. *Emerg. Infect. Dis.* **14**, 626–631 (2008).

18. Annan, A. et al. Human betacoronavirus 2c EMC/2012-related viruses in bats, Ghana and Europe. *Emerg. Infect. Dis.* **19**, 456–459 (2013).

19. Anthony, S. J. et al. Global patterns in coronavirus diversity. *Virus Evol.* **3**, 25 (2017).

20. Montecino-Latorre, D. et al. Reproduction of East-African bats may guide risk mitigation for coronavirus spillover. *One Health Outlook* **2**, 2 (2020).

21. Maganga, G. D. et al. Genetic diversity and ecology of coronaviruses hosted by cave-dwelling bats in Gabon. *Sci. Rep.* **10**, 7314 (2020).

22. Hu, B. et al. Discovery of a rich gene pool of bat SARS-related coronaviruses provides new insights into the origin of SARS coronavirus. *PLoS Pathog.* **13**, e1006698 (2017).

23. Wacharapluesadee, S. et al. A longitudinal study of the prevalence of Nipah virus in *Pteropus lylei* bats in Thailand: Evidence for seasonal preference in disease transmission. *Vector-Borne Zoonot. Dis.* **10**, 183–190 (2010).

24. Cappelle, J. et al. Nipah virus circulation at human-bat interfaces, Cambodia. *Bull. World Health Organ.* **98**, 539–547 (2020).

25. Thavry, H., Cappelle, J., Bumrungsri, S., Thona, L. & Furey, N. M. The diet of the cave nectar bat (#Eonycteris spelaea# Dobson) suggests it pollinates economically and ecologically significant plants in Southern Cambodia. *Zool. Stud.* **56**, 25 (2017).

26. Lim, T., Cappelle, J., Hoem, T. & Furey, N. Insectivorous bat reproduction and human cave visitation in Cambodia: A perfect conservation storm?. *PLoS One* **13**, e0196554 (2018).

27. Sikes, R. S. 2016 Guidelines of the American Society of Mammalogists for the use of wild mammals in research and education. *J. Mammal.* **97**, 663–688 (2016).

28. Anthony, E. L. P. Age determination in bats. In *Ecological and Behavioral Methods for the Study of Bats* 47–58 (Smithsonian Press, 1988).

29. Racey, P. A. Reproductive assessment. In *Behavioural and Ecological Methods for the Study of Bats* 249–264 (Johns Hopkins University Press, 2009).

30. Watanabe, S. et al. Bat coronaviruses and experimental infection of bats, the Philippines. *Emerg. Infect. Dis.* **16**, 1217–1223 (2010).

31. Quan, P.-L. et al. Identification of a severe acute respiratory syndrome coronavirus-like virus in a leaf-nosed bat in Nigeria. *MBio* **1**, 25 (2010).

32. Tamura, K., Stecher, G., Peterson, D., Filipski, A. & Kumar, S. MEGA6: Molecular evolutionary genetics analysis version 6.0. *Mol. Biol. Evol.* **30**, 2725–2729 (2013).

33. Bates, D., Mächler, M., Bolker, B. & Walker, S. Fitting linear mixed-effects models using lme4. *J. Stat. Softw.* **67**, 1–48 (2015).

34. Burgin, C. Family Rhinolophidae, horseshoe bats. In *Handbook of the Mammals of the World* Vol. 9 260–332 (Lynx Edicions, 2019).

35. Martín-Martín, A., Orduna-Malea, E., Thelwall, M. & Delgado López-Cózar, E. Google Scholar, Web of Science, and Scopus: A systematic comparison of citations in 252 subject categories. *J. Informetr.* **12**, 1160–1177 (2018).

36. Racey, P. A. & Entwistle, E. Life history and reproductive strategies of bats. *Reprod. Biol. Bats* **20**, 363–468 (2000).

37. Furey, N. M., Mackie, I. J. & Racey, P. A. Reproductive phenology of bat assemblages in Vietnamese karst and its conservation implications. *Acta Chiropterol.* **13**, 341–354 (2011).

38. Sterling, E. J., Hurley, M. M. & Le, D. M. *Vietnam: A Natural History* (Yale University Press, 2006).

39. Van, N. K., Hzien, N. T., Loc, P. K. & Hiep, N. T. *Bioclimatic Diagrams of Vietnam* (Vietnam National University Publishing House, 2000).

40. Plowright, R. K. et al. Transmission or within-host dynamics driving pulses of zoonotic viruses in reservoir-host populations. *PLoS Negl. Trop. Dis.* **10**, 0004796 (2016).

41. Peel, A. J. et al. Support for viral persistence in bats from age-specific serology and models of maternal immunity. *Sci. Rep.* **8**, 3859 (2018).

42. Wanger, T. C., Darras, K., Bumrungsri, S., Tscharntke, T. & Klein, A.-M. Bat pest control contributes to food security in Thailand. *Biol. Conserv.* **171**, 220–223 (2014).

43. Furey, N. M., Racey, P. A., Ith, S., Touch, V. & Cappelle, J. Reproductive ecology of wrinkle-lipped free-tailed bats *Chaerephon plicatus* (Buchannan, 1800) in relation to Guano production in Cambodia. *Diversity* **10**, 91 (2018).

44. Ades, G. W. J. & Dudgeon, D. Insect seasonality in Hong Kong: A monsoonal environment in the northern tropics (1999).

45. Kai, K. H. & Corlett, R. T. Seasonality of forest invertebrates in Hong Kong, South China. *J. Trop. Ecol.* **18**, 637–644 (2002).

46. Kingston, T., Lim, B. L. & Zubaid, A. *Bats of Krau Wildlife Reserve* (Universiti Kebangsaan Malaysia, 2006).

47. Nurul-Ain, E., Rosli, H. & Kingston, T. Resource availability and roosting ecology shape reproductive phenology of rain forest insectivorous bats. *Biotropica* **49**, 382–394 (2017).

48. Fleming, T. H., Hooper, E. T. & Wilson, D. E. Three Central American Bat Communities: Structure, reproductive cycles, and movement patterns. *Ecology* **53**, 555–569 (1972).

49. Bernard, R. T. & Cumming, G. S. African bats: Evolution of reproductive patterns and delays. *Q. Rev. Biol.* **72**, 253–274 (1997).

50. Nguyen, S. T. et al. Bats (Chiroptera) of Bidoup Nui Ba National Park, Dalat Plateau, Vietnam. *Mammal Stud.* **46**, 53–68 (2021).

51. Plowright, R. K. et al. Urban habituation, ecological connectivity and epidemic dampening: The emergence of Hendra virus from flying foxes (*Pteropus* spp.). *Proc. R. Soc. B Biol. Sci.* **278**, 3703–3712 (2011).

52. Peel, A. J. et al. Synchronous shedding of multiple bat paramyxoviruses coincides with peak periods of Hendra virus spillover. *Emerg. Microbes Infect.* **8**, 1314–1323 (2019).

## Acknowledgements

We thank all the participants of Kandal and Kampot provinces for their involvement in the study. Particular gratitude is due to Dr Keo Omaliss of the Forestry Administration of the Cambodia Ministry of Agriculture, Forest and Fisheries for arranging study approval and permissions by the national authorities of Cambodia. This study was supported by the European Commission Innovate program (ComAcross project, Grant no. DCI-ASIE/2013/315-047) and the French ANR (ZooCov project, Flash COVID-19 call).

## Author contributions

J.C., N.F., V.C., P.D. and V.D. designed and supervised the implementation of the study. J.C., N.F., T.H. and T.L. coordinated the field work and collected biological samples. V.H., P.D. and D.V. supervised the biological analyses.
J.C. and N.F. designed and performed the literature review. T.P., V.H. and O.H. performed the biological analyses.
All the authors discussed the results and commented on the manuscript.

**Competing interests**

The authors declare no competing interests.

**Additional information**

**Supplementary Information** The online version contains supplementary material available at https://doi.org/10.1038/s41598-021-03169-z.

**Correspondence** and requests for materials should be addressed to J.C.

**Reprints and permissions information** is available at www.nature.com/reprints.

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

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

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2021