Next Article in Journal
Assessment of the Impact of Trace Essential Metals on Cancer Development
Next Article in Special Issue
Analysis of the Correlation between Persimmon Fruit-Sugar Components and Taste Traits from Germplasm Evaluation
Previous Article in Journal
Redefining the Role of Ornithine Aspartate and Vitamin E in Metabolic-Dysfunction-Associated Steatotic Liver Disease through Its Biochemical Properties
Previous Article in Special Issue
Exogenous Melatonin Application Accelerated the Healing Process of Oriental Melon Grafted onto Squash by Promoting Lignin Accumulation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Comparative Transcriptome Analysis Reveals Key Functions of MiMYB Gene Family in Macadamia Nut Pericarp Formation

1
Guangxi South Subtropical Agricultural Research Institute, Longzhou 532415, China
2
Rubber Research Institute, Chinese Academy of Tropical Agricultural Sciences, Haikou 571101, China
3
Guangxi Academy of Agricultural Sciences, Nanning 530007, China
*
Authors to whom correspondence should be addressed.
Int. J. Mol. Sci. 2024, 25(13), 6840; https://doi.org/10.3390/ijms25136840
Submission received: 26 April 2024 / Revised: 15 June 2024 / Accepted: 19 June 2024 / Published: 21 June 2024
(This article belongs to the Special Issue Advances in Research on Fruit Crop Breeding and Genetics: 3rd Edition)

Abstract

:
Macadamia nuts are one of the most important economic food items in the world. Pericarp thickness and flavonoid composition are the key quality traits of Macadamia nuts, but the underlying mechanism of pericarp formation is still unknown. In this study, three varieties with significantly different pericarp thicknesses, namely, A38, Guire No.1, and HAES 900, at the same stage of maturity, were used for transcriptome analysis, and the results showed that there were significant differences in their gene expression profile. A total of 3837 new genes were discovered, of which 1532 were functionally annotated. The GO, COG, and KEGG analysis showed that the main categories in which there were significant differences were flavonoid biosynthesis, phenylpropanoid biosynthesis, and the cutin, suberine, and wax biosynthesis pathways. Furthermore, 63 MiMYB transcription factors were identified, and 56 R2R3-MYB transcription factors were clustered into different subgroups compared with those in Arabidopsis R2R3-MYB. Among them, the S4, S6, and S7 subgroups were involved in flavonoid biosynthesis and pericarp formation. A total of 14 MiMYBs’ gene expression were verified by RT-qPCR analysis. These results provide fundamental knowledge of the pericarp formation regulatory mechanism in macadamia nuts.

Graphical Abstract

1. Introduction

Macadamia (Macadamia integrifolia) is a perennial evergreen fruit tree of the family Macadamia, which is native to the subtropical rainforest areas of southeastern Queensland and the northeastern coast of New South Wales, Australia [1]. Known as the “Queen of Dried Fruits” and the “King of Nuts”, macadamias are the world’s most popular edible nuts and have become one of the world’s emerging nut industries. China began to commercially introduce macadamias in the 1970s, and they are now mainly distributed in Yunnan, Guangxi, and Guangdong provinces and regions. At present, China is the country with the fastest growth and the largest planting area of macadamias in the world, with more than 330,000 hectares, ranking first in the world, and the output of unhulled fruit is 110,000 tons, ranking fourth in the world. Macadamias comprise three parts: the pericarp, the shell, and the kernel. The pericarp contains a variety of nutrients, such as crude fat, tannins, crude protein, and phenolic acids [2], which have various properties, such as disease resistance, antioxidant effects, and a protective effect on the shell and kernel. However, since the pericarp is prone to oxidation and fermentation when stored in a stack, it is necessary to peel it on the day of harvesting. Machinery suitable for macadamia fruit processing has been developed. Significant progress has been made in the reuse of macadamia pericarp, as large quantities of macadamia pericarp have been separated [3]. Our previous studies found that the nutrient composition of different macadamia germplasm pericarp was significantly different, and the coefficient of variation of inclusions was between 15.75% and 48.33%. The tannin contents range from 1.06 to 2.16% in macadamia nut pericarps. These pericarps can be fermented into organic fertilizer to increase soil nutrients and organic matter. Recent studies have shown that macadamia pericarp color is associated with changes in kernel quality and can be used as a predictor of harvest time. A novel ultrasound-assisted extraction method can effectively extract anti-aging components from the pericarp and improve its economic value. With advances in macadamia separation and processing technology, macadamia pericarp is widely used in high value-added industries, such as cosmetics and organic fertilizers [3]. The main problem facing China’s macadamia industry is the urgent need for new varieties that are easy to process and rich in active ingredients. It is an effective way to solve the problem of industrial demand by screening and creating germplasm with significant differences in macadamia kernel and pericarp through grafting, hybridization, and transgenic technology, then cultivating high-quality varieties.
Recently, research on plant pericarp or peel development and nutrient composition has made rapid progress. For instance, to measure pericarp thickness, a new non-destructive method using extended-focus optical coherence microscopy was developed [4]. In sweetcorn, it was found that auxin (AUX), gibberellin (GA), and brassinosteroid (BR) signal transduction may indirectly mediate PCD to regulate pericarp thickness in the corn variety with a thin pericarp [5]. By using QTL mapping and transcriptome analysis, it was found that AUX/IAA transcription factor and ZIM transcription factor regulate pericarp thickness [6]. In pomegranate, naringenin, pelargonidin, kaempferol, and CHI (Pgr25966.1), F3′5′H (Pgr26644.1), and CHS (Pgr005566.1) genes were related to pericarp’s response to sunburn, and MYB93 and MYB111 may be involved in phenylpropanoid and flavonoid biosynthesis by regulating these genes [7]. Combining transcriptome and metabolome analyses, the transcription factors NAC and MYB were found to possibly be involved in major transcriptional regulatory mechanisms [8]. In rice, it was found that R2R3 MYB TF OsKala3 is a common key player for black rice pericarp [9]. In pepper, CaR2R3-MYB genes are involved in regulating the synthesis of capsaicin and dihydrocapsaicin [10,11]. In kiwifruit, AcMYB10 with AcbHLH42 may strongly activate anthocyanin biosynthesis by activating the transcription of AcLDOX and AcF3GT [12]. AcMYB123 and AcbHLH42 are the components involved in the spatiotemporal regulation of anthocyanin biosynthesis, specifically in the inner pericarp of kiwifruit [13]. Cytokinin is involved in the regulation of tomato pericarp thickness and fruit size [14]. The highest concentration of lignin was detected in the pericarp cell wall of the GP12 inbred line, with extremely high popping expansion [15]. Based on these results, we proposed the scientific hypothesis that the regulatory mechanism of pericarp thickness and components is related to specific transcription factors. Hence, screening new transcription factors related to pericarp thickness, creating new germplasm with high-quality pericarp traits and cultivating new varieties of macadamia nuts are effective ways to address the needs of the industry both in China and globally.

