PRGdb 2.0: towards a community-based database model for the analysis of R-genes in plants

De novo transcriptome assembly of Dalbergia sissoo Roxb. (Fabaceae) under Botryodiplodia theobromae-induced dieback disease

Dalbergia sissoo Roxb. (Shisham) is a timber-producing species of economic, cultural, and medicinal importance in the Indian subcontinent. In the past few decades, Shisham's dieback disease caused by the fungus Botryodiplodia theobromae has become an evolving issue in the subcontinent endangering its survival. To gain insights into this issue, a standard transcriptome assembly was deployed to assess the response of D. sissoo at the transcriptomic level under the stress of B. theobromae infection. For RNA isolation, the control and infected leaf tissue samples were taken from 1-year-old greenhouse-grown D. sissoo plants after 20 days of stem-base spore inoculation. cDNA synthesis was performed from these freshly isolated RNA samples that were then sent for sequencing. About 18.14 Gb (Giga base) of data was generated using the BGISEQ-500 sequencing platform. In terms of Unigenes, 513,821 were identified after a combined assembly of all samples and then filtering the abundance. The total length of Unigenes, their average length, N50, and GC-content were 310,523,693 bp, 604 bp, 1,101 bp, and 39.95% respectively. The Unigenes were annotated using 7 functional databases i.e., 200,355 (NR: 38.99%), 164,973 (NT: 32.11%), 123,733 (Swissprot: 24.08%), 142,580 (KOG: 27.75%), 139,588 (KEGG: 27.17%), 99,752 (GO: 19.41%), and 137,281 (InterPro: 26.72%). Furthermore, the Transdecoder detected 115,762 CDS. In terms of SSR (Simple Sequence Repeat) markers, 62,863 of them were distributed on 51,508 Unigenes and on the predicted 4673 TF (Transcription Factor) coding Unigenes. A total of 16,018 up- and 19,530 down-regulated Differentially Expressed Genes (DEGs) were also identified. Moreover, the Plant Resistance Genes (PRGs) had a count of 9230. We are hopeful that in the future, these identified Unigenes, SSR markers, DEGs and PRGs will provide the prerequisites for managing Shisham dieback disease, its breeding, and in tree improvement programs.

Transcriptome Analysis Reveals the Defense Mechanism of Cotton against Verticillium dahliae Induced by Hypovirulent Fungus Gibellulopsis nigrescens CEF08111

Verticillium wilt is a kind of plant vascular disease caused by the soilborne fungus Verticillium dahliae, which severely limits cotton production. Our previous studies showed that the endophytic fungus Gibellulopsis nigrescens CEF08111 can effectively control Verticillium wilt and induce a defense response in cotton plants. However, the comprehensive molecular mechanism governing this response is not yet clear. To study the signaling mechanism induced by strain CEF08111, the transcriptome of cotton seedlings pretreated with CEF08111 was sequenced. The results revealed 249, 3559 and 33 differentially expressed genes (DEGs) at 3, 12 and 48 h post inoculation with CEF08111, respectively. At 12 h post inoculation with CEF08111, Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis indicated that the DEGs were enriched mainly in the plant−pathogen interaction, mitogen-activated protein kinase (MAPK) signaling pathway-plant, and plant hormone signal transduction pathways. Gene ontology (GO) analysis revealed that these DEGs were enriched mainly in the following terms: response to external stimulus, systemic acquired resistance, kinase activity, phosphotransferase activity, xyloglucan: xyloglucosyl transferase activity, xyloglucan metabolic process, cell wall polysaccharide metabolic process and hemicellulose metabolic process. Moreover, many genes, such as calcium-dependent protein kinase (CDPK), flagellin-sensing 2 (FLS2), resistance to Pseudomonas syringae pv. maculicola 1(RPM1) and myelocytomatosis protein 2 (MYC2), that regulate crucial points in defense-related pathways were identified and may contribute to V. dahliae resistance in cotton. Seven DEGs of the pathway phenylpropanoid biosynthesis were identified by weighted gene co-expression network analysis (WGCNA), and these genes are related to lignin synthesis. The above genes were compared and analyzed, a total of 710 candidate genes that may be related to the resistance of cotton to Verticillium wilt were identified. These results provide a basis for understanding the molecular mechanism by which the biocontrol fungus CEF08111 increases the resistance of cotton to Verticillium wilt.

Bioinformatics approaches and applications in plant biotechnology

BACKGROUND

In recent years, major advance in molecular biology and genomic technologies have led to an exponential growth in biological information. As the deluge of genomic information, there is a parallel growth in the demands of tools in the storage and management of data, and the development of software for analysis, visualization, modelling, and prediction of large data set.

MAIN BODY

Particularly in plant biotechnology, the amount of information has multiplied exponentially with a large number of databases available from many individual plant species. Efficient bioinformatics tools and methodologies are also developed to allow rapid genome sequence and the study of plant genome in the 'omics' approach. This review focuses on the various bioinformatic applications in plant biotechnology, and their advantages in improving the outcome in agriculture. The challenges or limitations faced in plant biotechnology in the aspect of bioinformatics approach that explained the low progression in plant genomics than in animal genomics are also reviewed and assessed.

CONCLUSION

There is a critical need for effective bioinformatic tools, which are able to provide longer reads with unbiased coverage in order to overcome the complexity of the plant's genome. The advancement in bioinformatics is not only beneficial to the field of plant biotechnology and agriculture sectors, but will also contribute enormously to the future of humanity.

ResCap: plant resistance gene prediction and probe generation pipeline for resistance gene sequence capture

Summary

The discovery of novel resistance genes (R-genes) is an important component in disease resistance breeding. Nevertheless, R-gene identification from wild species and close relatives of plants is not only a difficult but also a cumbersome process. In this study, ResCap, a support vector machine-based high-throughput R-gene prediction and probe generation pipeline has been developed to generate probes from genomic datasets. ResCap contains two integral modules. The first module identifies the R-genes and R-gene like sequences under four categories containing different domains such as TIR-NBS-LRR (TNL), CC-NBS-LRR (CNL), Receptor-like kinase (RLK) and Receptor-like proteins (RLPs). The second module generates probes from extracted nucleotide sequences of resistance genes to conduct sequence capture (SeqCap) experiments. For the validation of ResCap pipeline, ResCap generated probes were synthesized and a sequence capture experiment was performed to capture expressed resistance genes among six spring barley genotypes. The developed ResCap pipeline in combination with the performed sequence capture experiment has shown to increase precision of R-gene identification while simultaneously allowing rapid gene validation including non-sequenced plants.

Availability and implementation

The ResCap pipeline is available at http://rescap.ltj.slu.se/ResCap/.

Contact

sandeep.kushwaha@slu.se or sandeep@niab.org.in.

Supplementary information

Supplementary materials are available at Bioinformatics Advances online.

Transcriptome Analysis of a Cotton Cultivar Provides Insights into the Differentially Expressed Genes Underlying Heightened Resistance to the Devastating Verticillium Wilt

Cotton is an important economic crop worldwide. Verticillium wilt (VW) caused by Verticillium dahliae (V. dahliae) is a serious disease in cotton, resulting in massive yield losses and decline of fiber quality. Breeding resistant cotton cultivars is an efficient but elaborate method to improve the resistance of cotton against V. dahliae infection. However, the functional mechanism of several excellent VW resistant cotton cultivars is poorly understood at present. In our current study, we carried out RNA-seq to discover the differentially expressed genes (DEGs) in the roots of susceptible cotton Gossypium hirsutum cultivar Junmian 1 (J1) and resistant cotton G.hirsutum cultivar Liaomian 38 (L38) upon Vd991 inoculation at two time points compared with the mock inoculated control plants. The potential function of DEGs uniquely expressed in J1 and L38 was also analyzed by GO enrichment and KEGG pathway associations. Most DEGs were assigned to resistance-related functions. In addition, resistance gene analogues (RGAs) were identified and analyzed for their role in the heightened resistance of the L38 cultivar against the devastating Vd991. In summary, we analyzed the regulatory network of genes in the resistant cotton cultivar L38 during V. dahliae infection, providing a novel and comprehensive insight into VW resistance in cotton.

