Francismar Corrêa Marcelino-Guimarães

Identification of novel soybean microRNAs involved in abiotic and biotic stresses

BackgroundSmall 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.ResultsThe 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.ConclusionsValidation 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.

The use of microRNAs as reference genes for quantitative polymerase chain reaction in soybean

Reverse transcription quantitative real-time polymerase chain reaction (RT-qPCR) is a robust and widely applied technique used to investigate gene expression. However, for correct analysis and interpretation of results, the choice of a suitable gene to use as an internal control is a crucial factor. These genes, such as housekeeping genes, should have a constant expression level in different tissues and across different conditions. The advances in genome sequencing have provided high-throughput gene expression analysis and have contributed to the identification of new genes, including microRNAs (miRNAs). The miRNAs are fundamental regulatory genes of eukaryotic genomes, acting on several biological functions. In this study, miRNA expression stability was investigated in different soybean tissues and genotypes as well as after abiotic or biotic stress treatments. The present study represents the first investigation into the suitability of miRNAs as housekeeping genes in plants. The transcript stability of 10 miRNAs was compared to those of six previously reported housekeeping genes for the soybean. In this study, we provide evidence that the expression stabilities of miR156b and miR1520d were the highest across the soybean experiments. Furthermore, these miRNAs genes were more stable than the most commonly protein-coding genes used in soybean gene expression studies involving RT-qPCR.

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

BackgroundThe Hsp 20 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 GmHsp 20 gene family members was performed, and the genes were characterized and subjected to in vivo expression analysis under biotic and abiotic stresses.ResultsA search of the available soybean genome databases revealed 51 gene models as potential GmHsp 20 candidates. The 51 GmHsp 20 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 Hsp 20 genes were responsive to heat shock stress. Among the GmHsp 20 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 GmHsp 20 was shown to be responsive only to nematode infection; no expression change was detected under other stress conditions. Some of the biotic stress-responsive GmHsp 20 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 GmHsp 20 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.ConclusionsThe evolution of Hsp 20 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 GmHsp 20 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 GmHsp 20 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.

Genome-wide annotation of the soybean WRKY family and functional characterization of genes involved in response to Phakopsora pachyrhizi infection

BackgroundMany previous studies have shown that soybean WRKY transcription factors are involved in the plant response to biotic and abiotic stresses. Phakopsora pachyrhizi is the causal agent of Asian Soybean Rust, one of the most important soybean diseases. There are evidences that WRKYs are involved in the resistance of some soybean genotypes against that fungus. The number of WRKY genes already annotated in soybean genome was underrepresented. In the present study, a genome-wide annotation of the soybean WRKY family was carried out and members involved in the response to P. pachyrhizi were identified.ResultsAs a result of a soybean genomic databases search, 182 WRKY-encoding genes were annotated and 33 putative pseudogenes identified. Genes involved in the response to P. pachyrhizi infection were identified using superSAGE, RNA-Seq of microdissected lesions and microarray experiments. Seventy-five genes were differentially expressed during fungal infection. The expression of eight WRKY genes was validated by RT-qPCR. The expression of these genes in a resistant genotype was earlier and/or stronger compared with a susceptible genotype in response to P. pachyrhizi infection. Soybean somatic embryos were transformed in order to overexpress or silence WRKY genes. Embryos overexpressing a WRKY gene were obtained, but they were unable to convert into plants. When infected with P. pachyrhizi, the leaves of the silenced transgenic line showed a higher number of lesions than the wild-type plants.ConclusionsThe present study reports a genome-wide annotation of soybean WRKY family. The participation of some members in response to P. pachyrhizi infection was demonstrated. The results contribute to the elucidation of gene function and suggest the manipulation of WRKYs as a strategy to increase fungal resistance in soybean plants.

Identification of reference genes for expression analysis by real-time quantitative PCR in drought-stressed soybean