2. Results

DEGs in the Three Macadamia Varieties with Different Pericarp Thicknesses

A transcriptome analysis of nine samples from the three macadamia varieties with different pericarp thicknesses (A38, Guire No.1, and HAES900), and three repeats, was completed. A total of 61.30 Gb of clean data was obtained, and the clean data of each sample reached 6.07 Gb, while the Q30 base percentage was 94.00% or above. The clean reads of each sample were compared with the specified reference genome, and the comparison efficiency ranged from 87.52% to 90.94% (Supplementary Table S1). Based on the comparative results, an alternative splicing prediction analysis, a gene structure optimization analysis, and new gene discovery were carried out, and 3837 new genes were discovered, of which 1532 were functionally annotated. As can be seen from Figure 1, the number of differentially expressed genes (DEGs) in the transcriptome analysis of pericarp between the three varieties reached a significant level. The total number of differential genes between A38 and Guire No.1 was 4349, of which 2590 were upregulated and 1759 downregulated. The total number of DEGs between A38 and HAES900 was 3615, of which 2059 were upregulated and 1556 were downregulated. The total number of DEGs between Guire No.1 and HAES900 was 4094, of which 1931 were upregulated and 2163 downregulated (Figure 1).
As shown in Figure 2, the upregulated differential expression genes between A38 and Guire No.1 were mainly in flavonoid biosynthesis, phenylpropanoid biosynthesis, and phenylalanine, tyrosine, and tryptophan biosynthesis, and the gene ratio was 4.53%, 5.58%, and 1.79%, respectively. The downregulated differential expression genes between A38 and Guire No.1 were mainly in photosynthesis and cutin, suberine, and wax biosynthesis, and the gene ratio was 3.23% and 2.04%, respectively. The upregulated differential expression genes between A38 and HAES900 were mainly in phenylpropanoid biosynthesis, phenylalanine, tyrosine, and tryptophan biosynthesis, and the flavonoid biosynthesis pathways, and the gene ratio was 5.33%; 1.47%; and 2.27%, respectively. The downregulated differential expression genes between A38 and HAES900 were mainly in fatty acid degradation, flavonoid biosynthesis, and cutin, suberine, and wax biosynthesis pathways, and the gene ratio was 1.98%, 2.58%, and 1.79%, respectively. The upregulated differential expression genes between Guire No.1 and HAES900 were mainly in photosynthesis—antenna proteins, photosynthesis, and carbon fixation in photosynthetic organisms, and the gene ratio was 2.40%, 4.38%, and 2.26%, respectively. The downregulated differential expression genes between Guire No.1 and HAES900 were mainly in phenylpropanoid biosynthesis and phenylalanine, tyrosine, and tryptophan biosynthesis, and the gene ratio was 4.41% and 1.13%, respectively (Supplementary Table S2).
The COG (Cluster of Orthologous Groups of proteins) database enables the orthologous classification of gene products. As shown in Figure 3 and Supplementary Table S3, the most significant differences in the pathway between A38, Guire No.1, and HAES900 are within three main branches: biological process, molecular function, and cellular component. For molecular function, transcription factor activity, protein binding, and nucleic acid binding transcription factor activity have a high DEG ratio. To further explore the key transcription factors in pericarp thickness and nutrition composition, the differentially expressed transcription factors were analyzed (Supplementary Table S4), and the MYB, bHLH, and NAC transcription factor families were the main differentially expressed gene families. Furthermore, all upregulated and downregulated genes were comprehensively analyzed to find potential target genes of MYB, bHLH, and NAC transcription factor (Supplementary Tables S5 and S6).
As shown in Figure 4, for flavonoid biosynthesis-related differentially expressed genes, most of which were highly expressed in Guire No.1, the FPKM of the most highly expressed PAL, C4H, F3H, and DFR was 106.62, 444.66, 231.89, and 134.45, respectively (Supplementary Table S7). For phenylpropanoid biosynthesis-related differentially expressed enzyme genes, most of which were highly expressed in Guire No.1, the FPKM of LOC122092787, LOC122089414, LOC122069164, LOC122078387, LOC122067873, LOC122088272, LOC122073309, and LOC122063458 was 291.94, 245.63, 239.09, 204.56, 181.53, 172.22, 170.63, and 140.55, significantly higher than that in A38 and HAES900 (Supplementary Table S8). Most of the cutin, suberine, and wax biosynthesis-related differentially expressed genes were highly expressed in Guire No.1. The FPKM of LOC122066392, LOC122081968, and LOC122071897 was 216.37, 24.35, and 11.67, respectively (Supplementary Table S9).
As shown in Figure 5, a total of 63 MiMYB TFs were identified in this study. According to their motif and protein structures, four 1R-MYB, fifty-six 2R-MYB, one 3R-MYB, and two 4R-MYB transcription factors were classified (Supplementary Table S10). Three specific motifs were found and used as the main indices for classification.
Furthermore, a phylogenetic analysis of 56 R2R3-MiMYBs compared with those of Arabidopsis was conducted; as shown in Figure 6, the 56 R2R3-MiMYBs can be divided into nine subgroups, namely, S4, S5, S6, S7, S13, S14, S15, S22, and S25 (Supplementary Table S10). Of these, S4, S6, and S7 participate in flavonoid biosynthesis. In the S4 subgroup, LOC122066430 and LOC122083708 were the homologous genes of AtMYB4 and AtMYB6 [16]. In the S6 subgroup, LOC122085905 and LOC122080976 were the homologous genes of AtMYB114 and AtMYB90 [17,18]. In the S7 subgroup, LOC122089129 and LOC122089330 were the homologous genes of AtMYB11 and AtMYB12 [19,20].
Furthermore, 14 MiMYBs’ gene expression level from S4, S6, S7, and S22 subgroups were verified by RT-qPCR method (Supplement Table S12). As shown in Figure 7, the fluorescence quantitative results were consistent with those of transcriptome analysis.

