Characterization of intronic structures and alternative splicing in Phytophthora sojae by comparative analysis of expressed sequence tags and genomic sequences

Phytophthora effector PSR1 hijacks the host pre-mRNA splicing machinery to modulate small RNA biogenesis and plant immunity

Phytophthora effector PSR1 suppresses small RNA (sRNA)-mediated immunity in plants, but the underlying mechanism remains unknown. Here, we show that Phytophthora suppressor of RNA silencing 1 (PSR1) contributes to the pathogenicity of Phytophthora sojae and specifically binds to three conserved C-terminal domains of the eukaryotic PSR1-Interacting Protein 1 (PINP1). PINP1 encodes PRP16, a core pre-mRNA splicing factor that unwinds RNA duplexes and binds to primary microRNA transcripts and general RNAs. Intriguingly, PSR1 decreased both RNA helicase and RNA-binding activity of PINP1, thereby dampening sRNA biogenesis and RNA metabolism. The PSR1-PINP1 interaction caused global changes in alternative splicing (AS). A total of 5,135 genes simultaneously exhibited mis-splicing in both PSR1-overexpressing and PINP1-silenced plants. AS upregulated many mRNA transcripts that had their introns retained. The high occurrence of intron retention in AS-induced transcripts significantly promoted Phytophthora pathogen infection in Nicotiana benthamiana, and this might be caused by the production of truncated proteins. Taken together, our findings reveal a key role for PINP1 in regulating sRNA biogenesis and plant immunity.

Animal, Fungi, and Plant Genome Sequences Harbor Different Non-Canonical Splice Sites

Most protein-encoding genes in eukaryotes contain introns, which are interwoven with exons. Introns need to be removed from initial transcripts in order to generate the final messenger RNA (mRNA), which can be translated into an amino acid sequence. Precise excision of introns by the spliceosome requires conserved dinucleotides, which mark the splice sites. However, there are variations of the highly conserved combination of GT at the 5' end and AG at the 3' end of an intron in the genome. GC-AG and AT-AC are two major non-canonical splice site combinations, which have been known for years. Recently, various minor non-canonical splice site combinations were detected with numerous dinucleotide permutations. Here, we expand systematic investigations of non-canonical splice site combinations in plants across eukaryotes by analyzing fungal and animal genome sequences. Comparisons of splice site combinations between these three kingdoms revealed several differences, such as an apparently increased CT-AC frequency in fungal genome sequences. Canonical GT-AG splice site combinations in antisense transcripts are a likely explanation for this observation, thus indicating annotation errors. In addition, high numbers of GA-AG splice site combinations were observed in Eurytemoraaffinis and Oikopleuradioica. A variant in one U1 small nuclear RNA (snRNA) isoform might allow the recognition of GA as a 5' splice site. In depth investigation of splice site usage based on RNA-Seq read mappings indicates a generally higher flexibility of the 3' splice site compared to the 5' splice site across animals, fungi, and plants.

Splice Sites Seldom Slide: Intron Evolution in Oomycetes

We examine exon junctions near apparent amino acid insertions and deletions in alignments of orthologous protein-coding genes. In 1,917 ortholog families across nine oomycete genomes, 10–20% of introns are near an alignment gap, indicating at first sight that splice-site displacements are frequent. We designed a robust algorithmic procedure for the delineation of intron-containing homologous regions, and combined it with a parsimony-based reconstruction of intron loss, gain, and splice-site shift events on a phylogeny. The reconstruction implies that 12% of introns underwent an acceptor-site shift, and 10% underwent a donor-site shift. In order to offset gene annotation problems, we amended the procedure with the reannotation of intron boundaries using alignment evidence. The corresponding reconstruction involves much fewer intron gain and splice-site shift events. The frequency of acceptor- and donor-side shifts drops to 4% and 3%, respectively, which are not much different from what one would expect by random codon insertions and deletions. In other words, gaps near exon junctions are mostly artifacts of gene annotation rather than evidence of sliding intron boundaries. Our study underscores the importance of using well-supported gene structure annotations in comparative studies. When transcription evidence is not available, we propose a robust ancestral reconstruction procedure that corrects misannotated intron boundaries using sequence alignments. The results corroborate the view that boundary shifts and complete intron sliding are only accidental in eukaryotic genome evolution and have a negligible impact on protein diversity.

The Pectin Methylesterase Gene Complement of Phytophthora sojae: Structural and Functional Analyses, and the Evolutionary Relationships with Its Oomycete Homologs

Phytophthora sojae is an oomycete pathogen that causes the disease known as root and stem rot in soybean plants, frequently leading to massive economic damage. Additionally, P. sojae is increasingly being utilized as a model for phytopathogenic oomycete research. Despite the economic and scientific importance of P. sojae, the mechanism by which it penetrates the host roots is not yet fully understood. It has been found that oomycetes are not capable of penetrating the cell wall solely through mechanical force, suggesting that alternative factors facilitate breakdown of the host cell wall. Pectin methylesterases have been suggested to be important for Phytophthora pathogenicity, but no data exist on their role in the P. sojae infection process. We have scanned the newly revised version of the annotated P. sojae genome for the presence of putative pectin methylesterases genes and conducted a sequence analysis of all gene models found. We also searched for potential regulatory motifs in the promoter region of the proposed P. sojae models, and investigated the gene expression levels throughout the early course of infection on soybean plants. We found that P. sojae contains a large repertoire of pectin methylesterase-coding genes and that most of these genes display similar motifs in the promoter region, indicating the possibility of a shared regulatory mechanism. Phylogenetic analyses confirmed the evolutionary relatedness of the pectin methylesterase-coding genes within and across Phytophthora spp. In addition, the gene duplication events that led to the emergence of this gene family appear to have occurred prior to many speciation events in the genus Phytophthora. Our results also indicate that the highest levels of expression occurred in the first 24 hours post inoculation, with expression falling after this time. Our study provides evidence that pectin methylesterases may be important for the early action of the P. sojae infection process.

Alternative Splicing in the Obligate Biotrophic Oomycete Pathogen Pseudoperonospora cubensis

Pseudoperonospora cubensis is an obligate pathogen and causative agent of cucurbit downy mildew. To help advance our understanding of the pathogenicity of P. cubensis, we used RNA-Seq to improve the quality of its reference genome sequence. We also characterized the RNA-Seq dataset to inventory transcript isoforms and infer alternative splicing during different stages of its development. Almost half of the original gene annotations were improved and nearly 4,000 previously unannotated genes were identified. We also demonstrated that approximately 24% of the expressed genome and nearly 55% of the intron-containing genes from P. cubensis had evidence for alternative splicing. Our analyses revealed that intron retention is the predominant alternative splicing type in P. cubensis, with alternative 5'- and alternative 3'-splice sites occurring at lower frequencies. Representatives of the newly identified genes and predicted alternatively spliced transcripts were experimentally validated. The results presented herein highlight the utility of RNA-Seq for improving draft genome annotations and, through this approach, we demonstrate that alternative splicing occurs more frequently than previously predicted. In total, the current study provides evidence that alternative splicing plays a key role in transcriptome regulation and proteome diversification in plant-pathogenic oomycetes.

Origin of spliceosomal introns and alternative splicing

In this work we review the current knowledge on the prehistory, origins, and evolution of spliceosomal introns. First, we briefly outline the major features of the different types of introns, with particular emphasis on the nonspliceosomal self-splicing group II introns, which are widely thought to be the ancestors of spliceosomal introns. Next, we discuss the main scenarios proposed for the origin and proliferation of spliceosomal introns, an event intimately linked to eukaryogenesis. We then summarize the evidence that suggests that the last eukaryotic common ancestor (LECA) had remarkably high intron densities and many associated characteristics resembling modern intron-rich genomes. From this intron-rich LECA, the different eukaryotic lineages have taken very distinct evolutionary paths leading to profoundly diverged modern genome structures. Finally, we discuss the origins of alternative splicing and the qualitative differences in alternative splicing forms and functions across lineages.

Comparative structural and functional analysis of genes encoding pectin methylesterases in Phytophthora spp

We have scanned the Phytophthora infestans, P. ramorum, and P. sojae genomes for the presence of putative pectin methylesterase genes and conducted a sequence analysis of all gene models found. We also searched for potential regulatory motifs in the promoter region of the proposed P. infestans models, and investigated the gene expression levels throughout the course of P. infestans infection on potato plants, using in planta and detached leaf assays. We found that genes located on contiguous chromosomal regions contain similar motifs in the promoter region, indicating the possibility of a shared regulatory mechanism. Results of our investigations also suggest that, during the pathogenicity process, the expression levels of some of the analyzed genes vary considerably when compared to basal expression observed in in vitro cultures of non-sporulating mycelium. These results were observed both in planta and in detached leaf assays.

Phylogenetic and transcriptional analysis of an expanded bZIP transcription factor family in Phytophthora sojae

BACKGROUND

Basic leucine zipper (bZIP) transcription factors are present exclusively in eukaryotes and constitute one of the largest and most diverse transcription factor families. The proteins are responsible for central developmental and physiological processes in plants, animals, and fungi, including the pathogenicity of fungal plant pathogens. However, there is limited understanding of bZIPs in oomycetes, which are fungus-like organisms in the kingdom Stramenopila. Oomycetes include many destructive plant pathogens, including the well-studied species Phytophthora sojae, which causes soybean stem and root rot.

RESULTS

Candidate bZIPs encoded in the genomes of P. sojae and four other oomycetes, two diatoms, and two fungal species were predicted using bioinformatic methods. Comparative analysis revealed expanded numbers of bZIP candidates in oomycetes, especially the Phytophthora species, due to the expansion of several novel bZIP classes whose highly conserved asparagines in basic DNA-binding regions were substituted by other residues such as cysteine. The majority of these novel bZIP classes were mostly restricted to oomycetes. The large number of novel bZIPs appears to be the result of widespread gene duplications during oomycete evolution. The majority of P. sojae bZIP candidates, including both conventional and novel bZIP classes, were predicted to contain canonical protein secondary structures. Detection of gene transcripts using digital gene expression profiling and qRT-PCR suggested that most of the candidates were not pseudogenes. The major transcriptional shifts of bZIPs occurred during the zoosporangia/zoospore/cyst and host infection stages. Several infection-associated bZIP genes were identified that were positively regulated by H2O2 exposure.

CONCLUSIONS

The identification of large classes of bZIP proteins in oomycetes with novel bZIP motif variants, that are conserved and developmentally regulated and thus presumably functional, extends our knowledge of this important family of eukaryotic transcription factors. It also lays the foundation for detailed studies of the roles of these proteins in development and infection in P. sojae and other oomycetes.

Effector identification in the lettuce downy mildew Bremia lactucae by massively parallel transcriptome sequencing

Lettuce downy mildew (Bremia lactucae) is a rapidly adapting oomycete pathogen affecting commercial lettuce cultivation. Oomycetes are known to use a diverse arsenal of secreted proteins (effectors) to manipulate their hosts. Two classes of effector are known to be translocated by the host: the RXLRs and Crinklers. To gain insight into the repertoire of effectors used by B. lactucae to manipulate its host, we performed massively parallel sequencing of cDNA derived from B. lactucae spores and infected lettuce (Lactuca sativa) seedlings. From over 2.3 million 454 GS FLX reads, 59 618 contigs were assembled representing both plant and pathogen transcripts. Of these, 19 663 contigs were determined to be of B. lactucae origin as they matched pathogen genome sequences (SOLiD) that were obtained from >270 million reads of spore-derived genomic DNA. After correction of cDNA sequencing errors with SOLiD data, translation into protein models and filtering, 16 372 protein models remained, 1023 of which were predicted to be secreted. This secretome included elicitins, necrosis and ethylene-inducing peptide 1-like proteins, glucanase inhibitors and lectins, and was enriched in cysteine-rich proteins. Candidate host-translocated effectors included 78 protein models with RXLR effector features. In addition, we found indications for an unknown number of Crinkler-like sequences. Similarity clustering of secreted proteins revealed additional effector candidates. We provide a first look at the transcriptome of B. lactucae and its encoded effector arsenal.

Dynamics and innovations within oomycete genomes: insights into biology, pathology, and evolution

The eukaryotic microbes known as oomycetes are common inhabitants of terrestrial and aquatic environments and include saprophytes and pathogens. Lifestyles of the pathogens extend from biotrophy to necrotrophy, obligate to facultative pathogenesis, and narrow to broad host ranges on plants or animals. Sequencing of several pathogens has revealed striking variation in genome size and content, a plastic set of genes related to pathogenesis, and adaptations associated with obligate biotrophy. Features of genome evolution include repeat-driven expansions, deletions, gene fusions, and horizontal gene transfer in a landscape organized into gene-dense and gene-sparse sectors and influenced by transposable elements. Gene expression profiles are also highly dynamic throughout oomycete life cycles, with transcriptional polymorphisms as well as differences in protein sequence contributing to variation. The genome projects have set the foundation for functional studies and should spur the sequencing of additional species, including more diverse pathogens and nonpathogens.

Alternative splicing of a multi-drug transporter from Pseudoperonospora cubensis generates an RXLR effector protein that elicits a rapid cell death

Pseudoperonospora cubensis, an obligate oomycete pathogen, is the causal agent of cucurbit downy mildew, a foliar disease of global economic importance. Similar to other oomycete plant pathogens, Ps. cubensis has a suite of RXLR and RXLR-like effector proteins, which likely function as virulence or avirulence determinants during the course of host infection. Using in silico analyses, we identified 271 candidate effector proteins within the Ps. cubensis genome with variable RXLR motifs. In extending this analysis, we present the functional characterization of one Ps. cubensis effector protein, RXLR protein 1 (PscRXLR1), and its closest Phytophthora infestans ortholog, PITG_17484, a member of the Drug/Metabolite Transporter (DMT) superfamily. To assess if such effector-non-effector pairs are common among oomycete plant pathogens, we examined the relationship(s) among putative ortholog pairs in Ps. cubensis and P. infestans. Of 271 predicted Ps. cubensis effector proteins, only 109 (41%) had a putative ortholog in P. infestans and evolutionary rate analysis of these orthologs shows that they are evolving significantly faster than most other genes. We found that PscRXLR1 was up-regulated during the early stages of infection of plants, and, moreover, that heterologous expression of PscRXLR1 in Nicotiana benthamiana elicits a rapid necrosis. More interestingly, we also demonstrate that PscRXLR1 arises as a product of alternative splicing, making this the first example of an alternative splicing event in plant pathogenic oomycetes transforming a non-effector gene to a functional effector protein. Taken together, these data suggest a role for PscRXLR1 in pathogenicity, and, in total, our data provide a basis for comparative analysis of candidate effector proteins and their non-effector orthologs as a means of understanding function and evolutionary history of pathogen effectors.

Transcriptome analysis of the model protozoan, Tetrahymena thermophila, using Deep RNA sequencing

BACKGROUND

The ciliated protozoan Tetrahymena thermophila is a well-studied single-celled eukaryote model organism for cellular and molecular biology. However, the lack of extensive T. thermophila cDNA libraries or a large expressed sequence tag (EST) database limited the quality of the original genome annotation.

METHODOLOGY/PRINCIPAL FINDINGS

This RNA-seq study describes the first deep sequencing analysis of the T. thermophila transcriptome during the three major stages of the life cycle: growth, starvation and conjugation. Uniquely mapped reads covered more than 96% of the 24,725 predicted gene models in the somatic genome. More than 1,000 new transcribed regions were identified. The great dynamic range of RNA-seq allowed detection of a nearly six order-of-magnitude range of measurable gene expression orchestrated by this cell. RNA-seq also allowed the first prediction of transcript untranslated regions (UTRs) and an updated (larger) size estimate of the T. thermophila transcriptome: 57 Mb, or about 55% of the somatic genome. Our study identified nearly 1,500 alternative splicing (AS) events distributed over 5.2% of T. thermophila genes. This percentage represents a two order-of-magnitude increase over previous EST-based estimates in Tetrahymena. Evidence of stage-specific regulation of alternative splicing was also obtained. Finally, our study allowed us to completely confirm about 26.8% of the genes originally predicted by the gene finder, to correct coding sequence boundaries and intron-exon junctions for about a third, and to reassign microarray probes and correct earlier microarray data.

CONCLUSIONS/SIGNIFICANCE

RNA-seq data significantly improve the genome annotation and provide a fully comprehensive view of the global transcriptome of T. thermophila. To our knowledge, 5.2% of T. thermophila genes with AS is the highest percentage of genes showing AS reported in a unicellular eukaryote. Tetrahymena thus becomes an excellent unicellular model eukaryote in which to investigate mechanisms of alternative splicing.