Summary Cleistogenes songorica (2n = 4x = 40) is a desert grass with a unique dimorphic flowering mechanism and an ability to survive extreme drought. Little is known about the genetics underlying drought tolerance and its reproductive adaptability. Here, we sequenced and assembled a high‐quality chromosome‐level C. songorica genome (contig N50 = 21.28 Mb). Complete assemblies of all telomeres, and of ten chromosomes were derived. C. songorica underwent a recent tetraploidization (~19 million years ago) and four major chromosomal rearrangements. Expanded genes were significantly enriched in fatty acid elongation, phenylpropanoid biosynthesis, starch and sucrose metabolism, and circadian rhythm pathways. By comparative transcriptomic analysis we found that conserved drought tolerance related genes were expanded. Transcription of CsMYB genes was associated with differential development of chasmogamous and cleistogamous flowers, as well as drought tolerance. Furthermore, we found that regulation modules encompassing miRNA, transcription factors and target genes are involved in dimorphic flower development, validated by overexpression of CsAP2_9 and its targeted miR172 in rice. Our findings enable further understanding of the mechanisms of drought tolerance and flowering in C. songorica, and provide new insights into the adaptability of native grass species in evolution, along with potential resources for trait improvement in agronomically important species.

Summary Melilotus species are used as green manure and rotation crops worldwide and contain abundant pharmacologically active coumarins. However, there is a paucity of information on its genome and coumarin production and function. Here, we reported a chromosome‐scale assembly of Melilotus albus genome with 1.04 Gb in eight chromosomes, containing 71.42% repetitive elements. Long terminal repeat retrotransposon bursts coincided with declining of population sizes during the Quaternary glaciation. Resequencing of 94 accessions enabled insights into genetic diversity, population structure, and introgression. Melilotus officinalis had relatively larger genetic diversity than that of M. albus. The introgression existed between M. officinalis group and M. albus group, and gene flows was from M. albus to M. officinalis. Selection sweep analysis identified candidate genes associated with flower colour and coumarin biosynthesis. Combining genomics, BSA, transcriptomics, metabolomics, and biochemistry, we identified a β‐glucosidase (BGLU) gene cluster contributing to coumarin biosynthesis. MaBGLU1 function was verified by overexpression in M. albus, heterologous expression in Escherichia coli, and substrate feeding, revealing its role in scopoletin (coumarin derivative) production and showing that nonsynonymous variation drives BGLU enzyme activity divergence in Melilotus. Our work will accelerate the understanding of biologically active coumarins and their biosynthetic pathways, and contribute to genomics‐enabled Melilotus breeding.

Coumarins, natural products abundant in Melilotus albus, confer features in response to abiotic stresses, and are mainly present as glycoconjugates. UGTs (UDP-glycosyltransferases) are responsible for glycosylation modification of coumarins. However, information regarding the relationship between coumarin biosynthesis and stress-responsive UGTs remains limited. Here, a total of 189 MaUGT genes were identified from the M. albus genome, which were distributed differentially among its eight chromosomes. According to the phylogenetic relationship, MaUGTs can be classified into 13 major groups. Sixteen MaUGT genes were differentially expressed between genotypes of Ma46 (low coumarin content) and Ma49 (high coumarin content), suggesting that these genes are likely involved in coumarin biosynthesis. About 73.55% and 66.67% of the MaUGT genes were differentially expressed under ABA or abiotic stress in the shoots and roots, respectively. Furthermore, the functions of MaUGT68 and MaUGT186, which were upregulated under stress and potentially involved in coumarin glycosylation, were characterized by heterologous expression in yeast and Escherichia coli. These results extend our knowledge of the UGT gene family along with MaUGT gene functions, and provide valuable findings for future studies on developmental regulation and comprehensive data on UGT genes in M. albus.