De novo characterization of fall dormant and nondormant alfalfa (Medicago sativa L.) leaf transcriptome and identification of candidate genes related to fall dormancy

Identification of Cold Tolerance Transcriptional Regulatory Genes in Seedlings of Medicago sativa L. and Medicago falcata L

Alfalfa species Medicago sativa L. (MS) and Medicago falcata L. (MF), globally prominent perennial leguminous forages, hold substantial economic value. However, our comprehension of the molecular mechanisms governing their resistance to cold stress remains limited. To address this knowledge gap, we scrutinized and compared MS and MF cold-stress responses at the molecular level following 24 h and 120 h low-temperature exposure (4 °C). Our study revealed that MF had superior physiological resilience to cold stress compared with MS, and its morphology was healthier under cold stress, and its malondialdehyde content and superoxide dismutase activity increased, first, and then decreased, while the soluble sugar content continued to accumulate. Transcriptome analysis showed that after 120 h of exposure, there were different gene-expression patterns between MS and MF, including 1274 and 2983 genes that were continuously up-regulated, respectively, and a total of 923 genes were included, including star cold-resistant genes such as ICE1 and SIP1. Gene ontology (GO) enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses revealed numerous inter-species differences in sustained cold-stress responses. Notably, MS-exclusive genes included a single transcription factor (TF) gene and several genes associated with a single DNA repair-related pathway, whereas MF-exclusive genes comprised nine TF genes and genes associated with 14 pathways. Both species exhibited high-level expression of genes encoding TFs belonging to AP2-EREBP, ARR-B, and bHLH TF families, indicating their potential roles in sustaining cold resistance in alfalfa-related species. These findings provide insights into the molecular mechanisms governing cold-stress responses in MS and MF, which could inform breeding programs aimed at enhancing cold-stress resistance in alfalfa cultivars.

MODMS: a multi-omics database for facilitating biological studies on alfalfa (Medicago sativa L.)

Alfalfa (Medicago sativa L.) is a globally important forage crop. It also serves as a vegetable and medicinal herb because of its excellent nutritional quality and significant economic value. Multi-omics data on alfalfa continue to accumulate owing to recent advances in high-throughput techniques, and integrating this information holds great potential for expediting genetic research and facilitating advances in alfalfa agronomic traits. Therefore, we developed a comprehensive database named MODMS (multi-omics database of M. sativa) that incorporates multiple reference genomes, annotations, comparative genomics, transcriptomes, high-quality genomic variants, proteomics, and metabolomics. This report describes our continuously evolving database, which provides researchers with several convenient tools and extensive omics data resources, facilitating the expansion of alfalfa research. Further details regarding the MODMS database are available at https://modms.lzu.edu.cn/.

Identification of genetic loci associated with five agronomic traits in alfalfa using multi-environment trials

The use of multi-environment trials to test yield-related traits in a diverse alfalfa panel allowed to find multiple molecular markers associated with complex agronomic traits. Yield is one of the most important target traits in alfalfa breeding; however, yield is a complex trait affected by genetic and environmental factors. In this study, we used multi-environment trials to test yield-related traits in a diverse panel composed of 200 alfalfa accessions and varieties. Phenotypic data of maturity stage measured as mean stage by count (MSC), dry matter content, plant height (PH), biomass yield (Yi), and fall dormancy (FD) were collected in three locations in Idaho, Oregon, and Washington from 2018 to 2020. Single-trial and stagewise analyses were used to obtain estimated trait means of entries by environment. The plants were genotyped using a genotyping by sequencing approach and obtained a genotypic matrix with 97,345 single nucleotide polymorphisms. Genome-wide association studies identified a total of 84 markers associated with the traits analyzed. Of those, 29 markers were in noncoding regions and 55 markers were in coding regions. Ten significant SNPs at the same locus were associated with FD and they were linked to a gene annotated as a nuclear fusion defective 4-like (NFD4). Additional SNPs associated with MSC, PH, and Yi were annotated as transcription factors such as Cysteine3Histidine (C3H), Hap3/NF-YB family, and serine/threonine-protein phosphatase 7 proteins, respectively. Our results provide insight into the genetic factors that influence alfalfa maturity, yield, and dormancy, which is helpful to speed up the genetic gain toward alfalfa yield improvement.

Transcriptomics and metabolomics revealed that phosphate improves the cold tolerance of alfalfa

INTRODUCTION

Alfalfa (Medicago sativa L.) is a highly nutritious leguminous forage that plays an essential role in animal husbandry. In the middle and high latitudes of the northern hemisphere, there are problems with its low rates of overwintering and production. The application of phosphate (P) is an important measure to improve the cold resistance and production of alfalfa, but little is known about the mechanism of P in improving the cold resistance of alfalfa.

METHODS

This study integrated the transcriptome and metabolome to explain the mechanism of alfalfa in response to low-temperature stress under two applications of P (50 and 200 mg kg^-1) and a control of none applied.

RESULTS

