Chromosome-Level Reference Genome and Population Genomic Analysis Provide Insights into the Evolution and Improvement of Domesticated Mulberry (Morus alba)

Safe utilization and remediation potential of the mulberry-silkworm system in heavy metal-contaminated lands: A review

Mulberry cultivation and silkworm rearing hold a prominent position in the agricultural industries of many Asian countries, contributing to economic growth, sustainable development, and cultural heritage preservation. Applying the soil-mulberry-silkworm system (SMSS) to heavy metal (HM)-contaminated areas is significant economically, environmentally, and socially. The ultimate goal of this paper is to review the main research progress of SMSS under HM stress, examining factors affecting its safe utilization and remediation potential for HM-contaminated soils. HM tolerance of mulberry and silkworms relates to their growth stages. Based on the standards for HM contaminants in various mulberry and silkworm products and the bioconcentration factor of HMs at different parts of SMSS, we calculated maximum safe Cd and Pb levels for SMSS application on contaminated lands. Several remediation practices demonstrated mulberry's ability to grow on barren lands, absorb various HMs, while silkworm excreta can adsorb HMs and improve soil fertility. Considering multiple factors influencing HM tolerance and accumulation, we propose a decision model to guide SMSS application in polluted areas. Finally, we discussed the potential of using molecular breeding techniques to screen or develop varieties better suited for HM-contaminated regions. However, actual pollution scenarios are often complex, requiring consideration of multiple factors. More large-scale applications are crucial to enhance the theoretical foundation for applying SMSS in HM pollution risk areas.

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.

Genome-wide analysis of the mulberry (Morus abla L.) GH9 gene family and the functional characterization of MaGH9B6 during the development of the abscission zone

Plant glycoside hydrolase family 9 genes (GH9s) are widely distributed in plants and involved in a variety of cellular and physiological processes. In the current study, nine GH9 genes were identified in the mulberry and were divided into two subfamilies based on the phylogenetic analysis. Conserved motifs and gene structure analysis suggested that the evolution of the two subfamilies is relatively conserved and the glycoside hydrolase domain almost occupy the entire coding region of the GH9s gene. Only segmental duplication has played a role in the expansion of gene family. Collinearity analysis showed that mulberry GH9s had the closest relationship with poplar GH9s. MaGH9B1, MaGH9B6, MaGH9B5, and MaGH9B3 were detected to have transcript accumulation in the stalk of easy-to drop mature fruit drop, suggesting that these could play a role in mulberry fruit drop. Multiple cis-acting elements related to plant hormones and abiotic stress responses were found in the mulberry GH9 promoter regions and showed different activities under exogenous abscisic acid (ABA) and 2,4- dichlorophenoxyacetic acid (2,4-D) stresses. We found that the lignin content in the fruit stalk decreased with the formation of the abscission zone (AZ), which could indirectly reflect the formation process of the AZ. These results provide a theoretical basis for further research on the role of GH9s in mulberry abscission.

Grafting seedling rootstock strengthens tolerance to drought stress in polyploid mulberry (Morus alba L.)

The economically adaptable mulberry (Morus alba L.) has a long history of grafting in China, yet the physiological mechanisms and advantages in drought tolerance remain unexplored. In our study, we investigated the responses of self-rooted 2X (diploid), 3X (triploid), and 4X (tetraploid) plants, as well as polyploid plants grafted onto diploid seedling rootstocks (2X/2X, 3X/2X, and 4X/2X) under drought stress. We found that self-rooted diploid plants exhibited the most severe phenotypic damage, lowest water retention, photosynthetic capacity, and the least effective osmotic stress adjustment compared to tetraploid and triploid plants. However, grafted diploid and triploid plants showed effective mitigation of drought-induced damage, with higher relative water content and improved soil water retention. Grafted plants also improved the photosystem response to drought stress through elevated photosynthetic potential, closed stomatal aperture, and faster recovery of chlorophyll biosynthesis in the leaves. Additionally, grafted plants altered osmotic protective compound levels, including starch, soluble sugar, and proline content, thereby enhancing drought resistance. Absolute quantification PCR indicated that the expression levels of proline synthesis-related genes in grafted plants were not influenced after drought stress, whereas they were significantly increased in self-rooted plants. Consequently, our findings support that self-rooted triploid and tetraploid mulberries exhibited superior drought resistance compared to diploid plants. Moreover, grafting onto seedling rootstocks enhanced tolerance against drought stress in diploid and triploid mulberry, but not in tetraploid. Our study provides valuable insights for a comprehensive analysis of physiological effects in response to drought stress between stem-roots and seedling rootstocks.

Antibacterial Activity and Mechanism of Three Root Exudates from Mulberry Seedlings against Ralstonia pseudosolanacearum

Bacterial wilt is a significant soil-borne disease that poses a threat to mulberry production yield and quality of agricultural production worldwide. However, the disease resistance mechanisms dependent on root exudates are not well understood. In this present study, we investigated the antibacterial mechanisms of the main active substances (erucamide, oleamide, and camphor bromide) present in mulberry root exudates (MRE) against Ralstonia pseudosolanacearum (Rp), the causal agent of bacterial wilt. Our findings revealed that these three active substances inhibited the growth activity of Rp by affecting the cell morphology and extracellular polysaccharide content, as well as triggering a burst of reactive oxygen species. The active substances induced oxidative stress, leading to a decrease in Rp growth. Additionally, the expression levels of key genes in the hrp gene cluster (hrpB, hrpX, and hrpF) and other virulence-related genes (such as ripAW, ripAE, Rs5-4819, Rs5-4374, ace, egl3, and pehB) were significantly reduced upon treatment with the active substances. Further pathogenicity experiments demonstrated that root exudates (at a concentration of 1.5 mg·mL-1) delayed or slowed down the occurrence of bacterial wilt in mulberry. These findings provide valuable insight into the antimicrobial mechanisms of MRE against Rp and lay a theoretical foundation for the development and application of biocontrol agents to control mulberry bacterial wilt.

Whole-genome sequencing in medicinal plants: current progress and prospect

Advancements in genomics have dramatically accelerated the research on medicinal plants, and the development of herbgenomics has promoted the "Project of 1K Medicinal Plant Genome" to decipher their genetic code. However, it is difficult to obtain their high-quality whole genomes because of the prevalence of polyploidy and/or high genomic heterozygosity. Whole genomes of 123 medicinal plants were published until September 2022. These published genome sequences were investigated in this review, covering their classification, research teams, ploidy, medicinal functions, and sequencing strategies. More than 1,000 institutes or universities around the world and 50 countries are conducting research on medicinal plant genomes. Diploid species account for a majority of sequenced medicinal plants. The whole genomes of plants in the Poaceae family are the most studied. Almost 40% of the published papers studied species with tonifying, replenishing, and heat-cleaning medicinal effects. Medicinal plants are still in the process of domestication as compared with crops, thereby resulting in unclear genetic backgrounds and the lack of pure lines, thus making their genomes more difficult to complete. In addition, there is still no clear routine framework for a medicinal plant to obtain a high-quality whole genome. Herein, a clear and complete strategy has been originally proposed for creating a high-quality whole genome of medicinal plants. Moreover, whole genome-based biological studies of medicinal plants, including breeding and biosynthesis, were reviewed. We also advocate that a research platform of model medicinal plants should be established to promote the genomics research of medicinal plants.

Functional characterization of MaEXPA11 and its roles in response to biotic and abiotic stresses in mulberry

