Soybean-Phakopsora pachyrhizi interactions: towards the development of next-generation disease-resistant plants.

Soybean rust (SBR), caused by the obligate biotrophic fungus Phakopsora pachyrhizi, is a devastating foliar disease threatening soybean production. To date, no commercial cultivars conferring durable resistance to SBR are available. The development of long-lasting SBR resistance has been hindered by the lack of understanding of this complex pathosystem, encompassing challenges posed by intricate genetic structures in both the host and pathogen, leading to a gap in the knowledge of gene-for-gene interactions between soybean and P. pachyrhizi. In this review, we focus on recent advancements and emerging technologies that can be used to improve our understanding of the P. pachyrhizi-soybean molecular interactions. We further explore approaches used to combat SBR, including conventional breeding, transgenic approaches and RNA interference, and how advances in our understanding of plant immune networks, the availability of new molecular tools, and the recent sequencing of the P. pachyrhizi genome could be used to aid in the development of better genetic resistance against SBR. Lastly, we discuss the research gaps of this pathosystem and how new technologies can be used to shed light on these questions and to develop durable next-generation SBR-resistant soybean plants.

Mapping of a soybean rust resistance in PI 594756 at the Rpp1 locus.

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 F2 and F2: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 F2 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.

Major proliferation of transposable elements shaped the genome of the soybean rust pathogen Phakopsora pachyrhizi.

With >7000 species the order of rust fungi has a disproportionately large impact on agriculture, horticulture, forestry and foreign ecosystems. The infectious spores are typically dikaryotic, a feature unique to fungi in which two haploid nuclei reside in the same cell. A key example is Phakopsora pachyrhizi, the causal agent of Asian soybean rust disease, one of the world's most economically damaging agricultural diseases. Despite P. pachyrhizi's impact, the exceptional size and complexity of its genome prevented generation of an accurate genome assembly. Here, we sequence three independent P. pachyrhizi genomes and uncover a genome up to 1.25 Gb comprising two haplotypes with a transposable element (TE) content of ~93%. We study the incursion and dominant impact of these TEs on the genome and show how they have a key impact on various processes such as host range adaptation, stress responses and genetic plasticity.

Overexpression of the GmEXPA1 gene reduces plant susceptibility to Meloidogyne incognita.

KEY MESSAGE: The overexpression of the soybean GmEXPA1 gene reduces plant susceptibility to M. incognita by the increase of root lignification. Plant expansins are enzymes that act in a pH-dependent manner in the plant cell wall loosening and are associated with improved tolerance or resistance to abiotic or biotic stresses. Plant-parasitic nematodes (PPN) can alter the expression profile of several expansin genes in infected root cells. Studies have shown that overexpression or downregulation of particular expansin genes can reduce plant susceptibility to PPNs. Root-knot nematodes (RKN) are obligate sedentary endoparasites of the genus Meloidogyne spp. of which M. incognita is one of the most reported species. Herein, using a transcriptome dataset and real-time PCR assays were identified an expansin A gene (GmEXPA1; Glyma.02G109100) that is upregulated in the soybean nematode-resistant genotype PI595099 compared to the susceptible cultivar BRS133 during plant parasitism by M. incognita. To understand the role of the GmEXPA1 gene during the interaction between soybean plant and M. incognita were generated stable A. thaliana and N. tabacum transgenic lines. Remarkably, both A. thaliana and N. tabacum transgenic lines overexpressing the GmEXPA1 gene showed reduced susceptibility to M. incognita. Furthermore, plant growth, biomass accumulation, and seed yield were not affected in these transgenic lines. Interestingly, significant upregulation of the NtACC oxidase and NtEFE26 genes, involved in ethylene biosynthesis, and NtCCR and Nt4CL genes, involved in lignin biosynthesis, was observed in roots of the N. tabacum transgenic lines, which also showed higher lignin content. These data suggested a possible link between GmEXPA1 gene expression and increased lignification of the root cell wall. Therefore, these data support that engineering of the GmEXPA1 gene in soybean offers a powerful biotechnology tool to assist in RKN management.

Time Course RNA-seq Reveals Soybean Responses against Root-Lesion Nematode and Resistance Players.

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.

Overexpression of a soybean Globin (GmGlb1-1) gene reduces plant susceptibility to Meloidogyne incognita.

