Effector Cs02526 from Ciboria shiraiana induces cell death and modulates plant immunity.

Sclerotinia disease is one of the most devastating fungal diseases worldwide, as it reduces the yields of many economically important crops. Pathogen-secreted effectors play crucial roles in infection processes. However, key effectors of Ciboria shiraiana, the pathogen primarily responsible for sclerotinia disease in mulberry (Morus spp.), remain poorly understood. In this study, we identified and functionally characterized the effector Cs02526 in C. shiraiana and found that Cs02526 could induce cell death in a variety of plants. Moreover, Cs02526-induced cell death was mediated by the central immune regulator brassinosteroid insensitive 1-associated receptor kinase 1, dependent on a 67-amino acid fragment. Notably, Cs02526 homologs were widely distributed in hemibiotrophic and necrotrophic phytopathogenic fungi, but the homologs failed to induce cell death in plants. Pretreatment of plants with recombinant Cs02526 protein enhanced resistance against both C. shiraiana and Sclerotinia sclerotiorum. Furthermore, the pathogenicity of C. shiraiana was diminished upon spraying plants with synthetic dsRNA-Cs02526. In conclusion, our findings highlight the cell death-inducing effector Cs02526 as a potential target for future biological control strategies against plant diseases.

Transcriptomic profiling and discovery of key transcription factors involved in adventitious roots formation from root cuttings of mulberry.

ARs plays a crucial role in plant morphogenesis and development. The limited and inefficient rooting of scions poses a significant challenge to the efficiency and quality of clonal propagation of forest trees in silvicultural practices. Building on previous research conducted by our team, we found that applying IBA at a concentration of 1000 mg/L significantly enhanced mulberry rooting. This study aims to uncover the molecular mechanisms underlying this effect by analyzing RNA sequencing data from mulberry phloem before and after treatment with IBA over time intervals of 10, 20, 30, and 40 days. We identified 5226 DEGs, which were then classified into GO terms and KEGG pathways, showing significant enrichment in hormone signaling processes. Using WGCNA, we identified eight co-expression modules, two of which were significantly correlated with the IBA treatment. Additionally, 18 transcription factors that potentially facilitate ARs formation in mulberry were identified, and an exploratory analysis on the cis-regulatory elements associated with these transcription factors was conducted. The findings of this study provide a comprehensive understanding of the mechanisms of ARs in mulberry and offer theoretical support for the discovery and utilization of exceptional genetic resources within the species.

Molecular and metabolic insights into purplish leaf coloration through the investigation of two mulberry (Morus alba) genotypes.

BACKGROUND: Leaf coloration in plants, attributed to anthocyanin compounds, plays a crucial role in various physiological functions, and also for pharmaceutical and horticultural uses. However, the molecular mechanisms governing leaf coloration and the physiological significance of anthocyanins in leaves remain poorly understood. RESULTS: In this study, we investigated leaf color variation in two closely related mulberry genotypes, one with purplish-red young leaves (EP) and another with normal leaf color (EW). We integrated transcriptomic and metabolomic approaches to gain insights into the metabolic and genetic basis of purplish-red leaf development in mulberry. Our results revealed that flavonoid biosynthesis, particularly the accumulation of delphinidin-3-O-glucoside, is a key determinant of leaf color. Additionally, the up-regulation of CHS genes and transcription factors, including MYB family members, likely contributes to the increased flavonoid content in purplish-red leaves. CONCLUSION: These findings enhance our understanding of the molecular mechanisms responsible for the purplish coloration observed in mulberry leaves and also offer supporting evidence for the hypothesis that anthocyanins serve a protective function in plant tissues until the processes of light absorption and carbon fixation reach maturity, thereby ensuring a balanced equilibrium between energy capture and utilization.

Integrated transcriptomic and metabolomic analyses elucidate the mechanism of flavonoid biosynthesis in the regulation of mulberry seed germination under salt stress.

