PRGdb: a bioinformatics platform for plant resistance gene analysis

GhCNGC31 is critical for conferring resistance to Verticillium wilt in cotton

In the past decades, cyclic nucleotide-gated ion channels (CNGCs) have been extensively studied in diploid species Arabidopsis thaliana. However, the functional diversification of CNGCs in crop plants, mostly polyploid, remains poorly understood. In allotetraploid Upland cotton (Gossypium hirsutum), GhCNGC31 is one of the multiple orthologs of AtCNGC2, being present in the plasma membrane, capable of interacting with itself and binding to calmodulins and cyclic nucleotides. GhCNGC31 knockdown plants exhibited slight growth inhibition, and became more susceptible to Verticillium dahliae infection, which was associated with the reduced lignin and flavonoid accumulation, impaired ROS (reactive oxygen species) burst, and down-regulation of defense-related genes PR1, JAZ2, LOX2, and RBOH10. RNA-Seq analysis identified 1817 differentially expressed genes from GhCNGC31 knockdown, of which 1184 (65%) were responsive to V. dahliae infection and accounted for 57% among a total of 2065 V. dahliae-responsive genes identified in this study. These GhCNGC31-regulated genes mainly function with cell wall organization and biogenesis, cellular carbohydrate metabolic or biosynthetic process, cellular component macromolecule biosynthetic process, and rhythmic process. They are significantly enriched in the pathways of plant MAPK signaling, plant-pathogen interaction, phenylpropanoid biosynthesis, and plant hormone signal transduction. A set of transcription factors (TFs) and resistance (R) genes are among the GhCNGC31-regulated genes, which are significantly over-represented with the TCP and WRKY TFs families, as well as with the R genes of T (TIR) and TNL (TIR-NB-LRR) classes. Together, our results unraveled a critical role of GhCNGC31 for conferring resistance to Verticillium wilt in cotton.

Comparative genomic and transcriptomic analyses provide new insight into symbiotic host specificity

Host specificity plays important roles in expanding the host range of rhizobia, while the genetic information responsible for host specificity remains largely unexplored. In this report, the roots of four symbiotic systems with notable different symbiotic phenotypes and the control were studied at four different post-inoculation time points by RNA sequencning (RNA-seq). The differentially expressed genes (DEGs) were divided into "found only in soybean or Lotus," "only expressed in soybean or Lotus," and "expressed in both hosts" according to the comparative genomic analysis. The distributions of enriched function ontology (GO) terms and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways vary significantly in different symbiotic systems. Host specific genes account for the majority of the DEGs involved in response to stimulus, associated with plant-pathogen interaction pathways, and encoding resistance (R) proteins, the symbiotic nitrogen fixation (SNF) proteins and the target proteins in the SNF-related modules. Our findings provided molecular candidates for better understanding the mechanisms of symbiotic host-specificity.

Mining of disease-resistance genes in Crocus sativus based on transcriptome sequencing

Introduction: Crocus sativus L. has an important medicinal and economic value in traditional perennial Chinese medicine. However, due to its unique growth characteristics, during cultivation it is highly susceptible to disease. The absence of effective resistance genes restricts us to breed new resistant varieties of C. sativus. Methods: In present study, comprehensive transcriptome sequencing was introduced to explore the disease resistance of the candidate gene in healthy and corm rot-infected C. sativus. Results and discussion: Totally, 43.72 Gb of clean data was obtained from the assembly to generate 65,337 unigenes. By comparing the gene expression levels, 7,575 differentially expressed genes (DEGs) were primarily screened. A majority of the DEGs were completely in charge of defense and metabolism, and 152 of them were annotated as pathogen recognition genes (PRGs) based on the PGRdb dataset. The expression of some transcription factors including NAC, MYB, and WRKY members, changed significantly based on the dataset of transcriptome sequencing. Therefore, this study provides us some valuable information for exploring candidate genes involved in the disease resistance in C. sativus.

The spatiotemporal profile of Dendrobium huoshanense and functional identification of bHLH genes under exogenous MeJA using comparative transcriptomics and genomics

Introduction

Alkaloids are one of the main medicinal components of Dendrobium species. Dendrobium alkaloids are mainly composed of terpene alkaloids. Jasmonic acid (JA) induce the biosynthesis of such alkaloids, mainly by enhancing the expression of JA-responsive genes to increase plant resistance and increase the content of alkaloids. Many JA-responsive genes are the target genes of bHLH transcription factors (TFs), especially the MYC2 transcription factor.

Methods

In this study, the differentially expressed genes involved in the JA signaling pathway were screened out from Dendrobium huoshanense using comparative transcriptomics approaches, revealing the critical roles of basic helix-loop-helix (bHLH) family, particularly the MYC2 subfamily.

Results and discussion

Microsynteny-based comparative genomics demonstrated that whole genome duplication (WGD) and segmental duplication events drove bHLH genes expansion and functional divergence. Tandem duplication accelerated the generation of bHLH paralogs. Multiple sequence alignments showed that all bHLH proteins included bHLH-zip and ACT-like conserved domains. The MYC2 subfamily had a typical bHLH-MYC_N domain. The phylogenetic tree revealed the classification and putative roles of bHLHs. The analysis of cis-acting elements revealed that promoter of the majority of bHLH genes contain multiple regulatory elements relevant to light response, hormone responses, and abiotic stresses, and the bHLH genes could be activated by binding these elements. The expression profiling and qRT-PCR results indicated that bHLH subgroups IIIe and IIId may have an antagonistic role in JA-mediated expression of stress-related genes. DhbHLH20 and DhbHLH21 were considered to be the positive regulators in the early response of JA signaling, while DhbHLH24 and DhbHLH25 might be the negative regulators. Our findings may provide a practical reference for the functional study of DhbHLH genes and the regulation of secondary metabolites.

PtrVCS2 Regulates Drought Resistance by Changing Vessel Morphology and Stomatal Closure in Populus trichocarpa

Drought has severe effects on plant growth, forest productivity, and survival throughout the world. Understanding the molecular regulation of drought resistance in forest trees can enable effective strategic engineering of novel drought-resistant genotypes of tree species. In this study, we identified a gene, PtrVCS2, encoding a zinc finger (ZF) protein of the ZF-homeodomain transcription factor in Populus trichocarpa (Black Cottonwood) Torr. & A. Gray. ex Hook. Overexpression of PtrVCS2 (OE-PtrVCS2) in P. trichocarpa resulted in reduced growth, a higher proportion of smaller stem vessels, and strong drought-resistance phenotypes. Stomatal movement experiments revealed that the OE-PtrVCS2 transgenics showed lower stomata apertures than wild-type plants under drought conditions. RNA-seq analysis of the OE-PtrVCS2 transgenics showed that PtrVCS2 regulates the expression of multiple genes involved in regulation of stomatal opening and closing, particularly the PtrSULTR3;1-1 gene, and several genes related to cell wall biosynthesis, such as PtrFLA11-12 and PtrPR3-3. Moreover, we found that the water use efficiency of the OE-PtrVCS2 transgenic plants was consistently higher than that of wild type plants when subjected to chronic drought stress. Taken together, our results suggest that PtrVCS2 plays a positive role in improving drought adaptability and resistance in P. trichocarpa.

Effect of Shading on the Morphological, Physiological, and Biochemical Characteristics as Well as the Transcriptome of Matcha Green Tea

High-quality tea leaves are required for matcha production. Shading is one of the key agronomic practices that can increase the quality of green tea. The objectives among matcha tea producers include increasing the ammonia and chlorophyll contents of tea buds, decreasing tea polyphenol contents, and enhancing tea aroma formation. In this study, Fuding white tea plants were cultivated under open-air conditions (control) as well as under 85% (S85) and 95% (S95) shade. The chlorophyll contents were highest for the S85 treatment, followed by the S95 and control treatments. Moreover, shading increased the theanine and caffeine contents, while decreasing the polyphenol (epicatechin and epigallocatechin) contents, thereby optimizing matcha tea flavors. A total of 2788 differentially expressed genes (DEGs) were identified, of which 1151 and 1637 were respectively upregulated and downregulated in response to shading. The GO and KEGG enrichment analyses indicated that most of the DEGs were associated with metabolic processes (e.g., MAPK signaling, plant-pathogen interactions, and phenylpropanoid biosynthesis). Therefore, shading may modulate tea plant metabolism, signaling, biosynthetic activities, and environment-related changes to gene transcription. The expression of amino acid permeases (APP) encoding genes was downregulated in tea plants. Thus, shading influences theanine biosynthesis and the AAP-mediated distribution of theanine in tea plants.

