Comparative Analyses of Four Complete Genomes in Pseudomonas amygdali Revealed Differential Adaptation to Hostile Environments and Secretion Systems

Pseudomonas amygdali is a hemibiotrophic phytopathogen that causes disease in woody and herbaceous plants. Complete genomes of four P. amygdali pathovars were comparatively analyzed to decipher the impact of genomic diversity on host colonization. The pan-genome indicated that 3,928 core genes are conserved among pathovars, while 504–1,009 are unique to specific pathovars. The unique genome contained many mobile elements and exhibited a functional distribution different from the core genome. Genes involved in O-antigen biosynthesis and antimicrobial peptide resistance were significantly enriched for adaptation to hostile environments. While the type III secretion system was distributed in the core genome, unique genomes revealed a different organization of secretion systems as follows: type I in pv. tabaci, type II in pv. japonicus, type IV in pv. morsprunorum, and type VI in pv. lachrymans. These findings provide genetic insight into the dynamic interactions of the bacteria with plant hosts.     

Comparative genomic analysis of Pseudomonas amygdali pv. lachrymans NM002: Insights into its potential virulence genes and putative invasion determinants

Pseudomonas amygdali pv. lachrymans is currently of important plant pathogenic bacteria that causes cucumber angular leaf spot worldwide. The pathogen has been studied for its roles in pathogenicity and plant inheritance resistance. To further delineate traits critical to virulence, invasion and survival in the phyllosphere, we reported the first complete genome of P. amygdali pv. lachrymans NM002. Analysis of the whole genome in comparison with three closely-related representative pathovars of P. syringae identified the conservation of virulence genes, including flagella and chemotaxis, quorum-sensing systems, two-component systems, and lipopolysaccharide and antiphagocytosis. It also revealed differences of invasion determinants, such as type III effectors, phytotoxin (coronatine, syringomycin and phaseolotoxin) and cell wall-degrading enzyme, which may contribute to infectivity. The aim of this study was to derive genomic information that would reveal the probable molecular mechanisms underlying the virulence, infectivity and provide a better understanding of the pathogenesis of the P. syringae pathovars.     

Global genomic analysis of Pseudomonas savastanoi pv. savastanoi plasmids

Pseudomonas savastanoi pv. savastanoi strains harbor native plasmids belonging to the pPT23A plasmid family (PFPs) which are detected in all pathovars of the related species Pseudomonas syringae examined and contribute to the ecological and pathogenic fitness of their host. However, there is a general lack of information about the gene content of P. savastanoi pv. savastanoi plasmids and their role in the interaction of this pathogen with olive plants. We designed a DNA macroarray containing 135 plasmid-borne P. syringae genes to conduct a global genetic analysis of 32 plasmids obtained from 10 P. savastanoi pv. savastanoi strains. Hybridization results revealed that the number of PFPs per strain varied from one to four. Additionally, most strains contained at least one plasmid (designated non-PFP) that did not hybridize to the repA gene of pPT23A. Only three PFPs contained genes involved in the biosynthesis of the virulence factor indole-3-acetic acid (iaaM, iaaH, and iaaL). In contrast, ptz, a gene involved in the biosynthesis of cytokinins, was found in five PFPs and one non-PFP. Genes encoding a type IV secretion system (T4SS), type IVA, were found in both PFPs and non-PFPs; however, type IVB genes were found only on PFPs. Nine plasmids encoded both T4SSs, whereas seven other plasmids carried none of these genes. Most PFPs and non-PFPs hybridized to at least one putative type III secretion system effector gene and to a variety of additional genes encoding known P. syringae virulence factors and one or more insertion sequence transposase genes. These results indicate that non-PFPs may contribute to the virulence and fitness of the P. savastanoi pv. savastanoi host. The overall gene content of P. savastanoi pv. savastanoi plasmids, with their repeated information, mosaic arrangement, and insertion sequences, suggests a possible role in adaptation to a changing environment.     