Seed propagation is the main method of mulberry expansion in China, an important economic forest species. However, seed germination is the most sensitive stage to various abiotic stresses, especially salinity stress. To reveal the molecular regulatory mechanism of mulberry seed germination under salt stress, flavonoid metabolomics and transcriptomics analyses were performed on mulberry seeds germinated under 50 and 100 mmol/L NaCl stress. Analysis of the flavonoid metabolome revealed that a total of 145 differential flavonoid metabolites (DFMs) were classified into 9 groups, 40 flavonols, 32 flavones, 16 chalcones and 14 flavanones. Among them, 61.4% (89) of the DFMs accumulated continuously with increasing salt concentration, reaching the highest level at a 100 mmol/L salt concentration; these DFMs included quercetin-3-O-glucoside (isoquercitrin), kaempferol (3,5,7,4'-tetrahydroxyflavone), quercetin-7-O-glucoside, taxifolin (dihydroquercetin) and apigenin (4',5,7-trihydroxyflavone), indicating that these flavonoids may be key metabolites involved in the response to salt stress. Transcriptional analysis identified a total of 3055 differentially expressed genes (DEGs), most of which were enriched in flavonoid biosynthesis (ko00941), phenylpropanoid biosynthesis (ko00940) and biosynthesis of secondary metabolites (ko01110). Combined analysis of flavonoid metabolomic and transcriptomic data indicated that phenylalanine ammonia-lyase (PAL), 4-coumarate-CoA ligase (4CL), chalcone synthase (CHS), flavonol synthase (FLS), bifunctional dihydroflavonol 4-reductase/flavanone 4-reductase (DFR) and anthocyanidin reductase (ANR) were the key genes involved in flavonoid accumulation during mulberry seed germination under 50 and 100 mmol/L NaCl stress. In addition, three transcription factors, MYB, bHLH and NAC, were involved in the regulation of flavonoid accumulation under salt stress. The results of quantitative real-time PCR (qRT‒PCR) validation showed that the expression levels of 11 DEGs, including 7 genes involved in flavonoid biosynthesis, under different salt concentrations were consistent with the transcriptomic data, and parallel reaction monitoring (PRM) results showed that the expression levels of 6 key enzymes (proteins) involved in flavonoid synthesis were consistent with the accumulation of flavonoids. This study provides a new perspective for investigating the regulatory role of flavonoid biosynthesis in the regulation of mulberry seed germination under salt stress at different concentrations.

Mining and functional characterization of NADPH-cytochrome P450 reductases of the DNJ biosynthetic pathway in mulberry leaves.

BACKGROUND: 1-Deoxynojirimycin (DNJ), the main active ingredient in mulberry leaves, with wide applications in the medicine and food industries due to its significant functions in lowering blood sugar, and lipids, and combating viral infections. Cytochrome P450 is a key enzyme for DNJ biosynthesis, its activity depends on the electron supply of NADPH-cytochrome P450 reductases (CPRs). However, the gene for MaCPRs in mulberry leaves remains unknown. RESULTS: In this study, we successfully cloned and functionally characterized two key genes, MaCPR1 and MaCPR2, based on the transcriptional profile of mulberry leaves. The MaCPR1 gene comprised 2064 bp, with its open reading frame (ORF) encoding 687 amino acids. The MaCPR2 gene comprised 2148 bp, and its ORF encoding 715 amino acids. The phylogenetic tree indicates that MaCPR1 and MaCPR2 belong to Class I and Class II, respectively. In vitro, we found that the recombinant enzymes MaCPR2 protein could reduce cytochrome c and ferricyanide using NADPH as an electron donor, while MaCPR1 did not. In yeast, heterologous co-expression indicates that MaCPR2 delivers electrons to MaC3'H hydroxylase, a key enzyme catalyzing the production of chlorogenic acid from 3-O-p-coumaroylquinic acid. CONCLUSIONS: These findings highlight the orchestration of hydroxylation process mediated by MaCPR2 during the biosynthesis of secondary metabolite biosynthesis in mulberry leaves. These results provided a foundational understanding for fully elucidating the DNJ biosynthetic pathway within mulberry leaves.

