Reference Gene Selection for qPCR Analysis in Tomato-Bipartite Begomovirus Interaction and Validation in Additional Tomato-Virus Pathosystems

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.

Genome wide inherited modifications of the tomato epigenome by trans-activated bacterial CG methyltransferase

BACKGROUND

Epigenetic variation is mediated by epigenetic marks such as DNA methylation occurring in all cytosine contexts in plants. CG methylation plays a critical role in silencing transposable elements and regulating gene expression. The establishment of CG methylation occurs via the RNA-directed DNA methylation pathway and CG methylation maintenance relies on METHYLTRANSFERASE1, the homologue of the mammalian DNMT1.

PURPOSE

Here, we examined the capacity to stably alter the tomato genome methylome by a bacterial CG-specific M.SssI methyltransferase expressed through the LhG4/pOP transactivation system.

RESULTS

Methylome analysis of M.SssI expressing plants revealed that their euchromatic genome regions are specifically hypermethylated in the CG context, and so are most of their genes. However, changes in gene expression were observed only with a set of genes exhibiting a greater susceptibility to CG hypermethylation near their transcription start site. Unlike gene rich genomic regions, our analysis revealed that heterochromatic regions are slightly hypomethylated at CGs only. Notably, some M.SssI-induced hypermethylation persisted even without the methylase or transgenes, indicating inheritable epigenetic modification.

CONCLUSION

Collectively our findings suggest that heterologous expression of M.SssI can create new inherited epigenetic variations and changes in the methylation profiles on a genome wide scale. This open avenues for the conception of epigenetic recombinant inbred line populations with the potential to unveil agriculturally valuable tomato epialleles.

The Role of CLE Peptides in the Suppression of Mycorrhizal Colonization of Tomato

Symbioses with beneficial microbes are widespread in plants, but these relationships must balance the energy invested by the plants with the nutrients acquired. Symbiosis with arbuscular mycorrhizal (AM) fungi occurs throughout land plants, but our understanding of the genes and signals that regulate colonization levels is limited, especially in non-legumes. Here, we demonstrate that in tomato, two CLV3/EMBRYO-SURROUNDING REGION (CLE) peptides, SlCLE10 and SlCLE11, act to suppress AM colonization of roots. Mutant studies and overexpression via hairy transformation indicate that SlCLE11 acts locally in the root to limit AM colonization. Indeed, SlCLE11 expression is strongly induced in AM-colonized roots, but SlCLE11 is not required for phosphate suppression of AM colonization. SlCLE11 requires the FIN gene that encodes an enzyme required for CLE peptide arabinosylation to suppress mycorrhizal colonization. However, SlCLE11 suppression of AM does not require two CLE receptors with roles in regulating AM colonization, SlFAB (CLAVATA1 ortholog) or SlCLV2. Indeed, multiple parallel pathways appear to suppress mycorrhizal colonization in tomato, as double mutant studies indicate that SlCLV2 and FIN have an additive influence on mycorrhizal colonization. SlCLE10 appears to play a more minor or redundant role, as cle10 mutants did not influence intraradical AM colonization. However, the fact that cle10 mutants had an elevated number of hyphopodia and that ectopic overexpression of SlCLE10 did suppress mycorrhizal colonization suggests that SlCLE10 may also play a role in suppressing AM colonization. Our findings show that CLE peptides regulate AM colonization in tomato and at least SlCLE11 likely requires arabinosylation for activity.

Immunity priming uncouples the growth-defense trade-off in tomato

Plants have developed an array of mechanisms to protect themselves against pathogen invasion. The deployment of defense mechanisms is imperative for plant survival, but can come at the expense of plant growth, leading to the 'growth-defense trade-off' phenomenon. Following pathogen exposure, plants can develop resistance to further attack. This is known as induced resistance, or priming. Here, we investigated the growth-defense trade-off, examining how defense priming via systemic acquired resistance (SAR), or induced systemic resistance (ISR), affects tomato development and growth. We found that defense priming can promote, rather than inhibit, plant development, and that defense priming and growth trade-offs can be uncoupled. Cytokinin response was activated during induced resistance, and found to be required for the observed growth and disease resistance resulting from ISR activation. ISR was found to have a stronger effect than SAR on plant development. Our results suggest that growth promotion and induced resistance can be co-dependent, and that, in certain cases, defense priming can drive developmental processes and promote plant yield.

The selection and validation of reference genes for quantitative real-time PCR studies in near-isogenic susceptible and resistant tomato lines, infected with the geminivirus tomato curly stunt virus

Quantitative real-time PCR (qPCR) is a sensitive and commonly used technique for gene expression profiling and provides insight into biological systems. Successful qPCR requires the use of appropriate reference genes for the normalization of data. In the present study, we aimed to identify and assess the best-suited reference genes in near-isogenic resistant (R) and susceptible (S) tomato lines infected with begomovirus Tomato curly stunt virus (ToCSV). Ten candidate reference genes namely, Actin7 (ACT), β-6 Tubulin (TUB), Ubiquitin 3 (UBI), Clathrin adaptor complexes medium subunit (CAC), Phytoene desaturase (PDS), Expressed protein (EXP), Glyceraldehyde-3-phosphate dehydrogenase (GAPDH), Adenine phosphoribosyl transferase-like protein (APT1), TAP42-interacting protein (TIP41) and Elongation factor 1-alpha (EF1α) were selected and evaluated for their expression stability in resistant and susceptible tomato leaves using the analytical tools geNorm, NormFinder, BestKeeper, and RefFinder. After ranking the reference genes from most to least stable, the results suggested that a combination of ACT, EXP, and EF1α in the S lines and a combination of TIP41, APT1, and ACT in the R line is appropriate for qPCR normalization. Furthermore, to validate the identified reference genes, iron superoxide dismutase (SOD), heat shock protein 70 (HSP70) and Glutathione-S-transferase (GST) were selected as targets for normalization. The relative expression of the target genes varied when normalized against the most stable reference genes in comparison to the least stable genes. These results highlight the importance of careful selection of reference genes for accurate normalization in qPCR studies.

Transcriptome analysis during ToLCBaV disease development in contrasting tomato genotypes

Tomato leaf curl Bangalore virus (ToLCBaV) is one of the most important plant viruses. The infection causes substantial yield losses in tomato crop. The current viral disease management is based mainly on introgression of Ty locus into new tomato cultivars. Unfortunately, strains of the leaf curl virus have been evolving and are breaking Ty based tolerance in tomato. In this study, the defence response to ToLCBaV infection has been compared between contrasting tomato genotypes, resistant line (IIHR 2611; without any known Ty markers) and the susceptible line (IIHR 2843). We carried out comparative transcriptome profiling, and gene expression analysis in an effort to identify gene networks that are associated with a novel ToLCBaV resistance. A total of 22,320 genes were examined to identify differentially expressed genes (DEGs). We found that 329 genes of them were expressed significantly and differentially between ToLBaV-infected samples of both IIHR 2611 and IIHR 2843. A good number of DEGs were related to defence response, photosynthesis, response to wounding, toxin catabolic process, glutathione metabolic process, regulation of transcription DNA-template, transcription factor activity, and sequence-specific DNA binding. A few selected genes such as, nudix hydrolase 8, MIK 2-like, RING-H2 finger protein ATL2-like, MAPKKK 18-like, EDR-2, SAG 21 wound-induced basic protein, GRXC6 and P4 were validated using qPCR. The pattern of gene expression was significantly different in resistant and susceptible plants during disease progression. Both positive and negative regulators of virus resistance were found in the present study. These findings will facilitate breeding and genetic engineering efforts to incorporate novel sources of ToLCBaV resistance in tomatoes.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-023-03629-5.

Evaluation of reference genes for miRNA and mRNA normalization in tobacco infected with PVYNTN, PVYN-Wi and PVYZ-NTN strains

This is the first report on identification of the most suitable reference genes for RT-qPCR quantification of miRNA and mRNA in tobacco response to the prevalent recombinant potato virus Y (PVY) strains PVY^NTN, PVY^N-Wi and the newly identified PVY^Z-NTN. Of 10 tested genes, the expression levels of neIF5C, nU2af and nPP2A were the most stable in samples taken from non-inoculated, mock-inoculated, and infected plants at three days post-inoculation (dpi) and 14 dpi. While the homologues of eIF5 were most stably expressed in tobacco in this study and in potato in our previous study (Yin et al., 2021) following inoculation with the same three PVY strains, the homologues of other two genes PP2A and U2af were stably expressed only in tobacco but unstable in potato. The tobacco homologue of PP2A, which was the most stably expressed one in tobacco interaction with PVY^NTN, PVY^N-Wi and PVY^Z-NTN strains in this study, was the least stable one in tobacco interaction with the non-recombinant PVY^O strain in a previous study (Baek et al., 2017). This study provides evidence on the influence of host species on expression of housekeeping genes and points out virus strain as a new factor influencing expression stability of reference gene. Caution should be taken when choosing reference genes in gene expression study in Solanaceae hosts response to different PVY strains.

The Impact of Tobamovirus Infection on Root Development Involves Induction of Auxin Response Factor 10a in Tomato

Plant viruses cause systemic diseases that severely impair plant growth and development. While the accumulation of viruses in the root system has long been established, little is known as to how viruses affect root architecture. Here, we examined how the emerging tobamovirus, tomato brown rugose fruit virus (ToBRFV), alters root development in tomato. We found that ToBRFV and tobacco mosaic virus both invaded root systems during the first week of infection. ToBRFV infection of tomato plants resulted in a significant decrease in root biomass and elongation and root-to-shoot ratio and a marked suppression of root branching. Mutation in RNA-dependent RNA polymerase 6 increased the susceptibility of tomato plants to ToBRFV, resulting in severe reduction of various root growth parameters including root branching. Viral root symptoms were associated with the accumulation of auxin response factor 10a (SlARF10a) transcript, a homolog of Arabidopsis ARF10, a known suppressor of lateral root development. Interestingly, loss-of-function mutation in SlARF10a moderated the effect of ToBRFV on root branching. In contrast, downregulation of sly-miR160a, which targets SlARF10a, was associated with constitutive suppression root branching independent of viral infection. In addition, overexpression of a microRNA-insensitive mutant of SlARF10a mimicked the effect of ToBRFV on root development, suggesting a specific role for SlARF10a in ToBRFV-mediated suppression of root branching. Taken together, our results provide new insights into the impact of tobamoviruses on root development and the role of ARF10a in the suppression of root branching in tomato.

The tomato resistance gene Bs4 suppresses leaf watersoaking phenotypes induced by AvrHah1, a transcription activator-like effector from tomato-pathogenic xanthomonads

The Xanthomonas transcription activator-like effector (TALE) protein AvrBs3 transcriptionally activates the executor-type resistance (R) gene Bs3 from pepper (Capsicum annuum), thereby triggering a hypersensitive cell death reaction (HR). AvrBs3 also triggers an HR in tomato (Solanum lycopersicum) upon recognition by the nucleotide-binding leucine-rich repeat (NLR) R protein Bs4. Whether the executor-type R protein Bs3 and the NLR-type R protein Bs4 use common or distinct signalling components to trigger an HR remains unclear. CRISPR/Cas9-mutagenesis revealed, that the immune signalling node EDS1 is required for Bs4- but not for Bs3-dependent HR, suggesting that NLR- and executor-type R proteins trigger an HR via distinct signalling pathways. CRISPR/Cas9-mutagenesis also revealed that tomato Bs4 suppresses the virulence function of both TALEs, the HR-inducing AvrBs3 protein and of AvrHah1, a TALE that does not trigger an HR in tomato. Analysis of AvrBs3- and AvrHah1-induced host transcripts and disease phenotypes in CRISPR/Cas9-induced bs4 mutant plants indicates that both TALEs target orthologous transcription factor genes to promote disease in tomato and pepper host plants. Our studies display that tomato mutants lacking the TALE-sensing Bs4 protein provide a novel platform to either uncover TALE-induced disease phenotypes or genetically dissect components of executor-triggered HR.

Molecular and Physiological Effects of Magnesium-Polyphenolic Compound as Biostimulant in Drought Stress Mitigation in Tomato

Plant biostimulants are being recognized as innovative tools to improve sustainable agricultural practices to mitigate the drastic effects of climate change, which is leading to a severe reduction in agricultural yields. In this work, a new biostimulant (EnNuVi^® ALPAN^®) was evaluated for its effectiveness on tomato (Solanum lycopersicum Mill. cv. Rio Grande) plants subjected to water deficit conditions. The molecular effects were elucidated through transcriptomic RNA-seq and gene expression qPCR analysis and the physiological responses were evaluated through qualitative analysis of pigments and proline content, membrane stability, and lipid peroxidation. ALPAN^® was shown to adjust the transcriptional response by upregulating genes involved in source to sink carbohydrate metabolism and translocation, stomatal closure, and cell homeostasis. ALPAN^® was shown to mitigate the deteriorating effects of water deficit on the physiological status of the plants by stabilizing the levels of the photosynthetic pigments, regulating the accumulation of osmo-protectants, and preserving the cell wall lipid bilayer from oxidation. In conclusion, transcriptomic and physiological analysis provided insightful information on the biostimulant effects, indicating a positive role of ALPAN^® foliar application in alleviating the negative costs of water deficit.

Deciphering fruit sugar transport and metabolism from tolerant and sensitive tomato plants subjected to simulated field conditions

In the current state of climate change, we must assume that abiotic stresses act together under natural field conditions, these will increase in the coming years. Therefore, in this report we investigated how sugar metabolism was affected under simulated field conditions, where plants faced high ambient temperatures and a low-quality water irrigation. Our studies were carried out on fruits of two tomato recombinant lines, a tolerant and a sensitive one exposed to the combination of heat and salinity. Two ripening stages (mature green and red ripe fruits) were used in our analyzes, where the gene expression levels of the main biosynthetic genes and transporters, enzymatic activities and compounds related to the synthesis, accumulation, and degradation of sugars in plants were analyzed. The tolerant line showed highly significant differences in red ripe fruits in comparison to the sensitive one under the simulated field conditions (35°C + 60 mM NaCl), with an overexpression of the genes SlFBP, SlSPS, SlSUS3, and SlNi. These expression patterns correlated with a higher activity of the enzymes FBP, SPS, SUS3, AI, and G6PDH, which resulted in the accumulation of fructose, glucose and UDP-glucose. Our results showed the advantage of using tomato recombinant lines for rescuing important traits, such as the resistance to some abiotic stresses, and for the identification of important molecular and metabolic markers that could be used to determine fruit quality in green or red maturity stages under detrimental environmental field conditions.

Stability analysis of reference genes for RT-qPCR assays involving compatible and incompatible Ralstonia solanacearum-tomato 'Hawaii 7996' interactions

Reverse transcription-quantitative PCR (RT-qPCR) is an analytical tool for gene expression quantification. Reference genes are not yet available for gene expression analysis during interactions of Ralstonia solanacearum with ‘Hawaii 7996’ (the most stable source of resistance in tomato). Here, we carried out a multi-algorithm stability analysis of eight candidate reference genes during interactions of ‘Hawaii 7996’ with one incompatible/avirulent and two compatible/virulent (= resistance-breaking) bacterial isolates. Samples were taken at 24- and 96-h post-inoculation (HPI). Analyses were performed using the ∆∆Ct method and expression stability was estimated using BestKeeper, NormFinder, and geNorm algorithms. TIP41 and EF1α (with geNorm), TIP41 and ACT (with NormFinder), and UBI3 and TIP41 (with BestKeeper), were the best combinations for mRNA normalization in incompatible interactions at 24 HPI and 96 HPI. The most stable genes in global compatible and incompatible interactions at 24 HPI and 96 HPI were PDS and TIP41 (with geNorm), TIP41 and ACT (with NormFinder), and UBI3 and PDS/EXP (with BestKeeper). Global analyses on the basis of the three algorithms across 20 R. solanacearum-tomato experimental conditions identified UBI3, TIP41 and ACT as the best choices as reference tomato genes in this important pathosystem.

The Tomato Brown Rugose Fruit Virus Movement Protein Overcomes Tm-22 Resistance in Tomato While Attenuating Viral Transport

Tomato brown rugose fruit virus is a new virus species in the Tobamovirus genus, causing substantial damage to tomato crops. Reports of recent tomato brown rugose fruit virus (ToBRFV) outbreaks from around the world indicate an emerging global epidemic. ToBRFV overcomes all tobamovirus resistances in tomato, including the durable Tm-2² resistance gene, which had been effective against multiple tobamoviruses. Here, we show that the ToBRFV movement protein (MP^ToBRFV) enables the virus to evade Tm-2² resistance. Transient expression of MP^ToBRFV failed to activate the Tm-2² resistance response. Replacement of the original MP sequence of tomato mosaic virus (ToMV) with MP^ToBRFV enabled this recombinant virus to infect Tm-2²-resistant plants. Using hybrid protein analysis, we show that the elements required to evade Tm-2² are located between MP^ToBRFV amino acids 1 and 216 and not the C terminus, as previously assumed. Analysis of ToBRFV systemic infection in tomato revealed that ToBRFV spreads more slowly compared with ToMV. Interestingly, replacement of tobacco mosaic virus (TMV) and ToMV MPs with MP^ToBRFV caused an attenuation of systemic infection of both viruses. Cell-to-cell movement analysis showed that MP^ToBRFV moves less effectively compared with the TMV MP (MP^TMV). These findings suggest that overcoming Tm-2² is associated with attenuated MP function. This may explain the high durability of Tm-2² resistance, which had remained unbroken for over 60 years.[Formula: see text] Copyright © 2021 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Gene expression analysis for selection and validation of suitable housekeeping gene(s) in cadmium exposed pigeonpea plants inoculated with arbuscular mycorrhizae

The expression stability of six commonly used housekeeping genes (18S rRNA-18S ribosomal RNA, EF1α-elongation factor 1α, ACT1-Actin 1, GAPDH-Glyceraldehyde-3-phosphate dehydrogenase, TUB6-Tubulin/FtsZ family and UBC-Ubiquitin-conjugating enzyme) were scrutinized in leaves and roots of Cd stressed pigeonpea plants inoculated with arbuscular mycorrhizal (AM) species- Rhizoglomus intraradices (Ri), Funneliformis mosseae (Fm), Claroideoglomus etunicatum (Ce), C. claroideum (Cc). The stability profile of each gene was assessed using ΔC_t, BestKeeper, NormFinder, RefFinder and geNorm algorithmic programs, which ranked different genes as most and least stable according to the tissues analysed. All the statistical algorithms ranked TUB6 as most stable and EF1α least stable housekeeping (HK) genes in both the plant tissues. The selected HK genes were verified using metallothionein (CcMT1) i.e. a stress responsive gene, whose expression altered under conditions of metal stress and AM inoculation. The expression pattern of CcMT1 varied highly when least stable reference gene was used for normalization as compared to most stable gene, under different treatments. Thus, there is a need of selecting suitable reference gene to achieve reliable results in gene expression studies using quantitative real time PCR (qRT-PCR). The study conducted will help future gene expression analysis in pigeonpea under specific stress.

Selection and optimization of reference genes for RT-qPCR normalization: A case study in Solanum lycopersicum exposed to UV-B

Quantitative RT- PCR is one of the most common methods to study gene expression in response to stress. Therefore, it is crucial to have suitable reference genes (RGs) for result normalization. Although several reports describe UV-B-modulated gene expression in Solanum lycopersicum, there are no suitable RGs identified until now. The aim of this work was to evaluate the suitability of seven traditional genes: actin (ACT), tubulin (TUB), ubiquitin (UBI), glyceraldehyde- 3 phosphate dehydrogenase (GAPDH), elongation factor 1α (EF1α), phosphatase 2A catalytic subunit (PP2A) and GAGA binding transcriptional activator (GAGA); and two non-traditional genes: thioredoxin h1 (TRX h1) and UV-B RESISTANCE LOCUS 8 (UVR8), as candidate RGs for their potential use as reliable internal controls in leaves, stems and roots of tomato seedlings exposed to acute and chronic UV-B. The stability of these genes expression was evaluated using five statistical algorithms: geNorm, NormFinder, BestKeeper, Delta Ct and ANOVA. Considering the comprehensive stability ranking, we recommend ACT+TUB as the best pair of RGs for leaves, PP2A+GAPDH+TRX h1 for stems and TUB+UVR8 for roots. The reliability of the selected RGs for each tissue was verified amplifying tomato chalcone synthase 1 (CHS1) and cyclobutane pyrimidine dimer (CPD) photolyase (PHR1-LIKE). Under UV-B treatment, CHS1 was upregulated in leaves, stems and roots whereas PHR1-LIKE was only upregulated in leaves and stems. This interpretation differs when the most and least stable RGs are chosen. This is the first report regarding suitable RGs selection for accurate normalization of gene expression in tomato seedlings exposed to UV-B irradiation.

Reference gene selection for miRNA and mRNA normalization in potato in response to potato virus Y

This was the first report on evaluating candidate reference genes for quantifying the expression profiles of both coding (e.g., mRNA) and non-coding (e.g., miRNA) genes in potato response to potato virus Y (PVY) inoculation. The reverse transcription-quantitative real-time polymerase chain reaction (RT-qPCR) method was employed to quantify the expression profiles of eight selected candidate reference genes; their expression stability was analyzed by four statistical algorithms, i.e., geNorm, BestKeeper, NormFinder and RefFinder. The most stable reference genes were sEF1a, sTUBb and seIF5 with a high stability. The least stable ones were sPP2A, sSUI1 and sGAPDH. The same reference gene allows for normalization of both miRNA and mRNA levels from a single RNA sample using cDNAs synthesized in a single RT reaction, in which a stem-loop primer was used for miRNAs and the oligo (dT) for mRNAs.

Selection and validation of optimal endogenous reference genes for analysis of quantitative PCR in four tissues pathologically associated with Kidney-yang deficiency syndrome following influenza A infection

Kidney-yang deficiency syndrome (KYDS) infected with the influenza virus is a suitable model to imitate a population at high-risk to influenza infection with a high rate of morbidity and mortality. However, the specific molecular mechanisms underlying this disease remain unclear. A stable reference gene is essential as an internal control for molecular biology research of this condition. Reverse-transcription-quantitative PCR (RT-qPCR) is considered an extremely sensitive technique used for absolute and relative quantification of target genes transcript levels. To accurately estimate the relative expression of genes in cells from mice with KYDS in response to infection with influenza A virus subtype H1N1 (A/H1N1) virus using RT-qPCR, it is necessary to identify suitable reference genes. In the present study, analysis of 10 reference genes (Act-β, β2m, GAPDH, Gusβ, Tubα, Grcc10, Eif4h, Rnf187, Nedd8 and Ywhae) was performed across a set of 4 tissue types: Lung; heart; liver; and kidney. KYDS mice were inoculated with A/H1N1 virus or a mock control. For analysis, geNorm, BestKeeper, NormFinder, and Bio-Rad Maestro^™ statistical programs were used for the estimation of the stability of the reference genes. The results were authenticated through extended experimental settings using a group of 10 samples, parallel to 3 additional innate immune system-associated genes of the host, TLR3, TLR7 and RIG-I, which were also analyzed using the same algorithms. From the 4 algorithms, taking into account the joint analyses of the ranking order outputs, the 2 genes Ywhae and Nedd8 were identified to be the most stable for mice with KYDS following infection with A/H1N1 virus. In contrast, the least stable genes in all 4 tissues were GAPDH and β2m. These results may affect the choice of reference genes in future studies that use RT-qPCR analysis of target genes in experimental conditions, such as mice with KYDS infected with influenza A virus.

Molecular insights into biochar-mediated plant growth promotion and systemic resistance in tomato against Fusarium crown and root rot disease

Molecular mechanisms associated with biochar-elicited suppression of soilborne plant diseases and improved plant performance are not well understood. A stem base inoculation approach was used to explore the ability of biochar to induce systemic resistance in tomato plants against crown rot caused by a soilborne pathogen, Fusarium oxysporum f. sp. radicis lycopersici. RNA-seq transcriptome profiling of tomato, and experiments with jasmonic and salycilic acid deficient tomato mutants, were performed to elucidate the in planta molecular mechanisms involved in induced resistance. Biochar (produced from greenhouse plant wastes) was found to mediate systemic resistance against Fusarium crown rot and to simultaneously improve tomato plant growth and physiological parameters by up to 63%. Transcriptomic analysis (RNA-seq) of tomato demonstrated that biochar had a priming effect on gene expression and upregulated the pathways and genes associated with plant defense and growth such as jasmonic acid, brassinosteroids, cytokinins, auxin and synthesis of flavonoid, phenylpropanoids and cell wall. In contrast, biosynthesis and signaling of the salicylic acid pathway was downregulated. Upregulation of genes and pathways involved in plant defense and plant growth may partially explain the significant disease suppression and improvement in plant performance observed in the presence of biochar.

Tomato Brown Rugose Fruit Virus Contributes to Enhanced Pepino Mosaic Virus Titers in Tomato Plants

The tobamovirus tomato brown rugose fruit virus (ToBRFV), a major threat to tomato production worldwide, has recently been documented in mixed infections with the potexvirus pepino mosaic virus (PepMV) CH2 strain in traded tomatoes in Israel. A study of greenhouse tomato plants in Israel revealed severe new viral disease symptoms including open unripe fruits and yellow patched leaves. PepMV was only detected in mixed infections with ToBRFV in all 104 tested sites, using serological and molecular analyses. Six PepMV isolates were identified, all had predicted amino acids characteristic of CH2 mild strains excluding an isoleucine at amino acid position 995 of the replicase. High-throughput sequencing of viral RNA extracted from four selected symptomatic plants showed solely the ToBRFV and PepMV, with total aligned read ratios of 40.61% and 11.73%, respectively, indicating prevalence of the viruses. Analyses of interactions between the co-infecting viruses by sequential and mixed viral inoculations of tomato plants, at various temperatures, showed a prominent increase in PepMV titers in ToBRFV pre-inoculated plants and in mixed-infected plants at 18–25 °C, compared to PepMV-single inoculations, as analyzed by Western blot and quantitative RT-PCR tests. These results suggest that Israeli mild PepMV isolate infections, preceded by ToBRFV, could induce symptoms characteristic of PepMV aggressive strains.

Real-time PCR data for reference candidate gene selection in tomato infected with Tomato curly stunt virus

Real-time PCR (qPCR) is a useful and robust method of quantifying gene expression, provided that suitable reference genes are used to normalize the data. To date, suitable reference genes have not been validated for tomato gene expression changes in response to Tomato curly stunt virus (ToCSV). RT-qPCR was conducted on resistent (R) and susceptible (S) tomato leave tissue infected with ToCSV at 35 days post infection. Ten candidate reference genes were selected and validated using SYBR green. Here, we report a set of primers designed for the ten candidate genes and the data for the melt curve analysis and standard curves generated for each candidate reference gene. This data provides a useful resourse in reference gene selection for future use in the normalization of qPCR data investigating tomato-virus interactions. To our knowledge, this data provides the first selection and testing of candidate reference genes in a tomato-ToCSV pathosystem.

HEBE, a novel positive regulator of senescence in Solanum lycopersicum

Leaf senescence and plant aging are traits of great interest for breeders. Senescing cells undergo important physiological and biochemical changes, while cellular structures such as chloroplasts are degraded with dramatic metabolic consequences for the whole plant. The possibility of prolonging the photosynthetic ability of leaves could positively impact the plant's life span with benefits for biomass production and metabolite accumulation; plants with these characteristics display a stay-green phenotype. A group of plant transcription factors known as NAC play a pivotal role in controlling senescence: here we describe the involvement of the tomato NAC transcription factor Solyc12g036480, which transcript is present in leaves and floral buds. Since its silencing delays leaf senescence and prevents plants from ageing, we renamed Solyc12g0364 HḖBĒ, for the Greek goddess of youth. In this manuscript we describe how HEB downregulation negatively affects the progression of senescence, resulting in changes in transcription of senescence-promoting genes, as well as the activity of enzymes involved in chlorophyll degradation, thereby explaining the stay-green phenotype.

Selection of reference genes for normalization of RT-qPCR data in gene expression studies in Anthonomus eugenii Cano (Coleoptera: Curculionidae)

The pepper weevil, Anthonomus eugenii Cano (Coleoptera: Curculionidae), is the main insect pest of peppers (Capsicum spp.) throughout the southern U.S. and a potential target for novel control methods that may require gene expression analyses. Careful selection of adequate reference genes to normalize RT-qPCR data is an important prerequisite for gene expression studies since the expression stability of reference genes can be affected by the experimental conditions leading to biased or erroneous results. The lack of studies on validation of reference genes for RT-qPCR analysis in A. eugenii limits the investigation of gene expression, therefore it is needed a systematic selection of suitable reference genes for data normalization. In the present study, three programs (BestKeeper, geNorm and NormFinder) were used to analyze the expression stability of candidate reference genes (β-ACT, ArgK, EF1-α, GAPDH, RPL12, RPS23, α-TUB, 18S and 28S) in A. eugenii under different experimental conditions. Our results revealed that the most stably expressed reference genes in A. eugenii varied according to the experimental condition evaluated: developmental stages (EF1-α, 18S and RPL12), sex (RPS23 and RPL12), low temperature (GAPDH and α-TUB), high temperature (α-TUB and RPS23), all temperatures (α-TUB and GAPDH), starvation (RPL12 and α-TUB), and dsRNA exposure (α-TUB and RPL12). Our study provides for the first time valuable information on appropriate reference genes that can be used in the analysis of gene expression by RT-qPCR in biological experiments involving A. eugenii.

Identification of two new trichome-specific promoters of Nicotiana tabacum

MAIN CONCLUSION: pRbcS-T1 and pMALD1, two new trichome-specific promoters of Nicotiana tabacum, were identified and their strength and specificity were compared to those of previously described promoters in this species. Nicotiana tabacum has emerged as a suitable host for metabolic engineering of terpenoids and derivatives in tall glandular trichomes, which actively synthesize and secrete specialized metabolites. However, implementation of an entire biosynthetic pathway in glandular trichomes requires the identification of trichome-specific promoters to appropriately drive the expression of the transgenes needed to set up the desired pathway. In this context, RT-qPCR analysis was carried out on wild-type N. tabacum plants to compare the expression pattern and gene expression level of NtRbcS-T1 and NtMALD1, two newly identified genes expressed in glandular trichomes, with those of NtCYP71D16, NtCBTS2α, NtCPS2, and NtLTP1, which were reported in the literature to be specifically expressed in glandular trichomes. We show that NtRbcS-T1 and NtMALD1 are specifically expressed in glandular trichomes like NtCYP71D16, NtCBTS2α, and NtCPS2, while NtLTP1 is also expressed in other leaf tissues as well as in the stem. Transcriptional fusions of each of the six promoters to the GUS-VENUS reporter gene were introduced in N. tabacum by Agrobacterium-mediated transformation. Almost all transgenic lines displayed GUS activity in tall glandular trichomes, indicating that the appropriate cis regulatory elements were included in the selected promoter regions. However, unlike for the other promoters, no trichome-specific line was obtained for pNtLTP1:GUS-VENUS, in agreement with the RT-qPCR data. These data thus provide two new transcription promoters that could be used in metabolic engineering of glandular trichomes.

Using Tomato Recombinant Lines to Improve Plant Tolerance to Stress Combination Through a More Efficient Nitrogen Metabolism

The development of plant varieties with a better nitrogen use efficiency (NUE) is a means for modern agriculture to decrease environmental pollution due to an excess of nitrate and to maintain a sufficient net income. However, the optimum environmental conditions for agriculture will tend to be more adverse in the coming years, with increases in temperatures, water scarcity, and salinity being the most important productivity constrains for plants. NUE is inherently a complex trait, as each step, including N uptake, translocation, assimilation, and remobilization, is governed by multiple interacting genetic and environmental factors. In this study, two recombinant inbred lines (RIL-66 and RIL-76) from a cross between Solanum lycopersicum and Solanum pimpinellifoilum with different degree of tolerance to the combination of salinity and heat were subjected to a physiological, ionomic, amino acid profile, and gene expression study to better understand how nitrogen metabolism is affected in tolerant plants as compared to sensitive ones. The ionomics results showed a different profile between the two RILs, with K+ and Mg2+ being significantly lower in RIL-66 (low tolerant) as compared to RIL-76 (high tolerant) under salinity and heat combination. No differences were shown between the two RILs in N total content; however, N-NO3 - was significantly higher in RIL-66, whereas N-Norg was lower as compared to the other genotype, which could be correlated with its tolerance to the combination of salinity and heat. Total proteins and total amino acid concentration were significantly higher in RIL-76 as compared to the sensitive recombinant line under these conditions. Glutamate, but more importantly glutamine, was also highly synthesized and accumulated in RIL-76 under the combination of salinity and heat, which was in agreement with the upregulation of the nitrogen metabolism related transcripts studied (SlNR, SlNiR, SlGDH, SlGLT1, SlNRT1.2, SlAMT1, and SlAMT2). This study emphasized the importance of studying abiotic stress in combination and how recombinant material with different degrees of tolerance can be highly important for the improvement of nitrogen use efficiency in horticultural plants through the targeting of N-related markers.

Expression of Putative Defense Responses in Cannabis Primed by Pseudomonas and/or Bacillus Strains and Infected by Botrytis cinerea

Cannabis (Cannabis sativa L.) offers many industrial, agricultural, and medicinal applications, but is commonly threatened by the gray mold disease caused by the fungus Botrytis cinerea. With few effective control measures currently available, the use of beneficial rhizobacteria represents a promising biocontrol avenue for cannabis. To counter disease development, plants rely on a complex network of inducible defense pathways, allowing them to respond locally and systemically to pathogens attacks. In this study, we present the first attempt to control gray mold in cannabis using beneficial rhizobacteria, and the first investigation of cannabis defense responses at the molecular level. Four promising Pseudomonas (LBUM223 and WCS417r) and Bacillus strains (LBUM279 and LBUM979) were applied as single or combined root treatments to cannabis seedlings, which were subsequently infected by B. cinerea. Symptoms were recorded and the expression of eight putative defense genes was monitored in leaves by reverse transcription quantitative polymerase chain reaction. The rhizobacteria did not significantly control gray mold and all infected leaves were necrotic after a week, regardless of the treatment. Similarly, no systemic activation of putative cannabis defense genes was reported, neither triggered by the pathogen nor by the rhizobacteria. However, this work identified five putative defense genes (ERF1, HEL, PAL, PR1, and PR2) that were strongly and sustainably induced locally at B. cinerea's infection sites, as well as two stably expressed reference genes (TIP41 and APT1) in cannabis. These markers will be useful in future researches exploring cannabis defense pathways.

Reference genes for gene expression studies targeting sugarcane infected with Sugarcane mosaic virus (SCMV)

Objective

The selection of reference genes in sugarcane under Sugarcane mosaic virus (SCMV) infection has not been reported and is indispensable to get reliable reverse transcription quantitative PCR (RT-qPCR) results for validation of transcriptome analysis. In this regard, seven potential reference genes were tested by RT-qPCR and ranked according to their stability using BestKeeper, NormFinder and GeNorm algorithms, and RefFinder WEB-based software in an experiment performed with samples from two sugarcane cultivars contrasting for SCMV resistance, when mechanically inoculated with a severe SCMV strain and using mock inoculated plant controls.

Results

The genes Uridylate kinase (UK) and Ubiquitin-conjugating enzyme 18 (UBC18) were the most stable according to GeNorm algorithm and the Pearson correlation coefficients with the BestKeeper index. On the other hand, ribosomal protein L35-4 (RPL1), Actin (ACT) and Ubiquitin1 (UBQ1) were the least stable genes for all algorithms tested.

Electronic supplementary material

The online version of this article (10.1186/s13104-019-4168-5) contains supplementary material, which is available to authorized users.

Transcriptome-based identification and validation of reference genes for plant-bacteria interaction studies using Nicotiana benthamiana

RT-qPCR is a widely used technique for the analysis of gene expression. Accurate estimation of transcript abundance relies strongly on a normalization that requires the use of reference genes that are stably expressed in the conditions analyzed. Initially, they were adopted from those used in Northern blot experiments, but an increasing number of publications highlight the need to find and validate alternative reference genes for the particular system under study. The development of high-throughput sequencing techniques has facilitated the identification of such stably expressed genes. Nicotiana benthamiana has been extensively used as a model in the plant research field. In spite of this, there is scarce information regarding suitable RT-qPCR reference genes for this species. Employing RNA-seq data previously generated from tomato plants, combined with newly generated data from N. benthamiana leaves infiltrated with Pseudomonas fluorescens, we identified and tested a set of 9 candidate reference genes. Using three different algorithms, we found that NbUbe35, NbNQO and NbErpA exhibit less variable gene expression in our pathosystem than previously used genes. Furthermore, the combined use of the first two is sufficient for robust gene expression analysis. We encourage employing these novel reference genes in future RT-qPCR experiments involving N. benthamiana and Pseudomonas spp.

Evaluation of internal control genes for quantitative realtime PCR analyses for studying fruit development of dwarf tomato cultivar 'Micro-Tom'

Quantitative real-time PCR (qRT-PCR) is widely used to analyze the expression profiles of the genes of interest. In order to obtain accurate quantification data, normalization by using reliable internal control genes is essential. In this study, we evaluated the stability and applicability of eight internal control gene candidates for analyzing gene expression during fruit development in dwarf tomato cultivar Micro-Tom. We collected seventeen different samples from flowers and fruits at different developmental stages, and estimated the expression stability of the candidate genes by two statistical algorithms, geNorm and NormFinder. The combined ranking order and qRT-PCR analyses for expression profiles of SlYABBY2a, SlYABBY1a, FRUITFULL1 and APETALA2c suggested that EXPRESSED was the most stable and reliable internal control gene among the candidates. Our analysis also suggested that RPL8 was also suitable if the sample group is limited to fruits at different maturation stages. In addition to EXPRESSED, GAPDH was also applicable for relative quantitation to monitor gene expression profiles through fruit development from pistil to pericarp.

Selection of reference genes for expression analysis in mouse models of acute alcoholic liver injury

Investigations of hepatic gene expression are crucial for determining the molecular factors involved in acute alcoholic liver injury. The results of liver molecular investigations may reveal etiologically important genomic alterations. Therefore, it is necessary to normalize gene expression data to identify stable genes, which may be used as a reference under different experimental conditions. The aim of the present study was to apply reverse transcription‑quantitative polymerase chain reaction analysis and use analysis software to investigate the expression stability of candidate reference genes in hepatic tissues from mice with acute alcoholic liver injury. The acute alcoholic liver injury models were established by the intragastric administration of alcohol (5 mg/kg) in Imprinting Control Region mice. Total RNA was isolated from the mouse livers, following which the expression levels of seven reference genes, β-actin, glyceraldehyde 3-phosphate dehydrogenase (Gadph), glucuronidase β, hypoxanthine phosphoribosyltransferase 1 (Hprt1), 18S ribosomal RNA, TATA binding protein and β‑2 microglobulin, were examined, and gene expression stability was assessed using the geNorm, NormFinder and BestKeeper tools. The geNorm analysis revealed that the gene with the lowest variability was Hprt1. Hprt1 and Gapdh were validated as the optimal reference gene pair in all samples from all groups. The NormFinder and BestKeeper results showed that Hprt1 was the most stable gene in all samples. Alcohol induces endoplasmic reticulum (ER) stress, causing changes in the expression levels of ER stress‑associated genes. The stability of Hprt1 was verified by the expression analysis of ER stress‑associated genes, and gene expression levels in the ethanol groups were upregulated, with a significant difference in expression, compared with those in the control group. Therefore, Hprt1 was selected as the most stable gene, and Hprt1 and Gapdh were determined to be the optimum gene pair in mouse models of acute alcoholic liver injury. The reliability of the Hprt1 gene was confirmed by expression analysis of ER stress‑associated genes.

Genome-wide identification and validation of new reference genes for transcript normalization in developmental and post-harvested fruits of Actinidia chinensis

The appropriate reference genes are important and essential for reliable results of transcript normalization in real-time qRT-PCR. In the current study, we identified 1203 stably expressed genes from 35,286 genes' expression profiles in developmental fruits of Actinidia chinensis. We manually selected six candidate genes and assessed their expression levels, using two sets of fruit samples of A. chinensis: flesh fruits at four developmental stages and post-harvested fruits. The expression stability of these six genes was assessed by three independent algorithms: geNorm, NormFinder, and BestKeeper. Statistical results indicated these six genes can serve as internal control in both developmental and post-harvested fruits. Among these genes, UBQ_CONJ_E2 (Ubiquitin-conjugating enzyme E2 36) and TUB_FCB (Tubulin folding cofactor B) were the two best reference genes identified in this study. The identification and validation of these reference genes can be helpful for elucidating the studies of fruit development and post-harvested fruits' storage in A. chinensis and other fruit crops of Actinidiaceae.

Use of RNA-seq data to identify and validate RT-qPCR reference genes for studying the tomato-Pseudomonas pathosystem

The agronomical relevant tomato-Pseudomonas syringae pv. tomato pathosystem is widely used to explore and understand the underlying mechanisms of the plant immune response. Transcript abundance estimation, mainly through reverse transcription-quantitative PCR (RT-qPCR), is a common approach employed to investigate the possible role of a candidate gene in certain biological process under study. The accuracy of this technique relies heavily on the selection of adequate reference genes. Initially, genes derived from other techniques (such as Northern blots) were used as reference genes in RT-qPCR experiments, but recent studies in different systems suggest that many of these genes are not stably expressed. The development of high throughput transcriptomic techniques, such as RNA-seq, provides an opportunity for the identification of transcriptionally stable genes that can be adopted as novel and robust reference genes. Here we take advantage of a large set of RNA-seq data originating from tomato leaves infiltrated with different immunity inducers and bacterial strains. We assessed and validated 9 genes that are much more stable than two traditional reference genes. Specifically, ARD2 and VIN3 were the most stably expressed genes and consequently we propose they be adopted for RT-qPCR experiments involving this pathosystem.

Selection of Reference Genes for qRT-PCR Analysis of Gene Expression in Stipa grandis during Environmental Stresses

Stipa grandis P. Smirn. is a dominant plant species in the typical steppe of the Xilingole Plateau of Inner Mongolia. Selection of suitable reference genes for the quantitative real-time reverse transcription polymerase chain reaction (qRT-PCR) is important for gene expression analysis and research into the molecular mechanisms underlying the stress responses of S. grandis. In the present study, 15 candidate reference genes (EF1 beta, ACT, GAPDH, SamDC, CUL4, CAP, SNF2, SKIP1, SKIP5, SKIP11, UBC2, UBC15, UBC17, UCH, and HERC2) were evaluated for their stability as potential reference genes for qRT-PCR under different stresses. Four algorithms were used: GeNorm, NormFinder, BestKeeper, and RefFinder. The results showed that the most stable reference genes were different under different stress conditions: EF1beta and UBC15 during drought and salt stresses; ACT and GAPDH under heat stress; SKIP5 and UBC17 under cold stress; UBC15 and HERC2 under high pH stress; UBC2 and UBC15 under wounding stress; EF1beta and UBC17 under jasmonic acid treatment; UBC15 and CUL4 under abscisic acid treatment; and HERC2 and UBC17 under salicylic acid treatment. EF1beta and HERC2 were the most suitable genes for the global analysis of all samples. Furthermore, six target genes, SgPOD, SgPAL, SgLEA, SgLOX, SgHSP90 and SgPR1, were selected to validate the most and least stable reference genes under different treatments. Our results provide guidelines for reference gene selection for more accurate qRT-PCR quantification and will promote studies of gene expression in S. grandis subjected to environmental stress.

Inducible Expression of the De-Novo Designed Antimicrobial Peptide SP1-1 in Tomato Confers Resistance to Xanthomonas campestris pv. vesicatoria

Antimicrobial peptides (AMPs) are small peptides with less than 50 amino acids and are part of the innate immune response in almost all organisms, including bacteria, vertebrates, invertebrates and plants. AMPs are active against a broad-spectrum of pathogens. The inducible expression of AMPs in plants is a promising approach to combat plant pathogens with minimal negative side effects, such as phytotoxicity or infertility. In this study, inducible expression of the de-novo designed AMP SP1-1 in Micro Tom tomato protected tomato fruits against bacterial spot disease caused by Xanthomonas campestris pv. vesicatoria. The peptide SP1-1 was targeted to the apoplast which is the primary infection site for plant pathogens, by fusing SP1-1 peptide to the signal peptide RsAFP1 of radish (Raphanus sativus). The pathogen inducibility of the expression was enabled by using an optimized inducible 4XW2/4XS promoter. As a result, the tomato fruits of independently generated SP1-1 transgenic lines were significantly more resistant to X. campestris pv. vesicatoria than WT tomato fruits. In transgenic lines, bacterial infection was reduced up to 65% in comparison to the infection of WT plants. Our study demonstrates that the combination of the 4XW2/4XS cis-element from parsley with the synthetic antimicrobial peptide SP1-1 is a good alternative to protect tomato fruits against infections with X. campestris pv. vesicatoria.

Selection of Reference Genes for Quantitative Real-Time RT-PCR Studies in Tomato Fruit of the Genotype MT-Rg1

Quantitative real-time RT-PCR (qRT-PCR) has become one of the most widely used methods for accurate quantification of gene expression. Since there are no universal reference genes for normalization, the optimal strategy to normalize raw qRT-PCR data is to perform an initial comparison of a set of independent reference genes to assess the most stable ones in each biological model. Normalization of a qRT-PCR experiment helps to ensure that the results are both statistically significant and biologically meaningful. Tomato is the model of choice to study fleshy fruit development. The miniature tomato (Solanum lycopersicum L.) cultivar Micro-Tom (MT) is considered a model system for tomato genetics and functional genomics. A new genotype, containing the Rg1 allele, improves tomato in vitro regeneration. In this work, we evaluated the expression stability of four tomato reference genes, namely CAC, SAND, Expressed, and ACTIN2. We showed that the genes CAC and Exp are the best reference genes of the four we tested during fruit development in the MT-Rg1 genotype. Furthermore, we validated the reference genes by showing that the expression profiles of the transcription factors FRUITFULL1 and APETALA2c during fruit development are comparable to previous reports using other tomato cultivars.