Clubroot-Induced Changes in the Root and Rhizosphere Microbiome of Susceptible and Resistant Canola.

Clubroot is a soilborne disease of canola (Brassica napus) and other crucifers caused by the obligate parasite Plasmodiophora brassicae. In western Canada, clubroot is usually managed by planting-resistant cultivars, but the emergence of resistance-breaking pathotypes of P. brassicae represents a major threat to sustainable canola production. The rhizosphere and root contain beneficial microorganisms that can improve plant health. In this study, we evaluated the effect of two P. brassicae isolates (termed A and B) with different levels of virulence on the root and rhizosphere microbiomes of clubroot-resistant and clubroot-susceptible canola. Additionally, potential biocontrol microorganisms were identified based on taxa antagonistic to clubroot. Although both P. brassicae isolates were classified as pathotype 3A, isolate A caused a higher disease severity index in the resistant canola genotype compared with isolate B. Metabarcoding analysis indicated a shift in the bacterial and fungal communities in response to inoculation with either field isolate. Root endophytic bacterial and fungal communities responded to changes in inoculation, isolate type, sampling time, and canola genotype. In contrast, fungal communities associated with the rhizosphere exhibited significant differences between sampling times, while bacterial communities associated with the rhizosphere exhibited low variability.

Seed protein biotyping in Amaranthus species: a tool for rapid identification of weedy amaranths of concern.

