Skip to main content

The Mesorhizobium huakuii transcriptional regulator AbiEi plays a critical role in nodulation and is important for bacterial stress response



Bacterial abortive infection (Abi) systems are type IV toxin–antitoxin (TA) system, which could elicit programmed cell death and constitute a native survival strategy of pathogenic bacteria under various stress conditions. However, no rhizobial AbiE family TA system has been reported so far. Here, a M. huakuii AbiE TA system was identified and characterized.


A mutation in M. huakuii abiEi gene, encoding an adjacent GntR-type transcriptional regulator, was generated by homologous recombination. The abiEi mutant strain grew less well in rich TY medium, and displayed increased antioxidative capacity and enhanced gentamicin resistance, indicating the abiEi operon was negatively regulated by the antitoxin AbiEi in response to the oxidative stress and a particular antibiotic. The mRNA expression of abiEi gene was significantly up-regulated during Astragalus sinicus nodule development. The abiEi mutant was severely impaired in its competitive ability in rhizosphere colonization, and was defective in nodulation with 97% reduction in nitrogen-fixing capacity. The mutant infected nodule cells contained vacuolation and a small number of abnormal bacteroids with senescence character. RNA-seq experiment revealed it had 5 up-regulated and 111 down-regulated genes relative to wild type. Of these down-regulated genes, 21 are related to symbiosis nitrogen fixation and nitrogen mechanism, 16 are involved in the electron transport chain and antioxidant responses, and 12 belong to type VI secretion system (T6SS).


M. huakuii AbiEi behaves as a key transcriptional regulator mediating root nodule symbiosis.

Peer Review reports


Toxin-antitoxin (TA) systems are ubiquitous constituents found on plasmids or chromosomes of countless bacteria, archaea, and possibly also unicellular fungi [1]. Bacterial TA systems are composed of bicistronic operons encoding a stable toxin that can harm the host cell and its cognate labile antitoxin, which protects the host from the toxin’s deleterious effect [2]. The product of the toxin gene is a protein, while the product of the antitoxin is either a protein or a non-coding RNA [3]. TA systems have been shown to play various physiological roles in the formation of dormancy and persister cells, survival during infection, adaptation to hostile environments, programmed cell death and biofilm formation [4,5,6,7]. Based on the biochemical nature and mode of action of the antitoxin gene product, bacterial TA systems have been divided into six types: small regulatory RNAs in types I and III, and antitoxins of proteinaceous nature in type II, IV, V and VI TA systems [8, 9].

The proteinaceous antitoxin of type IV system neutralizes its cognate toxin by forming toxin–antitoxin complexes instead of a direct protein-protein interaction [10]. Several common families of type IV system have been identified on the chromosomes of bacteria and archaea: CbtA/CbeA, YkfI/YafW, YpjF/YfjZ, and AbiEii/AbiEi. The pair CbtA/CbeA was the first type IV system found in E. coli. The toxin CbtA alters cell shape by inhibiting the polymerization of cytoskeletal proteins FtsZ and MreB through direct interaction, without showing direct interaction with its cognate antitoxin [3]. The first AbiE system was reported in Lactococcus lactis where it acts at the post-transcriptional level in the lytic cycle to abort phage development [11]. AbiE was found to be a type IV TA system as there is no apparently direct interaction between the antitoxin and toxin, and the antitoxin antagonizes toxin activity by stabilizing its targets [12]. The AbiEii toxin is a putative nucleotidyltransferase containing a C-terminal domain involved in toxin neutralization. The AbiEi antitoxin, a transcriptional regulator, contains an N-terminal domain required for repression of abiE transcription, and a bi-functional C-terminal domain required for transcriptional repression and sufficient for toxin neutralization [13]. Previous work on the Streptococcus agalactiae AbiE system revealed that AbiEi negatively autoregulates the abiE promoter and is sufficient for antitoxicity, yet it is unknown whether all systems provide dual resistance functions [3, 13].

Bacteria of the family Rhizobiaceae are able to induce symbiotic nodules on the roots of leguminous plants where bacteroids convert atmospheric nitrogen to ammonia [14]. This highly complex process requires a specific signal exchanges between the partners. Compatible rhizobia sense the flavonoids released from their host roots through the transcriptional activator NodD, which trigger nodule organogenesis [15]. Earlier reports described the presence of type II TA modules in Sinorhizobium meliloti and Bradyrhizobium japonicum. The bat/bto module is classified as a type II TA system belonging to the vapBC-family, and deletion of the B. japonicum bat/bto operon resulted in alterations of several metabolic pathways and defective symbiotic performance due to the changes in lipopolysaccharide (LPS) [16]. Milunovic et al. reported that loss of 4 S. meliloti type II TA systems results in growth inhibition, but does affect the S. meliloti-alfalfa symbiosis [17]. However, the S. meliloti ntrR gene, a member of the type II TA system, was controlled negatively by its own product and positively by the symbiotic regulator syrM. Mutation in the ntrR gene induced nodules with enhanced nitrogen fixation capacity [18].

Rhizobial genomes frequently contain type IV AbiE toxin-antitoxin operator, but unlike with the type II TA system, knowledge of regulation of the type IV AbiE system in Rhizobium species is still poorly documented. Here, we identified a type IV TA antitoxin gene abiEi in M. huakuii 7653R, and the roles of M. huakuii abiEi in free-living bacteria and during N2-fixing symbiosis with A. sinicus were investigated by analyzing the phenotypes of abiEi mutant strain. A transcriptome analysis was also carried out to discover clues that might explain the differences in the nodules induced by the abiEi mutant and the wild type strain. To our knowledge, this work represents the first transcriptome analysis of TA gene in symbiotic root nodules reported to date.


Determination of AbiE system activity

In order to elucidate the function of the TA system, the AbiE gene cassette and the toxin gene abiEii was overexpressed in E. coli, and the effect of expressed protein product was ruled out by comparing the growth of E. coli BL21(DE3) cells in the presence versus the absence of the plasmid (Fig. 1). Under IPTG induction, both E. coli strains harboring either recombination vector pETAbiE or pETAbiEii grew more slowly than control cells with the empty vector pET-28a, which indicated that induced expression of AbiE system or AbiEii toxin protein has a negative effect on the growth of host cells. However, E. coli strains harboring both toxin and antitoxin genes (pETAbiE) grew slightly, but significantly faster than that harboring toxin gene abiEii (pETAbiEii), suggesting antitoxin AbiEi is a transcriptional regulator that can repress abiEii expression to influence toxin.

Fig. 1

The growth curve of the E.coli BL21 strains after IPTG induction. E. coli strain BL21(pET-28a), the recombinant strains BL21 (pETAbiE) and BL21 (pETAbiEii) were cultured in LB at OD600 of 0.4, then induced with 1 mM IPTG. Data are from three biological samples plus and minus the standard error of the mean (± SEM)