The objective of this work was to validate, by quantitative PCR in real time (RT-qPCR), genes to be used as reference in studies of gene expression in soybean in drought-stressed trials. Four genes commonly used in soybean were evaluated: Gmβ-actin, GmGAPDH, GmLectin and GmRNAr18S. Total RNA was extracted from six samples: three from roots in a hydroponic system with different drought intensities (0, 25, 50, 75 and 100 minutes of water stress), and three from leaves of plants grown in sand with different soil moistures (15, 5 and 2.5% gravimetric humidity). The raw cycle threshold (Ct) data were analyzed, and the efficiency of each primer was calculated for an overall analysis of the Ct range among the different samples. The GeNorm application was used to evaluate the best reference gene, according to its stability. The GmGAPDH was the least stable gene, with the highest mean values of expression stability (M), and the most stable genes, with the lowest M values, were the Gmβ-actin and GmRNAr18S, when both root and leaves samples were tested. These genes can be used in RT-qPCR as reference gene for expression analysis.

Positive and Negative Roles for Soybean MPK6 in Regulating Defense Responses

It has been well established that MPK6 is a positive regulator of defense responses in model plants such as Arabidopsis and tobacco. However, the functional importance of soybean MPK6 in disease resistance has not been investigated. Here, we showed that silencing of GmMPK6 in soybean using virus-induced gene silencing mediated by Bean pod mottle virus (BPMV) caused stunted growth and spontaneous cell death on the leaves, a typical phenotype of activated defense responses. Consistent with this phenotype, expression of pathogenesis-related (PR) genes and the conjugated form of salicylic acid were significantly increased in GmMPK6-silenced plants. As expected, GmMPK6-silenced plants were more resistant to downy mildew and Soybean mosaic virus compared with vector control plants, indicating a negative role of GmMPK6 in disease resistance. Interestingly, overexpression of GmMPK6, either transiently in Nicotiana benthamiana or stably in Arabidopsis, resulted in hypersensitive response (HR)-like cell death. The HR-like cell death was accompanied by increased PR gene expression, suggesting that GmMPK6, like its counterpart in other plant species, also plays a positive role in cell death induction and defense response. Using bimolecular fluorescence complementation analysis, we determined that GmMKK4 might function upstream of GmMPK6 and GmMKK4 could interact with GmMPK6 independent of its phosphorylation status. Taken together, our results indicate that GmMPK6 functions as both repressor and activator in defense responses of soybean.

Expression patterns of GmAP2/EREB-like transcription factors involved in soybean responses to water deficit.

Soybean farming has faced several losses in productivity due to drought events in the last few decades. However, plants have molecular mechanisms to prevent and protect against water deficit injuries, and transcription factors play an important role in triggering different defense mechanisms. Understanding the expression patterns of transcription factors in response to water deficit and to environmental diurnal changes is very important for unveiling water deficit stress tolerance mechanisms. Here, we analyzed the expression patterns of ten APETALA2/Ethylene Responsive Element Binding-like (AP2/EREB-like) transcription factors in two soybean genotypes (BR16: drought-sensitive; and Embrapa 48: drought-tolerant). According to phylogenetic and domain analyses, these genes can be included in the DREB and ERF subfamilies. We also analyzed a GmDRIP-like gene that encodes a DREB negative regulator. We detected the up-regulation of 9 GmAP2/EREB-like genes and identified transcriptional differences that were dependent on the levels of the stress applied and the tissue type analyzed (the expression of the GmDREB1F-like gene, for example, was four times higher in roots than in leaves). The GmDRIP-like gene was not induced by water deficit in BR16 during the longest periods of stress, but was significantly induced in Embrapa 48; this suggests a possible genetic/molecular difference between the responses of these cultivars to water deficit stress. Additionally, RNAseq gene expression analysis over a 24-h time course indicates that the expression patterns of several GmDREB-like genes are subject to oscillation over the course of the day, indicating a possible circadian regulation.

Introduction of the rd29A: AtDREB2A CA gene into soybean (Glycine max L. Merril) and its molecular characterization in leaves and roots during dehydration

