The Soybean Genome Database (SoyGD): a browser for display of duplicated, polyploid, regions and sequence tagged sites on the integrated physical and genetic maps of Glycine max

TML1 AND TML2 SYNERGISTICALLY REGULATE NODULATION BUT NOT ARBUSCULAR MYCORRHIZA IN MEDICAGO TRUNCATULA

Two symbiotic processes, nodulation and arbuscular mycorrhiza, are primarily controlled by the plant's need for nitrogen (N) and phosphorus (P), respectively. Autoregulation of Nodulation (AON) and Autoregulation of Mycorrhization (AOM) share multiple components - plants that make too many nodules usually have higher arbuscule density. The protein TML (TOO MUCH LOVE) was shown to function in roots to maintain susceptibly to rhizobial infection under low N conditions and control nodule number through AON in Lotus japonicus. M. truncatula has two sequence homologs: MtTML1 and MtTML2. We report the generation of stable single and double mutants harboring multiple allelic variations in MtTML1 and MtTML2 using CRISPR-Cas9 targeted mutagenesis and screening of a transposon mutagenesis library. Plants containing single mutations in either gene produced twice the nodules of wild type plants whereas plants containing mutations in both genes displayed a synergistic effect, forming 20x more nodules and short roots compared to wild type plants. The synergistic effect on nodulation was maintained in the presence of 10mM nitrogen, but not observed in root length phenotypes. Examination of expression and heterozygote effects suggest genetic compensation may play a role in the observed synergy. However, plants with mutations in both TMLs had no detectable change in arbuscular mycorrhizal associations, suggesting that MtTMLs are specific to nodulation and nitrate signaling. The mutants created will be useful tools to dissect the mechanism of synergistic action of MtTML1 and MtTML2 in M. truncatula nodulation as well as the separation of AON from AOM.

Bridging the Gap: Combining Genomics and Transcriptomics Approaches to Understand Stylosanthes scabra, an Orphan Legume from the Brazilian Caatinga

Stylosanthes scabra is a scientifically orphaned legume found in the Brazilian Caatinga biome (a semi-arid environment). This work utilized omics approaches to investigate some ecophysiological aspects of stress tolerance/resistance in S. scabra, study its genomic landscape, and predict potential metabolic pathways. Considering its high-confidence conceptual proteome, 1694 (~2.6%) proteins were associated with resistance proteins, some of which were found in soybean QTL regions that confer resistance to Asian soybean rust. S. scabra was also found to be a potential source of terpenes, as biosynthetic gene clusters associated with terpene biosynthesis were identified in its genome. The analysis revealed that mobile elements comprised approximately 59% of the sequenced genome. In the remaining 41% of the sections, some of the 22,681 protein-coding gene families were categorized into two informational groups: those that were specific to S. scabra and those that expanded significantly compared to their immediate ancestor. Biological process enrichment analyses indicated that these gene families play fundamental roles in the adaptation of S. scabra to extreme environments. Additionally, phylogenomic analysis indicated a close evolutionary relationship between the genera Stylosanthes and Arachis. Finally, this study found a high number (57) of aquaporin-encoding loci in the S. scabra genome. RNA-Seq and qPCR data suggested that the PIP subfamily may play a key role in the species' adaptation to water deficit conditions. Overall, these results provide valuable insights into S. scabra biology and a wealth of gene/transcript information for future legume omics studies.

A soybean non-coding RNA mining and co-expression resource based on 1,596 RNA-seq and small RNA-seq libraries

The SoyNcRNAExp soybean non-coding RNA expression/co-expression resource can be used for ncRNA expression, mining, and co-expression analysis.

High-Resolution Translatome Analysis Reveals Cortical Cell Programs During Early Soybean Nodulation

Nodule organogenesis in legumes is regulated temporally and spatially through gene networks. Genome-wide transcriptome, proteomic, and metabolomic analyses have been used previously to define the functional role of various plant genes in the nodulation process. However, while significant progress has been made, most of these studies have suffered from tissue dilution since only a few cells/root regions respond to rhizobial infection, with much of the root non-responsive. To partially overcome this issue, we adopted translating ribosome affinity purification (TRAP) to specifically monitor the response of the root cortex to rhizobial inoculation using a cortex-specific promoter. While previous studies have largely focused on the plant response within the root epidermis (e.g., root hairs) or within developing nodules, much less is known about the early responses within the root cortex, such as in relation to the development of the nodule primordium or growth of the infection thread. We focused on identifying genes specifically regulated during early nodule organogenesis using roots inoculated with Bradyrhizobium japonicum. A number of novel nodulation gene candidates were discovered, as well as soybean orthologs of nodulation genes previously reported in other legumes. The differential cortex expression of several genes was confirmed using a promoter-GUS analysis, and RNAi was used to investigate gene function. Notably, a number of differentially regulated genes involved in phytohormone signaling, including auxin, cytokinin, and gibberellic acid (GA), were also discovered, providing deep insight into phytohormone signaling during early nodule development.

Gene duplication and rate variation in the evolution of plastid ACCase and Clp genes in angiosperms

While the chloroplast (plastid) is known for its role in photosynthesis, it is also involved in many other metabolic pathways essential for plant survival. As such, plastids contain an extensive suite of enzymes required for non-photosynthetic processes. The evolution of the associated genes has been especially dynamic in flowering plants (angiosperms), including examples of gene duplication and extensive rate variation. We examined the role of ongoing gene duplication in two key plastid enzymes, the acetyl-CoA carboxylase (ACCase) and the caseinolytic protease (Clp), responsible for fatty acid biosynthesis and protein turnover, respectively. In plants, there are two ACCase complexes-a homomeric version present in the cytosol and a heteromeric version present in the plastid. Duplications of the nuclear-encoded homomeric ACCase gene and retargeting of one resultant protein to the plastid have been previously reported in multiple species. We find that these retargeted homomeric ACCase proteins exhibit elevated rates of sequence evolution, consistent with neofunctionalization and/or relaxation of selection. The plastid Clp complex catalytic core is composed of nine paralogous proteins that arose via ancient gene duplication in the cyanobacterial/plastid lineage. We show that further gene duplication occurred more recently in the nuclear-encoded core subunits of this complex, yielding additional paralogs in many species of angiosperms. Moreover, in six of eight cases, subunits that have undergone recent duplication display increased rates of sequence evolution relative to those that have remained single copy. We also compared substitution patterns between pairs of Clp core paralogs to gain insight into post-duplication evolutionary routes. These results show that gene duplication and rate variation continue to shape the plastid proteome.

Genome-Wide Association Study of Soybean Germplasm Derived From Canadian × Chinese Crosses to Mine for Novel Alleles to Improve Seed Yield and Seed Quality Traits

Genome-wide association study (GWAS) has emerged in the past decade as a viable tool for identifying beneficial alleles from a genomic diversity panel. In an ongoing effort to improve soybean [Glycine max (L.) Merr.], which is the third largest field crop in Canada, a GWAS was conducted to identify novel alleles underlying seed yield and seed quality and agronomic traits. The genomic panel consisted of 200 genotypes including lines derived from several generations of bi-parental crosses between modern Canadian × Chinese cultivars (CD-CH). The genomic diversity panel was field evaluated at two field locations in Ontario in 2019 and 2020. Genotyping-by-sequencing (GBS) was conducted and yielded almost 32 K high-quality SNPs. GWAS was conducted using Fixed and random model Circulating Probability Unification (FarmCPU) model on the following traits: seed yield, seed protein concentration, seed oil concentration, plant height, 100 seed weight, days to maturity, and lodging score that allowed to identify five QTL regions controlling seed yield and seed oil and protein content. A candidate gene search identified a putative gene for each of the three traits. The results of this GWAS study provide insight into potentially valuable genetic resources residing in Chinese modern cultivars that breeders may use to further improve soybean seed yield and seed quality traits.

Genome-wide identification and characterization of legume T2 Ribonuclease gene family and analysis of GmaRNS9, a soybean T2 Ribonuclease gene, function in nodulation

T2 ribonuclease family (RNaseT2) proteins are secretory and nonspecific endoribonucleases that have a large conserved biological role. Family members of RNaseT2 are found in every organism and carry out important biological functions. However, little is known about the functions of these proteins in legumes, including potential roles in symbiotic nodulation. This study aimed to characterize and perform bioinformatic analysis of RNaseT2 genes in four legume species that their genome was sequenced. In total, 60 RNaseT2 genes were identified and characterized. By analyzing their phylogeny, we divided these RNaseT2 into five clades. Expression analysis of RNaseT2 genes indicated that these genes are expressed in various tissues, and the most expression level was related to the pod, flower, and root. Moreover, GmaRNS9 expression analysis in soybean was consistent with in silico studies and demonstrated that this gene usually has high root tip expression. GmaRNS9 expression was reduced by Bradyrhizobium japonicum inoculation and nodule formation. Reduced expression of this gene was possibly controlled by the GmNARK gene either directly or pleiotropically through increased phosphorus requirements during increased nodulation. However, the nutrient stress (phosphate and nitrate starvation) led to an increase in the expression level of GmRNS9. In silico and quantitative gene expression analyses showed that RNaseT2 genes could play important roles in the growth and development of legumes as well as nodulation.

