Characterization of the subtilase gene family in tomato (Lycopersicon esculentum Mill.)

The gene family of subtilisin-like serine proteases (subtilases, SBTs) in tomato (Lycopersicon esculentum Mill.) comprises at least 15 members, 12 of which have been characterized in this study. Sequence comparison revealed that tomato subtilases fall into 5 distinct subfamilies. Single genes were shown to exist for LeSBT1, LeSBT2 and tmp, while 5 and 6 genes were found in the LeSBT3/4 and P69 subfamilies, respectively. With the exception of tmp, tomato subtilase genes were found to lack introns. Expression of subtilase genes was confirmed at the mRNA level by northern blot analysis and/or by primer extension experiments. For each of the 5 subtilase subfamilies, a distinctive pattern of expression was observed in tomato organs. At least one of the subtilases was found to be expressed in each organ analysed. Structural features evident from deduced amino acid sequences are discussed with reference to the related mammalian proprotein convertases.     

Disruption of the Subtilase Gene, albin1, in Ophiostoma piliferum

Wood sapstaining fungi produce multiple proteases that break down wood protein. Three groups of subtilases have been identified in sapstaining fungi; however, it is not known if these groups have distinct physiological roles (B. Hoffman and C. Breuil, Curr. Genet. 41:168-175, 2002). In this work we examined the role of the subtilase Albin1 from Ophiostoma piliferum. Reamplification of cDNA ends PCR was used to obtain the albin1 gene sequence. The encoded subtilase is probably extracellular and involved in nutrient acquisition. This gene was disrupted with an Agrobacterium tumefaciens-mediated transformation system. Two of the disruptants obtained had significantly lower levels of proteolytic activity, slower growth in bovine serum albumin, and significantly reduced growth on wood. Thus, albin1 plays an important role in O. piliferum's ability to acquire nitrogen from wood proteins.     

Structural features of plant subtilases