Characterization of a Unique Chromosomal Copper Resistance Gene Cluster from Xanthomonas campestris pv. vesicatoria

We characterized the copper resistance genes in strain XvP26 of Xanthomonas campestris pv. vesicatoria, which was originally isolated from a pepper plant in Taiwan. The copper resistance genes were localized to a 7,652-bp region which, based on pulsed-field gel electrophoresis and Southern hybridization, was determined to be located on the chromosome. These genes hybridized only weakly, as determined by Southern analysis, to other copper resistance genes in Xanthomonas and Pseudomonas strains. We identified five open reading frames (ORFs) whose products exhibited high levels of amino acid sequence identity to the products of previously reported copper genes. Mutations in ORF1, ORF3, and ORF4 removed copper resistance, whereas mutations in ORF5 resulted in an intermediate copper resistance phenotype and insertions in ORF2 had no effect on resistance conferred to a copper-sensitive recipient in transconjugant tests. Based on sequence analysis, ORF1 was determined to have high levels of identity with the CopR (66%) and PcoR (63%) genes in Pseudomonas syringae pv. tomato and Escherichia coli, respectively. ORF2 and ORF5 had high levels of identity with the PcoS gene in E. coli and the gene encoding a putative copper-containing oxidoreductase signal peptide protein in Sinorhizobium meliloti, respectively. ORF3 and ORF4 exhibited 23% identity to the gene encoding a cation efflux system membrane protein, CzcC, and 62% identity to the gene encoding a putative copper-containing oxidoreductase protein, respectively. The latter two ORFs were determined to be induced following exposure to low concentrations of copper, while addition of Co, Cd, or Zn resulted in no significant induction. PCR analysis of 51 pepper and 34 tomato copper-resistant X. campestris pv. vesicatoria strains collected from several regions in Taiwan between 1987 and 2000 and nine copper-resistant strains from the United States and South America showed that successful amplification of DNA was obtained only for strain XvP26. The organization of this set of copper resistance genes appears to be uncommon, and the set appears to occur rarely in X. campestris pv. vesicatoria.     

The genome sequence of the biocontrol fungus Metarhizium anisopliae and comparative genomics of Metarhizium species

Background

Metarhizium anisopliae is an important fungal biocontrol agent of insect pests of agricultural crops. Genomics can aid the successful commercialization of biopesticides by identification of key genes differentiating closely related species, selection of virulent microbial isolates which are amenable to industrial scale production and formulation and through the reduction of phenotypic variability. The genome of Metarhizium isolate ARSEF23 was recently published as a model for M. anisopliae, however phylogenetic analysis has since re-classified this isolate as M. robertsii. We present a new annotated genome sequence of M. anisopliae (isolate Ma69) and whole genome comparison to M. robertsii (ARSEF23) and M. acridum (CQMa 102).

Results

Whole genome analysis of M. anisopliae indicates significant macrosynteny with M. robertsii but with some large genomic inversions. In comparison to M. acridum, the genome of M. anisopliae shares lower sequence homology. While alignments overall are co-linear, the genome of M. acridum is not contiguous enough to conclusively observe macrosynteny. Mating type gene analysis revealed both MAT1-1 and MAT1-2 genes present in M. anisopliae suggesting putative homothallism, despite having no known teleomorph, in contrast with the putatively heterothallic M. acridum isolate CQMa 102 (MAT1-2) and M. robertsii isolate ARSEF23 (altered MAT1-1). Repetitive DNA and RIP analysis revealed M. acridum to have twice the repetitive content of the other two species and M. anisopliae to be five times more RIP affected than M. robertsii. We also present an initial bioinformatic survey of candidate pathogenicity genes in M. anisopliae.

Conclusions

The annotated genome of M. anisopliae is an important resource for the identification of virulence genes specific to M. anisopliae and development of species- and strain- specific assays. New insight into the possibility of homothallism and RIP affectedness has important implications for the development of M. anisopliae as a biopesticide as it may indicate the potential for greater inherent diversity in this species than the other species. This could present opportunities to select isolates with unique combinations of pathogenicity factors, or it may point to instability in the species, a negative attribute in a biopesticide.