Construction of M. huakuii abiEi mutant

The gene MCHK_RS33180 in M. huakuii encodes a putative type IV TA system transcriptional regulator AbiEi. The abiEi gene is predicted to encode a 221-amino acid polypeptide with an expected molecular mass of 24.94 kDa and a pI value of 9.95. To confirm the function of the abiEi gene in growth, environment stress and symbiotic nitrogen fixation, a mutant HKabiEi strain of this gene was constructed by single crossover homologous recombination. qRT-PCR was firstly employed to examine the relative mRNA levels of antitoxin and toxin in free-living cells from M. huakuii. The expression of antitoxin gene abiEi was almost not detected in the mutant HKabiEi, while the expression of toxin gene abiEii was 0.86 ± 0.11, and there was no significant difference in the mutant HKabiEi compared to the wild-type strain 7653R. The results confirmed that abiEi was disrupted by the insertion of vector pK19mob in HKabiEi, and also indicated that antitoxin AbiEi didnot directly regulate toxin AbiEi expression, and may antagonize toxin activity by stabilizing its targets. The mutant showed no significant difference in AMS minimal medium with NH4Cl as nitrogen source and glucose as carbon source, but exhibited defective growth in rich TY medium (Fig. 2).

Fig. 2

Growth of M. huakuii in AMS or TY media. Wild-type 7653R and mutant HKabiEi were grown in AMS supplemented with Glc/NH4+ (A) and TY (B). Data are from three biological samples plus and minus the standard error of the mean (± SEM)

Role of the abiEi gene in the regulation of antibiotic resistance

It has been reported that TA systems constitute a native survival strategy of pathogenic bacteria and thus are potential targets of antibiotic drugs [19]. To examine the possibility that abiEi affects the antibiotic-resistance of M. huakuii, sensitivity to gentamicin (Gm) and chloramphenicol (Cm) at low concentration was assayed in AMS minimal medium. In the presence of 2 μg mL− 1 Gm, the mutant HKabiEi was grown to the early logarithmic phase at 24 h, and entered the stationary at 48 h postinoculation, whereas the parent strain 7653R was grown to the early logarithmic phase at 48 h, and entered the stationary phase at 72 h postinoculation (Fig. 3A). The resistance of the M. huakuii strain to 2 μg mL− 1 Cm stress was also estimated, the abiEi mutant HKabiEi showed no difference in growth as compared with the wild-type (Fig. 3B). These data indicated that mutantion in abiEi gene can display different antibiotic susceptibilities of M. huakuii.

Fig. 3

The effect of antibiotics on bacterial growth. Wild-type 7653R and mutant HKabiEi were grown in AMS Glc/NH4+ supplemented with 2 μg ml−1 Gm (A) or Cm (B). Data are from three biological samples plus and minus the standard error of the mean (± SEM)

Role of the abiEi gene in the regulation of redox homeostasis

It has been reported that the bacterial TA system is a module that plays a role in cell survival under stress conditions [20]. To determine the function of the M. huakuii abiEi gene in the regulation of redox homeostasis, the sensitivity of the abiEi mutant strain to the inhibition of growth by H2O2 or SNP, which could be due to changes in the cellular redox status [21], was examined. The zone of inhibition induced by H2O2 for mutant HKabiEi significantly decreased as compared with the wild-type 7653R strain (Table 1, Additional file 1). When abiEi on plasmid pBBR1MCS-5 was introduced into mutant HKabiEi, the resulting strain HKabiEi(pBBRabiEi) could rescue the variation of the inhibition zone, and showed hypersensitive to H2O2 as it made significantly bigger (p < 0.05) inhibition zone than the wild-type 7653R strain, whereas HKabiEi harboring an empty plasmid showed no significant difference with regard to mutant strain (Table 1, Additional file 1). The resistant to SNP stress was also estimated, and mutant HKabiEi was more resistant to 10 μg mL− 1 SNP compared with the wild-type 7653R strain (Fig. 4). These results indicated that abiEi gene plays a negative regulatory role in the oxidative stress response.

Table 1 Tolerance of strains to different concentrations of H2O2
Fig. 4

The effect of SNP on bacterial growth. Wild-type 7653R and mutant HKabiEi were grown in AMS Glc/NH4+ supplemented with 10 μg ml− 1 SNP. Data are from three biological samples plus and minus the standard error of the mean (± SEM)

Effect of mutation of abiEi on rhizosphere competition

Competition between the abiEi mutant HKabiEi and parent strain 7653R for growth in the plant rhizosphere was measured by inoculating a low number of M. huakuii strains into the A. sinicus rhizosphere (103 to 104 CFU per seedling), and determining the total amount of bacteria after 7 days. When the mutant HKabiEi and the parent 7653R were inoculated alone into short-term colonization of sterile plant rhizosphere, the population density of the mutant strain was a little (but significantly) higher than that of wild-type (Fig. 5). When these strains were inoculated in equal proportions, mutant HKabiEi was at a significant disadvantage (32.90% of bacteria recovered). When strain HKabiEi was inoculated at a 10-fold excess over wild-type 7653R, it still accounted for only 48.68% of bacteria recovered (t-test; P ≤ 0.01). TA systems are ubiquitous bacterial systems that may function in metabolic stress management, but are also thought to play a role in virulence by helping pathogens survive stress [22]. The decreased ability of the abiEi mutant to complete in host plant rhizosphere shows that transcriptional regulator AbiEi is essential for colonization of the host plant rhizosphere by M. huakuii.

Fig. 5

Bacteria recovered (7 dpi) from the rhizosphere of Astragalus sinicus plants following inoculation with wild-type 7653R and HKabiEi, both individually and together. Inoculation ratios are given on the x axis, with 1 corresponding to 103 CFU. Number of bacteria (per plant) recovered from at least 5 plants (mean ± SEM) are shown. a,b Values with different letters are signifificantly different between mutant and wild-type control (two-way ANOVA, P < 0.05)

Effect of mutation of abiEi on nodulation

In order to assess the nodulation and nitrogen fixing capacity of the abiEi mutant, A. sinicus seedlings were inoculated with M. huakuii strains. Four weeks post-inoculation, the number, structure, and acetylene reduction activity (ARA) values of the nodules were analyzed (Table 2 and Fig. 6). No statistically significant difference in the number of nodules per plant was observed between plants inoculated with strain HKabiEi and plants inoculated with wild-type 7653R (Table 2). In contrast, the abiEi mutant HKabiEi elicited small white nodules on A. sinicus, and the plants inoculated with the strain HKabiEi were similar to that of plants without inoculation. A notable feature of our study was that the nitrogen fixation capacity was severely affected in the abiEi mutant HKabiEi, with a reduction from 35.20 nmol of C2H4 plant− 1 h− 1 in 7653R to 0.90 nmol of C2H4 plant− 1 h− 1 in the abiEi mutant (Table 2). Plants nodulated with the complemented strain HKabiEi(pBBRabiEi) had somewhat wild-type properties; approximately the same plant fresh weight, red nodules and reduce acetylene per nodule at the same rate as 7653R-inoculated.

