Michelle A. Graham

RNA-Seq Atlas of Glycine max : A guide to the soybean transcriptome

BackgroundNext generation sequencing is transforming our understanding of transcriptomes. It can determine the expression level of transcripts with a dynamic range of over six orders of magnitude from multiple tissues, developmental stages or conditions. Patterns of gene expression provide insight into functions of genes with unknown annotation.ResultsThe RNA Seq-Atlas presented here provides a record of high-resolution gene expression in a set of fourteen diverse tissues. Hierarchical clustering of transcriptional profiles for these tissues suggests three clades with similar profiles: aerial, underground and seed tissues. We also investigate the relationship between gene structure and gene expression and find a correlation between gene length and expression. Additionally, we find dramatic tissue-specific gene expression of both the most highly-expressed genes and the genes specific to legumes in seed development and nodule tissues. Analysis of the gene expression profiles of over 2,000 genes with preferential gene expression in seed suggests there are more than 177 genes with functional roles that are involved in the economically important seed filling process. Finally, the Seq-atlas also provides a means of evaluating existing gene model annotations for the Glycine max genome.ConclusionsThis RNA-Seq atlas extends the analyses of previous gene expression atlases performed using Affymetrix GeneChip technology and provides an example of new methods to accommodate the increase in transcriptome data obtained from next generation sequencing. Data contained within this RNA-Seq atlas of Glycine max can be explored at http://www.soybase.org/soyseq.

Phosphorus Stress in Common Bean: Root Transcript and Metabolic Responses

Phosphorus (P) is an essential element for plant growth. Crop production of common bean (Phaseolus vulgaris), the most important legume for human consumption, is often limited by low P in the soil. Functional genomics were used to investigate global gene expression and metabolic responses of bean plants grown under P-deficient and P-sufficient conditions. P-deficient plants showed enhanced root to shoot ratio accompanied by reduced leaf area and net photosynthesis rates. Transcript profiling was performed through hybridization of nylon filter arrays spotted with cDNAs of 2,212 unigenes from a P deficiency root cDNA library. A total of 126 genes, representing different functional categories, showed significant differential expression in response to P: 62% of these were induced in P-deficient roots. A set of 372 bean transcription factor (TF) genes, coding for proteins with Inter-Pro domains characteristic or diagnostic for TF, were identified from The Institute of Genomic Research/Dana Farber Cancer Institute Common Bean Gene Index. Using real-time reverse transcription-polymerase chain reaction analysis, 17 TF genes were differentially expressed in P-deficient roots; four TF genes, including MYB TFs, were induced. Nonbiased metabolite profiling was used to assess the degree to which changes in gene expression in P-deficient roots affect overall metabolism. Stress-related metabolites such as polyols accumulated in P-deficient roots as well as sugars, which are known to be essential for P stress gene induction. Candidate genes have been identified that may contribute to root adaptation to P deficiency and be useful for improvement of common bean.

Computational Identification and Characterization of Novel Genes from Legumes

The Fabaceae, the third largest family of plants and the source of many crops, has been the target of many genomic studies. Currently, only the grasses surpass the legumes for the number of publicly available expressed sequence tags (ESTs). The quantity of sequences from diverse plants enables the use of computational approaches to identify novel genes in specific taxa. We used BLAST algorithms to compare unigene sets from Medicago truncatula, Lotus japonicus, and soybean (Glycine max and Glycine soja) to nonlegume unigene sets, to GenBanks nonredundant and EST databases, and to the genomic sequences of rice (Oryza sativa) and Arabidopsis. As a working definition, putatively legume-specific genes had no sequence homology, below a specified threshold, to publicly available sequences of nonlegumes. Using this approach, 2,525 legume-specific EST contigs were identified, of which less than three percent had clear homology to previously characterized legume genes. As a first step toward predicting function, related sequences were clustered to build motifs that could be searched against protein databases. Three families of interest were more deeply characterized: F-box related proteins, Pro-rich proteins, and Cys cluster proteins (CCPs). Of particular interest were the >300 CCPs, primarily from nodules or seeds, with predicted similarity to defensins. Motif searching also identified several previously unknown CCP-like open reading frames in Arabidopsis. Evolutionary analyses of the genomic sequences of several CCPs in M. truncatula suggest that this family has evolved by local duplications and divergent selection.

Sequencing and Analysis of Common Bean ESTs. Building a Foundation for Functional Genomics

