Search

Bacterial pustule (BP), caused by Xanthomonas citri pv. glycines , is an important disease that, under favorable conditions, can drastically affect soybean production. We performed a genome-wide association study (GWAS) with a panel containing Brazilian and American cultivars, which were screened qualitatively and quantitatively against two Brazilian X. citri isolates (IBS 333 and IBS 327). The panel was genotyped using a genotyping by sequencing (GBS) approach, and we identified two main new regions in soybeans associated with X. citri resistance on chromosomes 6 (IBS 333) and 18 (IBS 327), different from the traditional rxp gene located on chromosome 17. The region on chromosome 6 was also detected by QTL mapping using a biparental cross between Williams 82 (R) and PI 416937 (S), showing that Williams 82 has another recessive resistance gene besides rxp , which was also detected in nine BP-resistant ancestors of the Brazilian cultivars (including CNS, S-100), based on haplotype analysis. Furthermore, we identified additional SNPs in strong LD (0.8) with peak SNPs by exploring variation available in WGS (whole genome sequencing) data among 31 soybean accessions. In these regions in strong LD, two candidate resistance genes were identified (Glyma.06g311000 and Glyma.18g025100) for chromosomes 6 and 18, respectively. Therefore, our results allowed the identification of new chromosomal regions in soybeans associated with BP disease, which could be useful for marker-assisted selection and will enable a reduction in time and cost for the development of resistant cultivars.

Soybean is a crucial crop for the Brazilian economy, but it faces challenges from the biotrophic fungus Phakopsora pachyrhizi, which causes Asian Soybean Rust (ASR). In this study, we aimed to identify SNPs associated with resistance within the Rpp1 locus, which is effective against Brazilian ASR populations. We employed GWAS and re-sequencing analyzes to pinpoint SNP markers capable of differentiating between soybean accessions harboring the Rpp1, Rpp1-b and other alternative alleles in the Rpp1 locus and from susceptible soybean cultivars. Seven SNP markers were found to be associated with ASR resistance through GWAS, with three of them defining haplotypes that efficiently distinguished the accessions based on their ASR resistance and source of the Rpp gene. These haplotypes were subsequently validated using a bi-parental population and a diverse set of Rpp sources, demonstrating that the GWAS markers co-segregate with ASR resistance. We then examined the presence of these haplotypes in a diverse set of soybean genomes worldwide, finding a few new potential sources of Rpp1/Rpp1-b. Further genomic sequence analysis revealed nucleotide differences within the genes present in the Rpp1 locus, including the ULP1-NBS-LRR genes, which are potential R gene candidates. These results provide valuable insights into ASR resistance in soybean, thus helping the development of resistant soybean varieties through genetic breeding programs., (© 2024. The Author(s).)

Asian soybean rust (ASR), caused by the fungus Phakopsora pachyrhizi , is the main disease affecting soybean in Brazil. This study aimed at investigating and mapping the resistance of the PI 594756 to P. pachyrhizi , by using Bulked Segregant Analysis (BSA). The PI 594756 and the susceptible PI 594891 were crossed and the resulting F 2 and F 2:3 populations (208 and 1770 plants, respectively) were tested against ASR. Also, these PIs and differential varieties were tested against a panel of monosporic isolates. Plants presenting tan lesions were classified as susceptible ( S ) while plants presenting reddish-brown (RB) lesions were classified as resistant. DNA bulks were genotyped with Infinium BeadChips and the genomic region identified was further analyzed in the F 2 individuals with target GBS (tGBS). PI 594,56 presented a unique resistance profile compared to the differential varieties. The resistance was monogenic dominant; however, it was classified as incompletely dominant when quantitatively studied. Genetic and QTL mapping placed the PI 594756 gene between the genomic region located at 55,863,741 and 56,123,516 bp of chromosome 18. This position is slightly upstream mapping positions of Rpp1 (PI 200492) and Rpp1-b (PI 594538A). Finally, we performed a haplotype analysis in a whole genomic sequencing-SNP database composed of Brazilian historical germplasm and sources of Rpp genes. We found SNPs that successfully differentiated the new PI 594756 allele from Rpp1 and Rpp1-b sources. The haplotype identified can be used as a tool for marker-assisted selection (MAS)., Supplementary Information: The online version contains supplementary material available at 10.1007/s11032-023-01358-4., Competing Interests: Conflict of interestThe authors declare no competing interests., (© The Author(s), under exclusive licence to Springer Nature B.V. 2023. Springer Nature or its licensor (e.g. a society or other partner) holds exclusive rights to this article under a publishing agreement with the author(s) or other rightsholder(s); author self-archiving of the accepted manuscript version of this article is solely governed by the terms of such publishing agreement and applicable law.)

