Insights into organ-specific pathogen defense responses in plants: RNA-seq analysis of potato tuber-Phytophthora infestans interactions

The host-pathogen interaction between wheat and yellow rust induces temporally coordinated waves of gene expression

Background

Understanding how plants and pathogens modulate gene expression during the host-pathogen interaction is key to uncovering the molecular mechanisms that regulate disease progression. Recent advances in sequencing technologies have provided new opportunities to decode the complexity of such interactions. In this study, we used an RNA-based sequencing approach (RNA-seq) to assess the global expression profiles of the wheat yellow rust pathogen Puccinia striiformis f. sp. tritici (PST) and its host during infection.

Results

We performed a detailed RNA-seq time-course for a susceptible and a resistant wheat host infected with PST. This study (i) defined the global gene expression profiles for PST and its wheat host, (ii) substantially improved the gene models for PST, (iii) evaluated the utility of several programmes for quantification of global gene expression for PST and wheat, and (iv) identified clusters of differentially expressed genes in the host and pathogen. By focusing on components of the defence response in susceptible and resistant hosts, we were able to visualise the effect of PST infection on the expression of various defence components and host immune receptors.

Conclusions

Our data showed sequential, temporally coordinated activation and suppression of expression of a suite of immune-response regulators that varied between compatible and incompatible interactions. These findings provide the framework for a better understanding of how PST causes disease and support the idea that PST can suppress the expression of defence components in wheat to successfully colonize a susceptible host.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2684-4) contains supplementary material, which is available to authorized users.

Dual RNA-Sequencing of Eucalyptus nitens during Phytophthora cinnamomi Challenge Reveals Pathogen and Host Factors Influencing Compatibility

Damage caused by Phytophthora cinnamomi Rands remains an important concern on forest tree species. The pathogen causes root and collar rot, stem cankers, and dieback of various economically important Eucalyptus spp. In South Africa, susceptible cold tolerant Eucalyptus plantations have been affected by various Phytophthora spp. with P. cinnamomi considered one of the most virulent. The molecular basis of this compatible interaction is poorly understood. In this study, susceptible Eucalyptus nitens plants were stem inoculated with P. cinnamomi and tissue was harvested five days post inoculation. Dual RNA-sequencing, a technique which allows the concurrent detection of both pathogen and host transcripts during infection, was performed. Approximately 1% of the reads mapped to the draft genome of P. cinnamomi while 78% of the reads mapped to the Eucalyptus grandis genome. The highest expressed P. cinnamomi gene in planta was a putative crinkler effector (CRN1). Phylogenetic analysis indicated the high similarity of this P. cinnamomi CRN1 to that of Phytophthora infestans. Some CRN effectors are known to target host nuclei to suppress defense. In the host, over 1400 genes were significantly differentially expressed in comparison to mock inoculated trees, including suites of pathogenesis related (PR) genes. In particular, a PR-9 peroxidase gene with a high similarity to a Carica papaya PR-9 ortholog previously shown to be suppressed upon infection by Phytophthora palmivora was down-regulated two-fold. This PR-9 gene may represent a cross-species effector target during P. cinnamomi infection. This study identified pathogenicity factors, potential manipulation targets, and attempted host defense mechanisms activated by E. nitens that contributed to the susceptible outcome of the interaction.

Host Transcriptional Profiling at Early and Later Stages of the Compatible Interaction Between Phaseolus vulgaris and Meloidogyne incognita

The root-knot nematode (Meloidogyne incognita) is one of most devastating pathogens that attack the common bean crop. Although there is evidence that some cultivars have race-specific resistance against M. incognita, these resistance sources have not proved effective, and nematodes are able to circumvent the host's defense system. We constructed RNA-seq based libraries and used a high-throughput sequencing platform to analyze the plant responses to M. incognita. Assessments were performed at 4 and 10 days after inoculation corresponding to the stages of nematode penetration and giant cell development, respectively. Large-scale transcript mapping to the common bean reference genome (G19833) resulted in the identification of 27,195 unigenes. Of these, 797 host genes were found to be differentially expressed. The functional annotation results confirm the complex interplay between abiotic and biotic stress signaling pathways. High expression levels of the wounding-responsive genes were observed over the interaction. At early response, an overexpression of the N gene, a TIR-NBS-LRR resistance gene, was understood as a host attempt to overcome the pathogen attack. However, the repression of heat shock proteins resulted in a lack of reactive oxygen species accumulation and absence of a hypersensitive response. Furthermore, the host basal response was broken by the repression of the ethylene/jasmonate pathway later in the response, resulting in a continuous compatible process with consequent plant susceptibility.

