The rhomboids: a nearly ubiquitous family of intramembrane serine proteases that probably evolved by multiple ancient horizontal gene transfers

Background

The rhomboid family of polytopic membrane proteins shows a level of evolutionary conservation unique among membrane proteins. They are present in nearly all the sequenced genomes of archaea, bacteria and eukaryotes, with the exception of several species with small genomes. On the basis of experimental studies with the developmental regulator rhomboid from Drosophila and the AarA protein from the bacterium Providencia stuartii, the rhomboids are thought to be intramembrane serine proteases whose signaling function is conserved in eukaryotes and prokaryotes.

Results

Phylogenetic tree analysis carried out using several independent methods for tree constructions and the corresponding statistical tests suggests that, despite its broad distribution in all three superkingdoms, the rhomboid family was not present in the last universal common ancestor of extant life forms. Instead, we propose that rhomboids evolved in bacteria and have been acquired by archaea and eukaryotes through several independent horizontal gene transfers. In eukaryotes, two distinct, ancient acquisitions apparently gave rise to the two major subfamilies, typified by rhomboid and PARL (presenilins-associated rhomboid-like protein), respectively. Subsequent evolution of the rhomboid family in eukaryotes proceeded by multiple duplications and functional diversification through the addition of extra transmembrane helices and other domains in different orientations relative to the conserved core that harbors the protease activity.

Conclusions

Although the near-universal presence of the rhomboid family in bacteria, archaea and eukaryotes appears to suggest that this protein is part of the heritage of the last universal common ancestor, phylogenetic tree analysis indicates a likely bacterial origin with subsequent dissemination by horizontal gene transfer. This emphasizes the importance of explicit phylogenetic analysis for the reconstruction of ancestral life forms. A hypothetical scenario for the origin of intracellular membrane proteases from membrane transporters is proposed.

     

Drosophila rhomboid-1 defines a family of putative intramembrane serine proteases.

The polytopic membrane protein Rhomboid-1 promotes the cleavage of the membrane-anchored TGFalpha-like growth factor Spitz, allowing it to activate the Drosophila EGF receptor. Until now, the mechanism of this key signaling regulator has been obscure, but our analysis suggests that Rhomboid-1 is a novel intramembrane serine protease that directly cleaves Spitz. In accordance with the putative Rhomboid active site being in the membrane bilayer, Spitz is cleaved within its transmembrane domain, and thus is, to our knowledge, the first example of a growth factor activated by regulated intramembrane proteolysis. Rhomboid-1 is conserved throughout evolution from archaea to humans, and our results show that a human Rhomboid promotes Spitz cleavage by a similar mechanism. This growth factor activation mechanism may therefore be widespread.     

The HtrA family of serine proteases

HtrA, also known as DegP and probably identical to the Do protease, is a heat shock-induced serine protease that is active in the periplasm of Escherichia coli . Homologues of HtrA have been described in a wide range of bacteria and in eukaryotes. Its chief role is to degrade misfolded proteins in the periplasm. Substrate recognition probably involves the recently described PDZ domains in the C-terminal half of HtrA and, we suspect, has much in common with the substrate recognition system of the tail-specific protease, Prc (which also possesses a PDZ domain). The expression of htrA is regulated by a complex set of signal transduction pathways, which includes an alternative sigma factor, RpoE, an anti-sigma factor, RseA, a two-component regulatory system, CpxRA, and two phosphoprotein phosphatases, PrpA and PrpB. Mutations in the htrA genes of Salmonella , Brucella and Yersinia cause decreased survival in mice and/or macrophages, and htrA mutants can act as vaccines, as cloning hosts and as carriers of heterologous antigens.     

Granzymes: a family of lymphocyte granule serine proteases

