Identification of expressed resistance gene-like sequences by data mining in 454-derived transcriptomic sequences of common bean (Phaseolus vulgaris L.)

Exploring Chickpea Germplasm Diversity for Broadening the Genetic Base Utilizing Genomic Resourses

Legume crops provide significant nutrition to humans as a source of protein, omega-3 fatty acids as well as specific macro and micronutrients. Additionally, legumes improve the cropping environment by replenishing the soil nitrogen content. Chickpeas are the second most significant staple legume food crop worldwide behind dry bean which contains 17%–24% protein, 41%–51% carbohydrate, and other important essential minerals, vitamins, dietary fiber, folate, β-carotene, anti-oxidants, micronutrients (phosphorus, calcium, magnesium, iron, and zinc) as well as linoleic and oleic unsaturated fatty acids. Despite these advantages, legumes are far behind cereals in terms of genetic improvement mainly due to far less effort, the bottlenecks of the narrow genetic base, and several biotic and abiotic factors in the scenario of changing climatic conditions. Measures are now called for beyond conventional breeding practices to strategically broadening of narrow genetic base utilizing chickpea wild relatives and improvement of cultivars through advanced breeding approaches with a focus on high yield productivity, biotic and abiotic stresses including climate resilience, and enhanced nutritional values. Desirable donors having such multiple traits have been identified using core and mini core collections from the cultivated gene pool and wild relatives of Chickpea. Several methods have been developed to address cross-species fertilization obstacles and to aid in inter-specific hybridization and introgression of the target gene sequences from wild Cicer species. Additionally, recent advances in “Omics” sciences along with high-throughput and precise phenotyping tools have made it easier to identify genes that regulate traits of interest. Next-generation sequencing technologies, whole-genome sequencing, transcriptomics, and differential genes expression profiling along with a plethora of novel techniques like single nucleotide polymorphism exploiting high-density genotyping by sequencing assays, simple sequence repeat markers, diversity array technology platform, and whole-genome re-sequencing technique led to the identification and development of QTLs and high-density trait mapping of the global chickpea germplasm. These altogether have helped in broadening the narrow genetic base of chickpeas.

Genome-wide identification and characterization of NB-ARC resistant genes in wheat (Triticum aestivum L.) and their expression during leaf rust infection

NB-ARC domain-containing resistance genes from the wheat genome were identified, characterized and localized on chromosome arms that displayed differential yet positive response during incompatible and compatible leaf rust interactions. Wheat (Triticum aestivum L.) is an important cereal crop; however, its production is affected severely by numerous diseases including rusts. An efficient, cost-effective and ecologically viable approach to control pathogens is through host resistance. In wheat, high numbers of resistance loci are present but only few have been identified and cloned. A comprehensive analysis of the NB-ARC-containing genes in complete wheat genome was accomplished in this study. Complete NB-ARC encoding genes were mined from the Ensembl Plants database to predict 604 NB-ARC containing sequences using the HMM approach. Genome-wide analysis of orthologous clusters in the NB-ARC-containing sequences of wheat and other members of the Poaceae family revealed maximum homology with Oryza sativa indica and Brachypodium distachyon. The identification of overlap between orthologous clusters enabled the elucidation of the function and evolution of resistance proteins. The distributions of the NB-ARC domain-containing sequences were found to be balanced among the three wheat sub-genomes. Wheat chromosome arms 4AL and 7BL had the most NB-ARC domain-containing contigs. The spatio-temporal expression profiling studies exemplified the positive role of these genes in resistant and susceptible wheat plants during incompatible and compatible interaction in response to the leaf rust pathogen Puccinia triticina. Two NB-ARC domain-containing sequences were modelled in silico, cloned and sequenced to analyze their fine structures. The data obtained in this study will augment isolation, characterization and application NB-ARC resistance genes in marker-assisted selection based breeding programs for improving rust resistance in wheat.

Sequencing and characterization of leaf transcriptomes of six diploid Nicotiana species

BACKGROUND