Molecular candidates for early-stage flower-to-fruit transition in stenospermocarpic table grape (Vitis vinifera L.) inflorescences ascribed by differential transcriptome and metabolome profiles

Flower-to-fruit transition depends of nutrient availability and regulation at the molecular level by sugar and hormone signalling crosstalk. However, in most species, the identities of fruit initiation regulators and their targets are largely unknown. To ascertain the main pathways involved in stenospermocarpic table grape fruit set, comprehensive transcriptional and metabolomic analyses were conducted specifically targeting the early phase of this developmental stage in 'Thompson Seedless'. The high-throughput analyses performed disclosed the involvement of 496 differentially expressed genes and 28 differently accumulated metabolites in the sampled inflorescences. Our data show broad transcriptome reprogramming of molecule transporters, globally down-regulating gene expression, and suggest that regulation of sugar- and hormone-mediated pathways determines the downstream activation of berry development. The most affected gene was the SWEET14 sugar transporter. Hormone-related transcription changes were observed associated with increased indole-3-acetic acid, stimulation of ethylene and gibberellin metabolisms and cytokinin degradation, and regulation of MADS-box and AP2-like ethylene-responsive transcription factor expression. Secondary metabolism, the most representative biological process at transcriptome level, was predominantly repressed. The results add to the knowledge of molecular events occurring in grapevine inflorescence fruit set and provide a list of candidates, paving the way for genetic manipulation aimed at model research and plant breeding.

Transcriptomic Crosstalk between Fungal Invasive Pathogens and Their Host Cells: Opportunities and Challenges for Next-Generation Sequencing Methods

Fungal invasive infections are an increasing health problem. The intrinsic complexity of pathogenic fungi and the unmet clinical need for new and more effective treatments requires a detailed knowledge of the infection process. During infection, fungal pathogens are able to trigger a specific transcriptional program in their host cells. The detailed knowledge of this transcriptional program will allow for a better understanding of the infection process and consequently will help in the future design of more efficient therapeutic strategies. Simultaneous transcriptomic studies of pathogen and host by high-throughput sequencing (dual RNA-seq) is an unbiased protocol to understand the intricate regulatory networks underlying the infectious process. This protocol is starting to be applied to the study of the interactions between fungal pathogens and their hosts. To date, our knowledge of the molecular basis of infection for fungal pathogens is still very limited, and the putative role of regulatory players such as non-coding RNAs or epigenetic factors remains elusive. The wider application of high-throughput transcriptomics in the near future will help to understand the fungal mechanisms for colonization and survival, as well as to characterize the molecular responses of the host cell against a fungal infection.

Genomic limitations to RNA sequencing expression profiling

The field of genomics has grown rapidly with the advent of massively parallel sequencing technologies, allowing for novel biological insights with regards to genomic, transcriptomic, and epigenomic variation. One widely utilized application of high-throughput sequencing is transcriptional profiling using RNA sequencing (RNAseq). Understanding the limitations of a technology is critical for accurate biological interpretations, and clear interpretation of RNAseq data can be difficult in species with complex genomes. To understand the limitations of accurate profiling of expression levels we simulated RNAseq reads from annotated gene models in several plant species including Arabidopsis, brachypodium, maize, potato, rice, soybean, and tomato. The simulated reads were aligned using various parameters such as unique versus multiple read alignments. This allowed the identification of genes recalcitrant to RNAseq analyses by having over- and/or under-estimated expression levels. In maize, over 25% of genes deviated by more than 20% from the expected count values, suggesting the need for cautious interpretation of RNAseq data for certain genes. The reasons identified for deviation from expected expression varied between species due to differences in genome structure including, but not limited to, genes encoding short transcripts, overlapping gene models, and gene family size. Utilizing existing empirical datasets we demonstrate the potential for biological misinterpretation resulting from inclusion of 'flagged genes' in analyses. While RNAseq is a powerful tool for understanding biology, there are limitations to this technology that need to be understood in order to improve our biological interpretations.