Granzymes, a family of serine proteases, are expressed exclusively by cytotoxic T lymphocytes and natural killer (NK) cells, components of the immune system that protect higher organisms against viral infection and cellular transformation. Following receptor-mediated conjugate formation between a granzyme-containing cell and an infected or transformed target cell, granzymes enter the target cell via endocytosis and induce apoptosis. Granzyme B is the most powerful pro-apoptotic member of the granzyme family. Like caspases, cysteine proteases that play an important role in apoptosis, it can cleave proteins after acidic residues, especially aspartic acid. Other granzymes may serve additional functions, and some may not induce apoptosis. Granzymes have been well characterized only in human and rodents, and can be grouped into three subfamilies according to substrate specificity: members of the granzyme family that have enzymatic activity similar to the serine protease chymotrypsin are encoded by a gene cluster termed the 'chymase locus'; granzymes with trypsin-like specificities are encoded by the 'tryptase locus'; and a third subfamily cleaves after unbranched hydrophobic residues, especially methionine, and is encoded by the 'Met-ase locus'. All granzymes are synthesized as zymogens and, after clipping of the leader peptide, maximal enzymatic activity is achieved by removal of an amino-terminal dipeptide. They can all be blocked by serine protease inhibitors, and a new group of inhibitors has recently been identified - serpins, some of which are specific for granzymes. Future studies of serpins may bring insights into how cells that synthesize granzymes are protected from inadvertent cell suicide.     

ClpP: a distinctive family of cylindrical energy-dependent serine proteases

Processes maintaining protein homeostasis in the cell are governed by the activities of molecular chaperones that mainly assist in the folding of polypeptide chains and by a large class of proteases that regulate protein levels through degradation. ClpP proteases define a distinctive family of cylindrical, energy-dependent serine proteases that are highly conserved throughout bacteria and eukaryota. They typically interact with ATP-dependent AAA+ chaperones that bind and unfold target substrates and then translocate them into ClpP for degradation. Structural and functional studies have provided a detailed view of the mechanism of function of this class of proteases.     

Clp P represents a unique family of serine proteases

The amino acid sequence of Clp P, the proteolytic subunit of the ATP-dependent Clp protease of Escherichia coli, closely resembles a protein encoded by chloroplast DNA, which is well conserved between chloroplasts of different plant species. The homology extends over almost the full length of the sequences of both proteins and consists of approximately 46% identical and approximately 70% similar amino acids. Antibodies against E. coli Clp P cross-reacted with proteins with Mr of 20,000-30,000 in bacteria, lower eukaryotes, plants, and animal cells. Since the regulatory subunit of Clp protease, Clp A, also has a homolog in plants, as well as in other bacteria and in lower eukaryotes, it is likely that ATP-dependent proteolysis in chloroplasts is catalyzed in part by a Clp-like protease and that both components of Clp-like proteases are widespread in living cells. We have identified Ser-111 as the active site serine in E. coli Clp P modified by diisopropyl fluorophosphate. Mutational alteration of Ser-111 or His-136 eliminates proteolytic activity of Clp P. Both residues are found in highly conserved regions of the protein. The sequences around the active site residues suggest that Clp P represents a unique class of serine protease. Amino-terminal processing of cloned Clp P mutated at either Ser-111 or His-136 occurs efficiently when wild-type clpP is present in the chromosome but is blocked in clpP- hosts. Processing of Clp P appears, therefore, to involve an intermolecular autocatalytic cleavage reaction. Since processing of Clp P occurs in clpA- cells, the autoprocessing activity of Clp P is independent of Clp A.     

Identification of multiple independent horizontal gene transfers into poxviruses using a comparative genomics approach

Background  

Poxviruses are important pathogens of humans, livestock and wild animals. These large dsDNA viruses have a set of core orthologs whose gene order is extremely well conserved throughout poxvirus genera. They also contain many genes with sequence and functional similarity to host genes which were probably acquired by horizontal gene transfer.     

The mycosins of Mycobacterium tuberculosis H37Rv: a family of subtilisin-like serine proteases