3. Discussion

With the continuous innovation and expansion of genetics, breeding, and food science and technology, studies on the content and antioxidant activity of different solvent extracts of macadamia pericarp have been carried out at home and abroad, and obvious differences have been identified in the contents of total phenols, total flavonoids, tannins, and total antioxidant capacity in macadamia pericarps [21,22]. A genomic analysis of Guire No.1 showed that the macadamia genome has 14 chromosomes, encoding 37,728 genes [23].
In this study, the main innovative finding is that the S4, S6, and S7 MiMYB subgroups participate in flavonoid biosynthesis, which is related to pericarp thickness and nutrition composition. The MYB protein family is a large and functionally diverse family that is expressed in all eukaryotes [24]. The MYB protein has a highly conserved DNA-binding domain, which typically consists of four incomplete amino acid repeats (R), each of which has approximately 52 amino acids, and each repeat forms three α helixes. The second and third helixes of each repeat form a helix–turn–helix (HTH) structure with three regularly spaced tryptophan (or hydrophobic) residues, forming a hydrophobic core in a three-dimensional HTH structure [25]. MYB proteins can be classified into different classes based on the number of adjacent domain repeats (1, 2, 3, or 4): R1/R2-MYB, R2R3-MYB, R1R2R3-MYB, and 4R-MYB [26,27]. MYB plays an important role in regulating plant cell growth and development, abiotic stress and biotic stress, and metabolic reactions. The regulatory role of MYB family members in plant fruit traits mainly includes the direct regulation of pericarp thickness and the regulation of the active ingredients in fruit [28,29,30]. Based on a genome study of the Guire No.1 variety [23], a transcriptome analysis was performed of macadamia nut varieties A38, Guire No.1, and HAES900, with significant differences in pericarp thickness. A total of 63 MYB family members were found in the transcriptome (Figure 5 and Figure 6). We further found that the S4, S6, and S7 MiMYB subgroups participate in flavonoid biosynthesis. LOC122066430 and LOC122083708 are homologous genes of AtMYB4 and AtMYB6, etc. These are new findings which can be supported by previous studies. In a joint study of transcriptome and miRNA in maize pericarp thickness, it was found that miRNA164, miRNA167, and miRNA156 regulated miRNA–mRNA pairs that were involved in regulating pericarp thickness. miRNA164 regulates MYB transcription factors [5]. A transcriptome analysis revealed a maximum of 64 MYB differential genes, suggesting that MYB may directly regulate peel thickness [5]. On the other hand, MYB is a key regulatory gene for active ingredients in the peel. For example, MYB93 and MYB11 prevent damage to the peel from sunburn by modulating the biosynthesis of phenylpropionic acid and flavonoids [7]. AcMYB123 and AcbHLH42 in kiwifruit regulated the anthocyanin synthesis genes AcF3GT1 and AcANS during fruit and peel color change [13]. In the study of fruit acidity regulation, it was found that LcMYB5 in litchi regulates pH by regulating the acidification gene LcPH1 [31]. The spraying of gibberellin can brown the lychee peel, which is associated with LcMYB1 [32]. Phenols alter the thickness of corn peels and are associated with the c1/pl1/p1 Myb gene family of fumaric acid toxins [33]. A combined nuclear magnetic resonance (NMR) and transcriptome analysis showed that MYB was associated with proximal axial development [34]. In the study of the mechanism of peel protection, it was proved that rice OsMYB102 and tomato SlMYB102 transgenic Arabidopsis thaliana can improve antioxidant and senescence properties, and the mechanism is the MYB transcription factors directly regulating the expression of downstream SOS1 and CYP707A3 genes [35,36]. The expression patterns of MiMYBs and their target genes in macadamia nuts with different pericarp thickness, light intensity, gibberellin, and other hormones also required analysis. This study reveals the transcriptional regulatory mechanism of macadamia pericarp flavonoid synthesis and provides a solid foundation for the screening of high-quality and high-yield macadamia nut germplasm resources.

