LncRNA regulates tomato fruit cracking by coordinating gene expression via a  hormone-redox-cell wall network

doi:10.21203/rs.2.14406/v5

Download PDF

Research article

LncRNA regulates tomato fruit cracking by coordinating gene expression via a hormone-redox-cell wall network

https://doi.org/10.21203/rs.2.14406/v5

This work is licensed under a CC BY 4.0 License

Journal Publication

published 15 Apr, 2020

Read the published version in BMC Plant Biology →

You are reading this older preprint version

Read the latest preprint version →

Background Fruit cracking occurs easily under unsuitable environmental conditions and is one of the main types of damage that occurs in fruit production. It is widely accepted that plants have developed defence mechanisms and regulatory networks that respond to abiotic stress, which involves perceiving, integrating and responding to stress signals by modulating the expression of related genes. Fruit cracking is also a physiological disease caused by abiotic stress. It has been reported that a single or several genes may regulate fruit cracking. However, almost none of these reports have involved cracking regulatory networks.

Results Here, RNA expression in 0 h, 8 h and 30 h saturated irrigation-treated fruits from two contrasting tomato genotypes, ‘LA1698’ (cracking-resistant, CR) and ‘LA2683’ (cracking-susceptible, CS), was analysed by mRNA and lncRNA sequencing. The GO pathways of the differentially expressed mRNAs were mainly enriched in the ‘hormone metabolic process’, ‘cell wall organization’, ‘oxidoreductase activity’ and ‘catalytic activity’ categories. According to the gene expression analysis, significantly differentially expressed genes included Solyc02g080530.3 ( Peroxide, POD ), Solyc01g008710.3 ( Mannan endo-1,4-beta-mannosidase, MAN ), Solyc08g077910.3 ( Expanded, EXP ), Solyc09g075330.3 ( Pectinesterase , PE ), Solyc07g055990.3 ( Xyloglucan endotransglucosylase-hydrolase 7, XTH7 ), Solyc12g011030.2 ( X yloglucan endotransglucosylase-hydrolase 9 , XTH9 ), Solyc10g080210.2 ( Polygalacturonase-2, PG2 ), Solyc08g081010.2 ( Gamma-glutamylcysteine synthetase, gamma-GCS ), Solyc09g008720.2 ( Ethylene receptor , ER ), Solyc11g042560.2 ( Ethylene-responsive transcription factor 4, ERF4 ) etc. In addition, the lncRNAs (XLOC_16662 and XLOC_033910, etc) regulated the expression of their neighbouring genes, and genes related to tomato cracking were selected to construct a lncRNA-mRNA network influencing tomato cracking.

Conclusions This study provides insight into the responsive network for water-induced cracking in tomato fruit. Specifically, lncRNAs regulate the hormone-redox-cell wall network, including plant hormone (auxin, ethylene) and ROS (H 2 O 2 ) signal transduction and many cell wall-related mRNAs ( EXP, PG, XTH ), as well as some lncRNAs ( XLOC_16662 and XLOC_033910, etc.).

Plant Molecular Biology and Genetics

Tomato

LncRNA

mRNA

Transcriptome

Network

Fruit cracking

Fruit cracking, one of the main disorders in fruit production, can easily cause adverse impacts in fruit marketability such as reducing fruit quality due to a poor appearance, decreasing shelf life, and even making the fruit unmarketable because of fungal infection [1].

Fruit cracking occurs easily under unsuitable environmental conditions. For instance, under abiotic stress caused by dry to very wet conditions, there will be a rapid flow into the fruit, and if the skin loses strength and elasticity due to factors such as maturation, cracking is most likely to occur [2]. It is widely accepted that plants have developed defence mechanisms and regulatory networks to respond to abiotic stress, which involve perceiving, integrating and responding to stress signals by modulating the expression of related genes [3-8]. Fruit cracking is also a kind of physiological disease caused by abiotic stress. Is there a regulatory network involved in fruit cracking?

Since the 1930s, researchers have performed many theoretical and practical studies on cracking [2, 9, 10]. Cracking is the result of a combination of internal and external factors. The internal factors are the fruit's own characteristics (fruit size, shape, firmness, deposition of cutin, wax, strength of the pericarp, arrangement of cells in the pericarp, quantity and status of stomata, accumulation of osmoregulatory substances such as soluble sugars, growth stage of the fruit, etc.), while the external factors mainly include environmental factors (humidity, light, temperature, wind, etc.) and cultivation management measures (irrigation, mineral nutrition, plant regulation, etc.) [10-13]. Cortes [14] comprehensively analysed 62 genotypes and found that the correlation of cracked fruit with heredity was significantly greater than that with the environment, indicating that the cracking characteristic can show stable heritability and be regulated by certain genes.

Notably, cell wall components and modifications appear to be correlated with the strength of the skin and fruit cracking [15-18]. As ripening proceeds, cell wall degradation gradually occurs, and the fruit cracking rate increases [15-16]. The cell wall is composed of a cellulose-hemicellulose (Cel-Hem) network and pectin, which is essential to maintain the mechanical strength of the cell wall. As the fruit matures, enzymes and proteins that degrade the polysaccharide components of the cell wall are produced, such as polygalacturonase (PG), extended protein (EXP), pectin methylesterase (PME), pectate lyases (PL), pectinase, β-galactosidase. (β-gal) and cellulase (Cx) [19-22]. The synergistic action of these enzymes leads to the degradation of cell wall polysaccharides and softening of the mature fruit peel [23]. Previous research has shown that genes such as EXP, PG, β-gal and XET are associated with fruit cracking [24-29]. Inhibition of β-gal gene expression increases the rate of fruit cracking [27]. In tomato, inhibition of LePG expression slightly reduces the rate of fruit cracking [28]. Simultaneous suppression of SlPG and SlEXP1 expression in ripening fruits reduces cell wall disassembly and thereby reduces the fruit cracking rate by approximately 12% [30]. It is not a single gene but many genes working together that regulate fruit cracking [30-33]. It remains unclear whether there are other genes related to fruit cracking and which one is the major gene.

NcRNAs (Non-coding RNA) are involved in a lot of life processes, such as cell growth, differentiation, proliferation, and apoptosis [34-38]. In contrast to the approximately 2% of protein-coding genes, more than 90% of genes do not have the ability to encode proteins and are transcribed into ncRNAs [34]. These ncRNAs were originally thought to represent "expression noise" or "expression waste", but they have now been proven to be strictly regulated to play important roles in the biological processes of organisms and exhibit extremely complex biological functions [36-38]. While much of the published work on ncRNAs has been conducted in humans and animals, the studies on plants are limited to certain model plants, such as Arabidopsis, maize, and wheat [39-40]. Xin [39] identified 125 stress-related lncRNAs (long non-coding RNA) in wheat, among which 71 responded to powdery mildew, and 77 responded to heat stress. Swiezewski [40] discovered that cold-induced long antisense intragenic RNA (COOLAIR) is involved in the vernalization process and regulates the expression of the plant flowering suppressor “gene flowering locus C” (FLC). Wang [41] reported the expression and evolution of lncRNAs in Solanaceae. Cui’s [42] results provided insights into the WRKY1−lncRNA33732−RBOH module involved in the regulation of H₂O₂ accumulation and resistance to pathogens in tomato. Wang [43] identified several lnRNAs that are involved in TYLCV infection by virus-induced gene silencing (VIGS) and genome-wide analysis. However, none of these studies have focused on lncRNAs and fruit cracking. Is it possible that lncRNAs play important roles in fruit cracking too?

This study aimed to obtain a global view of the transcriptional regulation (mRNAs and lncRNAs) of fruit cracking induced by irrigation in tomato. Differentially expressed mRNAs and lncRNAs related to fruit cracking were identified through transcriptome profiling and bioinformatic analysis. Finally, we determined a lncRNA-regulated hormone-redox-cell wall network for water-induced cracking in tomato. The findings reported here can increase our understanding of the transcriptional regulatory mechanisms of fruit cracking.

RNA sequencing and identification of lncRNA and mRNA

In total, we obtained 0.81 to 1.14 billion raw reads and 0.79 to 1.14 billion clean reads from CR and CS tomatoes at various time points (0 h, 8 h and 30 h of saturated irrigation treatment) (Table 1). Through genomic comparison, Cufflinks splicing, and CPC2 and PFAM analysis, we identified 1 annotated lncRNA, 2508 putative lncRNAs , 33784 annotated mRNAs and 409 novel mRNAs ( Additional file 1).

Feature analysis of lncRNAs and identification of lncRNA-mRNA pairs

The average length of the obtained lncRNAs was 1470 nt, which was similar to that of the mRNAs (1221 nt); the average number of exons and average ORF length of the identified lncRNAs were 2.6 and 88.5 bp, which were much lower values than those for the mRNAs (4.7 and 347 bp) (Fig. 1), consistent with previous studies [44, 45]. At the same time, we used phyloP to separately score the lncRNAs and mRNAs, and the sequence conservation of the lncRNAs was lower than that of mRNAs, which was consistent with previous studies [46]. We identified 21048 lncRNA-mRNA pairs with target relationships upstream and downstream of 2508 lncRNAs (Additional file 2).

Differential expression analysis

Differentially expressed mRNAs and lncRNAs were analysed in CR and CS tomatoes by using edge R software (Fig. 2) and the number of differentially expressed genes was listed in the Additional file 3. mRNAs and lncRNAs with a Q-value<0.05 and |log2 fold-change|> 1 were selected as differentially expressed genes.

Functional prediction of DEGs

To investigate the trends in gene functions and enrichment for DEGs, we performed GO (Gene Ontology) analysis of the selected mRNAs (Fig. 3; Additional file 4). The results showed that DEGs in the CR tomato were involved in a series of biological processes, such as regulation of biological process, biological regulation and regulation of cellular process, as well as catalytic activity. For the CS tomato, DEGs were mainly involved in single-organism metabolic process, biological process and catalytic activity. Between the CR and CS tomatoes, before irrigation treatment (0 h), DEGs were significantly enriched in oxidoreductase activity; after 8 h of irrigation treatment, there were some DEGs enriched in fruit ripening, anatomical structure maturation, and ageing; after 30 h of irrigation treatment, the number of DEGs enriched in the catalytic category was the highest, followed by the single-organism metabolic process and oxidation-reduction process categories, and cell components.

To better understand the function of the DEGs, significantly enriched KEGG pathways were analysed (Fig. 4; Additional file 5). The results showed that the DEGs were mainly enriched in the ‘biosynthesis of secondary metabolites’, ‘cysteine and methionine metabolism’, ‘metabolic pathways’, ‘plant-pathogen interaction’, ‘photosynthesis-antenna protein’, ‘photosynthesis’, ‘histidine metabolism’ and ‘circadian rhythm-plant’ categories.

DEGs (cell wall, redox, hormone related) regulate tomato fruit cracking

According to the gene expression analysis, 16 significantly differentially expressed genes (Additional file 6) were predicted to be related to fruit cracking in tomato, such as Solyc02g080530.3 (Peroxide, POD), Solyc01g008710.3 (Mannan endo-1,4-beta-mannosidase, MAN), Solyc08g077910.3 (Expanded, EXP), Solyc09g075330.3 (Pectinesterase, PE), Solyc07g055990.3 (Xyloglucan endotransglucosylase-hydrolase 7, XTH7), Solyc12g011030.2 (Xyloglucan endotransglucosylase-hydrolase 9, XTH9), Solyc10g080210.2 (Polygalacturonase-2, PG2), Solyc08g081010.2 (Gamma-glutamylcysteine synthetase, gamma-GCS), Solyc09g008720.2 (Ethylene receptor, ER), Solyc11g042560.2 (Ethylene-responsive transcription factor 4, ERF4) etc. (Table. 2 [47-70] ). Hierarchical clustering analysis showed that the expression trends or levels of these genes in the two varieties were completely different after the irrigation treatment (Fig. 5a). For instance, the expression of XTH7, XTH9, PE and POD in the CR tomato showed a downward trend, while the expression in the CS tomato presented an upward trend. These genes play important roles in cell wall loosing and expansion. As disassembly of the fruit cell wall can influence fruit cracking [30], These plant cell-wall loosing genes may also play a key regulatory role in tomato fruit cracking. At the same time, we used Tomato Gene Expression Atlas (http://tea.solgenomics.net/expression_viewer/input) to verify the gene expression, and found that most differentially expressed genes in this experiment were expressed in tomato pericarp in red ripe stage. Among them, GCS, MAN and PG, the antioxidative genes and cell-wall degrading enzyme-associated genes showed the highest expression. Whereas high mRNA levels were never present in ERF. These might be because that the ERF is an upstream regulator [71], so it never presented a high expression in red ripe stage. The differences in gene expression in this experiment and Tomato Gene Expression Atlas might be due to different varieties, treatment and detection standards (Fig. 5b).

Finally, we mapped a pathway diagram (Fig. 5c) of fruit cracking based on these differentially expressed lncRNAs, mRNAs and previous studies [71-77]. Within this pathway, ERF, POD, PG and PE play important roles. Previous researches suggests that ethylene influences fruit development and ripening (regulating cell wall-related PG and EXP gene expression) [71] and promotes programmed cell death of epithelial cells under ROS signalling [72]. Li et al. [73] showed that ARFs represent a point of cross-talk between ethylene and auxin signalling. Furthermore, auxin induces the production of ROS, and H₂O₂ decomposes polymers at the cell wall by producing ·OH [74]. Programmed cell death leads to a reduction in or loss of permeability of the plasma membrane, which in turn influences fruit cell activity, water absorption and cracking [75]. Simultaneously, the increase of auxin can promote the accumulation of H₂O₂ and the elongation of cells [76]. Furthermore, Rayle and Cleland [77] proposed the acid growth theory indicating that hydrogen ions may exert a purely chemical or physical effect, such as cleavage of acid-labile bonds on the wall, or they may activate normal enzymatic processes directly or indirectly, potentially leading to wall loosening. We analyzed the regulatory element in the promoter sequence of PG, PE, EXP and XTH7 and found that there are Ethylene-responsive element and auxin-responsive element. Based on these findings, we speculate that the regulatory network of fruit cracking, especially the coexpression of cell wall-, redox-, and hormone-related mRNAs and their corresponding lncRNAs, influences fruit cracking.

qRT-PCR validation of DEGs

Genes showing upregulated and the downregulated expression were randomly selected from the DEGs for qRT-PCR verification. The results of qRT-PCR revealed that most of these mRNAs shared similar expression tendencies to those indicated by the mRNA-Seq data, which can validate the reliability of our sequence data and our research results from the present study (Additional file 7 ). The expression levels detected by the two methods were slightly different, which might have been due to the different detection ranges and sensitivities of the two detection methods. The comparison of the relative expression measured by qRT-PCR and RNA-seq was showed in Additional file 8, R²＞0.7 confirmed the reliability of the RNA-Seq analysis results.

Tomato is one of the most popular commercial vegetables [78], however, its fruit shows high susceptibility to cracking [16, 31]. Cracks can occur throughout the fruit development stage during the ripening and post-harvest period [79, 80], which may cause serious economic losses. Different hypotheses have been presented to explain the occurrence of tomato fruit cracking. Previous studies have shown that rapid fruit swelling and fruit cracking are closely related [81]. Irregular temperatures or watering, especially a shift from a lower temperature to a much higher temperature or from extremely dry to very humid conditions, will lead to rapid swelling. The pressure of the rapidly expanding pulp on the peel may lead to fruit cracking [82]. Cell senescence and apoptosis also influence skin strength and water absorption, which can in turn affect fruit cracking [10]. In addition, a large differences between day and night temperatures can lead to the accumulation of carbohydrates [83]. Fruits with high levels of carbohydrates absorb more water, grow much faster and are more likely to crack [84]. In general, fruit cracking is a complex problem involving a mixture of genetics and the environment. Previous studies have suggested it is not a single gene but many genes that work together to regulate fruit cracking [31, 32].

Cell wall polysaccharide metabolic

The DEG Solyc08g077910.3 encodes an Expansin-like protein that breaks down the hydrogen bonds between molecules in the cell wall macromolecular network to promote the depolymerization of the network, which can lead to relaxation of the cell wall [52]. In this experiment, the expression level of Solyc08g077910.3 was increased significantly after 8 h of irrigation (log2 fold-change=7.13395) in CS tomato. The increased expression of the expansin-like gene can relax the cell wall and may influence fruit cracking.

Solyc07g055990.3 and Solyc12g011030.2 encode xyloglucan endotransglucosylase/hydrolases 7 and 9, respectively, which mediate the cleavage and polymerization of β-1,4-xyloglucan in the primary cell wall. Xyloglucan is usually fused to the cell wall, and its oligosaccharides determine tissue tension [47]. Jan [48] found that OsXTH8 is involved in the cell wall modification process in rice and is highly expressed in the vascular bundle of the sheath and the young roots, in which the cells are rapidly elongated and differentiated. In addition, it can respond to gibberellin. He [49] found that OsXTH5, OsXTH19, OsXTH20, OsXTH24 and OsXTH28 play important roles in the elongation of rice peduncles and can respond to drought stress. These studies indicate that the OsXTH gene family plays an important role in the regulation of the structural function of rice cell walls. In this experiment, the expression levels of Solyc12g011030.2 and Solyc07g055990.3 in CS tomato showed an upward trend, while they showed a downward trend in the CR tomato (Fig. 5a). Simultaneously, the expression level in CS tomato was significantly higher than that in CR tomato. This illustrates that the CR tomato may exhibit a greater osmotic stress resistance ability with downregulation of the XTH gene that can strengthen the cell wall upon encountering water stress.

The DEG Solyc10g080210.2 (Polygalacturonase-2) can remove the methyl group from polygalacturonic acid; during tomato maturation, the degree of methylation decreases from 90% in the green ripen period to 35% in the red ripen period [50], which accelerates the degradation of the cell wall. In an antisense PaPG1 transgenic study of strawberry, the expression level of PG was significantly inhibited, and the degree of fruit softening was significantly delayed [51].

Redox processes

Previous studies have shown that peroxidase in the cell wall leads to cell wall sclerosis by causing cross-linking of cell wall components, thereby inhibiting cell elongation [56-58]. Peroxidase can also directly regulate plant cell elongation by controlling H₂O₂ levels [59]. Solyc02g080530.3 encodes peroxide, whose levels are significantly higher in CS tomato than in CR tomato. The expression of these genes in CS tomato fruits may increase cell wall hardness and hinder the elongation of the cell wall, which will lead to fruit cracking when water absorption swelling occurs. Solyc01g081250.3 encodes glutathione-S-transferase (GST). The GST superfamily enzymes exhibit multiple functions in plants. They are not only involved in primary metabolism and secondary metabolism [60], but they can also protect plants from oxidative damage and heterogeneous substances [61-63]. According to the data analysis, the gene expression of Solyc01g081250.3 in the CR tomato was significantly higher than that in the CS tomato after 0 h, 8 h and 30 h of irrigation treatment. Higher expression of GST in CR tomato can better maintain cell vigour and be beneficial to tomato fruits when coupled with water stress.

Hormone-related

Previous research has shown that hormones can regulate the expression of cell wall-related genes. Trainotti [64] studied the expression of 32 genes related to cell wall synthesis and degradation. Their research showed that the expression of these genes in fleshy fruits can be inhibited by ethylene, while ethylene promotes the expression of these genes during fruit ripening and softening. At the same time, ethylene inhibits and promotes dual regulatory effects on the formation of plant secondary metabolites [65-67]; TAPG1, encoding a cell wall-degrading enzyme, can be induced by ethylene at the transcriptional level in tomato [68]; Rose [69] showed that ethylene regulates LeEXP1, which is specifically expressed only during fruit ripening. The pathway of ethylene biosynthesis in plants is the methionine cycle[85-87]. In this study, KEGG functional analysis of DEGs revealed significant enrichment in the methionine metabolic pathway. Solyc11g042560.2 encodes an ethylene receptor, while Solyc09g008720.2 encodes an ethylene-responsive transcription factor, and their expression levels are significantly upregulated after irrigation and are higher in CS tomato than in CR tomato.

LncRNAs regulate tomato fruit cracking by coordinating gene expression in the hormone-redox-cell wall network

Liu’s research suggests that plants have gradually developed complex signalling pathways to cope with adverse environmental stimuli [88]. That is, plants perceive different stress signals from the circumstances in which they occur and then integrate these signals and respond to these different stresses by modulating the expression of related genes. Is it possible that cracking is also regulated by a complex network?

LncRNAs play important roles in epigenetic regulation, cell cycle regulation and many other activities. Here, we identified several lncRNAs that are involved in fruit cracking. Most lncRNAs are not annotated, and we do not know their functions. To predict the functions of these lncRNAs, we performed functional analysis of lncRNA-targeted mRNAs and constructed an lncRNA-mRNA network (Fig. 6; Additional file 9). The results showed that the mRNAs in the network (Fig. 6a) were mainly enriched in the ‘oxidation-reduction process’, ‘oxidoreductase activity’, ‘hormone metabolic process’, ‘response to hormone stimulus’, ‘catalytic activity’, ‘cell wall organization’ and ‘external encapsulating structure’ categories. We classified the target genes into four categories (cell wall polysaccharide metabolism, oxidation-reduction processes, hormones and others) based on their functions and amounts.

Some lncRNAs specifically target functional mRNAs, and we can assume that the lncRNAs perform similar functions to their target mRNAs. For example, many of the target genes of XLOC_033910, XLOC_007053, and XLOC_008464 (Fig. 6b) are enriched in categories such as ‘dioxygenase activity’, ‘oxidation-reduction process’ and ‘oxidoreductase activity’, so we predicted their gene function as ‘redox regulation’.

Some lncRNAs are targets of significantly differentially expressed mRNAs with various functions. For example, for XLOC_033910 (Fig. 6c), the target genes (Solyc06g060070.3,Solyc07g014670.3, Solyc05g055590.2, etc.) are enriched in the 'oxidoreductase activity' and 'dioxygenase activity' terms. XLOC_008464 also has other target genes, such as Solyc03g123550.1 ('response to hormone stimulus'), and Solyc11g005750.2 ('cell wall'). Previous studies [58-80] have shown that redox, hormone and cell wall terms are all very important factors that can influence fruit cracking, so we speculate that lncRNA XLOC_7053 may plays an important role in regulating tomato fruit cracking.

Liao’s research strongly supported Ethylene-responsive factor 4 (ERF4) was associated with rind hardness and cracking resistance in watermelon fruits [89]. In our study, DEGs also included ethylene, as well as cell wall and redox related genes. Based on our research and previous research [68-77], we mapped a “hormone-redox-cell wall” pathway diagram including lncRNA associated with fruit cracking.

In this study, the key genes involved in the response to tomato cracking were identified by high-throughput sequencing, which has important significance for guiding the selection of new tomato varieties. We have also established an lncRNA-mRNA (hormone-redox-cell wall) network to learn about the precise regulation of fruit cracking by lncRNAs. To the best of our knowledge, this is the first discovery of the lncRNA-mRNA network involved in tomato fruit cracking.

Plant materials and sample collection

RNA expression in 0 h, 8 h and 30 h irrigation-treated fruits from two contrasting tomato genotypes, ‘LA1698’ (cracking-resistant, CR) and ‘LA2683’ (cracking-susceptible, CS), was analysed by mRNA and lncRNA sequencing. ‘LA2683’ and ‘LA1698’ (Fig. 7) were both introduced by the Tomato Genetics Resource Centre (TGRC, University of California, Davis). The fruit cracking rates of ‘LA2683’ and ‘LA1698’ are 77.53% and 20.17% (Additional file 10), respectively. Both lines were selected and self-pollinated for more than 6 generations. All the seedlings were grown in 72-plug trays on 18 March 2016. On 28 April 2016, they were transplanted to the same greenhouse of the Kunshan Yuye Leaf Vegetable Base (31°95ʹE, 119°16ʹN), Suzhou, Jiangsu Province, China. The climate of this area belongs to the north subtropical south monsoon climate zone, with four distinct seasons and plenty of rainfall. The average annual temperature is 15.7 ℃, and the annual average precipitation is 1094 mm according to China Weather Network (http://www.weather.com.cn/jiangsu). These two genotypes were planted side by side. Plant spacing followed a 30×50×100 cm pattern. Drip irrigation were adopted to guarantee the consistent of irrigation amount. Fertilization practices and pest control were those usually used by local growers. When fruits of the third cluster ripen, saturated irrigation based on previous researches [2, 16, 31] with slightly modificationwas adopted to induce fruit cracking. Specifically, plants were firstly irrigated using flood irrigation for 2 h. The field was essentially flooded with water which was allowed to totally soak into the soil. We then adopted drip irrigation to keep the soil saturated with water. After 0 h, 8 h and 30 h of irrigation treatment, for each genotype, treatment, one fruit were cut into four pieces from peduncle to blossom-end including exocarp and mesocarp without pulp. Each group has two repetitions. Samples of both genetypes were immediately frozen in liquid nitrogen and stored at -80℃. Sampling time is based on morphological observation. At 8 h, cracked fruits began to appear in cracking susceptible tomato ‘LA2683’ but not in cracking resistant tomato ‘LA1698’. At 30 h, the fruits of ‘LA2683’ severely cracked and fruits of ‘LA1698’ slightly cracked .

RNA-seq analysis

Twelve fruit samples (2 genotypes, 3 time points, each with two replications) were collected and sequenced by Novogene, Beijing, China. RNA was extracted from the tomato pericarp, and qualified RNA samples were used to construct a cDNA library. Transcriptome sequencing was carried out on the Illumina HiSeq 2500 platform. Sequences showing low quality, linker contamination or a high unknown base (N) content were filtered from the raw reads obtained after sequencing. The filtered high-quality clean data were aligned to the ITAG4.0 reference genome using Hisat2 [90], and the transcriptome was assembled by using Cufflinks [91].

Identification of mRNA and lncRNA

The transcripts that could be compared to known transcript data were identified as annotated mRNAs.

Then, the transcripts were screened according to the following criteria: (1) exon number ≥2 and (2) transcript length ≥200 bp. Additionally, (3) Transcripts that overlapped with the database annotation of the exon area according to Cuffcompare software were screened out. The lncRNAs overlapping with the exon region of a spliced transcript were included in the subsequent analysis. (4) The expression level of each transcript was calculated with Cuffquant, and transcripts with an FPKM≥0.5 were selected. (5) Two algorithms for evaluating protein-coding potential (CPC2 [92] and PFAM [93]) were used to predict the protein-coding potential of the remaining transcripts. Only when these two algorithms simultaneously indicated no protein-coding potential were the sequences considered to be predicted lncRNAs. Finally, the predicted lncRNAs were obtained. The transcripts showing the potential to encode proteins by “CPC2” and “PFAM” were identified as novel mRNAs (Fig. 8).

Analysis of mRNAs and lncRNAs

The mRNA and lncRNA abundance of the unigenes was normalized via the fragments per kilobase of exon model per million mapped reads (FPKM) approach. The log2 fold changes between two samples were tested statistically to determine whether the expression of an individual gene was significantly altered. A Q-value <0.05 and |log2 fold-change| > 1 was considered to indicate a differentially expressed gene (DEG). Analyses of mRNA and lncRNA expression trends in tomato fruit after irrigation were then carried out. To understand the function of differentially expressed mRNAs, these mRNAs were further subjected to GO (Gene Ontology) and KEGG (Kyoto Encyclopedia of Genes and Genomes) analysis by using GOseq [94] and KOBAS [95] software, respectively.

Feature analysis of lncRNAs and identification of lncRNA-mRNA pairs

Genomic characterization of the predicted lncRNAs was performed and compared with the mRNA results. The parameters for comparison to understand the genomic characteristics of the lncRNAs included the number of exons, ORF length, transcript nucleic acid length and sequence conservation between species.

LncRNAs mainly acts on their protein-coding target genes by cis-regulation or trans-regulation to achieve their regulatory functions. One of the functions of lncRNAs is cis-regulation of their neighbouring genes of the same allele; lncRNAs, which were located in less than 100 kb up/down stream of a gene, probably acted as cis-regulators. The up-stream lncRNAs having intersection of promoter or other cis-elements may regulate gene expression during or after transcription. The 3’UTR region or down-stream lncRNAs may perform other regulatory functions [96]. Another function of lncRNAs is the trans-regulation of co-expressed genes that are not adjacent to lncRNAs. For trans-regulation analysis, multiple repetitions are required for the result to be accurate and persuasive, samples with a value < 6 are not recommended [97, 98]. So in this study, lncRNAs by cis-regulation were analyzed. LncRNA within 100kb upstream and downstream of the mRNA was identified as the lncRNA-mRNA pair and the function of lncRNAs was predicted through functional enrichment analysis of the cis target gene.

qRT-PCR verification

Total RNA was extracted from tomato fruit by using the RNAprep Pure Plant Kit (Tiangen Biotech Co., Ltd. (Beijing, China)) following the manufacturer’s instructions. cDNA was generated using 5 µl of total RNA and abm’s 5X All-In-One MasterMix. qRT-PCR was performed in a CFX 96 Touch RT-PCR detection system (Bio-Rad, USA) with abm’s EvaGreen 2X qPCR MasterMix-Low ROX. Then, 9 DEGs were randomly selected from the DEGs to verify the RNA-seq results. Gene-specific primers were designed using Premier 5.0 [99] and actin as a reference gene [30]. The relative levels of gene expression were calculated using the 2^-ΔΔCTmethod [100] . The sequences of the primers are listed in Additional file 11.

CS: Cracking-susceptible

CR: Cracking-resistant

CS 0 h: Cracking-susceptible tomato fruit before irrigation (0 hours)

CS 8 h: Cracking-susceptible tomato fruit after 8 hours of irrigation

CS 30 h: Cracking-susceptible tomato fruit after 30 hours of irrigation

CR 0 h: Cracking-resistant tomato fruit before irrigation (0 hours)

CR 8 h: Cracking-resistant tomato fruit after 8 hours of irrigation

CR 30 h: Cracking-resistant tomato fruit after 30 hours of irrigation

DEGs: Differentially expressed genes

GO: Gene Ontology

KEGG: Kyoto Encyclopedia of Genes and Genomes

qRT-PCR: Quantitative real-time PCR

FPKM: Expected number of fragments per kilobase of transcript sequence per millions base pairs sequenced

Ethics approval and consent to participate

Not applicable.

Consent to publish

Not applicable.

Availability of data and materials

The datasets generated during and/or analysed during the current study are included in this published article (and its supplementary information files). All RNA sequencing data from this study are available in the NCBI sequence read archive (SRA) under accession numbers: SAMN14409497 , SAMN14409498 , SAMN14409499 , SAMN14409500 , SAMN14409501 , SAMN14409502 , SAMN14409503 , SAMN14409504 , SAMN14409505 , SAMN14409506 , SAMN14409507 and SAMN14409508 (https://www.ncbi.nlm.nih.gov/sra/PRJNA613785).

Competing interests

The authors declare that there are no competing interests in the reported research.

Funding

The work was supported by grants from the National Natural Science Foundation of China (No. 31701924), the Fundamental Research Funds for the Central Universities (No. KJQN201814, KYZZ201909) and Priority Academic Program Development of Jiangsu Higher Education Institutions. The funders provided the financial but played no role in the study design, analysis of the data and in writing the manuscript.

Authors' Contributions

JFL led and coordinated the project. XLZ and SMT collected the plant materials. YL and XLZ isolated the RNA. SMT and TYP conducted real-time quantitative PCR. XLZ, WZ and YQH conducted the bioinformatics analysis. XLZ, SMT, and JFL wrote the paper. JFL is the corresponding author and is responsible for all contact and correspondence. All authors have read and approved the final manuscript.

Acknowledgements

Not applicable.

1.Khadivikhub A, Zamani Z, Bouzari N. Evaluation of genetic diversity in some Iranian and foreign sweet cherry cultivars by using RAPD molecular markers and morphological traits. Horticulture Environment & Biotechnology. 2008;49(3):188-196.

Balbontín, Cristián, Ayala, et al. Cracking in sweet cherries：A comprehensive review from a physiological molecular and genomic perspective. Chilean journal of agricultural research. 2013;73(1):66-72.
Mukhopadhyay A, Vij S, Tyagi AK. Overexpression of a Zinc finger protein gene from rice confers tolerance to cold, dehydration, and salt stress in transgenic tobacco. Proceedings of the National Academy of Sciences of the USA. 2004;101: 6309-6314.
Kang M Y, Fokar M, Abdelmageed H, Allen RD. Arabidopsis SAP5 functions as a positive regulator of stress responses and exhibits E3 ubiquitin ligase activity. Plant Molecular Biology. 2011;75: 451-466.
Wu X, Shiroto Y, Kishitani S, Ito Y, Toriyama K. Enhanced heat and drought tolerance in transgenic rice seedlings overexpressing OsWRKY11 under the control of HSP101 promoter. Plant Cell Reports. 2009;28( 1) : 21-30.
Sakuma Y，Maruyama K，Osakabe Y，et al. Functional analysis of an Arabidopsis transcription factor，DREB2A，involved in drought responsive geneexpression. Plant Cell. 2006;18( 5) : 1292- 1309.
Yamaguchi-Shinozaki, K. and Shinozaki K. Organization of cis-acting regulatory elements in osmotic- and cold-stress-responsive promoters. Trends in plant science. 2005;10:88-94.
Liu JH，Peng T，Dai WS. Critical cis-acting elements and interacting transcription factors: key players associated with abiotic stress responses in plants. Plant molecular biology reporter. 2014; 32:303-317.
Sawada E. Studies on the cracking of cherries. I. The cause of cracking. Agriculture and Horticulture. 1931;6: 964-892.
Khadivi-Khub A. Physiological and genetic factors influencing fruit cracking. Acta Physiologiae Plantarum. 2015;37:1718.
Cristián B, Héctor A, Bastías RM, et al. Cracking in sweet cherries: A comprehensive review from a physiological, molecular, and genomic perspective. Chilean journal of agricultural research. 2013;73(1):66-72.
Beyer M, Hahn R, Peschel S, Harz M, Knoche M. Analysing fruit shape in sweet cherry (Prunus avium L.). Scientia Horticulturae. 2002;96:139-150.
Considine J, Brown K. Physical aspects of fruit growth. Theoretical analysis of distribution of surface growth forces in relation to cracking and splitting. Plant Physiology. 1981; 68:371-376.
Cortes C, Ayuso MC, Palomares G, Cuartero J, Nuez F. Relationship between radial and concentric cracking of tomato fruit. Scientia Horticulturae. 1983;21:323-328.
Domínguez E, Ferna´ndez MD, Herna´ndez JCL, et al. Tomato fruit continues growing while ripening, affecting cuticle properties and cracking. Physiol Plantarum. 2012;146(4):473-486.

16.Emmons CL W, Scott JW. Environmental and physiological effects on cuticle cracking in tomato. Journal of the American Society for Horticultural Science. 1997;122(6):797-801.

López-Casado G, Matas AJ, Domínguez E, Cuartero J, Heredia A. Biomechanics of isolated tomato (Solanum lycopersicumL.) fruit cuticles: the role of the cutin matrix and polysaccharides. Journal of Experimental Botany. 2007;58:3875-3883.
Domínguez E, España L, López-Casado G, Cuartero J, Heredia A. Biomechanics of isolated tomato (Solanum lycopersicum) fruit cuticles during ripening: the role of flavonoids. Functional Plant Biology . 2009;36:613–620.
Ruiz-May E, Rose JKC. Progress toward the tomato fruit cell wall proteome. Functional Plant Biology. 2013;29:7.
Carpita NC, Campbell M, Tierney M. Pectin: cell biology and prospects for functional analysis. Plant Cell Walls. 2001; 2:9-27.
Marin-Rodriguez MC . Pectate lyases, cell wall degradation and fruit softening. Journal of Experimental Botany. 2002;53(377):2115-2119.
Geitmann PA. Pectin and the role of the physical properties of the cell wall in pollen tube growth of Solanum chacoense. Planta. 2005;220(4):582-592.
Cantu D, Vicente A, Greve L, et al. The intersection between cell wall disassembly, ripening, and fruit susceptibility to Botrytis cinerea. Proceedings of The National Academy of Sciences of The United States of America. 2008;105(3):859-864.
Wang Y, Lu WJ, Li JG, Jiang Y. Differential expression of two expansion genes in developing fruit of cracking-susceptible and resistant litchi cultivars. Journal of the American Society for Horticultural Science. 2006;131(1):118-121.
Kasai S, Hayama H, Kashimura Y, Kudo S, Osanai Y. Relationship between fruit cracking and expression of the expansin gene MdEXPA3 in ‘Fuji’ apples (Malus domesticaBorkh.). Scientia Horticulturae. 2008;116(2):194-198.
Huang XM, Wang HC, Lu XJ, et al. Cell wall modifications in the pericarp of litchi (Litchi chinensis Sonn.) cultivars that differ in their resistance tocracking. Journal of Horticultural Science & Biotechnology. 2006;81(2):231-237.
Moctezuma E, Smith DL, Gross KC. Antisense suppression of a ß-galactosidase gene (TBG6) in tomato increases fruit cracking. Journal of Experimental Botany. 2003;54:2025-2033
Schuch W, Kanczler J, Robertson D, et al. Fruit quality characteristics of transgenic tomato fruit with altered polygalacturonase activity. Hortscience. 1991;26(12):1517-1520.
Lu W, Wang Y, Jiang Y, et al. Differential expression of litchi XET genes in relation to fruit growth. Plant Physiology and Bochemistry. 2006;44(11):707-713
Jiang FL, Lopez A, Jeon S, et al. Disassembly of the fruit cell wall by the ripening-associated polygalacturonase and expansin influences tomato cracking. Horticulture Research. 2019; 6 :17.
Peet M. Fruit cracking in tomato. Hort Technology. 1992; 2(2):216-223.
Cuartero J, Palomares G, Balasch S. Tomato fruit cracking under plastic-house and in the open air. II. General and specific combining abilities. Bmc Infectious Diseases. 1981; 15(1):1-13.
Yang ZE, Wu Z, Zhang C, Hu EM, Zhou R, Jiang FL. The composition of pericarp, cell aging, and changes in water absorption in two tomato genotypes: mechanism, factors, and potential role in fruit cracking. Acta Physiologiae Plantarum. 2016;38: 215.
Birney E, Stamato yannopoulos JA, Dutta A, et al. Indentification and analysis of function elements in 1% of the human genome by the Encode pilot project. Nature. 2007;447(7146):799-816.
Derrien T, Johnson R, Bussotti G, et al. The GENCODE v7 catalog of human long noncoding RNAs: Analysis of their gene structure, evolution, and expression. Genome Research. 2012;22(9):1775-1789.
Ponting CP, Belgard TG. Transcribed dark matter: meaning or myth? Human Molecular Genetics. 2010;19(R2):R162-8.
Huttenhofer A, Schattner P, Polacek N. Non-coding RNAs: hope or hype? Trends in Genetics. 2005; 21(5):289-97.

38.Wang JQ, Zhang YJ, Ren JM, et al. Advances in research on biological functions and significance of Long non-coding RNA. Life science. 2012;24(06):543-548.

Xin MM, Wang Y, Yao YY, et al. Identification and characterization of wheat long non-protein coding RNAs responsive to powdery mildew infection and heat stress by using microarray analysis and SBS sequencing. BMC Plant Biology. 2011;11: 61.
Swiezewski S, Liu F, Magusin A, Dean C. Cold-induced silencing by long antisense transcripts of an Arabidopsis Polycomb target. Nature. 2009;462(7274):799-802.

41.Wang X, Ai G, Zhang C, et al. Expression and diversification analysis reveals transposable elements play important roles in the origin of Lycopersicon‐specific lncRNAs in tomato. New Phytologist. 2016;209.

Cui J , Jiang N , Meng J , et al. LncRNA33732‐respiratory burst oxidase module associated with WRKY1 in tomato‐ Phytophthora infestans interactions. The Plant Journal. 2019;97(5):933-946.
Wang J , Yu W , Yang Y , et al. Genome-wide analysis of tomato long non-coding RNAs and identification as endogenous target mimic for microRNA in response to TYLCV infection. Scientific Reports. 2015, 5:16946.
Consortium TF, Carninci P, Kasukawa T, et al. The transcriptional landscape of the mammalian genome. Science. 2005; 309 (5740):1559-1563
Luo H , Yang H , Lin Y , et al. LncRNA and mRNA profiling during activation of tilapia macrophages by HSP70 and Streptococcus agalactiae antigen. Oncotarget. 2017; 8(58):98455-98470.
Johnsson P, Lipovich L, Dan G, Morris KV. Evolutionary conservation of long non-coding RNAs; sequence, structure, function. Biochimica et Biophysica Acta (BBA)-General Subjects. 2014;1840:1063-1071.
Park YB, Cosgrove DJ. A revised architecture of primary cell walls based on biomechanical changes induced by substrate-specific endoglucanases. Plant Physiology. 2012;158, 1933-1943.
Jan A, Yang GX, Nakamura H, et al. Characterization of a xyloglucan endotransglucosylase gene that is up-regulated by gibberellin in rice. Plant Physiology. 2004;136(3): 3670-3681.
He H, Serraj R, Yang Q. Changes in OsXTH gene expression, ABA content, and peduncle elongation in rice subjected to drought at the reproductive stage. Acta Physiology Plant. 2009;31(4): 749-756.
Tieman DM, Harriman RW, Ramamohan G, Handa AK. An antisense pectin methylesterase fene alters pectin chemistry and soluble solids in tomato fruit. The Plant Cell. 1992;4(6): 667-679.
Quesada MA, Blanco-Portales R, Posé S, et al. Antisense down-regulation of the Fa PG1 gene reveals an unexpected central role for polygalacturonase in strawberry fruit softening. Plant Physiology. 2009;(150): 1022-1032.
McQueen-Mason S, Cosgrove DJ. Disruption of hydrogen bonding between plant cell wall polymers by proteins that induce wall extension. Proceedings of The National Academy of Sciences of The United States of America. 1994;91(14): 6574-6578.
Mo B, Bewley JD. The relationship between beta-mannosidase and endo-beta-mannanase activities in tomato seeds during and following germination: a comparison of seed populations and individual seeds. Journal of Experimental Botany. 2003; 54(392):2503-2510.
Stålbrand, H, Saloheimo A , Vehmaanperä, J, et al. Cloning and expression in Saccharomyces cerevisiae of a Trichoderma reesei beta-mannanase gene containing a cellulose binding domain. Applied & Environmental Microbiology. 1995;61(3):1090-1097.
Mohammad I, Sumit G, Meli VS, et al. Fruit Ripening Regulation of α-Mannosidase Expression by the MADS Box Transcription Factor RIPENING INHIBITOR and Ethylene. Frontiers in Plant Science. 2016;7.
Andrews J, Adams SR, Burton KS, Edmondson RN. Partial purification of tomato fruit peroxidase and its effect on mechanical properties of tomato fruit skin. Journal of Experimental Botany. 2002;53:2393-2399.
Cordoba-Pedregosa C, Cordoba F, Villalba JM, Gonzaez-Reyes JA. Zonal changes in ascorbate and hydrogen peroxide contents, peroxidase, and ascorbate-related enzyme activities in onion roots. Plant Physiol. 2003;131:.697-706.
Lin CC, Kao CH. Cell wall peroxidase against ferulic acid lignin, and NaCl-reddued root growth of rice seeding. Journal of plant physiology. 2001;158:667-671.
Passardi F, Cosio C, Penel C. Dunand C. Peroxidases have more functions than a Swiss army knife. Plant Cell Reports. 2005;24:255-265.
Jablonkai I，Hatzios KK. In vitro conjugation of chloroacetanilide herbicides and atrazine with thiols and contribution of nonenzymatic conjugation to their glutathione mediated metabolism in corn. Journal of Agricultural and Food Chemistry. 1993;41( 10) : 1736-1742．
Kampranis SC，Damianova R，Atallah M , et al. A novel plant glutathione S-transferase-peroxidase suppresses bax lethality in yeast. Journal of biological chemistry. 2000;275( 38) : 29207-29216．
Marrs KA，Alfenito M，Lloyd AM, Walbot V. A glutathioneS-transferase involved in vacuolar transfer encoded bythe maize gene Bronze-2. Nature. 1995;375( 6530) :397-400．
Bartling D, Radzio R, Steiner U, Weiler EW. A glutathione S-transferase with glutathione-peroxidase activity from Ara-bidopsis thaliana. European journal of biochemistry. 1993;216( 2) : 579-586．
Trainotti L, Zanin D, Casadoro G. A cell wall-oriented genomic approach reveals a new and unexpected complexity of the softening in peaches. Journal of experimental botany. 2003;(54): 1821-1832.
Ruther J, Kleier S. Plant-plant signaling: ethylene synergizes volatile emission in Zea mays induced by exposure to (Z)-3-hexen-1-ol. Journal of chemical ecology. 2005;31:2217-2222.
Hossain MA, Kim S, Kim KH, Lee SJ, Lee H. Flavonoid compounds are enriched in lemon balm (Melissa officinalis) leaves by a high level of sucrose and confer increased antioxidant activity. Hortscience. 2009;44:1907- 1913.
Jeong SW, Das PK, Jeoung SC, et al. Ethylene suppression of sugar-induced anthocyanin pigmentation in Arabidopsis. Plant Physiology. 2010;154:1514-1531.
Kalaitzis P, Koehler SM, Tucker ML. Cloning of a tomato polygalacturonase expressed in abscission. Plant Molecular Biology. 1995;(28):647-656.
Rose JK, Cosgrove DJ, Albersheim P, Darvill AG, Bennett AB. Detection of expansin proteins and activity duringtomato fruit ontogeny. Plant Physiology. 2000;123(4):1583-1592.
Yoshida S, Tamaoki M, Ioki M, et al. Ethylene and salicylic acid control glutathione biosynthesis in ozone-exposed Arabidopsis thaliana. Physiologia Plantarum. 2009;136(3):284-298.
Gao J, Zhang YX, Li ZG, Liu MC. Role of ethylene response factors (ERFs) in fruit ripening. Food Quality and Safety. 2020;275.
Steffens B , Kovalev A , Gorb SN , Sauter M. Emerging Roots Alter Epidermal Cell Fate through Mechanical and Reactive Oxygen Species Signaling. The Plant Cell. 2012;24(8):3296-3306.
Li JS, Dai XH, Zhao YD. A role for auxin response factor 19 in auxin and ethylene signaling in Arabidopsis. Plant physiology. 2006;140(3):899-908.
Liszkay, A. Production of Reactive Oxygen Intermediates (O2{middle dot-, H₂O₂, and {middle dotOH) by Maize Roots and Their Role in Wall Loosening and Elongation Growth. Plant Physiology. 2004;136(2):3114-3123.
Wang BM, Ding GX, Wang XY et al. Changes of histological structure and water potential of huping jujube fruit cracking. Scientia Agricultura Sinica, 2013;46(21): 4558-4568.
Ivanchenko MG, den Os D, Monshausen GB, Dubrovsky JG, Bednárová A, Krishnan N. Auxin increases the hydrogen peroxide (H₂O₂) concentration in tomato (Solanum lycopersicum) root tips while inhibiting root growth. Annals of Botany. 2013;112(6):1107-1116.
Rayle DL, Cleland R. Enhancement of wall loosening and elongation by Acid solutions. Plant physiology. 1970;46(2):250-3.
Nuruddin MM，Madramootoo CA，Dodds GT. Effects of water stress at different growth stages on greenhouse tomato yield and quality．Hort Science，2003，38( 7) : 1389-1393．
Walter JM. Hereditary resistance to disease in tomato. Annual Review Phytopathology. 1967;5:131-162.
Cotner SD, Burns EE, Leeper PW. Pericarp anatomy of crackresistant and susceptible tomato fruits. Journal of the American Society for Horticultural Science. 1969;94:136-137.
Capel C, Yuste-Lisbona FJ, López-Casado G, et al. Multi-environment QTL mapping reveals genetic architecture of fruit cracking in a tomato RIL Solanum lycopersicum×S. pimpinellifolium population. Tag.theoretical & Applied Genetics. theoretische Und Angewandte Genetik. 2016;130(1):1-10.
Garcia-Luis A, Duarte A MM, Porras I, et al. Fruit splitting in “Nova” hybrid mandarin in relation to the anatomy of the fruit and fruit set treatments. Scientia Horticulturae. 1994;57: 215-231.

83.Wang SY, Camp MJ. Temperatures after bloom affect plant growth and fruit quality of strawberry. Scientia Horticulturae. 2000;85:183-199

Young HW. Inheritance of fruit cracking. Tomato Genet Coop. 1958;8:38.
Lieberman M, Kunishi A, Mapson LW, Wardale DA. Stimulation of ethylene production in apple tissue slices by methionine. Plant Physiol. 1966;41: 76–82.

86.Yang SF, Hoffman NE. Ethylene biosynthesis and its regulation in higher plants. Annu Rev Plant Physiol. 1984;35: 155-189.

Wang KL-C, Li H, Ecker JR. Ethylene Biosynthesis and Signaling Networks. The Plant Cell. 2002;14:S131-S151.
Liu JH, Peng T, Dai WS. Critical cis-Acting Elements and Interacting Transcription Factors: Key Players Associated with Abiotic Stress Responses in Plants. Plant Molecular Biology Reporter. 2014;32, 303-317.
Liao NQ, Hu ZY, Li YY, et al. Ethylene-responsive factor 4 is associated with the desirable rind hardness trait conferring cracking resistance in fresh fruits of watermelon. Plant Biotechnology Journal. 2019; 1-12.
Kim D, Langmead B, Salzberg SL. HISAT: a fast spliced aligner with low memory requirements. Nature Methods. 2015;12(4):357-360.
Trapnell C, Williams BA, Pertea G, et al. Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nature biotechnology. 2010;28:511-515.
Yu-Jian K , De-Chang Y , Lei K , et al. CPC2: a fast and accurate coding potential calculator based on sequence intrinsic features. Nucleic Acids Research.(W1):W1.
Finn RD, Tate J, Mistry J, et al. The Pfam protein families database. Nucleic Acids Research. 2008;32(1):D138.
Young MD, Wakefield MJ, Smyth GK. goseq: Gene Ontology testing for RNA-seq datasets. 2012.
Mao XZ, Cai T, Olyarchuk JG, Wei LP. Automated genome annotation and pathway identification using the KEGG Orthology(KO) as a controlled vocabulary. Bioinformatics. 2005;21: 3787-3793.
Guli S, Esteller M. Cis-acting noncoding RNAs: friends and foes. Nature structural ＆ molecular biology. 2012;19:1068-1075.
Zhan SY, Dong Y, Zhao W, et al. Genome-wide identification and characterization of long non-coding RNAs in developmental skeletal muscle of fetal goat. BMC Genomics. 2016, 17(1):666.
Yang B, Xia ZA , Zhong B, et al. Distinct Hippocampal Expression Profiles of Long Non-coding RNAs in an Alzheimer’s Disease Model. Molecular Neurobiology. 2016.
Lalitha S. Primer Premier 5. Biotech Software & Internet Report. 2000, 1(6):270-272.
Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(−Delta Delta C(T)) method. Methods. 2001;25(4):402–8.

Additional file 1: Table S1. LncRNAs in 12 tomato fruit libraries

Additional file 2: Table S2. LncRNA-mRNA pairs were identified by cis-regulation

Additional file 3: Table S3. Number of differentially expressed mRNAs and lncRNAs in different libraries

Additional file 4: Table S4. GO analysis of DEGs between groups

Additional file 5: Table S5. KEGG analysis of DEGs between groups

Additional file 6: Table S6. Key genes related to tomato fruit cracking

Additional file 7: Figure S1. Real-time PCR validation of high-throughput sequencing data. The x-axis represents the different time points of sampling, the left y-axis represents relative expression levels, and the right y-axis represents FPKM values. Blue bars represent data yielded by qRT-PCR, and red points represent data obtained by RNA sequencing. Different letters indicate significant differences (P < 0.05). (a) LA2683, (b) LA1698

Additional file 8: Figure S2. The comparison of the relative expression measured by qRT-PCR and RNA-seq

Additional file 9: Table S7. Data for the lncRNA-mRNA network

Additional file 10: Table S8. Statistics for tomato cracking rate of ‘LA2683’ and ‘LA1698’.

Additional file 11: Table S9. Detailed primer sequences for qRT-PCR.

Table 1 The overall assessment of the sequencing data.

Sample_name	Raw_reads	Clean_reads	Raw_bases(G)	Clean_bases(G)	Error rate(%)	Q20(%)	Q30(%)	GC_content(%)
CR_0h_1	102656082	102656082	15.4	15.4	0.01	98.13	95.03	42.82
CR_30h_1	80567758	80567758	12.09	12.09	0.02	96.32	91.05	42.29
CR_8h_1	79478876	79478876	11.92	11.92	0.02	96.41	91.23	42.94
CR_0h_2	97730610	97730610	14.66	14.66	0.01	97.99	94.69	42.52
CR_30h_2	100192938	100192938	15.03	15.03	0.01	97.98	94.69	42.44
CR_8h_2	94324174	94324174	14.15	14.15	0.01	98.1	94.97	42.77
CS_0h_1	95783116	95783116	14.37	14.37	0.01	97.58	93.81	42.95
CS_30h_1	90916088	90916088	13.64	13.64	0.02	96.86	92.29	42.84
CS_8h_1	86534460	86534460	12.98	12.98	0.01	97.79	94.28	42.96
CS_0h_2	103713024	103713024	15.56	15.56	0.01	98.02	94.78	42.79
CS_30h_2	98559578	98559578	14.78	14.78	0.01	97.97	94.71	42.39
CS_8h_2	114298442	114298442	17.14	17.14	0.01	97.89	94.49	43

Table 2 Key genes related to tomato fruit cracking.

Gene ID

Gene Function

Previous Studies

Solyc07g055990.3

Solyc12g011030.2

Solyc01g081060.2

xyloglucan endotransglucosylase-hydrolase 7

xyloglucan endotransglucosylase-hydrolase 9

xyloglucan endotransglucosylase-hydrolase 14 (XTH 7, XTH 9, XTH 14)

Park Y B[47]

Jan A [48]

He H [49]

Solyc10g080210.2

Solyc09g075330.3

Polygalacturonase-2 (PG 2)

Pectinesterase（PE）

Tieman DM [50]

Quesada MA [51]

Solyc08g077910.3

Expanded protein (EXP)

McQueen-Mason S [52]

Solyc01g008710.3

Mannan endo-1,4-beta-mannosidase (MAN)

Mo B [53]

Stålbrand H [54]

Mohammad I [55]

Solyc02g080530.3

Peroxide (POD)

Andrews J [56]

Cordoba-Pedregosa C [57]

Lin CC [58]

Passardi F [59]

Solyc01g081250.3

Solyc08g081010.3

Glutathione-S-transferase (GST)

Gamma-glutamylcysteine synthetase (gamma-GCS)

Jablonkai I [60]

Kampranis SC [61]

Marrs KA [62]

Bartling D [63]

Solyc09g008720.2

Solyc11g042560.2

Ethylene receptor（ER）

Ethylene-responsive transcription factor（ERF4）

Trainotti L [64]

Ruther J [65]

Hossain MA [66]

Jeong SW [67]

Kalaitzis P [68]

Rose JK [69]

Yoshida S [70]

Download PDF

Journal Publication

published 15 Apr, 2020

Read the published version in BMC Plant Biology →

Editorial decision: Minor revision
26 Mar, 2020

You are reading this older preprint version

Read the latest preprint version →

LncRNA regulates tomato fruit cracking by coordinating gene expression via a hormone-redox-cell wall network

Status:

Journal Publication

Version 5

Abstract

Figures

Background

Results

Discussion

Conclusions

Methods

Abbreviations

Declarations

References

Additional File Information

Tables

Supplementary Files

Status:

Journal Publication

Version 5