Cereal cyst nematode Heterodera avenae, an important soil-borne pathogen in wheat, causes numerous annual yield losses worldwide, and use of resistant cultivars is the best strategy for control. However, target genes are not readily available for breeding resistant cultivars. Therefore, comparative transcriptomic analyses were performed to identify more applicable resistance genes for cultivar breeding. The developing nematodes within roots were stained with acid fuchsin solution. Transcriptome assemblies and redundancy filteration were obtained by Trinity, TGI Clustering Tool and BLASTN, respectively. Gene Ontology annotation was yielded by Blast2GO program, and metabolic pathways of transcripts were analyzed by Path_finder. The ROS levels were determined by luminol-chemiluminescence assay. The transcriptional gene expression profiles were obtained by quantitative RT-PCR. The RNA-sequencing was performed using an incompatible wheat cultivar VP1620 and a compatible control cultivar WEN19 infected with H. avenae at 24 h, 3 d and 8 d. Infection assays showed that VP1620 failed to block penetration of H. avenae but disturbed the transition of developmental stages, leading to a significant reduction in cyst formation. Two types of expression profiles were established to predict candidate resistance genes after developing a novel strategy to generate clean RNA-seq data by removing the transcripts of H. avenae within the raw data before assembly. Using the uncoordinated expression profiles with transcript abundance as a standard, 424 candidate resistance genes were identified, including 302 overlapping genes and 122 VP1620-specific genes. Genes with similar expression patterns were further classified according to the scales of changed transcript abundances, and 182 genes were rescued as supplementary candidate resistance genes. Functional characterizations revealed that diverse defense-related pathways were responsible for wheat resistance against H. avenae. Moreover, phospholipase was involved in many defense-related pathways and localized in the connection position. Furthermore, strong bursts of reactive oxygen species (ROS) within VP1620 roots infected with H. avenae were induced at 24 h and 3 d, and eight ROS-producing genes were significantly upregulated, including three class III peroxidase and five lipoxygenase genes. Large-scale identification of wheat resistance genes were processed by comparative transcriptomic analysis. Functional characterization showed that phospholipases associated with ROS production played vital roles in early defense responses to H. avenae via involvement in diverse defense-related pathways as a hub switch. This study is the first to investigate the early defense responses of wheat against H. avenae, not only provides applicable candidate resistance genes for breeding novel wheat cultivars, but also enables a better understanding of the defense mechanisms of wheat against H. avenae.

Rice is a staple crop for more than half of the world’s population, and its sustainable production is vital to ensure global food security. However, rice is susceptible to several devastating fungal diseases¹, including blast disease caused by Magnaporthe oryzae, sheath blight by Rhizoctonia solani, false smut by Ustilaginoidea virens, brown spot by Bipolaris oryzae, bakanae by Fusarium fujikuroi and head blight by Fusarium graminearum. The mechanisms underlying the susceptibility to these fungal diseases remain unclear. Here we report that the β subunit of SnRK1, SnRK1β1A, confers broad-spectrum susceptibility to these fungal diseases. Our findings show that diverse rice fungal pathogens have convergently evolved an effector-like protein, Gas2, which interacts with SnRK1β1A to prevent its ubiquitination-mediated degradation and promotes its nuclear translocation. SnRK1β1A is markedly induced on fungal infection, promoting susceptibility by inhibiting SnRK1α1, an α subunit of SnRK1 known to positively regulate broad-spectrum resistance in rice². Notably, rice lines with disrupted SnRK1β1A are resistant to several fungal diseases without compromising growth and yield in the field under normal farming conditions. This study demonstrates that broad-spectrum disease resistance in crops can be achieved by disrupting inducible susceptibility genes whose encoded proteins are targeted by effectors conserved across several pathogens.

Previous studies have shown that the high-osmolarity glycerol mitogen-activated protein kinase (HOG MAPK) signaling pathway and its downstream transcription factor CsAtf1 are involved in the regulation of fludioxonil sensitivity in C. siamense. However, the downstream target genes of CsAtf1 related to the fludioxonil stress response remain unclear. Here, we performed chromatin immunoprecipitation sequencing (ChIP-Seq) and high-throughput RNA-sequencing (RNA-Seq) to identify genome-wide potential CsAtf1 target genes. A total of 3809 significantly differentially expressed genes were predicted to be directly regulated by CsAtf1, including 24 cytochrome oxidase-related genes. Among them, a cytochrome P450-encoding gene, designated CsCyp51G1, was confirmed to be a target gene, and its transcriptional expression was negatively regulated by CsAtf1, as determined using an electrophoretic mobility shift assay (EMSA), a yeast one-hybrid (Y1H) assay, and quantitative real-time PCR (qRT-PCR). Moreover, the overexpression mutant CsCYP51G1 of C. siamense exhibited increased fludioxonil tolerance, and the CsCYP51G1 deletion mutant exhibited decreased fludioxonil resistance, which revealed that CsCyp51G1 is involved in fludioxonil sensitivity regulation in C. siamense. However, the cellular ergosterol content of the mutants was not consistent with the phenotype of fludioxonil sensitivity, which indicated that CsCyp51G1 regulates fludioxonil sensitivity by affecting factors other than the ergosterol level in C. siamense. In conclusion, our data indicate that the transcription factor CsAtf1 negatively regulates the cytochrome P450 gene CsCyp51G1 to increase fludioxonil sensitivity in C. siamense.

Powdery mildew severely affects several important crops and cash plants. Disruption of mildew resistance locus O ( MLO) genes elevates resistance against powdery mildew in several plants. However, whether rubber tree ( Heveae brasiliensis) MLO proteins are linked to susceptibility remains unknown, owing to technical limitations in the genetic manipulation of this woody plant. A previous study showed that the H. brasiliensis MLO-like protein HbMLO12 demonstrates high amino acid sequence similarity with the known Arabidopsis MLO protein AtMLO12. In this study, we investigated whether HbMLO12 regulates susceptibility to powdery mildew. H. brasiliensis leaves take up exogenously synthesized double-stranded RNAs (dsRNAs), and foliar application of dsRNA homologous to HbMLO12 gene specifically induces HbMLO12 silencing in H. brasiliensis leaf tissues. Notably, HbMLO12 silencing inhibited fungal infection and elevated the immune response during interaction with the rubber tree powdery mildew fungus. Furthermore, the heterologous expression of HbMLO12 suppressed bacterial flg22– and fungal chitin–induced immune responses and enhanced bacterial infection in Arabidopsis. Our study provides evidence that HbMLO12 contributes to susceptibility to powdery mildew. [ ] Copyright © 2023 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license .

Rice blast caused by Pyricularia oryzae (syn., Magnaporthe oryzae) was one of the most destructive diseases of rice throughout the world. Genome assembly was fundamental to genetic variation identification and critically impacted the understanding of its ability to overcome host resistance. Here, we report a gapless genome assembly of rice blast fungus P. oryzae strain P131 using PacBio, Illumina and high throughput chromatin conformation capture (Hi-C) sequencing data. This assembly contained seven complete chromosomes (43,237,743 bp) and a circular mitochondrial genome (34,866 bp). Approximately 14.31% of this assembly carried repeat sequences, significantly greater than its previous assembled version. This assembly had a 99.9% complement in BUSCO evaluation. A total of 14,982 genes protein-coding genes were predicted. In summary, we assembled the first telomere-to-telomere gapless genome of P. oryzae, which would be a valuable genome resource for future research on the genome evolution and host adaptation.