Summary: Drought and Verticillium wilt disease are two main factors that limit cotton production, which necessitates the identification of key molecular switch to simultaneously improve cotton resistance to Verticillium dahliae and tolerance to drought stress. R2R3‐type MYB proteins could play such a role because of their conserved functions in plant development, growth, and metabolism regulation, however, till date a MYB gene conferring the desired resistance to both biotic and abiotic stresses has not been found in cotton. Here, we describe the identification of GhMYB36, a gene encoding a R2R3‐type MYB protein in Gossypium hirsutum, which confers drought tolerance and Verticilium wilt resistance in both Arabidopsis and cotton. GhMYB36 was highly induced by PEG‐simulated drought stress in G. hirsutum. GhMYB36‐silenced cotton plants were more sensitive to both drought stress and Verticillium wilt. GhMYB36 overexpression in transgenic Arabidopsis and cotton plants gave rise to improved drought tolerance and Verticillium wilt resistance. Transient expression of fused GhMYB36‐GFP in tobacco cells was able to localize GhMYB36 in the cell nucleus. In addition, RNA‐seq analysis together with qRT‐PCR validation in transgenic Arabidopsis overexpressing GhMYB36 revealed significantly enhanced PR1 expression. Luciferase interaction assays indicated that GhMYB36 are probably bound to the promoter of PR1 to activate its expression and the interaction, which was further verified by Yeast one hybrid assay. Taken together, our results suggest that GhMYB36 functions as a transcription factor that is involved in drought tolerance and Verticillium wilt resistance in Arabidopsis and cotton by enhancing PR1 expression.

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. 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. 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.

Long Noncoding-RNAs (LncRNAs) are known to be involved in some biological processes, but their roles in plant-virus interactions remain largely unexplored. While circular RNAs (circRNAs) have been studied in animals, there has yet to be extensive research on them in a plant system, especially in tomato-tomato yellow leaf curl virus (TYLCV) interaction. In this study, RNA transcripts from the susceptible tomato line JS-CT-9210 either infected with TYLCV or untreated, were sequenced in a pair-end strand-specific manner using ribo-zero rRNA removal library method. A total of 2056 lncRNAs including 1767 long intergenic non-coding RNA (lincRNAs) and 289 long non-coding natural antisense transcripts (lncNATs) were obtained. The expression patterns in lncRNAs were similar in susceptible tomato plants between control check (CK) and TYLCV infected samples. Our analysis suggested that lncRNAs likely played a role in a variety of functions, including plant hormone signaling, protein processing in the endoplasmic reticulum, RNA transport, ribosome function, photosynthesis, glulathione metabolism, and plant-pathogen interactions. Using virus-induced gene silencing (VIGS) analysis, we found that reduced expression of the lncRNA S-slylnc0957 resulted in enhanced resistance to TYLCV in susceptible tomato plants. Moreover, we identified 184 circRNAs candidates using the CircRNA Identifier (CIRI) software, of which 32 circRNAs were specifically expressed in untreated samples and 83 circRNAs in TYLCV samples. Approximately 62% of these circRNAs were derived from exons. We validated the circRNAs by both PCR and Sanger sequencing using divergent primers, and found that most of circRNAs were derived from the exons of protein coding genes. The silencing of these circRNAs parent genes resulted in decreased TYLCV virus accumulation. In this study, we identified novel lncRNAs and circRNAs using bioinformatic approaches and showed that these RNAs function as negative regulators of TYLCV infection. Moreover, the expression patterns of lncRNAs in susceptible tomato plants were different from that of resistant tomato plants, while exonic circRNAs expression positively associated with their respective protein coding genes. This work provides a foundation for elaborating the novel roles of lncRNAs and circRNAs in susceptible tomatoes following TYLCV infection.

Cold stress is known to influence tomato growth, development, and yield. In this study, we analyzed the germination of tomato seeds treated with exogenous glycine betaine (GB) at a low temperature (14 °C). The results showed that cold stress inhibited tomato seed germination, and pretreatment with exogenous GB reduced this inhibition and enhanced the germination rate (GR), germination index (GI), and viability of tomato seeds at low temperatures. Analysis of gene expression and metabolism revealed that GB positively regulated endogenous hormone gibberellin (GA) content and negatively regulated abscisic acid (ABA) content, while GB reduced the starch content in the seeds by up-regulating the amylase gene expression. Gene expression analysis showed that the key genes (SlSOD, SlPOD, and SlchlAPX) involved in reactive oxygen species (ROS) scavenging systems were up-regulated in GB-pretreated tomato seeds compared with the control. At the same time, levels of malondialdehyde and hydrogen peroxide were significantly lower, while the proline content and peroxidase (POD), superoxide dismutase (SOD), and catalase (CAT) levels were elevated compared with those in the control. These results demonstrate that exogenous GB as a positive regulator effectively alleviated the inhibition of tomato seed germination under cold stress by different signal pathways.

The basic leucine zipper (bZIP) is a transcription factor family that plays critical roles in abiotic and biotic stress responses as well as plant development and growth. A comprehensive genome-wide study in Liriodendron chinense was conducted to identify 45 bZIP transcription factors (LchibZIPs), which were divided into 13 subgroups according the phylogenetic analysis. Proteins in the same subgroup shared similar gene structures and conserved domains, and a total of 20 conserved motifs were revealed in LchibZIP proteins. Gene localization analysis revealed that LchibZIP genes were unequally distributed across 16 chromosomes, and that 4 pairs of tandem and 9 segmental gene duplications existed. Concluding that segmental duplication events may be strongly associated with the amplification of the L. chinense bZIP gene family. We also assessed the collinearity of LchibZIPs between the Arabidopsis and Oryza and showed that the LchibZIP is evolutionarily closer to O. sativa as compared to the A. thaliana. The cis-regulatory element analysis showed that LchibZIPs clustered in one subfamily are involved in several functions. In addition, we gathered novel research suggestions for further exploration of the new roles of LchibZIPs from protein-protein interactions and gene ontology annotations of the LchibZIP proteins. Using the RNA-seq data and qRT-PCR we analyzed the gene expression patterns of LchibZIP genes, and showed that LchibZIP genes regulate cold stress, especially LchibZIP4 and LchibZIP7; and LchibZIP2 and LchibZIP28 which were up-regulated and down-regulated by cold stress, respectively. Studies of genetic engineering and gene function in L. chinense can benefit greatly from the thorough investigation and characterization of the L. chinense bZIP gene family.