Genetic Loci Associated with Resistance to Zucchini Yellow Mosaic Virus in Squash

Zucchini Yellow Mosaic Virus (ZYMV) is an aphid-transmitted potyvirus that causes severe yield losses in squash (Cucurbita moschata) production worldwide. Development of resistant cultivars using traditional breeding approaches relies on rigorous and resource-intensive phenotypic assays. QTL-seq, a whole genome re-sequencing based bulked segregant analysis, is a powerful tool for mapping quantitative trait loci (QTL) in crop plants. In the current study, the QTL-seq approach was used to identify genetic loci associated with ZYMV resistance in an F₂ population (n = 174) derived from a cross between Nigerian Local (resistant) and Butterbush (susceptible). Whole genome re-sequencing of the parents and bulks of resistant and susceptible F₂ progeny revealed a mapping rate between 94.04% and 98.76%, and a final effective mapping depth ranging from 81.77 to 101.73 across samples. QTL-seq analysis identified four QTLs significantly (p < 0.05) associated with ZYMV resistance on chromosome 2 (QtlZYMV-C02), 4 (QtlZYMV-C04), 8 (QtlZYMV-C08) and 20 (QtlZYMV-C20). Seven markers within the QTL intervals were tested for association with ZYMV resistance in the entire F₂ population. For QtlZYMV-C08, one single nucleotide polymorphism (SNP) marker (KASP-6) was found to be significantly (p < 0.05) associated with ZYMV resistance, while two SNPs (KASP-1 and KASP-3) and an indel (Indel-2) marker were linked to resistance within QtlZYMV-C20. KASP-3 and KASP-6 are non-synonymous SNPs leading to amino acid substitutions in candidate disease resistant gene homologs on chromosomes 20 (CmoCh20G003040.1) and 8 (CmoCh08G007140.1), respectively. Identification of QTL and SNP markers associated with ZYMV resistance will facilitate marker-assisted selection for ZYMV resistance in squash.

Bioinformatic tools support decision-making in plant disease management

Food loss due to pathogens is a major concern in agriculture, requiring the need for advanced disease detection and prevention measures to minimize pathogen damage to plants. Novel bioinformatic tools have opened doors for the low-cost rapid identification of pathogens and prevention of disease. The number of these tools is growing fast and a comprehensive and comparative summary of these resources is currently lacking. Here, we review all current bioinformatic tools used to identify the mechanisms of pathogen pathogenicity, plant resistance protein identification, and the detection and treatment of plant disease. We compare functionality, data volume, data sources, performance, and applicability of all tools to provide a comprehensive toolbox for researchers in plant disease management.

Structural and Functional Genomics of the Resistance of Cacao to Phytophthora palmivora

Black pod disease, caused by Phytophthora spp., is one of the main diseases that attack cocoa plantations. This study validated, by association mapping, 29 SSR molecular markers flanking to QTL (Quantitative Trait Loci) associated with Phytophthora palmivora Butler (Butler) (PP) resistance, in three local ancient varieties of the Bahia (Comum, Pará, and Maranhão), varieties that have a high potential in the production of gourmet chocolate. Four SSR loci associated with resistance to PP were detected, two on chromosome 8, explaining 7.43% and 3.72% of the Phenotypic Variation (%PV), one on chromosome 2 explaining 2.71%PV and one on chromosome 3 explaining 1.93%PV. A functional domains-based annotation was carried out, in two Theobroma cacao (CRIOLLO and MATINA) reference genomes, of 20 QTL regions associated with cocoa resistance to the pathogen. It was identified 164 (genome CRIOLLO) and 160 (genome MATINA) candidate genes, hypothetically involved in the recognition and activation of responses in the interaction with the pathogen. Genomic regions rich in genes with Coiled-coils (CC), nucleotide binding sites (NBS) and Leucine-rich repeat (LRR) domains were identified on chromosomes 1, 3, 6, 8, and 10, likewise, regions rich in Receptor-like Kinase domain (RLK) and Ginkbilobin2 (GNK2) domains were identified in chromosomes 4 and 6.

Evolution of NLR resistance genes with noncanonical N-terminal domains in wild tomato species

Nucleotide-binding and leucine-rich repeat immune receptors (NLRs) provide resistance against diverse pathogens. To create comparative NLR resources, we conducted resistance gene enrichment sequencing (RenSeq) with single-molecule real-time sequencing of PacBio for 18 accessions in Solanaceae, including 15 accessions of five wild tomato species. We investigated the evolution of a class of NLRs, CNLs with extended N-terminal sequences previously named Solanaceae Domain. Through comparative genomic analysis, we revealed that the extended CNLs (exCNLs) anciently emerged in the most recent common ancestor between Asterids and Amaranthaceae, far predating the Solanaceae family. In tomatoes, the exCNLs display exceptional modes of evolution in a clade-specific manner. In the clade G3, exCNLs have substantially elongated their N-termini through tandem duplications of exon segments. In the clade G1, exCNLs have evolved through recent proliferation and sequence diversification. In the clade G6, an ancestral exCNL has lost its N-terminal domains in the course of evolution. Our study provides high-quality NLR gene models for close relatives of domesticated tomatoes that can serve as a useful resource for breeding and molecular engineering for disease resistance. Our findings regarding the exCNLs offer unique backgrounds and insights for future functional studies of the NLRs.

Gene expression profiling identifies pathways involved in seed maturation of Jatropha curcas

Background

Jatropha curcas, a tropical shrub, is a promising biofuel crop, which produces seeds with high content of oil and protein. To better understand the maturation process of J. curcas seeds and to improve its agronomic performance, a two-step approach was performed in six different maturation stages of seeds: 1) generation of the entire transcriptome of J. curcas seeds using 454-Roche sequencing of a cDNA library, 2) comparison of transcriptional expression levels using a custom Agilent 8x60K oligonucleotide microarray.

Results

A total of 793,875 high-quality reads were assembled into 19,382 unique full-length contigs, of which 13,507 could be annotated with Gene Ontology (GO) terms. Microarray data analysis identified 9111 probes (out of 57,842 probes), which were differentially expressed between the six maturation stages. The expression results were validated for 75 selected transcripts based on expression levels, predicted function, pathway, and length.

Result from cluster analyses showed that transcripts associated with fatty acid, flavonoid, and phenylpropanoid biosynthesis were over-represented in the early stages, while those of lipid storage were over-represented in the late stages. Expression analyses of different maturation stages of J. curcas seed showed that most changes in transcript abundance occurred between the two last stages, suggesting that the timing of metabolic pathways during seed maturation in J. curcas occurs in late stages. The co-expression results showed that the hubs (CB5-D, CDR1, TT8, DFR, HVA22) with the highest number of edges, associated with fatty acid and flavonoid biosynthesis, are showing a decrease in their expression during seed maturation. Furthermore, seed development and hormone pathways are significantly well connected.

Conclusion

The obtained results revealed differentially expressed sequences (DESs) regulating important pathways related to seed maturation, which could contribute to the understanding of the complex regulatory network during seed maturation with the focus on lipid, flavonoid and phenylpropanoid biosynthesis. This study provides detailed information on transcriptional changes during J. curcas seed maturation and provides a starting point for a genomic survey of seed quality traits. The results highlighted specific genes and processes relevant to the molecular mechanisms involved in Jatropha seed maturation. These data can also be utilized regarding other Euphorbiaceae species.

Transcriptome analysis reveals the defense mechanism of cotton against Verticillium dahliae in the presence of the biocontrol fungus Chaetomium globosum CEF-082

BACKGROUND