Although common bean (Phaseolus vulgaris) is the most important grain legume in the developing world for human consumption, few genomic resources exist for this species. The objectives of this research were to develop expressed sequence tag (EST) resources for common bean and assess nodule gene expression through high-density macroarrays. We sequenced a total of 21,026 ESTs derived from 5 different cDNA libraries, including nitrogen-fixing root nodules, phosphorus-deficient roots, developing pods, and leaves of the Mesoamerican genotype, Negro Jamapa 81. The fifth source of ESTs was a leaf cDNA library derived from the Andean genotype, G19833. Of the total high-quality sequences, 5,703 ESTs were classified as singletons, while 10,078 were assembled into 2,226 contigs producing a nonredundant set of 7,969 different transcripts. Sequences were grouped according to 4 main categories, metabolism (34%), cell cycle and plant development (11%), interaction with the environment (19%), and unknown function (36%), and further subdivided into 15 subcategories. Comparisons to other legume EST projects suggest that an entirely different repertoire of genes is expressed in common bean nodules. Phaseolus-specific contigs, gene families, and single nucleotide polymorphisms were also identified from the EST collection. Functional aspects of individual bean organs were reflected by the 20 contigs from each library composed of the most redundant ESTs. The abundance of transcripts corresponding to selected contigs was evaluated by RNA blots to determine whether gene expression determined by laboratory methods correlated with in silico expression. Evaluation of root nodule gene expression by macroarrays and RNA blots showed that genes related to nitrogen and carbon metabolism are integrated for ureide production. Resources developed in this project provide genetic and genomic tools for an international consortium devoted to bean improvement.

Distinct Biphasic mRNA Changes in Response to Asian Soybean Rust Infection

Asian soybean rust (ASR), caused by Phakopsora pachyrhizi, is now established in all major soybean-producing countries. Currently, there is little information about the molecular basis of ASR-soybean interactions, which will be needed to assist future efforts to develop effective resistance. Toward this end, abundance changes of soybean mRNAs were measured over a 7-day ASR infection time course in mock-inoculated and infected leaves of a soybean accession (PI230970) carrying the Rpp2 resistance gene and a susceptible genotype (Embrapa-48). The expression profiles of differentially expressed genes (ASR-infected compared with the mock-inoculated control) revealed a biphasic response to ASR in each genotype. Within the first 12 h after inoculation (hai), which corresponds to fungal germination and penetration of the epidermal cells, differential gene expression changes were evident in both genotypes. mRNA expression of these genes mostly returned to levels found in mock-inoculated plants by 24 hai. In the susceptible genotype, gene expression remained unaffected by rust infection until 96 hai, a time period when rapid fungal growth began. In contrast, gene expression in the resistant genotype diverged from the mock-inoculated control earlier, at 72 h, demonstrating that Rpp2-mediated defenses were initiated prior to this time. These data suggest that ASR initially induces a nonspecific response that is transient or is suppressed when early steps in colonization are completed in both soybean genotypes. The race-specific resistance phenotype of Rpp2 is manifested in massive gene expression changes after the initial response prior to the onset of rapid fungal growth that occurs in the susceptible genotype.

Identification and Analyses of Candidate Genes for Rpp4-Mediated Resistance to Asian Soybean Rust in Soybean

Asian soybean rust is a formidable threat to soybean (Glycine max) production in many areas of the world, including the United States. Only five sources of resistance have been identified (Resistance to Phakopsora pachyrhizi1 [Rpp1], Rpp2, Rpp3, Rpp4, and Rpp5). Rpp4 was previously identified in the resistant genotype PI459025B and mapped within 2 centimorgans of Satt288 on soybean chromosome 18 (linkage group G). Using simple sequence repeat markers, we developed a bacterial artificial chromosome contig for the Rpp4 locus in the susceptible cv Williams82 (Wm82). Sequencing within this region identified three Rpp4 candidate disease resistance genes (Rpp4C1–Rpp4C3 [Wm82]) with greatest similarity to the lettuce (Lactuca sativa) RGC2 family of coiled coil-nucleotide binding site-leucine rich repeat disease resistance genes. Constructs containing regions of the Wm82 Rpp4 candidate genes were used for virus-induced gene silencing experiments to silence resistance in PI459025B, confirming that orthologous genes confer resistance. Using primers developed from conserved sequences in the Wm82 Rpp4 candidate genes, we identified five Rpp4 candidate genes (Rpp4C1–Rpp4C5 [PI459025B]) from the resistant genotype. Additional markers developed from the Wm82 Rpp4 bacterial artificial chromosome contig further defined the region containing Rpp4 and eliminated Rpp4C1 (PI459025B) and Rpp4C3 (PI459025B) as candidate genes. Sequencing of reverse transcription-polymerase chain reaction products revealed that Rpp4C4 (PI459025B) was highly expressed in the resistant genotype, while expression of the other candidate genes was nearly undetectable. These data support Rpp4C4 (PI459025B) as the single candidate gene for Rpp4-mediated resistance to Asian soybean rust.

