Genome-wide mapping of NBS-LRR genes and their association with disease resistance in soybean

Genome-Wide Characterization of NBS-Encoding Genes in Watermelon and Their Potential Association with Gummy Stem Blight Resistance

Watermelon (Citrullus lanatus) is a nutritionally rich and economically important horticultural crop of the Cucurbitaceae family. Gummy stem blight (GSB) is a major disease of watermelon, which is caused by the fungus Didymella bryoniae, and results in substantial economic losses in terms of yield and quality. However, only a few molecular studies have focused on GSB resistance in watermelon. Nucleotide binding site (NBS)-encoding resistance (R) genes play important roles in plant defense responses to several pathogens, but little is known about the role of NBS-encoding genes in disease resistance in watermelon. The analyzed NBS-encoding R genes comprises several domains, including Toll/interleukin-1 receptor(TIR), NBS, leucine-rich repeat (LRR), resistance to powdery mildew8(RPW8) and coiled coil (CC), which are known to be involved in disease resistance. We determined the expression patterns of these R genes in resistant and susceptible watermelon lines at different time points after D. bryoniae infection by quantitative RT-PCR. The R genes exhibited various expression patterns in the resistant watermelon compared to the susceptible watermelon. Only six R genes exhibited consistent expression patterns (Cla001821, Cla019863, Cla020705, Cla012430, Cla012433 and Cla012439), which were higher in the resistant line compared to the susceptible line. Our study provides fundamental insights into the NBS-LRR gene family in watermelon in response to D. bryoniae infection. Further functional studies of these six candidate resistance genes should help to advance breeding programs aimed at improving disease resistance in watermelons.

Identification of expressed R-genes associated with leaf spot diseases in cultivated peanut

Peanut (Arachis hypogaea L.) is an important food and oilseed crop worldwide. Yield and quality can be significantly reduced by foliar fungal diseases, such as early and late leaf spot diseases. Acceptable levels of leaf spot resistance in cultivated peanut have been elusive due to environmental interactions and the proper combination of QTLs in any particular peanut genotype. Resistance gene analogs, as potential resistance (R)-genes, have unique roles in the recognition and activation of disease resistance responses. Novel R-genes can be identified by searches for conserved domains such as nucleotide binding site, leucine rich repeat, receptor like kinase, and receptor like protein from expressed genes or through genomic sequences. Expressed R-genes represent necessary plant signals in a disease response. The goals of this research are to identify expressed R-genes from cultivated peanuts that are naturally infected by early and late spot pathogens, compare these to the closest diploid progenitors, and evaluate specific gene expression in cultivated peanuts. Putative peanut R-genes (381) were available from a public database (NCBI). Primers were designed and PCR products were sequenced. A total of 214 sequences were produced which matched to proteins with the corresponding R-gene motifs. These R-genes were mapped to the genome sequences of Arachis duranensis and Arachis ipaensis, which are the closest diploid progenitors for tetraploid cultivated peanut, A. hypogaea. Identification and association of specific gene-expression will elucidate potential disease resistance mechanism in peanut and may facilitate the selection of breeding lines with high levels of leaf spot resistance.

Evolutionary Divergence of TNL Disease-Resistant Proteins in Soybean (Glycine max) and Common Bean (Phaseolus vulgaris)

Disease-resistant genes (R genes) encode proteins that are involved in protecting plants from their pathogens and pests. Availability of complete genome sequences from soybean and common bean allowed us to perform a genome-wide identification and analysis of the Toll interleukin-1 receptor-like nucleotide-binding site leucine-rich repeat (TNL) proteins. Hidden Markov model (HMM) profiling of all protein sequences resulted in the identification of 117 and 77 regular TNL genes in soybean and common bean, respectively. We also identified TNL gene homologs with unique domains, and signal peptides as well as nuclear localization signals. The TNL genes in soybean formed 28 clusters located on 10 of the 20 chromosomes, with the majority found on chromosome 3, 6 and 16. Similarly, the TNL genes in common bean formed 14 clusters located on five of the 11 chromosomes, with the majority found on chromosome 10. Phylogenetic analyses of the TNL genes from Arabidopsis, soybean and common bean revealed less divergence within legumes relative to the divergence between legumes and Arabidopsis. Syntenic blocks were found between chromosomes Pv10 and Gm03, Pv07 and Gm10, as well as Pv01 and Gm14. The gene expression data revealed basal level expression and tissue specificity, while analysis of available microRNA data showed 37 predicted microRNA families involved in targeting the identified TNL genes in soybean and common bean.

Genome-Wide Identification of NBS-Encoding Resistance Genes in Sunflower (Helianthus annuus L.)

Nucleotide Binding Site—Leucine-Rich Repeat (NBS-LRR) genes encode disease resistance proteins involved in plants’ defense against their pathogens. Although sunflower is affected by many diseases, only a few molecular details have been uncovered regarding pathogenesis and resistance mechanisms. Recent availability of sunflower whole genome sequences in publicly accessible databases allowed us to accomplish a genome-wide identification of Toll-interleukin-1 receptor-like Nucleotide-binding site Leucine-rich repeat (TNL), Coiled Coil (CC)-NBS-LRR (CNL), Resistance to powdery mildew8 (RPW8)-NBS-LRR (RNL) and NBS-LRR (NL) protein encoding genes. Hidden Markov Model (HMM) profiling of 52,243 putative protein sequences from sunflower resulted in 352 NBS-encoding genes, among which 100 genes belong to CNL group including 64 genes with RX_CC like domain, 77 to TNL, 13 to RNL, and 162 belong to NL group. We also identified signal peptides and nuclear localization signals present in the identified genes and their homologs. We found that NBS genes were located on all chromosomes and formed 75 gene clusters, one-third of which were located on chromosome 13. Phylogenetic analyses between sunflower and Arabidopsis NBS genes revealed a clade-specific nesting pattern in CNLs, with RNLs nested in the CNL-A clade, and species-specific nesting pattern for TNLs. Surprisingly, we found a moderate bootstrap support (BS = 50%) for CNL-A clade being nested within TNL clade making both the CNL and TNL clades paraphyletic. Arabidopsis and sunflower showed 87 syntenic blocks with 1049 high synteny hits between chromosome 5 of Arabidopsis and chromosome 6 of sunflower. Expression data revealed functional divergence of the NBS genes with basal level tissue-specific expression. This study represents the first genome-wide identification of NBS genes in sunflower paving avenues for functional characterization and potential crop improvement.

What Slows Down Phytoplasma Proliferation? Speculations on the Involvement of AtSEOR2 Protein in Plant Defence Signalling

Considering the crude methods used to control phytoplasma diseases, a deeper knowledge on the defence mechanisms recruited by the plant to face phytoplasma invasion is required. Recently, we demonstrated that Arabidopsis mutants lacking AtSEOR1 gene showed a low phytoplasma titre. In wild type plants AtSEOR1 and AtSEOR2 are tied in filamentous proteins. Knockout of the AtSEOR1 gene may pave the way for an involvement of free AtSEOR2 proteins in defence mechanisms. Among the proteins conferring resistance against pathogenic bacteria, AtRPM1-interacting protein has been found to interact with AtSEOR2 in a high-quality, matrix-based yeast-two hybrid assay. For this reason, we investigated the expression levels of Arabidopsis AtRIN4, and the associated AtRPM1 and AtRPS2 genes in healthy and Chrysanthemum yellows-infected wild-type and Atseor1ko lines.

Assembly and annotation of a draft genome sequence for Glycine latifolia, a perennial wild relative of soybean

Glycine latifolia (Benth.) Newell & Hymowitz (2n = 40), one of the 27 wild perennial relatives of soybean, possesses genetic diversity and agronomically favorable traits that are lacking in soybean. Here, we report the 939-Mb draft genome assembly of G. latifolia (PI 559298) using exclusively linked-reads sequenced from a single Chromium library. We organized scaffolds into 20 chromosome-scale pseudomolecules utilizing two genetic maps and the Glycine max (L.) Merr. genome sequence. High copy numbers of putative 91-bp centromere-specific tandem repeats were observed in consecutive blocks within predicted pericentromeric regions on several pseudomolecules. No 92-bp putative centromeric repeats, which are abundant in G. max, were detected in G. latifolia or Glycine tomentella. Annotation of the assembled genome and subsequent filtering yielded a high confidence gene set of 54 475 protein-coding loci. In comparative analysis with five legume species, genes related to defense responses were significantly overrepresented in Glycine-specific orthologous gene families. A total of 304 putative nucleotide-binding site (NBS)-leucine-rich-repeat (LRR) genes were identified in this genome assembly. Different from other legume species, we observed a scarcity of TIR-NBS-LRR genes in G. latifolia. The G. latifolia genome was also predicted to contain genes encoding 367 LRR-receptor-like kinases, a family of proteins involved in basal defense responses and responses to abiotic stress. The genome sequence and annotation of G. latifolia provides a valuable source of alternative alleles and novel genes to facilitate soybean improvement. This study also highlights the efficacy and cost-effectiveness of the application of Chromium linked-reads in diploid plant genome de novo assembly.

Regulation and Evolution of NLR Genes: A Close Interconnection for Plant Immunity