Verticillium wilt of cotton is a serious soil-borne disease that causes a substantial reduction in cotton yields. A previous study showed that the endophytic fungus Chaetomium globosum CEF-082 could control Verticillium wilt of cotton, and induce a defense response in cotton plants. However, the comprehensive molecular mechanism governing this response is not yet clear.

RESULTS

To study the signalling mechanism induced by CEF-082, the transcriptome of cotton seedlings pretreated with CEF-082 was sequenced. The results revealed 5638 DEGs at 24 h post inoculation with CEF-082, and 2921 and 2153 DEGs at 12 and 48 h post inoculation with Verticillium dahliae, respectively. At 24 h post inoculation with CEF-082, KEGG enrichment analysis indicated that the DEGs were enriched mainly in the plant-pathogen interaction, MAPK signalling pathway-plant, flavonoid biosynthesis, and phenylpropanoid biosynthesis pathways. There were 1209 DEGs specifically induced only in cotton plants inoculated with V. dahliae in the presence of the biocontrol fungus CEF-082, and not when cotton plants were only inoculated with V. dahliae. GO analysis revealed that these DEGs were enriched mainly in the following terms: ROS metabolic process, H2O2 metabolic process, defense response, superoxide dismutase activity, and antioxidant activity. Moreover, many genes, such as ERF, CNGC, FLS2, MYB, GST and CML, that regulate crucial points in defense-related pathways were identified and may contribute to V. dahliae resistance in cotton. These results provide a basis for understanding the molecular mechanism by which the biocontrol fungus CEF-082 increases the resistance of cotton to Verticillium wilt.

CONCLUSIONS

The results of this study showed that CEF-082 could regulate multiple metabolic pathways in cotton. After treatment with V. dahliae, the defense response of cotton plants preinoculated with CEF-082 was strengthened.

Transposons played a major role in the diversification between the closely related almond and peach genomes: results from the almond genome sequence

We sequenced the genome of the highly heterozygous almond Prunus dulcis cv. Texas combining short- and long-read sequencing. We obtained a genome assembly totaling 227.6 Mb of the estimated almond genome size of 238 Mb, of which 91% is anchored to eight pseudomolecules corresponding to its haploid chromosome complement, and annotated 27 969 protein-coding genes and 6747 non-coding transcripts. By phylogenomic comparison with the genomes of 16 additional close and distant species we estimated that almond and peach (Prunus persica) diverged around 5.88 million years ago. These two genomes are highly syntenic and show a high degree of sequence conservation (20 nucleotide substitutions per kb). However, they also exhibit a high number of presence/absence variants, many attributable to the movement of transposable elements (TEs). Transposable elements have generated an important number of presence/absence variants between almond and peach, and we show that the recent history of TE movement seems markedly different between them. Transposable elements may also be at the origin of important phenotypic differences between both species, and in particular for the sweet kernel phenotype, a key agronomic and domestication character for almond. Here we show that in sweet almond cultivars, highly methylated TE insertions surround a gene involved in the biosynthesis of amygdalin, whose reduced expression has been correlated with the sweet almond phenotype. Altogether, our results suggest a key role of TEs in the recent history and diversification of almond and its close relative peach.

The PTI to ETI Continuum in Phytophthora-Plant Interactions

Phytophthora species are notorious pathogens of several economically important crop plants. Several general elicitors, commonly referred to as Pathogen-Associated Molecular Patterns (PAMPs), from Phytophthora spp. have been identified that are recognized by the plant receptors to trigger induced defense responses in a process termed PAMP-triggered Immunity (PTI). Adapted Phytophthora pathogens have evolved multiple strategies to evade PTI. They can either modify or suppress their elicitors to avoid recognition by host and modulate host defense responses by deploying hundreds of effectors, which suppress host defense and physiological processes by modulating components involved in calcium and MAPK signaling, alternative splicing, RNA interference, vesicle trafficking, cell-to-cell trafficking, proteolysis and phytohormone signaling pathways. In incompatible interactions, resistant host plants perceive effector-induced modulations through resistance proteins and activate downstream components of defense responses in a quicker and more robust manner called effector-triggered-immunity (ETI). When pathogens overcome PTI-usually through effectors in the absence of R proteins-effectors-triggered susceptibility (ETS) ensues. Qualitatively, many of the downstream defense responses overlap between PTI and ETI. In general, these multiple phases of Phytophthora-plant interactions follow the PTI-ETS-ETI paradigm, initially proposed in the zigzag model of plant immunity. However, based on several examples, in Phytophthora-plant interactions, boundaries between these phases are not distinct but are rather blended pointing to a PTI-ETI continuum.

Transcriptome characterization and expression profiling in chestnut cultivars resistant or susceptible to the gall wasp Dryocosmus kuriphilus

The oriental gall wasp Dryocosmus kuriphilus represents a limiting pest for the European Chestnut (Castanea sativa, Fagaceae) as it creates severe yield losses. The European Chestnut is a deciduous tree, having major social, economic and environmental importance in Southern Europe, covering an area of 2.53 million hectares, including 75,000 ha devoted to fruit production. Cultivars show different susceptibility and very few are resistant to gall wasp. To deeply investigate the plant response and understand which factors can lead the plant to develop or not the gall, the study of transcriptome is basic (fundamental). To date, little transcriptomic information are available for C. sativa species. Hence, we present a de novo assembly of the chestnut transcriptome of the resistant Euro-Japanese hybrid 'Bouche de Bétizac' (BB) and the susceptible cultivar 'Madonna' (M), collecting RNA from buds at different stages of budburst. The two transcriptomes were assembled into 34,081 (BB) and 30,605 (M) unigenes, respectively. The former was used as a reference sequence for further characterization analyses, highlighting the presence of 1444 putative resistance gene analogs (RGAs) and about 1135 unigenes, as putative MiRNA targets. A global quantitative transcriptome profiling comparing the resistant and the susceptible cultivars, in the presence or not of the gall wasp, revealed some GO enrichments as "response to stimulus" (GO:0050896), and "developmental processes" (e.g., post-embryonic development, GO:0009791). Many up-regulated genes appeared to be transcription factors (e.g., RAV1, AP2/ERF, WRKY33) or protein regulators (e.g., RAPTOR1B) and storage proteins (e.g., LEA D29) involved in "post-embryonic development". Our analysis was able to provide a large amount of information, including 7k simple sequence repeat (SSR) and 335k single-nucleotide polymorphism (SNP)/INDEL markers, and generated the first reference unigene catalog for the European Chestnut. The transcriptome data for C. sativa will contribute to understand the genetic basis of the resistance to gall wasp and will provide useful information for next molecular genetic studies of this species and its relatives.

PRGdb 3.0: a comprehensive platform for prediction and analysis of plant disease resistance genes

