Silvas J. Prince

An Atlas of Soybean Small RNAs Identifies Phased siRNAs from Hundreds of Coding Genes

An extensive analysis of small RNAs in soybean identified many miRNAs and phased, secondary siRNA (phasiRNA) loci; some of these miRNAs were the triggers of the phasiRNA loci. Small RNAs are ubiquitous, versatile repressors and include (1) microRNAs (miRNAs), processed from mRNA forming stem-loops; and (2) small interfering RNAs (siRNAs), the latter derived in plants by a process typically requiring an RNA-dependent RNA polymerase. We constructed and analyzed an expression atlas of soybean (Glycine max) small RNAs, identifying over 500 loci generating 21-nucleotide phased siRNAs (phasiRNAs; from PHAS loci), of which 483 overlapped annotated protein-coding genes. Via the integration of miRNAs with parallel analysis of RNA end (PARE) data, 20 miRNA triggers of 127 PHAS loci were detected. The primary class of PHAS loci (208 or 41% of the total) corresponded to NB-LRR genes; some of these small RNAs preferentially accumulate in nodules. Among the PHAS loci, novel representatives of TAS3 and noncanonical phasing patterns were also observed. A noncoding PHAS locus, triggered by miR4392, accumulated preferentially in anthers; the phasiRNAs are predicted to target transposable elements, with their peak abundance during soybean reproductive development. Thus, phasiRNAs show tremendous diversity in dicots. We identified novel miRNAs and assessed the veracity of soybean miRNAs registered in miRBase, substantially improving the soybean miRNA annotation, facilitating an improvement of miRBase annotations and identifying at high stringency novel miRNAs and their targets.

Integrating omic approaches for abiotic stress tolerance in soybean

Soybean production is greatly influenced by abiotic stresses imposed by environmental factors such as drought, water submergence, salt, and heavy metals. A thorough understanding of plant response to abiotic stress at the molecular level is a prerequisite for its effective management. The molecular mechanism of stress tolerance is complex and requires information at the omic level to understand it effectively. In this regard, enormous progress has been made in the omics field in the areas of genomics, transcriptomics, and proteomics. The emerging field of ionomics is also being employed for investigating abiotic stress tolerance in soybean. Omic approaches generate a huge amount of data, and adequate advancements in computational tools have been achieved for effective analysis. However, the integration of omic-scale information to address complex genetics and physiological questions is still a challenge. In this review, we have described advances in omic tools in the view of conventional and modern approaches being used to dissect abiotic stress tolerance in soybean. Emphasis was given to approaches such as quantitative trait loci (QTL) mapping, genome-wide association studies (GWAS), and genomic selection (GS). Comparative genomics and candidate gene approaches are also discussed considering identification of potential genomic loci, genes, and biochemical pathways involved in stress tolerance mechanism in soybean. This review also provides a comprehensive catalog of available online omic resources for soybean and its effective utilization. We have also addressed the significance of phenomics in the integrated approaches and recognized high-throughput multi-dimensional phenotyping as a major limiting factor for the improvement of abiotic stress tolerance in soybean.

Soybean (Glycine max) SWEET gene family: insights through comparative genomics, transcriptome profiling and whole genome re-sequence analysis

BackgroundSWEET (MtN3_saliva) domain proteins, a recently identified group of efflux transporters, play an indispensable role in sugar efflux, phloem loading, plant-pathogen interaction and reproductive tissue development. The SWEET gene family is predominantly studied in Arabidopsis and members of the family are being investigated in rice. To date, no transcriptome or genomics analysis of soybean SWEET genes has been reported.ResultsIn the present investigation, we explored the evolutionary aspect of the SWEET gene family in diverse plant species including primitive single cell algae to angiosperms with a major emphasis on Glycine max. Evolutionary features showed expansion and duplication of the SWEET gene family in land plants. Homology searches with BLAST tools and Hidden Markov Model-directed sequence alignments identified 52 SWEET genes that were mapped to 15 chromosomes in the soybean genome as tandem duplication events. Soybean SWEET (GmSWEET) genes showed a wide range of expression profiles in different tissues and developmental stages. Analysis of public transcriptome data and expression profiling using quantitative real time PCR (qRT-PCR) showed that a majority of the GmSWEET genes were confined to reproductive tissue development. Several natural genetic variants (non-synonymous SNPs, premature stop codons and haplotype) were identified in the GmSWEET genes using whole genome re-sequencing data analysis of 106 soybean genotypes. A significant association was observed between SNP-haplogroup and seed sucrose content in three gene clusters on chromosome 6.ConclusionPresent investigation utilized comparative genomics, transcriptome profiling and whole genome re-sequencing approaches and provided a systematic description of soybean SWEET genes and identified putative candidates with probable roles in the reproductive tissue development. Gene expression profiling at different developmental stages and genomic variation data will aid as an important resource for the soybean research community and can be extremely valuable for understanding sink unloading and enhancing carbohydrate delivery to developing seeds for improving yield.

Understanding abiotic stress tolerance mechanisms in soybean: A comparative evaluation of soybean response to drought and flooding stress

Many sources of drought and flooding tolerance have been identified in soybean, however underlying molecular and physiological mechanisms are poorly understood. Therefore, it is important to illuminate different plant responses to these abiotic stresses and understand the mechanisms that confer tolerance. Towards this goal we used four contrasting soybean (Glycine max) genotypes (PI 567690--drought tolerant, Pana--drought susceptible, PI 408105A--flooding tolerant, S99-2281--flooding susceptible) grown under greenhouse conditions and compared genotypic responses to drought and flooding at the physiological, biochemical, and cellular level. We also quantified these variations and tried to infer their role in drought and flooding tolerance in soybean. Our results revealed that different mechanisms contribute to reduction in net photosynthesis under drought and flooding stress. Under drought stress, ABA and stomatal conductance are responsible for reduced photosynthetic rate; while under flooding stress, accumulation of starch granules played a major role. Drought tolerant genotypes PI 567690 and PI 408105A had higher plastoglobule numbers than the susceptible Pana and S99-2281. Drought stress increased the number and size of plastoglobules in most of the genotypes pointing to a possible role in stress tolerance. Interestingly, there were seven fibrillin proteins localized within the plastoglobules that were up-regulated in the drought and flooding tolerant genotypes PI 567690 and PI 408105A, respectively, but down-regulated in the drought susceptible genotype Pana. These results suggest a potential role of Fibrillin proteins, FBN1a, 1b and 7a in soybean response to drought and flooding stress.

Identification and Comparative Analysis of Differential Gene Expression in Soybean Leaf Tissue under Drought and Flooding Stress Revealed by RNA-Seq

Drought and flooding are two major causes of severe yield loss in soybean worldwide. A lack of knowledge of the molecular mechanisms involved in drought and flood stress has been a limiting factor for the effective management of soybeans; therefore, it is imperative to assess the expression of genes involved in response to flood and drought stress. In this study, differentially expressed genes (DEGs) under drought and flooding conditions were investigated using Illumina RNA-Seq transcriptome profiling. A total of 2724 and 3498 DEGs were identified under drought and flooding treatments, respectively. These genes comprise 289 Transcription Factors (TFs) representing Basic Helix-loop Helix (bHLH), Ethylene Response Factors (ERFs), myeloblastosis (MYB), No apical meristem (NAC), and WRKY amino acid motif (WRKY) type major families known to be involved in the mechanism of stress tolerance. The expression of photosynthesis and chlorophyll synthesis related genes were significantly reduced under both types of stresses, which limit the metabolic processes and thus help prolong survival under extreme conditions. However, cell wall synthesis related genes were up-regulated under drought stress and down-regulated under flooding stress. Transcript profiles involved in the starch and sugar metabolism pathways were also affected under both stress conditions. The changes in expression of genes involved in regulating the flux of cell wall precursors and starch/sugar content can serve as an adaptive mechanism for soybean survival under stress conditions. This study has revealed the involvement of TFs, transporters, and photosynthetic genes, and has also given a glimpse of hormonal cross talk under the extreme water regimes, which will aid as an important resource for soybean crop improvement.

Genome-wide transcriptome analysis of soybean primary root under varying water- deficit conditions

BackgroundSoybean is a major crop that provides an important source of protein and oil to humans and animals, but its production can be dramatically decreased by the occurrence of drought stress. Soybeans can survive drought stress if there is a robust and deep root system at the early vegetative growth stage. However, little is known about the genome-wide molecular mechanisms contributing to soybean root system architecture. This study was performed to gain knowledge on transcriptome changes and related molecular mechanisms contributing to soybean root development under water limited conditions.ResultsThe soybean Williams 82 genotype was subjected to very mild stress (VMS), mild stress (MS) and severe stress (SS) conditions, as well as recovery from the severe stress after re-watering (SR). In total, 6,609 genes in the roots showed differential expression patterns in response to different water-deficit stress levels. Genes involved in hormone (Auxin/Ethylene), carbohydrate, and cell wall-related metabolism (XTH/lipid/flavonoids/lignin) pathways were differentially regulated in the soybean root system. Several transcription factors (TFs) regulating root growth and responses under varying water-deficit conditions were identified and the expression patterns of six TFs were found to be common across the stress levels. Further analysis on the whole plant level led to the finding of tissue-specific or water-deficit levels specific regulation of transcription factors. Analysis of the over-represented motif of different gene groups revealed several new cis-elements associated with different levels of water deficit. The expression patterns of 18 genes were confirmed byquantitative reverse transcription polymerase chain reaction method and demonstrated the accuracy and effectiveness of RNA-Seq.ConclusionsThe primary root specific transcriptome in soybean can enable a better understanding of the root response to water deficit conditions. The genes detected in root tissues that were associated with key hormones, carbohydrates, and cell wall-related metabolism could play a vital role in achieving drought tolerance and could be promising candidates for future functional characterization. TFs involved in the soybean root and at the whole plant level could be used for future network analysis between TFs and cis-elements. All of these findings will be helpful in elucidating the molecular mechanisms associated with water stress responses in soybean roots.

Genetic variants in root architecture-related genes in a Glycine soja accession, a potential resource to improve cultivated soybean

BackgroundRoot system architecture is important for water acquisition and nutrient acquisition for all crops. In soybean breeding programs, wild soybean alleles have been used successfully to enhance yield and seed composition traits, but have never been investigated to improve root system architecture. Therefore, in this study, high-density single-feature polymorphic markers and simple sequence repeats were used to map quantitative trait loci (QTLs) governing root system architecture in an inter-specific soybean mapping population developed from a cross between Glycine max and Glycine soja.ResultsWild and cultivated soybean both contributed alleles towards significant additive large effect QTLs on chromosome 6 and 7 for a longer total root length and root distribution, respectively. Epistatic effect QTLs were also identified for taproot length, average diameter, and root distribution. These root traits will influence the water and nutrient uptake in soybean. Two cell division-related genes (D type cyclin and auxin efflux carrier protein) with insertion/deletion variations might contribute to the shorter root phenotypes observed in G. soja compared with cultivated soybean. Based on the location of the QTLs and sequence information from a second G. soja accession, three genes (slow anion channel associated 1 like, Auxin responsive NEDD8-activating complex and peroxidase), each with a non-synonymous single nucleotide polymorphism mutation were identified, which may also contribute to changes in root architecture in the cultivated soybean. In addition, Apoptosis inhibitor 5-like on chromosome 7 and slow anion channel associated 1-like on chromosome 15 had epistatic interactions for taproot length QTLs in soybean.ConclusionRare alleles from a G. soja accession are expected to enhance our understanding of the genetic components involved in root architecture traits, and could be combined to improve root system and drought adaptation in soybean.

Identification of novel QTL governing root architectural traits in an interspecific soybean population.

