Next Article in Journal
A Novel Interaction Network Used by Potyviruses in Virus–Host Interactions at the Protein Level
Previous Article in Journal
Pathogenic Characteristics of a Porcine Astrovirus Strain Isolated in China
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Viral Diversity of Microbats within the South West Botanical Province of Western Australia

1
School of Veterinary Medicine, Murdoch University, Perth, WA 6150, Australia
2
Health and Biosecurity Business Unit, Australian Animal Health Laboratories, CSIRO, Geelong, VIC 3220, Australia
*
Author to whom correspondence should be addressed.
These authors contributed equally.
Viruses 2019, 11(12), 1157; https://doi.org/10.3390/v11121157
Submission received: 26 November 2019 / Revised: 6 December 2019 / Accepted: 6 December 2019 / Published: 13 December 2019
(This article belongs to the Section Animal Viruses)

Abstract

:
Bats are known reservoirs of a wide variety of viruses that rarely result in overt clinical disease in the bat host. However, anthropogenic influences on the landscape and climate can change species assemblages and interactions, as well as undermine host-resilience. The cumulative result is a disturbance of bat–pathogen dynamics, which facilitate spillover events to sympatric species, and may threaten bat communities already facing synergistic stressors through ecological change. Therefore, characterisation of viral pathogens in bat communities provides important basal information to monitor and predict the emergence of diseases relevant to conservation and public health. This study used targeted molecular techniques, serological assays and next generation sequencing to characterise adenoviruses, coronaviruses and paramyxoviruses from 11 species of insectivorous bats within the South West Botanical Province of Western Australia. Phylogenetic analysis indicated complex ecological interactions including virus–host associations, cross-species infections, and multiple viral strains circulating concurrently within selected bat populations. Additionally, we describe the entire coding sequences for five alphacoronaviruses (representing four putative new species), and one novel adenovirus. Results indicate that viral burden (both prevalence and richness) is not homogeneous among species, with Chalinolobus gouldii identified as a key epidemiological element within the studied communities.

1. Introduction

In the past three decades, viral surveillance in wild bat populations has accelerated due to human fatalities and the socio-economic impacts of emerging infectious diseases with bat origins, including respiratory syndrome coronaviruses, and paramyxoviruses such as Nipah and Hendra viruses [1,2]. Describing viral diversity within bat communities facilitates early detection of direct and indirect zoonotic disease risks [3], which are intensified by the anthropogenic modification of the landscape, climate change and increased human–domestic animal–wildlife interactions [4,5,6]. In bats, this link has been clearly illustrated by the loss of suitable habitat for flying foxes, and the risks of spillover events either through increased viral shedding [7], or movements of bat communities into bushland–urban interfaces [8].
Viral surveillance of wild bat populations is also necessary for proactive conservation management of this diverse taxonomic group. Viral shedding can act as an indicator of broader environmental change driving host population stress [9,10], which can influence disease dynamics through behavioural and physiological changes in the host [10]. Recent findings suggest that higher loads of coronavirus may hinder the survival of bats infected with Pseudogymnoascus destructans, the causative agent of White Nose Syndrome which has had a devastating impact on the conservation of a range of North American microbat species [11]. Further, a relationship between reduced nectar-based resources in winter, and Hendra virus spillover has been described [12], arguably mediated by impacts on the host energy budget and therefore immune system, as well as dispersed foraging patterns of the host species. Given the plausible environmental and biological drivers of viral shedding in bats, characterising their viral diversity is a crucial first step to develop strategies to detect changes in pathogen dynamics with implications for public health and conservation [13].
Among the variety of viral groups harboured by microbats [14], Adenoviridae, Coronaviridae, and Paramyxoviridae are diverse viral families infecting a wide range of vertebrates, and include species of known zoonotic potential, with key examples being the precursors to the Severe Acute Respiratory Syndrome and the Middle East Respiratory Syndrome coronaviruses which originated in insectivorous bats [15,16,17]. Even though adenoviruses of bat origin have not been associated with human disease outbreaks, in-vitro experiments have demonstrated their capacity to infect other vertebrate cell lines, including human cell lines [18]. Additionally, cross species infections of adenoviruses from New World monkeys to humans and consequent human to human transmission have also been reported [19]. Even though microbats are not known to harbour paramyxoviruses of zoonotic or conservation concern, this viral group is of interest as Hendra and Nipah viruses, two paramyxoviruses of flying fox origin, are able to infect humans with fatal disease progression [20,21].
To date, only a few viral groups have been detected and characterised in Australian bats, with surveillance predominantly focusing on the role of Australian flying foxes as reservoirs of Australian bat lyssavirus (ABLV) and Hendra virus [22,23,24]. ABLV has been isolated from a single Australian microbat species [25], with evidence of previous exposure in several other microbat and megabat species [26,27]. Coronavirus surveillance in northern Australia resulted in the detection of novel Australian alpha- and betacoronavirus strains [28], and transmission dynamics for coronaviruses have been modelled for the microbat, Myotis macropus [29]. Viral characterisation of populations of Miniopterus orianae bassanii and oceanensis showed a high prevalence of herpesviruses but neither coronaviruses, filoviruses, henipaviruses, nor lyssaviruses were detected [30]. Currently, the pathogen–host associations, and potential human risks or bat conservation implications of viral diversity within Australian insectivorous bats remain relatively unknown.
Here we describe the viral diversity of microbat communities within the South West Botanical Province (SWBP) of Western Australia, which is home to 13 species of microbats from two families, Vespertilionidae and Molossidae, including the locally endemic Falsistrellus mackenziei. We screened faecal samples for the presence of adenoviruses, coronaviruses and paramyxoviruses using targeted molecular approaches, and employed next generation sequencing (NGS) and serological protocols to enable a broader characterisation of viral diversity.

2. Materials and Methods

2.1. Study Area

The South West Botanical Province (SWBP), a global biodiversity hotspot [31], contains nine bioregions covering approximately 44 million hectares [32] and has a predominantly Mediterranean-type climate. Landscapes in the area are floristically diverse and have significantly changed due to anthropogenic activities [32]. Remnant natural cover on the east and west of the region is separated by the extensive monoculture known as the Western Australian wheatbelt (Figure 1).
Sampling took place to the south-west and the north-east of the province, with the 14 sampling sites falling on the east and west boundaries of the wheatbelt. The south-western sites were distributed across four bioregions: Swan Coastal Plain, Jarrah Forest, Warren and the Esperance Plains. The north-eastern sites were predominantly within the Avon bioregion (Figure 1).

2.2. Sample Collection

All sampling was approved by the Department of Biodiversity Conservation and Attractions, permits 08-001359-1, and CE005517. Capture, handling and sampling procedures were approved by the Murdoch University Animal Ethics Committee (R2882/16). Sampling occurred over two summers between November and April (2016–2018) using a cross-sectional sample design. The sites within the Avon bioregion were predominately sampled during the first year and the southwest sites during the second.
Bats were captured using harp traps and mist nets set up near water bodies or across forest tracks. Captured individuals were placed inside clean calico bags for at least ten minutes to allow time for defecation. Each individual was weighed and identified to species or genus level, using morphological features [33].
Collected faecal pellets were placed into 250 µL of RNAlater® (Ambion, Life Technologies, Carlsbad, CA, USA) and kept at either room temperature or 4 °C for a maximum of three days. Upon arrival to the laboratory, samples were stored at −80 °C until processing. A total of 10 µL of blood was collected, which was diluted 1:10 in phosphate buffered saline (PBS) and kept at 4 °C before processing

2.3. Molecular Analysis

Faecal samples were homogenised in the preservation solution by vigorous vortexing and centrifuged at 17,000 × g. A total of 50 µL of the supernatant was used as starting material for all extractions, using a Magmax viral RNA extraction kit (Ambion, Life Technologies, Vilnius, Lithuania) according to the manufacturer’s instructions.
Screening of coronaviruses and paramyxoviruses consisted of a nested and semi-nested reverse transcription PCR (RT-PCR) amplifying a segment of the RNA-dependent RNA polymerase (RdRp).
Coronavirus screening was performed with two primer sets that amplify different sections of the RdRp gene. These assays were run in parallel in order to increase sensitivity [34]. Optimisation of the paramyxovirus assay involved testing two previously published primer sets [35] on a sample subset (n= 48). The first set (pan-paramyxovirus) targeted all genera within the family and the second was specific to Morbillivirus, Respirovirus, and Henipavirus (Mor-Res-Hen primer set). Following an assessment of the results, subsequent screening was carried out using the group-specific primer set (Mor-Res-Hen primer set). Screening of adenoviruses involved the amplification of a section of the DNA polymerase gene (DNApol) (Table 1).
Reverse transcription PCR amplifications were performed using SuperScript III One-Step with Platinum Taq (Invitrogen, Thermofisher, Carlsbad, CA, USA). Additional PCRs were carried out with MyTaq PCR mix (Bioline, London, UK). Annealing temperatures are reported in Table 1.
Amplicons were visualised in 2% agarose gels stained with SYBR™ safe (Invitrogen, Thermofisher, Carlsbad, CA, USA). Products of the expected size were purified using Agencourt AMpure XP beads (Beckman Coulter, Brea, CA, USA). All fragments underwent Sanger sequencing in forward and reverse directions with the BigDye Terminator v3.1 Cycle Sequencing Kit (Applied Biosystems, Carlsbad, CA, USA) on an ABI PRISM 3130 DNA Analyser (Applied Biosystems, Carlsbad, CA, USA).
For full-length reconstruction of viral genomes, faecal nucleic acid samples were subject to sequence independent single primer amplification (SISPA) based PCR to capture RNA and DNA viral material [36]. Library preparation was performed using a Nextera XT library kit (Illumina, San Diego, CA, USA) as per the manufacturer’s instructions, and sequencing was performed on a NextSeq 500 using a 2x150 mid-output flowcell (Illumina, San Diego, CA, USA). De novo assembly of genomes was performed using SPAdes [37] with the –careful parameter enabled. Contigs were searched against the GenBank non-redundant protein database using DIAMOND v0.9 [38], and those returning significant similarity to viral sequences (E value ≤10-5) were retained. Contigs matching coronavirus and adenovirus sequences were extracted and examined manually by ORF prediction followed by BLASTn and BLASTp search. For contigs in which entire coding sequences were present, manual genome annotation was performed using Geneious v10.2.6 (BioMatters, Auckland, New Zealand).