Serine proteases of the subtilase family are present in Archaea, Bacteria and Eukarya. Many more subtilases are found in plants as compared to other organisms, implying adaptive significance for the expansion of the subtilase gene family in plants. Structural data, however, were hitherto available only for non-plant subtilases. We recently solved the first structure of a plant subtilase, SlSBT3 from tomato (Solanum lycopersicum). SlSBT3 is a multidomain enzyme displaying a subtilisin, a Protease-Associated (PA) and a fibronectin (Fn) III-like domain. Two prominent features set SlSBT3 apart from other structurally elucidated subtilases: (i) activation by PA domain-mediated homo-dimerization and (ii) calcium-independent activity and thermostability. To address the question whether these characteristics are unique features of SlSBT3, or else, general properties of plant subtilases, homology models were calculated for representative proteases from tomato and Arabidopsis using the SlSBT3 structure as template. We found the major structural elements required for the stabilization of the subtilisin domain to be conserved among all enzymes analyzed. PA domain-mediated dimerization as an auto-regulatory mechanism of enzyme activation, on the other hand, appears to be operating in only a subset of the analyzed subtilasesKey words: Arabidopsis thaliana, enzyme structure, homology modeling, proprotein convertase, protease-associated domain, proteolysis, Solanum lycopersicum, SBT3, subtilisin, thermostabilitySubtilases constitute the second largest family of serine peptidases, both in terms of number of sequences and characterized enzymes (merops.sanger.ac.uk/).¹ The function of subtilases ranges from the non-selective degradation of proteins by e.g., subtilisin Carlsberg in Bacillus licheniformis, to the highly specific maturation of peptide hormones and processing of protein precursors by e.g., kexin in Saccharomyces cerevisiae and proprotein convertases in mammals. In higher plants, subtilases are represented by large gene families comprising 56 members in Arabidopsis thaliana and 63 in rice.^2,3 The expansion of the subtilase family in plants as compared to animals has apparently been accompanied by the acquisition of novel physiological roles that are plant-specific. Plant subtilases were shown to be involved in stomata and seed development,^4,5 in the maintenance of the shoot apical meristem and the cell wall,^6,7 in the processing of peptide growth factors,^8,9 and in responses to the biotic and abiotic environment.^10,11 To address the question whether the adoption of specific roles in plant physiology is reflected in unique structural or biochemical features that distinguish subtilases in plants from those in other organisms, we recently characterized the subtilase SlSBT3 from tomato¹² and solved its structure by X-ray crystallography.^13,14SlSBT3 is an extracellular 79 kDa glycoprotein that exhibits a remarkable level of stability at elevated temperatures and alkaline pH.¹² Like most other subtilases, SlSBT3 is synthesized as a pre-pro-protein and targeted for secretion by an N-terminal signal peptide. Maturation of SlSBT3 involves cleavage of its prodomain, which is a prerequisite for passage through the secretory pathway.¹² Mature SlSBT3 features a protease-associated (PA) domain as a large insertion between the His and Ser active site residues of the protease domain and a fibronectin (Fn) III-like domain as C-terminal extension.¹⁴ This domain architecture (Fig. 1A) is shared with the majority of plant subtilases (e.g., 54 of the 56 subtilases in Arabidopsis). Unlike other structurally elucidated subtilases from bacteria, fungi and animals, SlSBT3 was found to be free of Ca²⁺ in its native state and independent of Ca²⁺ with respect to activity and thermostability. The ability of SlSBT3 to form homodimers is also unique and appears to be critical for enzyme activity and stability.¹⁴ For this addendum, we calculated homology models for representative subtilases from Arabidopsis and tomato to investigate whether PA domain-mediated homo-dimerization and calcium-independent thermostability are unique features of the SlSBT3 structure, or whether they can serve as a first paradigm for the structural biology of plant subtilases in general.Open in a separate windowFigure 1Structural comparison of plant subtilases. (A) Domain architecture of SlSBT3. In addition to the domain borders the three residues constituting the active site are displayed. (B) Structure of the SlSBT3 monomer. Color coding of the domains is like in (A). The bound chloromethylketone (cmk)-inhibitor is shown as ball model in yellow (carbon), red (oxygen) and blue (nitrogen). (C) Functional homodimer of SlSBT3. The region of the direct contact between the two PA domains (gold, purple) is highlighted. In this and all the following panels, a sequence alignment of the relevant regions in SlSBT3 and the modeled subtilases is shown on the right. The sequences highlighted in color were included in the structural alignment on the left. (D) Structure of the partially conserved β-hairpin. (E) Structure of the region corresponding to the conserved calcium-binding site 1 (Ca-1) in thermitase (yellow sticks, PDB code: 1THM). (F) Functional substitution of the conserved Ca-2 (white sticks, PDB code 1S2N) site by a lysine side chain in plant subtilases. (G) Structure of the region corresponding to the less conserved Ca-3 site in thermitase (yellow sticks, PDB code: 1THM). For details, see text.     

Genome-wide and molecular evolution analysis of the subtilase gene family in Vitis vinifera

Background

Vitis vinifera (grape) is one of the most economically significant fruit crops in the world. The availability of the recently released grape genome sequence offers an opportunity to identify and analyze some important gene families in this species. Subtilases are a group of subtilisin-like serine proteases that are involved in many biological processes in plants. However, no comprehensive study incorporating phylogeny, chromosomal location and gene duplication, gene organization, functional divergence, selective pressure and expression profiling has been reported so far for the grape.

Results

In the present study, a comprehensive analysis of the subtilase gene family in V. vinifera was performed. Eighty subtilase genes were identified. Phylogenetic analyses indicated that these subtilase genes comprised eight groups. The gene organization is considerably conserved among the groups. Distribution of the subtilase genes is non-random across the chromosomes. A high proportion of these genes are preferentially clustered, indicating that tandem duplications may have contributed significantly to the expansion of the subtilase gene family. Analyses of divergence and adaptive evolution show that while purifying selection may have been the main force driving the evolution of grape subtilases, some of the critical sites responsible for the divergence may have been under positive selection. Further analyses of real-time PCR data suggested that many subtilase genes might be important in the stress response and functional development of plants.

