Functional and Transcriptome Analysis Reveal Specific Roles of Dimocarpus longan DlRan3A and DlRan3B in Root Hair Development, Reproductive Growth, and Stress Tolerance.

Ran GTPases play essential roles in plant growth and development. Our previous studies revealed the nuclear localization of DlRan3A and DlRan3B proteins and proposed their functional redundancy and distinction in Dimocarpus longan somatic embryogenesis, hormone, and abiotic stress responses. To further explore the possible roles of DlRan3A and DlRan3B, gene expression analysis by qPCR showed that their transcripts were both more abundant in the early embryo and pulp in longan. Heterologous expression of DlRan3A driven by its own previously cloned promoter led to stunted growth, increased root hair density, abnormal fruits, bigger seeds, and enhanced abiotic stress tolerance. Conversely, constitutive promoter CaMV 35S (35S)-driven expression of DlRan3A, 35S, or DlRan3B promoter-controlled expression of DlRan3B did not induce the alterations in growth phenotype, while they rendered different hypersensitivities to abiotic stresses. Based on the transcriptome profiling of longan Ran overexpression in tobacco plants, we propose new mechanisms of the Ran-mediated regulation of genes associated with cell wall biosynthesis and expansion. Also, the transgenic plants expressing DlRan3A or DlRan3B genes controlled by 35S or by their own promoter all exhibited altered mRNA levels of stress-related and transcription factor genes. Moreover, DlRan3A overexpressors were more tolerant to salinity, osmotic, and heat stresses, accompanied by upregulation of oxidation-related genes, possibly involving the Ran-RBOH-CIPK network. Analysis of a subset of selected genes from the Ran transcriptome identified possible cold stress-related roles of brassinosteroid (BR)-responsive genes. The marked presence of genes related to cell wall biosynthesis and expansion, hormone, and defense responses highlighted their close regulatory association with Ran.

Integrated Metabolomic and Transcriptomic Analyses Reveal the Potential Regulation of Flavonoids in the Production of Embryogenic Cultures during Early Somatic Embryogenesis of Longan (Dimocarpus longan Lour.).

Embryogenic cultures of longan (Dimocarpus longan Lour.) contain various metabolites with pharmacological properties that may function in the regulation of somatic embryogenesis (SE). In this study, based on widely targeted metabolomics, 501 metabolites were obtained from the embryogenic calli, incomplete compact proembryogenic cultures, and globular embryos during early SE of longan, among which 41 flavonoids were differentially accumulated during the SE. Using RNA sequencing, 36 flavonoid-biosynthesis-related genes and 43 MYB and 52 bHLH transcription factors were identified as differentially expressed genes. Furthermore, Kyoto Encyclopedia of Genes and Genomes enrichment analysis revealed that the flavonoid metabolism-related pathways were significantly enriched during the early SE. These results suggested that the changes in flavonoid levels in the embryogenic cultures of longan were mediated by MYBs and bHLHs via regulating flavonoid-biosynthesis-related genes, thus potentially regulating early SE. The identified metabolites in the embryogenic cultures of longan can be used to develop pharmaceutical ingredients.

Genome-wide high-throughput chromosome conformation capture analysis reveals hierarchical chromatin interactions during early somatic embryogenesis.

Somatic embryogenesis (SE), like zygotic embryo development, is a progressive process. Early SE is the beginning of a switch from a somatic to an embryogenic state and is an important stage for initiating chromatin reprogramming of SE. Previous studies suggest that changes in chromatin accessibility occur during early SE, although information on the 3D structure of chromatin is not yet available. Here, we present a chromosome-level genome assembly of longan (Dimocarpus longan) using PacBio combined with high-through chromosome conformation capture scaffolding, which resulted in a 446 Mb genome assembly anchored onto 15 scaffolds. During early SE, chromatin was concentrated and then decondensed, and a large number of long terminal repeat retrotransposons (LTR-RTs) were enriched in the local chromatin interaction region, suggesting LTR-RTs were involved in chromatin reorganization. Early SE was accompanied by the transformation from A to B compartments, and the interactions between B compartments were enhanced. Results from chromatin accessibility, monomethylation of histone H3 at lysine 4 (H3K4me1) modification, and transcription analyses further revealed a gene regulatory network for cell wall thickening during SE. Particularly, we found that the H3K4me1 differential peak binding motif showed abnormal activation of ethylene response factor transcription factors and participation in SE. The chromosome-level genomic and multiomics analyses revealed the 3D conformation of chromatin during early SE, providing insight into the molecular mechanisms underlying cell wall thickening and the potential regulatory networks of TFs during early SE in D. longan. These results provide additional clues for revealing the molecular mechanisms of plant SE.