MAIN CONCLUSION: The overexpression of the GmGlb1-1 gene reduces plant susceptibility to Meloidogyne incognita. Non-symbiotic globin class #1 (Glb1) genes are expressed in different plant organs, have a high affinity for oxygen, and are related to nitric oxide (NO) turnover. Previous studies showed that soybean Glb1 genes are upregulated in soybean plants under flooding conditions. Herein, the GmGlb1-1 gene was identified in soybean as being upregulated in the nematode-resistant genotype PI595099 compared to the nematode-susceptible cultivar BRS133 during plant parasitism by Meloidogyne incognita. The Arabidopsis thaliana and Nicotiana tabacum transgenic lines overexpressing the GmGlb1-1 gene showed reduced susceptibility to M. incognita. Consistently, gall morphology data indicated that pJ2 nematodes that infected the transgenic lines showed developmental alterations and delayed parasitism progress. Although no significant changes in biomass and seed yield were detected, the transgenic lines showed an elongated, etiolation-like growth under well-irrigation, and also developed more axillary roots under flooding conditions. In addition, transgenic lines showed upregulation of some important genes involved in plant defense response to oxidative stress. In agreement, higher hydrogen peroxide accumulation and reduced activity of reactive oxygen species (ROS) detoxification enzymes were also observed in these transgenic lines. Thus, based on our data and previous studies, it was hypothesized that constitutive overexpression of the GmGlb1-1 gene can interfere in the dynamics of ROS production and NO scavenging, enhancing the acquired systemic acclimation to biotic and abiotic stresses, and improving the cellular homeostasis. Therefore, these collective data suggest that ectopic or nematode-induced overexpression, or enhanced expression of the GmGlb1-1 gene using CRISPR/dCas9 offers great potential for application in commercial soybean cultivars aiming to reduce plant susceptibility to M. incognita.

A Phakopsora pachyrhizi Effector Suppresses PAMP-Triggered Immunity and Interacts with a Soybean Glucan Endo-1,3-ß-Glucosidase to Promote Virulence.

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.

Mapping Major Disease Resistance Genes in Soybean by Genome-Wide Association Studies.

Soybean is one of the most valuable agricultural crops in the world. Besides, this legume is constantly attacked by a wide range of pathogens (fungi, bacteria, viruses, and nematodes) compromising yield and increasing production costs. One of the major disease management strategies is the genetic resistance provided by single genes and quantitative trait loci (QTL). Identifying the genomic regions underlying the resistance against these pathogens on soybean is one of the first steps performed by molecular breeders. In the past, genetic mapping studies have been widely used to discover these genomic regions. However, over the last decade, advances in next-generation sequencing technologies and their subsequent cost decreasing led to the development of cost-effective approaches to high-throughput genotyping. Thus, genome-wide association studies applying thousands of SNPs in large sets composed of diverse soybean accessions have been successfully done. In this chapter, a comprehensive review of the majority of GWAS for soybean diseases published since this approach was developed is provided. Important diseases caused by Heterodera glycines, Phytophthora sojae, and Sclerotinia sclerotiorum have been the focus of the several GWAS. However, other bacterial and fungi diseases also have been targets of GWAS. As such, this GWAS summary can serve as a guide for future studies of these diseases. The protocol begins by describing several considerations about the pathogens and bringing different procedures of molecular characterization of them. Advice to choose the best isolate/race to maximize the discovery of multiple R genes or to directly map an effective R gene is provided. A summary of protocols, methods, and tools to phenotyping the soybean panel is given to several diseases. We also give details of options of DNA extraction protocols and genotyping methods, and we describe parameters of SNP quality to soybean data. Websites and their online tools to obtain genotypic and phenotypic data for thousands of soybean accessions are highlighted. Finally, we report several tricks and tips in Subheading 4, especially related to composing the soybean panel as well as generating and analyzing the phenotype data. We hope this protocol will be helpful to achieve GWAS success in identifying resistance genes on soybean.

A Target-Site Mutation Confers Cross-Resistance to ALS-Inhibiting Herbicides in Erigeron sumatrensis from Brazil.

Cases of weed resistant to herbicides have changed the dynamics of agricultural areas in Brazil, and in recent years, Erigeron species have caused major problems to farmers in the country, mainly in relation to the ineffectiveness of herbicide treatments used. The objective of this study was to confirm the cross-resistance to ALS inhibitors in populations of Erigeron sumatrensis as well as to investigate the existence of mutations in the site of action of ALS-inhibiting herbicides. To do this, 30 populations collected in the 2016/2017 crop season were grown in a greenhouse. Dose-response (chlorimuron-ethyl and cloransulam-methyl), inhibition of cytochrome P-450 with malathion, and ALS gene sequencing experiments were carried out in the F1 generations of two fleabane populations. The results proved the cross-resistance to chlorimuron-ethyl and cloransulam-methyl herbicides applied in the post-emergence of the resistant population of E. sumatrensis. The higher activity of P450 enzymes is unlikely responsible for the resistance of the population studied. The resistance mechanism found in R was the target site mutation Pro197Ser at the ALS gene. This is the first study in Brazil to identify a target-site change as a survival mechanism in E. sumatrensis for the resistance to ALS-inhibiting herbicides.

