Samba virus: a novel mimivirus from a giant rain forest, the Brazilian Amazon
© Campos et al.; licensee BioMed Central Ltd. 2014
Received: 19 February 2014
Accepted: 1 May 2014
Published: 14 May 2014
The identification of novel giant viruses from the nucleocytoplasmic large DNA viruses group and their virophages has increased in the last decade and has helped to shed light on viral evolution. This study describe the discovery, isolation and characterization of Samba virus (SMBV), a novel giant virus belonging to the Mimivirus genus, which was isolated from the Negro River in the Brazilian Amazon. We also report the isolation of an SMBV-associated virophage named Rio Negro (RNV), which is the first Mimivirus virophage to be isolated in the Americas.
Based on a phylogenetic analysis, SMBV belongs to group A of the putative Megavirales order, possibly a new virus related to Acanthamoeba polyphaga mimivirus (APMV). SMBV is the largest virus isolated in Brazil, with an average particle diameter about 574 nm. The SMBV genome contains 938 ORFs, of which nine are ORFans. The 1,213.6 kb SMBV genome is one of the largest genome of any group A Mimivirus described to date. Electron microscopy showed RNV particle accumulation near SMBV and APMV factories resulting in the production of defective SMBV and APMV particles and decreasing the infectivity of these two viruses by several logs.
This discovery expands our knowledge of Mimiviridae evolution and ecology.
KeywordsMimiviridae DNA virus Giant virus NCLDV Virophage Amazon Brazil
The discovery of the Acanthamoeba polyphaga mimivirus (APMV), arguably the most elusive member of the nucleocytoplasmic large DNA virus (NCLDV) group and the first discovered member of the Mimiviridae family, revived discussions regarding the evolution and origin of the viruses, as well as the differentiation between viruses and living organisms . The complexity of NCLDVs in terms of genome size, particle size and metabolic capabilities (such as their role in photosynthesis and apoptosis) has challenged many concepts in virology . However, it was the discovery of APMV that spotlighted NCLDVs . This novel member of the NCLDV group, which belongs to the proposed Megavirales order, is extremely large and complex and contains genes related to translational activity, which were hitherto considered to be exclusive to cellular organisms .
The family Mimiviridae is comprised of double stranded DNA viruses up to 750 nm of diameter with genomes containing up to 1.2 Mb. The mimiviruses are some of the most complex viruses known to date and are important members of the NCLDV group . One of this family’s key members is the Mimivirus genus, whose type species is APMV. APMV was isolated in 1992 from a water cooling tower at a hospital in Bradford, England and was investigated as a putative etiological agent of pneumonia . The APMV particle is composed of a core, internal membrane, a capsid and external fibrils. The capsid has semi-icosahedral pentagonal symmetry and with a star-shaped structure called the star gate [7, 8].
To date, members of the Mimiviridae family have been isolated in England, France, Tunisia, Chile and a few other countries [3, 4, 9]. Additionally, DNA from these viruses was identified in the Sargasso Sea and other ocean samples using metagenomic approaches [10–12]. While APMV was isolated from a cooling tower at a hospital in England, Megavirus chilensis (MCHV), a putative new species of the Mimiviridae family, was isolated on the coast of Chile, indicating that members of this family can be found in a range of environmental conditions [3, 4].
Acanthamoeba castellanii mamavirus (ACMV), a strain of APMV, was isolated from a water cooling tower at a hospital in France, together with another virus, Sputnik virus (SNV) . SNV decreased the infectivity of ACMV in cultured Acanthamoeba castellanii, leading to its classification as the first virophage . Since then, other viruses with similar biological activity to SNV have been isolated from other NCLDVs, consolidating the emerging concept of virophages [14, 15].
Mimiviridae family viruses have been isolated in different countries from aquatic environments that display as wide range of temperatures and salinity [3, 4, 8, 9]. However, to date, only one Mimivirus has been isolated in the Americas . In this work we report the isolation, biological characterization, genome sequencing and annotation of the giant Samba virus (SMBV) and its associated virophage Rio Negro virophage (RNV), both isolated from the Brazilian Amazon, known to be one of the most biodiverse ecosystems in the world.
Results and Discussion
Collection area data
SMBV: isolation of a giant virus from the Negro River in the Amazon, Brazil
To characterize this potential new giant virus, SMBV was grown and purified in A. castellanii as described previously , and A. castellanii cells infected with SMBV at a TCID50 per cell rate of 10 were analyzed by electron microscopy (EM). Uninfected amoebae were used as controls. After seven hours of infection, EM images revealed the presence of giant viruses with multi-layered capsids covered with fibrils (Figure 2B). The capsids averaged 352 nm in diameter, the fibrils averaged 112 nm in length, and the average diameter of the particles was 574 nm (Figure 2C), making SMBV the largest virus ever isolated in Brazil (a total of 50 particles images were analyzed) and the first mimivirus isolated in this country. Actually, the size of the particles is likely to be significantly larger since chemical preparation might be related to particles shrinkage. In some of the images we detected a hypothetical star-gate structure, which has been described in other giant viruses (Figure 2D). We next observed purified virus using a light microscope. Remarkably, we were able to detect several particles at 1000× magnification on agarose surface, similarly to APMV.
Isolation and characterization of RNV, the virophage associated with SMBV
To analyze the influence of the virophage on giant virus replication, purified SMBV were diluted and filtered. The aliquoted virophage solution was stored at -80°C until use. An A. castellanii monolayer was co-infected with APMV (TCID50 per cell rate of = 10) and 100 μl of the solution containing the undiluted virophage isolated from SMBV. After 16 hours, we analyzed the cells by EM. Remarkably, virophage particles derived from SMBV were observed in association with APMV particles during viral assembly (Figure 4F). We also observed defective APMV particles similar to those seen in SMBV and virophage co-infected amoebas. Large areas of virophage accumulation covering more area than the APMV factories themselves were observed in some of the co-infected amoebas (Figure 4G).
To determine which viruses SMBV is most closely related to, we performed the following analyses:  we constructed phylogenetic trees by aligning the ribonucleotide reductase (Figure 6B) and helicase gene (data not shown) sequences from SMBV with GenBank sequences and  we constructed a Venn diagram showing a presence-absence analysis of ORFs from SMBV and related viruses (Additional file 3: Figure S3). Both phylogenetic trees show that SMBV clusters with members of Megavirales order group A (MGA), which includes APMV (the prototype of the family) and ACMV (Figure 6B). MGA was most closely related to Megavirales group B (MGB), which includes Acanthamoeba polyphaga moumouvirus (APMOUV), and was distinct from Megavirales group C, which includes MCHV and other viruses. Other NCLDV family members were also included in the analysis, and their position in the phylogenetic tree generated in this study corroborates previous phylogenetic studies. The phylogenetic trees were constructed by MEGA 4.1 using neighbor joining, maximum parsimony and other methods with 1,000 bootstrap replicates (n. differences model) based on the ribonucleotide reductase predicted amino acid from the SMBV and other nucleocytoplasmic large DNA viruses. The same tree topology was observed for all phylogenetic approaches.
Next, a Venn diagram (Additional file 3: Figure S3) was constructed using the predicted ORFs from SMBV and the ORF data set from its closest relatives, including APMV and ACMV (from the MGA group), APMOUV (from the MGB group) and MCHV (from the MGC group). The SMBV ORFs were aligned with the ORF data set of each viral counterpart using the Blastn all-against-all method feature of Blastall 2.2.9. Only hits with an e-value ≤ 10-5 were considered valid. The diagram (Additional file 3: Figure S3) shows that nine SMBV ORFs did not appear in other viral genomes. Of these, ORF-L331 appears to contain a signal-peptide cleavage site, and ORF-R518 exhibits a transmembrane domain. The remaining ORFans had no known function or domain. SMBV shares 925 ORFs with APMV, 909 with ACMV, 346 with MCHV and 312 with APMOUV. SMBV, APMV and ACMV (MGA strains) share 503 ORFs. SMBV, MGA and APMOUV (MGB) share 61 ORFs, while SMBV, MGA and MCHV (MGC) shared 93 ORFs. All of the viruses analyzed share a core repertoire of 249 ORFs (Additional file 3: Figure S3).
RNV virophage sequence analysis
The RNV capsid gene was also sequenced and analyzed (primers: 5’ATGTCTAATTCAGCTATTCCTCTTA3’ and 5’TCACATTTTAAGTTCTTTTCTCAAT3”). We then manually aligned the sequence with conserved virophage sequences from GenBank and used Modeltest software to determine which model of evolution was most appropriate for our analysis. Phylogenetic trees based on the capsid gene sequence were constructed using MEGA 4.1 with neighbor joining, maximum parsimony and other methods with 1,000 bootstrap replicates. The RNV capsid sequence was deposited in GenBank (KJ183141). Our results show that the RNV capsid gene shares 100% identity with SNV and high identity with other SNV genes. Therefore, RNV clusters with SNV-like viruses and does not cluster with Mavirus or the Organic Lake Virophage (Additional file 4: Figure S4). Unfortunately, it was not possible to recover RNV contigs from SMBV sequencing data.
Giant viruses exhibit strikingly large genome and particle sizes, and have shed light on the evolution of large DNA viruses [3, 4, 17–25]. This study describes the discovery, isolation and characterization of Samba virus, a novel mimivirus, from the Negro River in the Brazilian Amazon. SMBV is phylogenetically related to Megavirales order group A, which comprises APMV and ACMV, and contains one of the the largest genome described for this group.
However, despite its close relationship to APMV, the SMBV genome is 50,000 bp longer than APMV and is unique for its ORF content. The SMBV and APMV genome alignment showed a high degree of synteny among the conserved genes, although we did observe some inversions (Additional file 2: Figure S2). Many of the SMBV ORFans exhibited homology to APMV sequences, although they did not appear to correspond to ORFs in the APMV genome. This may be explained by the use of different bioinformatics tools for ORF prediction during the genome annotation. The difference of 50,000 bp between the SMBV and APMV genomes is most likely due to intergenic regions and/or ORF size variation. Regarding the SBMV cycle of life, our data showed viral entry by phagocytosis and morphogenesis in association with the viral factories. All steps of viral morphogenesis resembled those described previously to APMV and other mimiviruses [3, 4].
The presence of a Mimiviridae family virus and virophage in the Negro River correlates well with previous studies that indicate the presence of these viruses in aquatic ecosystems [10, 11, 24, 25]. SMBV appeared to tolerate the virophage better than APMV, although we have not yet been able to generate virophage-free cultures. Nevertheless, these results suggest SMBV has developed mechanisms to circumvent the virophage. However, future studies are necessary to confirm this hypothesis. Interestingly, recent data suggested that virophages may be important for control of giant viruses and amoeba populations [26, 27].
The discovery and characterization of SMBV and its virophage raise new questions regarding the role of these viral agents in microbial ecology. The Amazon rain forest contains a striking diversity of flora and fauna, although little is known regarding the virosphere in this environment. The discovery of SMBV in the Amazon corroborates that these fascinating viruses are ubiquitous and that their isolation and characterization can yield important insights into their life cycle and complexity.
Acanthamoeba castellanii (ATCC 30010) was kindly provided by the Laboratório de Amebíases (Departamento de Parasitologia, ICB/UFMG), and cultivated in PYG medium (Visvesvara & Balamuth, 1975) supplemented with 7% fetal bovine serum (FBS) (Cultilab, Brazil), 200 U/mL penicillin (Cristália, Brazil), 50 μg/mL gentamycin (Sigma, USA) and 2.5 μg/mL amphotericin B (Sigma, USA) at 28°C.
APMV was obtained from a cooling tower at a hospital in Bradford (England) in 1993 and characterized in 2003 . APMV was used as a control for the molecular and biological assays. For viral replication, A. castellanii cultures were inoculated at a TCID50 per cell rate of 0.1 and incubated in sealed bottles at 32°C.
Sample collection and virus isolation
The water samples were collected from Negro River, Amazonas and stored at 4°C overnight. The sample collections were performed with permission of Instituto Chico Mendes (ICM) – protocol numbers: 34293-1 and 33326-2. The field studies did not involve endangered or protected species. Then, 500 μl of each sample was added to 4.5 mL of autoclaved rice and water medium made with 40 rice grains in 1 liter of water . The samples were stored for 20 days in the dark at room temperature , then 5 × 103A. castellanii trophozoites were added, and the samples were re-incubated under the same conditions for 10 days. After the enrichment process, samples were pooled in groups of five (totaling seven pools), and filtered through a 1.2 μm membrane to remove impurities and a 0.2 μm membrane to retain giant viruses. The samples were then subjected in parallel to real-time PCR and to viral isolation in A. castellanii.
The virus titration was performed in 96-well plates containing approximately 4 × 104A. castellanii per well in 200 μL of PYG medium supplemented with 7% FCS. The viral samples were serially diluted (from 10-1 to 10-9) in 50 μL of PBS, 150 μL of PYG medium with 10% FCS was added to each well, and the plates were incubated. CPEs such as rounding, loss of motility and trophozoite lysis were monitored daily in each well. After five days of incubation, the titer (TCID50) was calculated as described by Reed and Muench .
To purify the giant viruses, cell lysates were centrifuged at 900 g for 5 minutes at 4°C. The supernatant was transferred to a fresh tube, and the pellet was subjected to three cycles of freezing and thawing to release virus trapped in the trophozoites. The lysate was homogenized in 10 mL of PBS and subjected to an additional two rounds of 50 homogenization cycles in a Dounce (Wheaton, USA). The supernatant and cell lysates were then filtered through a 2 μm filter (Millipore, USA). This filtrate was slowly dripped over 10 mL of a 22% sucrose solution (Merck, Germany) and ultra-centrifuged in a Sorvall Combi at 14,000 rpm for 30 minutes at 4°C. The pellet was homogenized in 500 μL of PBS. Aliquots of the virus were stored at -80°C and then titrated. To purify the virophages, 5 vials containing 50 μl each of previously purified infectious SMBV particles were diluted in 20 mL of PBS and filtered through 0.2 μm filters, which retain the giant viruses but not the virophages. The flow-through was collected and used in the biological assays.
The infectivity assays were performed by inoculating A. castellanii with the viral samples at an m.o.i of 10. The viruses were allowed to adsorb for 1 hour, and then the cells were washed with PBS and incubated at 32°C. At 1, 2, 4, 6, 8, 12, 24 and 48 hours, the cultures were frozen and thawed three times and titrated.
Inhibition of APMV growth after infection with isolated virophages
A. castellanii was infected with APMV at a TCID50 per cell rate of 1 and then superinfected with 100 μl of virophages, undiluted or diluted to 10-9. After 48 hours of infection, the infectivity (TCID50) was determined by observing APMV-induced CPEs and by titration (TCID50) in A. castellanii for five days.
The PCR assays were performed using primers constructed based on the helicase gene from APMV gene (primers: 5’ACCTGATCCACATCCCATAACTAAA3’ and 5’GGCCTCATCAACAAATGGTTTCT3’). The PCR contained 2.0 mM MgCl2, 10 mM nucleotides (dATP, dCTP, dGTP, dTTP), 2 U of Taq DNA polymerase (Promega, USA), 2.0 μl 10× Taq polymerase buffer, 4 mM primers, and 2 μl of the sample in a 20 μl total reaction volume. The amplification was performed according to the conditions recommended for StepOne (Applied Biosystems, USA) with an annealing temperature of 60°C. The PCR-amplified DNA was fractionated on a 1% agarose gel at 100 V and stained with Gel Red (Biotium, USA). The real-time PCR was performed using a commercial mix (Applied Biosystems, USA), primers (4 mM each) and 1 μl of sample in each 10 μl reaction.
Transmission electron microscopy
A. castellanii were infected at a TCID50 per cell rate of 10. Uninfected amoebae were used as controls. After 7 hours of infection, the amoebae were washed twice with 0.1 M phosphate buffer (pH 7.4) and fixed with 2.5% glutaraldehyde (grade I) in 0.1 M phosphate buffer (pH 7.4) (Electron Microscopy Sciences, Germany) for one hour at room temperature. The amoeba monolayer was scraped from the plates and recovered by centrifugation at 900 g for 5 minutes. The amoebae were postfixed with 2% osmium tetroxide and embedded in EPON resin. Ultrathin sections were stained with 2% uranyl acetate and examined using a Tecnai G2-Spirit FEI 2006 transmission electron microscope operating at 80 kV at the Microscopy Center, UFMG, Brazil.
The SMBV genome was sequenced using a 454 platform (Roche). For genome assembly, we used the read mapping approach implemented by the CLC Genomics Workbench 5.5.1 program to generate contig sequences, and the Perl-based genome assembly tool ABACAS.1.3.1 (algorithm-based automatic contiguation of assembled sequences) that yielded a partial genome (scaffold) of 1,213,607 bp (KF959826) (draft). The open-reading frames (ORF) were predicted using the Markov-based methods employed by Glimmer3 and FGENES. We also transferred annotations from a closely related genome using the Rapid Annotation Transfer Tool. The predictions were manually curated, and the ORFs were assigned a final identity. ORFs smaller than 150 bp were ruled out. Finally, 938 ORFs were annotated as putative genes. A gene similarity search was conducted using Blast2GO.
We thank João Rodrigues dos Santos and Gisele Cirilo dos Santos, colleagues from Gepvig and the Laboratório de Vírus, for their excellent technical support. We would also like to thank CNPq, CAPES, FAPEMIG, Pro-Reitoria de Pesquisa da Universidade Federal de Minas Gerais (PRPq-UFMG), and Centro de Microscopia da UFMG.
- Yamada T: Giant viruses in the environment: their origins and evolution. Curr Opin Virol. 2011, 1: 58-62. 10.1016/j.coviro.2011.05.008.PubMedView Article
- Culley AI: Virophages to viromes: a report from the frontier of viral oceanography. Curr Opin Virol. 2011, 1: 52-57. 10.1016/j.coviro.2011.05.003.PubMedView Article
- La Scola B, Audic S, Robert C, Jungang L, de Lamballerie X, Drancourt M, Birtles R, Claverie JM, Raoult D: A giant virus in amoebae. Science. 2003, 299: 2033-10.1126/science.1081867.PubMedView Article
- Arslan D, Legendre M, Seltzer V, Abergel C, Claverie J: Distant Mimivirus relative with a larger genome highlights the fundamental features of Megaviridae. Proc Natl Acad Sci U S A. 2011, 108: 17486-17491. 10.1073/pnas.1110889108.PubMedPubMed CentralView Article
- Legendre M, Arslan D, Abergel C, Claverie JM: Genomics of Megavirus and the elusive fourth domain of Life. Commun Integr Biol. 2012, 5: 102-106. 10.4161/cib.18624.PubMedPubMed CentralView Article
- La Scola B, Marrie TJ, Auffray JP, Raoult D: Mimivirus in pneumonia patients. Emerg Infect Dis. 2005, 11: 449-452. 10.3201/eid1103.040538.PubMedPubMed CentralView Article
- Xiao C, Chipman PR, Battisti AJ, Bowman VD, Renesto P, Raoult D, Rossmann MG: Cryo-electron microscopy of the giant Mimivirus. J Mol Biol. 2005, 28 (3): 493-496. 353View Article
- Zauberman N, Mutsafi Y, Halevy DB, Shimoni E, Klein E, Xiao C, Sun S, Minsky A: Distinct DNA exit and packaging portals in the virus Acanthamoeba polyphaga mimivirus. PLoS Biol. 2008, 13;6 (5): e114-View Article
- Boughalmi M, Saadi H, Pagnier I, Colson P, Fournous G, Raoult D, La Scola B: High-throughput isolation of giant viruses of the Mimiviridae and Marseilleviridae families in the Tunisian environment. Environ Microbiol. 2012, 53: 344-353.
- Ghedin E, Claverie JM: Mimivirus relatives in the Sargasso sea. Virol J. 2005, 2: 62-10.1186/1743-422X-2-62.PubMedPubMed CentralView Article
- Monier A, Larsen JB, Sandaa RA, Bratbak G, Claverie JM, Ogata H: Marine mimivirus relatives are probably large algal viruses. Virol J. 2008, 5: 12-10.1186/1743-422X-5-12.PubMedPubMed CentralView Article
- Williamson SJ, Allen LZ, Lorenzi HA, Fadrosh DW, Brami D, Thiagarajan M, McCrow JP, Tovchigrechko A, Yooseph S, Venter JC: Metagenomic exploration of viruses throughout the Indian Ocean. PLoS One. 2012, 7: e42047-10.1371/journal.pone.0042047.PubMedPubMed CentralView Article
- La Scola B, Desnues C, Pagnier I, Robert C, Barrassi L, Fournous G, Merchat M, Suzan-Monti M, Forterre P, Koonin E, Raoult D: The virophage as a unique parasite of the giant mimivirus. Nature. 2008, 455: 100-104. 10.1038/nature07218.PubMedView Article
- Yau S, Lauro FM, DeMaere MZ, Brown MV, Thomas T, Raftery MJ, Andrews-Pfannkoch C, Lewis M, Hoffman JM, Gibson JA, Cavicchioli R: Virophage control of antarctic algal host-virus dynamics. Proc Natl Acad Sci U S A. 2011, 108: 6163-6168. 10.1073/pnas.1018221108.PubMedPubMed CentralView Article
- Fischer MG, Suttle CA: A virophage at the origin of large DNA transposons. Science. 2011, 332: 231-234. 10.1126/science.1199412.PubMedView Article
- Dare RK, Chittaganpitch M, Erdman DD: Screening pneumonia patients for mimivirus. Emerg Infect Dis. 2008, 14: 465-467. 10.3201/eid1403.071027.PubMedPubMed CentralView Article
- Philippe N, Legendre M, Doutre G, Couté Y, Poirot O, Lescot M, Arslan D, Seltzer V, Bertaux L, Bruley C, Garin J, Claverie JM, Abergel C: Pandoraviruses: amoeba viruses with genomes up to 2.5 Mb reaching that of parasitic eukaryotes. Science. 2013, 341: 281-286. 10.1126/science.1239181.PubMedView Article
- Xiao C, Kuznetsov YG, Sun S, Hafenstein SL, Kostyuchenko VA, Chipman PR, Suzan-Monti M, Raoult D, McPherson A, Rossmann MG: Structural studies of the giant mimivirus. PLoS Biol. 2009, 7: e92-10.1371/journal.pbio.1000092.PubMedView Article
- Claverie JM, Abergel C: Mimivirus: the emerging paradox of quasi-autonomous viruses. Trends Genet. 2010, 26: 431-437. 10.1016/j.tig.2010.07.003.PubMedView Article
- Boyer M, Azza S, Barrassi L, Klose T, Campocasso A, Pagnier I, Fournous G, Borg A, Robert C, Zhang X, Desnues C, Henrissat B, Rossmann MG, La Scola B, Raoult D: Mimivirus shows dramatic genome reduction after intraamoebal culture. Proc Natl Acad Sci U S A. 2011, 108: 10296-10301. 10.1073/pnas.1101118108.PubMedPubMed CentralView Article
- Suzan-Monti M, La Scola B, Barrassi L, Espinosa L, Raoult D: Ultrastructural characterization of the giant volcano-like virus factory of Acanthamoeba polyphaga Mimivirus. PLoS One. 2007, 2: e328-10.1371/journal.pone.0000328.PubMedPubMed CentralView Article
- Krupovic M, Cvirkaite-Krupovic V: Virophages or satellite viruses?. Nat Rev Microbiol. 2011, 9: 762-763. 10.1038/nrmicro2676.PubMedView Article
- Desnues C, Raoult D: Inside the lifestyle of the virophage. Intervirology. 2010, 53: 293-303. 10.1159/000312914.PubMedView Article
- Claverie JM, Grzela R, Lartigue A, Bernadac A, Nitsche S, Vacelet J, Ogata H, Abergel C: Mimivirus and Mimiviridae: giant viruses with an increasing number of potential hosts, including corals and sponges. J Invertebr Pathol. 2009, 101: 172-180. 10.1016/j.jip.2009.03.011.PubMedView Article
- Fischer MG, Allen MJ, Wilson WH, Suttle CA: Giant virus with a remarkable complement of genes infects marine zooplankton. Proc Natl Acad Sci U S A. 2010, 107: 19508-19513. 10.1073/pnas.1007615107.PubMedPubMed CentralView Article
- Desnues C, La Scola B, Yutin N, Fournous G, Robert C, Azza S, Jardot P, Monteil S, Campocasso A, Koonin EV, Raoult D: Provirophages and transpovirons as the diverse mobilome of giant viruses. Proc Natl Acad Sci U S A. 2012, 30: 18078-18083.View Article
- Slimani M, Pagnier I, Raoult D, La Scola B: Amoebae as battlefields for bacteria, giant viruses, and virophages. J Virol. 2013, 87 (8): 4783-4785. 10.1128/JVI.02948-12.PubMedPubMed CentralView Article
- Reed LJ, Muench H: A simple method of estimating fifty percent endpoints. Am J Hyg. 1938, 27: 493-497.
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.