Transcriptome analysis of the Brassica napus-Leptosphaeria maculans pathosystem identifies receptor, signaling and structural genes underlying plant resistance

Leptosphaeria maculans isolates with variations in AvrLm1 and AvrLm4 effector genes induce differences in defence responses but not in resistance phenotypes in cultivars carrying the Rlm7 gene

BACKGROUND

The phoma stem canker pathogen Leptosphaeria maculans is one of the most widespread and devastating pathogens of oilseed rape (Brassica napus) in the world. Pathogen colonization is stopped by an interaction of a pathogen Avr effector gene with the corresponding host resistance (R) gene. While molecular mechanisms of this gene-for-gene interaction are being elucidated, understanding of effector function remains limited. The purpose of this study was to determine the action of L. maculans effector (AvrLm) genes on incompatible interactions triggered by B. napus noncorresponding R (Rlm) genes. Specifically, effects of AvrLm4-7 and AvrLm1 on Rlm7-mediated resistance were studied.

RESULTS

Although there was no major effect on symptom expression, induction of defence genes (e.g. PR1) and accumulation of reactive oxygen species was reduced when B. napus cv. Excel carrying Rlm7 was challenged with a L. maculans isolate containing AvrLm1 and a point mutation in AvrLm4-7 (AvrLm1, avrLm4-AvrLm7) compared to an isolate lacking AvrLm1 (avrLm1, AvrLm4-AvrLm7). AvrLm7-containing isolates, isogenic for presence or absence of AvrLm1, elicited similar symptoms on hosts with or without Rlm7, confirming results obtained with more genetically diverse isolates.

CONCLUSION

Careful phenotypic examination of isogenic L. maculans isolates and B. napus introgression lines demonstrated a lack of effect of AvrLm1 on Rlm7-mediated resistance despite an apparent alteration of the Rlm7-dependent defence response using more diverse fungal isolates with differences in AvrLm1 and AvrLm4. As deployment of Rlm7 resistance in crop cultivars increases, other effectors need to be monitored because they may alter the predominance of AvrLm7. © 2023 The Authors. Pest Management Science published by John Wiley & Sons Ltd on behalf of Society of Chemical Industry.

Transcriptomics of temperature-sensitive R gene-mediated resistance identifies a WAKL10 protein interaction network

Understanding temperature-sensitivity of R gene-mediated resistance against apoplastic pathogens is important for sustainable food production in the face of global warming. Here, we show that resistance of Brassica napus cotyledons against Leptosphaeria maculans was temperature-sensitive in introgression line Topas-Rlm7 but temperature-resilient in Topas-Rlm4. A set of 1,646 host genes was differentially expressed in Topas-Rlm4 and Topas-Rlm7 in response to temperature. Amongst these were three WAKL10 genes, including BnaA07g20220D, representing the temperature-sensitive Rlm7-1 allele and Rlm4. Network analysis identified a WAKL10 protein interaction cluster specifically for Topas-Rlm7 at 25 °C. Diffusion analysis of the Topas-Rlm4 network identified WRKY22 as a putative regulatory target of the ESCRT-III complex-associated protein VPS60.1, which belongs to the WAKL10 protein interaction community. Combined enrichment analysis of gene ontology terms considering gene expression and network data linked vesicle-mediated transport to defence. Thus, dysregulation of effector-triggered defence in Topas-Rlm7 disrupts vesicle-associated resistance against the apoplastic pathogen L. maculans.

Comparative transcriptome profiling reveals differential defense responses among Alternaria brassicicola resistant Sinapis alba and susceptible Brassica rapa

Alternaria blight is a devastating disease that causes significant crop losses in oilseed Brassicas every year. Adoption of conventional breeding to generate disease-resistant varieties has so far been unsuccessful due to the lack of suitable resistant source germplasms of cultivated Brassica spp. A thorough understanding of the molecular basis of resistance, as well as the identification of defense-related genes involved in resistance responses in closely related wild germplasms, would substantially aid in disease management. In the current study, a comparative transcriptome profiling was performed using Illumina based RNA-seq to detect differentially expressed genes (DEGs) specifically modulated in response to Alternaria brassicicola infection in resistant Sinapis alba, a close relative of Brassicas, and the highly susceptible Brassica rapa. The analysis revealed that, at 48 hpi (hours post inoculation), 3396 genes were upregulated and 23239 were downregulated, whereas at 72 hpi, 4023 genes were upregulated and 21116 were downregulated. Furthermore, a large number of defense response genes were detected to be specifically regulated as a result of Alternaria infection. The transcriptome data was validated using qPCR-based expression profiling for selected defense-related DEGs, that revealed significantly higher fold change in gene expression in S. alba when compared to B. rapa. Expression of most of the selected genes was elevated across all the time points under study with significantly higher expression towards the later time point of 72 hpi in the resistant germplasm. S. alba activates a stronger defense response reaction against the disease by deploying an array of genes and transcription factors involved in a wide range of biological processes such as pathogen recognition, signal transduction, cell wall modification, antioxidation, transcription regulation, etc. Overall, the study provides new insights on resistance of S. alba against A. brassicicola, which will aid in devising strategies for breeding resistant varieties of oilseed Brassica.

Two-Component System Genes in Brassica napus: Identification, Analysis, and Expression Patterns in Response to Abiotic and Biotic Stresses

The two-component system (TCS), consisting of histidine kinases (HKs), histidine phosphotransfer proteins (HPs) and response regulators (RRs) in eukaryotes, plays pivotal roles in regulating plant growth, development, and responses to environment stimuli. However, the TCS genes were poorly characterized in rapeseed, which is an important tetraploid crop in Brassicaceae. In this work, a total of 182 BnaTCS genes were identified, including 43 HKs, 16 HPs, and 123 RRs, which was more than that in other crops due to segmental duplications during the process of polyploidization. It was significantly different in genetic diversity between the three subfamilies, and some members showed substantial genetic differentiation among the three rapeseed ecotypes. Several hormone- and stress-responsive cis-elements were identified in the putative promoter regions of BnaTCS genes. Furthermore, the expression of BnaTCS genes under abiotic stresses, exogenous phytohormone, and biotic stresses was analyzed, and numerous candidate stress-responsive genes were screened out. Meanwhile, using a natural population with 505 B. napus accessions, we explored the genetic effects of BnaTCS genes on salt tolerance by association mapping analysis and detected some significant association SNPs/genes. The result will help to further understand the functions of TCS genes in the developmental and stress tolerance improvement in B. napus.

Novel gene loci associated with susceptibility or cryptic quantitative resistance to Pyrenopeziza brassicae in Brassica napus

Quantitative disease resistance (QDR) controls the association of the light leaf spot pathogen with Brassica napus; four QDR loci that were in linkage disequilibrium and eight gene expression markers were identified. Quantitative disease resistance (QDR) can provide durable control of pathogens in crops in contrast to resistance (R) gene-mediated resistance which can break down due to pathogen evolution. QDR is therefore a desirable trait in crop improvement, but little is known about the causative genes, and so it is difficult to incorporate into breeding programmes. Light leaf spot, caused by Pyrenopeziza brassicae, is an important disease of oilseed rape (canola, Brassica napus). To identify new QDR gene loci, we used a high-throughput screening pathosystem with P. brassicae on 195 lines of B. napus combined with an association transcriptomics platform. We show that all resistance against P. brassicae was associated with QDR and not R gene-mediated. We used genome-wide association analysis with an improved B. napus population structure to reveal four gene loci significantly (P = 0.0001) associated with QDR in regions showing linkage disequilibrium. On chromosome A09, enhanced resistance was associated with heterozygosity for a cytochrome P450 gene co-localising with a previously described locus for seed glucosinolate content. In addition, eight significant gene expression markers with a false discovery rate of 0.001 were associated with QDR against P. brassicae. For seven of these, expression was positively correlated with resistance, whereas for one, a HXXXD-type acyl-transferase, negative correlation indicated a potential susceptibility gene. The study identifies novel QDR loci for susceptibility and resistance, including novel cryptic QDR genes associated with heterozygosity, that will inform future crop improvement.

Identification of candidate genes for LepR1 resistance against Leptosphaeria maculans in Brassica napus