Mulberry is a traditional economic tree with various values in sericulture, ecology, food industry and medicine. Expansins (EXPs) are known as cell wall expansion related proteins and have been characterized to involve in plant development and responses to diverse stresses. In present study, twenty EXP and expansin-like (EXL) genes were identified in mulberry. RNA-seq results indicated that three EXP and EXL genes showed up-regulated expression level under sclerotiniose pathogen infection in three independent RNA-seq datasets. The most significant upregulated EXPA11 was selected as key EXP involving in response to sclerotiniose pathogen infection in mulberry. Furthermore, a comprehensive functional analysis was performed to reveal subcellular location, tissue expression profile of MaEXPA11 in mulberry. Down-regulation of MaEXPA11 using virus induced gene silence (VIGS) was performed to explore the function of MaEXPA11 in Morus alba. Results showed that MaEXPA11 can positively regulate mulberry resistance to Ciboria shiraiana infection and negatively regulate mulberry resistance to cold or drought stress.

Novel Tri-Segmented Rhabdoviruses: A Data Mining Expedition Unveils the Cryptic Diversity of Cytorhabdoviruses

Cytorhabdoviruses (genus Cytorhabdovirus, family Rhabdoviridae) are plant-infecting viruses with enveloped, bacilliform virions. Established members of the genus Cytorhabdovirus have unsegmented single-stranded negative-sense RNA genomes (ca. 10-16 kb) which encode four to ten proteins. Here, by exploring large publicly available metatranscriptomics datasets, we report the identification and genomic characterization of 93 novel viruses with genetic and evolutionary cues of cytorhabdoviruses. Strikingly, five unprecedented viruses with tri-segmented genomes were also identified. This finding represents the first tri-segmented viruses in the family Rhabdoviridae, and they should be classified in a novel genus within this family for which we suggest the name "Trirhavirus". Interestingly, the nucleocapsid and polymerase were the only typical rhabdoviral proteins encoded by those tri-segmented viruses, whereas in three of them, a protein similar to the emaravirus (family Fimoviridae) silencing suppressor was found, while the other predicted proteins had no matches in any sequence databases. Genetic distance and evolutionary insights suggest that all these novel viruses may represent members of novel species. Phylogenetic analyses, of both novel and previously classified plant rhabdoviruses, provide compelling support for the division of the genus Cytorhabdovirus into three distinct genera. This proposed reclassification not only enhances our understanding of the evolutionary dynamics within this group of plant rhabdoviruses but also illuminates the remarkable genomic diversity they encompass. This study not only represents a significant expansion of the genomics of cytorhabdoviruses that will enable future research on the evolutionary peculiarity of this genus but also shows the plasticity in the rhabdovirus genome organization with the discovery of tri-segmented members with a unique evolutionary trajectory.

Phylogenomic and syntenic data demonstrate complex evolutionary processes in early radiation of the rosids

Some of the most vexing problems of deep level relationship that remain in angiosperms involve the superrosids. The superrosid clade contains a quarter of all angiosperm species, with 18 orders in three subclades (Vitales, Saxifragales and core rosids) exhibiting remarkable morphological and ecological diversity. To help resolve deep-level relationships, we constructed a high-quality chromosome-level genome assembly for Tiarella polyphylla (Saxifragaceae) thus providing broader genomic representation of Saxifragales. Whole genome microsynteny analysis of superrosids showed that Saxifragales shared more synteny clusters with core rosids than Vitales, further supporting Saxifragales as more closely related with core rosids. To resolve the ordinal phylogeny of superrosids, we screened 122 single copy nuclear genes from genomes of 36 species, representing all 18 superrosid orders. Vitales were recovered as sister to all other superrosids (Saxifragales + core rosids). Our data suggest dramatic differences in relationships compared to earlier studies within core rosids. Fabids should be restricted to the nitrogen-fixing clade, while Picramniales, the Celastrales-Malpighiales (CM) clade, Huerteales, Oxalidales, Sapindales, Malvales and Brassicales formed an "expanded" malvid clade. The Celastrales-Oxalidales-Malpighiales (COM) clade (sensu APG IV) was not monophyletic. Crossosomatales, Geraniales, Myrtales and Zygophyllales did not belong to either of our well-supported malvids or fabids. There is strong discordance between nuclear and plastid phylogenetic hypotheses for superrosid relationships; we show that this is best explained by a combination of incomplete lineage sorting and ancient reticulation.

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.

S. A, H. B. M, Singh P, M. R. B, P S, Burji SM, T. M, V. S. Genome size, genetic diversity, and phenotypic variability imply the effect of genetic variation instead of ploidy on trait plasticity in the cross-pollinated tree species of mulberry

Elucidation of genome size (GS), genetic and phenotypic variation is the fundamental aspect of crop improvement programs. Mulberry is a cross-pollinated, highly heterozygous tree eudicot, and comprised of huge ploidy variation with great adaptability across the world. However, because of inadequate information on GS, ploidy-associated traits, as well as the correlation between genetic and phenotypic variation hinder the further improvement of mulberry. In this present research, a core set of 157 germplasm accessions belonging to eight accepted species of Morus including promising functional varieties were chosen to represent the genetic spectrum from the whole germplasm collection. To estimate the GS, accessions were subjected to flow cytometry (FCM) analysis and the result suggested that four different ploidies (2n = 2x, 3x, 4x, and 6x) with GS ranging from 0.72±0.005pg (S-30) to 2.89±0.015pg (M. serrata), accounting~4.01 fold difference. The predicted polyploidy was further confirmed with metaphase chromosome count. In addition, the genetic variation was estimated by selecting a representative morphologically, diverse population of 82 accessions comprised of all ploidy variations using simple sequence repeats (SSR). The estimated average Polymorphism Information Content (PIC) and expected heterozygosity showed high levels of genetic diversity. Additionally, three populations were identified by the model-based population structure (k = 3) with a moderate level of correlation between the populations and different species of mulberry, which imply the effect of genetic variation instead of ploidy on trait plasticity that could be a consequence of the high level of heterozygosity imposed by natural cross-pollination. Further, the correlation between ploidies, especially diploid and triploid with selected phenotypic traits was identified, however, consistency could not be defined with higher ploidy levels (>3x). Moreover, incite gained here can serve as a platform for future omics approaches to the improvement of mulberry traits.

Genomic Resequencing Unravels the Genetic Basis of Domestication, Expansion, and Trait Improvement in Morus Atropurpurea

Mulberry is an economically important plant in the sericulture industry and traditional medicine. However, the genetic and evolutionary history of mulberry remains largely unknown. Here, this work presents the chromosome-level genome assembly of Morus atropurpurea (M. atropurpurea), originating from south China. Population genomic analysis using 425 mulberry accessions reveal that cultivated mulberry is classified into two species, M. atropurpurea and M. alba, which may have originated from two different mulberry progenitors and have independent and parallel domestication in north and south China, respectively. Extensive gene flow is revealed between different mulberry populations, contributing to genetic diversity in modern hybrid cultivars. This work also identifies the genetic architecture of the flowering time and leaf size. In addition, the genomic structure and evolution of sex-determining regions are identified. This study significantly advances the understanding of the genetic basis and domestication history of mulberry in the north and south, and provides valuable molecular markers of desirable traits for mulberry breeding.

Genome size, chromosome number variation and its correlation with stomatal characters for assessment of ploidy levels in a core subset of mulberry (Morus spp.) germplasm