Functional genomics of tomato: opportunities and challenges in post-genome NGS era

The Tomato Genome Sequencing Project represented a landmark venture in the history of sequencing projects where both Sanger's and next-generation sequencing (NGS) technologies were employed, and a highly accurate and one of the best assembled plant genomes along with a draft of the wild relative, Solanum pimpinellifolium, were released in 2012. However, the functional potential of the major portion of this newly generated resource is still undefined. The very first challenge before scientists working on tomato functional biology is to exploit this high-quality reference sequence for tapping of the wealth of genetic variants for improving agronomic traits in cultivated tomatoes. The sequence data generated recently by 150 Tomato Genome Consortium would further uncover the natural alleles present in different tomato genotypes. Therefore, we found it relevant to have a fresh outlook on tomato functional genomics in the context of application of NGS technologies in its post-genome sequencing phase. Herein, we provide an overview how NGS technologies vis-a-vis available reference sequence have assisted each other for their mutual improvement and how their combined use could further facilitate the development of other 'omics' tools, required to propel the Solanaceae research. Additionally, we highlight the challenges associated with the application of these cutting-edge technologies.

De novo assembly and characterisation of the field pea transcriptome using RNA-Seq

BACKGROUND

Field pea (Pisum sativum L.) is a cool-season grain legume that is cultivated world-wide for both human consumption and stock-feed purposes. Enhancement of genetic and genomic resources for field pea will permit improved understanding of the control of traits relevant to crop productivity and quality. Advances in second-generation sequencing and associated bioinformatics analysis now provide unprecedented opportunities for the development of such resources. The objective of this study was to perform transcriptome sequencing and characterisation from two genotypes of field pea that differ in terms of seed and plant morphological characteristics.

RESULTS

Transcriptome sequencing was performed with RNA templates from multiple tissues of the field pea genotypes Kaspa and Parafield. Tissue samples were collected at various growth stages, and a total of 23 cDNA libraries were sequenced using Illumina high-throughput sequencing platforms. A total of 407 and 352 million paired-end reads from the Kaspa and Parafield transcriptomes, respectively were assembled into 129,282 and 149,272 contigs, which were filtered on the basis of known gene annotations, presence of open reading frames (ORFs), reciprocal matches and degree of coverage. Totals of 126,335 contigs from Kaspa and 145,730 from Parafield were subsequently selected as the reference set. Reciprocal sequence analysis revealed that c. 87% of contigs were expressed in both cultivars, while a small proportion were unique to each genotype. Reads from different libraries were aligned to the genotype-specific assemblies in order to identify and characterise expression of contigs on a tissue-specific basis, of which 87% were expressed in more than one tissue, while others showed distinct expression patterns in specific tissues, providing unique transcriptome signatures.

CONCLUSION

This study provided a comprehensive assembled and annotated transcriptome set for field pea that can be used for development of genetic markers, in order to assess genetic diversity, construct linkage maps, perform trait-dissection and implement whole-genome selection strategies in varietal improvement programs, as well to identify target genes for genetic modification approaches on the basis of annotation and expression analysis. In addition, the reference field pea transcriptome will prove highly valuable for comparative genomics studies and construction of a finalised genome sequence.

Potato Tuber Blight Resistance Phenotypes Correlate with RB Transgene Transcript Levels in an Age-Dependent Manner

Plants have evolved strategies and mechanisms to detect and respond to pathogen attack. Different organs of the same plant may be subjected to different environments (e.g., aboveground versus belowground) and pathogens with different lifestyles. Accordingly, plants commonly need to tailor defense strategies in an organ-specific manner. Phytophthora infestans, causal agent of potato late blight disease, infects both aboveground foliage and belowground tubers. We examined the efficacy of transgene RB (known for conferring foliar late blight resistance) in defending against tuber late blight disease. Our results indicate that the presence of the transgene has a positive yet only marginally significant effect on tuber disease resistance on average. However, a significant association between transgene transcript levels and tuber resistance was established for specific transformed lines in an age-dependent manner, with higher transcript levels indicating enhanced tuber resistance. Thus, RB has potential to function in both foliage and tuber to impart late blight resistance. Our data suggest that organ-specific resistance might result directly from transcriptional regulation of the resistance gene itself.

