- Extended genome report
- Open Access
First draft genome sequence of a strain from the genus Fusibacter isolated from Salar de Ascotán in Northern Chile
Standards in Genomic Sciencesvolume 12, Article number: 43 (2017)
Fusibacter sp. 3D3 (ATCC BAA-2418) is an arsenate-reducing halotolerant strain within the Firmicutes phylum, isolated from the Salar de Ascotán, a hypersaline salt flat in Northern Chile. This high-Andean closed basin is an athalassohaline environment located at the bottom of a tectonic basin surrounded by mountain range, including some active volcanoes. This landscape can be an advantageous system to explore the effect of salinity on microorganisms that mediate biogeochemical reactions. Since 2000, microbial reduction of arsenic has been evidenced in the system, and the phylogenetic analysis of the original community plus the culture enrichments has revealed the predominance of Firmicutes phylum. Here, we describe the first whole draft genome sequence of an arsenic-reducing strain belonging to the Fusibacter genus showing the highest 16S rRNA gene sequence similarity (98%) with Fusibacter sp. strain Vns02. The draft genome consists of 57 contigs with 5,111,250 bp and an average G + C content of 37.6%. Out of 4780 total genes predicted, 4700 genes code for proteins and 80 genes for RNAs. Insights from the genome sequence and some microbiological features of the strain 3D3 are available under Bioproject accession PRJDB4973 and Biosample SAMD00055724. The release of the genome sequence of this strain could contribute to the understanding of the arsenic biogeochemistry in extreme environments.
Salt flats or salares are one of the most interesting biomes on earth . Based on their hypersaline conditions, prokaryotes have evolved to develop biochemical processes with potential applications in biotechnology , providing also opportunities for biosignature detection on other planets . Due to geological, climatic and geomorphological factors, dozens of endorheic basins are located in Northern Chile including evaporitic bodies and saline lakes. Brines and crusts of these saline deposits are enriched in arsenic .
The Salar de Ascotán  is an athalassohaline environment located at the bottom of a tectonic basin surrounded by volcanic systems in east-west direction, including some active volcanoes with altitudes from 5000 to 6000 m.a.s.l. . The saline crusts are mainly composed of chlorides (halite) and sulfates (gypsum) to economic boron-bearing minerals associated with significant amounts of arsenic sulfides , with the arsenic concentrations the highest found in the area .
In order to understand the bacterial role in the arsenic biogeochemical cycle at circumneutral pH, several sampling expeditions to Salar de Ascotán, in the Chilean High-Andes, have been taken since 2000. The microbial diversity of this salt flat was first analyzed , then enrichment , isolation and sequencing efforts  as well as geochemical in situ investigations wer performed . In addition, the distribution of genes for the As (III) oxidation (aioA), As (V) detoxifying respiration (arrA), As detoxification (arsC), and As (III) extrusion (acr3) was explored in Salar de Ascotán and other natural environments in Northern Chile with arsenic concentrations spanning six orders of magnitude. The abundance of Firmicutes -like arsC genes compared to the Enterobacterial-like arsC genes in these environments suggested an important role of thioredoxin and the Firmicutes phylum in the local As biogeochemistry .
Fusibacter is a minor genus into the Clostridiales order within the Firmicutes phylum. Currently, it comprises four Gram-positive species with validly published names. This group started with the discovery of the thiosulfate-reducing bacterium Fusibacter paucivorans , being the most studied, isolated from oil-producing wells . Fusibacter tunisiensis was isolated from an anaerobic reactor used to treat olive-mill wastewater . Recently, Fusibacter bizertensis was identified from a corroded kerosene storage tank , and more recently, Fusibacter fontis was the first species of this genus isolated from a natural environment . In general terms, the reported members of this genus are fermentative and halotolerant anaerobes. Moreover, these species share sulfur-reducing features capable of generating sulfide starting from elemental sulfur [13, 15] or thiosulfate [12, 13] sources. To date, a whole-genome sequence has not been reported for any species within this genus.
Here, we report the first draft genome of a strain of Fusibacter plus some microbiological properties of this halotolerant isolate, recovered from a saline environment in Northern Chile. The strain was deposited as Fusibacter sp. 3D3 as ATCC BAA-2418 because we are still running the necessary tests and deposits to describe the isolate as a new species and “ Fusibacter ascotence” is the proposed species name.
This report contributes to a better understanding of the ecophysiology of extreme halotolerant microorganisms inhabiting saline environments and their role in the arsenic biogeochemistry.
Classification and features
Fusibacter sp. 3D3 is an indigenous strain of the Salar de Ascotán hypersaline sediments isolated at the Centro de Biotecnología, Universidad Católica del Norte, Antofagasta, Chile. Enrichment, isolation, and growth experiments were performed in a fresh Newman-modified minimal medium  containing, 1% (w/v) NaCl, 0.1% (w/v) yeast extract, and 1 mM cysteine adjusted to pH 7.0. After autoclaving, 10 mM lactate as electron donor, and 20 mM sodium sulfate and 2 mM sodium arsenate as electron acceptors were added in order to complete 20 mL of medium into 50 mL-anaerobic-bottles (Supelco). The strain was incubated in an anaerobic chamber (Airlock, Coydrive), in dark, at 30 °C, under N2:CO2:H2 gas atmosphere (80:15:5, v/v) up to 10 days. The pure colonies were obtained by inclined tubes of agar prepared with the mentioned Newman modified medium plus 2% (w/v) agar incubated under anaerobic conditions at 30 °C. Single yellow colonies were restreaked several times to obtain pure isolates and then were transferred to the liquid medium. Transmission electron microscopy revealed rod-shaped cells (0.4 μm × 3-10 μm) (Fig. 1a).
Arsenate reduction was tested by inoculation of 1 × 10−6 cells mL−1 into 20 mL of fresh Newman-modified medium under incubation conditions described above. An abiotic control was carried out in sterile medium without inoculum. Growth curves were performed in triplicate and monitored by counting chambers (0.01 mm × 0.0025 mm2, Neubauer, Marienfeld). Samples were acquired periodically, then centrifuged (15,000×g; 10 min) to remove cells, and finally filtered through a 0.2 μm cellulose filter. The filtered supernatant was sealed and refrigerated at 4 °C to preserve arsenic speciation until analysis. As (V) and As (III) concentrations were measured with a mobile phase of 10 mM acid phosphate at 6.25 pH by Millennium Excalibur HPLC System (PS Analytical, Orpington, UK). To quantify lactate and acetate, each filtrate was injected in a Dionex IonPac AS11-HC column to run a high-performance liquid chromatography (Thermo Scientific model 3200) with an isocratic concentration of KOH. Arsenate reduction and simultaneous lactate consumption were evidenced (Fig. 1b). Arsenate reduction has not been reported in the other members of the Fusibacter genus.
The Initial identification of strain 3D3 was performed by 16S rRNA gene amplification using a previously described method . The 16S rRNA sequences of strain 3D3 clustered with type strains of Fusibacter species when those were aligned using Clustal W and manually corrected. A phylogenetic tree was constructed using neighbor-joining, maximum-parsimony, and maximum likelihood algorithms with bootstrap values of 500 replicates using the MEGA program version 6.22. Phylogenetic analysis of the 16S rRNA sequence indicated that the strain 3D3 belongs to the genus Fusibacter and exhibits a similarity of 98% with Fusibacter sp. Vns02, and 95% with both Fusibacter paucivorans and Fusibacter tunisiensis (Fig. 2).
RapID™ NF Plus and RapID™ One (Thermo Scientific), two qualitative micromethods employing conventional and chromogenic substrates for the biochemical features identification were performed following the manufacturer’s instructions (Table 1). A single colony was inoculated into the given fluid and incubated for 24 h at 30 °C. Then, the inoculation fluid was transferred to the corresponding panel. The interpreted results were entered into the electronic RapID code database (ERIC electronic compendium, version 1.0.771, patch 0613). Comparing with the ERIC™ database, our results failed to identify our isolate (Table 1).
Genome sequencing information
Genome project history
Based on its phylogenetic position and 16S rRNA similarity, Fusibacter sp. 3D3 (Taxonomy ID: 1,048,380) was previously submitted to NCBI in 2010 (Gene Bank 16S rRNA gene: FR873490.1) (Fig 2). Later, in 2013, it was deposited as Fusibacter sp. 3D3 in ATCC BAA-2418, being the first strain of this genus coming from an extreme arsenic bearing and saline biotope (Table 2). This organism was selected for genome sequencing based on its interesting phenotypic characteristics. Recently, in 2016, the submission of the whole shotgun project assembled as a draft genome was performed to the DNA Data Bank of Japan under the Bioproject accession number PRJDB4973 and Biosample number SAMD00055724 (ID 573014). This Whole Genome Shotgun project has been deposited at GenBank under the accession BDHH00000000. The version described in this paper is the first version, BDHH01000000 . Table 3 presents the project information and its association with MIGS version 2.0 compliance .
Growth conditions and genomic DNA preparation
Fusibacter sp. strain 3D3 was grown anaerobically on fresh Newman  modified medium and conditions of incubation described above. DNA was extracted using High Pure Template Preparation Kit (Roche, Germany), according to the manufacturer instructions. Both quantity and quality of the genomic DNA were measured using a NanoDrop ND-1000 spectrophotometer (Thermo-Fisher Scientific Inc.) and analyzed by DGGE (200 C. B. S. Scientific Company), respectively. The purity of strain 3D3 was confirmed by a single band in the DGGE profile.
Genome sequencing and assembly
The genome of Fusibacter sp. 3D3 was sequenced on an Illumina MiSeq platform at Molecular Research Laboratory (MR. DNA, Shallowater, TX). The library for each sample was prepared using a Nextera DNA Sample Preparation Kit (Illumina), following the manufacturer’s instructions. Sequencing of 2 × 300-bp paired-end reads allowed for an estimate of 20,000 output with an average coverage over 50 times fold. The assemblage of quality-filtered reads was executed by MR-DNA for the complete genome sequence. As a result, the draft genome of ~5.1 Mbp size was generated. Reads were assembled de novo using Newbler v2.0.01.14. The final draft assembly contained 57 contigs identifying 4780 genes using RAST .
Genes were predicted using Glimmer 3.02  as part of the RAST annotation pipeline using SEED platform for 4780 features identified. Whole RNA genes were also predicted by the same annotation platform . The predicted protein coding genes were analyzed for the presence of signal peptides using SignalP 4.1 Server . TMHMM Server v. 2.0 was utilized for prediction of transmembrane helices in proteins . Geneious 7.1.9 (Biomatters) software was used to analyze COGs utilizing the BLAST COG database (Table 4). Pfam domains were computed using InterProScan 5.19-58.0 . CRISPRs repeats were found submitting the contigs to the CRISPRs Finder web server .
The draft genome for Fusibacter sp. 3D3 contained 5,111,250 nucleotides with an average G + C content of 37.6% (Table 5). From 4780 genes, 4700 were predicted protein coding, and 80 RNA coding genes (12 rRNA, and 62 tRNA genes). The putative function was assigned to 63.1% of the genes, while the remaining genes were annotated as hypothetical proteins. The distribution of genes in COGs functional categories is presented in Table 4.
Insights from the genome sequence
Similarity analysis of genes involved in the arsenic metabolism indicated that the closest available genome of strain 3D3 in the database was Clostridium sticklandii , which belongs to the Clostridiales family as well. The subsystem information approach to genome annotation performed by RAST/SEED  confirmed the relation to other members in the Clostridiales order (Table 6). Arsenic detoxification genes are clearly present in Fusibacter sp. 3D3 genome, however, genes coding for arsenate respiratory reductases (arr) and arsenite oxidases (aio) have a very low percentage of similarity with genes coding for the enzymes evidenced at protein level. The arsC gene sequence identified in the Fusibacter sp. genome was clustered inside the Firmicutes -like arsC gene clade whose predominance has been reported in Salar de Ascotán .
In the vicinity of the ArsC coding gene is the gen F3D3_RS05420. This piqued our interest and, in a deeper analysis, we found that F3D3_RS05420 codify for pyridine nucleotide-disulfide oxidoreductase NADH dehydrogenase (accession number: WP_069871897). The preliminary information indicates that the protein encoded by the gen F3D3_RS05420 is part of a new family of proteins of unknown function. However, the genomic context shows us some clues to formulate a hypothesis. By means of comparative genomics we identify two common components accompanying genes like F3D3_RS05420: A) genes codifying for transcription regulators and, B) genes codifying for arsenical transporters (Fig. 3c). This could be an indicative of a possible role in the response to stress by As. The multiple sequence alignment carried out using MUSCLE application  in CLC Genome Workbench 8.0 (Qiagen) shows that the protein is distributed in the Firmicutes bacteria and it is strongly conserved (Additional file 1: Figure S1).
The protein architecture of WP_069871897 shows a CoA-disulfide reductase domain (TIGR03385) and a rhodanese domain (PFAM00581). A rhodanese domain is also present in the ACR2 protein of Saccharomyces cerevisiae which also has arsenate reductase activity . The catalytic loop of the rhodanese domain has two known configurations, a short version with four residues to accommodate sulfur or selenium atoms and an extended version with five residues to accommodate an arsenic or phosphorus atoms . The architecture of WP_069871897 and related proteins suggest a role similar to ACR2, namely arsenic reductase. However, the comparison between the catalytic loops of ACR2 (Q06597) “CTGSKNRG” with the “CNKGVTGN” of WP_069871897 does not show an apparent similitude, which makes it difficult to extrapolate the activity of ACR2 with WP_069871897. In addition, the presence of the arsC gene in the compared genomes (Fig. 3c), but not in the same genomic context, suggests that the WP_069871897 and related proteins does not substitute the ArsC activity. Then, it remains a challenge to the scientific community to answer if the proteins similar to WP_069871897 are a new kind of arsenic reductase or if they are in some way involved with the response to arsenic stress.
The 5.11 Mbp draft genome sequence of Fusibacter sp. 3D3 is arranged in 57 contigs, being the first Fusibacter draft genome published. It potentially includes 4700 protein-coding genes, 67.1% of which were assigned to function prediction. 80 RNA genes partitioned in 12 rRNA and 62 tRNAs were identified. The release of the genome sequence of this strain will provide new insights into arsenic reduction processes in hypersaline biomes and further understanding of the mechanisms used by halophile bacteria to endure high osmotic stress in natural and industrial saline environments.
Meters above sea level
The minimum information about a genome sequence
Rapid annotations using subsystems technology
Demergasso C, Dorador C, Meneses D, Blamey J, Cabrol N, Escudero L, et al. Prokaryotic diversity pattern in high-altitude ecosystems of the Chilean Altiplano. J Geophys Res Biogeosci. 2010;115:G00D09.
Bull AT. Asenjo J a. Microbiology of hyper-arid environments: recent insights from the Atacama Desert, Chile. Antonie Van Leeuwenhoek. 2013;103:1173–9.
Parro V, de Diego-Castilla G, Moreno-Paz M, Blanco Y, Cruz-Gil P, Rodríguez-Manfredi JA, et al. A microbial oasis in the hypersaline Atacama subsurface discovered by a life detector chip: implications for the search for life on Mars. Astrobiology. 2011;11:969–96.
Risacher F, Fritz B. Origin of salts and brine evolution of Bolivian and Chilean salars. Aquat Geochem. 2009;15:123–57.
Chong G, Pueyo JJ, Demergasso C. The borate deposits in Chile. Rev Geol Chile. 2000;27:99–119.
Chong G. Die Salare in Nordchile-Geologie. Struktur und geochimie Goetekt Forsch. 1984;67:1984.
Demergasso CS, Guillermo CD, Lorena EG, Mur JJP, Pedrós-Alió C. Microbial precipitation of arsenic sulfides in Andean salt flats. Geomicrobiol J. 2007;24:111–23.
Lara J, Escudero González L, Ferrero M, Chong Díaz G, Pedrós-Alió C, Demergasso C. Enrichment of arsenic transforming and resistant heterotrophic bacteria from sediments of two salt lakes in Northern Chile. Extremophiles. 2012;16:523–38.
Valdés N, Rivera-Araya J, Bijman J, Escudero L, Demergasso C, Fernández S. Resistant Gammaproteobacterium isolated from a salt flat. Genome Announc. 2014;2:2006–7.
Nuñez C. Estudio de los procesos Biogeoquímicos que controlan la distribución del arsénico en el Salar de Ascotán y Gorbea. Antofagasta: Master thesis, Universidad Catolica del Norte; 2011.
Escudero LV, Casamayor EO, Chong G, Pedrós-Alió C, Demergasso C. Distribution of microbial arsenic reduction, oxidation and extrusion genes along a wide range of environmental arsenic concentrations. PLoS One. 2013;8:e78890.
Ravot G, Magot M, Fardeau ML, Pateli BKC, Garcia JL, Ollivier B, et al. Fusibacter paucivorans gen. nov., sp. nov., an anaerobic, thiosulfate-reducing bacterium from an oil-producing well. Int J Syst Bacteriol. 1999;49 Pt 3:1141–7.
Ben Hania W, Fraj B, Postec A, Fadhlaoui K, Hamdi M, Ollivier B, et al. Fusibacter tunisiensis sp. nov., isolated from an anaerobic reactor used to treat olive-mill wastewater. Int J Syst Evol Microbiol. 2012;62 Pt 6:1365–8.
Smii L, Ben Hania W, Cayol J-L, Joseph M, Hamdi M, Ollivier B, et al. Fusibacter bizertensis sp. nov., isolated from a corroded kerosene storage tank. Int J Syst Evol Microbiol. 2015;65 Pt 1:117–21.
Fadhlaoui K, Ben Hania W, Postec A, Fauque G, Hamdi M, Ollivier B, et al. Fusibacter fontis sp. nov., a sulfur-reducing, anaerobic bacterium isolated from a mesothermic Tunisian spring. Int J Syst Evol Microbiol. 2015;65:3501–6.
Mashima J, Kodama Y, Kosuge T, Fujisawa T, Katayama T, Nagasaki H, et al. DNA data bank of Japan (DDBJ) progress report. Nucleic Acids Res. 2016;44:D51–7.
Field D, Garrity G, Gray T, Morrison N, Selengut J, Sterk P, et al. The minimum information about a genome sequence (MIGS) specification. Nat Biotechnol. 2008;26:541–7.
Newman DK, Beveridge TJ, Morel F. Precipitation of arsenic Trisulfide by Desulfotomaculum auripigmentum. Appl Environ Microbiol. 1997;63:2022–8.
Aziz RK, Bartels D, Best AA, DeJongh M, Disz T, Edwards RA, et al. The RAST server: rapid annotations using subsystems technology. BMC Genomics. 2008;9:75.
Delcher AL, Bratke KA, Powers EC, Salzberg SL. Identifying bacterial genes and endosymbiont DNA with glimmer. Bioinformatics. 2007;23:673–9.
Overbeek R, Olson R, Pusch GD, Olsen GJ, Davis JJ, Disz T, et al. The SEED and the rapid annotation of microbial genomes using subsystems technology (RAST). Nucleic Acids Res. 2014;42:D206–14.
Petersen TN, Brunak S, von Heijne G, Nielsen H. SignalP 4.0: discriminating signal peptides from transmembrane regions. Nat Methods. 2011;8:785–6.
Krogh A, Larsson B, von Heijne G, Sonnhammer ELL. Predicting transmembrane protein topology with a hidden Markov model: application to complete genomes. J Mol Biol. 2001;305:567–80.
Jones P, Binns D, Chang HY, Fraser M, Li W, McAnulla C, et al. InterProScan 5: genome-scale protein function classification. Bioinformatics. 2014;30:1236–40.
Grissa I, Vergnaud G, Pourcel C, Bland C, Ramsey TL, Sabree F, et al. CRISPRFinder: a web tool to identify clustered regularly interspaced short palindromic repeats. Nucleic Acids Res. 2007;35(Web Server issue):W52–7.
Fonknechten N, Chaussonnerie S, Tricot S, Lajus A, Andreesen JR, Perchat N, et al. Clostridium sticklandii, a specialist in amino acid degradation:revisiting its metabolism through its genome sequence. BMC Genomics. 2010;11:555.
Overbeek R, Begley T, Butler RM, Choudhuri JV, Chuang HY, Cohoon M, et al. The subsystems approach to genome annotation and its use in the project to annotate 1000 genomes. Nucleic Acids Res. 2005;33:5691–702.
Edgar RC. MUSCLE: a multiple sequence alignment method with reduced time and space complexity. BMC Bioinformatics. 2004;5:113.
Mukhopadhyay R, Rosen BP. Saccharomyces Cerevisiae ACR2 gene encodes an arsenate reductase. FEMS Microbiol Lett. 1998;168:127–36.
Bordo D, Bork P. The rhodanese/Cdc25 phosphatase superfamily. Sequence-structure-function relations. EMBO Rep. 2002;3:741–6.
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.
Ludwig W, Schleifer KH, Whitman WB. Revised road map to the phylum Firmicutes. In: De Vos P, Garrity G, Jones D, Krieg NR, Ludwig W, Rainey FA, Schleifer KH, Whitman WB, editors. Bergey's manual of systematic bacteriology, vol. 3. 2nd ed. New York: Springer-Verlag; 2009. p. 1–12.
Yutin N, Galperin MY. A genomic update on clostridial phylogeny: gram-negative spore formers and other misplaced clostridia. Environ Microbiol. 2013;15(10):2631–41.
Prévot AR.. Dictionnaire des Bactéries Pathogenes In: Hauduroy P, Ehringer G, Guillot G, Magrou J, Prevot AR, Rosset A, Urbain A, editors. 2nd ed. Masson: Paris; 1953.
Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, et al. Gene ontology: tool for the unification of biology. Nat Genet. 2000;25:25–9.
Thanks to Eng. Christian Castillo Herrera for all informatics support in Linux-based data analysis and MV Coalova for her contribution to improving our English.
This article was funded by FONDECYT Project 1,100,795 from the Chilean National Comission for Science and Technology (CONICYT), Fondo de Innovación para la Competitividad Regional de Antofagasta FIC-R. 2013/BIP-30136799-0 and FIC-R 2015/BIP 30412822-0 and the Research Support for Biodiversity from Minera Escondida Ltda. AES is a CONICYT Becas-Chile postdoctoral fellow N° 74,160,010.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
CLUSTAL multiple sequence alignment of proteins related to WP_069871897 of Fusibacter sp. strain 3D3. (ALN 76400 bytes). (DOCX 38 kb)