The large size of the germplasm collection along with scanty information on their cytological and genome constitution have hindered well-planned breeding schemes in mulberry. To address the issue, a study was undertaken to investigate the variability in DNA content and genome size, chromosome number, ploidy and its relation with important stomatal characteristics among 162 mulberry germplasm. These germplasm comprise a core subset of 150 collections along with a representative collection of different mulberry species including the wild. Among the germplasm belonging to 16 species, we identified 122 diploids (2n = 28), 4 aneuploids (2n = 30), 13 triploids (2n = 42), 15 tetraploids (2n = 56), 7 hexaploids (2n = 84) and 1 dodecosaploid (2n = 308) based on the chromosome count. Most of the cultivated mulberries are found to be diploids. The mean nuclear 2C DNA content estimated by Flow cytometry, varied from 0.723±0.006 pg (M. australis, 2n = 2x) to 7.732 pg (M. nigra, 2n = 22x). The 2C DNA content positively correlated with the ploidy status and stomatal length (r = 0.814, p<0.001). Based on the 1Cx value, the study also suggests that the majority of the polyploid species have experienced genome downsizing in relation to their diploid progenitors. This study provides the most essential information on chromosome number, ploidy and DNA content to facilitate the utilization of a core subset of germplasm in the mulberry breeding program.

The gap-free genome of mulberry elucidates the architecture and evolution of polycentric chromosomes

Mulberry is a fundamental component of the global sericulture industry, and its positive impact on our health and the environment cannot be overstated. However, the mulberry reference genomes reported previously remained unassembled or unplaced sequences. Here, we report the assembly and analysis of the telomere-to-telomere gap-free reference genome of the mulberry species, Morus notabilis, which has emerged as an important reference in mulberry gene function research and genetic improvement. The mulberry gap-free reference genome produced here provides an unprecedented opportunity for us to study the structure and function of centromeres. Our results revealed that all mulberry centromeric regions share conserved centromeric satellite repeats with different copies. Strikingly, we found that M. notabilis is a species with polycentric chromosomes and the only reported polycentric chromosome species up to now. We propose a compelling model that explains the formation mechanism of new centromeres and addresses the unsolved scientific question of the chromosome fusion-fission cycle in mulberry species. Our study sheds light on the functional genomics, chromosome evolution, and genetic improvement of mulberry species.

Genome sequencing and resequencing identified three horizontal gene transfers and uncovered the genetic mechanism on the intraspecies adaptive evolution of Gastrodia elata Blume

Horizontal gene transfer is a rare and useful genetic mechanism in higher plants. Gastrodia elata Blume (GE) (Orchidaceae), well known as traditional medicinal material in East Asia, adopts a heterotrophic lifestyle, thus being considered to be more prone to horizontal gene transfer (HGT). GE is a "polytypic species" that currently comprised of five recognized forms according to the plant morphology. G. elata Blume forma elata (GEE) and G. elata Bl.f.glauca (GEG) are two common forms that naturally grow in different habitats with difference in altitude and latitude. G. elata Bl.f.viridis (GEV) often occurs sporadically in cultivated populations of GEE and GEG. However, the genetic relationships and genetic mechanism underpinned the divergent ecological adaptations of GEE and GEG have not been revealed. Here, we assembled a chromosome-level draft genome of GEE with 1.04 Gb. Among predicted 17,895 protein coding genes, we identified three HGTs. Meanwhile, we resequenced 10 GEE accessions, nine GEG accessions, and 10 GEV accessions, and identified two independent genetic lineages: GEG_pedigree (GEG individuals and GEV individuals collected from GEG populations) and GEE_pedigree (GEE individuals and GEV individuals collected from GEE populations), which strongly support the taxonomic status of GEE and GEG as subspecies, not as different forms. In highly differentiated genomic regions of GEE_pedigree and GEG_pedigree, three chalcone synthase-encoding genes and one Phox/Bem1p (PB1) domain of encoding Auxin (AUX)/Indoleacetic acid (IAA) were identified in selection sweeping genome regions, which suggested that differentiation between GEE_pedigree and GEG_pedigree was promoted by the selection of genes related to photoresponse and growth and development. Overall, this new genome would be helpful for breeding and utilization of GE and the new findings would deepen the understanding about ecological adaptation and evolution of GE.

Exchanges of economic plants along the land silk road

BACKGROUD

The greatest contribution of the Silk Road is to communicate among different countries and nationalities, and promote two-way cultural exchanges between the East and the West. We now have clearer understanding about how material civilization and religious culture of Central Asia and West Asia spread eastward along the Land Silk Road. However, there is controversial about how crops migrate along the Land Silk Road.

RESULTS

We summarize archaeology, genetics, and genomics data to explore crop migration patterns. Of the 207 crops that were domesticated along the Land Silk Road, 19 for which genomic evidence was available were selected for discussion.

CONCLUSIONS

There were conflicting lines of evidence for the domestication of Tibetan barley, mustard, lettuce, buckwheat, and chickpea. The main reasons for the conflicting results may include incomplete early knowledge, record differences in different period, sample sizes, and data analysis techniques. There was strong evidence that Tibetan barley, barley, wheat, and jujube were introduced into China before the existence of the Land Silk Road; and mustard, lettuce, buckwheat, chickpea, alfalfa, walnut, cauliflower, grape, spinach, apple, cucumber, mulberry, and pea spread to China via trade and human migration along the Land Silk Road.

The Current Developments in Medicinal Plant Genomics Enabled the Diversification of Secondary Metabolites' Biosynthesis

Medicinal plants produce important substrates for their adaptation and defenses against environmental factors and, at the same time, are used for traditional medicine and industrial additives. Plants have relatively little in the way of secondary metabolites via biosynthesis. Recently, the whole-genome sequencing of medicinal plants and the identification of secondary metabolite production were revolutionized by the rapid development and cheap cost of sequencing technology. Advances in functional genomics, such as transcriptomics, proteomics, and metabolomics, pave the way for discoveries in secondary metabolites and related key genes. The multi-omics approaches can offer tremendous insight into the variety, distribution, and development of biosynthetic gene clusters (BGCs). Although many reviews have reported on the plant and medicinal plant genome, chemistry, and pharmacology, there is no review giving a comprehensive report about the medicinal plant genome and multi-omics approaches to study the biosynthesis pathway of secondary metabolites. Here, we introduce the medicinal plant genome and the application of multi-omics tools for identifying genes related to the biosynthesis pathway of secondary metabolites. Moreover, we explore comparative genomics and polyploidy for gene family analysis in medicinal plants. This study promotes medicinal plant genomics, which contributes to the biosynthesis and screening of plant substrates and plant-based drugs and prompts the research efficiency of traditional medicine.

Ectopic overexpression of mulberry MnT5H2 enhances melatonin production and salt tolerance in tobacco

Soil salinization severely inhibits plant growth and has become one of the major limiting factors for global agricultural production. Melatonin (N-acetyl-5-methoxytryptamine) plays an important role in regulating plant growth and development and in responding to abiotic stresses. Tryptamine-5-hydroxylase (T5H) is an enzyme essential for the biosynthesis of melatonin in plants. Previous studies have identified the gene MnT5H for melatonin synthesis in mulberry (Morus notabilis), but the role of this gene in response to salinity stress in mulberry is remain unclear. In this study, we ectopically overexpressed MnT5H2 in tobacco (Nicotiana tabacum L.) and treated it with NaCl solutions. Compared to wild-type (WT), melatonin content was significantly increased in the overexpression-MnT5H2 tobacco. Under salt stress, the expression of NtCAT, NtSOD, and NtERD10C and activity of catalase (CAT), peroxidase (POD), and the content of proline (Pro) in the transgenic lines were significantly higher than that in WT. The Malondialdehyde (MDA) content in transgenic tobacco was significantly lower than that of WT. Furthermore, transgenic tobacco seedlings exhibited faster growth in media with NaCl. This study reveals the changes of melatonin and related substance content in MnT5H2-overexpressing tobacco ultimately lead to improve the salt tolerance of transgenic tobacco, and also provides a new target gene for breeding plant resistance to salt.

Comparative and phylogenetic analyses of the chloroplast genome reveal the taxonomy of the Morus genus