NLR (NOD-like receptor) genes belong to one of the largest gene families in plants. Their role in plants' resistance to pathogens has been clearly described for many members of this gene family, and dysregulation or overexpression of some of these genes has been shown to induce an autoimmunity state that strongly affects plant growth and yield. For this reason, these genes have to be tightly regulated in their expression and activity, and several regulatory mechanisms are described here that tune their gene expression and protein levels. This gene family is subjected to rapid evolution, and to maintain diversity at NLRs, a plethora of genetic mechanisms have been identified as sources of variation. Interestingly, regulation of gene expression and evolution of this gene family are two strictly interconnected aspects. Indeed, some examples have been reported in which mechanisms of gene expression regulation have roles in promotion of the evolution of this gene family. Moreover, co-evolution of the NLR gene family and other gene families devoted to their control has been recently demonstrated, as in the case of miRNAs.

Salicylic acid and broad spectrum of NBS-LRR family genes are involved in SMV-soybean interactions

Soybean mosaic virus (SMV) is a severe pathogen reducing crop yield and seed quality of soybean. Although several resistance gene loci including Rsv1, Rsv3 and Rsv4 are identified in some soybean varieties, most of the soybean genes related to SMV infection are still not characterized. In order to reveal genome-wide gene expression profiles in response to SMV infection, we used transcriptome analysis to determine SMV-responsive genes in susceptible variety Hefeng25. Time course RNA-seq analysis at 1, 5 and 10 dpi identified many deregulated pathways and gene families. "Plant-pathogen interaction" pathway with KEGG No. of KO04626 was highly enriched and dozens of NBS-LRR family genes were significantly down-regulated at 5 dpi. qRT-PCR analyses were performed to verify expression patterns of these genes and most were in accordance with the RNA-seq data. As NBS-LRR family proteins are broadly involved in plant immunity responses, our results indicated the importance of this time point (5 dpi) for SMV-soybean interaction. Consistent with it, SMV titer was increased from 1 dpi to 10 dpi and peaked at 5 dpi. Expression of SA (salicylic acid) marker gene PR-1 was induced by SMV infection. Application of exogenous MeSA, an active form of SA, primed the plant resistant to virus infection and reduced SMV accumulation in soybean. Interestingly, MeSA treatment also significantly upregulated expressions of SMV-responsive NBS-LRR genes. Compared with susceptible line Hefeng25, endogenous SA level was higher and was consistently induced by SMV infection in resistant variety RV8143. Moreover, expressions of NBS-LRR family genes were up-regulated by SMV infection in RV8143, while they were down-regulated by SMV infection in Hefeng25. Our results implied that SA and NBS-LRR family genes were involved in SMV-soybean interaction. SMV could compromise soybean defense responses by repression of NBS-LRR family genes in Hefeng25, and SA was implicated in this interaction process.

Fine-mapping and identification of a novel locus Rsc15 underlying soybean resistance to Soybean mosaic virus

KEY MESSAGE

Rsc15, a novel locus underlying soybean resistance to SMV, was fine mapped to a 95-kb region on chromosome 6. The Rsc15- mediated resistance is likely attributed to the gene GmPEX14 , the relative expression of which was highly correlated with the accumulation of H 2 O 2 along with the activities of POD and CAT during the early stages of SMV infection in RN-9. Soybean mosaic virus (SMV) causes severe yield losses and seed quality deterioration in soybean [Glycine max (L.) Merr.] worldwide. A series of single dominant SMV resistance genes have been identified on respective soybean chromosomes 2, 13 and 14, while one novel locus, Rsc15, underlying resistance to the virulent SMV strain SC15 from soybean cultivar RN-9 has been recently mapped to a 14.6-cM region on chromosome 6. However, candidate gene has not yet been identified within this region. In the present study, we aimed to fine map the Rsc15 region and identify candidate gene(s) for this invaluable locus. High-resolution fine-mapping revealed that the Rsc15 gene was located in a 95-kb genomic region which was flanked by the two simple sequence repeat (SSR) markers SSR_06_17 and BARCSOYSSR_06_0835. Allelic sequence comparison and expression profile analysis of candidate genes inferred that the gene Glyma.06g182600 (designated as GmPEX14) was the best candidate gene attributing for the resistance of Rsc15, and that genes encoding receptor-like kinase (RLK) (i.e., Glyma.06g175100 and Glyma.06g184400) and serine/threonine kinase (STK) (i.e., Glyma.06g182900 and Glyma.06g183500) were also potential candidates. High correlations were established between the relative expression level of GmPEX14 and the hydrogen peroxide (H2O2) concentration and activities of catalase (CAT) and peroxidase (POD) during the early stages of SMV-SC15 infection in RN-9. The results of the present study will be useful in marker-assisted breeding for SMV resistance and will lead to further understanding of the molecular mechanisms of host resistance against SMV.

Comparative Genomics of Non-TNL Disease Resistance Genes from Six Plant Species

Disease resistance genes (R genes), as part of the plant defense system, have coevolved with corresponding pathogen molecules. The main objectives of this project were to identify non-Toll interleukin receptor, nucleotide-binding site, leucine-rich repeat (nTNL) genes and elucidate their evolutionary divergence across six plant genomes. Using reference sequences from Arabidopsis, we investigated nTNL orthologs in the genomes of common bean, Medicago, soybean, poplar, and rice. We used Hidden Markov Models for sequence identification, performed model-based phylogenetic analyses, visualized chromosomal positioning, inferred gene clustering, and assessed gene expression profiles. We analyzed 908 nTNL R genes in the genomes of the six plant species, and classified them into 12 subgroups based on the presence of coiled-coil (CC), nucleotide binding site (NBS), leucine rich repeat (LRR), resistance to Powdery mildew 8 (RPW8), and BED type zinc finger domains. Traditionally classified CC-NBS-LRR (CNL) genes were nested into four clades (CNL A-D) often with abundant, well-supported homogeneous subclades of Type-II R genes. CNL-D members were absent in rice, indicating a unique R gene retention pattern in the rice genome. Genomes from Arabidopsis, common bean, poplar and soybean had one chromosome without any CNL R genes. Medicago and Arabidopsis had the highest and lowest number of gene clusters, respectively. Gene expression analyses suggested unique patterns of expression for each of the CNL clades. Differential gene expression patterns of the nTNL genes were often found to correlate with number of introns and GC content, suggesting structural and functional divergence.

Plant STAND P-loop NTPases: a current perspective of genome distribution, evolution, and function : Plant STAND P-loop NTPases: genomic organization, evolution, and molecular mechanism models contribute broadly to plant pathogen defense

STAND P-loop NTPase is the common weapon used by plant and other organisms from all three kingdoms of life to defend themselves against pathogen invasion. The purpose of this study is to review comprehensively the latest finding of plant STAND P-loop NTPase related to their genomic distribution, evolution, and their mechanism of action. Earlier, the plant STAND P-loop NTPase known to be comprised of only NBS-LRRs/AP-ATPase/NB-ARC ATPase. However, recent finding suggests that genome of early green plants comprised of two types of STAND P-loop NTPases: (1) mammalian NACHT NTPases and (2) NBS-LRRs. Moreover, YchF (unconventional G protein and members of P-loop NTPase) subfamily has been reported to be exceptionally involved in biotic stress (in case of Oryza sativa), thereby a novel member of STAND P-loop NTPase in green plants. The lineage-specific expansion and genome duplication events are responsible for abundance of plant STAND P-loop NTPases; where "moderate tandem and low segmental duplication" trajectory followed in majority of plant species with few exception (equal contribution of tandem and segmental duplication). Since the past decades, systematic research is being investigated into NBS-LRR function supported the direct recognition of pathogen or pathogen effectors by the latest models proposed via 'integrated decoy' or 'sensor domains' model. Here, we integrate the recently published findings together with the previous literature on the genomic distribution, evolution, and distinct models proposed for functional molecular mechanism of plant STAND P-loop NTPases.

Genome-wide identification and tissue-specific expression analysis of nucleotide binding site-leucine rich repeat gene family in Cicer arietinum (kabuli chickpea)

The nucleotide binding site-leucine rich repeat (NBS-LRR) proteins play an important role in the defense mechanisms against pathogens. Using bioinformatics approach, we identified and annotated 104 NBS-LRR genes in chickpea. Phylogenetic analysis points to their diversification into two families namely TIR-NBS-LRR and non-TIR-NBS-LRR. Gene architecture revealed intron gain/loss events in this resistance gene family during their independent evolution into two families. Comparative genomics analysis elucidated its evolutionary relationship with other fabaceae species. Around 50% NBS-LRRs reside in macro-syntenic blocks underlining positional conservation along with sequence conservation of NBS-LRR genes in chickpea. Transcriptome sequencing data provided evidence for their transcription and tissue-specific expression. Four cis-regulatory elements namely WBOX, DRE, CBF, and GCC boxes, that commonly occur in resistance genes, were present in the promoter regions of these genes. Further, the findings will provide a strong background to use candidate disease resistance NBS-encoding genes and identify their specific roles in chickpea.

Genome-Wide Association Study Identifies NBS-LRR-Encoding Genes Related with Anthracnose and Common Bacterial Blight in the Common Bean