Utilising resistance (R) genes, such as LepR1, against Leptosphaeria maculans, the causal agent of blackleg in canola (Brassica napus), could help manage the disease in the field and increase crop yield. Here we present a genome wide association study (GWAS) in B. napus to identify LepR1 candidate genes. Disease phenotyping of 104 B. napus genotypes revealed 30 resistant and 74 susceptible lines. Whole genome re-sequencing of these cultivars yielded over 3 million high quality single nucleotide polymorphisms (SNPs). GWAS in mixed linear model (MLM) revealed a total of 2,166 significant SNPs associated with LepR1 resistance. Of these SNPs, 2108 (97%) were found on chromosome A02 of B. napus cv. Darmor bzh v9 with a delineated LepR1_mlm1 QTL at 15.11-26.08 Mb. In LepR1_mlm1, there are 30 resistance gene analogs (RGAs) (13 nucleotide-binding site-leucine rich repeats (NLRs), 12 receptor-like kinases (RLKs), and 5 transmembrane-coiled-coil (TM-CCs)). Sequence analysis of alleles in resistant and susceptible lines was undertaken to identify candidate genes. This research provides insights into blackleg resistance in B. napus and assists identification of the functional LepR1 blackleg resistance gene.

Tissue-specific mRNA profiling of the Brassica napus-Sclerotinia sclerotiorum interaction uncovers novel regulators of plant immunity

White mold is caused by the fungal pathogen Sclerotinia sclerotiorum and leads to rapid and significant loss in plant yield. Among its many brassicaceous hosts, including Brassica napus (canola) and Arabidopsis, the response of individual tissue layers directly at the site of infection has yet to be explored. Using laser microdissection coupled with RNA sequencing, we profiled the epidermis, mesophyll, and vascular leaf tissue layers of B. napus in response to S. sclerotiorum. High-throughput tissue-specific mRNA sequencing increased the total number of detected transcripts compared with whole-leaf assessments and provided novel insight into the conserved and specific roles of ontogenetically distinct leaf tissue layers in response to infection. When subjected to pathogen infection, the epidermis, mesophyll, and vasculature activate both specific and shared gene sets. Putative defense genes identified through transcription factor network analysis were then screened for susceptibility against necrotrophic, hemi-biotrophic, and biotrophic pathogens. Arabidopsis deficient in PR5-like RECEPTOR KINASE (PR5K) mRNA levels were universally susceptible to all pathogens tested and were further characterized to identify putative interacting partners involved in the PR5K signaling pathway. Together, these data provide insight into the complexity of the plant defense response directly at the site of infection.

Molecular Interactions Between Leptosphaeria maculans and Brassica Species

Canola is an important oilseed crop, providing food, feed, and fuel around the world. However, blackleg disease, caused by the ascomycete Leptosphaeria maculans, causes significant yield losses annually. With the recent advances in genomic technologies, the understanding of the Brassica napus-L. maculans interaction has rapidly increased, with numerous Avr and R genes cloned, setting this system up as a model organism for studying plant-pathogen associations. Although the B. napus-L. maculans interaction follows Flor's gene-for-gene hypothesis for qualitative resistance, it also puts some unique spins on the interaction. This review discusses the current status of the host-pathogen interaction and highlights some of the future gaps that need addressing moving forward.

Identification of receptor-like proteins induced by Sclerotinia sclerotiorum in Brassica napus

Heightening the resistance of plants to microbial infection is a widely concerned issue, especially for economical crops. Receptor-like proteins (RLPs), typically with tandem leucine-rich repeats (LRRs) domain, play a crucial role in mediating immune activation, being an indispensable constituent in the first layer of defense. Based on an analysis of orthologs among Brassica rapa, Brassica oleracea, and Brassica napus using Arabidopsis thaliana RLPs as a reference framework, we found that compared to A. thaliana, there were some obvious evolutionary diversities of RLPs among the three Brassicaceae species. BnRLP encoding genes were unevenly distributed on chromosomes, mainly on chrA01, chrA04, chrC03, chrC04, and chrC06. The orthologs of five AtRLPs (AtRLP3, AtRLP10, AtRLP17, AtRLP44, and AtRLP51) were highly conserved, but retrenchment and functional centralization occurred in Brassicaceae RLPs during evolution. The RLP proteins were clustered into 13 subgroups. Ten BnRLPs presented expression specificity between R and S when elicited by Sclerotinia sclerotiorum, which might be fabulous candidates for S. sclerotiorum resistance research.

Arabinogalactan Protein-Like Proteins From Ulva lactuca Activate Immune Responses and Plant Resistance in an Oilseed Crop

Natural compounds isolated from macroalgae are promising, ecofriendly, and multifunctional bioinoculants, which have been tested and used in agriculture. Ulvans, for instance, one of the major polysaccharides present in Ulva spp. cell walls, have been tested for their plant growth-promoting properties as well as their ability to activate plant immune defense, on a large variety of crops. Recently, we have characterized for the first time an arabinogalactan protein-like (AGP-like) from Ulva lactuca, which exhibits several features associated to land plant AGPs. In land plant, AGPs were shown to play a role in several plant biological functions, including cell morphogenesis, reproduction, and plant-microbe interactions. Thus, isolated AGP-like proteins may be good candidates for either the plant growth-promoting properties or the activation of plant immune defense. Here, we have isolated an AGP-like enriched fraction from Ulva lactuca and we have evaluated its ability to (i) protect oilseed rape (Brassica napus) cotyledons against Leptosphaeria maculans, and (ii) its ability to activate immune responses. Preventive application of the Ulva AGP-like enriched fraction on oilseed rape, followed by cotyledon inoculation with the fungal hemibiotroph L. maculans, resulted in a major reduction of infection propagation. The noticed reduction correlated with an accumulation of H₂O₂ in treated cotyledons and with the activation of SA and ET signaling pathways in oilseed rape cotyledons. In parallel, an ulvan was also isolated from Ulva lactuca. Preventive application of ulvan also enhanced plant resistance against L. maculans. Surprisingly, reduction of infection severity was only observed at high concentration of ulvan. Here, no such significant changes in gene expression and H₂O₂ production were observed. Together, this study indicates that U. lactuca AGP-like glycoproteins exhibit promising elicitor activity and that plant eliciting properties of Ulva extract, might result not only from an ulvan-originated eliciting activities, but also AGP-like originated.

Less Virulent Leptosphaeria biglobosa Immunizes the Canola Plant to Resist Highly Virulent L. maculans, the Blackleg Pathogen

Leptosphaeria biglobosa is a less virulent Leptosphaeria spp. that causes blackleg disease in canola. Previous studies from our lab have shown that inoculation with the less virulent L. biglobosa can boost the resistance of canola plants against the highly virulent L. maculans. The objective of this study was to confirm the effectiveness of L. biglobosa as a biocontrol agent against L. maculans utilizing morphology, fluorescence microscopy, gene quantification, and transcriptomic analysis. The in planta development of two Leptosphaeria species inoculated at different time points was assessed using fluorescent protein-tagged isolates which are GFP-tagged L. maculans and DsRed-tagged L. biglobosa. The growth inhibition of L. maculans by pre-and co-inoculated L. biglobosa was supported by no lesion development on cotyledons and no or weak fluorescence protein-tagged mycelia under the confocal microscope. The host defense-related genes, WRKY33, PR1, APX6, and CHI, were upregulated in L. biglobosa inoculated Westar cotyledons compared to L. maculans inoculated cotyledons. The quantification of each pathogen through qPCR assay and gene expressions analysis on host defense-related genes by RT-qPCR confirmed the potential of L. biglobosa “brassicae’ in the management of the blackleg disease pathogen, L. maculans ‘brassicae’, in canola.

Leptosphaeria maculans-Brassica napus Battle: A Comparison of Incompatible vs. Compatible Interactions Using Dual RNASeq

Leptosphaeria maculans causes blackleg disease, which is one of the most destructive diseases of canola (Brassica napus L.). Due to the erosion of the current resistance in B. napus, it is pivotal to introduce new resistant genotypes to the growers. This study evaluated the potential of Rlm7 gene as resistance to its corresponding avirulence AvrLm7 gene is abundant. The Rlm7 line was inoculated with L. maculans isolate with AvrLm7; UMAvr7; and the CRISPR/Cas9 knockout AvrLm7 mutant, umavr7, of the same isolate to cause incompatible and compatible interactions, respectively. Dual RNA-seq showed differential gene expressions in both interactions. High expressions of virulence-related pathogen genes-CAZymes, merops, and effector proteins after 7-dpi in compatible interactions but not in incompatible interaction—confirmed that the pathogen was actively virulent only in compatible interactions. Salicyclic and jasmonic acid biosynthesis and signaling-related genes, defense-related PR1 gene (GSBRNA2T00150001001), and GSBRNA2T00068522001 in the NLR gene family were upregulated starting as early as 1- and 3-dpi in the incompatible interaction and the high upregulation of those genes after 7-dpi in compatible interactions confirmed the early recognition of the pathogen by the host and control it by early activation of host defense mechanisms in the incompatible interaction.