The application of P fertilizer improved the root structure and increased the content of soluble sugar and soluble protein in the root crown. In addition, there were 49 differentially expressed genes (DEGs) with 23 upregulated and 24 metabolites with 12 upregulated when 50 mg kg^-1 of P was applied. In contrast, there were 224 DEGs with 173 upregulated and 12 metabolites with 6 upregulated in the plants treated with 200 mg kg^-1 of P compared with the Control Check (CK). These genes and metabolites were significantly enriched in the biosynthesis of other secondary metabolites and the metabolic pathways of carbohydrates and amino acids. The integration of the transcriptome and metabolome indicated that P affected the biosynthesis of N-acetyl-L-phenylalanine, L-serine, lactose, and isocitrate during the period of increasing cold. It could also affect the expression of related genes that regulate cold tolerance in alfalfa.

DISCUSSION

Our findings could contribute to a deeper understanding of the mechanism that alfalfa uses to tolerate cold and lay a theoretical foundation for breeding alfalfa that is highly efficient at utilizing phosphorus.

A Pan-Transcriptome Analysis Indicates Efficient Downregulation of the FIB Genes Plays a Critical Role in the Response of Alfalfa to Cold Stress

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.

Comparative Transcriptome and Anatomic Characteristics of Stems in Two Alfalfa Genotypes

Stems are more important to forage quality than leaves in alfalfa. To understand lignin formation at different stages in alfalfa, lignin distribution, anatomical characteristics and transcriptome profile were employed using two alfalfa cultivars. The results showed that the in vitro true digestibility (IVTD) of stems in WL168 was significantly higher than that of Zhungeer, along with the significantly lower neutral detergent fiber (NDF), acid detergent fiber (ADF) and lignin contents. In addition, Zhungeer exhibited increased staining of the xylem areas in the stems of different developmental stages compared to WL168. Interestingly, the stems of WL168 appeared intracellular space from the stage 3, while Zhungeer did not. The comparative transcriptome analysis showed that a total of 1993 genes were differentially expressed in the stem between the cultivars, with a higher number of expressed genes in the stage 4. Of the differentially expressed genes, starch and sucrose metabolism as well as phenylpropanoid biosynthesis pathways were the most significantly enriched pathways. Furthermore, expression of genes involved in lignin biosynthesis such as PAL, 4CL, HCT, CAD, COMT and POD coincides with the anatomic characteristics and lignin accumulation. These results may help elucidate the regulatory mechanisms of lignin biosynthesis and improve forage quality in alfalfa.

Genomic prediction for canopy height and dry matter yield in alfalfa using family bulks

Genomic selection (GS) has proven to be an effective method to increase genetic gain rates and accelerate breeding cycles in many crop species. However, its implementation requires large investments to phenotype of the training population and for routine genotyping. Alfalfa (Medicago sativa L.) is one of the major cultivated forage legumes, showing high-quality nutritional value. Alfalfa breeding is usually carried out by phenotypic recurrent selection and is commonly done at the family level. The application of GS in alfalfa could be simplified and less costly by genotyping and phenotyping families in bulks. For this study, an alfalfa reference population composed of 142 full-sib and 35 half-sib families was bulk-genotyped using target enrichment sequencing and phenotyped for dry matter yield (DMY) and canopy height (CH) in Florida, USA. Genotyping of the family bulks with 17,707 targeted probes resulted in 114,945 single-nucleotide polymorphisms. The markers revealed a population structure that matched the mating design, and the linkage disequilibrium slowly decayed in this breeding population. After exploring multiple prediction scenarios, a strategy was proposed including data from multiple harvests and accounting for the G×E in the training population, which led to a higher predictive ability of up to 38 and 24% for DMY and CH, respectively. Although this study focused on the implementation of GS in alfalfa families, the bulk methodology and the prediction schemes used herein could guide future studies in alfalfa and other crops bred in bulks.

Identification and Characterization of Regulatory Pathways Controlling Dormancy Under Lower Temperature in Alfalfa (Medicago sativa L.)

