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.

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.

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).

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.

Genome Sequencing Analysis of Scleromitrula shiraiana, a Causal Agent of Mulberry Sclerotial Disease With Narrow Host Range.

Scleromitrula shiraiana is a necrotrophic fungus with a narrow host range, and is one of the main causal pathogens of mulberry sclerotial disease. However, its molecular mechanisms and pathogenesis are unclear. Here, we report a 39.0 Mb high-quality genome sequence for S. shiraiana strain SX-001. The S. shiraiana genome contains 11,327 protein-coding genes. The number of genes and genome size of S. shiraiana are similar to most other Ascomycetes. The cross-similarities and differences of S. shiraiana with the closely related Sclerotinia sclerotiorum and Botrytis cinerea indicated that S. shiraiana differentiated earlier from their common ancestor. A comparative genomic analysis showed that S. shiraiana has fewer genes encoding cell wall-degrading enzymes (CWDEs) and effector proteins than that of S. sclerotiorum and B. cinerea, as well as many other Ascomycetes. This is probably a key factor in the weaker aggressiveness of S. shiraiana to other plants. S. shiraiana has many species-specific genes encoding secondary metabolism core enzymes. The diversity of secondary metabolites may be related to the adaptation of these pathogens to specific ecological niches. However, melanin and oxalic acid are conserved metabolites among many Sclerotiniaceae fungi, and may be essential for survival and infection. Our results provide insights into the narrow host range of S. shiraiana and its adaptation to mulberries.

Molecular phylogeny of mulberries reconstructed from ITS and two cpDNA sequences.

BACKGROUND: Species in the genus Morus (Moraceae) are deciduous woody plants of great economic importance. The classification and phylogenetic relationships of Morus, especially the abundant mulberry resources in China, is still undetermined. Internal transcribed spacer (ITS) regions are among the most widely used molecular markers in phylogenetic analyses of angiosperms. However, according to the previous phylogenetic analyses of ITS sequences, most of the mulberry accessions collected in China were grouped into the largest clade lacking for phylogenetic resolution. Compared with functional ITS sequences, ITS pseudogenes show higher sequence diversity, so they can provide useful phylogenetic information. METHODS: We sequenced the ITS regions and the chloroplast DNA regions TrnL-TrnF and TrnT-TrnL from 33 mulberry accessions, and performed phylogenetic analyses to explore the evolution of mulberry. RESULTS: We found ITS pseudogenes in 11 mulberry accessions. In the phylogenetic tree constructed from ITS sequences, clade B was separated into short-type sequence clades (clades 1 and 2), and a long-type sequence clade (clade 3). Pseudogene sequences were separately clustered into two pseudogroups, designated as pseudogroup 1 and pseudogroup 2. The phylogenetic tree generated from cpDNA sequences also separated clade B into two clades. CONCLUSIONS: Two species were separated in clade B. The existence of three connection patterns and incongruent distribution patterns between the phylogenetic trees generated from cpDNA and ITS sequences suggested that the ITS pseudogene sequences connect with genetic information from the female progenitor. Hybridization has played important roles in the evolution of mulberry, resulting in low resolution of the phylogenetic analysis based on ITS sequences. An evolutionary pattern illustrating the evolution history of mulberry is proposed. These findings have significance for the conservation of local mulberry resources. Polyploidy, hybridization, and concerted evolution have all played the roles in the evolution of ITS sequences in mulberry. This study will expand our understanding of mulberry evolution.

Cloning and Phylogenetic Analysis of Brassica napus L. Caffeic Acid O-Methyltransferase 1 Gene Family and Its Expression Pattern under Drought Stress.

For many plants, regulating lignin content and composition to improve lodging resistance is a crucial issue. Caffeic acid O-methyltransferase (COMT) is a lignin monomer-specific enzyme that controls S subunit synthesis in plant vascular cell walls. Here, we identified 12 BnCOMT1 gene homologues, namely BnCOMT1-1 to BnCOMT1-12. Ten of 12 genes were composed of four highly conserved exons and three weakly conserved introns. The length of intron I, in particular, showed enormous diversification. Intron I of homologous BnCOMT1 genes showed high identity with counterpart genes in Brassica rapa and Brassica oleracea, and intron I from positional close genes in the same chromosome were relatively highly conserved. A phylogenetic analysis suggested that COMT genes experience considerable diversification and conservation in Brassicaceae species, and some COMT1 genes are unique in the Brassica genus. Our expression studies indicated that BnCOMT1 genes were differentially expressed in different tissues, with BnCOMT1-4, BnCOMT1-5, BnCOMT1-8, and BnCOMT1-10 exhibiting stem specificity. These four BnCOMT1 genes were expressed at all developmental periods (the bud, early flowering, late flowering and mature stages) and their expression level peaked in the early flowering stage in the stem. Drought stress augmented and accelerated lignin accumulation in high-lignin plants but delayed it in low-lignin plants. The expression levels of BnCOMT1s were generally reduced in water deficit condition. The desynchrony of the accumulation processes of total lignin and BnCOMT1s transcripts in most growth stages indicated that BnCOMT1s could be responsible for the synthesis of a specific subunit of lignin or that they participate in other pathways such as the melatonin biosynthesis pathway.

Molecular Mapping and QTL for Expression Profiles of Flavonoid Genes in Brassica napus.

Flavonoids are secondary metabolites that are extensively distributed in the plant kingdom and contribute to seed coat color formation in rapeseed. To decipher the genetic networks underlying flavonoid biosynthesis in rapeseed, we constructed a high-density genetic linkage map with 1089 polymorphic loci (including 464 SSR loci, 97 RAPD loci, 451 SRAP loci, and 75 IBP loci) using recombinant inbred lines (RILs). The map consists of 19 linkage groups and covers 2775 cM of the B. napus genome with an average distance of 2.54 cM between adjacent markers. We then performed expression quantitative trait locus (eQTL) analysis to detect transcript-level variation of 18 flavonoid biosynthesis pathway genes in the seeds of the 94 RILs. In total, 72 eQTLs were detected and found to be distributed among 15 different linkage groups that account for 4.11% to 52.70% of the phenotypic variance atrributed to each eQTL. Using a genetical genomics approach, four eQTL hotspots together harboring 28 eQTLs associated with 18 genes were found on chromosomes A03, A09, and C08 and had high levels of synteny with genome sequences of A. thaliana and Brassica species. Associated with the trans-eQTL hotspots on chromosomes A03, A09, and C08 were 5, 17, and 1 genes encoding transcription factors, suggesting that these genes have essential roles in the flavonoid biosynthesis pathway. Importantly, bZIP25, which is expressed specifically in seeds, MYC1, which controls flavonoid biosynthesis, and the R2R3-type gene MYB51, which is involved in the synthesis of secondary metabolites, were associated with the eQTL hotspots, and these genes might thus be involved in different flavonoid biosynthesis pathways in rapeseed. Hence, further studies of the functions of these genes will provide insight into the regulatory mechanism underlying flavonoid biosynthesis, and lay the foundation for elaborating the molecular mechanism of seed coat color formation in B. napus.