First Report of Kiwifruit Root Rot Caused by Globisporangium spinosum in China.

Kiwifruit is widely cultivated for its high vitamin C content and nutritional value. In January 2022, root rot symptoms were found in about 30% of Actinidia chinensis cv. Jinyan plants grafted on A. deliciosa rootstocks in an orchard located in Sanming (26.32°N, 117.23°E), Fujian Province of China. The affected plants appeared stunted, with brown and decaying roots, some of which were covered with white hyphae. To isolate the pathogen, the surfaces of typical symptomatic roots were sterilized for 30 s using 75% ethanol, followed by four rinses in sterile water, placing on potato dextrose agar (PDA), and incubating away from light at 25°C for 7 days. 16 Globisporangium-like isolates were obtained through hyphal tip isolation, displaying a milky-white appearance with irregular protuberances on the surface, and yellow-white backs with radial fold lines. The isolates were then cultured on corn meal agar for 5 days at 25°C in dark for morphological characteristics. Under microscope, the hyphae appeared as long strips without septa and 4.1 to 8.2 µm wide (average 6.7 µm), containing irregularly sized spherical droplets. Both terminal and intercalary hyphae swellings were observed; these appeared either spherical or subspherical, with some having projections. Their dimensions were 12.3 to 27.6 µm (average 17.3 µm). The oospores were mostly spherical, either plerotic or aplerotic, 11.8 to 22.3 µm wide (average 18.9 µm), with occasional projections. The antheridia were rod-shaped and curved, with one end attached to the oogonia. Amplification of the sequences of internal transcribed spacer (ITS) regions and cytochrome c oxidase subunit I (COI) were conducted using the primers ITS1/ITS4 (White et al. 1990) and OomCoxI-Levlo/OomCoxI-Levup (Robideau et al. 2011), respectively. The sequencing results revealed identical ITS and COI sequences in all 16 isolates. BLASTn analysis of the 969-bp ITS sequence ON202808 showed 99.38-99.59% similarity (965/971bp, 967/971bp) with the KJ162353 and AY598701 sequences from Globisporangium spinosum isolates, while the 700-bp COI sequence ON075783 showed 100% and 99.41% identity (680/680bp, 676/680bp) with the GenBank sequences HQ708835 and HQ708832, respectively, from G. spinosum. Phylogenetic analysis also showed that the obtained isolate (termed MA16) clustered with isolates from G. spinosum on the same evolutionary branch. For pathogenicity testing, four-month-old healthy Jinyan (A. chinensis) plants grown in sterilized media were transferred to sterile petri dishes covered with wet filter paper, and their roots were inoculated with a 5-mm-wide disk of MA16 when cultivated on PDA medium for 5 days. Miliang-1 (A. deliciosa) and Hongyang (A. chinensis) plants were treated similarly. The control groups each included three plants that were inoculated with non-colonized PDA. The plants were kept at 25 °C with a 12-/12-h light/dark cycle for 10 days when the inoculated plants exhibited root rot symptoms similar to those seen in the field, together with rotting and browning of the leaves. The control plants appeared healthy with no symptoms. After re-isolated from infected tissues, the pathogen was verified to be G. spinosum according to its ITS sequence, thus fulfilling the Koch's postulates. Recently, Pythium spinosum has been classified as G. spinosum according to whole-genome sequencing and phylogenomic analysis (Nguyen et al. 2022). Based on the morphological features and pathogenicity results, MA16 was identified as G. spinosum (van der Plaats-Niterink 1981; Huo et al. 2023). This report appears to be the first description of kiwifruit root rots caused by G. spinosum in China, and its identification will assist the development of strategies to counteract the disease.

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 Identification of the MPK Gene Family and Expression Analysis under Low-Temperature Stress in the Banana.

