Differential expression of phosphorus acquisition genes in response to phosphorus stress in two Raphidiopsis raciborskii strains

The cyanobacterium Raphidiopsis raciborskii is a nuisance in freshwater ecosystems. Strains vary in their physiological responses to environmental drivers, thus a greater understanding of the magnitude of strain variation is required to characterize the species. In this study, two strains of R. raciborskii isolated from a tropical Australian water reservoir were grown with and without phosphorus (P) to determine any relative response to P stress. The strains had the same growth rates and under P free conditions, cells grew at the same rate as P replete conditions until day 9 when cell growth ceased. There was no difference in the alkaline phosphatase activity per cell for the P replete and P free conditions, but the level of activity per cell was greater in CS-505 than CS-506. P acquisition genes were identified from the sequenced genomes; both strains contained the same genes, but with differences in copy number of phoA (7 and 6), phnK (3 and 1) and phnH (2 and 1) between CS-505 and CS-506 (respectively). The expression of P acquisition genes under P stress was measured throughout the experiment and shown to vary in magnitude and timing across strains, and in P replete versus P free cultures. In strain CS-505, upregulation of the pstS1 and phoA genes occurred late in the growth phase and into senescence. These genes are involved in phosphate uptake and use of various forms of organic P. For strain CS-506, there was upregulation of the phosphate uptake gene, pit, and organic P utilization genes, phoA, phoU, phnD and phnK, commencing late in the growth phase. Our study shows that despite the fact that these two strains were isolated from the same waterbody, they differed markedly in their gene expression response to P free conditions. This capacity of R. raciborskii to vary in strain responses to P conditions gives the organism flexibility in responding to environmental change, particularly P stress conditions.     

Localized Tetrazolium Reduction in Relation to N2 Fixation, CO2 Fixation, and H2 Uptake in Aquatic Filamentous Cyanobacteria

The aquatic filamentous cyanobacteria Anabaena oscillarioides and Trichodesmium sp. reveal specific cellular regions of tetrazolium salt reduction. The effects of localized reduction of five tetrazolium salts on N₂ fixation (acetylene reduction), ¹⁴CO₂ fixation, and ³H₂ utilization were examined. During short-term (within 30 min) exposures in A. oscillarioides, salt reduction in heterocysts occurred simultaneously with inhibition of acetylene reduction. Conversely, when salts failed to either penetrate or be reduced in heterocysts, no inhibition of acetylene reduction occurred. When salts were rapidly reduced in vegetative cells, ¹⁴CO₂ fixation and ³H₂ utilization rates decreased, whereas salts exclusively reduced in heterocysts were not linked to blockage of these processes. In the nonheterocystous genus Trichodesmium, the deposition of reduced 2,3,5-triphenyl-2-tetrazolium chloride (TTC) in the internal cores of trichomes occurs simultaneously with a lowering of acetylene reduction rates. Since TTC deposition in heterocysts of A. oscillarioides occurs contemporaneously with inhibition of acetylene reduction, we conclude that the cellular reduction of this salt is of use in locating potential N₂-fixing sites in cyanobacteria. The possible applications and problems associated with interpreting localized reduction of tetrazolium salts in cyanobacteria are presented.     

Constitutive toxin production under various nitrogen and phosphorus regimes of three ecotypes of Cylindrospermopsis raciborskii ((Wołoszyńska) Seenayya et Subba Raju)

Cylindrospermopsis raciborskii is a global invasive cyanobacterium, with some ecotypes (i.e. strains) producing the toxin cylindrospermopsin, CYN. Multiple ecotypes can co-exist, complicating prediction of toxin concentrations based on cell concentrations. This study examined the growth response and toxin production of three Australian ecotypes of C. raciborskii, two toxic (CS-505, CS-506) and one non-toxic (CS-510), to a range of nitrogen (N) and phosphorus (P) concentrations. CYN cell quota was constant under all N:P ratios and concentration conditions, indicative of a constitutive response, yet the CYN cell quota was 6-fold higher in CS-506 compared to CS-505. The ecotypes differed in response to dissolved N depletion: there was a 4-fold difference in the number of cells heterocyst mL⁻¹ between CS-505 and CS-510, while CS-506 did not produce any heterocysts and was unable to grow in N deplete conditions. Growth rates were lower for all ecotypes as [P] increased, indicative of a species with a strategy of P storage rather than increased growth. Presumably this is an adaptation to low and fluctuating P conditions. However, the negative effect of increasing [P] on growth is surprising. In contrast, increasing [N] resulted in higher growth rates across ecotypes. This study highlights the importance of understanding differences in growth and toxin production between ecotypes in response to environmental conditions in order to more effectively predict blooms and toxin yields.     

