Deep sequencing of small RNAs in plants: applied bioinformatics

Small RNAs, including microRNA and short-interfering RNAs, play important roles in plants. In recent years, developments in sequencing technology have enabled the large-scale discovery of sRNAs in various cells, tissues and developmental stages and in response to various stresses. This review describes the bioinformatics challenges to analysing these large datasets of short-RNA sequences and some of the solutions to those challenges.     

The rulB gene of plasmid pWW0 is a hotspot for the site‐specific insertion of integron‐like elements found in the chromosomes of environmental Pseudomonas fluorescens group bacteria

The rulAB operon of Pseudomonas spp. confers fitness traits on the host and has been suggested to be a hotspot for insertion of mobile elements that carry avirulence genes. Here, for the first time, we show that rulB on plasmid pWW0 is a hotspot for the active site‐specific integration of related integron‐like elements (ILEs) found in six environmental pseudomonads (strains FH1–FH6). Integration into rulB on pWW0 occurred at position 6488 generating a 3 bp direct repeat. ILEs from FH1 and FH5 were 9403 bp in length and contained eight open reading frames (ORFs), while the ILE from FH4 was 16 233 bp in length and contained 16 ORFs. In all three ILEs, the first 5.1 kb (containing ORFs 1–4) were structurally conserved and contained three predicted site‐specific recombinases/integrases and a tetR homologue. Downstream of these resided ORFs of the ‘variable side’ with structural and sequence similarity to those encoding survival traits on the fitness enhancing plasmid pGRT1 (ILE_FH1 and ILE_FH5) and the NR‐II virulence region of genomic island PAGI‐5 (ILE_FH4). Collectively, these ILEs share features with the previously described type III protein secretion system effector ILEs and are considered important to host survival and transfer of fitness enhancing and (a)virulence genes between bacteria.     

Multidimensional protein identification technology (MudPIT) analysis of ubiquitinated proteins in plants

Protein conjugation with ubiquitin, known as ubiquitination, is a key regulatory mechanism to control protein abundance, localization, and activity in eukaryotic cells. To identify ubiquitin-dependent regulatory steps in plants, we developed a robust affinity purification/identification system for ubiquitinated proteins. Using GST-tagged ubiquitin binding domains, we performed a large scale affinity purification of ubiquitinated proteins from Arabidopsis cell suspension culture. High molecular weight ubiquitinated proteins were separated by SDS-PAGE, and the trypsin-digested samples were then analyzed by a multidimensional protein identification technology (MudPIT) system. A total of 294 proteins specifically bound by the GST-tagged ubiquitin binding domains were identified. From these we determined 85 ubiquitinated lysine residues in 56 proteins, confirming the enrichment of the target class of proteins. Our data provide the first view of the ubiquitinated proteome in plants. We also provide evidence that this technique can be broadly applied to the study of protein ubiquitination in diverse plant species.     

Integrated bioinformatic and phenotypic analysis of RpoN-dependent traits in the plant growth-promoting bacterium Pseudomonas fluorescens SBW25

The alternative sigma factor RpoN is a key regulator in the acclimation of Pseudomonas to complex natural environments. In this study we show that RpoN is required for efficient colonization of sugar beet seedlings by the plant growth-promoting bacterium Pseudomonas fluorescens SBW25, and use phenotypic and bioinformatic approaches to profile the RpoN-dependent traits and genes of P. fluorescens SBW25. RpoN is required for flagellar biosynthesis and for assimilation of a wide variety of nutrient sources including inorganic nitrogen, amino acids, sugar alcohols and dicarboxylic acids. Chemosensitivity assays indicate that RpoN-regulated genes contribute to acid tolerance and resistance to some antibiotics, including tetracyclines and aminoglycosides. Gain of function changes associated with loss of RpoN included increased tolerance to hydroxyurea and Guanazole. Bioinformatic predictions of RpoN-regulated genes show a close correspondence with phenotypic analyses of RpoN-regulated traits and suggest novel functions for RpoN in P. fluorescens, including regulation of poly(A) polymerase. The reduced plant colonization ability observed for an rpoN mutant of P. fluorescens is therefore likely to be due to defects in multiple traits including nutrient assimilation, protein secretion and stress tolerance.     

Draft genome sequence of Pseudomonas syringae pathovar syringae strain FF5, causal agent of stem tip dieback disease on ornamental pear

Pseudomonas syringae FF5 causes stem tip dieback disease on ornamental pear (Pyrus calleryana). Its genome encodes a complete type III secretion system (T3SS) and HopAC1, HopM1, AvrE1, HopI1, HopAA1, HopJ1, HopAH2, HopAH1, HopAG1, and HopAZ1. Lacking detectable homologues of other T3SS effectors, it may encode novel, undiscovered effectors.     

Comparative Genome Analysis Provides Insights into the Evolution and Adaptation of Pseudomonas syringae pv. aesculi on Aesculus hippocastanum

