MicroRNAs (miRNAs) are a class of non-coding, small RNAs recognized as important regulators of gene expression. Although plant miRNAs have been extensively studied in model systems, less is known in other plants with limited genome sequence data, including Paeonia ostii. In this work, we used high-throughput sequencing to identify conserved and nonconserved miRNAs and other short RNAs in Paeonia ostii under control and copper stressed condition. 102 previously known plant miRNAs were identified and classified into 89 families according to their gene sequence identity. Some miRNAs were highly conserved in the plant kingdom suggesting that these miRNA play important and conserved roles. Combined our transcriptome sequencing data of Paeonia ostii under same conditions, 34 novel potential miRNAs were identified. The potential targets of the identified known and novel miRNAs were also predicted based on sequence homology search. Comparing the two libraries, it was observed that 12 conserved miRNAs and 18 novel miRNAs showed significantly changes in response to copper stress. Some of the new identified potential miRNAs might be involved in Paeonia ostii-specific regulating mechanisms under copper stress. These results provide a framework for further analysis of miRNAs and their role in regulating Paeonia ostii response to copper stress.

Background: The basic leucine zipper (bZIP) family is a predominant group of transcription factors in plants, involved in regulating plant growth, development, and response to stressors. Additionally, the bZIP gene family has a key role in anthocyanin production. Despite the significant role of bZIP genes in plants, their potential contribution in lotus remains understudied. Results: A total of 124 bZIP genes (59 NnbZIPs and 65 NlbZIPs) were identified from genomes of two lotus species. These genes were classified into 13 groups according to the grouping principle of the Arabidopsis bZIP gene family. Analysis of promoter cis-acting elements indicated that most bZIP gene family members in lotus are associated with response to abiotic stresses. The promoters of some bZIP genes contain MYB binding sites that regulate anthocyanin synthesis. We examined the anthocyanin content of the petals from three different colored lotus, combined with transcriptome data analysis and qRT-PCR results, showing that the expression trends of NnbZIP36 and the homologous gene NlbZIP38 were significantly correlated with the anthocyanin content in lotus petals. Furthermore, we found that overexpression of NnbZIP36 in Arabidopsis promoted anthocyanin accumulation by upregulating the expression of genes (4CL, CHI, CHS, F3H, F3'H, DFR, ANS and UF3GT) related to anthocyanin synthesis. Conclusions: Our study enhances the understanding of the bZIP gene family in lotus and provides evidence for the role of NnbZIP36 in regulating anthocyanin synthesis. This study also sets the stage for future investigations into the mechanism by which the bZIP gene family regulates anthocyanin biosynthesis in lotus.

Trehalose-6-phosphate synthase (TPS) plays a key role in plant carbohydrate metabolism and the perception of carbohydrate availability. In the present work, the publicly available Nelumbo nucifera (lotus) genome sequence database was analyzed which led to identification of nine lotus TPS genes (NnTPS). It was found that at least two introns are included in the coding sequences of NnTPS genes. When the motif compositions were analyzed we found that NnTPS generally shared the similar motifs, implying that they have similar functions. The dN/dS ratios were always less than 1 for different domains and regions outside domains, suggesting purifying selection on the lotus TPS gene family. The regions outside TPS domain evolved relatively faster than NnTPS domains. A phylogenetic tree was constructed using all predicted coding sequences of lotus TPS genes, together with those from Arabidopsis, poplar, soybean and rice. The result indicated that those TPS genes could be clearly divided into two main subfamilies (I-II), where each subfamily could be further divided into 2 (I) and 5 (II) subgroups. Analyses of divergence and adaptive evolution show that purifying selection may have been the main force driving evolution of plant TPS genes. Some of the critical sites that contributed to divergence may have been under positive selection. Transcriptome data analysis revealed that most NnTPS genes were predominantly expressed in sink tissues. Expression pattern of NnTPS genes under copper and submergence stress indicated that NNU_014679 and NNU_022788 might play important roles in lotus energy metabolism and participate in stress response. Our results can facilitate further functional studies of TPS genes in lotus.