Comparative transcriptome analysis of resistant and cultivated tomato lines in response to Clavibacter michiganensis subsp. michiganensis

Unraveling the genetic basis of quantitative resistance to diseases in tomato: a meta-QTL analysis and mining of transcript profiles

KEY MESSAGE

Whole-genome QTL mining and meta-analysis in tomato for resistance to bacterial and fungal diseases identified 73 meta-QTL regions with significantly refined/reduced confidence intervals. Tomato production is affected by a range of biotic stressors, causing yield losses and quality reductions. While sources of genetic resistance to many tomato diseases have been identified and characterized, stability of the resistance genes or quantitative trait loci (QTLs) across the resources has not been determined. Here, we examined 491 QTLs previously reported for resistance to tomato diseases in 40 independent studies and 54 unique mapping populations. We identified 29 meta-QTLs (MQTLs) for resistance to bacterial pathogens and 44 MQTLs for resistance to fungal pathogens, and were able to reduce the average confidence interval (CI) of the QTLs by 4.1-fold and 6.7-fold, respectively, compared to the average CI of the original QTLs. The corresponding physical length of the CIs of MQTLs ranged from 56 kb to 6.37 Mb, with a median of 921 kb, of which 27% had a CI lower than 500 kb and 53% had a CI lower than 1 Mb. Comparison of defense responses between tomato and Arabidopsis highlighted 73 orthologous genes in the MQTL regions, which were putatively determined to be involved in defense against bacterial and fungal diseases. Intriguingly, multiple genes were identified in some MQTL regions that are implicated in plant defense responses, including PR-P2, NDR1, PDF1.2, Pip1, SNI1, PTI5, NSL1, DND1, CAD1, SlACO, DAD1, SlPAL, Ph-3, EDS5/SID1, CHI-B/PR-3, Ph-5, ETR1, WRKY29, and WRKY25. Further, we identified a number of candidate resistance genes in the MQTL regions that can be useful for both marker/gene-assisted breeding as well as cloning and genetic transformation.

Genotyping-by-Sequencing Analysis Reveals Associations between Agronomic and Oil Traits in Gamma Ray-Derived Mutant Rapeseed (Brassica napus L.)

Rapeseed (Brassica napus L.) holds significant commercial value as one of the leading oil crops, with its agronomic features and oil quality being crucial determinants. In this investigation, 73,226 single nucleotide polymorphisms (SNPs) across 95 rapeseed mutant lines induced by gamma rays, alongside the original cultivar ('Tamra'), using genotyping-by-sequencing (GBS) analysis were examined. This study encompassed gene ontology (GO) analysis and a genomewide association study (GWAS), thereby concentrating on agronomic traits (e.g., plant height, ear length, thousand-seed weight, and seed yield) and oil traits (including fatty acid composition and crude fat content). The GO analysis unveiled a multitude of genes with SNP variations associated with cellular processes, intracellular anatomical structures, and organic cyclic compound binding. Through GWAS, we detected 320 significant SNPs linked to both agronomic (104 SNPs) and oil traits (216 SNPs). Notably, two novel candidate genes, Bna.A05p02350D (SFGH) and Bna.C02p22490D (MDN1), are implicated in thousand-seed weight regulation. Additionally, Bna.C03p14350D (EXO70) and Bna.A09p05630D (PI4Kα1) emerged as novel candidate genes associated with erucic acid and crude fat content, respectively. These findings carry implications for identifying superior genotypes for the development of new cultivars. Association studies offer a cost-effective means of screening mutants and selecting elite rapeseed breeding lines, thereby enhancing the commercial viability of this pivotal oil crop.

Transcriptome analysis of tomato plants following salicylic acid-induced immunity against Clavibacter michiganensis ssp. michiganensis