NLRscape: an atlas of plant NLR proteins

NLRscape is a webserver that curates a collection of over 80 000 plant protein sequences identified in UniProtKB to contain NOD-like receptor signatures, and hosts in addition a number of tools aimed at the exploration of the complex sequence landscape of this class of plant proteins. Each entry gathers sequence information, domain and motif annotations from multiple third-party sources but also in-house advanced annotations aimed at addressing caveats of the existing broad-based annotations. NLRscape provides a top-down perspective of the NLR sequence landscape but also services for assisting a bottom-up approach starting from a given input sequence. Sequences are clustered by their domain organization layout, global homology and taxonomic spread-in order to allow analysis of how particular traits of an NLR family are scattered within the plant kingdom. Tools are provided for users to locate their own protein of interest in the overall NLR landscape, generate custom clusters centered around it and perform a large number of sequence and structural analyses using included interactive online instruments. Amongst these, we mention: taxonomy distribution plots, homology cluster graphs, identity matrices and interactive MSA synchronizing secondary structure and motif predictions. NLRscape can be found at: https://nlrscape.biochim.ro/.

Transcriptome Meta-Analysis Identifies Candidate Hub Genes and Pathways of Pathogen Stress Responses in Arabidopsis thaliana

Following a pathogen attack, plants defend themselves using multiple defense mechanisms to prevent infections. We used a meta-analysis and systems-biology analysis to search for general molecular plant defense responses from transcriptomic data reported from different pathogen attacks in Arabidopsis thaliana. Data from seven studies were subjected to meta-analysis, which revealed a total of 3694 differentially expressed genes (DEGs), where both healthy and infected plants were considered. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analysis further suggested that the DEGs were involved in several biosynthetic metabolic pathways, including those responsible for the biosynthesis of secondary metabolites and pathways central to photosynthesis and plant-pathogen interactions. Using network analysis, we highlight the importance of WRKY40, WRKY46 and STZ, and suggest that they serve as major points in protein-protein interactions. This is especially true regarding networks of composite-metabolic responses by pathogens. In summary, this research provides a new approach that illuminates how different mechanisms of transcriptome responses can be activated in plants under pathogen infection and indicates that common genes vary in their ability to regulate plant responses to the pathogens studied herein.

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.

Efficiency of microbial bio-agents as elicitors in plant defense mechanism under biotic stress: A review

Numerous harmful microorganisms and insect pests have the ability to cause plant infections or damage, which is mostly controlled by toxic chemical agents. These chemical compounds and their derivatives exhibit hazardous effects on habitats and human life too. Hence, there's a need to develop novel, more effective and safe bio-control agents. A variety of microbes such as viruses, bacteria, and fungi possess a great potential to fight against phytopathogens and thus can be used as bio-control agents instead of harmful chemical compounds. These naturally occurring microorganisms are applied to the plants in order to control phytopathogens. Moreover, practicing them appropriately for agriculture management can be a way towards a sustainable approach. The MBCAs follow various modes of action and act as elicitors where they induce a signal to activate plant defense mechanisms against a variety of pathogens. MBCAs control phytopathogens and help in disease suppression through the production of enzymes, antimicrobial compounds, antagonist activity involving hyper-parasitism, induced resistance, competitive inhibition, etc. Efficient recognition of pathogens and prompt defensive response are key factors of induced resistance in plants. This resistance phenomenon is pertaining to a complex cascade that involves an increased amount of defensive proteins, salicylic acid (SA), or induction of signaling pathways dependent on plant hormones. Although, there's a dearth of information about the exact mechanism of plant-induced resistance, the studies conducted at the physiological, biochemical and genetic levels. These studies tried to explain a series of plant defensive responses triggered by bio-control agents that may enhance the defensive capacity of plants. Several natural and recombinant microorganisms are commercially available as bio-control agents that mainly include strains of Bacillus, Pseudomonads and Trichoderma. However, the complete understanding of microbial bio-control agents and their interactions at cellular and molecular levels will facilitate the screening of effective and eco-friendly bio-agents, thereby increasing the scope of MBCAs. This article is a comprehensive review that highlights the importance of microbial agents as elicitors in the activation and regulation of plant defense mechanisms in response to a variety of pathogens.

PRGdb 4.0: an updated database dedicated to genes involved in plant disease resistance process