In-Depth Transcriptome Sequencing of Mexican Lime Trees Infected with Candidatus Phytoplasma aurantifolia

Witches' broom disease of acid lime greatly affects the production of Mexican lime in Iran. It is caused by a phytoplasma (Candidatus Phytoplasma aurantifolia). However, the molecular mechanisms that underlie phytoplasma pathogenicity and the mode of interactions with host plants are largely unknown. Here, high-throughput transcriptome sequencing was conducted to explore gene expression signatures associated with phytoplasma infection in Mexican lime trees. We assembled 78,185 unique transcript sequences (unigenes) with an average length of 530 nt. Of these, 41,805 (53.4%) were annotated against the NCBI non-redundant (nr) protein database using a BLASTx search (e-value ≤ 1e-5). When the abundances of unigenes in healthy and infected plants were compared, 2,805 transcripts showed significant differences (false discovery rate ≤ 0.001 and log2 ratio ≥ 1.5). These differentially expressed genes (DEGs) were significantly enriched in 43 KEGG metabolic and regulatory pathways. The up-regulated DEGs were mainly categorized into pathways with possible implication in plant-pathogen interaction, including cell wall biogenesis and degradation, sucrose metabolism, secondary metabolism, hormone biosynthesis and signalling, amino acid and lipid metabolism, while down-regulated DEGs were predominantly enriched in ubiquitin proteolysis and oxidative phosphorylation pathways. Our analysis provides novel insight into the molecular pathways that are deregulated during the host-pathogen interaction in Mexican lime trees infected by phytoplasma. The findings can be valuable for unravelling the molecular mechanisms of plant-phytoplasma interactions and can pave the way for engineering lime trees with resistance to witches' broom disease.

RNA-Seq analysis of resistant and susceptible potato varieties during the early stages of potato virus Y infection

Background

Potato virus Y (PVY) is one of the most important plant viruses affecting potato production. The interactions between potato and PVY are complex and the outcome of the interactions depends on the potato genotype, the PVY strain, and the environmental conditions. A potato cultivar can induce resistance to a specific PVY strain, yet be susceptible to another. How a single potato cultivar responds to PVY in both compatible and incompatible interactions is not clear.

Results

In this study, we used RNA-sequencing (RNA-Seq) to investigate and compare the transcriptional changes in leaves of potato upon inoculation with PVY. We used two potato varieties: Premier Russet, which is resistant to the PVY strain O (PVY^O) but susceptible to the strain NTN (PVY^NTN), and Russet Burbank, which is susceptible to all PVY strains that have been tested. Leaves were inoculated with PVY^O or PVY^NTN, and samples were collected 4 and 10 h post inoculation (hpi). A larger number of differentially expressed (DE) genes were found in the compatible reactions compared to the incompatible reaction. For all treatments, the majority of DE genes were down-regulated at 4 hpi and up-regulated at 10 hpi. Gene Ontology enrichment analysis showed enrichment of the biological process GO term “Photosynthesis, light harvesting” specifically in PVY^O-inoculated Premier Russet leaves, while the GO term “nucleosome assembly” was largely overrepresented in PVY^NTN-inoculated Premier Russet leaves and PVY^O-inoculated Russet Burbank leaves but not in PVY^O-inoculated Premier Russet leaves. Fewer genes were DE over 4-fold in the incompatible reaction compared to the compatible reactions. Amongst these, five genes were DE only in PVY^O-inoculated Premier Russet leaves, and all five were down-regulated. These genes are predicted to encode for a putative ABC transporter, a MYC2 transcription factor, a VQ-motif containing protein, a non-specific lipid-transfer protein, and a xyloglucan endotransglucosylase-hydroxylase.

Conclusions

Our results show that the incompatible and compatible reactions in Premier Russet shared more similarities, in particular during the initial response, than the compatible reactions in the two different hosts. Our results identify potential key processes and genes that determine the fate of the reaction, compatible or incompatible, between PVY and its host.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1666-2) contains supplementary material, which is available to authorized users.

Transcriptome responses to Ralstonia solanacearum infection in the roots of the wild potato Solanum commersonii

BACKGROUND

Solanum commersonii is a wild potato species that exhibits high tolerance to both biotic and abiotic stresses and has been used as a source of genes for introgression into cultivated potato. Among the interesting features of S. commersonii is resistance to the bacterial wilt caused by Ralstonia solanacearum, one of the most devastating bacterial diseases of crops.