Conclusions

Tandem duplications as well as purifying and positive selections have contributed to the functional divergence of subtilase genes in V. vinifera. The data may contribute to a better understanding of the grape subtilase gene family.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1116) contains supplementary material, which is available to authorized users.     

A Surface-Exposed Region of a Novel Outer Membrane Protein (P66) of Borrelia spp. Is Variable in Size and Sequence

A model of the 66-kDa outer membrane protein (P66) of Lyme disease Borrelia spp. predicts a surface-exposed loop near the C terminus. This region contains an antigen commonly recognized by sera from Lyme disease patients. In the present study, this region of P66 and homologous proteins of other Borrelia spp. were further investigated by using monoclonal antibodies, epitope mapping of P66 of Borrelia burgdorferi, and DNA sequencing. A monoclonal antibody specific for B. burgdorferi bound to the portion of P66 that was accessible to proteolysis in situ. The linear epitope for the antibody was mapped within a variable segment of the surface-exposed region. To further study this protein, the complete gene of Borrelia hermsii for a protein homologous to P66 was cloned. The deduced protein was 589 amino acids in length and 58% identical to P66 of B. burgdorferi. The B. hermsii P66 protein was predicted to have a surface-exposed region in the same location as that of B. burgdorferi’s P66 protein. With primers designed on the basis of conserved sequences and PCR, we identified and cloned the same regions of P66 proteins of Borrelia turicatae, Borrelia parkeri, Borrelia coriaceae, and Borrelia anserina. The deduced protein sequences from all species demonstrated two conserved hydrophobic regions flanking a surface-exposed loop. The loop sequences were highly variable between different Borrelia spp. in both sequence and size, varying between 35 and 45 amino acids. Although the actual function of P66 of Borrelia spp. is unknown, the results suggest that its surface-exposed region is subject to selective pressure.     

Characterization of Chemotactic Responses and Flagella of Hyphomicrobium Strain W1-1B

Motile swarmer cells of Hyphomicrobium strain W1-1B displayed positive chemotactic responses toward methylamine, dimethylamine, and trimethylamine but did not display significant chemotactic responses towards methanol and arginine. Electron micrographs of negatively stained intact flagellar filaments indicated a novel striated surface pattern. The flagella were composed of two proteins of 39 and 41 kDa. Neither protein was a glycoprotein as determined by Schiff’s staining and by enzyme immunoassay. Protein fingerprints visualized from silver-stained polyacrylamide gels and Western blots of protease-digested samples indicated that the two proteins were similar but not identical. Monoclonal antibodies prepared to the complex flagella of Rhizobium meliloti cross-reacted with the striated flagella of Hyphomicrobium strain W1-1B; however, these antibodies did not cross-react with smooth-surface flagella. These results suggest that complex and striated flagella possess homologous epitope regions.     

Subtilase Genes Diversity in the Biogas Digester Microbiota

Biogas digesters contain microbial assemblages that process a mass of extracellular polymeric substances from animal manure and domestic wastewater; however, due to the limitation of available technology in cultivation of majority of the micro-organisms in biogas digesters, the enzymatic potential of these microbial communities remains largely unexplored. In this study, to evaluate subtilase gene diversity in a biogas digester, the partial sequences of the gene were directly amplified from the metagenomic DNA by using consensus-degenerate primers. The desired PCR products were cloned into pGEM-T Easy vector, and thirty positive clones were chose for Polymerase chain reaction–restriction fragment length polymorphism (PCR–RFLP) analysis, from which thirteen distinguished patterns were obtained and then sequenced. Phylogenetic analysis showed that ten out of the thirteen sequences were related to the subtilase genes in GenBank and were grouped into three families of the subtilases superfamily. The nucleotide sequences analysis through BLAST search revealed that none of the partial genes the authors isolated showed significant similarity against the non-redundant Nucleotide database of NCBI. Meanwhile, the deduced amino acid sequences of ten partial subtilase genes showed moderate identities to the previously identified sequences in GenBank, with a range from 39 to 61%. Collectively, the data indicate that there is a great diversity of subtilase genes in the biogas digester; and may be a rich reservoir for novel subtilase genes.     