Nucleotide-binding site and leucine-rich repeat (NBS-LRR) genes represent the largest and most important disease resistance genes in plants. The genome sequence of the common bean (Phaseolus vulgaris L.) provides valuable data for determining the genomic organization of NBS-LRR genes. However, data on the NBS-LRR genes in the common bean are limited. In total, 178 NBS-LRR-type genes and 145 partial genes (with or without a NBS) located on 11 common bean chromosomes were identified from genome sequences database. Furthermore, 30 NBS-LRR genes were classified into Toll/interleukin-1 receptor (TIR)-NBS-LRR (TNL) types, and 148 NBS-LRR genes were classified into coiled-coil (CC)-NBS-LRR (CNL) types. Moreover, the phylogenetic tree supported the division of these PvNBS genes into two obvious groups, TNL types and CNL types. We also built expression profiles of NBS genes in response to anthracnose and common bacterial blight using qRT-PCR. Finally, we detected nine disease resistance loci for anthracnose (ANT) and seven for common bacterial blight (CBB) using the developed NBS-SSR markers. Among these loci, NSSR24, NSSR73, and NSSR265 may be located at new regions for ANT resistance, while NSSR65 and NSSR260 may be located at new regions for CBB resistance. Furthermore, we validated NSSR24, NSSR65, NSSR73, NSSR260, and NSSR265 using a new natural population. Our results provide useful information regarding the function of the NBS-LRR proteins and will accelerate the functional genomics and evolutionary studies of NBS-LRR genes in food legumes. NBS-SSR markers represent a wide-reaching resource for molecular breeding in the common bean and other food legumes. Collectively, our results should be of broad interest to bean scientists and breeders.

Comparative analysis of NBS-LRR genes and their response to Aspergillus flavus in Arachis

Studies have demonstrated that nucleotide-binding site-leucine-rich repeat (NBS-LRR) genes respond to pathogen attack in plants. Characterization of NBS-LRR genes in peanut is not well documented. The newly released whole genome sequences of Arachis duranensis and Arachis ipaënsis have allowed a global analysis of this important gene family in peanut to be conducted. In this study, we identified 393 (AdNBS) and 437 (AiNBS) NBS-LRR genes from A. duranensis and A. ipaënsis, respectively, using bioinformatics approaches. Full-length sequences of 278 AdNBS and 303 AiNBS were identified. Fifty-one orthologous, four AdNBS paralogous, and six AiNBS paralogous gene pairs were predicted. All paralogous gene pairs were located in the same chromosomes, indicating that tandem duplication was the most likely mechanism forming these paralogs. The paralogs mainly underwent purifying selection, but most LRR 8 domains underwent positive selection. More gene clusters were found in A. ipaënsis than in A. duranensis, possibly owing to tandem duplication events occurring more frequently in A. ipaënsis. The expression profile of NBS-LRR genes was different between A. duranensis and A. hypogaea after Aspergillus flavus infection. The up-regulated expression of NBS-LRR in A. duranensis was continuous, while these genes responded to the pathogen temporally in A. hypogaea.

Genetic Analysis of NBS-LRR Gene Family in Chickpea and Their Expression Profiles in Response to Ascochyta Blight Infection

Ascochyta blight is one of the major diseases of chickpea worldwide. The genetic resistance to ascochyta blight in chickpea is complex and governed by multiple QTLs. However, the molecular mechanism of quantitative disease resistance to ascochyta blight and the genes underlying these QTLs are still unknown. Most often disease resistance is determined by resistance (R) genes. The most predominant R-genes contain nucleotide binding site and leucine rich repeat (NBS-LRR) domains. A total of 121 NBS-LRR genes were identified in the chickpea genome. Ninety-eight of these genes contained all essential conserved domains while 23 genes were truncated. The NBS-LRR genes were grouped into eight distinct classes based on their domain architecture. Phylogenetic analysis grouped these genes into two major clusters based on their structural variation, the first cluster with toll or interleukin-1 like receptor (TIR) domain and the second cluster either with or without a coiled-coil domain. The NBS-LRR genes are distributed unevenly across the eight chickpea chromosomes and nearly 50% of the genes are present in clusters. Thirty of the NBS-LRR genes were co-localized with nine of the previously reported ascochyta blight QTLs and were tested as potential candidate genes for ascochyta blight resistance. Expression pattern of these genes was studied in two resistant (CDC Corinne and CDC Luna) and one susceptible (ICCV 96029) genotypes at different time points after ascochyta blight infection using real-time quantitative PCR. Twenty-seven NBS-LRR genes showed differential expression in response to ascochyta blight infection in at least one genotype at one time point. Among these 27 genes, the majority of the NBS-LRR genes showed differential expression after inoculation in both resistant and susceptible genotypes which indicates the involvement of these genes in response to ascochyta blight infection. Five NBS-LRR genes showed genotype specific expression. Our study provides a new insight of NBS-LRR gene family in chickpea and the potential involvement of NBS-LRR genes in response to ascochyta blight infection.

Identification, characterization, and gene expression analysis of nucleotide binding site (NB)-type resistance gene homologues in switchgrass

BACKGROUND

Switchgrass (Panicum virgatum L.) is a warm-season perennial grass that can be used as a second generation bioenergy crop. However, foliar fungal pathogens, like switchgrass rust, have the potential to significantly reduce switchgrass biomass yield. Despite its importance as a prominent bioenergy crop, a genome-wide comprehensive analysis of NB-LRR disease resistance genes has yet to be performed in switchgrass.

RESULTS

In this study, we used a homology-based computational approach to identify 1011 potential NB-LRR resistance gene homologs (RGHs) in the switchgrass genome (v 1.1). In addition, we identified 40 RGHs that potentially contain unique domains including major sperm protein domain, jacalin-like binding domain, calmodulin-like binding, and thioredoxin. RNA-sequencing analysis of leaf tissue from 'Alamo', a rust-resistant switchgrass cultivar, and 'Dacotah', a rust-susceptible switchgrass cultivar, identified 2634 high quality variants in the RGHs between the two cultivars. RNA-sequencing data from field-grown cultivar 'Summer' plants indicated that the expression of some of these RGHs was developmentally regulated.

CONCLUSIONS

Our results provide useful insight into the molecular structure, distribution, and expression patterns of members of the NB-LRR gene family in switchgrass. These results also provide a foundation for future work aimed at elucidating the molecular mechanisms underlying disease resistance in this important bioenergy crop.

Characterization of Disease Resistance Loci in the USDA Soybean Germplasm Collection Using Genome-Wide Association Studies

Genetic resistance is a key strategy for disease management in soybean. Over the last 50 years, soybean germplasm has been phenotyped for resistance to many pathogens, resulting in the development of disease-resistant elite breeding lines and commercial cultivars. While biparental linkage mapping has been used to identify disease resistance loci, genome-wide association studies (GWAS) using high-density and high-quality markers such as single nucleotide polymorphisms (SNPs) has become a powerful tool to associate molecular markers and phenotypes. The objective of our study was to provide a comprehensive understanding of disease resistance in the United States Department of Agriculture Agricultural Research Service Soybean Germplasm Collection by using phenotypic data in the public Germplasm Resources Information Network and public SNP data (SoySNP50K). We identified SNPs significantly associated with disease ratings from one bacterial disease, five fungal diseases, two diseases caused by nematodes, and three viral diseases. We show that leucine-rich repeat (LRR) receptor-like kinases and nucleotide-binding site-LRR candidate resistance genes were enriched within the linkage disequilibrium regions of the significant SNPs. We review and present a global view of soybean resistance loci against multiple diseases and discuss the power and the challenges of using GWAS to discover disease resistance in soybean.

Genome-wide association mapping of soybean chlorophyll traits based on canopy spectral reflectance and leaf extracts

BACKGROUND

Chlorophyll is a major component of chloroplasts and a better understanding of the genetic basis of chlorophyll in soybean [Glycine max (L.) Merr.] might contribute to improving photosynthetic capacity and yield in regions with adverse environmental conditions. A collection of 332 diverse soybean genotypes were grown in 2 years (2009 and 2010) and chlorophyll a (eChl_A), chlorophyll b (eChl_B), and total chlorophyll (eChl_T) content as well as chlorophyll a/b ratio (eChl_R) in leaf tissues were determined by extraction and spectrometric determination. Total chlorophyll was also derived from canopy spectral reflectance measurements using a model of wavelet transformed spectra (tChl_T) as well as with a spectral reflectance index (iChl_T).

RESULTS

A genome-wide associating mapping approach was employed using 31,253 single nucleotide polymorphisms (SNPs) to identify loci associated with the extract based eChl_A, eChl_B, eChl_R and eChl_T measurements and the two canopy spectral reflectance-based methods (tChl_T and iChl_T). A total of 23 (14 loci), 15 (7 loci) and 14 SNPs (10 loci) showed significant association with eChl_A, eChl_B and eChl_R respectively. A total of 52 unique SNPs were significantly associated with total chlorophyll content based on at least one of the three approaches (eChl_T, tChl_T and iChl_T) and likely tagged 27 putative loci for total chlorophyll content, four of which were indicated by all three approaches.

CONCLUSIONS

