Skip to main content

Genome-wide identification of the longan R2R3-MYB gene family and its role in primary and lateral root

Abstract

R2R3-MYB is an important transcription factor family that regulates plant growth and development. Root development directly affects the absorption of water and nutrients by plants. Therefore, to understand the regulatory role of R2R3-MYB transcription factor family in root development of longan, this study identified the R2R3-MYB gene family members at the genome-wide level, and analyzed their phylogenetic characteristics, physical and chemical properties, gene structure, chromosome location and tissue expression. The analysis identified 124 R2R3-MYB family members in the longan genome. Phylogenetic analysis divided these members into 22 subfamilies, and the members of the unified subfamily had similar motifs and gene structures. The result of qRT-PCR showed that expression levels of DlMYB33, DlMYB34, DlMYB59, and DlMYB77 were significantly higher in main roots than in lateral as opposed to those of DlMYB35, DlMYB69, DlMYB70, and DlMYB83, which were significantly lower. SapBase database prediction and miRNAs sequencing results showed that 34 longan miRNAs could cleave R2R3-MYB, including 17 novel miRNAs unique to longan. The qRT-PCR and subcellular localization experiments of DlMYB92 and DlMYB98 showed that DlMYB92 is a key factor that regulates transcription in the nucleus and participates in the regulation of longan lateral root development. Longan also has a conserved miRNA-MYB-lateral root development regulation mechanism. This study provides a reference for further research on the transcriptional regulation of the miRNA-R2R3-MYB module in the root development of longan.

Peer Review reports

Introduction

Gene expression in plants is regulated by several factors at the transcriptional and post-transcriptional levels. Transcription factors (TFs) are a class of trans-acting factors that regulate gene expression at the transcriptional level [1, 2]. The MYB (v-Myb avian myeloblastosis viral oncogene homolog) transcription factor family is a large and functionally diverse family of TFs unique to plants [3], and the TFs are classified into four types based on the number of structural domains in the amino acid sequence: 1R-MYB, R2R3-MYB, 3R-MYB, and 4R-MYB. R2R3-MYB contains two tandem MYB structural domains, R2 and R3. [4]. MicroRNAs (miRNAs), a class of single-stranded endogenous non-coding RNAs ranging in size from 20 ~ 24 nt, play an important role as post-transcriptional regulators [5]. Recent research has shown that miRNAs play a critical role in plant growth and resistance to stress, and have become key molecules in understanding the molecular basis of biological activity [6,7,8,9]. Studies show that miRNAs exert their regulatory roles by targeting TFs [10, 11]. Therefore, understanding MYB regulation of target genes at the transcriptional level and miRNA regulation of MYB at the post-transcriptional level is important for insights into plant growth and development.

Root growth and development directly affect plant uptake of minerals and water [12]. While plant roots exhibit high phenotypic plasticity [13], their growth and development are regulated by multiple factors. It has been found that miRNAs and TFs can form miRNA-TFs regulatory modules to regulate plant root development in a stepwise manner. For example, MYB interacts with abscisic acid (ABA) signaling in the endodermis of Arabidopsis roots to form a cascading regulatory network of lignification and corkiness that affects lateral root development [14]. In Arabidopsis AtmmiR159ab double mutants, the expression of MYB33, MYB65, and MYB101 increased, the primary root became significantly longer, and the meristematic tissue became larger. The number of cells in the meristematic and elongation zones increased significantly compared with the wild type, indicating that miR159 may inhibit primary root growth by cleaving MYB33, MYB65, and MYB101 [15]. LAC2 (LACCASE2), a negative regulator of lignin deposition in root vascular tissue, is post-transcriptionally regulated by miR397b under water-deficit conditions, resulting in increased root length and reduced lignin content in Arabidopsis root vascular tissue [16]. Additionally, It was also demonstrated that Arabidopsis miR156 represses adventitious root development by regulating the target gene SPL (SQUAMOSA promoter binding protein-like) [17]. The miR156-SPL regulatory module was found to play a role in root development in rice [18], maize [19], alfalfa [20], and apple [21]. In summary, miRNA and MYB TFs are important regulators of plant root development. However, the role of most miRNAs in regulating root development has not been directly validated in transgenic plants due to the limitations of transgenic technology.

Longan (Dimocarpus longan) is widely cultivated in many Southeast Asian countries, including China [22, 23]. Despite the large cultivation area of longan in China, yields are low and still need improvement [24]. Irrational fertilization of land during cultivation has been shown to cause the degradation of longan roots, resulting in weaker root nutrient uptake and limiting yield improvement. In order to effectively enhance the quantity of lateral roots and their accompanying nutrient absorption capacity, a comprehensive understanding of how miRNA and MYB regulate the development of longan roots is imperative. Consequently, members of the R2R3-MYB transcription factor family in longan were screened and identified through whole genome sequencing. MicroRNA sequencing was used to analyze the post-transcriptional regulation of the R2R3-MYB family members in longan [25, 26]. This work will facilitate our further understanding of the potential of MYB TFs and miRNAs in regulating root development in longan.

Results

Identification and physicochemical characterization of the R2R3-MYB gene family in longan

The identification of R2R3-MYB gene family members in longan was carried out through an HMM search and Batch-CDD-search, targeting protein sequences with the R2R3-MYB structural domain. A total of 124 R2R3-MYB gene family members were successfully identified in the ‘Shi Xia’ longan genome and designated as DlMYB1-DlMYB124 (Table 1). The lengths of the protein sequences of the124 longan R2R3-MYB family members ranged from 134 to 661 aa, with an average length of 307 aa. The molecular masses ranged from 15.69 to 775.81 kDa, with an average molecular mass of 34.68 kDa. Theoretical isoelectric points ranged from 4.91 to 10.6. The instability coefficients ranged from 35.62 to 72.31, and only DlMYB101, DlMYB21, DlMYB91, DlMYB58, DlMYB119, DlMYB20 and DlMYB114 proteins had instability coefficients less than 40, while the rest of the 117 R2R3-MYB proteins had instability coefficients greater than 40, indicating that these proteins were unstable. The fat coefficients ranged from 50.78 to 84.16, indicating that the proteins encoded by the R2R3-MYB transcription factor gene of D. longan were thermally stable. The average coefficients of hydrophilicity ranged from -1.22 to -0.5, meaning that all longan R2R3-MYBs were hydrophilic proteins.

Table 1 Information of Dimocarpus longan R2R3-MYB gene family

Phylogenetic analysis of longan R2R3-MYB gene family members

A phylogenetic tree was constructed using the conserved amino acid sequences of R2R3-MYB gene family members from both Longan and Arabidopsis. The 124 members of the R2R3-MYB gene family in Arabidopsis were classified into 22 subfamilies, denoted by C1 to C22 (Fig. 1) [4]. Eighteen of the R2R3-MYB gene subfamilies in Arabidopsis contained some known Arabidopsis subfamilies whereas six were inconsistent with the reported classification results of Arabidopsis subfamilies. The C9 subfamily, unique to Arabidopsis, did not contain any longan R2R3-MYB genes among its members. All subfamilies had more than six members, and the largest number of members was the C8 subfamily, with 24 family members. The phylogenetic analysis results suggest a potential correlation between the subfamily classification of longan R2R3-MYB genes and their functional similarities.

Fig. 1
figure 1

Phylogenetic tree of R2R3-MYB genes in Dimocarpus longan and Arabidopsis thaliana. The inner ring black star and black circle represent the members of R2R3-MYB gene family in Arabidopsis and longan, respectively. The middle blue leaf label represents the known Arabidopsis subgroup members, and the pink leaf label represents the longan subgroup members. The outer circle is represented as different subfamilies, C1 ~ C22 are the names of different subgroups of longan R2R3-MYB, and the S group in brackets is the known Arabidopsis subfamily name, this figure was created by the first author, Xinmin Lv

Analysis of gene structures, conserved motifs and domains of longan R2R3-MYB gene family members

The analysis of 124 longan R2R3-MYB members included the examination of conserved motifs, conserved domains, and gene structures. The results of conserved motif analysis showed that all R2R3-MYB family members contained four highly conserved motifs; motif1, motif2, motif3, and motif5 (Fig. 2B). The positions of the conserved motifs were consistent with those of the conserved MYB domains, indicating that both the R2 and R3 domains of the R2R3-MYB proteins in longan are highly conserved (Fig. 2C). Additionally, it was observed that R2R3-MYB members within the same subfamily exhibited similar composition and distribution of motifs, as depicted in Fig. 2A. The conserved domain prediction results showed that all Longan R2R3-MYB members contained SANT and SNAT superfamily conserved domains (Fig. 2C). The number of introns of the longan R2R3-MYB gene family members ranged from 0 to 13, and the number of exons from 1 to 14. Moreover, a noteworthy observation was the substantial variation in the number of introns and exons among members belonging to distinct subfamilies. Conversely, members within the same subfamily exhibited greater similarity in terms of gene structure, as illustrated in Fig. 2D.

Fig. 2
figure 2

Gene structures of R2R3-MYB genes in Dimocarpus longan. A phylogenetic analysis; B conserved motif analysis; C protein conserved domain analysis; D gene structure display

Chromosomal distribution and colinearity analysis of R2R3-MYB gene family members in longan

