EST-based in silico identification and in vitro test of antimicrobial peptides in Brassica napus

Background

Brassica napus is the third leading source of vegetable oil in the world after soybean and oil palm. The accumulation of gene sequences, especially expressed sequence tags (ESTs) from plant cDNA libraries, has provided a rich resource for genes discovery including potential antimicrobial peptides (AMPs). In this study, we used ESTs including those generated from B. napus cDNA libraries of seeds, pathogen-challenged leaves and deposited in the public databases, as a model, to perform in silico identification and consequently in vitro confirmation of putative AMP activities through a highly efficient system of recombinant AMP prokaryotic expression.

Results

In total, 35,788 were generated from cDNA libraries of pathogen-challenged leaves and 187,272 ESTs from seeds of B. napus, and the 644,998 ESTs of B. napus were downloaded from the EST database of PlantGDB. They formed 201,200 unigenes. First, all the known AMPs from the AMP databank (APD2 database) were individually queried against all the unigenes using the BLASTX program. A total of 972 unigenes that matched the 27 known AMP sequences in APD2 database were extracted and annotated using Blast2GO program. Among these unigenes, 237 unigenes from B. napus pathogen-challenged leaves had the highest ratio (1.15 %) in this unigene dataset, which is 13 times that of the unigene datasets of B. napus seeds (0.09 %) and 2.3 times that of the public EST dataset. About 87 % of each EST library was lipid-transfer protein (LTP) (32 % of total unigenes), defensin, histone, endochitinase, and gibberellin-regulated proteins. The most abundant unigenes in the leaf library were endochitinase and defensin, and LTP and histone in the pub EST library. After masking of the repeat sequence, 606 peptides that were orthologous matched to different AMP families were found. The phylogeny and conserved structural motifs of seven AMPs families were also analysed. To investigate the antimicrobial activities of the predicted peptides, 31 potential AMP genes belonging to different AMP families were selected to test their antimicrobial activities after bioinformatics identification. The AMP genes were all optimized according to Escherichia coli codon usage and synthetized through one-step polymerase chain reaction method. The results showed that 28 recombinant AMPs displayed expected antimicrobial activities against E. coli and Micrococcus luteus and Sclerotinia sclerotiorum strains.

Conclusion

The study not only significantly expanded the number of known/predicted peptides, but also contributed to long-term plant genetic improvement for increased resistance to diverse pathogens of B.napus. These results proved that the high-throughput method developed that combined an in silico procedure with a recombinant AMP prokaryotic expression system is considerably efficient for identification of new AMPs from genome or EST sequence databases.

Electronic supplementary material

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

Construction of EST Database for Comparative Gene Studies of Acanthamoeba

