Zhiyong Li

Background: Medicago ruthenica, a perennial legume forage valuable for ecological restoration and improved breeding, suffers significant harvest losses due to pod shattering. Pod shattering is a trait not only linked to not only pod ventral suture, but also pericarp properties. In this study, we aimed to (1) elucidate the role of pericarp in explosive pod shattering by comparing shattering-susceptible (SPD) and shattering-resistant (RPD) M. ruthenica genotypes, and (2) identify key regulatory genes and pathways underlying this mechanism. Methods: We conducted comparative analyses of pericarp anatomy and physiological traits (pericarp components such as water content, cellulose, hemicellulose, pectin, and lignin; and the activities of enzymes such as cellulose synthase A (CesA), phenylalanine ammonia-lyase (PAL), 4-coumarate: CoA ligase (4CL), cinnamyl alcohol dehydrogenase (CAD), and peroxidase (POD) in SPD and RPD pods). Transcriptome of pod pericarps identified differentially expressed genes (DEGs) for the selection of candidates functional genes. Promoter analysis was performed on candidate functional genes to identify specific regulated factors. The functional role of auxin signaling was validated through exogenous auxin application and the assessment of pod shattering rates and gene expression. Results: SPD pod pericarps exhibited significantly higher lignification of endocarp, lignin, cellulose, hemicellulose and pectin content, but lower water content than RPD. Principal component analysis identified that lignin contributes the highest loading value (0.727) contributor to pod shattering. The activities of five cell wall biosynthesis enzymes were higher in SPD pod pericarps than RPD. Transcriptome analysis identified more than 3419 DEGs in SPD pericarps. KEGG enrichment highlighted “phenylpropanoid biosynthesis” as the most significant pathway. A total of 57 lignin-biosynthesis-related DEGs were upregulated in SPD, including 15 PODs. Promoters of 11 POD genes contained MYB-binding motifs and 8 contained auxin-responsive elements, a total of 76 MYB transcription factors (mostly upregulated) and 9 auxin biosynthesis genes (mostly downregulated) were differentially expressed in SPD. Exogenous auxin application significantly reduced SPD pod shattering to 23.6% and concurrently downregulated PODs expression. Conclusions: This study establishes that enhanced lignification within the pericarp endocarp by the upregulation of lignin biosynthetic genes (particularly PODs), coupled with upregulation by MYB transcription factors and downregulation by auxin, is a core mechanism of explosive pod shattering in M. ruthenica. The identified DEGs, especially MYBs, PODs, and auxin pathway genes, provide gene information for breeding shattering-resistant M. ruthenica varieties through molecular design or marker-assisted selection.

Alfalfa (Medicago sativa) is a globally distributed economic legume crop used for forage and ecological restoration. We aimed to explore the mechanisms underlying the cold acclimation observed in this species. Our results for fall plant growth showed that non-dormant alfalfa (SD) maintained a vigorous growth rate compared to that of fall-dormant alfalfa (ZD); however, the winter survival rate of ZD was higher than that of SD. Among the ZD samples, the starch content first accumulated and then decreased; the sucrose content was consumed first along with simultaneous raffinose accumulation, which was followed by sucrose content accumulation, with consistent changes in the corresponding related synthase and hydrolase activity. SD exhibited the opposite trend. The transcriptome data showed that most of the differentially expressed genes were involved in carbon metabolism (ko01200), amino acid biosynthesis (ko01230), and starch and sucrose metabolism (ko00500). Our data clearly show that alfalfa’s cold acclimation mechanism is a complex process, with the establishment of stable carbon homeostasis; sucrose is first converted into starch and raffinose, and then, starch is converted into sucrose, which enables alfalfa’s cold resistance. The process is accompanied by CBF/DREB1A TF regulation. This study provides important insights into the cold acclimation mechanisms of alfalfa.