Single-cell RNA sequencing analysis of the embryogenic callus clarifies the spatiotemporal developmental trajectories of the early somatic embryo in Dimocarpus longan.

Plant embryogenic calli (ECs) can undergo somatic embryogenesis to regenerate plants. This process is mediated by regulatory factors, such as transcription factors and specifically expressed genes, but the precise molecular mechanisms underlying somatic embryogenesis at the single-cell level remain unclear. In this study, we performed high-resolution single-cell RNA sequencing analysis to determine the cellular changes in the EC of the woody plant species Dimocarpus longan (longan) and clarify the continuous cell differentiation trajectories at the transcriptome level. The highly heterogeneous cells in the EC were divided into 12 putative clusters (e.g., proliferating, meristematic, vascular, and epidermal cell clusters). We determined cluster-enriched expression marker genes and found that overexpression of the epidermal cell marker gene GDSL ESTERASE/LIPASE-1 inhibited the hydrolysis of triacylglycerol. In addition, the stability of autophagy was critical for the somatic embryogenesis of longan. The pseudo-timeline analysis elucidated the continuous cell differentiation trajectories from early embryonic cell division to vascular and epidermal cell differentiation during the somatic embryogenesis of longan. Moreover, key transcriptional regulators associated with cell fates were revealed. We found that ETHYLENE RESPONSIVE FACTOR 6 was characterized as a heat-sensitive factor that negatively regulates longan somatic embryogenesis under high-temperature stress conditions. The results of this study provide new spatiotemporal insights into cell division and differentiation during longan somatic embryogenesis at single-cell resolution.

The function of BoTCP25 in the regulation of leaf development of Chinese kale.

XG Chinese kale (Brassica oleracea cv. 'XiangGu') is a variety of Chinese kale and has metamorphic leaves attached to the true leaves. Metamorphic leaves are secondary leaves emerging from the veins of true leaves. However, it remains unknown how the formation of metamorphic leaves is regulated and whether it differs from normal leaves. BoTCP25 is differentially expressed in different parts of XG leaves and respond to auxin signals. To clarify the function of BoTCP25 in XG Chinese kale leaves, we overexpressed BoTCP25 in XG and Arabidopsis, and interestingly, its overexpression caused Chinese kale leaves to curl and changed the location of metamorphic leaves, whereas heterologous expression of BoTCP25 in Arabidopsis did not show metamorphic leaves, but only an increase in leaf number and leaf area. Further analysis of the expression of related genes in Chinese kale and Arabidopsis overexpressing BoTCP25 revealed that BoTCP25 could directly bind the promoter of BoNGA3, a transcription factor related to leaf development, and induce a significant expression of BoNGA3 in transgenic Chinese kale plants, whereas this induction of NGA3 did not occur in transgenic Arabidopsis. This suggests that the regulation of Chinese kale metamorphic leaves by BoTCP25 is dependent on a regulatory pathway or elements specific to XG and that this regulatory element may be repressed or absent from Arabidopsis. In addition, the expression of miR319's precursor, a negative regulator of BoTCP25, also differed in transgenic Chinese kale and Arabidopsis. miR319's transcrips were significantly up-regulated in transgenic Chinese kale mature leaves, while in transgenic Arabidopsis, the expression of miR319 in mature leaves was kept low. In conclusion, the differential expression of BoNGA3 and miR319 in the two species may be related to the exertion of BoTCP25 function, thus partially contributing to the differences in leaf phenotypes between overexpressed BoTCP25 in Arabidopsis and Chinese kale.

Riboflavin mediates m6A modification targeted by miR408, promoting early somatic embryogenesis in longan.