Pratylenchus brachyurus causes serious damage to soybean production and other crops worldwide. Plant molecular responses to RLN infection remain largely unknown and no resistance genes have been identified in soybean. In this study, we analyzed molecular responses to RLN infection in moderately resistant BRSGO (Chapadões-BRS) and susceptible TMG115 RR (TMG) Glycine max genotypes. Differential expression analysis revealed two stages of response to RLN infection and a set of differentially expressed genes (DEGs) in the first stage suggested a pattern-triggered immunity (PTI) in both genotypes. The divergent time-point of DEGs between genotypes was observed four days post-infection, which included the activation of mitogen-activated protein kinase (MAPK) and plant-pathogen interaction genes in the BRS, suggesting the occurrence of an effector-triggered immunity response (ETI) in BRS. The co-expression analyses combined with single nucleotide polymorphism (SNP) uncovered a key element, a transcription factor phytochrome-interacting factor (PIF7) that is a potential regulator of moderate resistance to RLN infection. Two genes for resistance-related leucine-rich repeat (LRR) proteins were found as BRS-specific expressed genes. In addition, alternative splicing analysis revealed an intron retention in a myo-inositol oxygenase (MIOX) transcript, a gene related to susceptibility, may cause a loss of function in BRS., Competing Interests: The authors declare that they have no competing interests.

Asian soybean rust, caused by the fungus Phakopsora pachyrhizi , is one of the most important diseases affecting soybean production in tropical areas. During infection, P. pachyrhizi secretes proteins from haustoria that are transferred into plant cells to promote virulence. To date, only one candidate P. pachyrhizi effector protein has been characterized in detail to understand the mechanism by which it suppresses plant defenses to enhance infection. Here, we aimed to extend understanding of the pathogenic mechanisms of P. pachyrhizi based on the discovery of host proteins that interact with the effector candidate Phapa-7431740. We demonstrated that Phapa-7431740 suppresses pathogen-associated molecular pattern-triggered immunity (PTI) and that it interacts with a soybean glucan endo-1,3-β-glucosidase (GmβGLU), a pathogenesis-related (PR) protein belonging to the PR-2 family. Structural and phylogenetic characterization of the PR-2 protein family predicted in the soybean genome and comparison to PR-2 family members in Arabidopsis thaliana and cotton, demonstrated that GmβGLU is a type IV β-1,3-glucanase. Transcriptional profiling during an infection time course showed that the GmβGLU mRNA is highly induced during the initial hours after infection, coinciding with peak of expression of Phapa-7431740. The effector was able to interfere with the activity of GmβGLU in vitro, with a dose-dependent inhibition. Our results suggest that Phapa-7431740 may suppress PTI by interfering with glucan endo-1,3-β-glucosidase activity. [Formula: see text] The author(s) have dedicated the work to the public domain under the Creative Commons CC0 "No Rights Reserved" license by waiving all of his or her rights to the work worldwide under copyright law, including all related and neighboring rights, to the extent allowed by law, 2022.

Key Message: A locus on chromosome 13, containing multiple TIR-NB-LRR genes and SNPs associated with M. javanica resistance, was identified using a combination of GWAS, resequencing, genetic mapping and expression profiling. Meloidogyne javanica, a root-knot nematode, is an important problem in soybean-growing areas, leading to severe yield losses. Some accessions have been identified carrying resistance loci to this nematode. In this study, a set of 317 soybean accessions was characterized for resistance to M. javanica. A genome-wide association study was performed using SNPs from genotyping-by-sequencing, and a region of 29.2 kb on chromosome 13 was identified. An analysis of haplotypes showed that SNPs were able to discriminate between susceptible and resistant accessions, with 25 accessions sharing the haplotype associated with resistance. Furthermore, five accessions that exhibited resistance without carrying this haplotype may carry different loci conferring resistance to M. javanica. We also conducted the screening of the SNPs in the USDA soybean germplasm, revealing that several soybean accessions previously reported as resistant to other nematodes also shared the resistance haplotype on chromosome 13. Two SNP-based TaqMan® assays were developed and validated in two panels of soybean cultivars and in biparental populations. In silico analysis of the region associated with resistance identified the occurrence of genes with structural similarity with classical major resistance genes (NBS-LRR genes). Specifically, several nonsynonymous SNPs were observed in Glyma.13g194800 and Glyma.13g194900. The expression profile of these candidate genes demonstrated that the two gene models were up-regulated in the resistance source PI 505,099 after nematode infection. Overall, the SNPs associated with resistance and the genes identified constitute an important tool for introgression of resistance to the root-knot nematode by marker-assisted selection in soybean breeding programs.