Phylogenomics of mulberries (Morus, Moraceae) inferred from plastomes and single copy nuclear genes.

Mulberries (genus Morus), belonging to the order Rosales, family Moraceae, are important woody plants due to their economic values in sericulture, as well as for nutritional benefits and medicinal values. However, the taxonomy and phylogeny of Morus, especially for the Asian species, remains challenging due to its wide geographical distribution, morphological plasticity, and interspecific hybridization. To better understand the evolutionary history of Morus, we combined plastomes and a large-scale nuclear gene analyses to investigate their phylogenetic relationships. We assembled the plastomes and screened 211 single-copy nuclear genes from 13 Morus species and related taxa. The plastomes of Morus species were relatively conserved in terms of genome size, gene content, synteny, IR boundary and codon usage. Using nuclear data, our results elucidated identical topologies based on coalescent and concatenation methods. The genus Morus was supported as monophyletic, with M. notabilis as the first diverging lineage and the two North American Morus species, M. microphylla and M. rubra, as sister to the other Asian species. In the Asian Morus species, interspecific relationships were completely resolved. However, cyto-nuclear discordances and gene tree-species tree conflicts were detected in the phylogenies of Morus, with multiple evidences supporting hybridization/introgression as the main cause of discordances between nuclear and plastid phylogenies, while gene tree-species tree conflicts were mainly caused by ILS.

Mno-miR164a and MnNAC100 regulate the resistance of mulberry to Botrytis cinerea.

Although microRNAs (miRNAs) regulate the defense response of a variety of plant species against a variety of pathogenic fungi, the involvement of miRNAs in mulberry's defense against Botrytis cinerea has not yet been documented. In this study, we identified responsive B. cinerea miRNA mno-miR164a in mulberry trees. After infection with B. cinerea, the expression of mno-miR164a was reduced, which was fully correlated with the upregulation of its target gene, MnNAC100, responsible for encoding a transcription factor. By using transient infiltration/VIGS mulberry that overexpressed mno-miR164a or knocked-down MnNAC100, our study revealed a substantial enhancement in mulberry's resistance to B. cinerea when mno-miR164a was overexpressed or MnNAC100 expression was suppressed. This enhancement was accompanied by increased catalase (CAT) activity and reduced malondialdehyde (MDA) content. In addition, mno-miR164a-mediated inhibition of MnNAC100 enhanced the expression of a cluster of defense-related genes in transgenic plants upon exposure to B. cinerea. Meanwhile, MnNAC100 acts as a transcriptional repressor, directly suppressing the expression of MnPDF1.2. Our study indicated that the mno-miR164a-MnNAC100 regulatory module manipulates the defense response of mulberry to B. cinerea infection. This discovery has great potential in breeding of resistant varieties and disease control.

Rhizophagus Irregularis regulates flavonoids metabolism in paper mulberry roots under cadmium stress.

