Transcriptomic Profiling of Tomato Leaves Identifies Novel Transcription Factors Responding to Dehydration Stress.

Drought is among the most challenging environmental restrictions to tomatoes (Solanum lycopersi-cum), which causes dehydration of the tissues and results in massive loss of yield. Breeding for dehydration-tolerant tomatoes is a pressing issue as a result of global climate change that leads to increased duration and frequency of droughts. However, the key genes involved in dehydration response and tolerance in tomato are not widely known, and genes that can be targeted for dehydration-tolerant tomato breeding remains to be discovered. Here, we compared phenotypes and transcriptomic profiles of tomato leaves between control and dehydration conditions. We show that dehydration decreased the relative water content of tomato leaves after 2 h of dehydration treatment; however, it promoted the malondialdehyde (MDA) content and ion leakage ratio after 4 h and 12 h of dehydration, respectively. Moreover, dehydration stress triggered oxidative stress as we detected significant increases in H2O2 and O2- levels. Simultaneously, dehydration enhanced the activities of antioxidant enzymes including peroxidase (POD), superoxide dismutase (SOD), catalase (CAT), and phenylalanine ammonia-lyase (PAL). Genome-wide RNA sequencing of tomato leaves treated with or without dehydration (control) identified 8116 and 5670 differentially expressed genes (DEGs) after 2 h and 4 h of dehydration, respectively. These DEGs included genes involved in translation, photosynthesis, stress response, and cytoplasmic translation. We then focused specifically on DEGs annotated as transcription factors (TFs). RNA-seq analysis identified 742 TFs as DEGs by comparing samples dehydrated for 2 h with 0 h control, while among all the DEGs detected after 4 h of dehydration, only 499 of them were TFs. Furthermore, we performed real-time quantitative PCR analyses and validated expression patterns of 31 differentially expressed TFs of NAC, AP2/ERF, MYB, bHLH, bZIP, WRKY, and HB families. In addition, the transcriptomic data revealed that expression levels of six drought-responsive marker genes were upregulated by de-hydration treatment. Collectively, our findings not only provide a solid foundation for further functional characterization of dehydration-responsive TFs in tomatoes but may also benefit the improvement of dehydration/drought tolerance in tomatoes in the future.

Differences in Physiological Responses of Two Tomato Genotypes to Combined Waterlogging and Cadmium Stresses.

Waterlogging and heavy mental (e.g., cadmium) stress are two primary threats to crop growth. The combination of abiotic stresses was common and frequent, especially in the field condition. Even though the effects of individual waterlogging and cadmium on tomato plants have been widely investigated, the response of tomatoes under combined waterlogging and cadmium stress remains unclear. This study aimed to clarify and compare physiological, biochemical characteristics and plant growth of two tomato genotypes under individual and combined stress. Two tomato genotypes ('MIX-002' and 'LA4440') were treated under control, waterlogging, cadmium stress and their combination. The results showed that chloroplast ultrastructure of tomatoes under individual and combined stress was damaged with disordered stroma and grana lamellae. The H2O2 (hydrogen peroxide) content and O2·- (superoxide anion radical) production rate of plants under all the three stresses was not significantly higher than the control except for 'LA4440' under the combined stress. Antioxidant enzymes actively responded in the two tomato genotypes, as shown by significant increase in SOD activity from 'MIX-002' under waterlogging and combined stress and from 'LA4440' under cadmium. Meanwhile, CAT activity of 'MIX-002' under waterlogging and 'LA4440' under combined stress significantly decreased, and the POD activity of 'MIX-002' under combined stress significantly increased as compared with the respective control. The APX activity of 'MIX-002' and 'LA4440' under combined stress was significantly lower and higher than the respective controls. This indicated that tomato plants were able to secure redox homeostasis and protect plants from oxidative damage through the synergetic regulation of antioxidant enzymes. Plant height and biomass of the two genotypes under individual and combined stress significantly decreased, which could be a direct result from the chloroplast alteration and resource re-allocation. Overall, the effects of combined waterlogging and cadmium stress were not simply the sum of individual effects on two tomato genotypes. Distinct ROS (reactive oxygen species) scavenging systems of two tomato genotypes under stresses suggest a genotype-dependent antioxidant enzymes regulation.