Mulberry (genus Morus) is an economically important woody plant with an altered ploidy level. The variable number of Morus species recognized by different studies indicates that the genus is in need of revision. In this study, the chloroplast (CP) genomes of 123 Morus varieties were de novo assembled and systematically analyzed. The 123 varieties represented six Morus species, namely, Morus alba, Morus nigra, Morus notabilis, Morus rubra, Morus celtidifolia, and Morus serrata. The Morus CP genome was found to be 158,969~159,548 bp in size with 125 genes, including 81 protein coding, 36 tRNA, and 8 rRNA genes. The 87 out of 123 mulberry accessions were assigned to 14 diverse groups with identical CP genome, which indicated that they are maternally inherited and share 14 common ancestors. Then 50 diverse CP genomes occurred in 123 mulberry accessions for further study. The CP genomes of the Morus genus with a quadripartite structure have two inverted repeat (IR) regions (25,654~25,702 bp) dividing the circular genome into a large single-copy (LSC) region (87,873~88,243 bp) and small single-copy (SSC) region (19,740~19,994 bp). Analysis of the phylogenetic tree constructed using the complete CP genome sequences of Morus revealed a monophyletic genus and that M. alba consisted of two clades, M. alba var. alba and M. alba var. multicaulis. The Japanese cultivated germplasms were derived from M. alba var. multicaulis. We propose that the Morus genus be classified into six species, M. nigra, M. notabilis, M. serrata, M. celtidifolia, M. rubra, and M. alba with two subspecies, M. alba var. alba and M. alba var. multicaulis. Our findings provide a valuable resource for the classification, domestication, and breeding improvement of mulberry.

The complete chloroplast genome sequences of three Broussonetia species and comparative analysis within the Moraceae

Background

Species of Broussonetia (family Moraceae) are commonly used to make textiles and high-grade paper. The distribution of Broussonetia papyrifera L. is considered to be related to the spread and location of humans. The complete chloroplast (cp) genomes of B. papyrifera, Broussonetia kazinoki Sieb., and Broussonetia kaempferi Sieb. were analyzed to better understand the status and evolutionary biology of the genus Broussonetia.

Methods

The cp genomes were assembled and characterized using SOAPdenovo2 and DOGMA. Phylogenetic and molecular dating analysis were performed using the concatenated nucleotide sequences of 35 species in the Moraceae family and were based on 66 protein-coding genes (PCGs). An analysis of the sequence divergence (pi) of each PCG among the 35 cp genomes was conducted using DnaSP v6. Codon usage indices were calculated using the CodonW program.

Results

All three cp genomes had the typical land plant quadripartite structure, ranging in size from 160,239 bp to 160,841 bp. The ribosomal protein L22 gene (RPL22) was either incomplete or missing in all three Broussonetia species. Phylogenetic analysis revealed two clades. Clade 1 included Morus and Artocarpus, whereas clade 2 included the other seven genera. Malaisia scandens Lour. was clustered within the genus Broussonetia. The differentiation of Broussonetia was estimated to have taken place 26 million years ago. The PCGs' pi values ranged from 0.0005 to 0.0419, indicating small differences within the Moraceae family. The distribution of most of the genes in the effective number of codons plot (ENc-plot) fell on or near the trend line; the slopes of the trend line of neutrality plots were within the range of 0.0363-0.171. These results will facilitate the identification, taxonomy, and utilization of the Broussonetia species and further the evolutionary studies of the Moraceae family.

Systematic functional characterization of cinnamyl alcohol dehydrogenase family members revealed their functional divergence in lignin biosynthesis and stress responses in mulberry

Mulberry (Morus) is used as a feed additive and biofuel materials. Cinnamyl alcohol dehydrogenase (CAD; EC 1.1.1.95) catalyzes the final step of monolignol biosynthesis and is responsible for various monolignols. Five MaCADs from Morus alba were cloned and functionally characterized in the present study. These MaCADs encoded proteins with 357-364 amino acids, and the putative protein sequences conservatively possessed two Zn²⁺ binding motifs and an NADP(H) cofactor binding motif. However, MaCAD1, 2, and 5 shared similar amino acids at substrate binding positions that differed from those possessed by bona fide CADs. MaCAD3 and 4 had conservative substrate binding sites, and both phylogenetic and expression profile analysis indicated they were bona fide CADs involved in lignin biosynthesis. The enzymatic assay showed that MaCAD1 and 5 had a high affinity to p-coumaryl aldehyde. MaCAD4 preferentially used coniferyl aldehyde and sinapyl aldehyde as substrates. His-72 and Tyr-124 in MaCAD1 stabilized p-coumaryl aldehyde, and may have resulted in the substrate preference for p-coumaryl aldehyde. Down-regulation of MaCADs in mulberry showed that MaCAD3/4 were dominant CADs that functioned in monolignol biosynthesis, and decreased MaCAD3/4 resulted in significant decreases of lignin content in both stems and leaves. MaCADs exhibited different expression patterns in response to various stresses, indicating their possible diverse roles. MaCAD2 and MaCAD5 may play positive roles in response to drought and cold stresses, respectively. These results provide a systematic functional analysis of MaCADs in mulberry and an important foundation for the genetic modification of the monolignol pathway in mulberry.

Multi-omics analyses reveal MdMYB10 hypermethylation being responsible for a bud sport of apple fruit color

Apple bud sports offer a rich resource for clonal selection of numerous elite cultivars. The accumulation of somatic mutations as plants develop may potentially impact the emergence of bud sports. Previous studies focused on somatic mutation in the essential genes associated with bud sports. However, the rate and function of genome-wide somatic mutations that accumulate when a bud sport arises remain unclear. In this study, we identified a branch from a 10-year-old tree of the apple cultivar 'Oregon Spur II' as a bud sport. The mutant branch showed reduced red coloration on fruit skin. Using this plant material, we assembled a high-quality haplotype reference genome consisting of 649.61 Mb sequences with a contig N50 value of 2.04 Mb. We then estimated the somatic mutation rate of the apple tree to be 4.56 × 10 -8 per base per year, and further identified 253 somatic single-nucleotide polymorphisms (SNPs), including five non-synonymous SNPs, between the original type and mutant samples. Transcriptome analyses showed that 69 differentially expressed genes between the original type and mutant fruit skin were highly correlated with anthocyanin content. DNA methylation in the promoter of five anthocyanin-associated genes was increased in the mutant compared with the original type as determined using DNA methylation profiling. Among the genetic and epigenetic factors that directly and indirectly influence anthocyanin content in the mutant apple fruit skin, the hypermethylated promoter of MdMYB10 is important. This study indicated that numerous somatic mutations accumulated at the emergence of a bud sport from a genome-wide perspective, some of which contribute to the low coloration of the bud sport.

Screening and Verification of Photosynthesis and Chloroplast-Related Genes in Mulberry by Comparative RNA-Seq and Virus-Induced Gene Silencing