Broussonetia papyrifera is widely found in cadmium (Cd) contaminated areas, with an inherent enhanced flavonoids metabolism and inhibited lignin biosynthesis, colonized by lots of symbiotic fungi, such as arbuscular mycorrhizal fungi (AMF). However, the physiological and molecular mechanisms by which Rhizophagus irregularis, an AM fungus, regulates flavonoids and lignin in B. papyrifera under Cd stress remain unclear. Here, a pot experiment of B. papyrifera inoculated and non-inoculated with R. irregularis under Cd stress was carried out. We determined flavonoids and lignin concentrations in B. papyrifera roots by LC-MS and GC-MS, respectively, and measured the transcriptional levels of flavonoids- or lignin-related genes in B. papyrifera roots, aiming to ascertain the key components of flavonoids or lignin, and key genes regulated by R. irregularis in response to Cd stress. Without R. irregularis, the concentrations of eriodictyol, quercetin and myricetin were significantly increased under Cd stress. The concentrations of eriodictyol and genistein were significantly increased by R. irregularis, while the concentration of rutin was significantly decreased. Total lignin and lignin monomer had no alteration under Cd stress or with R. irregularis inoculation. As for flavonoids- or lignin-related genes, 26 genes were co-regulated by Cd stress and R. irregularis. Among these genes, BpC4H2, BpCHS8 and BpCHI5 were strongly positively associated with eriodictyol, indicating that these three genes participate in eriodictyol biosynthesis and were involved in R. irregularis assisting B. papyrifera to cope with Cd stress. This lays a foundation for further research revealing molecular mechanisms by which R. irregularis regulates flavonoids synthesis to enhance tolerance of B. papyrifera to Cd stress.

Unraveling the Contribution of MulSOS2 in Conferring Salinity Tolerance in Mulberry (Morus atropurpurea Roxb).

Salinity is one of the most serious threats to sustainable agriculture. The Salt Overly Sensitive (SOS) signaling pathway plays an important role in salinity tolerance in plants, and the SOS2 gene plays a critical role in this pathway. Mulberry not only has important economic value but also is an important ecological tree species; however, the roles of the SOS2 gene associated with salt stress have not been reported in mulberry. To gain insight into the response of mulberry to salt stress, SOS2 (designated MulSOS2) was cloned from mulberry (Morus atropurpurea Roxb), and sequence analysis of the amino acids of MulSOS2 showed that it shares some conserved domains with its homologs from other plant species. Our data showed that the MulSOS2 gene was expressed at different levels in different tissues of mulberry, and its expression was induced substantially not only by NaCl but also by ABA. In addition, MulSOS2 was exogenously expressed in Arabidopsis, and the results showed that under salt stress, transgenic MulSOS2 plants accumulated more proline and less malondialdehyde than the wild-type plants and exhibited increased tolerance to salt stress. Moreover, the MulSOS2 gene was transiently overexpressed in mulberry leaves and stably overexpressed in the hairy roots, and similar results were obtained for resistance to salt stress in transgenic mulberry plants. Taken together, the results of this study are helpful to further explore the function of the MulSOS2 gene, which provides a valuable gene for the genetic breeding of salt tolerance in mulberry.

Metabolic and Transcriptional Analysis Reveals Flavonoid Involvement in the Drought Stress Response of Mulberry Leaves.

Abiotic stress, especially drought stress, poses a significant threat to terrestrial plant growth, development, and productivity. Although mulberry has great genetic diversity and extensive stress-tolerant traits in agroforestry systems, only a few reports offer preliminary insight into the biochemical responses of mulberry leaves under drought conditions. In this study, we performed a comparative metabolomic and transcriptomic analysis on the "drooping mulberry" (Morus alba var. pendula Dippel) under PEG-6000-simulated drought stress. Our research revealed that drought stress significantly enhanced flavonoid accumulation and upregulated the expression of phenylpropanoid biosynthetic genes. Furthermore, the activities of superoxide dismutase (SOD), catalase (CAT) and malondialdehyde (MDA) content were elevated. In vitro enzyme assays and fermentation tests indicated the involvement of flavonol synthase/flavanone 3-hydroxylase (XM_010098126.2) and anthocyanidin 3-O-glucosyltransferase 5 (XM_010101521.2) in the biosynthesis of flavonol aglycones and glycosides, respectively. The recombinant MaF3GT5 protein was found to recognize kaempferol, quercetin, and UDP-glucose as substrates but not 3-/7-O-glucosylated flavonols and UDP-rhamnose. MaF3GT5 is capable of forming 3-O- and 7-O-monoglucoside, but not di-O-glucosides, from kaempferol. This implies its role as a flavonol 3, 7-O-glucosyltransferase. The findings from this study provided insights into the biosynthesis of flavonoids and could have substantial implications for the future diversified utilization of mulberry.