Cultivated soybean (Glycine max L.) cv. Dunbar (PI 552538) and wild G. soja (PI 326582A) exhibited significant differences in root architecture and root-related traits. In this study, phenotypic variability for root traits among 251 BC2F5 backcross inbred lines (BILs) developed from the cross Dunbar/PI 326582A were identified. The root systems of the parents and BILs were evaluated in controlled environmental conditions using a cone system at seedling stage. The G. max parent Dunbar contributed phenotypically favorable alleles at a major quantitative trait locus on chromosome 8 (Satt315-I locus) that governed root traits (tap root length and lateral root number) and shoot length. This QTL accounted for >10% of the phenotypic variation of both tap root and shoot length. This QTL region was found to control various shoot- and root-related traits across soybean genetic backgrounds. Within the confidence interval of this region, eleven transcription factors (TFs) were identified. Based on RNA sequencing and Affymetrix expression data, key TFs including MYB, AP2-EREBP and bZIP TFs were identified in this QTL interval with high expression in roots and nodules. The backcross inbred lines with different parental allelic combination showed different expression pattern for six transcription factors selected based on their expression pattern in root tissues. It appears that the marker interval Satt315–I locus on chromosome 8 contain an essential QTL contributing to early root and shoot growth in soybean.

Core clock, SUB1, and ABAR genes mediate flooding and drought responses via alternative splicing in soybean

Circadian clocks are a great evolutionary innovation and provide competitive advantage during the day/night cycle and under changing environmental conditions. The circadian clock mediates expression of a large proportion of genes in plants, achieving a harmonious relationship between energy metabolism, photosynthesis, and biotic and abiotic stress responses. Here it is shown that multiple paralogues of clock genes are present in soybean (Glycine max) and mediate flooding and drought responses. Differential expression of many clock and SUB1 genes was found under flooding and drought conditions. Furthermore, natural variation in the amplitude and phase shifts in PRR7 and TOC1 genes was also discovered under drought and flooding conditions, respectively. PRR3 exhibited flooding- and drought-specific splicing patterns and may work in concert with PRR7 and TOC1 to achieve energy homeostasis under flooding and drought conditions. Higher expression of TOC1 also coincides with elevated levels of abscisic acid (ABA) and variation in glucose levels in the morning and afternoon, indicating that this response to abiotic stress is mediated by ABA, endogenous sugar levels, and the circadian clock to fine-tune photosynthesis and energy utilization under stress conditions. It is proposed that the presence of multiple clock gene paralogues with variation in DNA sequence, phase, and period could be used to screen exotic germplasm to find sources for drought and flooding tolerance. Furthermore, fine tuning of multiple clock gene paralogues (via a genetic engineering approach) should also facilitate the development of flooding- and drought-tolerant soybean varieties.

Comparative analysis of the drought-responsive transcriptome in soybean lines contrasting for canopy wilting

Drought stress causes significant yield losses in major oil seed crops, such as soybean [Glycine max (L.) Merr]. Few soybean lines have been identified as canopy-wilting tolerant; however, the molecular mechanism conferring tolerance is not fully understood. To understand the biological process, a whole genome transcriptome analysis was performed for leaf tissues of two contrasting soybean lines: drought-susceptible (DS) Pana and drought-tolerant (DT) PI 567690. A pairwise comparison of the DS and DT lines under drought and control conditions detected 1914 and 670 genes with a greater than two-fold change in expression under drought conditions. Pairwise treatment comparison and gene enrichment analysis on the DT line showed the down-regulation of genes associated with protein binding, hydrolase activity, carbohydrate/lipid metabolism, xyloglucan endo-transglycosylases associated with cell-wall, apoplast, and chlorophyll a/b binding proteins. On the other hand, genes that were associated with the biotic stress response, ion binding and transport, the oxido-reductive process and electron carrier activity were up-regulated. Gene enrichment analysis detected UDP glucuronosyl transferase activity-encoding genes to be differentially expressed in PI 567690 under drought stress conditions. We found valuable SNPs variation in aquaporin genes of the DT line that are conserved in known slower canopy-wilting lines, this should facilitate marker-assisted selection in soybeans with improved drought tolerance.