The prediction of a pathogenesis-related secretome of Puccinia helianthi through high-throughput transcriptome analysis

Codon usage bias of secretory protein in Fusarium oxysporum f. sp. cubense tropical race 4

Banana Fusarium oxysporum f. sp. cubense tropical race 4 (Foc-TR4) is a highly destructive pathogen that infects nearly all major banana cultivars and has a tendency to spread further. Secreted proteins play a crucial role in the process of Fusarium wilt infection in bananas. In this study, we analyzed the codon usage bias (CUB) of the Foc-TR4 classical secretory protein genome for the first time and observed a strong bias toward codons ending with C. We found that 572 out of the 14,543 amino acid sequences in the Foc-TR4 genome exhibited characteristics of classical secretory proteins. The CUB was largely influenced by selection optimization pressure, as indicated by the ENC value and neutral plot analysis. Among the identified codons, such as UCC and CCC, 11 were found to be optimal for Foc-TR4 gene expression. Codons with higher GC content and a C base in the third position showed greater selectivity. The CUB in the secretory proteins encoded by Foc-TR4 provides insights into their evolutionary patterns, contributing to the development and screening of novel and effective antifungal drugs.

In silico secretome analyses of the polyphagous root-knot nematode Meloidogyne javanica: a resource for studying M. javanica secreted proteins

BACKGROUND

Plant-parasitic nematodes (PPNs) that cause most damage include root-knot nematodes (RKNs) which are a major impediment to crop production. Root-knot nematodes, like other parasites, secrete proteins which are required for parasite proliferation and survival within the host during the infection process.

RESULTS

Here, we used various computational tools to predict and identify classically and non-classically secreted proteins encoded in the Meloidogyne javanica genome. Furthermore, functional annotation analysis was performed using various integrated bioinformatic tools to determine the biological significance of the predicted secretome. In total, 7,458 proteins were identified as secreted ones. A large percentage of this secretome is comprised of small proteins of ≤ 300 aa sequence length. Functional analyses showed that M. javanica secretome comprises cell wall degrading enzymes for facilitating nematode invasion, and migration by disintegrating the complex plant cell wall components. In addition, peptidases and peptidase inhibitors are an important category of M. javanica secretome involved in compatible host-nematode interactions.

CONCLUSION

This study identifies the putative secretome encoded in the M. javanica genome. Future experimental validation analyses can greatly benefit from this global analysis of M. javanica secretome. Equally, our analyses will advance knowledge of the interaction between plants and nematodes.

Genome-Wide Characterization of Effector Protein-Encoding Genes in Sclerospora graminicola and Its Validation in Response to Pearl Millet Downy Mildew Disease Stress

Pearl millet [Pennisetum glaucum (L.) R. Br.] is the essential food crop for over ninety million people living in drier parts of India and South Africa. Pearl millet crop production is harshly hindered by numerous biotic stresses. Sclerospora graminicola causes downy mildew disease in pearl millet. Effectors are the proteins secreted by several fungi and bacteria that manipulate the host cell structure and function. This current study aims to identify genes encoding effector proteins from the S. graminicola genome and validate them through molecular techniques. In silico analyses were employed for candidate effector prediction. A total of 845 secretory transmembrane proteins were predicted, out of which 35 proteins carrying LxLFLAK (Leucine–any amino acid–Phenylalanine–Leucine–Alanine–Lysine) motif were crinkler, 52 RxLR (Arginine, any amino acid, Leucine, Arginine), and 17 RxLR-dEER putative effector proteins. Gene validation analysis of 17 RxLR-dEER effector protein-producing genes was carried out, of which 5genes were amplified on the gel. These novel gene sequences were submitted to NCBI. This study is the first report on the identification and characterization of effector genes in Sclerospora graminicola. This dataset will aid in the integration of effector classes that act independently, paving the way to investigate how pearl millet responds to effector protein interactions. These results will assist in identifying functional effector proteins involving the omic approach using newer bioinformatics tools to protect pearl millet plants against downy mildew stress. Considered together, the identified effector protein-encoding functional genes can be utilized in screening oomycetes downy mildew diseases in other crops across the globe.