Selective recognition of oncogene promoter G-quadruplexes by Mg2+

The epitope sequences within the hemagglutinin (HA) of influenza A virus H3N2 at amino acid residues 173-181 and 227-239 that forms anti-parallel β-sheet structure are similarly recognized by human monoclonal antibodies (HuMAbs), B-1 and D-1 that we recently obtained using the peripheral blood lymphocytes from two influenza-vaccinated volunteers. Both HuMAbs showed strong global neutralization of H3N2 strains. Here we show the significant conservation of the β-sheet region consisting of the above-mentioned two epitope regions in H3N2. In addition, we also identified the corresponding regions with similar structure in other subtypes such as H1N1 and H5N1. These two regions are similarly located underneath the receptor-binding sites of individual subtypes. Analysis of those regions using sequences available from the Influenza Virus Resource at the National Center for Biotechnology Information revealed that compared with those in the known neutralizing epitopes A-E, those sequences were fairly conserved in human H3N2 (n = 7955), swine H1N1 (n = 360) and swine H3N2 (n = 235); and highly conserved in human H1N1 (n = 2722), swine-origin pandemic H1N1 (n = 1474), human H5N1 (n = 319) and avian H5N1 (n = 2349). Phylogenetic tree for these regions formed clearly separable clusters for H1N1, H3N2 and H5N1, irrespective of different host origin. These data may suggest a possible significance of those regions for development of alternative vaccine that could induce neutralizing antibodies reactive against wide-range of influenza virus strains.     

Chimeric fab fragments of antibodies to recombinant Ebola virus glycoprotein

Previously, we have determined the nucleotide and amino acid sequences of the variable domains of three mouse monoclonal antibodies specific to the individual epitopes of the Ebola virus glycoprotein: GPE118 (IgG), GPE325 (IgM) and GPE534 (IgG) [1]. In the present paper, chimeric Fab fragments of Fab118, Fab325, and Fab534 antibodies were obtained based on the variable domains of murine antibodies by attaching CH1 and CL constant regions of human kappa-IgG1 to them. The recombinant chimeric Fab fragments were synthesized in the heterologous expression system Escherichia coli, isolated and purified using metal chelate affinity chromatography. The immunochemical properties of the obtained Fab fragments were studied by immunoblotting techniques as well as indirect and competitive ELISA using recombinant Ebola virus proteins: EBOV rGPdTM (recombinant glycoprotein of Ebola hemorrhagic fever virus without the transmembrane domain), NP (nucleoprotein) and VP40 (structural protein). The identity of recombinant chimeric Fab fragments, as well as their specificity to the recombinant glycoprotein of Ebola hemorrhagic fever virus (EBOV GP) was proved. The results of indirect ELISA evidence the absence of immunological cross-reactivity to NP and VP40 proteins of Ebola virus. The dissociation constants of the antigen-antibody complex K _d equal to 5.0, 1.0 and 1.0 nM for Fab118, Fab325 and Fab534, respectively, were determined; they indicate high affinity of the obtained experimental samples to EBOV GP. The epitope specificity of Fab fragments was studied using a panel of commercial neutralizing antibodies. It was found that all studied antibodies to EBOV GP are targeted to different epitopes, while the epitopes of the recombinant chimeric Fab fragments and original murine monoclonal antibodies (mAbs) coincide. All the obtained and studied mAbs to EBOV GP are specific to epitopes that coincide or overlap the epitopes of three commercial neutralizing mAbs to Ebola virus: epitopes Fab118 and Fab325 overlap the epitope of the known commercial mAb h13F6; Fab325 epitope also overlaps mAb c6D8 epitope; Fab534 epitope is located near mAb KZ52 conformational epitope, in the formation of which amino acid residues of GP1 and GP2 domains of EBOV GP are involved.     

Conservation of proteins involved in oocyst wall formation in Eimeria maxima, Eimeria tenella and Eimeria acervulina