Alfalfa (Medicago sativa L.), a kind of high-quality perennial legume forage, is widely distributed in the northern regions of China. In recent years, low temperatures have frequently occurred and limited alfalfa productivity and survival in early spring and late fall. However, the underlying molecular mechanisms of alfalfa response to cold tolerance are not well-documented. In this study, dormancy and non-dormancy alfalfa standard varieties were characterized under low-temperature stress. Our analysis revealed that plant height of the dormancy genotype was strongly inhibited by low temperature; flavonoids content, and higher expression of flavonoids biosynthesis genes (chalcone synthase, leucoanthocyanidin dioxygenase, and flavonoid 3'-monooxygenase) may play essential roles in response to low-temperature stress in dormancy genotype alfalfa. Further analyses revealed that receptor-like kinase family genes (such as cysteine-rich RLK10, lectin protein kinase, and S-locus glycoprotein like kinase), RNA and protein synthesis genes (RNA polymerases, ribosomal protein, and protein phosphatase 2C family protein), and proteasome degradation pathway genes (such as F-box family protein, RING/U-box superfamily protein, and zinc finger family protein) also highly upregulated and contributed to cold tolerance phenotype in dormancy genotype alfalfa. This will provide new insights into future studies for cold tolerance in alfalfa and offer new target genes for further functional characterization and genetic improvement of alfalfa.

Reference nodule transcriptomes for Melilotus officinalis and Medicago sativa cv. Algonquin

Host/symbiont compatibility is a hallmark of the symbiotic nitrogen-fixing interaction between rhizobia and legumes, mediated in part by plant-produced nodule-specific cysteine-rich (NCR) peptides and the bacterial BacA membrane protein that can act as a NCR peptide transporter. In addition, the genetic and metabolic properties supporting symbiotic nitrogen fixation often differ between compatible partners, including those sharing a common partner, highlighting the need for multiple study systems. Here, we report high-quality nodule transcriptome assemblies for Medicago sativa cv. Algonquin and Melilotus officinalis, two legumes able to form compatible symbioses with Sinorhizobium meliloti. The compressed M. sativa and M. officinalis assemblies consisted of 79,978 and 64,593 contigs, respectively, of which 33,341 and 28,278 were assigned putative annotations, respectively. As expected, the two transcriptomes showed broad similarity at a global level. We were particularly interested in the NCR peptide profiles of these plants, as these peptides drive bacterial differentiation during the symbiosis. A total of 412 and 308 NCR peptides were predicted from the M. sativa and M. officinalis transcriptomes, respectively, with approximately 9% of the transcriptome of both species consisting of NCR transcripts. Notably, transcripts encoding highly cationic NCR peptides (isoelectric point > 9.5), which are known to have antimicrobial properties, were ∼2-fold more abundant in M. sativa than in M. officinalis, and ∼27-fold more abundant when considering only NCR peptides in the six-cysteine class. We hypothesize that the difference in abundance of highly cationic NCR peptides explains our previous observation that some rhizobial bacA alleles which can support symbiosis with M. officinalis are unable to support symbiosis with M. sativa.

A Review of Unreduced Gametes and Neopolyploids in Alfalfa: How to Fill the Gap between Well-Established Meiotic Mutants and Next-Generation Genomic Resources

The gene flow mediated by unreduced gametes between diploid and tetraploid plants of the Medicagosativa-coerulea-falcata complex is pivotal for alfalfa breeding. Sexually tetraploidized hybrids could represent the best way to exploit progressive heterosis simultaneously derived from gene diversity, heterozygosity, and polyploidy. Moreover, unreduced gametes combined with parthenogenesis (i.e., apomixis) would enable the cloning of plants through seeds, providing a unique opportunity for the selection of superior genotypes with permanently fixed heterosis. This reproductive strategy has never been detected in the genus Medicago, but features of apomixis, such as restitutional apomeiosis and haploid parthenogenesis, have been reported. By means of an original case study, we demonstrated that sexually tetraploidized plants maintain apomeiosis, but this trait is developmentally independent from parthenogenesis. Alfalfa meiotic mutants producing unreduced egg cells revealed a null or very low capacity for parthenogenesis. The overall achievements reached so far are reviewed and discussed along with the efforts and strategies made for exploiting reproductive mutants that express apomictic elements in alfalfa breeding programs. Although several studies have investigated the cytological mechanisms responsible for 2n gamete formation and the inheritance of this trait, only a very small number of molecular markers and candidate genes putatively linked to unreduced gamete formation have been identified. Furthermore, this scenario has remained almost unchanged over the last two decades. Here, we propose a reverse genetics approach, by exploiting the genomic and transcriptomic resources available in alfalfa. Through a comparison with 9 proteins belonging to Arabidopsis thaliana known for their involvement in 2n gamete production, we identified 47 orthologous genes and evaluated their expression in several tissues, paving the way for novel candidate gene characterization studies. An overall view on strategies suitable to fill the gap between well-established meiotic mutants and next-generation genomic resources is presented and discussed.

Phenotypic Diversity and Productivity of Medicago sativa Subspecies from Drought-Prone Environments in Mediterranean Type Climates

The phenotypic diversity and productivity of a diverse alfalfa (M. sativa subspp.) panel of cultivars, landraces and wild relatives with putative drought tolerance were evaluated in two Mediterranean environments (central Chile and Southern Australia). In Chile, 70 accessions were evaluated in rainfed conditions and in Australia 30 accessions under rainfed and irrigated conditions, during three growing seasons. Large phenotypic variation was observed among and within subspecies for NDVI, stem length, intercepted PAR and forage yield. Principal component analysis indicated that the first two principal components (PC) accounted for 84.2% of total variance; fall dormancy, taxa, and breeding status were closely related to the agronomical performance of alfalfa accessions. Forage yield varied largely among accessions across years and locations. A linear relationship was found between annual forage yield and annual water added to the experiments (R2 = 0.60, p < 0.001). The GxE analysis for forage yield allowed the detection of the highest yielding accessions for each of the two mega-environments identified. The accessions CTA002 and CTA003 showed greater forage yield in both Chile and Australia environments. It is concluded that new breeding lines derived from crosses between cultivated alfalfa (M. sativa subsp. sativa) and wild relatives belonging to the primary (M. sativa subsp. falcata) and tertiary (M. arborea) gene pool, achieve outstanding agronomical performance in drought-prone environments.

High-Throughput Sequencing for Deciphering the Virome of Alfalfa (Medicago sativa L.)

Alfalfa (Medicago sativa L.), also known as lucerne, is a major forage crop worldwide. In the United States, it has recently become the third most valuable field crop, with an estimated value of over $9.3 billion. Alfalfa is naturally infected by many different pathogens, including viruses, obligate parasites that reproduce only inside living host cells. Traditionally, viral infections of alfalfa have been considered by breeders, growers, producers and researchers to be diseases of limited importance, although they are widespread in all major cultivation areas. However, over the past few years, due to the rapid development of high-throughput sequencing (HTS), viral metagenomics, bioinformatics tools for interpreting massive amounts of HTS data and the increasing accessibility of public data repositories for transcriptomic discoveries, several emerging viruses of alfalfa with the potential to cause serious yield losses have been described. They include alfalfa leaf curl virus (family Geminiviridae), alfalfa dwarf virus (family Rhabdoviridae), alfalfa enamovirus 1 (family Luteoviridae), alfalfa virus S (family Alphaflexiviridae) and others. These discoveries have called into question the assumed low economic impact of viral diseases in alfalfa and further suggested their possible contribution to the severity of complex infections involving multiple pathogens. In this review, we will focus on viruses of alfalfa recently described in different laboratories on the basis of the above research methodologies.

Biotechnological Perspectives of Omics and Genetic Engineering Methods in Alfalfa

For several decades, researchers are working to develop improved major crops with better adaptability and tolerance to environmental stresses. Forage legumes have been widely spread in the world due to their great ecological and economic values. Abiotic and biotic stresses are main factors limiting legume production, however, alfalfa (Medicago sativa L.) shows relatively high level of tolerance to drought and salt stress. Efforts focused on alfalfa improvements have led to the release of cultivars with new traits of agronomic importance such as high yield, better stress tolerance or forage quality. Alfalfa has very high nutritional value due to its efficient symbiotic association with nitrogen-fixing bacteria, while deep root system can help to prevent soil water loss in dry lands. The use of modern biotechnology tools is challenging in alfalfa since full genome, unlike to its close relative barrel medic (Medicago truncatula Gaertn.), was not released yet. Identification, isolation, and improvement of genes involved in abiotic or biotic stress response significantly contributed to the progress of our understanding how crop plants cope with these environmental challenges. In this review, we provide an overview of the progress that has been made in high-throughput sequencing, characterization of genes for abiotic or biotic stress tolerance, gene editing, as well as proteomic and metabolomics techniques bearing biotechnological potential for alfalfa improvement.

Identification of emerging viral genomes in transcriptomic datasets of alfalfa (Medicago sativa L.)

BACKGROUND

Publicly available transcriptomic datasets have become a valuable tool for the discovery of new pathogens, particularly viruses. In this study, several coding-complete viral genomes previously not found or experimentally confirmed in alfalfa were identified in the plant datasets retrieved from the NCBI Sequence Read Archive.

METHODS

Publicly available Medicago spp. transcriptomic datasets were retrieved from the NCBI SRA database. The raw reads were first mapped to the reference genomes of Medicago sativa and Medigago truncatula followed by the alignment of the unmapped reads to the NCBI viral genome database and de novo assembly using the SPAdes tool. When possible, assemblies were experimentally confirmed using 5'/3' RACE and RT-PCRs.

RESULTS

Twenty three different viruses were identified in the analyzed datasets, of which several represented emerging viruses not reported in alfalfa prior to this study. Among them were two strains of cnidium vein yellowing virus, lychnis mottle virus and Cactus virus X, for which coding-complete genomic sequences were obtained by a de novo assembly.

CONCLUSIONS

The results improve our knowledge of the diversity and host range of viruses infecting alfalfa, provide essential tools for their diagnostics and characterization and demonstrate the utility of transcriptomic datasets for the discovery of new pathogens.

How fall dormancy benefits alfalfa winter-survival? Physiologic and transcriptomic analyses of dormancy process

BACKGROUND

Fall dormancy and freezing tolerance characterized as two important phenotypic traits, have great effects on productivity and persistence of alfalfa (Medicago sativa L.). Despite the fact that one of the most limiting traits for alfalfa freezing tolerance in winter is fall dormancy, the interplay between fall dormancy and cold acclimation processes of alfalfa remains largely unknown. We compared the plant regrowth, winter survival, raffinose and amino acids accumulation, and genome-wide differentially expressed genes of fall-dormant cultivar with non-dormant cultivar under cold acclimation.

RESULTS

Averaged over both years, the non-dormant alfalfa exhibited largely rapid regrowth compared with fall dormant alfalfa after last cutting in autumn, but the winter survival rate of fall dormant alfalfa was about 34-fold higher than that of non-dormant alfalfa. The accumulation of raffinose and amino acids were significantly increased in fall dormant alfalfa, whereas were decreased in non-dormant alfalfa under cold acclimation. Expressions of candidate genes encoding raffinose biosynthesis genes were highly up-regulated in fall dormant alfalfa, but down-regulated in non-dormant alfalfa under cold acclimation. In fall dormant alfalfa, there was a significantly down-regulated expression of candidate genes encoding the glutamine synthase, which is indirectly involved in the proline metabolism. A total of eight significantly differentially expressed transcription factors (TFs) related to CBF and ABRE-BFs were identified. The most up-regulated TFs in fall dormant alfalfa cultivar were ABF4 and DREB1C.

CONCLUSIONS

Fall dormant alfalfa drastically increased raffinose and amino acids accumulation under cold acclimation. Raffinose-associated and amino acid-associated genes involved in metabolic pathways were more highly expressed in fall dormant alfalfa than non-dormant alfalfa under cold acclimation. This global survey of transcriptome profiles provides new insights into the interplay between fall dormancy and cold acclimation in alfalfa.

Computational Methods for Mapping, Assembly and Quantification for Coding and Non-coding Transcripts

The measurement of gene expression has long provided significant insight into biological functions. The development of high-throughput short-read sequencing technology has revealed transcriptional complexity at an unprecedented scale, and informed almost all areas of biology. However, as researchers have sought to gather more insights from the data, these new technologies have also increased the computational analysis burden. In this review, we describe typical computational pipelines for RNA-Seq analysis and discuss their strengths and weaknesses for the assembly, quantification and analysis of coding and non-coding RNAs. We also discuss the assembly of transposable elements into transcripts, and the difficulty these repetitive elements pose. In summary, RNA-Seq is a powerful technology that is likely to remain a key asset in the biologist's toolkit.

Analysis of transcripts and splice isoforms in Medicago sativa L. by single-molecule long-read sequencing

The full-length transcriptome of alfalfa was analyzed with PacBio single-molecule long-read sequencing technology. The transcriptome data provided full-length sequences and gene isoforms of transcripts in alfalfa, which will improve genome annotation and enhance our understanding of the gene structure of alfalfa. As an important forage, alfalfa (Medicago sativa L.) is world-wide planted. For its complexity of genome and unfinished whole genome sequencing, the sequences and complete structure of mRNA transcripts remain unclear in alfalfa. In this study, single-molecule long-read sequencing was applied to investigate the alfalfa transcriptome using the Pacific Biosciences platform, and a total of 113,321 transcripts were obtained from young, mature and senescent leaves. We identified 72,606 open reading frames including 46,616 full-length ORFs, 1670 transcription factors from 54 TF families and 44,040 simple sequence repeats from 30,797 sequences. A total of 7568 alternative splicing events was identified and the majority of alternative splicing events in alfalfa was intron retention. In addition, we identified 17,740 long non-coding RNAs. Our results show the feasibility of deep sequencing full-length RNA from alfalfa transcriptome on a single-molecule level.

First Report of Alphapartitiviruses Infecting Alfalfa (Medicago sativa L.) in the United States

Evidence is presented demonstrating that alfalfa cultivars in the United States could be widely infected with cryptic virus of the genus Alphapartitivirus. The nucleotide sequences of several U.S. isolates were obtained. The biological significance or negative effects of the virus on alfalfa are unknown and require further investigation.

Multiple Regulatory Networks Are Activated during Cold Stress in Medicago sativa L

Cultivated alfalfa (Medicago sativa L.) is one of the most important perennial legume forages in the world, and it has considerable potential as a valuable forage crop for livestock. However, the molecular mechanisms underlying alfalfa responses to cold stress are largely unknown. In this study, the transcriptome changes in alfalfa under cold stress at 4 °C for 2, 6, 24, and 48 h (three replicates for each time point) were analyzed using the high-throughput sequencing platform, BGISEQ-500, resulting in the identification of 50,809 annotated unigenes and 5283 differentially expressed genes (DEGs). Metabolic pathway enrichment analysis demonstrated that the DEGs were involved in carbohydrate metabolism, photosynthesis, plant hormone signal transduction, and the biosynthesis of amino acids. Moreover, the physiological changes of glutathione and proline content, catalase, and peroxidase activity were in accordance with dynamic transcript profiles of the relevant genes. Additionally, some transcription factors might play important roles in the alfalfa response to cold stress, as determined by the expression pattern of the related genes during 48 h of cold stress treatment. These findings provide valuable information for identifying and characterizing important components in the cold signaling network in alfalfa and enhancing the understanding of the molecular mechanisms underlying alfalfa responses to cold stress.

Characterization of alfalfa virus F, a new member of the genus Marafivirus

Viral infections of alfalfa are widespread in major cultivation areas and their impact on alfalfa production may be underestimated. A new viral species, provisionally named alfalfa virus F (AVF), was identified using a virion-associated nucleic acid (VANA) metagenomics-based approach in alfalfa (Medicago sativa L.) samples collected in Southern France. The nucleotide sequence of the viral genome was determined by de-novo assembly of VANA reads and by 5'/3' RACE with viral RNA extracted from enriched viral particles or with total RNA, respectively. The virus shares the greatest degree of overall sequence identity (~78%) with Medicago sativa marafivirus 1 (MsMV1) recently deduced from alfalfa transcriptomic data. The tentative nucleotide sequence of the AVF coat protein shares ~83% identity with the corresponding region of MsMV1. A sequence search of the predicted single large ORF encoding a polyprotein of 235kDa in the Pfam database resulted in identification of five domains, characteristic of the genus Marafivirus, family Tymoviridae. The AVF genome also contains a conserved "marafibox", a 16-nt consensus sequence present in all known marafiviruses. Phylogenetic analysis of the complete nucleotide sequences of AVF and other viruses of the family Tymoviridae grouped AVF in the same cluster with MsMV1. In addition to 5' and 3' terminal extensions, the identity of the virus was confirmed by RT-PCRs with primers derived from VANA-contigs, transmission electron microscopy with virus-infected tissues and transient expression of the viral coat protein gene using a heterologous virus-based vector. Based on the criteria demarcating species in the genus Marafivirus that include overall sequence identity less than 80% and coat protein identity less than 90%, we propose that AVF represents a distinct viral species in the genus Marafivirus, family Tymoviridae.

Dissecting Key Adaptation Traits in the Polyploid Perennial Medicago sativa Using GBS-SNP Mapping

Understanding key adaptation traits is crucial to developing new cultivars with broad adaptations. The main objective of this research is to understand the genetic basis of winter hardiness (WH) and fall dormancy (FD) in alfalfa and the association between the two traits. QTL analysis was conducted in a pseudo-testcross F1 population developed from two cultivars contrasting in FD (3010 with FD = 2 and CW 1010 with FD = 10). The mapping population was evaluated in three replications at two locations (Watkinsville and Blairsville, GA). FD levels showed low to moderate correlations with WH (0.22-0.57). Assessing dormancy in winter is more reliable than in the fall in southern regions with warm winters. The mapping population was genotyped using Genotyping-by-sequencing (GBS). Single dose allele SNPs (SDA) were used for constructing linkage maps. The parental map (CW 1010) consisted of 32 linkage groups spanning 2127.5 cM with 1377 markers and an average marker density of 1.5 cM/SNP. The maternal map (3010) had 32 linkage groups spanning 2788.4 cM with 1837 SDA SNPs with an average marker density of 1.5 cM/SNP. Forty-five significant (P < 0.05) QTLs for FD and 35 QTLs for WH were detected on both male and female linkage maps. More than 75% (22/28) of the dormancy QTL detected from the 3010 parent did not share genomic regions with WH QTLs and more than 70% (12/17) dormancy QTLs detected from CW 1010 parent were localized in different genomic regions than WH QTLs. These results suggest that the two traits have independent inheritance and therefore can be improved separately in breeding programs.

Transcriptome analysis of microRNA156 overexpression alfalfa roots under drought stress

Drought is one of the major abiotic stresses that negatively impact alfalfa growth and productivity. The role of microRNA156 (miR156) in drought has been demonstrated in plants. To date, there are no published studies investigating the role of miR156 in regulating global gene expression in alfalfa under drought. In our study, alfalfa genotypes overexpressing miR156 (miR156OE) exhibited reduced water loss, and enhanced root growth under drought. Our RNA-seq data showed that in response to drought, a total of 415 genes were upregulated and 169 genes were downregulated specifically in miR156OE genotypes. Genotypic comparison revealed that 285 genes were upregulated and 253 genes were downregulated in miR156OE genotypes relative to corresponding WT under drought. Gene Ontology enrichment analysis revealed that the number of differentially expressed genes belonging to biological process, molecular function and cell component functional groups was decreased in miR156OE genotypes under drought. Furthermore, RNA-Seq data showed downregulation of a gene encoding WD40 repeat in a miR156-specific manner. 5' RACE experiments verified cleavage of WD40-2 transcript under drought. Moreover, alfalfa plants overexpressing WD40-2 showed drought sensitive, whereas those with silenced WD40-2 exhibited drought tolerant phenotypes. These findings suggest that miR156 improves drought tolerance in alfalfa by targeting WD40-2.

Screening and identification of key genes regulating fall dormancy in alfalfa leaves

Fall dormancy (FD) determines the adaptation of an alfalfa variety and affects alfalfa production and quality. However, the molecular mechanism underlying FD remains poorly understood. Here, 44 genes regulating FD were identified by comparison of the transcriptomes from leaves of Maverick (fall-dormant alfalfa) and CUF101(non-fall-dormant), during FD and non-FD and were classified them depending on their function. The transcription of IAA-amino acid hydrolase ILR1-like 1, abscisic acid receptor PYL8, and monogalactosyldiacylglycerol synthase-3 in Maverick leaves was regulated by daylength and temperature, and the transcription of the abscisic acid receptor PYL8 was mainly affected by daylength. The changes in the expression of these genes and the abundance of their messenger RNA (mRNA) in Maverick leaves differed from those in CUF101 leaves, as evidenced by the correlation analysis of their mRNA abundance profiles obtained from April to October. The present findings suggested that these genes are involved in regulating FD in alfalfa.

Effect of photoperiod prior to cold acclimation on freezing tolerance and carbohydrate metabolism in alfalfa (Medicago sativa L.)

Cold acclimation proceeds sequentially in response to decreases in photoperiod and temperature. This study aimed at assessing the impact of photoperiod prior to cold acclimation on freezing tolerance and related biochemical and molecular responses in two alfalfa cultivars. The fall dormant cultivar Evolution and semi-dormant cultivar 6010 were grown in growth chambers under different photoperiods (8, 10, 12, 14 or 16h) prior to cold acclimation. Freezing tolerance was evaluated as well as carbohydrate concentrations, levels of transcripts encoding enzymes of carbohydrate metabolism as well as a K-3dehydrin, before and after cold acclimation. The fall dormant cultivar Evolution had a better freezing tolerance than the semi-dormant cultivar 6010. The effect of photoperiod prior to cold acclimation on the level of freezing tolerance differed between the two cultivars: an 8h-photoperiod induced the highest level of freezing tolerance in Evolution and the lowest in 6010. In Evolution, the 8h-induced superior freezing tolerance was associated with higher concentration of raffinose-family oligosaccharides (RFO). The transcript levels of sucrose synthase (SuSy) decreased whereas those of sucrose phosphatase synthase (SPS) and galactinol synthase (GaS) increased in response to cold acclimation in both cultivars. Our results indicate that RFO metabolism could be involved in short photoperiod-induced freezing tolerance in dormant alfalfa cultivars.

Transcriptomic and physiological analyses of Medicago sativa L. roots in response to lead stress

Lead (Pb) is one of the nonessential and toxic metals that threaten the environment and human health. Medicago sativa L. is a legume with high salt tolerance and high biomass production. It is not only a globally important forage crop but is also an ideal plant for phytoremediation. However, the biological and molecular mechanisms that respond to heavy metals are still not well defined in M. sativa. In this study, de novo and strand-specific RNA-sequencing was performed to identify genes involved in the Pb stress response in M. sativa roots. A total of 415,350 unigenes were obtained from the assembled cDNA libraries, among which 5,416 were identified as significantly differentially expressed genes (DEGs) (false discovery rate < 0.005) between cDNA libraries from control and Pb-treated plants. Gene ontology and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses of the DEGs showed they mainly clustered with terms associated with binding, transport, membranes, and the pathways related to signal and energy metabolism. Moreover, a number of candidate genes included antioxidant enzymes, metal transporters, and transcription factors involved in heavy metal response were upregulated under Pb stress. Quantitative real-time PCR(qRT-PCR) validation of the expression patterns of 10 randomly selected candidate DEGs were consistent with the transcriptome analysis results. Thus, this study offers new information towards the investigation of biological changes and molecular mechanisms related to Pb stress response in plants.

A Combined Comparative Transcriptomic, Metabolomic, and Anatomical Analyses of Two Key Domestication Traits: Pod Dehiscence and Seed Dormancy in Pea (Pisum sp.)

The origin of the agriculture was one of the turning points in human history, and a central part of this was the evolution of new plant forms, domesticated crops. Seed dispersal and germination are two key traits which have been selected to facilitate cultivation and harvesting of crops. The objective of this study was to analyze anatomical structure of seed coat and pod, identify metabolic compounds associated with water-impermeable seed coat and differentially expressed genes involved in pea seed dormancy and pod dehiscence. Comparative anatomical, metabolomics, and transcriptomic analyses were carried out on wild dormant, dehiscent Pisum elatius (JI64, VIR320) and cultivated, indehiscent Pisum sativum non-dormant (JI92, Cameor) and recombinant inbred lines (RILs). Considerable differences were found in texture of testa surface, length of macrosclereids, and seed coat thickness. Histochemical and biochemical analyses indicated genotype related variation in composition and heterogeneity of seed coat cell walls within macrosclereids. Liquid chromatography-electrospray ionization/mass spectrometry and Laser desorption/ionization-mass spectrometry of separated seed coats revealed significantly higher contents of proanthocyanidins (dimer and trimer of gallocatechin), quercetin, and myricetin rhamnosides and hydroxylated fatty acids in dormant compared to non-dormant genotypes. Bulk Segregant Analysis coupled to high throughput RNA sequencing resulted in identification of 770 and 148 differentially expressed genes between dormant and non-dormant seeds or dehiscent and indehiscent pods, respectively. The expression of 14 selected dormancy-related genes was studied by qRT-PCR. Of these, expression pattern of four genes: porin (MACE-S082), peroxisomal membrane PEX14-like protein (MACE-S108), 4-coumarate CoA ligase (MACE-S131), and UDP-glucosyl transferase (MACE-S139) was in agreement in all four genotypes with Massive analysis of cDNA Ends (MACE) data. In case of pod dehiscence, the analysis of two candidate genes (SHATTERING and SHATTERPROOF) and three out of 20 MACE identified genes (MACE-P004, MACE-P013, MACE-P015) showed down-expression in dorsal and ventral pod suture of indehiscent genotypes. Moreover, MACE-P015, the homolog of peptidoglycan-binding domain or proline-rich extensin-like protein mapped correctly to predicted Dpo1 locus on PsLGIII. This integrated analysis of the seed coat in wild and cultivated pea provides new insight as well as raises new questions associated with domestication and seed dormancy and pod dehiscence.

Comparative transcriptome investigation of global gene expression changes caused by miR156 overexpression in Medicago sativa

Background

Medicago sativa (alfalfa) is a low-input forage and potential bioenergy crop, and improving its yield and quality has always been a focus of the alfalfa breeding industry. Transgenic alfalfa plants overexpressing a precursor of alfalfa microRNA156 (MsmiR156) were recently generated by our group. These plants (miR156OE) showed enhanced biomass yield, reduced internodal length, increased shoot branching and trichome density, and a delay in flowering time. Transcripts of three SQUAMOSA-PROMOTER BINDING PROTEIN-LIKE (SPL) genes (MsSPL6, MsSPL12, and MsSPL13) were found to be targeted for cleavage by MsmiR156 in alfalfa.

Results

To further illustrate the molecular mechanisms underlying the effects of miR156 in alfalfa, two miR156OE genotypes (A11a and A17) were subjected to Next Generation RNA Sequencing with Illumina HiSeq. More than 1.11 billion clean reads were obtained from our available sequenced samples. A total of 160,472 transcripts were generated using Trinity de novo assembly and 4,985 significantly differentially expressed genes were detected in miR156OE plants A11a and A17 using the Medicago truncatula genome as reference. A total of 17 genes (including upregulated, downregulated, and unchanged) were selected for quantitative real-time PCR (qRT-PCR) validation, which showed that gene expression levels were largely consistent between qRT-PCR and RNA-Seq data. In addition to the established SPL genes MsSPL6, MsSPL12 and MsSPL13, four new SPLs; MsSPL2, MsSPL3, MsSPL4 and MsSPL9 were also down-regulated significantly in both miR156OE plants. These seven SPL genes belong to genes phylogeny clades VI, IV, VIII, V and VII, which have been reported to be targeted by miR156 in Arabidopsis thaliana. The gene ontology terms characterized electron transporter, starch synthase activity, sucrose transport, sucrose-phosphate synthase activity, chitin binding, sexual reproduction, flavonoid biosynthesis and lignin catabolism correlate well to the phenotypes of miR156OE alfalfa plants.

Conclusions

This is the first report of changes in global gene expression in response to miR156 overexpression in alfalfa. The discovered miR156-targeted SPL genes belonging to different clades indicate miR156 plays fundamental and multifunctional roles in regulating alfalfa plant development.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3014-6) contains supplementary material, which is available to authorized users.

