Deconstructing the genetic architecture of iron deficiency chlorosis in soybean using genome-wide approaches.

BACKGROUND: Iron (Fe) is an essential micronutrient for plant growth and development. Iron deficiency chlorosis (IDC), caused by calcareous soils or high soil pH, can limit iron availability, negatively affecting soybean (Glycine max) yield. This study leverages genome-wide association study (GWAS) and a genome-wide epistatic study (GWES) with previous gene expression studies to identify regions of the soybean genome important in iron deficiency tolerance. RESULTS: A GWAS and a GWES were performed using 460 diverse soybean PI lines from 27 countries, in field and hydroponic iron stress conditions, using more than 36,000 single nucleotide polymorphism (SNP) markers. Combining this approach with available RNA-sequencing data identified significant markers, genomic regions, and novel genes associated with or responding to iron deficiency. Sixty-nine genomic regions associated with IDC tolerance were identified across 19 chromosomes via the GWAS, including the major-effect quantitative trait locus (QTL) on chromosome Gm03. Cluster analysis of significant SNPs in this region deconstructed this historically prominent QTL into four distinct linkage blocks, enabling the identification of multiple candidate genes for iron chlorosis tolerance. The complementary GWES identified SNPs in this region interacting with nine other genomic regions, providing the first evidence of epistatic interactions impacting iron deficiency tolerance. CONCLUSIONS: This study demonstrates that integrating cutting edge genome wide association (GWA), genome wide epistasis (GWE), and gene expression studies is a powerful strategy to identify novel iron tolerance QTL and candidate loci from diverse germplasm. Crops, unlike model species, have undergone selection for thousands of years, constraining and/or enhancing stress responses. Leveraging genomics-enabled approaches to study these adaptations is essential for future crop improvement.

Glycerol-3-phosphate mediates rhizobia-induced systemic signaling in soybean.

Glycerol-3-phosphate (G3P) is a well-known mobile regulator of systemic acquired resistance (SAR), which provides broad spectrum systemic immunity in response to localized foliar pathogenic infections. We show that G3P-derived foliar immunity is also activated in response to genetically-regulated incompatible interactions with nitrogen-fixing bacteria. Using gene knock-down we show that G3P is essential for strain-specific exclusion of non-desirable root-nodulating bacteria and the associated foliar pathogen immunity in soybean. Grafting studies show that while recognition of rhizobium incompatibility is root driven, bacterial exclusion requires G3P biosynthesis in the shoot. Biochemical analyses support shoot-to-root transport of G3P during incompatible rhizobia interaction. We describe a root-shoot-root signaling mechanism which simultaneously enables the plant to exclude non-desirable nitrogen-fixing rhizobia in the root and pathogenic microbes in the shoot.

Transcriptome divergence during leaf development in two contrasting switchgrass (Panicum virgatum L.) cultivars.

The genetics and responses to biotic stressors of tetraploid switchgrass (Panicum virgatum L.) lowland cultivar 'Kanlow' and upland cultivar Summer are distinct and can be exploited for trait improvement. In general, there is a paucity of data on the basal differences in transcription across tissue developmental times for switchgrass cultivars. Here, the changes in basal and temporal expression of genes related to leaf functions were evaluated for greenhouse grown 'Kanlow', and 'Summer' plants. Three biological replicates of the 4th leaf pooled from 15 plants per replicate were harvested at regular intervals beginning from leaf emergence through senescence. Increases and decreases in leaf chlorophyll and N content were similar for both cultivars. Likewise, multidimensional scaling (MDS) analysis indicated both cultivar-independent and cultivar-specific gene expression. Cultivar-independent genes and gene-networks included those associated with leaf function, such as growth/senescence, carbon/nitrogen assimilation, photosynthesis, chlorophyll biosynthesis, and chlorophyll degradation. However, many genes encoding nucleotide-binding leucine rich repeat (NB-LRRs) proteins and wall-bound kinases associated with detecting and responding to environmental signals were differentially expressed. Several of these belonged to unique cultivar-specific gene co-expression networks. Analysis of genomic resequencing data provided several examples of NB-LRRs genes that were not expressed and/or apparently absent in the genomes of Summer plants. It is plausible that cultivar (ecotype)-specific genes and gene-networks could be one of the drivers for the documented differences in responses to leaf-borne pathogens between these two cultivars. Incorporating broad resistance to plant pathogens in elite switchgrass germplasm could improve sustainability of biomass production under low-input conditions.

A Two-Amino Acid Difference in the Coat Protein of Satellite panicum mosaic virus Isolates Is Responsible for Differential Synergistic Interactions with Panicum mosaic virus.

