Randy C. Shoemaker

Sequencing crop genomes: approaches and applications

Many challenges face plant scientists, in particular those working on crop production, such as a projected increase in population, decrease in water and arable land, changes in weather patterns and predictability. Advances in genome sequencing and resequencing can and should play a role in our response to meeting these challenges. However, several barriers prevent rapid and effective deployment of these tools to a wide variety of crops. Because of the complexity of crop genomes, de novo sequencing with next-generation sequencing technologies is a process fraught with difficulties that then create roadblocks to the utilization of these genome sequences for crop improvement. Collecting rapid and accurate phenotypes in crop plants is a hindrance to integrating genomics with crop improvement, and advances in informatics are needed to put these tools in the hands of the scientists on the ground.

Bifurcation and Enhancement of Autonomous-Nonautonomous Retrotransposon Partnership through LTR Swapping in Soybean

This work demonstrates that region-specific interelement recombinational exchange, behind natural selection, plays a primary role in maintaining preexisting partnership and establishing new partnership between nonautonomous and autonomous long terminal repeat retrotransposons in soybean. Long terminal repeat (LTR) retrotransposons, the most abundant genomic components in flowering plants, are classifiable into autonomous and nonautonomous elements based on their structural completeness and transposition capacity. It has been proposed that selection is the major force for maintaining sequence (e.g., LTR) conservation between nonautonomous elements and their autonomous counterparts. Here, we report the structural, evolutionary, and expression characterization of a giant retrovirus-like soybean (Glycine max) LTR retrotransposon family, SNARE. This family contains two autonomous subfamilies, SAREA and SAREB, that appear to have evolved independently since the soybean genome tetraploidization event ∼13 million years ago, and a nonautonomous subfamily, SNRE, that originated from SAREA. Unexpectedly, a subset of the SNRE elements, which amplified from a single founding SNRE element within the last ∼3 million years, have been dramatically homogenized with either SAREA or SAREB primarily in the LTR regions and bifurcated into distinct subgroups corresponding to the two autonomous subfamilies. We uncovered evidence of region-specific swapping of nonautonomous elements with autonomous elements that primarily generated various nonautonomous recombinants with LTR sequences from autonomous elements of different evolutionary lineages, thus revealing a molecular mechanism for the enhancement of preexisting partnership and the establishment of new partnership between autonomous and nonautonomous elements.

Transcriptome analyses and virus induced gene silencing identify genes in the Rpp4-mediated Asian soybean rust resistance pathway

Rpp4 (Resistance to Phakopsora pachyrhizi 4) confers resistance to Phakopsora pachyrhizi Sydow, the causal agent of Asian soybean rust (ASR). By combining expression profiling and virus induced gene silencing (VIGS), we are developing a genetic framework for Rpp4-mediated resistance. We measured gene expression in mock-inoculated and P. pachyrhizi-infected leaves of resistant soybean accession PI459025B (Rpp4) and the susceptible cultivar (Williams 82) across a 12-day time course. Unexpectedly, two biphasic responses were identified. In the incompatible reaction, genes induced at 12h after infection (hai) were not differentially expressed at 24 hai, but were induced at 72 hai. In contrast, genes repressed at 12 hai were not differentially expressed from 24 to 144 hai, but were repressed 216 hai and later. To differentiate between basal and resistance-gene (R-gene) mediated defence responses, we compared gene expression in Rpp4-silenced and empty vector-treated PI459025B plants 14 days after infection (dai) with P. pachyrhizi. This identified genes, including transcription factors, whose differential expression is dependent upon Rpp4. To identify differentially expressed genes conserved across multiple P. pachyrhizi resistance pathways, Rpp4 expression datasets were compared with microarray data previously generated for Rpp2 and Rpp3-mediated defence responses. Fourteen transcription factors common to all resistant and susceptible responses were identified, as well as fourteen transcription factors unique to R-gene-mediated resistance responses. These genes are targets for future P. pachyrhizi resistance research.

MAPPING OF IRON AND ZINC QUANTITATIVE TRAIT LOCI IN SOYBEAN FOR ASSOCIATION TO IRON DEFICIENCY CHLOROSIS RESISTANCE

