Mulberry fruits with high concentrations of anthocyanins are favored by consumers because of their good taste, bright color, and high nutritional value. However, neither the regulatory mechanism controlling flavonoid biosynthesis in mulberry nor the molecular basis of different mulberry fruit colors is fully understood. Here, we report that a flavonoid homeostasis network comprising activation and feedback regulation mechanisms determines mulberry fruit color. In vitro and in vivo assays showed that MYBA-bHLH3-TTG1 regulates the biosynthesis of anthocyanins, while TT2L1 and TT2L2 work with bHLH3 or GL3 and form a MYB-bHLH-WD40 (MBW) complex with TTG1 to regulate proanthocyanidin (PA) synthesis. Functional and expression analyses showed that bHLH3 is a key regulator of the regulatory network controlling mulberry fruit coloration and that MYB4 is activated by MBW complexes and participates in negative feedback control of the regulatory network to balance the accumulation of anthocyanins and proanthocyanidins. Our research demonstrates that the interaction between bHLH3 and MYB4 in the homeostasis regulatory network ensures that the fruits accumulate desirable flavonoids and that this network is stable in pigmentrich mulberry fruits. However, the abnormal expression of bHLH3 disrupts the balance of the network and redirects flavonoid metabolic flux in pale-colored fruits, resulting in differences in the levels and proportions of anthocyanins, flavones, and flavonols among differently colored mulberry fruits (red, yellow, and white). The results of our study reveal the molecular basis of the diversity of mulberry fruit colors.

No abstract

Background Miniature inverted-repeat transposable elements (MITEs) are common in eukaryotic genomes, and are important for genomic evolution. Results In the present study, the identification of MITEs in the mulberry genome revealed 286,122 MITE-related sequences, including 90,789 full-length elements. The amplification of mulberry MITEs and the influence of MITEs on the evolution of the mulberry genome were analyzed. The timing of MITE amplifications varied considerably among the various MITE families. Fifty-one MITE families have undergone a single round of amplification, while the other families developed from multiple amplifications. Most mulberry MITEs were inserted near genes and some could regulate gene expression through small RNAs. An analysis of transgenic plants indicated that MITE insertions can upregulate the expression of a target gene. Moreover, MITEs are frequently associated with alternative splicing events (exonizations). Conclusion The data presented herein provide insights into the generation of MITEs as well as their impact on gene regulation and evolution in mulberry. Electronic supplementary material The online version of this article (10.1186/s13100-019-0169-0) contains supplementary material, which is available to authorized users.

DNA methylation has been proposed to regulate plant stress resistance. However, the dynamic changes in DNA methylation in woody plants and their correlations with pathogenic responses are not fully understood. Here, we present single-base maps of the DNA methylomes of mulberry (Morus notabilis) leaves that were subjected to a mock treatment or inoculation with Botrytis cinerea. Compared with the former, the latter showed decreased mCG and mCHG levels and increased mCHH levels. DNA methylation inhibitors reduced resistance gene methylation levels and enhanced mulberry resistance, suggesting that the hypomethylation of resistance genes affects mulberry resistance to B. cinerea. Virus-induced gene silencing of MnMET1 enhanced the expression of mulberry-resistance genes, thereby increasing the plant’s resistance to B. cinerea. We also found that MITEs play a dominant role in controlling DNA methylation levels. MITEs appear to be the main sources of 24-nt siRNAs that regulate gene expression through the RNA-directed DNA methylation pathway.

Core Ideas Percentage coverage of TEs is significantly negatively correlated with that of genes. Different TE superfamilies exhibit distinct distribution patterns in mulberry genome. Copia elements may have a dominant influence on the regulation of mulberry genes. TE‐containing genes assigned to pathways were mainly in metabolism‐related pathways. TEs are a driving force in the formation of alternatively spliced genes. Mulberry (Morus notabilis C. K. Schneid) leaves have been used as the food for the domesticated silkworm, Bombyx mori, for more than 5000 yr, and the mulberry–silkworm relationship is one of the best‐known and oldest models of plant defense–insect adaptation. The availability of a genome assembly of mulberry provides us with an opportunity to mine the characteristics and distribution of transposable elements (TEs) in this species and to examine their relationship to genes and gene expression. In this study, a significantly correlated inverse relationship between the percentage coverage of genes and TEs was observed. The TE‐rich regions appeared to have a lower percentage of putatively expressed genes. Distribution patterns between different TE superfamilies were detected in the mulberry genome. The Copia elements (the TE making up the greatest proportion of the mulberry genome) were significantly overrepresented within genes in the mulberry genome, and they may have a dominant influence on evolution of the mulberry genome. Approximately 96.93% (330/344) of the TE‐containing genes assigned to pathways were assigned to metabolism‐related pathways. The TE‐related alternative splicing events accounted for 7.58% (402/5,302) of all alternative splicing types in the mulberry genome, suggesting that TEs are one of the driving forces in the formation of the alternatively spliced genes. The results will be valuable in improving our understanding of the important roles of TEs in mulberry genome evolution.

