It has been confirmed that the plant-specific Teosinte-branched 1/Cycloidea/Proliferating (TCP) gene family plays a pivotal role during plant growth and development. M. candidum is a native ornamental species and has a wide range of pharmacodynamic effects. However, there is still a lack of research on TCP’s role in controlling M. candidum’s development, abiotic stress responses and hormone metabolism. A comprehensive description of the TCP gene family in M. candidum is urgently needed. In this study, we used the HMMER search method in conjunction with the BLASTp method to identify the members of the TCP gene family, and a total of 35 TCP genes were identified. A domain analysis further confirmed that all 35 TCPs contained a TCP superfamily, a characteristic involved in dimerization and DNA binding that can be found in most genes from this gene family, suggesting that our identification was effective. As a result of the domain conservation analysis, the 35 TCP genes could be classified into two classes, TCP-P and TCP-C, based on the conservative regions of 55 and 59 amino acids, respectively. Gene-duplication analysis revealed that most TCP genes were present in duplication events that eventually led to TCP gene expansion in M. candidum. All the detected gene pairs had a Ka/Ks value of less than one, suggesting that purification selection is the most important factor that influences the evolution of TCP genes. Phylogenetic analysis of three species displayed the evolutionary relationship of TCP genes across different species and further confirmed our results. The real-time quantitative PCR (qRT-PCR) results showed that McTCP2a, McTCP7a, McTCP10, McTCP11, McTCP12a, McTCP13, McTCP16, McTCP17, McTCP18, McTCP20 and McTCP21 may be involved in leaf development; McTCP4a, McTCP1, McTCP14, McTCP17, McTCP18, McTCP20, McTCP22 and McTCP24 may be involved in flower development; and McTCP2a, McTCP3, McTCP5a, McTCP6, McTCP7a, McTCP9, McTCP11, McTCP14 and McTCP16 may be involved in seed development. Our results dissect the TCP gene family across the genome of M. candidum and provide valuable information for exploring TCP genes to promote molecular breeding and property improvement of M. candidum in the future.

M. candidum, an evergreen shrubby flower known for its superior adaptation ability in South China, has gained increased attention in garden applications. However, scant attention has been paid to its flower development and color formation process at the non-coding RNA level. To fill this gap, we conducted a comprehensive analysis based on long non-coding RNA sequencing (lncRNA-seq), RNA-seq, small RNA sequencing (sRNA-seq), and widely targeted metabolome detection of three different flower developmental stages of M. candidum. After differentially expressed lncRNAs (DElncRNAs), differentially expressed mRNAs (DEmRNAs), differentially expressed microRNAs (DEmiRNAs), and differentially synthesized metabolites (DSmets) analyses between the different flower developmental stages, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) were conducted to identify some key genes and metabolites in flavonoid, flavone, anthocyanin, carotenoid, and alkaloid-related GO terms and biosynthetic pathways. Three direct-acting models, including antisense-acting, cis-acting, and trans-acting between lncRNAs and mRNAs, were detected to illustrate the direct function of lncRNAs on target genes during flower development and color formation. Based on the competitive endogenous RNA (ceRNA) regulatory theory, we constructed a lncRNA-mediated regulatory network composed of DElncRNAs, DEmiRNAs, DEmRNAs, and DSmets to elucidate the indirect role of lncRNAs in the flower development and color formation of M. candidum. By utilizing correlation analyses between DERNAs and DSmets within the ceRNA regulatory network, alongside verification trials of the ceRNA regulatory mechanism, the study successfully illustrated the significance of lncRNAs in flower development and color formation process. This research provides a foundation for improving and regulating flower color at the lncRNA level in M. candidum, and sheds light on the potential applications of non-coding RNA in studies of flower development.

The mitogenome of Bauhinia variegate was assembled and characterized in this study. The mitogenome size was 437,271 bp, and its GC content was 45.5%. 36 protein-coding genes, 17 tRNAs and 3 rRNAs were annotated in the mitogenome. A total of 12 MTPTs, ranging from 71 bp to 3562 bp, were identified in the mitogenome and covered 1.46% (6373 bp) of the mitogenome. Phylogenetic analysis of 15 species of Leguminosae based on 23 core protein-coding genes showed that B. variegata was sister to Tylosema esculentum, another member from the subfamily Cercidoideae. The mitogenome of B. variegata provides a valuable genetic resource for further phylogenetic studies of this family.

