Complete genome sequence of Pseudomonas fluorescens strain PICF7, an indigenous root endophyte from olive (Olea europaea L.) and effective biocontrol agent against Verticillium dahliae

Pseudomonas fluorescens strain PICF7 is a native endophyte of olive roots. Previous studies have shown this motile, Gram-negative, non-sporulating bacterium is an effective biocontrol agent against the soil-borne fungus Verticillium dahliae, the causal agent of one of the most devastating diseases for olive (Olea europaea L.) cultivation. Here, we announce and describe the complete genome sequence of Pseudomonas fluorescens strain PICF7 consisting of a circular chromosome of 6,136,735 bp that encodes 5,567 protein-coding genes and 88 RNA-only encoding genes. Genome analysis revealed genes predicting factors such as secretion systems, siderophores, detoxifying compounds or volatile components. Further analysis of the genome sequence of PICF7 will help in gaining insights into biocontrol and endophytism.


Introduction
Pseudomonas fluorescens PICF7 is a native colonizer of olive (Olea europaea L.) roots and an in vitro antagonist of the soil-borne fungal phytopathogen Verticillium dahliae Kleb. [1], the causal agent of Verticillium wilts in a large number of plant species [2]. This strain has been demonstrated to be an effective BCA against Verticillium wilt of olive [1,3], one of the most important biotic constraints for olive cultivation [4]. Moreover, strain PICF7 is able to display an endophytic lifestyle within olive root tissues under different experimental conditions [3,5,6] and induces a broad range of defence responses at both local (roots) and systemic (above-ground organs) level, as well as to activate diverse transcription factors known to be involved in systemic defence responses [7,8]. Accordingly, a recent study has shown the ability of PICF7 to influence the establishment of the pathogen Pseudomonas savastanoi pv. savastanoi in olive stems and to affect the normal development of olive knots [9], its associated disease [10].
In this report, we summarize the complete genome sequence and annotation of PICF7. We also describe its genomic properties, highlighting genes encoding plantassociated factors, colonization abilities and well-known bacterial biocontrol traits. The genome sequencing of PICF7 and its comparison with related published genomes will provide a framework for further functional studies of its rhizosphere competence, biocontrol effectiveness and endophytic lifestyle.
Classification and features P. fluorescens PICF7 is a motile, Gram-negative, nonsporulating rod in the order Pseudomonadales of the class Gammaproteobacteria. Rod-shaped cells are approximately 0.5 μm in width and 2.0-2.5 μm in length ( Figure 1 Left and Centre). The strain is moderately fast-growing, forming 2 mm colonies within 2-3 days at 28°C. Colonies formed on King's B (KB) [11] agar plates are yellow-green opaque, domed and moderately mucoid with smooth margins (Figure 1 Right).
PICF7 was isolated from the roots of healthy nurseryproduced olive plants cv. Picual in Córdoba province (Southern Spain) [1]. It grows in complex media such as LB [12] or KB, as well as in minimal media such as Standard Succinate Medium (SSM; pH 7.0) [13]. Even though the optimal growth temperature is 28°C, PICF7 can also slightly replicate at 5°C in liquid LB and KB. However, growth at 37°C was not observed in these culturing media after 24 h. The bacterium is an efficient colonizer of the olive rhizosphere [1] and displays an endophytic lifestyle [3,5,6]. It does not cause any deleterious effect on its original host (olive) [1,5,9]. Strain PICF7 has natural resistance to kanamycin (50 mg/L) and nalidixic acid (25 mg/L), and it is possible to develop spontaneous rifampicin-resistant mutants [1].
Minimun Information about the Genome Sequence (MIGS) of P. fluorescens PICF7 is summarized in Table 1, and its phylogenetic position is shown in Figure 2.

