Comparative transcriptome analysis of tomato (Solanum lycopersicum) in response to exogenous abscisic acid

Characterization of tomato autophagy-related SlCOST family genes

Autophagy is a eukaryote-specific cellular process that can engulf unwanted targets with double-membrane autophagosomes and subject them to the vacuole or lysosome for breaking down and recycling, playing dual roles in plant growth and environmental adaptions. However, perception of specific environmental signals for autophagy induction is largely unknown, limiting its application in agricultural usage. Identification of plant-unique DUF641 family COST1 (Constitutively Stressed 1) protein directly links drought perception and autophagy induction, shedding light on manipulating autophagy for breeding stress tolerant crops. In this study, we performed a genome-wide analysis of DUF641/COST family in tomato, and identified five SlCOST genes SlCOST1, -2, -3, -4, and -5. SlCOST genes show both overlapping and distinct expression patterns in plant growth and stress responding. In addition, SlCOST1, -3, -4, -5 proteins demonstrate co-localization with autophagy adaptor protein ATG8e, and all five SlCOST proteins show interactions ATG8e in planta. However, only SlCOST1, the closest ortholog of Arabidopsis AtCOST1, can restore cost1 mutant to WT level, suggesting conserved role of COST1 and functional diversification of SlCOST family in tomato. Our study provides clues for future investigation of autophagy-related COST family and its promising implementations in breeding crops with robust environmental plasticity.

Tetraclinis articulata essential oil emulsion use as alternative to chemical fungicide to control tomato grey mould disease

Tetraclinis articulata essential oil proved to be effective in controlling tomato grey mould, so we would investigate its effect on some tomato defense mechanisms. The pretreatment of Botrytis cinerea infected tomato plants with TAEO emulsion enhanced the activity of antioxidant enzymes activity SOD, CAT, APX, and GPX, and total polyphenols content and it decreased IC50 of free radical-scavenging activity and H2O2 content. Results showed amelioration in antioxidant status in TAEO emulsion treated and B. cinerea infected plants indicating that treatment decreased infection in tomato plants. The qRT-PCR analysis of defense genes expression Chitinase SlChi, transcription factors SlWRKY and SlAP2/ERF, Lipoxygenase SlLOX, and Thioredoxin SlTRX showed that they were up-regulated as early as 12 hpi sustained with a second increase at 48 hpi in TAEO emulsion pretreated and infected plants. These results suggest the potential use of TAEO emulsion as natural product to induce tomato antioxidant status and activate defense genes.

Exogenous abscisic acid treatment regulates protein secretion in sorghum cell suspension cultures

Drought stress adversely affects plant growth, often leading to total crop failure. Upon sensing soil water deficits, plants switch on biosynthesis of abscisic acid (ABA), a stress hormone for drought adaptation. Here, we used exogenous ABA application to dark-grown sorghum cell suspension cultures as an experimental system to understand how a drought-tolerant crop responds to ABA. We evaluated intracellular and secreted proteins using isobaric tags for relative and absolute quantification. While the abundance of only ~ 7% (46 proteins) intracellular proteins changed in response to ABA, ~32% (82 proteins) of secreted proteins identified in this study were ABA responsive. This shows that the extracellular matrix is disproportionately targeted and suggests it plays a vital role in sorghum adaptation to drought. Extracellular proteins responsive to ABA were predominantly defense/detoxification and cell wall-modifying enzymes. We confirmed that sorghum plants exposed to drought stress activate genes encoding the same proteins identified in the in vitro cell culture system with ABA. Our results suggest that ABA activates defense and cell wall remodeling systems during stress response. This could underpin the success of sorghum adaptation to drought stress.

Uncovering the Role of Hormones in Enhancing Antioxidant Defense Systems in Stressed Tomato (Solanum lycopersicum) Plants

Tomato is one of the most important fruits worldwide. It is widely consumed due to its sensory and nutritional attributes. However, like many other industrial crops, it is affected by biotic and abiotic stress factors, reducing its metabolic and physiological processes. Tomato plants possess different mechanisms of stress responses in which hormones have a pivotal role. They are responsible for a complex signaling network, where the antioxidant system (enzymatic and non-enzymatic antioxidants) is crucial for avoiding the excessive damage caused by stress factors. In this sense, it seems that hormones such as ethylene, auxins, brassinosteroids, and salicylic, jasmonic, abscisic, and gibberellic acids, play important roles in increasing antioxidant system and reducing oxidative damage caused by different stressors. Although several studies have been conducted on the stress factors, hormones, and primary metabolites of tomato plants, the effect of endogenous and/or exogenous hormones on the secondary metabolism is still poorly studied, which is paramount for tomato growing management and secondary metabolites production. Thus, this review offers an updated overview of both endogenous biosynthesis and exogenous hormone application in the antioxidant system of tomato plants as a response to biotic and abiotic stress factors.

Abscisic Acid Is Involved in the Resistance Response of Arabidopsis thaliana Against Meloidogyne paranaensis

Abscisic acid (ABA) is a classical hormone involved in the plant defense against abiotic stresses, especially drought. However, its role in the defense response against biotic stresses is controversial: it can induce resistance to some pathogens but can also increase the susceptibility to other pathogens. Information regarding the effect of ABA on the relationship between plants and sedentary phytonematodes, such as Meloidogyne paranaensis, is scarce. In this study, we found that ABA changed the susceptibility level of Arabidopsis thaliana against M. paranaensis. The population of M. paranaensis was reduced by 58.3% with the exogenous application of ABA 24 h before the nematode inoculation, which demonstrated that ABA plays an important role in the preinfectional defense of A. thaliana against M. paranaensis. The increase in the nematode population density in the ABA biosynthesis mutant, aba2-1, corroborated the results observed with the exogenous application of ABA. The phytohormone did not show nematicide or nematostatic effects on M. paranaensis juveniles in in vitro tests, indicating that the response is linked to intrinsic plant factors, which was corroborated by the decrease in the number of nematodes in the abi4-1 mutant. This reduction indicates that the gene expression regulation by transcript factors is possibly related to regulatory cascades mediated by ABA in the response of A. thaliana against M. paranaensis.

Tomato Biodiversity and Drought Tolerance: A Multilevel Review

Ongoing global climate change suggests that crops will be exposed to environmental stresses that may affect their productivity, leading to possible global food shortages. Among these stresses, drought is the most important contributor to yield loss in global agriculture. Drought stress negatively affects various physiological, genetic, biochemical, and morphological characteristics of plants. Drought also causes pollen sterility and affects flower development, resulting in reduced seed production and fruit quality. Tomato (Solanum lycopersicum L.) is one of the most economically important crops in different parts of the world, including the Mediterranean region, and it is known that drought limits crop productivity, with economic consequences. Many different tomato cultivars are currently cultivated, and they differ in terms of genetic, biochemical, and physiological traits; as such, they represent a reservoir of potential candidates for coping with drought stress. This review aims to summarize the contribution of specific physio-molecular traits to drought tolerance and how they vary among tomato cultivars. At the genetic and proteomic level, genes encoding osmotins, dehydrins, aquaporins, and MAP kinases seem to improve the drought tolerance of tomato varieties. Genes encoding ROS-scavenging enzymes and chaperone proteins are also critical. In addition, proteins involved in sucrose and CO₂ metabolism may increase tolerance. At the physiological level, plants improve drought tolerance by adjusting photosynthesis, modulating ABA, and pigment levels, and altering sugar metabolism. As a result, we underline that drought tolerance depends on the interaction of several mechanisms operating at different levels. Therefore, the selection of drought-tolerant cultivars must consider all these characteristics. In addition, we underline that cultivars may exhibit distinct, albeit overlapping, multilevel responses that allow differentiation of individual cultivars. Consequently, this review highlights the importance of tomato biodiversity for an efficient response to drought and for preserving fruit quality levels.

Transcriptome analysis reveals the regulation of cyclic nucleotide-gated ion channels in response to exogenous abscisic acid and calcium treatment under drought stress in tomato

Background: Drought stress can limit the growth and development of tomato seedlings and cause considerable loss of tomato yield. Exogenous abscisic acid (ABA) and calcium (Ca²⁺) can effectively alleviate the damage of drought stress to plants in part because Ca²⁺ acts as a second messenger in the drought resistance pathway. Although cyclic nucleotide-gated ion channels (CNGCs) are common non-specific Ca²⁺ osmotic channels in cell membranes, a thorough understanding of the transcriptome characteristics of tomato treated with exogenous ABA and Ca²⁺ under drought stress is necessary to characterize the molecular mechanism of CNGC involved in tomato drought resistance. Results: There were 12,896 differentially expressed genes in tomato under drought stress, as well as 11,406 and 12,502 differentially expressed genes after exogenous ABA and Ca²⁺ application, respectively. According to functional annotations and reports, the 19 SlCNGC genes related to Ca²⁺ transport were initially screened, with 11 SlCNGC genes that were upregulated under drought stress and downregulated after exogenous ABA application. After exogenous Ca²⁺ application, the data showed that two of these genes were upregulated, while nine genes were downregulated. Based on these expression patterns, we predicted the role of SlCNGC genes in the drought resistance pathway and their regulation by exogenous ABA and Ca²⁺ in tomato. Conclusion: The results of this study provide foundational data for further study of the function of SlCNGC genes and a more comprehensive understanding of drought resistance mechanisms in tomato.

Comprehensive In Silico Analysis of RNA Silencing-Related Genes and Their Regulatory Elements in Wheat (Triticum aestivum L.)

