Transcriptomic, metabolomic, and ATAC-seq analysis reveal the regulatory mechanism of senescence of post-harvest tomato fruit.

Several physiological changes occur during fruit storage, which include the regulation of genes, metabolisms and transcription factors. In this study, we compared 'JF308' (a normal tomato cultivar) and 'YS006' (a storable tomato cultivar) to determine the difference in accumulated metabolites, gene expression, and accessible chromatin regions through metabolome, transcriptome, and ATAC-seq analysis. A total of 1006 metabolites were identified in two cultivars. During storage time, sugars, alcohols and flavonoids were found to be more abundant in 'YS006' compared to 'JF308' on day 7, 14, and 21, respectively. Differentially expressed genes, which involved in starch and sucrose biosynthesis were observed higher in 'YS006'. 'YS006' had lower expression levels of CesA (cellulose synthase), PL (pectate lyase), EXPA (expansin) and XTH (xyglucan endoglutransglucosylase/hydrolase) than 'JF308'. The results showed that phenylpropanoid pathway, carbohydrate metabolism and cell wall metabolism play important roles in prolonging the shelf life of tomato (Solanum lycopersicum) fruit. The ATAC-seq analysis revealed that the most significantly up-regulated transcription factors during storage were TCP 2,3,4,5, and 24 in 'YS006' compared to 'JF308' on day 21. This information on the molecular regulatory mechanisms and metabolic pathways of post-harvest quality changes in tomato fruit provides a theoretical foundation for slowing post-harvest decay and loss, and has theoretical importance and application value in breeding for longer shelf life cultivars.

Physiological, transcriptomic, and metabolic analyses reveal that mild salinity improves the growth, nutrition, and flavor properties of hydroponic Chinese chive (Allium tuberosum Rottler ex Spr).

Environmental stressors such as salinity have pronounced impacts on the growth, productivity, nutrition, and flavor of horticultural crops, though yield loss sometimes is inevitable. In this study, the salinity influences were evaluated using hydroponic Chinese chive (Allium tuberosum) treated with different concentrations of sodium chloride. The results demonstrated that lower salinity could stimulate plant growth and yield. Accordingly, the contents of soluble sugar, ascorbic acid, and soluble protein in leaf tissues increased, following the decrease of the nitrate content, under mild salinity (6.25 or 12.5 mM NaCl). However, a higher level of salinity (25 or 50 mM NaCl) resulted in growth inhibition, yield reduction, and leaf quality deterioration of hydroponic chive plants. Intriguingly, the chive flavor was boosted by the salinity, as evidenced by pungency analysis of salinity-treated leaf tissues. UPLC-MS/MS analysis reveals that mild salinity promoted the accumulation of glutamic acid, serine, glycine, and proline in leaf tissues, and thereby enhanced the umami and sweet flavors of Chinese chive upon salinity stress. Considering the balance between yield and flavor, mild salinity could conduce to hydroponic Chinese chive cultivation. Transcriptome analysis revealed that enhanced pungency could be ascribed to a salt stress-inducible gene, AtuFMO1, associated with the biosynthesis of S-alk(en)yl cysteine sulphoxides (CSOs). Furthermore, correlation analysis suggested that two transcription factors, AtubHLH and AtuB3, were potential regulators of AtuFMO1 expressions under salinity. Thus, these results revealed the molecular mechanism underlying mild salinity-induced CSO biosynthesis, as well as a practical possibility for producing high-quality Chinese chive hydroponically.

Identifying Quantitative Trait Loci for Thousand Grain Weight in Eggplant by Genome Re-Sequencing Analysis.

Eggplant (Solanum melongena L.; 2n = 24) is one of the most important Solanaceae vegetables and is primarily cultivated in China (approximately 42% of world production) and India (approximately 39%). Thousand-grain weight (TGW) is an important trait that affects eggplant breeding cost and variety promotion. This trait is controlled by quantitative trait loci (QTLs); however, no quantitative trait loci (QTL) has been reported for TGW in eggplant so far, and its potential genetic basis remain unclear. In this study, two eggplant lines, 17C01 (P1, wild resource, small seed) and 17C02 (P2, cultivar, large seed), were crossed to develop F1, F2 (308 lines), BC1P1 (44 lines), and BC1P2 (44 lines) populations for quantitative trait association analysis. The TGWs of P1, P2 and F1 were determined as 3.00, 3.98 and 3.77 g, respectively. The PG-ADI (polygene-controlled additive-dominance-epistasis) genetic model was identified as the optimal model for TGW and the polygene heritability value in the F2 generation was as high as 80.87%. A high-quality genetic linkage bin map was constructed with resequencing analysis. The map contained 3,918 recombination bins on 12 chromosomes, and the total length was 1,384.62 cM. A major QTL (named as TGW9.1) located on chromosome 9 was identified to be strongly associated with eggplant TGW, with a phenotypic variance explanation of 20.51%. A total of 45 annotated genes were identified in the genetic region of TGW9.1. Based on the annotation of Eggplant genome V3 and orthologous genes in Arabidopsis thaliana, one candidate gene SMEL_009g329850 (SmGTS1, encoding a putative ubiquitin ligase) contains 4 SNPs and 2 Indels consecutive intron mutations in the flank of the same exon in P1. SmGTS1 displayed significantly higher expression in P1 and was selected as a potential candidate gene controlling TGW in eggplant. The present results contribute to shed light on the genetic basis of the traits exploitable in future eggplant marker-assisted selection (MAS) breeding.

First report of tomato spotted wilt virus infecting water dropwort in China.

Water dropwort (Oenanthe javanica) is an aquatic perennial plant that has been cultivated in many regions in Asia for thousands of years. In China, it is an economically important vegetable that has been consumed as food, while also being used as a folk remedy to alleviate diseases (Liu et al., 2021). In 2021, during a disease survey of a greenhouse in Beijing, China, chlorotic spots were detected on many water dropwort plants (Fig. S1A). Twenty-seven water dropwort samples were collected for the extraction of total RNA using the TRIzol reagent (Invitrogen, USA). High-quality RNA samples from three water dropwort plants were combined and used as the template for constructing a single small RNA library (BGI-Shenzhen Company, China). The Velvet 1.0.5 software was used to assemble the clean reads (18 to 28 nt) into larger contigs, which were then compared with the nucleotide sequences in the National Center database using the BLASTn algorithm. Thirty-eight contigs matched sequences in the tomato spotted wilt virus (TSWV) genome. No other viruses were detected. Twenty-seven leaf samples were analyzed in an enzyme-linked immunosorbent assay (ELISA) with anti-TSWV antibody (Agdia, USA), which revealed 17 positive reaction. Two sets of primer pairs targeting different parts of the S RNA (Table S1) was used to verify the TSWV infection on water dropwort by reverse transcription (RT)-PCR followed by Sanger sequencing (BGI-Shenzhen, China). The TSWV target sequences were amplified from 17 samples, which was consistent with the ELISA results. The sequenced 861-bp PCR product shared 99.8% nucleotide sequence identity with TSWV isolate MR-01 (MG593199), while the 441-bp amplicon shared a 99.2% nucleotide sequence identity with MR-01 (MG593199). To obtain the whole genome sequence of TSWV (S, M, and L RNA sequences), specific RT-PCR primers were designed (Table S1) and used to amplify their respective fragments from one representative sample (TSWV-water dropwort). The amplified products were inserted into PCE2TA/Blunt-Zero vector (Vazyme Biotech Co., Ltd, China) and then sequenced (BGI-Shenzhen, China). The S, M, and L RNA sequences were determined to be 2,952 nt (accession no. OM154969), 4,776 nt (accession no. OM154970), and 8,914 nt (accession no. OM154971), respectively. BLASTn analysis demonstrated that the whole genome sequence was highly conserved. The nucleotide identities between this isolate and other TSWV isolates ranged from 98.6% to 99.6% (S RNA), 98.9% to 99.2% (M RNA), and 97.3% to 98.7% (L RNA). Using MEGA 7.0, the phylogenetic relationships of TSWV were determined on the basis of the S, M, and L RNA full-length sequences (Kumar et al., 2016). In the S RNA derived phylogenetic tree, the water dropwort isolate was closely related to the MR-01 isolate from the USA (MG593199). In the M RNA and L RNA derived phylogenetic trees, the water dropwort isolate formed a branch with only a TSWV isolate from eggplant. Additionally, the M and L RNA sequences were most similar to sequences in TSWV isolates from China and Korea, respectively (Fig. S1B). To the best of our knowledge, this is the first report of water dropwort as a natural host for TSWV in China and the second report worldwide since the first finding in the Korea (Kil et al. 2020). TSWV has caused serious problems on many crops in the world, and the infection of TSWV on water dropwort in a greenhouse should not be looked lightly. Firstly, the virus can be passed on from generation to generation in infected water dropwort due to the vegetative propagation mode of the plant in production, thus threaten the production of this vegetable crop. In addition, infected water dropwort may serve as a reservoir for the virus, thus potentially posing a threat for causing TSWV spread in the affected greenhouses. The author(s) declare no conflict of interest. Funding: This research was supported by the Beijing Academy of Agriculture and Forestry Foundation, China (QNJJ202131, KJCX20200212, and KJCX20200113). References: Kil et al. 2020. Plant Pathol. J. 36: 67-75 Kumar et al. 2016. Mol Biol Evol, 33: 1870-1874 Liu et al. 2021. Horticulture Research. 8:1-17.

Determination of reliable reference genes for gene expression studies in Chinese chive (Allium tuberosum) based on the transcriptome profiling.

Chinese chive (Allium tuberosum) is widely cultivated around the world for its unique flavor, nutrient, and medicinal values, yet its molecular mechanism on flavor formation and other metabolic pathways remains intangible. The elucidation of these complex processes begins with investigating the expression of the genes of interest, however the appropriate reference genes (RGs) for normalizing the gene expression are still unavailable in A. tuberosum. To fill this lacuna, transcriptome-wide screening was undertaken to identify the most stable genes according to the analysis of their FPKM values. The expression stability of the RGs was further evaluated using geNorm, NormFinder, BestKeeper, and RefFinder algorithms. The comprehensive analysis showed that GLY1 and SKP1, instead of two traditionally used RGs (eIF1α and ACT2), were the most stable genes across diverse A. tuberosum tissues, indicating the necessity to carefully validate the stability of RGs prior to their use for normalizations. As indicated by geNorm, the normalizations with at least two RGs could give more accurate results. qRT-PCR experiments were conducted with randomly selected genes, demonstrating that normalization with a combination of GLY1 and SKP1 resulted in reliable normalization results. Our finding represents the first attempt toward establishing a standardized qRT-PCR analysis in this economically important vegetable.

Transcriptome landscapes of multiple tissues highlight the genes involved in the flavor metabolic pathway in Chinese chive (Allium tuberosum).

The unique flavor of Allium tuberosum is primarily associated with the hydrolysis of a series of organosulfur compounds, S-alk(en)yl cysteine sulphoxides (CSOs), upon tissue bruising or maceration. To obtain the tissue-specific transcriptomes, 18 RNA-Seq libraries representing leaf, root, stem, mature flower, inflorescence, and seed tissues of A. tuberosum were sequenced, finally yielding 133.7 Gb clean reads. The de novo assembled transcriptomes enabled the identification of 223,529 unigenes, which were functionally annotated and analyzed for the gene ontology and metabolic pathways. Furthermore, to reveal the flavor metabolic pathways, a total of 205 unigenes involved in the sulfur assimilation and CSO biosynthesis were identified, and their expression profiles were analyzed by RNA-Seq and qRT-PCR. Collectively, this study provides a valuable resource for in-depth molecular and functional researches especially on flavor formation, as well as for the development of molecular markers, and other genetic studies in A. tuberosum.

Comprehensive Transcriptome Reveals an Opposite Regulatory Effect of Plant Growth Retardants in Controlling Seedling Overgrowth between Roots and Shoots.

Seedling overgrowth always develops in undernourished plants due to biotic or abiotic stresses, which significantly decrease the yield of crops and vegetables. It is known that the plant growth retardants paclobutrazol (PBZ) and chlormequat chloride (CCC) are the most commonly used chemicals in controlling seedling height in plants by regulating the gibberellin (GA) biosynthesis pathway. However, the exact molecular regulation mechanism remains largely unknown. This study performed a comprehensive transcriptome profile to identify significantly differentially expressed genes after adding CCC and PBZ to the water culture seedling raising system for the first time. According to the obviously restrained shoots and roots, the GA biosynthesis genes were significantly decreased, as well as the endogenous GA content being reduced. Intriguingly, the GA signaling pathway genes were affected in opposite ways, increasing in roots but decreasing in shoots, especially regarding the phytochrome interacting factor SlPIF1 and the downstream genes expansins (SlEXPs), which promote cell wall remodeling. Further study found that the most down-regulated genes SlEXPA5 and SlEXPA15 were expressed specifically in shoot tissue, performing the function of repressing elongation, while the up-regulated genes SlEXPB2 and SlEXPB8 were proven to be root-specific expressed genes, which may promote horizontal elongation in roots. This research reported the comprehensive transcriptome profiling of plant growth retardants in controlling seedling overgrowth and restraining GA biosynthesis through the regulation of the GA signaling-related genes SlPIF1 and SlEXPs, with an opposite expression pattern between roots and shoots.