RNA-seq analyses on gametogenic tissues of alfalfa (Medicago sativa) revealed plant reproduction- and ploidy-related genes

Background In alfalfa (Medicago sativa), the coexistence of interfertile subspecies (i.e. sativa, falcata and coerulea) characterized by different ploidy levels (diploidy and tetraploidy) and the occurrence of meiotic mutants capable of producing unreduced (2n) gametes, have been efficiently combined for the establishment of new polyploids. The wealth of agronomic data concerning forage quality and yield provides a thorough insight into the practical benefits of polyploidization. However, many of the underlying molecular mechanisms regarding gene expression and regulation remained completely unexplored. In this study, we aimed to address this gap by examining the transcriptome profiles of leaves and reproductive tissues, corresponding to anthers and pistils, sampled at different time points from diploid and tetraploid Medicago sativa individuals belonging to progenies produced by bilateral sexual polyploidization (dBSP and tBSP, respectively) and tetraploid individuals stemmed from unilateral sexual polyploidization (tUSP). Results Considering the crucial role played by anthers and pistils in the reduced and unreduced gametes formation, we firstly analyzed the transcriptional profiles of the reproductive tissues at different stages, regardless of the ploidy level and the origin of the samples. By using and combining three different analytical methodologies, namely weighted-gene co-expression network analysis (WGCNA), tau (τ) analysis, and differentially expressed genes (DEGs) analysis, we identified a robust set of genes and transcription factors potentially involved in both male sporogenesis and gametogenesis processes, particularly in crossing-over, callose synthesis, and exine formation. Subsequently, we assessed at the same floral stage, the differences attributable to the ploidy level (tBSP vs. dBSP) or the origin (tBSP vs. tUSP) of the samples, leading to the identification of ploidy and parent-specific genes. In this way, we identified, for example, genes that are specifically upregulated and downregulated in flower buds in the comparison between tBSP and dBSP, which could explain the reduced fertility of the former compared to the latter materials. Conclusions While this study primarily functions as an extensive investigation at the transcriptomic level, the data provided could represent not only a valuable original asset for the scientific community but also a fully exploitable genomic resource for functional analyses in alfalfa.