Photosynthesis is one of the most important factors in mulberry growth and production. To study the photosynthetic regulatory network of mulberry we sequenced the transcriptomes of two high-yielding (E1 and E2) and one low-yielding (H32) mulberry genotypes at two-time points (10:00 and 12:00). Re-annotation of the mulberry genome based on the transcriptome sequencing data identified 22,664 high-quality protein-coding genes with a BUSCO-assessed completeness of 93.4%. A total of 6587 differentially expressed genes (DEGs) were obtained in the transcriptome analysis. Functional annotation and enrichment revealed 142 out of 6587 genes involved in the photosynthetic pathway and chloroplast development. Moreover, 3 out of 142 genes were further examined using the VIGS technique; the leaves of MaCLA1- and MaTHIC-silenced plants were markedly yellowed or even white, and the leaves of MaPKP2-silenced plants showed a wrinkled appearance. The expression levels of the ensiled plants were reduced, and the levels of chlorophyll b and total chlorophyll were lower than those of the control plants. Co-expression analysis showed that MaCLA1 was co-expressed with CHUP1 and YSL3; MaTHIC was co-expressed with MaHSP70, MaFLN1, and MaEMB2794; MaPKP2 was mainly co-expressed with GH9B7, GH3.1, and EDA9. Protein interaction network prediction revealed that MaCLA1 was associated with RPE, TRA2, GPS1, and DXR proteins; MaTHIC was associated with TH1, PUR5, BIO2, and THI1; MaPKP2 was associated with ENOC, LOS2, and PGI1. This study offers a useful resource for further investigation of the molecular mechanisms involved in mulberry photosynthesis and preliminary insight into the regulatory network of photosynthesis.

Chromosome-level Genomes Reveal the Genetic Basis of Descending Dysploidy and Sex Determination in Morus Plants

Multiple plant lineages have independently evolved sex chromosomes and variable karyotypes to maintain their sessile lifestyles through constant biological innovation. Morus notabilis, a dioecious mulberry species, has the fewest chromosomes among Morus spp., but the genetic basis of sex determination and karyotype evolution in this species has not been identified. In this study, three high-quality genome assemblies were generated for Morus spp. [including dioecious M. notabilis (male and female) and Morus yunnanensis (female)] with genome sizes of 301-329 Mb and were grouped into six pseudochromosomes. Using a combination of genomic approaches, we found that the putative ancestral karyotype of Morus species was close to 14 protochromosomes, and that several chromosome fusion events resulted in descending dysploidy (2n = 2x = 12). We also characterized a ∼ 6.2-Mb sex-determining region on chromosome 3. Four potential male-specific genes, a partially duplicatedDNA helicase gene (named MSDH) and three Ty3_Gypsy long terminal repeat retrotransposons (named MSTG1/2/3), were identified in the Y-linked area and considered to be strong candidate genes for sex determination or differentiation. Population genomic analysis showed that Guangdong accessions in China were genetically similar to Japanese accessions of mulberry. In addition, genomic areas containing selective sweeps that distinguish domesticated mulberry from wild populations in terms of flowering and disease resistance were identified. Our findings provide an important genetic resource for sex identification research and molecular breeding in mulberry.

First Report of Colletotrichum aenigma Causing Anthracnose on Mulberry Leaves in China

Mulberry (Morus alba L.) has been grown worldwide as a crop for silkworm rearing for over five thousand years (Jiao et al. 2020). In July 2021, a leaf spot disease was observed on mulberry leaves in Wuhan city (114°33'E, 30°48'N), Hubei province, China, with approximately 40% of leaves (about 300 trees) affected. Early symptoms were light brown, with small lesions subsequently expanding to larger sometimes irregular dark brown or black spots surrounded by yellow-brown margins, with easily perforated necrotic lesions. Leaf tissues (5 mm×5 mm) were excised from the border between diseased and healthy tissues, surface sterilized with 75% ethanol solution for 30 s and 2.5% sodium hypochlorite for 2 min, washed thrice in sterile distilled water, and then placed on potato dextrose agar (PDA), and incubated at 25°C in darkness. Four isolates (C1, C9, CHS2, and CHS6) were subcultured using the single-spore method. On PDA, colonies were cottony, pale white from above, and white to grayish-green on the reverse side. Conidia were aseptate, hyaline, subcylindrical with broadly rounded ends, 8.4 to 18.3×4.1 to 7.7 μm (mean = 13.9×5.5 μm, n = 30). Appressoria were typically elliptic or irregular with a few lobes, dark brown, 5.9 to 9.6×4.2 to 8.1 μm (mean = 7.9 ×5.7 μm, n = 30). The morphological characteristics of the isolates matched the descriptions of Colletotrichum gloeosporioides species complex (Weir et al. 2012). The isolates were further identified by analysis of the ribosomal internal transcribed spacers (ITS), glyceraldehyde-3-phosphate dehydrogenase (GAPDH), calmodulin (CAL), actin (ACT), chitin synthase (CHS-1), glutamine synthetase (GS), and β-tubulin 2 (TUB2) genes, amplified respectively with ITS1/ITS4, GDF/GDR, CL1C/CL2C, ACT-512F/ACT-783R, CHS-79F/CHS-345R, GSF1/GSR, and Bt2a/Bt2b (Glass and Donaldson 1995; Weir et al. 2012; White et al. 1990). The sequences were deposited in GenBank (ON492187-ON492214). Concatenated sequences of the seven genes in addition to Colletotrichum species sequences from GenBank were used to conduct a phylogenetic analysis using Maximum-Likelihood (ML) method in MEGA7. The four isolates were grouped into a clade with Colletotrichum aenigma supported by a high bootstrap value (89%), and hence, they were identified as C. aenigma based on the morphological and molecular analyses. To confirm Koch's postulates, wounded leaves of six healthy 2-month-old seedlings made by a sterile needle were inoculated with each isolate by spraying 10 ml of conidial suspensions (105 conidia/ml) on each plant, and the control plants were treated with sterile distilled water. All the treated plants were kept in a plastic box containing sterile water and incubated at 28°C in a 12 h/12 h light/dark cycle. The test was performed three times. After 7 days, typical anthracnose lesions appeared on all inoculated leaves, whereas control plants remained asymptotic. Furthermore, C. aenigma was only reisolated from the symptomatic leaves. Previous studies reported five Colletotrichum species (C. morifolium, C. fioriniae, C. brevisporum, C. karstii, and C. kahawae subsp. ciggaro) to cause this disease on mulberry in China (Tian, 1981; Xue et al. 2019). To our knowledge, this is the first report of C. aenigma causing anthracnose on mulberry in China. The finding will facilitate epidemiological studies and the development of effective control strategies for the disease.

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.

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.

Systematic Identification of Bioactive Compositions in Leaves of Morus Cultivars Using UHPLC-ESI-QTOF-MS/MS and Comprehensive Screening of High-Quality Resources

Morus spp. leaves (MSLs) show various beneficial effects in the treatment of metabolic-related diseases, which have created a growing interest in MSL development as dietary supplements and functional foods. The illustration of chemical compositions and screening of high-quality MSL resources are therefore necessary for further application. This study developed a new UHPLC-ESI-QTOF-MS/MS strategy of in-source collision-induced dissociation (IS-CID) and target collision-cell CID (TCC-CID) to quickly capture analogues with consistent skeleton, and combined global natural product social molecular networking (GNPS) to efficiently annotate bioactive phytochemicals in MSLs. For the results, 49 bioactive ingredients, including quercetin-type flavonoids, kaempferol-type flavonoids, chlorogenic acid isomers, 1-deoxynojirimycin, γ-aminobutyric acid, amino acids, and unsaturated fatty acids, were systematically identified in MSLs for the first time. Quantification for the typical components was simultaneously carried out in MSLs of 90 Morus resources collected from different locations. Partial least squares discriminant analysis (PLS-DA) indicated that quercetin-3-O-(6″-O-malonyl)-glucoside, rutin, kaempferol-3-O-(6″-O-malonyl)-glucoside, kaempferol-3-O-rutinoside, and chlorogenic acid showed high variable importance in the project (VIP > 1) that were significant constituents for the differences between MSL species. Then, high-quality MSLs were comprehensively screened in multiple Morus cultivars based on the criteria importance through intercriteria correlation (CRITIC) method. This study presented an efficient strategy to annotate bioactive compounds, revealed the difference of bioactive components in MSLs, and provided important information for the high-value production of Morus cultivars in food and supplement fields.