De novo transcriptome assembly and analysis of differentially expressed genes of two barley genotypes reveal root-zone-specific responses to salt exposure

Plant roots are the first organs sensing and responding to salinity stress, manifested differentially between different root types, and also at the individual tissue and cellular level. High genetic diversity and the current lack of an assembled map-based sequence of the barley genome severely limit barley research potential. We used over 580 and 600 million paired-end reads, respectively, to create two de novo assemblies of a barley landrace (Sahara) and a malting cultivar (Clipper) with known contrasting responses to salinity. Generalized linear models were used to statistically access spatial, treatment-related, and genotype-specific responses. This revealed a spatial gene expression gradient along the barley root, with more differentially expressed transcripts detected between different root zones than between treatments. The root transcriptome also showed a gradual transition from transcripts related to sugar-mediated signaling at the root meristematic zone to those involved in cell wall metabolism in the elongation zone, and defense response-related pathways toward the maturation zone, with significant differences between the two genotypes. The availability of these additional transcriptome reference sets will serve as a valuable resource to the cereal research community, and may identify valuable traits to assist in breeding programmes.

Deep-sequencing transcriptome analysis of field-grown Medicago sativa L. crown buds acclimated to freezing stress

Medicago sativa L. (alfalfa) 'Zhaodong' is an important forage legume that can safely survive in northern China where winter temperatures reach as low as -30 °C. Survival of alfalfa following freezing stress depends on the amount and revival ability of crown buds. In order to investigate the molecular mechanisms of frost tolerance in alfalfa, we used transcriptome sequencing technology and bioinformatics strategies to analyze crown buds of field-grown alfalfa during winter. We statistically identified a total of 5605 differentially expressed genes (DEGs) involved in freezing stress including 1900 upregulated and 3705 downregulated DEGs. We validated 36 candidate DEGs using qPCR to confirm the accuracy of the RNA-seq data. Unlike other recent studies, this study employed alfalfa plants grown in the natural environment. Our results indicate that not only the CBF orthologs but also membrane proteins, hormone signal transduction pathways, and ubiquitin-mediated proteolysis pathways indicate the presence of a special freezing adaptation mechanism in alfalfa. The antioxidant defense system may rapidly confer freezing tolerance to alfalfa. Importantly, biosynthesis of secondary metabolites and phenylalanine metabolism, which is of potential importance in coordinating freezing tolerance with growth and development, were downregulated in subzero temperatures. The adaptive mechanism for frost tolerance is a complex multigenic process that is not well understood. This systematic analysis provided an in-depth view of stress tolerance mechanisms in alfalfa.