Iron deficiency chlorosis (IDC) in soybean results in yield losses or in extreme cases death. Breeding for resistance has shown limited success with no cultivar having complete resistance. Mineral content of the soybean could be an indicator of the ability of the plant to withstand the effects of IDC. Iron (Fe) and zinc (Zn) concentration was examined in soybean seed and leaves. SSR, RFLP, and BARCSOYSSR markers were used to construct a linkage map used for mapping of Fe and Zn concentrations. The QTL analysis for the combined data identified one major QTL for seed Fe accumulation on chromosome 20 that explained 21.5% of the variation. This QTL was in the marker interval pa_515-1-Satt239, with marker pa_515-1 previously being used to map an Fe-efficiency QTL. This provides the first evidence of a potential genetic link between Fe-efficiency and Fe accumulation in the soybean seed.

Applications and methods utilizing the Simple Semantic Web Architecture and Protocol (SSWAP) for bioinformatics resource discovery and disparate data and service integration

BackgroundScientific data integration and computational service discovery are challenges for the bioinformatic community. This process is made more difficult by the separate and independent construction of biological databases, which makes the exchange of data between information resources difficult and labor intensive. A recently described semantic web protocol, the Simple Semantic Web Architecture and Protocol (SSWAP; pronounced swap) offers the ability to describe data and services in a semantically meaningful way. We report how three major information resources (Gramene, SoyBase and the Legume Information System [LIS]) used SSWAP to semantically describe selected data and web services.MethodsWe selected high-priority Quantitative Trait Locus (QTL), genomic mapping, trait, phenotypic, and sequence data and associated services such as BLAST for publication, data retrieval, and service invocation via semantic web services. Data and services were mapped to concepts and categories as implemented in legacy and de novo community ontologies. We used SSWAP to express these offerings in OWL Web Ontology Language (OWL), Resource Description Framework (RDF) and eXtensible Markup Language (XML) documents, which are appropriate for their semantic discovery and retrieval. We implemented SSWAP services to respond to web queries and return data. These services are registered with the SSWAP Discovery Server and are available for semantic discovery at http://sswap.info.ResultsA total of ten services delivering QTL information from Gramene were created. From SoyBase, we created six services delivering information about soybean QTLs, and seven services delivering genetic locus information. For LIS we constructed three services, two of which allow the retrieval of DNA and RNA FASTA sequences with the third service providing nucleic acid sequence comparison capability (BLAST).ConclusionsThe need for semantic integration technologies has preceded available solutions. We report the feasibility of mapping high priority data from local, independent, idiosyncratic data schemas to common shared concepts as implemented in web-accessible ontologies. These mappings are then amenable for use in semantic web services. Our implementation of approximately two dozen services means that biological data at three large information resources (Gramene, SoyBase, and LIS) is available for programmatic access, semantic searching, and enhanced interaction between the separate missions of these resources.

ESTminer: a suite of programs for gene and allele identification

UNLABELLEDnESTminer is a collection of programs that use expressed sequence tag (EST) data from inbred genomes to identify unique genes within gene families. The algorithm utilizes Cap3 to perform an initial clustering of related EST sequences to produce a consensus sequence of a gene family. These consensus sequences are then used to collect all ESTs in the original EST library that are related using BLAST. A redundancy based criterion is applied to each EST to identify reliable unique gene-sequences. Using a highly inbred genome as a source of ESTs eliminates the necessity of computing covariance on each polymorphism to identify alleles of the same gene, thus making this algorithm more streamlined than other alternatives which must computationally attempt to distinguish genes from alleles.nnnAVAILABILITYnThe programs were written in PERL and are freely available at http://www.soybase.org/publication_data/Nelson/ESTminer/ESTminer.htmlnnnCONTACTnnelsonrt@iastate.edunnnSUPPLEMENTARY INFORMATIONnFigures and dataset can be obtained from: http://www.soybase.org/publication_data/Nelson/ESTminer/ESTminer.html.

Evaluation and QTL mapping of phosphorus concentration in soybean seed

Phosphorus (P) is an essential macronutrient required for many biological and metabolic plant functions. Phosphate is the form of P used by plants. Reducing the amount of P in the seed can further aide breeding efforts to reduce the amount of P released into the environment due to the indigestibility of the complexes formed with metal ions. Analysis of the variation of phosphorus concentration in soybean seed under non-stressed conditions revealed that phosphorus ranged from 3,948.1 to 5,695.8xa0μg/g total phosphorous (TP) in combined years. The averages for independent years were significantly different from one another. Quantitative trait loci (QTL) analysis of TP was performed to identify candidate gene(s) that is (are) involved in P accumulation in soybean seed. One putative QTL region was identified on chromosome 12 in the combined data that contained a phosphate transporter gene. Two additional suggestive QTL were identified on chromosomes 7 and 17 with chromosome 7 having both a phosphate transport gene and a ZIP transporter gene in the region of the QTL. There were additional genes in these regions that are involved in phosphate metabolism and transport.