RNA-Seq and Gene Regulatory Network Analyses Uncover Candidate Genes in the Early Defense to Two Hemibiotrophic Colletorichum spp. in Strawberry

Two hemibiotrophic pathogens, Colletotrichum acutatum (Ca) and C. gloeosporioides (Cg), cause anthracnose fruit rot and anthracnose crown rot in strawberry (Fragaria × ananassa Duchesne), respectively. Both Ca and Cg can initially infect through a brief biotrophic phase, which is associated with the production of intracellular primary hyphae that can infect host cells without causing cell death and establishing hemibiotrophic infection (HBI) or quiescent (latent infections) in leaf tissues. The Ca and Cg HBI in nurseries and subsequent distribution of asymptomatic infected transplants to fruit production fields is the major source of anthracnose epidemics in North Carolina. In the absence of complete resistance, strawberry varieties with good fruit quality showing rate-reducing resistance have frequently been used as a source of resistance to Ca and Cg. However, the molecular mechanisms underlying the rate-reducing resistance or susceptibility to Ca and Cg are still unknown. We performed comparative transcriptome analyses to examine how rate-reducing resistant genotype NCS 10-147 and susceptible genotype ‘Chandler’ respond to Ca and Cg and identify molecular events between 0 and 48 h after the pathogen-inoculated and mock-inoculated leaf tissues. Although plant response to both Ca and Cg at the same timepoint was not similar, more genes in the resistant interaction were upregulated at 24 hpi with Ca compared with those at 48 hpi. In contrast, a few genes were upregulated in the resistant interaction at 48 hpi with Cg. Resistance response to both Ca and Cg was associated with upregulation of MLP-like protein 44, LRR receptor-like serine/threonine-protein kinase, and auxin signaling pathway, whereas susceptibility was linked to modulation of the phenylpropanoid pathway. Gene regulatory network inference analysis revealed candidate transcription factors (TFs) such as GATA5 and MYB-10, and their downstream targets were upregulated in resistant interactions. Our results provide valuable insights into transcriptional changes during resistant and susceptible interactions, which can further facilitate assessing candidate genes necessary for resistance to two hemibiotrophic Colletotrichum spp. in strawberry.

Gene expression profiling reveals transcription factor networks and subgenome bias during Brassica napus seed development

We profiled the global gene expression landscape across the reproductive lifecycle of Brassica napus. Comparative analysis of this nascent amphidiploid revealed the contribution of each subgenome to plant reproduction. Whole-genome transcription factor networks identified BZIP11 as a transcriptional regulator of early B. napus seed development. Knockdown of BZIP11 using RNA interference resulted in a similar reduction in gene activity of predicted gene targets, and a reproductive-lethal phenotype. Global mRNA profiling revealed lower accumulation of Cⁿ subgenome transcripts relative to the Aⁿ subgenome. Subgenome-specific transcription factor networks identified distinct transcription factor families enriched in each of the Aⁿ and Cⁿ subgenomes early in seed development. Analysis of laser-microdissected seed subregions further reveal subgenome expression dynamics in the embryo, endosperm and seed coat of early stage seeds. Transcription factors predicted to be regulators encoded by the Aⁿ subgenome are expressed primarily in the seed coat, whereas regulators encoded by the Cⁿ subgenome were expressed primarily in the embryo. Data suggest subgenome bias are characteristic features of the B. napus seed throughout development, and that such bias might not be universal across the embryo, endosperm and seed coat of the developing seed. Transcriptional networks spanning both the Aⁿ and Cⁿ genomes of the whole B. napus seed can identify valuable targets for seed development research and that -omics level approaches to studying gene regulation in B. napus can benefit from both broad and high-resolution analyses.

Constitutive expression of transcription factor SirZ blocks pathogenicity in Leptosphaeria maculans independently of sirodesmin production

Sirodesmin, the major secondary metabolite produced by the plant pathogenic fungus Leptosphaeria maculans in vitro, has been linked to disease on Brassica species since the 1970s, and yet its role has remained ambiguous. Re-examination of gene expression data revealed that all previously described genes and two newly identified genes within the sir gene cluster in the genome are down-regulated during the crucial early establishment stages of blackleg disease on Brassica napus. To test if this is a strategy employed by the fungus to avoid damage to and then detection by the host plant during the L. maculans asymptomatic biotrophic phase, sirodesmin was produced constitutively by overexpressing the sirZ gene encoding the transcription factor that coordinates the regulation of the other genes in the sir cluster. The sirZ over-expression strains had a major reduction in pathogenicity. Mutation of the over-expression construct restored pathogenicity. However, mutation of two genes, sirP and sirG, required for specific steps in the sirodesmin biosynthesis pathway, in the sirZ over-expression background resulted in strains that were unable to synthesize sirodesmin, yet were still non-pathogenic. Elucidating the basis for this pathogenicity defect or finding ways to overexpress sirZ during disease may provide new strategies for the control of blackleg disease.

Analysis of the Oxidative Burst and Its Relevant Signaling Pathways in Leptosphaeria maculans-Brassica napus Pathosystem

An oxidative burst is an early response of plants to various biotic/abiotic stresses. In plant-microbe interactions, the plant body can induce oxidative burst to activate various defense mechanisms to combat phytopathogens. A localized oxidative burst is also one of the typical behaviors during hypersensitive response (HR) caused by gene-for-gene interaction. In this study, the occurrence of oxidative burst and its signaling pathways was studied from different levels of disease severity (i.e., susceptible, intermediate, and resistant) in the B. napus-L. maculans pathosystem. Canola cotyledons with distinct levels of resistance exhibited differential regulation of the genes involved in reactive oxygen species (ROS) accumulation and responses. Histochemical assays were carried out to understand the patterns of H2O2 accumulation and cell death. Intermediate and resistant genotypes exhibited earlier accumulation of H2O2 and emergence of cell death around the inoculation origins. The observations also suggested that the cotyledons with stronger resistance were able to form a protective region of intensive oxidative bursts between the areas with and without hyphal intrusions to block further fungal advancement to the uninfected regions. The qPCR analysis suggested that different onset patterns of some marker genes in ROS accumulation/programmed cell death (PCD) such as RBOHD, MPK3 were associated with distinct levels of resistance from B. napus cultivars against L. maculans. The observations and datasets from this article indicated the distinct differences in ROS-related cellular behaviors and signaling between compatible and incompatible interactions.

Hormonal Responses to Susceptible, Intermediate, and Resistant Interactions in the Brassica napus-Leptosphaeria maculans Pathosystem

Hormone signaling plays a pivotal role in plant-microbe interactions. There are three major phytohormones in plant defense: salicylic acid (SA), jasmonic acid (JA), and ethylene (ET). The activation and trade-off of signaling between these three hormones likely determines the strength of plant defense in response to pathogens. Here, we describe the allocation of hormonal signaling in Brassica napus against the fungal pathogen Leptosphaeria maculans. Three B. napus genotypes (Westar, Surpass400, and 01-23-2-1) were inoculated with two L. maculans isolates (H75 8-1 and H77 7-2), subsequently exhibiting three levels of resistance: susceptible, intermediate, and resistant. Quantitative analyses suggest that the early activation of some SA-responsive genes, including WRKY70 and NPR1, contribute to an effective defense against L. maculans. The co-expression among factors responding to SA/ET/JA was also observed in the late stage of infection. The results of conjugated SA measurement also support that early SA activation plays a crucial role in durable resistance. Our results demonstrate the relationship between the onset patterns of certain hormone regulators and the effectiveness of the defense of B. napus against L. maculans.

The Effect of Temperature on the Hypersensitive Response (HR) in the Brassica napus-Leptosphaeria maculans Pathosystem

Temperature is considered one of the crucial environmental elements in plant pathological interactions, and previous studies have indicated that there is a relationship between temperature change and host-pathogen interactions. The objective of this research is to investigate the link between temperature and the incompatible interactions of the host and pathogen. In this study, two Leptosphaeria maculans isolates (HCRT75 8-1 and HCRT77 7-2) and two Brassica napus genotypes (Surpass400 and 01-23-2-1) were selected. The selected B. napus genotypes displayed intermediate and resistant phenotypes. The inoculated seedlings were tested under three temperature conditions: 16 °C/10 °C, 22 °C/16 °C and 28 °C/22 °C (day/night: 16 h/8 h). Lesion measurements demonstrated that the necrotic lesions from the 28 °C/22 °C treatment were enlarged compared with the other two temperature treatments (i.e., 16 °C/10 °C and 22 °C/16 °C). The results of expression analysis indicated that the three temperature treatments displayed distinct differences in two marker genes (PATHOGENESIS-RELATED (PR) 1 and 2) for plant defense and one temperature-sensitive gene BONZAI 1 (BON1). Additionally, seven dpi at 22 °C/16 °C appeared to be the optimal pre-condition for the induction of PR1 and 2. These findings suggest that B. napus responds to temperature changes when infected with L. maculans.