Low light intensity elongates period and defers peak time of photosynthesis: a computational approach to circadian-clock-controlled photosynthesis in tomato.

Photosynthesis is involved in the essential process of transforming light energy into chemical energy. Although the interaction between photosynthesis and the circadian clock has been confirmed, the mechanism of how light intensity affects photosynthesis through the circadian clock remains unclear. Here, we propose a first computational model for circadian-clock-controlled photosynthesis, which consists of the light-sensitive protein P, the core oscillator, photosynthetic genes, and parameters involved in the process of photosynthesis. The model parameters were determined by minimizing the cost function ( [Formula: see text]), which is defined by the errors of expression levels, periods, and phases of the clock genes (CCA1, PRR9, TOC1, ELF4, GI, and RVE8). The model recapitulates the expression pattern of the core oscillator under moderate light intensity (100 µmol m -2 s-1). Further simulation validated the dynamic behaviors of the circadian clock and photosynthetic outputs under low (62.5 µmol m-2 s-1) and normal (187.5 µmol m-2 s-1) intensities. When exposed to low light intensity, the peak times of clock and photosynthetic genes were shifted backward by 1-2 hours, the period was elongated by approximately the same length, and the photosynthetic parameters attained low values and showed delayed peak times, which confirmed our model predictions. Our study reveals a potential mechanism underlying the circadian regulation of photosynthesis by the clock under different light intensities in tomato.

FIS1 encodes a GA2-oxidase that regulates fruit firmness in tomato.

Fruit firmness is a target trait in tomato breeding because it facilitates transportation and storage. However, it is also a complex trait and uncovering the molecular genetic mechanisms controlling fruit firmness has proven challenging. Here, we report the map-based cloning and functional characterization of qFIRM SKIN 1 (qFIS1), a major quantitative trait locus that partially determines the difference in compression resistance between cultivated and wild tomato accessions. FIS1 encodes a GA2-oxidase, and its mutation leads to increased bioactive gibberellin content, enhanced cutin and wax biosynthesis, and increased fruit firmness and shelf life. Importantly, FIS1 has no unfavorable effect on fruit weight or taste, making it an ideal target for breeders. Our study demonstrates that FIS1 mediates gibberellin catabolism and regulates fruit firmness, and it offers a potential strategy for tomato breeders to produce firmer fruit.

Unique miRNAs and their targets in tomato leaf responding to combined drought and heat stress.

BACKGROUND: Both drought and heat stress are serious global problems, leading to agricultural production loss. MicroRNAs (miRNAs) play important roles in plant species responding to individual drought and heat stress. However, the miRNAs and mRNAs in association with combined drought and heat in crops like tomato remains unclear. RESULTS: We studied the crosstalk of miRNAs and their target genes in tomato plants grown under simultaneous drought and heat stress that frequently happen in field conditions. In total, 335 known miRNAs representing 55 miRNA families and 430 potential novel miRNAs were identified in Solanum lycopersicum L. using small RNA deep sequencing. Through expression analysis, miRNAs in association with drought, heat and the combination of these were investigated. In total, 61, 74 and 37 miRNAs were differentially regulated for combination (of both stresses) vs control, combination vs drought and combination vs heat, respectively. Target genes with different expression levels were found using degradome sequencing, which were mainly involved in transcription factor activity, sequence-specific DNA binding, transcription, regulation of transcription, nucleus, DNA binding etc. The quantitative real-time polymerase chain reaction (qRT-PCR) results confirmed the accuracy of sequencing. CONCLUSIONS: Our study serves as valuable knowledge on how crop adapted to combined drought and heat stress by regulating miRNAs and mRNAs, which provide information for crop improvement to deal with future climate changes.

Physiological analysis and transcriptome sequencing reveal the effects of combined cold and drought on tomato leaf.

BACKGROUND: Co-occurrence of cold and drought stress can alter the response of plants at morphological, physiological and molecular levels, which finally affect crop production, more than individual stress. Understanding the responses of crop to combined stress is necessary to improve tolerance and maintain crop production especially in the field where combined stress frequently occurs. We aimed to clarify the underlying leaf physiological and molecular mechanisms of tomato by imposing combining cold and drought on one popular tomato cultivar 'Jinlingmeiyu' as an example. RESULTS: The physiological and genetic responses were identified in tomatoes after 42 h exposure to control, cold, drought and combined treatments. As compared with control, water loss rate at the three stresses including cold, drought and combined stress significantly decreased until 40 min after taking samples from the plants. The content of H2O2, zeatin riboside (ZR) and melatonin in all stress treatments were significantly higher than the control. Drought stress alone and combined stress induced the accumulation of abscisic acid (ABA) and auxin (IAA) as compared with control. The individual cold and combined stress significantly decreased the maximum quantum efficiency of PSII (Fv/Fm), quantum yield of PSII (Fq'/Fm') and electron transport rate (ETR). In total, 7141, 1850 and 7841 genes were involved in the stress response to cold, drought and their combination. Functional analysis of the stress-inducible genes provided more insights concerning the complex regulatory mechanisms that were involved in combined stress. The expression level of 12 genes were validated by quantitative real-time PCR (qRT-PCR). CONCLUSIONS: We found that the expression of stress-specific genes changed with physiological variation, indicating the close crosstalk between physiological and genetic response especially under combined stress. This study provides new knowledge on the complex regulatory mechanism genes in tomato ('Jinlingmeiyu') leaf to abiotic stresses.

Physiological response of tomatoes at drought, heat and their combination followed by recovery.

In nature, crops encounter a combination of abiotic stresses that severely limit yield. Our aim was to dynamically expose the changes of tomatoes' physiological parameters to drought, heat and their combination and thereby clarify the relationship between the responses to single and combined stress. We studied the effect of single and combined drought and heat stresses on the shoot and root of two tomato cultivars (Sufen No.14 as CV1; Jinlingmeiyu as CV2). After being exposed to combined stress for 6 days, the dry weight of shoot and root significantly decreased. The Fq '/Fm ' (quantum yield of photosystem II) was significantly lower in CV1 upon drought and combined stress and in CV2 subjected to combined stress (between days 4 and 6) compared to control. The relative water content during combined stress was significantly lower than control from day 4 to recovery day 2. On days 3 and 6, the water loss rate significantly increased under heat stress and decreased at drought and combined stress, respectively. The combined stress caused severe damages on photosystem II and chloroplast ultrastructure. The root activity after stress recovered even though drought significantly increased the activity from day 2 to day 6. Combined stress result in complex responses during tomato growth. The CV1 was more heat tolerant than CV2, but there was no varietal difference at drought and combined stress. This study contributes to the understanding of the underlying physiological response mechanism of plant to combined stress and crop improvement by providing valuable information for abiotic stress-tolerant tomato breeding.

Application of Whole Genome Resequencing in Mapping of a Tomato Yellow Leaf Curl Virus Resistance Gene.

Tomato yellow leaf curl virus (TYLCV) has significantly impacted the tomato industry around the world, and the use of insecticides and insect nets have not effectively controlled the spread of this pathogen. The tomato line AVTO1227 is highly resistant to TYLCV. In this study, F2 and BC1 populations derived from AVTO1227 and the susceptible line Money maker were used to assess the genetic mechanism underlying TYLCV resistance. We have identified a recessive TYLCV resistance gene, hereby designated as ty-5, which is linked to SlNACI. Genomic DNA pools from resistant and susceptible groups were constructed, and their genomes were resequenced. The ty-5 gene was identified on an interval encompassing the genomic positions 2.22 Mb to 3.19 Mb on tomato chromosome 4. Genotyping using linkage markers further mapped ty-5 within the interval between markers ty5-25 and ty5-29, where only the pelota gene is located. Consequently, pelota was considered as the candidate gene corresponding to ty-5. Two nucleotide transversions within the promoter region and one transversion in exon region of the pelota gene were detected in the parental lines. However, the relative transcript levels of pelota did not significantly differ among the three tomato lines, regardless of TYLCV infection. This study will facilitate marker-assisted breeding for resistance to TYLCV and lay a foundation for the research of the resistance mechanism of ty-5 in tomato.

Genome-wide identification of circRNAs involved in tomato fruit coloration.

Circular RNAs (circRNAs) are an emerging class of non-coding RNAs in plants. Our aim is to identify the circRNAs with different expression levels between tomato fruits from two cultivars ('Jinling Fenyu' and 'Jinling Moyu') at two different developmental stages. Following high-throughput sequencing, 3796 circRNAs were identified, and 243 circRNAs were shared in the four samples. As compared with the fruit at mature green stage, the expression levels of 273 and 89 circRNAs were significantly altered in the fruit at turning stage from 'Jinling Fenyu' and 'Jinling Moyu', respectively. Moreover, the parental genes of the circRNAs with significantly different expression level were mainly involved in metabolic, cellular and single-organism process and played roles in catalytic activity and binding based on GO (Gene Ontology) analysis. The results suggested that circRNAs were widespread in tomato and were generated from different chromosomes and diverse genomic regions. Some circRNAs were specific in tomato fruits at different developmental stages, which enriches the number of circRNAs in plants involved in fruit coloration and ripening. This study provides the first genome-wide profile of circRNAs involved in tomato fruit coloration and lays a foundation for studying the potential biological functions of circRNAs involved in fruit ripening.

Drought stress had a predominant effect over heat stress on three tomato cultivars subjected to combined stress.

BACKGROUND: Abiotic stresses due to environmental factors could adversely affect the growth and development of crops. Among the abiotic stresses, drought and heat stress are two critical threats to crop growth and sustainable agriculture worldwide. Considering global climate change, incidence of combined drought and heat stress is likely to increase. The aim of this study was to shed light on plant growth performance and leaf physiology of three tomatoes cultivars ('Arvento', 'LA1994' and 'LA2093') under control, drought, heat and combined stress. RESULTS: Shoot fresh and dry weight, leaf area and relative water content of all cultivars significantly decreased under drought and combined stress as compared to control. The net photosynthesis and starch content were significantly lower under drought and combined stress than control in the three cultivars. Stomata and pore length of the three cultivars significantly decreased under drought and combined stress as compared to control. The tomato 'Arvento' was more affected by heat stress than 'LA1994' and 'LA2093' due to significant decreases in shoot dry weight, chlorophyll a and carotenoid content, starch content and NPQ (non-photochemical quenching) only in 'Arvento' under heat treatment. By comparison, the two heat-tolerant tomatoes were more affected by drought stress compared to 'Arvento' as shown by small stomatal and pore area, decreased sucrose content, ΦPSII (quantum yield of photosystem II), ETR (electron transport rate) and qL (fraction of open PSII centers) in 'LA1994' and 'LA2093'. The three cultivars showed similar response when subjected to the combination of drought and heat stress as shown by most physiological parameters, even though only 'LA1994' and 'LA2093' showed decreased Fv/Fm (maximum potential quantum efficiency of photosystem II), ΦPSII, ETR and qL under combined stress. CONCLUSIONS: The cultivars differing in heat sensitivity did not show difference in the combined stress sensitivity, indicating that selection for tomatoes with combined stress tolerance might not be correlated with the single stress tolerance. In this study, drought stress had a predominant effect on tomato over heat stress, which explained why simultaneous application of heat and drought revealed similar physiological responses to the drought stress. These results will uncover the difference and linkage between the physiological response of tomatoes to drought, heat and combined stress and be important for the selection and breeding of tolerant tomato cultivars under single and combine stress.

Genome-wide analysis of bHLH transcription factor and involvement in the infection by yellow leaf curl virus in tomato (Solanum lycopersicum).

BACKGROUND: The basic helix-loop-helix (bHLH) proteins are a superfamily of transcription factors that can bind to specific DNA target sites. They have been well characterized in model plants such as Arabidopsis and rice and have been shown to be important regulatory components in many different biological processes. However, no systemic analysis of the bHLH transcription factor family has yet been reported in tomatoes. Tomato yellow leaf curl virus (TYLCV) threatens tomato production worldwide by causing leaf yellowing, leaf curling, plant stunting and flower abscission. RESULTS: A total of 152 bHLH transcription factors were identified from the entire tomato genome. Phylogenetic analysis of bHLH domain sequences from Arabidopsis and tomato facilitated classification of these genes into 26 subfamilies. The evolutionary and possible functional relationships revealed during this analysis are supported by other criteria, including the chromosomal distribution of these genes, the conservation of motifs and exon/intron structural patterns, and the predicted DNA binding activities within subfamilies. Distribution mapping results showed bHLH genes were localized on the 12 tomato chromosomes. Among the 152 bHLH genes from the tomato genome, 96 bHLH genes were detected in the TYLCV-susceptible and resistant tomato breeding line before (0 dpi) and after TYLCV (357 dpi) infection. As anticipated, gene ontology (GO) analysis indicated that most bHLH genes are related to the regulation of macromolecule metabolic processes and gene expression. Only four bHLH genes were differentially expressed between 0 and 357 dpi. Virus-induced gene silencing (VIGS) of one bHLH genes SlybHLH131 in resistant lines can lead to the cell death. CONCLUSION: In the present study, 152 bHLH transcription factor genes were identified. One of which bHLH genes, SlybHLH131, was found to be involved in the TYLCV infection through qRT-PCR expression analysis and VIGS validation. The isolation and identification of these bHLH transcription factors facilitated clarification of the molecular genetic basis for the genetic improvement of tomatoes and the development of functional gene resources for transgenic research. In addition, these findings may aid in uncovering an unexplored mechanism during the TYLCV infection in tomatoes.

Functional conservation of the glycosyltransferase gene GT47A in the monocot rice.

Glucuronoarabinoxylan is the major hemicellulose in grass cell walls, yet the mechanism of xylan synthesis in monocot plants is still unclear. Unraveling the genes involved in the biosynthesis of xylan in rice will be very important for the utilization of rice straw as a source of bioenergy in the future. In this report, we investigated the functional role of a rice gene homologous to Arabidopsis IRREGULAR XYLEM10 (IRX10), belonging to the glycosyl transferase (GT) gene family 47 (GT47), in the biosynthesis of xylan. The protein sequence of OsGT47A from rice exhibits a 93.49% similarity to IRX10, which is involved in the biosynthesis of glucuronoxylan in Arabidopsis. Phylogenetic analysis of the GT47 glycosyl transferase family in the rice genome revealed that OsGT47A is a closely related homolog of IRX10 and IRX10L. Expression pattern analysis showed that the OsGT47A gene is highly expressed in the rice stem. Overexpression of OsGT47A in the irx10 irx10L double mutant rescued the plant growth phenotype and restored secondary wall thickness. Analysis of monosaccharides indicated that the rescued plants had levels of xylose identical to those of the wild type plants, and the fluorescence signals were restored in the complementation plants by xylan immunolocalization. The OsGT47A complementation under the native promoter of Arabidopsis IRX10L (ProIRX10L) partially rescued the double mutant, indicating that OsGT47A is functionally equivalent to IRX10L. Together, these results suggest that the IRX10 homolog OsGT47A exhibits functional conservation and is most likely involved in xylan synthesis in rice.