RESULTS

In this study, we used deep sequencing of S. commersonii RNA (RNA-seq) to analyze the below-ground plant transcriptional responses to R. solanacearum. While a majority of S. commersonii RNA-seq reads could be aligned to the Solanum tuberosum Group Phureja DM reference genome sequence, we identified 2,978 S. commersonii novel transcripts through assembly of unaligned S. commersonii RNA-seq reads. We also used RNA-seq to study gene expression in pathogen-challenged roots of S. commersonii accessions resistant (F118) and susceptible (F97) to the pathogen. Expression profiles obtained from read mapping to the S. tuberosum reference genome and the S. commersonii novel transcripts revealed a differential response to the pathogen in the two accessions, with 221 (F118) and 644 (F97) differentially expressed genes including S. commersonii novel transcripts in the resistant and susceptible genotypes. Interestingly, 22.6% of the F118 and 12.8% of the F97 differentially expressed genes had been previously identified as responsive to biotic stresses and half of those up-regulated in both accessions had been involved in plant pathogen responses. Finally, we compared two different methods to eliminate ribosomal RNA from the plant RNA samples in order to allow dual mapping of RNAseq reads to the host and pathogen genomes and provide insights on the advantages and limitations of each technique.

CONCLUSIONS

Our work catalogues the S. commersonii transcriptome and strengthens the notion that this species encodes specific genes that are differentially expressed to respond to bacterial wilt. In addition, a high proportion of S. commersonii-specific transcripts were altered by R. solanacearum only in F118 accession, while phythormone-related genes were highly induced in F97, suggesting a markedly different response to the pathogen in the two plant accessions studied.

Uncovering the genetic architecture of Colletotrichum lindemuthianum resistance through QTL mapping and epistatic interaction analysis in common bean

Colletotrichum lindemuthianum is a hemibiotrophic fungal pathogen that causes anthracnose disease in common bean. Despite the genetics of anthracnose resistance has been studied for a long time, few quantitative trait loci (QTLs) studies have been conducted on this species. The present work examines the genetic basis of quantitative resistance to races 23 and 1545 of C. lindemuthianum in different organs (stem, leaf and petiole). A population of 185 recombinant inbred lines (RIL) derived from the cross PMB0225 × PHA1037 was evaluated for anthracnose resistance under natural and artificial photoperiod growth conditions. Using multi-environment QTL mapping approach, 10 and 16 main effect QTLs were identified for resistance to anthracnose races 23 and 1545, respectively. The homologous genomic regions corresponding to 17 of the 26 main effect QTLs detected were positive for the presence of resistance-associated gene cluster encoding nucleotide-binding and leucine-rich repeat (NL) proteins. Among them, it is worth noting that the main effect QTLs detected on linkage group 05 for resistance to race 1545 in stem, petiole and leaf were located within a 1.2 Mb region. The NL gene Phvul.005G117900 is located in this region, which can be considered an important candidate gene for the non-organ-specific QTL identified here. Furthermore, a total of 39 epistatic QTL (E-QTLs) (21 for resistance to race 23 and 18 for resistance to race 1545) involved in 20 epistatic interactions (eleven and nine interactions for resistance to races 23 and 1545, respectively) were identified. None of the main and epistatic QTLs detected displayed significant environment interaction effects. The present research provides essential information not only for the better understanding of the plant-pathogen interaction but also for the application of genomic assisted breeding for anthracnose resistance improvement in common bean through application of marker-assisted selection (MAS).

Analysis of gene expression changes in peach leaves in response to Plum pox virus infection using RNA-Seq