Comparison of the terrestrial cyanobacterium Leptolyngbya sp. NIES-2104 and the freshwater Leptolyngbya boryana PCC 6306 genomes

The cyanobacterial genus Leptolyngbya is widely distributed throughout terrestrial environments and freshwater. Because environmental factors, such as oxygen level, available water content, and light intensity, vary between soil surface and water bodies, terrestrial Leptolyngbya should have genomic differences with freshwater species to adapt to a land habitat. To study the genomic features of Leptolyngbya species, we determined the complete genome sequence of the terrestrial strain Leptolyngbya sp. NIES-2104 and compared it with that of the near-complete sequence of the freshwater Leptolyngbya boryana PCC 6306. The greatest differences between these two strains were the presence or absence of a nitrogen fixation gene cluster for anaerobic nitrogen fixation and several genes for tetrapyrrole synthesis, which can operate under micro-oxic conditions. These differences might reflect differences in oxygen levels where these strains live. Both strains have the genes for trehalose biosynthesis, but only Leptolyngbya sp. NIES-2104 has genetic capacity to produce a mycosporine-like amino acid, mycosporine-glycine. Mycosporine-glycine has an antioxidant action, which may contribute to adaptation to terrestrial conditions. These features of the genomes yielded additional insights into the classification and physiological characteristics of these strains.     

Dynamics and Cell-Type Specificity of the DNA Double-Strand Break Repair Protein RecN in the Developmental Cyanobacterium Anabaena sp. Strain PCC 7120

DNA replication and repair are two fundamental processes required in life proliferation and cellular defense and some common proteins are involved in both processes. The filamentous cyanobacterium Anabaena sp. strain PCC 7120 is capable of forming heterocysts for N₂ fixation in the absence of a combined-nitrogen source. This developmental process is intimately linked to cell cycle control. In this study, we investigated the localization of the DNA double-strand break repair protein RecN during key cellular events, such as chromosome damaging, cell division, and heterocyst differentiation. Treatment by a drug causing DNA double-strand breaks (DSBs) induced reorganization of the RecN focus preferentially towards the mid-cell position. RecN-GFP was absent in most mature heterocysts. Furthermore, our results showed that HetR, a central player in heterocyst development, was involved in the proper positioning and distribution of RecN-GFP. These results showed the dynamics of RecN in DSB repair and suggested a differential regulation of DNA DSB repair in vegetative cell and heterocysts. The absence of RecN in mature heterocysts is compatible with the terminal nature of these cells.     

A meiotically reproducible chromosome length polymorphism in the ascomycete fungus Ophiostoma ulmi (sensu lato)

We have followed the transmission of Ophiostoma ulmis.l. chromosome length polymorphisms (CLPs) into the F₂ generation to determine the reproducibility of a genome rearrangement culminating in the conversion of a 1.0 Mb chromosome into a 800 kb chromosome. The 1.0 Mb chromosome in strain CESS16K is thus far unique among O. ulmi s.l. wild-type strains, as no other wild-type strains have been observed with chromosomes smaller than 2.3 Mb. It has been previously shown that the 1.0 Mb chromosome is mitotically stable, carries at least one normally expressed gene, and is transmitted through meiosis. In this study, a series of crosses were performed to further elucidate the pattern of inheritance of the 1.0 Mb chromosome and the process of conversion of the 1.0 Mb species to 800 kb. In crosses where the 1.0 Mb chromosome was allowed to pair with itself or with the 800 kb chromosome, all progeny inherited a copy of the 1.0 Mb or 800 kb form, further demonstrating the A-type nature of these small chromosomes. When a cross was repeated between the strains CESS16K (1.0 Mb chromosome) and FG245B^r-O (no 1.0 Mb or 800 kb chromosome), the occurrence of a 800 kb chromosome was observed in 9% of the progeny. A reciprocal cross between an 800 kb strain and a strain with no 800 kb or 1.0 Mb chromosome was conducted, and a progeny strain containing a 1.0 Mb chromosome was recovered. The reproducibility and reciprocality of the 1.0 Mb to 800 kb chromosome conversion demonstrates that meiotic processes are responsible for this CLP, and that O. ulmi s.l. strains with various divergent genome architectures can remain sexually compatible.     