Genome-Wide Identification, Evolutionary and Functional Analyses of WRKY Family Members in Ginkgo biloba

WRKY transcription factors (TFs) are one of the largest families in plants which play essential roles in plant growth and stress response. Ginkgo biloba is a living fossil that has remained essentially unchanged for more than 200 million years, and now has become widespread worldwide due to the medicinal active ingredients in its leaves. Here, 37 WRKY genes were identified, which were distributed randomly in nine chromosomes of G. biloba. Results of the phylogenetic analysis indicated that the GbWRKY could be divided into three groups. Furthermore, the expression patterns of GbWRKY genes were analyzed. Gene expression profiling and qRT-PCR revealed that different members of GbWRKY have different spatiotemporal expression patterns in different abiotic stresses. Most of the GbWRKY genes can respond to UV-B radiation, drought, high temperature and salt treatment. Meanwhile, all GbWRKY members performed phylogenetic tree analyses with the WRKY proteins of other species which were known to be associated with abiotic stress. The result suggested that GbWRKY may play a crucial role in regulating multiple stress tolerances. Additionally, GbWRKY13 and GbWRKY37 were all located in the nucleus, while GbWRKY15 was located in the nucleus and cytomembrane.

Cloning and sequence analysis of a serine protease gene from Rhizoctonia solani Kühn AG5

The present study aimed to identify the subtilisin-like proteases (SLPs) of Rhizoctonia solani Kühn potentially involved in the virulence of this phytopathogenic fungus, which has 14 anastomosis groups (AGs) responsible for many crop diseases. Through mycelial microscope observation and strain identification of pathogenic fungus MS-3, it was determined to be R. solani AG-5. Both 5' and 3' rapid amplification of cDNA ends were used to clone the serine protease gene RsSLP from R. solani AG-5. The full-length obtained for RsSLP was 1714 bp with an open reading frame of 1587 bp, encoding a protein of 528 amino acids with a molecular mass of 55.8 kDa. This protein contained a predicted signal peptide for secretion but lacked a transmembrane domain or membrane anchor site. Bioinformatics analysis identified this protein as a serine protease with the Peptidase_S8 and Inhibitor_I9 characteristic domains of SLPs. Phylogenetic analysis suggested that frequent gene duplications of the SLPs occurred in R. solani (RsSLP), and RsSLP shares characteristic sequence features with virulence factors of other phytopathogenic fungi. Because the secretory serine protease RsSLP from R. solani AG5 is similar to the virulence factors of other phytopathogenic fungi, its identification will be helpful in studies considering the roles of these proteases in pathogen virulence.

Genome-wide identification of the pectin methylesterase inhibitor genes in Brassica napus and expression analysis of selected members

Pectin methylesterase inhibitors (PMEIs) modulate the status of pectin methylesterification by inhibiting the activity of pectin methylesterase (PME). Recent advances indicate PMEIs play an important role in regulating plant cell wall properties and defense responses. In this study, a genome-wide analysis of PMEI gene family in Brassica napus (B. napus) was conducted and the expression patterns of PMEI genes in response to Sclerotinia sclerotiorum (S. sclerotiorum) was investigated. A total of 190 PMEI proteins were identified from the genome of B. napus. Chromosomal location, gene structure and properties of the PMEI family were analyzed, and these features were compared with Arabidopsis thaliana (A. thaliana). A total of 123 syntenic ortholog pairs were detected from BnPMEI family by synteny analysis. Results showed the expansion of BnPMEI genes was likely predominately from whole-genome duplication (WGD) or segmental duplications. Multiple cis-elements related to plant growth and development, environmental stress responses, hormone responses were detected in the promoters of BnPMEI genes, implying they were regulated by both internal and external factors. Furthermore, expression analysis of transcriptome data combined with quantitative RT-PCR (qRT-PCR) validation identified several candidates that were strongly responsive to S. sclerotiorum infection. These BnPMEI genes are candidates for manipulation to breed novel and improved genotypes that are more resistant to sclerotinia stem rot (SSR). Extensive interactions were detected among 30 BnPMEI proteins, forming complex protein-protein interaction networks. Besides, 48 BnPMEIs showed interactions with other proteins including a range of cell wall structure-related enzymes. This study provides new insights into the evolution and function of PMEIs in B. napus and lays a foundation for breeding novel genotypes for crop improvement.