4. Materials and Methods

4.1. Plant Materials

Three macadamia nut varieties, namely, A38, Guire No.1, and HAES900, with significant differences in pericarp size, were selected as study materials to determine their pericarp thickness, weight, and total flavonoid content (Figure 8). The fruit samples of the above cultivars were obtained from the macadamia germplasm resource nursery (E106°79′85″, N22°34′13″) of the Institute of Subtropical Agricultural Sciences in Guangxi, China. Usually, the full-blossom stage of the three macadamia nut varieties are from 24 to 27 March 2023. The fruit samples were collected at the same time, i.e., 120 days after pollination (dpp). The samples were collected in liquid nitrogen and stored in a −80 °C ultra-low-temperature freezer (Thermo Scientific, Waltham, MA, USA).

4.2. Methods

4.2.1. Determination of Total Flavonoids

The total flavonoid content was determined using the method described by Benamar et al. (2010). An aliquot of 500 μL of pericarp extracts in methanol was mixed with 1500 μL of distilled water, followed by the addition of 150 μL of 5% NaNO2 (in water), and the mixture was allowed to react for 5 min. Following this, 150 μL of 10% AlCl3 (in water) was added and the mixture stood for a further 6 min. Finally, the reaction mixture was treated with 500 μL of 1 M NaOH (in water) and the absorbance at 510 nm was obtained against a blank prepared similarly, by replacing the extract with methanol using a spectrophotometer (YoukeL6, Shanghai, China). Total flavonoid content was calculated from a calibration curve, using catechin (5–25 μg mL−1) as the standard, and expressed as mg of catechin equivalents g−1 of lyophilized pericarp extracts. The experiment was performed in triplicate [37].

4.2.2. RNA Extraction and RNA-Seq Analysis

The total RNA in the macadamia nut pericarp was extracted using an RNAprep Pure Plant Plus Kit (Tiangen Biotech, Beijing, China). The pericarp samples were subjected to RNA-seq with Biomarker technologies. Three biological replicates were carried out. RNA was extracted from each sample, and the concentration and integrity of each RNA sample were examined using NanoDrop (Thermo Scientific, Waltham, MA, USA), Qubit 2.0 (Invitrogen, Carlsbad, CA, USA), Agilent 2100 (Agilent, Palo Alto, CA, USA), etc. Only RNA with good quality could move on to the following procedures. Qualified RNA samples were processed for library construction. In order to ensure the quality of the library, Qubit 2.0 and Agilent 2100 were used to examine the concentration of cDNA and insert size. qPCR was performed to obtain a more accurate library concentration. A library with a concentration larger than 2 nM was acceptable. The qualified library was pooled based on pre-designed target data volume and then sequenced on an Illumina sequencing platform. Clean data with high quality were obtained by filtering raw data, which removed adapter sequences and reads with low quality. These clean data were further mapped to pre-defined reference genome-generating mapped data. Assessments of insert size and sequencing randomness were performed on the mapped data as library quality control. A basic analysis of the mapped data included gene expression quantification, alternative splicing analysis, novel gene prediction, and gene structure optimization. DEGs with a fold change ≥ 1.5 and a p-value < 0.05 were selected to analyze GO and KEGG enrichment. In order to understand the GO entries that were significantly enriched compared with the whole genomic background, clusterProfiler was used to conduct an enrichment analysis of biological processes, molecular functions, and cell components by using hypergeometric testing methods for the differential gene sets of each group. The term obtained from the enrichment results was visualized with a histogram. The function of different genes in this group could be predicted based on the information of GO functional enrichment (Supplementary Table S11) [38]. The filtered sequences were aligned to the reference genome and the transcriptome sequence data were assembled using StringTie software v2.2.0 [39]. The transcription factors were identified based on gene annotation in the database and the PlantTFDB platform v5.0 https://planttfdb.gao-lab.org/tf.php?sp=Mtr&did=Medtr1g086510.1 (20 March 2024) [40]. Using Muscle v5 [41] for multiple sequence alignment, IQ-TREE software v2.3.4 [42] was employed to construct a phylogenetic tree based on the maximum likelihood method. FigTree v1.4.4 was used to adjust the constructed phylogenetic tree.

4.2.3. cDNA Synthesis and RT-qPCR Verification

First, 1 µg of the sample was reverse-transcribed to cDNA with RevertAid Reverse Transcriptase (Thermo Scientific, Waltham, MA, USA) as follows: (1) Oligo(dT)18 primer at 1 µL and total RNA at 1 ng were added with water to 12 µL, and the reaction was run at 65 °C for 5 min. (2) The following components were added in the indicated order: 5× Reaction Buffer at 4 µL, RiboLock RNase Inhibitor (20 U/µL) at 1 µL, 10 mM dNTP Mix at 2 µL, and RevertAid M-MuLV RT (200 U/µL) at 1 µL for a total volume of 20 µL; the reaction was run at 42 °C for 60 min, and then at 70 °C for 5 min. RT-qPCR was conducted using CFX96 (Bio-Rad, Hercules, CA, USA) according to the following steps: 8.2 µL of ddH2O, 0.4 µL of forward primer, 0.4 µL of reverse primer, 10.4 µL of cDNA, and 2 × 10 µL of ChamQ Universal SYBR qPCR Master Mix (Vazyme, Nanjing, China) were added into one reaction solution with real-time PCR. The following reaction steps were used: 95 °C for 30 s to denature; 5 °C for 10 s and 60 °C for 30 s, up to 40 cycles; 95 °C for 15 s, 60 °C for 60 s, and 95 °C for 15 s. MiActin was used as a reference gene (HQ260674.1) (Supplement Table S12). Three biological replicates were performed for each sample. Expression levels were calculated using 2−∆∆Ct. Three biological repeats and three experimental replicates were carried out for each sample [38].