Background: Southern stem canker (SSC), caused by Diaporthe aspalathi (E. Jansen, Castl. & Crous), is an important soybean disease that has been responsible for severe losses in the past. The main strategy for controlling this fungus involves the introgression of resistance genes. Thus far, five main loci have been associated with resistance to SSC. However, there is a lack of information about useful allelic variation at these loci. In this work, a genome-wide association study (GWAS) was performed to identify allelic variation associated with resistance against Diaporthe aspalathi and to provide molecular markers that will be useful in breeding programs., Results: We characterized the response to SSC infection in a panel of 295 accessions from different regions of the world, including important Brazilian elite cultivars. Using a GBS approach, the panel was genotyped, and we identified marker loci associated with Diaporthe aspalathi resistance through GWAS. We identified 19 SNPs associated with southern stem canker resistance, all on chromosome 14. The peak SNP showed an extremely high degree of association (p-value = 6.35E-27) and explained a large amount of the observed phenotypic variance (R 2  = 70%). This strongly suggests that a single major gene is responsible for resistance to D. aspalathi in most of the lines constituting this panel. In resequenced soybean materials, we identified other SNPs in the region identified through GWAS in the same LD block that clearly differentiate resistant and susceptible accessions. The peak SNP was selected and used to develop a cost-effective molecular marker assay, which was validated in a subset of the initial panel. In an accuracy test, this SNP assay demonstrated 98% selection efficiency., Conclusions: Our results suggest relevance of this locus to SSC resistance in soybean cultivars and accessions from different countries, and the SNP marker assay developed in this study can be directly applied in MAS studies in breeding programs to select materials that are resistant against this pathogen and support its introgression.

Asian soybean rust (ASR) is one of the most destructive diseases affecting soybeans. The causative agent of ASR, the fungus Phakopsora pachyrhizi, presents characteristics that make it difficult to study in vitro, limiting our knowledge of plant-pathogen dynamics. Therefore, this work used leaf lesion laser microdissection associated with deep sequencing to determine the pathogen transcriptome during compatible and incompatible interactions with soybean. The 36,350 generated unisequences provided an overview of the main genes and biological pathways that were active in the fungus during the infection cycle. We also identified the most expressed transcripts, including sequences similar to other fungal virulence and signaling proteins. Enriched P. pachyrhizi transcripts in the resistant (PI561356) soybean genotype were related to extracellular matrix organization and metabolic signaling pathways and, among infection structures, in amino acid metabolism and intracellular transport. Unisequences were further grouped into gene families along predicted sequences from 15 other fungi and oomycetes, including rust fungi, allowing the identification of conserved multigenic families, as well as being specific to P. pachyrhizi. The results revealed important biological processes observed in P. pachyrhizi, contributing with information related to fungal biology and, consequently, a better understanding of ASR.

Asian soybean rust (ASR), caused by the obligate biotrophic fungus Phakopsora pachyrhizi, can cause losses greater than 80%. Despite its economic importance, there is no soybean cultivar with durable ASR resistance. In addition, the P. pachyrhizi genome is not yet available. However, the availability of other rust genomes, as well as the development of sample enrichment strategies and bioinformatics tools, has improved our knowledge of the ASR secretome and its potential effectors. In this context, we used a combination of laser capture microdissection (LCM), RNAseq and a bioinformatics pipeline to identify a total of 36 350 P. pachyrhizi contigs expressed in planta and a predicted secretome of 851 proteins. Some of the predicted secreted proteins had characteristics of candidate effectors: small size, cysteine rich, do not contain PFAM domains (except those associated with pathogenicity) and strongly expressed in planta. A comparative analysis of the predicted secreted proteins present in Pucciniales species identified new members of soybean rust and new Pucciniales- or P. pachyrhizi-specific families (tribes). Members of some families were strongly up-regulated during early infection, starting with initial infection through haustorium formation. Effector candidates selected from two of these families were able to suppress immunity in transient assays, and were localized in the plant cytoplasm and nuclei. These experiments support our bioinformatics predictions and show that these families contain members that have functions consistent with P. pachyrhizi effectors., (© 2016 BSPP AND JOHN WILEY & SONS LTD.)