2.4. Phylogenetic Analysis

Sequence read data were manually checked and trimmed in Geneious v10.2.6 (Biomatters, Auckland, New Zealand), and subjected to BLASTn search to verify the identity of the sequenced amplicons. Reference sequences, as recognised by the International Committee on Taxonomy of Viruses (ICTV) and other publicly available sequences, were retrieved from GenBank and aligned to the sequences generated by this study using MUSCLE [42] for the Adenoviridae and Coronaviridae alignments and MAFFT [43] for Paramyxoviridae.
Nucleotide-based phylogenetic relationships were inferred from maximum likelihood and bayesian analysis, based on the general time-reversible evolution model with gamma-distributed rate variation among sites and a proportion of invariable sites. Maximum likelihood trees were constructed in Mega X [44]; all gaps and missing data were subject to total deletion and trees bootstrapped 1000 times. Bayesian trees were constructed in MrBayes v.3.2 [45]. The default settings were employed to set the prior probabilities on the model parameters. Two runs were performed for 1,000,000 generations and convergence of the posterior probabilities were assessed by checking the standard deviation of split frequencies and the resulting PSRF statistics. Additional generations were added until the standard deviation reached a value between 0.01 and 0.03.
Phylogenetic analyses for five whole coronavirus genomes constructed from NGS data were performed on the polyprotein 1ab and spike gene amino acid sequences. Representative coronavirus sequences were retrieved from GenBank and aligned using MUSCLE [42]. Phylogenetic relationships were inferred by constructing maximum likelihood trees using the PROTGAMMAWAG model with 1000 bootstraps in RAxML [46].

2.5. Serology

Samples were tested for reactivity to four paramyxovirus and two betacoronavirus antigens in an indirect-binding Luminex® assay at a final working dilution of 1:50. All serological testing was performed at the CSIRO Australian Animal Health Laboratory (AAHL; Victoria, Australia). The assay coupled soluble G (sG) glycoproteins from Hendra virus (HeV), Nipah virus (NiV) and Cedar virus (CedV), and nucleoproteins (N) from Tioman virus (TioV), SARS, and MERS coronaviruses to individual microsphere sets as previously described [47].
A total of 662 samples were obtained for serological analysis. Those collected during the first season were pooled one in three (n= 246) or one in four (n= 24). Positive pools were then tested at the individual sample level where possible. Pooling has shown not to have a diluting or detrimental effect on the assay result [48]. All samples collected during the second season (n= 392) were tested individually. Median Fluorescence Intensity (MFI) was read using a Bio-Plex 200 instrument (Bio-Rad laboratories Hercules, CA, USA). Positive and non-template controls were included in all runs. Positive sera were of non-bat origin from naturally or experimentally infected animals. Serology results from this study are referred to as paramyxovirus-like (PaV), and betacoronavirus-like (βCoV).
As the Luminex assay has not been validated for microbat species due to the lack of species-specific positive controls, the MFI threshold to differentiate positive and negative samples was set up at 1000 MFI, which is at least three times the mean MFI of negative sera. It is based on previous studies published by the AAHL and elsewhere using the same Bio-Plex platform on sera from other bat species with values below 250 MFI considered negative [49,50,51,52].

2.6. Prevalence Estimates

The Wilson’s method [53] was used to estimate prevalence values and 95% confidence intervals in R package epitools [54]. Estimates for the molecular data were stratified by viral family, and host species. Seroprevalence was calculated using all samples within the serology dataset (n= 662). Positive results in pooled samples were assigned to a single individual within each pool.
Species-specific prevalence estimates were derived from samples tested individually and calculated for sample sizes greater than 20 individuals as estimates with a small denominator provide meaningless values. Differences between prevalence values were assessed using the Fisher Exact test as implemented in R.

3. Results

A total of 571 faecal samples were collected from 11 species of microbats, with Chalinolobus gouldii (n= 232) and Vespadelus regulus (n= 141) having the largest representation in the dataset. Details on number of faecal and blood samples per species and collection sites are shown in Table 2.

3.1. Serology

Resulting overall antibody prevalence estimates were generally low for all serology assays. The highest values were observed for SARS-CoV at 5.8% (95% CI 4.3–7.9%) and CedV at 5.3% (95% CI 3.8–7.3%), and the lowest for NiV at 1.3% (95% CI 0.7–2.6%) and MERS-CoV at 1.5% (95% CI 0.8–2.7%). Sample sizes, overall seroprevalance estimates and associated confidence intervals are shown in Table 3.
PaV-like and βCoV-like antibody responses were detected in eight and six species, respectively, (Table 4) but only C. gouldii and V. regulus reacted to all antigens. Species-specific antibody prevalence values were also low, with the highest value recorded for CeV antibodies in C. morio at 10.3% (95% CI 4.8–20.7%) and βCoV antibodies in N. gouldi at 10% (95% CI 5.1–20%).

3.2. Faecal PCR Analysis

Overall, the highest prevalence of viral nucleic acid shedding was observed for coronaviruses at 18% (95% CI 15–22%), followed by adenoviruses at 11% (95% CI 8–14%) and paramyxoviruses at 3% (95% CI 2–5%). All viral families were predominantly amplified from four species, C. gouldii, C. morio, N. geoffroyi, and V. regulus, with the former two species showing the highest prevalence estimates for adenovirus DNA and coronavirus RNA (CoV-RNA) (Table 5).
Prevalence values across species were significantly different for adenoviruses and coronaviruses but not for paramyxoviruses. This differentiation was driven by the high prevalence in Chalinolobus spp, compared to all other species.

3.3. Summary of Virus Prevalence in Relation to Serology

Seroprevalence values for βCoV antibodies markedly differed from CoV-RNA prevalence estimates, especially for C. gouldii, V. regulus and C. morio. These three species showed the greatest levels of CoV-RNA shedding with values ranging from 12% to 36% (Table 5), while the corresponding βCoV seoprevalence was only between 1.6% and 6.7% (Table 4). To the contrary, PaV-RNA prevalence rates for these three species were comparable to their observed seroprevalence ranges. Thus, C. morio, which showed the highest PaV-RNA prevalence at 9% (95% CI 3–21%), had a seroprevalence range between 3.3% and 10.3%. C. gouldii showed a PaV-RNA shedding prevalence of 4% (95% CI 2–7%) and a PaV-seroprevalence ranging between 3% to 6%, and in V. regulus PaV antibody prevalence varied between 0.6% to 6.1% and the corresponding PaV-RNA shedding was 1.4% (95% CI 0.4–5%). Individuals that were seropositive in either coronavirus or paramyxovirus assays yielded a negative PCR result for either family, except for a single C. morio which was RNA and antibody positive for coronaviruses.
The detection rate of co-infections was low and limited to Chalinolobus spp. CoV and PaV- RNA were found in three individuals, CoV-RNA and AdV-DNA were detected in 17 individuals, and all three viral families were amplified from a single sample.

3.4. Phylogenetic Analysis and Whole Genome Sequencing

Phylogenetic analysis showed a diverse number of coronavirus, adenovirus and paramyxovirus strains circulating within microbat populations of the SWBP. Two different strains of coronavirus were detected in C. gouldii, V. regulus and N. geoffroyi. For adenoviruses, three strains were observed in C. gouldii and two in N. geoffroyi, and two paramyxovirus strains were amplified in C. gouldii and V. regulus. GeneBank accession numbers of representative sequences generated by this study are reported in Table S1. Resulting tree topologies were similar under the maximum likelihood and Bayesian analyses. Additionally, a total of six full genomes were described: five coronaviruses and one adenovirus. Details are provided in the sections below.
In general, a host species–strain association independent of the sampling region was observed for all viral families. In the majority of cases, individuals captured between 50 to 600 km apart, shared the same viral strain. However, for C. morio, a geographic clustering pattern was observed in the coronavirus analysis, with cluster delineation appearing to be based on the south-west and north-east sampling sites.

3.4.1. Adenovirus

The final working alignment for adenovirus amplicons contained 61 sequences of 228 bp. All amplified strains grouped into five clades (I-V) within the genus Mastadenovirus (Figure 2).
Most clades were dominated by strains detected from individuals of the same species sampled at different locations. WA AdV clade I consisted of seven C. morio, clade III and IV of 26 and 12 C. gouldii respectively, and clade V of 11 V. regulus. Clade WA AdV II contained three highly similar sequences, amplified from two different species, N. geoffroyi, and N. gouldi.
Additionally, two separate sequences which did not cluster directly within the WA AdV clades, were detected from single C. gouldii and N. geoffroyi individuals. The former was most similar to those within WA AdV V and the latter, grouped with Bovine AdV 3 (U57334).
WA AdV clades IV, V and sequence 3206-Cg-Dry clustered with sequences amplified from other Vespertilionidae species, Bat AdV 2 (JN252129) and Bat AdV 205A (KX871230), a pattern not observed for the remaining sequences. Under the Bayesian phylogenetic analyses, the remaining clades grouped with sequences amplified from Rhinolophus sinicus (KT698853, NC029899, NC09902) a member of Rhinolophidae, while under the maximum likelihood approach there was not a clear association to any other sequence of bat origin.
A single full-length adenovirus genome was obtained from a C. gouldii (Bat AdV WA 3301 Cg Dwe, clade IV, GenBank accession MK472072). The genome was 37,617 nucleotides in length, with inverted terminal repeats of 288 bp and 293 bp at the 5’ and 3’ termini, respectively. A BLASTn search on the full-length genome returned most similarity to Human adenovirus D isolate (KF268332.1) with 72% identity over a 32% query cover. Of the 27 putative open reading frames identified, the conserved DNA polymerase protein showed the greatest identity (62%) to canine (AP000050.1), equine (ANG08548.1) and bat (YP_004782100.1) adenoviruses. The hexon protein appeared to be more evolutionarily conserved with the highest similarity to Bat mastadenovirus G (YP009325345.1) (query coverage 100%, identity 75%), identified from Corynorhinus rafinesquii, a member of Vespertilionidae.