Electronic supplementary material

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

Genomic comparative analysis of the environmental Enterococcus mundtii against enterococcal representative species

Background

Enterococcus mundtii is a yellow-pigmented microorganism rarely found in human infections. The draft genome sequence of E. mundtii was recently announced. Its genome encodes at least 2,589 genes and 57 RNAs, and 4 putative genomic islands have been detected. The objective of this study was to compare the genetic content of E. mundtii with respect to other enterococcal species and, more specifically, to identify genes coding for putative virulence traits present in enterococcal opportunistic pathogens.

Results

An in-depth mining of the annotated genome was performed in order to uncover the unique properties of this microorganism, which allowed us to detect a gene encoding the antimicrobial peptide mundticin among other relevant features. Moreover, in this study a comparative genomic analysis against commensal and pathogenic enterococcal species, for which genomic sequences have been released, was conducted for the first time. Furthermore, our study reveals significant similarities in gene content between this environmental isolate and the selected enterococci strains (sharing an “enterococcal gene core” of 805 CDS), which contributes to understand the persistence of this genus in different niches and also improves our knowledge about the genetics of this diverse group of microorganisms that includes environmental, commensal and opportunistic pathogens.

Conclusion

Although E. mundtii CRL1656 is phylogenetically closer to E. faecium, frequently responsible of nosocomial infections, this strain does not encode the most relevant relevant virulence factors found in the enterococcal clinical isolates and bioinformatic predictions indicate that it possesses the lowest number of putative pathogenic genes among the most representative enterococcal species. Accordingly, infection assays using the Galleria mellonella model confirmed its low virulence.

Electronic supplementary material

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

Genomic Survey of Pathogenicity Determinants and VNTR Markers in the Cassava Bacterial Pathogen Xanthomonas axonopodis pv. Manihotis Strain CIO151

Xanthomonas axonopodis pv. manihotis (Xam) is the causal agent of bacterial blight of cassava, which is among the main components of human diet in Africa and South America. Current information about the molecular pathogenicity factors involved in the infection process of this organism is limited. Previous studies in other bacteria in this genus suggest that advanced draft genome sequences are valuable resources for molecular studies on their interaction with plants and could provide valuable tools for diagnostics and detection. Here we have generated the first manually annotated high-quality draft genome sequence of Xam strain CIO151. Its genomic structure is similar to that of other xanthomonads, especially Xanthomonas euvesicatoria and Xanthomonas citri pv. citri species. Several putative pathogenicity factors were identified, including type III effectors, cell wall-degrading enzymes and clusters encoding protein secretion systems. Specific characteristics in this genome include changes in the xanthomonadin cluster that could explain the lack of typical yellow color in all strains of this pathovar and the presence of 50 regions in the genome with atypical nucleotide composition. The genome sequence was used to predict and evaluate 22 variable number of tandem repeat (VNTR) loci that were subsequently demonstrated as polymorphic in representative Xam strains. Our results demonstrate that Xanthomonas axonopodis pv. manihotis strain CIO151 possesses ten clusters of pathogenicity factors conserved within the genus Xanthomonas. We report 126 genes that are potentially unique to Xam, as well as potential horizontal transfer events in the history of the genome. The relation of these regions with virulence and pathogenicity could explain several aspects of the biology of this pathogen, including its ability to colonize both vascular and non-vascular tissues of cassava plants. A set of 16 robust, polymorphic VNTR loci will be useful to develop a multi-locus VNTR analysis scheme for epidemiological surveillance of this disease.     

HopAZ1, a type III effector of Pseudomonas amygdali pv. tabaci,induces a hypersensitive response in tobacco wildfire‐resistant Nicotiana tabacum ‘N509’