Soybean homologs of MPK4 negatively regulate defense responses and positively regulate growth and development.

Mitogen-activated protein kinase (MAPK) cascades play important roles in disease resistance in model plant species such as Arabidopsis (Arabidopsis thaliana) and tobacco (Nicotiana tabacum). However, the importance of MAPK signaling pathways in the disease resistance of crops is still largely uninvestigated. To better understand the role of MAPK signaling pathways in disease resistance in soybean (Glycine max), 13, nine, and 10 genes encoding distinct MAPKs, MAPKKs, and MAPKKKs, respectively, were silenced using virus-induced gene silencing mediated by Bean pod mottle virus. Among the plants silenced for various MAPKs, MAPKKs, and MAPKKKs, those in which GmMAPK4 homologs (GmMPK4s) were silenced displayed strong phenotypes including stunted stature and spontaneous cell death on the leaves and stems, the characteristic hallmarks of activated defense responses. Microarray analysis showed that genes involved in defense responses, such as those in salicylic acid (SA) signaling pathways, were significantly up-regulated in GmMPK4-silenced plants, whereas genes involved in growth and development, such as those in auxin signaling pathways and in cell cycle and proliferation, were significantly down-regulated. As expected, SA and hydrogen peroxide accumulation was significantly increased in GmMPK4-silenced plants. Accordingly, GmMPK4-silenced plants were more resistant to downy mildew and Soybean mosaic virus compared with vector control plants. Using bimolecular fluorescence complementation analysis and in vitro kinase assays, we determined that GmMKK1 and GmMKK2 might function upstream of GmMPK4. Taken together, our results indicate that GmMPK4s negatively regulate SA accumulation and defense response but positively regulate plant growth and development, and their functions are conserved across plant species.

Complementary genetic and genomic approaches help characterize the linkage group I seed protein QTL in soybean

BackgroundThe nutritional and economic value of many crops is effectively a function of seed protein and oil content. Insight into the genetic and molecular control mechanisms involved in the deposition of these constituents in the developing seed is needed to guide crop improvement. A quantitative trait locus (QTL) on Linkage Group I (LG I) of soybean (Glycine max (L.) Merrill) has a striking effect on seed protein content.ResultsA soybean near-isogenic line (NIL) pair contrasting in seed protein and differing in an introgressed genomic segment containing the LG I protein QTL was used as a resource to demarcate the QTL region and to study variation in transcript abundance in developing seed. The LG I QTL region was delineated to less than 8.4 Mbp of genomic sequence on chromosome 20. Using Affymetrix® Soy GeneChip and high-throughput Illumina® whole transcriptome sequencing platforms, 13 genes displaying significant seed transcript accumulation differences between NILs were identified that mapped to the 8.4 Mbp LG I protein QTL region.ConclusionsThis study identifies gene candidates at the LG I protein QTL for potential involvement in the regulation of protein content in the soybean seed. The results demonstrate the power of complementary approaches to characterize contrasting NILs and provide genome-wide transcriptome insight towards understanding seed biology and the soybean genome.

Recent Advances in Legume-Microbe Interactions: Recognition, Defense Response, and Symbiosis from a Genomic Perspective

The ability of legumes to form symbiotic mutualistic relationships with certain bacteria in the Rhizobiales (collectively called rhizobia) and harness the ability of the bacteria to fix atmospheric N2 into ammonia has a tremendous impact on natural and agricultural ecosystems. The interaction

Defensin gene family in Medicago truncatula: structure, expression and induction by signal molecules.

A large gene family encoding the putative cysteine-rich defensins was discovered in Medicago truncatula. Sixteen members of the family were identified by screening a cloned seed defensin from M. sativa (Gao et al. 2000) against the Institute for Genomic Research’s (TIGR) M. truncatula gene index (MtGI version 7). Based on the comparison of their amino acid sequences, M. truncatula defensins fell arbitrarily into three classes displaying extensive sequence divergence outside of the eight canonical cysteine residues. The presence of Class II defensins is reported for the first time in a legume plant. In silico as well as Northern blot and RT-PCR analyses indicated these genes were expressed in a variety of tissues including leaves, flowers, developing pods, mature seed and roots. The expression of these genes was differentially induced in response to a variety of biotic and abiotic stimuli. For the first time, a defensin gene (TC77480) was shown to be induced in roots in response to infection by the mycorrhizal fungus, Glomus versiforme. Northern blot analysis indicated that the tissue-specific expression patterns of the cloned Def1 and Def2 genes differed substantially between M. truncatula and M. sativa. Furthermore, the induction profiles of the Def1 and Def2 genes in response to the signaling molecules methyl jasmonate, ethylene and salicylic acid differed markedly between these two legumes.