Panicum mosaic virus (PMV) (genus Panicovirus, family Tombusviridae) and its molecular parasite, Satellite panicum mosaic virus (SPMV), synergistically interact in coinfected proso and pearl millet (Panicum miliaceum L.) plants resulting in a severe symptom phenotype. In this study, we examined synergistic interactions between the isolates of PMV and SPMV by using PMV-NE, PMV85, SPMV-KS, and SPMV-Type as interacting partner viruses in different combinations. Coinfection of proso millet plants by PMV-NE and SPMV-KS elicited severe mosaic, chlorosis, stunting, and eventual plant death compared with moderate mosaic, chlorotic streaks, and stunting by PMV85 and SPMV-Type. In reciprocal combinations, coinfection of proso millet by either isolate of PMV with SPMV-KS but not with SPMV-Type elicited severe disease synergism, suggesting that SPMV-KS was the main contributor for efficient synergistic interaction with PMV isolates. Coinfection of proso millet plants by either isolate of PMV and SPMV-KS or SPMV-Type caused increased accumulation of coat protein (CP) and genomic RNA copies of PMV, compared with infections by individual PMV isolates. Additionally, CP and genomic RNA copies of SPMV-KS accumulated at substantially higher levels, compared with SMPV-Type in coinfected proso millet plants with either isolate of PMV. Hybrid viruses between SPMV-KS and SPMV-Type revealed that SPMV isolates harboring a CP fragment with four differing amino acids at positions 18, 35, 59, and 98 were responsible for differential synergistic interactions with PMV in proso millet plants. Mutation of amino acid residues at these positions in different combinations in SPMV-KS, similar to those as in SPMV-Type or vice-versa, revealed that A35 and R98 in SPMV-KS CP play critical roles in enhanced synergistic interactions with PMV isolates. Taken together, these data suggest that the two distinct amino acids at positions 35 and 98 in the CP of SPMV-KS and SPMV-Type are involved in the differential synergistic interactions with the helper viruses.

Genetic Architecture of Charcoal Rot (Macrophomina phaseolina) Resistance in Soybean Revealed Using a Diverse Panel.

Charcoal rot (CR) disease caused by Macrophomina phaseolina is responsible for significant yield losses in soybean production. Among the methods available for controlling this disease, breeding for resistance is the most promising. Progress in breeding efforts has been slow due to the insufficient information available on the genetic mechanisms related to resistance. Genome-wide association studies (GWAS) enable unraveling the genetic architecture of resistance and identification of causal genes. The aims of this study were to identify new sources of resistance to CR in a collection of 459 diverse plant introductions from the USDA Soybean Germplasm Core Collection using field and greenhouse screenings, and to conduct GWAS to identify candidate genes and associated molecular markers. New sources for CR resistance were identified from both field and greenhouse screening from maturity groups I, II, and III. Five significant single nucleotide polymorphism (SNP) and putative candidate genes related to abiotic and biotic stress responses are reported from the field screening; while greenhouse screening revealed eight loci associated with eight candidate gene families, all associated with functions controlling plant defense response. No overlap of markers or genes was observed between field and greenhouse screenings suggesting a complex molecular mechanism underlying resistance to CR in soybean with varied response to different environments; but our findings provide useful information for advancing breeding for CR resistance as well as the genetic mechanism of resistance.

Main and epistatic loci studies in soybean for Sclerotinia sclerotiorum resistance reveal multiple modes of resistance in multi-environments.

Genome-wide association (GWAS) and epistatic (GWES) studies along with expression studies in soybean [Glycine max (L.) Merr.] were leveraged to dissect the genetics of Sclerotinia stem rot (SSR) [caused by Sclerotinia sclerotiorum (Lib.) de Bary], a significant fungal disease causing yield and quality losses. A large association panel of 466 diverse plant introduction accessions were phenotyped in multiple field and controlled environments to: (1) discover sources of resistance, (2) identify SNPs associated with resistance, and (3) determine putative candidate genes to elucidate the mode of resistance. We report 58 significant main effect loci and 24 significant epistatic interactions associated with SSR resistance, with candidate genes involved in a wide range of processes including cell wall structure, hormone signaling, and sugar allocation related to plant immunity, revealing the complex nature of SSR resistance. Putative candidate genes [for example, PHYTOALEXIN DEFFICIENT 4 (PAD4), ETHYLENE-INSENSITIVE 3-LIKE 1 (EIL3), and ETHYLENE RESPONSE FACTOR 1 (ERF1)] clustered into salicylic acid (SA), jasmonic acid (JA), and ethylene (ET) pathways suggest the involvement of a complex hormonal network typically activated by both necrotrophic (ET/JA) and biotrophic (SA) pathogens supporting that S. sclerotiorum is a hemibiotrophic plant pathogen.

Virus-Induced Gene Silencing and Transient Gene Expression in Soybean (Glycine max) Using Bean Pod Mottle Virus Infectious Clones.