The Mulberry SPL Gene Family and the Response of MnSPL7 to Silkworm Herbivory through Activating the Transcription of MnTT2L2 in the Catechin Biosynthesis Pathway

SQUAMOSA PROMOTER BINDING PROTEIN-LIKE (SPL) genes, as unique plant transcription factors, play important roles in plant developmental regulation and stress response adaptation. Although mulberry is a commercially valuable tree species, there have been few systematic studies on SPL genes. In this work, we identified 15 full-length SPL genes in the mulberry genome, which were distributed on 4 Morus notabilis chromosomes. Phylogenetic analysis clustered the SPL genes from five plants (Malus × domestica Borkh, Populus trichocarpa, M. notabilis, Arabidopsis thaliana, and Oryza sativa) into five groups. Two zinc fingers (Zn1 and Zn2) were found in the conserved SBP domain in all of the MnSPLs. Comparative analyses of gene structures and conserved motifs revealed the conservation of MnSPLs within a group, whereas there were significant structure differences among groups. Gene quantitative analysis showed that the expression of MnSPLs had tissue specificity, and MnSPLs had much higher expression levels in older mulberry leaves. Furthermore, transcriptome data showed that the expression levels of MnSPL7 and MnSPL14 were significantly increased under silkworm herbivory. Molecular experiments revealed that MnSPL7 responded to herbivory treatment through promoting the transcription of MnTT2L2 and further upregulating the expression levels of catechin synthesis genes (F3′H, DFR, and LAR).

Chromosome restructuring and number change during the evolution of Morus notabilis and Morus alba

Mulberry (Morus spp.) is an economically important plant as the main food plant used for rearing domesticated silkworm and it has multiple uses in traditional Chinese medicine. Two basic chromosome numbers (Morus notabilis, n = 7, and Morus alba, n = 14) have been reported in the genus Morus, but the evolutionary history and relationship between them remain unclear. In the present study, a 335-Mb high-quality chromosome-scale genome was assembled for the wild mulberry species M. notabilis. Comparative genomic analyses indicated high chromosomal synteny between the 14 chromosomes of cultivated M. alba and the six chromosomes of wild M. notabilis. These results were successfully verified by fluorescence in situ hybridization. Chromosomal fission/fusion events played crucial roles in the chromosome restructuring process between M. notabilis and M. alba. The activity of the centromere was another key factor that ensured the stable inheritance of chromosomes. Our results also revealed that long terminal repeat retrotransposons were a major driver of the genome divergence and evolution of the mulberry genomes after they diverged from each other. This study provides important insights and a solid foundation for studying the evolution of mulberry, allowing the accelerated genetic improvement of cultivated mulberry species.

AGAMOUS Gene as a New Sex-Identification Marker in Fig (Ficus carica L.) Is More Efficient Than RAN1

Fig is an ancient gynodioecious fruit tree with females for commercial fruit production and hermaphrodites (males) sometimes used as pollen providers. An early sex-identification method would improve breeding efficiency. Three AGAMOUS (AG) genes were recruited from the Ficus carica genome using AG sequences from Ficus microcarpa and Ficus hispida. FcAG was 5230 bp in length, with 7 exons and 6 introns, and a 744-bp coding sequence. The gene was present in both female and male fig genomes, with a 15-bp deletion in the 7th exon. The other two AG genes (FcAG2-Gall_Stamen and FcAG3-Gall_Stamen) were male-specific, without the 15-bp deletion (759-bp coding sequence), and were only expressed in the gall and stamen of the male fig fruit. Using the deletion as the forward primer (AG-Marker), male plants were very efficiently identified by the presence of a 146-bp PCR product. The previously reported fig male and female polymorphism gene RESPONSIVE-TO-ANTAGONIST1 (RAN1) was also cloned and compared between male and female plants. Fifteen SNPs were found in the 3015-bp protein-coding sequence. Among them, 12 SNPs were identified as having sex-differentiating capacity by checking the sequences of 27 known male and 24 known female cultivars. A RAN1-Marker of 608 bp, including 6 SNPs, was designed, and a PCR and sequencing-based method was verified with 352 fig seedlings from two hybrid populations. Our results confirmed that the newly established AG-Marker is as accurate as the RAN1-Marker, and provide new clues to understanding Ficus sex determination.

An APSES Transcription Factor Xbp1 Is Required for Sclerotial Development, Appressoria Formation, and Pathogenicity in Ciboria shiraiana

Sclerotinia diseases are important plant fungal diseases that, causes huge economic worldwide losses every year. Ciboria shiraiana is the main pathogen that results in mulberry sclerotia diseases. Sclerotia and appressoria play important roles in long-term pathogen survival and in host infection during life and disease cycles. However, the molecular mechanisms of sclerotial development and appressoria formation in C. shiraiana have not been well studied. Here, an Asm1p, Phd1p, Sok2p, Efg1p and StuAp (APSES)-type transcription factor in C. shiraiana, CsXbp1, involved in sclerotial development and appressoria formation was functionally characterized. Bioinformatics analyses showed that CsXbp1 contained an APSES-type DNA binding domain. The expression levels of CsXbp1 were higher in sclerotia and during later stages of infection. Compared with wild-type strains, hyphal growth was slower, the number and weight of sclerotia were reduced significantly, and appressoria formation was obviously delayed in CsXbp1 RNA interference (RNAi) strains. Moreover, the CsXbp1 RNAi strains showed weakened pathogenicity owing to compound appressoria defects. Tobacco rattle virus-mediated host-induced gene silencing enabled Nicotiana benthamiana to increase its resistance to C. shiraiana by reducing the CsXbp1 transcripts level. Thus, CsXbp1 plays vital roles in sclerotial formation, appressoria formation, and pathogenicity in C. shiraiana. This study provides new insights into the infection mechanisms of C. shiraiana and plant resistance breeding.

Secondary Metabolism and Hormone Response Reveal the Molecular Mechanism of Triploid Mulberry (Morus Alba L.) Trees Against Drought

The improvement of a plant's tolerance to drought is a major endeavor in agriculture. Polyploid plants often exhibit enhanced stress tolerance relative to their diploid progenitor, but the matching stress tolerance is still little understood. Own-rooted stem cuttings of mulberry (Morus alba L.) cultivar Shinichinose (2n = 2x = 28) and Shaansang-305 (2n = 3x = 42) were used in this study, of which the latter (triploid) has more production and application purposes. The responses of triploid Shaansang-305 and diploid progenitor ShinIchinose under drought stress were compared through an investigation of their physiological traits, RNA-seq, and secondary metabolome analysis. The results showed that the triploid exhibited an augmented abscisic acid (ABA) content and a better stress tolerance phenotype under severe drought stress. Further, in the triploid plant some genes (TSPO, NCED3, and LOC21398866) and ATG gene related to ABA signaling showed significantly upregulated expression. Interestingly, the triploid accumulated higher levels of RWC and SOD activity, as well as more wax on the leaf surface, but with less reductive flavonoid than in diploid. Our results suggest triploid plants may better adapt to with drought events. Furthermore, the flavonoid metabolism involved in drought resistance identified here may be of great value to medicinal usage of mulberry. The findings presented here could have substantial implications for future studies of crop breeding.

Comparative transcriptome analysis provides insights into steviol glycoside synthesis in stevia (Stevia rebaudiana Bertoni) leaves under nitrogen deficiency