Genome sequencing and annotation
Genome project history P. fluorescens strain PICF7 was selected for sequencing due to its ability to exert biocontrol against Verticillium wilt of olive [1,3] and to develop an endophytic lifestyle within olive root tissues [5,6]. The genome project is deposited in the Genomes OnLine Database [31] and the NCBI BioProject database. The finished genome sequence is in GenBank. A summary of the project information is shown in Table 2.
Growth conditions and DNA isolation P. fluorescens strain PICF7 was grown in 50 ml of LB medium and incubated for 16 h at 28°C. After this period of time, the OD600 of the culture was 1.2. Serial dilutions from this culture and plating on LB plates yielded 2.8 × 10 8 CFU/mL of a pure bacterial culture (colonies showed uniform morphology and kanamycin resistance). The culture was divided into two 25-ml aliquots and total genomic DNA was extracted using the 'Jet-Flex genomic DNA purification' kit (Genomed GmbH, Löhne, Germany), according to the manufacturer's indications. DNA samples were further purified by extraction with phenol:chloroform and precipitation with ethanol. DNA quality and quantity were checked by agarose gel electrophoresis, spectrophotometry using a ND-1000 spectrophotometer (NanoDrop Technologies, Wilmington, DE), and digestion with different restriction enzymes. Two DNA aliquots (0.6 μg/μL,~20 μg each) were sent in a dry ice container to the sequencing service.

Genome sequencing and assembly
The genome of PICF7 was sequenced at the Beijing Genomics Institute (BGI) using Solexa paired-end sequencing. Draft assemblies were based on 3,482,351 reads with a length of 500 bp resulting in 1,200 Mb, 2,456,221 reads with a length of 2,000 bp resulting in 1,209 Mb and 1,924,515 reads with a length of 6,000 bp resulting in 1,309 Mb. The SOAPdenovo 1.05 software package [32][33][34] developed by BGI was used for sequence assembly and quality assessment.

Genome annotation
Automatic annotation was performed using the NCBI Prokaryotic Genome Annotation Pipeline. Identification of known type III effectors effectors was conducted by BLASTP searches of the effectors described in http:// pseudomonas-syringae.org/ against the proteome of PICF7. Functional annotation was performed by aligning the predicted protein sequences against the COG PSSM of the CDD using RPS-BLAST. Hits with an E-value < = 0.001 were first retained. Then, only the best hit was selected for each protein. Signal peptides and transmembrane helices were predicted using SignalP [35,36] and TMHMM [37,38], respectively. Figure 1 Image of P. fluorescens PICF7 cells using confocal laser scanning (Left) and phase-contrast (Centre) microscopy (CLSM and PCM, respectively) and the appearance of colony morphology after 24 h growing on KB agar medium at 28°C (Right). CLSM image was obtained using a PICF7 derivative carrying a plasmid with an enhanced green fluorescent protein (EGFP) [5].

Genome properties
The genome of PICF7 is composed of one circular chromosome of 6,136,735 bp with an average GC content of 60.4% (Table 3 and Figure 3), which is similar to that of other P. fluorescens strains. Among the 5,655 predicted genes, 5,567 were identified as protein coding genes. Of the last, 4,573 (82.1%) were assigned a putative function, while the other 994 (17.9%) were designated as hypothetical proteins. The classification of CDSs into functional categories according to the COG (Clusters of Orthologous Groups) [39,40] database is summarized in Table 4.

