- Short genome report
- Open Access
Complete genome sequence of the nitrogen-fixing bacterium Azospirillum humicireducens type strain SgZ-5T
Standards in Genomic Sciences volume 13, Article number: 28 (2018)
The Azospirillum humicireducens strain SgZ-5T, belonging to the Order Rhodospirillales and the Family Rhodospirillaceae, was isolated from a microbial fuel cell inoculated with paddy soil. A previous work has shown that strain SgZ-5T was able to fix atmospheric nitrogen involved in plant growth promotion. Here we present the complete genome of A. humicireducens SgZ-5T, which consists of a circular chromosome and six plasmids with the total genome size of 6,834,379 bp and the average GC content of 67.55%. Genome annotations predicted 5969 protein coding and 85 RNA genes including 14 rRNA and 67 tRNA genes. By genomic analysis, we identified a complete set of genes that is potentially involved in nitrogen fixation and its regulation. This genome also harbors numerous genes that are likely responsible for phytohormones production. We anticipate that the A. humicireducens SgZ-5T genome will contribute insights into plant growth promoting properties of Azospirillum strains.
Bacteria that live in the plant rhizosphere and possess a large array of potential mechanisms to enhance plant growth are considered as PGPR [1,2,3]. Azospirillum represents a well characterized genus of PGPR due to its capacity of fixing atmospheric nitrogen [4, 5]. Although the exact contribution of Azospirillum to biological nitrogen fixation in plant growth promotion is debated , agricultural applications of the genus Azospirillum have been still developed [6, 7]. Another main characteristic of Azospirillum proposed to explain plant growth promotion has been related to its ability to produce phytohormones [8, 9].
At present, there are 17 species within the genus Azospirillum , of which the nitrogen-fixing bacterium A. humicireducens SgZ-5T, the focus species of this study, was initially isolated from the anode biofilm of a MFC. A soil sample collected from paddy field in Guangzhou City, Guangdong Province, China (23.18o N 113.36o E) was used as inoculating source of the MFC. In a previous report , the nitrogen-fixing capability of strain SgZ-5T was confirmed by acetylene-reduction assay and identification of a nifH gene. Furthermore, this strain has the ability to grow under anaerobic conditions via the oxidation of various organic compounds coupled to the reduction of humus , showing its potential use in plant rhizosphere. Here, we describe the physiological features together with the whole genome sequence of A. humicireducens SgZ-5T.
Classification and features
A. humicireducens SgZ-5T is a Gram-negative, facultative anaerobic, motile, spiral, straight to slightly curved rod-shaped bacterium (Fig. 1), belonging to the Order Rhodospirillales and the Family Rhodospirillaceae . The strain grew optimally in the conditions of 30 °C, pH 7.2, and 1% NaCl . On NA, strain SgZ-5T formed cream-colored, round, smooth, convex and non-translucent colonies (Fig. 1). With AQDS as the sole terminal electron acceptor, strain SgZ-5T could utilize pyruvate, glucose and acetate as electron donors under anaerobic conditions . Strain SgZ-5T was able to use a range of carbon substrates including N-Acetyl-glucosamine, citrate, D-ribose, meso-inositol, D-saccharose, D-maltose, L-rhamnose, suberic acid, malonate, acetate, L-serine, salicin, L-lactate, L-alanine, gluconate, 2-keto-gluconate, glycogen, D-mannitol, D-glucose, D-melibiose, L-fucose, D-sorbierite, L-arabinose, L-histidine, 3-hydroxy-butyric acid, 4-hydroxy-benzoic acid, L-proline, capric acid, adipic acid and malic acid  (Table 1).
A phylogenetic tree was constructed from aligning the 16S rRNA gene sequences of strain SgZ-5T and type strains of the genus Azospirillum by MEGA 5 using the neighbour-joining method . The phylogenetic position of strain SgZ-5T is shown in Fig. 2, where A. humicireducens can be grouped as a Azospirillum species, forms a distinct subclade together with A. lipoferum that are known as a biofertilizer widely used for agricultural production [13, 14]. The 16S rRNA gene of strain SgZ-5T is 98% similar to that of A. lipoferum NCIMB 11861T. Since nifH gene is highly conserved among nitrogen-fixing Proteobacteria , a nifH-based phylogenetic tree was constructed to identify the relationship of A. humicireducens to other species within the genus Azospirillum and related genus (Additional file 1). The phylogenetic reconstruction indicated the close relationship of the A. humicireducens SgZ-5T nifH gene with that from Azospirillum sp. B510.
Genome sequencing information
Genome project history
A. humicireducens SgZ-5T was selected for genome sequencing on the basis of its biotechnological potential in agricultural applications as a PGPR likely harboring multiple PGPP . The complete genome sequences have been deposited at Gen-Bank/EMBL/DDBJ under the accession numbers CP015285.1, CP028902-CP028907. Project information is available from Genome Online database number Gp0150267 at Joint Genome Institute. In Table 2, we summarize the project information and its association with Minimum Information about a Genome Sequence (MIGS) .
Growth conditions and genomic DNA preparation
A. humicireducens SgZ-5T was routinely cultured in NB medium containing (L− 1) 5 g peptone, 3 g beef extract and 5 g NaCl at 30 °C. For genome sequencing, total genomic DNA was extracted from 10 mL overnight cultures using a DNA extraction kit following the manufacture’s instructions (Aidlab). Quantification and quality control of the genomic DNA were completed by using a Qubit fluorometer (Invitrogen, CA, USA) with Qubit dsDNA BR Assay kit and 0.7% agarose gel electrophoresis with λ-Hind III digest DNA marker.
Genome sequencing and assembly
Complete genome sequencing was performed on an Illumina HiSeq 2500 system by constructing three DNA libraries (a paired-end library with insert size of 491 bp, and two mate pair libraries with insert sizes of 2.5 and 6.9 kb). After filtering low quality and Illumina PCR adapter reads, a total of 1967 Mb clean data were obtained from 2052 Mb raw data. Subsequently, all reads data were denovo assembled into a circular contig with 259 folds of genomic coverage, using SOAPdenovo v.2.04 . Detailed genome sequencing project information is shown in Table 2.
Gene prediction was carried out by GeneMarkS v.4.6 . Function annotation of predicted ORFs was performed based on a BLASTP search against NCBI nonredundant protein database and COG database. Transfer RNAs, rRNAs and sRNA were predicted using tRNAscan-SE v.1.31 with the bacterial model, RNAmmer v.1.2 and Rfam database v.9.1, respectively [19,20,21]. The CRISPRs were identified by using the CRISPR database . The prediction of genes with signal peptides and transmembrane helices were performed by SignalP server v.4.1  and TMHMM server v.2.0 , respectively. The secondary metabolism gene cluster was predicted according to the antiSMASH v.3.0 procedure .
The genome of A. humicireducens SgZ-5T comprises a circular chromosome of 3,181,617 bp and six circular plasmids, designated as pYZ1 (715,112 bp), pYZ2 (1,008,603 bp), pYZ3 (252,411 bp), pYZ4 (338,445 bp), pYZ5 (626,509 bp) and pYZ6 (711,682 bp) (Table 3). The total size of the genome is 6,834,379 bp, and the average GC content is 67.55%. The genome contains 6054 genes with the total length of 5,902,731 bp, of which 5969 (98.6%) are protein coding genes. There are 85 RNA genes (1.4%), including 14 rRNA and 67 tRNA genes. A total of 4844 genes (80.0%) have been assigned a predicted function while the rest have been designated as hypothetical proteins. Genome statistics are summarized in Table 4 and a graphical map is represented in Fig. 3. Furthermore, 4550 (75.2%) genes were assigned to 21 COG functional categories. The distribution of genes into different COG functional categories is provided in Table 5. Six Azospirillum species genomes (including A. humicireducens ) of characterized strains are compared in Table 6. Almost all of these Azospirillum genomes consisting of 6–8 replicons have the total size of 6.5–7.6 Mb and the average GC content of 67.5–70.7%, and contain the total genes in the range of 5951 to 6982 [3, 6, 26, 27]. Furthermore, the main features of A. humicireducens SgZ-5T genome are close to those of A. lipoferum 4B genome.
Insights into the genome sequence
Nitrogen fixation is the major proposed mechanism, by which Azospirillum affects plant growth [2, 4]. A complete set of genes encoding enzymes involved in nitrogen fixation was found in the genomic analysis of A. humicireducens SgZ-5T (Table 7). The main genes involved in this process are nif genes, of which nifDK genes (A6A40_02900 and A6A40_02895) annotated as nitrogenase molybdenum-iron proteins and nifH gene (A6A40_02905) encoding dinitrogenase reductase protein have been identified. In the upstream region of the nifHDK operon, we have found that nifEN genes (A6A40_02875 and A6A40_02870) involved in synthesis of the molybdenum-iron cofactor of nitrogenase are clustered into a single operon together with nifX (A6A40_02865). Furthermore, the genome of A. humicireducens SgZ-5T has nifUSVW genes (A6A40_02235, A6A40_02230, A6A40_02225 and A6A40_02215), which are separated from the structural nifENX operon by about 160 kb.
Organization of the nitrogen fixation gene cluster in A. humicireducens SgZ-5T is presented in Fig. 4. Except for the separately transcribed nifA (A6A40_09040), nifB (A6A40_09050) and nifZ genes (A6A40_09070 and A6A40_09075), all the nif genes have resided in the nitrogen fixation gene cluster of 176.7 kb. Besides, an operon containing fixABCX genes (A6A40_02185, A6A40_02190, A6A40_02195 and A6A40_02220) responsible for electron transfer to nitrogenase is located upstream of this gene cluster. Nevertheless, the fixABCX operon is generally regulated by a transcriptional activator NifA protein for all nitrogen-fixing bacteria in the genus Azospirillum studied so far . Furthermore, draTG genes (A6A40_02920 and A6A40_02925) implicated in posttranslational regulatory process of nitrogenase activity were found in the downstream of and divergently oriented with respect to nifHDK genes. On the whole, the nitrogen fixation gene cluster of A. humicireducens SgZ-5T was in agreement with that in A. brasilense , A. lipoferum and Azospirillum sp. B510 [6, 26, 28, 29], suggesting that nitrogen fixation process demands the systematic action of various genes.
Since tryptophan is a main precursor for biosynthesis of IAA, a well-known phytohormone , the genes in A. humicireducens SgZ-5T related to the production of this amino acid have been analyzed (Additional file 2). The genome harbors three genes trpE, trpG and trpEG (A6A40_04380, A6A40_04655 and A6A40_05775), each encoding the key enzyme anthranilate synthase in tryptophan biosynthesis. Together with trpG, the genes trpD (A6A40_04650) and trpC (A6A40_04645) form a gene cluster of 2.4 kb. Except for anthranilate synthase, this trpGDC gene cluster encodes anthranilate phosphoribosyltransferase and indole-3-glycerol phosphate synthase, which plays a role in synthesis of tryptophan used in multiple biological processes including IAA biosynthesis . The same trpGDC cluster was previously found in A. brasilense . Although the ipdC gene, related to the indole-3-pyruvate pathway for the biosynthesis of IAA , was not discovered in the A. humicireducens SgZ-5T genome, alternative pathway might exist in SgZ-5T. In the genome, A6A40_22745 and A6A40_22755 were assigned as candidates for iaaM and iaaH genes, respectively. These two genes were also found in the Azospirillum sp. B510 genome, and are known to be involved in the IAM pathway for IAA biosynthesis by catalyzing the decarboxylation of tryptophan into IAM and the hydrolysis of IAM to produce IAA [6, 30].
The A. humicireducens SgZ-5T genome also contains a terpene gene cluster of 24.0 kb consisting of 23 genes (A6A40_04945, A6A40_04950, A6A40_04955, …, A6A40_05055) (Additional file 3). This gene cluster encodes a series of proteins, which are involved in the biosynthesis of secondary metabolite production of terpenoid. Thereinto, A6A40_05010 was indentified as the crtB gene, encoding phytoene synthase involved in the biosynthesis of carotenoid. Similar genes in this gene cluster were previously observed in the A. lipoferum 4B genome [7, 26]. Furthermore, some phytohormones including gibberellins and abscisic acid with over 120 types found in plants, fungi, and bacteria, are synthesized through the terpenoid pathway . Therefore, A. humicireducens SgZ-5 exhibits an attractive application as a PGPR likely harboring multiple PGPP in agriculture.
We report here an inventory of the genomic features of the nitrogen-fixing bacterium A. humicireducens SgZ-5T. The genome sequence of strain SgZ-5T revealed further genetic elements involved in nitrogen fixation and its regulation, as well as in the production of phytohormones. We anticipate that knowledge of this genome will contribute to new insights into the mechanisms of plant growth stimulation through genomic comparisons among available complete genomes of Azospirillum strains.
Microbial fuel cell
Plant growth promoting properties
Plant growth-promoting rhizobacteria
Okon Y, Vanderleyden J. Root-associated Azospirillum species can stimulate plants. ASM News. 1997;63:366–70.
Bashan Y, De-Bashan LE. How the plant growth-promoting bacterium Azospirillum promotes plant growth-a critical assessment. Adv Agron. 2010;108:77–136.
Kwak Y, Shin JH. First Azospirillum genome from aquatic environments: whole-genome sequence of Azospirillum thiophilum BV-ST, a novel diazotroph harboring a capacity of sulfur-chemolithotrophy from a sulfide spring. Mar Genomics. 2016;25:21–4.
Okon Y, Heytler PG, Hardy RWF. N2 fixation by Azospirillum brasilense and its incorporation into host Setaria italica. Appl Environ Microbiol. 1983;46:694–7.
Steenhoudt O, Vanderleyden J. Azospirillum, a free-living nitrogen-fixing bacterium closely associated with grasses: genetic, biochemical and ecological aspects. FEMS Microbiol Rev. 2000;24:487–506.
Kaneko T, Minamisawa K, Isawa T, Nakatsukasa H, Mitsui H, Kawaharada Y, et al. Complete genomic structure of the cultivated rice endophyte Azospirillum sp. B510. DNA Res. 2010;17:37–50.
Drogue B, Sanguin H, Borland S. Genome wide profiling of Azospirillum lipoferum 4B gene expression during interaction with rice roots. FEMS Microbiol Ecol. 2014;87:543–55.
Tsavkelova EA, Klimova SY, Cherdyntseva TA, Netrusov AI. Microbial producers of plant growth stimulators and their practical use: a review. Appl Biochem Microbiol. 2006;42:117–26.
Sant’Anna FH, Almeida LGP, Cecagno R, Reolon LA, Siqueira FM, Machado MRS, et al. Genomic insights into the versatility of the plant growth-promoting bacterium Azospirillum amazonense. BMC Genomics. 2011;12:409.
Lin SY, Liu YC, Hameed A, Hsu YH, Huang H, Lai WA, et al. Azospirillum agricola sp. nov., a nitrogen-fixing species isolated from cultivated soil. Int J Syst Evol Microbiol. 2016;66:1453–8.
Zhou SG, Han LC, Wang YQ, Yang GQ, Zhuang L, Hu P. Azospirillum humicireducens sp. nov., a nitrogen-fixing bacterium isolated from a microbial fuel cell. Int J Syst Evol Microbiol. 2013;63:2618–24.
Tamura K, Dudley J, Nei M, Kumar S. MEGA4: molecular evolutionary genetic analysis (MEGA) software version 4.0. Mol Biol Evol. 2007;24:1596–9.
Nakade DB. Halotolerent Azospirillum lipoferum N-29 as a biofertilizer for saline soils. J Pure Appl Microbiol. 2013;7:795–801.
Tamilselvi S, Dutta A, Sindhu M. Azospirillum lipoferum and Pseudomonas fluorescens as effective biological agents for enhanced agro-productivity. Biosci J. 2016;32:670–83.
Raymond J, Siefert JL, Staples CR, Blankenship RE. The natural history of nitrogen fixation. Mol Biol Evol. 2004;21:541–54.
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(5):541–7.
Luo R, Liu B, Xie Y, Li Z, Huang W, Yuan J, et al. SOAPdenovo2: an empirically improved memory-efficient short-read de novo assembler. GigaScience. 2012;1:18.
Besemer J, Lomsadze A, Borodovsky M. GeneMarkS: a self-training method for prediction of gene starts in microbial genomes. Implications for finding sequence motifs in regulatory regions. Nucleic Acids Res. 2001;29:2607–18.
Lowe TM, Eddy SR. tRNAscan-SE: a program for improved detection of transfer RNA genes in genomic sequence. Nucleic Acids Res. 1997;25:955–64.
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.
Gardner PP, Daub J, Tate JG, Nawrocki EP, Kolbe DL, Lindgreen S, et al. Rfam: updates to the RNA families database. Nucleic Acids Res. 2009;37:D136–40.
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(12):4576–9.
Petersen TN, Brunak S, von Heijne G, Nielsen H. Signal 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(3):567–80.
Weber T, Blin K, Duddela S, Krug D, Kim HU, Bruccoleri R, et al. antiSMASH 3.0-a comprehensive resource for the genome mining of biosynthetic gene clusters. Nucleic Acids Res. 2015;43:W237–43.
Wisniewski-Dyé F, Borziak K, Khalsa-Moyers G, Alexandre G, Sukharnikov LO, Wuichet K, et al. Azospirillum genomes reveal transition of Bacteria from aquatic to terrestrial environments. PLoS Genet. 2011;7(12):e1002430.
Rivera D, Revale S, Molina R, Gualpa J, Puente M, Maroniche G, et al. Complete genome sequence of the model rhizosphere strain Azospirillum brasilense az39, successfully applied in agriculture. Genome Announc. 2014;2(4):e00683–14.
Fu HA, Fitzmaurice WP, Roberts GP, Burris RH. Cloning and expression of draTG genes from Azospirillum lipoferum. Gene. 1990;86:95–8.
Zhang Y, Burris RH, Roberts GP. Cloning, sequencing, mutagenesis, and functional characterization of draT and draG genes from Azospirillum brasilense. J Bacteriol. 1992;174:3364–9.
Spaepen S, Vanderleyden J, Remans R. Indole-3-acetic acid in microbial and microorganism-plant signaling. FEMS Microbiol Rev. 2007;31:425–48.
Kang SM, Asaf S, Kim SJ, Yun BW, Lee IJ. Complete genome sequence of plant growth-promoting bacterium Leifsonia xyli SE134, a possible gibberellin and auxin producer. J Biotechnol. 2016;239:34–8.
Zimmer W, Aparicio C, Elmerich C. Relationship between tryptophan biosynthesis and indole-3-acetic acid production in Azospirillum: identification and sequencing of a trpGDC cluster. Mol Gen Genet. 1991;229:41–51.
Garrity GM, Bell JA, Liburn T. Phylum XIV. Proteobacteria phyl. nov. In: Garrity GM, Brenner DJ, Krieg NR, Staley JT, editors. Bergey’s manual of systematic bacteriology, volume 2, part B. 2nd ed. New York: Springer; 2005. p. 1.
Garrity GM, Bell JA, Liburn T. Class I. Alphaproteobacteria class. nov. In: Brenner DJ, Krieg NR, Staley JT, Garrity GM, editors. Bergey’s Manual of systematic bacteriology, volume 2, part C. 2nd ed. New York: Springer; 2005. p. 1.
Pfennig N, Trüper HG. Higher taxa of the phototrophic bacteria. Int J Syst Bacteriol. 1971;21:17–8.
Skerman VBD, McGowan V, Sneath PHA. Approved lists of bacterial names. Int J Syst Bacteriol. 1980;30:225–420.
Tarrand JJ, Krieg NR, Döbereiner J. A taxonomic study of the Spirillum lipoferum group, with descriptions of a new genus, Azospirillum gen. nov. and two species, Azospirillum lipoferum (Beijerinck) comb. nov. and Azospirillum brasilense sp. nov. Can J Microbiol. 1978;24:967–80.
List Editor. Validation list no. 2. Validation of the publication of new names and new combinations previously effectively published outside the IJSB. Int J Syst Bacteriol. 1979;29:79–80.
Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, et al. Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. Nat Genet. 2000;25:25–9.
This work was supported by the Guangdong Academy of Sciences Funds for Innovation Driven Development, China (2017GDASCX-0409), the National Natural Science Foundation of China (41501546), the Guangdong Natural Science Foundation, China (2016A030313779), and the Science and Technology Planning Project of Guangdong, China (2017A030303057).
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Phylogenetic tree based on the partial nifH gene sequences showing the position of A. humicireducens SgZ-5T relative to other species within the genus Azospirillum and related genus. The strains and their corresponding GenBank accession numbers of nifH gene were indicated in parentheses. The sequences were aligned using Clustal W and the neighbor-joining tree was constructed based on kimura 2-paramenter distance model by using MEGA 5. Bootstrap values above 50% were obtained from 1000 bootstrap replications. Bar, 0.01 substitutions per nucleotide position. Leptospirillum ferriphilum YSKT was used as an outgroup. (DOCX 64 kb)
Genes of A. humicireducens SgZ-5T involved in biosynthesis of tryptophan. (DOCX 16 kb)