3.4.2. Coronavirus

Coronavirus RNA was detected in 102 samples. Significant differences were observed in sensitivity between the two primer sets used for the screening. The CoV primers [40] detected 98 positives while the PQL primers [41] detected only nine. Since these two primer sets amplify different regions of the RdRp gene, only those sequences successfully amplified with the CoV primers were used to build the alignment, which consisted of 91 sequences of 396 bp.
Sequences grouped into four clades (I-IV) within the genus Alphacoronavirus (Figure 3). There was a clear host species dominance within each clade, WA CoV I and IV consisted predominantly of C. gouldii (n= 28, 24), WA CoV II of C. morio (n= 14), and WA CoV III of V. regulus (n= 16). Additionally, a single sequence was amplified from Ozimops sp. Despite the observed species-specific dominance, clades also included highly similar sequences detected from a minority of other species. Consequently, strains from two Nyctophilus spp, and one from V. regulus clustered with the 28 C. gouldii sequences within WA CoV I, two F. mackenziei, and two V. barvestocki clustered within the V. regulus clade (WA CoV III), and a sequence amplified from N. geoffroyi was identical to C. gouldii sequences within WA CoV IV.
Phylogenetic analysis also showed that all CoV groups were similar to sequences amplified from members of the same host family. Accordingly, the Ozimops sp sequence clustered with a Molossidae sequence (Chaerephon sp. HQ728486), and all the remaining clusters were most similar to sequences from members of Vespertilionidae. WA CoV clades I, II and III grouped with previously described strains from Australian and New Zealand vespertilionid bats (EU834951, KF545987-90).
Five full alphacoronavirus genomes were amplified (GenBank accession numbers MK472067–MK472071). Three genomes were obtained from C. gouldii individuals, 2028-Cg-CDR and 3301-Cg-Dwe, representing clade CoV I and Bat CoV WA 1087-Cg-MtG, within clade CoV IV. Additional genomes were amplified from a pooled faecal sample of V. regulus (Bat CoV WA Vrpool-Dwe) and one from Ozimops sp (Bat CoV WA 3607-Ozi-CDR). Genome sizes ranged from 27,405 nt to 28,171 nt in length, with all genomes containing the ORF1ab translated as a result of the presence of the ribosomal frameshift signal UUUAAAC. In addition to ORF1a, ORF1b and ORF1ab, all genomes contained coding regions for the spike protein, envelope protein, membrane glycoprotein, nucleocapsid and ORF3.
Phylogenetic analysis of the pp1ab amino acid sequences (Figure 4) and spike gene sequences (Figure S1) showed that genomes from clade CoV I clustered with the genome amplified from the pooled V. regulus faecal samples, while the Ozimops sp genome clustered with Rousettus bat coronavirus HKU10 (AFU92112). Even though these viruses predominantly clustered with other bat viruses, the genome representing clade CoV IV clustered with Porcine Epidemic Diarrhoea virus (PEDV) (ALB08471), with an amino acid identity of 74.1% for the pp1ab region and 67.7% for the spike protein.
Analysis of concatenated amino acid domains from ORF1ab (NSP3, NSP5, NSP12, NSP13, NSP14, NSP15, and NSP16) showed that the two samples from clade CoV I had over 99% identity to each other while levels of identity among the remaining three genomes ranged from 67.8% to 85.5%. Levels of identity to other coronaviruses also resulted in identity levels below 90%. Therefore, the five genomes reported here represent four putative new coronavirus species, as per the International Committee on the Taxonomy of Viruses guidelines.

3.4.3. Paramyxovirus

Paramyxovirus RNA was detected in 18 samples. The working alignment was 443 bp long and contained 17 sequences. Sequences clustered by host species, forming two main clades (Figure 5). WA PaV I consisted of nine strains amplified from C. gouldii, and WA PaV II consisted of four sequences from C. morio. These two clades and two additional sequences amplified from a single S. balstoni and a V. regulus clustered within the recently proposed genus Shaanvirus [55] and were phylogenetically related to other reported sequences amplified from host species of the same family. Two distinct sequences detected in V. regulus and C. gouldii clustered within the genus Morbillivirus. To further investigate the relationships to other bat morbilliviruses, these two samples were amplified with the Pan-paramyxovirus primer set. Only one sample was successfully sequenced, and the resulting sequence clustered within Shaanvirus instead of Morbillivirus.
No contigs with homology to known paramyxovirus species were found following next-generation sequencing analysis, inclusive of two samples paramyxovirus positive on RT-PCR (3197- Vr-Dry and 1087-Cg-MtG).

4. Discussion

Here we present a comprehensive characterisation of the viral diversity within Australian insectivorous bat communities, providing the first records of adenovirus and paramyxovirus shedding in Western Australian bats, as well as building on the current knowledge of coronaviruses in Australian microbats. Our results indicate that disease burden varies within bat species assemblages, and C. gouldii appears to be a key epidemiological element within the studied communities, sustaining the greatest viral richness and shedding. This study shows that WA microbat populations host a variety of viral families and sustain multiple viral strains concurrently, in line with other studies [56,57,58,59,60,61]. Our results are in contrast to the only additional Australian initiative to date aiming at characterizing viral diversity in microbats, where only herpesviruses were detected [30]. Difference in the discovery rate between these two studies may be due to the number of species surveyed (two subspecies vs. 11 species) as well as sample type (oral swabs vs. faeces), with oral swabs likely to carry lower levels of viral RNA/DNA, thus decreasing the sensitivity of molecular assays.
In general, phylogenetic analysis revealed an expected strain–host specificity independent of geographical distance, in line with previous findings for the targeted viral families [1,28,56,57,58,62,63,64,65,66,67,68]. At a finer scale, clustering patterns of coronavirus strains for C. gouldii and C. morio, may reflect the ecological connectivity between host populations. Host tropism over a wide spatial scale was particularly evident in the high similarity (99% identity) of the C. gouldii coronavirus genomes amplified from individuals captured at sites 350 km apart, while a partition between sites north and south of the wheatbelt was observed for coronavirus strains amplified from C. morio. Even though such partition was not reflected by either the adenovirus or paramyxovirus phylogenetic analysis, genetic or capture/recapture studies would clarify the current and historical connectivity of bat populations within the region, and their relationship to observed host–pathogen associations.
Phylogenetic analysis of coronavirus sequences showed that species-specific groups harboured a minority of identical or similar strains which were amplified from sympatric individuals of different genera, providing evidence for the occurrence of historic cross-species transmission events. For example, within clade WA CoV IV, a sequence amplified from N. geoffroyi was identical to those amplified from C. gouldii, with all individuals captured at the same location, while sequences amplified from F. mackenziei clustered with the clade dominated by V. regulus. Previous cross-species infections of bats have mainly been observed in cave-dwelling species that share roosting sites, with viral transmission plausible within this context [28]. However, bat species in the current study predominantly roost in tree hollows and decorticating bark, although cave colonies of C. morio are also known [33]. Given the predominant faecal–oral or respiratory transmission of CoV, the observed strain clustering patterns suggest that C. gouldii may either share roosts concurrently or use them alternatively with Nyctophilus species and V. regulus, and a similar interaction may occur between F. mackenziei and V. regulus. This is supported by studies showing that despite forest bats having strong intra and interspecific variation in roost requirements [69,70], some species also exhibit considerable flexibility in roost selection [71]. Additionally, unpublished records from the bat box monitoring program at Organ Pipes National Park (Victoria, Australia) document the occasional co-roosting of V. regulus with C. morio or C. gouldii [72].
Interestingly, the phylogenetic analysis of the spike protein from coronavirus genome Bat CoV 1087-Cg-MtG showed a higher level of similarity to Porcine Endemic Diarrhoea virus (PEDV) than to any other virus of bat origin. This raises the possibility of a spillover risk, as the spike protein determines cell tropism. For instance, the emergence of the alphacoronavirus Swine Acute Diarrhoea Syndrome (SADS) in Chinese swine, is a recent example of a cross-taxa infection likely originating from an insectivorous bat coronavirus [73]. Both SADS and HKU2 bat coronavirus, the likely strain of origin, share 86% identity across the spike protein gene. However, the nucleotide identity of only 67.7% between the spike protein of Bat CoV 1087-Cg- MtG and PEDV suggests that direct infection to swine would be unlikely, and intermediate host jumps would likely be required.
Cross-species infections were not generally supported by the paramyxovirus or adenovirus phylogenetic analysis. However, we detected an identical single adenovirus strain from N. geffroyi and N. gouldi collected at two separate regions. The small amplicon used for phylogenetic analysis may fail to discern a taxonomic separation between these two closely related species, although previous studies have reported the occurrence of divergent and cosmopolitan AdV strains [57,74,75]. However, evidence supporting this hypothesis has been based on a small amplicon and further studies on the co-evolution of adenoviruses and bats will benefit from whole genome analysis.
Host tropism was also supported by the full-length adenovirus amplified from C. gouldii. In general, the genome had low levels of identity to other available adenovirus genomes. A strong host–strain association was also supported by the significant difference of the fiber protein, which determines cellular tropism, to other sequences in GenBank, including those derived from bat adenoviruses. Additionally, identity analysis using the conserved hexon region was able to better identify this strain as of bat origin over DNApol or the fiber protein. This supports previous findings that tree topologies based on the hexon gene better reflect the phylogenetic relationships of the host than the DNApol gene [56,57,75].
Despite adenoviruses being detected in a smaller number of microbat species compared to coronaviruses (five versus eight species), they showed a similar level of strain richness, supporting the emergent consensus that this family encompasses a greater diversity than previously thought [56,68,76,77,78]. The relationship of the first three adenovirus clades to other published bat sequences was not resolved, likely due to the lack of comparable sequences. Despite the availability of several published Adenoviridae sequences, these represent alternative sections of the DNApol gene, highlighting the importance of unified protocols for the screening of wildlife populations to allow global comparisons [79].
Molecular detection of paramyxoviruses in faeces found a crude prevalence of 3% in line with previous reports [80,81]. We observed greater reactivity to antigens from Australian paramyxovirus species (Hendra virus and Cedar virus) than to those found in Asia (Nipah virus and Tioman virus). This pattern may be indicative of the closer relationship expected between Australian henipaviruses than to geographically distant viral species, and is consistent with cross-reactivity experiments on Hendra and Nipah antigens showing variability of the binding response despite the close relationship between these two viruses [47].
The serological results illustrate the challenges of sero-surveillance in novel species in the absence of a validated species-specific assay [82]. For example, RNA- and seroprevalence differed for the coronaviruses, in contrast to the paramyxoviruses where the RNA- and seroprevalence were consistent. A high prevalence of coronavirus RNA shedding suggests endemic infection of the studied communities, a finding not supported by the low seroprevalence. Such discrepancy is most likely caused by the lack of sensitivity of the serological assay at detecting alphacoronavirus antibodies, as it comprises SARS and MERS specific antigens. It could be argued that the use of alphacoronavirus antigens may have been more appropriate as a serological screening tool. However, given the number of betacornaviruses which are precursors of significant zoonotic viruses, we deemed it important to perform serological surveillance for this group of viruses. The serological assays employed here were developed for Australian flying foxes and have not been validated in insectivorous bats. Therefore, a true epidemiological picture cannot be resolved without further validation. Furthermore, unlike studies focusing on larger bats, the small amount of blood that may be ethically drawn from microbats makes it challenging to carry out confirmatory tests on seropositive samples.
Viral shedding, antibody reactivity, and viral richness varied across species, suggesting that viral burden is not homogeneous within bat communities, but harboured by a minority of key species. In this study, C. gouldii and C. morio consistently showed the highest levels of viral shedding across all viral families despite the marked differences in total sample size (232 vs. 45 individuals). Differential viral burden within bat communities has also been hypothesised for Australian pteropid bats as reservoirs of Hendra virus, with longitudinal studies indicating that spillover events are more likely associated with two of the four species of Australian flying foxes, despite all being known reservoirs of the virus [83].
While multiple studies discuss the life-traits that make bats competent viral reservoirs compared to other taxa [84,85,86], only a few have specifically looked into the drivers of pathogen persistence and richness within bat communities [87,88,89]. Such studies investigated the relationship between viral richness and geographic range, colony size, population fragmentation, body size, and threat status. In general, the role of home range appears to be an important driver of viral richness [88,89], which is congruent with the global pattern that widely distributed species at high densities harbour a greater diversity of pathogens [90,91]. Equally, phylogenetic affinity has been suggested as an underlying factor of viral richness, as pathogens are more likely to host-switch with closely related species [86,92]. Within this context, C. gouldii could be proposed as an important viral reservoir within bat communities of the SWBP. This species showed the greatest viral richness with a total of seven strains detected, it has a wide distribution being commonly found in a variety of habitats including natural, rural and urban environments [33], and it appears to have a higher level of tolerance to urbanisation than other Australian microbat species [93,94].
High viral prevalence and richness were expected from A. australis, as this species has a wide distribution [33], it is capable of covering long distances, and its suspected migratory behaviour [95] would facilitate contact with species with more limited ranges, thus making it a potential viral reservoir and super-spreader at large landscape scales. However, the small sample size (n= 9) precluded any in-depth characterization of the viral diversity within this group and its role as a super-spreader at large scales warrants further research.
Since temporal fluctuation of viral load has been documented in virus surveillance of north American and European microbats [96,97] as well as Hendra shedding in Australian peteropid bats [22,98], longitudinal studies would be required to fully understand the viral shedding and viral–host associations observed here. Furthermore, evidence of higher viral loads as indicators of individual and population scale stress [10] provides a good example of how long-term monitoring of viral prevalence can act as an indicator of environmental pressure on bat communities, particularly with regards to climate change where habitat contracture and distributional shifts may bring together populations harbouring novel viral strains. This approach would allow bat viral research to move beyond prediction of zoonotic disease emergence and be used as a management tool for monitoring the health of cryptic bat communities and the environment as a whole.