The evolutionary dynamics of long terminal repeat (LTR) retrotransposons in tree genomes has remained largely unknown. The availability of the complete genome sequences of the mulberry tree (Morus notabilis) has offered an unprecedented opportunity for us to characterize these retrotransposon elements. We investigated 202 and 114 families of Copia and Gypsy superfamilies, respectively, comprising 2916 intact elements in the mulberry genome. The tRNAMet was the most frequently used type of tRNA in both superfamilies. Phylogenetic analysis suggested that Copia and Gypsy from mulberry can be grouped into eight and six lineages, respectively. All previously characterized families of such elements could also be found in the mulberry genome. About 95% of the identified Copia and Gypsy full elements were estimated to have been inserted into the mulberry genome within the past 2–3 million years. Meanwhile, the estimated insertion times of members of the three most abundant families of the Copia superfamily (908 members from the three most abundant families) and Gypsy superfamily (783 members from the three most abundant families) revealed divergent life histories. Compared with the situation in Gypsy elements, three families of Copia elements are under positive selection pressure, which suggested that Copia elements may have a dominant influence in the evolution of mulberry genes. Analysis of insertion and deletion dynamics suggested that Copia and Gypsy elements exhibited a very long half-life in the mulberry genome. The present work provides new insights into the insertion and deletion dynamics of LTR retrotransposons, and it will greatly improve our understanding of the important roles transposable elements play in the architecture of the mulberry genome.

Chitinase is a hydrolase that uses chitin as a substrate. It plays an important role in plant resistance to fungal pathogens by degrading chitin. Here, we conducted bioinformatics analysis and transcriptome data analysis of the mulberry (Morus notabilis) chitinase gene family to determine its role in the resistance to Botrytis cinerea. A total of 26 chitinase genes were identified, belonging to the GH18 and GH19 families. Among them, six chitinase genes were differentially expressed under the infection of B. cinerea. MnChi18, which significantly responded to B. cinerea, was heterologously expressed in Arabidopsis (Arabidopsis thaliana). The resistance of MnChi18 transgenic Arabidopsis to B. cinerea was significantly enhanced, and after inoculation with B. cinerea, the activity of catalase (CAT) increased and the content of malondialdehyde (MDA) decreased. This shows that overexpression of MnChi18 can protect cells from damage. In addition, our study also indicated that MnChi18 may be involved in B. cinerea resistance through other resistance-related genes. This study provides an important basis for further understanding the function of mulberry chitinase.

Horizontal transposable element transfer (HTT) events have occurred among a large number of species and play important roles in the composition and evolution of eukaryotic genomes. HTTs are also regarded as effective forces in promoting genomic variation and biological innovation. In the present study, HTT events were identified and analyzed in seven sequenced species of Rosales using bioinformatics methods by comparing sequence conservation and K/K value of reverse transcriptase (RT) with 20 conserved genes, estimating the dating of HTTs, and analyzing the phylogenetic relationships. Seven HTT events involving long terminal repeat (LTR) retrotransposons, two HTTs between Morus notabilis and Ziziphus jujuba, and five between Malus domestica and Pyrus bretschneideri were identified. Further analysis revealed that these LTR retrotransposons had functional structures, and the copy insertion times were lower than the dating of HTTs, particularly in Mn.Zj.1 and Md.Pb.3. Altogether, the results demonstrate that LTR retrotransposons still have potential transposition activity in host genomes. These results indicate that HTT events are another strategy for exchanging genetic material among species and are important for the evolution of genomes.