Post-harvest ultraviolet irradiation induces changes in physical-chemical properties and levels of polycyclic aromatic hydrocarbons and gene expression in mulberry fruit.

BACKGROUND: Earlier studies reported that post-harvest ultraviolet (UV) irradiation could increase the health-promoting compounds in fruit but the effects of UV irradiation on the reduction of the polycyclic aromatic hydrocarbon (PAH) content in mulberries remain less known. Black mulberry fruit were exposed to two UV illumination dosages (3.5 and 7 kJ m-2 ) and were stored for 4, 8, and 12 days. RESULTS: Mulberries treated in this way displayed higher antioxidant enzyme activity and phenolic compound content in comparison with a control condition. The transcription factors (TFs) MdoMYB121, MdoMYB155, MdbZIP2, and MdbZIP48 were strongly expressed in two UV illumination dosages (about 45-95% higher than the control). The fluorine (Flu) and naphthalene (Nap) content in treated fruit decreased by 21-85% in comparison with the control condition. CONCLUSION: The findings of this study indicate that UV irradiation can be considered as a promising technique to remove some PAHs in black mulberries, to increase their health-promoting potential, and indirectly to improve their aesthetic quality due to the resulting desirable color parameters. © 2023 Society of Chemical Industry.

Chromosome constitution and genetic relationships of Morus spp. revealed by genomic in situ hybridization.

BACKGROUND: Mulberry (Morus spp.) is an economically important woody plant, which has been used for sericulture (silk farming) for thousands of years. The genetic background of mulberry is complex due to polyploidy and frequent hybridization events. RESULTS: Comparative genomic in situ hybridization (cGISH) and self-GISH were performed to illustrate the chromosome constitution and genetic relationships of 40 mulberry accessions belonging to 12 species and three varietas in the Morus genus and containing eight different ploidy levels. We identified six homozygous cGISH signal patterns and one heterozygous cGISH signal pattern using four genomic DNA probes. Using cGISH and self-GISH data, we defined five mulberry sections (Notabilis, Nigra, Wittiorum, and Cathayana, all contained only one species; and Alba, which contained seven closely related species and three varietas, was further divided into two subsections) and proposed the genetic relationships among them. Differential cGISH signal patterns detected in section Alba allowed us to refine the genetic relationships among the closely related members of this section. CONCLUSIONS: We propose that GISH is an efficient tool to investigate the chromosome constitution and genetic relationships in mulberry. The results obtained here can be used to guide outbreeding of heterozygous perennial crops like mulberry.

Morpho-cultural, pathological and molecular variability in Phloeospora maculans causing leaf spot of mulberry (Morus species) in India.

BACKGROUND: Leaf spot disease (LSD) of mulberry caused by Phloeospora maculans is a major threat to the silk industry of Jammu and Kashmir, India, therefore, it was necessary to study the population structure of the pathogen for successful management of the disease. METHODS AND RESULTS: To understand the diversity in the Phloeospora maculans, a combination of conventional (morphological, cultural and pathological) and molecular (ISSR markers) approaches were employed to discern the variability in 27 isolates collected from Srinagar, Bandipora, and Baramulla districts of Jammu and Kashmir, India. The studies revealed a high level of variability in the pathogen. Based on the morpho-cultural and pathological studies, the pathogen isolates were grouped into different categories based on colony growth, texture, margin and colour besides changes in colour of medium, incubation period, leaf area infected, etc.A high level of polymorphism was observed in different isolates of P. maculans using ISSR markers, which indicated that these markers are suitable for studying the genetic diversity in this pathogen. All the isolates (27) of P. maculans were clustered into two groups or populations as indicated by mean delta K value. Analysis of molecular variance revealed the low genetic variation among the populations (1.08%) and a high level of genetic variation within the populations (98.91%). Fst value was found to be 0.01 indicating smaller amount of genetic differentiation between the populations against calculated P-value of 0.29. CONCLUSION: A high level of diversity based on morphological, cultural, pathological and molecular levels was observed in Phloeospora maculans collected from various districts of Kashmir valley, which indicates that the study of population structure is necessary for successful management of the disease.