Genome-wide identification, evolutionary and functional analyses of KFB family members in potato

BACKGROUND

Kelch repeat F-box (KFB) proteins play vital roles in the regulation of multitudinous biochemical and physiological processes in plants, including growth and development, stress response and secondary metabolism. Multiple KFBs have been characterized in various plant species, but the family members and functions have not been systematically identified and analyzed in potato.

RESULTS

Genome and transcriptome analyses of StKFB gene family were conducted to dissect the structure, evolution and function of the StKFBs in Solanum tuberosum L. Totally, 44 StKFB members were identified and were classified into 5 groups. The chromosomal localization analysis showed that the 44 StKFB genes were located on 12 chromosomes of potato. Among these genes, two pairs of genes (StKFB15/16 and StKFB40/41) were predicted to be tandemly duplicated genes, and one pair of genes (StKFB15/29) was segmentally duplicated genes. The syntenic analysis showed that the KFBs in potato were closely related to the KFBs in tomato and pepper. Expression profiles of the StKFBs in 13 different tissues and in potato plants with different treatments uncovered distinct spatial expression patterns of these genes and their potential roles in response to various stresses, respectively. Multiple StKFB genes were differentially expressed in yellow- (cultivar 'Jin-16'), red- (cultivar 'Red rose-2') and purple-fleshed (cultivar 'Xisen-8') potato tubers, suggesting that they may play important roles in the regulation of anthocyanin biosynthesis in potato.

CONCLUSIONS

This study reports the structure, evolution and expression characteristics of the KFB family in potato. These findings pave the way for further investigation of functional mechanisms of StKFBs, and also provide candidate genes for potato genetic improvement.

Cross-Kingdom Gene Coexpression Analysis Using a Stemphylium botryosum-Lens ervoides System Revealed Plasticity of Intercommunication Between the Pathogen Secretome and the Host Immune Systems

Necrotrophic pathogens are responsible for significant declines in crop yield and quality worldwide. During the infection process, a pathogen releases a series of secretory proteins to counteract the plant immune system, and this interaction of pathogen and host molecules determines whether the pathogen will successfully invade the host plant tissues. In this study, we adopted co-transcriptomic approaches to analyze the Lens ervoides-Stemphylium botryosum system, with a focus on 1,216 fungal genes coding for secretory proteins and 8,810 disease-responsive genes of the host 48, 96, and 144 h postinoculation, captured in two F₉ recombinant inbred lines (RILs) displaying contrasting disease responses. By constructing in planta gene coexpression networks (GCNs) for S. botryosum, we found that the pathogen tended to co-upregulate genes regulating cell wall degradation enzymes, effectors, oxidoreductases, and peptidases to a much higher degree in the susceptible host LR-66-577 than in the resistant RIL LR-66-637, indicating that the promotion of these digestive enzymes and toxins increased S. botryosum virulence. Construction of cross-kingdom GCNs between pathogen and plant for the two RILs revealed that the co-upregulation of these fungal digestive enzymes and toxins simultaneously promoted a series of defense responses such as redox change, expression of membrane-related genes and serine/threonine kinase, and stress and disease responses in the susceptible RIL which was not observed in the resistant RIL, indicating that these activities exacerbated susceptibility to S. botryosum.[Formula: see text] Copyright © 2021 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Comparative transcriptomic analysis of races 1, 2, 5 and 6 of Fusarium oxysporum f.sp. pisi in a susceptible pea host identifies differential pathogenicity profiles

BACKGROUND