Morphological and habitat evolution in the Cyanobacteria using a compartmentalization approach

A phylogenomic approach was used to study the evolution of traits in the Cyanobacteria. A cyanobacterial backbone tree was constructed using multiple concatenated sequences from whole genome sequences. Additional taxa were added using a separate alignment that contained morphological characters, SSU (small subunit) and LSU (large subunit) rDNA, rpoC, rpoD, tufA, and gyrB genes. A compartmentalization approach was then used to construct a robust phylogeny with resolved deep branches. Additional morphological characters (e.g. unicellular or filamentous growth, presence or absence of heterocysts) were coded, mapped onto the backbone cyanobacterial tree, and the ancestral character states inferred. Our analyses show that the earliest cyanobacterial lineages were likely unicellular coccoid/ellipsoidal/short rods that lived in terrestrial/freshwater environments. Later cyanobacterial lineages independently gained the ability to colonize brackish, marine, and hypersaline environments while acquiring a large number of more complex traits: sheath, filamentous growth, nitrogen fixation, thermophily, motility, and use of sulphide as an electron donor. Many of these adaptations would have been important in the appearance of dense microbial mats early in Earth's history. Complex traits such as hormogonia, heterocysts, and akinetes had a single ancestor. Within the Nostocales, hormogonia and heterocysts arose before akinetes.     

3种水稻土中7株固氮蓝细菌的分离与特征

【背景】蓝细菌是水生和陆地生态系统中生物固氮的主要贡献者。【目的】增加对稻田土壤固氮蓝细菌的了解,获得用于进一步研究的可培养固氮蓝细菌菌株。【方法】选择3种具有不同固氮能力的水稻土,采用BG11-N培养基分离培养固氮蓝细菌菌株,对新分离菌株进行形态特征观察,通过基因组DNA的nifH基因扩增明确其固氮潜力,进一步采用乙炔还原法和~(15)N_2示踪法定量测定其固氮能力,通过基因组DNA的16SrRNA基因序列比对进行鉴定。【结果】在光照培养条件下,采用BG11-N培养基共分离纯化得到自养菌株7株,细胞呈圆形或椭圆形、单列、无分枝、丝状和念珠状,在固体培养基上形成团垫状菌落。新分离菌株在BG11-N培养基中生长状况良好,以基因组DNA为模板可扩增出nifH基因,乙炔还原法和~(15)N_2示踪法测定结果显示具有较高固氮能力,同时具有铁载体生成能力。结合16S rRNA基因序列比对和形态特征,7株菌被初步鉴定隶属于念珠藻科(Nostocaceae)。【结论】从水稻土中分离到在稻田生物固氮中发挥重要作用的蓝细菌(念珠藻科)菌株,可培养固氮蓝细菌菌株固氮能力较高,兼具铁载体生成能力,可作为进一步深入研究的微生物资源,具有潜在的研究应用价值。     

Pathogenicity,Phylogenetic relationship and NGS based identification and assembly of tumorigenic Agrobacterium radiabacter plasmidic and chromosomic reads isolated from Prunus duclcis

Although many Agrobacterium radiobacter strains have already been identified, only a few genomes of strains belonging to genomovar G4 have been sequenced so far. In this study, we report the first virulent genome sequence of Agrobacterium radiobacter strain tun 183, which is highly virulent to almond specie. The genome size was estimated to be 5.53 Mb, with 57.9%GC content. In total, 6486 genes encoding proteins and 61 genes encoding RNAs were identified in this genome. Comparisons with the available sequenced genomes of genomovar G4 as well as with other A. sp. were conducted, revealing a hexapartite genome containing circular and linear chromosomes in addition to two accessory plasmids and a tumor inducing plasmid (pTi) in strain tun 183. The phylogenetic analysis of recA gene clearly showed the clustering of tun 183 strain within genomovar G4, supporting the monophyly within this genomovar.     

Genomic architecture of three newly isolated unclassified Butyrivibrio species elucidate their potential role in the rumen ecosystem