Results presented here show that markers for chlorophyll traits can be identified in soybean using both extract-based and canopy spectral reflectance-based phenotypes, and confirm that high-throughput phenotyping-amenable canopy spectral reflectance measurements can be used for association mapping.

Comparative Transcriptome Analysis of Resistant and Susceptible Common Bean Genotypes in Response to Soybean Cyst Nematode Infection

Soybean cyst nematode (SCN; Heterodera glycines Ichinohe) reproduces on the roots of common bean (Phaseolus vulgaris L.) and can cause reductions in plant growth and seed yield. The molecular changes in common bean roots caused by SCN infection are unknown. Identification of genetic factors associated with SCN resistance could help in development of improved bean varieties with high SCN resistance. Gene expression profiling was conducted on common bean roots infected by SCN HG type 0 using next generation RNA sequencing technology. Two pinto bean genotypes, PI533561 and GTS-900, resistant and susceptible to SCN infection, respectively, were used as RNA sources eight days post inoculation. Total reads generated ranged between ~ 3.2 and 5.7 million per library and were mapped to the common bean reference genome. Approximately 70-90% of filtered RNA-seq reads uniquely mapped to the reference genome. In the inoculated roots of resistant genotype PI533561, a total of 353 genes were differentially expressed with 154 up-regulated genes and 199 down-regulated genes when compared to the transcriptome of non- inoculated roots. On the other hand, 990 genes were differentially expressed in SCN-inoculated roots of susceptible genotype GTS-900 with 406 up-regulated and 584 down-regulated genes when compared to non-inoculated roots. Genes encoding nucleotide-binding site leucine-rich repeat resistance (NLR) proteins, WRKY transcription factors, pathogenesis-related (PR) proteins and heat shock proteins involved in diverse biological processes were differentially expressed in both resistant and susceptible genotypes. Overall, suppression of the photosystem was observed in both the responses. Furthermore, RNA-seq results were validated through quantitative real time PCR. This is the first report describing genes/transcripts involved in SCN-common bean interaction and the results will have important implications for further characterization of SCN resistance genes in common bean.

Molecular phylogeny and dynamic evolution of disease resistance genes in the legume family

Background

Legumes are the second-most important crop family in agriculture for its economic and nutritional values. Disease resistance (R-) genes play an important role in responding to pathogen infections in plants. To further increase the yield of legume crops, we need a comprehensive understanding of the evolution of R-genes in the legume family.

Results

In this study, we developed a robust pipeline and identified a total of 4,217 R-genes in the genomes of seven sequenced legume species. A dramatic diversity of R-genes with structural variances indicated a rapid birth-and-death rate during the R-gene evolution in legumes. The number of R-genes transiently expanded and then quickly contracted after whole-genome duplications, which meant that R-genes were sensitive to subsequent diploidization. R proteins with the Coiled-coil (CC) domain are more conserved than others in legumes. Meanwhile, other types of legume R proteins with only one or two typical domains were subjected to higher rates of loss during evolution. Although R-genes evolved quickly in legumes, they tended to undergo purifying selection instead of positive selection during evolution. In addition, domestication events in some legume species preferentially selected for the genes directly involved in the plant-pathogen interaction pathway while suppressing those R-genes with low occurrence rates.

Conclusions

Our results provide insights into the dynamic evolution of R-genes in the legume family, which will be valuable for facilitating genetic improvements in the disease resistance of legume cultivars.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2736-9) contains supplementary material, which is available to authorized users.

Morphological features of different polyploids for adaptation and molecular characterization of CC-NBS-LRR and LEA gene families in Agave L

Polyploidy has been widely described in many Agave L. species, but its influence on environmental response to stress is still unknown. With the objective of knowing the morphological adaptations and regulation responses of genes related to biotic (LEA) and abiotic (NBS-LRR) stress in species of Agave with different levels of ploidy, and how these factors contribute to major response of Agave against environmental stresses, we analyzed 16 morphological trials on five accessions of three species (Agave tequilana Weber, Agave angustifolia Haw. and Agave fourcroydes Lem.) with different ploidy levels (2n=2x=60 2n=3x=90, 2n=5x=150, 2n=6x=180) and evaluated the expression of NBS-LRR and LEA genes regulated by biotic and abiotic stress. It was possible to associate some morphological traits (spines, nuclei, and stomata) to ploidy level. The genetic characterization of stress-related genes NBS-LRR induced by pathogenic infection and LEA by heat or saline stresses indicated that amino acid sequence analysis in these genes showed more substitutions in higher ploidy level accessions of A. fourcroydes Lem. 'Sac Ki' (2n=5x=150) and A. angustifolia Haw. 'Chelem Ki' (2n=6x=180), and a higher LEA and NBS-LRR representativeness when compared to their diploid and triploid counterparts. In all studied Agave accessions expression of LEA and NBS-LRR genes was induced by saline or heat stresses or by infection with Erwinia carotovora, respectively. The transcriptional activation was also higher in A. angustifolia Haw. 'Chelem Ki' (2n=6x=180) and A. fourcroydes 'Sac Ki' (2n=5x=150) than in their diploid and triploid counterparts, which suggests higher adaptation to stress. Finally, the diploid accession A. tequilana Weber 'Azul' showed a differentiated genetic profile relative to other Agave accessions. The differences include similar or higher genetic representativeness and transcript accumulation of LEA and NBS-LRR genes than in polyploid (2n=5x=150 and 2n=6x=180) Agave accessions, thus suggesting a differentiated selection pressure for overcoming the lower ploidy level of the diploid A. tequilana Weber 'Azul'.

TNL genes in peach: insights into the post-LRR domain

Background

Plants develop sustainable defence responses to pathogen attacks through resistance (R) genes contributing to effector-triggered immunity (ETI). TIR-NB-LRR genes (TNL genes) constitute a major family of ETI R genes in dicots. The putative functions or roles of the TIR, NB and LRR domains of the proteins they encode (TNLs) are well documented, but TNLs also have a poorly characterised C-terminal region, the function of which is unknown in most cases. We characterised this prevalent stress-response protein family in a perennial plant, using the genome of peach (Prunus persica), the model Prunus species. The first TNL gene from this genus to be cloned, the Ma gene, confers complete-spectrum resistance to root-knot nematodes (RKNs) and encodes a protein with a huge C-terminal region with five duplicated post-LRR (PL) domains. This gene was the cornerstone of this study.

Results

We investigated the role of this C-terminal region, by first describing the frequency, distribution and structural characteristics of i) TNL genes and ii) their PL domains in the peach genome, using the v1.0 Sanger sequence together with the v2.0 sequence, which has better genome annotation due to the incorporation of transcriptomic data. We detected 195 predicted TNL genes from the eight peach chromosomes: 85 % of these genes mapped to chromosomes 1, 2, 7 and 8. We reconstructed the putative structure of the predicted exons of all the TNL genes identified, and it was possible to retrieve the PL domains among two thirds of the TNL genes. We used our predicted TNL gene sequences to develop an annotation file for use with the Gbrowse tool in the v2.0 genome. The use of these annotation data made it possible to detect transcribed PL sequences in two Prunus species. We then used consensus sequences defined on the basis of 124 PL domains to design specific motifs, and we found that the use of these motifs significantly increased the numbers of PL domains and correlative TNL genes detected in diverse dicot genomes. Based on PL signatures, we showed that TNL genes with multiple PL domains were rare in peach and the other plants screened. The five-PL domain pattern is probably unique to Ma and its orthologues within Prunus and closely related genera from the Rosaceae and was probably inherited from the common ancestor of these plants in the subfamily Spiraeoideae.

Conclusions

The first physical TNL gene map for Prunus species can be used for the further investigation of R genes in this genus. The PL signature motifs are a complementary tool for the detection of TNL R genes in dicots. The low degree of similarity between PL domains and the neighbouring LRR exons and the specificity of PL signature motifs suggest that PL and LRR domains have different origins, with PL domains being specific to TNL genes, and possibly essential to the functioning of these genes in some cases. Investigations of the role of the oversized Ma PL region, in ligand binding or intramolecular interactions for example, may help to enrich our understanding of NB-LRR-mediated plant immunity to RKNs.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2635-0) contains supplementary material, which is available to authorized users.

Optimization of a Virus-Induced Gene Silencing System with Soybean yellow common mosaic virus for Gene Function Studies in Soybeans

Virus-induced gene silencing (VIGS) is an effective tool for the study of soybean gene function. Successful VIGS depends on the interaction between virus spread and plant growth, which can be influenced by environmental conditions. Recently, we developed a new VIGS system derived from the Soybean yellow common mosaic virus (SYCMV). Here, we investigated several environmental and developmental factors to improve the efficiency of a SYCMV-based VIGS system to optimize the functional analysis of the soybean. Following SYCMV: Glycine max-phytoene desaturase (GmPDS) infiltration, we investigated the effect of photoperiod, inoculation time, concentration of Agrobacterium inoculm, and growth temperature on VIGS efficiency. In addition, the relative expression of GmPDS between non-silenced and silenced plants was measured by qRT-PCR. We found that gene silencing efficiency was highest at a photoperiod of 16/8 h (light/dark) at a growth temperature of approximately 27°C following syringe infiltration to unrolled unifoliolate leaves in cotyledon stage with a final SYCMV:GmPDS optimal density (OD)600 of 2.0. Using this optimized protocol, we achieved high efficiency of GmPDS-silencing in various soybean germplasms including cultivated and wild soybeans. We also confirmed that VIGS occurred in the entire plant, including the root, stem, leaves, and flowers, and could transmit GmPDS to other soybean germplasms via mechanical inoculation. This optimized protocol using a SYCMV-based VIGS system in the soybean should provide a fast and effective method to elucidate gene functions and for use in large-scale screening experiments.