Diffusion and enrichment of high-risk antibiotic resistance genes (ARGs) via the transmission chain (mulberry leave, guts and feces of silkworm, and soil) in an ecological restoration area of manganese mining, China: Role of heavy metals.

This study investigated the diffusion and enrichment of antibiotic resistance genes (ARGs) and pathogens via the transmission chain (mulberry leaves - silkworm guts - silkworm feces - soil) near a manganese mine restoration area (RA) and control area (CA, away from RA). Horizontal gene transfer (HGT) of ARGs was testified by an IncP a-type broad host range plasmid RP4 harboring ARGs (tetA) and conjugative genes (e.g., korB, trbA, and trbB) as an indicator. Compared to leaves, the abundances of ARGs and pathogens in feces after silkworms ingested leaves from RA increased by 10.8% and 52.3%, respectively, whereas their abundance in feces from CA dropped by 17.1% and 97.7%, respectively. The predominant ARG types in feces involved the resistances to ß-lactam, quinolone, multidrug, peptide, and rifamycin. Therein, several high-risk ARGs (e.g., qnrB, oqxA, and rpoB) carried by pathogens were more enriched in feces. However, HGT mediated by plasmid RP4 in this transmission chain was not a main factor to promote the enrichment of ARGs due to the harsh survival environment of silkworm guts for the plasmid RP4 host E. coli. Notably, Zn, Mn, and As in feces and guts promoted the enrichment of qnrB and oqxA. Worriedly, the abundance of qnrB and oqxA in soil increased by over 4-fold after feces from RA were added into soil for 30 days regardless of feces with or without E. coli RP4. Overall, ARGs and pathogens could diffuse and enrich in environment via the sericulture transmission chain developed at RA, especially some high-risk ARGs carried by pathogens. Thus, greater attentions should be paid to dispel such high-risk ARGs to support benign development of sericulture industry in the safe utilization of some RAs.

Effects of artificial diet breeding on intestinal microbial populations at the young stage of silkworm (Bombyx mori).

The silkworm (Bombyx mori) is an economically important insect and serves as a model organism for Lepidoptera. To investigate the effects of the intestinal microbial population on the growth and development of larvae fed an artificial diet (AD) during the young stages, we analyzed the characteristics of the intestinal microbial population using 16S rRNA gene sequencing technology. Our results revealed that the intestinal flora of the AD group tended to be simple by the 3rd-instar, which Lactobacillus accounting for 14.85% and leading to a decreased pH in the intestinal fluid. In contrast, the intestinal flora of silkworms in the mulberry leaf (ML) group showed continuous growth of diversity, with Proteobacteria accounting for 37.10%, Firmicutes accounting for 21.44%, and Actinobacteria accounting for 17.36%. Additionally, we detected the activity of intestinal digestive enzymes at different instars and found that the activity of digestive enzymes in the AD group increased by larval instar. Protease activity in the AD group was lower during the 1st- to 3rd-instars compared to the ML group, while α-amylase and lipase activities were significantly higher in the AD group during the 2nd- and 3rd-instar compared to the ML group. Furthermore, our experimental results indicated that changes in the intestinal population decreased the pH and affected the activity of proteases, which might contribute to the slower growth and development of larvae in the AD group. In summary, this study provides a reference for investigating the relationship between artificial diet and intestinal flora balance.

Novel aldo-keto reductase AKR2E9 regulates aldehyde content in the midgut and antennae of the silkworm (Bombyx mori).