Transcriptome analysis and molecular marker discovery in Solanum incanum and S. aethiopicum, two close relatives of the common eggplant (Solanum melongena) with interest for breeding

Background

Solanum incanum is a close wild relative of S. melongena with high contents of bioactive phenolics and drought tolerance. S. aethiopicum is a cultivated African eggplant cross-compatible with S. melongena. Despite their great interest in S. melongena breeding programs, the genomic resources for these species are scarce.

Results

RNA-Seq was performed with NGS from pooled RNA of young leaf, floral bud and young fruit tissues, generating more than one hundred millions raw reads per species. The transcriptomes were assembled in 83,905 unigenes for S. incanum and in 87,084 unigenes for S. aethiopicum with an average length of 696 and 722 bp, respectively. The unigenes were structurally and functionally annotated based on comparison with public databases by using bioinformatic tools. The single nucleotide variant calling analysis (SNPs and INDELs) was performed by mapping our S. incanum and S. aethiopicum reads, as well as reads from S. melongena and S. torvum available on NCBI database (National Center for Biotechnology Information), against the eggplant genome. Both intraspecific and interspecific polymorphisms were identified and subsets of molecular markers were created for all species combinations. 36 SNVs were selected for validation in the S. incanum and S. aethiopicum accessions and 96 % were correctly amplified confirming the polymorphisms. In addition, 976 and 1,278 SSRs were identified in S. incanum and S. aethiopicum transcriptomes respectively, and a set of them were validated.

Conclusions

This work provides a broad insight into gene sequences and allelic variation in S. incanum and S. aethiopicum. This work is a first step toward better understanding of target genes involved in metabolic pathways relevant for eggplant breeding. The molecular markers detected in this study could be used across all the eggplant genepool, which is of interest for breeding programs as well as to perform marker-trait association and QTL analysis studies.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2631-4) contains supplementary material, which is available to authorized users.