Transcriptional Profile of Soybean Seeds with Contrasting Seed Coat Color.

Soybean is the primary source of vegetable protein and is used for various purposes, mainly to feed animals. This crop can have diverse seed coat colors, varying from yellow, black, brown, and green to bicolor. Black seed coat cultivars have already been assigned as favorable for both seed and grain production. Thus, this work aimed to identify genes associated with soybean seed quality by comparing the transcriptomes of soybean seeds with contrasting seed coat colors. The results from RNA-seq analyses were validated with real-time PCR using the cultivar BRS 715A (black seed coat) and the cultivars BRS 413 RR and DM 6563 IPRO (yellow seed coat). We found 318 genes differentially expressed in all cultivars (freshly harvested seeds and seeds stored in cold chamber). From the in silico analysis of the transcriptomes, the following genes were selected and validated with RT-qPCR: ACS1, ACSF3, CYP90A1, CYP710A1, HCT, CBL, and SAHH. These genes are genes induced in the black seed coat cultivar and are part of pathways responsible for ethylene, lipid, brassinosteroid, lignin, and sulfur amino acid biosynthesis. The BRSMG 715A gene has almost 4times more lignin than the yellow seed coat cultivars. These attributes are related to the BRSMG 715A cultivar's higher seed quality, which translates to more longevity and resistance to moisture and mechanical damage. Future silencing studies may evaluate the knockout of these genes to better understand the biology of soybean seeds with black seed coat.

Importance of inositols and their derivatives in cowpea under root dehydration: An omics perspective.

This work presents a robust analysis of the inositols (INSs) and raffinose family oligosaccharides (RFOs) pathways, using genomic and transcriptomic tools in cowpea under root dehydration. Nineteen (~70%) of the 26 scrutinized enzymes presented transcriptional up-regulation in at least one treatment time. The transcriptional orchestration allowed categorization of the analyzed enzymes as time-independent (those showing the same regulation throughout the assay) and time-dependent (those showing different transcriptional regulation over time). It is suggested that up-regulated time-independent enzymes (INSs: myo-inositol oxygenase, inositol-tetrakisphosphate 1-kinase 3, phosphatidylinositol 4-phosphate 5-kinase 4-like, 1-phosphatidylinositol-3-phosphate 5-kinase, phosphoinositide phospholipase C, and non-specific phospholipase C; RFOs: α-galactosidase, invertase, and raffinose synthase) actively participate in the reorganization of cowpea molecular physiology under the applied stress. In turn, time-dependent enzymes, especially those up-regulated in some of the treatment times (INSs: inositol-pentakisphosphate 2-kinase, phosphatidylinositol 4-kinase, phosphatidylinositol synthase, multiple inositol polyphosphate phosphatase 1, methylmalonate-semialdehyde dehydrogenase, triosephosphate isomerase, myo-inositol-3-phosphate synthase, phosphatidylinositol 3,4,5-trisphosphate 3-phosphatase and protein-tyrosine-phosphatase, and phosphatidylinositol 3-kinase; RFOs: galactinol synthase) seem to participate in fine-tuning of the molecular physiology, helping the cowpea plants to acclimatize under dehydration stress. Not all loci encoding the studied enzymes were expressed during the assay; most of the expressed ones exhibited a variable transcriptional profile in the different treatment times. Genes of the INSs and RFOs pathways showed high orthology with analyzed Phaseoleae members, suggesting a relevant role within this legume group. Regarding the promoter regions of INSs and RFOs genes, some bona fide cis-regulatory elements were identified in association with seven transcription factor families (AP2-EFR, Dof-type, MADS-box, bZIP, CPP, ZF-HD, and GATA-type). Members of INSs and RFOs pathways potentially participate in other processes regulated by these proteins in cowpea.

Insights into soybean transcriptome reconfiguration under hypoxic stress: Functional, regulatory, structural, and compositional characterization.