We studied the effects of green leaf volatiles (including reactive aldehydes) emitted by plants on insects that feed on these plants. The silkworm (Bombyx mori) is a model lepidopteran that eats mulberry leaves. Defense-related enzymes in silkworms can be targeted for developing new pest control methods. The aldo-keto reductase (AKR) superfamily catalyzes aldehyde reduction by converting a carbonyl group into an alcohol group. Here, we characterized a novel silkworm AKR, designated as AKR2E9. Recombinant AKR2E9 was overexpressed in Escherichia coli. The recombinant protein was used, along with nicotinamide adenine dinucleotide phosphate as a coenzyme, to reduce aldehydes present in mulberry (Morus alba) leaves. The catalytic efficiency of AKR2E9 toward various aldehyde substrates and its inhibitor sensitivity was lower than those of AKR2E8. High expression levels of akr2e9 messenger RNA (mRNA) were detected in the midgut and antennae of silkworms. In the antennae of adult silkworms, akr2e9 mRNA was more abundant than akr2e8 mRNA. The catalytic efficiency of AKR2E9 was low because of steric hindrance, due to which its active site is blocked. High expression levels of AKR2E9 in the midgut and antennae suggest that it may regulate the detoxification of toxic aldehydes in silkworms.

Draft genome sequence of Indian mulberry (Morus indica) provides a resource for functional and translational genomics.

Mulberry is an important crop plant for the sericulture industry. Here, we report high-quality genome sequence of a cultivated Indian mulberry (Morus indica cv K2) obtained by combining data from four different technologies, including Illumina, single-molecule real-time sequencing, chromosome conformation capture and optical mapping, with a gene completeness of 96.5%. Based on the genome sequence, we identified 49.2% of repetitive DNA and 27,435 high-confidence protein-coding genes with >90% of them supported by transcript evidence. A comparative analysis with other plant genomes identified 4.8% of species-specific genes in the M. indica genome. Transcriptome profiling revealed tissue-specific and differential expression across multiple accessions of ~4.7% and 2-5% of protein-coding genes, respectively, implicated in diverse biological processes. Whole genome resequencing of 21 accessions/species revealed ~2.5 million single nucleotide polymorphisms and ~ 0.2 million insertions/deletions. These data and results provide a comprehensive resource to accelerate the genomics research in mulberry for its improvement.

A Feature Selection Approach Guided an Early Prediction of Anthocyanin Accumulation Using Massive Untargeted Metabolomics Data in Mulberry.

Identifying the early predictive biomarkers or compounds represents a pivotal task for guiding a targeted agricultural practice. Despite the various available tools, it remains challenging to define the ideal compound combination and thereby elaborate an effective predictive model fitting that. Hence, we employed a stepwise feature selection approach followed by a maximum relevance and minimum redundancy (MRMR) on the untargeted metabolism in four mulberry genotypes at different fruit developmental stages (FDSs). Thus, we revealed that 7 out of 226 differentially abundant metabolites (DAMs) explained up to 80% variance of anthocyanin based on linear regression model and stepwise feature selection approach accompanied by an MRMR across the genotypes over the FDSs. Among them, the phosphoenolpyruvate, d-mannose and shikimate show the top 3 attribution indexes to the accumulation of anthocyanin in the fruits of these genotypes across the four FDSs. The obtained results were further validated by assessing the regulatory genes expression levels and the targeted metabolism approach. Taken together, our findings provide valuable evidences on the fact that the anthocyanin biosynthesis is somehow involved in the coordination between the carbon metabolism and secondary metabolic pathway. Our report highlights as well the importance of using the feature selection approach for the predictive biomarker identification issued from the untargeted metabolomics data.

Comprehensive analysis of the MYB transcription factor gene family in Morus alba.