Table 2 Symbiotic phenotype of 7653R and HKabiEi α,β
Fig. 6

Structure of 4-week-old Astragalus sinicus nodules and bacteroids. Nodules were induced by M. huakuii 7653R (A, D), HKabiEi (B, E), HKabiEi(pBBRabiEi) (C, F). Scale bars = 200 μm (A, B, C), 2 μm (D, E, F). i, incrassated membrane; dissociation; v, vacuolation; s, senescing bacteroid

To further investigate the symbiotic role of abiEi in more detail, sections of mature nodules were prepared and examined by an Olympus microscope and scanning electron microscopy (SEM). Microscopic analysis of HKabiEi nodules showed that they were small and spherical instead of large and elongated like the wild-type, and compared with wild-type nodules, relatively few nodule cells were infected in the mutant nodules (Fig. 6). SEM analysis demonstrated an aberrant shape of abiEi bacteroids, which were reduced compared to wild-type bacteroids (Fig. 6). A quantitative analysis of bacteroid lengths demonstrated that the mean length of wild-type bacteroids was more than three times those of the abiEi mutant bacteroids, but not different from those of the complemented strains HKabiEi(pBBRabiEi). Moreover, the mutant infected nodule cells contained vacuolation and a small number of abnormal bacteroids, with incrassated membrane, scarce content in bacteroids and signs of premature senescence of endosymbiotic bacteria. The results suggested that the mutant-induced nodules were functionally defective.

Activation of abiEi gene in 7653R-inoculated nodules

As toxin–antitoxin transcriptional regulator AbiEi plays a pivotal role in symbiotic nitrogen fixation, the expression of abiEi gene in root nodules was analyzed by quantitative RT-PCR (qRT-PCR). The abiEi gene expression was significantly up-regulated in the early stage of nodule formation (14 d), the nodule maturation stage (28 d) and the late stage (42 d) of nodule development and senescence, and the abiEi gene had the highest expression level (more than 5-fold) in nodules at 28 days post inoculation (Fig. 7). Therefore, abiEi gene expression was induced during the symbiotic interaction when compared with free-living cells growing in synthetic medium, and abiEi may play an important role in persistence of nodule bacteroids and prevention of premature nodule senescence.

Fig. 7

Expression patterns of abiEi gene in symbiotic nodules. Gene expression levels were examined by real-time RT-PCR. Nodules were collected on different days after inoculation with M. huakuii 7653R. Relative expression of abiEi gene involved in nodule bacteroids compared with 7653R cells growth in AMS medium. Data are the average of three independent biological samples (each with three technical replicates). Asterisk (*) indicates a significant difference (FC > 2, P < 0.05)

RNA-seq analyses of gene expression in the nodule bacteroids

RNA-Seq was performed to analyze the overall effect of abiEi on the transcription pattern of genes in M. huakuii bacteroids. A. sinicus root nodules induced by a M. huakuii wild-type or an abiEi gene mutant strain were collected at 28 dpi postinoculation. Total RNA was extracted, converted into cDNA libraries, and sequenced using Illumina paired-end sequencing technology. 52-million clean sequencing reads were gained from the RNA-seq transcriptomic analysis of the two samples, with an average of 26-million reads per sample. In total, 6590 expressed genes were detected by RNA-seq during the symbiosis between the strains and the A. sinicus. Gene expression level comparison analysis found that 116 genes were differentially expressed (p-value≤0.001, with log2 (FC) ≥ 3 and ≤ − 3), of which 5 were up-regulated and 111 were down-regulated (Tables 3 and 4). We validated the RNA-seq results by RT-PCR using several down-regulated genes as representatives. The datas showed that their expression was significantly decreased in mutant 28-day-old nodules compared to wildtype, which is consistent with the RNA-Seq data (Table 4). Among these differentially expressed genes, 23 (19.8%) were located on the chromosome, 54 (46.6%) were located on the plasmid pMHa, and 39 (33.6%) were located on the symbiotic plasmid pMHb, which contains nodulation genes (nod) as well as genes involved in nitrogen fixation (nif and fix). It is worth mentioning that all of the up-regulated genes were located on the chromosome, while more than 80% of the down-regulated genes were located on the extrachromosomal plasmids.

Table 3 List of 5 genes that showed significantly increased expression in Astragalus sinicus nodules, induced by the abiEi mutant strain compared to the wild-type 7653R
Table 4 List of 111 genes that showed significantly decreased expression in Astragalus sinicus nodules, induced by the abiEi mutant strain compared to the wild-type 7653R

Among the five up-regulated genes, one encodes type IVb pilin, one codes FMN transferase, one encodes cytochrome b561, and two encode proteins of unknown function. To categorize these differences into modules of biological relevance, the 111 down-regulated genes were annotated. They were functionally classified into 6 categories, which were involved in symbiotic N2-fixation and nitrogen mechanism (n = 21, 18.9%), type VI secretion system (n = 12, 10.8%), electron transport and antioxidant system (n = 16, 9.4%), nucleic acid metabolism (n = 27, 24.3%), Carbon mechanism (n = 11, 9.9%), and hypothetical protein (n = 16, 21.6%) (Table 4). In particular, one fix gene and seven nif genes were required for induction of nitrogen-fixing nodules on A. sinicus, and T6SS is rhizobial protein injection machinery with a positive role in Rhizobium-legume symbiosis [23]. The number of differentially expressed symbiosis-associated genes indicated that regulator AbiEi affects the transcription of a wide range of genes involved in the legume–rhizobium interactions. Further analysis of the differentially expressed genes identified a subset involved in electron transport and antioxidant system. The number of affected the expression of genes involved in the electron transport chain and antioxidant responses also suggested that regulator AbiEi plays an important role in antioxidant systems and the regulation of the electron transport chain. Moreover, 12 genes of T6SS were also found to be significantly over-represented. The genome of many rhizobia encodes T6SS but their role in symbiosis is mostly unknown. However, a functional T6SS, with an effect on symbiosis, has only been shown in R. leguminosarum [24]. Furthermore, 27 genes associated with nucleic acid metabolism were found among the 111 genes showing decreased expression in mutant nodules: 10 coding for transposases, six coding for transcriptional or response regulator. Previous studies suggested that transposases (or insertion sequence), transcriptional or response regulator are linked with functions to rhizobial nodulation [25,26,27,28,29]. In addition, 10 carbon mechanism genes and 24 genes of unknown function were found to be significantly down-regulated in the mutant nodules. qRT-PCR was further performed to confirm changes in gene expression determined by RNA-seq. Seven down-regulated genes in five different functional categories were significantly less in HKabiEi-induced nodules compared to nodules induced by wild-type. These results are largely consistent with the changes seen in the RNA-seq assay results.