There is little information regarding the role of proteolysis in Mycobacterium tuberculosis and no studies on the potential involvement of proteases in the pathogenesis of tuberculosis. We identified five M. tuberculosis genes (mycP1-5) that encode a family of serine proteases (mycosins-1 to 5), ranging from 36 to 47% identity. Each protein contains a catalytic triad (Asp, His, Ser) within highly conserved sequences, typical of proteases of the subtilisin family. These genes are also present in M. bovis BCG and other virulent mycobacteria, but only one homologue (mycP3) was detected in M. smegmatis. The mycosins have N-terminal signal sequences and C-terminal transmembrane anchors, and the localisation of the mycosins to the membrane/cell wall was verified by Western blot analysis of heterologously expressed proteins in cellular fractions of M. smegmatis. In M. tuberculosis, all the mycosins were expressed constitutively during growth in broth. Mycosins-2 and 3 were also expressed constitutively in M. bovis BCG, but no expression of mycosin-1 was detected. Mycosin-2 was modified by cleavage in all three mycobacterial species. The multiplicity and constitutive expression of these proteins suggests that they have an important role in the biology of M. tuberculosis.     

Core principles of intramembrane proteolysis: comparison of rhomboid and site-2 family proteases

Cleavage of proteins within their membrane-spanning segments is an ancient regulatory mechanism that has evolved to control a myriad of cellular processes in all forms of life. Although three mechanistic families of enzymes have been discovered that catalyze hydrolysis within the water-excluding environment of the membrane, how they achieve this improbable reaction has been both a point of controversy and skepticism. The crystal structures of rhomboid and site-2 protease, two different classes of intramembrane proteases, have been solved recently. Combined with current biochemical analyses, this advance provides an unprecedented view of how nature has solved the problem of facilitating hydrolysis within membranes in two independent instances. We focus on detailing the similarities between these unrelated enzymes to define core biochemical principles that govern this conserved regulatory mechanism.     

A despecialization step underlying evolution of a family of serine proteases

In the trypsin superfamily of serine proteases, non-trypsin-like primary specificities have arisen in only two monophyletic descendent subbranches. We have recreated an ancestor to one of these subbranches (granzyme) using phylogenetic inference, gene synthesis, and protein expression. This ancestor has two unusual properties. First, it has broad primary specificity encompassing the entire repertoire of novel primary specificities found in its descendents. Second, unlike extant members that have narrow primary specificities, the ancestor exhibits tolerance to mutational changes in primary specificity-conferring residues-that is, structural plasticity. Molecular modeling and mutagenesis studies indicate that these unusual properties are due to a particularly wide substrate binding pocket. These two crucial properties of the ancestor not only distinguish it from its extant descendents but also from the trypsin-like proteases that preceded it. This indicates that a despecialization step, characterized by broad specificity and structural plasticity, underlies evolution of new primary specificities in this protease superfamily.     

Comparative analysis of a genome fragment of an uncultivated mesopelagic crenarchaeote reveals multiple horizontal gene transfers

Marine planktonic crenarchaeota have escaped all cultivation attempts to date, all crenarchaeota growing in pure culture so far being hyperthermophiles. Here, we present a comparative genomic analysis of a 16S- plus 23S-rDNA-containing fragment of a crenarchaeote retrieved from an environmental genomic library constructed from picoplankton collected at 500 m depth in the Antarctic Polar Front. The clone DeepAnt-EC39 contained an insert of 33.3 kbp, which was completely sequenced. DeepAnt-EC39 appears to represent a lineage specific to deep-sea waters but widespread geographically, as revealed by the analysis of the 16S-23S-rDNA intergenic spacer region. A comparison with previously sequenced marine crenarchaeotal genomic clones also containing an rrn operon (74A4, 4B7 and Cenarchaeum symbiosum strains A and B) revealed a highly variable structure involving gene rearrangements and insertions/deletions. The surroundings of the rrn operon and the contiguous glutamate-1-semialdehyde aminotransferase gene appear hot spots for recombination. Phylogenetic analyses of all individual predicted proteins revealed the existence of several likely cases of horizontal gene transfer both, between the two archaeal kingdoms and between the two prokaryotic domains. The most frequent horizontal transfers appear to involve genes from mesophilic methanogenic euryarchaeota related to Methanosarcinales. We hypothesise that the acquisition of genes from mesophilic bacteria and euryarchaeota has played a major role in the adaptation of Group I crenarchaeota to life at lower temperatures.     