Plant somatic embryogenesis (SE) is an in vitro biological process wherein bipolar structures are induced to form somatic cells and regenerate into whole plants. MicroRNA (miRNA) is an essential player in plant SE. However, the mechanism of microRNA408 (miR408) in SE remains elusive. Here, we used stable transgenic technology in longan (Dimocarpus longan) embryogenic calli to verify the mechanism by which miR408 promotes cell division and differentiation of longan early SE. dlo-miR408-3p regulated riboflavin biosynthesis by targeting nudix hydrolase 23 (DlNUDT23), a previously unidentified gene mediating N6-methyladenosine (m6A) modification and influencing RNA homeostasis and cell cycle gene expression during longan early SE. We showed that DlMIR408 overexpression (DlMIR408-OE) promoted 21-nt miRNA biosynthesis. In DlMIR408-OE cell lines, dlo-miR408-3p targeted and downregulated DlNUDT23, promoted riboflavin biosynthesis, decreased flavin mononucleotide (FMN) accumulation, promoted m6A level, and influenced miRNA homeostasis. DNA replication, glycosylphosphatidylinositol (GPI)-anchor biosynthesis, the pentose phosphate pathway, and taurine and hypotaurine metabolism were also closely associated with riboflavin metabolism. In a riboflavin feeding assay, dlo-miR408-3p and pre-miR408 were upregulated and DlNUDT23 was downregulated, increasing the m6A level and cell division and differentiation in longan globular embryos. When riboflavin biosynthesis was inhibited, dlo-miR408-3p was downregulated and DlNUDT23 was upregulated, which decreased m6A modification and inhibited cell division but did not inhibit cell differentiation. FMN artificial demethylated m6A modification affected the homeostasis of precursor miRNA and miRNA. Our results revealed a mechanism underlying dlo-miR408-3p-activated riboflavin biosynthesis in which DlNUDT23 is targeted, m6A modification is dynamically mediated, and cell division is affected, promoting early SE in plants.

Genome-wide identification of miRNAs and targets associated with cell wall biosynthesis: Differential roles of dlo-miR397a and dlo-miR408-3p during early somatic embryogenesis in longan.

The dynamic alterations in cell wall (CW) biosynthesis play an essential role in physiological isolation during the plant somatic embryogenesis (SE). However, the mechanisms underlying the functions of cell wall-associated miRNAs (CW-miRNA) remain poorly understood in plant SE. Here, we have identified 36 distinct candidate miRNAs associated with CW biosynthesis from longan third-generation genome as well as miRNA transcriptome, and modified RLM-RACE validated four distinct miRNA, which specifically targeted four CW-related genes. More importantly, we found that the dlo-miR397a-antagomir significantly enhanced DlLAC7 expression and improved laccase activity. Interestingly, inhibition of dlo-miR397a increased CW lignin deposition and promoted the tightening of protodermal cell by miRNA-mimic technology during early SE. Moreover, overexpression of dlo-miR408-3p (dlo-miR408-3p-agomir) markedly decreased DlLAC12 expression. dlo-miR408-3p-agomir activated rapid cell division, thus promoting the globular embryo (GE) development, which might be due to high DNA synthesis activity in protoepidermal cells, rather than affecting lignin synthesis. The subcellular location also indicated that both DlLAC7 and DlLAC12 proteins were primarily localized in CW and regulated CW biosynthesis. Overall, our findings provided new insight on the molecular regulatory networks comprising various miRNAs associated with cell wall, and established that dlo-miR397a and dlo-miR408-3p played differential roles during early SE in longan. The findings also shed some light on the potential role of miRNA target DlLAC regulating in vivo embryonic development of plant.

Genome-wide identification, evolution analysis of cytochrome P450 monooxygenase multigene family and their expression patterns during the early somatic embryogenesis in Dimocarpus longan Lour.

Cytochrome P450 (CYP), a multi-gene superfamily, is involved in a broad range of physiological processes, including hormone responses and secondary metabolism throughout the plant life cycle. Longan (Dimocarpus longan), a subtropical and tropical evergreen fruit tree, its embryonic development is closely related to the yield and quality of fruits. And a large number of secondary metabolites, such as flavonoids and carotenoids, are also produced during the longan somatic embryogenesis (SE). It is important, therefore, to study potential functions of CYPs in longan. However, the knowledge of longan CYPs is still very limited. Here, a total of 327 DlCYPs were identified using the genome-search method, which could be classified into nine clans. The expansion of the DlCYP family was mainly caused by tandem duplication (TD) events. Promoter cis-acting elements analysis elucidated that DlCYPs played important roles in hormonal responses. A total of 246 DlCYPs exhibited six different expression patterns during the early SE based on longan transcriptomic data. Eight DlCYPs underwent alternative splicing (AS) events, and they might produce one to six isoforms. And the AS transcript of DlCYP97C1 might act as an alternative to the full-length transcript in ICpEC and GE stages. Finally, protein-protein interaction (PPI) networks and miRNA target prediction elucidated that DlCYPs might be involved in the phenylpropanoid metabolic pathway and primarily regulated and targeted by miR413. In summary, our results provided valuable inventory for understanding the classification and biological functions of DlCYPs and provided insight into further functional verification of DlCYPs during the longan early SE.

