Transcript profiling of two alfalfa genotypes with contrasting cell wall composition in stems using a cross-species platform: optimizing analysis by masking biased probes

Comparative Transcriptome Analysis Reveals New Insight of Alfalfa (Medicago sativa L.) Cultivars in Response to Abrupt Freezing Stress

Freezing stress is a major limiting environmental factor that affects the productivity and distribution of alfalfa (Medicago sativa L.). There is growing evidence that enhancing freezing tolerance through resistance-related genes is one of the most efficient methods for solving this problem, whereas little is known about the complex regulatory mechanism of freezing stress. Herein, we performed transcriptome profiling of the leaves from two genotypes of alfalfa, freezing tolerance "Gannong NO.3" and freezing-sensitive "WL326GZ" exposure to -10°C to investigate which resistance-related genes could improve the freezing tolerance. Our results showed that a total of 121,366 genes were identified, and there were 7,245 differentially expressed genes (DEGs) between the control and treated leaves. In particular, the DEGs in "Gannong NO.3" were mainly enriched in the metabolic pathways and biosynthesis of secondary metabolites, and most of the DEGs in "WL326GZ" were enriched in the metabolic pathways, the biosynthesis of secondary metabolites, and plant-pathogen interactions. Moreover, the weighted gene co-expression network analysis (WGCNA) showed that ATP-binding cassette (ABC) C subfamily genes were strongly impacted by freezing stress, indicating that ABCC8 and ABCC3 are critical to develop the freezing tolerance. Moreover, our data revealed that numerous Ca2+ signal transduction and CBF/DREB1 pathway-related genes were severely impacted by the freezing resistance, which is believed to alleviate the damage caused by freezing stress. Altogether, these findings contribute the comprehensive information to understand the molecular mechanism of alfalfa adaptation to freezing stress and further provide functional candidate genes that can adapt to abiotic stress.

Changes in phenotype and gene expression under lead stress revealed key genetic responses to lead tolerance in Medicago sativa L