The Plant Resistance Genes database (PRGdb; http://prgdb.org/prgdb4/) has been greatly expanded, keeping pace with the increasing amount of available knowledge and data (sequenced proteomes, cloned genes, public analysis data, etc.). The easy-to-use style of the database website has been maintained, while an updated prediction tool, more data and a new section have been added. This new section will contain plant resistance transcriptomic experiments, providing additional easy-to-access experimental information. DRAGO3, the tool for automatic annotation and prediction of plant resistance genes behind PRGdb, has been improved in both accuracy and sensitivity, leading to more reliable predictions. PRGdb offers 199 reference resistance genes and 586.652 putative resistance genes from 182 sequenced proteomes. Compared to the previous release, PRGdb 4.0 has increased the number of reference resistance genes from 153 to 199, the number of putative resistance genes from 177K from 76 proteomes to 586K from 182 sequenced proteomes. A new section has been created that collects plant-pathogen transcriptomic data for five species of agricultural interest. Thereby, with these improvements and data expansions, PRGdb 4.0 aims to serve as a reference to the plant scientific community and breeders worldwide, helping to further study plant resistance mechanisms that contribute to fighting pathogens.

Transcriptomic Analysis Reveals Candidate Genes Responding Maize Gray Leaf Spot Caused by Cercospora zeina

Gray leaf spot (GLS), caused by the fungal pathogen Cercospora zeina (C. zeina), is one of the most destructive soil-borne diseases in maize (Zea mays L.), and severely reduces maize production in Southwest China. However, the mechanism of resistance to GLS is not clear and few resistant alleles have been identified. Two maize inbred lines, which were shown to be resistant (R6) and susceptible (S8) to GLS, were injected by C. zeina spore suspensions. Transcriptome analysis was carried out with leaf tissue at 0, 6, 24, 144, and 240 h after inoculation. Compared with 0 h of inoculation, a total of 667 and 419 stable common differentially expressed genes (DEGs) were found in the resistant and susceptible lines across the four timepoints, respectively. The DEGs were usually enriched in 'response to stimulus' and 'response to stress' in GO term analysis, and 'plant-pathogen interaction', 'MAPK signaling pathways', and 'plant hormone signal transduction' pathways, which were related to maize's response to GLS, were enriched in KEGG analysis. Weighted-Genes Co-expression Network Analysis (WGCNA) identified two modules, while twenty hub genes identified from these indicated that plant hormone signaling, calcium signaling pathways, and transcription factors played a central role in GLS sensing and response. Combing DEGs and QTL mapping, five genes were identified as the consensus genes for the resistance of GLS. Two genes, were both putative Leucine-rich repeat protein kinase family proteins, specifically expressed in R6. In summary, our results can provide resources for gene mining and exploring the mechanism of resistance to GLS in maize.

Genome, Transcriptome, and Germplasm Sequencing Uncovers Functional Variation in the Warm-Season Grain Legume Horsegram Macrotyloma uniflorum (Lam.) Verdc

Horsegram is a grain legume with excellent nutritional and remedial properties and good climate resilience, able to adapt to harsh environmental conditions. Here, we used a combination of short- and long-read sequencing technologies to generate a genome sequence of 279.12Mb, covering 83.53% of the estimated total size of the horsegram genome, and we annotated 24,521 genes. De novo prediction of DNA repeats showed that approximately 25.04% of the horsegram genome was made up of repetitive sequences, the lowest among the legume genomes sequenced so far. The major transcription factors identified in the horsegram genome were bHLH, ERF, C2H2, WRKY, NAC, MYB, and bZIP, suggesting that horsegram is resistant to drought. Interestingly, the genome is abundant in Bowman-Birk protease inhibitors (BBIs), which can be used as a functional food ingredient. The results of maximum likelihood phylogenetic and estimated synonymous substitution analyses suggested that horsegram is closely related to the common bean and diverged approximately 10.17 million years ago. The double-digested restriction associated DNA (ddRAD) sequencing of 40 germplasms allowed us to identify 3,942 high-quality SNPs in the horsegram genome. A genome-wide association study with powdery mildew identified 10 significant associations similar to the MLO and RPW8.2 genes. The reference genome and other genomic information presented in this study will be of great value to horsegram breeding programs. In addition, keeping the increasing demand for food with nutraceutical values in view, these genomic data provide opportunities to explore the possibility of horsegram for use as a source of food and nutraceuticals.

The regulatory role of gibberellin related genes DKGA2ox1 and MIR171f_3 in persimmon dwarfism

'Nantongxiaofangshi' (Diospyros kaki Thunb., D. kaki Thunb.) is a local cultivar of persimmon with dwarf-like traits in Jiangsu, China. Closely spaced planting afforded by dwarfism is usually one of the most important ways to promote fruit cultivation and production. However, the understanding of dwarfism in D. kaki Thunb. is very limited at the molecular level, which hinders the further increase of the fruit production. In this work, a persimmon transgenic system was successfully established, and the field experiments of grafting phenotype were carried out. The results showed that D. kaki Thunb. could be used as an interstock to induce dwarfing in grafted scions, and the dwarf character was better when interstock lengths were between 20 and 25 cm. Furthermore, the key genes related to dwarfism in D. kaki Thunb. were screened and verified, and subsequently, the regulatory role of related genes in persimmon dwarfism was figured out. It was found that the gene encoding gibberellin 2-oxidase-1 (DkGA2ox1) involved in GA biosynthesis was associated with the dwarfing in D. kaki Thunb. Overexpression of DkGA2ox1 in Diospyros lotus resulted in a typical dwarf phenotype. Meanwhile, the microRNA data showed that the miR171f_3 demonstrated the active involvement in GA pathway response in persimmon dwarfism. DkGA2ox1 and MIR171f_3, as two highly expressed genes in D. kaki Thunb. interstock, could be used as stimulus signals to affect the content of GA in scion, however, the specific transmission mechanism still needs to be further explored. Ultimately, the bioactive GA level was decreased, resulting in the scion dwarfism.

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.

Effects of graphene on morphology, microstructure and transcriptomic profiling of Pinus tabuliformis Carr. roots

Graphene has shown great potential for improving growth of many plants, but its effect on woody plants remains essentially unstudied. In this work, Pinus tabuliformis Carr. bare-rooted seedlings grown outdoors in pots were irrigated with a graphene solution over a concentration range of 0-50 mg/L for six months. Graphene was found to stimulate root growth, with a maximal effect at 25 mg/L. We then investigated root microstructure and carried out transcript profiling of root materials treated with 0 and 25 mg/L graphene. Graphene treatment resulted in plasma-wall separation and destruction of membrane integrity in root cells. More than 50 thousand of differentially expressed genes (DEGs) were obtained by RNA sequencing, among which 6477 could be annotated using other plant databases. The GO enrichment analysis and KEGG pathway analysis of the annotated DEGs indicated that abiotic stress responses, which resemble salt stress, were induced by graphene treatment in roots, while responses to biotic stimuli were inhibited. Numerous metabolic processes and hormone signal transduction pathways were altered by the treatment. The growth promotion effects of graphene may be mediated by encouraging proline synthesis, and suppression of the expression of the auxin response gene SMALL AUXIN UP-REGULATED RNA 41 (SAUR41), PYL genes which encode ABA receptors, and GSK3 homologs.

Draft genome sequence of the pulse crop blackgram [Vigna mungo (L.) Hepper] reveals potential R-genes

Blackgram [Vigna mungo (L.) Hepper] (2n = 2x = 22), an important Asiatic legume crop, is a major source of dietary protein for the predominantly vegetarian population. Here we construct a draft genome sequence of blackgram, for the first time, by employing hybrid genome assembly with Illumina reads and third generation Oxford Nanopore sequencing technology. The final de novo whole genome of blackgram is ~ 475 Mb (82% of the genome) and has maximum scaffold length of 6.3 Mb with scaffold N50 of 1.42 Mb. Genome analysis identified 42,115 genes with mean coding sequence length of 1131 bp. Around 80.6% of predicted genes were annotated. Nearly half of the assembled sequence is composed of repetitive elements with retrotransposons as major (47.3% of genome) transposable elements, whereas, DNA transposons made up only 2.29% of the genome. A total of 166,014 SSRs, including 65,180 compound SSRs, were identified and primer pairs for 34,816 SSRs were designed. Out of the 33,959 proteins, 1659 proteins showed presence of R-gene related domains. KIN class was found in majority of the proteins (905) followed by RLK (239) and RLP (188). The genome sequence of blackgram will facilitate identification of agronomically important genes and accelerate the genetic improvement of blackgram.

The dual interplay of RAV5 in activating nitrate reductases and repressing catalase activity to improve disease resistance in cassava

Cassava bacterial blight (CBB) caused by Xanthomonas axonopodis pv. manihotis (Xam) seriously affects cassava yield. Nitrate reductase (NR) plays an important role in plant nitrogen metabolism in plants. However, the in vivo role of NR and the corresponding signalling pathway remain unclear in cassava. In this study, we isolated MeNR1/2 and revealed their novel upstream transcription factor MeRAV5. We also identified MeCatalase1 (MeCAT1) as the interacting protein of MeRAV5. In addition, we investigated the role of MeCatalase1 and MeRAV5-MeNR1/2 module in cassava defence response. MeNRs positively regulates cassava disease resistance against CBB through modulation of nitric oxide (NO) and extensive transcriptional reprogramming especially in mitogen-activated protein kinase (MAPK) signalling. Notably, MeRAV5 positively regulates cassava disease resistance through the coordination of NO and hydrogen peroxide (H2 O2 ) level. On the one hand, MeRAV5 directly activates the transcripts of MeNRs and NO level by binding to CAACA motif in the promoters of MeNRs. On the other hand, MeRAV5 interacts with MeCAT1 to inhibit its activity, so as to negatively regulate endogenous H2 O2 level. This study highlights the precise coordination of NR activity and CAT activity by MeRAV5 through directly activating MeNRs and interacting with MeCAT1 in plant immunity.

Comparative Transcriptome and Expression Profiling of Resistant and Susceptible Banana Cultivars during Infection by Fusarium oxysporum

Fusarium wilt caused by Fusarium oxysporum f. sp. cubense (Foc) is one of the most destructive diseases of banana. Methods to control the disease are still inadequate. The present investigation targeted expression of defense-related genes in tissue cultured banana plantlets of Fusarium resistant and susceptible cultivars after infection with biological control agents (BCAs) and Fusarium (Foc race 1). In total 3034 differentially expressed genes were identified which annotated to 58 transcriptional families (TF). TF families such as MYB, bHLH and NAC TFs were mostly up-regulated in response to pathogen stress, whereas AP2/EREBP were mostly down-regulated. Most genes were associated with plant-pathogen response, plant hormone signal transduction, starch and sucrose metabolism, cysteine and methionine metabolism, flavonoid biosynthesis, selenocompound metabolism, phenylpropanoid biosynthesis, mRNA surveillance pathway, mannose type O-glycan biosynthesis, amino acid and nucleotide sugar metabolism, cyanoamino acid metabolism, and hormone signal transduction. Our results showed that the defense mechanisms of resistant and susceptible banana cultivars treated with BCAs, were regulated by differentially expressed genes in various categories of defense pathways. Furthermore, the association with different resistant levels might serve as a strong foundation for the control of Fusarium wilt of banana.

Insertion of Badnaviral DNA in the Late Blight Resistance Gene (R1a) of Brinjal Eggplant (Solanum melongena)

Endogenous viral elements (EVEs) are widespread in plant genomes. They result from the random integration of viral sequences into host plant genomes by horizontal DNA transfer and have the potential to alter host gene expression. We performed a large-scale search for co-transcripts including caulimovirid and plant sequences in 1,678 plant and 230 algal species and characterized 50 co-transcripts in 45 distinct plant species belonging to lycophytes, ferns, gymnosperms and angiosperms. We found that insertion of badnavirus EVEs along with Ty-1 copia mobile elements occurred into a late blight resistance gene (R1) of brinjal eggplant (Solanum melongena) and wild relatives in genus Solanum and disrupted R1 orthologs. EVEs of two previously unreported badnaviruses were identified in the genome of S. melongena, whereas EVEs from an additional novel badnavirus were identified in the genome of S. aethiopicum, the cultivated scarlet eggplant. Insertion of these viruses in the ancestral lineages of the direct wild relatives of the eggplant would have occurred during the last 3 Myr, further supporting the distinctiveness of the group of the eggplant within the giant genus Solanum.

Transcriptome analysis of diploid and triploid Populus tomentosa

Triploid Chinese white poplar (Populus tomentosa Carr., Salicaceae) has stronger advantages in growth and better stress resistance and wood quality than diploid P. tomentosa. Using transcriptome sequencing technology to identify candidate transcriptome-based markers for growth vigor in young tree tissue is of great significance for the breeding of P. tomentosa varieties in the future. In this study, the cuttings of diploid and triploid P. tomentosa were used as plant materials, transcriptome sequencing was carried out, and their tissue culture materials were used for RT-qPCR verification of the expression of genes. The results showed that 12,240 differentially expressed genes in diploid and triploid P. tomentosa transcripts were annotated and enriched into 135 metabolic pathways. The top six pathways that enriched the most significantly different genes were plant-pathogen interaction, phenylpropanoid biosynthesis, MAPK signalling pathway-plant, ascorbate and aldarate metabolism, diterpenoid biosynthesis, and the betalain biosynthesis pathway. Ten growth-related genes were selected from pathways of plant hormone signal transduction and carbon fixation in photosynthetic organisms for RT-qPCR verification. The expression levels of MDH and CYCD3 in tissue-cultured and greenhouse planted triploid P. tomentosa were higher than those in tissue-cultured diploid P. tomentosa, which was consist ent with the TMM values calculated by transcriptome.

Transcriptome analysis provides insights into the mechanisms underlying wheat cultivar Shumai126 responding to stripe rust

Stripe rust, caused by Puccinia striiformis f. sp. tritici (Pst), is a destructive fungal disease of wheat globally. Breeding resistance cultivars is one of the most cost-effective methods to control Pst. Shumai126 (SM126), a high-yielding commercial wheat cultivar, showed strong stripe rust resistance for more than ten years. However, the molecular mechanisms and the responsive genes underlying the SM126 resistance to Pst have not been explored yet. In the present study, RNA-seq was used to analyze changes in the transcriptome at different time points of Pst infection in seedling leaves of SM126. In total, 520, 148 and 1439 differentially expressed genes (DEGs) were found to be up- or down-regulated after Pst infection at 1, 3, and 7 days post inoculation, respectively. The majority of DEGs exhibited transient expression patterns during Pst infection at different time points. GO and KEGG enrichment analysis revealed that many biological processes, such as photosynthesis, flavonoid biosynthesis, oxidative phosphorylation, MAPK signaling pathway, and phenylalanine metabolism are involved in SM126 response to Pst. Expression of genes involved in the plant-pathogen interaction pathway was detected and some key genes showed differential expression. DEGs encoded R proteins and transcription factors were also identified. Our study suggests the gene resources in SM126 related to stripe rust response could be valuable for understanding the mechanisms involved in stripe rust resistance and improvement of wheat resistance to Pst.

Resistance Genes Affect How Pathogens Maintain Plant Abundance and Diversity

AbstractSpecialized pathogens are thought to maintain plant community diversity; however, most ecological studies treat pathogens as a black box. Here we develop a theoretical model to test how the impact of specialized pathogens changes when plant resistance genes (R-genes) mediate susceptibility. This work synthesizes two major hypotheses: the gene-for-gene model of pathogen resistance and the Janzen-Connell hypothesis of pathogen-mediated coexistence. We examine three scenarios. First, R-genes do not affect seedling survival; in this case, pathogens promote diversity. Second, seedlings are protected from pathogens when their R-gene alleles and susceptibility differ from those of nearby conspecific adults, thereby reducing transmission. If resistance is not costly, pathogens are less able to promote diversity because populations with low R-gene diversity suffer higher mortality, putting those populations at a disadvantage and potentially causing their exclusion. R-gene diversity may also be reduced during population bottlenecks, creating a priority effect. Third, when R-genes affect survival but resistance is costly, populations can avoid extinction by losing resistance alleles, as they cease paying a cost that is unneeded. Thus, the impact pathogens can have on tree diversity depends on the mechanism of plant-pathogen interactions. Future empirical studies should examine which of these scenarios most closely reflects the real world.

Computational identification of receptor-like kinases "RLK" and receptor-like proteins "RLP" in legumes

Background

In plants, the plasma membrane is enclosed by the cell wall and anchors RLK and RLP proteins, which play a fundamental role in perception of developmental and environmental cues and are crucial in plant development and immunity. These plasma membrane receptors belong to large gene/protein families that are not easily classified computationally. This detailed analysis of these plasma membrane proteins brings a new source of information to the legume genetic, physiology and breeding research communities.

Results

A computational approach to identify and classify RLK and RLP proteins is presented. The strategy was evaluated using experimentally-validated RLK and RLP proteins and was determined to have a sensitivity of over 0.85, a specificity of 1.00, and a Matthews correlation coefficient of 0.91. The computational approach can be used to develop a detailed catalog of plasma membrane receptors (by type and domains) in several legume/crop species. The exclusive domains identified in legumes for RLKs are WaaY, APH Pkinase_C, LRR_2, and EGF, and for RLP are L-lectin LPRY and PAN_4. The RLK-nonRD and RLCK subclasses are also discovered by the methodology. In both classes, less than 20% of the total RLK predicted for each species belong to this class. Among the 10-species evaluated ~ 40% of the proteins in the kinome are RLKs. The exclusive legume domain combinations identified are B-Lectin/PR5K domains in G. max, M. truncatula, V. angularis, and V. unguiculata and a three-domain combination B-lectin/S-locus/WAK in C. cajan, M. truncatula, P. vulgaris, V. angularis. and V. unguiculata.

Conclusions

The analysis suggests that about 2% of the proteins of each genome belong to the RLK family and less than 1% belong to RLP family. Domain diversity combinations are greater for RLKs compared with the RLP proteins and LRR domains, and the dual domain combination LRR/Malectin were the most frequent domain for both groups of plasma membrane receptors among legume and non-legume species. Legumes exclusively show Pkinase extracellular domains, and atypical domain combinations in RLK and RLP compared with the non-legumes evaluated. The computational logic approach is statistically well supported and can be used with the proteomes of other plant species.

Discovery of two novel and adjacent QTLs on chromosome B02 controlling resistance against bacterial wilt in peanut variety Zhonghua 6

Two novel and adjacent genomics and candidate genes for bacterial wilt resistance were identified on chromosome B02 in peanut variety Zhonghua 6 using both traditional QTL mapping and QTL-seq methods. Peanut (Arachis hypogaea) is an important oilseed crop worldwide. Utilization of genetic resistance is the most economic and effective approach to control bacterial wilt, one of the most devastating plant diseases, in peanut production. To accelerate the genetic improvement of bacterial wilt resistance (BWR) in peanut breeding programs, quantitative trait locus (QTL) mapping has been conducted for two resistant varieties. In this context, we deployed linkage mapping as well as sequencing-based mapping approach, QTL-seq, to identify genomic regions and candidate genes for BWR in another highly resistant variety Zhonghua 6. The recombination inbred line population (268 progenies) from the cross Xuhua 13 × Zhonghua 6 was used in BWR evaluation across five environments. QTL mapping using both SSR- and SNP-based genetic maps identified a stable QTL (qBWRB02-1) on chromosome B02 with 37.79-78.86% phenotypic variation explained (PVE) across five environments. The QTL-seq facilitated further dissection of qBWRB02-1 into two adjacent genomic regions, qBWRB02-1-1 (2.81-4.24 Mb) and qBWRB02-1-2 (6.54-8.75 Mb). Mapping of newly developed Kompetitive allele-specific PCR (KASP) markers on the genetic map confirmed their stable expressions across five environments. The effects of qBWRB02-1-1 (49.43-68.86% PVE) were much higher than qBWRB02-1-2 (3.96-6.48% PVE) and other previously reported QTLs. Nineteen putative candidate genes affected by 49 non-synonymous SNPs were identified for qBWRB02-1-1, and ten of them were predicted to code for disease resistance proteins. The major and stable QTL qBWRB02-1-1 and validated KASP markers could be deployed in genomics-assisted breeding (GAB) to develop improved peanut varieties with enhanced BWR.

Comparison of leaf transcriptome in response to Rhizoctonia solani infection between resistant and susceptible rice cultivars

Background

Sheath blight (SB), caused by Rhizoctonia solani, is a common rice disease worldwide. Currently, rice cultivars with robust resistance to R. solani are still lacking. To provide theoretic basis for molecular breeding of R. solani-resistant rice cultivars, the changes of transcriptome profiles in response to R. solani infection were compared between a moderate resistant cultivar (Yanhui-888, YH) and a susceptible cultivar (Jingang-30, JG).

Results

In the present study, 3085 differentially express genes (DEGs) were detected between the infected leaves and the control in JG, with 2853 DEGs in YH. A total of 4091 unigenes were significantly upregulated in YH than in JG before infection, while 3192 were significantly upregulated after infection. Further analysis revealed that YH and JG showed similar molecular responses to R. solani infection, but the responses were earlier in JG than in YH. Expression levels of trans-cinnamate 4-monooxygenase (C4H), ethylene-insensitive protein 2 (EIN2), transcriptome factor WRKY33 and the KEGG pathway plant-pathogen interaction were significantly affected by R. solani infection. More importantly, these components were all over-represented in YH cultivar than in JG cultivar before and/or after infection.

Conclusions

These genes possibly contribute to the higher resistance of YH to R. solani than JG and were potential target genes to molecularly breed R. solani-resistant rice cultivar.

Transcriptome Analysis and Identification of Genes Associated with Starch Metabolism in Castanea henryi Seed (Fagaceae)

Starch is the most important form of carbohydrate storage and is the major energy reserve in some seeds, especially Castanea henryi. Seed germination is the beginning of the plant's life cycle, and starch metabolism is important for seed germination. As a complex metabolic pathway, the regulation of starch metabolism in C. henryi is still poorly understood. To explore the mechanism of starch metabolism during the germination of C. henryi, we conducted a comparative gene expression analysis at the transcriptional level using RNA-seq across four different germination stages, and analyzed the changes in the starch and soluble sugar contents. The results showed that the starch content increased in 0-10 days and decreased in 10-35 days, while the soluble sugar content continuously decreased in 0-30 days and increased in 30-35 days. We identified 49 candidate genes that may be associated with starch and sucrose metabolism. Three ADP-glucose pyrophosphorylase (AGPase) genes, two nucleotide pyrophosphatase/phosphodiesterases (NPPS) genes and three starch synthases (SS) genes may be related to starch accumulation. Quantitative real-time polymerase chain reaction (qRT-PCR) was used to validate the expression levels of these genes. Our study combined transcriptome data with physiological and biochemical data, revealing potential candidate genes that affect starch metabolism during seed germination, and provides important data about starch metabolism and seed germination in seed plants.

Deployment of Genetic and Genomic Tools Toward Gaining a Better Understanding of Rice-Xanthomonas oryzae pv. oryzae Interactions for Development of Durable Bacterial Blight Resistant Rice

Rice is the most important food crop worldwide and sustainable rice production is important for ensuring global food security. Biotic stresses limit rice production significantly and among them, bacterial blight (BB) disease caused by Xanthomonas oryzae pv. oryzae (Xoo) is very important. BB reduces rice yields severely in the highly productive irrigated and rainfed lowland ecosystems and in recent years; the disease is spreading fast to other rice growing ecosystems as well. Being a vascular pathogen, Xoo interferes with a range of physiological and biochemical exchange processes in rice. The response of rice to Xoo involves specific interactions between resistance (R) genes of rice and avirulence (Avr) genes of Xoo, covering most of the resistance genes except the recessive ones. The genetic basis of resistance to BB in rice has been studied intensively, and at least 44 genes conferring resistance to BB have been identified, and many resistant rice cultivars and hybrids have been developed and released worldwide. However, the existence and emergence of new virulent isolates of Xoo in the realm of a rapidly changing climate necessitates identification of novel broad-spectrum resistance genes and intensification of gene-deployment strategies. This review discusses about the origin and occurrence of BB in rice, interactions between Xoo and rice, the important roles of resistance genes in plant's defense response, the contribution of rice resistance genes toward development of disease resistance varieties, identification and characterization of novel, and broad-spectrum BB resistance genes from wild species of Oryza and also presents a perspective on potential strategies to achieve the goal of sustainable disease management.

Rice Stress-Resistant SNP Database

BACKGROUND

Rice (Oryza sativa L.) yield is limited inherently by environmental stresses, including biotic and abiotic stresses. Thus, it is of great importance to perform in-depth explorations on the genes that are closely associated with the stress-resistant traits in rice. The existing rice SNP databases have made considerable contributions to rice genomic variation information but none of them have a particular focus on integrating stress-resistant variation and related phenotype data into one web resource.

RESULTS

Rice Stress-Resistant SNP database (http://bioinformatics.fafu.edu.cn/RSRS) mainly focuses on SNPs specific to biotic and abiotic stress-resistant ability in rice, and presents them in a unified web resource platform. The Rice Stress-Resistant SNP (RSRS) database contains over 9.5 million stress-resistant SNPs and 797 stress-resistant candidate genes in rice, which were detected from more than 400 stress-resistant rice varieties. We incorporated the SNPs function, genome annotation and phenotype information into this database. Besides, the database has a user-friendly web interface for users to query, browse and visualize a specific SNP efficiently. RSRS database allows users to query the SNP information and their relevant annotations for individual variety or more varieties. The search results can be visualized graphically in a genome browser or displayed in formatted tables. Users can also align SNPs between two or more rice accessions.

CONCLUSION

RSRS database shows great utility for scientists to further characterize the function of variants related to environmental stress-resistant ability in rice.

Next-generation sequencing identified genomic region and diagnostic markers for resistance to bacterial wilt on chromosome B02 in peanut (Arachis hypogaea L.)

Bacterial wilt, caused by Ralstonia solanacearum, is a devastating disease affecting over 350 plant species. A few peanut cultivars were found to possess stable and durable bacterial wilt resistance (BWR). Genomics-assisted breeding can accelerate the process of developing resistant cultivars by using diagnostic markers. Here, we deployed sequencing-based trait mapping approach, QTL-seq, to discover genomic regions, candidate genes and diagnostic markers for BWR in a recombination inbred line population (195 progenies) of peanut. The QTL-seq analysis identified one candidate genomic region on chromosome B02 significantly associated with BWR. Mapping of newly developed single nucleotide polymorphism (SNP) markers narrowed down the region to 2.07 Mb and confirmed its major effects and stable expressions across three environments. This candidate genomic region had 49 nonsynonymous SNPs affecting 19 putative candidate genes including seven putative resistance genes (R-genes). Two diagnostic markers were successfully validated in diverse breeding lines and cultivars and could be deployed in genomics-assisted breeding of varieties with enhanced BWR.

Genome survey of resistance gene analogs in sugarcane: genomic features and differential expression of the innate immune system from a smut-resistant genotype

BACKGROUND

Resistance genes composing the two-layer immune system of plants are thought as important markers for breeding pathogen-resistant crops. Many have been the attempts to establish relationships between the genomic content of Resistance Gene Analogs (RGAs) of modern sugarcane cultivars to its degrees of resistance to diseases such as smut. However, due to the highly polyploid and heterozygous nature of sugarcane genome, large scale RGA predictions is challenging.

RESULTS

We predicted, searched for orthologs, and investigated the genomic features of RGAs within a recently released sugarcane elite cultivar genome, alongside the genomes of sorghum, one sugarcane ancestor (Saccharum spontaneum), and a collection of de novo transcripts generated for six modern cultivars. In addition, transcriptomes from two sugarcane genotypes were obtained to investigate the roles of RGAs differentially expressed (RGADE) in their distinct degrees of resistance to smut. Sugarcane references lack RGAs from the TNL class (Toll-Interleukin receptor (TIR) domain associated to nucleotide-binding site (NBS) and leucine-rich repeat (LRR) domains) and harbor elevated content of membrane-associated RGAs. Up to 39% of RGAs were organized in clusters, and 40% of those clusters shared synteny. Basically, 79% of predicted NBS-encoding genes are located in a few chromosomes. S. spontaneum chromosome 5 harbors most RGADE orthologs responsive to smut in modern sugarcane. Resistant sugarcane had an increased number of RGAs differentially expressed from both classes of RLK (receptor-like kinase) and RLP (receptor-like protein) as compared to the smut-susceptible. Tandem duplications have largely contributed to the expansion of both RGA clusters and the predicted clades of RGADEs.

CONCLUSIONS

Most of smut-responsive RGAs in modern sugarcane were potentially originated in chromosome 5 of the ancestral S. spontaneum genotype. Smut resistant and susceptible genotypes of sugarcane have a distinct pattern of RGADE. TM-LRR (transmembrane domains followed by LRR) family was the most responsive to the early moment of pathogen infection in the resistant genotype, suggesting the relevance of an innate immune system. This work can help to outline strategies for further understanding of allele and paralog expression of RGAs in sugarcane, and the results should help to develop a more applied procedure for the selection of resistant plants in sugarcane.

Machine learning approaches and their current application in plant molecular biology: A systematic review

Machine learning (ML) is a field of artificial intelligence that has rapidly emerged in molecular biology, thus allowing the exploitation of Big Data concepts in plant genomics. In this context, the main challenges are given in terms of how to analyze massive datasets and extract new knowledge in all levels of cellular systems research. In summary, ML techniques allow complex interactions to be inferred in several biological systems. Despite its potential, ML has been underused due to complex computational algorithms and definition terms. Therefore, a systematic review to disentangle ML approaches is relevant for plant scientists and has been considered in this study. We presented the main steps for ML development (from data selection to evaluation of classification/prediction models) with a respective discussion approaching functional genomics mainly in terms of pathogen effector genes in plant immunity. Additionally, we also considered how to access public source databases under an ML framework towards advancing plant molecular biology and introduced novel powerful tools, such as deep learning.

Genetic signatures of plant resistance genes with known function within and between species

Plant disease resistance (R) genes have undergone significant evolutionary divergence to cope with rapid changes in pathogens. These highly variable evolutionary patterns may have contributed to diversity in R gene protein families or structures. Here, the evolutionary patterns of 76 identified R genes and their homologs were investigated within and between plant species. Results demonstrated that nucleotide binding sites and leucine-rich-repeat genes located in loci with complex evolutionary histories tended to evolve rapidly, have high variation in copy numbers, exhibit high levels of nucleotide variation and frequent gene conversion events, and also exhibit high non-synonymous to synonymous substitution ratios in LRR regions. However, non-NBS-LRR R genes are relatively well conserved with constrained variation and are more likely to participate in the basic defense system of hosts. In addition, both conserved and highly divergent evolutionary patterns were observed for the same R genes and were consistent with inter- and intra-specific distributions of some R genes. These results thus indicate either continuous or altered evolutionary patterns between and within species. The present investigation is the first attempt to investigate evolutionary patterns among all clearly functional R genes. The results reported here thus provide a foundation for future plant disease studies.

Genome Sequencing and Analysis of the Peanut B-Genome Progenitor (Arachis ipaensis)

Peanut (Arachis hypogaea L.), an important leguminous crop, is widely cultivated in tropical and subtropical regions. Peanut is an allotetraploid, having A and B subgenomes that maybe have originated in its diploid progenitors Arachis duranensis (A-genome) and Arachis ipaensis (B-genome), respectively. We previously sequenced the former and here present the draft genome of the latter, expanding our knowledge of the unique biology of Arachis. The assembled genome of A. ipaensis is ~1.39 Gb with 39,704 predicted protein-encoding genes. A gene family analysis revealed that the FAR1 family may be involved in regulating peanut special fruit development. Genomic evolutionary analyses estimated that the two progenitors diverged ~3.3 million years ago and suggested that A. ipaensis experienced a whole-genome duplication event after the divergence of Glycine max. We identified a set of disease resistance-related genes and candidate genes for biological nitrogen fixation. In particular, two and four homologous genes that may be involved in the regulation of nodule development were obtained from A. ipaensis and A. duranensis, respectively. We outline a comprehensive network involved in drought adaptation. Additionally, we analyzed the metabolic pathways involved in oil biosynthesis and found genes related to fatty acid and triacylglycerol synthesis. Importantly, three new FAD2 homologous genes were identified from A. ipaensis and one was completely homologous at the amino acid level with FAD2 from A. hypogaea. The availability of the A. ipaensis and A. duranensis genomic assemblies will advance our knowledge of the peanut genome.

Red rot resistant gene characterization using RGAP markers among sugarcane cultivars resistant and susceptible to the red rot disease

Sugarcane is the major source of sugar in Asia and Europe, grown primarily in the tropical and sub-tropical zones of the world. The main disease responsible for its low yield is red rot. Therefore, in the present study, characterization of red rot disease was performed among 55 different sugarcane cultivars varying in red rot resistance level. 18 fragments were found to be associated with red rot resistance and were identified as resistant specific markers. The resistant specific fragments were amplified by RGA169, RGA396, RGA129, RGA231, RGA251, RGA057, RGA118, RGA152, RGA327, RGA542, RGA012, RGA173, RGA184, RGA275, RGA019, RGA267, RGA281 and RGA533. 7 fragments were found to be associated with red rot susceptibility and were considered as susceptible specific markers amplified by RGA088, RGA162, RGA396, RGA231, RGA251, RGA087 and RGA275. Sequencing of five resistant fragments, viz., RGA169, RGA231, RGA251, RGA267 and RGA533 was performed and the data thus obtained showed 80-99% similarity when compared with other resistant gene sequences previously submitted in NCBI database.

The pomegranate (Punica granatum L.) genome and the genomics of punicalagin biosynthesis

Pomegranate (Punica granatum L.) is a perennial fruit crop grown since ancient times that has been planted worldwide and is known for its functional metabolites, particularly punicalagins. We have sequenced and assembled the pomegranate genome with 328 Mb anchored into nine pseudo-chromosomes and annotated 29 229 gene models. A Myrtales lineage-specific whole-genome duplication event was detected that occurred in the common ancestor before the divergence of pomegranate and Eucalyptus. Repetitive sequences accounted for 46.1% of the assembled genome. We found that the integument development gene INNER NO OUTER (INO) was under positive selection and potentially contributed to the development of the fleshy outer layer of the seed coat, an edible part of pomegranate fruit. The genes encoding the enzymes for synthesis and degradation of lignin, hemicelluloses and cellulose were also differentially expressed between soft- and hard-seeded varieties, reflecting differences in their accumulation in cultivars differing in seed hardness. Candidate genes for punicalagin biosynthesis were identified and their expression patterns indicated that gallic acid synthesis in tissues could follow different biochemical pathways. The genome sequence of pomegranate provides a valuable resource for the dissection of many biological and biochemical traits and also provides important insights for the acceleration of breeding. Elucidation of the biochemical pathway(s) involved in punicalagin biosynthesis could assist breeding efforts to increase production of this bioactive compound.

Transcriptional Responses of Creeping Bentgrass to 2,3-Butanediol, a Bacterial Volatile Compound (BVC) Analogue

Bacterial volatile compounds (BVCs) have been reported to enhance plant growth and elicit plant defenses against fungal infection and insect damage. The objective of this study was to determine transcriptomic changes in response to synthetic BVC that could be associated with plant resistance to Rhizoctonia solani in creeping bentgrass. The 2,3-butanediol (BD) (250 µM) was sprayed on creeping bentgrass leaves grown in jam jars. The result showed that synthetic BD induced plant defense against R. solani for creeping bentgrass. Transcriptomic analysis demonstrated that more genes were repressed by BD while less showed up-regulation. BD suppressed the expression of some regular stress-related genes in creeping bentgrass, such as pheromone activity, calcium channel activity, photosystem II oxygen evolving complex, and hydrolase activity, while up-regulated defense related transcription factors (TFs), such as basic helix-loop-helix (bHLH) TFs, cysteine2-cysteine2-contans-like (C2C2-CO) and no apical meristem TFs (NAC). Other genes related to disease resistance, such as jasmonic acid (JA) signaling, leucine rich repeats (LRR)-transmembrane protein kinase, pathogen-related (PR) gene 5 receptor kinase and nucleotide binding site-leucine rich repeats (NBS-LRR) domain containing plant resistance gene (R-gene) were also significantly up-regulated. These results suggest that BD may induce changes to the plant transcriptome in induced systemic resistance (ISR) pathways.

Chickpea-Fusarium oxysporum interaction transcriptome reveals differential modulation of plant defense strategies

Fusarium wilt is one of the major biotic stresses reducing chickpea productivity. The use of wilt-resistant cultivars is the most appropriate means to combat the disease and secure productivity. As a step towards understanding the molecular basis of wilt resistance in chickpea, we investigated the transcriptomes of wilt-susceptible and wilt-resistant cultivars under both Fusarium oxysporum f.sp. ciceri (Foc) challenged and unchallenged conditions. Transcriptome profiling using LongSAGE provided a valuable insight into the molecular interactions between chickpea and Foc, which revealed several known as well as novel genes with differential or unique expression patterns in chickpea contributing to lignification, hormonal homeostasis, plant defense signaling, ROS homeostasis, R-gene mediated defense, etc. Similarly, several Foc genes characteristically required for survival and growth of the pathogen were expressed only in the susceptible cultivar with null expression of most of these genes in the resistant cultivar. This study provides a rich resource for functional characterization of the genes involved in resistance mechanism and their use in breeding for sustainable wilt-resistance. Additionally, it provides pathogen targets facilitating the development of novel control strategies.

Genome reconstruction in Cynara cardunculus taxa gains access to chromosome-scale DNA variation

The genome sequence of globe artichoke (Cynara cardunculus L. var. scolymus, 2n = 2x = 34) is now available for use. A survey of C. cardunculus genetic resources is essential for understanding the evolution of the species, carrying out genetic studies and for application of breeding strategies. We report on the resequencing analyses (~35×) of four globe artichoke genotypes, representative of the core varietal types, as well as a genotype of the related taxa cultivated cardoon. The genomes were reconstructed at a chromosomal scale and structurally/functionally annotated. Gene prediction indicated a similar number of genes, while distinctive variations in miRNAs and resistance gene analogues (RGAs) were detected. Overall, 23,5 M SNP/indel were discovered (range 6,34 M –14,50 M). The impact of some missense SNPs on the biological functions of genes involved in the biosynthesis of phenylpropanoid and sesquiterpene lactone secondary metabolites was predicted. The identified variants contribute to infer on globe artichoke domestication of the different varietal types, and represent key tools for dissecting the path from sequence variation to phenotype. The new genomic sequences are fully searchable through independent Jbrowse interfaces (www.artichokegenome.unito.it), which allow the analysis of collinearity and the discovery of genomic variants, thus representing a one-stop resource for C. cardunculus genomics.

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.

Variability among Cucurbitaceae species (melon, cucumber and watermelon) in a genomic region containing a cluster of NBS-LRR genes

BACKGROUND

Cucurbitaceae species contain a significantly lower number of genes coding for proteins with similarity to plant resistance genes belonging to the NBS-LRR family than other plant species of similar genome size. A large proportion of these genes are organized in clusters that appear to be hotspots of variability. The genomes of the Cucurbitaceae species measured until now are intermediate in size (between 350 and 450 Mb) and they apparently have not undergone any genome duplications beside those at the origin of eudicots. The cluster containing the largest number of NBS-LRR genes has previously been analyzed in melon and related species and showed a high degree of interspecific and intraspecific variability. It was of interest to study whether similar behavior occurred in other cluster of the same family of genes.

RESULTS

The cluster of NBS-LRR genes located in melon chromosome 9 was analyzed and compared with the syntenic regions in other cucurbit genomes. This is the second cluster in number within this species and it contains nine sequences with a NBS-LRR annotation including two genes, Fom1 and Prv, providing resistance against Fusarium and Ppapaya ring-spot virus (PRSV). The variability within the melon species appears to consist essentially of single nucleotide polymorphisms. Clusters of similar genes are present in the syntenic regions of the two species of Cucurbitaceae that were sequenced, cucumber and watermelon. Most of the genes in the syntenic clusters can be aligned between species and a hypothesis of generation of the cluster is proposed. The number of genes in the watermelon cluster is similar to that in melon while a higher number of genes (12) is present in cucumber, a species with a smaller genome than melon. After comparing genome resequencing data of 115 cucumber varieties, deletion of a group of genes is observed in a group of varieties of Indian origin.

CONCLUSIONS

Clusters of genes coding for NBS-LRR proteins in cucurbits appear to have specific variability in different regions of the genome and between different species. This observation is in favour of considering that the adaptation of plant species to changing environments is based upon the variability that may occur at any location in the genome and that has been produced by specific mechanisms of sequence variation acting on plant genomes. This information could be useful both to understand the evolution of species and for plant breeding.

RGAugury: a pipeline for genome-wide prediction of resistance gene analogs (RGAs) in plants

Background

Resistance gene analogs (RGAs), such as NBS-encoding proteins, receptor-like protein kinases (RLKs) and receptor-like proteins (RLPs), are potential R-genes that contain specific conserved domains and motifs. Thus, RGAs can be predicted based on their conserved structural features using bioinformatics tools. Computer programs have been developed for the identification of individual domains and motifs from the protein sequences of RGAs but none offer a systematic assessment of the different types of RGAs. A user-friendly and efficient pipeline is needed for large-scale genome-wide RGA predictions of the growing number of sequenced plant genomes.

Results

An integrative pipeline, named RGAugury, was developed to automate RGA prediction. The pipeline first identifies RGA-related protein domains and motifs, namely nucleotide binding site (NB-ARC), leucine rich repeat (LRR), transmembrane (TM), serine/threonine and tyrosine kinase (STTK), lysin motif (LysM), coiled-coil (CC) and Toll/Interleukin-1 receptor (TIR). RGA candidates are identified and classified into four major families based on the presence of combinations of these RGA domains and motifs: NBS-encoding, TM-CC, and membrane associated RLP and RLK. All time-consuming analyses of the pipeline are paralleled to improve performance. The pipeline was evaluated using the well-annotated Arabidopsis genome. A total of 98.5, 85.2, and 100 % of the reported NBS-encoding genes, membrane associated RLPs and RLKs were validated, respectively. The pipeline was also successfully applied to predict RGAs for 50 sequenced plant genomes. A user-friendly web interface was implemented to ease command line operations, facilitate visualization and simplify result management for multiple datasets.

Conclusions

RGAugury is an efficiently integrative bioinformatics tool for large scale genome-wide identification of RGAs. It is freely available at Bitbucket: https://bitbucket.org/yaanlpc/rgaugury.

Electronic supplementary material

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

De Novo Assembled Wheat Transcriptomes Delineate Differentially Expressed Host Genes in Response to Leaf Rust Infection

Pathogens like Puccinia triticina, the causal organism for leaf rust, extensively damages wheat production. The interaction at molecular level between wheat and the pathogen is complex and less explored. The pathogen induced response was characterized using mock- or pathogen inoculated near-isogenic wheat lines (with or without seedling leaf rust resistance gene Lr28). Four Serial Analysis of Gene Expression libraries were prepared from mock- and pathogen inoculated plants and were subjected to Sequencing by Oligonucleotide Ligation and Detection, which generated a total of 165,767,777 reads, each 35 bases long. The reads were processed and multiple k-mers were attempted for de novo transcript assembly; 22 k-mers showed the best results. Altogether 21,345 contigs were generated and functionally characterized by gene ontology annotation, mining for transcription factors and resistance genes. Expression analysis among the four libraries showed extensive alterations in the transcriptome in response to pathogen infection, reflecting reorganizations in major biological processes and metabolic pathways. Role of auxin in determining pathogenesis in susceptible and resistant lines were imperative. The qPCR expression study of four LRR-RLK (Leucine-rich repeat receptor-like protein kinases) genes showed higher expression at 24 hrs after inoculation with pathogen. In summary, the conceptual model of induced resistance in wheat contributes insights on defense responses and imparts knowledge of Puccinia triticina-induced defense transcripts in wheat plants.

Positive selection in the leucine-rich repeat domain of Gro1 genes in Solanum species

In pathogen resistant plants, solvent-exposed residues in the leucine-rich repeat (LRR) proteins are thought to mediate resistance by recognizing plant pathogen elicitors. In potato, the gene Gro1-4 confers resistance to Globodera rostochiensis. The investigation of variability in different copies of this gene represents a good model for the verification of positive selection mechanisms. Two datasets of Gro1 LRR sequences were constructed, one derived from the Gro1-4 gene, belonging to different cultivated and wild Solanum species, and the other belonging to paralogues of a resistant genotype. Analysis of nonsynonymous to synonymous substitution rates (K(a)/K(s)) highlighted 14 and six amino acids with K(a)/K(s) >1 in orthologue and paralogue datasets, respectively. Selection analysis revealed that the leucine-rich regions accumulate variability in a very specific way, and we found that some combinations of amino acids in these sites might be involved in pathogen recognition. The results confirm previous studies on positive selection in the LRR domain of R protein in Arabidopsis and other model plants and extend these to wild Solanum species. Moreover, positively selected sites in the Gro1 LRR domain show that coevolution mainly occurred in two regions on the internal surface of the three-dimensional horseshoe structure of the domain, albeit with different evolutionary forces between paralogues and orthologues.

Disease Resistance Gene Analogs (RGAs) in Plants

Plants have developed effective mechanisms to recognize and respond to infections caused by pathogens. Plant resistance gene analogs (RGAs), as resistance (R) gene candidates, have conserved domains and motifs that play specific roles in pathogens' resistance. Well-known RGAs are nucleotide binding site leucine rich repeats, receptor like kinases, and receptor like proteins. Others include pentatricopeptide repeats and apoplastic peroxidases. RGAs can be detected using bioinformatics tools based on their conserved structural features. Thousands of RGAs have been identified from sequenced plant genomes. High-density genome-wide RGA genetic maps are useful for designing diagnostic markers and identifying quantitative trait loci (QTL) or markers associated with plant disease resistance. This review focuses on recent advances in structures and mechanisms of RGAs, and their identification from sequenced genomes using bioinformatics tools. Applications in enhancing fine mapping and cloning of plant disease resistance genes are also discussed.

An advanced draft genome assembly of a desi type chickpea (Cicer arietinum L.)

Chickpea (Cicer arietinum L.) is an important pulse legume crop. We previously reported a draft genome assembly of the desi chickpea cultivar ICC 4958. Here we report an advanced version of the ICC 4958 genome assembly (version 2.0) generated using additional sequence data and an improved genetic map. This resulted in 2.7-fold increase in the length of the pseudomolecules and substantial reduction of sequence gaps. The genome assembly covered more than 94% of the estimated gene space and predicted the presence of 30,257 protein-coding genes including 2230 and 133 genes encoding potential transcription factors (TF) and resistance gene homologs, respectively. Gene expression analysis identified several TF and chickpea-specific genes with tissue-specific expression and displayed functional diversification of the paralogous genes. Pairwise comparison of pseudomolecules in the desi (ICC 4958) and the earlier reported kabuli (CDC Frontier) chickpea assemblies showed an extensive local collinearity with incongruity in the placement of large sequence blocks along the linkage groups, apparently due to use of different genetic maps. Single nucleotide polymorphism (SNP)-based mining of intra-specific polymorphism identified more than four thousand SNPs differentiating a desi group and a kabuli group of chickpea genotypes.

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.

Transcriptome analysis of resistant and susceptible alfalfa cultivars infected with root-knot nematode Meloidogyne incognita

Nematodes are one of the major limiting factors in alfalfa production. Root-knot nematodes (RKN, Meloidogyne spp.) are widely distributed and economically important sedentary endoparasites of agricultural crops and they may inflict significant damage to alfalfa fields. As of today, no studies have been published on global gene expression profiling in alfalfa infected with RKN or any other plant parasitic nematode. Very little information is available about molecular mechanisms that contribute to pathogenesis and defense responses in alfalfa against these pests and specifically against RKN. In this work, we performed root transcriptome analysis of resistant (cv. Moapa 69) and susceptible (cv. Lahontan) alfalfa cultivars infected with RKN Meloidogyne incognita, widespread root-knot nematode species and a major pest worldwide. A total of 1,701,622,580 pair-end reads were generated on an Illumina Hi-Seq 2000 platform from the roots of both cultivars and assembled into 45,595 and 47,590 transcripts in cvs Moapa 69 and Lahontan, respectively. Bioinformatic analysis revealed a number of common and unique genes that were differentially expressed in susceptible and resistant lines as a result of nematode infection. Although the susceptible cultivar showed a more pronounced defense response to the infection, feeding sites were successfully established in its roots. Characteristically, basal gene expression levels under normal conditions differed between the two cultivars as well, which may confer advantage to one of the genotypes toward resistance to nematodes. Differentially expressed genes were subsequently assigned to known Gene Ontology categories to predict their functional roles and associated biological processes. Real-time PCR validated expression changes in genes arbitrarily selected for experimental confirmation. Candidate genes that contribute to protection against M. incognita in alfalfa were proposed and alfalfa-nematode interactions with respect to resistance are discussed.

Transcriptome of Erysiphe necator-infected Vitis pseudoreticulata leaves provides insight into grapevine resistance to powdery mildew

Powdery mildew (PM), which is caused by the pathogen Erysiphe necator (Schw.) Burr., is the single most damaging disease of cultivated grapes (Vitis vinifera) worldwide. However, little is known about the transcriptional response of grapes to infection with PM. RNA-seq analysis was used for deep sequencing of the leaf transcriptome to study PM resistance in Chinese wild grapes (V. pseudoreticulata Baihe 35-1) to better understand the interaction between host and pathogen. Greater than 100 million (M) 90-nt cDNA reads were sequenced from a cDNA library derived from PM-infected leaves. Among the sequences obtained, 6541 genes were differentially expressed (DEG) and were annotated with Gene Ontology terms and by pathway enrichment. The significant categories that were identified included the following: defense, salicylic acid (SA) and jasmonic acid (JA) responses; systemic acquired resistance (SAR); hypersensitive response; plant-pathogen interaction; flavonoid biosynthesis; and plant hormone signal transduction. Various putative secretory proteins were identified, indicating potential defense responses to PM infection. In all, 318 putative R-genes and 183 putative secreted proteins were identified, including the defense-related R-genes BAK1, MRH1 and MLO3 and the defense-related secreted proteins GLP and PR5. The expression patterns of 16 genes were further illuminated by RT-qPCR. The present study identified several candidate genes and pathways that may contribute to PM resistance in grapes and illustrated that RNA-seq is a powerful tool for studying gene expression. The RT-qPCR results reveal that effective resistance responses of grapes to PM include enhancement of JA and SAR responses and accumulation of phytoalexins.