Background: Soybean [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., Results: We 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., Conclusions: The 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.

Asian soybean rust (ASR), caused by the obligate biotrophic fungus Phakopsora pachyrhizi, is one of most important diseases in the soybean (Glycine max (L.) Merr.) agribusiness. The identification and characterization of genes related to plant defense responses to fungal infection are essential to develop ASR-resistant plants. In this work, we describe four soybean genes, GmbZIP62, GmbZIP105, GmbZIPE1, and GmbZIPE2, which encode transcription factors containing a basic leucine zipper (bZIP) domain from two divergent classes, and that are responsive to P. pachyrhizi infection. Molecular phylogenetic analyses demonstrated that these genes encode proteins similar to bZIP factors responsive to pathogens. Yeast transactivation assays showed that only GmbZIP62 has strong transactivation activity in yeast. In addition, three of the bZIP transcription factors analyzed were also differentially expressed by plant defense hormones, and all were differentially expressed by fungal attack, indicating that these proteins might participate in response to ASR infection. The results suggested that these bZIP proteins are part of the plant defense response to P. pachyrhizi infection, by regulating the gene expression related to ASR infection responses. These bZIP genes are potential targets to obtain new soybean genotypes resistant to ASR.

Key Message: A 55% transformation efficiency was obtained by our optimized protocol; and we showed that GmELF1 - β and GmELF1 - α are the most stable reference genes for expression analyses under this specific condition. Gene functional analyses are essential to the validation of results obtained from in silico and/or gene-prospecting studies. Genetic transformation methods that yield tissues of transient expression quickly have been of considerable interest to researchers. Agrobacterium rhizogenes-mediated transformation methods, which are employed to generate plants with transformed roots, have proven useful for the study of stress caused by root phytopathogens via gene overexpression and/or silencing. While some protocols have been adapted to soybean plants, transformation efficiencies remain limited; thus, few viable plants are available for performing bioassays. Furthermore, mRNA analyses that employ reverse transcription quantitative polymerase chain reactions (RT-qPCR) require the use of reference genes with stable expression levels across different organs, development steps and treatments. In the present study, an A. rhizogenes-mediated soybean root transformation approach was optimized. The method delivers significantly higher transformation efficiency levels and rates of transformed plant recovery, thus enhancing studies of soybean abiotic conditions or interactions between phytopathogens, such as nematodes. A 55% transformation efficiency was obtained following the addition of an acclimation step that involves hydroponics and different selection processes. The present study also validated the reference genes GmELF1-β and GmELF1-α as the most stable to be used in RT-qPCR analysis in composite plants, mainly under nematode infection.

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., (Copyright © 2014 Elsevier Ireland Ltd. All rights reserved.)

Background: Many 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., Results: As 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., Conclusions: The 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.

Drought is a significant constraint to yield increase in soybean. The early perception of water deprivation is critical for recruitment of genes that promote plant tolerance. DeepSuperSAGE libraries, including one control and a bulk of six stress times imposed (from 25 to 150 min of root dehydration) for drought-tolerant and sensitive soybean accessions, allowed to identify new molecular targets for drought tolerance. The survey uncovered 120,770 unique transcripts expressed by the contrasting accessions. Of these, 57,610 aligned with known cDNA sequences, allowing the annotation of 32,373 unitags. A total of 1,127 unitags were up-regulated only in the tolerant accession, whereas 1,557 were up-regulated in both as compared to their controls. An expression profile concerning the most representative Gene Ontology (GO) categories for the tolerant accession revealed the expression "protein binding" as the most represented for "Molecular Function", whereas CDPK and CBL were the most up-regulated protein families in this category. Furthermore, particular genes expressed different isoforms according to the accession, showing the potential to operate in the distinction of physiological behaviors. Besides, heat maps comprising GO categories related to abiotic stress response and the unitags regulation observed in the expression contrasts covering tolerant and sensitive accessions, revealed the unitags potential for plant breeding. Candidate genes related to "hormone response" (LOX, ERF1b, XET), "water response" (PUB, BMY), "salt stress response" (WRKY, MYB) and "oxidative stress response" (PER) figured among the most promising molecular targets. Additionally, nine transcripts (HMGR, XET, WRKY20, RAP2-4, EREBP, NAC3, PER, GPX5 and BMY) validated by RT-qPCR (four different time points) confirmed their differential expression and pointed that already after 25 minutes a transcriptional reorganization started in response to the new condition, with important differences between both accessions.

Rpp4 (Resistance to Phakopsora pachyrhizi 4) confers resistance to Phakopsora pachyrhizi Sydow, the causal agent of Asian soybean rust (ASR). By combining expression profiling and virus induced gene silencing (VIGS), we are developing a genetic framework for Rpp4-mediated resistance. We measured gene expression in mock-inoculated and P. pachyrhizi-infected leaves of resistant soybean accession PI459025B (Rpp4) and the susceptible cultivar (Williams 82) across a 12-day time course. Unexpectedly, two biphasic responses were identified. In the incompatible reaction, genes induced at 12h after infection (hai) were not differentially expressed at 24 hai, but were induced at 72 hai. In contrast, genes repressed at 12 hai were not differentially expressed from 24 to 144 hai, but were repressed 216 hai and later. To differentiate between basal and resistance-gene (R-gene) mediated defence responses, we compared gene expression in Rpp4-silenced and empty vector-treated PI459025B plants 14 days after infection (dai) with P. pachyrhizi. This identified genes, including transcription factors, whose differential expression is dependent upon Rpp4. To identify differentially expressed genes conserved across multiple P. pachyrhizi resistance pathways, Rpp4 expression datasets were compared with microarray data previously generated for Rpp2 and Rpp3-mediated defence responses. Fourteen transcription factors common to all resistant and susceptible responses were identified, as well as fourteen transcription factors unique to R-gene-mediated resistance responses. These genes are targets for future P. pachyrhizi resistance research.

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.

Soybean [Glycine max (L.) Merril], one of the most important crop species in the world, is very susceptible to abiotic and biotic stress. Soybean plants have developed a variety of molecular mechanisms that help them survive stressful conditions. Hybrid proline-rich proteins (HyPRPs) constitute a family of cell-wall proteins with a variable N-terminal domain and conserved C-terminal domain that is phylogenetically related to non-specific lipid transfer proteins. Members of the HyPRP family are involved in basic cellular processes and their expression and activity are modulated by environmental factors. In this study, microarray analysis and real time RT-qPCR were used to identify putative HyPRP genes in the soybean genome and to assess their expression in different plant tissues. Some of the genes were also analyzed by time-course real time RT-qPCR in response to infection by Phakopsora pachyrhizi, the causal agent of Asian soybean rust disease. Our findings indicate that the time of induction of a defense pathway is crucial in triggering the soybean resistance response to P. pachyrhizi. This is the first study to identify the soybean HyPRP group B family and to analyze disease-responsive GmHyPRP during infection by P. pachyrhizi.

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.

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.

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.

Plants have the ability to recognize and respond to a multitude of pathogens, resulting in a massive reprogramming of the plant to activate defense responses including Resistance (R) and Pathogenesis-Related (PR) genes. Abiotic stresses can also activate PR genes and enhance pathogen resistance, representing valuable genes for breeding purposes. The present work offers an overview of soybean R and PR genes present in the GENOSOJA (Brazilian Soybean Genome Consortium) platform, regarding their structure, abundance, evolution and role in the plant-pathogen metabolic pathway, as compared with Medicago and Arabidopsis. Searches revealed 3,065 R candidates (756 in Soybean, 1,142 in Medicago and 1,167 in Arabidopsis), and PR candidates matching to 1,261 sequences (310, 585 and 366 for the three species, respectively). The identified transcripts were also evaluated regarding their expression pattern in 65 libraries, showing prevalence in seeds and developing tissues. Upon consulting the SuperSAGE libraries, 1,072 R and 481 PR tags were identified in association with the different libraries. Multiple alignments were generated for Xa21 and PR-2 genes, allowing inferences about their evolution. The results revealed interesting insights regarding the variability and complexity of defense genes in soybean, as compared with Medicago and Arabidopsis.

The soybean ubiquitous urease (encoded by GmEu4) is responsible for recycling metabolically derived urea. Additional biological roles have been demonstrated for plant ureases, notably in toxicity to other organisms. However, urease enzymatic activity is not related to its toxicity. The role of GmEu4 in soybean susceptibility to fungi was investigated in this study. A differential expression pattern of GmEu4 was observed in susceptible and resistant genotypes of soybeans over the course of a Phakopsora pachyrhizi infection, especially 24 h after infection. Twenty-nine adult, transgenic soybean plants, representing six independently transformed lines, were obtained. Although the initial aim of this study was to overexpress GmEu4, the transgenic plants exhibited GmEu4 co-suppression and decreased ureolytic activity. The growth of Rhizoctonia solani, Phomopsis sp., and Penicillium herguei in media containing a crude protein extract from either transgenic or non-transgenic leaves was evaluated. The fungal growth was higher in the protein extracts from transgenic urease-deprived plants than in extracts from non-transgenic controls. When infected by P. pachyrhizi uredospores, detached leaves of urease-deprived plants developed a significantly higher number of lesions, pustules and erupted pustules than leaves of non-transgenic plants containing normal levels of the enzyme. The results of the present work show that the soybean plants were more susceptible to fungi in the absence of urease. It was not possible to overexpress active GmEu4. For future work, overexpression of urease fungitoxic peptides could be attempted as an alternative approach.

We evaluated the molecular, anatomical and physiological properties of a soybean line transformed to improve drought tolerance with an rd29A:AtDREB1A construct. This construct expressed dehydration- responsive element binding protein DREB1A from the stress-inducible rd29A promoter. The greenhouse growth test included four randomized blocks of soybean plants, with each treatment performed in triplicate. Seeds from the non-transformed soybean cultivar BR16 and from the genetically modified soybean P58 line (T(2) generation) were grown at 15% gravimetric humidity for 31 days. To induce water deficit, the humidity was reduced to 5% gravimetric humidity (moderate stress) for 29 days and then to 2.5% gravimetric humidity (severe stress). AtDREB1A gene expression was higher in the genetically modified P58 plants during water deficit, demonstrating transgene stability in T(2) generations and induction of the rd29A promoter. Drought-response genes, including GmPI-PLC, GmSTP, GmGRP, and GmLEA14, were highly expressed in plants submitted to severe stress. Genetically modified plants had higher stomatal conductance and consequently higher photosynthetic and transpiration rates. In addition, they had more chlorophyll. Overexpression of AtDREB1A may contribute to a decrease in leaf thickness; however, a thicker abaxial epidermis was observed. Overexpression of AtDREB1A in soybean appears to enhance drought tolerance.

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.

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., (Copyright 2010 Elsevier Inc. All rights reserved.)

Soybean genotypes MG/BR46 (Conquista) and BR16, drought-tolerant and -sensitive, respectively, were compared in terms of morphophysiological and gene-expression responses to water stress during two stages of development. Gene-expression analysis showed differential responses in Gmdreb1a and Gmpip1b mRNA expression within 30 days of water-deficit initiation in MG/BR46 (Conquista) plants. Within 45 days of initiating stress, Gmp5cs and Gmpip1b had relatively higher expression. Initially, BR16 showed increased expression only for Gmdreb1a, and later (45 days) for Gmp5cs, Gmdefensin and Gmpip1b. Only BR16 presented down-regulated expression of genes, such as Gmp5cs and Gmpip1b, 30 days after the onset of moisture stress, and Gmgols after 45 days of stress. The faster perception of water stress in MG/BR46 (Conquista) and the better maintenance of up-regulated gene expression than in the sensitive BR16 genotype imply mechanisms by which the former is better adapted to tolerate moisture deficiency.

We determined the expression levels of DREB transcription factor (Gmdreb1) and of the genes Gmgols, Gmpip1b, Gmereb, and Gmdefensin in drought-tolerant (MG/BR46-Conquista) and drought-sensitive (BR16) genotypes of soybean, during drought. The trial was carried out in a controlled-environment chamber, set up to provide drought conditions. Sequences of Arabidopsis thaliana DREB-family proteins were used to build a phylogenetic tree through the alignment of the conserved regions near the AP2 domain. We found that Gmdreb1 is similar to Atrap2.1, which is located near the AtDREB1 and AtDREB2 families. The amplified fragment was cloned and sequenced; alignment with the sequence available at Genbank showed total similarity. Expression analysis showed that under drought: a) Gmdreb1 expression increased in leaves and roots of both genotypes and expression level changes occurred that were correlated with the length of the water-deficit period; b) there were increased expression levels of Gmdefensin in roots of MG/BR46; c) expression of Gmgols increased in leaves and roots of the two genotypes; d) Gmpip1b expression generally increased, except in roots of BR16, and e) the same was found for Gmereb, except in roots of MG/BR46.