Transcriptome analysis revealed the potential mechanism of nitrogen regulating steviol glycosides synthesis via shifting of leaf carbon metabolic flux or inducing certain transcription factors. Nitrogen (N) plays key regulatory roles in both stevia (Stevia rebaudiana) growth and the synthesis of its functional metabolite steviol glycosides (SGs), but the mechanism by which this nutrient regulates SGs synthesis remains to be elucidated. To address this question, a pot experiment was performed in a greenhouse where stevia plants fertilized with N (the control as CK plants) and compared with plants without the supply of N. Physiological and biochemical analyses were conducted to test the growth and metabolic responses of plants to N regimes. Our results showed that N deficiency significantly inhibited plant growth and leaf photosynthesis, while increased leaf SGs contents in stevia (49.97, 46.64 and 84.80% respectively for rebaudioside A, stevioside, and rebaudioside C), which may be partly due to "concentration effect". Then, transcriptome analysis was conducted to understand the underlying mechanisms. A total of 535 differentially expressed genes were identified, and carbon metabolism-related events were highlighted by Gene Ontology and Kyoto Encyclopedia of Genes and Genomes. Many of these genes were significantly upregulated by N-deficiency, including those involved in "phenylpropanoid biosynthesis", "flavonoid biosynthesis" and "starch and sucrose metabolism". Our study also analyzed the expression patterns of SGs synthesis-related genes under two N regimes and the potential transcription factors linking N nutrition and SG metabolism. N-deficiency may promote SGs synthesis by changing the carbon metabolism flux or inducing certain transcription factors. Our results provide deeper insight into the relationship between N nutrition and SGs synthesis in stevia plants.

Genome-wide analysis of the lignin toolbox for morus and the roles of lignin related genes in response to zinc stress

Mulberry (Morus, Moraceae) is an important economic plant with nutritional, medicinal, and ecological values. Lignin in mulberry can affect the quality of forage and the saccharification efficiency of mulberry twigs. The availability of the Morus notabilis genome makes it possible to perform a systematic analysis of the genes encoding the 11 protein families specific to the lignin branch of the phenylpropanoid pathway, providing the core genes for the lignin toolbox in mulberry. We performed genome-wide screening, which was combined with de novo transcriptome data for Morus notabilis and Morus alba variety Fengchi, to identify putative members of the lignin gene families followed by phylogenetic and expression profile analyses. We focused on bona fide clade genes and their response to zinc stress were further distinguished based on expression profiles using RNA-seq and RT-qPCR. We finally identified 31 bona fide genes in Morus notabilis and 25 bona fide genes in Fengchi. The putative function of these bona fide genes was proposed, and a lignin toolbox that comprised 19 genes in mulberry was provided, which will be convenient for researchers to explore and modify the monolignol biosynthesis pathway in mulberry. We also observed changes in the expression of some of these lignin biosynthetic genes in response to stress caused by excess zinc in Fengchi and proposed that the enhanced lignin biosynthesis in lignified organs and inhibition of lignin biosynthesis in leaf is an important response to zinc stress in mulberry.

Genomic analyses provide comprehensive insights into the domestication of bast fiber crop ramie (Boehmeria nivea)

Ramie (Boehmeria nivea) is an economically important natural fiber-producing crop that has been cultivated for thousands of years in China; however, the evolution of this crop remains largely unknown. Here, we report a ramie domestication analysis based on genome assembly and resequencing of cultivated and wild accessions. Two chromosome-level genomes representing wild and cultivated ramie were assembled de novo. Numerous structural variations between two assemblies, together with the genetic variations from population resequencing, constituted a comprehensive genomic variation map for ramie. Domestication analysis identified 71 high-confidence selective sweeps comprising 320 predicted genes, and 29 genes from sweeps were associated with fiber growth in the expression. In addition, we identified seven genetic loci associated with the fiber yield trait in the segregated population derived from the crossing of two assembled accessions, and two of which showed an overlap with the selective sweeps. These findings indicated that bast fiber traits were focused on during the domestication history of ramie. This study sheds light on the domestication of ramie and provides a valuable resource for biological and breeding studies of this important crop.

Phytoremediation potential of forage mulberry (Morus atropurpurea Roxb.) for cadmium contaminated paddy soils

Mulberry is an important material to utilize the Cd polluted farmland in China and planting forage mulberry is a new development direction. This study aimed to investigate the changes of annual biomass and Cd content in shoot of Guisangyou 62, Guisangyou 12 and Yuesang 11 in field XT-C1, XT-X1 and ZZ-M1 under the pressure of Cd in 3 years. The Cd extraction ability of forage mulberry was analyzed, and the safety of forage mulberry was also discussed. The results showed that the annual biomass of each forage mulberry shoot could reach 64.52 ∼ 86.61 t/hectare (ha). The total harvest biomass of Guisangyou 12 was the highest, followed by Guisangyou 62 and Yuesang 11. In the same test area, for different forage mulberry varieties, there were no significant differences in Cd content in the shoot at each sampling time, and the Cd concentrations in shoot were in the range of 0.05 ∼ 0.66 mg/kg, meeting the hygienical standard for feeds (GB 13078-2017, China). Without considering the test area, the average Cd removal amount of each forage mulberry in a year was about 18.52 g/ha. Planting forage mulberry may become a new ecological economic model to achieve the safe utilization of Cd polluted farmland.Novelty statement Mulberry is one of the most important plants for safe utilization the Cd polluted farmland in China and planting mulberry as animal feed is a new development direction. This study investigated the effects of Cd on the shoot biomass of 3 forage mulberry varieties at 3 experimental areas in 3 years. It also examined how much Cd could be removed from soil by harvesting forage mulberry for 4 times a year. The Cd content in the shoot of forage mulberry and its safety were also evaluated.

Comparison of De Novo Assembly Strategies for Bacterial Genomes

(1) Background: Short-read sequencing allows for the rapid and accurate analysis of the whole bacterial genome but does not usually enable complete genome assembly. Long-read sequencing greatly assists with the resolution of complex bacterial genomes, particularly when combined with short-read Illumina data. However, it is not clear how different assembly strategies affect genomic accuracy, completeness, and protein prediction. (2) Methods: we compare different assembly strategies for Haemophilus parasuis, which causes Glässer's disease, characterized by fibrinous polyserositis and arthritis, in swine by using Illumina sequencing and long reads from the sequencing platforms of either Oxford Nanopore Technologies (ONT) or SMRT Pacific Biosciences (PacBio). (3) Results: Assembly with either PacBio or ONT reads, followed by polishing with Illumina reads, facilitated high-quality genome reconstruction and was superior to the long-read-only assembly and hybrid-assembly strategies when evaluated in terms of accuracy and completeness. An equally excellent method was correction with Homopolish after the ONT-only assembly, which had the advantage of avoiding hybrid sequencing with Illumina. Furthermore, by aligning transcripts to assembled genomes and their predicted CDSs, the sequencing errors of the ONT assembly were mainly indels that were generated when homopolymer regions were sequenced, thus critically affecting protein prediction. Polishing can fill indels and correct mistakes. (4) Conclusions: The assembly of bacterial genomes can be directly achieved by using long-read sequencing techniques. To maximize assembly accuracy, it is essential to polish the assembly with homologous sequences of related genomes or sequencing data from short-read technology.

Hierarchical Action of Mulberry miR156 in the Vegetative Phase Transition

The vegetative phase transition is a prerequisite for flowering in angiosperm plants. Mulberry miR156 has been confirmed to be a crucial factor in the vegetative phase transition in Arabidopsis thaliana. The over-expression of miR156 in transgenic Populus × canadensis dramatically prolongs the juvenile phase. Here, we find that the expression of mno-miR156 decreases with age in all tissues in mulberry, which led us to study the hierarchical action of miR156 in mulberry. Utilizing degradome sequencing and dual-luciferase reporter assays, nine MnSPLs were shown to be directly regulated by miR156. The results of yeast one-hybrid and dual-luciferase reporter assays also revealed that six MnSPLs could recognize the promoter sequences of mno-miR172 and activate its expression. Our results demonstrate that mno-miR156 performs its role by repressing MnSPL/mno-miR172 pathway expression in mulberry. This work uncovered a miR156/SPLs/miR172 regulation pathway in the development of mulberry and fills a gap in our knowledge about the molecular mechanism of vegetative phase transition in perennial woody plants.