The genus Acanthamoeba can cause severe infections such as granulomatous amebic encephalitis and amebic keratitis in humans. However, little genomic information of Acanthamoeba has been reported. Here, we constructed Acanthamoeba expressed sequence tags (EST) database (Acanthamoeba EST DB) derived from our 4 kinds of Acanthamoeba cDNA library. The Acanthamoeba EST DB contains 3,897 EST generated from amebae under various conditions of long term in vitro culture, mouse brain passage, or encystation, and downloaded data of Acanthamoeba from National Center for Biotechnology Information (NCBI) and Taxonomically Broad EST Database (TBestDB). The almost reported cDNA/genomic sequences of Acanthamoeba provide stand alone BLAST system with nucleotide (BLAST NT) and amino acid (BLAST AA) sequence database. In BLAST results, each gene links for the significant information including sequence data, gene orthology annotations, relevant references, and a BlastX result. This is the first attempt for construction of Acanthamoeba database with genes expressed in diverse conditions. These data were integrated into a database (http://www.amoeba.or.kr).     

Psoroptes ovis: identification of vaccine candidates by immunoscreening

Serum from successful vaccine trials against the sheep scab mite, Psoroptes ovis, was used to immunoscreen a cDNA library constructed from mixed-stage and gender P. ovis to identify potential recombinant vaccine candidates. Immunodominant recombinant proteins recognised by IgG in these sera were selected for further analysis. Two candidates were identified in this way; a catchin-like protein (CLP) and a novel mu class glutathione S-transferase (GST). Both candidates were expressed in bacteria as recombinant proteins, the GST as an active enzyme, and combined with four other recombinant allergens in a multi-component recombinant vaccine. Strong serum IgG responses were induced in sheep against each of the components of the recombinant vaccine, however, the protective efficacy of the vaccine could not be determined because of variability in the establishment of a challenge infection.     

Rapid cloning and bioinformatic analysis of spinach Y chromosome-specific EST sequences

The genome of spinach single chromosome complement is about 1000 Mbp, which is the model material to study the molecular mechanisms of plant sex differentiation. The cytological study showed that the biggest spinach chromosome (chromosome 1) was taken as spinach sex chromosome. It had three alleles of sex-related X, X ^m and Y. Many researchers have been trying to clone the sex-determining genes and investigated the molecular mechanism of spinach sex differentiation. However, there are no successful cloned reports about these genes. A new technology combining chromosome microdissection with hybridization-specific amplification (HSA) was adopted. The spinach Y chromosome degenerate oligonucleotide primed-PCR (DOP-PCR) products were hybridized with cDNA of the male spinach flowers in florescence. The female spinach genome was taken as blocker and cDNA library specifically expressed in Y chromosome was constructed. Moreover, expressed sequence tag (EST) sequences in cDNA library were cloned, sequenced and bioinformatics was analysed. There were 63 valid EST sequences obtained in this study. The fragment size was between 53 and 486 bp. BLASTn homologous alignment indicated that 12 EST sequences had homologous sequences of nucleic acids, the rest were new sequences. BLASTx homologous alignment indicated that 16 EST sequences had homologous protein-encoding nucleic acid sequence. The spinach Y chromosome-specific EST sequences laid the foundation for cloning the functional genes, specifically expressed in spinach Y chromosome. Meanwhile, the establishment of the technology system in the research provided a reference for rapid cloning of other biological sex chromosome-specific EST sequences.     

Computational analysis of common bean (Phaseolus vulgaris L., genotype BAT93) lycopene β-cyclase and β-carotene hydroxylase gene's cDNA

The identification of genes and understanding of genes'' expression and regulation in common bean (Phaseolus vulgaris L.) is necessary in order to strategize its improvement using genetic engineering techniques. Generation of expressed sequence tags (ESTs) is useful in rapid isolation, identification and characterization of the genes. To study the gene expression in P. vulgaris pods tissue, ESTs generation work was initiated. Early stage and late stage bean-pod-tissues cDNA libraries were constructed using CloneMiner cDNA library construction kit. In total, 5972 EST clones were isolated using random method of gene isolation. While processing ESTs, we found lycopene β-cyclase (PvLCY-β) and β-carotene hydroxylase (PvCHY-β) gene''s cDNA. In carotenoid biosynthesis pathway, PvLCY-β catalyzes the production of carotene; and PvCHY-β is known to function as a catalyst in the production of lutein and zeaxanthin. To understand more about PvLCY-β and PvCHY-β, both strands of both cDNA clones were sequenced using M13 forward and reverse primers. Nucleotide and deduced protein sequences were analyzed and annotated using online bioinformatics tools. Results showed that PvLCY-β and PvCHY-β cDNAs are 1639 and 1107 bp in length, respectively. Analysis results showed that PvLCY-β and PvCHY-β gene''s cDNA contains an open reading frame (ORF) that encodes for 502 and 305 amino acid residues, respectively. The deduced protein sequence analysis results also showed the presence of conserved domains needed for PvLCY-β and PvCHY-β functions. The phylogenetic analysis of both PvLCY-β and PvCHY-β proteins showed it''s closeness with the LCY-β and CHY-β proteins from Glycine max, respectively. The nucleotide sequence of PvLCY-β and PvCHY-β gene''s cDNA and it''s annotation is reported in this paper.     

A novel Glycine soja tonoplast intrinsic protein gene responds to abiotic stress and depresses salt and dehydration tolerance in transgenic Arabidopsis thaliana

Tonoplast intrinsic protein (TIP) is a subfamily of the aquaporin (AQP), also known as major intrinsic protein (MIP) family, and regulates water movement across vacuolar membranes. Some reports have implied that TIP genes are associated with plant tolerance to some abiotic stresses that cause water loss, such as drought and high salinity. In our previous work, we found that an expressed sequence tag (EST) representing a TIP gene in our Glycine soja EST library was inducible by abiotic stresses. This TIP was subsequently isolated from G. soja with cDNA library screening, EST assembly and PCR, and named as GsTIP2;1. The expression patterns of GsTIP2;1 in G. soja under low temperature, salt and dehydration stress were different in leaves and roots. Though GsTIP2;1 is a stress-induced gene, overexpression of GsTIP2;1 in Arabidopsis thaliana depressed tolerance to salt and dehydration stress, but did not affect seedling growth under cold or favorable conditions. Higher dehydration speed was detected in Arabidopsis plants overexpressing GsTIP2;1, implying GsTIP2;1 might mediate stress sensitivity by enhancing water loss in the plant. Such a result is not identical to previous reports, providing some new information about the relationship between TIP and plant abiotic stress tolerance.     

Identification of a human immunodominant B-cell epitope within the GRA1 antigen of Toxoplasma gondii by phage display of cDNA libraries

Excreted secreted antigens of the protozoan parasite Toxoplasma gondii play a key role in stimulating the host immune system during acute and chronic infection. With the aim of identifying the immunodominant epitopes of T. gondii antigens involved in the human B-cell response against the parasite, we employed a novel immunological approach. A library of cDNA fragments from T. gondii tachyzoites was displayed as fusion proteins to the amino-terminus of lambda bacteriophage capsid protein D. The lambdaD-tachyzoite library was then affinity-selected by using a panel of sera of pregnant women, all infected with the parasite. Some of the clones identified through this procedure matched the sequence of the dense granule GRA1 protein (p24), allowing us to identify its antigenic regions. In particular, the analysis of human antibody response against the recombinant GRA1 antigen fragments revealed the existence of an immunodominant epitope (epi-24 peptide).     

Proteomics-Identified Bvg-Activated Autotransporters Protect against Bordetella pertussis in a Mouse Model

Pertussis is a highly infectious respiratory disease of humans caused by the bacterium Bordetella pertussis. Despite high vaccination coverage, pertussis has re-emerged globally. Causes for the re-emergence of pertussis include limited duration of protection conferred by acellular pertussis vaccines (aP) and pathogen adaptation. Pathogen adaptations involve antigenic divergence with vaccine strains, the emergence of strains which show enhanced in vitro expression of a number of virulence-associated genes and of strains that do not express pertactin, an important aP component. Clearly, the identification of more effective B. pertussis vaccine antigens is of utmost importance. To identify novel antigens, we used proteomics to identify B. pertussis proteins regulated by the master virulence regulatory system BvgAS in vitro. Five candidates proteins were selected and it was confirmed that they were also expressed in the lungs of naïve mice seven days after infection. The five proteins were expressed in recombinant form, adjuvanted with alum and used to immunize mice as stand-alone antigens. Subsequent respiratory challenge showed that immunization with the autotransporters Vag8 and SphB1 significantly reduced bacterial load in the lungs. Whilst these antigens induced strong opsonizing antibody responses, we found that none of the tested alum-adjuvanted vaccines - including a three-component aP - reduced bacterial load in the nasopharynx, suggesting that alternative immunological responses may be required for efficient bacterial clearance from the nasopharynx.     

Cloning of a cDNA fragment encoding part of the protein moiety of the 58-kDa fibrinogen-binding mannoprotein of Candida albicans

Immunoscreening of a Candida albicans expression library with antibodies against the 58 kDa fibrinogen-binding mannoprotein (mp58) of the fungus resulted in the isolation of clones encoding the protein moiety of this molecule. Sequence of the 0.9 kb cDNA of one of the clones selected for further analysis, revealed an open reading frame coding for 292 amino acids, which displays sequence similarity to proteins belonging to a family of immunodominant antigens of Aspergillus spp. The gene corresponding to this cDNA was named FBP1 (fibrinogen-binding protein). These results represent the first report on the identification of C. albicans genes encoding surface receptors for host proteins.     

Expressed sequence tags from the zhikong scallop (Chlamys farreri): Discovery and annotation of host-defense genes

A high-quality cDNA library was constructed from whole body tissues of the zhikong scallop, Chlamys farreri, challenged by Listonella anguillarum. A total of 5720 clones were sequenced, yielding 5123 expressed sequence tags (ESTs). Among the 3326 unique genes identified, 2289 (69%) genes had no significant (E-value < 1e?5) matches to known sequences in public databases and 194 (6%) matched proteins of unknown functions. The remaining 843 (25%) genes that exhibited homology with genes of known functions, showed broad involvement in metabolic processes (31%), cell structure and motility (20%), gene and protein expression (12%), cell signaling and cell communication (8%), cell division (4%), and notably, 25% of those genes were related to immune function. They included stress response genes, complement-like genes, proteinase and proteinase inhibitors, immune recognition receptors and immune effectors. The EST collection obtained in this study provides a useful resource for gene discovery and especially for the identification of host-defense genes and systems in scallops and other molluscs.     

Genes transcribed in the salivary glands of female Rhipicephalus appendiculatus ticks infected with Theileria parva

We describe the generation of an auto-annotated index of genes that are expressed in the salivary glands of four-day fed female adult Rhipicephalus appendiculatus ticks. A total of 9162 EST sequences were derived from an uninfected tick cDNA library and 9844 ESTs were from a cDNA library from ticks infected with Theileria parva, which develop in type III salivary gland acini. There were no major differences between abundantly expressed ESTs from the two cDNA libraries, although there was evidence for an up-regulation in the expression of some glycine-rich proteins in infected salivary glands. Gene ontology terms were also assigned to sequences in the index and those with potential enzyme function were linked to the Kyoto encyclopedia of genes and genomes database, allowing reconstruction of metabolic pathways. Several genes code for previously characterized tick proteins such as receptors for myokinin or ecdysteroid and an immunosuppressive protein. cDNAs coding for homologs of heme-lipoproteins which are major components of tick hemolymph were identified by searching the database with published N-terminal peptide sequence data derived from biochemically purified Boophilus microplus proteins. The EST data will be a useful resource for construction of microarrays to probe vector biology, vector-host and vector-pathogen interactions and to underpin gene identification via proteomics approaches.     

Expression and interaction of the CBLs and CIPKs from immature seeds of kidney bean (Phaseolus vulgaris L.)

Protein phosphorylation plays a key regulatory role in a variety of cellular processes. To better understand the function of protein phosphorylation in seed maturation, a PCR-based cloning method was employed and five cDNA clones (pvcipk1-5) for protein kinases were isolated from a cDNA library prepared from immature seeds of kidney bean (Phaseolus vulgaris L.). The deduced amino acid sequences showed that the five protein kinases (PvCIPK1-5) are members of the sucrose non-fermenting 1-related protein kinase type 3 (SnRK3) family, which interacts with calcineurin B-like proteins (CBLs). Two cDNA clones (pvcbl1 and 2) for CBLs were further isolated from the cDNA library. The predicted primary sequences of the proteins (PvCBL1 and 2) displayed significant identity (more than 90%) with those of other plant CBLs. Semi-quantitative RT-PCR analysis showed that the isolated genes, except pvcbl1, are expressed in leaves and early maturing seeds, whereas pvcbl1 is constitutively expressed during seed development. Yeast two-hybrid assay indicated that among the five PvCIPKs, only PvCIPK1 interacts with both PvCBL1 and PvCBL2. These results suggest that calcium-dependent protein phosphorylation-signaling via CBL-CIPK complexes occurs during seed development.