Defensin, thionin and lipid transfer protein (LTP) gene families, which antimicrobial activity has an attractive use in protein engineering and transgenic production of agronomical important plants, have been here functionally reviewed. Also, a transcriptional overview of a set of plant SuperSAGE libraries and analysis looking for 26 bp tags possibly annotated for those families is presented. Tags differentially expressed (p = 0.05) or constitutively transcribed were identified from leaves or roots SuperSAGE libraries from important Brazilian plant species [cowpea (Vigna unguiculata (L.) Walp.), soybean (Glycine max (L.) Merr.) and modern sugarcane hybrids (Saccharum spp.)] submitted to abiotic [salt (100 mM NaCl) or drought] or biotic stresses [fungus inoculation (Phakopsora pachyrhizi; Asiatic Soyben Rust phytopathogen)]. The diverse transcriptional patterns observed, probably related to the variable range of targets and functions involved, could be the first step to unravel the antimicrobial peptide world and the plant stress response relationship. Moreover, SuperSAGE opens the opportunity to find some SNPs or even rare transcript that could be important on plant stress resistance mechanisms. Putative defensin or LTP identified by SuperSAGE following a specific plant treatment or physiological condition could be useful for future use in genetic improvement of plants.

During earlier experiments, an SSR molecular marker (176 Soy HSP) showing high correlation (70%) with resistance/susceptibility to javanese root-knot nematode Meloidogyne javanica was identified in soybean. After being sequenced, results indicated that the SSR 176 Soy HSP marker was inserted in the promoter region of Gmhsp17.6-L gene. It was also detected in this region that resistant genotypes presented insertions between AT(31) and AT(33) in size and susceptible genotypes, AT(9). Gmhsp17.6-L gene coding region presented a perfect match in amino acid sequence in all soybean genotypes. A ribonuclease protection assay showed that Gmhsp17.6-L gene mRNA transcripts were present in all genotypes. A real-time relative quantification (qPCR) indicated in the resistant individuals higher mRNA transcripts levels, which presented in the sequencing more AT(n) insertions. These results suggest that the number of AT(n) insertions inside this promoter region could modulate up or down gene levels. Those findings can lead to the possibility of manipulating, between some limits, the mRNA transcripts levels using different sizes of AT(n) insertions.