Salicylic acid (SA) is known to be involved in the immunity against Clavibacter michiganensis ssp. michiganensis (Cmm) that causes bacterial canker in tomato. To identify the candidate genes associated with SA-inducible Cmm resistance, transcriptome analysis was conducted via RNA sequencing in tomato plants treated with SA. SA treatment upregulated various defense-associated genes, such as PR and GST genes, in tomato cotyledons. A comparison of SA- and Cmm-responsive genes revealed that both SA treatment and Cmm infection commonly upregulated a large number of genes. Gene Ontology (GO) analysis indicated that the GO terms associated with plant immunity were over-represented in both SA- and Cmm-induced genes. The genes commonly downregulated by both SA treatment and Cmm infection were associated with the cell cycle and may be involved in growth and immunity trade-off through cell division. After SA treatment, several proteins that were predicted to play a role in immune signaling, such as resistance gene analogs, Ca2+ sensors, and WRKY transcription factors, were transcriptionally upregulated. The W-box element, which was targeted by WRKYs, was over-represented in the promoter regions of genes upregulated by both SA treatment and Cmm infection, supporting the speculation that WRKYs are important for the SA-mediated immunity against Cmm. Prediction of protein-protein interactions suggested that genes encoding receptor-like kinases and EF-hand proteins play an important role in immune signaling. Thus, various candidate genes involved in SA-inducible Cmm resistance were identified.

Physiological responses and transcriptome analysis of soybean under gradual water deficit

Soybean is an important food and oil crop widely cultivated globally. However, water deficit can seriously affect the yield and quality of soybeans. In order to ensure the stability and increase of soybean yield and improve agricultural water use efficiency (WUE), research on improving drought tolerance and the efficiency of water utilization of soybeans under drought stress has become particularly important. This study utilized the drought-tolerant variety Heinong 44 (HN44) and the drought-sensitive variety Suinong 14 (SN14) to analyze physiological responses and transcriptome changes during the gradual water deficit at the early seed-filling stage. The results indicated that under drought conditions, HN44 had smaller stomata, higher stomatal density, and lower stomatal conductance (Gs) and transpiration rate as compared to SN14. Additionally, HN44 had a higher abscisic acid (ABA) content and faster changes in stomatal morphology and Gs to maintain a dynamic balance between net photosynthetic rate (Pn) and Gs. Additionally, drought-tolerant variety HN44 had high instantaneous WUE under water deficit. Further, HN44 retained a high level of superoxide dismutase (SOD) activity and proline content, mitigating malondialdehyde (MDA) accumulation and drought-induced damage. Comprehensive analysis of transcriptome data revealed that HN44 had fewer differentially expressed genes (DEGs) under light drought stress, reacting insensitivity to water deficit. At the initial stage of drought stress, both varieties had a large number of upregulated DEGs to cope with the drought stress. Under severe drought stress, HN44 had fewer downregulated genes enriched in the photosynthesis pathway than SN14, while it had more upregulated genes enriched in the ABA-mediated signaling and glutathione metabolism pathways than SN14. During gradual water deficit, HN44 demonstrated better drought-tolerant physiological characteristics and water use efficiency than SN14 through key DEGs such as GmbZIP4, LOC100810474, and LOC100819313 in the major pathways. Key transcription factors were screened and identified, providing further clarity on the molecular regulatory pathways responsible for the physiological differences in drought tolerance among these varieties. This study deepened the understanding of the drought resistance mechanisms in soybeans, providing valuable references for drought-resistant soybean breeding.

Deletion of pbpC Enhances Bacterial Pathogenicity on Tomato by Affecting Biofilm Formation, Exopolysaccharides Production, and Exoenzyme Activities in Clavibacter michiganensis

Penicillin-binding proteins (PBPs) are considered essential for bacterial peptidoglycan biosynthesis and cell wall assembly. Clavibacter michiganensis is a representative Gram-positive bacterial species that causes bacterial canker in tomato. pbpC plays a significant role in maintaining cell morphological characteristics and stress responses in C. michiganensis. The current study demonstrated that the deletion of pbpC commonly enhances bacterial pathogenicity in C. michiganensis and revealed the mechanisms through which this occurs. The expression of interrelated virulence genes, including celA, xysA, xysB, and pelA, were significantly upregulated in △pbpC mutants. Compared with those in wild-type strains, exoenzyme activities, the formation of biofilm, and the production of exopolysaccharides (EPS) were significantly increased in △pbpC mutants. It is noteworthy that EPS were responsible for the enhancement in bacterial pathogenicity, with the degree of necrotic tomato stem cankers intensifying with the injection of a gradient of EPS from C. michiganensis. These findings highlight new insights into the role of pbpC affecting bacterial pathogenicity, with an emphasis on EPS, advancing the current understanding of phytopathogenic infection strategies for Gram-positive bacteria.

Inactivation efficacy of low-pressure plasma treatment against seed-borne tomato pathogen Clavibacter michiganensis and effect of seed setting position and mesh sheet usage

Clavibacter michiganensis, a gram-positive actinomycete, is a major seed-borne tomato pathogen. We investigated the inactivation efficacy of low-pressure plasma treatment against C. michiganensis inoculated on tomato seeds by placing them on a mesh sheet above the bottom dielectric glass plate. The 2- and 5-minute plasma treatment reduced C. michiganensis populations on the tomato seeds by 0.8 and 1.8 log cfu/seed, respectively. The reduction rates were similar to those of C. michiganensis on shirona (cruciferous) seeds, which have different shapes and surface structures. In contrast, the inactivation of C. michiganensis cells using plasma was more difficult than that of X. campestris cells. Additionally, it was found that placing seeds on a mesh sheet laid on the dielectric glass plate was remarkably effective in inactivating the pathogens on tomato seeds. Since the tomato seeds were susceptible to damage from plasma treatment, methods to reduce its damage need to be investigated.

Transcriptome analysis reveals key drought-stress-responsive genes in soybean

Drought is the most common environmental stress and has had dramatic impacts on soybean (Glycine max L.) growth and yield worldwide. Therefore, to investigate the response mechanism underlying soybean resistance to drought stress, the drought-sensitive cultivar "Liaodou 15" was exposed to 7 (mild drought stress, LD), 17 (moderate drought stress, MD) and 27 (severe drought stress, SD) days of drought stress at the flowering stage followed by rehydration until harvest. A total of 2214, 3684 and 2985 differentially expressed genes (DEGs) in LD/CK1, MD/CK2, and SD/CK3, respectively, were identified by RNA-seq. Weighted gene co-expression network analysis (WGCNA) revealed the drought-response TFs such as WRKY (Glyma.15G021900, Glyma.15G006800), MYB (Glyma.15G190100, Glyma.15G237900), and bZIP (Glyma.15G114800), which may be regulated soybean drought resistance. Second, Glyma.08G176300 (NCED1), Glyma.03G222600 (SDR), Glyma.02G048400 (F3H), Glyma.14G221200 (CAD), Glyma.14G205200 (C4H), Glyma.19G105100 (CHS), Glyma.07G266200 (VTC) and Glyma.15G251500 (GST), which are involved in ABA and flavonoid biosynthesis and ascorbic acid and glutathione metabolism, were identified, suggesting that these metabolic pathways play key roles in the soybean response to drought. Finally, the soybean yield after rehydration was reduced by 50% under severe drought stress. Collectively, our study deepens the understanding of soybean drought resistance mechanisms and provides a theoretical basis for the soybean drought resistance molecular breeding and effectively adjusts water-saving irrigation for soybean under field production.

Comprehensive Genome-Wide Analysis and Expression Pattern Profiling of PLATZ Gene Family Members in Solanum Lycopersicum L. under Multiple Abiotic Stresses

PLATZ (plant AT-rich sequence and zinc-binding) family proteins with two conserved zinc-dependent DNA-binding motifs are transcription factors specific to the plant kingdom. The functions of PLATZ proteins in growth, development, and adaptation to multiple abiotic stresses have been investigated in various plant species, but their role in tomato has not been explored yet. In the present work, 20 non-redundant Solanum lycopersicum PLATZ (SlPLATZ) genes with three segmentally duplicated gene pairs and four tandemly duplicated gene pairs were identified on eight tomato chromosomes. The comparative modeling and gene ontology (GO) annotations of tomato PLATZ proteins indicated their probable roles in defense response, transcriptional regulation, and protein metabolic processes as well as their binding affinity for various ligands, including nucleic acids, peptides, and zinc. SlPLATZ10 and SlPLATZ17 were only expressed in 1 cm fruits and flowers, respectively, indicating their preferential involvement in the development of these organs. The expression of SlPLATZ1, SlPLATZ12, and SlPLATZ19 was up- or down-regulated following exposure to various abiotic stresses, whereas that of SlPLATZ11 was induced under temperature stresses (i.e., cold and heat stress), revealing their probable function in the abiotic stress tolerance of tomato. Weighted gene co-expression network analysis corroborated the aforementioned findings by spotlighting the co-expression of several stress-associated genes with SlPLATZ genes. Confocal fluorescence microscopy revealed the localization of SlPLATZ−GFP fusion proteins in the nucleus, hinting at their functions as transcription factors. These findings provide a foundation for a better understanding of the structure and function of PLATZ genes and should assist in the selection of potential candidate genes involved in the development and abiotic stress adaptation in tomato.

Transcriptomic analysis reveals the significant effects of fertilization on the biosynthesis of sesquiterpenes in Phoebe bournei

Phoebe bournei is a potential medicinal plant. Its essential oils (Eos) are mainly composed of sesquiterpenes that has potential activities of anti-bacteria and anti-tumors. In this study, we evaluated the effects of compost and compound fertilizer on the total amount and main components of Eos in P. bournei, we also studied the molecular mechanism undergoing this process by deep sequencing the genes involved in the biosynthesis of sesquiterpenes. Fertilization enhanced the total amount of main components in Eos from both leaves and twigs. Bicyclogermacrene, the primary sesquiterpene in the leaf EO, was significantly increased under compost treatment, while bicyclogermacrene and δ-cadinene (the second most abundant sesquiterpene) were decreased under compound fertilizer treatment. The two fertilizers had no significant effect on the abundance of the primary (+) - δ-cadinene in the twig EO, but had a positive effect on the second most abundant sesquiterpene copaene. Significant differences were observed in the number of differentially expressed genes (DEGs) with the leaves showing greater number of DEGs as compared to the twigs after compost treatment. Terpenoid backbone biosynthesis (TBB) is a key pathway of sesquiterpenes synthesis. The expression of genes regulating several important enzymes in TBB was altered after fertilization. After the compost treatment, the expression of the leaf DXS gene (ACQ66107.1), being closely related to the sesquiterpene biosynthesis in P. bournei leaves, was decreased. Compost and compound fertilizer altered the expression of the two important branch-point enzymes (FPPS and GGPPS) genes (ART33314.1 and ATT59265.1), which contributed to the changes of the total amount and components of P. bournei sesquiterpenes. This study provides a new insight into the future use of P. bournei for Eos.

QTL-Seq Analysis for Identification of Resistance Loci to Bacterial Canker in Tomato

Bacterial canker caused by Clavibacter michiganensis (Cm) is one of the most economically important vascular diseases causing unilateral leaf wilting, stem canker, a bird's-eye lesion on fruit, and whole plant wilting in tomato. There is no commercially available cultivar with bacterial canker resistance, and genomics-assisted breeding can accelerate the development of cultivars with enhanced resistance. Solanum lycopersicum "Hawaii 7998" was found to show bacterial canker resistance. A Quantitative trait loci (QTL)-seq was performed to identify the resistance loci using 909 F₂ individuals derived from a cross between S. lycopersicum "E6203" (susceptible) and "Hawaii 7998," and a genomic region (37.24-41.15 Mb) associated with bacterial canker resistance on chromosome 6 (Rcm6) was found. To dissect the Rcm6 region, 12 markers were developed and several markers were associated with the resistance phenotypes. Among the markers, the Rcm6-9 genotype completely matched with the phenotype in the 47 cultivars. To further validate the Rcm6 as a resistance locus and the Rcm6-9 efficiency, subsequent analysis using F₂ and F₃ progenies was conducted. The progeny individuals with homozygous resistance allele at the Rcm6-9 showed significantly lower disease severity than those possessing homozygous susceptibility alleles. Genomes of five susceptible and two resistant cultivars were analyzed and previously known R-genes were selected to find candidate genes for Rcm6. Nucleotide-binding leucine-rich repeat, receptor-like kinase, and receptor-like protein were identified to have putative functional mutations and show differential expression upon the Cm infection. The DNA markers and candidate genes will facilitate marker-assisted breeding and provide genetic insight of bacterial canker resistance in tomato.

Comparative RNA-Seq Analysis Reveals Potentially Resistance-Related Genes in Response to Bacterial Canker of Tomato

Tomato is one of the most important crops for human consumption. Its production is affected by the actinomycete Clavibacter michiganensis subsp. michiganensis (Cmm), one of the most devastating bacterial pathogens of this crop. Several wild tomato species represent a source of natural resistance to Cmm. Here, we contrasted the transcriptomes of the resistant wild tomato species Solanum arcanum LA2157 and the susceptible species Solanum lycopersicum cv. Ailsa Craig, during the first 24 h of challenge with Cmm. We used three analyses approaches which demonstrated to be complementary: mapping to S. lycopersicum reference genome SL3.0; semi de novo transcriptome assembly; and de novo transcriptome assembly. In a global context, transcriptional changes seem to be similar between both species, although there are some specific genes only upregulated in S. arcanum during Cmm interaction, suggesting that the resistance regulatory mechanism probably diverged during the domestication process. Although S. lycopersicum showed enriched functional groups related to defense, S. arcanum displayed a higher number of induced genes related to bacterial, oomycete, and fungal defense at the first few hours of interaction. This study revealed genes that may contribute to the resistance phenotype in the wild tomato species, such as those that encode for a polyphenol oxidase E, diacyl glycerol kinase, TOM1-like protein 6, and an ankyrin repeat-containing protein, among others. This work will contribute to a better understanding of the defense mechanism against Cmm, and the development of new control methods.

Transcriptome analysis of Clavibacter michiganensis subsp. michiganensis-infected tomatoes: a role of salicylic acid in the host response

Bacterial canker of tomato (Solanum lycopersicon) caused by the Gram-positive bacterium Clavibacter michiganensis subsp. michiganensis (Cmm) is an economically important disease. To understand the host defense response to Cmm infection, transcriptome sequences in tomato cotyledons were analyzed by RNA-seq. Overall, 1788 and 540 genes were upregulated and downregulated upon infection, respectively. Gene Ontology enrichment analysis revealed that genes involved in the defense response, phosphorylation, and hormone signaling were over-represented by the infection. Induced expression of defense-associated genes suggested that the tomato response to Cmm showed similarities to common plant disease responses. After infection, many resistance gene analogs (RGAs) were transcriptionally upregulated, including the expressions of some receptor-like kinases (RLKs) involved in pattern-triggered immunity. The expressions of WRKYs, NACs, HSFs, and CBP60s encoding transcription factors (TFs) reported to regulate defense-associated genes were induced after infection with Cmm. Tomato genes orthologous to Arabidopsis EDS1, EDS5/SID1, and PAD4/EDS9, which are causal genes of salicylic acid (SA)-deficient mutants, were upregulated after infection with Cmm. Furthermore, Cmm infection drastically stimulated SA accumulation in tomato cotyledons. Genes involved in the phenylalanine ammonia lyase pathway were upregulated, whereas metabolic enzyme gene expression in the isochorismate synthase pathway remained unchanged. Exogenously applied SA suppressed bacterial growth and induced the expression of WRKYs, suggesting that some Cmm-responsive genes are regulated by SA signaling, and SA signaling activation should improve tomato immunity against Cmm.

Clavibacter michiganensis Downregulates Photosynthesis and Modifies Monolignols Metabolism Revealing a Crosstalk with Tomato Immune Responses

The gram-positive pathogenic bacterium Clavibacter michiganensis subsp. michiganensis (Cmm) causes bacterial canker disease in tomato, affecting crop yield and fruit quality. To understand how tomato plants respond, the dynamic expression profile of host genes was analyzed upon Cmm infection. Symptoms of bacterial canker became evident from the third day. As the disease progressed, the bacterial population increased in planta, reaching the highest level at six days and remained constant till the twelfth day post inoculation. These two time points were selected for transcriptomics. A progressive down-regulation of key genes encoding for components of the photosynthetic apparatus was observed. Two temporally separated defense responses were observed, which were to an extent interdependent. During the primary response, genes of the phenylpropanoid pathway were diverted towards the synthesis of monolignols away from S-lignin. In dicots, lignin polymers mainly consist of G- and S-units, playing an important role in defense. The twist towards G-lignin enrichment is consistent with previous findings, highlighting a response to generate an early protective barrier and to achieve a tight interplay between lignin recomposition and the primary defense response mechanism. Upon progression of Cmm infection, the temporal deactivation of phenylpropanoids coincided with the upregulation of genes that belong in a secondary response mechanism, supporting an elegant reprogramming of the host transcriptome to establish a robust defense apparatus and suppress pathogen invasion. This high-throughput analysis reveals a dynamic reorganization of plant defense mechanisms upon bacterial infection to implement an array of barriers preventing pathogen invasion and spread.