RNA-Seq profile of flavescence dorée phytoplasma in grapevine.

BACKGROUND: The phytoplasma-borne disease flavescence dorée is still a threat to European viticulture, despite mandatory control measures and prophylaxis against the leafhopper vector. Given the economic importance of grapevine, it is essential to find alternative strategies to contain the spread, in order to possibly reduce the current use of harmful insecticides. Further studies of the pathogen, the vector and the mechanisms of phytoplasma-host interactions could improve our understanding of the disease. In this work, RNA-Seq technology followed by three de novo assembly strategies was used to provide the first comprehensive transcriptomics landscape of flavescence dorée phytoplasma (FD) infecting field-grown Vitis vinifera leaves. RESULTS: With an average of 8300 FD-mapped reads per library, we assembled 347 sequences, corresponding to 215 annotated genes, and identified 10 previously unannotated genes, 15 polycistronic transcripts and three genes supposedly localized in the gaps of the FD92 draft genome. Furthermore, we improved the annotation of 44 genes with the addition of 5'/3' untranslated regions. Functional classification revealed that the most expressed genes were either related to translation and protein biosynthesis or hypothetical proteins with unknown function. Some of these hypothetical proteins were predicted to be secreted, so they could be bacterial effectors with a potential role in modulating the interaction with the host plant. Interestingly, qRT-PCR validation of the RNA-Seq expression values confirmed that a group II intron represented the FD genomic region with the highest expression during grapevine infection. This mobile element may contribute to the genomic plasticity that is necessary for the phytoplasma to increase its fitness and endorse host-adaptive strategies. CONCLUSIONS: The RNA-Seq technology was successfully applied for the first time to analyse the FD global transcriptome profile during grapevine infection. Our results provided new insights into the transcriptional organization and gene structure of FD. This may represent the starting point for the application of high-throughput sequencing technologies to study differential expression in FD and in other phytoplasmas with an unprecedented resolution.

Genome wide sequence analysis grants unbiased definition of species boundaries in "Candidatus Phytoplasma".

The phytoplasmas are currently named using the Candidatus category, as the inability to grow them in vitro prevented (i) the performance of tests, such as DNA-DNA hybridization, that are regarded as necessary to establish species boundaries, and (ii) the deposition of type strains in culture collections. The recent accession to complete or nearly complete genome sequence information disclosed the opportunity to apply to the uncultivable phytoplasmas the same taxonomic approaches used for other bacteria. In this work, the genomes of 14 strains, belonging to the 16SrI, 16SrIII, 16SrV and 16SrX groups, including the species "Ca. P. asteris", "Ca. P. mali", "Ca. P. pyri", "Ca. P. pruni", and "Ca. P. australiense" were analyzed along with Acholeplasma laidlawi, to determine their taxonomic relatedness. Average nucleotide index (ANIm), tetranucleotide signature frequency correlation index (Tetra), and multilocus sequence analysis of 107 shared genes using both phylogenetic inference of concatenated (DNA and amino acid) sequences and consensus networks, were carried out. The results were in large agreement with the previously established 16S rDNA based classification schemes. Moreover, the taxonomic relationships within the 16SrI, 16SrIII and 16SrX groups, that represent clusters of strains whose relatedness could not be determined by 16SrDNA analysis, could be comparatively evaluated with non-subjective criteria. "Ca. P. mali" and "Ca. P. pyri" were found to meet the genome characteristics for the retention into two different, yet strictly related species; representatives of subgroups 16SrI-A and 16SrI-B were also found to meet the standards used in other bacteria to distinguish separate species; the genomes of the strains belonging to 16SrIII were found more closely related, suggesting that their subdivision into Candidatus species should be approached with caution.

Mapping the phytoplasma chromosome.

Physical and genetic mapping of the phytoplasma chromosome can be a useful tool in a genome sequencing project in order to assemble the in silico-predicted contigs robustly. Mapping consists of four distinct steps: preparation of phytoplasma chromosomes from infected plants, single- and double-digestion of chromosomes with rare-cutting restriction enzymes, separation of large DNA fragments by pulsed-field gel electrophoresis, and hybridization with various genetic markers. Materials and methods needed for each step are described and the technique is illustrated using the flavescence dorée phytoplasma genome map as an example.

Genetic diversity of European phytoplasmas of the 16SrV taxonomic group and proposal of 'Candidatus Phytoplasma rubi'.

In addition to the grapevine flavescence dorée phytoplasmas, other members of taxonomic group 16SrV phytoplasmas infect grapevines, alders and species of the genera Clematis and Rubus in Europe. In order to investigate which phytoplasmas constitute discrete, species-level taxa, several strains were analysed by comparing their 16S rRNA gene sequences and a set of five housekeeping genes. Whereas 16S rRNA gene sequence similarity values were >97.5 %, the proposed threshold to distinguish two 'Candidatus Phytoplasma' taxa, phylogenetic analysis of the combined sequences of the tuf, rplV-rpsC, rplF-rplR, map and uvrB-degV genetic loci showed that two discrete phylogenetic clusters could be clearly distinguished. The first cluster grouped flavescence dorée (FD) phytoplasmas, alder yellows (AldY) phytoplasmas, Clematis (CL) phytoplasmas and the Palatinate grapevine yellows (PGY) phytoplasmas. The second cluster comprised Rubus stunt (RS) phytoplasmas. In addition to the specificity of the insect vector, the Rubus stunt phytoplasma contained specific sequences in the 16S rRNA gene. Hence, the Rubus stunt phytoplasma 16S rRNA gene was sufficiently differentiated to represent a novel putative taxon: 'Candidatus Phytoplasma rubi'.