Asian soybean rust is a formidable threat to soybean (Glycine max) production in many areas of the world, including the United States. Only five sources of resistance have been identified (Resistance to Phakopsora pachyrhizi1 [Rpp1], Rpp2, Rpp3, Rpp4, and Rpp5). Rpp4 was previously identified in the resistant genotype PI459025B and mapped within 2 centimorgans of Satt288 on soybean chromosome 18 (linkage group G). Using simple sequence repeat markers, we developed a bacterial artificial chromosome contig for the Rpp4 locus in the susceptible cv Williams82 (Wm82). Sequencing within this region identified three Rpp4 candidate disease resistance genes (Rpp4C1-Rpp4C3 [Wm82]) with greatest similarity to the lettuce (Lactuca sativa) RGC2 family of coiled coil-nucleotide binding site-leucine rich repeat disease resistance genes. Constructs containing regions of the Wm82 Rpp4 candidate genes were used for virus-induced gene silencing experiments to silence resistance in PI459025B, confirming that orthologous genes confer resistance. Using primers developed from conserved sequences in the Wm82 Rpp4 candidate genes, we identified five Rpp4 candidate genes (Rpp4C1-Rpp4C5 [PI459025B]) from the resistant genotype. Additional markers developed from the Wm82 Rpp4 bacterial artificial chromosome contig further defined the region containing Rpp4 and eliminated Rpp4C1 (PI459025B) and Rpp4C3 (PI459025B) as candidate genes. Sequencing of reverse transcription-polymerase chain reaction products revealed that Rpp4C4 (PI459025B) was highly expressed in the resistant genotype, while expression of the other candidate genes was nearly undetectable. These data support Rpp4C4 (PI459025B) as the single candidate gene for Rpp4-mediated resistance to Asian soybean rust.

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.