Soybean (Glycine max) is one of the major crops worldwide and flooding stress affects the production and expansion of cultivated areas. Oxygen is essential for mitochondrial aerobic respiration to supply the energy demand of plant cells. Because oxygen diffusion in water is 10,000 times lower than in air, partial (hypoxic) or total (anoxic) oxygen deficiency is important component of flooding. Even when oxygen is externally available, oxygen deficiency frequently occurs in bulky, dense or metabolically active tissues such as phloem, meristems, seeds, and fruits. In this study, we analyzed conserved and divergent root transcriptional responses between flood-tolerant Embrapa 45 and flood-sensitive BR 4 soybean cultivars under hypoxic stress conditions with RNA-seq. To understand how soybean genes evolve and respond to hypoxia, stable and differentially expressed genes were characterized structurally and compositionally comparing its mechanistic relationship. Between cultivars, Embrapa 45 showed less up- and more down-regulated genes, and stronger induction of phosphoglucomutase (Glyma05g34790), unknown protein related to N-terminal protein myristoylation (Glyma06g03430), protein suppressor of phyA-105 (Glyma06g37080), and fibrillin (Glyma10g32620). RNA-seq and qRT-PCR analysis of non-symbiotic hemoglobin (Glyma11g12980) indicated divergence in gene structure between cultivars. Transcriptional changes for genes in amino acids and derivative metabolic process suggest involvement of amino acids metabolism in tRNA modifications, translation accuracy/efficiency, and endoplasmic reticulum stress in both cultivars under hypoxia. Gene groups differed in promoter TATA box, ABREs (ABA-responsive elements), and CRT/DREs (C-repeat/dehydration-responsive elements) frequency. Gene groups also differed in structure, composition, and codon usage, indicating biological significances. Additional data suggests that cis-acting ABRE elements can mediate gene expression independent of ABA in soybean roots under hypoxia.

Functional Characterization of a Putative Glycine max ELF4 in Transgenic Arabidopsis and Its Role during Flowering Control.

Flowering is an important trait in major crops like soybean due to its direct relation to grain production. The circadian clock mediates the perception of seasonal changes in day length and temperature to modulate flowering time. The circadian clock gene EARLY FLOWERING 4 (ELF4) was identified in Arabidopsis thaliana and is believed to play a key role in the integration of photoperiod, circadian regulation, and flowering. The molecular circuitry that comprises the circadian clock and flowering control in soybeans is just beginning to be understood. To date, insufficient information regarding the soybean negative flowering regulators exist, and the biological function of the soybean ELF4 (GmELF4) remains unknown. Here, we investigate the ELF4 family members in soybean and functionally characterize a GmELF4 homologous gene. The constitutive overexpression of GmELF4 delayed flowering in Arabidopsis, showing the ELF4 functional conservation among plants as part of the flowering control machinery. We also show that GmELF4 alters the expression of Arabidopsis key flowering time genes (AtCO and AtFT), and this down-regulation is the likely cause of flowering delay phenotypes. Furthermore, we identified the GmELF4 network genes to infer the participation of GmELF4 in soybeans. The data generated in this study provide original insights for comprehending the role of the soybean circadian clock ELF4 gene as a negative flowering controller.

Characterization of Soybean Genetically Modified for Drought Tolerance in Field Conditions.

Drought is one of the most stressful environmental factor causing yield and economic losses in many soybean-producing regions. In the last decades, transcription factors (TFs) are being used to develop genetically modified plants more tolerant to abiotic stresses. Dehydration responsive element binding (DREB) and ABA-responsive element-binding (AREB) TFs were introduced in soybean showing improved drought tolerance, under controlled conditions. However, these results may not be representative of the way in which plants behave over the entire season in the real field situation. Thus, the objectives of this study were to analyze agronomical traits and physiological parameters of AtDREB1A (1Ab58), AtDREB2CA (1Bb2193), and AtAREB1 (1Ea2939) GM lines under irrigated (IRR) and non-irrigated (NIRR) conditions in a field experiment, over two crop seasons and quantify transgene and drought-responsive genes expression. Results from season 2013/2014 revealed that line 1Ea2939 showed higher intrinsic water use and leaf area index. Lines 1Ab58 and 1Bb2193 showed a similar behavior to wild-type plants in relation to chlorophyll content. Oil and protein contents were not affected in transgenic lines in NIRR conditions. Lodging, due to plentiful rain, impaired yield from the 1Ea2939 line in IRR conditions. qPCR results confirmed the expression of the inserted TFs and drought-responsive endogenous genes. No differences were identified in the field experiment performed in crop season 2014/2015, probably due to the optimum rainfall volume during the cycle. These field screenings showed promising results for drought tolerance. However, additional studies are needed in further crop seasons and other sites to better characterize how these plants may outperform the WT under field water deficit.