Vaccination with proteins from gametocytes of Eimeria maxima protects chickens, via transfer of maternal antibodies, against infection with several species of Eimeria. Antibodies to E. maxima gametocyte proteins recognise proteins in the wall forming bodies of macrogametocytes and oocyst walls of E. maxima, Eimeria tenella and Eimeria acervulina. Homologous genes for two major gametocyte proteins - GAM56 and GAM82 - were found in E. maxima, E. tenella and E. acervulina. Alignment of the predicted protein sequences of these genes reveals that, as well as sharing regions of tyrosine richness, strong homology exists in their amino-terminal regions, where protective antibodies bind. This study confirms the conservation of the roles of GAM56 and GAM82 in oocyst wall formation and shows that antibodies to gametocyte antigens of E. maxima cross-react with homologous proteins in other species, helping to explain cross-species maternal immunity.     

Evolution of prokaryotic subtilases: genome-wide analysis reveals novel subfamilies with different catalytic residues

Subtilisin-like serine proteases (subtilases) are a very diverse family of serine proteases with low sequence homology, often limited to regions surrounding the three catalytic residues. Starting with different Hidden Markov Models (HMM), based on sequence alignments around the catalytic residues of the S8 family (subtilisins) and S53 family (sedolisins), we iteratively searched all ORFs in the complete genomes of 313 eubacteria and archaea. In 164 genomes we identified a total of 567 ORFs with one or more of the conserved regions with a catalytic residue. The large majority of these contained all three regions around the "classical" catalytic residues of the S8 family (Asp-His-Ser), while 63 proteins were identified as S53 (sedolisin) family members (Glu-Asp-Ser). More than 30 proteins were found to belong to two novel subsets with other evolutionary variations in catalytic residues, and new HMMs were generated to search for them. In one subset the catalytic Asp is replaced by an equivalent Glu (i.e. Glu-His-Ser family). The other subset resembles sedolisins, but the conserved catalytic Asp is not located on the same helix as the nucleophile Glu, but rather on a beta-sheet strand in a topologically similar position, as suggested by homology modeling. The Prokaryotic Subtilase Database (www.cmbi.ru.nl/subtilases) provides access to all information on the identified subtilases, the conserved sequence regions, the proposed family subdivision, and the appropriate HMMs to search for them. Over 100 proteins were predicted to be subtilases for the first time by our improved searching methods, thereby improving genome annotation.     

Pre-clinical assessment of novel multivalent MSP3 malaria vaccine constructs

Background

MSP3 has been shown to induce protection against malaria in African children. The characterization of a family of Plasmodium falciparum merozoite surface protein 3 (MSP3) antigens sharing a similar structural organization, simultaneously expressed on the merozoite surface and targeted by a cross-reactive network of protective antibodies, is intriguing and offers new perspectives for the development of subunit vaccines against malaria.

Methods

Eight recombinant polyproteins containing carefully selected regions of this family covalently linked in different combinations were all efficiently produced in Escherichia coli. The polyproteins consisted of one monovalent, one bivalent, one trivalent, two tetravalents, one hexavalent construct, and two tetravalents incorporating coiled-coil repeats regions from LSA3 and p27 vaccine candidates.

Results

All eight polyproteins induced a strong and homogeneous antibody response in mice of three distinct genotypes, with a dominance of cytophilic IgG subclasses, lasting up to six months after the last immunization. Vaccine-induced antibodies exerted a strong monocyte-mediated in vitro inhibition of P. falciparum growth. Naturally acquired antibodies from individuals living in an endemic area of Senegal recognized the polyproteins with a reactivity mainly constituted of cytophilic IgG subclasses.

Conclusions

Combination of genetically conserved and antigenically related MSP3 proteins provides promising subunit vaccine constructs, with improved features as compared to the first generation construct employed in clinical trials (MSP3-LSP). These multivalent MSP3 vaccine constructs expand the epitope display of MSP3 family proteins, and lead to the efficient induction of a wider range of antibody subclasses, even in genetically different mice. These findings are promising for future immunization of genetically diverse human populations.     

