Alfalfa (Medicago sativa L.) is a perennial forage legume that is widely distributed throughout the world, and cold stress is an important environmental factor limiting the growth and production of alfalfa in cold regions. However, little is known of the molecular mechanisms regarding cold tolerance in alfalfa. Here, we conducted physiological metabolism assays and pan-transcriptome sequencing on eight cultivars of alfalfa under cold stress conditions. The results of the RNA-seq analysis showed that the genes are “oxidoreductase activity” and “transcription regulator activity”, suggesting that genes with such functions are more likely to play important roles in the response to cold stress by alfalfa. In addition, to identify specific gene modules and hub genes in response to alfalfa cold stress, we applied weighted gene co-expression network (WGCNA) analyses to the RNA-seq data. Our results indicate that the modules of genes that focus on the ATPase complex, ribosome biogenesis, are more likely to be involved in the alfalfa response to cold stress. It is important to note that we identified two fibronectin (FIB) genes as hub genes in alfalfa in response to cold stress and that they negatively regulate alfalfa response to chilling stress, and it is possible that dormant alfalfa is more effective at down-regulating FIB expression and therefore more resistant to cold stress.

A perennial leguminous forage, Medicago ruthenica has outstanding tolerance to abiotic stresses. The genome of Medicago ruthenica is large and has a complex genetic background, making it challenging to accurately determine genetic information. However, the chloroplast genome is widely used for researching issues related to evolution, genetic diversity, and other studies. To better understand its chloroplast characteristics and adaptive evolution, chloroplast genomes of 61 Medicago ruthenica were assembled (including 16 cultivated Medicago ruthenica germplasm and 45 wild Medicago ruthenica germplasm). These were used to construct the pan-chloroplast genome of Medicago ruthenica, and the chloroplast genomes of cultivated and wild Medicago ruthenica were compared and analyzed. Phylogenetic and haplotype analyses revealed two main clades of 61 Medicago ruthenica germplasm chloroplast genomes, distributed in eastern and western regions. Meanwhile, based on chloroplast variation information, 61 Medicago ruthenica germplasm can be divided into three genetic groups. Unlike the phylogenetic tree constructed from the chloroplast genome, a new intermediate group has been identified, mainly consisting of samples from the eastern region of Inner Mongolia, Shanxi Province, and Hebei Province. Transcriptomic analysis showed that 29 genes were upregulated and three genes were downregulated. The analysis of these genes mainly focuses on enhancing plant resilience and adapting adversity by stabilizing the photosystem structure and promoting protein synthesis. Additionally, in the analysis of adaptive evolution, the accD, clpP and ycf1 genes showed higher average Ka/Ks ratios and exhibited significant nucleotide diversity, indicating that these genes are strongly positively selected. The editing efficiency of the ycf1 and clpP genes significantly increases under abiotic stress, which may positively contribute to plant adaptation to the environment. In conclusion, the construction and comparative analysis of the complete chloroplast genomes of 61 Medicago ruthenica germplasm from different regions not only revealed new insights into the genetic variation and phylogenetic relationships of Medicago ruthenica germplasm, but also highlighted the importance of chloroplast transcriptome analysis in elucidating the model of chloroplast responses to abiotic stress. These provide valuable information for further research on the adaptive evolution of Medicago ruthenica.

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.

White clover (Trifolium repens L.) is an allopolyploid plant and an excellent perennial legume forage. However, white clover is subjected to various stresses during its growth, with cold stress being one of the major limiting factors affecting its growth and development. Beta-amylase (BAM) is an important starch-hydrolyzing enzyme that plays a significant role in starch degradation and responses to environmental stress. In this study, 21 members of the BAM gene family were identified in the white clover genome. A phylogenetic analysis using BAMs from Arabidopsis divided TrBAMs into four groups based on sequence similarity. Through analysis of conserved motifs, gene duplication, synteny analysis, and cis-acting elements, a deeper understanding of the structure and evolution of TrBAMs in white clover was gained. Additionally, a gene regulatory network (GRN) containing TrBAMs was constructed; gene ontology (GO) annotation analysis revealed close interactions between TrBAMs and AMY (α-amylase) and DPE (4-alpha-glucanotransferase). To determine the function of TrBAMs under various tissues and stresses, RNA-seq datasets were analyzed, showing that most TrBAMs were significantly upregulated in response to biotic and abiotic stresses and the highest expression in leaves. These results were validated through qRT-PCR experiments, indicating their involvement in multiple gene regulatory pathways responding to cold stress. This study provides new insights into the structure, evolution, and function of the white clover BAM gene family, laying the foundation for further exploration of the functional mechanisms through which TrBAMs respond to cold stress.

Plants are frequently exposed to various abiotic stresses, among which low-temperature stress markedly impairs growth and physiological functions. WRKY transcription factors are key regulators in plant responses to abiotic stress. In this study, a novel WRKY transcription factor gene, TrWRKY79, was cloned from white clover. Functional characterization revealed that the full-length TrWRKY79 protein possesses typical features of transcription factors, including transcriptional activation activity located at its C-terminal domain. Heterologous expression of TrWRKY79 in Arabidopsis thaliana significantly enhanced cold tolerance under low-temperature stress. Physiological assays showed that the transgenic lines exhibited higher chlorophyll content and elevated activities of superoxide dismutase (SOD), peroxidase (POD), and catalase (CAT) compared to wild-type plants. Furthermore, Protenix was employed to predict the potential target genes of TrWRKY transcription factors, and their expression profiles were analyzed to help elucidate the regulatory network underlying cold tolerance. qRT-PCR analysis confirmed that several cold-responsive genes, such as COR47 and ABI5, were significantly upregulated in the transgenic lines. Collectively, these findings indicate that TrWRKY79 plays a positive regulatory role in enhancing cold tolerance, providing valuable insights into the molecular mechanisms of cold resistance in white clover and offering promising candidate genes for improving stress resilience in forage crops.

The brown frog (Rana dybowskii) exhibits a wide distribution and is extensively cultured in northeast China. This species undergoes a prolonged period of hibernation lasting from several months to half a year. The frog’s fallopian tubes are considered a valuable tonic medicine known as “Oviducts Ranae” in traditional Chinese medicine. To enhance comprehension of the molecular mechanism underlying the process of emergence from hibernation, a transcriptome-based comparative analysis was performed on differentially expressed genes (DEGs) across various organs of female frogs during hibernation and upon emergence. The organs examined encompassed the brain, liver, spleen, fallopian tube, and ovary. Subsequently, GO and KEGG enrichment analyses were performed to gain further insights into these DEGs. A total of 51,634 transcripts were identified in all five tissues. The spleen exhibited the highest number of DEGs, with 3651 members, while the liver, brain, and fallopian tube had 3182, 3115, and 3186 DEGs, respectively. The ovary had the lowest number of DEGs, with only 1661. Interestingly, only 45 genes were found to be differentially expressed in all five tissues, and GO enrichment analysis revealed distinct functional differences among the DEGs in the various tissues. Only two meaningful DEG-enriched KEGG pathways, “00310 Lysine degradation” and “05202 Transcriptional misregulation in cancer”, were present in all five tissues, and the pathway “B 09182 Protein families: genetic information processing” was significantly enriched in four tissues except the ovary, and it had the most enriched DEGs. Our findings offer insights to grasp the factors that need to be controlled in the design of overwintering sites and offer a novel perspective for the conservation and management of the overwintering population of female R. dybowskii.