Supplementary Materials

The following are available online at https://www.mdpi.com/1999-4915/11/12/1157/s1, Table S1: GenBank accession numbers for representative sequences produced by this study, Table S2: GenBank accession numbers for outgroup sequences used in the construction of the adenovirus, coronavirus and paramyxovirus phylogenetic trees, Figure S1: Maximum likelihood phylogenetic analysis of the spike protein amino acid sequences derived from five coronavirus genomes.

Author Contributions

M.O., B.J. and D.P. designed the study and collected the data; V.B. and M.L.B. carried out all the serology laboratory work; D.P. performed all the molecular laboratory work and analyzed the data from the molecular and serology datasets; D.P. drafted the manuscript; M.O. and B.J., V.B. and M.L.B. reviewed and edited the manuscript; M.O. and B.J. supervised the research.

Funding

This research was funded by the Murdoch University Small Grants Scheme, Gunduwa Regional Conservation Association, the Holsworth Wildlife Research Endowment, The Australian Wildlife Society University Grants (2016), the Graduate Women (WA) Inc. Mary Walters Bursary, and the Alistair Bursary (2017).

Acknowledgments

We would like to acknowledge Christopher C. Broder and Lianying Yan from USU for the development of the antigens used in this study. We will also like to thank all the volunteers involved in the data collection. This project would not have been possible without their help and good disposition in the field. We would also like to thank Robert Bullen, Nicholas Dunlop, Terry Reardon, Andrew Grigg and Cameron Richardson from Alcoa, and members of the Department of Biodiversity, Conservation and Attractions in particular Sarah Comer, Janine Liddelow, Mark True and Peter Lacey for their assistance in identifying adequate bat trapping sites, which was fundamental to this project. We also thank the private landowners who provided access to their properties, as well as John and Lisa Lawson from the Lions Dryandra Village for facilitating our stay there. We would also like to thank the Australian Wildlife Conservancy and Bush Heritage Australia for supporting this project.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Drexler, J.F.; Corman, V.M.; Drosten, C. Ecology, evolution and classification of bat coronaviruses in the aftermath of SARS. Antiviral Res. 2014, 101, 45–56. [Google Scholar] [CrossRef] [PubMed]
  2. Marsh, G.A.; de Jong, C.; Barr, J.A.; Tachedjian, M.; Smith, C.; Middleton, D.; Yu, M.; Todd, S.; Foord, A.J.; Haring, V.; et al. Cedar Virus: A novel Henipavirus isolated from Australian bats. PLoS Pathog. 2012, 8, e1002836. [Google Scholar] [CrossRef] [PubMed]
  3. Carroll, D.; Daszak, P.; Wolfe, N.D.; Gao, G.F.; Morel, C.M.; Morzaria, S.; Pablos-Méndez, A.; Tomori, O.; Mazet, J.A.K. The Global Virome Project. Science (80-. ) 2018, 359, 872–874. [Google Scholar] [CrossRef] [PubMed]
  4. Bradley, C.A.; Altizer, S. Urbanization and the ecology of wildlife diseases. Trends Ecol. Evol. 2007, 22, 95–102. [Google Scholar] [CrossRef]
  5. Altizer, S.; Ostfeld, R.S.; Johnson, P.T.J.; Kutz, S.; Harvell, C.D. Climate change and infectious diseases: from evidence to a predictive framework. Science (80-. ) 2013, 341, 514–519. [Google Scholar] [CrossRef] [Green Version]
  6. Plowright, R.K.; Foley, P.; Field, H.E.; Dobson, A.P.; Foley, J.E.; Eby, P.; Daszak, P. Urban habituation, ecological connectivity and epidemic dampening: the emergence of Hendra virus from flying foxes (Pteropus spp.). Proceedings. Biol. Sci. 2011, 278, 3703–3712. [Google Scholar] [CrossRef] [Green Version]
  7. Walsh, M.G.; Wiethoelter, A.; Haseeb, M.A. The impact of human population pressure on flying fox niches and the potential consequences for Hendra virus spillover. Sci. Rep. 2017, 7, 8226. [Google Scholar] [CrossRef]
  8. Jung, K.; Threlfall, C.G. Urbanisation and Its Effects on Bats—A Global Meta-Analysis. In Bats in the Anthropocene: Conservation of Bats in a Changing World; Voigt, C.C., Kingston, T., Eds.; Springer International Publishing: London, UK, 2016; pp. 13–33. [Google Scholar]
  9. Roberts, K.E.; Hadfield, J.D.; Sharma, M.D.; Longdon, B. Changes in temperature alter the potential outcomes of virus host shifts. PLOS Pathog. 2018, 14, e1007185. [Google Scholar] [CrossRef] [Green Version]
  10. Gervasi, S.S.; Burgan, S.C.; Hofmeister, E.; Unnasch, T.R.; Martin, L.B. Stress hormones predict a host superspreader phenotype in the West Nile virus system. Proc. R. Soc. B Biol. Sci. 2017, 284, 20171090. [Google Scholar] [CrossRef]
  11. Davy, C.M.; Donaldson, M.E.; Subudhi, S.; Rapin, N.; Warnecke, L.; Turner, J.M.; Bollinger, T.K.; Kyle, C.J.; Dorville, N.A.S.-Y.; Kunkel, E.L.; et al. White-nose syndrome is associated with increased replication of a naturally persisting coronaviruses in bats. Sci. Rep. 2018, 8, 15508. [Google Scholar] [CrossRef]
  12. 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, 1–19. [Google Scholar] [CrossRef] [PubMed]
  13. Tompkins, D.M.; Carver, S.; Jones, M.E.; Krkošek, M.; Skerratt, L.F. Emerging infectious diseases of wildlife: A critical perspective. Trends Parasitol. 2015, 31, 149–159. [Google Scholar] [CrossRef] [PubMed]
  14. Hayman, D.T.S. Bats as viral reservoirs. Annu. Rev. Virol. 2016, 3, 77–99. [Google Scholar] [CrossRef] [PubMed]
  15. Memish, Z.A.; Mishra, N.; Olival, K.J.; Fagbo, S.F.; Kapoor, V.; Epstein, J.H.; AlHakeem, R.; Durosinloun, A.; Al Asmari, M.; Islam, A.; et al. Middle East Respiratory Syndrome Coronavirus in bats, Saudi Arabia. Emerg. Infect. Dis. 2013, 19, 1819–1824. [Google Scholar] [CrossRef] [Green Version]
  16. Lau, S.K.P.; Woo, P.C.Y.; Li, K.S.M.; Huang, Y.; Tsoi, H.-W.; Wong, B.H.L.; Wong, S.S.Y.; Leung, S.-Y.; Chan, K.-H.; Yuen, K.-Y. Severe acute respiratory syndrome coronavirus-like virus in Chinese horseshoe bats. Proc. Natl. Acad. Sci. USA 2005, 102, 14040–14045. [Google Scholar] [CrossRef] [Green Version]
  17. Reusken, C.B.; Haagmans, B.L.; Müller, M.A.; Gutierrez, C.; Godeke, G.-J.; Meyer, B.; Muth, D.; Raj, V.S.; Vries, L.S.-D.; Corman, V.M.; et al. Middle East respiratory syndrome coronavirus neutralising serum antibodies in dromedary camels: a comparative serological study. Lancet Infect. Dis. 2013, 13, 859–866. [Google Scholar] [CrossRef] [Green Version]
  18. Li, Y.; Ge, X.; Zhang, H.; Zhou, P.; Zhu, Y.; Zhang, Y.; Yuan, J.; Wang, L.-F.; Shi, Z. Host range, prevalence, and genetic diversity of adenoviruses in bats. J. Virol. 2010, 84, 3889–3897. [Google Scholar] [CrossRef] [Green Version]
  19. Chen, E.C.; Yagi, S.; Kelly, K.R.; Mendoza, S.P.; Maninger, N.; Rosenthal, A.; Spinner, A.; Bales, K.L.; Schnurr, D.P.; Lerche, N.W.; et al. Cross-species transmission of a novel adenovirus associated with a fulminant pneumonia outbreak in a new world monkey colony. PLoS Pathog. 2011, 7, e1002155. [Google Scholar] [CrossRef]
  20. Mahalingam, S.; Herrero, L.J.; Playford, E.G.; Spann, K.; Herring, B.; Rolph, M.S.; Middleton, D.; McCall, B.; Field, H.; Wang, L.-F. Hendra virus: an emerging paramyxovirus in Australia. Lancet Infect. Dis. 2012, 12, 799–807. [Google Scholar] [CrossRef]
  21. Mohd Nor, M.N.; Gan, C.H.; Ong, B.L. Nipah virus infection of pigs in peninsular Malaysia. Rev. Sci. Tech. 2000, 19, 160–165. [Google Scholar] [CrossRef]
  22. Field, H.; Jordan, D.; Edson, D.; Morris, S.; Melville, D.; Parry-Jones, K.; Broos, A.; Divljan, A.; McMichael, L.; Davis, R.; et al. Spatiotemporal aspects of hendra virus infection in pteropid bats (Flying-Foxes) in Eastern Australia. PLoS ONE 2015, 10, 1–14. [Google Scholar] [CrossRef] [PubMed]
  23. Mackenzie, J.S.; Field, H.E.; Guyatt, K.J. Managing emerging diseases borne by fruit bats (flying foxes), with particular reference to henipaviruses and Australian bat lyssavirus. J. Appl. Microbiol. 2003, 94, 59–69. [Google Scholar] [CrossRef] [PubMed]
  24. Ng, J.H.J.; Baker, M.L. Bats and bat-borne diseases: a perspective on Australian megabats. Aust. J. Zool. 2013, 61, 48. [Google Scholar] [CrossRef]
  25. Gould, A.R.; Kattenbelt, J.A.; Gumley, S.G.; Lunt, R.A. Characterisation of an Australian bat lyssavirus variant isolated from an insectivorous bat. Virus Res. 2002, 89, 1–28. [Google Scholar] [CrossRef]
  26. Field, H.E. Evidence of Australian bat lyssavirus infection in diverse Australian bat taxa. Zoonoses Public Health 2018, 65, 742–748. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Prada, D.; Boyd, V.; Baker, M.; Jackson, B.; O’Dea, M. Insights into Australian Bat Lyssavirus in insectivorous bats of Western Australia. Trop. Med. Infect. Dis. 2019, 4, 46. [Google Scholar] [CrossRef] [Green Version]
  28. Smith, C.S.; De Jong, C.E.; Meers, J.; Henning, J.; Wang, L.-F.; Field, H.E. Coronavirus infection and diversity in bats in the Australasian region. Ecohealth 2016, 13, 72–82. [Google Scholar] [CrossRef] [Green Version]
  29. Jeong, J.; Smith, C.S.; Peel, A.J.; Plowright, R.K.; Kerlin, D.H.; McBroom, J.; McCallum, H. Persistent infections support maintenance of a coronavirus in a population of Australian bats (Myotis macropus). Epidemiol. Infect. 2017, 145, 2053–2061. [Google Scholar] [CrossRef] [Green Version]
  30. Holz, P.H.; Lumsden, L.F.; Druce, J.; Legione, A.R.; Vaz, P.; Devlin, J.M.; Hufschmid, J. Virus survey in populations of two subspecies of bent-winged bats (Miniopterus orianae bassanii and oceanensis) in south-eastern Australia reveals a high prevalence of diverse herpesviruses. PLoS ONE 2018, 13, e0197625. [Google Scholar] [CrossRef]
  31. Myers, N.; Mittermeier, R.A.; Mittermeier, C.G.; da Fonseca, G.A.; Kent, J. Biodiversity hotspots for conservation priorities. Nature 2000, 403, 853–858. [Google Scholar] [CrossRef]
  32. Wardell-Johnson, G.; Wardell-Johnson, A.; Bradby, K.; Robinson, T.; Bateman, P.W.; Williams, K.; Keesing, A.; Braun, K.; Beckerling, J.; Burbridge, M. Application of a Gondwanan perspective to restore ecological integrity in the south-western Australian global biodiversity hotspot. Restor. Ecol. 2016, 1–11. [Google Scholar] [CrossRef]
  33. Churchill, S. Australian bats, 2nd ed.; Allen & Unwin: Crows Nest, Australia, 2008. [Google Scholar]
  34. Goldstein, T.; University of California, Davis, CA, USA. Personal communication, 2017.
  35. Tong, S.; Chern, S.-W.W.; Li, Y.; Pallansch, M.A.; Anderson, L.J. Sensitive and broadly reactive reverse transcription-PCR assays to detect novel paramyxoviruses. J. Clin. Microbiol. 2008, 46, 2652–2658. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. O’Dea, M.A.; Jackson, B.; Jackson, C.; Xavier, P.; Warren, K. Discovery and partial genomic characterisation of a novel nidovirus associated with respiratory disease in wild shingleback lizards (Tiliqua rugosa). PLoS ONE 2016, 11. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Bankevich, A.; Nurk, S.; Antipov, D.; Gurevich, A.A.; Dvorkin, M.; Kulikov, A.S.; Lesin, V.M.; Nikolenko, S.I.; Pham, S.; Prjibelski, A.D.; et al. SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing. J. Comput. Biol. 2012, 19, 455–477. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Buchfink, B.; Xie, C.; Huson, D.H. Fast and sensitive protein alignment using DIAMOND. Nat. Methods 2015, 12, 59–60. [Google Scholar] [CrossRef] [PubMed]
  39. Roy, S.; Vandenberghe, L.H.; Kryazhimskiy, S.; Grant, R.; Calcedo, R.; Yuan, X.; Keough, M.; Sandhu, A.; Wang, Q.; Medina-Jaszek, C.A.; et al. Isolation and characterization of adenoviruses persistently shed from the gastrointestinal tract of non-human primates. PLoS Pathog. 2009, 5, e1000503. [Google Scholar] [CrossRef]
  40. Lacroix, A.; Duong, V.; Hul, V.; San, S.; Davun, H.; Omaliss, K.; Chea, S.; Hassanin, A.; Theppangna, W.; Silithammavong, S.; et al. Genetic diversity of coronaviruses in bats in Lao PDR and Cambodia. Infect. Genet. Evol. 2017, 48, 10–18. [Google Scholar] [CrossRef]
  41. Quan, P.-L.; Firth, C.; Street, C.; Henriquez, J.A.; Petrosov, A.; Tashmukhamedova, A.; Hutchison, S.K.; Egholm, M.; Osinubi, M.O.V.; Niezgoda, M.; et al. Identification of a severe acute respiratory syndrome coronavirus-like virus in a leaf-nosed bat in Nigeria. MBio 2010, 1, e00208-10. [Google Scholar] [CrossRef] [Green Version]
  42. Edgar, R.C. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004, 32, 1792–1797. [Google Scholar] [CrossRef] [Green Version]
  43. Katoh, K.; Misawa, K.; Kuma, K.; Miyata, T. MAFFT: a novel method for rapid multiple sequence alignment based on fast Fourier transform. Nucleic Acids Res. 2002, 30, 3059–3066. [Google Scholar] [CrossRef] [Green Version]
  44. Kumar, S.; Stecher, G.; Li, M.; Knyaz, C.; Tamura, K. MEGA X: Molecular evolutionary genetics analysis across computing platforms. Mol. Biol. Evol. 2018, 35, 1547–1549. [Google Scholar] [CrossRef]
  45. Ronquist, F.; Teslenko, M.; Van Der Mark, P.; Ayres, D.L.; Darling, A.; Höhna, S.; Larget, B.; Liu, L.; Suchard, M.A.; Huelsenbeck, J.P. Mrbayes 3.2: Efficient bayesian phylogenetic inference and model choice across a large model space. Syst. Biol. 2012, 61, 539–542. [Google Scholar] [CrossRef] [Green Version]
  46. Stamatakis, A. RAxML version 8: A tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics 2014, 30, 1312–1313. [Google Scholar] [CrossRef]
  47. Bossart, K.N.; McEachern, J.A.; Hickey, A.C.; Choudhry, V.; Dimitrov, D.S.; Eaton, B.T.; Wang, L.-F. Neutralization assays for differential henipavirus serology using Bio-Plex Protein Array Systems. J. Virol. Methods 2007, 142, 29–40. [Google Scholar] [CrossRef]
  48. Boyd, V.; Smith, I.; Crameri, G.; Burroughs, A.L.; Durr, P.a.; White, J.; Cowled, C.; Marsh, G.a.; Wang, L.F. Development of multiplexed bead arrays for the simultaneous detection of nucleic acid from multiple viruses in bat samples. J. Virol. Methods 2015, 223, 5–12. [Google Scholar] [CrossRef]
  49. Hayman, D.T.S.; Wang, L.F.; Barr, J.; Baker, K.S.; Suu-Ire, R.; Broder, C.C.; Cunningham, A.A.; Wood, J.L.N. Antibodies to henipavirus or henipa-like viruses in domestic pigs in Ghana, West Africa. PLoS ONE 2011, 6, e25256. [Google Scholar] [CrossRef] [Green Version]
  50. Hayman, D.T.S.; Suu-Ire, R.; Breed, A.C.; McEachern, J.A.; Wang, L.; Wood, J.L.N.; Cunningham, A.A. Evidence of henipavirus infection in West African fruit bats. PLoS ONE 2008, 3, e2739. [Google Scholar] [CrossRef]
  51. 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 Biol. Sci. 2008, 275, 861–869. [Google Scholar] [CrossRef] [Green Version]
  52. Breed, A.C.; Yu, M.; Barr, J.A.; Crameri, G.; Thalmann, C.M.; Wang, L.F. Prevalence of henipavirus and rubulavirus antibodies in pteropid bats, Papua New Guinea. Emerg. Infect. Dis. 2010, 16, 1997–1999. [Google Scholar] [CrossRef]
  53. Brown, L.D.; Cai, T.T.; Dasgupta, A. Interval estimation for a binomial proportion. Stat. Sci. 2001, 16, 101–133. [Google Scholar] [CrossRef]
  54. Aragon, T. epitools: Epidemiology Tools. R package version 0.5-10. Published 26 October 2017. Available online: https://CRAN.R-project.org/package=epitools (accessed on 20 April 2019).
  55. Wu, Z.; Yang, L.; Ren, X.; He, G.; Zhang, J.; Yang, J.; Qian, Z.; Dong, J.; Sun, L.; Zhu, Y.; et al. Deciphering the bat virome catalog to better understand the ecological diversity of bat viruses and the bat origin of emerging infectious diseases. ISME J. 2016, 10, 609–620. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  56. Iglesias-Caballero, M.; Juste, J.; Vázquez-Morón, S.; Falcon, A.; Aznar-Lopez, C.; Ibáñez, C.; Pozo, F.; Ruiz, G.; Berciano, J.M.; Garin, I.; et al. New adenovirus groups in western palaearctic bats. Viruses 2018, 10, 443. [Google Scholar] [CrossRef] [Green Version]
  57. Vidovszky, M.Z.; Kohl, C.; Boldogh, S.; Görföl, T.; Wibbelt, G.; Kurth, A.; Harrach, B. Random sampling of the Central European bat fauna reveals the existence of numerous hitherto unknown adenoviruses. Acta Vet. Hung. 2015, 63, 508–525. [Google Scholar] [CrossRef] [Green Version]
  58. Wray, A.K.; Olival, K.J.; Morán, D.; Lopez, M.R.; Alvarez, D.; Navarrete-Macias, I.; Liang, E.; Simmons, N.B.; Lipkin, W.I.; Daszak, P.; et al. Viral diversity, prey preference, and Bartonella prevalence in Desmodus rotundus in Guatemala. Ecohealth 2016, 13, 761–774. [Google Scholar] [CrossRef]
  59. Chu, D.K.W.; Poon, L.L.M.; Chan, K.H.; Chen, H.; Guan, Y.; Yuen, K.Y.; Peiris, J.S.M. Coronaviruses in bent-winged bats (Miniopterus spp.). J. Gen. Virol. 2006, 87, 2461–2466. [Google Scholar] [CrossRef]
  60. Dominguez, S.R.; O’Shea, T.J.; Oko, L.M.; Holmes, K.V. Detection of group 1 coronaviruses in bats in North America. Emerg. Infect. Dis. 2007, 13, 1295–1300. [Google Scholar] [CrossRef]
  61. Moreira-Soto, A.; Taylor-Castillo, L.; Vargas-Vargas, N.; Rodríguez-Herrera, B.; Jiménez, C.; Corrales-Aguilar, E. Neotropical bats from Costa Rica harbour diverse coronaviruses. Zoonoses Public Health 2015, 62, 501–505. [Google Scholar] [CrossRef]
  62. Wilkinson, D.A.; Temmam, S.; Lebarbenchon, C.; Lagadec, E.; Chotte, J.; Guillebaud, J.; Ramasindrazana, B.; Héraud, J.-M.; de Lamballerie, X.; Goodman, S.M.; et al. Identification of novel paramyxoviruses in insectivorous bats of the Southwest Indian Ocean. Virus Res. 2012, 170, 159–163. [Google Scholar] [CrossRef]
  63. Yuan, L.; Li, M.; Li, L.; Monagin, C.; Chmura, A.A.; Schneider, B.S.; Epstein, J.H.; Mei, X.; Shi, Z.; Daszak, P.; et al. Evidence for retrovirus and paramyxovirus infection of multiple bat species in china. Viruses 2014, 6, 2138–2154. [Google Scholar] [CrossRef]
  64. Gloza-Rausch, F.; Ipsen, A.; Seebens, A.; Göttsche, M.; Panning, M.; Drexler, J.F.; Petersen, N.; Annan, A.; Grywna, K.; Müller, M.; et al. Detection and prevalence patterns of group I coronaviruses in bats, northern Germany. Emerg. Infect. Dis. 2008, 14, 626–631. [Google Scholar] [CrossRef]
  65. Rizzo, F.; Edenborough, K.M.; Toffoli, R.; Culasso, P.; Zoppi, S.; Dondo, A.; Robetto, S.; Rosati, S.; Lander, A.; Kurth, A.; et al. Coronavirus and paramyxovirus in bats from Northwest Italy. BMC Vet. Res. 2017, 13, 396. [Google Scholar] [CrossRef] [Green Version]
  66. Asano, K.M.; Hora, A.S.; Scheffer, K.C.; Fahl, W.O.; Iamamoto, K.; Mori, E.; Brandão, P.E. Alphacoronavirus in urban Molossidae and Phyllostomidae bats, Brazil. Virol. J. 2016, 13, 110. [Google Scholar] [CrossRef] [Green Version]
  67. Zheng, X.Y.; Qiu, M.; Chen, H.F.; Chen, S.W.; Xiao, J.P.; Jiang, L.N.; Huo, S.T.; Shi, T.L.; Ma, L.Z.; Liu, S.; et al. Molecular detection and phylogenetic characterization of bat and human adenoviruses in Southern China. Vector-Borne Zoonotic Dis. 2016, 16, 423–427. [Google Scholar] [CrossRef]
  68. Conrardy, C.; Tao, Y.; Kuzmin, I.V.; Niezgoda, M.; Agwanda, B.; Breiman, R.F.; Anderson, L.J.; Rupprecht, C.E.; Tong, S. Short report: Molecular detection of adenoviruses, rhabdoviruses, and paramyxoviruses in bats from Kenya. Am. J. Trop. Med. Hyg. 2014, 91, 258–266. [Google Scholar] [CrossRef] [Green Version]
  69. Lumsden, L.F.; Bennett, A.F.; Silins, J.E. Selection of roost sites by the lesser long-eared bat (Nyctophilus geoffroyi) and Gould’s wattled bat (Chalinolobus gouldii) in south-eastern Australia. J. Zool. 2002, 257, S095283690200081X. [Google Scholar] [CrossRef]
  70. Rueegger, N.; Law, B.; Goldingay, R. Interspecific differences and commonalities in maternity roosting by tree cavity-roosting bats over a maternity season in a timber production landscape. PLoS ONE 2018, 13, e0194429. [Google Scholar] [CrossRef] [Green Version]
  71. Burgar, J.M.; Craig, M.D.; Stokes, V.L. The importance of mature forest as bat roosting habitat within a production landscape. For. Ecol. Manage. 2015, 356, 112–123. [Google Scholar] [CrossRef] [Green Version]
  72. Bender, R.; The Melbourne bat box monitoring program, Melbourne, Australia. Personal communication, 2018.
  73. Zhou, P.; Fan, H.; Lan, T.; Yang, X.-L.; Shi, W.-F.; Zhang, W.; Zhu, Y.; Zhang, Y.-W.; Xie, Q.-M.; Mani, S.; et al. Fatal swine acute diarrhoea syndrome caused by an HKU2-related coronavirus of bat origin. Nature 2018, 556, 255–258. [Google Scholar] [CrossRef]
  74. Benkö, M.; Harrach, B. Molecular evolution of adenoviruses. In Adenoviruses: Model and Vectors in Virus-Host Interactions; Doerfler, W., Böhm, P., Eds.; Springer: Berlin/Heidelberg, Germany, 2003; pp. 3–35. [Google Scholar]
  75. Hackenbrack, N.; Rogers, M.B.; Ashley, R.E.; Keel, M.K.; Kubiski, S.V.; Bryan, J.A.; Ghedin, E.; Holmes, E.C.; Hafenstein, S.L.; Allison, A.B. Evolution and Cryo-electron Microscopy Capsid Structure of a North American Bat Adenovirus and Its Relationship to Other Mastadenoviruses. J. Virol. 2017, 91, e01504-16. [Google Scholar] [CrossRef] [Green Version]
  76. Sonntag, M.; Mühldorfer, K.; Speck, S.; Wibbelt, G.; Kurth, A. New adenovirus in bats, Germany. Emerg. Infect. Dis. 2009, 15, 2052–2055. [Google Scholar] [CrossRef]
  77. Jánoska, M.; Vidovszky, M.; Molnár, V.; Liptovszky, M.; Harrach, B.; Benkő, M. Novel adenoviruses and herpesviruses detected in bats. Vet. J. 2011, 189, 118–121. [Google Scholar] [CrossRef] [PubMed]
  78. Wu, Z.; Ren, X.; Yang, L.; Hu, Y.; Yang, J.; He, G.; Zhang, J.; Dong, J.; Sun, L.; Du, J.; et al. Virome analysis for identification of novel mammalian viruses in bat species from Chinese provinces. J. Virol. 2012, 86, 10999–11012. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  79. Bensch, S.; Hellgren, O.; Pérez-Tris, J. MalAvi: a public database of malaria parasites and related haemosporidians in avian hosts based on mitochondrial cytochrome b lineages. Mol. Ecol. Resour. 2009, 9, 1353–1358. [Google Scholar] [CrossRef] [PubMed]
  80. Drexler, J.F.; Corman, V.M.; Müller, M.A.; Maganga, G.D.; Vallo, P.; Binger, T.; Gloza-Rausch, F.; Cottontail, V.M.; Rasche, A.; Yordanov, S.; et al. Bats host major mammalian paramyxoviruses. Nat. Commun. 2012, 3, 796. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  81. Fischer, K.; Zeus, V.; Kwasnitschka, L.; Kerth, G.; Haase, M.; Groschup, M.H.; Balkema-Buschmann, A. Insectivorous bats carry host specific astroviruses and coronaviruses across different regions in Germany. Infect. Genet. Evol. 2016, 37, 108–116. [Google Scholar] [CrossRef]
  82. Peel, A.J.; McKinley, T.J.; Baker, K.S.; Barr, J.A.; Crameri, G.; Hayman, D.T.S.; Feng, Y.-R.; Broder, C.C.; Wang, L.-F.; Cunningham, A.A.; et al. Use of cross-reactive serological assays for detecting novel pathogens in wildlife: Assessing an appropriate cutoff for henipavirus assays in African bats. J. Virol. Methods 2013, 193, 295–303. [Google Scholar] [CrossRef] [Green Version]
  83. Burroughs, A.L.; Durr, P.A.; Boyd, V.; Graham, K.; White, J.R.; Todd, S.; Barr, J.; Smith, I.; Baverstock, G.; Meers, J.; et al. Hendra virus infection dynamics in the grey-headed flying fox (Pteropus poliocephalus) at the southern-most extent of its range: Further evidence this species does not readily transmit the virus to horses. PLoS ONE 2016, 11, e0155252. [Google Scholar] [CrossRef]
  84. Wong, S.; Lau, S.; Woo, P.; Yuen, K.-Y. Bats as a continuing source of emerging infections in humans. Rev. Med. Virol. 2007, 17, 67–91. [Google Scholar] [CrossRef]
  85. Luis, A.D.; Hayman, D.T.S.; O’Shea, T.J.; Cryan, P.M.; Gilbert, A.T.; Pulliam, J.R.C.; Mills, J.N.; Timonin, M.E.; Willis, C.K.R.; Cunningham, A.A.; et al. A comparison of bats and rodents as reservoirs of zoonotic viruses: are bats special? Proc. Biol. Sci. 2013, 280, 1–9. [Google Scholar] [CrossRef] [Green Version]
  86. 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. [Google Scholar] [CrossRef]
  87. Turmelle, A.S.; Olival, K.J. Correlates of Viral Richness in Bats (Order Chiroptera). Ecohealth 2009, 6, 522–539. [Google Scholar] [CrossRef] [PubMed]
  88. Maganga, G.D.; Bourgarel, M.; Vallo, P.; Dallo, T.D.; Ngoagouni, C.; Drexler, J.F.; Drosten, C.; Nakouné, E.R.; Leroy, E.M.; Morand, S. Bat distribution size or shape as determinant of viral richness in African bats. PLoS ONE 2014, 9, e100172. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  89. Gay, N.; Olival, K.J.; Bumrungsri, S.; Siriaroonrat, B.; Bourgarel, M.; Morand, S. Parasite and viral species richness of Southeast Asian bats: Fragmentation of area distribution matters. Int. J. Parasitol. Parasites Wildl. 2014, 3, 161–170. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  90. Bordes, F.; Morand, S. The impact of multiple infections on wild animal hosts: a review. Infect. Ecol. Epidemiol. 2011, 1. [Google Scholar] [CrossRef] [Green Version]
  91. Woolhouse, M.E.J.; Gowtage-Sequeria, S. Host range and emerging and reemerging pathogens. Emerg. Infect. Dis. 2005, 11, 1842–1847. [Google Scholar] [CrossRef]
  92. Streicker, D.G.; Turmelle, A.S.; Vonhof, M.J.; Kuzmin, I.V.; McCracken, G.F.; Rupprecht, C.E. Host phylogeny constrains cross-species emergence and establishment of rabies virus in bats. Science (80-. ) 2010, 329, 676–679. [Google Scholar] [CrossRef] [Green Version]
  93. Basham, R.; Law, B.; Banks, P. Microbats in a “leafy” urban landscape: Are they persisting, and what factors influence their presence? Austral Ecol. 2011, 36, 663–678. [Google Scholar] [CrossRef]
  94. Straka, T.M.; Lentini, P.E.; Lumsden, L.F.; Wintle, B.A.; van der Ree, R. Urban bat communities are affected by wetland size, quality, and pollution levels. Ecol. Evol. 2016, 6, 4761–4774. [Google Scholar] [CrossRef] [Green Version]
  95. Bullen, R.D.; McKenzie, N.L. Seasonal range variation of Tadarida australis (Chiroptera : Molossidae) in Western Australia: The impact of enthalpy. Aust. J. Zool. 2005, 53, 145–156. [Google Scholar] [CrossRef]
  96. Hammarin, A.L.; Berndtsson, L.T.; Falk, K.; Nedinge, M.; Olsson, G.; Lundkvist, Å. Lyssavirus-reactive antibodies in Swedish bats. Infect. Ecol. Epidemiol. 2016, 6, 31262. [Google Scholar] [CrossRef]
  97. Amengual, B.; Bourhy, H.; López-Roig, M.; Serra-Cobo, J. Temporal dynamics of European bat lyssavirus type 1 and survival of Myotis myotis bats in natural colonies. PLoS ONE 2007, 2, e566. [Google Scholar] [CrossRef] [PubMed]
  98. Peel, A.J.; Wells, K.; Giles, J.; Boyd, V.; Burroughs, A.; Edson, D.; Crameri, G.; Baker, M.L.; Field, H.; Wang, L.-F.; et al. Synchronous shedding of multiple bat paramyxoviruses coincides with peak periods of Hendra virus spillover. Emerg. Microbes Infect. 2019, 8, 1314–1323. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. The South West Botanical province (SWBP), highlighted in yellow, sampling sites and extent of the Australian wheatbelt are displayed. The SWBP encompasses nine bioregions, Avon Wheatbelt (AVW), Coolgardie (COO), Esperance Plains (ESP), Geraldton Sandplains (GES), Hampton (HAM), Jarrah Forest (JAF), Mallee (MAL), Swan Coastal Plain (SWA), and Warren (WAR).
