Coronavirus Infection and Diversity in Bats in the Australasian Region.

Following the SARS outbreak, extensive surveillance was undertaken globally to detect and identify coronavirus diversity in bats. This study sought to identify the diversity and prevalence of coronaviruses in bats in the Australasian region. We identified four different genotypes of coronavirus, three of which (an alphacoronavirus and two betacoronaviruses) are potentially new species, having less than 90% nucleotide sequence identity with the most closely related described viruses. We did not detect any SARS-like betacoronaviruses, despite targeting rhinolophid bats, the putative natural host taxa. Our findings support the virus-host co-evolution hypothesis, with the detection of Miniopterus bat coronavirus HKU8 (previously reported in Miniopterus species in China, Hong Kong and Bulgaria) in Australian Miniopterus species. Similarly, we detected a novel betacoronavirus genotype from Pteropus alecto which is most closely related to Bat coronavirus HKU9 identified in other pteropodid bats in China, Kenya and the Philippines. We also detected possible cross-species transmission of bat coronaviruses, and the apparent enteric tropism of these viruses. Thus, our findings are consistent with a scenario wherein the current diversity and host specificity of coronaviruses reflects co-evolution with the occasional host shift.

Introduction

Novel coronaviruses were responsible for the global outbreak of severe acute respiratory syndrome (SARS) in 2003–2004, and for the ongoing sporadic outbreaks of Middle East respiratory syndrome (MERS), first described in 2012 (Drosten et al. 2003; Zaki et al. 2012). Bats have been identified as the natural hosts of the progenitor viruses of both SARS and MERS coronaviruses (Lau et al. 2005b; Li et al. 2005; Ge et al. 2013; Memish et al. 2013). Following the SARS outbreak, extensive surveillance was undertaken globally to detect and identify coronavirus diversity in bats, given their evident potential to host novel zoonotic coronaviruses (Woo et al. 2006). The majority of these studies were based on reverse-transcriptase PCR detection of coronavirus RNA in bat faeces and tissues, or serological detection of coronavirus antibodies in blood samples. Coronaviruses have been detected in 11 of the 18 extant families of bats in more than 25 countries in Asia, Africa, Europe and the Americas (Drexler et al. 2014).

Coronaviruses are enveloped positive-sense RNA viruses. Their taxonomy has undergone major revision over the past decade, primarily because of the identification of a suite of novel coronaviruses in bats. The family Coronaviridae is now divided into the two subfamilies of Coronavirinae and Torovirinae, with the former containing the genera Alpha-, Beta-, Gamma- and Deltacoronavirus (International Committee on Taxonomy of Viruses 2011). SARS and MERS coronaviruses belong to the genus Betacoronavirus, with all coronaviruses detected in bats being either alpha- or betacoronaviruses (Drexler et al. 2014).

Bat coronaviruses have a narrow host range and are typically species, genus or family-specific, independent of location (Poon et al. 2004; Chu et al. 2006; Tang et al. 2006; Woo et al. 2006; Gloza-Rausch et al. 2008; Pfefferle et al. 2009; Drexler et al. 2014). Drexler et al. (2010) hypothesised that this virus-host specificity could be used in a prospective manner to predict the geographic distribution of various bat coronavirus species, and showed that the same alphacoronavirus (Miniopterus bat coronavirus HKU8), which had previously been reported in Miniopterus species in China and Hong Kong, occurred in M. schreibersii in Bulgaria, over 8000 km distant. Further, related coronaviruses hosted by bats of the same family have been detected in Vespertilionidae (Cui et al. 2007) and Rhinolophidae bats (Lau et al. 2005b; Li et al. 2005).

The diversity and prevalence of coronaviruses in bats in the Australasian region has not yet been reported. Given the zoonotic potential of bat coronaviruses, a better understanding of the distribution and ecology of these viruses is essential to identify any potential threat to human health and to inform biosecurity preparedness. With this aim, we undertook a virological and serological survey of various bat populations in Australia and neighbouring countries for evidence of coronavirus infection. We focused primarily on insectivorous bats, which reportedly host the largest diversity of bat coronaviruses, and particularly targeted bats of the genus Rhinolophus, recognised as hosting SARS-like coronaviruses. We also compared relative detection frequency in various potential routes of excretion.