The type IV TA system AbiE was recently reported, and belonged to the poorly characterized but widespread abortive infection/TA family [12]. To date, there is almost no information of rhizobial AbiE TA system. The Mesorhizobium huakuii AbiE system consists of toxin AbiEii and antitoxin AbiEi. Here we examine the antitoxin AbiEi, which is essential for transcriptional repression of the abiE operon. To investigate the mode of action of this antitoxin, we focus on a abiEi mutant strain of M. huakuii that is affected with regard to its symbiotic capacity and stress response.

TA systems are key regulators of bacterial persistence, and are linked to many roles in cell physiology, such as plasmid maintenance, stress response, antibiotic resistance, virulence and programmed cell death [30]. Our experiments indicated that the M. huakuii AbiE system can work as a functional TA module in the growth, stress response and antibiotic resistance. Firstly, the abiEi mutant showed slight growth inhibition effect when the strains were cultured in nutrient-rich medium. Due to lack of the antitoxin AbiEi, the abiEi mutant can not neutralize the toxin’s activity, and therefore, causes growth stasis. This finding is in agreement with two previous studies reporting that overexpression of the toxin gene in Escherichia coli inhibited its growth [13, 31]. Interestingly, the abiEi mutant showed no significant difference in AMS minimal medium, suggesting that AbiE TA system does not perform growth inhibition when strains grow at low cell density in nutrient-limiting minimal medium. Secondly, the abiEi mutant displayed significant difference on susceptibility of antibiotics tested. Mutantion in abiEi did not affect M. huakuii chloramphenicol resistance, but displayed an enhanced capacity of M. huakuii to gentamicin resistance [32]. It has been reported that the activity of TA system may contribute to the maintenance of antibiotic resistance. Thirdly, compared to wild-type, the abiEi mutant showed enhanced resistance to the oxidizing and reducing agents, H2O2 and SNP. It was recently reported that the production of Reactive Oxygen Species (ROS) is a common mechanism of cell death induced by bactericidal antibiotics [33].

The symbiotic performance of the abiEi mutant was altered at large levels of nodule function and maintenance. We have observed that the M. huakuii abiEi mutant induced the same number of nodules, exhibited a dramatical decrease in the nitrogen-fixing activity of root nodules (reduced by approx. 97%). In fact, a S. meliloti strain with a mutant antitoxin gene also had a reduced nitrogen-fixing activity, but formed a higher number of root nodules, which suggested that the host plant was starved for nitrogen due to the inefficient nitrogen fixation, thus initiation of additional nodules was permitted on the roots [10]. The reason why the M. huakuii abiEi mutant did not form more nodules may be due to decreased competitive ability. The root nodules induced by antitoxin gene mutant had decreased size, small amounts of infected cells and bacteroid content, and a severe senescent phenotype. These results indicated that inactive antitoxin AbiEi may affect the formation of the symbiotic nodules on the host roots. In fact, it has been reported that mutant in a functional toxin NtrR increased the symbiotic efficiency of S. meliloti [18], whereas a S. meliloti strain with a mutant antitoxin gene had a reduced nitrogen-fixing activity [10]. This symbiotic defect of the antitoxin AbiEi insertion mutant can be explained as a result of the disruption of balance between toxin and antitoxin. The abiEi mutant was defective in transcriptional repression of the abiE operon, and therefore, overexpressed the toxin AbiEii. On the one hand, many of toxins may induce programmed cell death by directly interacting with bacteroids [34]. Actually, the antitoxin was required for maintenance of transcriptional repression throughout nodule development as its expression was consistently significantly up-regulated in the nodules. On the other hand, a toxin VapC from the Leptospiral VapBC toxin-antitoxin module was reported to display ribonuclease activity on the initiator tRNA [35]. As a result, overexpressing the toxin AbiEii may result in a high percentage of down-regulated genes such as those which are functionally symbiosis related (Table 4).

Previous reports indicated that the Toxin-Antitoxin is a posttranscriptional regulator of metabolic flux, and antitoxin functions as a transcriptional autoregulator [36, 37]. Henceforth, a RNA-seq experiment was performed to compare the transcript profiles between the root bacteroids infected by abiEi mutant and wild-type 7653R. Firstly, lack of M. sativa antitoxin protein has been reported to continuously decrease expression of nif genes in the nodules [10], similar result also happened in the abiEi mutant. Lack of the antitoxin AbiEi reduced transcription of 8 nitrogen fixation genes and 13 genes related to nitrogen transport and mechanism. These suggest that antitoxin AbiEi plays a crucial role in legume-rhizobial symbiosis by maintaining the level of nitrogen assimilation and mechanism activity. Secondly, Lack of the antitoxin AbiEi reduced transcription of T6SS, which comprises a series of proteins with structural homology to bacteriophage tail proteins and membrane proteins [38]. It has been reported that T6SS plays a major role in mediating interbacterial competition and might contribute to virulence in plant pathogenic bacteria [39], and Rhizobium etli mutants affected in T6SS structural genes produced plants with lower dry weight and smaller nodules [23]. Thirdly, mutation in AbiEi decreased the activity of redox enzymes such as peroxiredoxin, oxidoreductase, dehydrogenase, oxidase and reductase. Those enzymes were able to effectively reduced intracellular ROS production and involved in redox balance and respiration [40]. It has been reported that antioxidant proteins are critical for nodule activity, and a better antioxidant metabolism can lead to delayed senescence of nodules [41,42,43]. Fourthly, five down-regulation genes are involved in electron transfer. It has been reported that nitric oxide (NO) was also found to play a metabolic role in nodule energy metabolism, and electron transfer chains were found to significantly contribute to NO production in N2-fixing nodules [44]. Fifthly, inactivation of the AbiEi response regulator directly reduced the activation of four transcriptional regulators (MCHK_RS31390, MCHK_RS32640, MCHK_RS31250 and MCHK_RS34740) and two response regulators (MCHK_RS31480 and MCHK_RS34845). Previous studies suggested that many transcriptional or response regulators are linked with functions to rhizobial nodulation [25,26,27,28]. Taken together, M. huakuii AbiEi plays an important role in root nodule symbiosis by regulation of nitrogen fixation gene expression, interbacterial competition, redox balance and respiration, bacteroid formation and senescence.