Genomic and transcriptomic characterization of the transcription factor family R2R3-MYB in soybean and its involvement in the resistance responses to Phakopsora pachyrhizi.

Myb genes constitute one of the largest transcription factor families in the plant kingdom. Soybean MYB transcription factors have been related to the plant response to biotic stresses. Their involvement in response to Phakopsora pachyrhizi infection has been reported by several transcriptional studies. Due to their apparently highly diverse functions, these genes are promising targets for developing crop varieties resistant to diseases. In the present study, the identification and phylogenetic analysis of the soybean R2R3-MYB (GmMYB) transcription factor family was performed and the expression profiles of these genes under biotic stress were determined. GmMYBs were identified from the soybean genome using bioinformatic tools, and their putative functions were determined based on the phylogenetic tree and classified into subfamilies using guides AtMYBs describing known functions. The transcriptional profiles of GmMYBs upon infection with different pathogen were revealed by in vivo and in silico analyses. Selected target genes potentially involved in disease responses were assessed by RT-qPCR after different times of inoculation with P. pachyrhizi using different genetic backgrounds related to resistance genes (Rpp2 and Rpp5). R2R3-MYB transcription factors related to lignin synthesis and genes responsive to chitin were significantly induced in the resistant genotypes.

Genome-wide analysis of the Hsp20 gene family in soybean: comprehensive sequence, genomic organization and expression profile analysis under abiotic and biotic stresses.

BACKGROUND: The Hsp20 genes are associated with stress caused by HS and other abiotic factors, but have recently been found to be associated with the response to biotic stresses. These genes represent the most abundant class among the HSPs in plants, but little is known about this gene family in soybean. Because of their apparent multifunctionality, these proteins are promising targets for developing crop varieties that are better adapted to biotic and abiotic stresses. Thus, in the present study an in silico identification of GmHsp20 gene family members was performed, and the genes were characterized and subjected to in vivo expression analysis under biotic and abiotic stresses. RESULTS: A search of the available soybean genome databases revealed 51 gene models as potential GmHsp20 candidates. The 51 GmHsp20 genes were distributed across a total of 15 subfamilies where a specific predicted secondary structure was identified. Based on in vivo analysis, only 47 soybean Hsp20 genes were responsive to heat shock stress. Among the GmHsp20 genes that were potentials HSR, five were also cold-induced, and another five, in addition to one GmAcd gene, were responsive to Meloidogyne javanica infection. Furthermore, one predicted GmHsp20 was shown to be responsive only to nematode infection; no expression change was detected under other stress conditions. Some of the biotic stress-responsive GmHsp20 genes exhibited a divergent expression pattern between resistant and susceptible soybean genotypes under M. javanica infection. The putative regulatory elements presenting some conservation level in the GmHsp20 promoters included HSE, W-box, CAAT box, and TA-rich elements. Some of these putative elements showed a unique occurrence pattern among genes responsive to nematode infection. CONCLUSIONS: The evolution of Hsp20 family in soybean genome has most likely involved a total of 23 gene duplications. The obtained expression profiles revealed that the majority of the 51 GmHsp20 candidates are induced under HT, but other members of this family could also be involved in normal cellular functions, unrelated to HT. Some of the GmHsp20 genes might be specialized to respond to nematode stress, and the predicted promoter structure of these genes seems to have a particular conserved pattern related to their biological function.

The Lesion Simulating Disease (LSD) gene family as a variable in soybean response to Phakopsora pachyrhizi infection and dehydration.