BACKGROUND: The V-myb myeloblastosis viral oncogene homolog (MYB) family of proteins is large, containing functionally diverse transcription factors. However, MYBs in Morus are still poorly annotated and a comprehensive functional analysis of these transcription factors is lacking. RESULTS: In the present study, a genome-wide identification of MYBs in Morus alba was performed. In total 166 MaMYBs were identified, including 103 R2R3-MYBs and four 3R-MaMYBs. Comprehensive analyses, including the phylogenetic analysis with putative functional annotation, motif and structure analysis, gene structure organization, promoter analysis, chromosomal localization, and syntenic relationships of R2R3-MaMYBs and 3R-MaMYBs, provided primary characterization for these MaMYBs. R2R3-MaMYBs covered the subgroups reported for R2R3-MYBs in Arabidopsis and Populus, and had two Morus-specific subgroups, indicating the high retention of MYBs in Morus. Motif analysis revealed high conservative residues at the start and end of each helix and residues consisting of the third helix in R2 and R3 repeats. Thirteen intron/exon patterns (a-m) were summarized, and the intron/exon pattern of two introns with phase numbers of 0 and 2 was the prevalent pattern for R2R3-MaMYBs. Various cis-elements in promoter regions were identified, and were mainly related to light response, development, phytohormone response, and abiotic and biotic stress response and secondary metabolite production. Expression patterns of R2R3-MaMYBs in different organs showed that MaMYBs involved in secondary cell wall components and stress responsiveness were preferentially expressed in roots or stems. R2R3-MaMYBs involved in flavonoid biosynthesis and anthocyanin accumulation were identified and characterized based on functional annotation and correlation of their expression levels with anthocyanin contents. CONCLUSION: Based on a comprehensive analysis, this work provided functional annotation for R2R3-MYBs and an informative reference for further functional dissection of MYBs in Morus.

Molecular mechanism of mulberry response to drought stress revealed by complementary transcriptomic and iTRAQ analyses.

BACKGROUND: The use of mulberry leaves has long been limited to raising silkworms, but with the continuous improvement of mulberry (Morus alba) resource development and utilization, various mulberry leaf extension products have emerged. However, the fresh leaves of mulberry trees have a specific window of time for picking and are susceptible to adverse factors, such as drought stress. Therefore, exploring the molecular mechanism by which mulberry trees resist drought stress and clarifying the regulatory network of the mulberry drought response is the focus of the current work. RESULTS: In this study, natural and drought-treated mulberry grafted seedlings were used for transcriptomic and proteomic analyses (CK vs. DS9), aiming to clarify the molecular mechanism of the mulberry drought stress response. Through transcriptome and proteome sequencing, we identified 9889 DEGs and 1893 DEPs enriched in stress-responsive GO functional categories, such as signal transducer activity, antioxidant activity, and transcription regulator activity. KEGG enrichment analysis showed that a large number of codifferentially expressed genes were enriched in flavonoid biosynthesis pathways, hormone signalling pathways, lignin metabolism and other pathways. Through subsequent cooperation analysis, we identified 818 codifferentially expressed genes in the CK vs. DS9 comparison group, including peroxidase (POD), superoxide dismutase (SOD), aldehyde dehydrogenase (ALDHs), glutathione s-transferase (GST) and other genes closely related to the stress response. In addition, we determined that the mulberry gene MaWRKYIII8 (XP_010104968.1) underwent drought- and abscisic acid (ABA)-induced expression, indicating that it may play an important role in the mulberry response to drought stress. CONCLUSIONS: Our research shows that mulberry can activate proline and ABA biosynthesis pathways and produce a large amount of proline and ABA, which improves the drought resistance of mulberry. MaWRKYIII8 was up-regulated and induced by drought and exogenous ABA, indicating that MaWRKYIII8 may be involved in the mulberry response to drought stress. These studies will help us to analyse the molecular mechanism underlying mulberry drought tolerance and provide important gene information and a theoretical basis for improving mulberry drought tolerance through molecular breeding in the future.