Zea mays RNA-seq estimated transcript abundances are strongly affected by read mapping bias

Background

Genetic variation for gene expression is a source of phenotypic variation for natural and agricultural species. The common approach to map and to quantify gene expression from genetically distinct individuals is to assign their RNA-seq reads to a single reference genome. However, RNA-seq reads from alleles dissimilar to this reference genome may fail to map correctly, causing transcript levels to be underestimated. Presently, the extent of this mapping problem is not clear, particularly in highly diverse species. We investigated if mapping bias occurred and if chromosomal features associated with mapping bias. Zea mays presents a model species to assess these questions, given it has genotypically distinct and well-studied genetic lines.

Results

In Zea mays, the inbred B73 genome is the standard reference genome and template for RNA-seq read assignments. In the absence of mapping bias, B73 and a second inbred line, Mo17, would each have an approximately equal number of regulatory alleles that increase gene expression. Remarkably, Mo17 had 2–4 times fewer such positively acting alleles than did B73 when RNA-seq reads were aligned to the B73 reference genome. Reciprocally, over one-half of the B73 alleles that increased gene expression were not detected when reads were aligned to the Mo17 genome template. Genes at dissimilar chromosomal ends were strongly affected by mapping bias, and genes at more similar pericentromeric regions were less affected. Biased transcript estimates were higher in untranslated regions and lower in splice junctions. Bias occurred across software and alignment parameters.

Conclusions

Mapping bias very strongly affects gene transcript abundance estimates in maize, and bias varies across chromosomal features. Individual genome or transcriptome templates are likely necessary for accurate transcript estimation across genetically variable individuals in maize and other species.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07577-3.

Genome-wide transcriptome reveals mechanisms underlying Rlm1-mediated blackleg resistance on canola

Genetic resistance to blackleg (Leptosphaeria maculans, Lm) of canola (Brassica napus, Bn) has been extensively studied, but the mechanisms underlying the host-pathogen interaction are still not well understood. Here, a comparative transcriptome analysis was performed on a resistant doubled haploid Bn line carrying the resistance gene Rlm1 following inoculation with a virulent (avrLm1) or avirulent (AvrLm1) Lm isolate on cotyledons. A total of 6999 and 3015 differentially expressed genes (DEGs) were identified, respectively, in inoculated local tissues with compatible (susceptible) and incompatible (resistant) interactions. Functional enrichment analysis found several biological processes, including protein targeting to membrane, ribosome and negative regulation of programmed cell death, were over-represented exclusively among up-regulated DEGs in the resistant reaction, whereas significant enrichment of salicylic acid (SA) and jasmonic acid (JA) pathways observed for down-regulated DEGs occurred only in the susceptible reaction. A heat-map analysis showed that both biosynthesis and signaling of SA and JA were induced more significantly in the resistant reaction, implying that a threshold level of SA and JA signaling is required for the activation of Rlm1-mediated resistance. Co-expression network analysis revealed close correlation of a gene module with the resistance, involving DEGs regulating pathogen-associated molecular pattern recognition, JA signaling and transcriptional reprogramming. Substantially fewer DEGs were identified in mock-inoculated (control) cotyledons, relative to those in inoculated local tissues, including those involved in SA pathways potentially contributing to systemic acquired resistance (SAR). Pre-inoculation of cotyledon with either an avirulent or virulent Lm isolate, however, failed to induce SAR on remote tissues of same plant despite elevated SA and PR1 protein. This study provides insights into the molecular mechanism of Rlm1-mediated resistance to blackleg.

Recent Findings Unravel Genes and Genetic Factors Underlying Leptosphaeria maculans Resistance in Brassica napus and Its Relatives

Among the Brassica oilseeds, canola (Brassica napus) is the most economically significant globally. However, its production can be limited by blackleg disease, caused by the fungal pathogen Lepstosphaeria maculans. The deployment of resistance genes has been implemented as one of the key strategies to manage the disease. Genetic resistance against blackleg comes in two forms: qualitative resistance, controlled by a single, major resistance gene (R gene), and quantitative resistance (QR), controlled by numerous, small effect loci. R-gene-mediated blackleg resistance has been extensively studied, wherein several genomic regions harbouring R genes against L. maculans have been identified and three of these genes were cloned. These studies advance our understanding of the mechanism of R gene and pathogen avirulence (Avr) gene interaction. Notably, these studies revealed a more complex interaction than originally thought. Advances in genomics help unravel these complexities, providing insights into the genes and genetic factors towards improving blackleg resistance. Here, we aim to discuss the existing R-gene-mediated resistance, make a summary of candidate R genes against the disease, and emphasise the role of players involved in the pathogenicity and resistance. The comprehensive result will allow breeders to improve resistance to L. maculans, thereby increasing yield.

Understanding Host-Pathogen Interactions in Brassica napus in the Omics Era

Brassica napus (canola/oilseed rape/rapeseed) is an economically important crop, mostly found in temperate and sub-tropical regions, that is cultivated widely for its edible oil. Major diseases of Brassica crops such as Blackleg, Clubroot, Sclerotinia Stem Rot, Downy Mildew, Alternaria Leaf Spot and White Rust have caused significant yield and economic losses in rapeseed-producing countries worldwide, exacerbated by global climate change, and, if not remedied effectively, will threaten global food security. To gain further insights into the host-pathogen interactions in relation to Brassica diseases, it is critical that we review current knowledge in this area and discuss how omics technologies can offer promising results and help to push boundaries in our understanding of the resistance mechanisms. Omics technologies, such as genomics, proteomics, transcriptomics and metabolomics approaches, allow us to understand the host and pathogen, as well as the interaction between the two species at a deeper level. With these integrated data in multi-omics and systems biology, we are able to breed high-quality disease-resistant Brassica crops in a more holistic, targeted and accurate way.

BrassicaEDB: A Gene Expression Database for Brassica Crops

The genus Brassica contains several economically important crops, including rapeseed (Brassica napus, 2n = 38, AACC), the second largest source of seed oil and protein meal worldwide. However, research in rapeseed is hampered because it is complicated and time-consuming for researchers to access different types of expression data. We therefore developed the Brassica Expression Database (BrassicaEDB) for the research community. In the current BrassicaEDB, we only focused on the transcriptome level in rapeseed. We conducted RNA sequencing (RNA-Seq) of 103 tissues from rapeseed cultivar ZhongShuang11 (ZS11) at seven developmental stages (seed germination, seedling, bolting, initial flowering, full-bloom, podding, and maturation). We determined the expression patterns of 101,040 genes via FPKM analysis and displayed the results using the eFP browser. We also analyzed transcriptome data for rapeseed from 70 BioProjects in the SRA database and obtained three types of expression level data (FPKM, TPM, and read counts). We used this information to develop the BrassicaEDB, including “eFP”, “Treatment”, “Coexpression”, and “SRA Project” modules based on gene expression profiles and “Gene Feature”, “qPCR Primer”, and “BLAST” modules based on gene sequences. The BrassicaEDB provides comprehensive gene expression profile information and a user-friendly visualization interface for rapeseed researchers. Using this database, researchers can quickly retrieve the expression level data for target genes in different tissues and in response to different treatments to elucidate gene functions and explore the biology of rapeseed at the transcriptome level.

Integrated transcriptomic and metabolomic analyses reveal the effects of callose deposition and multihormone signal transduction pathways on the tea plant-Colletotrichum camelliae interaction

Colletotrichum infects diverse hosts, including tea plants, and can lead to crop failure. Numerous studies have reported that biological processes are involved in the resistance of tea plants to Colletotrichum spp. However, the molecular and biochemical responses in the host during this interaction are unclear. Cuttings of the tea cultivar Longjing 43 (LJ43) were inoculated with a conidial suspension of Colletotrichum camelliae, and water-sprayed cuttings were used as controls. In total, 10,592 differentially expressed genes (DEGs) were identified from the transcriptomic data of the tea plants and were significantly enriched in callose deposition and the biosynthesis of various phytohormones. Subsequently, 3,555 mass spectra peaks were obtained by LC-MS detection in the negative ion mode, and 27, 18 and 81 differentially expressed metabolites (DEMs) were identified in the tea leaves at 12 hpi, 24 hpi and 72 hpi, respectively. The metabolomic analysis also revealed that the levels of the precursors and intermediate products of jasmonic acid (JA) and indole-3-acetate (IAA) biosynthesis were significantly increased during the interaction, especially when the symptoms became apparent. In conclusion, we suggest that callose deposition and various phytohormone signaling systems play important roles in the tea plant-C. camelliae interaction.

Exploring Mechanisms of Quantitative Resistance to Leptosphaeria maculans (Blackleg) in the Cotyledons of Canola (Brassica napus) Based on Transcriptomic and Microscopic Analyses

Using resistant cultivars is a common approach to managing blackleg of canola/rapeseed caused by Leptosphaeria maculans (Lm). Quantitative resistance (QR), as opposed to major-gene resistance, is of interest because it is generally more durable, due to its multi-genetic basis. However, the mechanisms and genes underlying QR are mostly unknown. In this study, potential QR modes of action in “74-44 BL” was explored. This Canadian canola cultivar showed moderate but consistent race-nonspecific resistance at the cotyledon and adult-plant stages. A susceptible cultivar, “Westar”, was used as a control. After inoculation, the lesions developed more slowly on the cotyledons of 74-44 BL than those of Westar. We used RNA sequencing (-RNA-seq) to identify genes and their functions, putatively related to this resistance, and found that genes involved in programmed cell death (PCD), reactive oxygen species (ROS), signal transduction or intracellular endomembrane transport were most differentially expressed. ROS production was assessed in relation to Lm hyphal growth and lesion size; it occurred beyond the tissue colonized by Lm in 74-44 BL and appeared to trigger rapid cell death, limiting cotyledon colonization by Lm. In contrast, Lm grew more rapidly in Westar, often catching up with the ring of ROS and surpassing lesion boundaries. It appears that QR in 74-44 BL cotyledons is associated with limited colonization by Lm possibly mediated via ROS. The RNA-seq data also showed a link between ROS, signal transduction, and endomembrane vesicle trafficking, as well as PCD in the resistance. These results provide a starting point for a better understanding of the mechanisms behind QR against Lm in canola.

Interactive Gene Expression Patterns of Susceptible and Resistant Lens ervoides Recombinant Inbred Lines and the Necrotroph Ascochyta lentis

Ascochyta lentis is a foliar pathogen of Lens species and is of worldwide importance in cultivated lentil production. High levels of resistance were identified in the wild species Lens ervoides. This resistance was explored through histopathology, qPCR estimation of fungal biomass and transcriptome sequencing in a susceptible and a resistant recombinant inbred line (RIL) of L. ervoides infected with an aggressive isolate of A. lentis. Necrotrophic growth was delayed in the resistant RIL compared to accelerated necrotrophy of A. lentis in the susceptible RIL. Analysis of the fungal secretome indicated that the early activation of cell wall-degrading enzymes contributed to increased virulence of A. lentis. On the host side, gene co-expression analysis revealed that the invasion by A. lentis caused mRNA, DNA and protein decay in infected plants regardless of the level of resistance in the host. The resistant RIL exhibited a stronger gene co-expression in lipid localization and sulfur processes, and cellular responses to nutrients and stimuli than the susceptible RIL. In addition, differential gene analysis revealed that the repression of both, gibberellin signaling and cell death associated with the hypersensitive response (HR), were associated with enhanced A. lentis resistance.

Genome-wide mining and comparative analysis of fatty acid elongase gene family in Brassica napus and its progenitors

Very long chain fatty acids (VLCFAs) that are structural components of cell membrane lipid, cuticular waxes and seed oil, play crucial roles in plant growth, development and stress response. Fatty acid elongases (FAEs) comprising KCS and ELO, are key enzymes for VLCFA biosynthesis in plants. Although reference genomes of Brassica napus and its parental speices both have been sequenced, whole-genome analysis of FAE gene family in these Brassica speices is not reported. Here, 58, 33 and 30 KCS genes were identified in B. napus, B. rapa and B. oleracea genomes, respectively, whereas 14, 6 and 8 members were obtained for ELO genes. These KCS genes were unevenly located in 37 chromosomes and 3 scaffolds of 3 Brassica species, while these ELO genes were mapped to 19 chromosomes. The KCS and ELO proteins were divided into 8 and 4 subclasses, respectively. Gene structure and protein motifs remained highly conserved in each KCS or ELO subclass. Most promoters of KCS and ELO genes harbored various plant growth-, phytohormone-, and stress response-related cis-acting elements. 20 SSR loci existed in the KCS and ELO genes/promoters. The whole-genome duplication and segmental duplication mainly contributed to expansion of KCS and ELO genes in these genomes. Transcriptome analysis showed that KCS and ELO genes in 3 Brassica species were expressed in various tissues/organs with different levels, whereas 1 BnELO gene and 6 BnKCS genes might be pathogen-responsive genes. The qRT-PCR assay showed that BnKCS22 and BnELO04 responded to various phytohormone treatments and abiotic stresses. This work lays the foundation for further function identification of KCS and ELO genes in B. napus and its progenitors.

Pigment variation and transcriptional response of the pigment synthesis pathway in the S2309 triple-color ornamental kale (Brassica oleracea L. var. acephala) line

Ornamental kale is popular because of its colorful leaves and few studies have investigated the mechanism of color changes. In this study, an ornamental kale line (S2309) with three leaf colors was developed. Analysis of the anthocyanin, chlorophyll, and carotenoid contents and RNA-seq were performed on the three leaf color types. There was less chlorophyll in the white leaves and purple leaves than in the green leaves, and the anthocyanin content was greatest in the purple leaves. All the downregulated DEGs related to chlorophyll metabolism were detected only in the S2309_G vs. S2309_W comparison, which indicated that the decrease in chlorophyll content was caused mainly by the inhibition of chlorophyll biosynthesis during the leaf color change from green to white. Moreover, the expression of 19 DEGs involved in the anthocyanin biosynthesis pathway was upregulated. These results provide new insight into the mechanisms underlying the three-color formation.

Comparative transcriptome provides molecular insight into defense-associated mechanisms against spider mite in resistant and susceptible common bean cultivars

Common bean (Phaseolus vulgaris L.) is a major source of proteins and one of the most important edible foods for more than three hundred million people in the world. The common bean plants are frequently attacked by spider mite (Tetranychus urticae Koch), leading to a significant decrease in plant growth and economic performance. The use of resistant cultivars and the identification of the genes involved in plant-mite resistance are practical solutions to this problem. Hence, a comprehensive study of the molecular interactions between resistant and susceptible common bean cultivars and spider mite can shed light into the understanding of mechanisms and biological pathways of resistance. In this study, one resistant (Naz) and one susceptible (Akhtar) cultivars were selected for a transcriptome comparison at different time points (0, 1 and 5 days) after spider mite feeding. The comparison of cultivars in different time points revealed several key genes, which showed a change increase in transcript abundance via spider mite infestation. These included genes involved in flavonoid biosynthesis process; a conserved MYB-bHLH-WD40 (MBW) regulatory complex; transcription factors (TFs) TT2, TT8, TCP, Cys2/His2-type and C2H2-type zinc finger proteins; the ethylene response factors (ERFs) ERF1 and ERF9; genes related to metabolism of auxin and jasmonic acid (JA); pathogenesis-related (PR) proteins and heat shock proteins.

Mycovirus-Induced Hypervirulence of Leptosphaeria biglobosa Enhances Systemic Acquired Resistance to Leptosphaeria maculans in Brassica napus

Phoma stem canker (blackleg) is one of the most important diseases of winter oilseed rape (Brassica napus) worldwide and is caused by a complex that comprises at least two species: Leptosphaeria maculans and L. biglobosa. Screening a panel of field Leptosphaeria isolates from B. napus for the presence of mycoviruses revealed the presence of a novel double-stranded RNA quadrivirus in L. biglobosa and no viruses in L. maculans. Following elimination of the mycovirus, virus-infected and virus-free isogenic lines of L. biglobosa were created. A direct comparison of the growth and virulence of these isogenic lines illustrated that virus infection caused hypervirulence and resulted in induced systemic resistance toward L. maculans in B. napus following lower leaf preinoculation with the virus-infected isolate. Analysis of the plant transcriptome suggests that the presence of the virus leads to subtle alterations in metabolism and plant defenses. For instance, transcripts involved in carbohydrate and amino acid metabolism are enriched in plants treated with the virus-infected isolate, while pathogenesis-related proteins, chitinases and WRKY transcription factors are differentially expressed. These results illustrate the potential for deliberate inoculation of plants with hypervirulent L. biglobosa to decrease the severity of Phoma stem canker later in the growing season.[Formula: see text] Copyright © 2020 The Author(s). This is an open access article distributed under the CC BY 4.0 International license.

Directing Trophic Divergence in Plant-Pathogen Interactions: Antagonistic Phytohormones With NO Doubt?

A fundamental process culminating in the mechanisms of plant-pathogen interactions is the regulation of trophic divergence into biotrophic, hemibiotrophic, and necrotrophic interactions. Plant hormones, of almost all types, play significant roles in this regulatory apparatus. In plant-pathogen interactions, two classical mechanisms underlying hormone-dependent trophic divergence are long recognized. While salicylic acid dominates in the execution of host defense response against biotrophic and early-stage hemibiotrophic pathogens, jasmonic acid, and ethylene are key players facilitating host defense response against necrotrophic and later-stage hemibiotrophic pathogens. Evidence increasingly suggests that trophic divergence appears to be modulated by more complex signaling networks. Acting antagonistically or agonistically, other hormones such as auxins, cytokinins, abscisic acid, gibberellins, brassinosteroids, and strigolactones, as well as nitric oxide, are emerging candidates in the regulation of trophic divergence. In this review, the latest advances in the dynamic regulation of trophic divergence are summarized, emphasizing common and contrasting hormonal and nitric oxide signaling strategies deployed in plant-pathogen interactions.

Transcriptome Analysis of Rlm2-Mediated Host Immunity in the Brassica napus-Leptosphaeria maculans Pathosystem

Our study investigated disease resistance in the Brassica napus-Leptosphaeria maculans pathosystem using a combination of laser microdissection, dual RNA sequencing, and physiological validations of large-scale gene sets. The use of laser microdissection improved pathogen detection and identified putative L. maculans effectors and lytic enzymes operative during host colonization. Within 24 h of inoculation, we detected large shifts in gene activity in resistant cotyledons associated with jasmonic acid and calcium signaling pathways that accelerated the plant defense response. Sequencing data were validated through the direct quantification of endogenous jasmonic acid levels. Additionally, resistance against L. maculans was abolished when the calcium chelator EGTA was applied to the inoculation site, providing physiological evidence of the role of calcium in B. napus immunity against L. maculans. We integrated gene expression data with all available information on cis-regulatory elements and transcription factor binding affinities to better understand the gene regulatory networks underpinning plant resistance to hemibiotrophic pathogens. These in silico analyses point to early cellular reprogramming during host immunity that are coordinated by CAMTA, BZIP, and bHLH transcription factors. Together, we provide compelling genetic and physiological evidence into the programming of plant resistance against fungal pathogens.

Transcriptional Insight Into Brassica napus Resistance Genes LepR3 and Rlm2-Mediated Defense Response Against the Leptosphaeria maculans Infection

The phytopathogenic fungus Leptosphaeria maculans causes the blackleg disease on Brassica napus, resulting in severe loss of rapeseed production. Breeding of resistant cultivars containing race-specific resistance genes is provably effective to combat this disease. While two allelic resistance genes LepR3 and Rlm2 recognizing L. maculans avirulence genes AvrLm1 and AvrLm2 at plant apoplastic space have been cloned in B. napus, the downstream gene expression network underlying the resistance remains elusive. In this study, transgenic lines expressing LepR3 and Rlm2 were created in the susceptible "Westar" cultivar and inoculated with L. maculans isolates containing different sets of AvrLm1 and AvrLm2 for comparative transcriptomic analysis. Through grouping the RNA-seq data based on different levels of defense response, we find LepR3 and Rlm2 orchestrate a hierarchically regulated gene expression network, consisting of induced ABA acting independently of the disease reaction, activation of signal transduction pathways with gradually increasing intensity from compatible to incompatible interaction, and specifically induced enzymatic and chemical actions contributing to hypersensitive response with recognition of AvrLm1 and AvrLm2. This study provides an unconventional investigation into LepR3 and Rlm2-mediated plant defense machinery and adds novel insight into the interaction between surface-localized receptor-like proteins (RLPs) and apoplastic fungal pathogens.

Genome wide identification and comparative analysis of glutathione transferases (GST) family genes in Brassica napus

Glutathione transferases (GSTs) are multifunctional enzymes that play important roles in plant development and responses to biotic and abiotic stress. However, a systematic analysis of GST family members in Brassica napus has not yet been reported. In this study, we identified 179 full-length GST genes in B. napus, 44.2% of which are clustered on various chromosomes. In addition, we identified 141 duplicated GST gene pairs in B. napus. Molecular evolutionary analysis showed that speciation and whole-genome triplication played important roles in the divergence of the B. napus GST duplicated genes. Transcriptome analysis of 21 tissues at different developmental stages showed that 47.6% of duplicated GST gene pairs have divergent expression patterns, perhaps due to structural divergence. We constructed a GST gene coexpression network with genes encoding various transcription factors (NAC, MYB, WRKY and bZIP) and identified six modules, including genes expressed during late seed development (after 40 days; BnGSTU19, BnGSTU20 and BnGSTZ1) and in the seed coat (BnGSTF6 and BnGSTF12), stamen and anther (BnGSTF8), root and stem (BnGSTU21), leaves and funiculus, as well as during the late stage of pericarp development (after 40 days; BnGSTU12 and BnGSTF2) and in the radicle during seed germination (BnGSTF14, BnGSTU1, BnGSTU28, and BnGSTZ1). These findings lay the foundation for elucidating the roles of GSTs in B. napus.

Spatial identification of transcripts and biological processes in laser micro-dissected sub-regions of waterlogged corn roots with altered expression of phytoglobin

Over-expression of the corn phytoglobin ZmPgb1.2 increases tolerance to waterlogging, while suppression of ZmPgb1.2 compromises plant growth. To unravel compartment-specific transcriptional changes evoked by ZmPgb1.2 during hypoxia, laser micro-dissected sub-regions from waterlogged roots of WT and ZmPgb1.2 overexpressing [ZmPgb1.2(S)] plants were probed for global transcriptional analysis using next generation RNA sequencing. These sub-regions included compartments within the meristematic, elongation, and maturation zone. Of the 149 genes differentially expressed by the up-regulation of ZmPgb1.2, 78 occurred within the meristematic region and included genes involved in jasmonic acid synthesis and response, ascorbic acid metabolism, and ethylene signalling. The ZmPgb1.2 regulation of these genes, discussed in the context of known functions of Pgbs, was further validated by monitoring their expression in meristematic cells of waterlogged roots suppressing ZmPgb1.2. Of the 27 genes differentially expressed by the over-expression of ZmPgb1.2 in the elongation zone, pyruvate kinase and alcohol dehydrogenase showed an expression pattern correlated to the level of ZmPgb1.2 in the tissue. The transcriptional induction of these two enzymes in hypoxic domains of the elongation zone over-expressing ZmPgb1.2 suggests the activation of the fermentation pathway which might be required to sustain metabolic flux and production of ATP in support of cell elongation.

Dissecting R gene and host genetic background effect on the Brassica napus defense response to Leptosphaeria maculans

While our understanding of the genetics underlying the Brassica-Leptosphaeria pathosystem has advanced greatly in the last decade, differences in molecular responses due to interaction between resistance genes and host genetic background has not been studied. We applied RNAseq technology to monitor the transcriptome profiles of Brassica napus (Bn) lines carrying one of four blackleg R genes (Rlm2, Rlm3, LepR1 & LepR2) in Topas or Westar background, during the early stages of infection by a Leptosphaeria maculans (Lm) isolate carrying the corresponding Avr genes. We observed upregulation of host genes involved in hormone signalling, cell wall thickening, response to chitin and glucosinolate production in all R gene lines at 3 day after inoculation (dai) albeit having higher level of expression in LepR1 and Rlm2 than in Rlm3 and LepR2 lines. Bn-SOBIR1 (Suppressor Of BIR1-1), a receptor like kinase (RLK) that forms complex receptor like proteins (RLPs) was highly expressed in LepR1 and Rlm2 at 3 dai. In contrast Bn-SOBIR1 induction was low in Rlm3 line, which could indicate that Rlm3 may function independent of SOBIR1. Expression of Salicylic acid (SA) related defense was enhanced in LepR1 and Rlm2 at 3 dai. In contrast to SA, expression of Bn genes with homology to PDF1.2, a jasmonic acid (JA) pathway marker, were increased in all Rlm and LepR lines at 6 and 9 dai. Effect of host genetic background on induction of defense, was determined by comparison of LepR1 and LepR2 in Topas vs Westar genotype (i.e. T-LepR1 vs W-LepR1 and T-LepR2 vs W-LepR2). In both cases (regardless of R gene) overall number of defense related genes at the earliest time point (3 dai) was higher in Tops compared to Westar. SA and JA markers genes such as PR1 and PDF1.2 were more induced in Topas compared to Westar introgression lines at this time point. Even in the absence of any R gene, effect of Topas genotype in enhanced defense, was also evident by the induction of PDF1.2 that started at a low level at 3 dai and peaked at 6 and 9 dai, while no induction in Westar genotype was observed at any of these time points. Overall, variation in time and intensity of expression of genes related to defense, was clearly dependent on both R gene and the host genotype.

Evolutionary and ecological functional genomics, from lab to the wild

Plant phenotypes are the result of both genetic and environmental forces that act to modulate trait expression. Over the last few years, numerous approaches in functional genomics and systems biology have led to a greater understanding of plant phenotypic variation and plant responses to the environment. These approaches, and the questions that they can address, have been loosely termed evolutionary and ecological functional genomics (EEFG), and have been providing key insights on how plants adapt and evolve. In particular, by bringing these studies from the laboratory to the field, EEFG studies allow us to gain greater knowledge of how plants function in their natural contexts.

Tissue-specific localization of polyketide synthase and other associated genes in the lichen, Cladonia rangiferina, using laser microdissection

The biosynthesis of two polyketides, atranorin and fumarprotocetraric acid, produced from a lichen-forming fungus, Cladonia rangiferina (L.) F. H. Wigg. was correlated with the expression of eight fungal genes (CrPKS1, CrPKS3, CrPKS16, Catalase (CAT), Sugar Transporter (MFsug), Dioxygenase (YQE1), C₂H₂ Transcription factor (C₂H₂), Transcription Factor PacC (PacC), which are thought to be involved in polyketide biosynthesis, and one algal gene, NAD-dependent deacetylase sirtuin 2 (AsNAD)), using laser microdissection (LMD). The differential gene expression levels within the thallus tissue layers demonstrate that the most active region for potential polyketide biosynthesis within the lichen is the outer apical region proximal to the photobiont but some expression also occurs in reproductive tissue. This is the first study using laser microdissection to explore gene expression of these nine genes and their location of expression; it provides a proof-of-concept for future experiments exploring tissue-specific gene expression within lichens; and it highlights the utility of LMD for use in lichen systems.

Characterization of Callose Deposition and Analysis of the Callose Synthase Gene Family of Brassica napus in Response to Leptosphaeria maculans

Callose plays a critical role in different biological processes including development as well as in the response to multiple biotic and abiotic stresses. In this study, we characterized the callose deposition in cotyledons of different Brassica napus varieties post-inoculated with different Leptosphaeria maculans isolates. Further, members of the callose synthase gene were identified from the whole genome of B. napus using the 12 Arabidopsis thaniana callose synthase protein sequences, and were then classified into three groups based on their phylogenetic relationships. Chromosomal location and duplication patterns indicated uneven distribution and segmental duplication patterns of BnCalS genes in the B. napus genome. Subsequently, gene structures, conserved domains analysis, and protein properties were analyzed for BnCalS genes. In addition, 12 B. napus orthologs of the AtCalS were selected for investigating the tissue expression pattern, indicating diverse expression patterns for these BnCalS genes. Responses of the selected 12 orthologs and all the BnCalS genes were characterized in the different types (AvrLm1-Rlm1, AvrLm4-Rlm4, AvrLepR1-LepR1) of B. napus–L. maculans interactions and B. napus-Leptosphaeria biglobosa interactions, implying their potential roles in response to Leptosphaeria infection.

Stable Quantitative Resistance Loci to Blackleg Disease in Canola (Brassica napus L.) Over Continents

The hemibiotrophic fungus, Leptosphaeria maculans is the most devastating pathogen, causing blackleg disease in canola (Brassica napus L). To study the genomic regions involved in quantitative resistance (QR), 259-276 DH lines from Darmor-bzh/Yudal (DYDH) population were assessed for resistance to blackleg under shade house and field conditions across 3 years. In different experiments, the broad sense heritability varied from 43 to 95%. A total of 27 significant quantitative trait loci (QTL) for QR were detected on 12 chromosomes and explained between 2.14 and 10.13% of the genotypic variance. Of the significant QTL, at least seven were repeatedly detected across different experiments on chromosomes A02, A07, A09, A10, C01, and C09. Resistance alleles were mainly contributed by 'Darmor-bzh' but 'Yudal' also contributed few of them. Our results suggest that plant maturity and plant height may have a pleiotropic effect on QR in our conditions. We confirmed that Rlm9 which is present in 'Darmor-bzh' is not effective to confer resistance in our Australian field conditions. Comparative mapping showed that several R genes coding for nucleotide-binding leucine-rich repeat (LRR) receptors map in close proximity (within 200 Kb) of the significant trait-marker associations on the reference 'Darmor-bzh' genome assembly. More importantly, eight significant QTL regions were detected across diverse growing environments: Australia, France, and United Kingdom. These stable QTL identified herein can be utilized for enhancing QR in elite canola germplasm via marker- assisted or genomic selection strategies.

The initial deficiency of protein processing and flavonoids biosynthesis were the main mechanisms for the male sterility induced by SX-1 in Brassica napus

Background

Rapeseed (Brassica napus) is an important oil seed crop in the Brassicaceae family. Chemical induced male sterility (CIMS) is one of the widely used method to produce the hybrids in B. napus. Identification of the key genes and pathways that involved in CIMS were important to understand the underlying molecular mechanism. In the present report, a multi-omics integrative analysis, including of the proteomic, transcriptomic and miRNAs, combined with morphological and physiological analysis were conducted.

Results

Earlier degeneration of the tapetosomes and elaioplasts, aberrantly stacking in tapetal cells and incompletely deposition in tryphine of pollen wall were observed in chemical hybridization agent (CHA) of SX-1 treated B. napus through SEM and TEM analysis. It was revealed that the deficiencies in protein processing in endoplasmic reticulum (ER) and flavonoids biosynthesis were occurred at early stage in the SX-1 treated materials. Subsequently, plant hormone signal transduction, biosynthesis of amino acids, fatty acids and steroid in anther at later stages were identified down-regulated after SX-1 treatment. 144 transcript factors (TFs) were also indentified to down-regulated at early stage, which suggested the early regulation in anther and pollen wall development were disordered in CHA treated B. napus. In addition, 7 important miRNAs were identified and 2 of the predicted target genes of miRNAs were Rf-like genes.

Conclusions

Taken together, an interaction network of candidate genes and the putative metabolism pathways were constructed based on the multi-omics integrative analysis, it provided a new insight into the male sterility induced by CHA of SX-1 in B. napus.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5203-y) contains supplementary material, which is available to authorized users.

Multi-year linkage and association mapping confirm the high number of genomic regions involved in oilseed rape quantitative resistance to blackleg

A repertoire of the genomic regions involved in quantitative resistance to Leptosphaeria maculans in winter oilseed rape was established from combined linkage-based QTL and genome-wide association (GWA) mapping. Linkage-based mapping of quantitative trait loci (QTL) and genome-wide association studies are complementary approaches for deciphering the genomic architecture of complex agronomical traits. In oilseed rape, quantitative resistance to blackleg disease, caused by L. maculans, is highly polygenic and is greatly influenced by the environment. In this study, we took advantage of multi-year data available on three segregating populations derived from the resistant cv Darmor and multi-year data available on oilseed rape panels to obtain a wide overview of the genomic regions involved in quantitative resistance to this pathogen in oilseed rape. Sixteen QTL regions were common to at least two biparental populations, of which nine were the same as previously detected regions in a multi-parental design derived from different resistant parents. Eight regions were significantly associated with quantitative resistance, of which five on A06, A08, A09, C01 and C04 were located within QTL support intervals. Homoeologous Brassica napus genes were found in eight homoeologous QTL regions, which corresponded to 657 pairs of homoeologous genes. Potential candidate genes underlying this quantitative resistance were identified. Genomic predictions and breeding are also discussed, taking into account the highly polygenic nature of this resistance.

Genomic evidence for genes encoding leucine-rich repeat receptors linked to resistance against the eukaryotic extra- and intracellular Brassica napus pathogens Leptosphaeria maculans and Plasmodiophora brassicae

