Solanum lycopersicum and Solanum habrochaites are closely related plant species; however, their cold tolerance capacities are different. The wild species S. habrochaites is more cold tolerant than the cultivated species S. lycopersicum. The transcriptomes of S. lycopersicum and S. habrochaites leaf tissues under cold stress were studied using Illumina high-throughput RNA sequencing. The results showed that more than 200 million reads could be mapped to identify genes, microRNAs (miRNAs), and alternative splicing (AS) events to confirm the transcript abundance under cold stress. The results indicated that 21 % and 23 % of genes were differentially expressed in the cultivated and wild tomato species, respectively, and a series of changes in S. lycopersicum and S. habrochaites transcriptomes occur when plants are moved from warm to cold conditions. Moreover, the gene expression patterns for S. lycopersicum and S. habrochaites were dissimilar; however, there were some overlapping genes that were regulated by low temperature in both tomato species. An AS analysis identified 75,885 novel splice junctions among 172,910 total splice junctions, which suggested that the relative abundance of alternative intron isoforms in S. lycopersicum and S. habrochaites shifted significantly under cold stress. In addition, we identified 89 miRNA sequences that may regulate relevant target genes. Our data indicated that some miRNAs (e.g., miR159, miR319, and miR6022) play roles in the response to cold stress. Differences in gene expression, AS events, and miRNAs under cold stress may contribute to the observed differences in cold tolerance of these two tomato species.

Background: APETALA2/ethylene responsive factor (AP2/ERF) transcription factors are a plant-specific family of transcription factors and one of the largest families of transcription factors. Ethylene response factors (ERF) regulate plant growth, development, and responses to biotic and abiotic stress. In a previous study, the ERF2 gene was significantly upregulated in both resistant and susceptible tomato cultivars in response to Stemphylium lycopersici. The main purpose of this study was to systematically analyze the ERF family and to explore the mechanism of ERF2 in tomato plants resisting pathogen infection by the Virus-induced Gene Silencing technique. Results: In this experiment, 134 ERF genes were explored and subjected to bioinformatic analysis and divided into twelve groups. The spatiotemporal expression characteristics of ERF transcription factor gene family in tomato were diverse. Combined with RNA-seq, we found that the expression of 18 ERF transcription factors increased after inoculation with S. lycopersici. In ERF2-silenced plants, the susceptible phenotype was observed after inoculation with S. lycopersici. The hypersensitive response and ROS production were decreased in the ERF2-silenced plants. Physiological analyses showed that the superoxide dismutase, peroxidase and catalase activities were lower in ERF2-silenced plants than in control plants, and the SA and JA contents were lower in ERF2-silenced plants than in control plants after inoculation with S. lycopersici. Furthermore, the results indicated that ERF2 may directly or indirectly regulate Pto, PR1b1 and PR-P2 expression and enhance tomato resistance. Conclusions: In this study, we identified and analyzed members of the tomato ERF family by bioinformatics methods and classified, described and analyzed these genes. Subsequently, we used VIGS technology to significantly reduce the expression of ERF2 in tomatoes. The results showed that ERF2 had a positive effect on tomato resistance to S. lycopersici. Interestingly, ERF2 played a key role in multiple SA, JA and ROS signaling pathways to confer resistance to invasion by S. lycopersici. In addition, ERF2 may directly or indirectly regulate Pto, PR1b1 and PR-P2 expression and enhance tomato resistance to S. lycopersici. In summary, this study provides gene resources for breeding for disease resistance in tomato.

Heterosis and combining ability are two important considerations in the utilization of heterosis, which can be used to generate excellent hybrid resource candidates and is very important in conventional hybrid breeding. In this study, the combining ability, heritability, and heterosis of eight major agronomic traits were analyzed in 10 tomato parents and 45 crosses between them. As well as TY-301, a recognized and official excellent variety that is currently selling well on the market was used as a control to conduct a control heterosis analysis, with the goal of selecting ideal parents with high combining ability and new hybrids with commodity value, high yield, early maturity, and high quality. The results showed that both additive and nonadditive genetic effects are involved in the expression of the traits and that the additive genetic effect is dominant in trait inheritance. Although general combining ability (GCA) and specific combining ability (SCA) were not correlated, and the strength of heterosis depends on SCA, the sum of the parental GCA values (GCAsum) did predict heterosis for some traits with higher predictive accuracy than did SCA. Compared with heterosis, GCAsum can better predict hybrid performance. Finally, the parent 17,969 was the breeding material with the best comprehensive trait performance, especially in yield. We screened a high-yielding candidate combination 17,927 × 17,969 and a precocious and good taste candidate combination 17,666 × 17,927. This information may play an important role in the selection of superior parents and hybrid combinations based on combining ability and heterosis analysis.

Cf-12 is an effective gene for resisting tomato leaf mold disease caused by Cladosporium fulvum (C. fulvum). Unlike many other Cf genes such as Cf-2, Cf-4, Cf-5 and Cf-9, no physiological races of C. fulvum that are virulent to Cf-12 carrying plant lines have been identified. In order to better understand the molecular mechanism of Cf-12 gene resistance response, RNA-Seq was used to analyze the transcriptome changes at three different stages of C. fulvum infection (0, 4 and 8 days post infection [dpi]). A total of 9,100 differentially expressed genes (DEGs) between 4 and 0 dpi, 8,643 DEGs between 8 and 0 dpi and 2,547 DEGs between 8 and 4 dpi were identified. In addition, we found that 736 DEGs shared among the above three groups, suggesting the presence of a common core of DEGs in response to C. fulvum infection. These DEGs were significantly enriched in defense-signaling pathways such as the calcium dependent protein kinases pathway and the jasmonic acid signaling pathway. Additionally, we found that many transcription factor genes were among the DEGs, indicating that transcription factors play an important role in C. fulvum defense response. Our study provides new insight on the molecular mechanism of Cf resistance to C. fulvum, especially the unique features of Cf-12 in responding to C. fulvum infection.

Golden 2-Like (G2-like) transcription factors play an important role in plant development. However, the roles of these G2-like regulatory genes in response to abiotic stresses in tomato are not well understood. In this study, we identified 66 putative G2-like genes in tomato (Solanum lycopersicum) and classified them into 5 groups (I to V) according to gene structure, motif composition and phylogenetic analysis. The G2-like genes were unevenly distributed across all 12 chromosomes. There were nine pairs of duplicated gene segments and four tandem duplicated SlGlk genes. Analysis of the cis-regulatory elements (CREs) showed that the promoter regions of SlGlks contain many kinds of stress- and hormone-related CREs. Based on RNA-seq, SlGlks were expressed in response to three abiotic stresses. Thirty-six differentially expressed SlGlks were identified; these genes have multiple functions according to Gene Ontology (GO) analysis and are enriched mainly in the zeatin biosynthesis pathway. Further studies exhibited that silencing SlGlk16 in tomato would reduce drought stress tolerance by earlier wilted, lower superoxide dismutase (SOD), peroxidase (POD) activities, less Pro contents and more MDA contents. Overall, the results of this study provide comprehensive information on G2-like transcription factors and G2-like genes that may be expressed in response to abiotic stresses.