The Lesion Simulating Disease (LSD) genes encode a family of zinc finger proteins that are reported to play an important role in the hypersensitive response and programmed cell death (PCD) that are caused by biotic and abiotic stresses. In the present study, 117 putative LSD family members were identified in Viridiplantae. Genes with one, two, or three conserved LSD domains were identified. Proteins with three LSD domains were highly represented in the species analyzed and were present in basal organisms. Proteins with two LSD domains were identified only in the Embryophyte clade, and proteins possessing one LSD domain were highly represented in grass species. Expression analyses of Glycine max LSD (GmLSD) genes were performed by real-time quantitative polymerase chain reaction. The results indicated that GmLSD genes are not ubiquitously expressed in soybean organs and that their expression patterns are instead organ-dependent. The expression of the majority of GmLSD genes is modulated in soybean during Phakopsora pachyrhizi infection. In addition, the expression of some GmLSD genes is modulated in plants under dehydration stress. These results suggest the involvement of GmLSD genes in the response of soybean to both biotic and abiotic stresses.

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.

In silico identification of known osmotic stress responsive genes from Arabidopsis in soybean and Medicago.

Plants experience various environmental stresses, but tolerance to these adverse conditions is a very complex phenomenon. The present research aimed to evaluate a set of genes involved in osmotic response, comparing soybean and medicago with the well-described Arabidopsis thaliana model plant. Based on 103 Arabidopsis proteins from 27 categories of osmotic stress response, comparative analyses against Genosoja and Medicago truncatula databases allowed the identification of 1,088 soybean and 1,210 Medicago sequences. The analysis showed a high number of sequences and high diversity, comprising genes from all categories in both organisms. Genes with unknown function were among the most representative, followed by transcription factors, ion transport proteins, water channel, plant defense, protein degradation, cellular structure, organization & biogenesis and senescence. An analysis of sequences with unknown function allowed the annotation of 174 soybean and 217 Medicago sequences, most of them concerning transcription factors. However, for about 30% of the sequences no function could be attributed using in silico procedures. The establishment of a gene set involved in osmotic stress responses in soybean and barrel medic will help to better understand the survival mechanisms for this type of stress condition in legumes.

Identification of novel soybean microRNAs involved in abiotic and biotic stresses.

BACKGROUND: Small RNAs (19-24 nt) are key regulators of gene expression that guide both transcriptional and post-transcriptional silencing mechanisms in eukaryotes. Current studies have demonstrated that microRNAs (miRNAs) act in several plant pathways associated with tissue proliferation, differentiation, and development and in response to abiotic and biotic stresses. In order to identify new miRNAs in soybean and to verify those that are possibly water deficit and rust-stress regulated, eight libraries of small RNAs were constructed and submitted to Solexa sequencing. RESULTS: The libraries were developed from drought-sensitive and tolerant seedlings and rust-susceptible and resistant soybeans with or without stressors. Sequencing the library and subsequent analyses detected 256 miRNAs. From this total, we identified 24 families of novel miRNAs that had not been reported before, six families of conserved miRNAs that exist in other plants species, and 22 families previously reported in soybean. We also observed the presence of several isomiRNAs during our analyses. To validate novel miRNAs, we performed RT-qPCR across the eight different libraries. Among the 11 miRNAs analyzed, all showed different expression profiles during biotic and abiotic stresses to soybean. The majority of miRNAs were up-regulated during water deficit stress in the sensitive plants. However, for the tolerant genotype, most of the miRNAs were down regulated. The pattern of miRNAs expression was also different for the distinct genotypes submitted to the pathogen stress. Most miRNAs were down regulated during the fungus infection in the susceptible genotype; however, in the resistant genotype, most miRNAs did not vary during rust attack. A prediction of the putative targets was carried out for conserved and novel miRNAs families. CONCLUSIONS: Validation of our results with quantitative RT-qPCR revealed that Solexa sequencing is a powerful tool for miRNA discovery. The identification of differentially expressed plant miRNAs provides molecular evidence for the possible involvement of miRNAs in the process of water deficit- and rust-stress responses.

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.

Differential gene expression and mitotic cell analysis of the drought tolerant soybean (Glycine max L. Merrill Fabales, Fabaceae) cultivar MG/BR46 (Conquista) under two water deficit induction systems