The Soybean High Density 'Forrest' by 'Williams 82' SNP-Based Genetic Linkage Map Identifies QTL and Candidate Genes for Seed Isoflavone Content

Isoflavones are secondary metabolites that are abundant in soybean and other legume seeds providing health and nutrition benefits for both humans and animals. The objectives of this study were to construct a single nucleotide polymorphism (SNP)-based genetic linkage map using the ‘Forrest’ by ‘Williams 82’ (F×W82) recombinant inbred line (RIL) population (n = 306); map quantitative trait loci (QTL) for seed daidzein, genistein, glycitein, and total isoflavone contents in two environments over two years (NC-2018 and IL-2020); identify candidate genes for seed isoflavone. The FXW82 SNP-based map was composed of 2075 SNPs and covered 4029.9 cM. A total of 27 QTL that control various seed isoflavone traits have been identified and mapped on chromosomes (Chrs.) 2, 4, 5, 6, 10, 12, 15, 19, and 20 in both NC-2018 (13 QTL) and IL-2020 (14 QTL). The six QTL regions on Chrs. 2, 4, 5, 12, 15, and 19 are novel regions while the other 21 QTL have been identified by other studies using different biparental mapping populations or genome-wide association studies (GWAS). A total of 130 candidate genes involved in isoflavone biosynthetic pathways have been identified on all 20 Chrs. And among them 16 have been identified and located within or close to the QTL identified in this study. Moreover, transcripts from four genes (Glyma.10G058200, Glyma.06G143000, Glyma.06G137100, and Glyma.06G137300) were highly abundant in Forrest and Williams 82 seeds. The identified QTL and four candidate genes will be useful in breeding programs to develop soybean cultivars with high beneficial isoflavone contents.

Genome-wide association analyses reveal the genetic basis of biomass accumulation under symbiotic nitrogen fixation in African soybean

KEY MESSAGE

We explored the genetic basis of SNF-related traits through GWAS and identified 40 candidate genes. This study provides fundamental insights into SNF-related traits and will accelerate efforts for SNF breeding. Symbiotic nitrogen fixation (SNF) increases sustainability by supplying biological nitrogen for crops to enhance yields without damaging the ecosystem. A better understanding of this complex biological process is critical for addressing the triple challenges of food security, environmental degradation, and climate change. Soybean plants, the most important legume worldwide, can form a mutualistic interaction with specialized soil bacteria, bradyrhizobia, to fix atmospheric nitrogen. Here we report a comprehensive genome-wide association study of 11 SNF-related traits using 79K GBS-derived SNPs in 297 African soybeans. We identified 25 QTL regions encompassing 40 putative candidate genes for SNF-related traits including 20 genes with no prior known role in SNF. A line with a large deletion (164 kb), encompassing a QTL region containing a strong candidate gene (CASTOR), exhibited a marked decrease in SNF. This study performed on African soybean lines provides fundamental insights into SNF-related traits and yielded a rich catalog of candidate genes for SNF-related traits that might accelerate future efforts aimed at sustainable agriculture.

Characterization of the Spatial and Temporal Expression of Two Soybean miRNAs Identifies SCL6 as a Novel Regulator of Soybean Nodulation

MicroRNAs (miRNAs) control expression of endogenous target genes through transcript cleavage or translational inhibition. Legume plants can form a specialized organ, the nodule, through interaction with nitrogen fixing soil bacteria. To understand the regulatory roles of miRNAs in the nodulation process, we functionally validated gma-miR171o and gma-miR171q and their target genes in soybean. These two miRNA sequences are identical in sequence but their miRNA genes are divergent and show unique, tissue-specific expression patterns. The expression levels of the two miRNAs are negatively correlated with that of their target genes. Ectopic expression of these miRNAs in transgenic hairy roots resulted in a significant reduction in nodule formation. Both gma-miR171o and gma-miR171q target members of the GRAS transcription factor superfamily, namely GmSCL-6 and GmNSP2. Transient interaction of miRNAs and their target genes in tobacco cells further confirmed their cleavage activity. The results suggest that gma-miR171o and gma-miR171q regulate GmSCL-6 and GmNSP2, which in turn, influence expression of the Nodule inception (NIN), Early Nodulin 40 (ENOD40), and Ethylene Response Factor Required for Nodulation (ERN) genes during the Bradyrhizobium japonicum-soybean nodulation process. Collectively, our data suggest a role for two miRNAs, gma-miR171o and gma-miR171q, in regulating the spatial and temporal aspects of soybean nodulation.

In Silico Genome-Wide Analysis of the ATP-Binding Cassette Transporter Gene Family in Soybean (Glycine max L.) and Their Expression Profiling

ATP-binding cassette (ABC) transporters constitute one of the largest gene families in all living organisms, most of which mediate transport across biological membranes by hydrolyzing ATP. However, detailed studies of ABC transporter genes in the important oil crop, soybean, are still lacking. In the present study, we carried out genome-wide identification and phylogenetic and transcriptional analyses of the ABC gene family in G. max. A total of 261 G. max ABC (GmABCs) genes were identified and unevenly localized onto 20 chromosomes. Referring to protein-domain orientation and phylogeny, the GmABC family could be classified into eight (ABCA-ABCG and ABCI) subfamilies and ABCG were the most abundantly present. Further, investigation of whole genome duplication (WGD) signifies the role of segmental duplication in the expansion of the ABC transporter gene family in soybean. The Ka/Ks ratio indicates that several duplicated genes are governed by intense purifying selection during evolution. In addition, in silico expression analysis based on RNA-sequence using publicly available database revealed that ABC transporters are differentially expressed in tissues and developmental stages and in dehydration. Overall, we provide an extensive overview of the GmABC transporter gene family and it promises the primary basis for the study in development and response to dehydration tolerance.

Regulated Disorder: Posttranslational Modifications Control the RIN4 Plant Immune Signaling Hub

RIN4 is an intensively studied immune regulator in Arabidopsis and is involved in perception of microbial features outside and bacterial effectors inside plant cells. Furthermore, RIN4 is conserved in land plants and is targeted for posttranslational modifications by several virulence proteins from the bacterial pathogen Pseudomonas syringae. Despite the important roles of RIN4 in plant immune responses, its molecular function is not known. RIN4 is an intrinsically disordered protein (IDP), except at regions where pathogen-induced posttranslational modifications take place. IDP act as hubs for protein complex formation due to their ability to bind to multiple client proteins and, thus, are important players in signal transduction pathways. RIN4 is known to associate with multiple proteins involved in immunity, likely acting as an immune-signaling hub for the formation of distinct protein complexes. Genetically, RIN4 is a negative regulator of immunity, but diverse posttranslational modifications can either enhance its negative regulatory function or, on the contrary, render it a potent immune activator. In this review, we describe the structural domains of RIN4 proteins, their intrinsically disordered regions, posttranslational modifications, and highlight the implications that these features have on RIN4 function. In addition, we will discuss the potential role of plasma membrane subdomains in mediating RIN4 protein complex formations.

The Evolution of Soybean Knowledge Base (SoyKB)

Soybean Knowledge Base (SoyKB) is a comprehensive all-inclusive web resource for bridging the gap between soybean translational genomics and molecular breeding. It provides information for six entities including genes/proteins, microRNAs (miRNAs)/small interfering RNAs (sRNA), metabolites, single nucleotide polymorphisms (SNPs), and plant introduction lines and traits. It has a user-friendly web interface publicly available at http://soykb.org , which integrates and presents data in an intuitive manner to the soybean researchers, breeders, and consumers. It incorporates several informatics and analytical tools for integrating and merging various multi-omics datasets.

VAMP721a and VAMP721d are important for pectin dynamics and release of bacteria in soybean nodules

In root nodules rhizobia enter host cells via infection threads. The release of bacteria to a host cell is possible from cell wall-free regions of the infection thread. We hypothesized that the VAMP721d and VAMP721e exocytotic pathway, identified before in Medicago truncatula, has a role in the local modification of cell wall during the release of rhizobia. To clarify the role of VAMP721d and VAMP721e we used Glycine max, a plant with a determinate type of nodule. The localization of the main polysaccharide compounds of primary cell walls was analysed in control vs nodules with partially silenced GmVAMP721d. The silencing of GmVAMP721d blocked the release of rhizobia. Instead of rhizobia-containing membrane compartments - symbiosomes - the infected cells contained big clusters of bacteria embedded in a matrix of methyl-esterified and de-methyl-esterified pectin. These clusters were surrounded by a membrane. We found that GmVAMP721d-positive vesicles were not transporting methyl-esterified pectin. We hypothesized that they may deliver the enzymes involved in pectin turnover. Subsequently, we found that GmVAMP721d is partly co-localized with pectate lyase. Therefore, the biological role of VAMP721d may be explained by its action in delivering pectin-modifying enzymes to the site of release.

A toolbox of genes, proteins, metabolites and promoters for improving drought tolerance in soybean includes the metabolite coumestrol and stomatal development genes

BACKGROUND

The purpose of this project was to identify metabolites, proteins, genes, and promoters associated with water stress responses in soybean. A number of these may serve as new targets for the biotechnological improvement of drought responses in soybean (Glycine max).

RESULTS

We identified metabolites, proteins, and genes that are strongly up or down regulated during rapid water stress following removal from a hydroponics system. 163 metabolites showed significant changes during water stress in roots and 93 in leaves. The largest change was a root-specific 160-fold increase in the coumestan coumestrol making it a potential biomarker for drought and a promising target for improving drought responses. Previous reports suggest that coumestrol stimulates mycorrhizal colonization and under certain conditions mycorrhizal plants have improved drought tolerance. This suggests that coumestrol may be part of a call for help to the rhizobiome during stress. About 3,000 genes were strongly up-regulated by drought and we identified regulators such as ERF, MYB, NAC, bHLH, and WRKY transcription factors, receptor-like kinases, and calcium signaling components as potential targets for soybean improvement as well as the jasmonate and abscisic acid biosynthetic genes JMT, LOX1, and ABA1. Drought stressed soybean leaves show reduced mRNA levels of stomatal development genes including FAMA-like, MUTE-like and SPEECHLESS-like bHLH transcription factors and leaves formed after drought stress had a reduction in stomatal density of 22.34 % and stomatal index of 17.56 %. This suggests that reducing stomatal density may improve drought tolerance. MEME analyses suggest that ABRE (CACGT/CG), CRT/DRE (CCGAC) and a novel GTGCnTGC/G element play roles in transcriptional activation and these could form components of synthetic promoters to drive expression of transgenes. Using transformed hairy roots, we validated the increase in promoter activity of GmWRKY17 and GmWRKY67 during dehydration and after 20 μM ABA treatment.

CONCLUSIONS

Our toolbox provides new targets and strategies for improving soybean drought tolerance and includes the coumestan coumestrol, transcription factors that regulate stomatal density, water stress-responsive WRKY gene promoters and a novel DNA element that appears to be enriched in water stress responsive promoters.

Evolution of plant δ(1)-pyrroline-5-carboxylate reductases from phylogenetic and structural perspectives

Proline plays a crucial role in cell growth and stress responses, and its accumulation is essential for the tolerance of adverse environmental conditions in plants. Two routes are used to biosynthesize proline in plants. The main route uses glutamate as a precursor, while in the other route proline is derived from ornithine. The terminal step of both pathways, the conversion of δ(1)-pyrroline-5-carboxylate (P5C) to L-proline, is catalyzed by P5C reductase (P5CR) using NADH or NADPH as a cofactor. Since P5CRs are important housekeeping enzymes, they are conserved across all domains of life and appear to be relatively unaffected throughout evolution. However, global analysis of these enzymes unveiled significant functional diversity in the preference for cofactors (NADPH vs. NADH), variation in metal dependence and the differences in the oligomeric state. In our study we investigated evolutionary patterns through phylogenetic and structural analysis of P5CR representatives from all kingdoms of life, with emphasis on the plant species. We also attempted to correlate local sequence/structure variation among the functionally and structurally characterized members of the family.

Evolution of proline biosynthesis: enzymology, bioinformatics, genetics, and transcriptional regulation

Proline is not only an essential component of proteins but it also has important roles in adaptation to osmotic and dehydration stresses, redox control, and apoptosis. Here, we review pathways of proline biosynthesis in the three domains of life. Pathway reconstruction from genome data for hundreds of eubacterial and dozens of archaeal and eukaryotic organisms revealed evolutionary conservation and variations of this pathway across different taxa. In the most prevalent pathway of proline synthesis, glutamate is phosphorylated to γ-glutamyl phosphate by γ-glutamyl kinase, reduced to γ-glutamyl semialdehyde by γ-glutamyl phosphate reductase, cyclized spontaneously to Δ(1)-pyrroline-5-carboxylate and reduced to proline by Δ(1)-pyrroline-5-carboxylate reductase. In higher plants and animals the first two steps are catalysed by a bi-functional Δ(1) -pyrroline-5-carboxylate synthase. Alternative pathways of proline formation use the initial steps of the arginine biosynthetic pathway to ornithine, which can be converted to Δ(1)-pyrroline-5-carboxylate by ornithine aminotransferase and then reduced to proline or converted directly to proline by ornithine cyclodeaminase. In some organisms, the latter pathways contribute to or could be fully responsible for the synthesis of proline. The conservation of proline biosynthetic enzymes and significance of specific residues for catalytic activity and allosteric regulation are analysed on the basis of protein structural data, multiple sequence alignments, and mutant studies, providing novel insights into proline biosynthesis in organisms. We also discuss the transcriptional control of the proline biosynthetic genes in bacteria and plants.

SFGD: a comprehensive platform for mining functional information from soybean transcriptome data and its use in identifying acyl-lipid metabolism pathways

BACKGROUND

Soybean (Glycine max L.) is one of the world's most important leguminous crops producing high-quality protein and oil. Increasing the relative oil concentration in soybean seeds is many researchers' goal, but a complete analysis platform of functional annotation for the genes involved in the soybean acyl-lipid pathway is still lacking. Following the success of soybean whole-genome sequencing, functional annotation has become a major challenge for the scientific community. Whole-genome transcriptome analysis is a powerful way to predict genes with biological functions. It is essential to build a comprehensive analysis platform for integrating soybean whole-genome sequencing data, the available transcriptome data and protein information. This platform could also be used to identify acyl-lipid metabolism pathways.

DESCRIPTION

In this study, we describe our construction of the Soybean Functional Genomics Database (SFGD) using Generic Genome Browser (Gbrowse) as the core platform. We integrated microarray expression profiling with 255 samples from 14 groups' experiments and mRNA-seq data with 30 samples from four groups' experiments, including spatial and temporal transcriptome data for different soybean development stages and environmental stresses. The SFGD includes a gene co-expression regulatory network containing 23,267 genes and 1873 miRNA-target pairs, and a group of acyl-lipid pathways containing 221 enzymes and more than 1550 genes. The SFGD also provides some key analysis tools, i.e. BLAST search, expression pattern search and cis-element significance analysis, as well as gene ontology information search and single nucleotide polymorphism display.

CONCLUSION

The SFGD is a comprehensive database integrating genome and transcriptome data, and also for soybean acyl-lipid metabolism pathways. It provides useful toolboxes for biologists to improve the accuracy and robustness of soybean functional genomics analysis, further improving understanding of gene regulatory networks for effective crop improvement. The SFGD is publically accessible at http://bioinformatics.cau.edu.cn/SFGD/, with all data available for downloading.

SoyFN: a knowledge database of soybean functional networks

Many databases for soybean genomic analysis have been built and made publicly available, but few of them contain knowledge specifically targeting the omics-level gene-gene, gene-microRNA (miRNA) and miRNA-miRNA interactions. Here, we present SoyFN, a knowledge database of soybean functional gene networks and miRNA functional networks. SoyFN provides user-friendly interfaces to retrieve, visualize, analyze and download the functional networks of soybean genes and miRNAs. In addition, it incorporates much information about KEGG pathways, gene ontology annotations and 3'-UTR sequences as well as many useful tools including SoySearch, ID mapping, Genome Browser, eFP Browser and promoter motif scan. SoyFN is a schema-free database that can be accessed as a Web service from any modern programming language using a simple Hypertext Transfer Protocol call. The Web site is implemented in Java, JavaScript, PHP, HTML and Apache, with all major browsers supported. We anticipate that this database will be useful for members of research communities both in soybean experimental science and bioinformatics. Database URL: http://nclab.hit.edu.cn/SoyFN.

Soybean knowledge base (SoyKB): a web resource for integration of soybean translational genomics and molecular breeding