Bacterial Toxin-antitoxin systems are composed of bicistronic operons encoding a stable toxin that can harm the host cell and its cognate labile antitoxin. The contribution of toxin-antitoxin to symbiosis and stress response was investigated using the M. huakuii abiEi mutant. The results showed the abiEi mutant strain displayed decreased antioxidative capacity and enhanced gentamicin resistance, and was severely impaired in symbiotic nitrogen-fixing capacity. A quantitative RNA-Seq based transcriptomics approach was also applied to reveal the global transcriptomic responses to AbiEi defect in M. huakuii bacteroids from A. sinicus root nodules. Compared to the 7653R bacteroids, there were 5 genes were up-regulated and 111 genes were down-regulated in HKabiEi bacteroids. This study provided majority of these differentially expressed genes were grouped into 6 categories and a valuable insight into AbiEi-mediated mechanisms during M. huakuii-A. sinicus symbiosis. Furthermore, this study has generated an abundant list of transcript from M. huakuii which will provide a fundamental basis for future functional genomic research in M. huakuii and other closely related species.


Bacterial growth and media

The host plant is Astragalus sinicus, which is distributed in all provinces of the Yangtze River Basin in China and cultivated all over the country. The strains, plasmids and primers used in this study are listed in Table 5. M. huakuii strains were grown at 28 °C in either Tryptone Yeast extract (TY) [49] or Acid Minimal Salts medium (AMS) [50] with D-glucose (10 mM) as a carbon source and NH4Cl (10 mM) as a nitrogen source. For growth and qRT-PCR experiments, cells were grown in AMS. When required, the following antibiotics were used at the following final concentrations (μg mL− 1): Streptomycin (Str), 500; Ampicillin (Amp), 50; Kanamycin (Km), 20, or 50 (for E. coli growth); Neomycin (Neo), 80, or 250 (for making abiEi mutant); Gentamicin (Gm), 20; Spectinomycin (Spe), 100; Tetracycline (Tc), 5. To monitor culture growth, strains were grown at 28 °C with shaking (200 rpm) in liquid AMS or TY, and culture optical density at 600 nm (OD600) was measured during the culture period. For antibiotic and sodium nitroprusside (SNP) sensitivity assay, gentamicin, chloramphenicol and SNP were added to each AMS medium at a final concentration of 2 μg mL− 1, 2 μgmL− 1, and 10 μg mL− 1, respectively.

Table 5 Strains, plasmids and primers used in this experiment

Expression of AbiE system and AbiEii protein in E. coli

The AbiE system coding sequence of 1.42 kb was amplified by PCR from 7653R genomic DNA by using primers AbiE_F and AbiE_R, and the AbiEii coding sequence of 0.7 kb was amplified by using primers AbiE_F and AbiEii_R. PCR products were digested with BamHI and HindIII and cloned into the pET-28a(+), and the resulting plasmids were designated pETAbiE for AbiE system and pETAbiEii for toxin gene abiEii. The recombinant plasmids were further transformed in expression host E. coli BL21(DE3). Transformants obtained were grown in LB supplemented with appropriate antibiotics. The cultures were induced at OD600 0.4 with 1 mM isopropyl b-thiogalactopyranoside (IPTG) and samples were collected every hour for 6 h. The experiment was repeated three times.

Construction and complementation of abiEi gene mutant strain of M. huakuii 7653R

A single-crossover integration mutation in abiEi was made in 7653R. Primers abiEiUP and abiEiLW were used to PCR amplify the abiEi region from 7653R genomic DNA, and the 650 bp internal fragment of the abiEi gene was cloned into the HindIII and XbaI sites of pK19mob, giving plasmid pKabiEi. Plasmid pKabiEi was transferred from E. coli to 7653R and recombined into the genomic abiEi region via single crossover to give strain HKabiEi. Insertions into the abiEi gene of strain 7653R were confirmed by colony PCR using the abiEimap primer and a pK19mob-specific primer pK19A or pK19B [45, 49].

To complement the abiEi mutant, primers cabiEi_F and cabiEi_R were used to amplify the complete abiEi gene from M. huakuii 7653R genomic DNA. The PCR product was digested with HindIII and XbaI and cloned into the broad-host-range vector pBBR1MCS-5, resulting in plasmid pBBRabiEi. Plasmid pBBRabiEi was mated into the mutant strain HKabiEi using the triparental mating method as previously described, while conjugation of HKabiEi with empty vector (pBBR1MCS-5) was used as negative control.

Cellular sensitivity to H2O2

Rhizobial cultures were grown in TY medium up to an OD600 of 1, washed and resuspended in sterilized phosphate buffered saline (PBS) 1 × (136 mM NaCl, 2.6 mM KCl, 8.0 mM Na2HPO4, 1.5 mM KH2PO4). Antioxidation activity were determined by agar diffusion assay using filter paper disc diffusion method. Agar plates made from TY were spread with 100 μL of bacterial cultures containing 1 × 109 cfu mL − 1. Sterile and dried filter paper discs 6 mm in diameter, impregnated with different concentrations of H2O2, were placed on the TY media plate on which bacteria had been spread. Plates were incubated at 28 °C for 96 h. Results were determined by measuring the diameter of the zones of growth inhibition surrounding the disc. The presence of clear inhibition zones around the bacterial growth indicated the antimicrobial activity. The experiment was repeated at least three times to ensure the reproducibility of the results. The data were analyzed with two-way ANOVA (P < 0.05).

Plant experiment and cytological study of nodules

Astragalus sinicus L. was used as a host plant to test nodulation of the M. huakuii strains. Seeds were surface-sterilized, placed in 500 mL pots at two seeds per pot filled with sterile black vermiculite containing nitrogen-free Fahraeus solution. Inoculation with a bacterial suspension of about 108 cells/seed was performed on 7-day-old seedlings. The cultivation was carried out in a controlled environment chamber with 16 h light/8 h dark period (day/night temperature, 22 °C and 20 °C). Acetylene reduction rate per plant was determined at 31 days postinoculation (dpi) as previously described. The experiment consisted of two independent experiments, each of which had five repeats, and statistical differences were analyzed with one-way ANOVA (P < 0.05).

Nodules at 31 dpi were fixed for 12 h at 4 °C with 2.5% glutaraldehyde, rinsed, and post-fixed in 1.5% phosphate-buffered osmium tetroxide. Ultra-thin sections stained with lead citrate were examined using a Hitachi H-7100 transmission electron microscope. Sections were cut with a microtome and stained with toluidine blue for light microscopy.

Rhizosphere colonization

Rhizosphere colonization was performed as described previously. Astragalus sinicus seedlings were germinated and grown for 7 days as described above for acetylene reduction activity, and inoculated with M. huakuii 7653R and HKabiEi in the cfu ratios 1000:0, 0:1000, 1000:1000 and 1000:10000. Shoots were cut-off after 7 days (14 days after plant), and 10 mL of sterile phosphate-buffered saline (PBS) buffer (pH 7.4) was added to the roots and vortexed for 15 mins. The samples were further serially diluted and plated on TY agar plates containing either streptomycin and neomycin (only HKabiEi will grow) or streptomycin (both 7653R and HKabiEi will grow), giving the total number of viable rhizosphere- and root-associated bacteria. The ratio of the number of mutant strains to total number of mutant and wild-type strains was calculated for each inoculation test. Each treatment consisted of 10 replications, each of which consisted of a single plant, and statistical differences were analyzed with one-way ANOVA (P < 0.05).