Partial chromosome sequence of Spiroplasma citri reveals extensive viral invasion and important gene decay.

The assembly of 20,000 sequencing reads obtained from shotgun and chromosome-specific libraries of the Spiroplasma citri genome yielded 77 chromosomal contigs totaling 1,674 kbp (92%) of the 1,820-kbp chromosome. The largest chromosomal contigs were positioned on the physical and genetic maps constructed from pulsed-field gel electrophoresis and Southern blot hybridizations. Thirty-eight contigs were annotated, resulting in 1,908 predicted coding sequences (CDS) representing an overall coding density of only 74%. Cellular processes, cell metabolism, and structural-element CDS account for 29% of the coding capacity, CDS of external origin such as viruses and mobile elements account for 24% of the coding capacity, and CDS of unknown function account for 47% of the coding capacity. Among these, 21% of the CDS group into 63 paralog families. The organization of these paralogs into conserved blocks suggests that they represent potential mobile units. Phage-related sequences were particularly abundant and include plectrovirus SpV1 and SVGII3 and lambda-like SpV2 sequences. Sixty-nine copies of transposases belonging to four insertion sequence (IS) families (IS30, IS481, IS3, and ISNCY) were detected. Similarity analyses showed that 21% of chromosomal CDS were truncated compared to their bacterial orthologs. Transmembrane domains, including signal peptides, were predicted for 599 CDS, of which 58 were putative lipoproteins. S. citri has a Sec-dependent protein export pathway. Eighty-four CDS were assigned to transport, such as phosphoenolpyruvate phosphotransferase systems (PTS), the ATP binding cassette (ABC), and other transporters. Besides glycolytic and ATP synthesis pathways, it is noteworthy that S. citri possesses a nearly complete pathway for the biosynthesis of a terpenoid.

A chromosome map of the Flavescence doree phytoplasma.

The Flavescence dorée phytoplasma (FD-P), a non-cultivable, plant-pathogenic bacterium of the class Mollicutes, is the causal agent of a quarantine disease affecting vineyards of southern Europe, mainly in southern France and northern Italy. To investigate FD-P diversity and phytoplasma genetic determinants governing the FD-P life cycle, a genome project has been initiated. A physical map of the chromosome of FD-P strain FD92, purified from infected broad beans, was constructed by performing restriction digests of the chromosome and resolving the fragments by PFGE. Single and double digestions of the chromosome with the enzymes SalI, BssHII, MluI and EagI were performed and used to map 13 restriction sites on the FD-P chromosome. The size of the chromosome was calculated to be 671 kbp. Southern blot analyses using cloned phytoplasma probes were carried out to assist in the arrangement of contiguous restriction fragments and to map eight genetic loci, including the two rRNA operons, the tuf, uvrB-degV and secY-map (FD9) genes, the FD2 marker and two orphan sequences (FDDH29 and FDSH05) isolated through subtractive suppression hybridization.

The abundant extrachromosomal DNA content of the Spiroplasma citri GII3-3X genome.

BACKGROUND: Spiroplama citri, the causal agent of citrus stubborn disease, is a bacterium of the class Mollicutes and is transmitted by phloem-feeding leafhopper vectors. In order to characterize candidate genes potentially involved in spiroplasma transmission and pathogenicity, the genome of S. citri strain GII3-3X is currently being deciphered. RESULTS: Assembling 20,000 sequencing reads generated seven circular contigs, none of which fit the 1.8 Mb chromosome map or carried chromosomal markers. These contigs correspond to seven plasmids: pSci1 to pSci6, with sizes ranging from 12.9 to 35.3 kbp and pSciA of 7.8 kbp. Plasmids pSci were detected as multiple copies in strain GII3-3X. Plasmid copy numbers of pSci1-6, as deduced from sequencing coverage, were estimated at 10 to 14 copies per spiroplasma cell, representing 1.6 Mb of extrachromosomal DNA. Genes encoding proteins of the TrsE-TraE, Mob, TraD-TraG, and Soj-ParA protein families were predicted in most of the pSci sequences, in addition to members of 14 protein families of unknown function. Plasmid pSci6 encodes protein P32, a marker of insect transmissibility. Plasmids pSci1-5 code for eight different S. citri adhesion-related proteins (ScARPs) that are homologous to the previously described protein P89 and the S. kunkelii SkARP1. Conserved signal peptides and C-terminal transmembrane alpha helices were predicted in all ScARPs. The predicted surface-exposed N-terminal region possesses the following elements: (i) 6 to 8 repeats of 39 to 42 amino acids each (sarpin repeats), (ii) a central conserved region of 330 amino acids followed by (iii) a more variable domain of about 110 amino acids. The C-terminus, predicted to be cytoplasmic, consists of a 27 amino acid stretch enriched in arginine and lysine (KR) and an optional 23 amino acid stretch enriched in lysine, aspartate and glutamate (KDE). Plasmids pSci mainly present a linear increase of cumulative GC skew except in regions presenting conserved hairpin structures. CONCLUSION: The genome of S. citri GII3-3X is characterized by abundant extrachromosomal elements. The pSci plasmids could not only be vertically inherited but also horizontally transmitted, as they encode proteins usually involved in DNA element partitioning and cell to cell DNA transfer. Because plasmids pSci1-5 encode surface proteins of the ScARP family and pSci6 was recently shown to confer insect transmissibility, diversity and abundance of S. citri plasmids may essentially aid the rapid adaptation of S. citri to more efficient transmission by different insect vectors and to various plant hosts.