GlyGly-CTERM and rhombosortase: a C-terminal protein processing signal in a many-to-one pairing with a rhomboid family intramembrane serine protease

The rhomboid family of serine proteases occurs in all domains of life. Its members contain at least six hydrophobic membrane-spanning helices, with an active site serine located deep within the hydrophobic interior of the plasma membrane. The model member GlpG from Escherichia coli is heavily studied through engineered mutant forms, varied model substrates, and multiple X-ray crystal studies, yet its relationship to endogenous substrates is not well understood. Here we describe an apparent membrane anchoring C-terminal homology domain that appears in numerous genera including Shewanella, Vibrio, Acinetobacter, and Ralstonia, but excluding Escherichia and Haemophilus. Individual genomes encode up to thirteen members, usually homologous to each other only in this C-terminal region. The domain's tripartite architecture consists of motif, transmembrane helix, and cluster of basic residues at the protein C-terminus, as also seen with the LPXTG recognition sequence for sortase A and the PEP-CTERM recognition sequence for exosortase. Partial Phylogenetic Profiling identifies a distinctive rhomboid-like protease subfamily almost perfectly co-distributed with this recognition sequence. This protease subfamily and its putative target domain are hereby renamed rhombosortase and GlyGly-CTERM, respectively. The protease and target are encoded by consecutive genes in most genomes with just a single target, but far apart otherwise. The signature motif of the Rhombo-CTERM domain, often SGGS, only partially resembles known cleavage sites of rhomboid protease family model substrates. Some protein families that have several members with C-terminal GlyGly-CTERM domains also have additional members with LPXTG or PEP-CTERM domains instead, suggesting there may be common themes to the post-translational processing of these proteins by three different membrane protein superfamilies.     

Comparative modeling methods: application to the family of the mammalian serine proteases

Comparative modeling methods are described that can be used to construct a three-dimensional model structure of a new protein from knowledge of its sequence and of the experimental structures and sequences of other members of its homology family. The methods are illustrated with the mammalian serine protease family, for which seven experimental structures have been reported in the literature, and the sequences for over 35 different protein members of the family are available. The strategy for modeling these proteins is presented, and criteria are developed for determining and assigning the reliability of the modeled structure. Criteria are described that are specially designed to help detect cases in which it is likely that the local structure diverges significantly from the usual conformation of the family.     

The structural basis of mode of activation and functional diversity: a case study with HtrA family of serine proteases

HtrA (High temperature requirement protease A) proteins that are primarily involved in protein quality control belong to a family of serine proteases conserved from bacteria to humans. HtrAs are oligomeric proteins that share a common trimeric pyramidal architecture where each monomer comprises a serine protease domain and one or two PDZ domains. Although the overall structural integrity is well maintained and they exhibit similar mechanism of activation, subtle conformational changes and structural plasticity especially in the flexible loop regions and domain interfaces lead to differences in their active site conformation and hence in their specificity and functions.     

Characterization of a ubiquitous expressed gene family encoding polygalacturonase in Arabidopsis thaliana