The fungal pathogen Fusarium oxysporum f.sp. pisi (Fop) causes Fusarium wilt in peas. There are four races globally: 1, 2, 5 and 6 and all of these races are present in Australia. Molecular infection mechanisms have been studied in a few other F. oxysporum formae speciales; however, there has been no transcriptomic Fop-pea pathosystem study.

RESULTS

A transcriptomic study was carried out to understand the molecular pathogenicity differences between the races. Transcriptome analysis at 20 days post-inoculation revealed differences in the differentially expressed genes (DEGs) in the Fop races potentially involved in fungal pathogenicity variations. Most of the DEGs in all the races were engaged in transportation, metabolism, oxidation-reduction, translation, biosynthetic processes, signal transduction, proteolysis, among others. Race 5 expressed the most virulence-associated genes. Most genes encoding for plant cell wall degrading enzymes, CAZymes and effector-like proteins were expressed in race 2. Race 6 expressed the least number of genes at this time point.

CONCLUSION

Fop races deploy various factors and complex strategies to mitigate host defences to facilitate colonisation. This investigation provides an overview of the putative pathogenicity genes in different Fop races during the necrotrophic stage of infection. These genes need to be functionally characterised to confirm their pathogenicity/virulence roles and the race-specific genes can be further explored for molecular characterisation.

Genomic Characterization Provides an Insight into the Pathogenicity of the Poplar Canker Bacterium Lonsdalea populi

An emerging poplar canker caused by the gram-negative bacterium, Lonsdalea populi, has led to high mortality of hybrid poplars Populus × euramericana in China and Europe. The molecular bases of pathogenicity and bark adaptation of L. populi have become a focus of recent research. This study revealed the whole genome sequence and identified putative virulence factors of L. populi. A high-quality L. populi genome sequence was assembled de novo, with a genome size of 3,859,707 bp, containing approximately 3434 genes and 107 RNAs (75 tRNA, 22 rRNA, and 10 ncRNA). The L. populi genome contained 380 virulence-associated genes, mainly encoding for adhesion, extracellular enzymes, secretory systems, and two-component transduction systems. The genome had 110 carbohydrate-active enzyme (CAZy)-coding genes and putative secreted proteins. The antibiotic-resistance database annotation listed that L. populi was resistant to penicillin, fluoroquinolone, and kasugamycin. Analysis of comparative genomics found that L. populi exhibited the highest homology with the L. britannica genome and L. populi encompassed 1905 specific genes, 1769 dispensable genes, and 1381 conserved genes, suggesting high evolutionary diversity and genomic plasticity. Moreover, the pan genome analysis revealed that the N-5-1 genome is an open genome. These findings provide important resources for understanding the molecular basis of the pathogenicity and biology of L. populi and the poplar-bacterium interaction.

Smut fungi as a stratagem to characterize rust effectors: opportunities and challenges

The rust pathogens are one of the most complex fungi in the Basidiomycetes. The development of genomic resources for rust and other plant pathogens has opened the opportunities for functional genomics of fungal genes. Despite significant progress in the field of fungal genomics, functional characterization of the genome components has lacked, especially for the rust pathogens. Their obligate nature and lack of standard stable transformation protocol are the primary reasons for rusts to be one of the least explored genera despite its significance. In the recently sequenced rust genomes, a vast catalogue of predicted effectors and pathogenicity genes have been reported. However, most of these candidate genes remained unexplored due to the lack of suitable characterization methods. The heterologous expression of putative effectors in Nicotiana benthamiana and Arabidopsis thaliana has proved to be a rapid screening method for identifying the role of these effectors in virulence. However, no fungal system has been used for the functional validation of these candidate genes. The smuts, from the evolutionary point of view, are closely related to the rust pathogens. Moreover, they have been widely studied and hence could be a suitable model system for expressing rust fungal genes heterologously. The genetic manipulation methods for smuts are also well standardized. Complementation assays can be used for functional validation of the homologous genes present in rust and smut fungal pathogens, while the species-specific proteins can be expressed in the mutant strains of smut pathogens having reduced or no virulence for virulence analysis. We propose that smuts, especially Ustilago maydis, may prove to be a good model system to characterize rust effector proteins in the absence of methods to manipulate the rust genomes directly.