Nicotiana belongs to the Solanaceae family that includes important crops such as tomato, potato, eggplant, and pepper. Nicotiana species are of worldwide economic importance and are important model plants for scientific research. Here we present the comparative analysis of the transcriptomes of six wild diploid Nicotiana species. Wild relatives provide an excellent study system for the analysis of the genetic basis for various traits, especially disease resistance.

RESULTS

Whole transcriptome sequencing (RNA-seq) was performed for leaves of six diploid Nicotiana species, i.e. Nicotiana glauca, Nicotiana noctiflora, Nicotiana cordifolia, Nicotiana knightiana, Nicotiana setchellii and Nicotiana tomentosiformis. For each species, 9.0-22.3 Gb high-quality clean data were generated, and 67,073-182,046 transcripts were assembled with lengths greater than 100 bp. Over 90 % of the ORFs in each species had significant similarity with proteins in the NCBI non-redundant protein sequence (NR) database. A total of 2491 homologs were identified and used to construct a phylogenetic tree from the respective transcriptomes in Nicotiana. Bioinformatic analysis identified resistance gene analogs, major transcription factor families, and alkaloid transporter genes linked to plant defense.

CONCLUSIONS

This is the first report on the leaf transcriptomes of six wild Nicotiana species by Illumina paired-end sequencing and de novo assembly without a reference genome. These sequence resources hopefully will provide an opportunity for identifying genes involved in plant defense and several important quality traits in wild Nicotiana and will accelerate functional genomic studies and genetic improvement efforts of Nicotiana or other important Solanaceae crops in the future.

Genome-wide investigation and expression analysis of AP2-ERF gene family in salt tolerant common bean

Apetala2-ethylene-responsive element binding factor (AP2-ERF) superfamily with common AP2-DNA binding domain have developmentally and physiologically important roles in plants. Since common bean genome project has been completed recently, it is possible to identify all of the AP2-ERF genes in the common bean genome. In this study, a comprehensive genome-wide in silico analysis identified 180 AP2-ERF superfamily genes in common bean (Phaseolus vulgaris). Based on the amino acid alignment and phylogenetic analyses, superfamily members were classified into four subfamilies: DREB (54), ERF (95), AP2 (27) and RAV (3), as well as one soloist. The physical and chemical characteristics of amino acids, interaction between AP2-ERF proteins, cis elements of promoter region of AP2-ERF genes and phylogenetic trees were predicted and analyzed. Additionally, expression levels of AP2-ERF genes were evaluated by in silico and qRT-PCR analyses. In silico micro-RNA target transcript analyses identified nearly all PvAP2-ERF genes as targets of by 44 different plant species' miRNAs were identified in this study. The most abundant target genes were PvAP2/ERF-20-25-62-78-113-173. miR156, miR172 and miR838 were the most important miRNAs found in targeting and BLAST analyses. Interactome analysis revealed that the transcription factor PvAP2-ERF78, an ortholog of Arabidopsis At2G28550, was potentially interacted with at least 15 proteins, indicating that it was very important in transcriptional regulation. Here we present the first study to identify and characterize the AP2-ERF transcription factors in common bean using whole-genome analysis, and the findings may serve as a references for future functional research on the transcription factors in common bean.

Genome-wide identification of the Phaseolus vulgaris sRNAome using small RNA and degradome sequencing

BACKGROUND

MiRNAs and phasiRNAs are negative regulators of gene expression. These small RNAs have been extensively studied in plant model species but only 10 mature microRNAs are present in miRBase version 21, the most used miRNA database, and no phasiRNAs have been identified for the model legume Phaseolus vulgaris. Thanks to the recent availability of the first version of the common bean genome, degradome data and small RNA libraries, we are able to present here a catalog of the microRNAs and phasiRNAs for this organism and, particularly, we suggest new protagonists in the symbiotic nodulation events.

RESULTS

We identified a set of 185 mature miRNAs, including 121 previously unpublished sequences, encoded by 307 precursors and distributed in 98 families. Degradome data allowed us to identify a total of 181 targets for these miRNAs. We reveal two regulatory networks involving conserved miRNAs: those known to play crucial roles in the establishment of nodules, and novel miRNAs present only in common bean, suggesting a specific role for these sequences. In addition, we identified 125 loci that potentially produce phased small RNAs, with 47 of them having all the characteristics of being triggered by a total of 31 miRNAs, including 14 new miRNAs identified in this study.