Insights from the genome sequence
The genome contains a complete canonical type III secretion system and two known effector proteins, namely, AvrE1 and HopB1. In addition, two complete type VI secretion system (T6SS) clusters were identified. T6SS has been described to promote antibacterial activity against a wide range of competitor bacteria [41]. PICF7 genome also encodes gene clusters for the synthesis of the siderophores pyochelin and pyoverdine and the hemophore HasAp. A repertoire of cell adhesion proteins has been also identified, including two filamentous hemagglutinin proteins and several fimbrial proteins clustered together with a number of pilus assembly proteins. Notably, two genes have been found to show high similarity with attC and attG genes from Agrobacterium, whose mutation leads to lack of attachment on tomato, carrot, and Bryophyllum daigremontiana [42]. It is worth mentioning the presence of genome components presumably involved in the synthesis of detoxifying compounds. Such is the case of two clusters containing genes for copper resistance and for production of a cbb (3)-type cytochrome C oxidase, respectively. An ortholog of the gene that codes for Dps, a ferritin-like protein reported to protect plant-associated bacteria against oxidative stress [43], has also been found. Additional identified traits involved in detoxification are orthologs of catalase KatB and hydroperoxidase KatG, which detoxify plant- Figure 2 Phylogenetic tree highlighting the position of P. fluorescens strain PICF7 relative to its closest Pseudomonas strains for which complete genomes are available. P. entomophila strain L48 was used as an outgroup. For the construction of the tree, five protein-coding house-keeping genes were first aligned, namely: argF, atpA, nusA, pyrH and rpoH. Then, Maximum Likelihood method based on the JTT (Jones-Taylor-Thornton) matrix-based model [29] was used. The percentage of trees in which the associated taxa clustered in the bootstrap test (1000 replicates) is shown next to the branches [30].   [44], and a gene coding for a proline iminopeptidase, which has been shown to have dealanylating activity toward the antibiotic ascamycin [45]. A gene predicting a salycilic hydroxylase has been also identified in PICF7 genome. This gene could be involved in the degradation of the plant defence hormone salicylic acid, thus disrupting the systemic response against colonizing bacteria. In addition, all genes required for biosynthesis of the exopolysaccharide alginate [46] are present in a gene cluster.
Genes predicting volatile components are present in PICF7 genome as well. Volatile components have been shown to act as antibiotics and to induce plant growth [47,48]. An example is hydrogen cyanide (HCN), an inorganic compound with antagonistic effects against soil microbes [49]. Orthologs of genes required for Figure 3 Graphical map of the chromosome. From outside to the centre: genes on forward strand (coloured by COG categories), genes on reverse strand (coloured by COG categories), RNA genes: tRNAs -blue, rRNAspink, G + C in relation to the mean G + C in 2 kb windows and trinucleotide distribution in 2 kb windows. The latter was defined as the χ2 statistic on the difference between the trinucleotide composition of 2 kb windows and that of the whole chromosome.
the biosynthesis of other volatile components such as 2,3-butanediol and acetoin were also found. Further genome analysis revealed other factors presumably involved in the endophytic fitness of PICF7. Such is the case of enzymes like a cellulase and a phytase, as well as the gene coding for aminocyclopropane-1-carboxylate deaminase suggested to be key in the modulation of ethylene levels in plants by bacteria [50].

Conclusions
In this report we describe the complete genome sequence of Pseudomonas fluorescens strain PICF7, a "Pseudomonadales" in the order Gammaproteobacteria that was originally isolated from the roots of healthy nursery-produced olive plants cv. Picual in Córdoba province, Spain. This strain was selected for sequencing based on its ability to exert biocontrol against Verticillium wilt of olive and to develop an endophytic lifestyle within olive root tissues. Such properties likely have origins in a repertoire of genes including a putative T3SS, two putative T6SS, and several genes presumably implicated in siderophore production. It also has a collection of genes predicting adhesion proteins, detoxifying compounds, volatile components and enzymes such as a cellulase, aphytase and a deaminase. Further functional studies and comparative genomics with related isolates will provide insights into biocontrol and endophytism.

Competing interests
The authors declare that they have no competing interests.
Authors' contributions JMB conceived and designed the experiments. DRR, ES and PP performed the experiments. PMMG performed the annotation and sequence homology searches. PMMG, CR, JMB and PRP wrote the manuscript. All authors commented on the manuscript before submission. All authors read and approved the final manuscript. The total is based on the total number of protein coding genes in the annotated genome.