The loss of soybean yield to Brazilian producers because of a water deficit in the 2011–2012 season was 12.9%. To reduce such losses, molecular biology techniques, including plant transformation, can be used to insert genes of interest into conventional soybean cultivars to produce lines that are more tolerant to drought. The abscisic acid (ABA)-independent Dehydration Responsive Element Binding (DREB) gene family has been used to obtain plants with increased tolerance to abiotic stresses. In the present study, the rd29A:AtDREB2A CA gene from Arabidopsis thaliana was inserted into soybean using biolistics. Seventy-eight genetically modified (GM) soybean lines containing 2–17 copies of the AtDREB2A CA gene were produced. Two GM soybean lines (P1397 and P2193) were analyzed to assess the differential expression of the AtDREB2A CA transgene in leaves and roots submitted to various dehydration treatments. Both GM lines exhibited high expression of the transgene, with the roots of P2193 showing the highest expression levels during water deficit. Physiological parameters examined during water deficit confirmed the induction of stress. This analysis of AtDREB2A CA expression in GM soybean indicated that line P2193 had the greatest stability and highest expression in roots during water deficit-induced stress.

Overall picture of expressed Heat Shock Factors in Glycine max, Lotus japonicus and Medicago truncatula

Heat shock (HS) leads to the activation of molecular mechanisms, known as HS-response, that prevent damage and enhance survival under stress. Plants have a flexible and specialized network of Heat Shock Factors (HSFs), which are transcription factors that induce the expression of heat shock proteins. The present work aimed to identify and characterize the Glycine max HSF repertory in the Soybean Genome Project (GENOSOJA platform), comparing them with other legumes (Medicago truncatula and Lotus japonicus) in view of current knowledge of Arabidopsis thaliana. The HSF characterization in leguminous plants led to the identification of 25, 19 and 21 candidate ESTs in soybean, Lotus and Medicago, respectively. A search in the SuperSAGE libraries revealed 68 tags distributed in seven HSF gene types. From the total number of obtained tags, more than 70% were related to root tissues (water deficit stress libraries vs. controls), indicating their role in abiotic stress responses, since the root is the first tissue to sense and respond to abiotic stress. Moreover, as heat stress is related to the pressure of dryness, a higher HSF expression was expected at the water deficit libraries. On the other hand, expressive HSF candidates were obtained from the library inoculated with Asian Soybean Rust, inferring crosstalk among genes associated with abiotic and biotic stresses. Evolutionary relationships among sequences were consistent with different HSF classes and subclasses. Expression profiling indicated that regulation of specific genes is associated with the stage of plant development and also with stimuli from other abiotic stresses pointing to the maintenance of HSF expression at a basal level in soybean, favoring its activation under heat-stress conditions.

Evaluation of genetic variation among Brazilian soybean cultivars through genome resequencing

BackgroundSoybean [Glycine max (L.) Merrill] is one of the most important legumes cultivated worldwide, and Brazil is one of the main producers of this crop. Since the sequencing of its reference genome, interest in structural and allelic variations of cultivated and wild soybean germplasm has grown. To investigate the genetics of the Brazilian soybean germplasm, we selected soybean cultivars based on the year of commercialization, geographical region and maturity group and resequenced their genomes.ResultsWe resequenced the genomes of 28 Brazilian soybean cultivars with an average genome coverage of 14.8X. A total of 5,835,185 single nucleotide polymorphisms (SNPs) and 1,329,844 InDels were identified across the 20 soybean chromosomes, with 541,762 SNPs, 98,922 InDels and 1,093 CNVs that were exclusive to the 28 Brazilian cultivars. In addition, 668 allelic variations of 327 genes were shared among all of the Brazilian cultivars, including genes related to DNA-dependent transcription-elongation, photosynthesis, ATP synthesis-coupled electron transport, cellular respiration, and precursors of metabolite generation and energy. A very homogeneous structure was also observed for the Brazilian soybean germplasm, and we observed 41 regions putatively influenced by positive selection. Finally, we detected 3,880 regions with copy-number variations (CNVs) that could help to explain the divergence among the accessions evaluated.ConclusionsThe large number of allelic and structural variations identified in this study can be used in marker-assisted selection programs to detect unique SNPs for cultivar fingerprinting. The results presented here suggest that despite the diversification of modern Brazilian cultivars, the soybean germplasm remains very narrow because of the large number of genome regions that exhibit low diversity. These results emphasize the need to introduce new alleles to increase the genetic diversity of the Brazilian germplasm.