Genome-wide circular RNA profiling and competing endogenous RNA regulatory network analysis provide new insights into the molecular mechanisms underlying early somatic embryogenesis in Dimocarpus longan Lour.

Circular RNAs (circRNAs) are widely involved in plant growth and development. However, the function of circRNAs in plant somatic embryogenesis (SE) remains elusive. Here, by using high-throughput sequencing, a total of 5029 circRNAs were identified in the three stages of longan (Dimocarpus longan Lour.) early SE. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis revealed that differentially expressed (DE) circRNA host genes were enriched in the 'non-homologous end-joining' (NHEJ) and 'butanoate metabolism' pathways. In addition, the reactive oxygen species (ROS) content during longan early SE was determined. The results indicated that ROS-induced DNA double-strand breaks may not depend on the NHEJ repair pathway. Correlation analyses of the levels of related metabolites (glutamate, γ-aminobutyrate and pyruvate) and the expression levels of circRNAs and their host genes involved in butanoate metabolism were performed. The results suggested that circRNAs may act as regulators of the expression of cognate mRNAs, thereby affecting the accumulation of related compounds. A competing endogenous RNA (ceRNA) network of DE circRNAs, DE mRNAs, DE long noncoding RNAs (lncRNAs) and DE microRNAs (miRNAs) was constructed. The results showed that the putative targets of the noncoding RNA (ncRNAs) were significantly enriched in the KEGG pathways 'mitogen-activated protein kinase signaling' and 'nitrogen metabolism'. Furthermore, the expression patterns of the candidate circRNAs, lncRNAs, miRNAs and mRNAs confirmed the negative correlation between miRNAs and ceRNAs. In addition, two circRNA overexpression vectors were constructed to further verify the ceRNA network correlations in longan early SE. Our study revealed the potential role of circRNAs in longan early SE, providing new insights into the intricate regulatory mechanism underlying plant SE.

Genome-wide identification and comprehensive analyses of NAC transcription factor gene family and expression patterns during somatic embryogenesis in Dimocarpus longan Lour.

The NAM, ATAF1/2, and CUC2 form a huge plant-specific gene family of NAC TFs that are involved in the growth, development, and regulation of biotic and abiotic stress responses. Although the draft genome of longan (Dimocarpus longan Lour.) has been published, however the comprehensive data regarding the functions, evolution, and expression patterns of the NAC family are still unavailable. In this study, a comprehensive analysis of the NAC transcription factor family in longan was performed, and a total of 114 NAC genes were found. We investigated the NAC gene family exploring the phylogeny, domain conservation, intron/exon, motifs, cis-regulatory elements, protein-protein interaction, and expression profiles of RNA-seq samples in different tissues and early somatic embryogenesis of longan. Phylogenetic analysis showed that the genes with similar gene structure and motif distribution were clustered in the same group. Cis-element identification indicates the possible role of NAC genes in biological and physiological processes. Protein-protein interaction identified the DlNACs homologous with Arabidopsis proteins. We further investigated the expression pattern of DlNAC genes in different tissues (pulp, stem, large fruit, young fruit, and flower) during somatic embryogenesis at embryogenic callus (EC), incomplete compact pro-embryogenic cultures (ICpEC), and globular embryos (GE) stages. The qRT-PCR results showed that the DlNAC genes were expressed higher at EC and GE stage compared with ICpEC stage. In conclusion, our results provide insight into the evolution, diversity, and characterization of NAC genes in the longan and provide a base for understanding their biological roles and molecular mechanisms in plants.

Genome-wide identification and expression analysis of the SR gene family in longan (Dimocarpus longan Lour.).