MAPPING OF HIV-1 GAG EPITOPES RECOGNIZED BY POLYCLONAL ANTIBODIES USING GENE-FRAGMENT PHAGE DISPLAY SYSTEM

Phage display has emerged as a powerful technique for mapping epitopes recognised by monoclonal and polyclonal antibodies. We have recently developed a simple gene-fragment phage display system and have shown its utility in mapping epitope recognised by a monoclonal antibody. In the present study, we have employed this system in mapping epitopes recognised by polyclonal antibodies raised against HIV-1 capsid protein, p24 which is derived from proteolytic cleavage of Gag polyprotein. HIV-1 gag DNA was fragmented by DNase I and the fragments (50–250 bp) were cloned into gene-fragment phage display vector to construct a library of phages displaying peptides. This phage library was used for affinity selection of phages displaying epitopes recognised by rabbit anti-p24 polyclonal antibodies. Selected phages contained sequences from two discrete regions of p24, demonstrating the presence of two antigenic regions.

The DNA sequences encoding these regions were also cloned and expressed as GST fusion proteins. The immunoreactivity of these epitopes as GST fusion proteins, or as phage-displayed peptides, was comparable in ELISA system using same anti-p24 polyclonal antibodies. The results indicate that the gene-fragment based phage display system can be used efficiently to identify epitopes recognised by polyclonal antibodies, and phage displayed epitopes can be directly employed in ELISA to detect antibodies.     

Chromosomal Mapping of Repetitive DNAs in the Grasshopper Abracris flavolineata Reveal Possible Ancestry of the B Chromosome and H3 Histone Spreading

Supernumerary chromosomes (B chromosomes) occur in approximately 15% of eukaryote species. Although these chromosomes have been extensively studied, knowledge concerning their specific molecular composition is lacking in most cases. The accumulation of repetitive DNAs is one remarkable characteristic of B chromosomes, and the occurrence of distinct types of multigene families, satellite DNAs and some transposable elements have been reported. Here, we describe the organization of repetitive DNAs in the A complement and B chromosome system in the grasshopper species Abracris flavolineata using classical cytogenetic techniques and FISH analysis using probes for five multigene families, telomeric repeats and repetitive C₀t-1 DNA fractions. The 18S rRNA and H3 histone multigene families are highly variable and well distributed in A. flavolineata chromosomes, which contrasts with the conservation of U snRNA genes and less variable distribution of 5S rDNA sequences. The H3 histone gene was an extensively distributed with clusters occurring in all chromosomes. Repetitive DNAs were concentrated in C-positive regions, including the pericentromeric region and small chromosomal arms, with some occurrence in C-negative regions, but abundance was low in the B chromosome. Finally, the first demonstration of the U2 snRNA gene in B chromosomes in A. flavolineata may shed light on its possible origin. These results provide new information regarding chromosomal variability for repetitive DNAs in grasshoppers and the specific molecular composition of B chromosomes.     

Antipeptide Antibodies That Can Distinguish Specific Subunit Polypeptides of Glutamine Synthetase from Bean (Phaseolus vulgaris L.)

The amino acid sequences of the β and γ subunit polypeptides of glutamine synthetase from bean (Phaseolus vulgaris L.) root nodules are very similar. However, there are small regions within the sequences that are significantly different between the two polypeptides. The sequences between amino acids 2 and 9 and between 264 and 274 are examples. Three peptides (γ_2-9, γ_264-274, and β_264-274) corresponding to these sequences were synthesized. Antibodies against these peptides were raised in rabbits and purified with corresponding peptide-Sepharose affinity chromatography. Western blot analysis of polyacrylamide gel electrophoresis of bean nodule proteins demonstrated that the anti-β_264-274 antibodies reacted specifically with the β polypeptide and the anti-γ_264-274 and anti-γ_2-9 antibodies reacted specifically with the γ polypeptide of the native and denatured glutamine synthetase. These results showed the feasibility of using synthetic peptides in developing antibodies that are capable of distinguishing proteins with similar primary structures.