CONCLUSIONS

We provide here a set of new small RNAs that contribute to the broader knowledge of the sRNAome of Phaseolus vulgaris. Thanks to the identification of the miRNA targets from degradome analysis and the construction of regulatory networks between the mature microRNAs, we present here the probable functional regulation associated with the sRNAome and, particularly, in N2-fixing symbiotic nodules.

De novo transcriptome assembly of the wild relative of tea tree (Camellia taliensis) and comparative analysis with tea transcriptome identified putative genes associated with tea quality and stress response

Background

Camellia taliensis is one of the most important wild relatives of cultivated tea tree, C. sinensis. The species extensively occupies mountainous habitats representing a wide-range abiotic tolerance and biotic resistance and thus harbors valuable gene resources that may greatly benefit genetic improvement of cultivated tea tree. However, owning to a large genome size of ~3 Gb and structurally complex genome, there are fairly limited genetic information and particularly few genomic resources publicly available for this species. To better understand the key pathways determining tea flavor and enhance tea tree breeding programs, we performed a high-throughput transcriptome sequencing for C. taliensis.

Results

In this study, approximate 241.5 million high-quality paired-end reads, accounting for ~24 Gb of sequence data, were generated from tender shoots, young leaves, flower buds and flowers using Illumina HiSeq 2000 platform. De novo assembly with further processing and filtering yielded a set of 67,923 transcripts with an average length of 685 bp and an N50 of 995 bp. Based on sequence similarity searches against public databases, a total of 39,475 transcripts were annotated with gene descriptions, conserved protein domains or gene ontology (GO) terms. Candidate genes for major metabolic pathways involved in tea quality were identified and experimentally validated using RT-qPCR. Further gene expression profiles showed that they are differentially regulated at different developmental stages. To gain insights into the evolution of these genes, we aligned them to the previously cloned orthologous genes in C. sinensis, and found that considerable nucleotide variation within several genes involved in important secondary metabolic biosynthesis pathways, of which flavone synthase II gene (FNSII) is the most variable between these two species. Moreover, comparative analyses revealed that C. taliensis shows a remarkable expansion of LEA genes, compared to C. sinensis, which might contribute to the observed stronger stress resistance of C. taliensis.

Conclusion

We reported the first large-coverage transcriptome datasets for C. taliensis using the next-generation sequencing technology. Such comprehensive EST datasets provide an unprecedented opportunity for identifying genes involved in several major metabolic pathways and will accelerate functional genomic studies and genetic improvement efforts of tea trees in the future.

Electronic supplementary material

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

The effect on the transcriptome of Anemone coronaria following infection with rust (Tranzschelia discolor)

In order to understand plant/pathogen interaction, the transcriptome of uninfected (1S) and infected (2I) plant was sequenced at 3'end by the GS FLX 454 platform. De novo assembly of high-quality reads generated 27,231 contigs leaving 37,191 singletons in the 1S and 38,393 in the 2I libraries. ESTcalc tool suggested that 71% of the transcriptome had been captured, with 99% of the genes present being represented by at least one read. Unigene annotation showed that 50.5% of the predicted translation products shared significant homology with protein sequences in GenBank. In all 253 differential transcript abundance (DTAs) were in higher abundance and 52 in lower abundance in the 2I library. 128 higher abundance DTA genes were of fungal origin and 49 were clearly plant sequences. A tBLASTn-based search of the sequences using as query the full length predicted polypeptide product of 50 R genes identified 16 R gene products. Only one R gene (PGIP) was up-regulated. The response of the plant to fungal invasion included the up-regulation of several pathogenesis related protein (PR) genes involved in JA signaling and other genes associated with defense response and down regulation of cell wall associated genes, non-race-specific disease resistance1 (NDR1) and other genes like myb, presqualene diphosphate phosphatase (PSDPase), a UDP-glycosyltransferase 74E2-like (UGT). The DTA genes identified here should provide a basis for understanding the A. coronaria/T. discolor interaction and leads for biotechnology-based disease resistance breeding.