We visualized the chromosomal distribution and colinearity of the R2R3-MYB gene family members in longan. The results showed that the members were unevenly distributed on the 15 longan chromosomes (Fig. 3). The largest number of R2R3-MYB gene family members, comprising of 15, was found on Chr15 and the smallest number, comprising of 2, was on Chr11 (Fig. 3). This indicates that the number of R2R3-MYB gene family members distributed on each chromosome did not correlate with the length of the chromosome. The colinearity analysis showed several genes duplicated in tandem among the members of the R2R3-MYB gene family in longan, such as DlMYB1, DlMYB2, DlMYB3, DlMYB5, and DlMYB6, in addition to 15 pairs of fragmented duplicated genes (Fig. 3). The above results suggest that longan R2R3-MYB genes are mainly generated by chromosomal fragment duplication, which is distant on the chromosome. An unequal interchange of alleles generates a small number of tandem duplicated genes, resulting in homologous sequence clusters.

Fig. 3
figure 3

Chromosomal location and collinearity analysis of longan R2R3-MYB family genes. Jasmine boxes represent chromosomes. Segmental duplication genes are connected with black lines

Tissue expression analysis of R2R3-MYB gene family members in longan

Multiple Experiment Comparison types of Gene Expression analysis were conducted using SapBase (Sapindaceae Genomic DataBase) (http://www.sapindaceae.com/index.html). The obtained expression values of certain longan R2R3-MYB gene family members encompassed various plant tissues, including roots, stems, leaves, flowers, fruits, and seeds. The results showed that eight members, including DlMYB69, DlMYB72, and DlMYB74, were the most highly expressed in longan flowers (Fig. 4). Thirteen members, including DlMYB6, were the most highly expressed in leaves. Thirteen members, including DlMYB7, were the most highly expressed in roots. Fourteen members, including DlMYB33, were the most highly expressed in stems. Sixteen members, including DlMYB69, were the most highly expressed in fruits. DlMYB69 had the highest expression in all tissues.

Fig. 4
figure 4

The expression of R2R3-MYB gene in different tissues of longan. The expression level differences of the longan R2R3-MYB gene family members in the roots, stems, leaves, flowers, fruits, and seeds of longan are shown

In order to gain a deeper understanding of the potential role of R2R3-MYB transcription factors (TFs) in the development of longan roots, an analysis was conducted to examine the expression of R2R3-MYB genes in both primary and lateral roots of longan. This analysis utilized root RNA-seq data from the PRJNA554213 dataset. As shown in Fig. 5A, DlMYB33 and DlMYB83 were more expressed in lateral roots than in primary roots, while DlMYB69 and DlMYB70 had the highest expression in both primary and lateral roots. To screen for R2R3-MYB that may regulate lateral root development in longan, we further differentially showed the top 20 MYB TFs with the greatest variation in the expression on the top (Fig. 5B), middle (Fig. 5C), bottom (Fig. 5D), and overall (Fig. 5E) of the main and lateral roots. The results showed that the largest difference in expression was in DlMYB33, followed by DlMYB69, DlMYB83, and DlMYB13.

Fig. 5
figure 5

The expression of R2R3-MYB gene in primary root and lateral root of longan. A Heat map of R2R3-MYB gene expression in primary roots and measured roots of longan; the bar chart of the top 20 R2R3-MYB genes with the largest difference in expression between the primary root and the upper B middle C, lower D and overall E root of longan

Prediction of miRNAs regulating R2R3-MYB in longan

In order to investigate the potential regulatory role of miRNAs on R2R3-MYB gene expression in longan, miRNA prediction was carried out using the SapBase website and miRNA sequencing data. A total of 34 miRNAs were predicted, which are likely to exert post-transcriptional regulatory effects on R2R3-MYB genes in longan. The details of these predicted miRNAs are provided in Table 2. Among all predicted miRNAs, the number of novel miRNAs was 17. Different miRNAs were found to regulate the same R2R3-MYB, and majorly DlMYB98, DlMYB92, DlMYB103, and DlMYB115 (background markers in different colors in Table 2). The corresponding homologous genes of these R2R3-MYB family members in Arabidopsis were determined through sequence comparison and analysis. The homologs of DlMYB98 and DlMYB92 were found to be GAMYB and AtMYB33, respectively, two MYB TFs reported to be closely associated with root development in Arabidopsis [15].

Table 2 Prediction of miRNAs regulating longan R2R3-MYB gene

qRT-PCR analysis of differential R2R3-MYB gene expression in longan primary and lateral roots

The expression of eleven R2R3-MYBs in the primary and lateral roots was randomly detected using qRT-PCR, as shown in Fig. 5E. It was observed that DlMYB33, DlMYB34, DlMYB59, and DlMYB77 exhibited significant upregulation in the primary roots compared to the lateral roots of longan (Fig. 6). Additionally, DlMYB45 and DlMYB52 showed higher expression in the primary roots than in the lateral roots, although the differences were not statistically significant. DlMYB35, DlMYB69, DlMYB70 and DlMYB83 were less expressed in primary roots, i.e., the present residence, than in lateral roots. Inevitably, variations in the expression levels of the same genes were observed between qRT-PCR and RNA-seq results due to inconsistencies in the test materials and normalization methods employed (Fig. 6).

Fig. 6
figure 6

Analysis of relative expression of 11 R2R3-MYB genes from Dimocarpus longan. Relative expression of A DlMYB33, B DlMYB34, C DlMYB35, D DlMYB47, E DlMYB48, F DlMYB52, G DlMYB59, H DlMYB69, I DlMYB70, J DlMYB77 and K DlMYB83 in longan primary and lateral root

Analysis of subcellular localization and expression analysis of DlMYB92

MYB33 and AtGAMYB are widely studied TFs involved in plant root development in Arabidopsis. Our results showed that in longan, DlMYB92 (AtMYB33 homolog) and DlMYB98 (AtGAMYB homolog) were regulated by multiple miRNAs and differed significantly in expression in primary and lateral roots. The expression of DlMYB92 and DlMYB98 in the primary and lateral roots of longan was further examined using qRT-PCR. The results showed that the expression of DlMYB92 and DlMYB98 in the primary roots was significantly lower than that in the lateral roots, consistent with the results in Fig. 7A, B. The subcellular localization of DlMYB92 and AtMYB33 was also predicted. As shown in Fig. 7C, most predictions showed that both DlMYB92 and its homolog in Arabidopsis, AtMYB33, were localized in the nucleus. Next, we performed a subcellular localization assay for DlMYB33 in tobacco leaves. The results showed that DlMYB92 was localized in the nucleus of tobacco cells compared with the negative control and NLS-nuclear localization signal (Fig. 7D). It is speculated that DlMYB33 plays a transcriptional regulatory role in the nucleus and may be involved in the regulation of longan lateral root development.

Fig. 7
figure 7

Expression and subcellular localization analysis of DlMYB92 and DlMYB98. Expression analysis of DlMYB92 (A) and DlMYB98 (B) in root and lateral roots of longan. C Prediction of subcellular localization of DlMYB92 and DlMYB98 in longan. D Subcellular localization of DlMYB92 in tobacco leaves

Discussion

Longan R2R3-MYB family members display evolutionary conservation along with interspecies diversity

In this study, we identified a total of 124 members belonging to the R2R3-MYB transcription factor family in longan. By comparing the number of R2R3-MYB genes in longan with other species such as jujube [27], pear [28], Arabidopsis [3], soybean [29], and pea [30], we observed quantitative variation in R2R3-MYB gene numbers among different species. This variability can be attributed to differences in the occurrence of whole genome duplication (WGD) events in plants across species. Some species may have a lower number of R2R3-MYB genes due to the loss of genes that expanded after the WGD event. Phylogenetic analysis of R2R3-MYB family members in longan and Arabidopsis revealed that 18 out of the 22 longan R2R3-MYB gene subfamilies matched with known Arabidopsis subfamilies. However, six subfamilies did not align with the reported classification of Arabidopsis subfamilies [3]. These findings suggest that functional adaptive differentiation of R2R3-MYB genes may have occurred in longan and Arabidopsis following the WGD event. Furthermore, our phylogenetic and genetic structure analyses demonstrated that members of the same R2R3-MYB subfamily in longan share similar motifs and gene structures. This pattern is consistent with previous studies conducted in other species, such as potato [31] and pepper [32]. In summary, although longan exhibits quantitative diversity in R2R3-MYB gene numbers compared to other plant species, these genes remain functionally conserved, indicating their evolutionary significance.

Multiple miRNAs orchestrate the post-transcriptional regulation of Longan R2R3-MYB family members

miRNAs are a class of single-stranded endogenous non-coding RNAs that post-transcriptionally regulate gene expression [5]. They play an important role in plant growth and development, participating in developmental processes, signal transduction, protein degradation, response to adversity stress, and pathogen invasion [6]. Numerous studies have shown that miRNAs play a leading role in post-transcriptional regulation, especially in targeting transcription factors [10, 11]. MiRNAs that potentially regulate the R2R3-MYB members in longan were predicted using online prediction tools and miRNA sequencing. The results showed that miRNA abundance obtained from online database prediction was significantly less than that from miRNA sequencing. Therefore, a single online prediction approach may lead to the loss of the obtained miRNA abundance. Multiple miRNAs were found to have the potential to regulate the same R2R3-MYB gene, and a single miRNA was also found to have the potential to regulate different R2R3-MYB genes in longan (Supplemental Fig. 1). The longan R2R3-MYBs regulated by the largest number of miRNAs were DlMYB92, DlMYB98, DlMYB115 and DlMYB103. Among them, both DlMYB98 and DlMYB92 were cleaved by different members of the miR159 and miR319 families. The homologs of DlMYB92 and DlMYB98 in Arabidopsis were identified as AtMYB33 and AtGAMYB, respectively. These genes share a close evolutionary relationship in both longan and Arabidopsis, as reported by Dubos et al. [4]. This finding suggests that the same miRNA family can regulate members of the R2R3-MYB family with similar gene structures. It was observed that the R2R3-MYB genes regulated by the same miRNA did not show overlap in covariance results. This suggests that R2R3-MYBs regulated by the same miRNA are not functionally redundant and may exhibit differences in their regulatory roles.

Conserved miRNA-MYB modules involved in root development are present in longan>

Out of the 124 identified R2R3-MYB family members in longan, 26 were found to be regulated by 34 miRNAs. Notably, 17 of these miRNAs were specific to longan and considered novel. The regulatory pathways governed by these miRNA-MYB interactions were directly associated with plant root development, highlighting their significance in this biological process. In Arabidopsis, miR159 has been shown to promote lateral root growth by targeting MYB transcription factors. Double mutants of miR159ab exhibit larger root tip meristem tissues, higher cell numbers, and longer primary roots compared to the wild type. miR159 negatively regulates the expression of AtMYB33, AtMYB65, and AtMYB101 genes, leading to increased cell division and primary root growth in the root tip meristematic tissue [15]. In our study, we predicted the involvement of the miR159 gene family (including miR159a, miR159b, and miR159c) in the regulation of DlMYB92 (homolog of AtMYB33) in longan.

Another regulator of longan R2R3-MYB is miR319. Previous research has demonstrated the existence of differentially expressed mRNAs and miRNAs during the development of radish (Raphanus sativus) tubers. Among them, miR319 candidate target genes, such as RSG11844.t1, RSG42419.t1, and RSG49768.t1, are involved in tuber formation and development [33]. In citrus, up-regulation of miR319 expression was observed in roots under 400 μmol/ L H3BO3 stress, resulting in the down-regulation of the target MYB gene and subsequent alteration of root morphological structure, including a decrease in root tip number [34]. We also discovered that DlMYB98 exhibited significantly higher expression in lateral roots compared to primary roots and was regulated by miR319 in longan. In summary, our findings suggest the presence of a conserved mechanism involving miRNA-DlMYB regulation in root development in longan. However, further in-depth investigation is required to elucidate the precise mechanisms by which miR159 and miR319 regulate the expression of DlMYB92 and DlMYB98, respectively, and subsequently influence longan root development.

Methods

Material

The materials used in this study were ‘Chu Liang’ longan, all from the Guangdong Longan Germplasm Resource Nursery. Longan seeds were sown in moist river sand in April 2023, and after about 30 days of growth, the resulting seedlings were used in the experiment. Primary and lateral root parts of the seedings were collected in triplicate. These samples were rapidly frozen in liquid nitrogen and stored in a -80℃ refrigerator for subsequent studies.

Identification and physicochemical property analysis of the R2R3-MYB gene of longan

The longan genome and gene annotation files were obtained through the SapBase database (http://www.sapindaceae.com/index.html). HMM files were downloaded from the Pfam database (http://pfam-legacy.xfam.org/). A sequence search of the longan genome was performed using Simple HMM Search in TBtools software, using PF00249 as a model [35]. Candidate genes were then submitted to the Batch CDD-search (https://www.ncbi.nlm.nih.gov/Structure/cdd/wrpsb.cgi) database for comparison to verify the R2R3-MYB conserved domains and remove sequences without MYB domains or with incomplete domains. The isoelectric point, molecular mass, hydrophilic mean coefficient, instability coefficient, and aliphatic coefficient of the R2R3-MYB protein were calculated through Expasy online (https://www.expasy.org/).

Phylogenetic characterization of longan R2R3-MYB family members

The ClustalW software compared the longan R2R3-MYB amino acid sequence with 124 Arabidopsis R2R3-MYB amino acid sequences. The Neighbor-Joining (NJ) method was used to construct a phylogenetic tree by MEGA software [36]. The members of the longan R2R3-MYB gene family were grouped with reference to the grouping of Arabidopsis in the phylogenetic tree [4].

Analysis of conserved motifs, conserved domains, and gene structures of longan R2R3-MYB

Conserved motif analysis of the longan R2R3-MYB protein sequence was performed using the MEME online tool (https://meme-suite.org/meme/doc/meme.html), and the number of conserved motifs was set at 10. Conserved domain analysis of the longan R2R3-MYB protein sequence was performed using the CDD-search tool on the NCBI website. The exon, intron, and untranslated (UTR) information corresponding to the R2R3-MYB gene were extracted from the longan genome using TBtools software for gene structure analysis. Finally, the conserved motifs, conserved domains, and gene structures of the genes were synchronized with the phylogenetic tree of R2R3-MYB in longan using TBtools software for display.

Chromosome distribution and colinearity analysis of R2R3-MYB gene in longan

Chromosomal location information of the R2R3-MYB genes was extracted from the longan genomic information using TBtools software and visualized for analysis. In addition, gene maps with duplicated fragments were selected and analyzed by the Multiple Colinearity Scanning Toolkit (MCSan X) for default parameters [37]. The synonymous relationship of the R2R3-MYB gene was determined by using Advanced Circos software in TBtools.

qRT-PCR analysis

Total RNA was extracted from the primary and lateral roots of ‘Chu Liang’ longan using RNAprep Pure Polysaccharide Polyphenol Plant Total RNA Extraction Kit (DP441). cDNA was synthesized using PC44-THERMOscript 1st Strand cDNA Synthesis Kit (PC4402). mRNA was used as template to design DlMYB33, DlMYB34, DlMYB35, DlMYB47, DlMYB48, DlMYB52, DlMYB59, DlMYB69, DlMYB70, DlMYB77, DlMYB83, DlMYB92 and DlMYB98 quantitative primers (Supplemental Table 1). qRT-PCR analysis was performed using the QuantStudio 3 Real-Time PCR System. The reaction system was performed according to the instructions of the PC60-2 x SYBR Green qPCR Mix (Low ROX) kit, using DlActin as the internal reference gene. Each sample run was repeated thrice, and the relative gene expression was calculated using the 2−ΔΔCT method [38].

Prediction of longan miRNA and miRNA sequencing analysis

The sRNA Target Prediction function in the online database SapBase was used to predict miRNAs capable of shearing the R2R3-MYB gene in longan. miRNA high-throughput sequencing was performed using ‘Chu Liang’ longan root as a sample, with three replicates set. miRNA library construction and high-throughput sequencing were co-performed with Biomarker Technologies Co and Beijing Biomarker Technologies Co The target gene analysis of miRNA was performed using BMKCloud (www.biocloud.net). The sequencing data are not yet available.

Subcellular Localization of DlMYB92 in tobacco leaf cells

The full-length coding sequence (CDS) of DlMYB92 was amplified by KOD FX using the cDNA of the longan root as a template. The target sequences were detected by gel electrophoresis and recovered, and the full-length CDS of DlMYB92 was ligated to the linearized pRI101-AN vector using NEB Builder HIFI DNA Assembly. The recombinant plasmid was then transferred into Trans1-T1 Phage Resistant Chemically Competent Cell, and PCR was used to confirm presence of the recombinant vector.

The recombinant vector was transformed into GV3101(pSoup-p19) Chemically Competent Cell, and PCR identified the positive clones. The identified positive clone was incubated in LB medium containing rifampicin (50 mg/L) and kanamycin (50 mg/L) with shaking (28 °C, 200 r/ min) until the OD600 value of the bacterial solution was 1. The bacterial solution was collected by centrifugation at 4000 g for 10 min. The bacteria were suspended in an MS liquid medium with 10 mM MgCl2, 15 μM acetosyringone, and 10 mM MES, and the OD600 was adjusted to 0.6. Wild-type tobacco (Nicotiana benthamiana) in good growth condition was selected for Agrobacterium infiltration. Each leaf was injected with 0.5 mL of bacterial solution, and a non-recombinant vector was injected as a control. This was performed using an OLYMPUS FV3000 laser confocal microscope.

Availability of data and materials

The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding author.

References

  1. Feng K, Hou XL, Xing GM, Liu JX, Duan AQ, Xu ZS, ..., Xiong AS. Advances in AP2/ERF superfamily transcription factors in plant. Crit Rev Biotechnol. 2020;40(6):750–776.

  2. Shoji T, Yuan L. ERF gene clusters: working together to regulate metabolism. Trends Plant Sci. 2021;26(1):23–32.

    Article  CAS  PubMed  Google Scholar 

  3. Wang X, Niu Y, Zheng Y. Multiple functions of myb transcription factors in abiotic stress responses. Int J Mol Sci. 2021;22(11):6125.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  4. Dubos C, Stracke R, Grotewold E, Weisshaar B, Martin C, Lepiniec L. MYB transcription factors in Arabidopsis. Trends Plant Sci. 2010;15(10):573–81.

    Article  CAS  PubMed  Google Scholar 

  5. Gangadhar BH, Venkidasamy B, Samynathan R, Saranya B, Chung IM, Thiruvengadam M. Overview of miRNA biogenesis and applications in plants. Biologia. 2021;76(8):2309–27.

    Article  CAS  Google Scholar 

  6. Choudhary A, Kumar A, Kaur H, Kaur N. MiRNA: the taskmaster of plant world. Biologia. 2021;76(5):1551–67.

    Article  CAS  Google Scholar 

  7. Curaba J, Singh MB, Bhalla PL. miRNAs in the crosstalk between phytohormone signalling pathways. J Exp Bot. 2014;65(6):1425–38.

    Article  CAS  PubMed  Google Scholar 

  8. Jin D, Wang Y, Zhao Y, Chen M. MicroRNAs and their cross-talks in plant development. J Genet Genomics. 2013;40(4):161–70.

    Article  CAS  PubMed  Google Scholar 

  9. Zhang B, Pan X, Cobb GP, Anderson TA. Plant microRNA: a small regulatory molecule with big impact. Dev Biol. 2006;289(1):3–16.

    Article  CAS  PubMed  Google Scholar 

  10. Jodder J. miRNA-mediated regulation of auxin signaling pathway during plant development and stress responses. J Biosci. 2020;45:91.

    Article  CAS  PubMed  Google Scholar 

  11. Li X, Xue C, Li J, Qiao X, Li L, Yu L, Huang Y, Wu J. Genome-wide identification, evolution and functional divergence of myb transcription factors in chinese white pear (Pyrus bretschneideri). Plant Cell Physiol. 2016;57(4):824–47.

    Article  CAS  PubMed  Google Scholar 

  12. Meng Y, Ma X, Chen D, Wu P, Chen M. MicroRNA-mediated signaling involved in plant root development. Biochem Biophys Res Commun. 2010;393(3):345–9.

    Article  CAS  PubMed  Google Scholar 

  13. Tang W, Tang AY. MicroRNAs associated with molecular mechanisms for plant root formation and growth. J Forestry Res. 2016;27(1):1–12.

    Article  CAS  Google Scholar 

  14. Xu H, Liu P, Wang C, Wu S, Dong C, Lin Q, Wu S. Transcriptional networks regulating suberin and lignin in endodermis link development and ABA response. Plant Physiol. 2022;190(2):1165–81.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Xue T, Liu ZH, Dai XH, Xiang FN. Primary root growth in Arabidopsis thaliana is inhibited by the miR159 mediated repression of MYB33, MYB65 and MYB101. Plant Sci. 2017;262:182–9.

    Article  CAS  PubMed  Google Scholar 

  16. Khandal H, Singh AP, Chattopadhyay D. The MicroRNA397b-LACCASE2 module regulates root lignification under water and phosphate deficiency. Plant Physiol. 2020;182(3):1387–403.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Xu M, Hu T, Zhao J, Park MY, Earley KW, Wu G, Yang L, Poethig RS. Developmental functions of miR156-regulated SQUAMOSA promoter binding protein-like (SPL) genes in Arabidopsis thaliana. PLoS Genet. 2016;12(8):e1006263.

    Article  PubMed  PubMed Central  Google Scholar 

  18. Shao YL, Zhou HZ, Wu YR, Zhang H, Lin J, Jiang XY, …, Mao CZ. OsSPL3, an SBP-domain protein, regulates crown root development in rice. Plant Cell. 2019;31(6):1257–75.

  19. Ding D, Zhang L, Wang H, Liu Z, Zhang Z, Zheng Y. Differential expression of miRNAs in response to salt stress in maize roots. Ann Bot. 2009;103(1):29–38.

    Article  CAS  PubMed  Google Scholar 

  20. Gao R, Austin RS, Amyot L, Hannoufa A. Comparative transcriptome investigation of global gene expression changes caused by miR156 overexpression in Medicago sativa. BMC Genomics. 2016;17(1):658.

    Article  PubMed  PubMed Central  Google Scholar 

  21. Li X, Shen F, Xu X, Zheng Q, Wang Y, Wu T, …, Zhang X. An HD-ZIP transcription factor, MxHB13, integrates auxin-regulated and juvenility-determined control of adventitious rooting in Malus xiaojinensis. Plant J. 2021;107(6):1663–80.

  22. Chen C, Chen H, Zhang Y, Thomas HR, Frank MH, He Y, Xia R. TBtools: an integrative toolkit developed for interactive analyses of big biological data. Mol Plant. 2020;13(8):1194–202.

    Article  CAS  PubMed  Google Scholar 

  23. Tang YY, He XM, Sun J, Li CB, Li L, Sheng JF, …, Ling DN. Polyphenols and alkaloids in byproducts of longan fruits (Dimocarpus Longan Lour.) And their bioactivities. Molecules. 2019;24(6):1186.

  24. Altendorf S. FAO Food Outlook 2018. Rome: Food and Agriculture Organization of the United Nations (FAO). p. 69–76. http://www.fao.org/.../ca0239en.pdf. ISSN: 0251-1959.

  25. Wang J, Li JG, Li ZY, Liu B, Zhang LL, Guo DL. …, Guo L.Genomic insights into longan evolution from a chromosome-level genome assembly and population genomics of longan accessions. Hortic Res. 2022;9:uhac021.

  26. Lin Y, Min J, Lai R, Wu Z, Chen Y, Yu L, …, Lai Z. Genome-wide sequencing of longan (Dimocarpus longan Lour.) Provides insights into molecular basis of its polyphenol-rich characteristics. GigaScience. 2017;6(5):1–14.

  27. Qing J, Dawei W, Jun Z, Yulan X, Bingqi S, Fan Z. Genome-wide characterization and expression analyses of the MYB superfamily genes during developmental stages in chinese jujube. PeerJ. 2019;7:e6353.

    Article  PubMed  PubMed Central  Google Scholar 

  28. Li C, Zhang BH. MicroRNAs in control of plant development. J Cell Physiol. 2016;231(2):303–13.

    Article  CAS  PubMed  Google Scholar 

  29. Du H, Yang SS, Liang Z, Feng BR, Liu L, Huang YB, Tang YX. Genome-wide analysis of the MYB transcription factor superfamily in soybean. BMC Plant Biol. 2012;12: 106.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Yang Y, Yuan Z, Ning C, Zhao B, Wang R, Zheng X, He L. The pea R2R3-MYB gene family and its role in anthocyanin biosynthesis in flowers. Front Genet. 2022;13: 936051.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Sun WJ, Ma ZT, Chen H, Liu MY. MYB gene family in potato (Solanum tuberosum L.): genome-wide identification of hormone-responsive reveals their potential functions in growth and development. Int J Mol Sci. 2019;20(19):4847.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Arce-Rodriguez ML, Martinez O, Ochoa-Alejo N. Genome-wide identification and analysis of the MYB transcription factor gene family in chili pepper (Capsicum spp.). Int J Mol Sci. 2021;22(5):2229.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. Liu C, Liu X, Xu W, Fu W, Wang F, Gao J, ..., Wang S. Identification of miRNAs and their targets in regulating tuberous root development in radish using small RNA and degradome analyses. 3 Biotech. 2018;8(7):311.

  34. Huang JH, Lin XJ, Zhang LY, Wang XD, Fan GC, Chen LS. MicroRNA sequencing revealed citrus adaptation to long-term boron toxicity through modulation of root development by miR319 and miR171. Int J Mol Sci. 2019;20(6):1422.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Chen YH, Sun JZ, Lin HT, Lin MS, Lin YF, Wang H, Hung YC. Salicylic acid reduces the incidence of Phomopsis longanae Chi infection in harvested longan fruit by affecting the energy status and respiratory metabolism. Postharvest Biol Tec. 2020;160:111035.

    Article  CAS  Google Scholar 

  36. Kumar S, Stecher G, Li M, Knyaz C, Tamura K. MEGA X: molecular evolutionary genetics analysis across computing platforms. Mol Biol Evol. 2018;35(6):1547–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Wang Y, Tang H, Debarry JD, Tan X, Li J, Wang X, Paterson AH. MCScanX: a toolkit for detection and evolutionary analysis of gene synteny and collinearity. Nucleic Acids Res. 2012;40(7):e49.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2 (-Delta Delta C(T)) method. Methods. 2001;25(4):402–8.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

Not applicable.

Research involving plants

The experimental research and field studies involving the ‘Chu Liang’ longan used in this study, including the collection of plant material, all methods comply with relevant guidelines, regulations, and legislations.

Funding

This study was supported by the Basic and Applied Basic Research Topics (SL2022A04J00407) and Key-Area Research and Development Program of Guangdong Province (2022B0202070002).

Author information

Authors and Affiliations

Authors

Contributions

XL and YZ designed the study. ST and JW performed the experiments and data collection. SH, and DH provided software and methods in data processing and analysis. JL and DG analyzed and discussed the data. All authors revised and discussed subsequent versions. All authors contributed to the article and approved the submitted version.

Corresponding authors

Correspondence to Dongliang Guo or Yan Zhou.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Supplemental Figure 1.

Regulatory network of miRNAs and Longan MYB family members.

Additional file 2: Supplemental Table 1.

Primers used in this study.

Additional file 3.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Lv, X., Tian, S., Huang, S. et al. Genome-wide identification of the longan R2R3-MYB gene family and its role in primary and lateral root. BMC Plant Biol 23, 448 (2023). https://doi.org/10.1186/s12870-023-04464-9

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12870-023-04464-9

Keywords