BACKGROUND: The Amaranthus genus contains at least 20 weedy and invasive species, including Amaranthus palmeri (palmer's amaranth) and Amaranthus tuberculatus (tall waterhemp), two species of regulatory concern in North America, impacting production and yield in crops like corn, soybean and cotton. Amaranthus tuberculatus is regulated in Canada with limited establishment, while current climate models predict a range expansion of A. palmeri impacting crop growing areas in Ontario, Quebec and Manitoba. Since many Amaranthus species are similar in their morphology, especially at the seed stage, this demands the development of additional methods that can efficiently aid in the detection and identification of these species. Protein biotyping using Matrix-Assisted Laser Desorption Ionization Time of Flight Mass Spectrometry (MALDI-TOF-MS) has been traditionally used to identify microorganism species, races and pathotypes. Major protein fractions extracted from an organism, ionized and run through a biotyper using mass spectrometry, result in protein spectra that represent a fingerprint at the species or lower taxonomic rank, providing an efficient molecular diagnostics method. Here we use a modified protein biotyping protocol to extract major protein fractions from seeds of the family Brassicaceae to test our protocol, and then implemented the standardized approach in seeds from Amaranthus species. We then created a database of Amaranthus protein spectra that can be used to test blind samples for a quick identification of species of concern. RESULTS: We generated a protein spectra database with 16 Amaranthus species and several accessions per species, spanning target species of regulatory concern and species which are phylogenetically related or easily confused at the seed stage due to phenotypic plasticity. Testing of two Amaranthus blind sample seed sets against this database showed accuracies of 100% and 87%, respectively. CONCLUSIONS: Our method is highly efficient in identifying Amaranthus species of regulatory concern. The mismatches between our protein biotyping approach and phenotypic identification of seeds are due to absence of the species in the database or close phylogenetic relationship between the species. While A. palmeri cannot be distinguished from A. watsonii, there is evidence these two species have the same native range and are closely related.

Application of the NanoString nCounter System as an Alternative Method to Investigate Molecular Mechanisms Involved in Host Plant Responses to Plasmodiophora brassicae.

Clubroot, caused by the soilborne pathogen Plasmodiophora brassicae, is an important disease of canola (Brassica napus) and other crucifers. The recent application of RNA sequencing (RNA-seq) technologies to study P. brassicae−host interactions has generated large amounts of gene expression data, improving knowledge of the molecular mechanisms of pathogenesis and host resistance. Quantitative PCR (qPCR) analysis has been widely applied to examine the expression of a limited number of genes and to validate the results of RNA-seq studies, but may not be ideal for analyzing larger suites of target genes or increased sample numbers. Moreover, the need for intermediate steps such as cDNA synthesis may introduce variability that could affect the accuracy of the data generated by qPCR. Here, we report the validation of gene expression data from a previous RNA-seq study of clubroot using the NanoString nCounter System, which achieves efficient gene expression quantification in a fast and simple manner. We first confirm the robustness of the NanoString system by comparing the results with those generated by qPCR and RNA-seq and then discuss the importance of some candidate genes for resistance or susceptibility to P. brassicae in the host. The results show that the expression of genes measured using NanoString have a high correlation with the values obtained using the other two technologies, with R > 0.90 and p < 0.01, and the same expression patterns for most genes. The three methods (qPCR, RNA-seq, and NanoString) were also compared in terms of laboratory procedures, time, and cost. We propose that the NanoString nCounter System is a robust, sensitive, highly reproducible, and simple technology for gene expression analysis. NanoString could become a common alternative to qPCR to validate RNA-seq data or to create panels of genes for use as markers of resistance/susceptibility when plants are challenged with different P. brassicae pathotypes.

Protocol: rhPCR and SNaPshot assays to distinguish Plasmodiophora brassicae pathotype clusters.

BACKGROUND: Clubroot of canola (Brassica napus), caused by the soilborne pathogen Plasmodiophora brassicae, has become a serious threat to canola production in Canada. The deployment of clubroot-resistant (CR) cultivars is the most commonly used management strategy; however, the widespread cultivation of CR canola has resulted in the emergence of new pathotypes of P. brassicae capable of overcoming resistance. Several host differential sets have been reported for pathotype identification, but such testing is time-consuming, labor-intensive, and based on phenotypic classifications. The development of rapid and objective methods that allow for efficient, cost-effective and convenient pathotyping would enable testing of a much larger number of samples in shorter times. The aim of this study was to develop two pathotyping assays, an RNase H2-dependent PCR (rhPCR) assay and a SNaPshot assay, which could quickly differentiate P. brassicae pathotypes. RESULTS: Both assays clearly distinguished between pathotype clusters in a collection of 38 single-spore isolates of P. brassicae. Additional isolates pathotyped from clubbed roots and samples from blind testing also were correctly clustered. The rhPCR assay generated clearly differentiating electrophoretic bands without non-specific amplification. The SNaPshot assay was able to detect down to a 10% relative allelic proportion in a 10:90 template mixture with both single-spore isolates and field isolates when evaluated in a relative abundance test. CONCLUSIONS: This study describes the development of two rapid and sensitive technologies for P. brassicae pathotyping. The high-throughput potential and accuracy of both assays makes them promising as SNP-based pathotype identification tools for clubroot diagnostics. rhPCR is a highly sensitive approach that can be optimized into a quantitative assay, while the main advantages of SNaPshot are its ability to multiplex samples and alleles in a single reaction and the detection of up to four allelic variants per target site.

Candidate Effectors of Plasmodiophora brassicae Pathotype 5X During Infection of Two Brassica napus Genotypes.

Clubroot, caused by Plasmodiophora brassicae, is one of the most important diseases of canola (Brassica napus) in Canada. Disease management relies heavily on planting clubroot resistant (CR) cultivars, but in recent years, new resistance-breaking pathotypes of P. brassicae have emerged. Current efforts against the disease are concentrated in developing host resistance using traditional genetic breeding, omics and molecular biology. However, because of its obligate biotrophic nature, limited resources have been dedicated to investigating molecular mechanisms of pathogenic infection. We previously performed a transcriptomic study with the cultivar resistance-breaking pathotype 5X on two B. napus hosts presenting contrasting resistance/susceptibility, where we evaluated the mechanisms of host response. Since cultivar-pathotype interactions are very specific, and pathotype 5X is one of the most relevant resistance-breaking pathotypes in Canada, in this study, we analyze the expression of genes encoding putative secreted proteins from this pathotype, predicted using a bioinformatics pipeline, protein modeling and orthologous comparisons with effectors from other pathosystems. While host responses were found to differ markedly in our previous study, many common effectors are found in the pathogen while infecting both hosts, and the gene response among biological pathogen replicates seems more consistent in the effectors associated with the susceptible interaction, especially at 21 days after inoculation. The predicted effectors indicate the predominance of proteins with interacting domains (e.g., ankyrin), and genes bearing kinase and NUDIX domains, but also proteins with protective action against reactive oxygen species from the host. Many of these genes confirm previous predictions from other clubroot studies. A benzoic acid/SA methyltransferase (BSMT), which methylates SA to render it inactive, showed high levels of expression in the interactions with both hosts. Interestingly, our data indicate that E3 ubiquitin proteasome elements are also potentially involved in pathogenesis. Finally, a gene with similarity to indole-3-acetaldehyde dehydrogenase is a promising candidate effector because of its involvement in indole acetic acid synthesis, since auxin is one of the major players in clubroot development.

Current and Future Pathotyping Platforms for Plasmodiophora brassicae in Canada.

Clubroot, caused by Plasmodiophora brassicae, is one of the most detrimental threats to crucifers worldwide and has emerged as an important disease of canola (Brassica napus) in Canada. At present, pathotypes are distinguished phenotypically by their virulence patterns on host differential sets, including the systems of Williams, Somé et al., the European Clubroot Differential set, and most recently the Canadian Clubroot Differential set and the Sinitic Clubroot Differential set. Although these are frequently used because of their simplicity of application, they are time-consuming, labor-intensive, and can lack sensitivity. Early, preventative pathotype detection is imperative to maximize productivity and promote sustainable crop production. The decreased turnaround time and increased sensitivity and specificity of genotypic pathotyping will be valuable for the development of integrated clubroot management plans, and interest in molecular techniques to complement phenotypic methods is increasing. This review provides a synopsis of current and future molecular pathotyping platforms for P. brassicae and aims to provide information on techniques that may be most suitable for the development of rapid, reliable, and cost-effective pathotyping assays.

Comparative Transcriptome Analysis of Rutabaga (Brassica napus) Cultivars Indicates Activation of Salicylic Acid and Ethylene-Mediated Defenses in Response to Plasmodiophora brassicae.

Clubroot, caused by Plasmodiophora brassicae Woronin, is an important soilborne disease of Brassica napus L. and other crucifers. To improve understanding of the mechanisms of resistance and pathogenesis in the clubroot pathosystem, the rutabaga (B. napus subsp. rapifera Metzg) cultivars 'Wilhelmsburger' (resistant) and 'Laurentian' (susceptible) were inoculated with P. brassicae pathotype 3A and their transcriptomes were analyzed at 7, 14, and 21 days after inoculation (dai) by RNA sequencing (RNA-seq). Thousands of transcripts with significant changes in expression were identified in each host at each time-point in inoculated vs. non-inoculated plants. Molecular responses at 7 and 14 dai supported clear differences in the clubroot response mechanisms of the two genotypes. Both the resistant and the susceptible cultivars activated receptor-like protein (RLP) genes, resistance (R) genes, and genes involved in salicylic acid (SA) signaling as clubroot defense mechanisms. In addition, genes related to calcium signaling and genes encoding leucine-rich repeat (LRR) receptor kinases, the respiratory burst oxidase homolog (RBOH) protein, and transcription factors such as WRKYs, ethylene responsive factors, and basic leucine zippers (bZIPs), appeared to be upregulated in 'Wilhelmsburger' to restrict P. brassicae development. Some of these genes are essential components of molecular defenses, including ethylene (ET) signaling and the oxidative burst. Our study highlights the importance of activation of genes associated with SA- and ET-mediated responses in the resistant cultivar. A set of candidate genes showing contrasting patterns of expression between the resistant and susceptible cultivars was identified and includes potential targets for further study and validation through approaches such as gene editing.

Response of Brassica napus to Plasmodiophora brassicae Involves Salicylic Acid-Mediated Immunity: An RNA-Seq-Based Study.

Clubroot, caused by the obligate parasite Plasmodiophora brassicae, is an important disease of the Brassicaceae and poses a significant threat to the $26.7 billion canola/oilseed rape (Brassica napus) industry in western Canada. While clubroot is managed most effectively by planting resistant host varieties, new pathotypes of P. brassicae have emerged recently that can overcome this resistance. Whole genome analyses provide both a toolbox and a systemic view of molecular mechanisms in host-pathogen interactions, which can be used to design new breeding strategies to increase P. brassicae resistance. We used RNA-seq to evaluate differential gene expression at 7, 14 and 21 days after inoculation (dai) of two B. napus genotypes with differential responses to P. brassicae pathotype 5X. Gall development was evident at 14 dai in the susceptible genotype (the oilseed rape 'Brutor'), while gall development in the resistant genotype (the rutabaga (B. napus) 'Laurentian') was limited and not visible until 21 dai. Immune responses were better sustained through the time-course in 'Laurentian', and numerous genes from immune-related functional categories were associated with salicylic acid (SA)-mediated responses. Jasmonic acid (JA)-mediated responses seemed to be mostly inhibited, especially in the resistant genotype. The upregulation of standard defense-related proteins, like chitinases and thaumatins, was evident in 'Laurentian'. The enrichment, in both host genotypes, of functional categories for syncytium formation and response to nematodes indicated that cell enlargement during P. brassicae infection, and the metabolic processes therein, share similarities with the response to infection by nematodes that produce similar anatomical symptoms. An analysis of shared genes between the two genotypes at different time-points, confirmed that the nematode-like responses occurred earlier for 'Brutor', along with cell metabolism and growth changes. Additionally, the susceptible cultivar turned off defense mechanisms earlier than 'Laurentian'. Collectively, this study showed the importance of SA in triggering immune responses and suggested some key resistance and susceptibility factors that can be used in future studies for resistance breeding through gene-editing approaches.

LTR-retrotransposons in plants: Engines of evolution.

LTR retrotransposons are the most abundant group of transposable elements (TEs) in plants. These elements can fall inside or close to genes, and therefore influence their expression and evolution. This review aims to examine how LTR retrotransposons, especially Ty1-copia elements, mediate gene regulation and evolution. Various stimuli, including polyploidization and biotic and abiotic elicitors, result in the transcription and movement of these retrotransposons, and can facilitate adaptation. The presence of cis-regulatory motifs in the LTRs are central to their stress-mediated responses and are shared with host stress-responsive genes, showing a complex evolutionary history in which TEs provide new regulatory units to genes. The presence of retrotransposon remnants in genes that are necessary for normal gene function, demonstrates the importance of exaptation and co-option, and is also a consequence of the abundance of these elements in plant genomes. Furthermore, insertions of LTR retrotransposons in and around genes provide potential for alternative splicing, epigenetic control, transduction, duplication and recombination. These characteristics can become an active part of the evolution of gene families as in the case of resistance genes (R-genes). The character of TEs as exclusively selfish is now being re-evaluated. Since genome-wide reprogramming via TEs is a long evolutionary process, the changes we can examine are case-specific and their fitness advantage may not be evident until TE-derived motifs and domains have been completely co-opted and fixed. Nevertheless, the presence of LTR retrotransposons inside genes and as part of gene promoter regions is consistent with their roles as engines of plant genome evolution.

Ty1-copia elements reveal diverse insertion sites linked to polymorphisms among flax (Linum usitatissimum L.) accessions.

BACKGROUND: Initial characterization of the flax genome showed that Ty1-copia retrotransposons are abundant, with several members being recently inserted, and in close association with genes. Recent insertions indicate a potential for ongoing transpositional activity that can create genomic diversity among accessions, cultivars or varieties. The polymorphisms generated constitute a good source of molecular markers that may be associated with phenotype if the insertions alter gene activity. Flax, where accessions are bred mainly for seed nutritional properties or for fibers, constitutes a good model for studying the relationship of transpositional activity with diversification and breeding. In this study, we estimated copy number and used a type of transposon display known as Sequence-Specific Amplification Polymorphisms (SSAPs), to characterize six families of Ty1-copia elements across 14 flax accessions. Polymorphic insertion sites were sequenced to find insertions that could potentially alter gene expression, and a preliminary test was performed with selected genes bearing transposable element (TE) insertions. RESULTS: Quantification of six families of Ty1-copia elements indicated different abundances among TE families and between flax accessions, which suggested diverse transpositional histories. SSAPs showed a high level of polymorphism in most of the evaluated retrotransposon families, with a trend towards higher levels of polymorphism in low-copy number families. Ty1-copia insertion polymorphisms among cultivars allowed a general distinction between oil and fiber types, and between spring and winter types, demonstrating their utility in diversity studies. Characterization of polymorphic insertions revealed an overwhelming association with genes, with insertions disrupting exons, introns or within 1 kb of coding regions. A preliminary test on the potential transcriptional disruption by TEs of four selected genes evaluated in three different tissues, showed one case of significant impact of the insertion on gene expression. CONCLUSIONS: We demonstrated that specific Ty1-copia families have been active since breeding commenced in flax. The retrotransposon-derived polymorphism can be used to separate flax types, and the close association of many insertions with genes defines a good source of potential mutations that could be associated with phenotypic changes, resulting in diversification processes.

RNA-seq Transcriptome Response of Flax (Linum usitatissimum L.) to the Pathogenic Fungus Fusarium oxysporum f. sp. lini.

Fusarium oxysporum f. sp. lini is a hemibiotrophic fungus that causes wilt in flax. Along with rust, fusarium wilt has become an important factor in flax production worldwide. Resistant flax cultivars have been used to manage the disease, but the resistance varies, depending on the interactions between specific cultivars and isolates of the pathogen. This interaction has a strong molecular basis, but no genomic information is available on how the plant responds to attempted infection, to inform breeding programs on potential candidate genes to evaluate or improve resistance across cultivars. In the current study, disease progression in two flax cultivars [Crop Development Center (CDC) Bethune and Lutea], showed earlier disease symptoms and higher susceptibility in the later cultivar. Chitinase gene expression was also divergent and demonstrated and earlier molecular response in Lutea. The most resistant cultivar (CDC Bethune) was used for a full RNA-seq transcriptome study through a time course at 2, 4, 8, and 18 days post-inoculation (DPI). While over 100 genes were significantly differentially expressed at both 4 and 8 DPI, the broadest deployment of plant defense responses was evident at 18 DPI with transcripts of more than 1,000 genes responding to the treatment. These genes evidenced a reception and transduction of pathogen signals, a large transcriptional reprogramming, induction of hormone signaling, activation of pathogenesis-related genes, and changes in secondary metabolism. Among these, several key genes that consistently appear in studies of plant-pathogen interactions, had increased transcript abundance in our study, and constitute suitable candidates for resistance breeding programs. These included: an induced RPMI-induced protein kinase; transcription factors WRKY3, WRKY70, WRKY75, MYB113, and MYB108; the ethylene response factors ERF1 and ERF14; two genes involved in auxin/glucosinolate precursor synthesis (CYP79B2 and CYP79B3); the flavonoid-related enzymes chalcone synthase, dihydroflavonol reductase and multiple anthocyanidin synthases; and a peroxidase implicated in lignin formation (PRX52). Additionally, regulation of some genes indicated potential pathogen manipulation to facilitate infection; these included four disease resistance proteins that were repressed, indole acetic acid amido/amino hydrolases which were upregulated, activated expansins and glucanases, amino acid transporters and aquaporins, and finally, repression of major latex proteins.

Ion Torrent sequencing as a tool for mutation discovery in the flax (Linum usitatissimum L.) genome.

BACKGROUND: Detection of induced mutations is valuable for inferring gene function and for developing novel germplasm for crop improvement. Many reverse genetics approaches have been developed to identify mutations in genes of interest within a mutagenized population, including some approaches that rely on next-generation sequencing (e.g. exome capture, whole genome resequencing). As an alternative to these genome or exome-scale methods, we sought to develop a scalable and efficient method for detection of induced mutations that could be applied to a small number of target genes, using Ion Torrent technology. We developed this method in flax (Linum usitatissimum), to demonstrate its utility in a crop species. RESULTS: We used an amplicon-based approach in which DNA samples from an ethyl methanesulfonate (EMS)-mutagenized population were pooled and used as template in PCR reactions to amplify a region of each gene of interest. Barcodes were incorporated during PCR, and the pooled amplicons were sequenced using an Ion Torrent PGM. A pilot experiment with known SNPs showed that they could be detected at a frequency > 0.3% within the pools. We then selected eight genes for which we wanted to discover novel mutations, and applied our approach to screen 768 individuals from the EMS population, using either the Ion 314 or Ion 316 chips. Out of 29 potential mutations identified after processing the NGS reads, 16 mutations were confirmed using Sanger sequencing. CONCLUSIONS: The methodology presented here demonstrates the utility of Ion Torrent technology in detecting mutation variants in specific genome regions for large populations of a species such as flax. The methodology could be scaled-up to test >100 genes using the higher capacity chips now available from Ion Torrent.

Integrated transcriptomic and proteomic profiling of white spruce stems during the transition from active growth to dormancy.

In the autumn, stems of woody perennials such as forest trees undergo a transition from active growth to dormancy. We used microarray transcriptomic profiling in combination with a proteomics analysis to elucidate processes that occur during this growth-to-dormancy transition in a conifer, white spruce (Picea glauca[Moench] Voss). Several differentially expressed genes were likely associated with the developmental transition that occurs during growth cessation in the cambial zone and the concomitant completion of cell maturation in vascular tissues. Genes encoding for cell wall and membrane biosynthetic enzymes showed transcript abundance patterns consistent with completion of cell maturation, and also of cell wall and membrane modifications potentially enabling cells to withstand the harsh conditions of winter. Several differentially expressed genes were identified that encoded putative regulators of cambial activity, cell development and of the photoperiodic pathway. Reconfiguration of carbon allocation figured centrally in the tree's overwintering preparations. For example, genes associated with carbon-based defences such as terpenoids were down-regulated, while many genes associated with protein-based defences and other stress mitigation mechanisms were up-regulated. Several of these correspond to proteins that were accumulated during the growth-to-dormancy transition, emphasizing the importance of stress protection in the tree's adaptive response to overwintering.