Common bean reaction to angular leaf spot comprises transcriptional modulation of genes in the ALS10.1 QTL

Genetic resistance of common bean (Phaseolus vulgaris L.) against angular leaf spot (ALS), caused by the fungus Pseudocercospora griseola, is conferred by quantitative trait loci (QTL). In this study, we determined the gene content of the major QTL ALS10.1 located at the end of chromosome Pv10, and identified those that are responsive to ALS infection in resistant (CAL 143) and susceptible (IAC-UNA) genotypes. Based on the current version of the common bean reference genome, the ALS10.1 core region contains 323 genes. Gene Ontology (GO) analysis of these coding sequences revealed the presence of genes involved in signal perception and transduction, programmed cell death (PCD), and defense responses. Two putative R gene clusters were found at ALS10.1 containing evolutionary related coding sequences. Among them, the Phvul.010G025700 was consistently up-regulated in the infected IAC-UNA suggesting its contribution to plant susceptibility to the fungus. We identified six other genes that were regulated during common bean response to P. griseola; three of them might be negative regulators of immunity as they showed opposite expression patterns during resistant and susceptible reactions at the initial phase of fungal infection. Taken together, these findings suggest that common bean reaction to P. griseola involves transcriptional modulation of defense genes in the ALS10.1 locus, contributing to resistance or susceptibility depending on the plant-pathogen interaction.

Common bean reaction to angular leaf spot comprises transcriptional modulation of genes in the ALS10.1 QTL

Genetic resistance of common bean (Phaseolus vulgaris L.) against angular leaf spot (ALS), caused by the fungus Pseudocercospora griseola, is conferred by quantitative trait loci (QTL). In this study, we determined the gene content of the major QTL ALS10.1 located at the end of chromosome Pv10, and identified those that are responsive to ALS infection in resistant (CAL 143) and susceptible (IAC-UNA) genotypes. Based on the current version of the common bean reference genome, the ALS10.1 core region contains 323 genes. Gene Ontology (GO) analysis of these coding sequences revealed the presence of genes involved in signal perception and transduction, programmed cell death (PCD), and defense responses. Two putative R gene clusters were found at ALS10.1 containing evolutionary related coding sequences. Among them, the Phvul.010G025700 was consistently up-regulated in the infected IAC-UNA suggesting its contribution to plant susceptibility to the fungus. We identified six other genes that were regulated during common bean response to P. griseola; three of them might be negative regulators of immunity as they showed opposite expression patterns during resistant and susceptible reactions at the initial phase of fungal infection. Taken together, these findings suggest that common bean reaction to P. griseola involves transcriptional modulation of defense genes in the ALS10.1 locus, contributing to resistance or susceptibility depending on the plant-pathogen interaction.

Tree immunity: growing old without antibodies

Perennial plants need to cope with changing environments and pathogens over their lifespan. Infections are compartmentalised by localised physiological responses, and multiple apical meristems enable repair and regrowth, but genes are another crucial component in the perception and response to pathogens. In this opinion article we suggest that the mechanism for dynamic pathogen-specific recognition in long-lived plants could be explained by extending our current understanding of plant defence genes. We propose that, in addition to physiological responses, tree defence uses a three-pronged genomic approach involving: (i) gene numbers, (ii) genomic architecture, and (iii) mutation loads accumulated over long lifespans.

Transcriptome analysis of salt tolerant common bean (Phaseolus vulgaris L.) under saline conditions