Soybean Knowledge Base (http://soykb.org) is a comprehensive web resource developed for bridging soybean translational genomics and molecular breeding research. It provides information for six entities including genes/proteins, microRNAs/sRNAs, metabolites, single nucleotide polymorphisms, plant introduction lines and traits. It also incorporates many multi-omics datasets including transcriptomics, proteomics, metabolomics and molecular breeding data, such as quantitative trait loci, traits and germplasm information. Soybean Knowledge Base has a new suite of tools such as In Silico Breeding Program for soybean breeding, which includes a graphical chromosome visualizer for ease of navigation. It integrates quantitative trait loci, traits and germplasm information along with genomic variation data, such as single nucleotide polymorphisms, insertions, deletions and genome-wide association studies data, from multiple soybean cultivars and Glycine soja.

Construction and characterization of a BAC library from the Coffea arabica genotype Timor Hybrid CIFC 832/2

Most of the world's coffee production originates from Coffea arabica, an allotetraploid species with low genetic diversity and for which few genomic resources are available. Genomic libraries with large DNA fragment inserts are useful tools for the study of plant genomes, including the production of physical maps, integration studies of physical and genetic maps, genome structure analysis and gene isolation by positional cloning. Here, we report the construction and characterization of a Bacterial Artificial Chromosome (BAC) library from C. arabica Timor Hybrid CIFC 832/2, a parental genotype for several modern coffee cultivars. The BAC library consists of 56,832 clones with an average insert size of 118 kb, which represents a dihaploid genome coverage of five to sixfold. The content of organellar DNA was estimated at 1.04 and 0.5 % for chloroplast and mitochondrial DNA, respectively. The BAC library was screened for the NADPH-dependent mannose-6-phosphate reductase gene (CaM6PR) with markers positioned on four linkage groups of a partial C. arabica genetic map. A mixed approach using PCR and membrane hybridization of BAC pools allowed for the discovery of nine BAC clones with the CaM6PR gene and 53 BAC clones that were anchored to the genetic map with simple sequence repeat markers. This library will be a useful tool for future studies on comparative genomics and the identification of genes and regulatory elements controlling major traits in this economically important crop species.

Genome-wide analysis and molecular characterization of heat shock transcription factor family in Glycine max

Heat shock transcription factors (Hsfs) play an essential role on the increased tolerance against heat stress by regulating the expression of heat-responsive genes. In this study, a genome-wide analysis was performed to identify all of the soybean (Glycine max) GmHsf genes based on the latest soybean genome sequence. Chromosomal location, protein domain, motif organization, and phylogenetic relationships of 26 non-redundant GmHsf genes were analyzed compared with AtHsfs (Arabidopsis thaliana Hsfs). According to their structural features, the predicted members were divided into the previously defined classes A-C, as described for AtHsfs. Transcript levels and subcellular localization of five GmHsfs responsive to abiotic stresses were analyzed by real-time RT-PCR. These results provide a fundamental clue for understanding the complexity of the soybean GmHsf gene family and cloning the functional genes in future studies.

Expression dynamics and genome distribution of osmoprotectants in soybean: identifying important components to face abiotic stress

BACKGROUND

Despite the importance of osmoprotectants, no previous in silico evaluation of high throughput data is available for higher plants. The present approach aimed at the identification and annotation of osmoprotectant-related sequences applied to short transcripts from a soybean HT-SuperSAGE (High Throughput Super Serial Analysis of Gene Expression; 26-bp tags) database, and also its comparison with other transcriptomic and genomic data available from different sources.

METHODS

A curated set of osmoprotectants related sequences was generated using text mining and selected seed sequences for identification of the respective transcripts and proteins in higher plants. To test the efficiency of the seed sequences, these were aligned against four HT-SuperSAGE contrasting libraries generated by our group using soybean tolerant and sensible plants against water deficit, considering only differentially expressed transcripts (p ≤ 0.05). Identified transcripts from soybean and their respective tags were aligned and anchored against the soybean virtual genome.

RESULTS

The workflow applied resulted in a set including 1,996 seed sequences that allowed the identification of 36 differentially expressed genes related to the biosynthesis of osmoprotectants [Proline (P5CS: 4, P5CR: 2), Trehalose (TPS1: 9, TPPB: 1), Glycine betaine (BADH: 4) and Myo-inositol (MIPS: 7, INPS1: 8)], also mapped in silico in the soybean genome (25 loci). Another approach considered matches using Arabidopsis full length sequences as seed sequences, and allowed the identification of 124 osmoprotectant-related sequences, matching ~10.500 tags anchored in the soybean virtual chromosomes. Osmoprotectant-related genes appeared clustered in all soybean chromosomes, with higher density in some subterminal regions and synteny among some chromosome pairs.

CONCLUSIONS

Soybean presents all searched osmoprotectant categories with some important members differentially expressed among the comparisons considered (drought tolerant or sensible vs. control; tolerant vs. sensible), allowing the identification of interesting candidates for biotechnological inferences. The identified tags aligned to corresponding genes that matched 19 soybean chromosomes. Osmoprotectant-related genes are not regularly distributed in the soybean genome, but clustered in some regions near the chromosome terminals, with some redundant clusters in different chromosomes indicating their involvement in previous duplication and rearrangements events. The seed sequences, transcripts and map represent the first transversal evaluation for osmoprotectant-related genes and may be easily applied to other plants of interest.

Molecular mapping and breeding with microsatellite markers

In genetics databases for crop plant species across the world, there are thousands of mapped loci that underlie quantitative traits, oligogenic traits, and simple traits recognized by association mapping in populations. The number of loci will increase as new phenotypes are measured in more diverse genotypes and genetic maps based on saturating numbers of markers are developed. A period of locus reevaluation will decrease the number of important loci as those underlying mega-environmental effects are recognized. A second wave of reevaluation of loci will follow from developmental series analysis, especially for harvest traits like seed yield and composition. Breeding methods to properly use the accurate maps of QTL are being developed. New methods to map, fine map, and isolate the genes underlying the loci will be critical to future advances in crop biotechnology. Microsatellite markers are the most useful tool for breeders. They are codominant, abundant in all genomes, highly polymorphic so useful in many populations, and both economical and technically easy to use. The selective genotyping approaches, including genotype ranking (indexing) based on partial phenotype data combined with favorable allele data and bulked segregation event (segregant) analysis (BSA), will be increasingly important uses for microsatellites. Examples of the methods for developing and using microsatellites derived from genomic sequences are presented for monogenic, oligogenic, and polygenic traits. Examples of successful mapping, fine mapping, and gene isolation are given. When combined with high-throughput methods for genotyping and a genome sequence, the use of association mapping with microsatellite markers will provide critical advances in the analysis of crop traits.

Mining plant genome browsers as a means for efficient connection of physical, genetic and cytogenetic mapping: An example using soybean

Physical maps are important tools to uncover general chromosome structure as well as to compare different plant lineages and species, helping to elucidate genome structure, evolution and possibilities regarding synteny and colinearity. The increasing production of sequence data has opened an opportunity to link information from mapping studies to the underlying sequences. Genome browsers are invaluable platforms that provide access to these sequences, including tools for genome analysis, allowing the integration of multivariate information, and thus aiding to explain the emergence of complex genomes. The present work presents a tutorial regarding the use of genome browsers to develop targeted physical mapping, providing also a general overview and examples about the possibilities regarding the use of Fluorescent In Situ Hybridization (FISH) using bacterial artificial chromosomes (BAC), simple sequence repeats (SSR) and rDNA probes, highlighting the potential of such studies for map integration and comparative genetics. As a case study, the available genome of soybean was accessed to show how the physical and in silico distribution of such sequences may be compared at different levels. Such evaluations may also be complemented by the identification of sequences beyond the detection level of cytological methods, here using members of the aquaporin gene family as an example. The proposed approach highlights the complementation power of the combination of molecular cytogenetics and computational approaches for the anchoring of coding or repetitive sequences in plant genomes using available genome browsers, helping in the determination of sequence location, arrangement and number of repeats, and also filling gaps found in computational pseudochromosome assemblies.

Genomic analysis of phospholipase D family and characterization of GmPLDαs in soybean (Glycine max)

Phospholipase D (PLD) and its product phosphatidic acid play important roles in the regulation of plant growth, development, and stress responses. The genome database analysis has revealed PLD family in Arabidopsis, rice, poplar and grape. In this study, we report a genomic analysis of 18 putative soybean (Glycine max) PLD genes (GmPLDs), which exist in the 14 of 20 chromosomes. GmPLDs were grouped into six types, α(3), β(4), γ, δ(5), ε(2), and ζ(3), based on gene architectures, protein domains, evolutionary relationship, and sequence identity. These GmPLDs contained two HKD domains, PX/PH domains (for GmPLDζs), and C2 domain (for the other GmPLDs). The expression patterns analyzed by quantitative reverse transcription PCR demonstrated that GmPLDs were expressed differentially in various tissues. GmPLDα1, α2, and β2 were highly expressed in most tissues, whereas GmPLDδ5 was only expressed in flowers and GmPLDζ1 was predominantly expressed in flowers and early pods. The expression of GmPLDα1 and α2 was increased and that of GmPLDγ was decreased by salt stress. GmPLDα1 protein expressed in E. coli was active under the reaction conditions for both PLDα and PLDδ, hydrolyzing the common membrane phospholipids phosphatidylcholine, phosphatidylethanolamine, and phosphatidylglycerol. The genomic analysis for soybean PLD family provides valuable data for further identity and characterization of their functions.

Recombination suppression at the dominant Rhg1/Rfs2 locus underlying soybean resistance to the cyst nematode

Host resistance to "yellow dwarf" or "moonlight" disease cause by any population (Hg type) of Heterodera glycines I., the soybean cyst nematode (SCN), requires a functional allele at rhg1. The host resistance encoded appears to mimic an apoptotic response in the giant cells formed at the nematode feeding site about 24-48 h after nematode feeding commences. Little is known about how the host response to infection is mediated but a linked set of 3 genes has been identified within the rhg1 locus. This study aimed to identify the role of the genes within the locus that includes a receptor-like kinase (RLK), a laccase and an ion antiporter. Used were near isogeneic lines (NILs) that contrasted at their rhg1 alleles, gene-based markers, and a new Hg type 0 and new recombination events. A syntenic gene cluster on Lg B1 was found. The effectiveness of SNP probes from the RLK for distinguishing homolog sequence variants on LgB1 from alleles at the rhg1 locus on LgG was shown. The resistant allele of the rhg1 locus was shown to be dominant in NILs. None of the recombination events were within the cluster of the three candidate genes. Finally, rhg1 was shown to reduce the plant root development. A model for rhg1 as a dominant multi-gene resistance locus based on the developmental control was inferred.

Soybean Knowledge Base (SoyKB): a web resource for soybean translational genomics

BACKGROUND

Soybean Knowledge Base (SoyKB) is a comprehensive all-inclusive web resource for soybean translational genomics. SoyKB is designed to handle the management and integration of soybean genomics, transcriptomics, proteomics and metabolomics data along with annotation of gene function and biological pathway. It contains information on four entities, namely genes, microRNAs, metabolites and single nucleotide polymorphisms (SNPs).

METHODS

SoyKB has many useful tools such as Affymetrix probe ID search, gene family search, multiple gene/metabolite search supporting co-expression analysis, and protein 3D structure viewer as well as download and upload capacity for experimental data and annotations. It has four tiers of registration, which control different levels of access to public and private data. It allows users of certain levels to share their expertise by adding comments to the data. It has a user-friendly web interface together with genome browser and pathway viewer, which display data in an intuitive manner to the soybean researchers, producers and consumers.

CONCLUSIONS

SoyKB addresses the increasing need of the soybean research community to have a one-stop-shop functional and translational omics web resource for information retrieval and analysis in a user-friendly way. SoyKB can be publicly accessed at http://soykb.org/.

Selection of a core set of RILs from Forrest x Williams 82 to develop a framework map in soybean

Soybean BAC-based physical maps provide a useful platform for gene and QTL map-based cloning, EST mapping, marker development, genome sequencing, and comparative genomic research. Soybean physical maps for "Forrest" and "Williams 82" representing the southern and northern US soybean germplasm base, respectively, have been constructed with different fingerprinting methods. These physical maps are complementary for coverage of gaps on the 20 soybean linkage groups. More than 5,000 genetic markers have been anchored onto the Williams 82 physical map, but only a limited number of markers have been anchored to the Forrest physical map. A mapping population of Forrest × Williams 82 made up of 1,025 F(8) recombinant inbred lines (RILs) was used to construct a reference genetic map. A framework map with almost 1,000 genetic markers was constructed using a core set of these RILs. The core set of the population was evaluated with the theoretical population using equality, symmetry and representativeness tests. A high-resolution genetic map will allow integration and utilization of the physical maps to target QTL regions of interest, and to place a larger number of markers into a map in a more efficient way using a core set of RILs.

Genome-wide identification of intron fragment insertion mutations and their potential use as SCAR molecular markers in the soybean

Introns often have a high probability of mutation as a result of DNA insertions and deletions (indels). In this study, 503 introns with exon-derived insertions were identified using a comprehensive search of the soybean genome. Of the 375 pairs of PCR primer sets designed for the loci in question, 161 primer sets amplified length polymorphism among nine soybean varieties and were identified as soybean gene-intron-driven functional sequence characterized amplified region (SCAR) markers. These SCAR markers are distributed among all 20 of the soybean chromosomes, and they developed from numerous genes involved in various physiological and biochemical processes that influence important agronomic traits of the soybean. The development of these novel gene-driven functional SCAR markers was fast and cost effective, and their use will facilitate molecular-assisted breeding of the soybean.

Ultrasensitive detection of non-amplified genomic DNA by nanoparticle-enhanced surface plasmon resonance imaging

Technologies today available for the DNA detection rely on a combination of labeled probes hybridized to target sequences which are amplified by polymerase chain reaction (PCR). Direct detection methods that eliminate the requirement for both PCR and labeling steps could afford faster, cheaper and simpler devices for the analysis of small amounts of unamplified DNA. In this work we describe the results obtained in the ultrasensitive detection of non-amplified genomic DNA. We analyzed certified reference materials containing different amounts of genetically modified DNA by using a detection method which combines the nanoparticle-enhanced surface plasmon resonance imaging (SPRI) biosensing to the peptide nucleic acids (PNAs) improved selectivity and sensitivity in targeting complementary DNA sequences. The method allowed us to obtain a 41 zM sensitivity in targeting genomic DNA even in the presence of a large excess of non-complementary DNA.

Iso-lines and inbred-lines confirmed loci that underlie resistance from cultivar 'Hartwig' to three soybean cyst nematode populations

Soybean [Glycine max (L.) Merr.] cultivars varied in their resistance to different populations of the soybean cyst nematode (SCN), Heterodera glycines, called HG Types. The rhg1 locus on linkage group G was necessary for resistance to all HG types. However, the loci for resistance to H. glycines HG Type 1.3- (race 14) and HG Type 1.2.5- (race 2) of the soybean cyst nematode have varied in their reported locations. The aims were to compare the inheritance of resistance to three nematode HG Types in a population segregating for resistance to SCN and to identify the underlying quantitative trait loci (QTL). 'Hartwig', a soybean cultivar resistant to most SCN HG Types, was crossed with the susceptible cultivar 'Flyer'. A total of 92 F5-derived recombinant inbred lines (RILs; or inbred lines) and 144 molecular markers were used for map development. The rhg1 associated QTL found in earlier studies were confirmed and shown to underlie resistance to all three HG Types in RILs (Satt309; HG Type 0, P = 0.0001 R (2) = 22%; Satt275; HG Type 1.3, P = 0.001, R (2) = 14%) and near isogeneic lines (NILs; or iso-lines; Satt309; HG Type 1.2.5-, P = 0.001 R (2) = 24%). A new QTL underlying resistance to HG Type 1.2.5- was detected on LG D2 (Satt574; P = 0.001, R (2) = 11%) among 14 RILs resistant to the other HG types. The locus was confirmed in a small NIL population consisting of 60 plants of ten genotypes (P = 0.04). This QTL (cqSCN-005) is located in an interval previously associated with resistance to both SDS leaf scorch from 'Pyramid' and 'Ripley' (cqSDS-001) and SCN HG Type 1.3- from Hartwig and Pyramid. The QTL detected will allow marker assisted selection for multigenic resistance to complex nematode populations in combination with sudden death syndrome resistance (SDS) and other agronomic traits.

SoyDB: a knowledge database of soybean transcription factors

BACKGROUND

Transcription factors play the crucial rule of regulating gene expression and influence almost all biological processes. Systematically identifying and annotating transcription factors can greatly aid further understanding their functions and mechanisms. In this article, we present SoyDB, a user friendly database containing comprehensive knowledge of soybean transcription factors.

DESCRIPTION

The soybean genome was recently sequenced by the Department of Energy-Joint Genome Institute (DOE-JGI) and is publicly available. Mining of this sequence identified 5,671 soybean genes as putative transcription factors. These genes were comprehensively annotated as an aid to the soybean research community. We developed SoyDB - a knowledge database for all the transcription factors in the soybean genome. The database contains protein sequences, predicted tertiary structures, putative DNA binding sites, domains, homologous templates in the Protein Data Bank (PDB), protein family classifications, multiple sequence alignments, consensus protein sequence motifs, web logo of each family, and web links to the soybean transcription factor database PlantTFDB, known EST sequences, and other general protein databases including Swiss-Prot, Gene Ontology, KEGG, EMBL, TAIR, InterPro, SMART, PROSITE, NCBI, and Pfam. The database can be accessed via an interactive and convenient web server, which supports full-text search, PSI-BLAST sequence search, database browsing by protein family, and automatic classification of a new protein sequence into one of 64 annotated transcription factor families by hidden Markov models.

CONCLUSIONS

A comprehensive soybean transcription factor database was constructed and made publicly accessible at http://casp.rnet.missouri.edu/soydb/.

Gene duplication, transfer, and evolution in the chloroplast genome

In addition to the nuclear genome, organisms have organelle genomes. Most of the DNA present in eukaryotic organisms is located in the cell nucleus. Chloroplasts have independent genomes which are inherited from the mother. Duplicated genes are common in the genomes of all organisms. It is believed that gene duplication is the most important step for the origin of genetic variation, leading to the creation of new genes and new gene functions. Despite the fact that extensive gene duplications are rare among the chloroplast genome, gene duplication in the chloroplast genome is an essential source of new genetic functions and a mechanism of neo-evolution. The events of gene transfer between the chloroplast genome and nuclear genome via duplication and subsequent recombination are important processes in evolution. The duplicated gene or genome in the nucleus has been the subject of several recent reviews. In this review, we will briefly summarize gene duplication and evolution in the chloroplast genome. Also, we will provide an overview of gene transfer events between chloroplast and nuclear genomes.

Map-based cloning of the gene associated with the soybean maturity locus E3

Photosensitivity plays an essential role in the response of plants to their changing environments throughout their life cycle. In soybean [Glycine max (L.) Merrill], several associations between photosensitivity and maturity loci are known, but only limited information at the molecular level is available. The FT3 locus is one of the quantitative trait loci (QTL) for flowering time that corresponds to the maturity locus E3. To identify the gene responsible for this QTL, a map-based cloning strategy was undertaken. One phytochrome A gene (GmPhyA3) was considered a strong candidate for the FT3 locus. Allelism tests and gene sequence comparisons showed that alleles of Misuzudaizu (FT3/FT3; JP28856) and Harosoy (E3/E3; PI548573) were identical. The GmPhyA3 alleles of Moshidou Gong 503 (ft3/ft3; JP27603) and L62-667 (e3/e3; PI547716) showed weak or complete loss of function, respectively. High red/far-red (R/FR) long-day conditions enhanced the effects of the E3/FT3 alleles in various genetic backgrounds. Moreover, a mutant line harboring the nonfunctional GmPhyA3 flowered earlier than the original Bay (E3/E3; PI553043) under similar conditions. These results suggest that the variation in phytochrome A may contribute to the complex systems of soybean flowering response and geographic adaptation.

Construction and characterization of twoCitrusBAC libraries and identification of clones containing the phytoene synthase gene

Two deep-coverage Bacterial Artificial Chromosome (BAC) libraries of Citrus sinensis (L.) Osbeck ‘Cara Cara’ navel orange and Citrus reticulata (L.) Blanco ‘Egan No. 1’ Ponkan mandarin, which belong to the two most important species of the Citrus genus, have been constructed and characterized to facilitate gene cloning and to analyze variety-specific genome composition. The C. sinensis BAC library consists of 36 000 clones with negligible false-positive clones and an estimated average insert size of 126 kb covering ~4.5 × 109 bp and thus providing an 11.8-fold coverage of haploid genome equivalents, whereas the C. reticulata library consists of 21 000 clones also with negligible false-positive clones and an estimated average of 120 kb covering ~2.5 × 109 bp representing a 6.6-fold coverage of haploid genome equivalents. Both libraries were evaluated for contamination with high-copy vector, empty pIndigoBAC536 vector, and organellar DNA sequences. Screening has been performed by Southern hybridization of BAC filters, which results in <0.5% chloroplast DNA contamination and no mitochondrial DNA contamination in both libraries. Eight and five positive clones harboring the gene encoding Phytoene synthase (Psy (EC 2.5.1.32)) were identified from the C. sinensis and C. reticulata libraries, respectively, using the filter hybridization procedure. These results suggest that the two BAC libraries are useful tools for the isolation of functional genes and advanced genomics research in the two important species C. sinensis and C. reticulata. Resources, high-density filters, individual clones, and whole libraries are available for public distribution and are accessible at the National Key Laboratory of Crop Genetic Improvement, Huazhong Agricultural University.

An overview of the apple genome through BAC end sequence analysis

The apple, Malus x domestica Borkh., is one of the most important fruit trees grown worldwide. A bacterial artificial chromosome (BAC)-based physical map of the apple genome has been recently constructed. Based on this physical map, a total of approximately 2,100 clones from different contigs (overlapping BAC clones) have been selected and sequenced at both ends, generating 3,744 high-quality BAC end sequences (BESs) including 1,717 BAC end pairs. Approximately 8.5% of BESs contain simple sequence repeats (SSRs), most of which are AT/TA dimer repeats. Potential transposable elements are identified in approximately 21% of BESs, and most of these elements are retrotransposons. About 11% of BESs have homology to the Arabidopsis protein database. The matched proteins cover a broad range of categories. The average GC content of the predicted coding regions of BESs is 42.4%; while, that of the whole BESs is 39%. A small number of BES pairs were mapped to neighboring chromosome regions of A. thaliana and Populus trichocarpa; whereas, no pairs are mapped to the Oryza sativa genome. The apple has a higher degree of synteny with the closely related Populus than with the distantly related Arabidopsis. BAC end sequencing can be used to anchor a small proportion of the apple genome to the Populus and possibly to the Arabidopsis genomes.

Re-annotation of the physical map of Glycine max for polyploid-like regions by BAC end sequence driven whole genome shotgun read assembly

BACKGROUND

Many of the world's most important food crops have either polyploid genomes or homeologous regions derived from segmental shuffling following polyploid formation. The soybean (Glycine max) genome has been shown to be composed of approximately four thousand short interspersed homeologous regions with 1, 2 or 4 copies per haploid genome by RFLP analysis, microsatellite anchors to BACs and by contigs formed from BAC fingerprints. Despite these similar regions,, the genome has been sequenced by whole genome shotgun sequence (WGS). Here the aim was to use BAC end sequences (BES) derived from three minimum tile paths (MTP) to examine the extent and homogeneity of polyploid-like regions within contigs and the extent of correlation between the polyploid-like regions inferred from fingerprinting and the polyploid-like sequences inferred from WGS matches.

RESULTS

Results show that when sequence divergence was 1-10%, the copy number of homeologous regions could be identified from sequence variation in WGS reads overlapping BES. Homeolog sequence variants (HSVs) were single nucleotide polymorphisms (SNPs; 89%) and single nucleotide indels (SNIs 10%). Larger indels were rare but present (1%). Simulations that had predicted fingerprints of homeologous regions could be separated when divergence exceeded 2% were shown to be false. We show that a 5-10% sequence divergence is necessary to separate homeologs by fingerprinting. BES compared to WGS traces showed polyploid-like regions with less than 1% sequence divergence exist at 2.3% of the locations assayed.

CONCLUSION

The use of HSVs like SNPs and SNIs to characterize BACs wil improve contig building methods. The implications for bioinformatic and functional annotation of polyploid and paleopolyploid genomes show that a combined approach of BAC fingerprint based physical maps, WGS sequence and HSV-based partitioning of BAC clones from homeologous regions to separate contigs will allow reliable de-convolution and positioning of sequence scaffolds (see BES_scaffolds section of SoyGD). This approach will assist genome annotation for paleopolyploid and true polyploid genomes such as soybean and many important cereal and fruit crops.

Mapping two genes in the purine metabolism pathway of soybean

Mapping genes in biochemical pathways allow study of the genomic organization of pathways and geneic relationships within these pathways. Additionally, molecular markers located within the boundaries of a specific gene sequence represent important marker assisted selection resources. We report map locations of two geneic markers from the purine synthesis pathway in soybean (Glycine max (L. merr.)), utilizing a 90 plant F(2) population created from the cross of "DT97-4290" x "DS97-84-1". Primers were designed based on sequences from annotated soybean complimentary DNA. A polymorphic, co-dominant, sequence-characterized amplified region marker was created for hypoxanthine phosphoribosyl transferase (EC 2.4.2.8). Linkage analysis placed this gene on linkage group (LG) O. In addition, a single-nucleotide polymorphism (SNP) marker was developed for a urate oxidase gene (EC 1.7.3.3). Linkage analysis of the SNP placed the urate oxidase gene on LG I. For both genes, amplicon sequence data confirmed the identification of the respective gene. Mapping these genes represents the first step in understanding the genomic organization of the purine biochemical pathway in soybean.

Molecular mapping of two loci that confer resistance to Asian rust in soybean

Asian soybean rust (ASR) is caused by the fungal pathogen Phakopsora pachyrhizi Sydow & Sydow. It was first identified in Brazil in 2001 and quickly infected soybean areas in several countries in South America. Primary efforts to combat this disease must involve the development of resistant cultivars. Four distinct genes that confer resistance against ASR have been reported: Rpp1, Rpp2, Rpp3, and Rpp4. However, no cultivar carrying any of those resistance loci has been released. The main objective of this study was to genetically map Rpp2 and Rpp4 resistance genes. Two F(2:3) populations, derived from the crosses between the resistant lines PI 230970 (Rpp2), PI 459025 (Rpp4) and the susceptible cultivar BRS 184, were used in this study. The mapping populations and parental lines were inoculated with a field isolate of P. pachyrhizi and evaluated for lesion type as resistant (RB lesions) or susceptible (TAN lesions). The mapping populations were screened with SSR markers, using the bulk segregant analysis (BSA) to expedite the identification of linked markers. Both resistance genes showed an expected segregation ratio for a dominant trait. This study allowed mapping Rpp2 and Rpp4 loci on the linkage groups J and G, respectively. The associated markers will be of great value on marker assisted selection for this trait.

Separate loci underlie resistance to root infection and leaf scorch during soybean sudden death syndrome

Soybean [Glycine max (L.) Merr.] cultivars show differences in their resistance to both the leaf scorch and root rot of sudden death syndrome (SDS). The syndrome is caused by root colonization by Fusarium virguliforme (ex. F. solani f. sp. glycines). Root susceptibility combined with reduced leaf scorch resistance has been associated with resistance to Heterodera glycines HG Type 1.3.6.7 (race 14) of the soybean cyst nematode (SCN). In contrast, the rhg1 locus underlying resistance to Hg Type 0 was found clustered with three loci for resistance to SDS leaf scorch and one for root infection. The aims of this study were to compare the inheritance of resistance to leaf scorch and root infection in a population that segregated for resistance to SCN and to identify the underlying quantitative trait loci (QTL). "Hartwig", a cultivar partially resistant to SDS leaf scorch, F. virguliforme root infection and SCN HG Type 1.3.6.7 was crossed with the partially susceptible cultivar "Flyer". Ninety-two F5-derived recombinant inbred lines and 144 markers were used for map development. Four QTL found in earlier studies were confirmed. One contributed resistance to leaf scorch on linkage group (LG) C2 (Satt277; P = 0.004, R2 = 15%). Two on LG G underlay root infection at R8 (Satt038; P = 0.0001 R2 = 28.1%; Satt115; P = 0.003, R2 = 12.9%). The marker Satt038 was linked to rhg1 underlying resistance to SCN Hg Type 0. The fourth QTL was on LG D2 underlying resistance to root infection at R6 (Satt574; P = 0.001, R2 = 10%). That QTL was in an interval previously associated with resistance to both SDS leaf scorch and SCN Hg Type 1.3.6.7. The QTL showed repulsion linkage with resistance to SCN that may explain the relative susceptibility to SDS of some SCN resistant cultivars. One additional QTL was discovered on LG G underlying resistance to SDS leaf scorch measured by disease index (Satt130; P = 0.003, R2 = 13%). The loci and markers will provide tagged alleles with which to improve the breeding of cultivars combining resistances to SDS leaf scorch, root infection and SCN HG Type 1.3.6.7.

Microsatellite discovery from BAC end sequences and genetic mapping to anchor the soybean physical and genetic maps

Whole-genome sequencing of the soybean (Glycine max (L.) Merr. 'Williams 82') has made it important to integrate its physical and genetic maps. To facilitate this integration of maps, we screened 3290 microsatellites (SSRs) identified from BAC end sequences of clones comprising the 'Williams 82' physical map. SSRs were screened against 3 mapping populations. We found the AAT and ACT motifs produced the greatest frequency of length polymorphisms, ranging from 17.2% to 32.3% and from 11.8% to 33.3%, respectively. Other useful motifs include the dinucleotide repeats AG, AT, and AG, with frequency of length polymorphisms ranging from 11.2% to 18.4% (AT), 12.4% to 20.6% (AG), and 11.3% to 16.4% (GT). Repeat lengths less than 16 bp were generally less useful than repeat lengths of 40-60 bp. Two hundred and sixty-five SSRs were genetically mapped in at least one population. Of the 265 mapped SSRs, 60 came from BAC singletons not yet placed into contigs of the physical map. One hundred and ten originated in BACs located in contigs for which no genetic map location was previously known. Ninety-five SSRs came from BACs within contigs for which one or more other BACs had already been mapped. For these fingerprinted contigs (FPC) a high percentage of the mapped markers showed inconsistent map locations. A strategy is introduced by which physical and genetic map inconsistencies can be resolved using the preliminary 4x assembly of the whole genome sequence of soybean.

The soybean-Phytophthora resistance locus Rps1-k encompasses coiled coil-nucleotide binding-leucine rich repeat-like genes and repetitive sequences

BACKGROUND

A series of Rps (resistance to Pytophthora sojae) genes have been protecting soybean from the root and stem rot disease caused by the Oomycete pathogen, Phytophthora sojae. Five Rps genes were mapped to the Rps1 locus located near the 28 cM map position on molecular linkage group N of the composite genetic soybean map. Among these five genes, Rps1-k was introgressed from the cultivar, Kingwa. Rps1-k has been providing stable and broad-spectrum Phytophthora resistance in the major soybean-producing regions of the United States. Rps1-k has been mapped and isolated. More than one functional Rps1-k gene was identified from the Rps1-k locus. The clustering feature at the Rps1-k locus might have facilitated the expansion of Rps1-k gene numbers and the generation of new recognition specificities. The Rps1-k region was sequenced to understand the possible evolutionary steps that shaped the generation of Phytophthora resistance genes in soybean.

RESULTS

Here the analyses of sequences of three overlapping BAC clones containing the 184,111 bp Rps1-k region are reported. A shotgun sequencing strategy was applied in sequencing the BAC contig. Sequence analysis predicted a few full-length genes including two Rps1-k genes, Rps1-k-1 and Rps1-k-2. Previously reported Rps1-k-3 from this genomic region 1 was evolved through intramolecular recombination between Rps1-k-1 and Rps1-k-2 in Escherichia coli. The majority of the predicted genes are truncated and therefore most likely they are nonfunctional. A member of a highly abundant retroelement, SIRE1, was identified from the Rps1-k region. The Rps1-k region is primarily composed of repetitive sequences. Sixteen simple repeat and 63 tandem repeat sequences were identified from the locus.

CONCLUSION

These data indicate that the Rps1 locus is located in a gene-poor region. The abundance of repetitive sequences in the Rps1-k region suggested that the location of this locus is in or near a heterochromatic region. Poor recombination frequencies combined with presence of two functional Rps genes at this locus has been providing stable Phytophthora resistance in soybean.

An integrated high-density linkage map of soybean with RFLP, SSR, STS, and AFLP markers using A single F2 population

Soybean [Glycine max (L.) Merrill] is the most important leguminous crop in the world due to its high contents of high-quality protein and oil for human and animal consumption as well as for industrial uses. An accurate and saturated genetic linkage map of soybean is an essential tool for studies on modern soybean genomics. In order to update the linkage map of a F2 population derived from a cross between Misuzudaizu and Moshidou Gong 503 and to make it more informative and useful to the soybean genome research community, a total of 318 AFLP, 121 SSR, 108 RFLP, and 126 STS markers were newly developed and integrated into the framework of the previously described linkage map. The updated genetic map is composed of 509 RFLP, 318 SSR, 318 AFLP, 97 AFLP-derived STS, 29 BAC-end or EST-derived STS, 1 RAPD, and five morphological markers, covering a map distance of 3080 cM (Kosambi function) in 20 linkage groups (LGs). To our knowledge, this is presently the densest linkage map developed from a single F2 population in soybean. The average intermarker distance was reduced to 2.41 from 5.78 cM in the earlier version of the linkage map. Most SSR and RFLP markers were relatively evenly distributed among different LGs in contrast to the moderately clustered AFLP markers. The number of gaps of more than 25 cM was reduced to 6 from 19 in the earlier version of the linkage map. The coverage of the linkage map was extended since 17 markers were mapped beyond the distal ends of the previous linkage map. In particular, 17 markers were tagged in a 5.7 cM interval between CE47M5a and Satt100 on LG C2, where several important QTLs were clustered. This newly updated soybean linkage map will enable to streamline positional cloning of agronomically important trait locus genes, and promote the development of physical maps, genome sequencing, and other genomic research activities.

End-sequencing and characterization of silkworm (Bombyx mori) bacterial artificial chromosome libraries

Background

We performed large-scale bacterial artificial chromosome (BAC) end-sequencing of two BAC libraries (an EcoRI- and a BamHI-digested library) and conducted an in silico analysis to characterize the obtained sequence data, to make them a useful resource for genomic research on the silkworm (Bombyx mori).

Results

More than 94000 BAC end sequences (BESs), comprising more than 55 Mbp and covering about 10.4% of the silkworm genome, were sequenced. Repeat-sequence analysis with known repeat sequences indicated that the long interspersed nuclear elements (LINEs) were abundant in BamHI BESs, whereas DNA-type elements were abundant in EcoRI BESs. Repeat-sequence analysis revealed that the abundance of LINEs might be due to a GC bias of the restriction sites and that the GC content of silkworm LINEs was higher than that of mammalian LINEs. In a BLAST-based sequence analysis of the BESs against two available whole-genome shotgun sequence data sets, more than 70% of the BESs had a BLAST hit with an identity of ≥ 99%. About 14% of EcoRI BESs and about 8% of BamHI BESs were paired-end clones with unique sequences at both ends. Cluster analysis of the BESs clarified the proportion of BESs containing protein-coding regions.

Conclusion

As a result of this characterization, the identified BESs will be a valuable resource for genomic research on Bombyx mori, for example, as a base for construction of a BAC-based physical map. The use of multiple complementary BAC libraries constructed with different restriction enzymes also makes the BESs a more valuable genomic resource. The GenBank accession numbers of the obtained end sequences are DE283657–DE378560.

The development of BAC-end sequence-based microsatellite markers and placement in the physical and genetic maps of soybean

The composite map of soybean shared among Soybase, LIS and SoyGD (March 2006) contained 3,073 DNA markers in the "Locus" class. Among the markers were 1,019 class I microsatellite markers with 2-3 bp simple sequence repeats (SSRs) of >10 iterations (BARC-SSR markers). However, there were few class II SSRs (2-5 bp repeats with <10 iterations; mostly SIUC-Satt markers). The aims here were to increase the number of classes I and II SSR markers and to integrate bacterial artificial chromosome (BAC) clones onto the soybean physical map using the markers. Used was 10 Mb of BAC-end sequence (BES) derived from 13,473 reads from 7,050 clones constituting minimum tile path 2 of the soybean physical map ( http://www.soybeangenome.siu.edu ; SoyGD). Identified were 1,053 1-6 bp motif, repeat sequences, 333 from class I (>10 repeats) and 720 from class II (<10 repeats). Potential markers were shown on the MTP_SSR track at Gbrowse. Primers were designed as 20-24 bp oligomers that had Tm of 55 +/- 1 C that would generate 100-500 bp amplicons. About 853 useful primer pairs were established. Motifs were not randomly distributed with biases toward AT rich motifs. Strong biases against the GC motif and all tetra-nucleotide repeats were found. The markers discovered were useful. Among the first 135 targeted for use in genetic map improvement about 60% of class II markers and 75% of class I markers were polymorphic among on the parents of four recombinant inbred line (RIL) populations. Many of the BES-based SSRs were located on the soybean genetic map in regions with few BARC-SSR markers. Therefore, BES-based SSRs represent useful tools for genetic map development in soybean. New members of a consortium to map the markers in additional populations are invited.

A sequence based synteny map between soybean and Arabidopsis thaliana

Background

Soybean (Glycine max, L. Merr.) is one of the world's most important crops, however, its complete genomic sequence has yet to be determined. Nonetheless, a large body of sequence information exists, particularly in the form of expressed sequence tags (ESTs). Herein, we report the use of the model organism Arabidopsis thaliana (thale cress) for which the entire genomic sequence is available as a framework to align thousands of short soybean sequences.

Results

A series of JAVA-based programs were created that processed and compared 341,619 soybean DNA sequences against A. thaliana chromosomal DNA. A. thaliana DNA was probed for short, exact matches (15 bp) to each soybean sequence, and then checked for the number of additional 7 bp matches in the adjacent 400 bp region. The position of these matches was used to order soybean sequences in relation to the A. thaliana genome.

Conclusion

Reported associations between soybean sequences and A. thaliana were within a 95% confidence interval of e^-30 – e^-100. In addition, the clustering of soybean expressed sequence tags (ESTs) based on A. thaliana sequence was accurate enough to identify potential single nucleotide polymorphisms (SNPs) within the soybean sequence clusters. An EST, bacterial artificial chromosome (BAC) end sequence and marker amplicon sequence synteny map of soybean and A. thaliana is presented. In addition, all JAVA programs used to create this map are available upon request and on the WEB.

Development of a pooled probe method for locating small gene families in a physical map of soybean using stress related paralogues and a BAC minimum tile path

BACKGROUND

Genome analysis of soybean (Glycine max L.) has been complicated by its paleo-autopolyploid nature and conserved homeologous regions. Landmarks of expressed sequence tags (ESTs) located within a minimum tile path (MTP) of contiguous (contig) bacterial artificial chromosome (BAC) clones or radiation hybrid set can identify stress and defense related gene rich regions in the genome. A physical map of about 2,800 contigs and MTPs of 8,064 BAC clones encompass the soybean genome. That genome is being sequenced by whole genome shotgun methods so that reliable estimates of gene family size and gene locations will provide a useful tool for finishing. The aims here were to develop methods to anchor plant defense- and stress-related gene paralogues on the MTP derived from the soybean physical map, to identify gene rich regions and to correlate those with QTL for disease resistance.

RESULTS

The probes included 143 ESTs from a root library selected by subtractive hybridization from a multiply disease resistant soybean cultivar 'Forrest' 14 days after inoculation with Fusarium solani f. sp. glycines (F. virguliforme). Another 166 probes were chosen from a root EST library (Gm-r1021) prepared from a non-inoculated soybean cultivar 'Williams 82' based on their homology to the known defense and stress related genes. Twelve and thirteen pooled EST probes were hybridized to high-density colony arrays of MTP BAC clones from the cv. 'Forrest' genome. The EST pools located 613 paralogues for 201 of the 309 probes used (range 1-13 per functional probe). One hundred BAC clones contained more than one kind of paralogue. Many more BACs (246) contained a single paralogue of one of the 201 probes detectable gene families. ESTs were anchored on soybean linkage groups A1, B1, C2, E, D1a+Q, G, I, M, H, and O.

CONCLUSION

Estimates of gene family sizes were more similar to those made by Southern hybridization than by bioinformatics inferences from EST collections. When compared to Arabidopsis thaliana there were more 2 and 4 member paralogue families reflecting the diploidized-tetraploid nature of the soybean genome. However there were fewer families with 5 or more genes and the same number of single genes. Therefore the method can identify evolutionary patterns such as massively extensive selective gene loss or rapid divergence to regenerate the unique genes in some families.

Genomic analysis of the rhg1 locus: candidate genes that underlie soybean resistance to the cyst nematode

The rhg1 gene or genes lie at a recessive or co-dominant locus, necessary for resistance to all Hg types of the soybean (Glycine max (L.) Merr.) cyst nematode (Heterodera glycines I.). The aim here was to identify nucleotide changes within a candidate gene found at the rhg1 locus that were capable of altering resistance to Hg types 0 (race 3). A 1.5 +/- 0.25 cM region of chromosome 18 (linkage group G) was shown to encompass rhg1 using recombination events from four near isogenic line populations and nine DNA markers. The DNA markers anchored two bacterial artificial chromosome (BAC) clones 21d9 and 73p6. A single receptor like kinase (RLK; leucine rich repeat-transmembrane-protein kinase) candidate resistance gene was amplified from both BACs using redundant primers. The DNA sequence showed nine alleles of the RLK at Rhg1 in the soybean germplasm. Markers designed to detect alleles showed perfect association between allele 1 and resistance to soybean cyst nematode Hg types 0 in three segregating populations, fifteen additional selected recombination events and twenty-two Plant Introductions. A quantitative trait nucleotide (QTN) [corrected] in the RLK at rhg1 was inferred that alters A87 to V87 in the context of H274 rather than N274. [corrected] Contiguous DNA sequence of 315 kbp of chromosome 18 (about 2 cM) contained additional gene candidates that may modulate resistance to other Hg-types including a variant laccase, a hydrogen-sodium ion antiport and two proteins of unknown function. A molecular basis for recessive and co-dominant resistance that involves interactions among paralagous disease-resistance genes was inferred that would improve methods for developing new nematode-resistant soybean cultivars.