Virus-induced gene silencing (VIGS) is a powerful and rapid approach for determining the functions of plant genes. The basis of VIGS is that a viral genome is engineered so that it can carry fragments of plant genes, typically in the 200 to 300 base pair size range. The recombinant viruses are used to infect experimental plants, and wherever the virus invades, the target gene or genes will be silenced. VIGS is thus transient, and in the span of a few weeks, it is possible to design VIGS constructs and then generate loss-of-function phenotypes through RNA silencing of the target genes. In soybean (Glycine max), the Bean pod mottle virus (BPMV) has been engineered to be valuable tool for silencing genes with diverse functions and also for over-expression of foreign genes. This protocol describes a method for designing BPMV constructs and using them to silence or transiently express genes in soybean. © 2016 by John Wiley & Sons, Inc.

Simultaneous mutations in multi-viral proteins are required for soybean mosaic virus to gain virulence on soybean genotypes carrying different R genes.

BACKGROUND: Genetic resistance is the most effective and sustainable approach to the control of plant pathogens that are a major constraint to agriculture worldwide. In soybean, three dominant R genes, i.e., Rsv1, Rsv3 and Rsv4, have been identified and deployed against Soybean mosaic virus (SMV) with strain-specificities. Molecular identification of virulent determinants of SMV on these resistance genes will provide essential information for the proper utilization of these resistance genes to protect soybean against SMV, and advance knowledge of virus-host interactions in general. METHODOLOGY/PRINCIPAL FINDINGS: To study the gain and loss of SMV virulence on all the three resistance loci, SMV strains G7 and two G2 isolates L and LRB were used as parental viruses. SMV chimeras and mutants were created by partial genome swapping and point mutagenesis and then assessed for virulence on soybean cultivars PI96983 (Rsv1), L-29 (Rsv3), V94-5152 (Rsv4) and Williams 82 (rsv). It was found that P3 played an essential role in virulence determination on all three resistance loci and CI was required for virulence on Rsv1- and Rsv3-genotype soybeans. In addition, essential mutations in HC-Pro were also required for the gain of virulence on Rsv1-genotype soybean. To our best knowledge, this is the first report that CI and P3 are involved in virulence on Rsv1- and Rsv3-mediated resistance, respectively. CONCLUSIONS/SIGNIFICANCE: Multiple viral proteins, i.e., HC-Pro, P3 and CI, are involved in virulence on the three resistance loci and simultaneous mutations at essential positions of different viral proteins are required for an avirulent SMV strain to gain virulence on all three resistance loci. The likelihood of such mutations occurring naturally and concurrently on multiple viral proteins is low. Thus, incorporation of all three resistance genes in a soybean cultivar through gene pyramiding may provide durable resistance to SMV.

Mutations in the P3 protein of Soybean mosaic virus G2 isolates determine virulence on Rsv4-genotype soybean.

Two Soybean mosaic virus (SMV) G2 isolates, L and L-RB, sharing high-sequence similarly but differing in ability to break Rsv4-mediated resistance in soybean, were investigated. Infectious clones corresponding to these two isolates and their chimeric clones resulting from swapping different regions of genomic cDNA between L and L-RB were constructed. Only L-RB or chimeras containing the middle fragment of L-RB cDNA showed virulence on Rsv4-genotype soybean. Sequence comparison analysis revealed that the middle genomic region of L and L-RB encodes four different amino acids. Point mutagenesis demonstrated that a single amino acid substitution (Q1033K) in the P3 protein determined virulence toward Rsv4 resistance. In addition, six new SMV Rsv4 resistance-breaking isolates, variants of the second passage on Williams 82 infected with the chimeras or mutants noninfectious on soybean carrying Rsv4, were obtained. Sequencing data indicated that these new isolates contain either the Q1033K mutation or a new substitution (G1054R) in P3. Site-directed mutagenesis confirmed the virulence role of the G1054R mutation on Rsv4-genotype soybean. Taken together, these data suggest that P3 of the SMV G2 strain is an avirulent determinant for Rsv4 and one single nucleotide mutation in P3 may be sufficient to compromise its elicitor function.

Molecular characterization and experimental host range of an isolate of Macroptilium golden mosaic virus that infects Wissadula amplissima in Jamaica.

Partial genome sequences for the tentative begomovirus Macroptilium golden mosaic virus (MGMV) have been previously reported and were originally obtained for an isolate that infected Macroptilium lathyroides in Jamaica. In this study, we PCR-amplified, cloned and determined the sequence for the complete genome of isolates of MGMV that we found infecting Wissadula amplissima collected from August Town and Spanish Town, Jamaica. Sequence analysis confirmed that MGMV is a distinct begomovirus species, based on the ICTV 89% rule for species demarcation. MGMV shared its highest nucleotide identity at 79% for DNA-A component and 66% for DNA-B component to Corchorus yellow spot virus [Mexico:Yucatan:2005]. The names Macroptilium golden mosaic virus [Jamaica1:Wissadula:AugustTown] (MGMV [JM1:Wd:AT]) and Macroptilium golden mosaic virus [Jamaica1:Wissadula:SpanishTown] (MGMV [JM1:Wd:ST]) are proposed herein for the MGMV isolates from August Town and Spanish Town, respectively. The genome organization of MGMV [JM1:Wd:AT] and MGMV [JM1:Wd:ST] is characteristic of Western Hemisphere bipartite begomoviruses. Excluding the replication enhancer protein (REn), all proteins encoded by the MGMV [JM1:Wd:AT] and MGMV [JM1:Wd:ST] genomes are most similar to their counterparts in Western Hemisphere begomoviruses. The REn proteins of MGMV [JM1:Wd:AT] and MGMV [JM1:Wd:ST], share greatest similarity to the REn protein of Corchorus yellow vein virus [Vietnam:Hoa Binh:2000], a New World-like begomovirus identified in Asia. Phylogenetic reconstruction places MGMV in a clade containing Potato yellow mosaic virus. Results of an experimental host range study indicated that MGMV [JM1:Wd:AT] can infect kidney bean, hot pepper and tomato.

The Tobacco etch virus P3 protein forms mobile inclusions via the early secretory pathway and traffics along actin microfilaments.

Plant potyviruses encode two membrane proteins, 6K and P3. The 6K protein has been shown to induce virus replication vesicles. However, the function of P3 remains unclear. In this study, subcellular localization of the Tobacco etch virus (TEV) P3 protein was investigated in Nicotiana benthamiana leaf cells. The TEV P3 protein localized on the endoplasmic reticulum (ER) membrane and formed punctate inclusions in association with the Golgi apparatus. The trafficking of P3 to the Golgi was mediated by the early secretory pathway. The Golgi-associated punctate structures originated from the ER exit site (ERES). Deletion analyses identified P3 domains required for the retention of P3 at the Golgi. Moreover, the P3 punctate structure was found to traffic along the actin filaments and colocalize with the 6K-containing replication vesicles. Taken together, these data support previous suggestions that P3 may play dual roles in virus movement and replication.

Characterization of the cassava geminivirus transcription activation protein putative nuclear localization signal.

The geminivirus transcription activation protein (TrAP) localizes to the nucleus and contains a putative nuclear localization signal (NLS) ((28)PRRRR(32)) on the N-terminus. The role of individual residues of this putative NLS on nuclear localization and symptom induction was investigated using TrAP of East African cassava mosaic Cameroon virus (EACMCV). Subcellular localization was conducted using the green fluorescent protein (GFP). Results showed that the proline residue at position 28 (Pro-28) is essential for nuclear localization whereas individually, none of the four contiguous arginines is necessary for nuclear targeting. The role of each of the five NLS amino acid residues on TrAP-mediated disease phenotype and gene silencing suppression was investigated by expressing these mutants in Nicotiana benthamiana from the PVX vector and under the control of the Cauliflower mosaic virus 35S promoter. Results showed that all five residues of the NLS play a role on disease phenotype production in N. benthamiana plants. Furthermore, each of the NLS residues appeared to be required for suppression of VIGS but appeared not to be required for the ability of TrAP to transactivate transcription and interact with adenosine kinase (ADK).

Role of a geminivirus AV2 protein putative protein kinase C motif on subcellular localization and pathogenicity.

Virus-derived genes or genome fragments are increasingly being used to generate transgenic plants with resistance to plant viruses. There is need to rapidly investigate these genes in plants using transient expression prior to using them as transgenes since they may be pathogenic to plants. In this study, we investigated the AV2 protein encoded by East African cassava mosaic Cameroon virus, a virus associated with a cassava disease epidemic in western Africa. For subcellular localization, AV2 was fused to the yellow fluorescent protein (YFP) and expressed in Nicotiana benthamiana. Confocal analyses showed that AV2-YFP localizes mainly in the cytoplasm. Because it overlaps with the coat protein gene and therefore could be used to generate transgenic plants for resistance to geminiviruses, we investigated its pathogenesis in N. benthamiana by using the Potato virus X (PVX) vector. The chimeric virus PVX-AV2 induced a mild mottling in infected plants and was shown to suppress virus-induced gene silencing (VIGS). Using point mutations, we show here that AV2 pathogenicity is dependent on a conserved putative protein kinase C (PKC) phosphorylation motif. Because of its pathogenicity and ability to suppress RNA silencing, AV2 transgenic plants will less likely provide a control to geminiviruses, indeed it may weaken the resistance of the plant. We therefore suggest the use of the AV2 putative PKC mutants to generate transgenic plants.