Secondary Metabolites of Pseudomonas aeruginosa LV Strain Decrease Asian Soybean Rust Severity in Experimentally Infected Plants.

Asian Soybean Rust (ASR), a disease caused by Phakopsora pachyrhizi, causing yield losses up to 90%. The control is based on the fungicides which may generate resistant fungi. The activation of the plant defense system, should help on ASR control. In this study, secondary metabolites of Pseudomonas aeruginosa LV strain were applied on spore germination and the expression of defense genes in infected soybean plants. The F4A fraction and the pure metabolites were used. In vitro, 10 µg mL-1 of F4A reduced spore germination by 54%, while 100 µg mL-1 completely inhibited. Overexpression of phenylalanine ammonia lyase (PAL), O-methyltransferase (OMT) and pathogenesis related protein-2 (PR-2; glucanases) defense-related genes were detected 24 and 72 h after soybean sprouts were sprayed with an organocopper antimicrobial compound (OAC). Under greenhouse conditions, the best control was observed in plants treated with 60 µg mL-1 of PCA, which reduced ASR severity and lesion frequency by 75% and 43%, respectively. Plants sprayed with 2 and 20 µg mL-1 of F4A also decreased severity (41%) and lesion frequency (32%). The significant reduction in spore germination ASR in plant suggested that the strain of these metabolites are effective against P. pachyrhizi, and they can be used for ASR control.

Soybean Metabolomics Based in Mass Spectrometry: Decoding the Plant's Signaling and Defense Responses under Biotic Stress.

Metabolomics is an omics technology that is extremely valuable to analyze all small-molecule metabolites in organisms. Recent advances in analytical instrumentation, such as mass spectrometry combined with data processing tools, chemometrics, and spectral data libraries, allow plant metabolomics studies to play a fundamental role in the agriculture field and food security. Few studies are found in the literature using the metabolomics approach in soybean plants on biotic stress. In this review, we provide a new perspective highlighting the potential of metabolomics-based mass spectrometry for soybean in response to biotic stress. Furthermore, we highlight the response and adaptation mechanisms of soybean on biotic stress about primary and secondary metabolism. Consequently, we provide subsidies for further studies of the resistance and improvement of the crop.

Untargeted Metabolomics Analysis by UHPLC-MS/MS of Soybean Plant in a Compatible Response to Phakopsora pachyrhizi Infection.

Phakopsora pachyrhizi is a biotrophic fungus, causer of the disease Asian Soybean Rust, a severe crop disease of soybean and one that demands greater investment from producers. Thus, research efforts to control this disease are still needed. We investigated the expression of metabolites in soybean plants presenting a resistant genotype inoculated with P. pachyrhizi through the untargeted metabolomics approach. The analysis was performed in control and inoculated plants with P. pachyrhizi using UHPLC-MS/MS. Principal component analysis (PCA) and the partial least squares discriminant analysis (PLS-DA), was applied to the data analysis. PCA and PLS-DA resulted in a clear separation and classification of groups between control and inoculated plants. The metabolites were putative classified and identified using the Global Natural Products Social Molecular Networking platform in flavonoids, isoflavonoids, lipids, fatty acyls, terpenes, and carboxylic acids. Flavonoids and isoflavonoids were up-regulation, while terpenes were down-regulated in response to the soybean-P. pachyrhizi interaction. Our data provide insights into the potential role of some metabolites as flavonoids and isoflavonoids in the plant resistance to ASR. This information could result in the development of resistant genotypes of soybean to P. pachyrhizi, and effective and specific products against the pathogen.

Genome-wide association study for resistance to the Meloidogyne javanica causing root-knot nematode in soybean.

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.

The soybean gene GmHsp22.4 is involved in the resistance response to Meloidogyne javanica in Arabidopsis thaliana.