QTL mapping for bacterial wilt resistance in peanut (Arachis hypogaea L.)

Bacterial wilt (BW) caused by Ralstonia solanacearum is a serious, global, disease of peanut (Arachis hypogaea L.), but it is especially destructive in China. Identification of DNA markers linked to the resistance to this disease will help peanut breeders efficiently develop resistant cultivars through molecular breeding. A F₂ population, from a cross between disease-resistant and disease-susceptible cultivars, was used to detect quantitative trait loci (QTL) associated with the resistance to this disease in the cultivated peanut. Genome-wide SNPs were identified from restriction-site-associated DNA sequencing tags using next-generation DNA sequencing technology. SNPs linked to disease resistance were determined in two bulks of 30 resistant and 30 susceptible plants along with two parental plants using bulk segregant analysis. Polymorphic SSR and SNP markers were utilized for construction of a linkage map and for performing the QTL analysis, and a moderately dense linkage map was constructed in the F₂ population. Two QTL (qBW-1 and qBW-2) detected for resistance to BW disease were located in the linkage groups LG1 and LG10 and account for 21 and 12 % of the bacterial wilt phenotypic variance. To confirm these QTL, the F₈ RIL population with 223 plants was utilized for genotyping and phenotyping plants by year and location as compared to the F₂ population. The QTL qBW-1 was consistent in the location of LG1 in the F₈ population though the QTL qBW-2 could not be clarified due to fewer markers used and mapped in LG10. The QTL qBW-1, including four linked SNP markers and one SSR marker within 14.4-cM interval in the F₈, was closely related to a disease resistance gene homolog and was considered as a candidate gene for resistance to BW. QTL identified in this study would be useful to conduct marker-assisted selection and may permit cloning of resistance genes. Our study shows that bulk segregant analysis of genome-wide SNPs is a useful approach to expedite the identification of genetic markers linked to disease resistance traits in the allotetraploidy species peanut.

The wheat homolog of putative nucleotide-binding site-leucine-rich repeat resistance gene TaRGA contributes to resistance against powdery mildew

Powdery mildew, one of the most destructive wheat diseases worldwide, is caused by Blumeria graminis f. sp. tritici (Bgt), a fungal species with a consistently high mutation rate that makes individual resistance (R) genes ineffective. Therefore, effective resistance-related gene cloning is vital for breeding and studying the resistance mechanisms of the disease. In this study, a putative nucleotide-binding site-leucine-rich repeat (NBS-LRR) R gene (TaRGA) was cloned using a homology-based cloning strategy and analyzed for its effect on powdery mildew disease and wheat defense responses. Real-time reverse transcription-PCR (RT-PCR) analyses revealed that a Bgt isolate 15 and salicylic acid stimulation significantly induced TaRGA in the resistant variety. Furthermore, the silencing of TaRGA in powdery mildew-resistant plants increased susceptibility to Bgt15 and prompted conidia propagation at the infection site. However, the expression of TaRGA in leaf segments after single-cell transient expression assay highly increased the defense responses to Bgt15 by enhancing callose deposition and phenolic autofluorogen accumulation at the pathogen invading sites. Meanwhile, the expression of pathogenesis-related genes decreased in the TaRGA-silenced plants and increased in the TaRGA-transient-overexpressing leaf segments. These results implied that the TaRGA gene positively regulates the defense response to powdery mildew disease in wheat.

Genome-wide association and epistasis studies unravel the genetic architecture of sudden death syndrome resistance in soybean

Soybean [Glycine max (L.) Merr.] is an economically important crop that is grown worldwide. Sudden death syndrome (SDS), caused by Fusarium virguliforme, is one of the top yield-limiting diseases in soybean. However, the genetic basis of SDS resistance, especially with respect to epistatic interactions, is still unclear. To better understand the genetic architecture of soybean SDS resistance, genome-wide association and epistasis studies were performed using a population of 214 germplasm accessions and 31,914 SNPs from the SoySNP50K Illumina Infinium BeadChip. Twelve loci and 12 SNP-SNP interactions associated with SDS resistance were identified at various time points after inoculation. These additive and epistatic loci together explained 24-52% of the phenotypic variance. Disease-resistant, pathogenesis-related and chitin- and wound-responsive genes were identified in the proximity of peak SNPs, including stress-induced receptor-like kinase gene 1 (SIK1), which is pinpointed by a trait-associated SNP and encodes a leucine-rich repeat-containing protein. We report that the proportion of phenotypic variance explained by identified loci may be considerably improved by taking epistatic effects into account. This study shows the necessity of considering epistatic effects in soybean SDS resistance breeding using marker-assisted and genomic selection approaches. Based on our findings, we propose a model for soybean root defense against the SDS pathogen. Our results facilitate identification of the molecular mechanism underlying SDS resistance in soybean, and provide a genetic basis for improvement of soybean SDS resistance through breeding strategies based on additive and epistatic effects.

Large-Scale Transcriptome Analysis in Faba Bean (Vicia faba L.) under Ascochyta fabae Infection

Faba bean is an important food crop worldwide. However, progress in faba bean genomics lags far behind that of model systems due to limited availability of genetic and genomic information. Using the Illumina platform the faba bean transcriptome from leaves of two lines (29H and Vf136) subjected to Ascochyta fabae infection have been characterized. De novo transcriptome assembly provided a total of 39,185 different transcripts that were functionally annotated, and among these, 13,266 were assigned to gene ontology against Arabidopsis. Quality of the assembly was validated by RT-qPCR amplification of selected transcripts differentially expressed. Comparison of faba bean transcripts with those of better-characterized plant genomes such as Arabidopsis thaliana, Medicago truncatula and Cicer arietinum revealed a sequence similarity of 68.3%, 72.8% and 81.27%, respectively. Moreover, 39,060 single nucleotide polymorphism (SNP) and 3,669 InDels were identified for genotyping applications. Mapping of the sequence reads generated onto the assembled transcripts showed that 393 and 457 transcripts were overexpressed in the resistant (29H) and susceptible genotype (Vf136), respectively. Transcripts involved in plant-pathogen interactions such as leucine rich proteins (LRR) or plant growth regulators involved in plant adaptation to abiotic and biotic stresses were found to be differently expressed in the resistant line. The results reported here represent the most comprehensive transcript database developed so far in faba bean, providing valuable information that could be used to gain insight into the pathways involved in the resistance mechanism against A. fabae and to identify potential resistance genes to be further used in marker assisted selection.

Detection of novel QTLs for foxglove aphid resistance in soybean

KEY MESSAGE

The Raso2 , novel QTL for Korea biotype foxglove aphid resistance in soybean from PI 366121 was identified on chromosome 7 using GoldenGate SNP microarray. Foxglove aphid, Aulacorthum solani (Kaltenbach), is a hemipteran insect that infects a wide variety of plants worldwide and causes serious yield losses in crops. The objective of this study was to identify the putative QTL for foxglove aphid resistance in wild soybean, PI 366121, (Glycine soja Sieb. and Zucc.). One hundred and forty-one F4-derived F8 recombinant inbred lines developed from a cross of susceptible Williams 82 and PI 366121 were used. The phenotyping of antibiosis and antixenosis resistance was done through choice and no-choice tests with total plant damage and primary infestation leaf damage; a genome-wide molecular linkage map was constructed with 504 single-nucleotide polymorphism markers utilizing a GoldenGate assay. Using inclusive composite interval mapping analysis for foxglove aphid resistance, one major candidate QTL on chromosome 7 and three minor QTL regions on chromosomes 3, 6 and 18 were identified. The major QTL on chromosome 7 showed both antixenosis and antibiosis resistance responses. However, the minor QTLs showed only antixenosis resistance response. The major QTL mapped to a different chromosome than the previously identified foxglove aphid resistance QTL, Raso1, from the cultivar Adams. Also, the responses to the Korea biotype foxglove aphid were different for Raso1, and the gene from PI 366121 against the Korea biotype foxglove aphid was different. Thus, the foxglove aphid resistance gene from PI 366121 was determined to be an independent gene from Raso1 and was designated as Raso2. This result could be useful in breeding for new foxglove aphid-resistant soybean cultivars.

Genome-Wide Profiling of Histone Modifications (H3K9me2 and H4K12ac) and Gene Expression in Rust (Uromyces appendiculatus) Inoculated Common Bean (Phaseolus vulgaris L.)

Histone modifications such as methylation and acetylation play a significant role in controlling gene expression in unstressed and stressed plants. Genome-wide analysis of such stress-responsive modifications and genes in non-model crops is limited. We report the genome-wide profiling of histone methylation (H3K9me2) and acetylation (H4K12ac) in common bean (Phaseolus vulgaris L.) under rust (Uromyces appendiculatus) stress using two high-throughput approaches, chromatin immunoprecipitation sequencing (ChIP-Seq) and RNA sequencing (RNA-Seq). ChIP-Seq analysis revealed 1,235 and 556 histone methylation and acetylation responsive genes from common bean leaves treated with the rust pathogen at 0, 12 and 84 hour-after-inoculation (hai), while RNA-Seq analysis identified 145 and 1,763 genes differentially expressed between mock-inoculated and inoculated plants. The combined ChIP-Seq and RNA-Seq analyses identified some key defense responsive genes (calmodulin, cytochrome p450, chitinase, DNA Pol II, and LRR) and transcription factors (WRKY, bZIP, MYB, HSFB3, GRAS, NAC, and NMRA) in bean-rust interaction. Differential methylation and acetylation affected a large proportion of stress-responsive genes including resistant (R) proteins, detoxifying enzymes, and genes involved in ion flux and cell death. The genes identified were functionally classified using Gene Ontology (GO) and EuKaryotic Orthologous Groups (KOGs). The Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis identified a putative pathway with ten key genes involved in plant-pathogen interactions. This first report of an integrated analysis of histone modifications and gene expression involved in the bean-rust interaction as reported here provides a comprehensive resource for other epigenomic regulation studies in non-model species under stress.

Gaining insight into soybean defense responses using functional genomics approaches

Soybean pathogens significantly impact yield, resulting in over $4 billion dollars in lost revenue annually in the United States. Despite the deployment of improved soybean cultivars, pathogens continue to evolve to evade plant defense responses. Thus, there is an urgent need to identify and characterize gene networks controlling defense responses to harmful pathogens. In this review, we focus on major advances that have been made in identifying the genes and gene networks regulating defense responses with an emphasis on soybean-pathogen interactions that have been amenable to gene function analyses using gene silencing technologies. Further we describe new research striving to identify genes involved in durable broad-spectrum resistance. Finally, we consider future prospects for functional genomic studies in soybean and demonstrate that understanding soybean disease and stress tolerance will be expedited at an unprecedented pace.

Genome-wide analysis of the gene families of resistance gene analogues in cotton and their response to Verticillium wilt

BACKGROUND

Gossypium raimondii is a Verticillium wilt-resistant cotton species whose genome encodes numerous disease resistance genes that play important roles in the defence against pathogens. However, the characteristics of resistance gene analogues (RGAs) and Verticillium dahliae response loci (VdRLs) have not been investigated on a global scale. In this study, the characteristics of RGA genes were systematically analysed using bioinformatics-driven methods. Moreover, the potential VdRLs involved in the defence response to Verticillium wilt were identified by RNA-seq and correlations with known resistance QTLs.

RESULTS

The G. raimondii genome encodes 1004 RGA genes, and most of these genes cluster in homology groups based on high levels of similarity. Interestingly, nearly half of the RGA genes occurred in 26 RGA-gene-rich clusters (Rgrcs). The homology analysis showed that sequence exchanges and tandem duplications frequently occurred within Rgrcs, and segmental duplications took place among the different Rgrcs. An RNA-seq analysis showed that the RGA genes play roles in cotton defence responses, forming 26 VdRLs inside in the Rgrcs after being inoculated with V. dahliae. A correlation analysis found that 12 VdRLs were adjacent to the known Verticillium wilt resistance QTLs, and that 5 were rich in NB-ARC domain-containing disease resistance genes.

CONCLUSIONS

The cotton genome contains numerous RGA genes, and nearly half of them are located in clusters, which evolved by sequence exchanges, tandem duplications and segmental duplications. In the Rgrcs, 26 loci were induced by the V. dahliae inoculation, and 12 are in the vicinity of known Verticillium wilt resistance QTLs.

Changes in the proteome of pad2-1, a glutathione depleted Arabidopsis mutant, during Pseudomonas syringae infection

The involvement of glutathione (GSH) in plant defense against pathogen invasion is an established fact. However, the molecular mechanism conferring this tolerance remains to be explored. Here, proteomic analysis of pad2-1, an Arabidopsis thaliana GSH-depleted mutant, in response to Pseudomonas syringae infection has been performed to explore the intricate position of GSH in defense against biotrophic pathogens. The pad2-1 mutant displayed severe susceptibility to P. syringae infection compared to the wild-type (Col-0) thus re-establishing a fundamental role of GSH in defense. Apart from general up-accumulation of energy metabolism-related protein-species in both infected Col-0 and pad2-1, several crucial defense-related protein-species were identified to be differentially accumulated. Leucine-rich repeat-receptor kinase (LRR-RK) and nucleotide-binding site-leucine-rich repeat resistance protein (NBS-LRR), known to play a pioneering role against pathogen attack, were only weakly up-accumulated in pad2-1 after infection. Transcriptional and post-transcriptional regulators like MYB-P1 and glycine-rich repeat RNA-binding protein (GRP) and several other stress-related protein-species like heat shock protein 17 (HSP17) and glutathione-S-transferase (GST) were also identified to be differentially regulated in pad2-1 and Col-0 in response to infection. Together, the present investigation reveals that the optimum GSH-level is essential for the efficient activation of plant defense signaling cascades thus conferring resistance to pathogen invasion.

Identification and distribution of the NBS-LRR gene family in the Cassava genome

BACKGROUND

Plant resistance genes (R genes) exist in large families and usually contain both a nucleotide-binding site domain and a leucine-rich repeat domain, denoted NBS-LRR. The genome sequence of cassava (Manihot esculenta) is a valuable resource for analysing the genomic organization of resistance genes in this crop.

RESULTS

With searches for Pfam domains and manual curation of the cassava gene annotations, we identified 228 NBS-LRR type genes and 99 partial NBS genes. These represent almost 1% of the total predicted genes and show high sequence similarity to proteins from other plant species. Furthermore, 34 contained an N-terminal toll/interleukin (TIR)-like domain, and 128 contained an N-terminal coiled-coil (CC) domain. 63% of the 327 R genes occurred in 39 clusters on the chromosomes. These clusters are mostly homogeneous, containing NBS-LRRs derived from a recent common ancestor.

CONCLUSIONS

This study provides insight into the evolution of NBS-LRR genes in the cassava genome; the phylogenetic and mapping information may aid efforts to further characterize the function of these predicted R genes.

Genome-wide SNP identification and characterization in two soybean cultivars with contrasting Mungbean Yellow Mosaic India Virus disease resistance traits

Mungbean yellow mosaic India virus (MYMIV) is a bipartite Geminivirus, which causes severe yield loss in soybean (Glycine max). Considering this, the present study was conducted to develop large-scale genome-wide single nucleotide polymorphism (SNP) markers and identify potential markers linked with known disease resistance loci for their effective use in genomics-assisted breeding to impart durable MYMIV tolerance. The whole-genome re-sequencing of MYMIV resistant cultivar 'UPSM-534' and susceptible Indian cultivar 'JS-335' was performed to identify high-quality SNPs and InDels (insertion and deletions). Approximately 234 and 255 million of 100-bp paired-end reads were generated from UPSM-534 and JS-335, respectively, which provided ~98% coverage of reference soybean genome. A total of 3083987 SNPs (1559556 in UPSM-534 and 1524431 in JS-335) and 562858 InDels (281958 in UPSM-534 and 280900 in JS-335) were identified. Of these, 1514 SNPs were found to be present in 564 candidate disease resistance genes. Among these, 829 non-synonymous and 671 synonymous SNPs were detected in 266 and 286 defence-related genes, respectively. Noteworthy, a non-synonymous SNP (in chromosome 18, named 18-1861613) at the 149th base-pair of LEUCINE-RICH REPEAT RECEPTOR-LIKE PROTEIN KINASE gene responsible for a G/C transversion [proline (CCC) to alanine(GCC)] was identified and validated in a set of 12 soybean cultivars. Taken together, the present study generated a large-scale genomic resource such as, SNPs and InDels at a genome-wide scale that will facilitate the dissection of various complex traits through construction of high-density linkage maps and fine mapping. In the present scenario, these markers can be effectively used to design high-density SNP arrays for their large-scale validation and high-throughput genotyping in diverse natural and mapping populations, which could accelerate genomics-assisted MYMIV disease resistance breeding in soybean.

CNL Disease Resistance Genes in Soybean and Their Evolutionary Divergence

Disease resistance genes (R-genes) encode proteins involved in detecting pathogen attack and activating downstream defense molecules. Recent availability of soybean genome sequences makes it possible to examine the diversity of gene families including disease-resistant genes. The objectives of this study were to identify coiled-coil NBS-LRR (= CNL) R-genes in soybean, infer their evolutionary relationships, and assess structural as well as functional divergence of the R-genes. Profile hidden Markov models were used for sequence identification and model-based maximum likelihood was used for phylogenetic analysis, and variation in chromosomal positioning, gene clustering, and functional divergence were assessed. We identified 188 soybean CNL genes nested into four clades consistent to their orthologs in Arabidopsis. Gene clustering analysis revealed the presence of 41 gene clusters located on 13 different chromosomes. Analyses of the K s-values and chromosomal positioning suggest duplication events occurring at varying timescales, and an extrapericentromeric positioning may have facilitated their rapid evolution. Each of the four CNL clades exhibited distinct patterns of gene expression. Phylogenetic analysis further supported the extrapericentromeric positioning effect on the divergence and retention of the CNL genes. The results are important for understanding the diversity and divergence of CNL genes in soybean, which would have implication in soybean crop improvement in future.

