Alfalfa, an important forage legume, is an ideal crop for sustainable agriculture and a potential crop for bioenergy resources. Drought, one of the most common environmental stresses, substantially affects plant growth, development, and productivity. MicroRNAs (miRNAs) are newly discovered gene expression regulators that have been linked to several plant stress responses. To elucidate the role of miRNAs in drought stress regulation of alfalfa, a high-throughput sequencing approach was used to analyze 12 small RNA libraries comprising of four samples, each with three biological replicates. From the 12 libraries, we identified 348 known miRNAs belonging to 80 miRNA families, and 281 novel miRNAs, using Mireap software. Eighteen known miRNAs in roots and 12 known miRNAs in leaves were screened as drought-responsive miRNAs. With the exception of miR319d and miR157a which were upregulated under drought stress, the expression pattern of drought-responsive miRNAs was different between roots and leaves in alfalfa. This is the first study that has identified miR3512, miR3630, miR5213, miR5294, miR5368 and miR6173 as drought-responsive miRNAs. Target transcripts of drought-responsive miRNAs were computationally predicted. All 447 target genes for the known miRNAs were predicted using an online tool. This study provides a significant insight on understanding drought-responsive mechanisms of alfalfa.

Anthraquinones constitute the largest group of natural quinones, which are used as safe natural dyes and have many pharmaceutical applications. In plants, anthraquinones are biosynthesized through two main routes: the polyketide pathway and the shikimate pathway. The latter primarily forms alizarin-type anthraquinones, and the prenylation of 1,4-dihydroxy-2-naphthoic acid is the first pathway-specific step. However, the prenyltransferase responsible for this key step remains uncharacterized. In this study, the cell suspension culture of Madder (Rubia cordifolia), a plant rich in alizarin-type anthraquinones, was discovered to be capable of prenylating 1,4-dihydroxy-2-naphthoic acid to form 2-carboxyl-3-prenyl-1,4-naphthoquinone and 3-prenyl-1,4-naphthoquinone. Then, a candidate gene belonging to the UbiA superfamily, R. cordifolia dimethylallyltransferase 1 (RcDT1), was shown to account for the prenylation activity. Substrate specificity studies revealed that the recombinant RcDT1 recognized naphthoic acids primarily, followed by 4-hydroxyl benzoic acids. The prenylation activity was strongly inhibited by 1,2- and 1,4-dihydroxynaphthalene. RcDT1 RNA interference significantly reduced the anthraquinones content in R. cordifolia callus cultures, demonstrating that RcDT1 is required for alizarin-type anthraquinones biosynthesis. The plastid localization and root-specific expression further confirmed the participation of RcDT1 in anthraquinone biosynthesis. The phylogenetic analyses of RcDT1 and functional validation of its rubiaceous homologs indicated that DHNA-prenylation activity evolved convergently in Rubiaceae via recruitment from the ubiquinone biosynthetic pathway. Our results demonstrate that RcDT1 catalyzes the first pathway-specific step of alizarin-type anthraquinones biosynthesis in R. cordifolia. These findings will have profound implications for understanding the biosynthetic process of the anthraquinone ring derived from the shikimate pathway.

The basic helix–loop–helix (bHLH) transcription factors possess DNA-binding and dimerization domains and are involved in various biological and physiological processes, such as growth and development, the regulation of secondary metabolites, and stress response. However, the bHLH gene family in C. tinctorius has not been investigated. In this study, we performed a genome-wide identification and analysis of bHLH transcription factors in C. tinctorius. A total of 120 CtbHLH genes were identified, distributed across all 12 chromosomes, and classified into 24 subfamilies based on their phylogenetic relationships. Moreover, the 120 CtbHLH genes were subjected to comprehensive analyses, including protein sequence alignment, evolutionary assessment, motif prediction, and the analysis of promoter cis-acting elements. The promoter region analysis revealed that CtbHLH genes encompass cis-acting elements and were associated with various aspects of plant growth and development, responses to phytohormones, as well as responses to both abiotic and biotic stresses. Expression profiles, sourced from transcriptome databases, indicated distinct expression patterns among these CtbHLH genes, which appeared to be either tissue-specific or specific to certain cultivars. To further explore their functionality, we determined the expression levels of fifteen CtbHLH genes known to harbor motifs related to abiotic and hormone responses. This investigation encompassed treatments with ABA, salt, drought, and MeJA. The results demonstrated substantial variations in the expression patterns of CtbHLH genes in response to these abiotic and hormonal treatments. In summary, our study establishes a solid foundation for future inquiries into the roles and regulatory mechanisms of the CtbHLH gene family.