Longan (Dimocarpus longan Lour.) is an important commercial fruit tree in southern China. The embryogenesis of longan affects the quality and yield of fruit. A large number of alternative splicing events occurs during somatic embryogenesis (SE), which is regulated by serine/arginine-rich (SR) proteins. However, the functions of SR proteins in longan are poorly understood. In this study, 21 Dlo-SR gene family members belonging to six subfamilies were identified, among which Dlo-RSZ20a, Dlo-SR30, Dlo-SR17, Dlo-SR53 and Dlo-SR32 were localized in the nucleus, Dlo-RSZ20b, Dlo-RSZ20c, Dlo-RSZ20d, Dlo-SC18, Dlo-RS2Z29, Dlo-SCL41, and Dlo-SR33 were localized in chloroplasts, and Dlo-RS43, Dlo-SC33, Dlo-SC37, Dlo-RS2Z33, Dlo-RS2Z16, Dlo-RS2Z24, Dlo-SCL43, Dlo-SR112, and Dlo-SR59 were localized in the nucleus and chloroplasts. The Dlo-SR genes exhibited differential expression patterns in different tissues of longan. The transcript levels of Dlo-RSZ20a, Dlo-SC18, Dlo-RS2Z29, DLo-SR59, Dlo-SR53, and Dlo-SR17 were low in all analyzed tissues, whereas Dlo-RS43, Dlo-RS2Z16, Dlo-RS2Z24, and Dlo-SR30 were highly expressed in all tissues. To clarify their function during SE, the transcript levels of Dlo-SR genes were analyzed at different four stages of SE, comprising non-embryonic callus (NEC), friable-embryogenic callus (EC), incomplete compact pro-embryogenic culture (ICpEC) and globular embryo (GE). Interestingly, the transcript levels of Dlo-RS2Z29 and Dlo-SR112 were increased in embryogenic cells compared with the NEC stage, whereas transcript levels of Dlo-RSZ20a, Dlo-RS43, Dlo-SC37, and Dlo-RS2Z16 were especially increased at the GE stage compared with the other stages. Alternative splicing events of Dlo-SR mRNA precursors (pre-mRNAs) was detected during SE, with totals of 41, 29, 35, and 44 events detected during NEC, EC, ICpEC, and GE respectively. Protein-protein interaction analysis showed that SR proteins were capable of interaction with each other. The results indicate that the alternative splicing of Dlo-SR pre-mRNAs occurs during SE and that Dlo-SR proteins may interact to regulate embryogenesis of longan.

Genome-wide investigation of DNA methylation dynamics reveals a critical role of DNA demethylation during the early somatic embryogenesis of Dimocarpus longan Lour.

DNA methylation plays essential roles in gene regulation, chromatin structure stability, gene imprinting, X chromosome inactivation and embryonic development. However, the dynamics and functions of DNA methylation during the early stage of longan (Dimocarpus longan) somatic embryogenesis (SE) are still unclear. In this study, we carried out whole genome bisulphite sequencing and transcriptome sequencing analyses for embryogenic callus (EC), incomplete compact pro-embryogenic cultures (ICpEC) and globular embryos (GE) in an early SE system. At a global level, the DNA 5-methylcytosine content in EC, ICpEC and GE was 24.59, 19.65 and 19.74%, respectively, suggesting a global decrease in DNA methylation from EC to ICpEC and then a slight increase from ICpEC to GE. Differentially methylated region (DMR) analysis showed that hypomethylation mainly occurred in CHH contexts. Gene ontology and Kyoto encyclopedia of genes and genomes analysis of hypomethylated-CHH-DMR-associated genes revealed that zein biosynthesis, fatty acid biosynthesis, circadian rhythm and mitophagy pathways were involved in longan early SE. Expression patterns of DNA methyltransferase and demethylase genes during longan early SE suggested that the decrease in DNA methylation was probably regulated by DNA methyltransferase genes and the DNA demethylase gene REPRESSOR OF SILENCING 1 (ROS1). The correlation between DNA hypomethylation and gene expression revealed that decreased DNA methylation did not cause extensive changes in gene expression during early longan SE and that gene expression may be affected by methylation changes in gene and downstream regions. Inhibiting DNA methylation with 5-azacytidine treatment in EC promoted the formation of GE and enhanced the capability of longan SE. Our results suggest that DNA demethylation has important roles in longan SE development.

Global scale transcriptome analysis reveals differentially expressed genes involve in early somatic embryogenesis in Dimocarpus longan Lour.