A recently emerging bleeding canker disease, caused by Pseudomonas syringae pathovar aesculi (Pae), is threatening European horse chestnut in northwest Europe. Very little is known about the origin and biology of this new disease. We used the nucleotide sequences of seven commonly used marker genes to investigate the phylogeny of three strains isolated recently from bleeding stem cankers on European horse chestnut in Britain (E-Pae). On the basis of these sequences alone, the E-Pae strains were identical to the Pae type-strain (I-Pae), isolated from leaf spots on Indian horse chestnut in India in 1969. The phylogenetic analyses also showed that Pae belongs to a distinct clade of P. syringae pathovars adapted to woody hosts. We generated genome-wide Illumina sequence data from the three E-Pae strains and one strain of I-Pae. Comparative genomic analyses revealed pathovar-specific genomic regions in Pae potentially implicated in virulence on a tree host, including genes for the catabolism of plant-derived aromatic compounds and enterobactin synthesis. Several gene clusters displayed intra-pathovar variation, including those encoding type IV secretion, a novel fatty acid biosynthesis pathway and a sucrose uptake pathway. Rates of single nucleotide polymorphisms in the four Pae genomes indicate that the three E-Pae strains diverged from each other much more recently than they diverged from I-Pae. The very low genetic diversity among the three geographically distinct E-Pae strains suggests that they originate from a single, recent introduction into Britain, thus highlighting the serious environmental risks posed by the spread of an exotic plant pathogenic bacterium to a new geographic location. The genomic regions in Pae that are absent from other P. syringae pathovars that infect herbaceous hosts may represent candidate genetic adaptations to infection of the woody parts of the tree.     

The plant pathogen Pseudomonas syringae pv. tomato is genetically monomorphic and under strong selection to evade tomato immunity

Recently, genome sequencing of many isolates of genetically monomorphic bacterial human pathogens has given new insights into pathogen microevolution and phylogeography. Here, we report a genome-based micro-evolutionary study of a bacterial plant pathogen, Pseudomonas syringae pv. tomato. Only 267 mutations were identified between five sequenced isolates in 3,543,009 nt of analyzed genome sequence, which suggests a recent evolutionary origin of this pathogen. Further analysis with genome-derived markers of 89 world-wide isolates showed that several genotypes exist in North America and in Europe indicating frequent pathogen movement between these world regions. Genome-derived markers and molecular analyses of key pathogen loci important for virulence and motility both suggest ongoing adaptation to the tomato host. A mutational hotspot was found in the type III-secreted effector gene hopM1. These mutations abolish the cell death triggering activity of the full-length protein indicating strong selection for loss of function of this effector, which was previously considered a virulence factor. Two non-synonymous mutations in the flagellin-encoding gene fliC allowed identifying a new microbe associated molecular pattern (MAMP) in a region distinct from the known MAMP flg22. Interestingly, the ancestral allele of this MAMP induces a stronger tomato immune response than the derived alleles. The ancestral allele has largely disappeared from today's Pto populations suggesting that flagellin-triggered immunity limits pathogen fitness even in highly virulent pathogens. An additional non-synonymous mutation was identified in flg22 in South American isolates. Therefore, MAMPs are more variable than expected differing even between otherwise almost identical isolates of the same pathogen strain.     

Coverage-based consensus calling (CbCC) of short sequence reads and comparison of CbCC results to identify SNPs in chickpea (Cicer arietinum; Fabaceae), a crop species without a reference genome

? Premise of the study: Next-generation sequencing (NGS) technologies are frequently used for resequencing and mining of single nucleotide polymorphisms (SNPs) by comparison to a reference genome. In crop species such as chickpea (Cicer arietinum) that lack a reference genome sequence, NGS-based SNP discovery is a challenge. Therefore, unlike probability-based statistical approaches for consensus calling and by comparison with a reference sequence, a coverage-based consensus calling (CbCC) approach was applied and two genotypes were compared for SNP identification. ? Methods: A CbCC approach is used in this study with four commonly used short read alignment tools (Maq, Bowtie, Novoalign, and SOAP2) and 15.7 and 22.1 million Illumina reads for chickpea genotypes ICC4958 and ICC1882, together with the chickpea trancriptome assembly (CaTA). ? Key results: A nonredundant set of 4543 SNPs was identified between two chickpea genotypes. Experimental validation of 224 randomly selected SNPs showed superiority of Maq among individual tools, as 50.0% of SNPs predicted by Maq were true SNPs. For combinations of two tools, greatest accuracy (55.7%) was reported for Maq and Bowtie, with a combination of Bowtie, Maq, and Novoalign identifying 61.5% true SNPs. SNP prediction accuracy generally increased with increasing reads depth. ? Conclusions: This study provides a benchmark comparison of tools as well as read depths for four commonly used tools for NGS SNP discovery in a crop species without a reference genome sequence. In addition, a large number of SNPs have been identified in chickpea that would be useful for molecular breeding.     

Genome Update. Let the consumer beware: Streptomyces genome sequence quality

A genome sequence assembly represents a model of a genome. This article explores some tools and methods for assessing the quality of an assembly, using publicly available data for Streptomyces species as the example. There is great variability in quality of assemblies deposited in GenBank. Only in a small minority of these assemblies are the raw data available, enabling full appraisal of the assembly quality.