We report a large-scale analysis of the patterns of genome-wide genetic variation in soybeans. We re-sequenced a total of 17 wild and 14 cultivated soybean genomes to an average of approximately ×5 depth and >90% coverage using the Illumina Genome Analyzer II platform. We compared the patterns of genetic variation between wild and cultivated soybeans and identified higher allelic diversity in wild soybeans. We identified a high level of linkage disequilibrium in the soybean genome, suggesting that marker-assisted breeding of soybean will be less challenging than map-based cloning. We report linkage disequilibrium block location and distribution, and we identified a set of 205,614 tag SNPs that may be useful for QTL mapping and association studies. The data here provide a valuable resource for the analysis of wild soybeans and to facilitate future breeding and quantitative trait analysis.

Using a whole-genome-sequencing approach to explore germplasm resources can serve as an important strategy for crop improvement, especially in investigating wild accessions that may contain useful genetic resources that have been lost during the domestication process. Here we sequence and assemble a draft genome of wild soybean and construct a recombinant inbred population for genotyping-by-sequencing and phenotypic analyses to identify multiple QTLs relevant to traits of interest in agriculture. We use a combination of de novo sequencing data from this work and our previous germplasm re-sequencing data to identify a novel ion transporter gene, GmCHX1, and relate its sequence alterations to salt tolerance. Rapid gain-of-function tests show the protective effects of GmCHX1 towards salt stress. This combination of whole-genome de novo sequencing, high-density-marker QTL mapping by re-sequencing and functional analyses can serve as an effective strategy to unveil novel genomic information in wild soybean to facilitate crop improvement.

Genes encoding ion transporters that regulate ion homeostasis in soybean have not been carefully investigated. Using degenerate primers, we cloned a putative chloride channel gene ( GmCLC1 ) and a putative Na + /H + antiporter gene ( GmNHX1 ) from soybean. Confocal microscopic studies using yellow fluorescent fusion proteins revealed that GmCLC1 and GmNHX1 were both localized on tonoplast. The expressions of GmCLC1 and GmNHX1 were both induced by NaCl or dehydration stress imposed by polyethylene glycol (PEG). Using mitochondrial integrity and cell death as the damage indicators, a clear alleviation under NaCl stress (but not PEG stress) was observed in both GmCLC1 and GmNHX1 transgenic cells. Using fluorescent dye staining and quenching, respectively, a higher concentration of chloride ion (Cl -) or sodium ion (Na + ) was observed in isolated vacuoles in the cells of GmCLC1 and of GmNHX1 transgenic lines. Our result suggested that these vacuolar-located ion transporters function to sequester ions from cytoplasm into vacuole to reduce its toxic effects.

To obtain a comprehensive understanding of transcriptomic reprogramming under salt stress, we performed whole-transcriptome sequencing on the leaf and root of soybean seedlings subjected to salt treatment in a time-course experiment (0, 1, 2, 4, 24, and 48 hr). This time series dataset enabled us to identify important hubs and connections of gene expressions. We highlighted the analysis on phytohormone signaling pathways and their possible crosstalks. Differential expressions were also found among those genes involved in carbon and nitrogen metabolism. In general, the salt-treated seedlings slowed down their photosynthetic functions and ramped up sugar catabolism to provide extra energy for survival. Primary nitrogen assimilation was shut down whereas nitrogen resources were redistributed. Overall, the results from the transcriptomic analyses indicate that the plant uses a multipronged approach to overcome salt stress, with both fast-acting, immediate physiological responses, and longer term reactions that may involve metabolic adjustment.

No abstract

The BURP-domain protein family comprises a diverse group of plant-specific proteins that share a conserved BURP domain at the C terminus. However, there have been only limited studies on the functions and subcellular localization of these proteins. Members of the RD22-like subfamily are postulated to associate with stress responses due to the stress-inducible nature of some RD22-like genes. In this report, we used different transgenic systems (cells and in planta) to show that the expression of a stress-inducible RD22-like protein from soybean (GmRD22) can alleviate salinity and osmotic stress. We also performed detailed microscopic studies using both fusion proteins and immuno-electron microscopic techniques to demonstrate the apoplast localization of GmRD22, for which the BURP domain is a critical determinant of the subcellular localization. The apoplastic GmRD22 interacts with a cell wall peroxidase and the ectopic expression of GmRD22 in both transgenic Arabidopsis thaliana and transgenic rice resulted in increased lignin production when subjected to salinity stress. It is possible that GmRD22 regulates cell wall peroxidases and hence strengthens cell wall integrity under such stress conditions.

BackgroundLegumes are the second-most important crop family in agriculture for its economic and nutritional values. Disease resistance (R-) genes play an important role in responding to pathogen infections in plants. To further increase the yield of legume crops, we need a comprehensive understanding of the evolution of R-genes in the legume family.ResultsIn this study, we developed a robust pipeline and identified a total of 4,217 R-genes in the genomes of seven sequenced legume species. A dramatic diversity of R-genes with structural variances indicated a rapid birth-and-death rate during the R-gene evolution in legumes. The number of R-genes transiently expanded and then quickly contracted after whole-genome duplications, which meant that R-genes were sensitive to subsequent diploidization. R proteins with the Coiled-coil (CC) domain are more conserved than others in legumes. Meanwhile, other types of legume R proteins with only one or two typical domains were subjected to higher rates of loss during evolution. Although R-genes evolved quickly in legumes, they tended to undergo purifying selection instead of positive selection during evolution. In addition, domestication events in some legume species preferentially selected for the genes directly involved in the plant-pathogen interaction pathway while suppressing those R-genes with low occurrence rates.ConclusionsOur results provide insights into the dynamic evolution of R-genes in the legume family, which will be valuable for facilitating genetic improvements in the disease resistance of legume cultivars.Electronic supplementary materialThe online version of this article (doi:10.1186/s12864-016-2736-9) contains supplementary material, which is available to authorized users.

Long non-coding RNAs (lncRNAs) are defined as non-protein-coding transcripts that are at least 200 nucleotides long. They are known to play pivotal roles in regulating gene expression, especially during stress responses in plants. We used a large collection of in-house transcriptome data from various soybean (Glycine max and Glycine soja) tissues treated under different conditions to perform a comprehensive identification of soybean lncRNAs. We also retrieved publicly available soybean transcriptome data that were of sufficient quality and sequencing depth to enrich our analysis. In total, RNA-sequencing data of 332 samples were used for this analysis. An integrated reference-based, de novo transcript assembly was developed that identified ;69,000 lncRNA gene loci. We showed that lncRNAs are distinct from both protein-coding transcripts and genomic background noise in terms of length, number of exons, transposable element composition, and sequence conservation level across legume species. The tissue-specific and time-specific transcriptional responses of the lncRNA genes under some stress conditions may suggest their biological relevance. The transcription start sites of lncRNA gene loci tend to be close to their nearest proteincoding genes, and they may be transcriptionally related to the protein-coding genes, particularly for antisense and intronic lncRNAs. A previously unreported subset of small peptide-coding transcripts was identified from these lncRNA loci via tandem mass spectrometry, which paved the way for investigating their functional roles. Our results also highlight the present inadequacy of the bioinformatic definition of lncRNA, which excludes those lncRNA gene loci with small open reading frames from being regarded as protein-coding.