Pseudomonas amygdali pv. tabaci (formerly Pseudomonas syringae pv. tabaci; Pta) is a gram‐negative bacterium that causes bacterial wildfire disease in Nicotiana tabacum. The pathogen establishes infections by using a type III secretion system to inject type III effector proteins (T3Es) into cells, thereby interfering with the host__s immune system. To counteract the effectors, plants have evolved disease‐resistance genes and mechanisms to induce strong resistance on effector recognition. By screening a series of Pta T3E‐deficient mutants, we have identified HopAZ1 as the T3E that induces disease resistance in N. tabacum ‘N509’. Inoculation with the Pta ∆hopAZ1 mutant did not induce resistance to Pta in N509. We also found that the Pta ∆hopAZ1 mutant did not induce a hypersensitive response and promoted severe disease symptoms in N509. Furthermore, a C‐terminal truncated HopAZ1 abolished HopAZ1‐dependent cell death in N509. These results indicate that HopAZ1 is the avirulence factor that induces resistance to Pta by N509.     

Insights into the genome of the xanthan-producing phytopathogen Xanthomonas arboricola pv. pruni 109 by comparative genomic hybridization

The phytopathogenic bacterium Xanthomonas arboricola pv. pruni is the causal agent of Prunus Bacterial Spot disease that infects cultivated Prunus species and their hybrids. Furthermore, X. arboricola pv. pruni (Xap) plays a role in biotechnology since it produces xanthan gum, an important biopolymer used mainly in the food, oil, and cosmetics industry. To gain first insights into the genome composition of this pathovar, genomic DNA of X. arboricola pv. pruni strains was compared to the genomes of reference strains X. campestris pv. campestris B100 (Xcc B100) and X. campestris pv. vesicatoria 85-10 (Xcv 85-10) applying microarray-based comparative genomic hybridizations (CGH). The results implied that X. arboricola pv. pruni 109 lacks 6.67% and 5.21% of the genes present in the reference strains Xcc B100 and Xcv 85-10, respectively. Most of the missing genes were found to be organized in clusters and do not belong to the core genome of the two reference strains. Often they encode mobile genetic elements. Furthermore, the absence of gene clusters coding for the lipopolysaccharide (LPS) O-antigens of Xcc B100 and Xcv 85-10 indicates that the structure of the O-antigen of X. arboricola pv. pruni 109 differs from that of Xcc B100 and Xcv 85-10.     

Genome-Wide Relatedness of Treponema pedis,from Gingiva and Necrotic Skin Lesions of Pigs,with the Human Oral Pathogen Treponema denticola

Treponema pedis and T. denticola are two genetically related species with different origins of isolation. Treponema denticola is part of the human oral microbiota and is associated with periodontitis while T. pedis has been isolated from skin lesions in animals, e.g., digital dermatitis in cattle and necrotic ulcers in pigs. Although multiple Treponema phylotypes may exist in ulcerative lesions in pigs, T. pedis appears to be a predominant spirochete in these lesions. Treponema pedis can also be present in pig gingiva. In this study, we determined the complete genome sequence of T. pedis strain T A4, isolated from a porcine necrotic ear lesion, and compared its genome with that of T. denticola. Most genes in T. pedis were homologous to those in T. denticola and the two species were similar in general genomic features such as size, G+C content, and number of genes. In addition, many homologues of specific virulence-related genes in T. denticola were found in T. pedis. Comparing a selected pair of strains will usually not give a complete picture of the relatedness between two species. We therefore complemented the analysis with draft genomes from six T. pedis isolates, originating from gingiva and necrotic ulcers in pigs, and from twelve T. denticola strains. Each strain carried a considerable amount of accessory genetic material, of which a large part was strain specific. There was also extensive sequence variability in putative virulence-related genes between strains belonging to the same species. Signs of lateral gene-transfer events from bacteria known to colonize oral environments were found. This suggests that the oral cavity is an important habitat for T. pedis. In summary, we found extensive genomic similarities between T. pedis and T. denticola but also large variability within each species.     