BACKGROUND: Somatic embryogenesis (SE) is a process of somatic cells that dedifferentiate to totipotent embryonic stem cells and generate embryos in vitro. Longan SE has been established and wildly used as model system for studying embryogenesis in woody plants, SE-related genes had been characterized. In spite of that, a comprehensive overview of SE at a molecular level is still absent. To understand the molecular mechanisms during longan SE, we examined the transcriptome changes by using Illumina HiSeq from the four distinct developmental stages, including non-embryogenic callus (NEC), embryogenic callus (EC), incomplete compact pro-embryogenic cultures (ICpEC), globular embryos (GE). RESULTS: RNA-seq of the four samples generated a total of 243.78 million high quality reads, approximately 81.5% of the data were mapped to longan genome. The cDNA libraries of NEC, EC, ICpEC and GE, generated 22,743, 19,745, 21,144, 21,102 expressed transcripts, 1935, 1710, 1816, 1732 novel transcripts, 2645, 366, 505, 588 unique genes, respectively. Comparative transcriptome analysis showed that a total of 10,642, 4180, 5846 and 1785 genes were differentially expressed in the pairwise comparisons of NEC_vs_EC, EC_vs_ICpEC, EC_vs_GE, ICpEC_vs_GE, respectively. Among them, plant hormones signalling related genes were significantly enriched, especially the auxin and cytokinin signalling components. The transcripts of flavonoid biosynthesis related genes were mainly expressed in NEC, while fatty acid biosynthesis related genes mainly accumulated in early SE. In addition, the extracelluar protein encoding genes LTP, CHI, GLP, AGP, EP1 were related to longan SE. Combined with the FPKM value of longan nine tissues transcription, 27 SE specific or preferential genes (LEC1, LEC1-like, PDF1.3, GH3.6, AGL80, PIN1, BBM, WOX9, WOX2, ABI3, et al.) and 28 NEC preferential genes (LEA5, CNOT3, DC2.15, PR1-1, NsLTP2, DIR1, PIP1, PIP2.1, TIP2-1, POD-P7 and POD5 et al.) were characterized as molecular markers for longan early SE. qRT-PCR validation of SE-related genes showed a high correlation between RNA-seq and qRT-PCR data. CONCLUSION: This study provides new insights into the role of the transcriptome during early SE in longan. Differentially expressed genes reveal that plant hormones signalling, flavonoid and fatty acid biosynthesis, and extracelluar protein related genes were involved in longan early SE. It could serve as a valuable platform resource for further functional studies addressing embryogenesis in woody plants.

Plant Cytosolic Ascorbate Peroxidase with Dual Catalytic Activity Modulates Abiotic Stress Tolerances.

Ascorbic acid-glutathione (AsA-GSH) cycle represents important antioxidant defense system in planta. Here we utilized Oncidium cytosolic ascorbate peroxidase (OgCytAPX) as a model to demonstrate that CytAPX of several plants possess dual catalytic activity of both AsA and GSH, compared with the monocatalytic activity of Arabidopsis APX (AtCytAPX). Structural modeling and site-directed mutagenesis identified that three amino acid residues, Pro63, Asp75, and Tyr97, are required for oxidization of GSH in dual substrate catalytic type. Enzyme kinetic study suggested that AsA and GSH active sites are distinctly located in cytosolic APX structure. Isothermal titration calorimetric and UV-visible analysis confirmed that cytosolic APX is a heme-containing protein, which catalyzes glutathione in addition to ascorbate. Biochemical and physiological evidences of transgenic Arabidopsis overexpressing OgCytAPX1 exhibits efficient reactive oxygen species-scavenging activity, salt and heat tolerances, and early flowering, compared with Arabidopsis overexpressing AtCytAPX. Thus results on dual activity CytAPX impose significant advantage on evolutionary adaptive mechanism in planta.

RNA-sequencing analysis reveals betalains metabolism in the leaf of Amaranthus tricolor L.

Amaranth plants contain large amounts of betalains, including betaxanthins and betacyanins. Amaranthin is a betacyanin, and its molecular structure and associated metabolic pathway differ from those of betanin in beet plants. The chlorophyll, carotenoid, betalain, and flavonoid contents in amaranth leaves were analyzed. The abundance of betalain, betacyanin, and betaxanthin was 2-5-fold higher in the red leaf sectors than in the green leaf sectors. Moreover, a transcriptome database was constructed for the red and green sectors of amaranth leaves harvested from 30-day-old seedlings. 22 unigenes were selected to analyze the expression profiles in the two leaf sectors. The RNA-sequencing data indicated that many unigenes are involved in betalain metabolic pathways. The potential relationships between diverse metabolic pathways and betalain metabolism were analyzed. The validation of the expression of 22 selected unigenes in a qRT-PCR assay revealed the genes that were differentially expressed in the two leaf sectors. Betalains were biosynthesized in specific tissues of the red sectors of amaranth leaves. Almost all of the genes related to betalain metabolism were identified in the transcriptome database, and the expression profiles were different between the red sectors and green sectors in the leaf. Amaranth plants consist of diverse metabolic pathways, and the betalain metabolic pathway is linked to a group of other metabolic pathways.