Cercidoideae, one of the six subfamilies of Leguminosae, contains one genus Cercis with its chromosome number 2n = 14 and all other genera with 2n = 28. An allotetraploid origin hypothesis for the common ancestor of non-Cercis genera in this subfamily has been proposed; however, no chromosome-level genomes from Cercidoideae have been available to test this hypothesis. Here, we conducted a chromosome-level genome assembly of Bauhinia variegata to test this hypothesis. The assembled genome is 326.4 Mb with the scaffold N50 of 22.1 Mb and contains 37,996 protein-coding genes. The Ks distribution between gene pairs in the syntenic regions indicates two whole-genome duplications (WGDs): one is B. variegata-specific, and the other is shared among core eudicots. Although Ks between gene pairs generated by the recent WGD in Bauhinia is greater than that between Bauhinia and Cercis, the WGD was not detected in Cercis, which can be explained by an accelerated evolutionary rate in Bauhinia after divergence from Cercis. Ks distribution and phylogenetic analysis for gene pairs generated by the recent WGD in Bauhinia and their corresponding orthologs in Cercis support the allopolyploidy origin hypothesis of Bauhinia. The genome of B. variegata also provides a genomic resource for dissecting genetic basis of its ornamental traits.

The MYB gene family is significant in plants, playing a role in numerous plant development processes, including metabolism, hormone signal transduction, cell identity, and biotic and abiotic stresses. Due to the recent availability of the Melastoma candidum genome, this is the first time that the MYB gene family has been identified in this species. This study identified 421 MYB gene members in the M. candidum genome using the HMMER search and BLASTp method. These MYBs were further divided into 10 sub-types, including R2R3, R-R, CPC-like, CCA1-like, TBP-like, R1R2R3, I-box, atypical MYB, MYB-CC, and MYB-like. Domain and conservation analyses revealed that each type of MYB was characterized by a different number and combination of SANTs/myb DNA-binding domains. Collinearity analysis revealed several gene duplication events within the MYB gene family. The Ka to Ks ratio suggested that most of the MYB genes underwent purifying selection during the evolution process. Phylogenetic analysis among three species confirmed our findings and displayed the evolutionary relationship of MYB genes in different species. RNA-seq of three developmental stages of flowers and WGCNA analysis identified McMYB113h, McMYB21b, and McGLK1c as playing a pivotal role during flower development in M. candidum. Finally, we conducted qRT-PCR experiments for 20 flower-development-related MYBs across 9 tissues to illustrate their expression patterns in M. candidum. This study establishes a foundation for exploring MYB gene resources and their potential applications in related industries of M. candidum.

Melastoma, consisting of ~100 species diversified in tropical Asia and Oceania in the past 1-2 million years, represents an excellent example of rapid speciation in flowering plants. Trichomes on hypanthia, twigs and leaves vary markedly among species of this genus and are the most important diagnostic traits for species identification. These traits also play critical roles in contributing to differential adaptation of these species to their own habitats. Here we sequenced the genome ofM. candidum, a common, erect-growing species from southern China, with the aim to provide genomic insights into trichome evolution in this genus. We generated a high-quality, chromosome-level genome assembly ofM. candidum, with the genome size of 256.2 Mb and protein-coding gene number of 40,938. The gene families specific to, and significantly expanded inMelastomaare enriched for GO terms related to trichome initiation and differentiation. We provide evidence thatMelastomaand its sister genusOsbeckiahave undergone two whole genome duplications (WGDs) after the triplication event (γ) shared by all core eudicots. Preferential retention of trichome development-related transcription factor genes such as C2H2, bHLH, HD-ZIP, WRKY, and MYB after both WGDs might provide raw materials for trichome evolution and thus contribute to rapid species diversification inMelastoma. Our study provides candidate transcription factor genes related to trichome evolution inMelastoma, which can be used to evolutionary and functional studies of trichome diversification among species of this genus.