Differences in gene expression were studied after Plum pox virus (PPV, sharka disease) infection in peach GF305 leaves with and without sharka symptoms using RNA-Seq. For each sample, more than 80% of 100-nucleotide paired-end (PE) Illumina reads were aligned on the peach reference genome. In the symptomatic sample, a significant proportion of reads were mapped to PPV reference genomes (1.04% compared with 0.00002% in non-symptomatic leaves), allowing for the ultra-deep assembly of the complete genome of the PPV isolate used (9775 nucleotides, missing only 11 nucleotides at the 5' genome end). In addition, significant alternative splicing events were detected in 359 genes and 12 990 single nucleotide polymorphisms (SNPs) were identified, 425 of which could be annotated. Gene ontology annotation revealed that the high-ranking mRNA target genes associated with the expression of sharka symptoms are mainly related to the response to biotic stimuli, to lipid and carbohydrate metabolism and to the negative regulation of catalytic activity. A greater number of differentially expressed genes were observed in the early asymptomatic phase of PPV infection in comparison with the symptomatic phase. These early infection events were associated with the induction of genes related to pathogen resistance, such as jasmonic acid, chitinases, cytokinin glucosyl transferases and Lys-M proteins. Once the virus had accumulated, the overexpression of Dicer protein 2a genes suggested a gene silencing plant response that was suppressed by the virus HCPro and P1 proteins. These results illustrate the dynamic nature of the peach-PPV interaction at the transcriptome level and confirm that sharka symptom expression is a complex process that can be understood on the basis of changes in plant gene expression.

Comparative Transcriptome Analysis of Resistant and Susceptible Tomato Lines in Response to Infection by Xanthomonas perforans Race T3

Bacterial spot, incited by several Xanthomonas sp., is a serious disease in tomato (Solanum lycopersicum L.). Although genetics of resistance has been widely investigated, the interactions between the pathogen and tomato plants remain unclear. In this study, tanscriptomes of X. perforans race T3 infected tomato lines were compared to those of controls. An average of 7 million reads were generated with approximately 21,526 genes mapped in each sample post-inoculation at 6 h (6 HPI) and 6 days (6 DPI) using RNA-sequencing technology. Overall, the numbers of differentially expressed genes (DEGs) were higher in the resistant tomato line PI 114490 than in the susceptible line OH 88119, and the numbers of DEGs were higher at 6 DPI than at 6 HPI. Fewer genes (78 in PI 114490 and 15 in OH 88119) were up-regulated and most DEGs were down-regulated, suggesting that the inducible defense response might not be fully activated at 6 HPI. Accumulation expression levels of 326 co-up regulated genes in both tomato lines at 6 DPI might be involved in basal defense, while the specific and strongly induced genes at 6 DPI might be correlated with the resistance in PI 114490. Most DEGs were involved in plant hormone signal transduction, plant-pathogen interaction and phenylalanine metabolism, and the genes significantly up-regulated in PI 114490 at 6 DPI were associated with defense response pathways. DEGs containing NBS-LRR domain or defense-related WRKY transcription factors were also identified. The results will provide a valuable resource for understanding the interactions between X. perforans and tomato plants.

AtROP1 negatively regulates potato resistance to Phytophthora infestans via NADPH oxidase-mediated accumulation of H2O2

BACKGROUND

Small GTPases are monomeric guanine nucleotide-binding proteins. In plants, ROPs regulate plant cell polarity, plant cell differentiation and development as well as biotic and abiotic stress signaling pathways.

RESULTS

We report the subcellular localization of the AtRop1 protein at the plasma membrane in tobacco epidermal cells using GFP fusions. Additionally, transient and stable expression of a dominant negative form (DN) of the Arabidopsis AtRop1 in potato led to H2O2 accumulation associated with the reduced development of Phytophthora infestans Montagne de Bary and smaller lesions on infected potato leaves. The expression of the Strboh-D gene, a NADPH oxidase homologue in potato, was analyzed by RT-PCR. Expression of this gene was maintained in DN-AtRop1 transgenic plants after infection with P. infestans. In transgenic potato lines, the transcript levels of salicylic acid (SA) and jasmonic acid (JA) marker genes (Npr1 and Lox, respectively) were analyzed. The Lox gene was induced dramatically whereas expression of Npr1, a gene up-regulated by SA, decreased slightly in DN-AtRop1 transgenic plants after infection with P. infestans.

CONCLUSIONS

In conclusion, our results indicate that DN-AtROP1 affects potato resistance to P. infestans. This is associated with increased NADPH oxidase-mediated H2O2 production and JA signaling.

Fine mapping of Rcr1 and analyses of its effect on transcriptome patterns during infection by Plasmodiophora brassicae

BACKGROUND

The protist Plasmodiophora brassicae is a biotrophic soil-borne pathogen that causes clubroot on Brassica crops worldwide. Clubroot disease is a serious threat to the 8 M ha of canola (Brassica napus) grown annually in western Canada. While host resistance is the key to clubroot management, sources of resistance are limited.

RESULTS

To identify new sources of clubroot resistance (CR), we fine mapped a CR gene (Rcr1) from B. rapa ssp. chinensis to the region between 24.26 Mb and 24.50 Mb on the linkage group A03, with several closely linked markers identified. Transcriptome analysis was conducted using RNA sequencing on a segregating F1 population inoculated with P. brassicae, with 2,212 differentially expressed genes (DEGs) identified between plants carrying and not carrying Rcr1. Functional annotation of these DEGs showed that several defense-related biological processes, including signaling and metabolism of jasmonate and ethylene, defensive deposition of callose and biosynthesis of indole-containing compounds, were up-regulated significantly in plants carrying Rcr1 while genes involved in salicylic acid metabolic and signaling pathways were generally not elevated. Several DEGs involved in metabolism potentially related to clubroot symptom development, including auxin biosynthesis and cell growth/development, showed significantly lower expression in plants carrying Rcr1.

CONCLUSION

The CR gene Rcr1 and closely linked markers will be highly useful for breeding new resistant canola cultivars. The identification of DEGs between inoculated plants carrying and not carrying Rcr1 is an important step towards understanding of specific metabolic/signaling pathways in clubroot resistance mediated by Rcr1. This information may help judicious use of CR genes with complementary resistance mechanisms for durable clubroot resistance.

The dihydrolipoyl acyltransferase gene BCE2 participates in basal resistance against Phytophthora infestans in potato and Nicotiana benthamiana

Dihydrolipoyl acyltransferase (EC 2.3.1.12), a branched-chain α-ketoacid dehydrogenase E2 subunit (BCE2), catalyzes the transfer of the acyl group from the lipoyl moiety to coenzyme A. However, the role of BCE2 responding to biotic stress in plant is not clear. In this study, we cloned and characterized a BCE2 gene from potato, namely StBCE2, which was previously suggested to be involved in Phytophthora infestans-potato interaction. We found that the expression of StBCE2 was strongly induced by both P. infestans isolate HB09-14-2 and salicylic acid. Besides, when the homolog of StBCE2 in Nicotiana benthamiana named NbBCE2 was silenced, plants showed increased susceptibility to P. infestans and reduced accumulation of hydrogen peroxide (H2O2). Furthermore, we found that a marker gene NbrbohB involved in the production of reactive oxygen species, was also suppressed in NbBCE2-silenced plants. However, silencing of NbBCE2 had no significant effect on the hypersensitive responses trigged by INF1, R3a-AVR3a(KI) pair or Rpi-vnt1.1-AVR-vnt1.1 pair. Our results suggest that BCE2 is associated with the basal resistance to P. infestans by regulating H2O2 production.

Characterization of Withania somnifera leaf transcriptome and expression analysis of pathogenesis-related genes during salicylic acid signaling

Withania somnifera (L.) Dunal is a valued medicinal plant with pharmaceutical applications. The present study was undertaken to analyze the salicylic acid induced leaf transcriptome of W. somnifera. A total of 45.6 million reads were generated and the de novo assembly yielded 73,523 transcript contig with average transcript contig length of 1620 bp. A total of 71,062 transcripts were annotated and 53,424 of them were assigned GO terms. Mapping of transcript contigs to biological pathways revealed presence of 182 pathways. Seventeen genes representing 12 pathogenesis-related (PR) families were mined from the transcriptome data and their pattern of expression post 17 and 36 hours of salicylic acid treatment was documented. The analysis revealed significant up-regulation of all families of PR genes by 36 hours post treatment except WsPR10. The relative fold expression of transcripts ranged from 1 fold to 6,532 fold. The two families of peroxidases including the lignin-forming anionic peroxidase (WsL-PRX) and suberization-associated anionic peroxidase (WsS-PRX) recorded maximum expression of 377 fold and 6532 fold respectively, while the expression of WsPR10 was down-regulated by 14 fold. Additionally, the most stable reference gene for normalization of qRT-PCR data was also identified. The effect of SA on the accumulation of major secondary metabolites of W. somnifera including withanoside V, withaferin A and withanolide A was also analyzed and an increase in content of all the three metabolites were detected. This is the first report on expression patterns of PR genes during salicylic acid signaling in W. somnifera.