Leguminosae is one of the three largest families of angiosperms after Compositae and Orchidaceae. It is widely distributed and grows in a variety of environments, including plains, mountains, deserts, forests, grasslands, and even waters where almost all legumes can be found. It is one of the most important sources of starch, protein and oil in the food of mankind and also an important source of high-quality forage material for animals, which has important economic significance. In our study, the codon usage patterns and variation sources of the chloroplast genome of nine important forage legumes were systematically analyzed. Meanwhile, we also constructed a phylogenetic tree based on the whole chloroplast genomes and protein coding sequences of these nine forage legumes. Our results showed that the chloroplast genomes of nine forage legumes end with A/T bases, and seven identical high-frequency (HF) codons were detected among the nine forage legumes. ENC-GC3s mapping, PR2 analysis, and neutral analysis showed that the codon bias of nine forage legumes was influenced by many factors, among which natural selection was the main influencing factor. The codon usage frequency showed that the Nicotiana tabacum and Saccharomyces cerevisiae can be considered as receptors for the exogenous expression of chloroplast genes of these nine forage legumes. The phylogenetic relationships of the chloroplast genomes and protein coding genes were highly similar, and the nine forage legumes were divided into three major clades. Among the clades Melilotus officinalis was more closely related to Medicago sativa, and Galega officinalis was more closely related to Galega orientalis. This study provides a scientific basis for the molecular markers research, species identification and phylogenetic studies of forage legumes.

Calcium-dependent protein kinases (CDPKs) are an important class of calcium-sensitive response proteins that play an important regulatory role in response to abiotic stresses. To date, little is known about the CDPK genes in white clover. White clover is a high-quality forage grass with high protein content, but it is susceptible to cold stress. Therefore, we performed a genome-wide analysis of the CDPK gene family in white clover and identified 50 members of the CDPK genes. Phylogenetic analysis using CDPKs from the model plant Arabidopsis divided the TrCDPK genes into four groups based on their sequence similarities. Motif analysis showed that TrCDPKs within the same group had similar motif compositions. Gene duplication analysis revealed the evolution and expansion of TrCDPK genes in white clover. Meanwhile, a genetic regulatory network (GRN) containing TrCDPK genes was reconstructed, and gene ontology (GO) annotation analysis of these functional genes showed that they contribute to signal transduction, cellular response to stimuli, and biological regulation, all of which are important processes in response to abiotic stresses. To determine the function of TrCDPK genes, we analyzed the RNA-seq dataset and found that most TrCDPK genes were highly up-regulated under cold stress, particularly in the early stages of cold stress. These results were validated by qRT-PCR experiments, implying that TrCDPK genes are involved in various gene regulatory pathways in response to cold stress. Our study may help to further investigate the function of TrCDPK genes and their role in response to cold stress, which is important for understanding the molecular mechanisms of cold tolerance in white clover and improving its cold tolerance.

Alfalfa (Medicago sativa) is a perennial forage legume that is widely distributed all over the world; therefore, it has an extremely complex genetic background. Though population structure and phylogenetic studies have been conducted on a large group of alfalfa nuclear genomes, information about the chloroplast genomes is still lacking. Chloroplast genomes are generally considered to be conservative and play an important role in population diversity analysis and species adaptation in plants. Here, 231 complete alfalfa chloroplast genomes were successfully assembled from 359 alfalfa resequencing data, on the basis of which the alfalfa chloroplast pan-genome was constructed. We investigated the genetic variations of the alfalfa chloroplast genome through comparative genomic, genetic diversity, phylogenetic, population genetic structure, and haplotype analysis. Meanwhile, the expression of alfalfa chloroplast genes under cold stress was explored through transcriptome analysis. As a result, chloroplast genomes of 231 alfalfa lack an IR region, and the size of the chloroplast genome ranges from 125,192 bp to 126,105 bp. Using population structure, haplotypes, and construction of a phylogenetic tree, it was found that alfalfa populations could be divided into four groups, and multiple highly variable regions were found in the alfalfa chloroplast genome. Transcriptome analysis showed that tRNA genes were significantly up-regulated in the cold-sensitive varieties, while rps7, rpl32, and ndhB were down-regulated, and the editing efficiency of ycf1, ycf2, and ndhF was decreased in the cold-tolerant varieties, which may be due to the fact that chloroplasts store nutrients through photosynthesis to resist cold. The huge number of genetic variants in this study provide powerful resources for molecular markers.