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.

Iron (Fe) is necessary for plant growth and development. The mechanism of uptake and translocation in Cadmium (Cd) is similar to iron, which shares iron transporters. Yellow stripe-like transporter (YSL) plays a pivotal role in transporting iron and other metal ions in plants. In this study, MsYSL6 and its promoter were cloned from leguminous forage alfalfa. The transient expression of MsYSL6-GFP indicated that MsYSL6 was localized to the plasma membrane and cytoplasm. The expression of MsYSL6 was induced in alfalfa by iron deficiency and Cd stress, which was further proved by GUS activity driven by the MsYSL6 promoter. To further identify the function of MsYSL6, it was heterologously overexpressed in tobacco. MsYSL6-overexpressed tobacco showed better growth and less oxidative damage than WT under Cd stress. MsYSL6 overexpression elevated Fe and Cd contents and induced a relatively high Fe translocation rate in tobacco under Cd stress. The results suggest that MsYSL6 might have a dual function in the absorption of Fe and Cd, playing a role in the competitive absorption between Fe and Cd. MsYSL6 might be a regulatory factor in plants to counter Cd stress. This study provides a novel gene for application in heavy metal enrichment or phytoremediation and new insights into plant tolerance to toxic metals.

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.

Gametocidal (Gc) chromosomes can ensure their preferential transmission by killing the gametes without themselves through causing chromosome breakage and therefore have been exploited as an effective tool for genetic breeding. However, to date very little is known about the molecular mechanism of Gc action. In this study, we used methylation-sensitive amplified polymorphism (MSAP) technique to assess the extent and pattern of cytosine methylation alterations at the whole genome level between two lines of wheat Gc addition line and their common wheat parent. The results indicated that the overall levels of cytosine methylation of two studied Gc addition lines (CS–3C and CS–3C3C, 48.68% and 48.65%, respectively) were significantly increased when compared to common wheat CS (41.31%) and no matter fully methylated or hemimethylated rates enhanced in Gc addition lines. A set of 30 isolated fragments that showed different DNA methylation or demethylation patterns between the three lines were sequenced and the results indicated that 8 fragments showed significant homology to known sequences, of which three were homologous to MITE transposon (Miniature inverted–repeat transposable elements), LTR-retrotransposon WIS-1p and retrotransposon Gypsy, respectively. Overall, our results showed that DNA methylation could play a role in the Gc action.

Polyploidization is a driving force in plant evolution. Chromosomal variation often occurs at early generations following polyploid formation due to meiotic pairing irregularity that may compromise segregation fidelity and cause homoeologous exchange (HE). The trends of chromosomal variation and especially factors affecting HE remain to be fully deciphered. Here, by whole-genome resequencing, we performed nuanced analyses of patterns of chromosomal number variation and explored genomic features that affect HE in two early generations of a synthetic rice segmental allotetraploid. We found a wide occurrence of whole-chromosome aneuploidy and, to a lesser extent, also large segment gains/losses in both generations (S2 and S4) of the tetraploids. However, while the number of chromosome gains was similar between S2 and S4, that of losses in S4 was lower than in S2. HEs were abundant across all chromosomes in both generations and showed variable correlations with different genomic features at chromosomal and/or local scales. Contents of genes and transposable elements (TEs) were positively and negatively correlated with HE frequencies, respectively. By dissecting TEs into different classes, retrotransposons were found to be negatively correlated with HE frequency to a stronger extent than DNA transposons, whereas miniature terminal inverted elements (MITEs) showed a strong positive correlation. Local HE frequencies in the tetraploids and homologous recombination (HR) rates in diploids within 1 Mb sliding windows were significantly correlated with each other and showed similar overall distribution profiles. Nonetheless, non-concordant trends between HE and HR rates were found at distal regions in some chromosomes. At local scale, both shared and polymorphic retrotransposons between parents were negatively correlated with HE frequency; in contrast, both shared and polymorphic MITEs showed positive correlations with HE frequency. Our results shed new light on the patterns of chromosomal number variation and reveal genomic features influencing HE frequency in early generations following plant polyploidization.