RNA isolation and quantitative RT-PCR analysis

Quantitative real-time reverse transcription PCR (qRT-PCR) was performed to determine the abiEi gene expression level in A. sinicus nodules and validate the results of RNA-seq. The total RNA was isolated using Trizol reagent from free-living M. huakuii 7653R or mutant HKabiEi cultivated in AMS liquid medium, or root nodules which were harvested from A. sinicus inoculated with strain 7653R or mutant HKabiEi. RNA were reverse transcribed into cDNA using the SuperScript II reverse transcriptase and random hexamers. qRT-PCR analysis was performed using a SYBR Premix ExTaq kit following the manufacturer’s instructions on the BIO-RAD CFX96 Real-Time PCR Detection System. The primer sequences are shown in Table 5. The 16S rRNA gene of M. huakuii 7653R was used as a calibrator gene, and the data were obtained. Three independent biological replicates were included in the analysis, and the relative expression levels of the target genes were normalized using the 2−ΔΔCT method.

RNA-seq library preparation and sequencing using the illumina genome analyzer

At 4 weeks post-inoculation, the nodules of plants inoculated with HKabiEi or 7653R were harvested, immediately frozen in liquid nitrogen and stored at − 80 °C. Total cellular RNA was isolated from frozen nodule tissues using Trizol Reagent (Invitrogen) and RNeasy Mini Kit (Qiagen). Total RNA of each nodule sample was treated with RNase-free DnaseI (TAKARA, Dalian, China), and the absence of contaminating DNA was confirmed by PCR-based assays. Total RNA was assessed using Agilent 2100 Bioanalyzer (Agilent Technologies, Palo Alto, CA, USA), and NanoDrop (Thermo Fisher Scientific Inc). 1 μg total RNA with RNA integrity number (RIN) value above 6.5 was used for following library preparation. The rRNA was depleted from total RNA using Ribo-Zero rRNA Removal Kit (Bacteria) (Illumina) [51]. The ribosomal depleted RNA was then fragmented and reverse-transcribed into cDNA with random primers. The purified double-stranded cDNA by beads was then treated with End Prep Enzyme Mix to repair both ends and add a dA-tailing in one reaction, followed by a T-A ligation to add adaptors to both ends. Next generation sequencing library preparations were constructed according to the manufacturer’s protocol. The Qsep100 (Bioptic, Taiwan, China) and Qubit 3.0 Fluorometer was used to determine the quality of the libraries.

The libraries with different indices were multiplexed and sequenced on an Illumina HiSeq instrument according to manufacturer’s instructions (Illumina, San Diego, CA, USA). Sequencing was carried out by Illumina paired-end configuration. The sequencing image processing and base calling were conducted following to Illumina’s protocol on the HiSeq instrument. The read length was 90 bp, and reads were mapped to the M. huakuii genome using Bowtie 2. The HTSeq-count tool was used to generate the total number of uniquely mapped reads for each gene. Three independent biological replicates per sample were processed and sequenced.

Data analysis

Differences between the average of gene expression for the control and experimental groups were analyzed by the Student’s t-test using SPSS software, version 18 (SPSS, Inc., Chicago, IL). For the RNA-seq study, the unique reads mapping to the M. huakuii genome were used for a differential gene expression analysis using the DESeq2 package. The P-values with false discovery rate were adjusted for multiple testing. The false discovery rate P-value < 0.001 and the absolute value of log2 (FC) ≥ 3 and ≤ − 3 were used to identify statistically significant changes in gene expression. For quantitative RT-PCR analysis, p < 0.05 was considered to be statistically significant.

Availability of data and materials

Raw sequence data from these RNA-seq studies can be accessed via the NCBI Sequence Read Archive, with accession number PRJNA718999.



Bacterial abortive infection




Type VI secretion system


Quantitative real-time PCR


Reactive oxygen species


Nitric oxide


Scanning electron microscopy


  1. 1.

    Schuster CF, Bertram R. Toxin-antitoxin systems are ubiquitous and versatile modulators of prokaryotic cell fate. FEMS Microbiol Lett. 2013;340(2):73–85.

    CAS  Article  PubMed  Google Scholar 

  2. 2.

    Melderen LV, Bast M. Bacterial toxin–antitoxin systems: more than selfish entities? PLoS Genet. 2009;5(3):e1000437.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  3. 3.

    Wen Z, Wang P, Sun C, Guo Y, Wang X. Interaction of type IV toxin/antitoxin systems in cryptic prophages of Escherichia coli K-12. Toxins (Basel). 2017;9(3):77.

    CAS  Article  PubMed Central  Google Scholar 

  4. 4.

    Maisonneuve E, Shakespeare LJ, Jørgensen MG, Gerdes K. Bacterial persistence by RNA endonucleases. Proc Natl Acad Sci U S A. 2011;108(32):13206–11.

    Article  PubMed  PubMed Central  Google Scholar 

  5. 5.

    Aizenman E, Engelberg-Kulka H, Glaser G. An Escherichia coli chromosomal "addiction module" regulated by guanosine [corrected] 3',5'-bispyrophosphate: a model for programmed bacterial cell death. Proc Natl Acad Sci U S A. 1996;93(12):6059–63.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  6. 6.

    Wang X, Wood TK. Toxin-antitoxin systems influence biofilm and persister cell formation and the general stress response. Appl Environ Microbiol. 2011;77(16):5577–83.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  7. 7.

    Magnuson RD. Hypothetical functions of toxin-antitoxin systems. J Bacteriol. 2007;189(17):6089–92.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  8. 8.

    Lobato-Márquez D, Díaz-Orejas R, García-Del Portillo F. Toxin-antitoxins and bacterial virulence. FEMS Microbiol Rev. 2016;40(5):592–609.

    CAS  Article  PubMed  Google Scholar 

  9. 9.

    Page R, Peti W. Toxin-antitoxin systems in bacterial growth arrest and persistence. Nat Chem Biol. 2016;12(4):208–14.

    CAS  Article  PubMed  Google Scholar 

  10. 10.

    Lipuma J, Cinege G, Bodogai M, Oláh B, Kiers A, Endre G, et al. A vapBC-type toxin-antitoxin module of Sinorhizobium meliloti influences symbiotic efficiency and nodule senescence of Medicago sativa. Environ Microbiol. 2014;16(12):3714–29.

    CAS  Article  PubMed  Google Scholar 

  11. 11.

    Deng YM, Liu CQ, Dunn NW. Genetic organization and functional analysis of a novel phage abortive infection system, AbiL, from Lactococcus lactis. J Biotechnol. 1999;67(2–3):135–49.

    CAS  Article  PubMed  Google Scholar 

  12. 12.

    Hampton HG, Jackson SA, Fagerlund RD, Vogel AIM, Dy RL, Blower TR, et al. AbiEi binds cooperatively to the type IV abiE toxin-antitoxin operator via a positively-charged surface and causes DNA bending and negative autoregulation. J Mol Biol. 2018;430(8):1141–56.

    CAS  Article  PubMed  Google Scholar 

  13. 13.

    Dy RL, Przybilski R, Semeijn K, Salmond GP, Fineran PC. A widespread bacteriophage abortive infection system functions through a Type IV toxin-antitoxin mechanism. Nucleic Acids Res. 2014;42(7):4590–605.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  14. 14.

    Miklashevichs E, Rohrig H, Schell J, Schmidt J. Perception and signal transduction of rhizobial NOD factors. Crit Rev Plant Sci. 2001;20(4):373–94.

    CAS  Article  Google Scholar 

  15. 15.

    Theunis M, Kobayashi H, Broughton WJ, Prinsen E. Flavonoids, NodD1, NodD2, and Nod-box NB15 modulate expression of the y4wEFG locus that is required for indole-3-acetic acid synthesis in Rhizobium sp. strain NGR234. Mol Plant-Microbe Interact. 2004;17(10):1153–61.

  16. 16.

    Miclea PS, Péter M, Végh G, Cinege G, Kiss E, Váró G, et al. Atypical transcriptional regulation and role of a new toxin-antitoxin-like module and its effect on the lipid composition of Bradyrhizobium japonicum. Mol Plant-Microbe Interact. 2010;23(5):638–50.

    CAS  Article  PubMed  Google Scholar 

  17. 17.

    Milunovic B, di Cenzo GC, Morton RA, Finan TM. Cell growth inhibition upon deletion of four toxin-antitoxin loci from the megaplasmids of Sinorhizobium meliloti. J Bacteriol. 2014;196(4):811–24.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  18. 18.

    Oláh B, Kiss E, Györgypál Z, Borzi J, Cinege G, Csanádi G, et al. Mutation in the ntrR gene, a member of the vap gene family, increases the symbiotic efficiency of Sinorhizobium meliloti. Mol Plant-Microbe Interact. 2001;14(7):887–94.

    Article  PubMed  Google Scholar 

  19. 19.

    Fernández-Bachiller MI, Brzozowska I, Odolczyk N, Zielenkiewicz U, Zielenkiewicz P, Rademann J. Mapping protein-protein interactions of the resistance-related bacterial zeta toxin-epsilon antitoxin complex (ε2ζ2) with high affinity peptide ligands using fluorescence polarization. Toxins. 2016;8(7):222.

    CAS  Article  PubMed Central  Google Scholar 

  20. 20.

    Park SJ, Son WS, Lee BJ. Structural overview of toxin-antitoxin systems in infectious bacteria: a target for developing antimicrobial agents. Biochim Biophys Acta. 2013;1834(6):1155–67.

    CAS  Article  PubMed  Google Scholar 

  21. 21.

    Tang G, Li N, Liu Y, Yu L, Yan J, Luo L. Sinorhizobium meliloti glutathione reductase is required for both redox homeostasis and symbiosis. Appl Environ Microbiol. 2018;84(3):e01937–17.

    Article  PubMed  PubMed Central  Google Scholar 

  22. 22.

    Triplett LR, Shidore T, Long J, Miao J, Wu S, Han Q, et al. AvrRxo1 is a bifunctional Type III secreted effector and toxin-antitoxin system component with homologs in diverse environmental contexts. PLoS One. 2016;11(7):e0158856.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  23. 23.

    Salinero-Lanzarote A, Pacheco-Moreno A, Domingo-Serrano L, Durán D, Ormeño-Orrillo E, Martínez-Romero E, et al. The Type VI secretion system of Rhizobium etli Mim1 has a positive effect in symbiosis. FEMS Microbiol Ecol. 2019;95(5):fiz054.

  24. 24.

    Nelson MS, Sadowsky MJ. Secretion systems and signal exchange between nitrogen-fixing rhizobia and legumes. Front Plant Sci. 2015;6:491.

  25. 25.

    Loh J, Lohar DP, Andersen B, Stacey G. A two-component regulator mediates population-density-dependent expression of the Bradyrhizobium japonicum nodulation genes. J Bacteriol. 2002;184(6):1759–66.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  26. 26.

    Janczarek M, Skorupska A. The Rhizobium leguminosarum bv. Trifolii RosR: transcriptional regulator involved in exopolysaccharide production. Mol Plant-Microbe Interact. 2007;20(7):867–81.

    CAS  Article  PubMed  Google Scholar 

  27. 27.

    Del Cerro P, Rolla-Santos AA, Valderrama-Fernández R, Gil-Serrano A, Bellogín RA, Gomes DF, et al. NrcR, a new transcriptional regulator of Rhizobium tropici CIAT 899 involved in the legume root-nodule symbiosis. PLoS One. 2016;11(4):e0154029.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  28. 28.

    Bolzan de Campos S, Deakin WJ, Broughton WJ, Passaglia LMP. Roles of flavonoids and the transcriptional regulator TtsI in the activation of the type III secretion system of Bradyrhizobium elkanii SEMIA587. Microbiology. 2011;157(3):627–35.

  29. 29.

    Selbitschka W, Zekri S, Schröder G, Pühler A, Toro N. The Sinorhizobium meliloti insertion sequence (IS) elements ISRm102F34-1/ISRm7 and ISRm220-13-5 belong to a new family of insertion sequence elements. FEMS Microbiol Lett. 1999;172(1):1–7.

    CAS  Article  PubMed  Google Scholar 

  30. 30.

    Otsuka Y. Prokaryotic toxin-antitoxin systems: novel regulations of the toxins. Curr Genet. 2016;62(2):379–82.

    CAS  Article  PubMed  Google Scholar 

  31. 31.

    Li Z, Song Q, Wang Y, Xiao X, Xu J. Identification of a functional toxin-antitoxin system located in the genomic island PYG1 of piezophilic hyperthermophilic archaeon Pyrococcus yayanosii. Extremophiles. 2018;22(3):347–57.

    CAS  Article  PubMed  Google Scholar 

  32. 32.

    Wozniak RA, Waldor MK. A toxin-antitoxin system promotes the maintenance of an integrative conjugative element. PLoS Genet. 2009;5(3):e1000439.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  33. 33.

    Kolodkin-Gal I, Sat B, Keshet A, Engelberg-Kulka H. The communication factor EDF and the toxin-antitoxin module mazEF determine the mode of action of antibiotics. PLoS Biol. 2008;6(12):e319.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  34. 34.

    Hazan R, Sat B, Engelberg-Kulka H. Escherichia coli mazEF-mediated cell death is triggered by various stressful conditions. J Bacteriol. 2004;186(11):3663–9.

  35. 35.

    Lopes AP, Lopes LM, Fraga TR, Chura-Chambi RM, Sanson AL, Cheng E, et al. VapC from the leptospiral VapBC toxin-antitoxin module displays ribonuclease activity on the initiator tRNA. PLoS One. 2014;9(7):e101678.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  36. 36.

    McKenzie JL, Robson J, Berney M, Smith TC, Ruthe A, Gardner PP, et al. A VapBC toxin-antitoxin module is a posttranscriptional regulator of metabolic flux in mycobacteria. J Bacteriol. 2012;194(9):2189–204.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  37. 37.

    Habib G, Zhu Q, Sun B. Bioinformatics and functional assessment of toxin-antitoxin systems in Staphylococcus aureus. Toxins. 2018;10(11):473.

    CAS  Article  PubMed Central  Google Scholar 

  38. 38.

    Gmiter D, Czerwonka G, Kaca WR. Type VB and VI secretion systems as competition agents of gram-negative bacteria[J]. Postep Mikrobiol. 2018;57(4):360–73.

  39. 39.

    Tian Y, Zhao Y, Shi L, Cui Z, Hu B, Zhao Y. Type VI secretion systems of Erwinia amylovora contribute to bacterial competition, virulence, and exopolysaccharide production. Phytopathology. 2017;107(6):654–61.

    CAS  Article  PubMed  Google Scholar 

  40. 40.

    Hu A, Chen X, Luo S, Zou Q, Xie J, He D, et al. Rhizobium leguminosarum glutathione peroxidase is essential for oxidative stress resistance and efficient nodulation. Front Microbiol. 2021;12:627562.

    Article  PubMed  PubMed Central  Google Scholar 

  41. 41.

    Yurgel SN, Simarov BV. Rhizobial electron transport chain and its genetic determinants. Russ J Genet. 1999;35(2):101–9.

    CAS  Google Scholar 

  42. 42.

    Becana M, Matamoros MA, Udvardi M, Dalton DA. Recent insights into antioxidant defenses of legume root nodules. New Phytol. 2010;188(4):960–76.

    CAS  Article  PubMed  Google Scholar 

  43. 43.

    Redondo FJ, de la Peña TC, Morcillo CN, Lucas MM, Pueyo JJ. Overexpression of flavodoxin in bacteroids induces changes in antioxidant metabolism leading to delayed senescence and starch accumulation in alfalfa root nodules. Plant Physiol. 2009;149(2):1166–78.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  44. 44.

    Puppo A, Pauly N, Boscari A, Mandon K, Brouquisse R. Hydrogen peroxide and nitric oxide: key regulators of the legume-Rhizobium and mycorrhizal symbioses. Antioxid Redox Signal. 2013;18(16):2202–19.

    CAS  Article  PubMed  Google Scholar 

  45. 45.

    Cheng GJ, Li YG, Zhou JC. Cloning and identification of opa22, a new gene involved in nodule formation by Mesorhizobium huakuii. FEMS Microbiol Lett. 2006;257(1):152–7.

    CAS  Article  PubMed  Google Scholar 

  46. 46.

    Schäfer A, Tauch A, Jäger W, Kalinowski J, Thierbach G, Pühler A. Small mobilizable multi-purpose cloning vectors derived from the Escherichia coli plasmids pK18 and pK19: selection of defined deletions in the chromosome of Corynebacterium glutamicum. Gene. 1994;145(1):69–73.

    Article  PubMed  Google Scholar 

  47. 47.

    Figurski DH, Helinski DR. Replication of an origin-containing derivative of plasmid RK2 dependent on a plasmid function provided in trans. Proc Natl Acad Sci U S A. 1979;76(4):1648–52.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  48. 48.

    Kovach ME, Phillips RW, Elzer PH, Roop RM 2nd, Peterson KM. pBBR1MCS: a broad-host-range cloning vector. Biotechniques. 1994;16(5):800–2.

    CAS  PubMed  Google Scholar 

  49. 49.

    Zou Q, Zhou Y, Cheng G, Peng Y, Luo S, Wu H, et al. Antioxidant ability of glutaredoxins and their role in symbiotic nitrogen fixation in Rhizobium leguminosarum bv. Viciae 3841. Appl Environ Microbiol. 2020;87(4):e01956–20.

    Article  Google Scholar 

  50. 50.

    Cheng G, Karunakaran R, East AK, Munoz-Azcarate O, Poole PS. Glutathione affects the transport activity of Rhizobium leguminosarum 3841 and is essential for efficient nodulation. FEMS Microbiol Lett. 2017;364(8):fnx045.

  51. 51.

    Owens NDL, De Domenico E, Gilchrist MJ. An RNA-Seq protocol for differential expression analysis. Cold Spring Harb Protoc. 2019;2019(6):514.