Supplementary Materials

The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/ijms25136840/s1.

Author Contributions

Conceptualization W.W. and L.W.; methodology, Q.T.; software, Q.T.; formal analysis, Q.T., X.H., Z.P., X.Y., Y.W. and C.Z.; investigation Q.T., X.H., Z.P., X.Y., Y.W. and C.Z.; resources W.W.; writing—original draft preparation, Q.T.; writing—review and editing, W.W. and L.W.; funding acquisition, W.W. and Q.T. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Science and Technology Development Fund Project of Guangxi Academy of Agricultural Sciences, grant number Guinongke2021YT156, 2023YM42 and Chongzuo Science and Technology Plan Project, grant number Chongke2022QN1220.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All raw data can be found on the website “http://www.ncbi.nlm.nih.gov/bioproject/1091847 (26 March 2024)”, the BioProject ID: PRJNA1091847.

Conflicts of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as potential conflicts of interest.

References

  1. Nock, C.J.; Baten, A.; Mauleon, R.; Langdon, K.S.; Topp, B.; Hardner, C.; Furtado, A.; Henry, R.J.; King, G.J. Chromosome-Scale Assembly and Annotation of the Macadamia Genome (Macadamia integrifolia HAES 741). G3 Genes Genomes Genet. 2020, 10, 3497–3504. [Google Scholar] [CrossRef] [PubMed]
  2. Abubaker, M.; Hawary, S.S.E.; Mahrous, E.A.; El-Kader, E.M.A. Study of Nutritional Contents of Macadamia integrifolia Maiden and Betche Leaves, Kernel and Pericarp Cultivated in Egypt. Int. J. Pharmacogn. Phytochem. Res. 2018, 9, 1442–1445. [Google Scholar] [CrossRef]
  3. Somwongin, S.; Sirilun, S.; Chantawannakul, P.; Anuchapreeda, S.; Yawootti, A.; Chaiyana, W. Ultrasound-assisted green extraction methods: An approach for cosmeceutical compounds isolation from Macadamia integrifolia pericarp. Ultrason. Sonochem. 2023, 92, 106266. [Google Scholar] [CrossRef] [PubMed]
  4. Sen, D.; Fernández, A.; Crozier, D.; Henrich, B.; Sokolov, A.V.; Scully, M.O.; Rooney, W.L.; Verhoef, A.J. Non-Destructive Direct Pericarp Thickness Measurement of Sorghum Kernels Using Extended-Focus Optical Coherence Microscopy. Sensors 2023, 23, 707. [Google Scholar] [CrossRef] [PubMed]
  5. Xiong, C.; Pei, H.; Zhang, Y.; Ren, W.; Ma, Z.; Tang, Y.; Huang, J. Integrative analysis of transcriptome and miRNAome reveals molecular mechanisms regulating pericarp thickness in sweet corn during kernel development. Front. Plant Sci. 2022, 13, 945379. [Google Scholar] [CrossRef] [PubMed]
  6. Wu, X.; Wang, B.; Xie, F.; Zhang, L.; Gong, J.; Zhu, W.; Li, X.; Feng, F.; Huang, J. QTL mapping and transcriptome analysis identify candidate genes regulating pericarp thickness in sweet corn. BMC Plant Biol. 2020, 20, 117. [Google Scholar] [CrossRef] [PubMed]
  7. Liu, C.; Su, Y.; Li, J.; Jia, B.; Cao, Z.; Qin, G. Physiological adjustment of pomegranate pericarp responding to sunburn and its underlying molecular mechanisms. BMC Plant Biol. 2022, 22, 169. [Google Scholar] [CrossRef] [PubMed]
  8. Li, Y.; Liao, B.; Wang, Y.; Luo, H.; Wang, S.; Li, C.; Song, W.; Zhang, K.; Yang, B.; Lu, S.; et al. Transcriptome and metabolome analyses provide insights into the relevance of pericarp thickness variations in Camellia drupifera and Camellia oleifera. Front. Plant Sci. 2022, 13, 1016475. [Google Scholar] [CrossRef] [PubMed]
  9. Kim, D.-H.; Yang, J.; Ha, S.-H.; Kim, J.-K.; Lee, J.-Y.; Lim, S.-H. An OsKala3, R2R3 MYB TF, Is a Common Key Player for Black Rice Pericarp as Main Partner of an OsKala4, bHLH TF. Front. Plant Sci. 2021, 12, 765049. [Google Scholar] [CrossRef]
  10. Liu, R.; Song, J.; Liu, S.; Chen, C.; Zhang, S.; Wang, J.; Xiao, Y.; Cao, B.; Lei, J.; Zhu, Z. Genome-wide identification of the Capsicum bHLH transcription factor family: Discovery of a candidate regulator involved in the regulation of species-specific bioactive metabolites. BMC Plant Biol. 2021, 21, 262. [Google Scholar] [CrossRef]
  11. Wang, J.; Liu, Y.; Tang, B.; Dai, X.; Xie, L.; Liu, F.; Zou, X. Genome-Wide Identification and Capsaicinoid Biosynthesis-Related Expression Analysis of the R2R3-MYB Gene Family in Capsicum annuum L. Front. Genet. 2020, 11, 598183. [Google Scholar] [CrossRef] [PubMed]
  12. Yu, M.; Man, Y.; Wang, Y. Light- and Temperature-Induced Expression of an R2R3-MYB Gene Regulates Anthocyanin Biosynthesis in Red-Fleshed Kiwifruit. Int. J. Mol. Sci. 2019, 20, 5228. [Google Scholar] [CrossRef]
  13. Wang, L.; Tang, W.; Hu, Y.; Zhang, Y.; Sun, J.; Guo, X.; Lu, H.; Yang, Y.; Fang, C.; Niu, X.; et al. A MYB/bHLH complex regulates tissue-specific anthocyanin biosynthesis in the inner pericarp of red-centered kiwifruit Actinidia chinensis cv. Hongyang. Plant J. 2019, 99, 359–378. [Google Scholar] [CrossRef]
  14. Gan, L.; Song, M.; Wang, X.; Yang, N.; Li, H.; Liu, X.; Li, Y. Cytokinins are involved in regulation of tomato pericarp thickness and fruit size. Hortic. Res. 2022, 9, uhab041. [Google Scholar] [CrossRef] [PubMed]
  15. Damasceno Junior, C.V.; Godoy, S.; Gonela, A.; Scapim, C.A.; Grandis, A.; Dos Santos, W.D.; Mangolin, C.A.; Buckeridge, M.S.; Machado, M.d.F.P.S. Biochemical composition of the pericarp cell wall of popcorn inbred lines with different popping expansion. Curr. Res. Food Sci. 2021, 5, 102–106. [Google Scholar] [CrossRef]
  16. Zhang, L.; Duan, Z.; Ma, S.; Sun, S.; Sun, M.; Xiao, Y.; Ni, N.; Irfan, M.; Chen, L.; Sun, Y. SlMYB7, an AtMYB4-Like R2R3-MYB Transcription Factor, Inhibits Anthocyanin Accumulation in Solanum lycopersicum Fruits. J. Agric. Food Chem. 2023, 71, 18758–18768. [Google Scholar] [CrossRef]
  17. Velten, J.; Cakir, C.; Youn, E.; Chen, J.; Cazzonelli, C.I. Transgene silencing and transgene-derived siRNA production in tobacco plants homozygous for an introduced AtMYB90 construct. PLoS ONE 2012, 7, e30141. [Google Scholar] [CrossRef] [PubMed]
  18. Velten, J.; Cakir, C.; Cazzonelli, C.I. A spontaneous dominant-negative mutation within a 35S::AtMYB90 transgene inhibits flower pigment production in tobacco. PLoS ONE 2010, 5, e9917. [Google Scholar] [CrossRef]
  19. Pandey, A.; Misra, P.; Trivedi, P.K. Constitutive expression of Arabidopsis MYB transcription factor, AtMYB11, in tobacco modulates flavonoid biosynthesis in favor of flavonol accumulation. Plant Cell Rep. 2015, 34, 1515–1528. [Google Scholar] [CrossRef]
  20. Petroni, K.; Falasca, G.; Calvenzani, V.; Allegra, D.; Stolfi, C.; Fabrizi, L.; Altamura, M.M.; Tonelli, C. The AtMYB11 gene from Arabidopsis is expressed in meristematic cells and modulates growth in planta and organogenesis in vitro. J. Exp. Bot. 2008, 59, 1201–1213. [Google Scholar] [CrossRef]
  21. Woźniak, M.; Waśkiewicz, A.; Ratajczak, I. The Content of Phenolic Compounds and Mineral Elements in Edible Nuts. Molecules 2022, 27, 4326. [Google Scholar] [CrossRef] [PubMed]
  22. Winkel-Shirley, B. Flavonoid biosynthesis. A colorful model for genetics, biochemistry, cell biology, and biotechnology. Plant Physiol. 2001, 126, 485–493. [Google Scholar] [CrossRef] [PubMed]
  23. Xia, C.; Jiang, S.; Tan, Q.; Wang, W.; Zhao, L.; Zhang, C.; Bao, Y.; Liu, Q.; Xiao, J.; Deng, K.; et al. Chromosomal-level genome of macadamia (Macadamia integrifolia). Trop. Plants 2022, 1, 3. [Google Scholar] [CrossRef]
  24. Dubos, C.; Stracke, R.; Grotewold, E.; Weisshaar, B.; Martin, C.; Lepiniec, L. MYB transcription factors in Arabidopsis. Trends Plant Sci. 2010, 15, 573–581. [Google Scholar] [CrossRef] [PubMed]
  25. Ogata, K.; Kanei-Ishii, C.; Sasaki, M.; Hatanaka, H.; Nagadoi, A.; Enari, M.; Nakamura, H.; Nishimura, Y.; Ishii, S.; Sarai, A. The cavity in the hydrophobic core of Myb DNA-binding domain is reserved for DNA recognition and trans-activation. Nat. Struct. Biol. 1996, 3, 178–187. [Google Scholar] [CrossRef] [PubMed]
  26. Rowan, D.D.; Cao, M.; Lin-Wang, K.; Cooney, J.M.; Jensen, D.J.; Austin, P.T.; Hunt, M.B.; Norling, C.; Hellens, R.P.; Schaffer, R.J.; et al. Environmental regulation of leaf colour in red 35S:PAP1 Arabidopsis thaliana. New Phytol. 2009, 182, 102–115. [Google Scholar] [CrossRef] [PubMed]
  27. Li, S.F.; Parish, R.W. Isolation of two novel myb-like genes from Arabidopsis and studies on the DNA-binding properties of their products. Plant J. 1995, 8, 963–972. [Google Scholar] [CrossRef] [PubMed]
  28. Zhuang, H.; Chong, S.L.; Priyanka, B.; Han, X.; Lin, E.; Tong, Z.; Huang, H. Full-length transcriptomic identification of R2R3-MYB family genes related to secondary cell wall development in Cunninghamia lanceolata (Chinese fir). BMC Plant Biol. 2021, 21, 581. [Google Scholar] [CrossRef] [PubMed]
  29. Huo, D.; Liu, X.; Zhang, Y.; Duan, J.; Zhang, Y.; Luo, J. A Novel R2R3-MYB Transcription Factor PqMYB4 Inhibited Anthocyanin Biosynthesis in Paeonia qiui. Int. J. Mol. Sci. 2020, 21, 5878. [Google Scholar] [CrossRef]
  30. Matsui, K.; Oshima, Y.; Mitsuda, N.; Sakamoto, S.; Nishiba, Y.; Walker, A.R.; Ohme-Takagi, M.; Robinson, S.P.; Yasui, Y.; Mori, M.; et al. Buckwheat R2R3 MYB transcription factor FeMYBF1 regulates flavonol biosynthesis. Plant Sci. 2018, 274, 466–475. [Google Scholar] [CrossRef]
  31. Lai, B.; Du, L.-N.; Hu, B.; Wang, D.; Huang, X.-M.; Zhao, J.-T.; Wang, H.-C.; Hu, G.-B. Characterization of a novel litchi R2R3-MYB transcription factor that involves in anthocyanin biosynthesis and tissue acidification. BMC Plant Biol. 2019, 19, 62. [Google Scholar] [CrossRef] [PubMed]
  32. Fahima, A.; Levinkron, S.; Maytal, Y.; Hugger, A.; Lax, I.; Huang, X.; Eyal, Y.; Lichter, A.; Goren, M.; Stern, R.A.; et al. Cytokinin treatment modifies litchi fruit pericarp anatomy leading to reduced susceptibility to post-harvest pericarp browning. Plant Sci. 2019, 283, 41–50. [Google Scholar] [CrossRef] [PubMed]
  33. Landoni, M.; Puglisi, D.; Cassani, E.; Borlini, G.; Brunoldi, G.; Comaschi, C.; Pilu, R. Phlobaphenes modify pericarp thickness in maize and accumulation of the fumonisin mycotoxins. Sci. Rep. 2020, 10, 1417. [Google Scholar] [CrossRef] [PubMed]
  34. Pielot, R.; Kohl, S.; Manz, B.; Rutten, T.; Weier, D.; Tarkowská, D.; Rolčík, J.; Strnad, M.; Volke, F.; Weber, H.; et al. Hormone-mediated growth dynamics of the barley pericarp as revealed by magnetic resonance imaging and transcript profiling. J. Exp. Bot. 2015, 66, 6927–6943. [Google Scholar] [CrossRef] [PubMed]
  35. Zhang, X.; Chen, L.; Shi, Q.; Ren, Z. SlMYB102, an R2R3-type MYB gene, confers salt tolerance in transgenic tomato. Plant Sci. 2020, 291, 110356. [Google Scholar] [CrossRef]
  36. Piao, W.; Sakuraba, Y.; Paek, N.-C. Transgenic expression of rice MYB102 (OsMYB102) delays leaf senescence and decreases abiotic stress tolerance in Arabidopsis thaliana. BMB Rep. 2019, 52, 653–658. [Google Scholar] [CrossRef]
  37. Benamar, H.; Rached, W.; Derdour, A.; Marouf, A. Screening of Algerian Medicinal Plants for Acetylcholinesterase Inhibitory Activity. J. Biol. Sci. 2010, 10, 1–9. [Google Scholar] [CrossRef]
  38. Guo, B.; Liu, M.; Yang, H.; Dai, L.; Wang, L. Brassinosteroids Regulate the Water Deficit and Latex Yield of Rubber Trees. Int. J. Mol. Sci. 2023, 24, 12857. [Google Scholar] [CrossRef]
  39. Pertea, M.; Pertea, G.M.; Antonescu, C.M.; Chang, T.-C.; Mendell, J.T.; Salzberg, S.L. StringTie enables improved reconstruction of a transcriptome from RNA-seq reads. Nat. Biotechnol. 2015, 33, 290–295. [Google Scholar] [CrossRef]
  40. Jin, J.; Tian, F.; Yang, D.-C.; Meng, Y.-Q.; Kong, L.; Luo, J.; Gao, G. PlantTFDB 4.0: Toward a central hub for transcription factors and regulatory interactions in plants. Nucleic Acids Res. 2017, 45, D1040–D1045. [Google Scholar] [CrossRef]
  41. Edgar, R.C. MUSCLE: Multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004, 32, 1792–1797. [Google Scholar] [CrossRef] [PubMed]
  42. Nguyen, L.-T.; Schmidt, H.A.; Von Haeseler, A.; Minh, B.Q. IQ-TREE: A fast and effective stochastic algorithm for estimating maximum-likelihood phylogenies. Mol. Biol. Evol. 2015, 32, 268–274. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Volcano map of DEGs in the 3 macadamia varieties with different pericarp thicknesses: A38, Guire No.1, and HAES900.