Dicer-like (DCL), Argonaute (AGO), and RNA-dependent RNA polymerase (RDR) are known as the three major gene families that act as the critical components of RNA interference or silencing mechanisms through the noncoding small RNA molecules (miRNA and siRNA) to regulate the expressions of protein-coding genes in eukaryotic organisms. However, most of their characteristics including structures, chromosomal location, subcellular locations, regulatory elements, and gene networking were not rigorously studied. Our analysis identified 7 TaDCL, 39 TaAGO, and 16 TaRDR genes as RNA interference (RNAi) genes from the wheat genome. Phylogenetic analysis of predicted RNAi proteins with the RNAi proteins of Arabidopsis and rice showed that the predicted proteins of TaDCL, TaAGO, and TaRDR groups are clustered into four, eight, and four subgroups, respectively. Domain, 3D protein structure, motif, and exon-intron structure analyses showed that these proteins conserve identical characteristics within groups and maintain differences between groups. The nonsynonymous/synonymous mutation ratio (Ka/Ks) < 1 suggested that these protein sequences conserve some purifying functions. RNAi genes networking with TFs revealed that ERF, MIKC-MADS, C2H2, BBR-BPC, MYB, and Dof are the key transcriptional regulators of the predicted RNAi-related genes. The cis-regulatory element (CREs) analysis detected some important CREs of RNAi genes that are significantly associated with light, stress, and hormone responses. Expression analysis based on an online database exhibited that almost all of the predicted RNAi genes are expressed in different tissues and organs. A case-control study from the gene expression level showed that some RNAi genes significantly responded to the drought and heat stresses. Overall results would therefore provide an excellent basis for in-depth molecular investigation of these genes and their regulatory elements for wheat crop improvement against different stressors.

Abscisic acid-polyacrylamide (ABA-PAM) treatment enhances forage grass growth and soil microbial diversity under drought stress

Drought restricts the growth of alpine grassland vegetation. This study aimed to explore a new technical system to improve the drought resistance of forage grass. Qinghai cold-land Poa pratensis seedlings were used in the drought stress experiment. A combination of abscisic acid (ABA) and polyacrylamide (PAM) were used to affect the growth, leaf physiology, soil enzyme activity, and rhizosphere microbial diversity of P. pratensis. The fresh leaf weight and root surface area were significantly increased after ABA-PAM combined treatment, while root length was significantly reduced. Besides, the leaf catalase (CAT) and superoxide dismutase (SOD) enzyme activity, proline and chlorophyll content, increased after the treatment, while malondialdehyde (MDA) content decreased. The treatment also increased sucrase, urease, and alkaline protease activities in rhizosphere soil, while decreasing acid phosphatase and neutral phosphatase enzyme activities. ABA-PAM combined treatment enhanced the rhizosphere microbial community and forage drought resistance by altering the abundance of various dominant microorganisms in the rhizosphere soil. The relative abundances of Actinobacteria, Chloroflexi, and Acidobacteria decreased, while Proteobacteria, Firmicutes, and Ascomycota increased. Unlike the relative abundance of Gibberella that decreased significantly, Komagataeibacter, Lactobacillus, Pichia, and Dekkera were significantly increased. Single-factor collinearity network analysis revealed a close relationship between the different rhizosphere microbial communities of forage grass, after ABA-PAM treatment. This study implies that ABA-PAM combined treatment can improve the drought resistance of forages. Therefore, it provides a theoretical and practical basis for restoring drought-induced grassland degradation.

Comparative Transcriptome Profiling Reveals Potential Candidate Genes, Transcription Factors, and Biosynthetic Pathways for Phosphite Response in Potato (Solanum tuberosum L.)

The study was conducted with C31 and C80 genotypes of the potato (Solanum tuberosum L.), which are tolerant and susceptible to phosphite (Phi, H₂PO₃), respectively. To decipher the molecular mechanisms underlying tolerance and susceptibility to Phi in the potato, RNA sequencing was used to study the global transcriptional patterns of the two genotypes. Media were prepared with 0.25 and 0.50 mM Phi, No-phosphorus (P), and 1.25 mM (phosphate, Pi as control). The values of fragments per kilobase of exon per million mapped fragments of the samples were also subjected to a principal component analysis, grouping the biological replicates of each sample. Using stringent criteria, a minimum of 819 differential (DEGs) were detected in both C80-Phi-0.25_vs_C80-Phi-0.50 (comprising 517 upregulated and 302 downregulated) and C80-Phi-0.50_vs_C80-Phi-0.25 (comprising 302 upregulated and 517 downregulated) and a maximum of 5214 DEGs in both C31-Con_vs_C31-Phi-0.25 (comprising 1947 upregulated and 3267 downregulated) and C31-Phi-0.25_vs_C31-Con (comprising 3267 upregulated and 1947 downregulated). DEGs related to the ribosome, plant hormone signal transduction, photosynthesis, and plant–pathogen interaction performed important functions under Phi stress, as shown by the Kyoto Encyclopedia of Genes and Genomes annotation. The expressions of transcription factors increased significantly in C31 compared with C80. For example, the expressions of Soltu.DM.01G047240, Soltu.DM.08G015900, Soltu.DM.06G012130, and Soltu.DM.08G012710 increased under P deficiency conditions (Phi-0.25, Phi-0.50, and No-P) relative to the control (P sufficiency) in C31. This study adds to the growing body of transcriptome data on Phi stress and provides important clues to the Phi tolerance response of the C31 genotype.

Drought tolerance improvement in Solanum lycopersicum: an insight into "OMICS" approaches and genome editing

Solanum lycopersicum (tomato) is an internationally acclaimed vegetable crop that is grown worldwide. However, drought stress is one of the most critical challenges for tomato production, and it is a crucial task for agricultural biotechnology to produce drought-resistant cultivars. Although breeders have done a lot of work on the tomato to boost quality and quantity of production and enhance resistance to biotic and abiotic stresses, conventional tomato breeding approaches have been limited to improving drought tolerance because of the intricacy of drought traits. Many efforts have been made to better understand the mechanisms involved in adaptation and tolerance to drought stress in tomatoes throughout the years. "Omics" techniques, such as genomics, transcriptomics, proteomics, and metabolomics in combination with modern sequencing technologies, have tremendously aided the discovery of drought-responsive genes. In addition, the availability of biotechnological tools, such as plant transformation and the recently developed genome editing system for tomatoes, has opened up wider opportunities for validating the function of drought-responsive genes and the generation of drought-tolerant varieties. This review highlighted the recent progresses for tomatoes improvement against drought stress through "omics" and "multi-omics" technologies including genetic engineering. We have also discussed the roles of non-coding RNAs and genome editing techniques for drought stress tolerance improvement in tomatoes.

Transcriptome Profiling Reveals Role of MicroRNAs and Their Targeted Genes during Adventitious Root Formation in Dark-Pretreated Micro-Shoot Cuttings of Tetraploid Robinia pseudoacacia L

Tetraploid Robinia pseudoacacia L. is a difficult-to-root species, and is vegetatively propagated through stem cuttings. Limited information is available regarding the adventitious root (AR) formation of dark-pretreated micro-shoot cuttings. Moreover, the role of specific miRNAs and their targeted genes during dark-pretreated AR formation under in vitro conditions has never been revealed. The dark pretreatment has successfully promoted and stimulated adventitious rooting signaling-related genes in tissue-cultured stem cuttings with the application of auxin (0.2 mg L-1 IBA). Histological analysis was performed for AR formation at 0, 12, 36, 48, and 72 h after excision (HAE) of the cuttings. The first histological events were observed at 36 HAE in the dark-pretreated cuttings; however, no cellular activities were observed in the control cuttings. In addition, the present study aimed to uncover the role of differentially expressed (DE) microRNAs (miRNAs) and their targeted genes during adventitious root formation using the lower portion (1-1.5 cm) of tetraploid R. pseudoacacia L. micro-shoot cuttings. The samples were analyzed using Illumina high-throughput sequencing technology for the identification of miRNAs at the mentioned time points. Seven DE miRNA libraries were constructed and sequenced. The DE number of 81, 162, 153, 154, 41, 9, and 77 miRNAs were upregulated, whereas 67, 98, 84, 116, 19, 16, and 93 miRNAs were downregulated in the following comparisons of the libraries: 0-vs-12, 0-vs-36, 0-vs-48, 0-vs-72, 12-vs-36, 36-vs-48, and 48-vs-72, respectively. Furthermore, we depicted an association between ten miRNAs (novel-m0778-3p, miR6135e.2-5p, miR477-3p, miR4416c-5p, miR946d, miR398b, miR389a-3p, novel m0068-5p, novel-m0650-3p, and novel-m0560-3p) and important target genes (auxin response factor-3, gretchen hagen-9, scarecrow-like-1, squamosa promoter-binding protein-like-12, small auxin upregulated RNA-70, binding protein-9, vacuolar invertase-1, starch synthase-3, sucrose synthase-3, probable starch synthase-3, cell wall invertase-4, and trehalose phosphatase synthase-5), all of which play a role in plant hormone signaling and starch and sucrose metabolism pathways. The quantitative polymerase chain reaction (qRT-PCR) was used to validate the relative expression of these miRNAs and their targeted genes. These results provide novel insights and a foundation for further studies to elucidate the molecular factors and processes controlling AR formation in woody plants.

β-Cyclocitral, a Master Regulator of Multiple Stress-Responsive Genes in Solanum lycopersicum L. Plants

β-cyclocitral (βCC), a major apocarotenoid of β-carotene, enhances plants' defense against environmental stresses. However, the knowledge of βCC's involvement in the complex stress-signaling network is limited. Here we demonstrate how βCC reprograms the transcriptional responses that enable Solanum lycopersicum L. (tomato) plants to endure a plethora of environmental stresses. Comparative transcriptome analysis of control and βCC-treated tomato plants was done by generating RNA sequences in the BGISEQ-500 platform. The trimmed sequences were mapped on the tomato reference genome that identifies 211 protein-coding differentially expressed genes. Gene ontology and Kyoto Encyclopedia of Genes and Genomes analysis and their enrichment uncovered that only upregulated genes are attributed to the stress response. Moreover, 80% of the upregulated genes are functionally related to abiotic and biotic stresses. Co-functional analysis of stress-responsive genes revealed a network of 18 genes that code for heat shock proteins, transcription factors (TFs), and calcium-binding proteins. The upregulation of jasmonic acid (JA)-dependent TFs (MYC2, MYB44, ERFs) but not the JA biosynthetic genes is surprising. However, the upregulation of DREB3, an abscisic acid (ABA)-independent TF, validates the unaltered expression of ABA biosynthetic genes. We conclude that βCC treatment upregulates multiple stress-responsive genes without eliciting JA and ABA biosynthesis.