Whole-genome gene expression profiling revealed genes and pathways potentially involved in regulating interactions of soybean with cyst nematode (Heterodera glycines Ichinohe)

BACKGROUND

Soybean cyst nematode (SCN, Heterodera glycines Ichinohe) is the most devastating pathogen of soybean. Many gene expression profiling studies have been conducted to investigate the responses of soybean to the infection by this pathogen using primarily the first-generation soybean genome array that covered approximately 37,500 soybean transcripts. However, no study has been reported yet using the second-generation Affymetrix soybean whole-genome transcript array (Soybean WT array) that represents approximately 66,000 predicted soybean transcripts.

RESULTS

In the present work, the gene expression profiles of two soybean plant introductions (PIs) PI 437654 and PI 567516C (both resistant to multiple SCN HG Types) and cultivar Magellan (susceptible to SCN) were compared in the presence or absence of the SCN inoculum at 3 and 8 days post-inoculation using the Soybean WT array. Data analysis revealed that the two resistant soybean lines showed distinctive gene expression profiles from each other and from Magellan not only in response to the SCN inoculation, but also in the absence of SCN. Overall, 1,413 genes and many pathways were revealed to be differentially regulated. Among them, 297 genes were constitutively regulated in the two resistant lines (compared with Magellan) and 1,146 genes were responsive to the SCN inoculation in the three lines, with 30 genes regulated both constitutively and by SCN. In addition to the findings similar to those in the published work, many genes involved in ethylene, protein degradation, and phenylpropanoid pathways were also revealed differentially regulated in the present study. GC-rich elements (e.g., GCATGC) were found over-represented in the promoter regions of certain groups of genes. These have not been observed before, and could be new defense-responsive regulatory elements.

CONCLUSIONS

Different soybean lines showed different gene expression profiles in the presence and absence of the SCN inoculum. Both inducible and constitutive gene expression may contribute to resistance to multiple SCN HG Types in the resistant soybean PI lines. Ethylene, protein degradation, and phenylpropanoid pathways, as well as many other pathways reported previously, may play important roles in mediating the soybean-SCN interactions. The revealed genes, pathways, and promoter elements can be further explored to regulate or engineer soybean for resistance to SCN.

Fine mapping of the Asian soybean rust resistance gene Rpp2 from soybean PI 230970

KEY MESSAGE

Asian soybean rust (ASR) resistance gene Rpp2 has been fine mapped into a 188.1 kb interval on Glyma.Wm82.a2, which contains a series of plant resistance ( R ) genes. Asian soybean rust (ASR), caused by the fungus Phakopsora pachyrihizi Syd. & P. Syd., is a serious disease in major soybean [Glycine max (L.) Merr.] production countries worldwide and causes yield losses up to 75 %. Defining the exact chromosomal position of ASR resistance genes is critical for improving the effectiveness of marker-assisted selection (MAS) for resistance and for cloning these genes. The objective of this study was to fine map the ASR resistance gene Rpp2 from the plant introduction (PI) 230970. Rpp2 was previously mapped within a 12.9-cM interval on soybean chromosome 16. The fine mapping was initiated by identifying recombination events in F2 and F3 plants using simple sequence repeat (SSR) and single nucleotide polymorphism (SNP) markers that flank the gene. Seventeen recombinant plants were identified and then tested with additional genetic markers saturating the gene region to localize the positions of each recombination. The progeny of these selected plants were tested for resistance to ASR and with SSR markers resulting in the mapping of Rpp2 to a 188.1 kb interval on the Williams 82 reference genome (Glyma.Wm82.a2). Twelve genes including ten toll/interleukin-1 receptor (TIR)-nucleotide-binding site (NBS)-leucine-rich repeat (LRR) genes were predicted to exist in this interval on the Glyma.Wm82.a2.v1 gene model map. Eight of these ten genes were homologous to the Arabidopsis TIR-NBS-LRR gene AT5G17680.1. The identified SSR and SNP markers close to Rpp2 and the candidate gene information presented in this study will be significant resources for MAS and gene cloning research.

An atlas of soybean small RNAs identifies phased siRNAs from hundreds of coding genes

Small RNAs are ubiquitous, versatile repressors and include (1) microRNAs (miRNAs), processed from mRNA forming stem-loops; and (2) small interfering RNAs (siRNAs), the latter derived in plants by a process typically requiring an RNA-dependent RNA polymerase. We constructed and analyzed an expression atlas of soybean (Glycine max) small RNAs, identifying over 500 loci generating 21-nucleotide phased siRNAs (phasiRNAs; from PHAS loci), of which 483 overlapped annotated protein-coding genes. Via the integration of miRNAs with parallel analysis of RNA end (PARE) data, 20 miRNA triggers of 127 PHAS loci were detected. The primary class of PHAS loci (208 or 41% of the total) corresponded to NB-LRR genes; some of these small RNAs preferentially accumulate in nodules. Among the PHAS loci, novel representatives of TAS3 and noncanonical phasing patterns were also observed. A noncoding PHAS locus, triggered by miR4392, accumulated preferentially in anthers; the phasiRNAs are predicted to target transposable elements, with their peak abundance during soybean reproductive development. Thus, phasiRNAs show tremendous diversity in dicots. We identified novel miRNAs and assessed the veracity of soybean miRNAs registered in miRBase, substantially improving the soybean miRNA annotation, facilitating an improvement of miRBase annotations and identifying at high stringency novel miRNAs and their targets.

De novo assembly of soybean wild relatives for pan-genome analysis of diversity and agronomic traits

Wild relatives of crops are an important source of genetic diversity for agriculture, but their gene repertoire remains largely unexplored. We report the establishment and analysis of a pan-genome of Glycine soja, the wild relative of cultivated soybean Glycine max, by sequencing and de novo assembly of seven phylogenetically and geographically representative accessions. Intergenomic comparisons identified lineage-specific genes and genes with copy number variation or large-effect mutations, some of which show evidence of positive selection and may contribute to variation of agronomic traits such as biotic resistance, seed composition, flowering and maturity time, organ size and final biomass. Approximately 80% of the pan-genome was present in all seven accessions (core), whereas the rest was dispensable and exhibited greater variation than the core genome, perhaps reflecting a role in adaptation to diverse environments. This work will facilitate the harnessing of untapped genetic diversity from wild soybean for enhancement of elite cultivars.

The plasticity of NBS resistance genes in sorghum is driven by multiple evolutionary processes

BACKGROUND

Increased disease resistance is a key target of cereal breeding programs, with disease outbreaks continuing to threaten global food production, particularly in Africa. Of the disease resistance gene families, the nucleotide-binding site plus leucine-rich repeat (NBS-LRR) family is the most prevalent and ancient and is also one of the largest gene families known in plants. The sequence diversity in NBS-encoding genes was explored in sorghum, a critical food staple in Africa, with comparisons to rice and maize and with comparisons to fungal pathogen resistance QTL.

RESULTS

In sorghum, NBS-encoding genes had significantly higher diversity in comparison to non NBS-encoding genes and were significantly enriched in regions of the genome under purifying and balancing selection, both through domestication and improvement. Ancestral genes, pre-dating species divergence, were more abundant in regions with signatures of selection than in regions not under selection. Sorghum NBS-encoding genes were also significantly enriched in the regions of the genome containing fungal pathogen disease resistance QTL; with the diversity of the NBS-encoding genes influenced by the type of co-locating biotic stress resistance QTL.

CONCLUSIONS

NBS-encoding genes are under strong selection pressure in sorghum, through the contrasting evolutionary processes of purifying and balancing selection. Such contrasting evolutionary processes have impacted ancestral genes more than species-specific genes. Fungal disease resistance hot-spots in the genome, with resistance against multiple pathogens, provides further insight into the mechanisms that cereals use in the "arms race" with rapidly evolving pathogens in addition to providing plant breeders with selection targets for fast-tracking the development of high performing varieties with more durable pathogen resistance.

Genome-wide identification and expression analysis of NBS-encoding genes in Malus x domestica and expansion of NBS genes family in Rosaceae

Nucleotide binding site leucine-rich repeats (NBS-LRR) disease resistance proteins play an important role in plant defense against pathogen attack. A number of recent studies have been carried out to identify and characterize NBS-LRR gene families in many important plant species. In this study, we identified NBS-LRR gene family comprising of 1015 NBS-LRRs using highly stringent computational methods. These NBS-LRRs were characterized on the basis of conserved protein motifs, gene duplication events, chromosomal locations, phylogenetic relationships and digital gene expression analysis. Surprisingly, equal distribution of Toll/interleukin-1 receptor (TIR) and coiled coil (CC) (1 ∶ 1) was detected in apple while the unequal distribution was reported in majority of all other known plant genome studies. Prediction of gene duplication events intriguingly revealed that not only tandem duplication but also segmental duplication may equally be responsible for the expansion of the apple NBS-LRR gene family. Gene expression profiling using expressed sequence tags database of apple and quantitative real-time PCR (qRT-PCR) revealed the expression of these genes in wide range of tissues and disease conditions, respectively. Taken together, this study will provide a blueprint for future efforts towards improvement of disease resistance in apple.