Figure 1. Volcano map of DEGs in the 3 macadamia varieties with different pericarp thicknesses: A38, Guire No.1, and HAES900.
Ijms 25 06840 g001
Figure 2. Upregulated and downregulated DEGs in the three macadamia varieties with different pericarp thicknesses: A38, Guire No.1, and HAES900.
Figure 2. Upregulated and downregulated DEGs in the three macadamia varieties with different pericarp thicknesses: A38, Guire No.1, and HAES900.
Ijms 25 06840 g002
Figure 3. GO analysis of DEGs in the three macadamia varieties with different pericarp thicknesses: A38, Guire No.1, and HAES900.
Figure 3. GO analysis of DEGs in the three macadamia varieties with different pericarp thicknesses: A38, Guire No.1, and HAES900.
Ijms 25 06840 g003
Figure 4. DEG expression of flavonoid, phenylpropanoid biosynthesis, and cutin, suberine, and wax biosynthesis in the three macadamia varieties with different pericarp thicknesses: A38, Guire No.1, and HAES900.
Figure 4. DEG expression of flavonoid, phenylpropanoid biosynthesis, and cutin, suberine, and wax biosynthesis in the three macadamia varieties with different pericarp thicknesses: A38, Guire No.1, and HAES900.
Ijms 25 06840 g004
Figure 5. Protein structures of 63 MiMYBs for transcriptome analysis in the three macadamia varieties with different pericarp thicknesses: A38, Guire No.1, and HAES900. Green, yellow and pink color box mean Motif 1, Motif 2 and motif 3, respectively.
Figure 5. Protein structures of 63 MiMYBs for transcriptome analysis in the three macadamia varieties with different pericarp thicknesses: A38, Guire No.1, and HAES900. Green, yellow and pink color box mean Motif 1, Motif 2 and motif 3, respectively.
Ijms 25 06840 g005
Figure 6. Phylogenetic and expression analysis of 56 R2R3-MiMYBs in the three macadamia varieties with different pericarp thicknesses: A38, Guire No.1, and HAES900. * means significant difference level (p < 0.05).
Figure 6. Phylogenetic and expression analysis of 56 R2R3-MiMYBs in the three macadamia varieties with different pericarp thicknesses: A38, Guire No.1, and HAES900. * means significant difference level (p < 0.05).
Ijms 25 06840 g006
Figure 7. RT-qPCR verification for 14 MiMYBs’ gene. In each bar graph, the left was the data from transcriptome analysis (FPKM), and the right was the RT-qPCR results (folds). Red, S4 subgroup; Green, S6 subgroup; Blue, S7 subgroup; Cyan, S22 subgroup.
Figure 7. RT-qPCR verification for 14 MiMYBs’ gene. In each bar graph, the left was the data from transcriptome analysis (FPKM), and the right was the RT-qPCR results (folds). Red, S4 subgroup; Green, S6 subgroup; Blue, S7 subgroup; Cyan, S22 subgroup.
Ijms 25 06840 g007
Figure 8. Phenotypes of the three macadamia varieties with different pericarp thicknesses—A38, Guire No.1, and HAES900: (A) phenotype; (B) pericarp thickness, weight, and total flavonoid content. Different uppercase letters mean extremely significant level (p < 0.01).
Figure 8. Phenotypes of the three macadamia varieties with different pericarp thicknesses—A38, Guire No.1, and HAES900: (A) phenotype; (B) pericarp thickness, weight, and total flavonoid content. Different uppercase letters mean extremely significant level (p < 0.01).
Ijms 25 06840 g008
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Tan, Q.; Huan, X.; Pan, Z.; Yang, X.; Wei, Y.; Zhou, C.; Wang, W.; Wang, L. Comparative Transcriptome Analysis Reveals Key Functions of MiMYB Gene Family in Macadamia Nut Pericarp Formation. Int. J. Mol. Sci. 2024, 25, 6840. https://doi.org/10.3390/ijms25136840

AMA Style

Tan Q, Huan X, Pan Z, Yang X, Wei Y, Zhou C, Wang W, Wang L. Comparative Transcriptome Analysis Reveals Key Functions of MiMYB Gene Family in Macadamia Nut Pericarp Formation. International Journal of Molecular Sciences. 2024; 25(13):6840. https://doi.org/10.3390/ijms25136840

Chicago/Turabian Style

Tan, Qiujin, Xiuju Huan, Zhenzhen Pan, Xiaozhou Yang, Yuanrong Wei, Chunheng Zhou, Wenlin Wang, and Lifeng Wang. 2024. "Comparative Transcriptome Analysis Reveals Key Functions of MiMYB Gene Family in Macadamia Nut Pericarp Formation" International Journal of Molecular Sciences 25, no. 13: 6840. https://doi.org/10.3390/ijms25136840

APA Style

Tan, Q., Huan, X., Pan, Z., Yang, X., Wei, Y., Zhou, C., Wang, W., & Wang, L. (2024). Comparative Transcriptome Analysis Reveals Key Functions of MiMYB Gene Family in Macadamia Nut Pericarp Formation. International Journal of Molecular Sciences, 25(13), 6840. https://doi.org/10.3390/ijms25136840

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop