Previous Article | Next Article ![]()
Applied and Environmental Microbiology, August 2009, p. 5037-5046, Vol. 75, No. 15
0099-2240/09/$08.00+0 doi:10.1128/AEM.00398-09
Copyright © 2009, American Society for Microbiology. All Rights Reserved.
,
Toshihiro Aono,
*
Noriko Akiba,
Kyung-Bum Lee,
Chi-Te Liu,¶
Hiroki Toyazaki, and
Hiroshi Oyaizu
Biotechnology Research Center, The University of Tokyo, 1-1-1 Yayoi, Bunkyo-ku, Tokyo 113-8657, Japan
Received 17 February 2009/ Accepted 8 June 2009
|
|
|---|
|
|
|---|
Azorhizobium caulinodans ORS571 is a microsymbiont of Sesbania rostrata (18-20). Nitrogen-fixing nodules are formed by A. caulinodans on the stems as well as on the roots of S. rostrata. Stem nodules occur at the site of adventitious root primordia located on the stems after crack-entry invasion by A. caulinodans (70). During crack-entry invasion, bacteria proliferate in the epidermal fissures at the adventitious root primordia on the stems (27). Cortical infection pockets are formed by Nod factor-dependent induction of local cell death and subsequent colonization by bacteria (13). From the infection pockets, infection threads guide bacteria toward the cells in the nodule primordia for symbiotic uptake (14).
A. caulinodans is taxonomically different from other rhizobia and belongs to the family Xanthobacteraceae, including Xanthobacter autotrophicus, although it is relatively close to Bradyrhizobium (44). Recently, the whole-genome sequence of A. caulinodans was determined (43), and it was revealed to consist of a single circular chromosome of 5.37 Mb, which is the smallest among the sequenced genomes of rhizobia. The size of the symbiosis island is only 87.6 kb and contains nod genes but not nif and fix genes. With the analysis of the genome sequence, we concurrently performed a large-scale screening of rhizobial factors involved in nodule development using 10,080 mutants of A. caulinodans created by random Tn5 mutagenesis. We analyzed a total of 108 Tn5 mutants and identified 86 genes involved in symbiosis (68). The list of mutants, which showed aberrant nodule development at various stages, shows the complexity of nodule development and maturation. Next, we investigated which genes in A. caulinodans were important for establishing symbiosis with host plants. To this end, we constructed a whole-genome microarray of A. caulinodans and analyzed the gene expression in free-living cells versus that in bacteroids isolated from stem nodules. In addition, we examined the effect of a flavonoid, naringenin, which is one of the most efficient inducers of A. caulinodans nod genes (26), on gene expression in free-living cells. In the present report, we discuss the symbiotic mechanism of A. caulinodans based on the results of the transcriptional profiling experiments described herein and our previous transposon mutagenesis study.
|
|
|---|
Plant growth, inoculation, and bacteroid isolation.
S. rostrata seeds were treated with concentrated sulfuric acid for 1 h, rinsed with sterile water, and soaked in sterile water on trays. The trays were stored for 3 days at 37°C under dark conditions. After germination, S. rostrata were grown for 2 weeks before bacterial inoculation at 35°C under a 24-h light regime at an intensity of 50,000 lx, as described in our previous report (68). Bacterial cultures grown overnight were inoculated onto the stems. At 7 days postinoculation, the stem nodules on the second stem internode of plants were harvested, immediately frozen in liquid nitrogen, and stored at –80°C. To isolate bacteroids, the stem nodules were vigorously homogenized with MMS buffer (40 mM MOPS [morpholinepropanesulfonic acid], 20 mM KOH, 2 mM MgSO4, and 0.3 M sucrose [pH 7.0]) (55) containing 4% polyvinylpyrrolidone. The resulting homogenate was passed through three layers of Miracloth (Calbiochem), and the filtrate was centrifuged at 500 x g for 5 min to remove plant cell fractions and other debris. Then, the supernatant was centrifuged at 5,000 x g for 10 min, and the pellet containing bacteroids was collected. To remove the green layer containing the contaminated chloroplasts from the surface of the white pellet, MMS buffer containing 25% polyvinylpyrrolidone was added onto the pellet, and the surface of the pellet was washed by vigorous pipetting; the buffer was then removed. After this washing procedure, the pellet was resuspended in MMS buffer without polyvinylpyrrolidone and centrifuged at 5,000 x g for 2 min, and the supernatant was removed. This procedure was performed twice. The white pellet was used as a bacteroid sample for RNA isolation.
A. caulinodans ORS571 DNA microarray.
The AzcNXa520550F whole-genome array (designed by Affymetrix) contains 7,096 probe pair sets (45 for target preparation controls, 4,725 for ORFs, and 2,326 for intergenic regions tiled in both directions) corresponding to 4,703 ORFs and 1,163 intergenic regions of A. caulinodans ORS571, based on the genome sequence (43). The nucleotide sequence of the entire genome of A. caulinodans ORS571 is available in the DDBJ/EMBL/GenBank databases under the accession number AP009384. The probe length was 25 mer, and for all ORFs, 16 probe pairs were chosen. Every probe pair set contained 16 perfect-match probes and mismatch probes. For the intergenic sequences, probes were selected at a density of 16 to 24 bp to fill the remaining space on the chip. As target preparation controls, the poly(A) controls (dap, lys, phe, thr, and trp from Bacillus subtilis) were included in the gene chip. Furthermore, the hybridization controls (bioB, bioC, and bioD from Escherichia coli and cre from the P1 bacteriophage) were also included.
RNA isolation, cDNA synthesis, and labeling.
Total RNA was isolated from cultured and bacteroid cells using Isogen (Nippon Gene, Japan) according to the manufacturer's instructions. The isolated RNA was treated with DNase I (Invitrogen) and purified by a phenol-chloroform-isoamyl alcohol (25:24:1) extraction and ethanol precipitation. The purified RNA was analyzed by agarose gel electrophoresis, and distinct bands of 16S and 23S rRNAs were detected. In the case of RNA samples from bacteroids, bands of 18S and 28S rRNAs of plants were not detected. To synthesize cDNA, 10 µg of total RNA was mixed with 2 µl of the 260-fold diluted poly(A) target preparation control containing dap, lys, phe, and thr mRNA (Affymetrix) and was reverse transcribed by SuperScript III transcriptase (Invitrogen). After the remaining RNA was digested with RNase H for 20 min at 37°C, cDNA was collected with the MinElute kit (Qiagen, Germany). Fragmentation was performed with DNase I with 1x One-Phor-All buffer (GE Healthcare, United Kingdom). Fragmented cDNA was end labeled with terminal deoxynucleotidyl transferase (Promega) and GeneChip DNA labeling reagent (Affymetrix).
Hybridization, scanning, and microarray data analyses.
Hybridization, washing, staining, and scanning were performed using the GeneChip hybridization, wash, and stain kit (Affymetrix), the GeneChip Hybridization Oven 640 (Affymetrix), the GeneChip Fluidics Station 450 (Affymetrix), and the GeneChip Scanner 3000 (Affymetrix) according to the Affymetrix manual. For hybridization, 1.0 to 1.5 µg of labeled cDNA was mixed with Control Oligo B2 (3 nM; Affymetrix) in a volume of 150 µl hybridization solution. Hybridization was performed overnight at 50°C for 16 h. Primary signal extraction and comparative analyses were performed using GeneChip Operating Software (GCOS) version 1.4 (Affymetrix). To normalize signals between arrays, global scaling was performed by Microarray Suite version 5.0 (Affymetrix) using all probe sets. The target intensity value was set to 500, and the scaling factors were between 1.5 and 4.9. Default statistical parameters of GCOS (Alpha1 = 0.04; Alpha2 = 0.06; Tau = 0.015; Gamma1L = 0.001111; Gamma1H = 0.001111; Gamma2L = 0.001333; Gamma2H = 0.001333; perturbation = 1.1) were applied. A probe pair set was called "present" (i.e., expressed) when the detection P value was lower than the Alpha1 value, i.e., when the intensity of the perfect-match probe cell was significantly greater than that of the corresponding mismatch probe cell. When the detection P value was greater than or equal to the Alpha2 value, i.e., when the intensity of the perfect-match probe cell was significantly less than that of the corresponding mismatch probe cell, the probe pair set was called "absent" (i.e., not expressed). When the detection P value was between Alpha1 and Alpha2, the probe pair set was called "marginal" (i.e., unable to call the transcript present or absent). Data were then processed using Microsoft Excel. For statistical comparison, a Student's t test was applied to the log2 signal with a P value threshold of
0.05. To evaluate the differentially expressed genes, the weighted average difference (WAD) method (36, 37) was also applied to the log2 signal. Genes were considered differentially expressed if |WAD| was
0.5.
Computational sequence analyses.
Homology searches based on amino acid sequences were performed using the BlastP programs on the National Center for Biotechnology Information (NCBI) server (www.ncbi.nlm.nih.gov/BLAST/). Searches for protein signatures were performed using the InterProScan program on the European Bioinformatics Institute (EBI) server (www.ebi.ac.uk/InterProScan/).
Microarray data accession number.
The entire set of microarray data was deposited in the Center for Information Biology Gene Expression Database (CIBEX) (http://cibex.nig.ac.jp/index.jsp) under the accession number CBX88.
|
|
|---|
0.05; |WAD|
0.5) in each functional category defined in our previous genome sequence study (43). Figure 1 shows the WAD values for each gene plotted against position within the genome for each pairwise comparison. |
View this table: [in a new window] |
TABLE 1. Average correlation of replicate microarray hybridizations for each experimental condition
|
|
View this table: [in a new window] |
TABLE 2. Overview of differential expression patterns based on gene functional categories
|
![]() View larger version (50K): [in a new window] |
FIG. 1. WAD values for each gene plotted against position within the genome, for each pairwise comparison. (A) Expressional difference between bacteroid and free-living cells grown in TY medium. (B) Expressional difference between bacteroid and free-living cells grown in MMO medium. (C) Expressional difference between free-living cells grown in MMO medium with 20 µM naringenin and MMO medium. (D) Expressional difference between free-living cells grown in MMO medium and TY medium.
|
The nutritional environment also drastically affected the gene expression of A. caulinodans (Table 2 and Fig. 1C). It was noted that 20% of the genes were differentially expressed in free-living cells grown in minimal and rich media; 466 genes were highly expressed, while 477 genes were expressed to a lesser extent in minimal medium than in rich medium. In minimal medium, the genes of the gene clusters involved in exopolysaccharide (EPS) production (AZC_3319 to AZC_3332) and organosulfur metabolism (AZC_3157 to AZC_3180) were expressed to a greater extent than in free-living cells grown in rich medium. The expression levels of these genes were also higher in bacteroids than in free-living cells grown in rich medium (Fig. 1A). On the other hand, a large gene cluster of flagellar biosynthesis genes (AZC_0619 to AZC_0666), a gene cluster of the gene transfer agent (GTA) genes (AZC_1104 to AZC_1127), and two clusters containing branched-chain amino acid transporter genes (AZC_4072 to AZC_4080 and AZC_4583 to AZC_4589) were expressed to a lesser extent in minimal medium (Fig. 1D). The expression levels of these gene clusters were also lower in bacteroids than in free-living cells grown in rich medium (Fig. 1A).
The addition of naringenin influenced the expression of 32 genes; 18 genes were highly expressed and 14 genes were expressed to a lesser extent in the cells supplied with naringenin (Table 2 and Fig. 1D). Of the 18 genes with increased expression, 13 were in the symbiosis island. Of the 14 genes with decreased expression, 7 were in the gene cluster related to organosulfur acquisition, and 5 were in the gene cluster related to EPS synthesis. Of the two remaining genes, one was a gene putatively encoding MarR-like transcriptional factor (AZC_1832) and the other was a gene putatively encoding molybdopterin oxidoreductase (AZC_1874). In the case of Bradyrhizobium japonicum USDA110, about 100 genes were induced in response to genistein, a flavonoid secreted by soybeans (42). Thus, it can be assumed that the regulatory mechanism responding to a flavonoid might be simpler in A. caulinodans than in B. japonicum.
The previous screening of Tn5 mutants showed that 86 genes, including nif, fix, and nod genes, were required for A. caulinodans to establish symbiosis (68). The expression profiling of these genes is listed in the supplemental material (see Table S2). Among these genes, 17 genes, such as nod and noe genes (AZC_3809 and AZC_3811), nif and fix genes (AZC_1038, AZC_1040, AZC_1049, AZC_3414, AZC_3446, AZC_3448, and AZC_3449), a putative transcriptional factor gene (phrR; AZC_0013), a trehalose-6-phosphate phosphatase gene (otsA; AZC_0483), a Lon protease gene (AZC_1610), a C4 dicarboxylate transporter gene (dctA; AZC_3014), EPS synthesis genes (AZC_3319 and AZC_3322), a cysteine desulfurase gene (sufS; AZC_3616), and a fatty acid coenzyme A (CoA) ligase gene (AZC_4536), showed higher expression in bacteroids than in free-living cells grown in both types of media. The expressions of 42 genes, including genes involved in lipopolysaccharide biosynthesis (rfaF [AZC_2251] and rfaD [AZC_2253]), were not significantly different. Five genes (AZC_1714, AZC_1894, AZC_2209, AZC_2370, and AZC_3038) showed lower expression in bacteroids than in free-living cells. Five genes (AZC_0859, AZC_1796, AZC_2131, AZC_2759, and AZC_3904) were called absent or marginal in all replicates applied in this study (data not shown), suggesting that these genes showed little or no expression under all experimental conditions. The nodule-defective phenotypes of the mutants of these 10 genes might have developed due to the polar effect of the Tn5 transposon or unlinked second mutations.
Hereafter, we discuss important features of the genome and the expression profiling of A. caulinodans. WAD values for the genes mentioned below are included in Table S3 in the supplemental material.
Symbiosis island.
As previously described (43), the symbiosis island of A. caulinodans is flanked by tRNA-Gly and interspersed by multiple transposases and integrases, which separate the symbiosis island into several clusters (see Fig. S1 in the supplemental material). Three of these clusters are nod loci (nodD, AZC_3792; the nodABCSUIJZ-noeCHO operon, AZC_3818 to AZC_3808; and nolK, AZC_3850). Two of the clusters (AZC_3826 to AZC_3844 and AZC_3856 to AZC_3877) contain conjugation-related genes. The cluster AZC_3856 to AZC_3877 also contains genes related to amino acid transport. Moreover, a cluster (AZC_3794 to AZC_3800) consists of genes related to nitrile hydratase (NHase).
The expression levels of the nodABCSUIJZ-noeCHO and nolK genes in free-living cells supplied with naringenin were much higher than in those without naringenin (WAD, 0.70 to 3.64) (Fig. 1D). Furthermore, the expression levels of the nodABCSZ-noeHO and nolK genes were higher in bacteroids than in free-living cells grown in both rich and minimal medium (WAD values were 0.73 to 1.25 in a comparison of expression levels in bacteroids and free-living cells grown in minimal medium). The expression of the genes related to conjugation and amino acid transport located in the symbiosis island did not show significant change under any experimental conditions in this study.
A. caulinodans harbors 45 transposase genes on the genome, 19 of which are located in the symbiosis island. Nine of the 45 transposase genes were highly expressed in bacteroids, compared to the levels in free-living cells grown in minimal medium, and 7 of these 9 genes were located in the symbiosis island (AZC_3806, AZC_3821, AZC_3824, AZC_3825, AZC_3846, AZC_3848, and AZC_3851; WAD, 0.63 to 3.86). These results suggest that the symbiosis island is the region where gene transfer has most frequently occurred. Although the high expression of transposase genes in bacteroids has been reported in the transcriptomic analysis of other rhizobia, the biological significance is unknown (9, 71). Interestingly, a transposase gene (AZC_3851) located in the symbiosis island was highly expressed in free-living cells supplied with naringenin. This result suggests that the transposase encoded by AZC_3851 could be activated in the rhizosphere. Such transcriptional behavior of transposase genes has not been reported in other transcriptomic research of rhizobia.
The gene cluster AZC_3794 to AZC_3800 contains the genes AZC_3797 and AZC_3798 encoding proteins with 53% and 30% identity, respectively, to the NHase
and β subunits of Agrobacterium tumefaciens. The expression levels of the genes within this cluster were higher in bacteroids than in free-living cells grown in both rich and minimal media (WAD values were 0.67 to 1.37 in a comparison of expression levels in bacteroids and free-living cells grown in rich medium, and WAD values were 0.56 to 1.62 in a comparison of expression levels in bacteroids and free-living cells grown in minimal medium). In A. tumefaciens and some rhizobia, it is known that NHases are able to convert indole-3-acetonitrile to indole-3-acetamide, and indole-3-acetamide is converted to indole-3-acetic acid by amidase (40). Accordingly, there is a possibility that auxin production may be enhanced in bacteroids.
nif, fix, and suf gene clusters.
Four gene clusters, including nif and fix genes, were found in the genome of A. caulinodans (see Fig. S2 in the supplemental material). The first cluster contained nif genes, along with sufE and sufBCDS genes (AZC_1042 and AZC_1044 to AZC_1047), which are involved in iron-sulfur (Fe-S) cluster assembly, and phn genes, which are involved in phosphonate metabolism. The second cluster contained nif and fix genes, along with phn genes. The first and second gene clusters were conserved among A. caulinodans, X. autotrophicus, and photosynthetic bradyrhizobia, although the syntenies were different (see Fig. S3 in the supplemental material). Most of the nif and fix genes within the first, second, and third clusters were strongly expressed in bacteroids compared to the expression levels in free-living cells (WAD, 0.79 to 5.01 in rich medium and 1.99 to 5.61 in minimal medium) (Fig. 1A and B).
The sufBCDS and sufE genes in the first cluster were strongly expressed in bacteroids (WAD, 1.48 to 3.70 in rich medium and 1.93 to 3.86 in minimal medium). In addition to these suf genes, A. caulinodans has other suf genes outside this cluster, namely, sufSDBCS (AZC_3612 to AZC_3616). The expression levels of these suf genes were also higher in bacteroids than in free-living cells (WAD, 0.90 to 1.82 in rich medium and 0.53 to 1.37 in minimal medium). SufS is cysteine desulfurase, similar to NifS, and supplies sulfur for Fe-S centers (47, 52). The SufBCD complex activates the cysteine desulfurase activity of SufS in conjunction with the SufE that is a sulfur acceptor protein (51). The high expression of these suf genes suggests that bacteroids require a large amount of Fe-S clusters and that Fe-S clusters provided by SufS proteins may be used to sustain nitrogen fixation as well as the NifS protein. In the previous screening, a Tn5 mutant (Ao15-F04) with disruption in the sufS gene (AZC_3616) was obtained, and this mutant formed small stem nodules lacking nitrogen-fixing ability (68). This result also supports the importance of suf genes in stem nodule formation.
Ammonium assimilation.
Unexpectedly, the expression levels of a putative glnK-amtB operon (AZC_4227 to AZC_4226) encoding a nitrogen regulatory protein PII and an ammonium transporter, respectively, were higher in bacteroids than in free-living cells (WAD values were 2.39 and 2.33 in a comparison of expression levels in bacteroids and free-living cells in rich medium, and WAD values were 2.16 and 1.79 in a comparison of expression levels in bacteroids and free-living cells grown in minimal medium) (Fig. 1A and B). In the free-living state, the fixed nitrogen is assimilated by glutathione synthase (17), whereas during symbiosis, the fixed nitrogen is exported from the bacteroids to the plant cells by passive diffusion and assimilated by the host. It has been reported that the expression levels of genes encoding the PII-like protein of Rhizobium etli and Sinorhizobium meliloti were lower in bacteroids than in free-living cells (1, 69). However, our result contradicted these findings. Furthermore, transcriptomic analysis of B. japonicum showed that the expression level of the glnK-amtB operon was high in nodules (54). Although the biochemical function of glnK-amtB in nodules is unknown, it is possible that the transcriptional regulation of glnK-amtB in A. caulinodans and B. japonicum is different from that of R. etli and S. meliloti.
Carbon uptake and metabolism.
A. caulinodans, as well as other rhizobia, uses C4-dicarboxylates as primary carbon sources and has multiple C4-dicarboxylate transport systems. Prokaryotes have three types of C4-dicarboxylate transport systems, that is, a major facilitator superfamily (MFS)-type transporter, an ABC-type transporter, and a tripartite ATP-independent periplasmic transporter. A. caulinodans has four genes (AZC_1325, AZC_2673, AZC_3014, and AZC_3353) putatively encoding an MFS-type transporter (DctA). Searches performed using BlastP revealed that the protein encoded by AZC_3014 was 81% identical to S. meliloti DctA and that the proteins encoded by the other three genes were about 50% identical to S. meliloti DctA. On the genome of S. meliloti, there is a single dctA gene, and a dctA mutant of S. meliloti forms ineffective nodules lacking nitrogen-fixing ability (23). A. caulinodans has 21 genes putatively encoding any of the DctP, -Q, and -M components of tripartite ATP-independent periplasmic C4-dicarboxylate transporters. Of the multiple C4-dicarboxylate transporter genes of A. caulinodans, only AZC_3014 was highly expressed in bacteroids compared to the expression level in free-living cells (WAD, 1.49 in rich medium and 2.38 in minimal medium). We have isolated a Tn5 mutant (Ao62-E02) with disruption in AZC_3014, and this mutant was observed to form inefficient stem nodules lacking nitrogen-fixing ability (68). These results suggest that the MFS transporter encoded by AZC_3014 is the most important C4-dicarboxylate transporter in bacteroids. Thus, C4-dicarboxylate uptake by bacteroids may depend mainly on this transporter.
In addition to C4-dicarboxylates, the array data suggest the possibility that bacteroids of A. caulinodans are supplied with other carbon compounds. A. caulinodans contains AZC_2923 to AZC_2921 genes putatively encoding acetone carboxylase subunits (AcxA, -B, and -C), which are 89, 88, and 79% identical to AcxA, -B, and -C of X. autotrophicus, respectively. It is known that X. autotrophicus has the ability to use acetone as a carbon source. X. autotrophicus produces acetoacetate by the carboxylation of acetone in an ATP-dependent manner (62, 74). The high similarities of amino acid sequences suggest that A. caulinodans may have the ability to use acetone as a carbon source. The gene AZC_2924 adjacent to the acxABC genes encodes a protein with 55% identity to the AcxR protein of X. autotrophicus, which is a transcriptional activator for acxABC genes (63). In the previous screening, a Tn5 mutant (Ao22-A06) with disruption in AZC_2924 was isolated, and this mutant was deficient in nodule formation. It is interesting that the expression levels of acetone carboxylase genes were higher in bacteroids than in free-living cells (WAD, 1.01 to 1.95 in rich medium and 1.48 to 2.40 in minimal medium) (Fig. 1A and B). Furthermore, it was reported that acetone was produced in soybean nodules (67). Although it is uncertain whether acetone is also produced in the nodules of S. rostrata, if acetone production is a general event in the nodules of legumes, A. caulinodans bacteroids may utilize acetone as a carbon source in addition to C4-dicarboxylates.
A. caulinodans contains genes (AZC_2937 to AZC_2939) putatively encoding carbon monoxide dehydrogenase (CODH) subunits (CoxM, -L, and -S) with 34, 38, and 61% identities, respectively, to the CoxM, -L, and -S of Oligotropha carboxidovorans, which is the most intensively studied representative of the carboxydotrophic bacteria. CODH is an O2-stable molybdenum-iron-sulfur-flavin hydroxylase that catalyzes the oxidation of CO to CO2. Carboxydotrophic bacteria constitute a small but diverse group of aerobes (48), and putative coxLMS genes have been identified in the genomes of bacteria that are not known to be CO oxidizers. Among the rhizobia, the B. japonicum strain 110 spc4 is known to have the ability to grow using CO as a sole carbon and energy source under aerobic conditions (45). The expression levels of putative coxLMS genes of A. caulinodans were higher in bacteroids than in free-living cells (WAD, 0.84 to 1.37 in rich medium and 1.50 to 2.34 in minimal medium) (Fig. 1A and B). If these genes really encode functional CODH, there is a possibility that bacteroids of A. caulinodans may oxidize CO in nodules to use CO as a carbon source.
Malic enzyme (ME) catalyzes the reversible oxidative decarboxylation of malate to pyruvate and CO2. The generation of acetyl-CoA via ME and pyruvate dehydrogenase is important for nitrogen-fixing bacteroids in order to metabolize C4-dicarboxylates via the trichloroacetic acid cycle. S. meliloti contains NAD(P)+- and NADP+-dependent ME encoded by dme and tme, respectively (21, 22). A dme mutant of S. meliloti forms nodules showing no nitrogen-fixing ability, while a tme mutant forms wild-type nodules (21, 22). On the genome of A. caulinodans, two genes (AZC_0119 and AZC_3656) putatively encoding ME were identified by BlastP analyses. The protein encoded by AZC_0119 was similar to plant or animal NADP+-dependent ME (44% identical to NADP-ME 1 of Arabidopsis thaliana). The protein encoded by AZC_3656 was 46% identical to TME of S. meliloti. Two Tn5 mutants (Ao30-B02 and Ao58-F11) with disruption in AZC_3656 formed nodules lacking nitrogen-fixing ability (68), suggesting that this putative ME similar to S. meliloti TME is essential for nodule formation by A. caulinodans, although TME is not essential for nodule formation by S. meliloti. The expression levels of AZC_3656 were not different between bacteroids and free-living cells. In contrast, the expression levels of AZC_0119 were much higher in bacteroids than in free-living cells (WAD, 3.68 in rich medium and 3.84 in minimal medium) (Fig. 1A and B). It is important to determine whether AZC_0119 is essential for nitrogen fixation or not. Therefore, the functional assignment of each ME encoded by AZC_0119 and AZC_3656 should be determined.
Sulfur uptake.
Sulfur is a primary component of sulfur-containing amino acids (cysteine and methionine) and plays an essential role as a substituent in certain enzyme cofactors. In rhizobia-legume symbiosis, sulfur has important roles. In S. meliloti, sulfation of the Nod factor is essential to elicit the nodule formation of host plants (56), and sulfated Nod factors are more resistant to degradation than nonsulfated ones (64). Furthermore, nitrogen fixation requires Fe-S clusters. Bacteria typically utilize inorganic sulfate or cysteine as their primary sulfur source, and they also utilize sulfonates or sulfate esters under sulfur-limited conditions (39). Sulfate uptake mechanisms have been well analyzed in E. coli (32, 59, 60). In E. coli, the ABC-type sulfate/thiosulfate transporters are encoded by the cysPTWA genes, with the unlinked sbp gene located separately on the genome. CysT and CysW are permease components, and CysA is an ATP-binding protein. CysP and Sbp are substrate-binding proteins with high affinity for thiosulfate and sulfate, respectively. However, CysP and Sbp possess an overlapping substrate range and share the CysTWA system. The genes related to sulfonate utilization by bacteria have been characterized in B. subtilis (72), E. coli (73), and Pseudomonas putida (38). E. coli contains two types of sulfonate transport and desulfonation gene clusters, ssuABCDE and tauABCD. SsuABC and TauABC are ABC transporters that import alkanesulfonates and taurine, respectively. Alkanesulfonates and taurine are desulfurized by a two-component reduced flavin mononucleotide-dependent monooxygenase (SsuDE) and an
-ketoglutarate-dependent dioxygenase (TauD), respectively. Sulfate ester transporter is characterized as an ABC-type transporter encoded by atsRBC genes in Pseudomonas aeruginosa (33).
BlastP analyses revealed that A. caulinodans contained two genes (AZC_0302 and AZC_0922) encoding proteins that were 50 to 72% identical to E. coli CysP/Sbp. The genes (AZC_0921 to AZC_0919) downstream of AZC_0922 encoded proteins that were 48%, 52%, and 57% identical to E. coli CysT, -W, and -A, respectively. It is noteworthy that a large gene cluster (AZC_3157 to AZC_3180) containing putative ssu genes (AZC_3157 and AZC_3160 to AZC_3164) and ats genes (AZC_3168 to AZC_3171, AZC_3173, and AZC_3177) was identified on the genome. Each protein encoded by these genes was 31 to 61% identical to the corresponding proteins of E. coli or P. aeruginosa (see Table S4 in the supplemental material). AZC_3159 located next to the putative ssu genes (AZC_3160 to AZC_3164) was previously annotated as actinorhodin polyketide dimerase. However, an InterProScan search revealed that the protein encoded by AZC_3159 contained a flavin mononucleotide-binding domain (data not shown), and a BlastP search revealed that this protein was 29% identical to E. coli RutF, a putative flavin, i.e., NADH component of oxidoreductase (see Table S4 in the supplemental material). From these results, it appears that the protein encoded by AZC_3159 might act as an NAD(P)H-flavin oxidoreductase, like SsuE, in cooperation with a putative SsuD encoded by AZC_3157.
The expression levels of putative cysP and sbp genes (AZC_0302 and AZC_0922) were higher in bacteroids (WAD, 0.54 and 1.09, respectively) and free-living cells grown in minimal medium (WAD, 1.04 and 1.61, respectively) than in free-living cells grown in rich medium. The expression levels of 12 genes (AZC_3157 to AZC_3161, AZC_3163, AZC_3165, AZC_3172 to AZC_3174, and AZC_3179 to AZC_3180) within the AZC_3157 to AZC_3180 cluster, including the ssu and ats genes, were also higher in bacteroids than in free-living cells grown in rich medium (WAD, 0.55 to 2.17). The expression levels of 14 genes (AZC_3157 to AZC_3161, AZC_3163, AZC_3167, AZC_3169 to AZC_3174, and AZC_3177) in free-living cells were higher in minimal medium than in rich medium (WAD, 0.57 to 2.25) (Fig. 1A and C). Furthermore, a Tn5 mutant (Ao61-B09) with disruption in AZC_3159 formed stem nodules lacking nitrogen-fixing ability (68). It is unclear whether these genes are related to sulfur acquisition; however, if they are, these data suggest that bacteroids as well as free-living cells grown in minimal medium might actively acquire organosulfur compounds in addition to inorganic sulfate/thiosulfate. In addition, organosulfur uptake may have advantages in obtaining carbon sources.
Vitamins.
Biotin is an essential cofactor in carboxylation, decarboxylation, and transcarboxylation reactions. Generally, biotin is synthesized via the conserved four enzymes encoded by bioF, bioA, bioD, and bioB starting from pimeloyl-CoA (see Fig. S4 in the supplemental material), although pimeloyl-CoA is synthesized in various ways (34, 57, 65). Bacteria have the ability not only to synthesize biotin but also to acquire biotin from the environment. It was proposed that bioYMN genes consisting of an operon are involved in biotin uptake in R. etli (30), and recently, it was found that BioY is a biotin transporter and that BioN and BioM are dedicated to being energy-coupling modules (31). The bioY genes have been identified on a variety of bacterial genomes, and bioY genes are not always adjacent to bioMN genes. A. caulinodans is biotin auxotrophic, and BlastP analyses revealed that A. caulinodans had a gene (AZC_0361) encoding a protein 57% identical to E. coli BioB, but no bioF, bioA, or bioD genes were identified by BlastP analyses. It was also revealed that A. caulinodans contained three genes (AZC_0362, AZC_2177, and AZC_3757) encoding proteins 34%, 31%, and 52% identical to R. etli BioY, respectively (see Fig. S4 in the supplemental material). AZC_0362 is located next to the putative bioB gene (AZC_0361), AZC_2177 is an orphan, and AZC_3757 is located next to bioMN genes (AZC_3755 and AZC_3756). The expression levels of AZC_0361 to AZC_0362 (bioBY cluster) were higher in bacteroids than in free-living cells (WAD values were 3.08 and 3.47 in a comparison of expression levels in bacteroids and free-living cells grown in rich medium, and WAD values were 1.70 and 1.96 in a comparison of expression levels in bacteroids and free-living cells grown in minimal medium), but those of AZC_2177 (orphan bioY) and AZC_3755 to AZC_3757 (bioMNY cluster) were not. This result suggests that the biotin uptake of bacteroids depends mainly on the biotin transporter encoded by AZC_0361, if this protein is really a functional biotin transporter. Why was the putative bioB gene highly expressed in bacteroids, even though A. caulinodans has only this gene among biotin synthesis genes? It is speculated that the BioY-like protein encoded by AZC_0362 may act as a dethiobiotin transporter, and dethiobiotin obtained by AZC_0362 may be converted to biotin by BioB encoded by AZC_0361.
Thiamine is an essential cofactor for several important enzymes of carbohydrate metabolism, including the trichloroacetic acid cycle. Thiamine synthesis and the related genes of bacteria are well studied in E. coli and B. subtilis, as reviewed in references 35 and 58. BlastP analyses revealed that A. caulinodans contained genes homologous to the thiamine-related genes of E. coli and B. subtilis (31 to 65% similar, based on the amino acid sequences [see Fig. S5 in the supplemental material]). The expression of putative thiCGSO genes (AZC_2607 and AZC_3552 to AZC_3554) was higher in bacteroids (WAD, 0.63 to 1.22) and in cells grown in minimal medium (WAD, 1.01 to 2.52) than in those grown in rich medium. We have isolated three Tn5 mutants (Ao5-E01, Ao55-D02, and Ao73-E10) with disruption in the putative thiC gene, and these mutants form stem nodules showing lower nitrogen-fixing activity than the wild-type strain (68). However, the nitrogen-fixing activities of these mutants grown in minimal medium were not significantly different from that of the wild-type strain (68). From these results, it can be deduced that highly activated de novo synthesis of thiamine is essential for bacteroids but not for free-living cells grown in minimal medium, although thiCGSO genes were highly expressed under both conditions.
EPS synthesis.
Rhizobial EPS is one of the important factors in the establishment of symbiosis. Rhizobial EPSs are linear or branched heteropolysaccharides that consist of common sugars (61). A. caulinodans produces a linear homopolysaccharide of
-1,3-linked 4,6-O-(1-carboxyethylidene)-D-galactopyranoside residues (15). S. meliloti produces two structurally distinct EPSs: EPS I (succinoglycan) and EPS II (galactoglycan). These EPSs are essential for the establishment of nitrogen-fixing symbiosis (3). S. meliloti contains two gene clusters (exo-exs and exp gene clusters), directing the biosynthesis of EPS I and EPS II, respectively (3).
BlastP searches suggested that no gene cluster corresponding to the S. meliloti exo-exs cluster existed on the genome of A. caulinodans (data not shown). On the contrary, two gene clusters putatively corresponding to the S. meliloti exp gene cluster, which is composed of five putative operons designated as expA (nine genes), expC (one gene), expD (two genes), expE (eight genes), and expG (one gene) (4, 25), were identified on the genome of A. caulinodans (see Table S5 in the supplemental material). One is composed of AZC_1831 to AZC_1834 encoding proteins that are 45 to 61% identical to ExpA7, -10, -9, and -8 of S. meliloti, respectively. The other is AZC_3319 to AZC_3332. AZC_3319 to AZC_3331 encodes proteins that are 21 to 64% identical to ExpA5, ExpE1 to -7, ExpD1 and -2, ExpG, ExpC, and ExpA4 of S. meliloti, respectively. AZC_3332 encodes a hypothetical protein. The former and latter clusters are designated as exp cluster I and II, respectively. Genes homologous to expA1, expA23, and expE8 of S. meliloti are not contained in these clusters.
There were no significant changes in the expression levels of the exp cluster I genes under any experimental conditions in this study. The expression levels of the exp cluster II genes were higher in free-living cells grown in minimal medium than in those grown in rich medium (WAD, 0.81 to 3.69) (Fig. 1C). In S. meliloti, the expression of the exo gene cluster increased in minimal medium, although the expression of the exp gene cluster did not (2), suggesting that the regulatory system of the exp genes of A. caulinodans is different from that of S. meliloti. Furthermore, the expression levels of the exp cluster II genes of A. caulinodans were higher in bacteroids than in free-living cells (WAD, 1.63 to 4.43 in rich medium and 0.58 to 2.73 in minimal medium) (Fig. 1A and B). In contrast, the expression of most exo genes of S. meliloti did not increase in bacteroids, and the expression of exoX, a posttranscriptional inhibitor of EPS synthesis, increased in bacteroids (2). These results suggest that EPS might be actively synthesized in bacteroids of A. caulinodans but not in bacteroids of S. meliloti. We have isolated two Tn5 mutants (Ao3-G08 and Ao11-A12) with disruption in the putative expA5 (AZC_3319) and expE5 (AZC_3322) genes, respectively (68). Ao3-G08 formed very small nodules called bumps, and Ao11-A12 formed small nodules showing nitrogen-fixing activity. Future studies on the EPS synthesis of these mutants will help to clarify the function of EPS in nodulation formation by A. caulinodans.
GTA genes.
It is noteworthy that a gene cluster related to an unusual class of a virus-like genetic exchange element, the GTA, which was discovered in Rhodobacter capsulatus (46), was identified on the genome of A. caulinodans. BlastP analyses revealed that 15 genes in the gene cluster of AZC_1104 to AZC_1127 encoded proteins 28 to 55% identical to the GTA proteins of R. capsulatus (see Fig. S6 in the supplemental material). The structure of GTA is similar to a small, tailed phage (77), but GTA seems to function only in the transfer of small, random pieces of genomic DNA between cells (41). Recent genome sequencing projects have revealed that GTA gene clusters are present exclusively in alphaproteobacteria (7, 41). More than half of the sequenced alphaproteobacteria genomes contain a GTA-like gene cluster, including Rhizobiales. However, among the symbiotic rhizobia, only A. caulinodans contains the GTA-like gene cluster.
The expression levels of the putative GTA genes of A. caulinodans were lower in bacteroids (WAD, –0.56 to –1.78) and free-living cells grown in minimal medium (WAD, –0.55 to –1.47) than in those grown in rich medium (Fig. 1A and C). As mentioned above, among the rhizobia, GTA genes were found only in A. caulinodans, and the expression of these genes was suppressed in bacteroids. From these results, GTA may have adverse effects on the establishment of symbiosis—or GTA may simply be unnecessary for symbiosis. To investigate which hypotheses are true, it will be necessary to introduce GTA genes into other rhizobia or to enhance artificially the expression of GTA genes in A. caulinodans.
Concluding remarks.
So far, the genome sequencing of A. caulinodans and the screening of genetic factors required for stem nodule formation have been carried out (43, 68). In the present study, we conducted the transcriptional profiling of A. caulinodans in free-living and symbiotic states. The data obtained from these studies enabled us to ascertain some aspects of the mechanism of stem nodule formation by A. caulinodans. Many of the genes we have mentioned in the present report, such as nod, nif, fix, hup, and hyp, were previously reported to be involved in nodule formation; however, transcriptional profiling also revealed that A. caulinodans adapts to the endosymbiotic state by uptaking a variety of compounds. The expression patterns of ssu genes and acxABC genes, especially, indicate that A. caulinodans uses a wide variety of compounds for the acquisition of sulfur and carbon sources. The active acquisition and metabolism of sulfur and carbon sources should reflect the higher ability of nitrogen fixation of stem nodules formed by A. caulinodans. Furthermore, the expression patterns of genes involved in EPS biosynthesis indicate that EPS has an extra function in stem nodules. In addition, A. caulinodans has some distinctive features among the rhizobia: the ability to induce stem nodules after crack-entry infection; the ability to fix nitrogen both in the symbiotic and free-living states even at relatively high concentrations of dissolved oxygen (up to 12 µM) (12); and the ability to colonize the xylem of wheat, rice, tomato, and Arabidopsis plants (11, 28, 29, 66, 75, 76). These traits of A. caulinodans make this bacterium one of the most qualified candidates for the application of rhizobia to form nitrogen-fixing nodules on plants other than legumes. The findings of the present study and further studies will accelerate not only the elucidation of the symbiotic mechanism of legumes and rhizobia but also their applications in agriculture in the future.
Published ahead of print on 19 June 2009. ![]()
Supplemental material for this article may be found at http://aem.asm.org/. ![]()
These authors contributed equally to the work. ![]()
Present address: Center for Information Biology and DNA Data Bank of Japan, National Institute of Genetics, Yata 1111, Mishima, Shizuoka 411-8540, Japan. ![]()
¶ Present address: Institute of Biotechnology, National Taiwan University, R412, No. 81, Chang-Xing St., Taipei 106, Taiwan (R.O.C.). ![]()
|
|
|---|
ucha, A. Böck, and D. Hulanicka. 1990. Sulfate and thiosulfate transport in Escherichia coli K-12: identification of a gene encoding a novel protein involved in thiosulfate binding. J. Bacteriol. 172:3358-3366.
-proteobacteria. Trends Microbiol. 15:54-62.[CrossRef][Medline]
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2010 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»