Lead (Pb) is a serious heavy metal soil pollutant. It can be absorbed and accumulated by plant roots and impact plant growth. Medicago sativa L. (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. Little is known about the mechanism by which alfalfa responds to Pb stress at the molecular level. In this study, three alfalfa genotypes (a lead-resistant type (LR), a lead-sensitive type (LS) and an intermediate type (IN)) with contrasting abilities to resist lead were exposed to different durations of Pb treatment. Next-generation sequencing (NGS)-based RNA-seq technology was employed to characterize the root transcriptomes of three genotypes of alfalfa and identify differentially expressed genes (DEGs) during Pb stress. Genotypes LR and LS displayed different mechanisms of tolerance. In LR, the accumulation of more resistant substances was induced by the upregulation of sucrose synthase, glucan endo-1,3-beta-glucosidase, beta-amylase 3, probable trehalose-phosphate phosphatase J, 6-phosphofructo-2-kinase delta-1-pyrroline-5-carboxylate synthase (P5CS) and δ-ornithine aminotransferase (δ-OAT). In addition, flavin monooxygenase (YUCCA), 4-coumarate:CoA ligase-like protein (4CL), cinnamoyl-CoA reductase-like protein (CCR), ferulate 5-hydroxylase (F5H) and caffeic acid O-methyltransferase (COMT) were upregulated, leading to root development in a short time under Pb stress. Further study of the expression levels of metal transport-related genes, such as NRAMP (metal transporter), MATE (multidrug and toxin extrusion), HIPPs (heavy metal-associated isoprenylated plant proteins), MTP (metal tolerance protein), and ABC transporter, suggested that these genes were differentially expressed after lead treatment in the three alfalfa genotypes. Our research provides useful information for further studies on the molecular mechanism of Pb resistance in Medicago sativa L.

iTRAQ-based comparative proteomic analysis of differences in the protein profiles of stems and leaves from two alfalfa genotypes

BACKGROUND

To explore the molecular regulatory mechanisms of early stem and leaf development, proteomic analysis was performed on leaves and stems of F genotype alfalfa, with thin stems and small leaves, and M genotype alfalfa, with thick stems and large leaves.

RESULTS

Based on fold-change thresholds of > 1.20 or < 0.83 (p < 0.05), a large number of proteins were identified as being differentially enriched between the M and F genotypes: 249 downregulated and 139 upregulated in stems and 164 downregulated and 134 upregulated in leaves. The differentially enriched proteins in stems were mainly involved in amino acid biosynthesis, phenylpropanoid biosynthesis, carbon fixation, and phenylalanine metabolism. The differentially enriched proteins in leaves were mainly involved in porphyrin and chlorophyll metabolism, phenylpropanoid biosynthesis, starch and sucrose metabolism, and carbon fixation in photosynthetic organisms. Six differentially enriched proteins were mapped onto the porphyrin and chlorophyll metabolism pathway in leaves of the M genotype, including five upregulated proteins involved in chlorophyll biosynthesis and one downregulated protein involved in chlorophyll degradation. Eleven differentially enriched proteins were mapped onto the phenylpropanoid pathway in stems of the M genotype, including two upregulated proteins and nine downregulated proteins.

CONCLUSION

Enhanced chlorophyll synthesis and decreased lignin synthesis provided a reasonable explanation for the larger leaves and lower levels of stem lignification in M genotype alfalfa. This proteomic study aimed to classify the functions of differentially enriched proteins and to provide information on the molecular regulatory networks involved in stem and leaf development.

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.

Secondary Wall Regulating NACs Differentially Bind at the Promoter at a CELLULOSE SYNTHASE A4 Cis-eQTL

Arabidopsis thaliana CELLULOSE SYNTHASE A4/7/8 (CESA4/7/8) are three non-redundant subunits of the secondary cell wall cellulose synthase complex. Transcript abundance of these genes can vary among genotypes and expression quantitative trait loci (eQTL) were identified in a recombinant population of the accessions Bay-0 and Shahdara. Genetic mapping and analysis of the transcript levels of CESAs between two distinct near isogenic lines (NILs) confirmed a change in CESA4 expression that segregates within that interval. We sequenced the promoters and identified 16 polymorphisms differentiating CESA4^Sha and CESA4^Bay . In order to determine which of these SNPs could be responsible for this eQTL, we screened for transcription factor protein affinity with promoter fragments of CESA4^Bay, CESA4^Sha , and the reference genome CESA4^Col . The wall thickening activator proteins NAC SECONDARY WALL THICKENING PROMOTING FACTOR2 (NST2) and NST3 exhibited a decrease in binding with the CESA4^Sha promoter with a tracheary element-regulating cis-element (TERE) polymorphism. While NILs harboring the TERE polymorphisms exhibited significantly different CESA4 expression, cellulose crystallinity and cell wall thickness were indistinguishable. These results suggest that the TERE polymorphism resulted in differential transcription factor binding and CESA4 expression; yet A. thaliana is able to tolerate this transcriptional variability without compromising the structural elements of the plant, providing insight into the elasticity of gene regulation as it pertains to cell wall biosynthesis and regulation. We also explored available DNA affinity purification sequencing data to resolve a core binding site, C(G/T)TNNNNNNNA(A/C)G, for secondary wall NACs referred to as the VNS element.

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.

Ups and downs in alfalfa: Proteomic and metabolic changes occurring in the growing stem

The expanding interest for using lignocellulosic biomass in industry spurred the study of the mechanisms underlying plant cell-wall synthesis. Efforts using genetic approaches allowed the disentanglement of major steps governing stem fibre synthesis. Nonetheless, little is known about the relations between the stem maturation and the evolution of its proteome. During Medicago sativa L. maturation, the different internodes grow asynchronously allowing the discrimination of various developmental stages on a same stem. In this study, the proteome of three selected regions of the stem of alfalfa (apical, intermediate and basal) was analyzed and combined with a compositional analysis of the different stem parts. Interestingly, the apical and the median regions share many similarities: high abundance of chloroplast- and mitochondrial-related proteins together with the accumulation of proteins acting in the early steps of fibre production. In the mature basal region, forisomes and stress-related proteins accumulate. The RT-qPCR assessment of the expression of genes coding for members of the cellulose synthase family likewise indicates that fibres and the machinery responsible for the deposition of secondary cell walls are predominantly formed in the apical section. Altogether, this study reflects the metabolic change from the fibre production in the upper stem regions to the acquisition of defence-related functions in the fibrous basal part.

Analysis of Cell Wall-Related Genes in Organs of Medicago sativa L. under Different Abiotic Stresses

Abiotic constraints are a source of concern in agriculture, because they can have a strong impact on plant growth and development, thereby affecting crop yield. The response of plants to abiotic constraints varies depending on the type of stress, on the species and on the organs. Although many studies have addressed different aspects of the plant response to abiotic stresses, only a handful has focused on the role of the cell wall. A targeted approach has been used here to study the expression of cell wall-related genes in different organs of alfalfa plants subjected for four days to three different abiotic stress treatments, namely salt, cold and heat stress. Genes involved in different steps of cell wall formation (cellulose biosynthesis, monolignol biosynthesis and polymerization) have been analyzed in different organs of Medicago sativa L. Prior to this analysis, an in silico classification of dirigent/dirigent-like proteins and class III peroxidases has been performed in Medicago truncatula and M. sativa. The final goal of this study is to infer and compare the expression patterns of cell wall-related genes in response to different abiotic stressors in the organs of an important legume crop.

Wood biosynthesis and typologies: a molecular rhapsody

Wood represents one of the most important renewable commodities for humanity and plays a crucial role in terrestrial ecosystem carbon-cycling. Wood formation is the result of a multitude of events that require the concerted action of endogenous and exogenous factors under the influence of photoperiod, for instance genes and plant growth regulators. Beyond providing mechanical support and being responsible for the increase in stem radial diameter, woody tissues constitute the vascular system of trees and are capable of reacting to environmental stimuli, and as such are therefore quite plastic and responsive. Despite the ecological and economic importance of wood, not all aspects of its formation have been unveiled. Many gaps in our knowledge are still present, which hinder the maximal exploitation of this precious bioresource. This review aims at surveying the current knowledge of wood formation and the available molecular data addressing the relationship between wood production and environmental factors, which have crucial influences on the rhythmic regulation of cambial activity and exert profound effects on tree stem growth, wood yield and properties. We will here go beyond wood sensu stricto, i.e., secondary xylem, and extend our survey to other tissues, namely vascular cambium, phloem and fibres. The purpose is to provide the reader with an overview of the complexity of the topic and to highlight the importance of progressing in the future towards an integrated knowledge on the subject.

Alfalfa Cellulose synthase gene expression under abiotic stress: a Hitchhiker's guide to RT-qPCR normalization

Abiotic stress represents a serious threat affecting both plant fitness and productivity. One of the promptest responses that plants trigger following abiotic stress is the differential expression of key genes, which enable to face the adverse conditions. It is accepted and shown that the cell wall senses and broadcasts the stress signal to the interior of the cell, by triggering a cascade of reactions leading to resistance. Therefore the study of wall-related genes is particularly relevant to understand the metabolic remodeling triggered by plants in response to exogenous stresses. Despite the agricultural and economical relevance of alfalfa (Medicago sativa L.), no study, to our knowledge, has addressed specifically the wall-related gene expression changes in response to exogenous stresses in this important crop, by monitoring the dynamics of wall biosynthetic gene expression. We here identify and analyze the expression profiles of nine cellulose synthases, together with other wall-related genes, in stems of alfalfa plants subjected to different abiotic stresses (cold, heat, salt stress) at various time points (e.g. 0, 24, 72 and 96 h). We identify 2 main responses for specific groups of genes, i.e. a salt/heat-induced and a cold/heat-repressed group of genes. Prior to this analysis we identified appropriate reference genes for expression analyses in alfalfa, by evaluating the stability of 10 candidates across different tissues (namely leaves, stems, roots), under the different abiotic stresses and time points chosen. The results obtained confirm an active role played by the cell wall in response to exogenous stimuli and constitute a step forward in delineating the complex pathways regulating the response of plants to abiotic stresses.

Comparison of stem morphology and anatomy of two alfalfa clonal lines exhibiting divergent cell wall composition

BACKGROUND

In previous research, two alfalfa clonal lines (252 and 1283) were identified that exhibited environmentally stable differences in stem cell walls. Compared with stems of 1283, stems of 252 have a higher cell wall concentration and greater amounts of lignin and cellulose but reduced levels of pectic sugar residues. These results suggest greater deposition of secondary xylem and a reduction in pith in stems of 252 compared with 1283.

RESULTS

The stem morphology and anatomy of first-cut and second-cut harvests of field-grown 1283 and 252 were examined. For both harvests, stems of 1283 were thicker and had a higher leaf/stem ratio compared with stems of 252. Stem cross-sections of both genotypes were stained for lignin, and the proportions of stem area that were pith and secondary xylem were measured using ImageJ. Stems of 252 exhibited greater deposition of secondary xylem and a reduction in pith proportion compared with stems of 1283 for the first-cut harvest, but this difference was not statistically significant for the second-cut harvest.

CONCLUSION

The results indicate that the proportions of secondary xylem and pith are not environmentally stable in these two genotypes and hence cannot be the sole basis for the differences in cell wall concentration/composition.

Modifying crops to increase cell wall digestibility

Improving digestibility of roughage cell walls will improve ruminant animal performance and reduce loss of nutrients to the environment. The main digestibility impediment for dicotyledonous plants is highly lignified secondary cell walls, notably in stem secondary xylem, which become almost non-digestible. Digestibility of grasses is slowed severely by lignification of most tissues, but these cell walls remain largely digestible. Cell wall lignification creates an access barrier to potentially digestible wall material by rumen bacteria if cells have not been physically ruptured. Traditional breeding has focused on increasing total dry matter digestibility rather than cell wall digestibility, which has resulted in minimal reductions in cell wall lignification. Brown midrib mutants in some annual grasses exhibit small reductions in lignin concentration and improved cell wall digestibility. Similarly, transgenic approaches down-regulating genes in monolignol synthesis have produced plants with reduced lignin content and improved cell wall digestibility. While major reductions in lignin concentration have been associated with poor plant fitness, smaller reductions in lignin provided measurable improvements in digestibility without significantly impacting agronomic fitness. Additional targets for genetic modification to enhance digestibility and improve roughages for use as biofuel feedstocks are discussed; including manipulating cell wall polysaccharide composition, novel lignin structures, reduced lignin/polysaccharide cross-linking, smaller lignin polymers, enhanced development of non-lignified tissues, and targeting specific cell types. Greater tissue specificity of transgene expression will be needed to maximize benefits while avoiding negative impacts on plant fitness.cauliflower mosiac virus (CaMV) 35S promoter.

Proteomic analysis of cell walls of two developmental stages of alfalfa stems

Cell walls are important for the growth and development of all plants. They are also valuable resources for feed and fiber, and more recently as a potential feedstock for bioenergy production. Cell wall proteins comprise only a fraction of the cell wall, but play important roles in establishing the walls and in the chemical interactions (e.g., crosslinking) of cell wall components. This crosslinking provides structure, but restricts digestibility of cell wall complex carbohydrates, limiting available energy in animal and bioenergy production systems. Manipulation of cell wall proteins could be a strategy to improve digestibility. An analysis of the cell wall proteome of apical alfalfa stems (less mature, more digestible) and basal alfalfa stems (more mature, less digestible) was conducted using a recently developed low-salt/density gradient method for the isolation of cell walls. Walls were subsequently subjected to a modified extraction utilizing EGTA to remove pectins, followed by a LiCl extraction to isolate more tightly bound proteins. Recovered proteins were identified using shotgun proteomics. We identified 272 proteins in the alfalfa stem cell wall proteome, 153 of which had not previously been identified in cell wall proteomic analyses. Nearly 70% of the identified proteins were predicted to be secreted, as would be expected for most cell wall proteins, an improvement over previously published studies using traditional cell wall isolation methods. A comparison of our and several other cell wall proteomic studies indicates little overlap in identified proteins among them, which may be largely due to differences in the tissues used as well as differences in experimental approach.

Using RNA-Seq for gene identification, polymorphism detection and transcript profiling in two alfalfa genotypes with divergent cell wall composition in stems

BACKGROUND

Alfalfa, [Medicago sativa (L.) sativa], a widely-grown perennial forage has potential for development as a cellulosic ethanol feedstock. However, the genomics of alfalfa, a non-model species, is still in its infancy. The recent advent of RNA-Seq, a massively parallel sequencing method for transcriptome analysis, provides an opportunity to expand the identification of alfalfa genes and polymorphisms, and conduct in-depth transcript profiling.

RESULTS

Cell walls in stems of alfalfa genotype 708 have higher cellulose and lower lignin concentrations compared to cell walls in stems of genotype 773. Using the Illumina GA-II platform, a total of 198,861,304 expression sequence tags (ESTs, 76 bp in length) were generated from cDNA libraries derived from elongating stem (ES) and post-elongation stem (PES) internodes of 708 and 773. In addition, 341,984 ESTs were generated from ES and PES internodes of genotype 773 using the GS FLX Titanium platform. The first alfalfa (Medicago sativa) gene index (MSGI 1.0) was assembled using the Sanger ESTs available from GenBank, the GS FLX Titanium EST sequences, and the de novo assembled Illumina sequences. MSGI 1.0 contains 124,025 unique sequences including 22,729 tentative consensus sequences (TCs), 22,315 singletons and 78,981 pseudo-singletons. We identified a total of 1,294 simple sequence repeats (SSR) among the sequences in MSGI 1.0. In addition, a total of 10,826 single nucleotide polymorphisms (SNPs) were predicted between the two genotypes. Out of 55 SNPs randomly selected for experimental validation, 47 (85%) were polymorphic between the two genotypes. We also identified numerous allelic variations within each genotype. Digital gene expression analysis identified numerous candidate genes that may play a role in stem development as well as candidate genes that may contribute to the differences in cell wall composition in stems of the two genotypes.

CONCLUSIONS

Our results demonstrate that RNA-Seq can be successfully used for gene identification, polymorphism detection and transcript profiling in alfalfa, a non-model, allogamous, autotetraploid species. The alfalfa gene index assembled in this study, and the SNPs, SSRs and candidate genes identified can be used to improve alfalfa as a forage crop and cellulosic feedstock.

Using RNA-Seq for gene identification, polymorphism detection and transcript profiling in two alfalfa genotypes with divergent cell wall composition in stems

Background

Alfalfa, [Medicago sativa (L.) sativa], a widely-grown perennial forage has potential for development as a cellulosic ethanol feedstock. However, the genomics of alfalfa, a non-model species, is still in its infancy. The recent advent of RNA-Seq, a massively parallel sequencing method for transcriptome analysis, provides an opportunity to expand the identification of alfalfa genes and polymorphisms, and conduct in-depth transcript profiling.

Results

Cell walls in stems of alfalfa genotype 708 have higher cellulose and lower lignin concentrations compared to cell walls in stems of genotype 773. Using the Illumina GA-II platform, a total of 198,861,304 expression sequence tags (ESTs, 76 bp in length) were generated from cDNA libraries derived from elongating stem (ES) and post-elongation stem (PES) internodes of 708 and 773. In addition, 341,984 ESTs were generated from ES and PES internodes of genotype 773 using the GS FLX Titanium platform. The first alfalfa (Medicago sativa) gene index (MSGI 1.0) was assembled using the Sanger ESTs available from GenBank, the GS FLX Titanium EST sequences, and the de novo assembled Illumina sequences. MSGI 1.0 contains 124,025 unique sequences including 22,729 tentative consensus sequences (TCs), 22,315 singletons and 78,981 pseudo-singletons. We identified a total of 1,294 simple sequence repeats (SSR) among the sequences in MSGI 1.0. In addition, a total of 10,826 single nucleotide polymorphisms (SNPs) were predicted between the two genotypes. Out of 55 SNPs randomly selected for experimental validation, 47 (85%) were polymorphic between the two genotypes. We also identified numerous allelic variations within each genotype. Digital gene expression analysis identified numerous candidate genes that may play a role in stem development as well as candidate genes that may contribute to the differences in cell wall composition in stems of the two genotypes.

Conclusions

Our results demonstrate that RNA-Seq can be successfully used for gene identification, polymorphism detection and transcript profiling in alfalfa, a non-model, allogamous, autotetraploid species. The alfalfa gene index assembled in this study, and the SNPs, SSRs and candidate genes identified can be used to improve alfalfa as a forage crop and cellulosic feedstock.

Transcript profiling of common bean (Phaseolus vulgaris L.) using the GeneChip Soybean Genome Array: optimizing analysis by masking biased probes

BACKGROUND

Common bean (Phaseolus vulgaris L.) and soybean (Glycine max) both belong to the Phaseoleae tribe and share significant coding sequence homology. This suggests that the GeneChip(R) Soybean Genome Array (soybean GeneChip) may be used for gene expression studies using common bean.

RESULTS

To evaluate the utility of the soybean GeneChip for transcript profiling of common bean, we hybridized cRNAs purified from nodule, leaf, and root of common bean and soybean in triplicate to the soybean GeneChip. Initial data analysis showed a decreased sensitivity and accuracy of measuring differential gene expression in common bean cross-species hybridization (CSH) GeneChip data compared to that of soybean. We employed a method that masked putative probes targeting inter-species variable (ISV) regions between common bean and soybean. A masking signal intensity threshold was selected that optimized both sensitivity and accuracy of measuring differential gene expression. After masking for ISV regions, the number of differentially-expressed genes identified in common bean was increased by 2.8-fold reflecting increased sensitivity. Quantitative RT-PCR (qRT-PCR) analysis of 20 randomly selected genes and purine-ureide pathway genes demonstrated an increased accuracy of measuring differential gene expression after masking for ISV regions. We also evaluated masked probe frequency per probe set to gain insight into the sequence divergence pattern between common bean and soybean. The sequence divergence pattern analysis suggested that the genes for basic cellular functions and metabolism were highly conserved between soybean and common bean. Additionally, our results show that some classes of genes, particularly those associated with environmental adaptation, are highly divergent.

CONCLUSIONS

The soybean GeneChip is a suitable cross-species platform for transcript profiling in common bean when used in combination with the masking protocol described. In addition to transcript profiling, CSH of the GeneChip in combination with masking probes in the ISV regions can be used for comparative ecological and/or evolutionary genomics studies.