Asian soybean rust (ASR), caused by Phakopsora pachyrhizi, is now established in all major soybean-producing countries. Currently, there is little information about the molecular basis of ASR-soybean interactions, which will be needed to assist future efforts to develop effective resistance. Toward this end, abundance changes of soybean mRNAs were measured over a 7-day ASR infection time course in mock-inoculated and infected leaves of a soybean accession (PI230970) carrying the Rpp2 resistance gene and a susceptible genotype (Embrapa-48). The expression profiles of differentially expressed genes (ASR-infected compared with the mock-inoculated control) revealed a biphasic response to ASR in each genotype. Within the first 12 h after inoculation (hai), which corresponds to fungal germination and penetration of the epidermal cells, differential gene expression changes were evident in both genotypes. mRNA expression of these genes mostly returned to levels found in mock-inoculated plants by 24 hai. In the susceptible genotype, gene expression remained unaffected by rust infection until 96 hai, a time period when rapid fungal growth began. In contrast, gene expression in the resistant genotype diverged from the mock-inoculated control earlier, at 72 h, demonstrating that Rpp2-mediated defenses were initiated prior to this time. These data suggest that ASR initially induces a nonspecific response that is transient or is suppressed when early steps in colonization are completed in both soybean genotypes. The race-specific resistance phenotype of Rpp2 is manifested in massive gene expression changes after the initial response prior to the onset of rapid fungal growth that occurs in the susceptible genotype.