Salinity is one of the important abiotic stress factors that limit crop production. Common bean, Phaseolus vulgaris L., a major protein source in developing countries, is highly affected by soil salinity and the information on genes that play a role in salt tolerance is scarce. We aimed to identify differentially expressed genes (DEGs) and related pathways by comprehensive analysis of transcriptomes of both root and leaf tissues of the tolerant genotype grown under saline and control conditions in hydroponic system. We have generated a total of 158 million high-quality reads which were assembled into 83,774 all-unigenes with a mean length of 813 bp and N50 of 1,449 bp. Among the all-unigenes, 58,171 were assigned with Nr annotations after homology analyses. It was revealed that 6,422 and 4,555 all-unigenes were differentially expressed upon salt stress in leaf and root tissues respectively. Validation of the RNA-seq quantifications (RPKM values) was performed by qRT-PCR (Quantitative Reverse Transcription PCR) analyses. Enrichment analyses of DEGs based on GO and KEGG databases have shown that both leaf and root tissues regulate energy metabolism, transmembrane transport activity, and secondary metabolites to cope with salinity. A total of 2,678 putative common bean transcription factors were identified and classified under 59 transcription factor families; among them 441 were salt responsive. The data generated in this study will help in understanding the fundamentals of salt tolerance in common bean and will provide resources for functional genomic studies.

Identification of expressed resistance gene analogs from peanut (Arachis hypogaea L.) expressed sequence tags

Low genetic diversity makes peanut (Arachis hypogaea L.) very vulnerable to plant pathogens, causing severe yield loss and reduced seed quality. Several hundred partial genomic DNA sequences as nucleotide-binding-site leucine-rich repeat (NBS-LRR) resistance genes (R) have been identified, but a small portion with expressed transcripts has been found. We aimed to identify resistance gene analogs (RGAs) from peanut expressed sequence tags (ESTs) and to develop polymorphic markers. The protein sequences of 54 known R genes were used to identify homologs from peanut ESTs from public databases. A total of 1,053 ESTs corresponding to six different classes of known R genes were recovered, and assembled 156 contigs and 229 singletons as peanut-expressed RGAs. There were 69 that encoded for NBS-LRR proteins, 191 that encoded for protein kinases, 82 that encoded for LRR-PK/transmembrane proteins, 28 that encoded for Toxin reductases, 11 that encoded for LRR-domain containing proteins and four that encoded for TM-domain containing proteins. Twenty-eight simple sequence repeats (SSRs) were identified from 25 peanut expressed RGAs. One SSR polymorphic marker (RGA121) was identified. Two polymerase chain reaction-based markers (Ahsw-1 and Ahsw-2) developed from RGA013 were homologous to the Tomato Spotted Wilt Virus (TSWV) resistance gene. All three markers were mapped on the same linkage group AhIV. These expressed RGAs are the source for RGA-tagged marker development and identification of peanut resistance genes.

Transcriptome Profiling Identifies Candidate Genes Associated with the Accumulation of Distinct Sulfur γ-Glutamyl Dipeptides in Phaseolus vulgaris and Vigna mungo Seeds

Common bean (Phaseolus vulgaris) and black gram (Vigna mungo) accumulate γ-Glutamyl-S-methylcysteine and γ-Glutamyl-methionine in seed, respectively. Transcripts were profiled by 454 pyrosequencing data at a similar developmental stage coinciding with the beginning of the accumulation of these metabolites. Expressed sequence tags were assembled into Unigenes, which were assigned to specific genes in the early release chromosomal assembly of the P. vulgaris genome. Genes involved in multiple sulfur metabolic processes were expressed in both species. Expression of Sultr3 members was predominant in P. vulgaris, whereas expression of Sultr5 members predominated in V. mungo. Expression of the cytosolic SERAT1;1 and -1;2 was approximately fourfold higher in P. vulgaris while expression of the plastidic SERAT2;1 was twofold higher in V. mungo. Among BSAS family members, BSAS4;1, encoding a cytosolic cysteine desulfhydrase, and BSAS1;1, encoding a cytosolic O-acetylserine sulphydrylase were most highly expressed in both species. This was followed by BSAS3;1 encoding a plastidic β-cyanoalanine synthase which was more highly expressed by 10-fold in P. vulgaris. The data identify BSAS3;1 as a candidate enzyme for the biosynthesis of S-methylcysteine through the use of methanethiol as substrate instead of cyanide. Expression of GLC1 would provide a complete sequence leading to the biosynthesis of γ-Glutamyl-S-methylcysteine in plastids. The detection of S-methylhomoglutathione in P. vulgaris suggested that homoglutathione synthetase may accept, to some extent, γ-Glutamyl-S-methylcysteine as substrate, which might lead to the formation of S-methylated phytochelatins. In conclusion, 454 sequencing was effective at revealing differences in the expression of sulfur metabolic genes, providing information on candidate genes for the biosynthesis of distinct sulfur amino acid γ-Glutamyl dipeptides between P. vulgaris and V. mungo.

