Skip directly to site content Skip directly to page options Skip directly to A-Z link Skip directly to A-Z link Skip directly to A-Z link
Volume 21, Number 4—April 2015
Letter

Bat Coronavirus in Brazil Related to Appalachian Ridge and Porcine Epidemic Diarrhea Viruses

On This Page
Article Metrics
5
citations of this article
EID Journal Metrics on Scopus

Cite This Article

To the Editor: Tadarida brasiliensis (I. Geoffroy, 1824) is a species of free-tailed bat that has resident and migratory populations in Brazil (1). This species has adapted to urban areas, occupying roofs, ceilings, and other human constructions, and often coexists with other bat species and humans (2), enabling epidemiologic risks (3). In recent studies, an alphacoronavirus has been detected in urban bat species Molossus molossus, M. rufus, and Tadarida brasiliensi in Brazil (4,5). Evidence suggests that alphacoronaviruses may use bats as hosts to spread human coronavirus (HCoV) NL63, which originated by evolution of Appalachian Ridge CoV strain 2 (ARCoV.2) (6).

In this study, a total of 20 anal and tracheal swab samples from 10 bats (T. brasiliensis) were collected at the Jequitibás Wood, in Campinas, São Paulo State, Brazil (22°54′31.34′′S 47°02′58.01′′W). We extracted viral genetic material using the RNA Extraction Mini Kit (QIAGEN, Hilden, Germany) and synthesized cDNA using random primers from the High Capacity cDNA Reverse Transcription Kit (Applied Biosystems, Foster City, CA, USA), following the manufacturer’s protocol.