Download references


Not applicable.


This study was supported by the National Natural Science Foundation of China (grant no. 31772399). The funders had no role in design of the study, analysis and interpretation of data and in the preparation of the manuscript.

Author information




GC conceived and designed the study. XC, AH, QZ, HW and SL performed the experiments. GC, CY, TL, DH and XL analyzed the results, GC, XC and QZ wrote the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Guojun Cheng.

Ethics declarations

Ethics approval and consent to participate

Ethical approval according to national legislation (the Law of the People’s Republic of China on the Protection of Wildlife) was not necessary as Astragalus sinicus is by far the wide spread green manure of rice field in China, and could improve soil fertility, increase gain yield, benefit nitrogen fixation and energy saving.

Consent for publication

Not applicable.

Competing interests

The authors declare that the research was conducted in the absence of any competing interest or commercial or financial relationships that could be construed as a potential conflict of interest.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Fig. S1.

Inhibition zones of the H2O2 in disc diffusion test of different strains. Filter paper discs were impregnated with a solution containing 20 (A-D), 100 (E-H) and 250 (I-L) mg L− 1 of H2O2, and placed on the TY media plate on which bacteria had been spread. Plates were incubated at 28 °C for 96 h. A, E, I, M. huakuii 7653R; B, F, J, M. huakuii HKabiEi; C, G, K, M. huakuii HKabiEi(pBBRabiEi); D, H, L, M. huakuii HKabiEi(pBBR1MCS-5).

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Chen, X., Hu, A., Zou, Q. et al. The Mesorhizobium huakuii transcriptional regulator AbiEi plays a critical role in nodulation and is important for bacterial stress response. BMC Microbiol 21, 245 (2021).

Download citation


  • Type IV toxin–antitoxin system
  • Mesorhizobium huakuii
  • Transcriptional regulator AbiEi
  • Symbiosis nitrogen fixation
  • RNA-Seq analysis of nodule bacteroids