Genome-wide identification and characterization of long non-coding RNAs involved in the early somatic embryogenesis in Dimocarpus longan Lour.

BACKGROUND: Long non-coding RNAs (lncRNAs) are involved in variable cleavage, transcriptional interference, regulation of DNA methylation and protein modification. However, the regulation of lncRNAs in plant somatic embryos remains unclear. The longan (Dimocarpus longan) somatic embryogenesis (SE) system is a good system for research on longan embryo development. RESULTS: In this study, 7643 lncRNAs obtained during early SE in D. longan were identified by high-throughput sequencing, among which 6005 lncRNAs were expressed. Of the expressed lncRNAs, 4790 were found in all samples and 160 were specifically expressed in embryogenic callus (EC), 154 in incomplete embryogenic compact structures (ICpECs), and 376 in globular embryos (GEs). We annotated the 6005 expressed lncRNAs, and 1404 lncRNAs belonged to 506 noncoding RNA (ncRNA) families and 4682 lncRNAs were predicted to target protein-coding genes. The target genes included 5051 cis-regulated target genes (5712 pairs) and 1605 trans-regulated target genes (3618 pairs). KEGG analysis revealed that most of the differentially expressed target genes (mRNAs) of the lncRNAs were enriched in the "plant-pathogen interaction" and "plant hormone signaling" pathways during early longan SE. Real-time quantitative PCR confirmed that 20 selected lncRNAs showed significant differences in expression and that five lncRNAs were related to auxin response factors. Compared with the FPKM expression trends, 16 lncRNA expression trends were the same in qPCR. In lncRNA-miRNA-mRNA relationship prediction, 40 lncRNAs were predicted to function as eTMs for 15 miRNAs and 7 lncRNAs were identified as potential miRNA precursors. In addition, we verified the lncRNA-miRNA-mRNA regulatory relationships by transient expression of miRNAs (miR172a, miR159a.1 and miR398a). CONCLUSION: Analyses of lncRNAs during early longan SE showed that differentially expressed lncRNAs were involved in expression regulation at each SE stage, and may form a regulatory network with miRNAs and mRNAs. These findings provide new insights into lncRNAs and lay a foundation for future functional analysis of lncRNAs during early longan SE.

Seed-Specific Gene MOTHER of FT and TFL1（MFT） Involved in Embryogenesis, Hormones and Stress Responses in Dimocarpus longan Lour.

Mother of FT and TFL1 (MFT) belongs to phosphatidylethanolamine-binding protein (PEBP) family, which plays an important role in flowering time regulation, seed development, and germination. To gain insight into the molecular function of DlMFT in Dimocarpus longan Lour., we isolated DlMFT and its promoter sequence from longan embryogenic callus (EC). Bioinformatic analysis indicated that the promoter contained multiphytohormones and light responsive regulatory elements. Subcellular localization showed that the given the DlMFT signal localized in the nucleus, expression profiling implied that DlMFT showed significant upregulation during somatic embryogenesis (SE) and zygotic embryogenesis (ZE), and particular highly expressed in late or maturation stages. The accumulation of DlMFT was mainly detected in mature fruit and seed, while it was undetected in abortive seeds, and notably decreased during seed germination. DlMFT responded differentially to exogenous hormones in longan EC. Auxins, salicylic acid (SA) and methyl jasmonate (MeJa) suppressed its expression, however, abscisic acid (ABA), brassinosteroids (BR) showed the opposite function. Meanwhile, DlMFT differentially responded to various abiotic stresses. Our study revealed that DlMFT might be a key regulator of longan somatic and zygotic embryo development, and in seed germination, it is involved in complex plant hormones and abiotic stress signaling pathways.

Genome-wide identification and characterization of mRNAs and lncRNAs involved in cold stress in the wild banana (Musa itinerans).