Identification of candidate signaling genes including regulators of chromosome condensation 1 protein family differentially expressed in the soybean–Phytophthora sojae interaction

Stem and root rot caused by the oomycete pathogen, Phytophthora sojae, is a serious soybean disease. Use of Phytophthora resistance genes (Rps) in soybean cultivars has been very effective in controlling this pathogen. Resistance encoded by Rps genes is manifested through activation of defense responses. In order to identify candidate signaling genes involved in the expression of Phytophthora resistance in soybean, a cDNA library was prepared from infected etiolated hypocotyl tissues of a Phytophthora resistant soybean cultivar harvested 2 and 4xa0h following P. sojae inoculation. In silico subtraction of 101,833 expressed sequence tags (ESTs) originating from unstressed cDNA libraries from 4,737 ESTs of this library resulted in identification of 204 genes that were absent in the unstressed libraries. Of the 204 identified genes, seven were P. sojae genes. Putative function of 91 of the 204 genes could not be assigned based on sequence comparison. Macroarray analyses of all 204 genes led to identification of 60 genes including 15 signaling-related soybean genes and three P. sojae genes, transcripts of which were induced twofold in P. sojae-infected tissues as compared to that in water controls. Eight soybean genes were down-regulated twofold following P. sojae infection as compared to water controls. Differential expression of a few selected genes was confirmed by conducting Northern and RT-PCR analyses. We have shown that two putative regulators of chromosome condensation 1 (RCC1) family proteins were down-regulated in the incompatible interaction. This observation suggested that the nucleocytoplasmic transport function for trafficking protein and non-coding RNA is suppressed during expression of race-specific Phytophthora resistance. Characterization of a cDNA library generated from tissues harvested almost immediately following P. sojae-infection of a resistant cultivar allowed us to identify many candidate signaling genes that are presumably involved in regulating the expression of defense-related pathways for expression of Phytophthora resistance in soybean.

Gene expression: sizing it all up.

Genomic architecture appears to be a largely unexplored component of gene expression. That architecture can be related to chromatin domains, transposable element neighborhoods, epigenetic modifications of the genome, and more. Although surely not the end of the story, we are learning that when it comes to gene expression, size is also important. We have been surprised to find that certain patterns of expression, tissue specific versus constitutive, or high expression versus low expression, are often associated with physical attributes of the gene and genome. Multiple studies have shown an inverse relationship between gene expression patterns and various physical parameters of the genome such as intron size, exon size, intron number, and size of intergenic regions. An increase in expression level and breadth often correlates with a decrease in the size of physical attributes of the gene. Three models have been proposed to explain these relationships. Contradictory results were found in several organisms when expression level and expression breadth were analyzed independently. However, when both factors were combined in a single study a novel relationship was revealed. At low levels of expression, an increase in expression breadth correlated with an increase in genic, intergenic, and intragenic sizes. Contrastingly, at high levels of expression, an increase in expression breadth inversely correlated with the size of the gene. In this article we explore the several hypotheses regarding genome physical parameters and gene expression.

Fertile progeny of a hybridization between soybean [Glycine max (L.) Merr.] and G. tomentella Hayata

SummaryA colchicine-doubled F1 hybrid (2n=118) of a cross between PI 360841 (Glycine max) (2n=40) x PI 378708 (G. tomentella) (2n=78), propagated by shoot cuttings since January 1984, produced approximately 100 F2 seed during October 1988. One-fourth of the F2 plants or their F3 progeny have been analyzed for chromosome number, pollen viability, pubescence tip morphology, seed coat color, and isoenzyme variation. Without exception, all plants evaluated possessed the chromosome number of the G. max parent (2n=40). Most F2 plants demonstrated a high level of fertility, although 2 of 24 plants had low pollen viability and had large numbers of fleshy pods. One F2 plant possessed sharp pubescence tip morphology, whereas all others were blunt-tipped. All evaluated F2 and F3 plants expressed the malate dehydrogenase and diaphorase isoenzyme patterns of the G. max parent and the endopeptidase isoenzyme pattern of the G. tomentella parent. Mobility variants were observed among progeny for the isoenzymes phosphoglucomutase, aconitase, and phosphoglucoisomerase. This study suggests that the G. Tomentella chromosome complement has been eliminated after genetic exchange and/or modification has taken place between the genomes.