Mitogen-activated protein kinases (MAPKs and MPKs) are important in the process of resisting plant stress. In this study, 21, 12, 18, 16, and 10 MPKs were identified from Musa acuminata, Musa balbisiana, Musa itinerans, Musa schizocarpa, and Musa textilis, respectively. These MPKs were divided into Group A, B, C, and D. Phylogenetic analysis revealed that this difference in number was due to the gene shrinkage of the Group B subfamily of Musa balbisiana and Musa textilis. KEGG annotations revealed that K14512, which is involved in plant hormone signal transduction and the plant-pathogen interaction, was the most conserved pathway of the MPKs. The results of promoter cis-acting element prediction and focTR4 (Fusarium oxysporum f. sp. cubense tropical race 4) transcriptome expression analysis preliminarily confirmed that MPKs were relevant to plant hormone and biotic stress, respectively. The expression of MPKs in Group A was significantly upregulated at 4 °C, and dramatically, the MPKs in the root were affected by low temperature. miR172, miR319, miR395, miR398, and miR399 may be the miRNAs that regulate MPKs during low-temperature stress, with miR172 being the most critical. miRNA prediction and qRT-PCR results indicated that miR172 may negatively regulate MPKs. Therefore, we deduced that MPKs might coordinate with miR172 to participate in the process of the resistance to low-temperature stress in the roots of the banana. This study will provide a theoretical basis for further analysis of the mechanism of MPKs under low-temperature stress of bananas, and this study could be applied to molecular breeding of bananas in the future.

Waterlogging Tolerance of Actinidia valvata Dunn Is Associated with High Activities of Pyruvate Decarboxylase, Alcohol Dehydrogenase and Antioxidant Enzymes.

Kiwifruit (Actinidia spp.) is susceptible to waterlogging stress. Although abundant wild germplasm resources exist among Actinidia plants for improving the waterlogging tolerance of kiwifruit cultivars, the underlying mechanisms remain largely unknown. Here, a comparative study was undertaken using one wild germplasm, Maorenshen (A. valvata Dunn, MRS), and one cultivar, Miliang-1 (A. chinensis var. deliciosa (A.Chev.) A.Chev. cv. Miliang-1, ML). Under stress, the ML plantlets were seriously damaged with wilted chlorotic leaves and blackened rotten roots, whereas the symptoms of injury in the MRS plantlets were much fewer, along with higher photosynthetic rates, chlorophyll fluorescence characteristics and root activity under stress conditions. However, neither aerenchyma in the root nor adventitious roots appeared in both germplasms upon stress exposure. The activities of pyruvate decarboxylase (PDC) and alcohol dehydrogenase (ADH), as well as their transcript levels, were constitutively higher in MRS than those in ML under both normal and stress conditions. Waterlogging stress significantly enhanced the PDC and ADH enzyme activities in both germplasms, which were 60.8% and 22.4% higher in the MRS roots than those in the ML roots under waterlogging stress, respectively. Moreover, MRS displayed higher activities of antioxidant enzymes, including SOD, CAT, and APX, as well as DPPH-radical scavenging ability, and decreased H2O2 and MDA accumulation under both normal and stress conditions. Our findings suggest that the waterlogging tolerance of the wild A. valvata germplasm was associated with high PDC and ADH, as well as antioxidant ability.

Identification and Characterization of Long Non-Coding RNAs: Implicating Insights into Their Regulatory Role in Kiwifruit Ripening and Softening during Low-Temperature Storage.

Long non-coding RNAs (lncRNAs) are crucial players regulating many biological processes in plants. However, limited knowledge is available regarding their roles in kiwifruit ripening and softening. In this study, using lncRNA-seq technology, 591 differentially expressed (DE) lncRNAs (DELs) and 3107 DE genes (DEGs) were identified from kiwifruit stored at 4 °C for 1, 2, and 3 weeks in comparison with non-treated control fruits. Of note, 645 DEGs were predicted to be targets of DELs (DEGTLs), including some DE protein-coding genes (such as ß-amylase and pectinesterase). DEGTL-based GO enrichment analysis revealed that these genes were significantly enriched in cell wall modification and pectinesterase activity in 1 W vs. CK and 3 W vs. CK, which might be closely related to the fruit softening during low-temperature storage. Moreover, KEGG enrichment analysis revealed that DEGTLs were significantly associated with starch and sucrose metabolism. Our study revealed that lncRNAs play critical regulatory roles in kiwifruit ripening and softening under low-temperature storage, mainly by mediating the expression of starch and sucrose metabolism and cell wall modification related genes.

Integrated Metabolomic and Transcriptomic Analysis Reveals Differential Flavonoid Accumulation and Its Underlying Mechanism in Fruits of Distinct Canarium album Cultivars.

