Alfalfa is the most widely cultivated forage legume and one of the most economically valuable crops in the world. The large size and complexity of the alfalfa genome has delayed the development of genomic resources for alfalfa research. Second-generation Illumina transcriptome sequencing is an efficient method for generating a global transcriptome sequence dataset for gene discovery and molecular marker development in alfalfa. More than 28 million sequencing reads (5.64 Gb of clean nucleotides) were generated by Illumina paired-end sequencing from 15 different alfalfa tissue samples. In total, 40,433 unigenes with an average length of 803 bp were obtained by de novo assembly. Based on a sequence similarity search of known proteins, a total of 36,684 (90.73%) unigenes were annotated. In addition, 1,649 potential EST-SSRs were identified as potential molecular markers from unigenes with lengths exceeding 1 kb. A total of 100 pairs of PCR primers were randomly selected to validate the assembly quality and develop EST-SSR markers from genomic DNA. Of these primer pairs, 82 were able to amplify sequences in initial screening tests, and 27 primer pairs successfully amplified DNA fragments and detected significant amounts of polymorphism among 10 alfalfa accessions. The present study provided global sequence data for autotetraploid alfalfa and demonstrates the Illumina platform is a fast and effective approach to EST-SSR markers development in alfalfa. The use of these transcriptome datasets will serve as a valuable public information platform to accelerate studies of the alfalfa genome.

Alfalfa is the most extensively cultivated forage legume. Salinity is a major environmental factor that impacts on alfalfa’s productivity. However, little is known about the molecular mechanisms underlying alfalfa responses to salinity, especially the relative contribution of the two important components of osmotic and ionic stress. In this study, we constructed the first full-length transcriptome database for alfalfa root tips under continuous NaCl and mannitol treatments for 1, 3, 6, 12, and 24 h (three biological replicates for each time points, including the control group) via PacBio Iso-Seq. This resulted in the identification of 52,787 full-length transcripts, with an average length of 2551 bp. Global transcriptional changes in the same 33 stressed samples were then analyzed via BGISEQ-500 RNA-Seq. Totals of 8861 NaCl-regulated and 8016 mannitol-regulated differentially expressed genes (DEGs) were identified. Metabolic analyses revealed that these DEGs overlapped or diverged in the cascades of molecular networks involved in signal perception, signal transduction, transcriptional regulation, and antioxidative defense. Notably, several well characterized signalling pathways, such as CDPK, MAPK, CIPK, and PYL-PP2C-SnRK2, were shown to be involved in osmotic stress, while the SOS core pathway was activated by ionic stress. Moreover, the physiological shifts of catalase and peroxidase activity, glutathione and proline content were in accordance with dynamic transcript profiles of the relevant genes, indicating that antioxidative defense system plays critical roles in response to salinity stress. Overall, our study provides evidence that the response to salinity stress in alfalfa includes both osmotic and ionic components. The key osmotic and ionic stress-related genes are candidates for future studies as potential targets to improve resistance to salinity stress via genetic engineering.

Background: Alfalfa, the “queen of forage”, is the most extensively cultivated forage legume in the world. The development and yield of alfalfa are seriously limited by abiotic stress. MADS-box transcription factors are one of the largest gene families and play a pivotal role in plant development and abiotic stress. However, little is known regarding the MADS-box transcription factors in autotetraploid cultivated alfalfa. Results: In the present study, we identified 120 MsMADS-box genes in the alfalfa genome. Phylogenetic analysis indicated that 75 type-I MsMADS-box genes were classified into the Mα, Mβ, and Mγ subgroups, and 45 type-II MsMADS-box genes were classified into 11 subgroups. The promoter region of MsMADS-box genes containing several hormone and stress related elements. Chromosomal location analysis revealed that 117 MsMADS-box genes were unevenly distributed on 32 chromosomes, and the remaining three genes were located on unmapped scaffolds. A total of nine pairs of segmental duplications and four groups of tandem duplications were found. Expression analysis showed that MsMADS-box genes were differentially expressed in various tissues and under abiotic stresses. qRT-PCR analysis revealed that the expression profiles of eight selected MsMADS-box genes were distinct under various stresses. Conclusions: In this study, MsMADS-box genes were identified in the cultivated alfalfa genome based on autotetraploid level, and further confirmed by Gene Ontology (GO) analysis, phylogenetic analysis, sequence features and expression analysis. Taken together, these findings will provide clues for further study of MsMADS-box functions and alfalfa molecular breeding. Our study is the first to systematically identify and characterize the MADS-box transcription factors in autotetraploid cultivated alfalfa (Medicago sativa L.), and eight MsMADS-box genes were significantly involved in response to various stresses.

Siberian wildrye (Elymus sibiricus L.) is a perennial, self-fertilizing grass that plays an important role in animal husbandry and environmental sustenance. However, the transcriptomic and genomic information on this species is very limited, which hinders genetic and breeding studies. In the present study, 76,686,804 clean reads were generated from 11 different tissue samples of E. sibiricus by Illumina paired-end sequencing and the reads were deposited into the NCBI SRA database (SRX574376). A total of 8,769 EST-SSRs were identified from 94,458 unigene sequences, which were obtained by de novo assembly. Moreover, 1,078 primer pairs were successfully designed and 500 pairs were randomly selected to assess polymorphisms in 15 E. sibiricus accessions. A total of 112 primer pairs were polymorphic and the polymorphism information content (PIC) values ranged from 0.39 to 0.81, indicating a high level of informativeness. Furthermore, these 112 polymorphic primer pairs were used to evaluate the transferability to 13 other related species and 55 EST-SSR markers were found to be polymorphic among these 13 Elymus species. This study collected the global sequence data for E. sibiricus, and the newly developed markers will prove valuable in facilitating genetic diversity in E. sibiricus and related Elymus species.

The seed dispersion caused by pod shattering is a form of propagation used by many wild species. Loss of seeds from pod shattering is frequent in the common vetch (Vicia sativa L.), an important self-pollinating annual forage legume. However, pod shattering is one of the most important defects that limits the reproduction of the vetch in the field and the usage as a leguminous forage crop. To better understand the vetch pod shattering mechanism, we used high-throughput RNA sequencing to assess the global changes in the transcriptomes of the pod ventral sutures of shattering-susceptible and shattering-resistant vetch accessions screened from 541 vetch germplasms. A total of 1,285 significantly differentially expressed unigenes (DEGs) were detected, including 575 up-regulated unigenes and 710 down-regulated unigenes. Analyses of Gene Ontology and KEGG metabolic enrichment pathways of 1,285 DEGs indicated that 22 DEGs encoding cell wall modifications and hydrolases associated with pod shattering were highly expressed in shattering-susceptible accessions. These genes were mainly enriched in "hydrolase activity", "cytoplasm" and "carbohydrate metabolic process" systems. These cell wall modifications and hydrolases genes included β-glucosidase and endo-polygalacturonase, which work together to break down the glycosidic bonds of pectin and cellulose, and to promote the dissolution and disappearance of the cell wall in the ventral suture of the pod and make the pod more susceptible to shattering. We demonstrated the differences in gene transcription levels between the shattering-susceptible and shattering-resistant vetch accessions for the first time and our results provided valuable information for the identifying and characterizing of pod shattering regulation networks in vetch. This information may facilitate the future identification of pod shattering-related genes and their underlying molecular mechanisms in the common vetch.