De novo Genome Assembly of the Fungal Plant Pathogen Pyrenophora semeniperda

Pyrenophora semeniperda (anamorph Drechslera campulata) is a necrotrophic fungal seed pathogen that has a wide host range within the Poaceae. One of its hosts is cheatgrass (Bromus tectorum), a species exotic to the United States that has invaded natural ecosystems of the Intermountain West. As a natural pathogen of cheatgrass, P. semeniperda has potential as a biocontrol agent due to its effectiveness at killing seeds within the seed bank; however, few genetic resources exist for the fungus. Here, the genome of P. semeniperda isolate assembled from sequence reads of 454 pyrosequencing is presented. The total assembly is 32.5 Mb and includes 11,453 gene models encoding putative proteins larger than 24 amino acids. The models represent a variety of putative genes that are involved in pathogenic pathways typically found in necrotrophic fungi. In addition, extensive rearrangements, including inter- and intrachromosomal rearrangements, were found when the P. semeniperda genome was compared to P. tritici-repentis, a related fungal species.     

Identifying pathogenicity genes in the rubber tree anthracnose fungus Colletotrichum gloeosporioides through random insertional mutagenesis

To gain more insight into the molecular mechanisms of Colletotrichum gloeosporioides pathogenesis, Agrobacterium tumefaciens-mediated transformation (ATMT) was used to identify mutants of C. gloeosporioides impaired in pathogenicity. An ATMT library of 4128 C. gloeosporioides transformants was generated. Transformants were screened for defects in pathogenicity with a detached copper brown leaf assay. 32 mutants showing reproducible pathogenicity defects were obtained. Southern blot analysis showed 60.4% of the transformants had single-site T-DNA integrations. 16 Genomic sequences flanking T-DNA were recovered from mutants by thermal asymmetric interlaced PCR, and were used to isolate the tagged genes from the genome sequence of wild-type C. gloeosporioides by Basic Local Alignment Search Tool searches against the local genome database of the wild-type C. gloeosporioides. One potential pathogenicity genes encoded calcium-translocating P-type ATPase. Six potential pathogenicity genes had no known homologs in filamentous fungi and were likely to be novel fungal virulence factors. Two putative genes encoded Glycosyltransferase family 28 domain-containing protein and Mov34/MPN/PAD-1 family protein, respectively. Five potential pathogenicity genes had putative function matched with putative protein of other Colletotrichum species. Two known C. gloeosporioides pathogenicity genes were also identified, the encoding Glomerella cingulata hard-surface induced protein and C. gloeosporioides regulatory subunit of protein kinase A gene involved in cAMP-dependent PKA signal transduction pathway.     

Genome characterization of a novel Burkholderia cepacia complex genomovar isolated from dieback affected mango orchards

We characterized the genome of the antibiotic resistant, caseinolytic and non-hemolytic Burkholderia sp. strain TJI49, isolated from mango trees (Mangifera indica L.) with dieback disease. This isolate produced severe disease symptoms on the indicator plants. Next generation DNA sequencing and short-read assembly generated the 60X deep 7,631,934 nucleotide draft genome of Burkholderia sp. TJI49 which comprised three chromosomes and at least one mega plasmid. Genome annotation studies revealed a total 8,992 genes, out of which 8,940 were protein coding genes. Comparative genomics and phylogenetics identified Burkholderia sp. TJI49 as a distinct species of Burkholderia cepacia complex (BCC), closely related to B. multivorans ATCC17616. Genome-wide sequence alignment of this isolate with replicons of BCC members showed conservation of core function genes but considerable variations in accessory genes. Subsystem-based gene annotation identified the active presence of wide spread colonization island and type VI secretion system in Burkholderia sp. TJI49. Sequence comparisons revealed (a) 28 novel ORFs that have no database matches and (b) 23 ORFs with orthologues in species other than Burkholderia, indicating horizontal gene transfer events. Fold recognition of novel ORFs identified genes encoding pertactin autotransporter-like proteins (a constituent of type V secretion system) and Hap adhesion-like proteins (involved in cell–cell adhesion) in the genome of Burkholderia sp. TJI49. The genomic characterization of this isolate provided additional information related to the ‘pan-genome’ of Burkholderia species.     