Materials and Methods

Study Samples

Faecal samples, anal swabs or intestinal content samples were obtained by either primary fieldwork or sourced from archived samples.

Fieldwork: Bat Capture and Sample Collection

Fieldwork was undertaken at multiple locations in Australia and Taiwan between 2006 and 2009, using harp traps to capture insectivorous bats. Particular effort was made to capture bats at locations where rhinolophid bats were previously known to be present. Captured bats were individually housed in clean cloth bags and a polystyrene cooler prior to sample collection. Bat identification was based on morphological characteristics (Churchill 2008) and was undertaken by experienced biologists. Occasionally, bats were only identified to genus level. All bats were released at their point of capture within 6 h. A single faecal pellet (collected directly from a defecating bat or from its clean cloth bag) was placed into 1 ml of sucrose potassium glutamate albumin (SPGA) with added penicillin, streptomycin and fungizone. When no faecal pellet was obtained, the anus was swabbed and the swab placed in 1 ml SPGA, and stored at 4°C prior to RNA extraction. Faecal samples collected from bats in Taiwan were placed into 1 ml of AVL buffer from the QIAamp® Viral RNA Mini Kit (QIAGEN) and stored at room temperature for 1 week prior to extraction. Faecal swabs from 33 Pteropus alecto and 27 P. poliocephalus opportunistically obtained from colleagues conducting concurrent henipavirus research in south-east Queensland in 2009 were placed in 1 ml SPGA and stored at 4°C prior to extraction.

Bats were bled as described by Smith et al. (2010), and blood was diluted 1:10 in phosphate-buffered saline. In addition to faecal and blood samples, oropharyngeal swabs were collected from 30 bats (14 Miniopterus australis, 16 M. schreibersii) from south-east Queensland.

Bat capture and sample collection was undertaken in direct accordance with animal ethics committee permits SA 2006/06/117 and SA 2007/005/194 issued by the (then) Queensland Department of Primary Industries and Fisheries, and scientific purposes permits WISP03887606 and WISP04906107 issued by the Queensland Parks and Wildlife Service of the (then) Queensland Environmental Protection Agency.

Archived Samples

Twenty archival samples of intestines (stored at −70°C) from bats caught in central Queensland in 1996 were homogenised (2 mm2 section) in 1 ml of SPGA. In addition, archived (−20°C) serum samples collected from various bat species between 1999 and 2006 by colleagues in the course of henipavirus and other surveillance were serologically screened (Western Australia, n = 119; Northern Territory, n = 40; East Timor, n = 36; Indonesia, n = 67; Malaysia, n = 101 and Papua New Guinea, n = 65). Archived samples were sourced from the Queensland Department of Agriculture and Fisheries Biological Sciences Laboratory (BSL) in Brisbane, Australia.

Coronavirus RNA Detection and Sequencing

Faecal pellets, intestinal samples, anal and oropharyngeal swabs in SPGA were centrifuged as required, and RNA was extracted from 560 μl of the supernatant using the QIAamp® Viral RNA Mini Kit (QIAGEN) following the manufacturer’s instructions. RNA was similarly extracted from 140 μl of diluted blood. Reverse transcription PCR (RT-PCR) targeting a conserved region of the coronavirus RNA-dependent RNA polymerase (RdRp) gene, as described by Poon et al. (2005), was performed using the Superscript III One-Step RT-PCR System with Platinum® Taq DNA Polymerase (Invitrogen). Amplicons consistent with the expected length of 440 bp were purified using the QIAquick® PCR Purification Kit (QIAGEN) as per the manufacturer’s instructions. Purified amplicons were directly sequenced in both directions using BigDye® Terminator v3.1 Cycle Sequencing Kit (Applied Biosystems) as per the manufacturer’s instructions. The extension products were purified using the ethanol/EDTA precipitation method (Applied Biosystems 2002) and analysed at the Griffith University DNA Sequencing Facility (Brisbane, Australia). Nucleotide sequence traces were edited using Sequence Scanner v1.0 (Applied Biosystems).

