- Research article
- Open Access
Characterisation of the mgo operon in Pseudomonas syringae pv. syringae UMAF0158 that is required for mangotoxin production
BMC Microbiology volume 12, Article number: 10 (2012)
Mangotoxin is an antimetabolite toxin that is produced by strains of Pseudomonas syringae pv. syringae; mangotoxin-producing strains are primarily isolated from mango tissues with symptoms of bacterial apical necrosis. The toxin is an oligopeptide that inhibits ornithine N-acetyl transferase (OAT), a key enzyme in the biosynthetic pathway of the essential amino acids ornithine and arginine. The involvement of a putative nonribosomal peptide synthetase gene (mgoA) in mangotoxin production and virulence has been reported.
In the present study, we performed a RT-PCR analysis, insertional inactivation mutagenesis, a promoter expression analysis and terminator localisation to study the gene cluster containing the mgoA gene. Additionally, we evaluated the importance of mgoC, mgoA and mgoD in mangotoxin production. A sequence analysis revealed an operon-like organisation. A promoter sequence was located upstream of the mgoB gene and was found to drive lacZ transcription. Two terminators were located downstream of the mgoD gene. RT-PCR experiments indicated that the four genes (mgoBCAD) constitute a transcriptional unit. This operon is similar in genetic organisation to those in the three other P. syringae pathovars for which complete genomes are available (P. syringae pv. syringae B728a, P. syringae pv. tomato DC3000 and P. syringae pv. phaseolicola 1448A). Interestingly, none of these three reference strains is capable of producing mangotoxin. Additionally, extract complementation resulted in a recovery of mangotoxin production when the defective mutant was complemented with wild-type extracts.
The results of this study confirm that mgoB, mgoC, mgoA and mgoD function as a transcriptional unit and operon. While this operon is composed of four genes, only the last three are directly involved in mangotoxin production.
Antimetabolite toxins are generally small metabolites that exhibit strong effects in plant cells by causing an increase in disease symptoms . Various toxic substances produced by pathovars of Pseudomonas syringae have been well characterised. Each antimetabolite toxin inhibits a specific step in the glutamine and arginine biosynthesis pathways of the host, enhancing disease symptoms and increasing the virulence of the bacterial pathogen. The most well-studied antimetabolite toxins are tabtoxin and phaseolotoxin .
Tabtoxin is a monocyclic β-lactam that specifically inhibits the enzyme glutamine synthetase (GS, EC 184.108.40.206). This toxin is produced by P. syringae pv. tabaci, pv. coronafaciens and pv. garcae . The biosynthetic pathway of tabtoxin is not well understood, and tabtoxin biosynthesis may diverge from the lysine biosynthetic pathway prior to the formation of diaminopimelate [4, 5]. A genetic analysis of tabtoxin production revealed the presence of biosynthetic genes at the att site adjacent to the lysC tRNA gene in Pseudomonas syringae BR2 . The various ORFs within this region include sequences similar to β-lactam synthase, clavaminic acid synthase and enzymes involved in amino acid synthesis. Additionally, novel ORFs were identified in a portion of the biosynthetic region that is known to be associated with a toxin hypersensitivity phenotype .
Phaseolotoxin is produced mainly by P. syringae pv. phaseolicola and pv. actinidae. The molecular structure of phaseolotoxin includes a sulphodiaminophosphinyl moiety linked to a tripeptide of ornithine, alanine and homoarginine . Phaseolotoxin inhibits ornithine carbamoyltransferase (OCT, EC 220.127.116.11) . The phaseolotoxin homoarginine and ornithine residues are synthesised by a transamidation reaction that requires arginine and lysine [8, 9].
Aguilera et al.  have reported a biosynthetic cluster, pht, which is composed of 23 genes flanked by insertion sequences and transposases, that is involved in the biosynthesis of phaseolotoxin. Mutations of 11 of the genes within the cluster led to a Tox- phenotype, and the mutation of three additional genes resulted in low levels of toxin production. Preliminary results also indicated that the product of phtL may be involved in the regulation of phaseolotoxin biosynthesis .
Pseudomonas syringae pv. syringae (Pss) is a pathogenic bacterium that can cause canker, blossom blights and leaf spots in more than 200 different plant species, many of which are of economic importance . Strains of this pathovar can cause bacterial apical necrosis on mango trees, limiting mango production in the Mediterranean area . More than 86% of the Pss strains isolated from mango tissues produce mangotoxin, an antimetabolite toxin that inhibits ornithine N-acetyl-transferase (OAT), a key enzyme in the biosynthesis of arginine . Mangotoxin also acts as a virulence factor that increases the necrotic symptoms of Pss strains during the infection of plant tissues . In a previous study, a DNA fragment from Pss, UMAF0158, was cloned into pCG2-6 and sequenced (DQ532441), revealing a cluster of 4 ORFs that included the mgoA gene. Our group identified mgoA as the first P. syringae pv. syringae gene known to be directly involved in mangotoxin production . This gene encodes a putative nonribosomal peptide synthetase (NRPS), and its inactivation by insertional mutagenesis abolishes mangotoxin production and drastically reduces virulence [14, 15]. The genetic organisation of the three remaining genes and their roles in the production of mangotoxin remain unknown. The goal of our current study is to determine the organisation of the four ORFs in this cluster (Figure 1) and their relative importance in the production of mangotoxin.
The gene cluster containing mgoA may constitute an operon composed of four ORFs. Our current study provides insight into the organisation of the operon and the involvement of the genes in the production of mangotoxin.
The construction and characterisation of insertion mutants derived from Pseudomonas syringae pv. syringaeUMAF0158
Each ORF that was cloned into plasmid pCG2-6 (Figure 1) was subjected to insertional inactivation mutagenesis in the P. syringae pv. syringae UMAF0158 chromosome by integration of the appropriately cloned PCR products. The ORFs were 92%-98% identical to the homologous genes in P. syringae pv. syringae strain B728a (accession no. CP000075, Table 1). The deduced ORF0 and ORF1 protein products are homologous to proteins of the HAD hydrolase family and aldo-keto oxidoreductases, respectively. The mutation of these ORFs by insertional inactivation did not affect mangotoxin production. ORF2 is located just upstream of the putative mgo operon (Figure 1) and contains a putative ribosomal binding site (RBS) at nucleotide -6 (AAGAAGT). This gene is 97% identical to Psyr_5008 from P. syringae pv. syringae B728a (Table 1), PSPTO_5454 from P. syringae pv. tomato DC3000 and PSPPH_5087 from P. syringae pv. phaseolicola 1448A. The protein products of the genes from each of these bacteria were annotated in the database as members of the GntR family of transcriptional regulators . When ORF2 was disrupted, the corresponding mutant UMAF0158::ORF2 still produced mangotoxin (Tables 1 and 2).
The four genes downstream of ORF2 exhibit a high degree of identity to four consecutive P. syringae pv. syringae B728a genes (Psyr_5009 to Psyr_5012) (Table 1). The mgoB gene, which contains a putative RBS at nucleotide -8 (AGGA), is 96% similar to Psyr_5009, which encodes a conserved hypothetical protein. The mgoB mutant UMAF0158::mgoB produced mangotoxin (Table 1), although the level of mangotoxin was decreased slightly (Table 2). A search of the Pfam database revealed a similarity to DUF3050, a protein of unknown function, between amino acids 15 and 244 with an e-value of 3.1e-97. Searches in the InterProScan (EMBL-EBI) database revealed that the theoretical MgoB protein product is similar to the haem oxygenase-like, multi-helical superfamily between amino acids 128 and 245 (e-value of 1.3e-8).
The inactivation of the mgoC, mgoA and mgoD genes yielded mutants (UMAF0158::mgoC, UMAF0158::mgoA and UMAF0158::mgoD) that were completely unable to produce mangotoxin (Tables 1 and 2). The mgoC gene, which contains a putative RBS at -7 (AAGGA), exhibits 95% similarity to the Psyr_5010 gene of P. syringae pv. syringae B728a, a conserved hypothetical protein (Table 1). Homology searches for the MgoC protein product in the Pfam database revealed a significant match with the p-aminobenzoate N-oxygenase AurF from Streptomyces thioluteus. The alignment was between amino acids 2 and 295 with an e-value of 7.2e-88. The disruption of mgoA by insertion (UMAF0158::mgoA) or miniTn5 (UMAF0158-6γF6) mutation resulted in the complete inactivation of a putative nonribosomal peptide synthetase (NRPS). The inactivation of mgoA has previously been shown to result in defects in mangotoxin production and considerably reduced virulence . However, a putative RBS for mgoA could not be located using the consensus sequences published to date. Finally, insertional mutagenesis of the mgoD gene, which contains a putative RBS at -6 (ATGGAG), resulted in the inactivation of a conserved hypothetical protein that is 94% identical to Psy_5012. A conserved-domain analysis of the hypothetical amino acid sequence of MgoD revealed sequence similarity to Polyketide_cyc2, a polyketide cyclase/dehydrase and lipid transporter domain, from amino acids 20 to 158. The e-values were 1e-17 (Specialized BLAST-NCBI) and 1.6e-23 (Pfam).
The genetic organisation of the mgooperon and complementation of insertional mutants
To define the mgo operon and determine its genetic organisation and co-transcription, reverse-transcription PCR (RT-PCR) experiments were performed (Figure 2). The total DNA and RNA from wild-type UMAF0158 grown in PMS minimal medium at 22°C were used, and the RT-PCR primers were designed to anneal between the ORFs. The total DNA was used as an amplification control, and the cDNA derived from the mRNA was used to detect the transcripts of genes belonging to the putative mgo operon. To confirm the co-transcription of mgoB, mgoC, mgoA and mgoD, we amplified the connecting areas between the sequential ORFs of the putative mgo operon (Figure 2A). Sequences within ORF2 and mgoB were also amplified to determine their mRNA transcripts (Figure 2A, B). Our results indicated that ORF2 and the upstream region and mgoB and the downstream region were amplified. However, there was no amplification of the inter-genetic region upstream of mgoB. These results suggest that the transcriptional unit is mgoB, mgoC, mgoA and mgoD (Figure 2B). The lack of amplification between ORF2 and mgoB supports the presence of a putative promoter in this DNA sequence.
Additional RT-PCR experiments showed that only the disrupted mgoB gene was not amplified in UMAF0158::mgoB while the transcripts of the disrupted mgoC gene as well as that of the downstream genes were absent in UMAF0158::mgoC (Figure 2C). A hybridisation analysis of the transcript of the mgo operon with the total mRNA from wild-type UMAF0158 and the insertional mutants UMAF0158::mgoB, and UMAF0158::mgoC showed that the transcript was present in the wild-type strain and reduced in the mgoB mutant strain (Figure 2D).
To confirm the role of these genes in mangotoxin production and to analyse the specific phenotype of each mutation, we performed a complementation analysis using plasmids containing all of the genes that were situated downstream of the mutations (Table 3). The mgo genes were cloned downstream of the PLAC promoter. Plasmid pLac36, which contains the structural genes of the operon (mgoB, mgoC, mgoA and mgoD), and a plasmid containing the genomic clone pCG2-6 were both able to restore mangotoxin production in all of the constructed mutants (Tables 3 and 2). These results demonstrate that the complemented plasmids were functional and rule out the possibility that secondary mutations influence mangotoxin production. Plasmid pLac56, which contains only mgoA and mgoD, was able to complement the phenotypes of the miniTn5 mutant UMAF0158-6γF6 and the insertional mutants UMAF0158::mgoA and UMAF0158::mgoD. Plasmid pLac6, however, was only able to complement UMAF0158::mgoD (Table 3). These complementation experiments show that the insertional mutants UMAF0158::mgoC, UMAF0158::mgoA and UMAF0158::mgoD were unable to produce mangotoxin even when the downstream genes were restored on a plasmid. The insertional mutation of the mgoC, mgoA and mgoD genes resulted in a loss of mangotoxin activity, which did not occur when mgoB was mutated (Tables 1 and 2). Therefore, we cannot eliminate the possibility that a polar effect of the insertional mutations affected the phenotypes of the mutants and downstream genes transcription. Apparently the insertional mutation in mgoB did not show polar effect on mgo genes located downstream (mgoC, mgoA and mgoD), in contrast with the insertional mutation in mgoC, which produce a polar effect on mgo downstream genes transcription (Figure 2, Table 3).
Mangotoxin production in mutants derived from Pseudomonas syringae pv. syringaeUMAF0158
To further support our results, we determined the amount of mangotoxin production in the insertional and miniTn5 mutants relative to wild-type UMAF0158 (Table 2). The production of the syringomycin complex by the insertional mutants confirmed that only mangotoxin production was affected (data not shown). The results obtained from the quantitative mangotoxin analysis indicated that the two miniTn5 mutants that were complemented with pCG2-6, UMAF2-6A and UMAF2-6-3H1, and the insertion mutant UMAF0158::ORF1 were able to produce mangotoxin at the same level as wild-type UMAF0158.
Upon complementation with pLac56 (mgoA and mgoD), mangotoxin production was restored in the mutants UMAF0158::ORF2 and UMAF0158::mgoB and the miniTn5 mutant UMAF0158-6γF6; however, the production was slightly lower and could be detected only until a 1:4 dilution (Table 2).
Promoter and terminator localisation in the mgo operon
Promoter expression and terminator localisation experiments were performed to characterise the structure of the operon.
The promoter prediction software BPROM (SoftBerry Inc.) was used to identify possible promoters in the putative mgo operon. The best candidates were found in the nucleotide sequence (814 bp) of the non-coding region located upstream of the mgoB gene. Two possible promoters were predicted and designated as P mgo . The first predicted promoter was located at position 134 from 5'-end with a linear discriminant function (LDF) of 0.59, a -10 box, CGTTTTTAT, at position 119 (score: 37) and a -35 box, TCGCCA, at position 95 (score: 24). The second predicted promoter was located at position 549 from the 5'-end of the sequence, with an LDF of 4.38, a -10 box, TGATAAATT, at position 534 (score: 55) and a -35 box, TTAAAA, at position 513 (score: 37) (Figure 3C). The scores of the first predicted promoter were lower than those of the second promoter. According to the in silica prediction, the 814 bp sequence containing both putative promoters was cloned into pMP220, and its activity was measured with a β-galactosidase assay (β-Gal) [17, 18]. The P mgo studies were performed in Pseudomonas fluorescens Pf-5, which contains no genomic sequences that are homologous to the mgo operon, and P. syringae pv. syringae B728a, which contains genomic sequences that are homologous to the mgo operon but cannot produce mangotoxin. We also used the insertional mutant UMAF0158::ORF2, which contains a disruption in the putative transcriptional regulator gene, and wild-type UMAF0158. P mgo activity was measured in three different culture media (LB, KB and PMS) and at two growth temperatures (28°C and 22°C). In the minimal medium PMS, the P mgo promoter was active in the wild-type strain at both temperatures and in the insertional mutant at 22°C (Figure 4). The β-Gal assays of the strains grown in rich LB and KB media did not indicate activity in any of the strains at either temperature (data not shown).
Once the presence of promoter activity in the analysed sequence was confirmed, the 5'RACE method was used to determine the transcript start point of the mgo operon (Figure 3A, B). With this method, we could determine which of the two putative promoters of the mgo operon was the functional promoter and also analyse the presence of an additional promoter between mgoB and mgoC, which was suggested by the results of the polarity and mangotoxin production experiments. The start point of the transcript (nucleotide +1) was located 18 nucleotides after the predicted -10 box of the second putative promoter (Figure 3C). Therefore, the second predicted promoter appears to be the functional promoter for the mgo operon. At this point using the known nucleotide sequence and the 5'RACE results, alternative -35 and -10 boxes were located in correct positions from nucleotide +1. The sequences of these alternative -35 and -10 boxes are more typical of Pseudomonas sigma70-dependent promoter sequences [19, 20] than the predicted boxes by BPROM software, which are similar to Escherichia coli sequences (Figure 3C). Additionally, the results do not support the presence of an alternative promoter at the end of mgoB, which could explain the previous results.
The location of the transcriptional terminator was then determined. A 118-bp sequence was located in the region downstream of the mgo operon (Figure 5A) and was compared with the equivalent DNA segment in Pss B728a by Blast (NCBI). A putative terminator (CCC CTC ATC GCG TAA GCG ATG AGG GG), which was 100% identical to the equivalent terminator in Pss B728a, was identified at position 79 from the mgoD stop codon. This terminator sequence was then analysed by FoldRNA software (SoftBerry Inc.), a program used to predict RNA secondary structure through energy minimisation, to calculate the free energy released during palindrome structure formation. A value of -24.4 kcal/mol was found in 84% of the helices. The entire sequence of 118 bp was also analysed by FindTerm software (SoftBerry Inc.) to locate putative Rho-independent bacterial terminators. Two putative terminators (T1 and T2) were found, the first (T1) of which contained more apparent poly-U tracts typical of Rho-independent terminators (Figure 5B, C). T1 was located at position 20-57 (-12.5 kcal/mol and 35% in helices), and T2 was located at position 75-108 (-24.9 kcal/mol and 40% in helices), which includes the sequence homologous to the B728a terminator. Both terminator sequences had negative free energy values, indicating that their folding would be favoured and spontaneous. Finally, to determine which putative terminator acted as the functional terminator, RT-PCR experiments were performed by amplifying the 3'-end of the transcript with primers designed to anneal before, in the middle of and after of the putative terminators (Figure 5D). The amplification test of the mgo transcript revealed that the T1 sequence but not the T2 sequence was included in the mgo transcript, indicating that T1 is the functional terminator of the mgo operon.
The hypothetical function of the mgooperon
Our study of the mgo locus demonstrates that the mgo operon is involved in the biosynthesis or regulation of mangotoxin. Recent studies of the pvf genes, which share high homology with the mgo operon, have indicated a possible regulatory function for those genes . Given these findings, it should be possible to isolate a signalling molecule that is required for virulence gene expression and use it to restore the virulence of an mgoA mutant (defective in the nonribosomal peptide synthetase ) by adding this molecule to the growth medium. Growing the UMAF0158 mutant, which possesses a deletion of mgoA (UMAF0158ΔmgoA) and is defective in mangotoxin production, in media supplemented with an extract from wild-type UMAF0158 restored mangotoxin production. An extract from the mgoA mutant did not restore toxin production. Strains that were defective in other regulatory genes were also used. Extracts from wild-type Pss UMAF0158 and the reference strain Pss B728a were used to complement UMAF0158-2βB7, which contains a miniTn5 disruption of the gacA gene, and UMAF0158-3αE10, which contains a miniTn5 disruption of the gacS gene (Table 4). Mangotoxin production was restored in the defective mutants when an extract from UMAF0158 was added. By contrast, an extract from Pss B728a only restored mangotoxin production in the gacS mutant (Table 4). These results suggest a possible regulatory role for the mgo operon.
The focus of the present study was to characterise the transcriptional organisation that is directly involved in mangotoxin production. We had previously identified the mgo operon (Mangotoxin-Generating Operon) . We determined which genes are involved in mangotoxin production by disrupting each chromosomal gene that was previously identified in pCG2-6 by mutagenesis. The disruption of ORF0 and ORF1 did not affect mangotoxin production. These two genes may belong to another independent gene cluster located close to the mgo operon that is not involved in mangotoxin production. ORF2 transcription was independent of the mgo operon, and ORF2 is homologous to the GntR family of transcriptional regulators. This family of regulatory proteins consists of the N-terminal HTH region of GntR-like bacterial transcription factors. An effector-binding/oligomerisation domain is usually located at the C-terminus . In the deposited genomes of other P. syringae pathovars, the genes in this family are often located close to gene clusters that are homologous to the mgo operon. The relationship between ORF2 and the regulation of the mgo operon remains unclear. In the present study, we observed promoter P mgo expression in the ORF2 mutant (UMAF0158::ORF2) when it was grown in minimal medium at 22°C but not at 28°C, in agreement with the production of mangotoxin by the ORF2 insertional mutant. These data suggest that ORF2 is not involved in mangotoxin production but provide no direct information on the possible influence of ORF2 on the mgo operon with respect to variations in temperature.
Our results demonstrate that the DNA sequence downstream of ORF2 constitutes an operon. Ma et al.  first established the correlation between the presence of a Shine-Dalgarno sequence, also known as a ribosomal binding site (RBS), and translational initiation, the expression levels of the predicted genes and operon structure . We found putative RBSs in almost all of the genes in the putative mgo operon. Only the mgoA gene, in which the start codon overlaps with the stop codon of mgoC, does contain a potential RBS sequence. mgoC and mgoA may share the same RBS, and post-translational changes may separate the two proteins; this situation could explain the absence of a putative RBS for the mgoA gene. The mutagenesis and bioinformatics analysis of each gene in the mgo operon provided insight into their relationship to mangotoxin production. The disruption of mgoB did not abrogate mangotoxin production; however, the production decreased noticeably compared with the wild-type strain. Protein domain searches indicated that mgoB is similar to haem oxygenase. This enzyme is a member of a superfamily represented by a multi-helical structural domain consisting of two structural repeats that is found in both eukaryotic and prokaryotic haem oxygenases and in proteins that enhance the expression of extracellular enzymes . The disruption mutants of the next three genes, mgoC, mgoA and mgoD, were unable to produce mangotoxin, indicating that these genes are essential for mangotoxin production. A similar conclusion was reached by Aguilera et al. , who obtained Tox- phenotypes when 11 different genes in the pht cluster, a region involved in phaseolotoxin production, were mutated by insertion, indicating that all of the genes located within this region encode proteins that are required at different stages of phaseolotoxin production, including synthesis, transport and regulation .
The sequence analysis of mgoC prompted us to search the superfamily protein domains, revealing a similarity to the N-oxygenase domain. This domain was identified in the protein PrnD, which is derived from the pyrrolnitrin biosynthesis gene cluster of Pseudomonas fluorescens. MgoC is also similar to AurF from Streptomyces thioluteus, which produces the starter unit p-nitrobenzoic acid (PNBA) for the polyketide synthase of the aureothin biosynthesis pathway . The gene mgoA, which is homologous to non-ribosomal peptide synthetases, is the largest gene in the mgo operon, and its disruption produces a mutant that is defective in mangotoxin production. Its structure, participation in mangotoxin production and influence on the virulence of the wild-type bacterium has been discussed previously . The final gene studied was mgoD; a domain localisation analysis indicated that mgoD could be a Polyketide_cyc2 belonging to the star-related lipid-transfer (START) domain superfamily. The START superfamily includes bacterial polyketide cyclase/aromatases and two families of previously uncharacterised proteins that are present only in plants and the cyanobacterium Prochlorococcus .
After analysing the elements that composed the putative mgo operon, we evaluated whether the four genes were transcribed together in a single transcript. RT-PCR experiments using the wild-type RNA showed that the four genes were connected in the single transcript (Figure 2). Moreover, the transcript size was analysed by hybridisation, which confirmed the presence of a single transcript with a sufficient size (about 6 kb) to contain the genes mgoBCAD; however, the exact size of the transcript could not be determined.
Following the identification of the mgo operon, the promoter and transcription terminator were identified and studied. The in silico analysis of the sequence identified two putative promoters. Promoter activity was detected only in a minimal medium, the same culture medium that is traditionally used for antimetabolite toxin assays [2, 13]. Promoter activity occurred in the wild-type strain at both temperatures and in the ORF2 insertion mutant at 22°C only. The other Pseudomonas spp. experimental strains, which do not produce mangotoxin, did not exhibit any β-Gal activity. The promoter activity in the wild-type strain was more intense at 28°C than 22°C. When the promoter activity was assayed at 22°C, the activity of the mutant UMAF0158::ORF2 was statistically comparable with that of the wild-type strain. These results suggest a possible influence of ORF2 on the mgo operon during its regulation in response to temperature variations. The promoter inactivity in the other two strains of Pseudomonas spp. may be due to the absence of genes homologous to the mgo operon in P. fluorescens Pf-5, but this explanation is not applicable to Pss B728a. The sequence in B728a that is homologous to the mgo operon is composed of genes that are orthologous to the mgo genes; theoretically, the promoter activity should have been similar to that of the wild-type strain, but it was not. This result suggests that there are additional genes that are necessary for mangotoxin production that are not present in B728a. In support of this explanation, additional genes involved in mangotoxin production have been identified in UMAF0158 and cloned into a different vector than pCG2-6 . The initial sequence analysis did not show any identity with the genome of B728a, and thus these additional genes may influence mgo promoter activity.
Finally, the functional promoter of the mgo operon was established by locating the start of the mgo transcript (Figure 4), which is located 18 nucleotides after the putative -10 box of the second promoter analysed in silico. Thus, the first putative promoter was eliminated as a functional promoter of the mgo operon. Once the +1 site was established, it was possible to locate additional -35 and -10 boxes, which were typical of sigma70 dependent promoters of Pseudomonas spp [19, 20] and were more closely related than the predicted -35 and -10 boxes by BPROM software developed for Escherichia coli, which are less accurate in the search for promoters of Pseudomonas spp. These results allowed us to determine the functional promoter of the mgo operon. The mgo operon terminator was found in a similar manner. The in silico analysis of the sequence identified two possible terminator sequences between the 3'-end of mgoD and the 5'-end of the 5S rRNA, both of which exhibited secondary structures typical of transcription terminators. We considered that the ribosomal transcript terminator is also likely present in the analysed sequence. RT-PCR was used to clarify which was the operon terminator, establishing T1 as the functional terminator of the mgo operon. This is a typical terminator with a stable hairpin having many GC pairs followed by a string of T's. So, it seems that the T1 terminator is a bifunctional terminator, serving this DNA region to terminate transcription of mgo operon in the sense strand and of the ribosomal operon in the antisense strand (Figure 5).
The results described above are sufficient to suggest that mgoBCAD is a transcriptional unit and therefore propose that mgo is an operon. If this argument is correct, mutations in each mgo gene should lead to the absence of a transcript for the downstream genes. A polar effect was demonstrated for UMAF0158::mgoC but not UMAF0158::mgoB. The mutation in mgoB did not prevent the transcription of the downstream genes, although the hybridisation experiments revealed that the transcription appeared to be less efficient. This reduction in transcription corresponds to the reduced production of mangotoxin by UMAF0158::mgoB relative to the wild-type strain. Therefore, the results obtained with wild-type UMAF0158 and the insertional mutants of mgoC, mgoA and mgoD support the hypothesis that the mgo genes form an operon in contrast, the results with the mutant UMAF0158::mgoB do not. We also evaluated the possible existence of an alternative promoter after the mgoB gene, which would explain the production of mangotoxin by the mutant UMAF0158::mgoB. However, during 5'RACE experiment (Figure 3) only a single transcription start site was located, eliminating the possibility of another promoter downstream of mgoB. Therefore there must be something different between the mutant and wild-type strain, which is probably the plasmid integration. In reviewing the process by which the mgo mutants were obtained, we observed that UMAF0158::mgoB was not easy to obtain. The size of mgoB is 777 bp, and the cloned sequence in pCR2.1 was 360 bp of mgoB. The integration of pCR::mgoB into mgoB occurred by single-crossover homologous recombination as it was confirmed. During this process, the plasmid could be integrated into mgoB sequence maintaining an important part of the gene. In this circumstances mgoB or sufficient fragment of it, and the remarkably other three genes of the mgo operon, could be under the influence of a promoter located in plasmid polylinker, lacZ promoter, allowing a reduced transcript expression (Figure 2) and mangotoxin production (Tables 1 and 2). To determine the insert position, a PCR was performed in which the forward primer annealed to the lacZ gene (M13F primer) and the reverse primer annealed to the 5'-end of the mgoC gene, with wild-type UMAF0158 used as the negative control. The amplicon obtained from the mutant UMAF0158::mgoB had a size of 1000 bp, confirming that the plasmid pCR::mgoB was integrated and the lacZ promoter is close to mgoB fragment (Additional file 1: Figure S1).
Because the chemical structure of mangotoxin is unknown , it is difficult to establish a hypothesis concerning the role of the mgo genes in mangotoxin biosynthesis or to determine whether they are related to the regulation of mangotoxin production. Recent studies in P. entomophila have focussed on the pvf gene cluster, which is homologous to the mgo operon, and suggest that the gene cluster serves as a regulator of certain virulence factors in pathogenic strains of Pseudomonas spp. The pvf gene cluster may be a new regulatory system that is specific to certain Pseudomonas species . In the present study, extract complementation restored mangotoxin production in the UMAF0158ΔmgoA mutant only when the culture medium was supplemented with an extract from wild-type UMAF0158. Polar effects of the deleted mgoA on mgoD expression were excluded because the construction of the deletion mutant preserved the reading phase of protein translation. Mangotoxin production was restored in the miniTn5 mutants, which contain disrupted regulatory genes, when their cultures were complemented with a wild-type extract. These results are in agreement with the results obtained by Vallet-Gely et al. , in which pvf and gac mutants were complemented by a wild-type extract. These results allow us to propose a putative regulatory role for the mgo operon in secondary metabolite production by P. syringae pv. syringae, in accordance with Vallet-Gely et al. .
To fully characterise the functions of the mgo operon, more data concerning the chemical structure of mangotoxin and a characterisation of the other genetic traits that regulate mangotoxin biosynthesis by P. syringae pv. syringae UMAF0158 are required.
In the present study, the organisation of the mgo operon in P. syringae pv. syringae UMAF0158 was characterised. The mgo operon is composed of four genes, mgoB, mgoC, mgoA and mgoD. Additionally, this operon possesses one active promoter and a terminator. The last three genes are essential for mangotoxin production, as insertional mutation of these genes results in a loss of mangotoxin production. This operon is only active in minimal medium, in agreement with the standard process for mangotoxin production. Moreover, experiments performed to determine the functional role of the mgo operon demonstrated a putative regulatory function in the production of mangotoxin.
Bacterial strains and plasmids used in this study
The strains of Escherichia coli, Pseudomonas fluorescens Pf-5 and Pseudomonas syringae pv. syringae as well as the vectors and plasmids used in this study are listed in Table 5. E. coli was grown in Luria-Bertani medium (LB) at 37°C for 24 h. The Pseudomonas strains were grown routinely in King's medium B (KB) at 28°C for 48 h. Derivative mutants of P. syringae pv. syringae UMAF0158 (Table 5) were grown and maintained in KB supplemented with the appropriate antibiotics (ampicillin, 50 μg/ml; streptomycin, 50 μg/ml; kanamycin, 50 μg/ml; and gentamicin, 20 μg/ml).
Detection of P. syringaetoxin production
Syringomycin complex production by strains of P. syringae strains was detected using growth inhibition tests performed on potato dextrose agar (PDA) against Geotrichum candidum  and nutrient agar against Rhodotorula pilimanae .
Mangotoxin production was assayed using the indicator technique, which has been described previously and involves growth inhibition of Escherichia coli on Pseudomonas Minimal Medium (PMS ). Briefly, a double layer of the indicator microorganism was made with E. coli CECT831. After solidification, the P. syringae wild-type strain and its derivatives mutants were stabbed, and the plates were incubated at 22°C for 24 h, followed by an additional 24 h at 37°C. To determine the identity of the biochemical step that is putatively targeted by mangotoxin, the same plate bioassay was performed in separate plates with the addition of 100 μl of a 6 mM solution of ornithine or N-acetyl-ornithine. To assess the production of mangotoxin in liquid cultures, we used a cell-free filtrate dilution as previously described .
Insertional inactivation mutagenesis
Insertional inactivation mutagenesis of P. syringae pv. syringae UMAF0158 was used to suppress the production of mangotoxin by inserting disruption vectors into the different ORFs of the mgo operon by single-crossover homologous recombination. To construct the integrative plasmids (Table 4), DNA fragments of the different ORFs within the gene cluster were obtained by PCR using primers specific to the sequence of the genomic clone pCG2-6 (accession number DQ532441) . PCR, cloning and plasmid purification were carried out following standard procedures. The plasmids were transformed into the wild-type strain UMAF0158 by standard electroporation. The mangotoxin-deficient phenotype of the mutants was evaluated by the mangotoxin production assay described previously. Additionally, the mutants were analysed by PCR and Southern blot analyses using the antibiotic resistance cassette or partial target gene sequences as probes to confirm gene disruption and select single-copy transformants.
To prevent potential polar effects from the mutations introduced into the mgo operon, a series of plasmids containing the mutated ORF in addition to each of the downstream ORFs located within the operon was constructed. To ensure expression, these constructs were fused to the PLAC promoter, which is constitutively activated in P. syringae. A fragment containing mgoB, mgoC, mgoA and mgoD (7808 bp) was amplified by PCR from UMAF0158 using the primers ORF3F (5'- CTG CAC AGC CGA CAC TTT TA -3') and ORF6R (5'- TCC GAG GAT CCT GTT GTG GTG CAG CAT CAG TC -3'). A fragment containing mgoA and mgoD (4107 bp) was amplified from P. syringae pv. syringae UMAF0158 using the primers ORF5F (5'- CCG CCG GAT CCC ACT GGT GGC TAA CAT CGT G -3') and ORF6R; both primers contained an artificial BamHI site at the 5' end to facilitate cloning. The amplifications were performed with a high-fidelity Taq polymerase (Expand High Fidelity PCR System, Roche, Basel, Switzerland), and the PCR products were cloned into the vector pGEM-T (Invitrogen, California, USA). The cloned amplicons were removed from the vector by digestion with BamHI and individually cloned into the BamHI site located within pBBR1MCS-5 . The amplicons were cloned in the direction of transcription downstream from the PLAC promoter, resulting in plasmids pLac36 (mgoB, mgoC, mgoA and mgoD) and pLac56 (mgoA and mgoD), which contained the 7.8-kb and the 4.1-kb amplicons, respectively. To obtain mgoD alone, pLac56 was digested with SalI, and the 0.8-kb fragment containing mgoD was recovered and cloned into pBBR1MCS-5, resulting in pLac6. The complementing plasmids were introduced into P. syringae by standard electroporation.
Preparation of RNA for RT-PCR and northern blot experiments
Pure cultures of the wild-type strain of P. syringae pv. syringae UMAF0158 were grown for 48 h at 28°C in KMB agar to prepare a bacterial suspension in PMS minimal medium that possessed an optical density of 1.0 at 600 nm (approximately 109 cfu/ml). One millilitre of this suspension was used to inoculate 50 ml of PMS minimal medium. The culture was incubated at 22°C for 48 h with orbital shaking. RNA was isolated from the bacterial culture with a commercial NucleoSpin RNA Plant kit (Macherey-Nagel GmbH & Co. KG, Germany). The RNA concentration was determined using a Nanodrop ND-1000 (NanoDrop Technologies Wilmington, DE) and was optimised up to 50 ng/μl for RT-PCR assays and 1 μg/μl for Northern blotting. The integrity of the RNA sample was assessed by agarose gel electrophoresis. RT-PCR was performed using 100 ng of RNA at a final volume of 50 μl using the Titan OneTube RT-PCR system, according to the manufacturer's instructions (Roche Diagnostics). The primers were designed by using sequences located between each gene (Additional file 2: Table S1). A 40-cycle amplification programme (94°C for 30 s, 58°C for 1 min, and 68°C for 1 min) was performed followed by a final extension cycle at 68°C for 7 min. Positive control reactions that contained DNA isolated from each corresponding bacterial strain were included in all assays.
Northern blotting was performed using a denaturing agarose gel (0.7%) and formaldehyde (2.2 M). The samples were prepared with 20 μg of total RNA in MOPS running buffer with 2.2. M formaldehyde and 50% formamide and denatured at 65°C for 10 min. The agarose gel was run for 90 min at 60 V. The RNA was transferred to a nylon membrane by capillary diffusion using 10× saline-sodium citrate buffer (SSC) and was immobilised by UV cross-linking. The hybridisation was performed with radioactively labelled probes (dCTP32).
Characterisation of the mgo operon promoter
We used pMP220  as the promoter-probe vector to measure transcriptional activity by β-galactosidase (β-Gal) expression. The amplicon (1008 bp), which included the putative promoter region upstream of mgoB, was cloned into the multicloning site using the EcoRI and PstI restriction sites, which were not present in the cloned sequence. The resulting plasmid, pMPmgo, was transformed into multiple bacterial species (Table 5), and β-Gal assays were performed [17, 18]. The protocol followed the assay described by J.H. Miller , except for the addition of an extra step. In our assay, the cells were pelleted and then resuspended in assay buffer to eliminate any error in the detection of β-galactosidase enzyme activity due to the effects of different carbon sources present in the growth medium.
Additionally, 5'-RACE (Rapid Amplification of cDNA Ends) experiments were performed to locate the +1 nucleotide in the mgo operon transcript and determine which putative promoter is active during mgo operon transcription. The commercial SMART™ RACE cDNA Amplification Kit (Conthech Laboratory, Inc.) was used. Moreover, mRNA from UMAF0158 was obtained by a commercial NucleoSpin RNA Plant kit (Macherey-Nagel GmbH & Co. KG), as described above.
Extracts from wild-type UMAF0158 and the mutant UMAF0158ΔmgoA were used in the complementation experiments. UMAF0158ΔmgoA, which lacks the mgoA gene, was created by deletion mutagenesis as described by Martinez-García and de Lorenzo , taking special care to maintain the reading phase downstream of the mgoD gene. Extractions from the culture supernatant were performed as described by Vallet-Gely et al. . Briefly, 200 ml of bacterial culture in PMS minimal medium was pelleted by centrifugation after 7 days of growth. The supernatants were passed through a 0.2-μm filter (Millipore Corporation, Bedford, MA); the pH was adjusted to 5.0 with HCl or NaOH, and the preparation was extracted three times with dichloromethane. Initially, the preparations were extracted with 100 ml of solvent, then again with 70 ml of solvent and finally with 50 ml of solvent. The extracts were pooled, dried with anhydrous Na2SO4, filtered through Whatman paper, evaporated to dryness and dissolved in 1 ml of methanol. To supplement the growth medium with extract, 150 μl of methanolic extract was added to a 15-ml PMS culture, which was subsequently allowed to grow for 24 h. The mangotoxin production was analysed as previously described, and cell-free filtrates of UMAF0158 and UMAF0158ΔmgoA supplemented with extracts from UMAF0158 and UMAF0158ΔmgoA were tested. Cell-free filtrates from P. syringae pv. syringae UMAF0158 and UMAF0158ΔmgoA grown in PMS supplemented with 150 μl of methanol were used as controls, as were cell-free filtrates of UMAF0158 and UMAF0158ΔmgoA that were grown in PMS under standard conditions.
Database searches were performed using the website of the National Center for Biotechnology Information (NCBI) (http://www.ncbi.nlm.nih.gov). Homology searches and the analysis of conserved protein domains were performed using the NCBI Specialized BLAST programme, the protein tools (InterProScan) of the EMBL European Bioinformatics Institute (http://www.ebi.ac.uk) and the Pfam database (http://pfam.sanger.ac.uk). The restriction maps were constructed and analysed using the JustBio website (http://www.justbio.com). The primers were designed using Primer3 online software (http://primer3.sourceforge.net). The annotation and general manipulation of sequences was performed using Artemis software (Sanger Institute, Cambridge, U.K.). The plasmid maps were constructed using the programme Plasmid Map Enhancer 3.1 (Scientific & Educational Software). The promoter prediction was performed by SoftBerry online software http://linux1.softberry.com/berry.phtml.
Mitchell RE: The relevance of non-host toxins in the expression of virulence by pathogens. Annu Rev Phytopathol. 1984, 22: 215-245. 10.1146/annurev.py.22.090184.001243.
Bender C, Alarcón-Chaidez F, Gross DC: Peudomonas syringa phytotoxins: mode of action, regulation, and biosynthesis by peptide and polyketide synthetases. Microbiol Mol Biol Rev. 1999, 63: 266-292.
Mitchell RE: Implications of toxins in the ecology and evolution of plant pathogenic microorganisms: bacteria. Experientia. 1991, 47: 791-803. 10.1007/BF01922459.
Roth P, Hädener A, Tamm C: Further studies on the biosynthesis of tabtoxin (wildfire toxin): incorporation of [2,3-13C2]pyruvate into the β-lactam moiety. Helv Chim Acta. 1990, 73: 476-482. 10.1002/hlca.19900730228.
Unkefer PJ: The biosynthesis of tabtoxinine-beta-lactam use of specially C-13-labeled glucose and C-13-NMR-spectroscopy to identify its biosynthetic precursors. J Biol Chem. 1987, 262: 4994-4999.
Kinscherf TG, Willis DK: The biosynthetic gene cluster for the b-lactam antibiotic tabtoxin in Pseudomonas syringa. J Antibiot. 2005, 58: 817-821. 10.1038/ja.2005.109.
Tamura K, Imamura M, Yoneyama K, Kohno Y, Takikawa Y, Yamaguchi I, Takahashi H: Role of phaseolotoxin production by Pseudomonas syringapv. actinidain the formation of halo lesions of kiwifruit canker disease. Physiol Mol Plant Pathol. 2002, 60: 207-214. 10.1006/pmpp.2002.0405.
Hernández-Guzmán G, Álvarez-Morales A: Isolation and characterization of the gene coding for the amidinotransferase involved in the biosynthesis of phaseolotoxin in Pseudomonas syringapv. phaseolicol. Mol Plant-Microbe Interact. 2001, 14: 1351-1363. 10.1094/MPMI.2001.14.12.1351.
Zhang YX, Patil SS: The phE locus in the phaseolotoxin gene cluster has ORFs with homologies to genes encoding amino acid transferase, the AraC family of transcriptional factors, and fatty acid desaturases. Mol Plant-Microbe Interact. 1997, 10: 947-960. 10.1094/MPMI.1918.104.22.1687.
Aguilera S, López-López K, Nieto Y, Garcidueñas-Piña R, Hernández-Guzmán G, Hernández-Flores JL, Murillo J, Álvarez-Morales A: Functional characterization of the gene cluster from Pseudomonas syringapv. phaseolicol NPS3121 involved in synthesis of phaseolotoxin. J Bacteriol. 2007, 189: 2834-2843. 10.1128/JB.01845-06.
Kennelly MM, Cazorla FM, de Vicente A, Ramos C, Sundin GM: Pseudomonas syringa diseases of fruit trees. Progress toward understanding and control. Plant Dis. 2007, 91: 4-17. 10.1094/PD-91-0004.
Cazorla FM, Torés JA, Olalla L, Pérez-García A, Farré JM, de Vicente A: Bacterial apical necrosis in mango in southern Spain: a disease produced by Pseudomonas syringapv. syringa. Phytopathology. 1998, 88: 614-620. 10.1094/PHYTO.1922.214.171.1244.
Arrebola E, Cazorla FM, Durán VE, Rivera E, Olea F, Codina JC, Pérez-García A, de Vicente A: Mangotoxin: a novel antimetabolite toxin produced by Pseudomonas syringainhibiting ornitine/arginine biosynthesis. Physiol Mol Plant Pathol. 2003, 63: 117-127. 10.1016/j.pmpp.2003.11.003.
Arrebola E, Cazorla FM, Codina JC, Gutierrez-Barranquero JA, Pérez-García A, de Vicente A: Contribution of mangotoxin to the virulence and epiphytic fitness of Pseudomonas syringapv. syringa. Int Microbiol. 2009, 12: 87-95.
Arrebola E, Cazorla FM, Romero D, Pérez-García A, de Vicente A: A nonribosomal peptide synthetase gene (mgA) of Pseudomonas syringapv. syringa is involved in mangotoxin biosynthesis and is required for full virulence. Mol Plant-Microbe Interact. 2007, 20: 500-509. 10.1094/MPMI-20-5-0500.
Feil H, Feil W, Chain P, Larimer F, DiBartolo G, Copeland A, Lykidis A, Trong S, Nolan M, Goltsman E, Thiel J, Malfatti S, Loper JE, Lapidus A, Detter JC, Land M, Richardson PM, Kyrpides NC, Ivanova N, Lindow SE: Comparison of the complete genome sequences of Pseudomonas syringapv. syringaB728a and pvtomatDC3000. PNAS. 2005, 102: 11064-11069. 10.1073/pnas.0504930102.
Solaiman DKY, Swigle BM: Isolation of novel Pseudomonas syringapromoters and functional characterization in polyhydroxyalkanoate-producing pseudomonads. New Biotechnol. 2010, 27: 1-9. 10.1016/j.nbt.2009.12.003.
Miller JH: Experiments in Molecular Genetics. 1972, NY: Cold Spring Harbor Laboratory, 352-355.
Ramos JL: Pseudomonas: Virulence and Gene Regulator. 2004, NY: Kluwer Academic/Plenum Publishers, ISBN 0-306-48376-9
Humair B, Wackwitz B, Haas D: GacA-controlled activation of promoters for small RNA genes in Pseudomonas fluorescen. Appl Environ Microb. 2010, 76: 1497-1506. 10.1128/AEM.02014-09.
Vallet-Gely I, Opota O, Boniface A, Novikov A, Lemaitre B: A secondary metabolite acting as a signalling molecule controls Pseudomonas entomophi virulence. Cell Microbiol. 2010, doi:10.1111/j.1462-5822.2010.0150.x
Hillerich B, Westpheling J: A new GntR family transcriptional regulator in Streptomyces coelicol is required for morphogenesis and antibiotic production and controls transcriptional of an ABC transporter in response to carbon source. J Bacteriol. 2006, 188: 7477-7487. 10.1128/JB.00898-06.
Ma J, Campbell A, Karlin S: Correlation between Shine-Dalgarno sequences and gene features such as predicted expression levels and operon structures. J Bacteriol. 2002, 184: 5733-5745. 10.1128/JB.184.20.5733-5745.2002.
Wegiel B, Otterbein L: Heme oxygenase 1. UCSD-Nature Mol Pages. 2001, doi:10.1038/mp.a004120.01
Zocher G, Winkler R, Hertweck C, Schulz GE: Structure and action of the N-oxygenase AurF from Streptomyces thiolute. J Mol Biol. 2007, 373: 65-74. 10.1016/j.jmb.2007.06.014.
Iyer LM, Koonin EV, Aravind L: Adaptations of the helix-grip fold for ligand binding and catalysis in the START domain superfamily. Proteins. 2001, 43: 134-144. 10.1002/1097-0134(20010501)43:2<134::AID-PROT1025>3.0.CO;2-I.
Hanahan D: Techniques for transformation on E. col. In DNA Cloning: a Practical Approach. Edited by: Glover DM. 1985, Oxford: IRL, 109-135.
Paulsen IT, Press CM, Ravel J, Kobayashi DY, Myers GS, Mavrodi DV, DeBoy RT, Seshadri R, Ren Q, Madupu R, Dodson RJ, Durkin AS, Brinkac LM, Daugherty SC, Sullivan SA, Rosovitz MJ, Gwinn ML, Zhou L, Schneider DJ, Cartinhour SW, Nelson WC, Weidman J, Watkins K, Tran K, Khouri H, Pierson EA, Pierson LS, Thomashow LS, Loper JE: Complete genome sequence of the plant commensal Pseudomonas fluorescen Pf-5. Nat Biotechnol. 2005, 23: 873-878. 10.1038/nbt1110.
Kovach ME, Elzer PH, Hill DS, Robertson GT, Farris MA, Roop RM, Peterson KM: Four news derivates of the broad-host-range cloning vector pBBR1MCS, carrying different antibiotics-resistance cassettes. Gene. 1995, 166: 175-176. 10.1016/0378-1119(95)00584-1.
Spaink HP, Okker RJH, Wijffelman CA, Pees E, Lugtenberg BJJ: Promoters in the nodulation region of the Rhizobium leguminosaru Sym plasmid pRL1JI. Plant Mol Biol. 1987, 9: 27-39. 10.1007/BF00017984.
Martínez-Garcia E, de Lorenzo V: Transposon-base and plasmid-based genetic tools for editing genomes of gram negatives bacteria. Methods Mol Biol. 2012, 813: 267-283. 10.1007/978-1-61779-412-4_16.
Gross DC, DeVay JE: Production and purification of syringomycin, a phytotoxins produced by Pseudomonas syringa. Physiol Plant Pathol. 1977, 11: 13-28.
Iacobellis NS, Lavermicocca P, Grgurina I, Simmaco M, Ballio A: Phytotoxic properties of Pseudomomas syringa pv. syringa toxins. Physiol Mol Plant Pathol. 1992, 40: 107-116. 10.1016/0885-5765(92)90038-W.
Cazorla FM, Olalla L, Torés JA, Codina JC, Pérez-García A, de Vicente A: Pseudomonas syringae pv. syringae as microorganism involved in apical necrosis of mango: characterization of some virulence factors. Pseudomonas syringae Pathovars and related Species. Edited by: Rudolph K, Burr TJ, Mansfield JW, Stead D, Vivian A, von Kietzell J. 1997, Dordrecht: Kluwer Academic Publishers, 82-87.
This study was supported by funding from Consejería de Innovación, Ciencia y Empresa, Secretaría General de Universidades, Investigación y Tecnología, Junta de Andalucía, Spain (Proyecto de Excelencia P07-AGR-2471), cofinanced by FEDER funds (EU). This work was developed during my hired by the CSIC in the program mode JAEDoc "Junta para la Ampliación de Estudios" cofinanced by ESF.
EA performed the RT-PCR assays, the promoter and terminator characterisations, the mutation experiments and the complementation experiments. EA also performed the mangotoxin test, the evaluation of mangotoxin production using the insertional, deletion and miniTn5 mutants and the Northern blot experiments. JM and EA designed the plasmids and created the constructs used for the complementation experiments. EA also drafted the manuscript. VJC performed the 5'-RACE experiments and the identification of the RBS sites and contributed to the mRNA extraction. FMC and AdV were responsible for initiating this study and participated in its design and coordination and the manuscript preparation. JM conceived the mutation strategy and participated in preparing the final manuscript. APG participated in helpful discussions and the creation of the final manuscript. All authors read and approved the final manuscript.