A comparative genomic analysis of Xanthomonas arboricola pv. juglandis strains reveal hallmarks of mobile genetic elements in the adaptation and accelerated evolution of virulence

Xanthomonas arboricola pv. juglandis (Xaj) is the most significant aboveground walnut bacterial pathogen. Disease management uses copper-based pesticides which induce pathogen resistance. We examined the genetic repertoire associated with adaptation and virulence evolution in Xaj. Comparative genomics of 32 Xaj strains reveal the possible acquisition and propagation of virulence factors via insertion sequences (IS). Fine-scale annotation revealed a Tn3 transposon (TnXaj417) encoding copper resistance genes acquired by horizontal gene transfer and associated with adaptation and tolerance to metal-based pesticides commonly used to manage pathogens in orchard ecosystems. Phylogenomic analysis reveals IS involvement in acquisition and diversification of type III effector proteins ranging from two to eight in non-pathogenic strains, 16 to 20 in pathogenic strains, besides six other putative effectors with a reduced identity degree found mostly among pathogenic strains. Yersiniabactin, xopK, xopAI, and antibiotic resistance genes are also located near ISs or inside genomic islands and structures resembling composite transposons.     

Characterization of the core and accessory genomes of Pseudomonas aeruginosa using bioinformatic tools Spine and AGEnt

Background

Pseudomonas aeruginosa is an important opportunistic pathogen responsible for many infections in hospitalized and immunocompromised patients. Previous reports estimated that approximately 10% of its 6.6 Mbp genome varies from strain to strain and is therefore referred to as “accessory genome”. Elements within the accessory genome of P. aeruginosa have been associated with differences in virulence and antibiotic resistance. As whole genome sequencing of bacterial strains becomes more widespread and cost-effective, methods to quickly and reliably identify accessory genomic elements in newly sequenced P. aeruginosa genomes will be needed.

Results

We developed a bioinformatic method for identifying the accessory genome of P. aeruginosa. First, the core genome was determined based on sequence conserved among the completed genomes of twelve reference strains using Spine, a software program developed for this purpose. The core genome was 5.84 Mbp in size and contained 5,316 coding sequences. We then developed an in silico genome subtraction program named AGEnt to filter out core genomic sequences from P. aeruginosa whole genomes to identify accessory genomic sequences of these reference strains. This analysis determined that the accessory genome of P. aeruginosa ranged from 6.9-18.0% of the total genome, was enriched for genes associated with mobile elements, and was comprised of a majority of genes with unknown or unclear function. Using these genomes, we showed that AGEnt performed well compared to other publically available programs designed to detect accessory genomic elements. We then demonstrated the utility of the AGEnt program by applying it to the draft genomes of two previously unsequenced P. aeruginosa strains, PA99 and PA103.

Conclusions

The P. aeruginosa genome is rich in accessory genetic material. The AGEnt program accurately identified the accessory genomes of newly sequenced P. aeruginosa strains, even when draft genomes were used. As P. aeruginosa genomes become available at an increasingly rapid pace, this program will be useful in cataloging the expanding accessory genome of this bacterium and in discerning correlations between phenotype and accessory genome makeup. The combination of Spine and AGEnt should be useful in defining the accessory genomes of other bacterial species as well.

Electronic supplementary material

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

Strain-specific genes of Helicobacter pylori: genome evolution driven by a novel type IV secretion system and genomic island transfer