Stability of the mitochondrial genome is controlled by nuclear loci. In plants, nuclear genes suppress mitochondrial DNA rearrangements during development. One nuclear gene involved in this process is Msh1. Msh1 appears to be involved in the suppression of illegitimate recombination in plant mitochondria. To test the hypothesis that Msh1 disruption leads to the type of mitochondrial DNA rearrangements associated with naturally occurring cytoplasmic male sterility in plants, a transgenic approach for RNAi was used to modulate expression of Msh1 in tobacco and tomato. In both species, these experiments resulted in reproducible mitochondrial DNA rearrangements and a condition of male (pollen) sterility. The male sterility was, in each case, heritable, associated with normal female fertility, and apparently maternal in its inheritance. Segregation of the transgene did not reverse the male sterile phenotype, producing stable, nontransgenic male sterility. The reproducible transgenic induction of mitochondrial rearrangements in plants is unprecedented, providing a means to develop novel cytoplasmic male sterile lines for release as non-GMO or transgenic materials.

Mitochondrial processes influence a broad spectrum of physiological and developmental events in higher eukaryotes, and their aberrant function can lead to several familiar disease phenotypes in mammals. In plants, mitochondrial genes directly influence pollen development and the occurrence of male sterility in natural plant populations. Likewise, in animal systems evidence accumulates to suggest important mitochondrial functions in spermatogenesis and reproduction. Here we present evidence for a convergent gene fusion involving a MutS-homologous gene functioning within the mitochondrion and designated Msh1. In only plants and soft corals, the MutS homologue has fused with a homing endonuclease sequence at the carboxy terminus of the protein. However, the endonuclease domains in the plants and the soft corals are members of different groups. In plants, Msh1 can influence mitochondrial genome organization and male sterility expression. Based on parallels in Msh1 gene structure shared by plants and corals, and their similarities in reproductive behavior, we postulate that this convergent gene fusion might have occurred in response to coincident adaptive pressures on reproduction.

The plant mitochondrial genome is retained in a multipartite structure that arises by a process of repeat-mediated homologous recombination. Low-frequency ectopic recombination also occurs, often producing sequence chimeras, aberrant ORFs, and novel subgenomic DNA molecules. This genomic plasticity may distinguish the plant mitochondrion from mammalian and fungal types. In plants, relative copy number of recombination-derived subgenomic DNA molecules within mitochondria is controlled by nuclear genes, and a genomic shifting process can result in their differential copy number suppression to nearly undetectable levels. We have cloned a nuclear gene that regulates mitochondrial substoichiometric shifting in Arabidopsis. The CHM gene was shown to encode a protein related to the MutS protein of Escherichia coli that is involved in mismatch repair and DNA recombination. We postulate that the process of substoichiometric shifting in plants may be a consequence of ectopic recombination suppression or replication stalling at ectopic recombination sites to effect molecule-specific copy number modulation.