Salinity is one of the major abiotic factors affecting alfalfa productivity. Identifying genes that control this complex trait will provide critical insights for alfalfa breeding programs. To date, no studies have been published on a deep sequencing-based profiling of the alfalfa transcriptome in response to salinity stress. Observations gathered through research on reference genomes may not always be applicable to alfalfa. In this work, Illumina RNA-sequencing was performed in two alfalfa genotypes contrasting in salt tolerance, in order to estimate a broad spectrum of genes affected by salt stress. A total of 367,619,586 short reads were generated from cDNA libraries originated from roots of both lines. More than 60,000 tentative consensus sequences (TCs) were obtained and, among them, 74.5% had a significant similarity to proteins in the NCBI database. Mining of simple sequence repeats (SSRs) from all TCs revealed 6,496 SSRs belonging to 3,183 annotated unigenes. Bioinformatics analysis showed that the expression of 1,165 genes, including 86 transcription factors (TFs), was significantly altered under salt stress. About 40% of differentially expressed genes were assigned to known gene ontology (GO) categories using Arabidopsis GO. A random check of differentially expressed genes by quantitative real-time PCR confirmed the bioinformatic analysis of the RNA-seq data. A number of salt-responsive genes in both tested genotypes were identified and assigned to functional classes, and gene candidates with roles in the adaptation to salinity were proposed. Alfalfa-specific data on salt-responsive genes obtained in this work will be useful in understanding the molecular mechanisms of salinity tolerance in alfalfa.

Nematodes are one of the major limiting factors in alfalfa production. Root-knot nematodes (RKN, Meloidogyne spp.) are widely distributed and economically important sedentary endoparasites of agricultural crops and they may inflict significant damage to alfalfa fields. As of today, no studies have been published on global gene expression profiling in alfalfa infected with RKN or any other plant parasitic nematode. Very little information is available about molecular mechanisms that contribute to pathogenesis and defense responses in alfalfa against these pests and specifically against RKN. In this work, we performed root transcriptome analysis of resistant (cv. Moapa 69) and susceptible (cv. Lahontan) alfalfa cultivars infected with RKN Meloidogyne incognita, widespread root-knot nematode species and a major pest worldwide. A total of 1,701,622,580 pair-end reads were generated on an Illumina Hi-Seq 2000 platform from the roots of both cultivars and assembled into 45,595 and 47,590 transcripts in cvs Moapa 69 and Lahontan, respectively. Bioinformatic analysis revealed a number of common and unique genes that were differentially expressed in susceptible and resistant lines as a result of nematode infection. Although the susceptible cultivar showed a more pronounced defense response to the infection, feeding sites were successfully established in its roots. Characteristically, basal gene expression levels under normal conditions differed between the two cultivars as well, which may confer advantage to one of the genotypes toward resistance to nematodes. Differentially expressed genes were subsequently assigned to known Gene Ontology categories to predict their functional roles and associated biological processes. Real-time PCR validated expression changes in genes arbitrarily selected for experimental confirmation. Candidate genes that contribute to protection against M. incognita in alfalfa were proposed and alfalfa-nematode interactions with respect to resistance are discussed.

At the moment, there are a number of publications describing gene expression profiling in virus-infected plants. Most of the data are limited to specific host-pathogen interactions involving a given virus and a model host plant – usually Arabidopsis thaliana. Even though several summarizing attempts have been made, a general picture of gene expression changes in susceptible virus-host interactions is lacking.

Pseudomonas syringae infects diverse crop plants and comprises at least 50 different pathovar strains with different host ranges. More information on the physiological and molecular effects of the host inhibitory environment on the pathogen is needed to develop resistant cultivars. Recently, we reported an in vitro model system that mimics the redox pulse associated with the oxidative burst in plant cells inoculated with Pseudomonas syringae pv. syringae. Using this system, we demonstrated that oxidation of acetosyringone, a major extracellular phenolic compound induced in some plants in response to bacteria, rendered Pseudomonas syringae. pv. syringae to a “viable but nonculturable” (VBNC) state. Here we performed a large scale transcriptome profiling of P.s.pv. syringae in the VBNC state induced by acetosyringone treatment and identified bacterial genes and pathways presumably associated with this condition. The findings offer insight into what events occur when bacterial pathogens are first encountered and host defense responses are triggered. The acquired knowledge will improve our understanding of the molecular mechanisms of stress tolerance. We believe that this is the first work on global gene expression profiling of VBNC cells in plant pathogenic bacteria.

Bacterial stem blight caused by Pseudomonas syringae pv. syringae is a common disease of alfalfa (Medicago sativa L). Little is known about host-pathogen interactions and host defense mechanisms. Here, individual resistant and susceptible plants were selected from cultivars Maverick and ZG9830 and used for transcript profiling at 24 and 72 hours after inoculation (hai) with the isolate PssALF3. Bioinformatic analysis revealed a number of differentially expressed genes (DEGs) in resistant and susceptible genotypes. Although resistant plants from each cultivar produced a hypersensitive response, transcriptome analyses indicated that they respond differently at the molecular level. The number of DEGs was higher in resistant plants of ZG9830 at 24 hai than in Maverick, suggesting that ZG9830 plants had a more rapid effector triggered immune response. Unique up-regulated genes in resistant ZG9830 plants included genes encoding putative nematode resistance HSPRO2-like proteins, orthologs for the rice Xa21 and soybean Rpg1-b resistance genes, and TIR-containing R genes lacking both NBS and LRR domains. The suite of R genes up-regulated in resistant Maverick plants had an over-representation of R genes in the CC-NBS-LRR family including two genes for atypical CC_R domains and a putative ortholog of the Arabidopsis RPM1 gene. Resistance in both cultivars appears to be mediated primarily by WRKY family transcription factors and expression of genes involved in protein phosphorylation, regulation of transcription, defense response including synthesis of isoflavonoids, and oxidation-reduction processes. These results will further the identification of mechanisms involved in resistance to facilitate selection of parent populations and development of commercial varieties.