Molecular Phylogenetic Analysis

Phylogenetic analysis of a 440 bp sequence from the RNA-dependent RNA polymerase (RdRp) gene from representative sequences from this study, and from 33 other sequences including the ICTV reference coronaviruses, was conducted using the Maximum Likelihood method based on the General Time Reversible model (Nei and Kumar 2000). Initial tree(s) for the heuristic search were obtained by applying the Neighbour-Joining method to a matrix of pairwise distances estimated using the Maximum Composite Likelihood (MCL) approach. A discrete Gamma distribution was used to model evolutionary rate differences amongst sites [5 categories (+G, parameter = 0.6965)]. The rate variation model allowed for some sites to be evolutionarily invariable ([+I], 0.0000% sites). Evolutionary analyses were conducted in MEGA6 (Tamura et al. 2013).

Anti-coronavirus Antibody Detection

The archived serum samples were screened for cross-reactive antibodies to SARS or SARS-like viruses using a SARS coronavirus crude antigen ELISA modified from the original protocol developed by Yu et al. (2008). Briefly, total viral antigen was used to coat a 96-well Nunc Max ELISA plate in 0.1 M carbonate coating buffer, followed by incubation for 30 min (all incubation was done at 37°C). After washing with PBST buffer, the plate was incubated for 30 min with blocking buffer (2% skimmed milk in PBST). Testing sera were diluted 1:50 in blocking buffer and incubated for 30 min, followed by washing with PBST. HRP-conjugated Protein-A/G (1:1 mix, ZYMED Laboratories) was used to detect the bound bat antibodies by incubation for 30 min with a 1:2000 diluted solution. After washing, TMB substrate was added for colour development and the reading was taken after 10 min at 450 nm using an ELISA plate reader (Multiskan Ascent, Thermo Labsystems).

Prevalence Estimation

Prevalence of coronavirus RNA and coronavirus antibodies was estimated for species and country of sample origin (in Australia, for species and state/territory of sample origin), with 95% confidence intervals calculated using the Jeffreys method (Brown et al. 2001).

Results

Field and archived samples were derived from 2195 bats from Australia, Papua New Guinea, East Timor, Indonesia, Malaysia and Taiwan (Fig. 1). The sampled bats belonged to 22 different genera within eight families (Table 1). In total, 1484 faecal/anal swab samples, 30 blood samples and 30 oropharyngeal swab samples from captured bats, and 20 archived bat intestinal samples were subjected to RT-PCR. Serology was performed on a total 1808 serum samples (Table 1).

Figure 1

Sampling locations. Locations of 2195 bats sampled between 1996 and 2009 and screened for evidence of coronavirus infection, comprising bats from south-east Queensland (SEQ, n = 1162), central Queensland (CQ, n = 42), far-north Queensland (FNQ, n = 222), the Northern Territory (NT, n = 333), Western Australia (WA, n = 119), Taiwan (n = 43), East Timor (n = 36), Indonesia (n = 67), Malaysia (n = 101) and Papua New Guinea (n = 65).

Table 1

Coronavirus RNA and Antibodies in Bats in Australasia.