Floral bud growth influences seed yield and quality; however, the molecular mechanism underlying the development of floral buds in alfalfa (Medicago sativa) is still unclear. Here, we comprehensively analyzed the transcriptome and targeted metabolome across the early, mid, and late bud developmental stages (D1, D2, and D3) in alfalfa. The metabolomic results revealed that gibberellin (GA), auxin (IAA), cytokinin (CK), and jasmonic acid (JA) might play an essential role in the developmental stages of floral bud in alfalfa. Moreover, we identified some key genes associated with GA, IAA, CK, and JA biosynthesis, including CPS, KS, GA20ox, GA3ox, GA2ox, YUCCA6, amid, ALDH, IPT, CYP735A, LOX, AOC, OPR, MFP2, and JMT. Additionally, many candidate genes were detected in the GA, IAA, CK, and JA signaling pathways, including GID1, DELLA, TF, AUX1, AUX/IAA, ARF, GH3, SAUR, AHP, B-ARR, A-ARR, JAR1, JAZ, and MYC2. Furthermore, some TFs related to flower growth were screened in three groups, such as AP2/ERF-ERF, MYB, MADS-M-type, bHLH, NAC, WRKY, HSF, and LFY. The findings of this study revealed the potential mechanism of floral bud differentiation and development in alfalfa and established a theoretical foundation for improving the seed yield of alfalfa.

The petals of Medicago sativa ssp. sativa and M. sativa ssp. falcata are purple and yellow, respectively. Free hybridization between M. sativa ssp. sativa and M. sativa ssp. falcata has created hybrids with various flower colors in nature. Moreover, the flower colors of alfalfa are closely correlated with yield, nutritional quality, stress tolerance and other agronomic characteristics. To elucidate the underlying mechanisms of flower color formation in M. sativa ssp. sativa and M. sativa ssp. falcata, we conducted an integrative analysis of the transcriptome and metabolome of alfalfa with three different petal colors (purple, yellow and cream). The metabolic profiles suggested that anthocyanins and carotenoids are the crucial pigments in purple and yellow flowers, respectively. A quantitative exploration of the anthocyanin and carotenoid components indicated that the accumulations of cyanidin, delphinidin, peonidin, malvidin, pelargonidin and petunidin derivatives are significantly higher in purple flowers than in cream flowers. In addition, the content of carotenes (phytoene, α-carotene and β-carotene) and xanthophylls (α-cryptoxanthin, lutein, β-cryptoxanthin, zeaxanthin, antheraxanthin and violaxanthin derivatives) was markedly higher in yellow flowers than in cream flowers. Furthermore, we found that delphinidin-3,5-O-diglucoside and lutein were the predominant pigments accumulated in purple and yellow flowers, respectively. The transcriptomic results revealed that twenty-five upregulated structural genes (one C4H, three 4CL, twelve CHS, two CHI, one F3H, one F3′H, one F3′5′H and four DFR) are involved in the accumulation of anthocyanins in purple flowers, and nine structural genes (two PSY, one ZDS, two CRTISO, two BCH, one ZEP and one ECH) exert an effect on the carotenoid biosynthesis pathway in yellow flowers. The findings of this study reveal the underlying mechanisms of anthocyanin and carotenoid biosynthesis in alfalfa with three classic flower colors.

Context: Red clover (Trifolium pratense L.) is an important legume forage in temperate agricultural zones. Evaluation of self- and cross-pollination fertility is important for setting up an effective breeding-program scale. However, the outcrossing rate of red clover under open-pollination conditions is not certain. Development of a reliable and time-saving marker system is needed to quantify and characterise outcrossing rates. Aim: We aimed to develop a duplex PCR-based protocol based on a genome-wide simple sequence repeat (SSR) screen, and to determine the outcrossing rate of red clover under open-pollination environments. Methods: We screened 209 SSR markers with pooled DNA samples of 60 plants from 20 red clover accessions, and selected 185 SSR markers that produced clear scorable bands for testing with 24 individual DNA samples to determine polymorphism. We selected 70 primer pairs, and then assembled a core set of 24 loci into 12 sets of duplex markers, which were used for outcrossing behaviour analysis of 60 maternal parents and their respective 22 half-sib progenies. Key results: Mean polymorphic information content (PIC) for the 70 markers was 0.490 (range 0.117–0.878). Minimum, mean and maximum PIC values for the 24 markers constituting the 12 duplexes were 0.226, 0.594 and 0.781, respectively. The outcrossing rate was identified as 99.4% for red clover in a natural environment. Conclusion: We successfully developed a duplex SSR-based PCR protocol consisting of 24 markers. This SSR system was applied to determine the outcrossing rate of red clover in a natural environment.