Cold stress seriously affects banana growth, yield and fruit quality. Long noncoding RNAs (lncRNAs) have been demonstrated as key regulators of biotic and abiotic stress in plants, but the identification and prediction of cold responsive mRNAs and lncRNAs in wild banana remains unexplored. In present study, a cold resistant wild banana line from China was used to profile the cold-responsive mRNAs and lncRNAs by RNA-seq under cold stress conditions, i.e. 13°C (critical growth temperature), 4°C (chilling temperature), 0°C (freezing temperature) and normal growing condition, i.e. 28°C (control group). A total of 12,462 lncRNAs were identified in cold-stressed wild banana. In mRNA, much more alternative splicing events occurred in wild banana under the cold stress conditions compared with that in the normal growing condition. The GO analysis of differential expression genes (DEGs) showed the biochemical processes and membrane related genes responded positively to the cold stress. The KEGG pathway enrichment analysis of the DEGs showed that the pathways of photosynthesis, photosynthesis-antenna proteins, circadian rhythm-plant, glutathione metabolism, starch and sucrose metabolism, cutin/suberine/biosynthesis were altered or affected by the cold stress conditions. Our analyses of the generated transcriptome and lncRNAs provide new insights into regulating expression of genes and lncRNAs that respond to cold stress in the wild banana.

Identification of miRNAs Affecting the Establishment of Brassica Alboglabra Seedling.

MicroRNAs (miRNAs) are important for plant development including seed formation, dormancy, and germination, as well as seedling establishment. The Brassica vegetable seedling establishment stage influences the development of high quality seedlings, but also affects the nutrient content of sprouts. Chinese kale (Brassica alboglabra) seedlings at different growth stages were used to construct two small-RNA (sRNA) libraries. We comprehensively analyzed the miRNAs in 2- and 9-day-old seedlings. An average of 11,722,490 clean reads were generated after removing low-quality reads and adapter contaminants. The results revealed that 37.65 and 26.69% of the sRNAs in 2- and 9-day-old seedlings, respectively, were 24 nt long. In total, 254 known mature miRNA sequences from 228 miRNA families and 343 novel miRNAs were identified. Of these miRNAs, 224 were differentially expressed between the two analyzed libraries. The most abundant miRNAs identified by sequence homology were miR156, miR167, and miR157, each with more than 100,000 sequenced reads. Compared with the expression levels in 2-day-old seedlings, MiR8154 and miR390 were the most up- and down-regulated miRNAs respectively in 9-day-old seedlings. Gene ontology enrichment analysis of the differentially expressed-miRNA target genes affecting biological processes revealed that most genes were in the "regulation of transcription" category. Additionally, the expression patterns of some miRNAs and target genes were validated by quantitative real-time polymerase chain reaction. We determined that development-associated miRNAs (e.g., bal-miR156/157/159/166/167/172/396), were highly-expressed during seedling-establishment stage, as were stress-related (bal-miR408) and metabolism-related (bal-miR826) miRNAs. Combined with the low level of targets SPL9 and AP2, it was concluded that miR156-SPL9 and miR172-AP modules play key roles during the B. alboglabra seedling establishment stage.

Comparative Transcriptome Analyses Reveal a Special Glucosinolate Metabolism Mechanism in Brassica alboglabra Sprouts.

Brassica sprouts contain abundant phytochemicals, especially glucosinolates (GSs). Various methods have been used to enhance GS content in sprouts. However, the molecular basis of GS metabolism in sprouts remains an open question. Here we employed RNA-seq analysis to compare the transcriptomes of high-GS (JL-08) and low-GS (JL-09) Brassica alboglabra sprouts. Paired-end Illumina RNA-seq reads were generated and mapped to the Brassica oleracea reference genome. The differentially expressed genes were analyzed between JL-08 and JL-09. Among these, 1477 genes were up-regulated and 1239 down-regulated in JL-09 compared with JL-08. Enrichment analysis of these differentially expressed genes showed that the GS biosynthesis had the smallest enrichment factor and the highest Q-value of all metabolic pathways in Kyoto Encyclopedia of Genes and Genomes database, indicating the main metabolic difference between JL-08 and JL-09 is the GS biosynthetic pathway. Thirty-seven genes of the sequenced data were annotated as putatively involved in GS biosynthesis, degradation, and regulation, of which 11 were differentially expressed in JL-08 and JL-09. The expression level of GS degradation enzyme myrosinase in high-GS JL-08 was lower compared with low-GS JL-09. Surprisingly, in high-GS JL-08, the expression levels of GS biosynthesis genes were also lower than those in low-GS JL-09. As the GS contents in sprouts are determined by dynamic equilibrium of seed stored GS mobilization, de novo synthesis, degradation, and extra transport, the result of this study leads us to suggest that efforts to increase GS content should focus on either raising GS content in seeds or decreasing myrosinase activity, rather than improving the expression level of GS biosynthesis genes in sprouts.