Drought cause serious yield losses in soybean (Glycine max), roots being the first plant organ to detect the water-stress signals triggering defense mechanisms. We used two drought induction systems to identify genes differentially expressed in the roots of the drought-tolerant soybean cultivar MG/BR46 (Conquista) and characterize their expression levels during water deficit. Soybean plants grown in nutrient solution hydroponically and in sand-pots were submitted to water stress and gene expression analysis was conducted using the differential display (DD) and real time polymerase chain reaction (PCR) techniques. Three differentially expressed mRNA transcripts showed homology to the Antirrhinum majus basic helix-loop-helix transcription factor bHLH, the Arabidopsis thaliana phosphatidylinositol transfer protein PITP and the auxin-independent growth regulator 1 (axi 1). The hydroponic experiments showed that after 100 min outside the nutrient solution photosynthesis completely stopped, stomata closed and leaf temperature rose. Both stress induction treatments produced significant decrease in the mitotic indices of root cells. Axi 1, PITP and bHLH were not only differentially expressed during dehydration in the hydroponics experiments but also during induced drought in the pot experiments. Although, there were differences between the two sets of experiments in the time at which up or down regulation occurred, the expression pattern of all three transcripts was related. Similar gene expression and cytological analysis results occurred in both systems, suggesting that hydroponics could be used to simulate drought detection by roots growing in soil and thus facilitate rapid and easy root sampling.

Population structure of the Brazilian southern green stink bug, Nezara viridula.

The Southern Green Stink Bug, Nezara viridula (L.) (Heteroptera: Pentatomidae), is a cosmopolitan and economically important pest to several crops. Studies on N. viridula migration and population structure have been neglected. We studied geographically distinct Brazilian N. viridula populations to assess their variability and to determine gene flow among them. DNA from specimens collected on soybean fields were subjected to RAPD analysis to determine genetic similarity and population structure parameters. All N. viridula populations studied were genetically distinct from the others. The maximum similarity occurred between populations from Londrina and Sertanópolis (Parana State). The Cruz Alta population was the most divergent from the others. Despite the short distance between Cambé and Londrina (ca. 29 km), and the absence of geographic barriers, both populations clustered in different groups and the estimated gene flow index (Nm) among them was 2.02, indicating relatively restricted migration. The estimated overall index, Nm was 1.41 suggesting that N. viridula is a better flier than the Neotropical Brown stink bug, Euschistus heros (Nm =0.83).

VSQual: a visual system to assist DNA sequencing quality control.

A lack of pliant software tools that support small- to medium-scale DNA sequencing efforts is a major hindrance for recording and using laboratory workflow information to monitor the overall quality of data production. Here we describe VSQual, a set of Perl programs intended to provide simple and powerful tools to check several quality features of the sequencing data generated by automated DNA sequencing machines. The core program of VSQual is a flexible Perl-based pipeline, designed to be accessible and useful for both programmers and non-programmers. This pipeline directs the processing steps and can be easily customized for laboratory needs. Basically, the raw DNA sequencing trace files are processed by Phred and Cross_match, then the outputs are parsed, reformatted into Web-based graphical reports, and added to a Web site structure. The result is a set of real time sequencing reports easily accessible and understood by common laboratory people. These reports facilitate the monitoring of DNA sequencing as well as the management of laboratory workflow, significantly reducing operational costs and ensuring high quality and scientifically reliable results.

VSQual: a visual system to assist DNA sequencing quality control

A lack of pliant software tools that support small- to medium-scale DNA sequencing efforts is a major hindrance for recording and using laboratory workflow information to monitor the overall quality of data production. Here we describe VSQual, a set of Perl programs intended to provide simple and powerful tools to check several quality features of the sequencing data generated by automated DNA sequencing machines. The core program of VSQual is a flexible Perl-based pipeline, designed to be accessible and useful for both programmers and non-programmers. This pipeline directs the processing steps and can be easily customized for laboratory needs. Basically, the raw DNA sequencing trace files are processed by Phred and Cross_match, then the outputs are parsed, reformatted into Web-based graphical reports, and added to a Web site structure. The result is a set of real time sequencing reports easily accessible and understood by common laboratory people. These reports facilitate the monitoring of DNA sequencing as well as the management of laboratory workflow, significantly reducing operational costs and ensuring high quality and scientifically reliable results.