Analysis of the leaf transcriptome of Musa acuminata during interaction with Mycosphaerella musicola: gene assembly, annotation and marker development

BACKGROUND

Although banana (Musa sp.) is an important edible crop, contributing towards poverty alleviation and food security, limited transcriptome datasets are available for use in accelerated molecular-based breeding in this genus. 454 GS-FLX Titanium technology was employed to determine the sequence of gene transcripts in genotypes of Musa acuminata ssp. burmannicoides Calcutta 4 and M. acuminata subgroup Cavendish cv. Grande Naine, contrasting in resistance to the fungal pathogen Mycosphaerella musicola, causal organism of Sigatoka leaf spot disease. To enrich for transcripts under biotic stress responses, full length-enriched cDNA libraries were prepared from whole plant leaf materials, both uninfected and artificially challenged with pathogen conidiospores.

RESULTS

The study generated 846,762 high quality sequence reads, with an average length of 334 bp and totalling 283 Mbp. De novo assembly generated 36,384 and 35,269 unigene sequences for M. acuminata Calcutta 4 and Cavendish Grande Naine, respectively. A total of 64.4% of the unigenes were annotated through Basic Local Alignment Search Tool (BLAST) similarity analyses against public databases.Assembled sequences were functionally mapped to Gene Ontology (GO) terms, with unigene functions covering a diverse range of molecular functions, biological processes and cellular components. Genes from a number of defense-related pathways were observed in transcripts from each cDNA library. Over 99% of contig unigenes mapped to exon regions in the reference M. acuminata DH Pahang whole genome sequence. A total of 4068 genic-SSR loci were identified in Calcutta 4 and 4095 in Cavendish Grande Naine. A subset of 95 potential defense-related gene-derived simple sequence repeat (SSR) loci were validated for specific amplification and polymorphism across M. acuminata accessions. Fourteen loci were polymorphic, with alleles per polymorphic locus ranging from 3 to 8 and polymorphism information content ranging from 0.34 to 0.82.

CONCLUSIONS

A large set of unigenes were characterized in this study for both M. acuminata Calcutta 4 and Cavendish Grande Naine, increasing the number of public domain Musa ESTs. This transcriptome is an invaluable resource for furthering our understanding of biological processes elicited during biotic stresses in Musa. Gene-based markers will facilitate molecular breeding strategies, forming the basis of genetic linkage mapping and analysis of quantitative trait loci.

Isolation and characterization of nucleotide-binding site resistance gene homologues in common bean (Phaseolus vulgaris)

Common bean production is constrained by many fungal, viral, and bacterial pathogens. Thus, the identification of resistance (R) genes is an important focal point of common bean research. The main goal of our study was to identify resistance gene homologues (RGH) in the crop, using degenerate primers designed from conserved sequences in the nucleotide-binding site (NBS) domains of R-genes from the model legume Medicago truncatula. Total DNA of the Andean common bean genotype G19833 was used for amplification of over 500 primer combinations. Sequencing of amplicons showed that 403 cloned fragments had uninterrupted open reading frames and were considered representative of functional RGH genes. The sequences were grouped at two levels of nucleotide identity (90 and 80%) and representative sequences of each group were used for phylogenetic analyses. The RGH sequence diversity of common bean was divided into TIR and non-TIR families, each with different clusters. The TIR sequences grouped into 14 clades while non-TIR sequences grouped into seven clades. Pairwise comparisons showed purifying selection, although some sequences may have been the result of diversifying selection. Knowledge about RGH genes in common bean can allow the design of molecular markers for pyramiding of resistance genes against various pathogens.