The availability of multiple bacterial genome sequences has revealed a surprising extent of variability among strains of the same species. The human gastric pathogen Helicobacter pylori is known as one of the most genetically diverse species. We have compared the genome sequence of the duodenal ulcer strain P12 and six other H. pylori genomes to elucidate the genetic repertoire and genome evolution mechanisms of this species. In agreement with previous findings, we estimate that the core genome comprises about 1200 genes and that H. pylori possesses an open pan-genome. Strain-specific genes are preferentially located at potential genome rearrangement sites or in distinct plasticity zones, suggesting two different mechanisms of genome evolution. The P12 genome contains three plasticity zones, two of which encode type IV secretion systems and have typical features of genomic islands. We demonstrate for the first time that one of these islands is capable of self-excision and horizontal transfer by a conjugative process. We also show that excision is mediated by a protein of the XerD family of tyrosine recombinases. Thus, in addition to its natural transformation competence, conjugative transfer of genomic islands has to be considered as an important source of genetic diversity in H. pylori.     

Pseudomonas syringae pv. actinidiae draft genomes comparison reveal strain-specific features involved in adaptation and virulence to Actinidia species

A recent re-emerging bacterial canker disease incited by Pseudomonas syringae pv. actinidiae (Psa) is causing severe economic losses to Actinidia chinensis and A. deliciosa cultivations in southern Europe, New Zealand, Chile and South Korea. Little is known about the genetic features of this pathovar. We generated genome-wide Illumina sequence data from two Psa strains causing outbreaks of bacterial canker on the A. deliciosa cv. Hayward in Japan (J-Psa, type-strain of the pathovar) and in Italy (I-Psa) in 1984 and 1992, respectively as well as from a Psa strain (I2-Psa) isolated at the beginning of the recent epidemic on A. chinensis cv. Hort16A in Italy. All strains were isolated from typical leaf spot symptoms. The phylogenetic relationships revealed that Psa is more closely related to P. s. pv. theae than to P. avellanae within genomospecies 8. Comparative genomic analyses revealed both relevant intrapathovar variations and putative pathovar-specific genomic regions in Psa. The genomic sequences of J-Psa and I-Psa were very similar. Conversely, the I2-Psa genome encodes four additional effector protein genes, lacks a 50 kb plasmid and the phaseolotoxin gene cluster, argK-tox but has acquired a 160 kb plasmid and putative prophage sequences. Several lines of evidence from the analysis of the genome sequences support the hypothesis that this strain did not evolve from the Psa population that caused the epidemics in 1984-1992 in Japan and Italy but rather is the product of a recent independent evolution of the pathovar actinidiae for infecting Actinidia spp. All Psa strains share the genetic potential for copper resistance, antibiotic detoxification, high affinity iron acquisition and detoxification of nitric oxide of plant origin. Similar to other sequenced phytopathogenic pseudomonads associated with woody plant species, the Psa strains isolated from leaves also display a set of genes involved in the catabolism of plant-derived aromatic compounds.     

High-quality genome assembly of Pseudocercospora ulei the main threat to natural rubber trees

Pseudocercospora ulei is the causal agent of South American Leaf Blight (SALB), the main disease affecting Hevea brasiliensis rubber tree, a native species to the Amazon. Rubber tree is a major crop in South American countries and SALB disease control strategies would benefit from the availability of genomic resources for the fungal pathogen. Here, we assembled and annotated the P. ulei genome. Shotgun sequencing was performed using second and third generation sequencing technologies. We present the first P. ulei high-quality genome assembly, the largest among Mycosphaerellaceae, with 93.8 Mbp, comprising 215 scaffolds, an N50 of 2.8 Mbp and a BUSCO gene completeness of 97.5%. We identified 12,745 protein-coding gene models in the P. ulei genome with 756 genes encoding secreted proteins and 113 genes encoding effector candidates. Most of the genome (80%) is composed of repetitive elements dominated by retrotransposons of the Gypsy superfamily. P. ulei has the largest genome size among Mycosphaerellaceae, with the highest TE content. In conclusion, we have established essential genomic resources for a wide range of studies on P. ulei and related species.