The Plant Resistance Genes database (PRGdb; http://prgdb.org) has been redesigned with a new user interface, new sections, new tools and new data for genetic improvement, allowing easy access not only to the plant science research community but also to breeders who want to improve plant disease resistance. The home page offers an overview of easy-to-read search boxes that streamline data queries and directly show plant species for which data from candidate or cloned genes have been collected. Bulk data files and curated resistance gene annotations are made available for each plant species hosted. The new Gene Model view offers detailed information on each cloned resistance gene structure to highlight shared attributes with other genes. PRGdb 3.0 offers 153 reference resistance genes and 177 072 annotated candidate Pathogen Receptor Genes (PRGs). Compared to the previous release, the number of putative genes has been increased from 106 to 177 K from 76 sequenced Viridiplantae and algae genomes. The DRAGO 2 tool, which automatically annotates and predicts (PRGs) from DNA and amino acid with high accuracy and sensitivity, has been added. BLAST search has been implemented to offer users the opportunity to annotate and compare their own sequences. The improved section on plant diseases displays useful information linked to genes and genomes to connect complementary data and better address specific needs. Through, a revised and enlarged collection of data, the development of new tools and a renewed portal, PRGdb 3.0 engages the plant science community in developing a consensus plan to improve knowledge and strategies to fight diseases that afflict main crops and other plants.

The Gossypium hirsutum TIR-NBS-LRR gene GhDSC1 mediates resistance against Verticillium wilt

Improving genetic resistance is a preferred method to manage Verticillium wilt of cotton and other hosts. Identifying host resistance is difficult because of the dearth of resistance genes against this pathogen. Previously, a novel candidate gene involved in Verticillium wilt resistance was identified by a genome-wide association study using a panel of Gossypium hirsutum accessions. In this study, we cloned the candidate resistance gene from cotton that encodes a protein sharing homology with the TIR-NBS-LRR receptor-like defence protein DSC1 in Arabidopsis thaliana (hereafter named GhDSC1). GhDSC1 expressed at higher levels in response to Verticillium wilt and jasmonic acid (JA) treatment in resistant cotton cultivars as compared to susceptible cultivars and its product was localized to nucleus. The transfer of GhDSC1 to Arabidopsis conferred Verticillium resistance in an A. thaliana dsc1 mutant. This resistance response was associated with reactive oxygen species (ROS) accumulation and increased expression of JA-signalling-related genes. Furthermore, the expression of GhDSC1 in response to Verticillium wilt and JA signalling in A. thaliana displayed expression patterns similar to GhCAMTA3 in cotton under identical conditions, suggesting a coordinated DSC1 and CAMTA3 response in A. thaliana to Verticillium wilt. Analyses of GhDSC1 sequence polymorphism revealed a single nucleotide polymorphism (SNP) difference between resistant and susceptible cotton accessions, within the P-loop motif encoded by GhDSC1. This SNP difference causes ineffective activation of defence response in susceptible cultivars. These results demonstrated that GhDSC1 confers Verticillium resistance in the model plant system of A. thaliana, and therefore represents a suitable candidate for the genetic engineering of Verticillium wilt resistance in cotton.

The Gossypium hirsutum TIR-NBS-LRR gene GhDSC1 mediates resistance against Verticillium wilt

Improving genetic resistance is a preferred method to manage Verticillium wilt of cotton and other hosts. Identifying host resistance is difficult because of the dearth of resistance genes against this pathogen. Previously, a novel candidate gene involved in Verticillium wilt resistance was identified by a genome-wide association study using a panel of Gossypium hirsutum accessions. In this study, we cloned the candidate resistance gene from cotton that encodes a protein sharing homology with the TIR-NBS-LRR receptor-like defence protein DSC1 in Arabidopsis thaliana (hereafter named GhDSC1). GhDSC1 expressed at higher levels in response to Verticillium wilt and jasmonic acid (JA) treatment in resistant cotton cultivars as compared to susceptible cultivars and its product was localized to nucleus. The transfer of GhDSC1 to Arabidopsis conferred Verticillium resistance in an A. thaliana dsc1 mutant. This resistance response was associated with reactive oxygen species (ROS) accumulation and increased expression of JA-signalling-related genes. Furthermore, the expression of GhDSC1 in response to Verticillium wilt and JA signalling in A. thaliana displayed expression patterns similar to GhCAMTA3 in cotton under identical conditions, suggesting a coordinated DSC1 and CAMTA3 response in A. thaliana to Verticillium wilt. Analyses of GhDSC1 sequence polymorphism revealed a single nucleotide polymorphism (SNP) difference between resistant and susceptible cotton accessions, within the P-loop motif encoded by GhDSC1. This SNP difference causes ineffective activation of defence response in susceptible cultivars. These results demonstrated that GhDSC1 confers Verticillium resistance in the model plant system of A. thaliana, and therefore represents a suitable candidate for the genetic engineering of Verticillium wilt resistance in cotton.

The island cotton NBS-LRR gene GbaNA1 confers resistance to the non-race 1 Verticillium dahliae isolate Vd991

Wilt caused by Verticillium dahliae significantly reduces cotton yields, as host resistance in commercially cultivated Gossypium species is lacking. Understanding the molecular basis of disease resistance in non-commercial Gossypium species could galvanize the development of Verticillium wilt resistance in cultivated species. Nucleotide-binding site leucine-rich repeat (NBS-LRR) proteins play a central role in plant defence against pathogens. In this study, we focused on the relationship between a locus enriched with eight NBS-LRR genes and Verticillium wilt resistance in G. barbadense. Independent virus-induced gene silencing of each of the eight NBS-LRR genes in G. barbadense cultivar Hai 7124 revealed that silencing of GbaNA1 alone compromised the resistance of G. barbadense to V. dahliae isolate Vd991. In cultivar Hai 7124, GbaNA1 could be induced by V. dahliae isolate Vd991 and by ethylene, jasmonic acid and salicylic acid. Nuclear protein localization of GbaNA1 was demonstrated by transient expression. Sequencing of the GbaNA1 orthologue in nine G. hirsutum accessions revealed that all carried a non-functional allele, caused by a premature peptide truncation. In addition, all 10 G. barbadense and nine G. hirsutum accessions tested carried a full-length (∼1140 amino acids) homologue of the V. dahliae race 1 resistance gene Gbve1, although some sequence polymorphisms were observed. Verticillium dahliae Vd991 is a non-race 1 isolate that lacks the Ave1 gene. Thus, the resistance imparted by GbaNA1 appears to be mediated by a mechanism distinct from recognition of the fungal effector Ave1.

Comparative genomic and transcriptomic analysis of selected fatty acid biosynthesis genes and CNL disease resistance genes in oil palm

Comparative genomics and transcriptomic analyses were performed on two agronomically important groups of genes from oil palm versus other major crop species and the model organism, Arabidopsis thaliana. The first analysis was of two gene families with key roles in regulation of oil quality and in particular the accumulation of oleic acid, namely stearoyl ACP desaturases (SAD) and acyl-acyl carrier protein (ACP) thioesterases (FAT). In both cases, these were found to be large gene families with complex expression profiles across a wide range of tissue types and developmental stages. The detailed classification of the oil palm SAD and FAT genes has enabled the updating of the latest version of the oil palm gene model. The second analysis focused on disease resistance (R) genes in order to elucidate possible candidates for breeding of pathogen tolerance/resistance. Ortholog analysis showed that 141 out of the 210 putative oil palm R genes had homologs in banana and rice. These genes formed 37 clusters with 634 orthologous genes. Classification of the 141 oil palm R genes showed that the genes belong to the Kinase (7), CNL (95), MLO-like (8), RLK (3) and Others (28) categories. The CNL R genes formed eight clusters. Expression data for selected R genes also identified potential candidates for breeding of disease resistance traits. Furthermore, these findings can provide information about the species evolution as well as the identification of agronomically important genes in oil palm and other major crops.

Identifying Quantitative Trait Loci (QTLs) and Developing Diagnostic Markers Linked to Orange Rust Resistance in Sugarcane (Saccharum spp.)

Sugarcane (Saccharum spp.) is an important economic crop, contributing up to 80% of table sugar used in the world and has become a promising feedstock for biofuel production. Sugarcane production has been threatened by many diseases, and fungicide applications for disease control have been opted out for sustainable agriculture. Orange rust is one of the major diseases impacting sugarcane production worldwide. Identifying quantitative trait loci (QTLs) and developing diagnostic markers are valuable for breeding programs to expedite release of superior sugarcane cultivars for disease control. In this study, an F1 segregating population derived from a cross between two hybrid sugarcane clones, CP95-1039 and CP88-1762, was evaluated for orange rust resistance in replicated trails. Three QTLs controlling orange rust resistance in sugarcane (qORR109, qORR4 and qORR102) were identified for the first time ever, which can explain 58, 12 and 8% of the phenotypic variation, separately. We also characterized 1,574 sugarcane putative resistance (R) genes. These sugarcane putative R genes and simple sequence repeats in the QTL intervals were further used to develop diagnostic markers for marker-assisted selection of orange rust resistance. A PCR-based Resistance gene-derived maker, G1 was developed, which showed significant association with orange rust resistance. The putative QTLs and marker developed in this study can be effectively utilized in sugarcane breeding programs to facilitate the selection process, thus contributing to the sustainable agriculture for orange rust disease control.

Heterologous Expression of the Cotton NBS-LRR Gene GbaNA1 Enhances Verticillium Wilt Resistance in Arabidopsis

Verticillium wilt caused by Verticillium dahliae results in severe losses in cotton, and is economically the most destructive disease of this crop. Improving genetic resistance is the cleanest and least expensive option to manage Verticillium wilt. Previously, we identified the island cotton NBS-LRR-encoding gene GbaNA1 that confers resistance to the highly virulent V. dahliae isolate Vd991. In this study, we expressed cotton GbaNA1 in the heterologous system of Arabidopsis thaliana and investigated the defense response mediated by GbaNA1 following inoculations with V. dahliae. Heterologous expression of GbaNA1 conferred Verticillium wilt resistance in A. thaliana. Moreover, overexpression of GbaNA1 enabled recovery of the resistance phenotype of A. thaliana mutants that had lost the function of GbaNA1 ortholog gene. Investigations of the defense response in A. thaliana showed that the reactive oxygen species (ROS) production and the expression of genes associated with the ethylene signaling pathway were enhanced significantly following overexpression of GbaNA1. Intriguingly, overexpression of the GbaNA1 ortholog from Gossypium hirsutum (GhNA1) in A. thaliana did not induce the defense response of ROS production due to the premature termination of GhNA1, which lacks the encoded NB-ARC and LRR motifs. GbaNA1 therefore confers Verticillium wilt resistance in A. thaliana by the activation of ROS production and ethylene signaling. These results demonstrate the functional conservation of the NBS-LRR-encoding GbaNA1 in a heterologous system, and the mechanism of this resistance, both of which may prove valuable in incorporating GbaNA1-mediated resistance into other plant species.

Transcriptomic Analysis Reveal the Molecular Mechanisms of Wheat Higher-Temperature Seedling-Plant Resistance to Puccinia striiformis f. sp. tritici

Stripe rust, caused by Puccinia striiformis f. sp. tritici (Pst), is a destructive disease of wheat worldwide. The disease is preferably controlled by growing resistant cultivars. Wheat cultivar Xiaoyan 6 (XY 6) has been resistant to stripe rust since its release. In the previous studies, XY 6 was found to have higher-temperature seedling-plant (HTSP) resistance. However, the molecular mechanisms of HTSP resistance were not clear. To identify differentially expressed genes (DEGs) involved in HTSP resistance, we sequenced 30 cDNA libraries constructed from XY 6 seedlings exposed to several temperature treatments. Compared to the constant normal (15°C) and higher (20°C) temperature treatments, 1395 DEGs were identified in seedlings exposed to 20°C for 24 h (to activate HTSP resistance) and then kept at 15°C. These DEGs were located on all 21 chromosomes, with 29.2% on A, 41.1% on B and 29.7% on D genomes, by mapping to the Chinese Spring wheat genome. The 1395 DEGs were enriched in ribosome, plant-pathogen interaction and glycerolipid metabolism pathways, and some of them were identified as hub proteins (phosphatase 2C10), resistance protein homologs, WRKY transcription factors and protein kinases. The majority of these genes were up-regulated in HTSP resistance. Based on the differential expression, we found that phosphatase 2C10 and LRR receptor-like serine/threonine protein kinases are particularly interesting as they may be important for HTSP resistance through interacting with different resistance proteins, leading to a hypersensitive response.

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.

Evidence-based gene models for structural and functional annotations of the oil palm genome

Background

Oil palm is an important source of edible oil. The importance of the crop, as well as its long breeding cycle (10-12 years) has led to the sequencing of its genome in 2013 to pave the way for genomics-guided breeding. Nevertheless, the first set of gene predictions, although useful, had many fragmented genes. Classification and characterization of genes associated with traits of interest, such as those for fatty acid biosynthesis and disease resistance, were also limited. Lipid-, especially fatty acid (FA)-related genes are of particular interest for the oil palm as they specify oil yields and quality. This paper presents the characterization of the oil palm genome using different gene prediction methods and comparative genomics analysis, identification of FA biosynthesis and disease resistance genes, and the development of an annotation database and bioinformatics tools.

Results

Using two independent gene-prediction pipelines, Fgenesh++ and Seqping, 26,059 oil palm genes with transcriptome and RefSeq support were identified from the oil palm genome. These coding regions of the genome have a characteristic broad distribution of GC₃ (fraction of cytosine and guanine in the third position of a codon) with over half the GC₃-rich genes (GC₃ ≥ 0.75286) being intronless. In comparison, only one-seventh of the oil palm genes identified are intronless. Using comparative genomics analysis, characterization of conserved domains and active sites, and expression analysis, 42 key genes involved in FA biosynthesis in oil palm were identified. For three of them, namely EgFABF, EgFABH and EgFAD3, segmental duplication events were detected. Our analysis also identified 210 candidate resistance genes in six classes, grouped by their protein domain structures.

Conclusions

We present an accurate and comprehensive annotation of the oil palm genome, focusing on analysis of important categories of genes (GC₃-rich and intronless), as well as those associated with important functions, such as FA biosynthesis and disease resistance. The study demonstrated the advantages of having an integrated approach to gene prediction and developed a computational framework for combining multiple genome annotations. These results, available in the oil palm annotation database (http://palmxplore.mpob.gov.my), will provide important resources for studies on the genomes of oil palm and related crops.

Reviewers

This article was reviewed by Alexander Kel, Igor Rogozin, and Vladimir A. Kuznetsov.

Electronic supplementary material

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

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.

Prediction and analysis of three gene families related to leaf rust (Puccinia triticina) resistance in wheat (Triticum aestivum L.)

BACKGROUND

The resistance to leaf rust (Lr) caused by Puccinia triticina in wheat (Triticum aestivum L.) has been well studied over the past decades with over 70 Lr genes being mapped on different chromosomes and numerous QTLs (quantitative trait loci) being detected or mapped using DNA markers. Such resistance is often divided into race-specific and race-nonspecific resistance. The race-nonspecific resistance can be further divided into resistance to most or all races of the same pathogen and resistance to multiple pathogens. At the molecular level, these three types of resistance may cover across the whole spectrum of pathogen specificities that are controlled by genes encoding different protein families in wheat. The objective of this study is to predict and analyze genes in three such families: NBS-LRR (nucleotide-binding sites and leucine-rich repeats or NLR), START (Steroidogenic Acute Regulatory protein [STaR] related lipid-transfer) and ABC (ATP-Binding Cassette) transporter. The focus of the analysis is on the patterns of relationships between these protein-coding genes within the gene families and QTLs detected for leaf rust resistance.

RESULTS

We predicted 526 ABC, 1117 NLR and 144 START genes in the hexaploid wheat genome through a domain analysis of wheat proteome. Of the 1809 SNPs from leaf rust resistance QTLs in seedling and adult stages of wheat, 126 SNPs were found within coding regions of these genes or their neighborhood (5 Kb upstream from transcription start site [TSS] or downstream from transcription termination site [TTS] of the genes). Forty-three of these SNPs for adult resistance and 18 SNPs for seedling resistance reside within coding or neighboring regions of the ABC genes whereas 14 SNPs for adult resistance and 29 SNPs for seedling resistance reside within coding or neighboring regions of the NLR gene. Moreover, we found 17 nonsynonymous SNPs for adult resistance and five SNPs for seedling resistance in the ABC genes, and five nonsynonymous SNPs for adult resistance and six SNPs for seedling resistance in the NLR genes. Most of these coding SNPs were predicted to alter encoded amino acids and such information may serve as a starting point towards more thorough molecular and functional characterization of the designated Lr genes. Using the primer sequences of 99 known non-SNP markers from leaf rust resistance QTLs, we found candidate genes closely linked to these markers, including Lr34 with distances to its two gene-specific markers being 1212 bases (to cssfr1) and 2189 bases (to cssfr2).

CONCLUSION

This study represents a comprehensive analysis of ABC, NLR and START genes in the hexaploid wheat genome and their physical relationships with QTLs for leaf rust resistance at seedling and adult stages. Our analysis suggests that the ABC (and START) genes are more likely to be co-located with QTLs for race-nonspecific, adult resistance whereas the NLR genes are more likely to be co-located with QTLs for race-specific resistance that would be often expressed at the seedling stage. Though our analysis was hampered by inaccurate or unknown physical positions of numerous QTLs due to the incomplete assembly of the complex hexaploid wheat genome that is currently available, the observed associations between (i) QTLs for race-specific resistance and NLR genes and (ii) QTLs for nonspecific resistance and ABC genes will help discover SNP variants for leaf rust resistance at seedling and adult stages. The genes containing nonsynonymous SNPs are promising candidates that can be investigated in future studies as potential new sources of leaf rust resistance in wheat breeding.

Genome-wide comparative analysis of NBS-encoding genes in four Gossypium species

Background

Nucleotide binding site (NBS) genes encode a large family of disease resistance (R) proteins in plants. The availability of genomic data of the two diploid cotton species, Gossypium arboreum and Gossypium raimondii, and the two allotetraploid cotton species, Gossypium hirsutum (TM-1) and Gossypium barbadense allow for a more comprehensive and systematic comparative study of NBS-encoding genes to elucidate the mechanisms of cotton disease resistance.

Results

Based on the genome assembly data, 246, 365, 588 and 682 NBS-encoding genes were identified in G. arboreum, G. raimondii, G. hirsutum and G. barbadense, respectively. The distribution of NBS-encoding genes among the chromosomes was nonrandom and uneven, and was tended to form clusters. Gene structure analysis showed that G. arboreum and G. hirsutum possessed a greater proportion of CN, CNL, and N genes and a lower proportion of NL, TN and TNL genes compared to that of G. raimondii and G. barbadense, while the percentages of RN and RNL genes remained relatively unchanged. The percentage changes among them were largest for TNL genes, about 7 times. Exon statistics showed that the average exon numbers per NBS gene in G. raimondii and G. barbadense were all greater than that in G. arboretum and G. hirsutum. Phylogenetic analysis revealed that the TIR-NBS genes of G. barbadense were closely related with that of G. raimondii. Sequence similarity analysis showed that diploid cotton G. arboreum possessed a larger proportion of NBS-encoding genes similar to that of allotetraploid cotton G. hirsutum, while diploid G. raimondii possessed a larger proportion of NBS-encoding genes similar to that of allotetraploid cotton G. barbadense. The synteny analysis showed that more NBS genes in G. raimondii and G. arboreum were syntenic with that in G. barbadense and G. hirsutum, respectively.

Conclusions

The structural architectures, amino acid sequence similarities and synteny of NBS-encoding genes between G. arboreum and G. hirsutum, and between G. raimondii and G. barbadense were the highest among comparisons between the diploid and allotetraploid genomes, indicating that G. hirsutum inherited more NBS-encoding genes from G. arboreum, while G. barbadense inherited more NBS-encoding genes from G. raimondii. This asymmetric evolution of NBS-encoding genes may help to explain why G. raimondii and G. barbadense are more resistant to Verticillium wilt, whereas G. arboreum and G. hirsutum are more susceptible to Verticillium wilt. The disease resistances of the allotetraploid cotton were related to their NBS-encoding genes especially in regard from which diploid progenitor they were derived, and the TNL genes may have a significant role in disease resistance to Verticillium wilt in G. raimondii and G. barbadense.

Electronic supplementary material

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

Quantitative Trait Loci and Resistance Gene Analogs Associated with Maize White Spot Resistance

Maize white spot (MWS), caused by the bacterium Pantoea ananatis, is one of the most important maize foliar diseases in tropical and subtropical regions, causing significant yield losses. Despite its economic importance, genetic studies of MWS are scarce. The aim of this study was to map quantitative trait loci (QTL) associated with MWS resistance and to identify resistance gene analogs (RGA) underlying these QTL. QTL mapping was performed in a tropical maize F_2:3 population, which was genotyped with simple-sequence repeat and RGA-tagged markers and phenotyped for the response to MWS in two Brazilian southeastern locations. Nine QTL explained approximately 70% of the phenotypic variance for MWS resistance at each location, with two of them consistently detected in both environments. Data mining using 112 resistance genes cloned from different plant species revealed 1,697 RGA distributed in clusters within the maize genome. The RGA Pto19, Pto20, Pto99, and Xa26.151.4 were genetically mapped within MWS resistance QTL on chromosomes 4 and 8 and were preferentially expressed in the resistant parental line at locations where their respective QTL occurred. The consistency of QTL mapping, in silico prediction, and gene expression analyses revealed RGA and genomic regions suitable for marker-assisted selection to improve MWS resistance.

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.

The Diversification of Plant NBS-LRR Defense Genes Directs the Evolution of MicroRNAs That Target Them

High expression of plant nucleotide binding site leucine-rich repeat (NBS-LRR) defense genes is often lethal to plant cells, a phenotype perhaps associated with fitness costs. Plants implement several mechanisms to control the transcript level of NBS-LRR defense genes. As negative transcriptional regulators, diverse miRNAs target NBS-LRRs in eudicots and gymnosperms. To understand the evolutionary benefits of this miRNA-NBS-LRR regulatory system, we investigated the NBS-LRRs of 70 land plants, coupling this analysis with extensive small RNA data. A tight association between the diversity of NBS-LRRs and miRNAs was found. The miRNAs typically target highly duplicated NBS-LRRs In comparison, families of heterogeneous NBS-LRRs were rarely targeted by miRNAs in Poaceae and Brassicaceae genomes. We observed that duplicated NBS-LRRs from different gene families periodically gave birth to new miRNAs. Most of these newly emerged miRNAs target the same conserved, encoded protein motif of NBS-LRRs, consistent with a model of convergent evolution for these miRNAs. By assessing the interactions between miRNAs and NBS-LRRs, we found nucleotide diversity in the wobble position of the codons in the target site drives the diversification of miRNAs. Taken together, we propose a co-evolutionary model of plant NBS-LRRs and miRNAs hypothesizing how plants balance the benefits and costs of NBS-LRR defense genes.

Using Genotyping by Sequencing to Map Two Novel Anthracnose Resistance Loci in Sorghum bicolor

Colletotrichum sublineola is an aggressive fungal pathogen that causes anthracnose in sorghum [Sorghum bicolor (L.) Moench]. The obvious symptoms of anthracnose are leaf blight and stem rot. Sorghum, the fifth most widely grown cereal crop in the world, can be highly susceptible to the disease, most notably in hot and humid environments. In the southeastern United States the acreage of sorghum has been increasing steadily in recent years, spurred by growing interest in producing biofuels, bio-based products, and animal feed. Resistance to anthracnose is, therefore, of paramount importance for successful sorghum production in this region. To identify anthracnose resistance loci present in the highly resistant cultivar ‘Bk7’, a biparental mapping population of F_3:4 and F_4:5 sorghum lines was generated by crossing ‘Bk7’ with the susceptible inbred ‘Early Hegari-Sart’. Lines were phenotyped in three environments and in two different years following natural infection. The population was genotyped by sequencing. Following a stringent custom filtering protocol, totals of 5186 and 2759 informative SNP markers were identified in the two populations. Segregation data and association analysis identified resistance loci on chromosomes 7 and 9, with the resistance alleles derived from ‘Bk7’. Both loci contain multiple classes of defense-related genes based on sequence similarity and gene ontologies. Genetic analysis following an independent selection experiment of lines derived from a cross between ‘Bk7’ and sweet sorghum ‘Mer81-4’ narrowed the resistance locus on chromosome 9 substantially, validating this QTL. As observed in other species, sorghum appears to have regions of clustered resistance genes. Further characterization of these regions will facilitate the development of novel germplasm with resistance to anthracnose and other diseases.

Whole Genome Re-Sequencing and Characterization of Powdery Mildew Disease-Associated Allelic Variation in Melon

Powdery mildew is one of the most common fungal diseases in the world. This disease frequently affects melon (Cucumis melo L.) and other Cucurbitaceous family crops in both open field and greenhouse cultivation. One of the goals of genomics is to identify the polymorphic loci responsible for variation in phenotypic traits. In this study, powdery mildew disease assessment scores were calculated for four melon accessions, 'SCNU1154', 'Edisto47', 'MR-1', and 'PMR5'. To investigate the genetic variation of these accessions, whole genome re-sequencing using the Illumina HiSeq 2000 platform was performed. A total of 754,759,704 quality-filtered reads were generated, with an average of 82.64% coverage relative to the reference genome. Comparisons of the sequences for the melon accessions revealed around 7.4 million single nucleotide polymorphisms (SNPs), 1.9 million InDels, and 182,398 putative structural variations (SVs). Functional enrichment analysis of detected variations classified them into biological process, cellular component and molecular function categories. Further, a disease-associated QTL map was constructed for 390 SNPs and 45 InDels identified as related to defense-response genes. Among them 112 SNPs and 12 InDels were observed in powdery mildew responsive chromosomes. Accordingly, this whole genome re-sequencing study identified SNPs and InDels associated with defense genes that will serve as candidate polymorphisms in the search for sources of resistance against powdery mildew disease and could accelerate marker-assisted breeding in melon.

A high-quality carrot genome assembly provides new insights into carotenoid accumulation and asterid genome evolution

We report a high-quality chromosome-scale assembly and analysis of the carrot (Daucus carota) genome, the first sequenced genome to include a comparative evolutionary analysis among members of the euasterid II clade. We characterized two new polyploidization events, both occurring after the divergence of carrot from members of the Asterales order, clarifying the evolutionary scenario before and after radiation of the two main asterid clades. Large- and small-scale lineage-specific duplications have contributed to the expansion of gene families, including those with roles in flowering time, defense response, flavor, and pigment accumulation. We identified a candidate gene, DCAR_032551, that conditions carotenoid accumulation (Y) in carrot taproot and is coexpressed with several isoprenoid biosynthetic genes. The primary mechanism regulating carotenoid accumulation in carrot taproot is not at the biosynthetic level. We hypothesize that DCAR_032551 regulates upstream photosystem development and functional processes, including photomorphogenesis and root de-etiolation.

Integration of decoy domains derived from protein targets of pathogen effectors into plant immune receptors is widespread

Plant immune receptors of the class of nucleotide-binding and leucine-rich repeat domain (NLR) proteins can contain additional domains besides canonical NB-ARC (nucleotide-binding adaptor shared by APAF-1, R proteins, and CED-4 (NB-ARC)) and leucine-rich repeat (LRR) domains. Recent research suggests that these additional domains act as integrated decoys recognizing effectors from pathogens. Proteins homologous to integrated decoys are suspected to be effector targets and involved in disease or resistance. Here, we scrutinized 31 entire plant genomes to identify putative integrated decoy domains in NLR proteins using the Interpro search. The involvement of the Zinc Finger-BED type (ZBED) protein containing a putative decoy domain, called BED, in rice (Oryza sativa) resistance was investigated by evaluating susceptibility to the blast fungus Magnaporthe oryzae in rice over-expression and knock-out mutants. This analysis showed that all plants tested had integrated various atypical protein domains into their NLR proteins (on average 3.5% of all NLR proteins). We also demonstrated that modifying the expression of the ZBED gene modified disease susceptibility. This study suggests that integration of decoy domains in NLR immune receptors is widespread and frequent in plants. The integrated decoy model is therefore a powerful concept to identify new proteins involved in disease resistance. Further in-depth examination of additional domains in NLR proteins promises to unravel many new proteins of the plant immune system.

Genome and transcriptome analysis of the Mesoamerican common bean and the role of gene duplications in establishing tissue and temporal specialization of genes

BACKGROUND

Legumes are the third largest family of angiosperms and the second most important crop class. Legume genomes have been shaped by extensive large-scale gene duplications, including an approximately 58 million year old whole genome duplication shared by most crop legumes.

RESULTS

We report the genome and the transcription atlas of coding and non-coding genes of a Mesoamerican genotype of common bean (Phaseolus vulgaris L., BAT93). Using a comprehensive phylogenomics analysis, we assessed the past and recent evolution of common bean, and traced the diversification of patterns of gene expression following duplication. We find that successive rounds of gene duplications in legumes have shaped tissue and developmental expression, leading to increased levels of specialization in larger gene families. We also find that many long non-coding RNAs are preferentially expressed in germ-line-related tissues (pods and seeds), suggesting that they play a significant role in fruit development. Our results also suggest that most bean-specific gene family expansions, including resistance gene clusters, predate the split of the Mesoamerican and Andean gene pools.

CONCLUSIONS

The genome and transcriptome data herein generated for a Mesoamerican genotype represent a counterpart to the genomic resources already available for the Andean gene pool. Altogether, this information will allow the genetic dissection of the characters involved in the domestication and adaptation of the crop, and their further implementation in breeding strategies for this important crop.

Genome-Wide Analysis of NBS-LRR Genes in Sorghum Genome Revealed Several Events Contributing to NBS-LRR Gene Evolution in Grass Species

The nucleotide-binding site (NBS)–leucine-rich repeat (LRR) gene family is crucially important for offering resistance to pathogens. To explore evolutionary conservation and variability of NBS-LRR genes across grass species, we identified 88, 107, 24, and 44 full-length NBS-LRR genes in sorghum, rice, maize, and Brachypodium, respectively. A comprehensive analysis was performed on classification, genome organization, evolution, expression, and regulation of these NBS-LRR genes using sorghum as a representative of grass species. In general, the full-length NBS-LRR genes are highly clustered and duplicated in sorghum genome mainly due to local duplications. NBS-LRR genes have basal expression levels and are highly potentially targeted by miRNA. The number of NBS-LRR genes in the four grass species is positively correlated with the gene clustering rate. The results provided a valuable genomic resource and insights for functional and evolutionary studies of NBS-LRR genes in grass species.

NBSPred: a support vector machine-based high-throughput pipeline for plant resistance protein NBSLRR prediction

The nucleotide binding site leucine-rich repeats (NBSLRRs) belong to one of the largest known families of disease resistance genes that encode resistance proteins (R-protein) against the pathogens of plants. Various defence mechanisms have explained the regulation of plant immunity, but still, we have limited understanding about plant defence against different pathogens. Identification of R-proteins and proteins having R-protein-like features across the genome, transcriptome and proteome would be highly useful to develop the global understanding of plant defence mechanisms, but it is laborious and time-consuming task. Therefore, we have developed a support vector machine-based high-throughput pipeline called NBSPred to differentiate NBSLRR and NBSLRR-like protein from Non-NBSLRR proteins from genome, transcriptome and protein sequences. The pipeline was tested and validated with input sequences from three dicot and two monocot plants including Arabidopsis thaliana, Boechera stricta, Brachypodium distachyon Solanum lycopersicum and Zea mays.

AVAILABILITY AND IMPLEMENTATION

The NBSPred pipeline is available at http://soilecology.biol.lu.se/nbs/

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

CONTACT

sandeep.kushwaha@biol.lu.se.

De novo assembly and comparative analysis of root transcriptomes from different varieties of Panax ginseng C. A. Meyer grown in different environments

Panax ginseng C. A. Meyer is an important traditional herb in eastern Asia. It contains ginsenosides, which are primary bioactive compounds with medicinal properties. Although ginseng has been cultivated since at least the Ming dynasty to increase production, cultivated ginseng has lower quantities of ginsenosides and lower disease resistance than ginseng grown under natural conditions. We extracted root RNA from six varieties of fifth-year P. ginseng cultivars representing four different growth conditions, and performed Illumina paired-end sequencing. In total, 163,165,706 raw reads were obtained and used to generate a de novo transcriptome that consisted of 151,763 contigs (76,336 unigenes), of which 100,648 contigs (66.3%) were successfully annotated. Differential expression analysis revealed that most differentially expressed genes (DEGs) were upregulated (246 out of 258, 95.3%) in ginseng grown under natural conditions compared with that grown under artificial conditions. These DEGs were enriched in gene ontology (GO) terms including response to stimuli and localization. In particular, some key ginsenoside biosynthesis-related genes, including HMG-CoA synthase (HMGS), mevalonate kinase (MVK), and squalene epoxidase (SE), were upregulated in wild-grown ginseng. Moreover, a high proportion of disease resistance-related genes were upregulated in wild-grown ginseng. This study is the first transcriptome analysis to compare wild-grown and cultivated ginseng, and identifies genes that may produce higher ginsenoside content and better disease resistance in the wild; these genes may have the potential to improve cultivated ginseng grown in artificial environments.

Whole genome re-sequencing of date palms yields insights into diversification of a fruit tree crop

Date palms (Phoenix dactylifera) are the most significant perennial crop in arid regions of the Middle East and North Africa. Here, we present a comprehensive catalogue of approximately seven million single nucleotide polymorphisms in date palms based on whole genome re-sequencing of a collection of 62 cultivars. Population structure analysis indicates a major genetic divide between North Africa and the Middle East/South Asian date palms, with evidence of admixture in cultivars from Egypt and Sudan. Genome-wide scans for selection suggest at least 56 genomic regions associated with selective sweeps that may underlie geographic adaptation. We report candidate mutations for trait variation, including nonsense polymorphisms and presence/absence variation in gene content in pathways for key agronomic traits. We also identify a copia-like retrotransposon insertion polymorphism in the R2R3 myb-like orthologue of the oil palm virescens gene associated with fruit colour variation. This analysis documents patterns of post-domestication diversification and provides a genomic resource for this economically important perennial tree crop.

A novel workflow correlating RNA-seq data to Phythophthora infestans resistance levels in wild Solanum species and potato clones

Comparative transcriptomics between species can provide valuable understanding of plant-pathogen interactions. Here, we focus on wild Solanum species and potato clones with varying degree of resistance against Phytophthora infestans, which causes the devastating late blight disease in potato. The transcriptomes of three wild Solanum species native to Southern Sweden, Solanum dulcamara, Solanum nigrum, and Solanum physalifolium were compared to three potato clones, Desiree (cv.), SW93-1015 and Sarpo Mira. Desiree and S. physalifolium are susceptible to P. infestans whereas the other four have different degrees of resistance. By building transcript families based on de novo assembled RNA-seq across species and clones and correlating these to resistance phenotypes, we created a novel workflow to identify families with expanded or depleted number of transcripts in relation to the P. infestans resistance level. Analysis was facilitated by inferring functional annotations based on the family structure and semantic clustering. More transcript families were expanded in the resistant clones and species and the enriched functions of these were associated to expected gene ontology (GO) terms for resistance mechanisms such as hypersensitive response, host programmed cell death and endopeptidase activity. However, a number of unexpected functions and transcripts were also identified, for example transmembrane transport and protein acylation expanded in the susceptible group and a cluster of Zinc knuckle family proteins expanded in the resistant group. Over 400 expressed putative resistance (R-)genes were identified and resistant clones Sarpo Mira and SW93-1015 had ca 25% more expressed putative R-genes than susceptible cultivar Desiree. However, no differences in numbers of susceptibility (S-)gene homologs were seen between species and clones. In addition, we identified P. infestans transcripts including effectors in the early stages of P. infestans-Solanum interactions.

Genetic analysis of Verticillium wilt resistance in a backcross inbred line population and a meta-analysis of quantitative trait loci for disease resistance in cotton

Background

Verticillium wilt (VW) and Fusarium wilt (FW), caused by the soil-borne fungi Verticillium dahliae and Fusarium oxysporum f. sp. vasinfectum, respectively, are two most destructive diseases in cotton production worldwide. Root-knot nematodes (Meloidogyne incognita, RKN) and reniform nematodes (Rotylenchulus reniformis, RN) cause the highest yield loss in the U.S. Planting disease resistant cultivars is the most cost effective control method. Numerous studies have reported mapping of quantitative trait loci (QTLs) for disease resistance in cotton; however, very few reliable QTLs were identified for use in genomic research and breeding.

Results

This study first performed a 4-year replicated test of a backcross inbred line (BIL) population for VW resistance, and 10 resistance QTLs were mapped based on a 2895 cM linkage map with 392 SSR markers. The 10 VW QTLs were then placed to a consensus linkage map with other 182 VW QTLs, 75 RKN QTLs, 27 FW QTLs, and 7 RN QTLs reported from 32 publications. A meta-analysis of QTLs identified 28 QTL clusters including 13, 8 and 3 QTL hotspots for resistance to VW, RKN and FW, respectively. The number of QTLs and QTL clusters on chromosomes especially in the A-subgenome was significantly correlated with the number of nucleotide-binding site (NBS) genes, and the distribution of QTLs between homeologous A- and D- subgenome chromosomes was also significantly correlated.

Conclusions

Ten VW resistance QTL identified in a 4-year replicated study have added useful information to the understanding of the genetic basis of VW resistance in cotton. Twenty-eight disease resistance QTL clusters and 24 hotspots identified from a total of 306 QTLs and linked SSR markers provide important information for marker-assisted selection and high resolution mapping of resistance QTLs and genes. The non-overlapping of most resistance QTL hotspots for different diseases indicates that their resistances are controlled by different genes.

Electronic supplementary material

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

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.

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.

Functional genomics of tomato for the study of plant immunity

Tomato (Solanum lycopersicum), along with many other economically valuable species, belongs to the Solanaceae family. Understanding how plants in this family defend themselves against pathogens offers the opportunity of improving yield and quality of their edible products. The use of functional genomics has contributed to this purpose through both traditional and recently developed techniques that allow determination of changes in transcript abundance during pathogen attack. Such changes can implicate the affected gene as participating in plant defense. Testing the involvement of these candidate genes in defense has relied largely on posttranscriptional gene silencing, particularly virus-induced gene silencing. We discuss how functional genomics has played a key role in our current understanding of the defense response in tomato and related species and what are the challenges and opportunities for the future.

Loss and retention of resistance genes in five species of the Brassicaceae family

BACKGROUND

Plants have evolved disease resistance (R) genes encoding for nucleotide-binding site (NB) and leucine-rich repeat (LRR) proteins with N-terminals represented by either Toll/Interleukin-1 receptor (TIR) or coiled-coil (CC) domains. Here, a genome-wide study of presence and diversification of CC-NB-LRR and TIR-NB-LRR encoding genes, and shorter domain combinations in 19 Arabidopsis thaliana accessions and Arabidopsis lyrata, Capsella rubella, Brassica rapa and Eutrema salsugineum are presented.

RESULTS

Out of 528 R genes analyzed, 12 CC-NB-LRR and 17 TIR-NB-LRR genes were conserved among the 19 A. thaliana genotypes, while only two CC-NB-LRRs, including ZAR1, and three TIR-NB-LRRs were conserved when comparing the five species. The RESISTANCE TO LEPTOSPHAERIA MACULANS 1 (RLM1) locus confers resistance to the Brassica pathogen L. maculans the causal agent of blackleg disease and has undergone conservation and diversification events particularly in B. rapa. On the contrary, the RLM3 locus important in the immune response towards Botrytis cinerea and Alternaria spp. has recently evolved in the Arabidopsis genus.

CONCLUSION

Our genome-wide analysis of the R gene repertoire revealed a large sequence variation in the 23 cruciferous genomes. The data provides further insights into evolutionary processes impacting this important gene family.