Suborder	Family	Genus	Species	Location1	Year	Coronavirus RNA2	Coronavirus antibodies3
Pteropodiformes	Hipposideridae	Hipposideros	ater	Australia (FNQ)	2006–2009	0, 0–8.2 (29)	0, 0–8.2 (29)
				Australia (NT)	2006–2009	0, 0–8.8 (27)	25, 2.8–71.6 (4)
				Australia (WA)	2000–2002		0, 0–7.7 (31)
			terasensis	Taiwan	2008	0, 0–66.7 (2)
		Rhinonicteris	aurantia	Australia (NT)4	2006–2009	1, 0.1–3.6 (126)	0, 0–2.4 (105)
	Megadermatidae	Macroderma	gigas	Australia (NT)	2006–2009	0, 0–4.3 (57)	1.6, 0.2–7.2 (63)
				Australia (WA)	2000–2002		81, 60.8–93.2 (21)
	Pteropodidae	Acerodon	celebensis	Indonesia	2000–2006		0, 0–15.2 (15)
		Cynopterus	spp.	Malaysia	1999		73.3, 48.3–90.3 (15)
		Dobsonia	anderseni	Papua New Guinea	1996–1999		5.6, 0.6–23.2 (18)
			peronii	Indonesia	1996–1999		0, 0–85.3 (1)
			praedatrix	Papua New Guinea	1996–1999		0, 0–11.1 (21)
		Eonycteris	spp.	Malaysia	1999		91.7, 67.2–91.1 (12)
		Macroglossus	minimus	Papua New Guinea	2000–2006		0, 0–66.7 (2)
			spp.	Indonesia	2000–2006		0, 0–53.6 (3)
		Pteropus	alecto	Australia (SEQ)5	2009	12, 4.2–26.3 (33)	26.5, 14–42.8 (34)
				Indonesia	2000–2006		0, 0–6.7 (36)
				Papua New Guinea	2000–2006		90.9, 64.7–99 (11)
			capistratus	Papua New Guinea	1996–1999		0, 0–29.2 (7)
			conspicillatus	Australia (FNQ)	2000–2006		15, 6.5–28.3 (40)
			griseus	East Timor	2000–2006		0, 0–85.3 (1)
			hypomelanus	Malaysia	1999		0, 0–7.1 (34)
			neohibernicus	Papua New Guinea	2000–2006		66.7, 28.6–92.3 (6)
			poliocephalus	Australia (SEQ)	2009	0, 0–8.8 (27)	16.4, 9.3–26.2 (73)
			scapulatus	Australia (NT)	2000–2006		7.5, 2.2–18.7 (40)
			vampyrus	East Timor	2000–2006		11.4, 4–24.9 (35)
				Malaysia	1999		37.5, 22.4–54.8 (32)
		Rousettus	amplexicaudatus	Indonesia	2000–2006	0, 0–33 (6)
			spp.	Indonesia	2000–2006		16.7, 1.9–55.8 (6)
	Rhinolophidae	Rhinolophus	megaphyllus	Australia (FNQ)6	2006–2009	1.7, 0.2–7.8 (58)	8.2, 3.2–17 (61)
				Australia (SEQ)	2006–2009	0, 0–0.6 (448)	3.3, 1.8–5.4 (399)
			monoceros	Taiwan	2008	0, 0–5.9 (41)
Vespertilioniformes	Emballonuridae	Saccolaimus	flaviventris	Australia (WA)	2000–2002		0, 0–12.9 (18)
		Taphozous	spp.	Australia (WA)	2000–2002		21.1, 10.5–35.8 (38)
				Malaysia	1999		25, 2.8–71.6 (4)
	Miniopteridae	Miniopterus	australis	Australia (CQ)6	1996	5, 0.5–21.1 (20)	50, 32.8–67.2 (30)
				Australia (FNQ)6	2006–2009	46.7, 29.8–64.1 (30)	53.3, 35.9–70.2 (30)
				Australia (SEQ)6	2006–2009	24.7, 18.4–31.9 (154)	64.5, 55.8–72.5 (124)
				Australia (WA)	2000–2002		100, 14.7–100 (1)
			schreibersii	Australia (NT)6	2006–2009	10.2, 4.4–19.8 (59)	46.4, 33.8–59.4 (56)
				Australia (SEQ)6	2006–2009	26.5, 21.2–32.3 (238)	68.7, 62.2–74.7 (211)
	Molossidae	Chaerephon	jobensis	Australia (WA)	2000–2002		50, 12.3–87.7 (4)
		Mormopterus	beccarii	Australia (SEQ)	2006–2009	0, 0–53.6 (3)	97.6, 89.2–99.7 (41)
			norfolkensis	Australia (SEQ)	2006–2009	0, 0–85.3 (1)
	Vespertilionidae	Chalinolobus	spp.	Australia (WA)	2000–2002		50, 12.3–87.7 (4)
		Myotis	macropus	Australia (FNQ)	2006–2009	0, 0–7.7 (31)	58.1, 40.6–74.1 (31)
				Australia (SEQ)	2006–2009	20.3, 11.9–31.3 (64)
		Nyctophilus	bifax	Australia (SEQ)	2006–2009	0, 0–33 (6)
			gouldi	Australia (SEQ)	2006–2009	0, 0–29.2 (7)
		Scotophilus	spp.	Malaysia	1999		100, 55.5–100 (4)
		Scotorepens	greyii	Australia (SEQ)	2006–2009	0, 0–85.3 (1)
			rueppellii	Australia (SEQ)	2006–2009	0, 0–85.3 (1)
			spp.	Australia (SEQ)	2006–2009		100, 90.2–100 (24)
				Australia (WA)	2000–2002		0, 0–85.3 (1)
		Vespadelus	finlaysoni	Australia (WA)	2000–2002		0, 0–85.3 (1)
			pumilus	Australia (SEQ)	2006–2009	25, 2.8–71.6 (4)
			troughtoni	Australia (FNQ)	2006–2009	0, 0–7.7 (31)	16.1, 6.4–31.8 (31)

1Locations within Australia: central Queensland (CQ), far-north Queensland (FNQ), south-east Queensland (SEQ), Northern Territory (NT) and Western Australia.

2Tested using universal coronavirus RT-PCR (Poon et al. 2005). Prevalence, 95% confidence interval (No. tested).

3Tested using SARS coronavirus crude antigen ELISA developed by Yu et al. (2008). Prevalence, 95% confidence interval (No. tested).

4 Rhinolophus aurantia from the Northern Territory, Australia yielded a novel betacoronavirus most closely related to a coronavirus identified in Hipposideros caffer ruber from Ghana.

5 Pteropus alecto from south-east Queensland, Australia yielded a second novel betacoronavirus genotype most closely related to a coronavirus identified in Rousettus aegyptiacus from Kenya and Cynopterus brachyotis from the Philippines.

6 Miniopterus australis in central Queensland, south-east Queensland and far-north Queensland, Australia; M. schreibersii in south-east Queensland and the Northern Territory, Australia and Rhinolophus megaphyllus from far-north Queensland, Australia yielded an alphacoronavirus that shared greater than 99% nucleotide sequence identity with the ICTV reference virus Miniopterus bat coronavirus HKU8.

7 Myotis macropus and Vespadelus pumilus from south-east Queensland, Australia yielded a novel alphacoronavirus genotype most closely related to a putative coronavirus species identified in Pipistrellus kuhlii from Italy and Spain.

Detection of Coronavirus RNA and Antibodies

Coronavirus RNA was detected in faecal, intestinal or anal swab samples of seven species of Australian bats from five families (Table 1). RNA-positive bats were found in the three broad sampling locations in Queensland and the Northern Territory, but not in Western Australia. Viral RNA prevalence in the seven positive species ranged from 1/126 (0.8%, 95% CI 0.1–3.6%) in Rhinonicteris aurantia to 14/30 (46.7%, 95% CI 29.8–64.1%) in Miniopterus australis. Limited sample size (and consequent wide 95% confidence intervals) constrained statistical confidence in prevalence estimates for some species in some locations.

In the 30 bats from which faecal samples, blood and oropharyngeal swabs were collected, coronavirus RNA was detected in 11 faecal samples and two oropharyngeal samples, the latter from bats with RNA-positive faecal samples. No blood samples yielded viral RNA.

Coronavirus antibodies were detected in 23 species of bats from Australia, East Timor, Indonesia, Malaysia and Papua New Guinea (Table 1). In Australia, antibody-positive bats were found in all five broad sampling locations, including remote Western Australia. Further, anti-coronavirus antibodies were detected in five of six species in which coronavirus RNA was detected. The seventh RNA-positive species did not have serum available for testing. As with viral RNA prevalence, seroprevalence varied markedly between species and locations: the highest seroprevalence estimate was 100% on three occasions, namely Scotorepens spp. from south-east Queensland (24/24, 100% 95% CI 90.2–100%), Scotophilus spp. from Malaysia (4/4, 95% CI 55.5–100%) and Miniopterus australis from Western Australia (1/1, 95% CI 14.7–100%); the lowest detected seroprevalence was 1/63 (1.6%, 95% CI 0.2–7.2%) in Macroderma gigas in the Northern Territory, Australia.

Sequencing and Phylogenetic Analysis

The phylogenetic analyses identified four different coronavirus genotypes amongst the sequenced samples. Two genotypes clustered in the genus Alphacoronavirus and two clustered in the genus Betacoronavirus. Significantly, three of the four genotypes shared less than 90% nucleotide sequence identity with the most closely related known coronaviruses. The tree with the highest log likelihood is shown in Fig. 2. The percentage of trees in which the associated taxa clustered together is shown next to the branches. The tree is drawn to scale, with branch lengths measured in the number of substitutions per site. The final consensus sequences, derived from sense and anti-sense primers, were deposited in GenBank under accessions numbers EU834950-EU834956.

Figure 2

Phylogenetic analysis of nucleotide sequences identified in Australian bats. Phylogenetic analysis of the 440 bp sequence from the RNA-dependent RNA polymerase (RdRp) gene from ten representative sequences from this study and 33 other sequences was conducted using the Maximum Likelihood method. The tree with the highest log likelihood (−16,168.6385) is shown. Values to the left of branches are the percentage of trees in which the associated taxa clustered. The tree is drawn to scale, with branch lengths indicating the number of substitutions per site. Coronavirus sequences identified in this study are shown in bold. Square brackets are used to identify coronavirus species and genus groups. Coronavirus nomenclature: Host species/country of origin/laboratory identification/year collected (GenBank accession).

The most widespread genotype, both in terms of host species and geographical location, was an alphacoronavirus that shared greater than 99% nucleotide sequence identity with the ICTV reference virus Miniopterus bat coronavirus HKU8. This genotype was identified in M. australis in central Queensland in 1996, and in south-east Queensland and far-north Queensland between 2006 and 2008; in M. schreibersii from south-east Queensland and the Northern Territory sampled between 2006 and 2008; and in a single Rhinolophus megaphyllus from far-north Queensland sampled in 2007 (Table 1).

A second and novel alphacoronavirus genotype was identified in Myotis macropus and Vespadelus pumilus from south-east Queensland (Table 1; Fig. 2). This nucleotide sequence shares only 89% identity with any reported coronavirus sequence, and is most closely related to a putative coronavirus species identified in Pipistrellus kuhlii from Italy and Spain (Lelli et al. 2013).

A novel betacoronavirus identified in a single Rhinonicteris aurantia from the Northern Territory was most closely related to one identified in Hipposideros caffer ruber from Ghana (Table 1; Fig. 2), although the two sequences had less than 81% identity. A second novel betacoronavirus genotype was identified in Pteropus alecto in south-east Queensland; the most closely related sequences are coronaviruses detected in Rousettus aegyptiacus from Kenya and Cynopterus brachyotis from the Philippines (Table 1; Fig. 2), although nucleotide identity was less than 87%.

Discussion

This is the first report of detection of coronavirus RNA and anti-coronavirus antibodies in bats in the Australasian region, and illustrates that bats in the region are host to a range of coronaviruses. The prevalence of each varied across species and location, with the prevalence of either viral RNA or antibodies being relatively high in some populations. This variability highlights the complexity of exposure, infection and immune dynamics over space and time. In two species with relatively large sample sizes (Miniopterus australis and M. schreibersii in south-east Queensland), 25% of 154 individuals and 27% of 238 individuals, respectively, were positive for coronavirus RNA. Studies elsewhere have shown similar prevalence rates (Tsuda et al. 2012; Suzuki et al. 2014); others have reported much lower prevalence (Anthony et al. 2013). Seroprevalence in our study was also high in some populations, including those that were negative for viral RNA. This antibody-positive, viral RNA-negative finding most plausibly reflects previous infection. Alternatively, virus excretion associated with current infection may have been occurring at a concentration below the detection level of the RT-PCR assay. In some cases, the limited sample size may have constrained detection of viral RNA in the sampled population. Although there have been fewer serological studies than molecular studies in bats, others have also found high coronavirus antibody prevalence in some bat species, albeit using a different serological assay (Tsuda et al. 2012). Indeed, given that our ELISA was based on SARS-CoV antigens, it may be that broader coronavirus seroprevalence is underestimated.

Our molecular analysis was based on highly conserved 440-bp sequence of the coronavirus RdRp gene, a common length and site for sequence data analyses conducted in many other bat coronavirus studies. Although constrained by the short length of sequence, we were able to identify four different genotypes of coronavirus, three of which (an alphacoronavirus and two betacoronaviruses) are potentially new species, having less than 90% nucleotide sequence identity with the most closely related described viruses. To obtain a surrogate estimation of taxonomy, Drexler et al. (2010) proposed a scheme based on 816-bp sequence of the RdRp gene to form RdRp grouping units (RGU), defining species separation in mammalian coronaviruses, with >4.8% amino acid sequence distance for alphacoronaviruses and >6.3% amino acid sequence distance for betacoronaviruses. While our field and laboratory work predated this publication and thus precluded use of the approach, based on nucleotide sequence identity, the three novel coronavirus genotypes identified in our study may represent new species. Regardless, our findings support the hypothesis of Woo et al. (2006) that bats host a large diversity of novel coronaviruses.

Despite targeting rhinolophid bats, the putative natural host taxa of SARS coronaviruses (Lau et al. 2005a; Li et al. 2005), and testing faecal samples from more than 500 individuals from multiple locations and detecting alphacoronavirus RNA, we did not detect any SARS-like betacoronaviruses. We did detect anti-coronavirus antibodies in a small percentage of R. megaphyllus, but we believe this more plausibly reflects alphacoronavirus exposure, and thus have no robust indication that SARS-like coronaviruses circulate in Australian bats.

Previous studies have reported that bat coronaviruses have a relatively narrow host range, independent of location (Poon et al. 2004; Chu et al. 2006; Tang et al. 2006; Woo et al. 2006; Gloza-Rausch et al. 2008; Pfefferle et al. 2009; Drexler et al. 2014). Our study supports this host species-virus association with the identification of Miniopterus bat coronavirus HKU8 (previously reported in Miniopterus species in China, Hong Kong and Bulgaria) in M. australis and M. schreibersii from Australia, some 7,000 km from Hong Kong and some 15,000 km from Bulgaria. Similarly, the novel betacoronavirus genotype from Rhinonicteris aurantia is most closely related to betacoronavirus sequence identified in Ghana from bats of the same family Hipposideridae (Pfefferle et al. 2009), although the lower sequence identity should be noted. Finally, our betacoronavirus genotype from P. alecto is most closely related to Bat coronavirus HKU9, identified in Rousettus spp from China and Kenya, and Cynopterus brachyotis from the Philippines, all of which belong to the family Pteropodidae. Notwithstanding the limited number of sequences and modest size of the amplicon evaluated, our findings support the contention that, regardless of location, bats belonging to the same species, genus or family are likely to host closely related coronaviruses. With the recent revision of the taxonomy of bats using comparative-method and molecular studies (Hutcheon and Kirsch 2006), the suborder Pteropodiformes now comprises, amongst others, bats from the families Hipposideridae, Rhinolophidae and Pteropodidae. To date, betacoronaviruses have been identified predominantly from bats of these families, so it appears that the host-virus association of bat coronaviruses may extend to bats of the same suborder. Thus, as yet unidentified betacoronaviruses may still be found in other species in other families (Craseonycteridae, Megadermatidae and Rhinopomatidae) in the suborder Pteropodiformes.

Our finding that an alphacoronavirus Miniopterus bat coronavirus HKU8 detected in a Rhinolophus megaphyllus bat was genetically identical to that identified in 50% of M. australis bats at the same roost site (Table 1) suggests that a cross-species transmission event has occurred. The detection of viral RNA in only one of the 58 R. megaphyllus sampled at that site may indicate that only sporadic spillover is occurring, and intra-species transmission in R. megaphyllus is not (yet) occurring. The genus Rhinolophus is hypothesised to be more likely to foster host shifts and foster the emergence of coronaviruses than other genera (Cui et al. 2007), and our observation is consistent with this. Similarly, our detection of coronavirus RNA in one of the 126 R. aurantia in the absence of any seropositivity in the population may also suggest transmission from another bat species. The sampled R. aurantia population was co-roosting with populations of Hipposideros ater and Macroderma gigas bats (Fig. 3), both of which were positive for anti-coronavirus antibodies, and although negative for coronavirus RNA at the time of sampling, constitute a plausible source of infection. Further investigation of these two species, beyond the scope of this study, is necessary to test this hypothesis. The alternative interpretation of the single RNA detection in both R. megaphyllus and R. aurantia is that the PCR results are spurious. However, we believe that the rigour of our field and laboratory techniques make cross-contamination of samples highly unlikely.

Two methods of evolution have been proposed to explain the genetic diversity of bat coronaviruses and their host specificity (Cui et al. 2007; Vijaykrishna et al. 2007). Divergent evolution requires the cross-species transmission of a common ancestor virus and subsequent adaption and establishment in the new host. To account for the identification of related coronaviruses in related species of bats in different locations throughout the world, divergent evolution requires the global distribution of each newly diverged coronavirus. This process is theoretically possible for those bat species that exhibit long-range nomadic movement (Breed et al. 2010). The second evolutionary method is the co-evolution of bat species and their coronaviruses (Cui et al. 2007), whereby the divergence of each bat species is mirrored by the divergence of the coronavirus it hosts. This method accounts for the diversity, relatedness and global distribution of bat coronaviruses, but would require the age of bat coronaviruses to parallel that of the most common bat ancestor, i.e. 65 million years (Churchill 2008). In addition, co-evolution alone does not explain the presence of different coronavirus genera in the same species or genus; for example, hipposiderid and rhinolophid bats hosting both alpha- and betacoronaviruses (Woo et al. 2006; Pfefferle et al. 2009), which intuitively appear to require episodes of host-shifting. Our findings are consistent with a scenario wherein the current diversity of coronaviruses and host specificity in Australasian bats reflects co-evolution with the occasional host shift, most plausibly by hipposiderid and/or rhinolophid bats (Cui et al. 2007). The significance of host geography and host ecology on virus evolution and host-virus interactions has recently been canvassed for bat rabies virus (Streicker et al. 2010, 2012). The concepts may be equally relevant for coronaviruses across the geographic and taxonomic spectrum, and warrant future investigation beyond the scope of this regional study.

The majority of coronaviruses previously reported in bats were detected in faecal samples or rectal swabs indicating a predominantly enteric tropism (Lau et al. 2005a; Poon et al. 2005; Chu et al. 2006; Tang et al. 2006; Dominguez et al. 2007; Lau et al. 2007). From the (albeit limited) number of bats from which additional sample types were tested in this study (n = 30), coronavirus RNA was detected in oropharyngeal swabs (n = 2) only from animals that were RNA-positive in faecal samples (n = 11), supporting a predominantly enteric tropism of bat coronaviruses. Corroborative studies are needed given our limited sample size. Lack of viral RNA detection in blood samples indicates that blood is not useful for the detection of bat coronaviruses.

Conclusion

This study identified coronavirus RNA and antibodies in a range of bat taxa in Australia and the neighbouring countries of East Timor, Indonesia, Malaysia and Papua New Guinea. Four different coronavirus genotypes were identified, including one alphacoronavirus and two betacoronaviruses that are evidently novel. We also detected possible cross-species transmission of bat coronaviruses, and provide further support for the apparent enteric tropism of these viruses. The findings advance our knowledge of the diversity of coronaviruses in bats. This diversity, the global distribution of bats, and the propensity of coronaviruses to successfully cross species barriers, suggests MERS and SARS-like coronaviruses may not be the only bat coronaviruses that can cause disease outbreaks in humans, and further highlights the need for ongoing investigation of bats and their coronaviruses.