- Extended genome report
- Open Access
Complete genome sequence analysis of Archaeoglobus fulgidus strain 7324 (DSM 8774), a hyperthermophilic archaeal sulfate reducer from a North Sea oil field
Standards in Genomic Sciencesvolume 12, Article number: 79 (2017)
Archaeoglobus fulgidus is the type species of genus Archaeoglobus Stetter 1998, a hyperthermophilic sulfate reducing group within the Archaeoglobi class of the euryarchaeota phylum. Members of this genus grow heterotrophically or chemolithoautotrophically with sulfate or thiosulfate as electron acceptors. Except for A. fulgidus strain 7324 and the candidate species “Archaeoglobus lithotrophicus”, which both originate from deep oil-fields, the other members of this genus have been recovered from marine hydrothermal systems. Here we describe the features of the A. fulgidus strain 7324 genome as compared to the A. fulgidus VC16 type strain. The 2.3 Mbp genome sequence of strain 7324 shares about 93.5% sequence identity with that of strain VC16T but is about 138 Kbp longer, which is mostly due to two large ‘insertions’ carrying one extra cdc6 (cell-cycle control protein 6) gene, extra CRISPR elements and mobile genetic elements, a high-GC ncRNA gene (hgcC) and a large number of hypothetical gene functions. A comparison with four other Archaeoglobus spp. genomes identified 1001 core Archaeoglobus genes and more than 2900 pan-genome orthologous genes.
Archaeoglobus fulgidus strain 7324 was recovered from hot oil-field water originating from a deep oil-well in the North Sea . It shares many features with the A. fulgidus type strain VC16, e.g. dissimilatory sulfate reduction, utilization of lactate and pyruvate as carbon sources, irregular coccoid to disc-shaped cells, and blue-green fluorescence under the UV microscope due to the presence of Coenzyme F420. Strain VC16T was isolated from a shallow marine hydrothermal system at Volcano island, Italy . The complete genome sequence of strain VC16T was reported in 1997 as the third archaeal genome to be fully sequenced  and A. fulgidus has since served as a prototype for studies of archaeal and hyperthermophilic sulfate reduction . Here we report a summary of the features of A. fulgidus strain 7324, together with the description of the complete genomic sequencing and annotation and comparison with the genome of the A. fulgidus type strain and other Archaeoglobus spp.
Classification and features
Genus Archaeoglobus comprises five validly published species; A. fulgidus , Archaeoglobus profundus , Archaeoglobus veneficus , Archaeoglobus infectus , Archaeoglobus sulfaticallidus , and one candidate species termed “ Archaeoglobus lithotrophicus” . All are hyperthermophilic sulfate-reducers capable of heterotrophic or chemolitoautotrophic growth on H2 and CO2. The Archaeoglobus ‘clade’ also encompasses a few non-sulfate reducing anaerobic hyperthermophiles; Geoglobus acetivorans  and “ Geoglobus ahangari ” [11, 12], which are both Fe(III) reducers, and Ferroglobus placidus , which is capable of using ferrous iron, H2 and sulfide as electron donors with nitrate as electron acceptor . Fig. 1 shows the phylogenetic affiliation of all current members of the Archaeoglobaceae family, including strain 7324. All Archaeoglobus species form small irregularly shaped cells. A scanning electron micrograph of A. fulgidus strain 7324 is shown in Fig. 2, revealing a similar cell shape as originally determined by transmission electron microscopy . Strain 7324 has not been phylogenetically characterized by 16S rRNA gene sequencing before, but a wet lab genomic DNA: DNA hybridization with A. fulgidus strain Z, which, like the type strain, was recovered from the Vulcano island , revealed a genome hybridization value of 100% . This close relationship was now confirmed via digital DNA-DNA hybridization  between strains VC16T and 7324 with a GLM-based DDH estimate of 93.9%. All three A. fulgidus strains share common physiological characteristics, like growth from 60 °C to above 84 °C, use of sulfate and thiosulfate as electron acceptors, optimal growth with lactate or pyruvate as carbon sources, and production of trace amounts of methane. Although the optimal growth temperature of strain 7324 was initially determined to 76 °C, we have routinely been cultivating it at 80 °C. In contrast to the other isolates, strain 7324 rapidly lyses after the stationary phase . The main features of the organism are listed in Table 1.
Genome sequencing information
Genome project history
A. fulgidus strain 7324 was chosen for whole genome sequencing because it was isolated from a deep and hot oil reservoir while the closely related type strain, VC-16, was isolated from a shallow marine hot vent. A genome comparison might reveal particular adaptations of strain 7324 to the deep biosphere. The genome project information is given in the Genomes OnLine Database (Gp0102124). The genome sequence is deposited in GenBank (CP006577.1). A summary of the project information is shown in Table 2.
Growth conditions and genomic DNA preparation
A. fulgidus strain 7324 was from our own collection at the University of Bergen. It was cultivated in anaerobic medium containing lactate and sulfate as described previously . The incubation temperature was 80 °C. Genomic DNA was isolated using a modification of the cetyl trimethylammonium bromide method as described .
Genome sequencing and assembly
The genome was sequenced using a combination of Illumina and 454 sequencing platforms. All general aspects of library construction and sequencing can be found at the JGI website . The initial assembly of 454 raw data suggested a contamination of the sequenced sample. Using blast search, all contigs (>500 nt in length) could be assigned either to A. fulgidus or Thermococcus litoralis , an archaeon that shares the same habitat . To overcome this issue, two additional blast searches including all contigs longer than 500 nt were performed against the previously sequenced genome of A. fulgidus VC-16T (NCBI/GenBank:AE000782) and all available genomic sequences of Thermococcus species in Genebank (Dec. 2010). Only sequences in length sharing more than 90% sequence identity with A. fulgidus VC-16T and having no hits in the Thermococcus blast database were kept. A total of 84 Newbler contigs could be assigned to A. fulgidus . Illumina raw reads were assembled to 223 contigs. Both draft assemblies were merged in a hybrid approach using the phred/phrap/consed pipeline . After manual curation, a total of 27 ordered gaps were closed by bridging PCRs at LGC Genomics (Berlin). The final consensus sequence represents a single circular chromosomal element (103× coverage).
Coding genes were predicted by GeneMark  as part of the genome annotation pipeline in the Integrated Microbial Genomes Expert Review system . The tRNAs were identified by tRNAScan-SE-1.23 , while ribosomal RNA genes within the genome were predicted using the tool RNAmmer . Other non-coding RNA genes were predicted using Infernal . CRISPR elements were identified by the program CRT . Manual functional annotation was performed within the IMG platform  and the Artemis Genome Browser .
The genome of A. fulgidus strain 7324 comprises one circular chromosome with a total size of 2,316,287 bp, which is 137,887 bp larger than A. fulgidus VC16T DSM 3404 . The mole percent G + C is 48.08, which is slightly higher than the 47% value estimated previously with thermal denaturation  and slightly lower than for the type strain DSM 3404 (48.6%); in any case within the 1% threshold with the species’ type strain VC16T sensu Meier-Kolthoff et al. . No plasmids were detected. The strain 7324 genome is the largest of the genome-sequenced Archaeoglobus species, the smallest one being the A. profundus genome with a total size of 1.56 Mbp . Out of the total 2615 genes annotated in the 7324 genome, 2558 were identified as protein coding genes and 56 as RNA genes (Table 3). Only 67.29% of the genes could be assigned to COG functional categories as listed in Table 4. Five CRISPR repeat regions were identified, as compared with only three in strain VC16T (AE000782). There is only one rRNA operon (Fig. 3). As for VC16T, there is no apparent GC skew in the genome, which could indicate the presence of multiple DNA replication origins and explain previous difficulties in precise mapping of replication origin(s) in this species using a marker rescue analysis approach [29, 30].
Insights from the genome sequence
Genes encoding central metabolic pathways like dissimilatory sulfate reduction, lactate oxidation, a complete TCA cycle and the Acetyl-CoA pathway were identified in strain 7324, confirming previous physiological characterization and similarity with strain VC16T . A genome alignment between strains VC16T and 7324 revealed a large degree of genomic similarity and sequence synteny, interrupted mostly by two large additional regions of about 64 and 109 Kbp (InsI and InsII, respectively) in strain 7324 as compared to VC16T (Fig. 4). Both these regions are flanked by a disrupted tRNA gene, which implies that InsI and II represent genomic insertions. They also possess a considerably lower G + C content (42%) as compared to the average of the genome (Fig. 3), indicating a result of recombination with an AT-rich organism. This is further supported by identification of a non-coding high-GC RNA gene in InsII belonging to the hgcC family of ncRNA (RFAM v12 accession code RF00062) typically found in AT-rich hyperthermophiles (Figs. 3 and 5b). This ncRNA family was originally identified in the genomes of Pyrococcus furiosus and Methanocaldococcus jannaschii  but its function is still unresolved.
InsII also contains a gene encoding an extra homologue of the Orc1/Cdc6 family of replication initiation control proteins in addition to the two other cdc6 homologous in the 7324 and VC16T genomes. The closest homologue identified by a blast search is from A. veneficus (62% amino acid sequence identity). The majority of the other genes are hypothetical or have a general function prediction only. InsI carries two CRISPR repeat regions and 14 genes encoding CRISPR-associated proteins including a Cas6 homologue (Fig. 5a). The rest of this insert mostly contains hypothetical genes.
A Venn diagram shows that A. fulgidus strains VC16T and 7324 share a large number of genes (2292) (Fig. 6a), reflecting the high degree of genome similarity. The 263 genes unique to strain 7324 include about 200 hypothetical genes/uncharacterized functions most of them belonging to the large insertions and the CRISPR-associated genes of InsII. The Venn diagram including all the five genome-sequenced strains revealed an Archaeoglobus core genome of 1001 genes (Fig. 6b), most of which encode energy-yielding, biosynthetic and regulatory functions. About 200 of the core genes belong to the hypothetical/uncharacterized category in the EggNog database . This is considerably lower than the 32% fraction of unassigned genes for the entire strain 7324, but underpins that a large part of central gene functions in this genus still remain to be disclosed. About 2900 genes belong to the Archaeoglobus pan-genome, being unique to one of the genomes or shared by 2 to 4 of the species.
Archaeoglobus fulgidus strain 7324 has been reported to degrade starch  and several enzymes involved in starch degradation have been purified from cells grown on starch. These enzymes include cyclodextrin gluconotransferase, cyclodextrinase, maltodextrin phosphorylase, and phosphoglucomutase, ADP-dependent glucokinase, ADP-dependent phosphofructokinase and pyruvate kinase [34,35,36,37]. The enzymes were biochemically characterized and the N-termini (with exception of cyclodextrinase) have been determined. They showed highest sequence identity to proteins from Thermococcales , e.g. from Thermococcus litoralis . However, in the genome of A. fulgidus presented here, genes encoding sugar degrading enzymes have not been identified. Rather, various analyses indicated that A. fulgidus strain 7324, deposited as DSM 8774, was contaminated with a Thermococcus strain. The genome of this Thermococcus strain has also been sequenced and all genes encoding the sugar degrading enzymes, originally attributed to A. fulgidus strain 7324, were found in the Thermococcus strain genome. Also, glutamate dehydrogenase from a DSM 8774 culture grown in the presence of small amounts of yeast extract (0.3 gL−) was previously purified and characterized [38, 39], but appears also to be encoded by the Thermococcus strain genome. The purity of the original A. fulgidus 7324 isolate was not assessed by 16S rRNA gene sequencing prior to deposition at DSMZ  and whether the Thermococcus contamination was present in the original culture or has been introduced at a different stage is not known. The genome analysis of this Thermococcus strain, which appears to represent a novel Thermococcus species, will be published separately.
The complete genome of A. fulgidus strain 7324, recovered from hot water produced from an oil well in the North Sea was sequenced and annotated. In addition to the A. fulgidus type strain, VC16, isolated from a shallow hot vent in the Mediterranean, this is the second A. fulgidus genome to be characterized. The two strains share 93.5% genome sequence similarity, and differ mostly by two large insertions of 64 and 109 Kbp in strain 7324 that seem to have originated from an AT-rich archaeon. The insertions carry two additional CRISPR elements, an extra cdc6 gene, a variety of mobile genetic elements and a large number of hypothetical and unassigned genes. Based on comparison with four other Archaeoglobus spp. genomes, the Archaeoglobus core genome was estimated to 1001 genes. No particular traits indicating adaptation to the petroleum reservoir subsurface environment could be identified.
Beeder J, Nilsen RK, Rosnes JT, Torsvik T, Lien T. Archaeoglobus fulgidus isolated from hot North Sea oil field waters. Appl Env Microbiol. 1994;60:1227–31.
Stetter KO. Archaeoglobus fulgidus gen. nov., sp. nov, a new taxon of extremely thermophilic Archaebacteria. Syst Appl Microbiol. 1988;10:172–3.
Klenk HP, Clayton RA, Tomb JF, White O, Nelson KE, Ketchum KA, Dodson RJ, Gwinn M, Hickey EK, Peterson JD, et al. The complete genome sequence of the hyperthermophilic, sulphate-reducing archaeon Archaeoglobus fulgidus. Nature. 1997;390:364–70.
Parey K, Fritz G, Ermler U, Kroneck PMH. Conserving energy with sulfate around 100 degrees C - structure and mechanism of key metal enzymes in hyperthermophilic Archaeoglobus fulgidus. Metallomics. 2013;5:302–17.
Burggraf S, Jannasch HW, Nicolaus B, Stetter KO. Archaeoglobus profundus sp. nov., represents a new species within the sulfate-reducing Archaebacteria. Syst Appl Microbiol. 1990;13:24–8.
Huber H, Jannasch H, Rachel R, Fuchs T, Stetter KO. Archaeoglobus veneficus sp. nov., a novel facultative chemolithoautotrophic hyperthermophilic sulfite reducer, isolated from abyssal black smokers. Syst Appl Microbiol. 1997;20:374–80.
Mori K, Maruyama A, Urabe T, Suzuki KI, Hanada S. Archaeoglobus infectus sp. nov., a novel thermophilic, chemolithoheterotrophic archaeon isolated from a deep-sea rock collected at Suiyo seamount, Izu-Bonin arc, western Pacific Ocean. Int J Syst Evol Micr. 2008;58:810–6.
Steinsbu BO, Thorseth IH, Nakagawa S, Inagaki F, Lever MA, Engelen B, Øvreas L, Pedersen RB. Archaeoglobus sulfaticallidus sp. nov., a thermophilic and facultatively lithoautotrophic sulfate-reducer isolated from black rust exposed to hot ridge flank crustal fluids. Int J Syst Evol Micr. 2010;60:2745–52.
Stetter KO, Huber R, Blochl E, Kurr M, Eden RD, Fielder M, Cash H, Vance I. Hyperthermophilic Archaea are thriving in deep North-Sea and Alaskan oil-reservoirs. Nature. 1993;365:743–5.
Slobodkina GB, Kolganova TV, Querellou J, Bonch-Osmolovskaya EA, Slobodkin AI. Geoglobus acetivorans sp. nov., an iron(III)-reducing archaeon from a deep-sea hydrothermal vent. Int J Syst Evol Micr. 2009;59:2880–3.
Kashefi K, Tor JM, Holmes DE, Van Praagh CVG, Reysenbach AL, Lovley DR. Geoglobus ahangari gen. nov., sp. nov., a novel hyperthermophilic archaeon capable of oxidizing organic acids and growing autotrophically on hydrogen with Fe(III) serving as the sole electron acceptor. Int J Syst Evol Micr. 2002;52:719–28.
Tindall BJ, Euzeby JP. Status of strains that contravene rules 27 (3) and 30 of the international code of nomenclature of bacteria. Opinion 81. Int J Syst Evol Microbiol. 2008;58:1755–63.
Hafenbradl D, Keller M, Dirmeier R, Rachel R, Rossnagel P, Burggraf S, Huber H, Stetter KO. Ferroglobus placidus gen. nov., sp. nov., a novel hyperthermophilic archaeum that oxidizes Fe2+ at neutral pH under anoxic conditions. Arch Microbiol. 1996;166:308–14.
Zellner G, Stackebrandt E, Kneifel H, Messner P, Sleytr UB, Demacario EC, Zabel HP, Stetter KO, Winter J. Isolation and characterization of a thermophilic, sulfate reducing Archaebacterium, Archaeoglobus fulgidus strain-Z. Syst Appl Microbiol. 1989;11:151–60.
Meier-Kolthoff JP, Auch AF, Klenk HP, Göker M. Genome sequence-based species delimitation with confidence intervals and improved distance functions. BMC Bioinf. 2013; doi:10.1186/1471-2105-14-60.
Lien T, Madsen M, Rainey FA, Birkeland NK. Petrotoga mobilis sp. nov., from a North Sea oil-production well. Int J Syst Bacteriol. 1998;48:1007–13.
The DOE Joint Genome Institute. http://www.jgi.doe.gov. Accessed 15 Feb 2017.
Kostyukova AS, Gongadze GM, Polosina YY, Bonch-Osmolovskaya EA, Miroshnichenko ML, Chernyh NA, Obraztsova MV, Svetlichny VA, Messner P, Sleytr UB, et al. Investigation of structure and antigenic capacities of Thermococcales cell envelopes and reclassification of “Caldococcus litoralis” Z-1301 as Thermococcus litoralis Z-1301. Extremophiles. 1999;3:239–45.
Phrap, Phred and Consed for Windows, MacOS, Linux, and Unix. www.phrap.com. Accessed 15 Feb 2017.
Besemer J, Borodovsky M. GeneMark: web software for gene finding in prokaryotes, eukaryotes and viruses. Nucleic Acids Res. 2005;33:451–4.
Markowitz VM, Mavromatis K, Ivanova NN, Chen IM, Chu K, Kyrpides NC. IMG ER: a system for microbial genome annotation expert review and curation. Bioinformatics. 2009;25:2271–8.
Hyatt D, Chen GL, Locascio PF, Land ML, Larimer FW, Hauser LJ. Prodigal: prokaryotic gene recognition and translation initiation site identification. BMC Bioinf. 2010; doi:10.1186/1471-2105-11-119.
Lagesen K, Hallin P, Rodland EA, Staerfeldt HH, Rognes T, Ussery DW. RNAmmer: consistent and rapid annotation of ribosomal RNA genes. Nucleic Acids Res. 2007;35:3100–8.
Nawrocki EP, Kolbe DL, Eddy SR. Infernal 1.0: inference of RNA alignments. Bioinformatics. 2009;25:1335–7.
Bland C, Ramsey TL, Sabree F, Lowe M, Brown K, Kyrpides NC, Hugenholtz P. CRISPR recognition tool (CRT): a tool for automatic detection of clustered regularly interspaced palindromic repeats. BMC Bioinf. 2007; doi:10.1186/1471-2105-8-209.
Rutherford K, Parkhill J, Crook J, Horsnell T, Rice P, Rajandream MA, Barrell B. Artemis: sequence visualization and annotation. Bioinformatics. 2000;16:944–5.
Meier-Kolthoff JP, Klenk HP, Göker M. Taxonomic use of DNA G+C content and DNA-DNA hybridization in the genomic age. Int J Syst Evol Microbiol. 2014;64:352–6.
von Jan M, Lapidus A, Del Rio TG, Copeland A, Tice H, Cheng JF, Lucas S, Chen F, Nolan M, Goodwin L, et al. Complete genome sequence of Archaeoglobus profundus type strain (AV18(T)). Stand Genomic Sci. 2010;2:327–46.
Maisnier-Patin S, Malandrin L, Birkeland NK, Bernander R. Chromosome replication patterns in the hyperthermophilic euryarchaea Archaeoglobus fulgidus and Methanocaldococcus (Methanococcus) jannaschiii (vol 45, pg 1443, 2001). Mol Microbiol. 2002;46:599.
Maisnier-Patin S, Malandrin L, Birkeland NK, Bernander R. Chromosome replication patternsin the hyperthermophilic euryarchaea Archaeoglobus fulgidus and Methanocaldococcus (Methanococcus) jannaschii. Mol Microbiol. 2002;45:1443–50.
Klein RJ, Misulovin Z, Eddy SR. Noncoding RNA genes identified in AT-rich hyperthermophiles. P Natl Acad Sci USA. 2002;99:7542–7.
Huerta-Cepas J, Szklarczyk D, Forslund K, Cook H, Heller D, Walter MC, Rattei T, Mende DR, Sunagawa S, Kuhn M, et al. eggNOG 4.5: a hierarchical orthology framework with improved functional annotations for eukaryotic, prokaryotic and viral sequences. Nucleic Acids Res. 2016;44:286–93.
Labes A, Schönheit P. Sugar utilization in the hyperthermophilic, sulfate-reducing archaeon Archaeoglobus fulgidus strain 7324: starch degradation to acetate and CO2 via a modified Embden-Meyerhof pathway and acetyl-CoA synthetase (ADP-forming). Arch Microbiol. 2001;176:329–38.
Hansen T, Schönheit P. ADP-dependent 6-phosphofructokinase, an extremely thermophilic, non-allosteric enzyme from the hyperthermophilic, sulfate-reducing archaeon Archaeoglobus fulgidus strain 7324. Extremophiles. 2004;8:29–35.
Johnsen U, Hansen T, Schönheit P. Comparative analysis of pyruvate kinases from the hyperthermophilic archaea Archaeoglobus fulgidus, Aeropyrum pernix, and Pyrobaculum aerophilum and the hyperthermophilic bacterium Thermotoga maritima - unusual regulatory properties in hyperthermophilic archaea. J Biol Chem. 2003;278:25417–27.
Labes A, Schönheit P. ADP-dependent glucokinase from the hyperthermophilic sulfate-reducing archaeon Archaeoglobus fulgidus strain 7324. Arch Microbiol. 2003;180:69–75.
Labes A, Schönhelt P. Unusual starch degradation pathway via cyclodextrins in the hyperthermophilic sulfate-reducing archaeon Archaeoglobus fulgidus strain 7324. J Bacteriol. 2007;189:8901–13.
Aalen N, Steen IH, Birkeland NK, Lien T. Purification and properties of an extremely thermostable NADP(+)-specific glutamate dehydrogenase from Archaeoglobus fulgidus. Arch Microbiol. 1997;168:536–9.
Steen IH, Hvoslef H, Lien T, Birkeland NK. Isocitrate dehydrogenase, malate dehydrogenase, and glutamate dehydrogenase from Archaeoglobus fulgidus. Method Enzymol. 2001;331:13–26.
Woese CR, Kandler O, Wheelis ML. Towards a natural system of organisms: proposal for the domains Archaea, bacteria, and Eucarya. Proc Natl Acad Sci U S A. 1990;87:4576–9.
Garrity GM, Holt JG. Phylum AII. Euryarchaeota phy. nov. In: Garrity GM, Boone DR, Castenholz RW, editors. Bergey’s manual of systematic bacteriology, second edition, volume 1. New York: Springer; 2001. p. 211–355.
Garrity GM, Holt JG. Class VI. Archaeoglobi class. nov. In: Garrity GM, Boone DR, Castenholz RW, editors. Bergey’s manual of systematic bacteriology, second edition, volume 1. The Archaea and the deeply branching and phototrophic Bacteria. New York: Springer; 2001. p. 349.
Huber H, Stetter KO. Order I. Archaeoglobales ord. nov. In: Garrity GM, Boone DR, Castenholz RW, editors. Bergey’s manual of systematic bacteriology, second edition, volume 1. The Archaea and the deeply branching and phototrophic Bacteria. New York: Springer; 2001. p. 349.
Huber H, Stetter KO. Family I. Archaeoglobaceae fam. nov. In: Garrity GM, Boone DR, Castenholz RW, editors. Bergey’s manual of systematic bacteriology, second edition, volume 1. The Archaea and the deeply branching and phototrophic Bacteria. New York: Springer; 2001. p. 349.
Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT, et al. Gene ontology: tool for the unification of biology. Nat Genet. 2000;25:25–9.
Larkin MA, Blackshields G, Brown NP, Chenna R, McGettigan PA, McWilliam H, Valentin F, Wallace IM, Wilm A, Lopez R, et al. Clustal W and clustal X version 2.0. Bioinformatics. 2007;23:2947–8.
Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. MEGA6: molecular evolutionary genetics analysis version 6.0. Mol Biol Evol. 2013;30:2725–9.
Darling ACE, Mau B, Blattner FR, Perna NT. Mauve: multiple alignment of conserved genomic sequence with rearrangements. Genome Res. 2004;14:1394–403.
Bardou P, Mariette J, Escudie F, Djemiel C, Klopp C. jvenn: an interactive Venn diagram viewer. BMV. Bioinformatics. 2014;15 doi:10.1186/1471-2105-15-293.
EzBioCloud’s Comparative Genomics Database. https://www.ezbiocloud.net/. Accessed 15 Mar 2017.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.