BACKGROUND: Small heat shock proteins (sHSPs) belong to the class of molecular chaperones that respond to biotic and abiotic stresses in plants. A previous study has showed strong induction of the gene GmHsp22.4 in response to the nematode Meloidogyne javanica in a resistant soybean genotype, while repression in a susceptible one. This study aimed to investigate the functional involvement of this small chaperone in response to M. javanica in Arabidopsis thaliana. First, it was evaluated the activation of the promoter region after the nematode inoculation, and the occurrence of polymorphisms between resistant and susceptible re-sequenced soybean accessions. Then functional analysis using A. thaliana lines overexpressing the soybean GmHsp22.4 gene, and knocked-out mutants were challenged with M. javanica infestation. RESULTS: High expression levels of the GFP gene marker in transformed A. thaliana plants revealed that the promoter region of GmHsp22.4 was strongly activated after nematode inoculation. Moreover, the multiplication of the nematode was significantly reduced in plants overexpressing GmHsp22.4 gene in A. thaliana compared to the wild type. Additionally, the multiplication of M. javanica in the A. thaliana mutants was significantly increased mainly in the event athsp22.0-2. This increase was not that evident in the event athsp22.0-1, the one that preserved a portion of the promoter region, including the HSEs in the region around - 83 bp. However, structural analysis at sequence level among soybean resistant and susceptible genotypes did not detect any polymorphisms in the whole gene model. CONCLUSIONS: The soybean chaperone GmHsp22.4 is involved in the defense response to root-knot nematode M. javanica in A. thaliana. Specifically, the promoter region covering until - 191 from the transcriptional start site (TSS) is necessary to promoter activation after nematode infection in Arabidopsis. No polymorphisms that could explain these differences in the defense response were detected in the GmHsp22.4 gene between resistant and susceptible soybean genotypes. Therefore, further investigation is needed to elucidate the triggering factor of the plant's defense mechanism, both at the sequence level of the soybean genotypes presenting contrasting reaction to root-knot nematode and by detecting cis-elements that are essential for the activation of the GmHsp22.4 gene promoter.

Transcriptional profile of genes involved in the production of terpenes and glyceollins in response to biotic stresses in soybean.

Terpenes produced by plants comprise a diverse range of secondary metabolites, including volatile organic compounds (VOCs). Terpene VOC production may be altered after damage or by biological stimuli such as bacterial, fungal and insects, and subsequent triggering of plant defense responses. These VOCs originate in plants from two independent pathways: the mevalonate and the methylerythritol phosphate pathways, which utilize dimethylallyl and isopentenyl diphosphates to form the terpenoidal precursors. Phakopsora pachyrhizi fungi causes Asian soybean rust, limiting soybean production and resulting in losses of up to 80% if no control strategies are applied. By using a transcriptome datasets, we investigated the regulation of genes of the mevalonate pathway under different biotic stresses. We studied the impact of P. pachyrhizi infection in vivo expression profile of genes involved in terpenoid and glyceollin biosynthesis in genotypes harboring different resistance genes (Rpp), and across the infection cycle. In addition, we used UPLC and UPGC analysis to evaluate glyceollin and VOC production, respectively, to identify metabolites associated with soybean responses to pathogen infection. The regulation of soybean genes involved in terpene production was influenced by genotypes, depending on the Rpp gene, while glyceollin was induced in all genotypes. Furthermore, a sesquiterpene was identified as a potential marker associated with rust symptoms on soybean.

Phakopsora pachyrhizi triggers the jasmonate signaling pathway during compatible interaction in soybean and GmbZIP89 plays a role of major component in the pathway.

The biotrophic fungus Phakopsora pachyrhizi is currently the major pathogen affecting soybean production worldwide. It has already been suggested for the non-host interaction between P. pachyrhizi and Arabidopsis thaliana that the fungus in early infection induces jasmonic acid (JA) pathway to the detriment of the salicylic acid (SA) pathway as a mechanism to the establishment of infection. In this study, we verified that this mechanism might also be occurring during the compatible interaction in soybean (Glycine max L. Merril). It was demonstrated that P. pachyrhizi triggers a JA pathway during the early and late stages of infection in a susceptible soybean cultivar. The expression of the GmbZIP89 was induced in a biphasic profile, similarly to other JA responsive genes, which indicates a new marker gene for this signaling pathway. Additionally, plants silenced for GmbZIP89 (iGmZIP89) by the virus-induced gene silencing (VIGS) approach present lower severity of infection and higher expression of pathogenesis related protein 1 (PR1). The lower disease severity showed that the iGmbZIP89 plants became more resistant to infection. These data corroborate the hypothesis that the GmbZIP89 may be a resistance negative regulator. In conclusion, we demonstrated that P. pachyrhizi mimics a necrotrophic fungus and activates the JA/ET pathway in soybean. It is possible to suppose that its direct penetration on epidermal cells or fungal effectors may modulate the expression of target genes aiming the activation of the JA pathway and inhibition of SA defense.

Association mapping of a locus that confers southern stem canker resistance in soybean and SNP marker development.

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 (R2 = 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.

New insights into Phakopsora pachyrhizi infection based on transcriptome analysis in planta.

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.

Prediction of the in planta Phakopsora pachyrhizi secretome and potential effector families.

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.

Evaluation of genetic variation among Brazilian soybean cultivars through genome resequencing.

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.