Genome-wide association studies identifies seven major regions responsible for iron deficiency chlorosis in soybean (Glycine max)

Iron deficiency chlorosis (IDC) is a yield limiting problem in soybean (Glycine max (L.) Merr) production regions with calcareous soils. Genome-wide association study (GWAS) was performed using a high density SNP map to discover significant markers, QTL and candidate genes associated with IDC trait variation. A stepwise regression model included eight markers after considering LD between markers, and identified seven major effect QTL on seven chromosomes. Twelve candidate genes known to be associated with iron metabolism mapped near these QTL supporting the polygenic nature of IDC. A non-synonymous substitution with the highest significance in a major QTL region suggests soybean orthologs of FRE1 on Gm03 is a major gene responsible for trait variation. NAS3, a gene that encodes the enzyme nicotianamine synthase which synthesizes the iron chelator nicotianamine also maps to the same QTL region. Disease resistant genes also map to the major QTL, supporting the hypothesis that pathogens compete with the plant for Fe and increase iron deficiency. The markers and the allelic combinations identified here can be further used for marker assisted selection.

Association mapping for partial resistance to Phytophthora sojae in soybean (Glycine max (L.) Merr.)

Association mapping is a powerful high-resolution mapping tool for complex traits. The objective of this study was to identify QTLs for partial resistance to Phytophthora sojae. In this study, we evaluated a total of 214 soybean accessions by the hypocotyl inoculation method, and 175 were susceptible. The 175 susceptible accessions were then evaluated for P. sojae partial resistance using slant board assays. The 175 accessions were screened with 138 SSR markers that generated 730 SSR alleles. A subset of 495 SSR loci with minor allele frequency (MAF) ≥ 0.05 was used for association mapping by the Tassel general linear model (GLM) and mixed linear model (MLM) programmes. This soybean population could be divided into two subpopulations and no or weak relatedness was detected between pairwise accessions. Four SSR alleles, Satt634-133, Satt634-149, Sat_222-168 and Satt301-190, associated with partial resistance to P. sojae were detected by both GLM and MLM methods. Of these identified markers, one marker, Satt301, was located in regions where P. sojae resistance QTL have been previously mapped using linkage analysis. The identified markers will help to understand the genetic basis of partial resistance, and facilitate future marker-assistant selection aimed to improve resistance to P. sojae and reduce disease-related mortality in soybean.

Characterization of resistance gene analogues (RGAs) in apple (Malus × domestica Borkh.) and their evolutionary history of the Rosaceae family

The family of resistance gene analogues (RGAs) with a nucleotide-binding site (NBS) domain accounts for the largest number of disease resistance genes and is one of the largest gene families in plants. We have identified 868 RGAs in the genome of the apple (Malus × domestica Borkh.) cultivar ‘Golden Delicious’. This represents 1.51% of the total number of predicted genes for this cultivar. Several evolutionary features are pronounced in M. domestica, including a high fraction (80%) of RGAs occurring in clusters. This suggests frequent tandem duplication and ectopic translocation events. Of the identified RGAs, 56% are located preferentially on six chromosomes (Chr 2, 7, 8, 10, 11, and 15), and 25% are located on Chr 2. TIR-NBS and non-TIR-NBS classes of RGAs are primarily exclusive of different chromosomes, and 99% of non-TIR-NBS RGAs are located on Chr 11. A phylogenetic reconstruction was conducted to study the evolution of RGAs in the Rosaceae family. More than 1400 RGAs were identified in six species based on their NBS domain, and a neighbor-joining analysis was used to reconstruct the phylogenetic relationships among the protein sequences. Specific phylogenetic clades were found for RGAs of Malus, Fragaria, and Rosa, indicating genus-specific evolution of resistance genes. However, strikingly similar RGAs were shared in Malus, Pyrus, and Prunus, indicating high conservation of specific RGAs and suggesting a monophyletic origin of these three genera.

Genome-wide association studies identifies seven major regions responsible for iron deficiency chlorosis in soybean (Glycine max)

Iron deficiency chlorosis (IDC) is a yield limiting problem in soybean (Glycine max (L.) Merr) production regions with calcareous soils. Genome-wide association study (GWAS) was performed using a high density SNP map to discover significant markers, QTL and candidate genes associated with IDC trait variation. A stepwise regression model included eight markers after considering LD between markers, and identified seven major effect QTL on seven chromosomes. Twelve candidate genes known to be associated with iron metabolism mapped near these QTL supporting the polygenic nature of IDC. A non-synonymous substitution with the highest significance in a major QTL region suggests soybean orthologs of FRE1 on Gm03 is a major gene responsible for trait variation. NAS3, a gene that encodes the enzyme nicotianamine synthase which synthesizes the iron chelator nicotianamine also maps to the same QTL region. Disease resistant genes also map to the major QTL, supporting the hypothesis that pathogens compete with the plant for Fe and increase iron deficiency. The markers and the allelic combinations identified here can be further used for marker assisted selection.

Evolution and Conservation of Plant NLR Functions

In plants and animals, nucleotide-binding domain and leucine-rich repeats (NLR)-containing proteins play pivotal roles in innate immunity. Despite their similar biological functions and protein architecture, comparative genome-wide analyses of NLRs and genes encoding NLR-like proteins suggest that plant and animal NLRs have independently arisen in evolution. Furthermore, the demonstration of interfamily transfer of plant NLR functions from their original species to phylogenetically distant species implies evolutionary conservation of the underlying immune principle across plant taxonomy. In this review we discuss plant NLR evolution and summarize recent insights into plant NLR-signaling mechanisms, which might constitute evolutionarily conserved NLR-mediated immune mechanisms.

Co-segregation analysis and mapping of the anthracnose Co-10 and angular leaf spot Phg-ON disease-resistance genes in the common bean cultivar Ouro Negro

Anthracnose (ANT) and angular leaf spot (ALS) are devastating diseases of common bean (Phaseolus vulgaris L.). Ouro Negro is a highly productive common bean cultivar, which contains the Co-10 and Phg-ON genes for resistance to ANT and ALS, respectively. In this study, we performed a genetic co-segregation analysis of resistance to ANT and ALS using an F2 population from the Rudá × Ouro Negro cross and the F2:3 families from the AND 277 × Ouro Negro cross. Ouro Negro is resistant to races 7 and 73 of the ANT and race 63-39 of the ALS pathogens. Conversely, cultivars AND 277 and Rudá are susceptible to races 7 and 73 of ANT, respectively. Both cultivars are susceptible to race 63-39 of ALS. Co-segregation analysis revealed that Co-10 and Phg-ON were inherited together, conferring resistance to races 7 and 73 of ANT and race 63-39 of ALS. The Co-10 and Phg-ON genes were co-segregated and were tightly linked at a distance of 0.0 cM on chromosome Pv04. The molecular marker g2303 was linked to Co-10 and Phg-ON at a distance of 0.0 cM. Because of their physical linkage in a cis configuration, the Co-10 and Phg-ON resistance alleles are inherited together and can be monitored with great efficiency using g2303. The close linkage between the Co-10 and Phg-ON genes and prior evidence are consistent with the existence of a resistance gene cluster at one end of chromosome Pv04, which also contains the Co-3 locus and ANT resistance quantitative trait loci. These results will be very useful for breeding programs aimed at developing bean cultivars with ANT and ALS resistance using marker-assisted selection.

Systematic analysis and comparison of nucleotide-binding site disease resistance genes in a diploid cotton Gossypium raimondii

Plant disease resistance genes are a key component of defending plants from a range of pathogens. The majority of these resistance genes belong to the super-family that harbors a Nucleotide-binding site (NBS). A number of studies have focused on NBS-encoding genes in disease resistant breeding programs for diverse plants. However, little information has been reported with an emphasis on systematic analysis and comparison of NBS-encoding genes in cotton. To fill this gap of knowledge, in this study, we identified and investigated the NBS-encoding resistance genes in cotton using the whole genome sequence information of Gossypium raimondii. Totally, 355 NBS-encoding resistance genes were identified. Analyses of the conserved motifs and structural diversity showed that the most two distinct features for these genes are the high proportion of non-regular NBS genes and the high diversity of N-termini domains. Analyses of the physical locations and duplications of NBS-encoding genes showed that gene duplication of disease resistance genes could play an important role in cotton by leading to an increase in the functional diversity of the cotton NBS-encoding genes. Analyses of phylogenetic comparisons indicated that, in cotton, the NBS-encoding genes with TIR domain not only have their own evolution pattern different from those of genes without TIR domain, but also have their own species-specific pattern that differs from those of TIR genes in other plants. Analyses of the correlation between disease resistance QTL and NBS-encoding resistance genes showed that there could be more than half of the disease resistance QTL associated to the NBS-encoding genes in cotton, which agrees with previous studies establishing that more than half of plant resistance genes are NBS-encoding genes.