Samples were analyzed by conventional reverse transcription PCR assays using panCoV primers targeting a 215-bp replicase fragment as previously described (7) with slight modifications to include more cycles and less extension time in order to obtain more PCR products. Sequencing reactions on a Pancoronavirus-positive anal swab sample (7) were performed at Central Laboratory of High Performance Technologies in Life Sciences (LaCTAD) at UNICAMP (http://www.lactad.unicamp.br) using an automated sequencer (3730xl DNA Analyzer; Applied Biosystems).

The chromatograms were edited using the program UGENE version 1.14 (UGENE, http://ugene.unipro.ru/forum/YaBB.pl?num=1407749393) and evaluated using Phred scores for base calling. Alignment was made with ClustalW v.2.1 software (http://www.clustal.org) implemented on Linux command line interface, and a similarity matrix was generated with sequences retrieved from the GenBank database. A 144-nt fragment of the replicase gene was obtained after editing, and phylogenetic analysis was performed after determining the best evolution model by using the jModelTest2 software (https://code.google.com/p/jmodeltest2/). Different CoV sequences were included to represent the genera Alpha-, Beta-, Gamma-, and Deltacoronavirus. Clustering with the ARCoV.2 and porcine epidemic diarrhea virus (PEDV) was obtained using the maximum-likelihood (ML) method after 1,000 Shimodaira-Hasegawa–like support values with the general time-reversible model and category approximation in 20 rates category in a gamma distribution (Technical Appendix Figure, panel A) and neighbor-joining methods under Kimura-2-parameter and 1,000 replicates of bootstrap (Technical Appendix Figure, panel B).

Subsequently, metagenomic analysis was made by creating a pool of the 10 bat samples. Samples were resuspended in Dulbecco Modified Eagle Medium (Life Technologies-GIBCO, Grand Island, NY, USA) and filtered through 0.22 μm. The recovered sample was then treated with DNase (Invitrogen, Carlsbad, CA, USA) to remove contaminating DNA and with Proteinase K (Invitrogen) to eliminate inhibitors and to disrupt viral capsids. Samples were then subjected to RNA extraction (QIAGEN) and sent to the sequencing core facility. Sequencing was performed on Illumina HiSeq2500 instrument by using the 2×100 bp kit according to manufacturer’s instructions.

Through these analyses, we obtained 34,409,110 reads, of which 76.47% had quality index ≥30. The contigs were assembled by de novo genome assembly (blastx E-value ≤1–5) (8) generating 10.742 scaffolds: 35 matches for coronaviruses (using the Coronavirus Database, http://covdb.microbiology.hku.hk), 3 matches for PEDV, and 2 matches for HCoV-NL63 (both using the UniProt database, http://www.uniprot.org) (Technical Appendix). The sequences obtained had 87.5% (126/144) nucleotide identity with ARCoV.2, an unclassified alphacoronavirus (GenBank accession no. JX537912) for which a zoonotic role has been suggested (6). Preliminary analysis indicated good coverage of the polymerase region of the ARCoV.2 reference sequence by the reads (quality index >40) by using reference assembly against CoV complete genomes. This finding reinforces the hypothesis of this viral agent in the specimens analyzed. Moreover, molecular assays are underway in our laboratory to elucidate the alternative hypothesis of PEDV presence in bats in Brazil.

In summary, we found that a CoV detected in T. brasiliensis bats in Brazil has close phylogenetic relationships to ARCoV.2 and PEDV. Considering the zoonotic impact of these viral agents on the emergence of new diseases in animal and human populations, we believe that both results may strongly contribute to a better understanding of the molecular eco-epidemiology of these alphacoronaviruses. The reconstruction of their evolutionary history to trace their occurrence in humans and in bat populations as well as in other animals is being conducted to clarify their evolutionary pathway.

Top

Acknowledgment

This work was supported by FAPESP (grant 2011/50919-5) and CNPq (grant 307738/2011-6).

Top

Paulo Vitor Marques Simas1, Ana Caroline de Souza Barnabé1, Ricardo Durães-Carvalho1, Daniel Ferreira de Lima Neto1, Leonardo Cardia Caserta1, Luiza Artacho, Fábio André Facco Jacomassa, Matheus Cavalheiro Martini, Márcia Mercês Aparecida Bianchi dos Santos, Paulo Anselmo Nunes Felippe, Helena Lage Ferreira, and Clarice Weis ArnsComments to Author 

Author affiliations: University of Campinas, Campinas, São Paulo, Brazil (P.V.M. Simas, A.C.S. Barnabé, R. Durães-Carvalho, D.F. Lima-Neto, L.C. Caserta, L. Artacho, M.C. Martini, P.A.N. Felippe); São Paulo State University, Rio Claro, São Paulo (F.A.F. Jacomassa); Federal University of Juiz de Fora, Juiz de Fora, Brazil (M.M.A. Bianchi dos Santos); São Paulo University, Pirassununga, Brazil (H.L. Ferreira).

Top

References

  1. Armstrong  K. Brazilian free-tailed bat (Tadarida brasiliensis). [cited 2014 Dec 12]. Mamm Species. 2008;4:16 http://www.cfr.msstate.edu/wildlife/mammals/pdf/Brazilianfree-tailedbat.pdf.
  2. Wilkins  KT. Tadarida brasiliensis. Mamm Species. Oxford University Press. 1989;331:1–10.
  3. Scheffer  KC, Carrieri  ML, Albas  A, Santos  HCP, Kotait  I, Ito  FH. Rabies virus in naturally infected bats in the state of São Paulo, southeastern Brazil. Rev Saude Publica. 2007;41:38995 . DOIPubMed
  4. Góes  LGB, Ruvalcaba  SG, Campos  AA, Queiroz  LH, De Carvalho  C, Jerez  JA, Novel Bat Coronaviruses, Brazil and Mexico. Emerg Infect Dis. 2013;19:17113 .PubMed
  5. Lima  FE, Campos  FS, Filho  H, Batista  HB, Junior  P, Cibulski  SP, Detection of Alphacoronavirus in velvety free-tailed bats (Molossus molossus) and Brazilian free-tailed bats (Tadarida brasiliensis) from urban areas of Southern Brazil. Virus Genes. 2013;47:1647 . DOIPubMed
  6. Huynh  J, Li  S, Yount  B, Smith  A, Sturges  L, Olsen  JC, Evidence Supporting a Zoonotic Origin of Human Coronavirus Strain NL63. J Virol. 2012;86:1281625 . DOIPubMed
  7. Vijgen  L, Moes  E, Keyaerts  E, Li  S, Van Ranst  M. A pancoronavirus RT-PCR assay for detection of all known coronaviruses. Methods Mol Biol. 2011;454:312. DOIPubMed
  8. Peng  Y, Henry  CM. Leung, S. M. Yiu, Francis Y. L. Chin. IDBA-UD: a de novo assembler for single-cell and metagenomic sequencing data with highly uneven depth. Bioinformatics. 2012;28:14208 . DOIPubMed

Top

Cite This Article

DOI: 10.3201/eid2104.141783

1These authors contributed equally to this article.

Related Links

Top

Table of Contents – Volume 21, Number 4—April 2015

Comments

Please use the form below to submit correspondence to the authors or contact them at the following address:

Clarice Weis Arns, State University of Campinas, Institute of Biology
R. Monteiro Lobato, 255 Cidade Zeferino Vaz, Campinas, Sao Paulo, 13083862
Brazil

Send To

character(s) remaining.

Comment submitted successfully, thank you for your feedback.

Top

Page created: March 17, 2015
Page updated: March 17, 2015
Page reviewed: March 17, 2015
The conclusions, findings, and opinions expressed by authors contributing to this journal do not necessarily reflect the official position of the U.S. Department of Health and Human Services, the Public Health Service, the Centers for Disease Control and Prevention, or the authors' affiliated institutions. Use of trade names is for identification only and does not imply endorsement by any of the groups named above.
file_external