Pectin, as one of the major components of plant cell wall, has been implicated in many developmental processes occurring during plant growth. Among the different enzymes known to participate in the pectin structure modifications, polygalacturonase (PG) activity has been shown to be associated with fruit ripening, organ abscission and pollen grain development. Until now, sequence analyses of the deduced polypeptides of the plant PG genes allowed their grouping into three clades corresponding to genes involved in one of these three activities. In this study, we report the sequence of three genomic clones encoding PG in Arabidopsis thaliana. These genes, together with 16 other genes present in the databases form a large gene family, ubiquitously expressed, present on the five chromosomes with at least two gene clusters on chromosomes II and V, respectively. Phylogenetic analyses suggest that the A. thaliana gene family contains five classes of genes, with three of them corresponding to the previously defined clades. Comparison of positions and numbers of introns among the A. thaliana genes reveals structural conservation between genes belonging to the same class. The pattern of intron losses that could have given rise to the PG gene family is consistent with a mechanism of intron loss by replacement of an ancestral intron-containing gene with a reverse-transcribed DNA copy of a spliced mRNA. Following this event of intron loss, the acquisition of introns in novel positions is consistent with a mechanism of intron gain at proto-splice sites.     

HGTector: an automated method facilitating genome-wide discovery of putative horizontal gene transfers

Background

First pass methods based on BLAST match are commonly used as an initial step to separate the different phylogenetic histories of genes in microbial genomes, and target putative horizontal gene transfer (HGT) events. This will continue to be necessary given the rapid growth of genomic data and the technical difficulties in conducting large-scale explicit phylogenetic analyses. However, these methods often produce misleading results due to their inability to resolve indirect phylogenetic links and their vulnerability to stochastic events.

Results

A new computational method of rapid, exhaustive and genome-wide detection of HGT was developed, featuring the systematic analysis of BLAST hit distribution patterns in the context of a priori defined hierarchical evolutionary categories. Genes that fall beyond a series of statistically determined thresholds are identified as not adhering to the typical vertical history of the organisms in question, but instead having a putative horizontal origin. Tests on simulated genomic data suggest that this approach effectively targets atypically distributed genes that are highly likely to be HGT-derived, and exhibits robust performance compared to conventional BLAST-based approaches. This method was further tested on real genomic datasets, including Rickettsia genomes, and was compared to previous studies. Results show consistency with currently employed categories of HGT prediction methods. In-depth analysis of both simulated and real genomic data suggests that the method is notably insensitive to stochastic events such as gene loss, rate variation and database error, which are common challenges to the current methodology. An automated pipeline was created to implement this approach and was made publicly available at: https://github.com/DittmarLab/HGTector. The program is versatile, easily deployed, has a low requirement for computational resources.

Conclusions

HGTector is an effective tool for initial or standalone large-scale discovery of candidate HGT-derived genes.

Electronic supplementary material

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

Enhanced function annotations for Drosophila serine proteases: a case study for systematic annotation of multi-member gene families

Systematically annotating function of enzymes that belong to large protein families encoded in a single eukaryotic genome is a very challenging task. We carried out such an exercise to annotate function for serine-protease family of the trypsin fold in Drosophila melanogaster, with an emphasis on annotating serine-protease homologues (SPHs) that may have lost their catalytic function. Our approach involves data mining and data integration to provide function annotations for 190 Drosophila gene products containing serine-protease-like domains, of which 35 are SPHs. This was accomplished by analysis of structure-function relationships, gene-expression profiles, large-scale protein-protein interaction data, literature mining and bioinformatic tools. We introduce functional residue clustering (FRC), a method that performs hierarchical clustering of sequences using properties of functionally important residues and utilizes correlation co-efficient as a quantitative similarity measure to transfer in vivo substrate specificities to proteases. We show that the efficiency of transfer of substrate-specificity information using this method is generally high. FRC was also applied on Drosophila proteases to assign putative competitive inhibitor relationships (CIRs). Microarray gene-expression data were utilized to uncover a large-scale and dual involvement of proteases in development and in immune response. We found specific recruitment of SPHs and proteases with CLIP domains in immune response, suggesting evolution of a new function for SPHs. We also suggest existence of separate downstream protease cascades for immune response against bacterial/fungal infections and parasite/parasitoid infections. We verify quality of our annotations using information from RNAi screens and other evidence types. Utilization of such multi-fold approaches results in 10-fold increase of function annotation for Drosophila serine proteases and demonstrates value in increasing annotations in multiple genomes.