One cellulose-degrading strain CB08 and two xylan-degrading strains XB500–5 and X503 were isolated from buffalo rumen. All the strains were designated as putative novel species of Butyrivibrio based on phylogeny, phylogenomy, digital DNA-DNA hybridization, and average nucleotide identity with their closest type strains. The draft genome length of CB08 was ~3.54 Mb, while X503 and XB500-5 genome sizes were ~3.24 Mb and ~3.27 Mb, respectively. Only 68.28% of total orthologous clusters were shared among three genomes, and 40–44% of genes were identified as hypothetical proteins. The presence of genes encoding diverse carbohydrate-active enzymes (CAZymes) exhibited the lignocellulolytic potential of these strains. Further, the genome annotations revealed the metabolic pathways for monosaccharide fermentation to acetate, butyrate, lactate, ethanol, and hydrogen. The presence of genes for chemotaxis, antibiotic resistance, antimicrobial activity, synthesis of vitamins, and essential fatty acid suggested the versatile metabolic nature of these Butyrivibrio strains in the rumen environment.     

Genome Analysis Coupled with Physiological Studies Reveals a Diverse Nitrogen Metabolism in Methylocystis sp. Strain SC2

Background

Methylocystis sp. strain SC2 can adapt to a wide range of methane concentrations. This is due to the presence of two isozymes of particulate methane monooxygenase exhibiting different methane oxidation kinetics. To gain insight into the underlying genetic information, its genome was sequenced and found to comprise a 3.77 Mb chromosome and two large plasmids.

Principal Findings

We report important features of the strain SC2 genome. Its sequence is compared with those of seven other methanotroph genomes, comprising members of the Alphaproteobacteria, Gammaproteobacteria, and Verrucomicrobia. While the pan-genome of all eight methanotroph genomes totals 19,358 CDS, only 154 CDS are shared. The number of core genes increased with phylogenetic relatedness: 328 CDS for proteobacterial methanotrophs and 1,853 CDS for the three alphaproteobacterial Methylocystaceae members, Methylocystis sp. strain SC2 and strain Rockwell, and Methylosinus trichosporium OB3b. The comparative study was coupled with physiological experiments to verify that strain SC2 has diverse nitrogen metabolism capabilities. In correspondence to a full complement of 34 genes involved in N₂ fixation, strain SC2 was found to grow with atmospheric N₂ as the sole nitrogen source, preferably at low oxygen concentrations. Denitrification-mediated accumulation of 0.7 nmol ³⁰N₂/hr/mg dry weight of cells under anoxic conditions was detected by tracer analysis. N₂ production is related to the activities of plasmid-borne nitric oxide and nitrous oxide reductases.

Conclusions/Perspectives

Presence of a complete denitrification pathway in strain SC2, including the plasmid-encoded nosRZDFYX operon, is unique among known methanotrophs. However, the exact ecophysiological role of this pathway still needs to be elucidated. Detoxification of toxic nitrogen compounds and energy conservation under oxygen-limiting conditions are among the possible roles. Relevant features that may stimulate further research are, for example, absence of CRISPR/Cas systems in strain SC2, high number of iron acquisition systems in strain OB3b, and large number of transposases in strain Rockwell.     

Observation of microplasmodesmata in both heterocyst-forming and non-heterocyst forming filamentous cyanobacteria by freeze-fracture electron microscopy

Structures which may establish cytoplasmic continuity between adjacent cells of filamentous cyanobacteria have been observed by freeze-fracture electron microscopy. They are visible in the septum region of the plasma membrane as pits on the E-face (EF) and corresponding protrusions on the P-face (PF). Between 100 and 250 of these structures, termed microplasmodesmata, were present between adjacent vegetative cells in all four strains of heterocyst-forming filamentous cyanobacteria, Anabaena cylindrica Lemm, A. variabilis (IUCC B377), A. variabilis Kütz. (ATCC 29413) and Nostoc muscorum, examined. Only 30–40 microplasmodesmata were observed between adjacent cells in two species, Phormidium luridum and Plectonema boryanum, that do not form heterocysts. The results suggest that in species that form heterocysts a greater degree of cytoplasmic continuity is established, presumably to facilitate the exchange of metabolites. In species capable of forming heterocysts, the number of microplasmodesmata per septum between two adjacent vegetative cells remained constant whether the filaments were grown in the presence of NH₄ and lacked heteroxysts or under N₂-fixing conditions and contained heterocysts. When a vegetative cell differentiates into a heterocyst, about 80% of the existing microplasmodesmata are destroyed as the poles of the cell become constricted into narrow necks leaving smaller areas of contact with the adjacent vegetative cells.     