An Update on Crop ABA Receptors

The hormone abscisic acid (ABA) orchestrates the plant stress response and regulates sophisticated metabolic and physiological mechanisms essential for survival in a changing environment. Plant ABA receptors were described more than 10 years ago, and a considerable amount of information is available for the model plant Arabidopsis thaliana. Unfortunately, this knowledge is still very limited in crops that hold the key to feeding a growing population. In this review, we summarize genomic, genetic and structural data obtained in crop ABA receptors. We also provide an update on ABA perception in major food crops, highlighting specific and common features of crop ABA receptors.

Exogenous 1',4'-trans-Diol-ABA Induces Stress Tolerance by Affecting the Level of Gene Expression in Tobacco (Nicotiana tabacum L.)

1',4'-trans-diol-ABA is a key precursor of the biosynthesis of abscisic acid (ABA) biosynthesis in fungi. We successfully obtained the pure compound from a mutant of Botrytis cinerea and explored its function and possible mechanism on plants by spraying 2 mg/L 1',4'-trans-diol-ABA on tobacco leaves. Our results showed that this compound enhanced the drought tolerance of tobacco seedlings. A comparative transcriptome analysis showed that a large number of genes responded to the compound, exhibiting 1523 genes that were differentially expressed at 12 h, which increased to 1993 at 24 h and 3074 at 48 h, respectively. The enrichment analysis demonstrated that the differentially expressed genes (DEGs) were primarily enriched in pathways related to hormones and resistance. The DEGs of transcription factors were generally up-regulated and included the bHLH, bZIP, ERF, MYB, NAC, WRKY and HSF families. Moreover, the levels of expression of PYL/PYR, PP2C, SnRK2, and ABF at the ABA signaling pathway responded positively to exogenous 1',4'-trans-diol-ABA. Among them, seven ABF transcripts that were detected were significantly up-regulated. In addition, the genes involved in salicylic acid, ethylene and jasmonic acid pathways, reactive oxygen species scavenging system, and other resistance related genes were primarily induced by 1',4'-trans-diol-ABA. These findings indicated that treatment with 1',4'-trans-diol-ABA could improve tolerance to plant abiotic stress and potential biotic resistance by regulating gene expression, similar to the effects of exogenous ABA.

Feline calicivirus strain 2280 p30 antagonizes type I interferon-mediated antiviral innate immunity through directly degrading IFNAR1 mRNA

Feline calicivirus (FCV) belongs to the Caliciviridae, which comprises small RNA viruses of both medical and veterinary importance. Once infection has occurred, FCV can persist in the cat population, but the molecular mechanism of how it escapes the innate immune response is still unknown. In this study, we found FCV strain 2280 to be relatively resistant to treatment with IFN-β. FCV 2280 infection inhibited IFN-induced activation of the ISRE (Interferon-stimulated response element) promoter and transcription of ISGs (Interferon-stimulated genes). The mechanistic analysis showed that the expression of IFNAR1, but not IFNAR2, was markedly reduced in FCV 2280-infected cells by inducing the degradation of IFNAR1 mRNA, which inhibited the phosphorylation of downstream adaptors. Further, overexpression of the FCV 2280 nonstructural protein p30, but not p30 of the attenuated strain F9, downregulated the expression of IFNAR1 mRNA. His-p30 fusion proteins were produced in Escherichia coli and purified, and an in vitro digestion assay was performed. The results showed that 2280 His-p30 could directly degrade IFNAR1 RNA but not IFNAR2 RNA. Moreover, the 5’UTR of IFNAR1 mRNA renders it directly susceptible to cleavage by 2280 p30. Next, we constructed two chimeric viruses: rFCV 2280-F9 p30 and rFCV F9-2280 p30. Compared to infection with the parental virus, rFCV 2280-F9 p30 infection displayed attenuated activities in reducing the level of IFNAR1 and inhibiting the phosphorylation of STAT1 and STAT2, whereas rFCV F9-2280 p30 displayed enhanced activities. Animal experiments showed that the virulence of rFCV 2280-F9 p30 infection was attenuated but that the virulence of rFCV F9-2280 p30 was increased compared to that of the parental viruses. Collectively, these data show that FCV 2280 p30 could directly and selectively degrade IFNAR1 mRNA, thus blocking the type I interferon-induced activation of the JAK-STAT signalling pathway, which may contribute to the pathogenesis of FCV infection.

Vaccination against FCV has been available for many years and has effectively reduced the incidence of clinical disease. However, vaccines cannot prevent infection, and vaccinated cats can still become persistently infected by FCV, suggesting that FCV has evolved several strategies for counteracting various components of the innate and adaptive immune systems. Here, we show that FCV strain 2280 is resistant to the antiviral effect of IFN. The molecular mechanism by which this occurs is that FCV 2280 infection blocks the JAK-STAT pathway through promoting the degradation of IFNAR1 mRNA by FCV p30 protein. An in vitro degradation assay demonstrated that 2280 p30, but not p30 of the vaccine strain F9, could directly and selectively decay IFNAR1 RNA. The exchange of p30 between 2280 and F9 strains using a reverse genetic system also showed that 2280 p30 is a key factor that contributes to the resistance to IFN and enhances virulence. Our findings reveal a new mechanism evolved by FCV to circumvent the host antiviral response.

SlEAD1, an EAR motif-containing ABA down-regulated novel transcription repressor regulates ABA response in tomato

EAR motif-containing proteins are able to repress gene expression, therefore play important roles in regulating plants growth and development, plant response to environmental stimuli, as well as plant hormone signal transduction. ABA is a plant hormone that regulates abiotic stress tolerance in plants via signal transduction. ABA signaling via the PYR1/PYLs/RCARs receptors, the PP2Cs phosphatases, and SnRK2s protein kinases activates the ABF/AREB/ABI5-type bZIP transcription factors, resulting in the activation/repression of ABA response genes. However, functions of many ABA response genes remained largely unknown. We report here the identification of the ABA-responsive gene SlEAD1 (Solanum lycopersicum EAR motif-containing ABA down-regulated 1) as a novel EAR motif-containing transcription repressor gene in tomato. We found that the expression of SlEAD1 was down-regulated by ABA treatment, and SlEAD1 repressed reporter gene expression in transfected protoplasts. By using CRISPR gene editing, we generated transgene-free slead1 mutants and found that the mutants produced short roots. By using seed germination and root elongation assays, we examined ABA response of the slead1 mutants and found that ABA sensitivity in the mutants was increased. By using qRT-PCR, we further show that the expression of some of the ABA biosynthesis and signaling component genes were increased in the slead1 mutants. Taken together, our results suggest that SlEAD1 is an ABA response gene, that SlEAD1 is a novel EAR motif-containing transcription repressor, and that SlEAD1 negatively regulates ABA responses in tomato possibly by repressing the expression of some ABA biosynthesis and signaling genes.

Identification of tRFs and phasiRNAs in tomato (Solanum lycopersicum) and their responses to exogenous abscisic acid

BACKGROUND

The non-coding small RNA tRFs (tRNA-derived fragments) and phasiRNAs (plant-specific) exert important roles in plant growth, development and stress resistances. However, whether the tRFs and phasiRNAs respond to the plant important stress hormone abscisic acid (ABA) remain enigma.

RESULTS

Here, the RNA-sequencing was implemented to decipher the landscape of tRFs and phasiRNAs in tomato (Solanum lycopersicum) leaves and their responses when foliar spraying exogenous ABA after 24 h. In total, 733 tRFs and 137 phasiRNAs were detected. The tRFs were mainly derived from the tRNA^Ala transporting alanine, which tended to be cleaved at the 5^'terminal guanine site and D loop uracil site to produce tRF^Ala with length of 20 nt. Most of phasiRNAs originated from NBS-LRR resistance genes. Expression analysis revealed that 156 tRFs and 68 phasiRNAs expressed differentially, respectively. Generally, exogenous ABA mainly inhibited the expression of tRFs and phasiRNAs. Furthermore, integrating analysis of target gene prediction and transcriptome data presented that ABA significantly downregulated the abundance of phsaiRNAs associated with biological and abiotic resistances. Correspondingly, their target genes such as AP2/ERF, WRKY and NBS-LRR, STK and RLK, were mainly up-regulated.

CONCLUSIONS

Combined with the previous analysis of ABA-response miRNAs, it was speculated that ABA can improve the plant resistances to various stresses by regulating the expression and interaction of small RNAs (such as miRNAs, tRFs, phasiRNAs) and their target genes. This study enriches the plant tRFs and phasiRNAs, providing a vital basis for further investigating ABA response-tRFs and phasiRNAs and their functions in biotic and abiotic stresses.

Fusarium oxysporum infection activates the plastidial branch of the terpenoid biosynthesis pathway in flax, leading to increased ABA synthesis

Upregulation of the terpenoid pathway and increased ABA content in flax upon Fusarium infection leads to activation of the early plant's response (PR genes, cell wall remodeling, and redox status). Plants have developed a number of defense strategies against the adverse effects of fungi such as Fusarium oxysporum. One such defense is the production of antioxidant secondary metabolites, which fall into two main groups: the phenylpropanoids and the terpenoids. While functions and biosynthesis of phenylpropanoids have been extensively studied, very little is known about the genes controlling the terpenoid synthesis pathway in flax. They can serve as antioxidants, but are also substrates for a plethora of different compounds, including those of regulatory functions, like ABA. ABA's function during pathogen attack remains obscure and often depends on the specific plant-pathogen interactions. In our study we showed that in flax the non-mevalonate pathway is strongly activated in the early hours of pathogen infection and that there is a redirection of metabolites towards ABA synthesis. The elevated synthesis of ABA correlates with flax resistance to F. oxysporum, thus we suggest ABA to be a positive regulator of the plant's early response to the infection.

Bacillus amyloliquefaciens MBI600 differentially induces tomato defense signaling pathways depending on plant part and dose of application

The success of Bacillus amyloliquefaciens as a biological control agent relies on its ability to outgrow plant pathogens. It is also thought to interact with its plant host by inducing systemic resistance. In this study, the ability of B. amyloliquefaciens MBI600 to elicit defense (or other) responses in tomato seedlings and plants was assessed upon the expression of marker genes and transcriptomic analysis. Spray application of Serifel, a commercial formulation of MBI600, induced responses in a dose-dependent manner. Low dosage primed plant defense by activation of SA-responsive genes. Suggested dosage induced defense by mediating synergistic cross-talk between JA/ET and SA-signaling. Saturation of tomato roots or leaves with MBI600 elicitors activated JA/ET signaling at the expense of SA-mediated responses. The complex signaling network that is implicated in MBI600-tomato seedling interactions was mapped. MBI600 and flg22 (a bacterial flagellin peptide) elicitors induced, in a similar manner, biotic and abiotic stress responses by the coordinated activation of genes involved in JA/ET biosynthesis as well as hormone and redox signaling. This is the first study to suggest the activation of plant defense following the application of a commercial microbial formulation under conditions of greenhouse crop production.

Advances in Omics Approaches for Abiotic Stress Tolerance in Tomato

Tomato, one of the most important crops worldwide, has a high demand in the fresh fruit market and processed food industries. Despite having considerably high productivity, continuous supply as per the market demand is hard to achieve, mostly because of periodic losses occurring due to biotic as well as abiotic stresses. Although tomato is a temperate crop, it is grown in almost all the climatic zones because of widespread demand, which makes it challenge to adapt in diverse conditions. Development of tomato cultivars with enhanced abiotic stress tolerance is one of the most sustainable approaches for its successful production. In this regard, efforts are being made to understand the stress tolerance mechanism, gene discovery, and interaction of genetic and environmental factors. Several omics approaches, tools, and resources have already been developed for tomato growing. Modern sequencing technologies have greatly accelerated genomics and transcriptomics studies in tomato. These advancements facilitate Quantitative trait loci (QTL) mapping, genome-wide association studies (GWAS), and genomic selection (GS). However, limited efforts have been made in other omics branches like proteomics, metabolomics, and ionomics. Extensive cataloging of omics resources made here has highlighted the need for integration of omics approaches for efficient utilization of resources and a better understanding of the molecular mechanism. The information provided here will be helpful to understand the plant responses and the genetic regulatory networks involved in abiotic stress tolerance and efficient utilization of omics resources for tomato crop improvement.

Transcriptome-Wide Identification and Characterization of Circular RNAs in Leaves of Chinese Cabbage (Brassica rapa L. ssp. pekinensis) in Response to Calcium Deficiency-Induced Tip-burn

Circular RNA (circRNA) is a newly discovered non-coding RNA, which play significant roles in the function and transcriptional regulation of microRNA. To date, in Chinese cabbage, the functional characteristic of circRNAs in response to calcium deficiency-induced tip-burn have not been reported. In this study, 730 circRNAs were isolated from Chinese cabbage leaves, of which 23 and 22 were differentially expressed in different calcium deficiency stages compared with the control. Forty-six host genes of the differentially expressed circRNAs were identified, and one circRNA was found to act as miRNAs sponges. Based on the functional analysis of host genes and target mRNAs of the corresponding miRNAs, the identified circRNAs might participated in response to stimulus, electron carrier activity, ATPase activity, cell wall metabolism, transcription factors and plant hormone signal transduction. ABF2, a positive regulator of the abiotic stress response in the abscisic acid (ABA) pathway, may play a role in calcium deficiency tolerance through a circRNA regulatory pathway. Correspondingly, the concentration of ABA is also increased during the Ca²⁺ deficiency stress. Our results suggest that circRNAs participate in a broad range of biological processes and physiological functions in the response to calcium deficiency-induced tip-burn and provide a basis for further studies of the biological roles that circRNAs play in the plant stress response.

Involvement of abscisic acid-responsive element-binding factors in cassava (Manihot esculenta) dehydration stress response

Cassava (Manihot esculenta) is a major staple food, animal feed and energy crop in the tropics and subtropics. It is one of the most drought-tolerant crops, however, the mechanisms of cassava drought tolerance remain unclear. Abscisic acid (ABA)-responsive element (ABRE)-binding factors (ABFs) are transcription factors that regulate expression of target genes involved in plant tolerance to drought, high salinity, and osmotic stress by binding ABRE cis-elements in the promoter regions of these genes. However, there is little information about ABF genes in cassava. A comprehensive analysis of Manihot esculenta ABFs (MeABFs) described the phylogeny, genome location, cis-acting elements, expression profiles, and regulatory relationship between these factors and Manihot esculenta betaine aldehyde dehydrogenase genes (MeBADHs). Here we conducted genome-wide searches and subsequent molecular cloning to identify seven MeABFs that are distributed unevenly across six chromosomes in cassava. These MeABFs can be clustered into three groups according to their phylogenetic relationships to their Arabidopsis (Arabidopsis thaliana) counterparts. Analysis of the 5′-upstream region of MeABFs revealed putative cis-acting elements related to hormone signaling, stress, light, and circadian clock. MeABF expression profiles displayed clear differences among leaf, stem, root, and tuberous root tissues under non-stress and drought, osmotic, or salt stress conditions. Drought stress in cassava leaves and roots, osmotic stress in tuberous roots, and salt stress in stems induced expression of the highest number of MeABFs showing significantly elevated expression. The glycine betaine (GB) content of cassava leaves also was elevated after drought, osmotic, or salt stress treatments. BADH1 is involved in GB synthesis. We show that MeBADH1 promoter sequences contained ABREs and that MeBADH1 expression correlated with MeABF expression profiles in cassava leaves after the three stress treatments. Taken together, these results suggest that in response to various dehydration stresses, MeABFs in cassava may activate transcriptional expression of MeBADH1 by binding the MeBADH1 promoter that in turn promotes GB biosynthesis and accumulation via an increase in MeBADH1 gene expression levels and MeBADH1 enzymatic activity. These responses protect cells against dehydration stresses by preserving an osmotic balance that enhances cassava tolerance to dehydration stresses.

Global transcriptome analysis provides new insights in Thellungiella salsuginea stress response

Thellungiella salsuginea is highly tolerant to abiotic stress, while its a close relative Arabidopsis thaliana is sensitive to stress. This characteristic makes T. salsuginea an excellent model for uncovering the mechanisms of abiotic stress tolerance. Abscisic acid (ABA) plays essential roles in plant abiotic and biotic stress tolerance. To test the changes in gene expression of T. salsuginea under ABA treatment, in this study, the transcriptomes of T. salsuginea roots and leaves were compared in response to exogenously application of ABA. The results showed that ABA treatment caused different expression of 2,200 and 3,305 genes in leaves and roots, respectively, compared with the untreated control. In particular, genes encoding transcription factors such as WRKY, MYB, NAC, GATA, ethylene-responsive factors (ERFs), heat stress transcription factors, basic helix-loop-helix, PLATZ and B3 domain-containing family members were enriched. In addition, 49 and 114 differentially expressed genes were identified as ABA-regulated genes, separately in leaves and roots, respectively, which were related to biotic and abiotic stresses. The expression levels of some genes were validated by qRT-PCR. Different responses of genes to ABA treatment were discovered in T. salsuginea and A. thaliana. This transcriptome analysis expands our understanding of the role of ABA in stress tolerance in T. salsuginea. Our study provides a wealth of information for improving stress tolerance in crop plants.

Using Transcriptome to Discover a Novel Melatonin-Induced Sodic Alkaline Stress Resistant Pathway in Solanum lycopersicum L

Melatonin plays important roles in multiple stress responses. However, the downstream signaling pathway and molecular mechanism are unclear until now. Here, we not only revealed the transcriptional control of melatonin-induced sodic alkaline stress tolerance, but also described a screen for key downstream transcriptional factors of melatonin through transcriptome analysis. The melatonin-induced transcriptional network of hormone, transcriptional factors and functional genes has been established under both control and stress conditions. Among these, six candidates of transcriptional factors have been identified via Gene Ontology and Kyoto Encyclopedia of Genes and Genomes analysis. Using the virus-induced gene silencing approach, we confirmed that DREB1α and IAA3 were key downstream transcriptional factors of melatonin-induced sodic alkaline stress tolerance at the genetic level. The transcriptions of DREB1α and IAA3 could be activated by melatonin or sodic alkaline treatment. Interestingly, we found that DREB1α could directly upregulate the expression of IAA3 by binding to its promoters. Moreover, several physiological processes of Na+ detoxification, dehydration resistance, high pH buffering and reactive oxygen species scavenging were confirmed to depend or partly depend on DREB1α and IAA3 pathway in melatonin-induced stress tolerance. Taken together, this study suggested that DREB1α and IAA3 are positive resistant modulators, and provided a direct link among melatonin, DREB1α and IAA3 in the sodic alkaline stress tolerance activating in tomato plants.

Transcriptome-based analysis of carotenoid accumulation-related gene expression in petals of Chinese cabbage (Brassica rapa L.)

To identify genes associated with carotenoid accumulation in petals of Chinese cabbage, the composition and content of carotenoids were analyzed, and comparative transcriptome sequencing was performed between the yellow flower line, 92S105, and the orange flower line, 94C9. High-performance liquid chromatography (HPLC) revealed that petals of 92S105 were high in violaxanthin as well as lutein, whereas petals of 94C9 showed considerable levels of lutein and β-carotene. Transcriptome analysis showed that 3534 and 3833 genes were up- and down-regulated in 94C9, respectively. Among these differentially expressed genes (DEGs), many related to carotenoid accumulation were identified, including 12 carotenoid biosynthesis pathway genes, 4 transcription factor genes, and 1028 specifically expressed genes. β-carotene hydroxylase 1 (BrBCH1), BrBCH2, zeaxanthin epoxidase (BrZEP), and MYB transcription factor gene (BrGAMYB) were down-regulated in petals of 94C9 when compared with petals of 92S105, which caused β-carotene accumulation and may lead to orange petal color in 94C9. Expression levels of 20 DEGs were verified by qPCR and the results were highly consistent with those of transcriptome sequencing. Moreover, Gene Ontology (GO) enrichment analysis revealed that membrane, binding, and metabolic processes were the most significantly enriched GO terms in cellular component, molecular function, and biological process ontologies, respectively. In conclusion, our study analyzed the differences in composition and content of carotenoids between 92S105 and 94C9 and identified potential candidate genes related to carotenoid accumulation in petals, thereby creating a solid foundation for future studies on the mechanism regulating carotenoid accumulation in petals of Chinese cabbage.

Transcript analyses reveal a comprehensive role of abscisic acid in modulating fruit ripening in Chinese jujube

BACKGROUND

Chinese jujube (Ziziphus jujuba Mill.) is a non-climacteric fruit; however, the underlying mechanism of ripening and the role of abscisic acid involved in this process are not yet understood for this species.

RESULTS

In the present study, a positive correlation between dynamic changes in endogenous ABA and the onset of jujube ripening was determined. Transcript analyses suggested that the expression balance among genes encoding nine-cis-epoxycarotenoid dioxygenase (ZjNCED3), ABA-8'-hydroxylase (ZjCYP707A2), and beta-glucosidase (ZjBG4, ZjBG5, ZjBG8, and ZjBG9) has an important role in maintaining ABA accumulation, while the expression of a receptor (ZjPYL8), protein phosphatase 2C (ZjPP2C4-8), and sucrose nonfermenting 1-related protein kinase 2 (ZjSnRK2-2 and ZjSnRK2-5) is important in regulating fruit sensitivity to ABA applications. In addition, white mature 'Dongzao' fruit were harvested and treated with 50 mg L- 1 ABA or 50 mg L- 1 nordihydroguaiaretic acid (NDGA) to explore the role of ABA in jujube fruit ripening. By comparative transcriptome analyses, 1103 and 505 genes were differentially expressed in response to ABA and NDGA applications on the 1st day after treatment, respectively. These DEGs were associated with photosynthesis, secondary, lipid, cell wall, and starch and sugar metabolic processes, suggesting the involvement of ABA in modulating jujube fruit ripening. Moreover, ABA also exhibited crosstalk with other phytohormones and transcription factors, indicating a regulatory network for jujube fruit ripening.

CONCLUSIONS

Our study further elucidated ABA-associated metabolic and regulatory processes. These findings are helpful for improving strategies for jujube fruit storage and for gaining insights into understand complex non-climacteric fruit ripening processes.

Expanding Alternative Splicing Identification by Integrating Multiple Sources of Transcription Data in Tomato

Tomato (Solanum lycopersicum) is an important vegetable and fruit crop. Its genome was completely sequenced and there are also a large amount of available expressed sequence tags (ESTs) and short reads generated by RNA sequencing (RNA-seq) technologies. Mapping transcripts including mRNA sequences, ESTs, and RNA-seq reads to the genome allows identifying pre-mRNA alternative splicing (AS), a post-transcriptional process generating two or more RNA isoforms from one pre-mRNA transcript. We comprehensively analyzed the AS landscape in tomato by integrating genome mapping information of all available mRNA and ESTs with mapping information of RNA-seq reads which were collected from 27 published projects. A total of 369,911 AS events were identified from 34,419 genomic loci involving 161,913 transcripts. Within the basic AS events, intron retention is the prevalent type (18.9%), followed by alternative acceptor site (12.9%) and alternative donor site (7.3%), with exon skipping as the least type (6.0%). Complex AS types having two or more basic event accounted for 54.9% of total AS events. Within 35,768 annotated protein-coding gene models, 23,233 gene models were found having pre-mRNAs generating AS isoform transcripts. Thus the estimated AS rate was 65.0% in tomato. The list of identified AS genes with their corresponding transcript isoforms serves as a catalog for further detailed examination of gene functions in tomato biology. The post-transcriptional information is also expected to be useful in improving the predicted gene models in tomato. The sequence and annotation information can be accessed at plant alternative splicing database (http://proteomics.ysu.edu/altsplice).

Comparative Transcriptomic and Physiological Analyses of Medicago sativa L. Indicates that Multiple Regulatory Networks Are Activated during Continuous ABA Treatment

Alfalfa is the most extensively cultivated forage legume worldwide. However, the molecular mechanisms underlying alfalfa responses to exogenous abscisic acid (ABA) are still unknown. In this study, the first global transcriptome profiles of alfalfa roots under ABA treatments for 1, 3 and 12 h (three biological replicates for each time point, including the control group) were constructed using a BGISEQ-500 sequencing platform. A total of 50,742 isoforms with a mean length of 2541 bp were generated, and 4944 differentially expressed isoforms (DEIs) were identified after ABA deposition. Metabolic analyses revealed that these DEIs were involved in plant hormone signal transduction, transcriptional regulation, antioxidative defense and pathogen immunity. Notably, several well characterized hormone signaling pathways, for example, the core ABA signaling pathway, was activated, while salicylic acid, jasmonate and ethylene signaling pathways were mainly suppressed by exogenous ABA. Moreover, the physiological work showed that catalase and peroxidase activity and glutathione and proline content were increased after ABA deposition, which is in accordance with the dynamic transcript profiles of the relevant genes in antioxidative defense system. These results indicate that ABA has the potential to improve abiotic stress tolerance, but that it may negatively regulate pathogen resistance in alfalfa.

Transcriptome Sequence Analysis Elaborates a Complex Defensive Mechanism of Grapevine (Vitis vinifera L.) in Response to Salt Stress

Salinity is ubiquitous abiotic stress factor limiting viticulture productivity worldwide. However, the grapevine is vulnerable to salt stress, which severely affects growth and development of the vine. Hence, it is crucial to delve into the salt resistance mechanism and screen out salt-resistance prediction marker genes; we implicated RNA-sequence (RNA-seq) technology to compare the grapevine transcriptome profile to salt stress. Results showed 2472 differentially-expressed genes (DEGs) in total in salt-responsive grapevine leaves, including 1067 up-regulated and 1405 down-regulated DEGs. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) annotations suggested that many DEGs were involved in various defense-related biological pathways, including ROS scavenging, ion transportation, heat shock proteins (HSPs), pathogenesis-related proteins (PRs) and hormone signaling. Furthermore, many DEGs were encoded transcription factors (TFs) and essential regulatory proteins involved in signal transduction by regulating the salt resistance-related genes in grapevine. The antioxidant enzyme analysis showed that salt stress significantly affected the superoxide dismutase (SOD), peroxidase (POD), catalase (CAT) and glutathione S-transferase (GST) activities in grapevine leaves. Moreover, the uptake and distribution of sodium (Na⁺), potassium (K⁺) and chlorine (Cl⁻) in source and sink tissues of grapevine was significantly affected by salt stress. Finally, the qRT-PCR analysis of DE validated the data and findings were significantly consistent with RNA-seq data, which further assisted in the selection of salt stress-responsive candidate genes in grapevine. This study contributes in new perspicacity into the underlying molecular mechanism of grapevine salt stress-tolerance at the transcriptome level and explore new approaches to applying the gene information in genetic engineering and breeding purposes.

Factors Influencing Gene Family Size Variation Among Related Species in a Plant Family, Solanaceae

Gene duplication and loss contribute to gene content differences as well as phenotypic divergence across species. However, the extent to which gene content varies among closely related plant species and the factors responsible for such variation remain unclear. Here, using the Solanaceae family as a model and Pfam domain families as a proxy for gene families, we investigated variation in gene family sizes across species and the likely factors contributing to the variation. We found that genes in highly variable families have high turnover rates and tend to be involved in processes that have diverged between Solanaceae species, whereas genes in low-variability families tend to have housekeeping roles. In addition, genes in high- and low-variability gene families tend to be duplicated by tandem and whole genome duplication, respectively. This finding together with the observation that genes duplicated by different mechanisms experience different selection pressures suggest that duplication mechanism impacts gene family turnover. We explored using pseudogene number as a proxy for gene loss but discovered that a substantial number of pseudogenes are actually products of pseudogene duplication, contrary to the expectation that most plant pseudogenes are remnants of once-functional duplicates. Our findings reveal complex relationships between variation in gene family size, gene functions, duplication mechanism, and evolutionary rate. The patterns of lineage-specific gene family expansion within the Solanaceae provide the foundation for a better understanding of the genetic basis underlying phenotypic diversity in this economically important family.

Analysis of the regulation networks in grapevine reveals response to waterlogging stress and candidate gene-marker selection for damage severity

Owing to the climate change impacts, waterlogging is one of the most hazardous abiotic stresses to crops, which also can result in a serious reduction in the quantity and quality of grape berry and wine production during the rainy season. Therefore, the exploration of the response mechanism of grape to waterlogging is necessary, for which the analysis of the transcriptomic regulation networks of grapevine leaves in response to waterlogging stress was carried out. In this study, 12 634 genes were detected in both waterlogging stress and control grapevine plants, out of which 6837 genes were differentially expressed. A comparative analysis revealed that genes functioning in the antioxidant system, glycolysis and fermentation pathway, chlorophyll metabolism, amino acid metabolism and hormones were activated to reduce injury to grapes under the waterlogging stress. Meanwhile, genes encoding class-2 non-symbiotic haemoglobin were determined as important in waterlogging acclimation. Additionally, the expression variations of three marker genes were found to be informative and can be used to predict the viability of the grapevines subjected to waterlogging. This research not only probes the molecular mechanism underlying grapevine waterlogging tolerance but also puts forward an idea about the application of gene expression information to practical management.

Genome-wide transcriptomic analysis of BR-deficient Micro-Tom reveals correlations between drought stress tolerance and brassinosteroid signaling in tomato

Brassinosteroids (BRs) are plant steroid hormones that play crucial roles in a range of growth and developmental processes. Although BR signal transduction and biosynthetic pathways have been well characterized in model plants, their biological roles in an important crop, tomato (Solanum lycopersicum), remain unknown. Here, cultivated tomato (WT) and a BR synthesis mutant, Micro-Tom (MT), were compared using physiological and transcriptomic approaches. The cultivated tomato showed higher tolerance to drought and osmotic stresses than the MT tomato. However, BR-defective phenotypes of MT, including plant growth and stomatal closure defects, were completely recovered by application of exogenous BR or complementation with a SlDWARF gene. Using genome-wide transcriptome analysis, 619 significantly differentially expressed genes (DEGs) were identified between WT and MT plants. Several DEGs were linked to known signaling networks, including those related to biotic/abiotic stress responses, lignification, cell wall development, and hormone responses. Consistent with the higher susceptibility of MT to drought stress, several gene sets involved in responses to drought and osmotic stress were differentially regulated between the WT and MT tomato plants. Our data suggest that BR signaling pathways are involved in mediating the response to abiotic stress via fine-tuning of abiotic stress-related gene networks in tomato plants.

Genetic mapping identifies loci that influence tomato resistance against Colorado potato beetles

The Colorado potato beetle (CPB; Leptinotarsa decemlineata Say), the most economically important insect pest on potato (Solanum tuberosum L.), also feeds on other Solanaceae, including cultivated tomato (Solanum lycopersicum L.). We used tomato genetic mapping populations to investigate natural variation in CPB resistance. CPB bioassays with 74 tomato lines carrying introgressions of Solanum pennellii in S. lycopersicum cv. M82 identified introgressions from S. pennellii on chromosomes 1 and 6 conferring CPB susceptibility, whereas introgressions on chromosomes 1, 8 and 10 conferred higher resistance. Mapping of CPB resistance using 113 recombinant inbred lines derived from a cross between S. lycopersicum cv UC-204B and Solanum galapagense identified significant quantitative trait loci on chromosomes 6 and 8. In each case, the S. galapagense alleles were associated with lower leaf damage and reduced larval growth. Results of both genetic mapping approaches converged on the same region of chromosome 6, which may have important functions in tomato defense against CPB herbivory. Although genetic mapping identified quantitative trait loci encompassing known genes for tomato acyl sugar and glycoalkaloid biosynthesis, experiments with acyl sugar near-isogenic lines and transgenic GAME9 glycoalkaloid-deficient and overproducing lines showed no significant effect of these otherwise insect-defensive metabolites on CPB performance.

Comparative Transcriptome Analyses Reveal Potential Mechanisms of Enhanced Drought Tolerance in Transgenic Salvia Miltiorrhiza Plants Expressing AtDREB1A from Arabidopsis

In our previous study, drought-resistant transgenic plants of Salvia miltiorrhiza were produced via overexpression of the transcription factor AtDREB1A. To unravel the molecular mechanisms underpinning elevated drought tolerance in transgenic plants, in the present study we compared the global transcriptional profiles of wild-type (WT) and AtDREB1A-expressing transgenic plants using RNA-sequencing (RNA-seq). Using cluster analysis, we identified 3904 differentially expressed genes (DEGs). Compared with WT plants, 423 unigenes were up-regulated in pRD29A::AtDREB1A-31 before drought treatment, while 936 were down-regulated and 1580 and 1313 unigenes were up- and down-regulated after six days of drought. COG analysis revealed that the 'signal transduction mechanisms' category was highly enriched among these DEGs both before and after drought stress. Based on the Kyoto Encyclopedia of Genes and Genomes (KEGG) annotation, DEGs associated with "ribosome", "plant hormone signal transduction", photosynthesis", "plant-pathogen interaction", "glycolysis/gluconeogenesis" and "carbon fixation" are hypothesized to perform major functions in drought resistance in AtDREB1A-expressing transgenic plants. Furthermore, the number of DEGs associated with different transcription factors increased significantly after drought stress, especially the AP2/ERF, bZIP and MYB protein families. Taken together, this study substantially expands the transcriptomic information for S. miltiorrhiza and provides valuable clues for elucidating the mechanism of AtDREB1A-mediated drought tolerance in transgenic plants.

Transcriptome profiling during mangrove viviparity in response to abscisic acid

Mangrove plants adapt to coastal tidal mudflats with specially evolved viviparity seed development. However, very little is known about the genetic and molecular mechanisms of mangrove viviparity. Here, we tested a hypothesis that plant hormone abscisic acid (ABA) plays a significant role in precocious germination of viviparous Kandelia obovata seeds by exogenous applications. Through transcriptome analysis of ABA treated seeds, it was found that ABA repressed mangrove fruit growth and development, and there were thousands of genes differentially expressed. As a result, dynamics of the pathways were dramatically altered. In particular, "Plant hormone signal transduction" and "MAPK signaling pathway" were represented significantly. Among differentially expressed genes, some key genes of ABA signal transduction were induced, while ABA biosynthesis genes were repressed. Take ABI1 and ABI2, key negative regulators in ABA signal pathway, as examples, homologous alignment and a phylogenetic tree in various species showed that ABI1 and ABI2 are highly conserved among various species. The functional similarity of these genes was confirmed by transgenic work in Arabidopsis. Taken together, ABA inhibited mangrove viviparity, but mangroves developed a mechanism to prevent accidently increase of ABA in the harsh environment for maintaining viviparous reproductive strategy.

Up-Regulation of HSFA2c and HSPs by ABA Contributing to Improved Heat Tolerance in Tall Fescue and Arabidopsis

Abscisic acid (ABA) is known to play roles in regulating plant tolerance to various abiotic stresses, but whether ABA’s effects on heat tolerance are associated with its regulation of heat stress transcription factors (HSFs) and heat shock proteins (HSPs) is not well documented. The objective of this study was to determine whether improved heat tolerance of tall fescue (Festuca arundinacea Schreb.) by ABA was through the regulation of HSFs and HSPs. ABA-responsive transcriptional factors, ABA-responsive element binding protein 3 (FaAREB3) and dehydration-responsive element binding protein 2A (FaDREB2A) of tall fescue, were able to bind to the cis-elements in the promoter of tall fescue heat stress transcription factor A2c (FaHSFA2c). Exogenous ABA (5 μM) application enhanced heat tolerance of tall fescue, as manifested by increased leaf photochemical efficiency and membrane stability under heat stress (37/32 °C, day/night). The expression levels of FaHSFA2c, several tall fescue HSPs (FaHSPs), and ABA-responsive transcriptional factors were up-regulated in plants treated with ABA. Deficiency of Arabidopsis heat stress transcription factor A2 (AtHSFA2) suppressed ABA-induction of AtHSPs expression and ABA-improved heat tolerance in Arabidopsis. These results suggested that HSFA2 plays an important role in ABA-mediated plant heat tolerance, and FaAREB3 and FaDREB2A may function as upstream trans-acting factors and regulate transcriptional activity of FaHSFA2c and the downstream FaHSPs, leading to improved heat tolerance.

Genomic and transcriptomic analyses of Agrobacterium tumefaciens S33 reveal the molecular mechanism of a novel hybrid nicotine-degrading pathway

Agrobacterium tumefaciens S33 is able to degrade nicotine via a novel hybrid of the pyridine and pyrrolidine pathways. It can be utilized to remove nicotine from tobacco wastes and transform nicotine into important functionalized pyridine precursors for some valuable drugs and insecticides. However, the molecular mechanism of the hybrid pathway is still not completely clear. Here we report the genome analysis of strain S33 and its transcriptomes grown in glucose-ammonium medium and nicotine medium. The complete gene cluster involved in nicotine catabolism was found to be located on a genomic island composed of genes functionally similar but not in sequences to those of the pyridine and pyrrolidine pathways, as well as genes encoding plasmid partitioning and replication initiation proteins, conjugal transfer proteins and transposases. This suggests that the evolution of this hybrid pathway is not a simple fusion of the genes involved in the two pathways, but the result of a complicated lateral gene transfer. In addition, other genes potentially involved in the hybrid pathway could include those responsible for substrate sensing and transport, transcription regulation and electron transfer during nicotine degradation. This study provides new insights into the molecular mechanism of the novel hybrid pathway for nicotine degradation.

Identification of TIFY/JAZ family genes in Solanum lycopersicum and their regulation in response to abiotic stresses

Plant phenotypic plasticity determines plant adaptation to changing environments and agricultural productivity. Phytohormones are essential plant signalling molecules regulating this plasticity through complex signalling networks. Jasmonates (JAs) are key phytohormones regulating many aspects of growth, development and defence responses. An important role of JAs in tolerance to abiotic stresses is also emerging. The expression of JAZ (JASMONATE-ZIM-DOMAIN PROTEIN) genes, encoding for the key repressors in the JA-pathway, is regulated by multiple abiotic stresses, suggesting a role for the JAZ proteins in response to these stresses. The JAZ proteins belong to the TIFY family, well described in many plant species. However, only the role of few tomato JAZ proteins in response to microbial infection has been analysed so far. Here, we identify the members of the tomato TIFY family, and characterize them phylogenetically. In addition, we analyse the transcriptional regulation of several SlJAZ in response to abiotic stresses and hormone treatments both in root and leaves to assess their specific expression in response to stresses. Most SlJAZ are JA-induced and responsive to one or more abiotic stresses, providing clues for functional analysis of JAZ genes in abiotic responses in tomato.

Transcriptome profiling of litchi leaves in response to low temperature reveals candidate regulatory genes and key metabolic events during floral induction

Background

Litchi (Litchi chinensis Sonn.) is an economically important evergreen fruit tree widely cultivated in subtropical areas. Low temperature is absolutely required for floral induction of litchi, but its molecular mechanism is not fully understood. Leaves of litchi played a key role during floral induction and could be the site of low temperature perception. Therefore, leaves were treated under different temperature (15 °C/25 °C), and high-throughput RNA sequencing (RNA-Seq) performed with leaf samples for the de novo assembly and digital gene expression (DGE) profiling analyses to investigate low temperature-induced gene expression changes.

Results

83,107 RNA-Seq unigenes were de novo assembled with a mean length of 1221 bp and approximately 61% of these unigenes (50,345) were annotated against public protein databases. Differentially-expressed genes (DEGs) under low temperature treatment in comparison with the control group were the main focus of our study. Hierarchical clustering analysis arranged 2755 DEGs into eight groups with three significant expression clusters (p-value ≤ 0.05) during floral induction. With the increasing contents of sugars and starch, the expression of genes involved in metabolism of sugars increased dramatically after low temperature induction. One FT gene (Unigene0025396, LcFT1) which produces a protein called ‘florigen’ was also detected among DEGs of litchi. LcFT1 exhibited an apparent specific tissue and its expression was highly increased after low temperature induction, GUS staining results also showed GUS activity driven by LcFT1 gene promoter can be induced by low temperature, which indicated LcFT1 probably played a pivotal role in the floral induction of litchi under low temperature.

Conclusions

Our study provides a global survey of transcriptomes to better understand the molecular mechanisms underlying changes of leaves in response to low temperature induction in litchi. The analyses of transcriptome profiles and physiological indicators will help us study the complicated metabolism of floral induction in the subtropic evergreen plants.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3747-x) contains supplementary material, which is available to authorized users.

De novo sequencing and comparative transcriptome analysis of adventitious root development induced by exogenous indole-3-butyric acid in cuttings of tetraploid black locust

Background

Indole-3-butyric acid (IBA) is applied to the cuttings of various plant species to induce formation of adventitious roots (ARs) in commercial settings. Tetraploid black locust is an attractive ornamental tree that is drought resistant, sand tolerant, can prevent sand erosion and has various commercial uses. To further elucidate the mechanisms of AR formation, we used Illumina sequencing to analyze transcriptome dynamics and differential gene expression at four developmental stages in control (CK) and IBA-treated groups.

Results

The short reads were assembled into 127,038 unitranscripts and 101,209 unigenes, with average lengths of 986 and 852 bp. In total, 10,181 and 14,924 differentially expressed genes (DEGs) were detected in the CK and IBA-treated groups, respectively. Comparison of the four consecutive developmental stages showed that 282 and 260 DEGs were shared between IBA-treated and CK, suggesting that IBA treatment increased the number of DEGs. We observed 1,721 up-regulated and 849 down-regulated genes in CI vs. II, 849 up-regulated and 836 down-regulated genes in CC vs. IC, 881 up-regulated and 631 down-regulated genes in CRP vs. IRP, and 5,626 up-regulated and 4,932 down-regulated genes in CAR vs. IAR, of which 25 up-regulated DEGs were common to four pairs, and these DEGs were significantly up-regulated at AR. These results suggest that substantial changes in gene expression are associated with adventitious rooting. GO functional category analysis indicated that IBA significantly up- or down-regulated processes associated with regulation of transcription, transcription of DNA dependent, integral to membrane and ATP binding during the development process. KEGG pathway enrichment indicated that glycolysis/gluconeogenesis, cysteine and methionine metabolism, photosynthesis, nucleotide sugar metabolism, and lysosome were the pathways most highly regulated by IBA. We identified a number of differentially regulated unigenes, including 12 methionine-related genes and 12 ethylene-related genes, associated with the KEGG pathway cysteine and methionine metabolism. The GO enrichment, pathway mapping, and gene expression profile analyses revealed molecular traits for root induction and initiation.

Conclusion

Our study presents a global view of the transcriptomic profiles of tetraploid black locust cuttings in response to IBA treatment and provides new insights into the fundamental mechanisms associated with auxin-induced adventitious rooting.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3554-4) contains supplementary material, which is available to authorized users.

Transcriptome Analyses Reveal Candidate Genes Potentially Involved in Al Stress Response in Alfalfa

Alfalfa is the most extensively cultivated forage legume, yet most alfalfa cultivars are not aluminum tolerant, and the molecular mechanisms underlying alfalfa responses to Al stress are largely unknown. In this study, we aimed to understand how alfalfa responds to Al stress by identifying and analyzing Al-stress-responsive genes in alfalfa roots at the whole-genome scale. The transcriptome changes in alfalfa roots under Al stress for 4, 8, or 24 h were analyzed using Illumina high-throughput sequencing platforms. A total of 2464 differentially expressed genes (DEGs) were identified, and most were up-regulated at early (4 h) and/or late (24 h) Al exposure time points rather than at the middle exposure time point (8 h). Metabolic pathway enrichment analysis demonstrated that the DEGs involved in ribosome, protein biosynthesis, and process, the citrate cycle, membrane transport, and hormonal regulation were preferentially enriched and regulated. Biosynthesis inhibition and signal transduction downstream of auxin- and ethylene-mediated signals occur during alfalfa responses to root growth inhibition. The internal Al detoxification mechanisms play important roles in alfalfa roots under Al stress. These findings provide valuable information for identifying and characterizing important components in the Al signaling network in alfalfa and enhance understanding of the molecular mechanisms underlying alfalfa responses to Al stress.

Comparative Transcriptomic Analysis of Grape Berry in Response to Root Restriction during Developmental Stages

Root restriction improved berry quality by being involved in diverse aspects of grapevine life. However, the molecular mechanism driving this process is not understood very well. In this study, the ‘Summer Black’ grape berry (Vitis vinifera × V. labrusca) under root restriction was investigated, which showed an increase of total soluble solids (TSS), color index of red grapes (CIRG) value, anthocyanins accumulation, total phenolics and total procyanidins contents during berry development compared with those in control berries. The transcriptomic changes induced by root restriction in ‘Summer Black’ grape over the course of berry development were analyzed by RNA-Seq method. A total of 29,971 genes were generated in ‘Summer Black’ grape berry during development, among which, 1606 genes were significantly responded to root restriction. Furthermore, 1264, 313, 141, 246 and 19 sequences were significantly changed at S1, S2, S3, S4 and S5 sample points, respectively. The gene (VIT_04s0023g02290) predicted as a salicylate O-methyltransferase was differentially expressed in all developmental stages. Gene Ontology (GO) enrichment showed that response to organic nitrogen, response to endogenous stimulus, flavonoid metabolic process, phenylpropanoid biosynthetic process and cell wall macromolecule metabolic process were the main significant differential categories. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment revealed plant–pathogen interaction, plant hormone signal transduction, flavone and flavonol biosynthesis, flavonoid biosynthesis and glucosinolate biosynthesis were the main significant differential pathways. The results of the present study provided a genetic base for the understanding of grape berry fruit quality improvement under root restriction.

Quantitative Resistance to Verticillium Wilt in Medicago truncatula Involves Eradication of the Fungus from Roots and Is Associated with Transcriptional Responses Related to Innate Immunity

Resistance mechanisms to Verticillium wilt are well-studied in tomato, cotton, and Arabidopsis, but much less in legume plants. Because legume plants establish nitrogen-fixing symbioses in their roots, resistance to root-attacking pathogens merits particular attention. The interaction between the soil-borne pathogen Verticillium alfalfae and the model legume Medicago truncatula was investigated using a resistant (A17) and a susceptible (F83005.5) line. As shown by histological analyses, colonization by the pathogen was initiated similarly in both lines. Later on, the resistant line A17 eliminated the fungus, whereas the susceptible F83005.5 became heavily colonized. Resistance in line A17 does not involve homologs of the well-characterized tomato Ve1 and V. dahliae Ave1 genes. A transcriptomic study of early root responses during initial colonization (i.e., until 24 h post-inoculation) similarly was performed. Compared to the susceptible line, line A17 displayed already a significantly higher basal expression of defense-related genes prior to inoculation, and responded to infection with up-regulation of only a small number of genes. Although fungal colonization was still low at this stage, the susceptible line F83005.5 exhibited a disorganized response involving a large number of genes from different functional classes. The involvement of distinct phytohormone signaling pathways in resistance as suggested by gene expression patterns was supported by experiments with plant hormone pretreatment before fungal inoculation. Gene co-expression network analysis highlighted five main modules in the resistant line, whereas no structured gene expression was found in the susceptible line. One module was particularly associated to the inoculation response in A17. It contains the majority of differentially expressed genes, genes associated with PAMP perception and hormone signaling, and transcription factors. An in silico analysis showed that a high number of these genes also respond to other soil-borne pathogens in M. truncatula, suggesting a core of transcriptional response to root pathogens. Taken together, the results suggest that resistance in M. truncatula line A17 might be due to innate immunity combining preformed defense and PAMP-triggered defense mechanisms, and putative involvement of abscisic acid.

Genomic profiling of exogenous abscisic acid-responsive microRNAs in tomato (Solanum lycopersicum)

BACKGROUND

Plant microRNAs (miRNAs) are involved in various biological pathways and stress responses as negative regulators at the posttranscriptional level. Abscisic acid (ABA) is a key signaling molecule that mediates plant stress response by activating many stress-related genes. Although some miRNAs in plants are previously identified to respond to ABA, a comprehensive profile of ABA-responsive miRNAs has not yet been elucidated.

RESULTS

Here, we identified miRNAs responding to exogenous application of ABA, and their predicted target genes in the model plant organism tomato (Solanum lycopersicum). Deep sequencing of small RNAs from ABA-treated and untreated tomatoes revealed that miRNAs can be up- or down-regulated upon treatment with ABA. A total of 1067 miRNAs were detected (including 365 known and 702 candidate novel miRNAs), of those, 416 miRNAs which had an abundance over two TPM (transcripts per million) were selected for differential expression analysis. We identified 269 (180 known and 89 novel) miRNAs that respond to exogenous ABA treatment with a change in expression level of |log2FC|≥0.25. 136 of these miRNAs (90 known and 46 novel) were expressed at significantly different levels |log2FC|≥1 between treatments. Furthermore, stem-loop RT-PCR was applied to validate the RNA-seq data. Target prediction and analysis of the corresponding ABA-responsive transcriptome data uncovered that differentially expressed miRNAs are involved in condition stress and pathogen resistance, growth and development. Among them, approximately 90 miRNAs were predicted to target transcription factors and pathogen resistance genes. Some miRNAs had functional overlap in biotic and abiotic stress. Most of these miRNAs were down-regulated following exposure to exogenous ABA, while their related target genes were inversely up-regulated, which is consistent with their negative regulatory role in gene expression.

CONCLUSIONS

Exogenous ABA application influences the composition and expression level of tomato miRNAs. ABA mainly down-regulates miRNAs that their target genes involve in abiotic stress adaption and disease resistance. ABA might increase expression of stress-related genes via miRNA-mediated posttranscriptional regulation, and our results indicate that ABA treatment has the potential to improve both abiotic stress tolerance and pathogen resistance. This study presents a comprehensive profile of ABA-regulated miRNAs in the tomato, and provides a robust database for further investigation of ABA regulatory mechanisms.

Salinity stress induces the production of 2-(2-phenylethyl)chromones and regulates novel classes of responsive genes involved in signal transduction in Aquilaria sinensis calli

BACKGROUND

Agarwood, is a resinous portion derived from Aquilaria sinensis, has been widely used in traditional medicine and incense. 2-(2-phenylethyl)chromones are principal components responsible for the quality of agarwood. However, the molecular basis of 2-(2-phenylethyl)chromones biosynthesis and regulation remains almost unknown. Our research indicated that salt stress induced production of several of 2-(2-phenylethyl)chromones in A. sinensis calli. Transcriptome analysis of A. sinensis calli treated with NaCl is required to further facilitate the multiple signal pathways in response to salt stress and to understand the mechanism of 2-(2-phenylethyl)chromones biosynthesis.

RESULTS

Forty one 2-(2-phenylethyl)chromones were identified from NaCl-treated A. sinensis calli. 93 041 unigenes with an average length of 1562 nt were generated from the control and salt-treated calli by Illmunina sequencing after assembly, and the unigenes were annotated by comparing with the public databases including NR, Swiss-Prot, KEGG, COG, and GO database. In total, 18 069 differentially expressed transcripts were identified by the transcriptome comparisons on the control calli and calli induced by 24 h or 120 h salinity stress. Numerous genes involved in signal transduction pathways including the genes responsible for hormone signal transduction, receptor-like kinases, MAPK cascades, Ca(2+) signal transduction, and transcription factors showed clear differences between the control calli and NaCl-treated calli. Furthermore, our data suggested that the genes annotated as chalcone synthases and O-methyltransferases may contribute to the biosynthesis of 2-(2-phenylethyl)chromones.

CONCLUSIONS

Salinity stress could induce the production of 41 2-(2-phenylethyl)chromones in A. sinensis calli. We conducted the first deep-sequencing transcriptome profiling of A. sinensis under salt stress and observed a large number of differentially expressed genes in response to salinity stress. Moreover, salt stress induced dynamic changes in transcript abundance for novel classes of responsive genes involved in signal transduction, including the genes responsible for hormone signal transduction, receptor-like kinases, MAPK cascades, Ca(2+) signal transduction, and transcription factors. This study will aid in selecting the target genes to genetically regulate A. sinensis salt-stress signal transduction and elucidating the biosynthesis of 2-(2-phenylethyl)chromones under salinity stress.

Comprehensive Analysis of ABA Effects on Ethylene Biosynthesis and Signaling during Tomato Fruit Ripening

ABA has been widely acknowledged to regulate ethylene biosynthesis and signaling during fruit ripening, but the molecular mechanism underlying the interaction between these two hormones are largely unexplored. In the present study, exogenous ABA treatment obviously promoted fruit ripening as well as ethylene emission, whereas NDGA (Nordihydroguaiaretic acid, an inhibitor of ABA biosynthesis) application showed the opposite biological effects. Combined RNA-seq with time-course RT-PCR analysis, our study not only helped to illustrate how ABA regulated itself at the transcription level, but also revealed that ABA can facilitate ethylene production and response probably by regulating some crucial genes such as LeACS4, LeACO1, GR and LeETR6. In addition, investigation on the fruits treated with 1-MCP immediately after ABA exposure revealed that ethylene might be essential for the induction of ABA biosynthesis and signaling at the onset of fruit ripening. Furthermore, some specific transcription factors (TFs) known as regulators of ethylene synthesis and sensibility (e.g. MADS-RIN, TAGL1, CNR and NOR) were also observed to be ABA responsive, which implied that ABA influenced ethylene action possibly through the regulation of these TFs expression. Our comprehensive physiological and molecular-level analysis shed light on the mechanism of cross-talk between ABA and ethylene during the process of tomato fruit ripening.

Tomato MBD5, a methyl CpG binding domain protein, physically interacting with UV-damaged DNA binding protein-1, functions in multiple processes

In tomato (Solanum lycopersicum), high pigment mutations (hp-1 and hp-2) were mapped to genes encoding UV-damaged DNA binding protein 1 (DDB1) and de-etiolated-1 (DET1), respectively. Here we characterized a tomato methyl-CpG-binding domain protein SlMBD5 identified by yeast two-hybrid screening using SlDDB1 as a bait. Yeast two-hybrid assay demonstrated that the physical interaction of SlMBD5 with SlDDB1 is mediated by the C-termini of SlMBD5 and the β-propeller-C (BPC) of SlDDB1. Co-immunoprecipitation analyses revealed that SlMBD5 associates with SlDDB1-interacting partners including SlDET1, SlCUL4, SlRBX1a and SlRBX1b in vivo. SlMBD5 was shown to target to nucleus and dimerizes via its MBD motif. Electrophoresis mobility shift analysis suggested that the MBD of SlMBD5 specifically binds to methylated CpG dinucleotides but not to methylated CpHpG or CpHpH dinucleotides. SlMBD5 expressed in protoplast is capable of activating transcription of CG islands, whereas CUL4/DDB1 antagonizes this effect. Overexpressing SlMBD5 resulted in diverse developmental alterations including darker green fruits with increased plastid level and elevated pigmentation, as well as enhanced expression of SlGLK2, a key regulator of plastid biogenesis. Taken together, we hypothesize that the physical interaction of SlMBD5 with the CUL4-DDB1-DET1 complex component may affect its binding activity to methylated DNA and subsequently attenuate its transcription activation of downstream genes.

Transcriptomic Analysis of Grapevine (cv. Summer Black) Leaf, Using the Illumina Platform

Proceeding to illumina sequencing, determining RNA integrity numbers for poly RNA were separated from each of the four developmental stages of cv. Summer Black leaves by using Illumina HiSeq^™ 2000. The sums of 272,941,656 reads were generated from vitis vinifera leaf at four different developmental stages, with more than 27 billion nucleotides of the sequence data. At each growth stage, RNA samples were indexed through unique nucleic acid identifiers and sequenced. KEGG annotation results depicted that the highest number of transcripts in 2,963 (2Avs4A) followed by 1Avs4A (2,920), and 3Avs4A (2,294) out of 15,614 (71%) transcripts were recorded. In comparison, a total of 1,532 transcripts were annotated in GOs, including Cellular component, with the highest number in “Cell part” 251 out of 353 transcripts (71.1%), followed by intracellular organelle 163 out of 353 transcripts (46.2%), while in molecular function and metabolic process 375 out of 525 (71.4%) transcripts, multicellular organism process 40 out of 525 (7.6%) transcripts in biological process were most common in 1Avs2A. While in case of 1Avs3A, cell part 476 out of 662 transcripts (71.9%), and membrane-bounded organelle 263 out of 662 transcripts (39.7%) were recorded in Cellular component. In the grapevine transcriptome, during the initial stages of leaf development 1Avs2A showed single transcript was down-regulated and none of them were up-regulated. While in comparison of 1A to 3A showed one up-regulated (photosystem II reaction center protein C) and one down regulated (conserved gene of unknown function) transcripts, during the hormone regulating pathway namely SAUR-like auxin-responsive protein family having 2 up-regulated and 7 down-regulated transcripts, phytochrome-associated protein showed 1 up-regulated and 9 down-regulated transcripts, whereas genes associated with the Leucine-rich repeat protein kinase family protein showed 7 up-regulated and 1 down-regulated transcript, meanwhile Auxin Resistant 2 has single up-regulated transcript in second developmental stage, although 3 were down-regulated at lateral growth stages (3A and 4A). In the present study, 489 secondary metabolic pathways related genes were identified during leaf growth, which mainly includes alkaloid (40), anthocyanins (21), Diterpenoid (144), Monoterpenoid (90) and Flavonoids (93). Quantitative real-time PCR was applied to validate 10 differentially expressed transcripts patterns from flower, leaf and fruit metabolic pathways at different growth stages.