Genes coding for nucleotide-binding leucine-rich repeat (LRR) receptors (NLRs) control resistance against intracellular (cell-penetrating) pathogens. However, evidence for a role of genes coding for proteins with LRR domains in resistance against extracellular (apoplastic) fungal pathogens is limited. Here, the distribution of genes coding for proteins with eLRR domains but lacking kinase domains was determined for the Brassica napus genome. Predictions of signal peptide and transmembrane regions divided these genes into 184 coding for receptor-like proteins (RLPs) and 121 coding for secreted proteins (SPs). Together with previously annotated NLRs, a total of 720 LRR genes were found. Leptosphaeria maculans-induced expression during a compatible interaction with cultivar Topas differed between RLP, SP and NLR gene families; NLR genes were induced relatively late, during the necrotrophic phase of pathogen colonization. Seven RLP, one SP and two NLR genes were found in Rlm1 and Rlm3/Rlm4/Rlm7/Rlm9 loci for resistance against L. maculans on chromosome A07 of B. napus. One NLR gene at the Rlm9 locus was positively selected, as was the RLP gene on chromosome A10 with LepR3 and Rlm2 alleles conferring resistance against L. maculans races with corresponding effectors AvrLm1 and AvrLm2, respectively. Known loci for resistance against L. maculans (extracellular hemi-biotrophic fungus), Sclerotinia sclerotiorum (necrotrophic fungus) and Plasmodiophora brassicae (intracellular, obligate biotrophic protist) were examined for presence of RLPs, SPs and NLRs in these regions. Whereas loci for resistance against P. brassicae were enriched for NLRs, no such signature was observed for the other pathogens. These findings demonstrate involvement of (i) NLR genes in resistance against the intracellular pathogen P. brassicae and a putative NLR gene in Rlm9-mediated resistance against the extracellular pathogen L. maculans.

Leptosphaeria maculans Effector Protein AvrLm1 Modulates Plant Immunity by Enhancing MAP Kinase 9 Phosphorylation

Leptosphaeria maculans, the causal agent of blackleg disease in canola (Brassica napus), secretes an array of effectors into the host to overcome host defense. Here we present evidence that the L. maculans effector protein AvrLm1 functions as a virulence factor by interacting with the B. napus mitogen-activated protein (MAP) kinase 9 (BnMPK9), resulting in increased accumulation and enhanced phosphorylation of the host protein. Transient expression of BnMPK9 in Nicotiana benthamiana induces cell death, and this phenotype is enhanced in the presence of AvrLm1, suggesting that induction of cell death due to enhanced accumulation and phosphorylation of BnMPK9 by AvrLm1 supports the initiation of necrotrophic phase of L. maculans infection. Stable expression of BnMPK9 in B. napus perturbs hormone signaling, notably salicylic acid response genes, to facilitate L. maculans infection. Our findings provide evidence that a MAP kinase is directly targeted by a fungal effector to modulate plant immunity.

•

Leptosphaeria maculans effector AvrLm1 interacts with the Brassica napus MPK9 (BnMPK9)

•

AvrLm1 increases the accumulation and enhances the phosphorylation of BnMPK9

•

AvrLm1 enhances BnMPK9-dependent cell death in Nicotiana benthamiana

•

Stable expression of BnMPK9 in B. napus facilitates L. maculans infection

Comparative genomic and transcriptomic analyses of Family-1 UDP glycosyltransferase in three Brassica species and Arabidopsis indicates stress-responsive regulation

In plants, UGTs (UDP-glycosyltransferases) glycosylate various phytohormones and metabolites in response to biotic and abiotic stresses. Little is known about stress-responsive glycosyltransferases in plants. Therefore, it is important to understand the genomic and transcriptomic portfolio of plants with regard to biotic and abiotic stresses. Here, we identified 140, 154, and 251 putative UGTs in Brassica rapa, Brassica oleracea, and Brassica napus, respectively, and clustered them into 14 major phylogenetic groups (A–N). Fourteen major KEGG pathways and 24 biological processes were associated with the UGTs, highlighting them as unique modulators against environmental stimuli. Putative UGTs from B. rapa and B. oleracea showed a negative selection pressure and biased gene fractionation pattern during their evolution. Polyploidization increased the intron proportion and number of UGT-containing introns among Brassica. The putative UGTs were preferentially expressed in developing tissues and at the senescence stage. Differential expression of up- and down-regulated UGTs in response to phytohormone treatments, pathogen responsiveness and abiotic stresses, inferred from microarray and RNA-Seq data in Arabidopsis and Brassica broaden the glycosylation impact at the molecular level. This study identifies unique candidate UGTs for the manipulation of biotic and abiotic stress pathways in Brassica and Arabidopsis.

Transcriptional analysis and histochemistry reveal that hypersensitive cell death and H2O2 have crucial roles in the resistance of tea plant (Camellia sinensis (L.) O. Kuntze) to anthracnose

Anthracnose causes severe losses of tea production in China. Although genes and biological processes involved in anthracnose resistance have been reported in other plants, the molecular response to anthracnose in tea plant is unknown. We used the susceptible tea cultivar Longjing 43 and the resistant cultivar Zhongcha 108 as materials and compared transcriptome changes in the leaves of both cultivars following Colletotrichum fructicola inoculation. In all, 9015 and 8624 genes were differentially expressed between the resistant and susceptible cultivars and their controls (0 h), respectively. In both cultivars, the differentially expressed genes (DEGs) were enriched in 215 pathways, including responses to sugar metabolism, phytohormones, reactive oxygen species (ROS), biotic stimuli and signalling, transmembrane transporter activity, protease activity and signalling receptor activity, but DEG expression levels were higher in Zhongcha 108 than in Longjing 43. Moreover, functional enrichment analysis of the DEGs showed that hydrogen peroxide (H₂O₂) metabolism, cell death, secondary metabolism, and carbohydrate metabolism are involved in the defence of Zhongcha 108, and 88 key genes were identified. Protein-protein interaction (PPI) network demonstrated that putative mitogen-activated protein kinase (MAPK) cascades are activated by resistance (R) genes and mediate downstream defence responses. Histochemical analysis subsequently validated the strong hypersensitive response (HR) and H₂O₂ accumulation that occurred around the hyphal infection sites in Zhongcha 108. Overall, our results indicate that the HR and H₂O₂ are critical mechanisms in tea plant defence against anthracnose and may be activated by R genes via MAPK cascades.

RNA sequencing of Brassica napus reveals cellular redox control of Sclerotinia infection

Brassica napus is one of the world's most valuable oilseeds and is under constant pressure by the necrotrophic fungal pathogen, Sclerotinia sclerotiorum, the causal agent of white stem rot. Despite our growing understanding of host pathogen interactions at the molecular level, we have yet to fully understand the biological processes and underlying gene regulatory networks responsible for determining disease outcomes. Using global RNA sequencing, we profiled gene activity at the first point of infection on the leaf surface 24 hours after pathogen exposure in susceptible (B. napus cv. Westar) and tolerant (B. napus cv. Zhongyou 821) plants. We identified a family of ethylene response factors that may contribute to host tolerance to S. sclerotiorum by activating genes associated with fungal recognition, subcellular organization, and redox homeostasis. Physiological investigation of redox homeostasis was further studied by quantifying cellular levels of the glutathione and ascorbate redox pathway and the cycling enzymes associated with host tolerance to S. sclerotiorum. Functional characterization of an Arabidopsis redox mutant challenged with the fungus provides compelling evidence into the role of the ascorbate-glutathione redox hub in the maintenance and enhancement of plant tolerance against fungal pathogens.

SeqEnrich: A tool to predict transcription factor networks from co-expressed Arabidopsis and Brassica napus gene sets

Transcription factors and their associated DNA binding sites are key regulatory elements of cellular differentiation, development, and environmental response. New tools that predict transcriptional regulation of biological processes are valuable to researchers studying both model and emerging-model plant systems. SeqEnrich predicts transcription factor networks from co-expressed Arabidopsis or Brassica napus gene sets. The networks produced by SeqEnrich are supported by existing literature and predicted transcription factor–DNA interactions that can be functionally validated at the laboratory bench. The program functions with gene sets of varying sizes and derived from diverse tissues and environmental treatments. SeqEnrich presents as a powerful predictive framework for the analysis of Arabidopsis and Brassica napus co-expression data, and is designed so that researchers at all levels can easily access and interpret predicted transcriptional circuits. The program outperformed its ancestral program ChipEnrich, and produced detailed transcription factor networks from Arabidopsis and Brassica napus gene expression data. The SeqEnrich program is ideal for generating new hypotheses and distilling biological information from large-scale expression data.