Comparative genomics of Bradyrhizobium japonicum CPAC 15 and Bradyrhizobium diazoefficiens CPAC 7: elite model strains for understanding symbiotic performance with soybean

Background

The soybean-Bradyrhizobium symbiosis can be highly efficient in fixing nitrogen, but few genomic sequences of elite inoculant strains are available. Here we contribute with information on the genomes of two commercial strains that are broadly applied to soybean crops in the tropics. B. japonicum CPAC 15 (=SEMIA 5079) is outstanding in its saprophytic capacity and competitiveness, whereas B. diazoefficiens CPAC 7 (=SEMIA 5080) is known for its high efficiency in fixing nitrogen. Both are well adapted to tropical soils. The genomes of CPAC 15 and CPAC 7 were compared to each other and also to those of B. japonicum USDA 6^T and B. diazoefficiens USDA 110^T.

Results

Differences in genome size were found between species, with B. japonicum having larger genomes than B. diazoefficiens. Although most of the four genomes were syntenic, genome rearrangements within and between species were observed, including events in the symbiosis island. In addition to the symbiotic region, several genomic islands were identified. Altogether, these features must confer high genomic plasticity that might explain adaptation and differences in symbiotic performance. It was not possible to attribute known functions to half of the predicted genes. About 10% of the genomes was composed of exclusive genes of each strain, but up to 98% of them were of unknown function or coded for mobile genetic elements. In CPAC 15, more genes were associated with secondary metabolites, nutrient transport, iron-acquisition and IAA metabolism, potentially correlated with higher saprophytic capacity and competitiveness than seen with CPAC 7. In CPAC 7, more genes were related to the metabolism of amino acids and hydrogen uptake, potentially correlated with higher efficiency of nitrogen fixation than seen with CPAC 15.

Conclusions

Several differences and similarities detected between the two elite soybean-inoculant strains and between the two species of Bradyrhizobium provide new insights into adaptation to tropical soils, efficiency of N₂ fixation, nodulation and competitiveness.

Electronic supplementary material

The online version of this article (doi: 10.1186/1471-2164-15-420) contains supplementary material, which is available to authorized users.     

Genome-Wide and Heterocyst-Specific Circadian Gene Expression in the Filamentous Cyanobacterium Anabaena sp. Strain PCC 7120

The filamentous, heterocystous cyanobacterium Anabaena sp. strain PCC 7120 is one of the simplest multicellular organisms that show both morphological pattern formation with cell differentiation (heterocyst formation) and circadian rhythms. Therefore, it potentially provides an excellent model in which to analyze the relationship between circadian functions and multicellularity. However, detailed cyanobacterial circadian regulation has been intensively analyzed only in the unicellular species Synechococcus elongatus. In contrast to the highest-amplitude cycle in Synechococcus, we found that none of the kai genes in Anabaena showed high-amplitude expression rhythms. Nevertheless, ∼80 clock-controlled genes were identified. We constructed luciferase reporter strains to monitor the expression of some high-amplitude genes. The bioluminescence rhythms satisfied the three criteria for circadian oscillations and were nullified by genetic disruption of the kai gene cluster. In heterocysts, in which photosystem II is turned off, the metabolic and redox states are different from those in vegetative cells, although these conditions are thought to be important for circadian entrainment and timekeeping processes. Here, we demonstrate that circadian regulation is active in heterocysts, as shown by the finding that heterocyst-specific genes, such as all1427 and hesAB, are expressed in a robust circadian fashion exclusively without combined nitrogen.     

Analysis of the hli gene family in marine and freshwater cyanobacteria

Certain cyanobacteria thrive in natural habitats in which light intensities can reach 2000 micromol photon m(-2) s(-1) and nutrient levels are extremely low. Recently, a family of genes designated hli was demonstrated to be important for survival of cyanobacteria during exposure to high light. In this study we have identified members of the hli gene family in seven cyanobacterial genomes, including those of a marine cyanobacterium adapted to high-light growth in surface waters of the open ocean (Prochlorococcus sp. strain Med4), three marine cyanobacteria adapted to growth in moderate- or low-light (Prochlorococcus sp. strain MIT9313, Prochlorococcus marinus SS120, and Synechococcus WH8102), and three freshwater strains (the unicellular Synechocystis sp. strain PCC6803 and the filamentous species Nostoc punctiforme strain ATCC29133 and Anabaena sp. [Nostoc] strain PCC7120). The high-light-adapted Prochlorococcus Med4 has the smallest genome (1.7 Mb), yet it has more than twice as many hli genes as any of the other six cyanobacterial species, some of which appear to have arisen from recent duplication events. Based on cluster analysis, some groups of hli genes appear to be specific to either marine or freshwater cyanobacteria. This information is discussed with respect to the role of hli genes in the acclimation of cyanobacteria to high light, and the possible relationships among members of this diverse gene family.     