Long-Read-Based de novo Genome Assembly and Comparative Genomics of the Wheat Leaf Rust Pathogen Puccinia triticina Identifies Candidates for Three Avirulence Genes

Leaf rust, caused by Puccinia triticina (Pt), is one of the most devastating diseases of wheat, affecting production in nearly all wheat-growing regions worldwide. Despite its economic importance, genomic resources for Pt are very limited. In the present study, we have used long-read sequencing (LRS) and the pipeline of FALCON and FALCON-Unzip (v4.1.0) to carry out the first LRS-based de novo genome assembly for Pt. Using 22.4-Gb data with an average read length of 11.6 kb and average coverage of 150-fold, we generated a genome assembly for Pt104 [strain 104-2,3,(6),(7),11; isolate S423], considered to be the founding isolate of a clonal lineage of Pt in Australia. The Pt104 genome contains 162 contigs with a total length of 140.5 Mb and N50 of 2 Mb, with the associated haplotigs providing haplotype information for 91% of the genome. This represents the best quality of Pt genome assembly to date, which reduces the contig number by 91-fold and improves the N50 by 4-fold as compared to the previous Pt race1 assembly. An annotation pipeline that combined multiple lines of evidence including the transcriptome assemblies derived from RNA-Seq, previously identified expressed sequence tags and Pt race 1 protein sequences predicted 29,043 genes for Pt104 genome. Based on the presence of a signal peptide, no transmembrane segment, and no target location to mitochondria, 2,178 genes were identified as secreted proteins (SPs). Whole-genome sequencing (Illumina paired-end) was performed for Pt104 and six additional strains with differential virulence profile on the wheat leaf rust resistance genes Lr26, Lr2a, and Lr3ka. To identify candidates for the corresponding avirulence genes AvrLr26, AvrLr2a, and AvrLr3ka, genetic variation within each strain was first identified by mapping to the Pt104 genome. Variants within predicted SP genes between the strains were then correlated to the virulence profiles, identifying 38, 31, and 37 candidates for AvrLr26, AvrLr2a, and AvrLr3ka, respectively. The identification of these candidate genes lays a good foundation for future studies on isolating these avirulence genes, investigating the molecular mechanisms underlying host-pathogen interactions, and the development of new diagnostic tools for pathogen monitoring.

Genome and secretome analysis of jute endophyte Grammothele lineata strain SDL-CO-2015-1: Insights into its lignocellulolytic structure and secondary metabolite profile

Grammothele lineata strain SDL-CO-2015-1, jute (Corchorus olitorius) endophyte has been reported to produce anti-cancer drug paclitaxel in culture condition. Here we investigated the genome using different bioinformatic tools to find its association with the production of commercially important compounds including taxol. Carbohydrate-active enzymes, proteases, and secretory proteins were annotated revealing a complex endophytic relationship with its plant host. The presences of a diverse range of CAZymes including numerous lignocellulolytic enzymes support its potentiality in biomass degradation. Genome annotation led to the identification of 28 clusters for secondary metabolite biosynthesis. Several biosynthesis gene clusters were identified for terpene biosynthesis from antiSMASH analysis but none could be specifically pinned to taxol synthesis. This study will direct us to understand the genomic organization of endophytic basidiomycetes with a potential for producing numerous commercially important enzymes and secondary metabolites taking G. lineata as a model.

Whole-genome and time-course dual RNA-Seq analyses reveal chronic pathogenicity-related gene dynamics in the ginseng rusty root rot pathogen Ilyonectria robusta

Ilyonectria robusta causes rusty root rot, the most devastating chronic disease of ginseng. Here, we for the first time report the high-quality genome of the I. robusta strain CD-56. Time-course (36 h, 72 h, and 144 h) dual RNA-Seq analysis of the infection process was performed, and many genes, including candidate effectors, were found to be associated with the progression and success of infection. The gene expression profile of CD-56 showed a trend of initial inhibition and then gradually returned to a profile similar to that of the control. Analyses of the gene expression patterns and functions of pathogenicity-related genes, especially candidate effector genes, indicated that the stress response changed to an adaptive response during the infection process. For ginseng, gene expression patterns were highly related to physiological conditions. Specifically, the results showed that ginseng defenses were activated by CD-56 infection and persisted for at least 144 h thereafter but that the mechanisms invoked were not effective in preventing CD-56 growth. Moreover, CD-56 did not appear to fully suppress plant defenses, even in late stages after infection. Our results provide new insight into the chronic pathogenesis of CD-56 and the comprehensive and complex inducible defense responses of ginseng root to I. robusta infection.

Genome based screening of epitope ensemble vaccine candidates against dreadful visceral leishmaniasis using immunoinformatics approach

Visceral leishmaniasis (VL) is a dreadful protozoan disease caused by Leishmania donovani that severely affects huge populations in tropical and sub-tropical regions. The present study reports an unbiased genome based screening of 4 potent vaccine antigens against 8023 L. donovani proteins by following the criteria of presence of signal peptides, GPI-anchors and ≤1 transmembrane helix using advanced bioinformatics tools viz. SignalP4.0, PredGPI and TMHMM2.0, respectively. They are designated as genome based predicted signal peptide antigens (GBPSPA). The antigenicity/immunogenicity of chosen vaccine antigens (GBPSPA) with 4 randomly selected known leishmanial antigens (RSKLA) was compared by simulation study employing C-ImmSim software for human immune responses. This revealed better immunological responses. These antigens were further evaluated for the presence of B- and T-cell epitopes using immune epitope database (IEDB) based recommended consensus method of MHC class I and II tools. It was found to forecast CD4⁺ and CD8⁺ T-cell responses in genetically diverse human population worldwide as well as different endemic regions through IEDB based predicted population coverage (PPC) analysis tool. The worldwide percent PPC value of combined (HLA class I and II) epitope ensemble forecast was found to be 99.98, 99.96 and 50.04, respectively for GBPSPA, RSKLA and experimentally known epitopes (EKE) of L. donovani. Therefore, these potential antigens/epitope ensembles could favor the design of prospective and novel vaccine constructs like self-assembled epitopes as nano vaccine formulations against VL. Overall, the present study will serve as a model framework that might improve the effectiveness of designed vaccine against L. donovani and other related pathogens.

In silico prediction of the secretome from the invasive neurotoxic marine dinoflagellate Alexandrium catenella

Alexandrium catenella, a marine dinoflagellate responsible for harmful algal blooms (HABs), proliferates with greater frequency, distribution and intensity, in disturbed marine coastal ecosystems. The proteins secreted into seawater may play a crucial role in maintaining this dinoflagellate in these ecosystems, but this possibility has never been investigated before. In this study, the A. catenella secretome was predicted from its transcriptome by combining several bioinformatics tools. Our results predict a secretome of 2 779 proteins, among which 79% contain less than 500 amino acids, suggesting that most secreted proteins are short in length. The predicted secretome includes 963 proteins (35%) with Pfam domains: 773 proteins with one Pfam domain and 190 proteins with two or more Pfam domains. Their functional annotation showed that they are mainly involved in (i) proteolysis, (ii) stress responses and (iii) primary metabolism. In addition, 47% of the secreted proteins appear to be enzymes, primarily peptidases, known to be biologically active in the extracellular medium during stress responses. Finally, this study provides a wealth of candidates of proteins secreted by A. catenella, which may interact with the marine environment and help this dinoflagellate develop in various environmental conditions.

Comparative Transcriptome Profiling of Gaeumannomyces graminis var. tritici in Wheat Roots in the Absence and Presence of Biocontrol Bacillus velezensis CC09