Figure 1. The South West Botanical province (SWBP), highlighted in yellow, sampling sites and extent of the Australian wheatbelt are displayed. The SWBP encompasses nine bioregions, Avon Wheatbelt (AVW), Coolgardie (COO), Esperance Plains (ESP), Geraldton Sandplains (GES), Hampton (HAM), Jarrah Forest (JAF), Mallee (MAL), Swan Coastal Plain (SWA), and Warren (WAR).
Viruses 11 01157 g001
Figure 2. Bayesian phylogenetic analysis of adenovirus sequences based on a fragment of the DNApol gene. Tree was constructed with Mr Bayes, posterior probabilities of mayor partitions are reported below branches. Clades identified in this study are coloured, naming convention represents the unique ID for each individual followed by the first letter of the genus, the first letter of the species and the geographical trapping site. A full adenovirus genome was amplified from the sample marked by an orange circle and samples that did not cluster within any of the identified WA clades are underlined. GeneBank accession numbers for representative sequences produced by this study are reported in Supplementary Table S1. Within Mastoadenovirus, clades of reference samples were collapsed for brevity, black circles, and accession numbers are shown following strain names. Genbank accession numbers for representatives of Siadavirus, Aviadvirus and Atadenovirus are reported in Supplementary Table S2.
Figure 2. Bayesian phylogenetic analysis of adenovirus sequences based on a fragment of the DNApol gene. Tree was constructed with Mr Bayes, posterior probabilities of mayor partitions are reported below branches. Clades identified in this study are coloured, naming convention represents the unique ID for each individual followed by the first letter of the genus, the first letter of the species and the geographical trapping site. A full adenovirus genome was amplified from the sample marked by an orange circle and samples that did not cluster within any of the identified WA clades are underlined. GeneBank accession numbers for representative sequences produced by this study are reported in Supplementary Table S1. Within Mastoadenovirus, clades of reference samples were collapsed for brevity, black circles, and accession numbers are shown following strain names. Genbank accession numbers for representatives of Siadavirus, Aviadvirus and Atadenovirus are reported in Supplementary Table S2.
Viruses 11 01157 g002
Figure 3. Bayesian phylogenetic analysis of coronavirus sequences based on a fragment of the RdRp gene. Tree was constructed with Mr Bayes, posterior probabilities of major partitions are reported below branches. Clades identified in this study are coloured, naming convention is based on the unique ID for each individual followed by the first letter of the genus, the first letter of the species and the geographical trapping site. Full coronavirus genomes were amplified from the samples marked by an orange circle. GenBank accession numbers for representative sequences produced by this study are reported in Supplementary Table S1. Within Alphacoronavirus, clades of reference samples were collapsed for brevity, black circles, and accession numbers are shown following strain names. GenBank accession numbers for representatives of Betacoronavirus, Deltacoronavirus and Gammacoronavirus are reported in Supplementary Table S2.
Figure 3. Bayesian phylogenetic analysis of coronavirus sequences based on a fragment of the RdRp gene. Tree was constructed with Mr Bayes, posterior probabilities of major partitions are reported below branches. Clades identified in this study are coloured, naming convention is based on the unique ID for each individual followed by the first letter of the genus, the first letter of the species and the geographical trapping site. Full coronavirus genomes were amplified from the samples marked by an orange circle. GenBank accession numbers for representative sequences produced by this study are reported in Supplementary Table S1. Within Alphacoronavirus, clades of reference samples were collapsed for brevity, black circles, and accession numbers are shown following strain names. GenBank accession numbers for representatives of Betacoronavirus, Deltacoronavirus and Gammacoronavirus are reported in Supplementary Table S2.
Viruses 11 01157 g003
Figure 4. Maximum likelihood phylogenetic analysis of the pp1ab amino acid sequences derived from five coronavirus genomes. The tree was constructed in RAxML using the PROTGAMMAWAG model with 1000 bootstraps; supports above 80% are shown below each branch. Naming convention for the sequences generated in this study is based on the unique ID for each individual followed by the first letter of the genus, the first letter of the species and the geographical trapping site.
Figure 4. Maximum likelihood phylogenetic analysis of the pp1ab amino acid sequences derived from five coronavirus genomes. The tree was constructed in RAxML using the PROTGAMMAWAG model with 1000 bootstraps; supports above 80% are shown below each branch. Naming convention for the sequences generated in this study is based on the unique ID for each individual followed by the first letter of the genus, the first letter of the species and the geographical trapping site.
Viruses 11 01157 g004
Figure 5. Bayesian phylogenetic analysis of paramyxovirus sequences based on a fragment of the RdRp gene. Tree was constructed with Mr Bayes, posterior probabilities of major partitions are reported below branches. Clades identified in this study are coloured, naming convention for sequences produced by this study is based on the unique ID for each individual, followed by the first letters of the genus, the species and the geographical trapping site. Samples that did not cluster within any the identified WA clades are underlined. GeneBank accession numbers for representative sequences produced by this study are reported in Supplementary Table S1. Clades of reference samples were collapsed for brevity, black circles, GenBank accession numbers for representatives of Rubulavirus, Avulavirus, Respirovirus, Henipavirus, Metapneuavirus and Pneumovirus are reported in Supplementary Table S2.
Figure 5. Bayesian phylogenetic analysis of paramyxovirus sequences based on a fragment of the RdRp gene. Tree was constructed with Mr Bayes, posterior probabilities of major partitions are reported below branches. Clades identified in this study are coloured, naming convention for sequences produced by this study is based on the unique ID for each individual, followed by the first letters of the genus, the species and the geographical trapping site. Samples that did not cluster within any the identified WA clades are underlined. GeneBank accession numbers for representative sequences produced by this study are reported in Supplementary Table S1. Clades of reference samples were collapsed for brevity, black circles, GenBank accession numbers for representatives of Rubulavirus, Avulavirus, Respirovirus, Henipavirus, Metapneuavirus and Pneumovirus are reported in Supplementary Table S2.
Viruses 11 01157 g005
Table 1. Primer sets used in this study for the screening of Adenoviruses, Coronaviruses, and Paramyxoviruses.
Table 1. Primer sets used in this study for the screening of Adenoviruses, Coronaviruses, and Paramyxoviruses.
Target GroupTarget RegionPrimer NameT °CAmplicon SizeReference
AdenovirusDNA polymerase gene (DNApol)F155205[39]
R1
F255
R2
CoronavirusRNA-dependent RNA polymerase gene (RdRp)CoV-Fwd153400[40]
CoV-Rvs2
CoV-Fwd250
PLQ-F150300[41]
PLQ-R1
PLQ-F250
PLQ-R2
ParamyxovirusRNA-dependent RNA polymerase gene (RdRp)PAR-F153500[35]
PAR-R
PAR-F250
RES-MOR-HEN-F153500
RES-MOR-HEN-R
RES-MOR-HEN-F250
Table 2. Details of species sampled, number of faecal and serum samples collected per species and collection sites.
Table 2. Details of species sampled, number of faecal and serum samples collected per species and collection sites.
FamilySpeciesN FaecalN SerumTrapping Locations
VespertilionidaeChalinolobus gouldii232264MtG, CDR, Coc, Dry, Dwe, Koj, Per
Chalinolobus morio4566MtG, Bla, CDR, Coc, Don, Dry, Dwe, Frk, Met, Twi, Two
Falsistrellus mackenziei117Bla, Dwe
Nyctophilus geoffroyi5150MtG, Bla, CDR, Coc, Don, Dry, Dwe, Ger, Per, Two
Nyctophilus gouldi5672Bla, Don, Dry, Dwe, Per, Two
Nyctophilus major106Bla, CDR, Dry, Dwe, Frk
Scotorepens balstoni98MtG, CDR
Vespadelus baverstocki45MtG, CDR
Vespadelus regulus141170Bla, Coc, Don, Dry, Dwe, Frk Met, Per, Twi, Two
MolossidaeAustronomus australis911CDR
Ozimops sp33CDR
Table 3. Overall seroprevalence to six antigens. Seroprevalence and confidence intervals (95% CI) presented as percentages.
Table 3. Overall seroprevalence to six antigens. Seroprevalence and confidence intervals (95% CI) presented as percentages.
FamilyGenusAntigenN 1Pos.2Seroprevalence %
ParamyxoviridaeHenipavirusHeV645335.1 (3.6–7.1)
CedV637345.3 (3.8–7.3)
NiV64891.3 (0.7–2.6)
RubulavirusTioV653172.6 (1.6–4.1)
CoronaviridaeBetacoronavirusSARS-CoV645385.8 (4.3–7.9)
MERS-CoV659101.5 (0.8–2.7)
1 Total number of samples that were successfully tested. 2 Positives: number of samples above the calculated MFI threshold (1000 MFI).
Table 4. Seroprevalence per antigen per species. Tested samples (N), followed by the number of positives () are shown. Seroprevalence (SeP) and 95% confident intervals (CI) are presented as percentages. Last two columns present the seroprevalence range for each species. Estimates were not calculated for species with positive results but total sample sizes under 20 individuals (NC). Antigen abbreviations are as follows: Hendra virus (HeV), Nipah virus (NiV), Cedar virus (CedV), Tioman virus (TioV), SARS coronavirus (SARS-CoV), and MERS coronavirus (MERS-CoV).
Table 4. Seroprevalence per antigen per species. Tested samples (N), followed by the number of positives () are shown. Seroprevalence (SeP) and 95% confident intervals (CI) are presented as percentages. Last two columns present the seroprevalence range for each species. Estimates were not calculated for species with positive results but total sample sizes under 20 individuals (NC). Antigen abbreviations are as follows: Hendra virus (HeV), Nipah virus (NiV), Cedar virus (CedV), Tioman virus (TioV), SARS coronavirus (SARS-CoV), and MERS coronavirus (MERS-CoV).
Paramyxovirus AntigensBetacoronavirus AntigensOverall Summary
HeVCedVNiVTioVSARS-CoVMERS-CoV
SpeciesNSePN SePNSePNSePNSePNSePPaV
Range
CoV
Range
C. gouldii66(2)3
(0.8–10.4)
66(4)6
(2.4–14.6)
68(2)3
(0.8–10)
68(4)6
(2.3–14.2)
66(2)3
(0.8–10.4)
68(2)3
(0.8–10.1)
3–63
C. morio60(5)8.3
(3.6–18.0)
58(6)10.3
(4.8–20.7)
59(0) 58(2)3.3
(0.9–11.5)
59(4)6.7
(2.6–16)
60(1)1.6
(0.2–8.7)
3.3–10.31.6–6.7
F. mackenziei7(7)NC7(7)NC7(3)NC7(0) 7(0) 7(0) NC
N. geoffroyi27(0) 26(0) 26(1)3.8
(0.6–1.8)
26(0) 26(0) 27(0) NC
N. gouldi69(2)3
(0.8–10)
63(2)1.6
(0.3–8.4)
67(0) 69(0) 67(7)10
(5.1–20)
70(0) 1.6–310
N. major5(0) 5(0) 5(0) 5(1)NC5(2)NC4(0) NCNC
S. balstoni2(0) 2(0) 2(0) 2(0) 2(0) 2(0)
V. baverstocki1(0) 1(0) 1(0) 1(0) 1(0) 1(0)
V. regulus147(7)4.8
(2.3–9.5)
146(9)6.1
(3.2–11.3)
150(1)0.6
(0.1–3.7)
154(2)1.3
(0.3–4.6)
150(9)6
(3–11)
155(4)2.5
(1.0–6.4)
0.6–6.12.5–6
A. australis11(0) 11(0) 11(0) 11(1)NC11(1)NC11(0) NCNC
O. sp1(0) 1(0) 1(0) 1(0) 1(0) 1(0)
Table 5. Viral shedding prevalence of three viral families in 11 species of microbats of the South West Botanical Province of Western Australia. Total number of samples tested (N), prevalence values (%), confidence interval (95% CI) and number of positive samples () are shown.
Table 5. Viral shedding prevalence of three viral families in 11 species of microbats of the South West Botanical Province of Western Australia. Total number of samples tested (N), prevalence values (%), confidence interval (95% CI) and number of positive samples () are shown.
FamilySpeciesNAdenoviridaeCoronaviridaeParamyxoviridae
VespertilionidaeC. gouldii23217 (13–23) (41)25 (20–31) (59)4 (2–7) (10)
C. morio4516 (8–29) (7)36 (24–51) (17)9 (3–21) (4)
F. mackenziei11 NC1 (2)
N. geoffroyi514 (1–13) (2)4 (1–13) (2)2 (0.3–10) (1)
N. gouldi561.8 (0.3–9) (1)3.5 (1–12) (2)
N. major10
S. balstoni9 NC1 (1)
V. baverstocki4 NC1 (2)
V. regulus1418 (4–13) (11)12 (8–18) (17)1.4 (0.4–5) (2)
MolossidaeA. australis9
O. sp3 NC1 (1)
1 Not calculated—calculations were not made for species with total sample sizes below 20 individuals.

Share and Cite

MDPI and ACS Style

Prada, D.; Boyd, V.; Baker, M.L.; O’Dea, M.; Jackson, B. Viral Diversity of Microbats within the South West Botanical Province of Western Australia. Viruses 2019, 11, 1157. https://doi.org/10.3390/v11121157

AMA Style

Prada D, Boyd V, Baker ML, O’Dea M, Jackson B. Viral Diversity of Microbats within the South West Botanical Province of Western Australia. Viruses. 2019; 11(12):1157. https://doi.org/10.3390/v11121157

Chicago/Turabian Style

Prada, Diana, Victoria Boyd, Michelle L. Baker, Mark O’Dea, and Bethany Jackson. 2019. "Viral Diversity of Microbats within the South West Botanical Province of Western Australia" Viruses 11, no. 12: 1157. https://doi.org/10.3390/v11121157

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop