Whole transcriptome sequencing reveals genes involved in plastid/chloroplast division and development are regulated by the HP1/DDB1 at an early stage of tomato fruit development

Recent progress on mechanisms that allocate cellular space to plastids

Mechanisms that allocate cellular space to organelles are of fundamental importance to biology but remain poorly understood. A detailed understanding of mechanisms that allocate cellular space to plastids, such as chloroplasts, will lead to high-yielding crops with enhanced nutritional value. The HIGH PIGMENT (HP) genes in tomato contribute to regulated proteolysis and abscisic acid metabolism. The HP1 gene was the first gene reported to influence the amount of cellular space occupied by chloroplasts and chromoplasts almost 20 years ago. Recently, our knowledge of mechanisms that allocate cellular space to plastids was enhanced by new information on the influence of cell type on the amount of cellular space occupied by plastids and the identification of new genes that help to allocate cellular space to plastids. These genes encode proteins with unknown and diverse biochemical functions. Several transcription factors were recently reported to regulate the numbers and sizes of chloroplasts in fleshy fruit. If these transcription factors do not induce compensating effects on cell size, they should affect the amount of cellular space occupied by plastids. Although we can now propose more detailed models for the network that allocates cellular space to plastids, many gaps remain in our knowledge of this network and the genes targeted by this network. Nonetheless, these recent breakthroughs provide optimism for future progress in this field.

Galactinol synthase 2 influences the metabolism of chlorophyll, carotenoids, and ethylene in tomato fruits

Galactinol synthase (GolS), which catalyses the synthesis of galactinol, is the first critical enzyme in the biosynthesis of raffinose family oligosaccharides (RFOs) and contributes to plant growth and development, and resistance mechanisms. However, its role in fruit development remains largely unknown. In this study, we used CRISPR/Cas9 gene-editing technology in tomato (Solanum lycopersicum) to create the gols2 mutant showing uniformly green fruits without dark-green shoulders, and promoting fruit ripening. Analysis indicated that galactinol was undetectable in the ovaries and fruits of the mutant, and the accumulation of chlorophyll and chloroplast development was suppressed in the fruits. RNA-sequencing analysis showed that genes related to chlorophyll accumulation and chloroplast development were down-regulated, including PROTOCHLOROPHYLLIDE OXIDOREDUCTASE, GOLDEN 2-LIKE 2, and CHLOROPHYLL A/B-BINDING PROTEINS. In addition, early color transformation and ethylene release was prompted in the gols2 lines by regulation of the expression of genes involved in carotenoid and ethylene metabolism (e.g. PHYTOENE SYNTHASE 1, CAROTENE CIS-TRANS ISOMERASE, and 1-AMINOCYCLOPROPANE-1-CARBOXYLIC ACID SYNTHASE2/4) and fruit ripening (e.g. RIPENING INHIBITOR, NON-RIPENING, and APETALA2a). Our results provide evidence for the involvement of GolS2 in pigment and ethylene metabolism of tomato fruits.

SlCK2α as a novel substrate for CRL4 E3 ligase regulates fruit size through maintenance of cell division homeostasis in tomato

This study unravels a novel regulatory module (CRL4-CK2α-CDK2) involving fruit size control by mediating cell division homeostasis (SlCK2α and SlCDK2) in tomato. Fruit size is one of the crucial agronomical traits for crop production. UV-damaged DNA binding protein 1 (DDB1), a core component of Cullin4-RING E3 ubiquitin ligase complex (CRL4), has been identified as a negative regulator of fruit size in tomato (Solanum lycopersicum). However, the underlying molecular mechanism remains largely unclear. Here, we report the identification and characterization of a SlDDB1-interacting protein putatively involving fruit size control through regulating cell proliferation in tomato. It is a tomato homolog SlCK2α, the catalytic subunit of the casein kinase 2 (CK2), identified by yeast two-hybrid (Y2H) assays. The interaction between SlDDB1 and SlCK2α was demonstrated by bimolecular fluorescence complementation (BiFC) and co-immunoprecipitation (Co-IP). RNA interference (RNAi) and CRISPR/Cas9-based mutant analyses showed that lack of SlCK2α resulted in reduction of fruit size with reduced cell number, suggesting it is a positive regulator on fruit size by promoting cell proliferation. We also showed SlDDB1 is required to ubiquitinate SlCK2α and negatively regulate its stability through 26S proteasome-mediated degradation. Furthermore, we found that a tomato homolog of cell division protein kinase 2 (SlCDK2) could interact with and specifically be phosphorylated by SlCK2α, resulting in an increase of SlCDK2 protein stability. CRISPR/Cas9-based genetic evidence showed that SlCDK2 is also a positive regulator of fruit size by influencing cell division in tomato. Taken together, our findings, thus, unravel a novel regulatory module CRL4-CK2α-CDK2 in finely modulating cell division homeostasis and the consequences on fruit size.

Progress on Understanding Transcriptional Regulation of Chloroplast Development in Fleshy Fruit

Edible fleshy fruits are important food sources in the human diet. Their yield and nutritional quality have long been considered as breeding targets for improvement. Various developing fleshy fruits with functional chloroplasts are capable of photosynthesis and contribute to fruit photosynthate, leading to the accumulation of metabolites associated with nutritional quality in ripe fruit. Although tomato high-pigment mutants with dark-green fruits have been isolated for more than 100 years, our understanding of the mechanism of chloroplast development in fleshy fruit remain poor. During the past few years, several transcription factors that regulate chloroplast development in fleshy fruit were identified through map-based cloning. In addition, substantial progress has been made in elucidating the mechanisms that how these transcription factors regulate chloroplast development. This review provides a summary and update on this progress, with a framework for further investigations of the multifaceted and hierarchical regulation of chloroplast development in fleshy fruit.

Identification and genetic analysis of alternative splicing of long non-coding RNAs in tomato initial flowering stage

Alternative splicing (AS) is a key modulator of development in many eukaryotic organisms. Long non-coding RNAs (lncRNAs) are a class of non-coding RNAs that play essential regulatory roles in various developmental processes and stress responses. However, the functions of AS lncRNAs during the initial flowering of tomato are largely unknown. This study was designed to investigate the AS pattern of lncRNAs in tomato flower, leaf, and root tissues at the initial flowering stage. Using RNA-Seq, we found that 72.55% of lncRNAs underwent AS in these tissues, yielding a total of 16,995 AS events. Among them, the main type of AS event is alternative first exon (AFE), followed by retained intron (RI). We performed candidate target genes analysis on tissue-specific AS lncRNA, and the results indicated that the candidate target genes of these lncRNAs may be involved in the regulation of circadian rhythm, plant immunity, cellulose synthesis and phosphate-containing compound metabolic process. Moreover, a total of 73,085 putative SNPs and 15,679 InDels were detected, and the potential relationship between the AS of lncRNAs and interesting SNP and InDel loci, as well as their numbers, revealed their effects on tomato genetic diversity and genomic stability. Our data provide new insights into the complexity of the transcriptome and the regulation of AS.

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

Transcriptome sequencing and analysis during seed growth and development in Euryale ferox Salisb

BACKGROUND

Euryale ferox Salisb., an annual aquatic plant, is the only species in the genus Euryale in the Nymphaeaceae. Seeds of E. ferox are a nutritious food and also used in traditional Chinese medicine (Qian Shi in Mandarin). The molecular events that occurred during seed development in E. ferox have not yet been characterized. In this study, we performed transcriptomic analysis of four developmental stages (T1, T2, T3, and T4) in E. ferox seeds with three biological replicates per developmental stage to understand the physiological and biochemical processes during E. ferox seeds development.

RESULTS

313,844,425 clean reads were assembled into 160,107 transcripts and 85,006 unigenes with N50 lengths of 2052 bp and 1399 bp, respectively. The unigenes were annotated using five public databases (NR, COG, Swiss-Prot, KEGG, and GO). In the KEGG database, all of the unigenes were assigned to 127 pathways, of which phenylpropanoid biosynthesis was associated with the synthesis of secondary metabolites during E. ferox seed growth and development. Phenylalanine ammonia-lyase (PAL) as the first key enzyme catalyzed the conversion of phenylalanine to trans-cinnamic acid, then was related to the synthesis of flavonoids, lignins and alkaloid. The expression of PAL1 reached its peak at T3 stage, followed by a slight decrease at T4 stage. Cytochrome P450 (P450), encoded by CYP84A1 (which also called ferulate-5-hydroxylase (F5H) in Arabidopsis), was mainly involved in the biosynthesis of lignins.

CONCLUSIONS

Our study provides a transcriptomic analysis to better understand the morphological changes and the accumulation of medicinal components during E. ferox seed development. The increasing expression of PAL and P450 encoded genes in phenylpropanoid biosynthesis may promote the maturation of E. ferox seed including size, color, hardness and accumulation of medicinal components.

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.

Comprehensive Transcriptome Profiling Reveals Long Noncoding RNA Expression and Alternative Splicing Regulation during Fruit Development and Ripening in Kiwifruit (Actinidia chinensis)

Genomic and transcriptomic data on kiwifruit (Actinidia chinensis) in public databases are very limited despite its nutritional and economic value. Previously, we have constructed and sequenced nine fruit RNA-Seq libraries of A. chinensis "Hongyang" at immature, mature, and postharvest ripening stages of fruit development, and generated over 66.2 million paired-end and 24.4 million single-end reads. From this dataset, here we have identified 7051 long noncoding RNAs (lncRNAs), 29,327 alternative splicing (AS) events and 2980 novel protein-coding genes that were not annotated in the draft genome of "Hongyang." AS events were demonstrated in genes involved in the synthesis of nutritional metabolites in fruit, such as ascorbic acids, carotenoids, anthocyanins, and chlorophylls, and also in genes in the ethylene signaling pathway, which plays an indispensable role in fruit ripening. Additionally, transcriptome profiles and the contents of sugars, organic and main amino acids were compared between immature, mature, and postharvest ripening stages in kiwifruits. A total of 5931 differentially expressed genes were identified, including those associated with the metabolism of sugar, organic acid, and main amino acids. The data generated in this study provide a foundation for further studies of fruit development and ripening in kiwifruit, and identify candidate genes and regulatory elements that could serve as targets for improving important agronomic traits through marker assisted breeding and biotechnology.

Ubiquitin-conjugated degradation of golden 2-like transcription factor is mediated by CUL4-DDB1-based E3 ligase complex in tomato

CULLIN4-RING ubiquitin ligases (CRL4s) as well as their targets are fundamental regulators functioning in many key developmental and stress responses in eukaryotes. In tomato (Solanum lycopersicum), molecular cloning has revealed that the underlying genes of natural spontaneous mutations high pigment 1 (hp1), high pigment 2 (hp2) and uniform ripening (u) encode UV-DAMAGED DNA BINDING PROTEIN 1 (DDB1), DE-ETIOLATED 1 (DET1) and GOLDEN 2-LIKE (GLK2), respectively. However, the molecular basis of the opposite actions of tomato GLK2 vs CUL4-DDB1-DET1 complex on regulating plastid level and fruit quality remains unknown. Here, we provide molecular evidence showing that the tomato GLK2 protein is a substrate of the CUL4-DDB1-DET1 ubiquitin ligase complex for the proteasome degradation. SlGLK2 is degraded by the ubiquitin-proteasome system, which is mainly determined by two lysine residues (K11 and K253). SlGLK2 associates with the CUL4-DDB1-DET1 E3 complex in plant cells. Genetically impairing CUL4, DDB1 or DET1 results in a retardation of SlGLK2 degradation by the 26S proteasome. These findings are relevant to the potential of nutrient accumulation in tomato fruit by mediating the plastid level and contribute to a deeper understanding of an important regulatory loop, linking protein turnover to gene regulation.

Genome-wide identification, sequence characterization, and protein-protein interaction properties of DDB1 (damaged DNA binding protein-1)-binding WD40-repeat family members in Solanum lycopersicum

MAIN CONCLUSIONS

One hundred DDB1 (damaged DNA binding protein-1)-binding WD40-repeat domain (DWD) family genes were identified in the S. lycopersicum genome. The DWD genes encode proteins presumably functioning as the substrate recognition subunits of the cullin4-ring ubiquitin E3 ligase complex. These findings provide candidate genes and a research platform for further gene functionality and molecular breeding study. A subclass of DDB1 (damaged DNA binding protein-1)-binding WD40-repeat domain (DWD) family proteins has been demonstrated to function as the substrate recognition subunits of the cullin4-ring ubiquitin E3 ligase complex. However, little information is available about the cognate subfamily genes in tomato (S. lycopersicum). In this study, based on the recently released tomato genome sequences, 100 tomato genes encoding DWD proteins that potentially interact with DDB1 were identified and characterized, including analyses of the detailed annotations, chromosome locations and compositions of conserved amino acid domains. In addition, a phylogenetic tree, which comprises of three main groups, of the subfamily genes was constructed. The physical interaction between tomato DDB1 and 14 representative DWD proteins was determined by yeast two-hybrid and co-immunoprecipitation assays. The subcellular localization of these 14 representative DWD proteins was determined. Six of them were localized in both nucleus and cytoplasm, seven proteins exclusively in cytoplasm, and one protein either in nucleus and cytoplasm, or exclusively in cytoplasm. Comparative genomic analysis demonstrated that the expansion of these subfamily members in tomato predominantly resulted from two whole-genome triplication events in the evolution history.

Comparative Transcriptome Analysis of Resistant and Susceptible Tomato Lines in Response to Infection by Xanthomonas perforans Race T3

Bacterial spot, incited by several Xanthomonas sp., is a serious disease in tomato (Solanum lycopersicum L.). Although genetics of resistance has been widely investigated, the interactions between the pathogen and tomato plants remain unclear. In this study, tanscriptomes of X. perforans race T3 infected tomato lines were compared to those of controls. An average of 7 million reads were generated with approximately 21,526 genes mapped in each sample post-inoculation at 6 h (6 HPI) and 6 days (6 DPI) using RNA-sequencing technology. Overall, the numbers of differentially expressed genes (DEGs) were higher in the resistant tomato line PI 114490 than in the susceptible line OH 88119, and the numbers of DEGs were higher at 6 DPI than at 6 HPI. Fewer genes (78 in PI 114490 and 15 in OH 88119) were up-regulated and most DEGs were down-regulated, suggesting that the inducible defense response might not be fully activated at 6 HPI. Accumulation expression levels of 326 co-up regulated genes in both tomato lines at 6 DPI might be involved in basal defense, while the specific and strongly induced genes at 6 DPI might be correlated with the resistance in PI 114490. Most DEGs were involved in plant hormone signal transduction, plant-pathogen interaction and phenylalanine metabolism, and the genes significantly up-regulated in PI 114490 at 6 DPI were associated with defense response pathways. DEGs containing NBS-LRR domain or defense-related WRKY transcription factors were also identified. The results will provide a valuable resource for understanding the interactions between X. perforans and tomato plants.

Transcriptome sequencing of a chimaera reveals coordinated expression of anthocyanin biosynthetic genes mediating yellow formation in herbaceous peony (Paeonia lactiflora Pall.)

BACKGROUND

Herbaceous peony (Paeonia lactiflora Pall.) is a traditional flower in China and a wedding attractive flower in worldwide. In its flower colour, yellow is the rarest which is ten times the price of the other colours. However, the breeding of new yellow P. lactiflora varieties using genetic engineering is severely limited due to the little-known biochemical and molecular mechanisms underlying its characteristic formation.

RESULTS

In this study, two cDNA libraries generated from P. lactiflora chimaera with red outer-petal and yellow inner-petal were sequenced using an Illumina HiSeq™ 2000 platform. 66,179,398 and 65,481,444 total raw reads from red outer-petal and yellow inner-petal cDNA libraries were generated, which were assembled into 61,431 and 70,359 Unigenes with an average length of 628 and 617 nt, respectively. Moreover, 61,408 non-redundant All-unigenes were obtained, with 37,511 All-unigenes (61.08%) annotated in public databases. In addition, 6,345 All-unigenes were differentially expressed between the red outer-petal and yellow inner-petal, with 3,899 up-regulated and 2,446 down-regulated All-unigenes, and the flavonoid metabolic pathway related to colour development was identified using the Kyoto encyclopedia of genes and genomes database (KEGG). Subsequently, the expression patterns of 10 candidate differentially expressed genes (DEGs) involved in the flavonoid metabolic pathway were examined, and flavonoids were qualitatively and quantitatively analysed. Numerous anthoxanthins (flavone and flavonol) and a few anthocyanins were detected in the yellow inner-petal, which were all lower than those in the red outer-petal due to the low expression levels of the phenylalanine ammonialyase gene (PlPAL), flavonol synthase gene (PlFLS), dihydroflavonol 4-reductase gene (PlDFR), anthocyanidin synthase gene (PlANS), anthocyanidin 3-O-glucosyltransferase gene (Pl3GT) and anthocyanidin 5-O-glucosyltransferase gene (Pl5GT).

CONCLUSION

Transcriptome sequencing (RNA-Seq) analysis based on the high throughput sequencing technology was an efficient approach to identify critical genes in P. lactiflora and other non-model plants. The flavonoid metabolic pathway and glucide metabolic pathway were identified as relatived yellow formation in P. lactiflora, PlPAL, PlFLS, PlDFR, PlANS, Pl3GT and Pl5GT were selected as potential candidates involved in flavonoid metabolic pathway, which inducing inhibition of anthocyanin biosynthesis mediated yellow formation in P. lactiflora. This study could lay a theoretical foundation for breeding new yellow P. lactiflora varieties.