Canarium album fruit has great potential to be consumed as a raw material not only for food but also medicine. The diverse active metabolites composition and content of C. album fruits greatly affect their pharmacological effects. However, up to now, there has been no report on the global metabolome differences among fruits from distinct C. album cultivars. In our present study, by using non-targeted metabolomics techniques, we identified 87 DAMs (differentially accumulated metabolites) including 17 types of flavonoids from fruits of four different C. album cultivars. KEGG (Kyoto Encyclopedia of Genes and Genomes) pathway enrichment analysis revealed that the flavone and flavonol biosynthesis- and flavonoid biosynthesis-related DAMs were major factors determining their metabolome differences. Comparative transcriptomic analysis revealed that 15 KEGG pathways were significantly enriched by genes of the identified 3655 DEGs (differentially expressed genes) among different C. album cultivars. Consistent with the metabolome data, flavonoid biosynthesis-related DEGs, including eight key structural genes (such as FLS, CCoAOMT, CHI, C4H, DFR, LAR, and C3'H, etc.) and several regulatory transcription factor (TF) genes (including 32 MYBs and 34 bHLHs, etc.), were found to be significantly enriched (p < 0.01). Our study indicated that the differential expression of flavonoid biosynthesis-related genes and accumulation of flavonoids played dominant roles in the various metabolome compositions of fruits from different C. album cultivars.

De novo transcriptome assembly and comparative transcriptomic analysis provide molecular insights into low temperature stress response of Canarium album.

A de novo transcriptome analysis was performed in C. album, a temperature sensitive fruit tree in China, after treatment with varied temperatures. A total number of 168,385 transcripts were assembled, comprising of 109,439 unigenes, of which 70,530 were successfully annotated. Compared with control check group (CK), which was treated under 25 °C, the chilling stress (4 °C) treated group (CT), showed about 2810 up-regulated and 2567 down-regulated genes. Whereas, group treated under freezing (- 3 °C) stress (FT) showed an up-regulation and a down-regulation of 1748 and 1459 genes, respectively. GO classification analysis revealed that DEGs related to metabolic processes, single-organism metabolic process, and catalytic activity are significantly enriched in both CT and FT conditions. KEGG pathway enrichment analysis for both CT and FT treatments showed an enrichment of genes encoding or related to glycine/serine and threonine metabolism, alpha-linolenic acid metabolism, carotenoid biosynthesis, photosynthesis-antenna proteins, and circadian rhythm. However, genes related to photosynthesis, carbon fixation in photosynthetic organisms, glutathione metabolism, pyruvate metabolism, nicotinate and nicotinamide metabolism were specifically enriched in CT condition. Nevertheless, FT treatment induced genes related to plant-pathogen interaction, linoleic acid metabolism, plant hormone signal transduction and pentose phosphate pathway. Many of the genes involved in plant hormone signal transduction showed significantly different expression in both FT and CT conditions. However, the change was more evident in FT. Here we present the first of the reports for a de novo transcriptomic analysis in C. album, suggesting that the plant shows differential responses in chilling and freezing temperatures, where the hormone signaling and transduction contribute greatly to FT responses. Our study thus paves way for future research regarding functions of these potentially identified genes.

Integrated analysis of lncRNA and mRNA transcriptomes reveals the potential regulatory role of lncRNA in kiwifruit ripening and softening.

Kiwifruit has gained increasing attention worldwide for its unique flavor and high nutritional value. Rapid softening after harvest greatly shortens its shelf-life and reduces the commercial value. Therefore, it is imperative and urgent to identify and clarify its softening mechanism. This study aimed to analyze and compare the long noncoding RNA (lncRNA) and mRNA expression patterns in ABA-treated (ABA) and room temperature (RT)-stored fruits with those in freshly harvested fruits (CK) as control. A total of 697 differentially expressed genes (DEGs) and 81 differentially expressed lncRNAs (DELs) were identified while comparing ABA with CK, and 458 DEGs and 143 DELs were detected while comparing RT with CK. The Kyoto Encyclopedia of Genes and Genomes analysis of the identified DEGs and the target genes of DELs revealed that genes involved in starch and sucrose metabolism, brassinosteroid biosynthesis, plant hormone signal transduction, and flavonoid biosynthesis accounted for a large part. The co-localization networks, including 38 DEGs and 31 DELs in ABA vs. CK, and 25 DEGs and 25 DELs in RT vs. CK, were also performed. Genes related to fruit ripening, such as genes encoding ß-galactosidase, mannan endo-1,4-ß-mannosidase, pectinesterase/pectinesterase inhibitor, and NAC transcription factor, were present in the co-localization network, suggesting that lncRNAs were involved in regulating kiwifruit ripening. Notably, several ethylene biosynthesis- and signaling-related genes, including one 1-aminocyclopropane-1-carboxylic acid oxidase gene and three ethylene response factor genes, were found in the co-localization network of ABA vs. CK, suggesting that the promoting effect of ABA on ethylene biosynthesis and fruit softening might be embodied by increasing the expression of these lncRNAs. These results may help understand the regulatory mechanism of lncRNAs in ripening and ABA-induced fruit softening of kiwifruit.

Characterization and complete genome sequences of two novel variants of the family Closteroviridae from Chinese kiwifruit.

Two distinct closterovirus-like genome sequences (termed AdV-1 v1 and v2) were identified in Actinidia chinensis var. deliciosa 'Miliang-1' that had no disease symptoms using high-throughput sequencing. Using overlapping reverse transcription-polymerase chain reaction and rapid amplification of cDNA ends, the genomic sequences of AdV-1 v1 and v2 were confirmed as 17,646 and 18,578 nucleotides in length, respectively. The two complete genomes contained 9 and 15 open reading frames, respectively, coding for proteins having domains typical of Closteroviridae, such as RNA-dependent RNA polymerase (RdRp), heat shock protein 70 homolog (HSP70h) and coat protein (CP). Sequence analysis showed that the amino acid sequences of RdRp, HSP70h, and CP of the two variants exhibited high similarity (> 80%), while their genomic organization was somewhat different. This suggested that the two viral genomes identified here are variants of the family Closteroviridae in a single kiwifruit host. Furthermore, phylogenetic relationship analysis revealed that the two variants had a closer relationship with the unclassified virus Persimmon virus B (PeVB) and Actinidia virus 1 (AcV-1) than with other members of the family Closteroviridae, as did their genomic organization. It is speculated that the two variants, together with PeVB and AcV-1 belong to a new subfamily of Closteroviridae.

The complete chloroplast genome sequence of actinidia valvata.

In this study, we first presented the complete chloroplast genome of Actinidia valvata by using Illumina Novaseq sequencing. Its complete chloroplast genome is 156,596 bp in length, containing a large single copy region of 88,477 bp and a small single copy region of 20,379 bp separated by a pair of inverted repeat regions of 23,870 bp. The chloroplast genome contains 112 unique genes, including 78 protein-coding genes, 30 tRNA, and four rRNA genes. Phylogenetic analysis based on chloroplast genome sequences of ten plants from the family Actinidiaceae showed that A. valvata is more closely related to A. polygama than other members.

The complete chloroplast genome characterization and phylogenetic analysis of Canarium album.

Canarium album is one of the precious and characteristic fruit trees of China. In this study, we first presented the complete chloroplast genome of C. album by using BGISEQ-500 sequencing. Its complete chloroplast genome is 163,140 bp in size, containing a pair of inverted repeat regions of 30,729 bp, a large single copy region of 87,748 bp and a small single copy region of 13,934 bp. The chloroplast genome contains 117 unique genes, including 83 protein coding genes, 30 tRNA and 4 rRNA genes. Phylogenetic analysis based on complete chloroplast genomes indicated that C. album was closest to Boswellia sacra.

Genome-wide sequencing of longan (Dimocarpus longan Lour.) provides insights into molecular basis of its polyphenol-rich characteristics.

Longan (Dimocarpus longan Lour.), an important subtropical fruit in the family Sapindaceae, is grown in more than 10 countries. Longan is an edible drupe fruit and a source of traditional medicine with polyphenol-rich traits. Tree size, alternate bearing, and witches' broom disease still pose serious problems. To gain insights into the genomic basis of longan traits, a draft genome sequence was assembled. The draft genome (about 471.88 Mb) of a Chinese longan cultivar, "Honghezi," was estimated to contain 31 007 genes and 261.88 Mb of repetitive sequences. No recent whole-genome-wide duplication event was detected in the genome. Whole-genome resequencing and analysis of 13 cultivated D. longan accessions revealed the extent of genetic diversity. Comparative transcriptome studies combined with genome-wide analysis revealed polyphenol-rich and pathogen resistance characteristics. Genes involved in secondary metabolism, especially those from significantly expanded (DHS, SDH, F3΄H, ANR, and UFGT) and contracted (PAL, CHS, and F3΄5΄H) gene families with tissue-specific expression, may be important contributors to the high accumulation levels of polyphenolic compounds observed in longan fruit. The high number of genes encoding nucleotide-binding site leucine-rich repeat (NBS-LRR) and leucine-rich repeat receptor-like kinase proteins, as well as the recent expansion and contraction of the NBS-LRR family, suggested a genomic basis for resistance to insects, fungus, and bacteria in this fruit tree. These data provide insights into the evolution and diversity of the longan genome. The comparative genomic and transcriptome analyses provided information about longan-specific traits, particularly genes involved in its polyphenol-rich and pathogen resistance characteristics.

Ras-Related Nuclear Protein Ran3B Gene Is Involved in Hormone Responses in the Embryogenic Callus of Dimocarpus longan Lour.

Ras-related guanosine triphosphate (GTP)-binding nuclear protein (Ran) GTPases function as molecular switches and regulate diverse cellular events in eukaryotes. Our previous work suggested that DlRan3B is active during longan (Dimocarpus longan Lour.) somatic embryogenesis (SE) processes. Herein, subcellular localization of DlRan3B was found to be localized in the nucleus and expression profiling of DlRan3B was performed during longan SE and after exposure to plant hormones (indoleacetic acid (IAA), gibberellin A3 (GA3), salicylic acid (SA), methyl jasmonte (MeJA), and abscisic acid (ABA)). We cloned and sequenced 1569 bp of 5'-flanking sequence of DlRan3B (GenBank: JQ279697). Bioinformatic analysis indicated that the promoter contained plant hormone-related regulatory elements. Deletion analysis and responses to hormones identified stimulative and repressive regulatory elements in the DlRan3B promoter. The key elements included those responding to auxin, gibberellin, SA, MeJA, and ABA. DlRan3B was located in the nucleus and accumulated in the late stage of longan SE. The expression of DlRan3B was significantly induced by IAA, GA3, and ABA, but suppressed by SA and MeJA. Promoter transcription was induced by IAA and GA3, but suppressed by SA. Thus, DlRan3B might participate in auxin, gibberellin, and ABA responses during longan late SE, and DlRan3B is involved in phytohormone responsiveness.

Endogenous target mimics down-regulate miR160 mediation of ARF10, -16, and -17 cleavage during somatic embryogenesis in Dimocarpus longan Lour.

MicroRNA160 plays a critical role in plant development by negatively regulating the auxin response factors ARF10, -16, and -17. However, the ways in which miR160 expression is regulated at the transcriptional level, and how miR160 interacts with its targets during plant embryo development, remain unknown. Here, we studied the regulatory relationships among endogenous target mimics (eTMs), and miR160 and its targets, and their involvement in hormone signaling and somatic embryogenesis (SE) in Dimocarpus longan. We identified miR160 family members and isolated the miR160 precursor, primary transcript, and promoter. The promoter contained cis-acting elements responsive to stimuli such as light, abscisic acid, salicylic acid (SA) and heat stress. The pri-miR160 was down-regulated in response to SA but up-regulated by gibberellic acid, ethylene, and methyl jasmonate treatment, suggesting that pri-miR160 was associated with hormone transduction. Dlo-miR160a, -a(∗) and -d(∗) reached expression peaks in torpedo-shaped embryos, globular embryos and cotyledonary embryos, respectively, but were barely detectable in friable-embryogenic callus. This suggests that they have expression-related and functional diversity, especially during the middle and later developmental stages of SE. Four potential eTMs for miR160 were identified. Two of them, glucan endo-1,3-beta- glucosidase-like protein 2-like and calpain-type cysteine protease DEK1, were confirmed to control the corresponding dlo-miR160a(∗) expression level. This suggests that they may function to abolish the binding between dlo-miR160a(∗) and its targets. These two eTMs also participated in 2,4-D and ABA signal transduction. DlARF10, -16, and -17 targeting by dlo-miR160a was confirmed; their expression levels were higher in friable-embryogenic callus and incomplete compact pro-embryogenic cultures and responded to 2,4-D, suggesting they may play a major role in the early stages of longan SE dependent on 2,4-D. The eTMs, miR160, and ARF10, -16, and -17 exhibited tissue specificity in 'Sijimi' longan vegetative and reproductive organs, but were not significant negatively correlated. These results provide insights into the possible role of the eTM-miR160-ARF10-16-17 pathway in longan somatic embryo development.

DlRan3A is involved in hormone, light, and abiotic stress responses in embryogenic callus of Dimocarpus longan Lour.

Ras-related nuclear protein (Ran) GTPase plays an important role in nucleo-cytoplasmic transportation of proteins and RNA, mitotic spindle assembly, microtubule assembly and nuclear envelope (NE) assembly. We previously identified the full-length cDNAs and a DNA of DlRan3A from longan (Dimocarpus longan Lour.) somatic embryos and demonstrated its possible roles in cell activities during longan somatic embryogenesis (SE). However, thus far little is known of how Ran functions in the signaling pathways in plant embryos in response to changing environmental stimuli. To discover more biological roles of DlRan3A, we observed DlRan3A located in the nucleus and detected abundant accumulation of DlRan3As in globular and cotyledon embryos during longan SE, which suggested its involvement in auxin signal pathways in longan early SE. The transcript level of the DlRan3A gene were also determined in longan embryogenic callus (EC) in response to different levels of exogenous phytohormones (indoleacetic acid (IAA), gibberellin A3 (GA3), salicylic acid (SA), methyl jasmonate (MeJA), and ethephon (Eth)), salt, sucrose, prolonged culturing period and light quality treatments. IAA (1mg/L), GA3 (12 mg/L), SA (75 µm), MeJA (50 and 100 µm), and Eth (50mg/L) modulate expression of DlRan3A to 1.3-1.4 folds, and the expression of DlRan3A was significantly affected by light quality, significantly induced to 2-fold by salt (10 g/L), 2.7-fold by sucrose (90 g/L) and was completely suppressed by prolonged cultivation (>40 days). Deletion analysis suggested that both activation and repression regulatory elements co-exist in the DlRan3A promoter sequence and that the key cis-acting elements included ones in response to auxin, SA, MeJA, and stress. Promoter activities were induced to the highest level by IAA followed by SA, GA3 and MeJA, while suppressed by ABA and Eth. Together, these results showed DlRan3A close involvement in phytohormones, light, and abiotic stress responsiveness during longan SE.

MicroRNA390-Directed TAS3 Cleavage Leads to the Production of tasiRNA-ARF3/4 During Somatic Embryogenesis in Dimocarpus longan Lour.

Trans-acting short-interfering RNAs (tasiRNAs) originate from TAS3 families through microRNA (miRNA) 390-guided cleavage of primary transcripts and target auxin response factors (ARF3/-4), which are involved in the normal development of lateral roots and flowers in plants. However, their roles in embryo development are still unclear. Here, the pathway miR390-TAS3-ARF3/-4 was identified systematically for the first time during somatic embryo development in Dimocarpus longan. We identified the miR390 primary transcript and promoter. The promoter contained cis-acting elements responsive to stimuli such as light, salicylic acid, anaerobic induction, fungal elicitor, circadian control, and heat stress. The longan TAS3 transcript, containing two miR390-binding sites, was isolated; the miR390- guided cleavage site located near the 3' end of the TAS3 transcript was verified. Eight TAS3-tasiRNAs with the 21-nucleotides phase were found among longan small RNA data, further confirming that miR390-directed TAS3 cleavage leads to the production of tasiRNA in longan. Among them, TAS3_5'D5+ and 5'D6+ tasiRNAs were highly abundant, and verified to target ARF3 and -4, implying that miR390-guided TAS3 cleavage with 21-nucleotides phase leading to the production of tasiRNA-ARF is conserved in plants. Pri-miR390 was highly expressed in friable-embryogenic callus (EC), and less expressed in incomplete compact pro-embryogenic cultures, while miR390 showed its lowest expression in EC and highest expression in torpedo-shaped embryos (TEs). DlTAS3 and DlARF4 both exhibited their lowest expressions in EC, and reached their peaks in the globular embryos stage, which were mainly inversely proportional to the expression of miR390, especially at the globular embryos to cotyledonary embryos (CEs) stages. While DlARF3 showed little variation from the EC to TEs stages, and exhibited its lowest expression in the CEs stage. There was a general lack of correlation between the expressions of DlARF3 and miR390. In addition, pri-miR390, DlTAS3, DlARF3 and -4 were up-regulated by 2,4-D in a concentration-dependent manner. They were also preferentially expressed in roots, pulp, and seeds of 'Sijimi' longan, implying their extended roles in the development of longan roots and fruit. This study provided insights into a possible role of miR390-tasiRNAs-ARF in plant somatic embryo development.