Medicago ruthenica, a wild legume forage widely distributed in the Eurasian steppe, demonstrates high genetic and phenotypic variation. M. ruthenica with a purely yellow flower (YFM), differing from the general phenotype of M. ruthenica with a purple flower (PFM), was recently discovered. The similar characteristics of YFM with Medicago falcata have led to conflicting opinions on its taxonomy using traditional morphological methods. The lack of chemotaxonomy information about M. ruthenica species and the unclear flower coloration mechanisms have hampered their study. Here, we investigated M. ruthenica using targeted metabolomics based on the chemotaxonomy method and elaborated the floral coloration mechanisms using transcriptomics. The identified flavonoids were the same types, but there were different contents in YFM and PFM, especially the contents of cyanidin-3-O-glucoside (C3G), an anthocyanin that causes the purple-reddish color of flowers. The over-accumulation of C3G in PFM was 1,770 times more than YFM. Nineteen anthocyanin-related genes were downregulated in YFM compared with their expression in PFM. Thus, YFM could be defined as a variety of M. ruthenica rather than a different species. The loss of purple flower coloration in YFM was attributed to the downregulation of these genes, resulting in reduced C3G accumulation. The taxonomic characteristics and molecular and physiological characteristics of this species will contribute to further research on other species with similar external morphologies.

Astragalus is the largest genus in Leguminosae. Several molecular studies have investigated the potential adulterants of the species within this genus; nonetheless, the evolutionary relationships among these species remain unclear. Herein, we sequenced and annotated the complete chloroplast genomes of three Astragalus species—Astragalus adsurgens, Astragalus mongholicus var. dahuricus, and Astragalus melilotoides using next-generation sequencing technology and plastid genome annotator (PGA) tool. All species belonged to the inverted repeat lacking clade (IRLC) and had similar sequences concerning gene contents and characteristics. Abundant simple sequence repeat (SSR) loci were detected, with single-nucleotide repeats accounting for the highest proportion of SSRs, most of which were A/T homopolymers. Using Astragalus membranaceus var. membranaceus as reference, the divergence was evident in most non-coding regions of the complete chloroplast genomes of these species. Seven genes (atpB, psbD, rpoB, rpoC1, trnV, rrn16, and rrn23) showed high nucleotide variability (Pi), and could be used as DNA barcodes for Astragalus sp. cemA and rpl33 were found undergoing positive selection by the section patterns in the coded protein. Phylogenetic analysis showed that Astragalus is a monophyletic group closely related to the genus Oxytropis within the tribe Galegeae. The newly sequenced chloroplast genomes provide insight into the unresolved evolutionary relationships within Astragalus spp. and are expected to contribute to species identification.

Red clover (Trifolium pratense L.) is a globally significant legume having economic value as a forage and green manure crop in temperate agricultural zones. Based on the gametophytic self-incompatibility system, and analysis of seed set rates between self and cross pollination, the clover has long been considered an outcrossing species. However, the outcrossing rates of red clover under open-pollination conditions are not definitive. Development of a reliable, timesaving, and easily usable marker system is needed to quantify and characterize rates of red clover crossing and selfing. Here, genome-wide screening of 209 mapped simple sequence repeat (SSR) markers was conducted, and 185 produced clear scorable bands on a pooled DNA sample of 20 red clover accessions. Seventy markers were selected based on polymerase chain reaction (PCR) amplification quality on 24 genotypes and their relatively even distribution in red clover genome. Mean polymorphic information content for the 70 markers was 0.490, ranging from 0.117 to 0.878. From the markers, a core set of 24 loci, which had been mapped on six linkage groups, was further tested to develop 12 sets of duplex markers. Using the established duplex PCR protocol, 10 out of 12 sets of duplex SSR markers were used to genotype 60 maternal parents and their respective 22 half-sib progenies. Eight plants among the 1,320 progenies were identified to be selfed, indicating that the outcrossing rate was 99.4% in a natural environment. The protocol and the finding of this rare selfing strategy under open-pollinating conditions contribute to red clover breeding efforts.