Stratification,nitrogen fixation,and cyanobacterial bloom stage regulate the planktonic food web structure

Changes in the complexity of planktonic food webs may be expected in future aquatic systems due to increases in sea surface temperature and an enhanced stratification of the water column. Under these conditions, the growth of unpalatable, filamentous, N₂‐fixing cyanobacterial blooms, and their effect on planktonic food webs will become increasingly important. The planktonic food web structure in aquatic ecosystems at times of filamentous cyanobacterial blooms is currently unresolved, with discordant lines of evidence suggesting that herbivores dominate the mesozooplankton or that mesozooplankton organisms are mainly carnivorous. Here, we use a set of proxies derived from amino acid nitrogen stable isotopes from two mesozooplankton size fractions to identify changes in the nitrogen source and the planktonic food web structure across different microplankton communities. A transition from herbivory to carnivory in mesozooplankton between more eutrophic, near‐coastal sites and more oligotrophic, offshore sites was accompanied by an increasing diversity of microplankton communities with aging filamentous cyanobacterial blooms. Our analyses of 124 biotic and abiotic variables using multivariate statistics confirmed salinity as a major driver for the biomass distribution of non‐N₂‐fixing microplankton species such as dinoflagellates. However, we provide strong evidence that stratification, N₂ fixation, and the stage of the cyanobacterial blooms regulated much of the microplankton diversity and the mean trophic position and size of the metabolic nitrogen pool in mesozooplankton. Our empirical, macroscale data set consistently unifies contrasting results of the dominant feeding mode in mesozooplankton during blooms of unpalatable, filamentous, N₂‐fixing cyanobacteria by identifying the at times important role of heterotrophic microbial food webs. Thus, carnivory, rather than herbivory, dominates in mesozooplankton during aging and decaying cyanobacterial blooms with hitherto uncharacterized consequences for the biogeochemical functions of mesozooplankton.     

Comparative Genomic Analysis of N2-Fixing and Non-N2-Fixing Paenibacillus spp.: Organization,Evolution and Expression of the Nitrogen Fixation Genes

We provide here a comparative genome analysis of 31 strains within the genus Paenibacillus including 11 new genomic sequences of N₂-fixing strains. The heterogeneity of the 31 genomes (15 N₂-fixing and 16 non-N₂-fixing Paenibacillus strains) was reflected in the large size of the shell genome, which makes up approximately 65.2% of the genes in pan genome. Large numbers of transposable elements might be related to the heterogeneity. We discovered that a minimal and compact nif cluster comprising nine genes nifB, nifH, nifD, nifK, nifE, nifN, nifX, hesA and nifV encoding Mo-nitrogenase is conserved in the 15 N₂-fixing strains. The nif cluster is under control of a σ⁷⁰-depedent promoter and possesses a GlnR/TnrA-binding site in the promoter. Suf system encoding [Fe–S] cluster is highly conserved in N₂-fixing and non-N₂-fixing strains. Furthermore, we demonstrate that the nif cluster enabled Escherichia coli JM109 to fix nitrogen. Phylogeny of the concatenated NifHDK sequences indicates that Paenibacillus and Frankia are sister groups. Phylogeny of the concatenated 275 single-copy core genes suggests that the ancestral Paenibacillus did not fix nitrogen. The N₂-fixing Paenibacillus strains were generated by acquiring the nif cluster via horizontal gene transfer (HGT) from a source related to Frankia. During the history of evolution, the nif cluster was lost, producing some non-N₂-fixing strains, and vnf encoding V-nitrogenase or anf encoding Fe-nitrogenase was acquired, causing further diversification of some strains. In addition, some N₂-fixing strains have additional nif and nif-like genes which may result from gene duplications. The evolution of nitrogen fixation in Paenibacillus involves a mix of gain, loss, HGT and duplication of nif/anf/vnf genes. This study not only reveals the organization and distribution of nitrogen fixation genes in Paenibacillus, but also provides insight into the complex evolutionary history of nitrogen fixation.