This study aimed to explore potential biocontrol mechanisms involved in the interference of antagonistic bacteria with fungal pathogenicity in planta. To do this, we conducted a comparative transcriptomic analysis of the “take-all” pathogenic fungus Gaeumannomyces graminis var. tritici (Ggt) by examining Ggt-infected wheat roots in the presence or absence of the biocontrol agent Bacillus velezensis CC09 (Bv) compared with Ggt grown on potato dextrose agar (PDA) plates. A total of 4,134 differentially expressed genes (DEGs) were identified in Ggt-infected wheat roots, while 2,011 DEGs were detected in Bv+Ggt-infected roots, relative to the Ggt grown on PDA plates. Moreover, 31 DEGs were identified between wheat roots, respectively infected with Ggt and Bv+Ggt, consisting of 29 downregulated genes coding for potential Ggt pathogenicity factors – e.g., para-nitrobenzyl esterase, cutinase 1 and catalase-3, and two upregulated genes coding for tyrosinase and a hypothetical protein in the Bv+Ggt-infected roots when compared with the Ggt-infected roots. In particular, the expression of one gene, encoding the ABA3 involved in the production of Ggt’s hormone abscisic acid, was 4.11-fold lower in Ggt-infected roots with Bv than without Bv. This is the first experimental study to analyze the activity of Ggt transcriptomes in wheat roots exposed or not to a biocontrol bacterium. Our results therefore suggest the presence of Bv directly and/or indirectly impairs the pathogenicity of Ggt in wheat roots through complex regulatory mechanisms, such as hyphopodia formation, cell wall hydrolase, and expression of a papain inhibitor, among others, all which merit further investigation.

Exploitation of reverse vaccinology and immunoinformatics as promising platform for genome-wide screening of new effective vaccine candidates against Plasmodium falciparum

Background

In the current scenario, designing of world-wide effective malaria vaccine against Plasmodium falciparum remain challenging despite the significant progress has been made in last few decades. Conventional vaccinology (isolate, inactivate and inject) approaches are time consuming, laborious and expensive; therefore, the use of computational vaccinology tools are imperative, which can facilitate the design of new and promising vaccine candidates.

Results

In current investigation, initially 5548 proteins of P. falciparum genome were carefully chosen for the incidence of signal peptide/ anchor using SignalP4.0 tool that resulted into 640 surface linked proteins (SLP). Out of these SLP, only 17 were predicted to contain GPI-anchors using PredGPI tool in which further 5 proteins were considered as malarial antigenic adhesins by MAAP and VaxiJen programs, respectively. In the subsequent step, T cell epitopes of 5 genome derived predicted antigenic adhesins (GDPAA) and 5 randomly selected known malarial adhesins (RSKMA) were analysed employing MHC class I and II tools of IEDB analysis resource. Finally, VaxiJen scored T cell epitopes from each antigen were considered for prediction of population coverage (PPC) analysis in the world-wide population including malaria endemic regions. The validation of the present in silico strategy was carried out by comparing the PPC of combined (MHC class I and II) predicted epitope ensemble among GDPAA (99.97%), RSKMA (99.90%) and experimentally known epitopes (EKE) of P. falciparum (97.72%) pertaining to world-wide human population.

Conclusions

The present study systematically screened 5 potential protective antigens from P. falciparum genome using bioinformatics tools. Interestingly, these GDPAA, RSKMA and EKE of P. falciparum epitope ensembles forecasted to contain highly promiscuous T cell epitopes, which are potentially effective for most of the world-wide human population with malaria endemic regions. Therefore, these epitope ensembles could be considered in near future for novel and significantly effective vaccine candidate against malaria.

Electronic supplementary material

The online version of this article (10.1186/s12859-018-2482-x) contains supplementary material, which is available to authorized users.

Growth medium and incubation temperature alter the Pseudogymnoascus destructans transcriptome: implications in identifying virulence factors