Heterotrimeric G-protein γ subunits regulate ABA signaling in response to drought through interacting with PP2Cs and SnRK2s in mulberry (Morus alba L.)

ABA signaling plays a central role in regulating plants respond to drought. Although much progress has been made in understanding the functions of ABA signaling in drought response, very little information is available regarding woody plants. In this study, the components of ABA signaling pathway were identified in mulberry which has excellent adaptation to drought, including three PYLs, two PP2Cs, two SnRK2s, four ABFs, and an ABA responsive gene MaRD29B. The gene expression of ABA signaling components exhibited significant response to ABA and drought, and their roles in drought response were revealed using a transient transformation system in mulberry seedlings. We discovered the ABA signaling components, MaABI1/2 and MaSnRK2.1/2.4, could directly interact with G-protein γ subunits, MaGγ1 and MaGγ2, which indicated that G-protein γ subunits may mediate the signal crosstalk between G-proteins and ABA signaling. Transient activation assay in tobacco and RNAi silencing assay in mulberry further demonstrated that MaGγ1 and MaGγ2 regulated drought response by enhancing ABA signaling. This study expands the repertoire of ABA signaling controlling drought responses in plants and provides the direct evidence about the crosstalk between ABA signaling and G-proteins for the first time.

Functional characterization of two chalcone isomerase (CHI) revealing their responsibility for anthocyanins accumulation in mulberry

Mulberry (Morus sp., Moraceae) is an important economic crop plant and mulberry fruits are rich in anthocyanidins. Chalcone isomerase (CHI) catalyzes the conversion of chalcones to flavanones providing precursors for biosynthesis of anthocyanidins. In this study, bona fide CHIs were cloned and characterized from different Morus species with differently colored fruits (Morus multicaulis, Mm and Morus alba variety LvShenZi, LSZ). Enzymatic assay of MmCHI1 and MmCHI2 showed that they can utilize naringenin chalcone as substrate. The catalytic efficiency of MmCHI2 and LSZCHI2 are approximately 200 and 120-fold greater than that of MmCHI1 respectively. Phylogenetic analysis showed the two mulberry CHIs belonged to different sub-clade of Type I CHI1 named type IA (CHI2) and type IB (CHI1). Type IB CHIs are mulberry specific. MmCHI1 and MmCHI2 had similar expression profiles and showed preferred expression in fruits. In addition, both mulberry CHI1 and CHI2 played roles in the response to excess zinc stress and sclerotiniose pathogen infection. Both MmCHI1 and MmCHI2 expression levels showed positive close relationship with anthocyanins content during fruit ripening process. The co-expression of MmCHI1 and MmCHI2 was observed during fruit ripening process and in transgenic mulberry. VIGS (virus induced gene silence) targeting on MmCHI1 and MmCHI2 showed significant down-regulation of MmCHI2 instead of MmCHI1 would result in significant (about 50%) decrease in anthocyanins content. MmCHI2 is the dominant CHI for anthocyanins accumulation in mulberry. The results presented in this work provided insight on bona fide CHIs in mulberry and reveal their roles in anthocyanins accumulation.

The Diversity of Plant Sex Chromosomes Highlighted through Advances in Genome Sequencing

For centuries, scientists have been intrigued by the origin of dioecy in plants, characterizing sex-specific development, uncovering cytological differences between the sexes, and developing theoretical models. Through the invention and continued improvements in genomic technologies, we have truly begun to unlock the genetic basis of dioecy in many species. Here we broadly review the advances in research on dioecy and sex chromosomes. We start by first discussing the early works that built the foundation for current studies and the advances in genome sequencing that have facilitated more-recent findings. We next discuss the analyses of sex chromosomes and sex-determination genes uncovered by genome sequencing. We synthesize these results to find some patterns are emerging, such as the role of duplications, the involvement of hormones in sex-determination, and support for the two-locus model for the origin of dioecy. Though across systems, there are also many novel insights into how sex chromosomes evolve, including different sex-determining genes and routes to suppressed recombination. We propose the future of research in plant sex chromosomes should involve interdisciplinary approaches, combining cutting-edge technologies with the classics to unravel the patterns that can be found across the hundreds of independent origins.

Molecular Mechanisms Underlying the Biosynthesis of Melatonin and Its Isomer in Mulberry

Mulberry (Morus alba L.) leaves and fruit are traditional Chinese medicinal materials with anti-inflammatory, immune regulatory, antiviral and anti-diabetic properties. Melatonin performs important roles in the regulation of circadian rhythms and immune activities. We detected, identified and quantitatively analyzed the melatonin contents in leaves and mature fruit from different mulberry varieties. Melatonin and three novel isoforms were found in the Morus plants. Therefore, we conducted an expression analysis of melatonin and its isomer biosynthetic genes and in vitro enzymatic synthesis of melatonin and its isomer to clarify their biosynthetic pathway in mulberry leaves. MaASMT4 and MaASMT20, belonging to class II of the ASMT gene family, were expressed selectively in mulberry leaves, and two recombinant proteins that they expressed catalyzed the conversion of N-acetylserotonin to melatonin and one of three isomers in vitro. Unlike the ASMTs of Arabidopsis and rice, members of the three ASMT gene families in mulberry can catalyze the conversion of N-acetylserotonin to melatonin. This study provides new insights into the molecular mechanisms underlying melatonin and its isomers biosynthesis and expands our knowledge of melatonin isomer biosynthesis.

Mulberry leaf attenuates atherosclerotic lesions in patients with coronary heart disease possibly via 1-Deoxynojirimycin: A placebo-controlled, double-blind clinical trial

Mulberry leaf tea reduces atherosclerotic lesions and its main component is Deoxynojirimycin (DNJ). We aimed to explore the effects of Mulberry leaf DNJ (MLD) on atherosclerotic lesions in the patients with coronary heart disease (CHD). MLD and serum DNJ was detected by high-performance liquid chromatography (HPLC). The CHD patients with low-density lipoprotein cholesterol >140 mg/dl were assigned into the EG (MLD treatment) and the CG (placebo treatment) groups. Serum biochemical indices, the cerebrovascular and cardiovascular events and carotid intima-media thickness (IMT) were measured before and after a 1-year intervention. Pearson Correlation Coefficient test was used to explore the relationship between serum levels of DNJ and IMT values. Serum DNJ was 70 ± 50 ng/ml in the EG group but no serum DNJ was detected in the CG group. The incidence of cerebrovascular and cardiovascular events in the EG was lower than that in the CG group (p < .05). MLD therapy improved antioxidant and anti-inflammatory properties and serum lipid profile (p < .05). The IMT values in the EG group were lower than those in the CG group (p < .05). Serum levels of DNJ had a strong negative relationship with IMT values. MLD treatment attenuates atherosclerotic lesions possibly via DNJ. PRACTICAL APPLICATIONS: Mulberry leaves Deoxynojirimycin (MLD) treatment improved antioxidant and anti-inflammatory properties and serum lipid profile in heart disease patients. MLD reduce carotid intima-media thickness (IMT) and serum levels of DNJ (the main ingredient of mulberry leaf) had a strong relationship with IMT values. MLD is an amazing treat used for reducing the incidence of cerebrovascular and cardiovascular events in heart disease patients.