four bat guano samples were collected from known roost sites on the remote offshore island of whenua hou ( codfish island ) ( 4647s , 16738e ) , which is situated at the southern coast of new zealand . it has special conservation status for the protection of endangered species , and public access is not permitted ( 9 ) . bat guano was held at 4c during transport ( < 48 h ) , resuspended in 2 ml phosphate - buffered saline , and subjected to centrifugation at 6,000  g for 5 min . rna was extracted from 400 l of supernatant by using the iprep purelink virus kit ( life technologies , carlsbad , ca , usa ) and eluted into 50 l reverse transcription pcr molecular - grade water ( ambion , austin , tx , usa ) . metagenomic sequencing was then conducted for 1 of the samples by using an illumina miseq instrument ( new zealand genomics ltd . , massey genome service , massey university , palmerston north , new zealand ) after a series of steps involving dnase i treatment , reverse transcription , multiple displacement amplification ( qiagen , valencia , ca , usa ) , and illumina truseq library preparation ( new zealand genomics , ltd . ) . a total of 10,749,878 paired - end sequence reads of 250 bp were generated and assembled into contigs by using velvet 1.2.07 ( 10 ) . assembled contigs were searched for viral sequence by comparison to the nonredundant nucleotide database in genbank ( downloaded march 2013 ) by using the nucleotide basic local alignment search tool ( national center for biotechnology information , bethesda , md , usa ) . forty - six contigs showed similarity to known genus alphacoronavirus sequences ( table ) . ( raw sequence data are available on request from the authors . ) * blastn , nucleotide basic local alignment search tool ( http://blast.ncbi.nlm.nih.gov/blast.cgi?program=blastn&page_type=blastsearch&link_loc=blasthome ) . virus sequence from the genus alphacoronavirus was confirmed in all 4 original guano samples by using a specific reverse transcription pcr based on the metagenomic data specific for 582 bp of the rna - dependent rna polymerase ( rdrp ) gene ( genbank accession nos . the rdrp sequence was identical in all 4 guano samples , and the closest relative was bat coronavirus hku8 ( genbank accession no . phylogenetic analysis was performed for rdrp ( figure 1 ) and for the spike protein , as derived from metagenomic data ( genbank accession no . we propose that these data support identification of a new alphacoronavirus , which has been designated as mystacina bat cov / new zealand/2013 . phylogenetic tree showing genetic relatedness of rna - dependent rna polymerase amino acid sequences for mystacina sp . bat coronavirus ( cov)/new zealand/2013 ( shown in boldface ) with those of known coronaviruses . evolutionary history was inferred for 183 informative amino acid sites by using the maximum - likelihood method based on the whilan and goldman model with gamma distribution in mega 5.05 software ( www.megasoftware.net ) . bootstrap values are calculated from 1,000 trees ( only bootstrap values > 50% are shown ) . tgev , transmissible gastroenteritis cov ; prcv , porcine respiratory cov ; sars , severe acute respiratory syndrome . bat coronavirus ( cov)/new zealand/2013 ( shown in boldface ) with those of known coronaviruses . evolutionary history was inferred for 492 informative amino acid sites by using the maximum - likelihood method based on the whilan and goldman + f model with gamma distribution and invariant sites in mega 5.05 software ( www.megasoftware.net ) . bootstrap values are calculated from 1,000 trees ( only bootstrap values > 50% are shown ) . tgev , transmissible gastroenteritis cov ; prcv , porcine respiratory cov ; sars , severe acute respiratory syndrome . the discovery of unknown coronaviruses provides information for a model of coronavirus evolution ( 11 ) and contributes to understanding the process of disease emergence , as in detection of middle east respiratory syndrome coronavirus ( 12 ) . the alphacoronavirus identified in this study from m. tuberculata bat guano is unique in respect to the extreme geographic and evolutionary isolation of the host bat species , which along with c. tuberculatus bats , has been separated from all other mammalian species for 1 million years . the current estimate for a common ancestor for all coronaviruses is 8,100 bce ( 13 ) . to be consistent with this estimate , mystacina bat coronavirus would need to have been introduced to bats on whenua hou within the past 800 years since humans first arrived on this island ( given that the island had no other mammals before this time ) ( 11 ) or is extant to modern alphacoronavirus phylogenetic radiation ( genesis and expansion ) . apart from humans , only 2 other terrestrial mammals have ever inhabited whenua hou : the brushtail possum ( trichosurus vulpecula ) and the polynesian rat ( rattus exulans ) , both of which were eliminated from the island in the late 1980s ( 11 ) ; neither mammal has been reported as a host of alphacoronaviruses . members of the genus alphacoronavirus infect only mammals . thus , an avian origin for this virus is unlikely . an alternative theory of an ancient origin for all coronaviruses has recently been proposed that involves an alternative evolutionary molecular clock analysis , which places the most recent common ancestor many millions of years ago ( 14 ) . the discovery of mystacina bat cov / new zealand/2013 virus could lend support to such a theory ; despite potentially millions of years of isolation , it has diverged relatively little from other extant alphacoronaviruses , as shown by the close relationship of the rdrp and spike protein genes to those of other extant alphacoronaviruses ( figures 1 , 2 ) . an expanded survey for mystacina bat coronavirus in mammals in new zealand and subsequent characterization of viral genomes would provide further insights into the origin of coronaviruses . only a small number of conservation staff handle these bats , these staff use standard personal protective equipment and work practice . staff are also offered prophylactic rabies vaccination as a precautionary measure , even though new zealand is free from rabies . the genus alphacoronavirus includes several human and animal pathogens , but on the basis of phylogenetic data in this study , it is not possible to estimate the risk posed by mystacina bat coronavirus to human or animal health . this country is free from many human and animal diseases , such as rabies and foot - and - mouth disease , and infections with human arboviruses because of recent colonization by humans and strict biosecurity border controls ( 15 ) . thus , indigenous wildlife in new zealand has generally been viewed as an almost sterile and unique biosphere . given detection of this coronavirus , more thorough characterization of the ecology of viruses and other microorganisms in native wildlife should be considered to fulfill conservation needs and further safeguard human and domestic animal health against cross - species transmission .
because of recent interest in bats as reservoirs of emerging diseases , we investigated the presence of viruses in mystacina tuberculata bats in new zealand . 
 a novel alphacoronavirus sequence was detected in guano from roosts of m. tuberculata bats in pristine indigenous forest on a remote offshore island ( codfish island ) .