Pseudogymnoascus destructans is the causal agent of bat white-nose syndrome (WNS), which is devastating some North American bat populations. Previous transcriptome studies provided insight regarding the molecular mechanisms involved in WNS; however, it is unclear how different environmental parameters could influence pathogenicity. This information could be useful in developing management strategies to mitigate the negative impacts of P. destructans on bats. We cultured three P. destructans isolates from Atlantic Canada on two growth media (potato dextrose agar and Sabouraud dextrose agar) that differ in their nitrogen source, and at two separate incubation temperatures (4 C and 15 C) that approximate the temperature range of bat hibernacula during the winter and a temperature within its optimal mycelial growth range. We conducted RNA sequencing to determine transcript levels in each sample and performed differential gene expression (DGE) analyses to test the influence of growth medium and incubation temperature on gene expression. We also compared our in vitro results with previous RNA-sequencing data sets generated from P. destructans growing on the wings of a susceptible host, Myotis lucifugus. Our findings point to a critical role for substrate and incubation temperature in influencing the P. destructans transcriptome. DGE analyses suggested that growth medium plays a larger role than temperature in determining P. destructans gene expression and that although the psychrophilic fungus responds to different nitrogen sources, it may have evolved for continued growth at a broad range of low temperatures. Further, our data suggest that down-regulation of the RNA-interference pathway and increased fatty acid metabolism are involved in the P. destructans-bat interaction. Finally, we speculate that to reduce the activation of host defense responses, P. destructans minimizes changes in the expression of genes encoding secreted proteins during bat colonization.

Gene Cloning, Expression, and Antifungal Activities of Permatin from Naked Oat (Avena nuda)

Thaumatin-like proteins (TLPs) are the products of a large, highly complex gene family involved in host defense. TLPs also belong to the pathogenesis-related family 5 (PR-5) of plant defense proteins. Most TLPs exhibit potential antifungal activities, and their accumulation in the plant is related to many physiological processes. In this study, a gene encoding TLP named permatin with an open reading frame of 678 bp encoding a protein of 225 amino acids with a calculated molecular mass of 23.5 kDa was cloned from naked oat leaves. Phylogenetic analysis revealed that permatin shares high homology with a number of other TLPs among diverse taxa. Model of structure by homology modeling showed that permatin consists of an acidic cleft region consistent with most TLPs. Recombinant NusA-permatin was overexpressed in Escherichia coli strain BL21 and purified by Heparin column combined with Sephacryl S-200 column. The protein exhibited antifungal activity to Fusarium oxysporum (half maximal inhibitory concentration, IC₅₀ = 21.42 μM). Morphological observation showed that NusA-permatin can induce mycelium deformation of F. oxysporum, the cell membrane is blurred, and the diaphragm is not obvious. NusA-permatin also causes membrane permeabilization and reactive oxygen species accumulation in the mycelium of F. oxysporum. Permatin may play an important role in the disease resistance responses of plants against pathogen attacks through its antifungal activity.

Cloning and Characterization of Two Iridoid Synthase Homologs from Swertia Mussotii

Swertia mussotii is an important medicinal plant found on the Qinghai Tibetan Plateau that has great economic and medicinal value. This plant has enjoyed a long history of use as a curative for hepatitis. The biological activity of secoiridoids, including gentiopicroside and swertiamarin, has been mainly tested for its anti-hepatitis effects. Here, we identify two candidate genes (SmIS1 and SmIS2) that are homologues of iridoid synthase and that are components of the secoiridoid pathway in S. mussotii. Using sequencing and phylogenetic analyses, we confirm that SmIS1 and SmIS2 contain six conserved short-chain dehydrogenases/reductase (SDR) motifs and thus belong to the P5βRs group. The two purified Escherichia coli-expressed proteins reduced 8-oxogeranial to both nepetalactol and iridodials. A comparison of the kinetic parameters of SmIS1 and SmIS2 recombinant proteins revealed that SmIS2 has a lower affinity than SmIS1 for 8-oxogeranial. Transcript levels of the two genes were analysed in three different tissues of S. mussotii using semi-quantitative RT-PCR and RT-qPCR. SmIS1 and SmIS2 expression levels were more abundant in leaves and stems. This investigation adds to our knowledge of P5βRs genes in the secoiridoid synthesis pathway and provides candidate genes for genetically improving S. mussotii by enhancing secondary metabolite production.