Mining the <it>Arabidopsis thaliana </it>genome for highly-divergent seven transmembrane receptors

gb-2006-7-10-r96 GBJ Method Mining the <it>Arabidopsis thaliana </it>genome for highly-divergent seven transmembrane receptors Moriyama N Etsuko emoriyama2@unl.edu Strope K Pooja poojastrope@gmail.com Opiyo O Stephen sopiyo@unlserve.unl.edu Chen Zhongying zchen@email.unc.edu Jones M Alan alanjones@bio.unc.edu

School of Biological Sciences and Plant Science Initiative, University of Nebraska-Lincoln, Lincoln, NE 68588-0660, USA

Department of Agronomy and Horticulture, University of Nebraska-Lincoln, Lincoln, NE 68583-0915, USA

Departments of Biology and Pharmacology, University of North Carolina at Chapel Hill, Chapel Hill, NC 27599, USA

Genome Biology 1465-6906 2006 7 10 R96 http://genomebiology.com/2006/7/10/R96 17064408 10.1186/gb-2006-7-10-r96 28 6 2006 24 8 2006 25 10 2006 25 10 2006 2006 Moriyama et al.; licensee BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Arabidopsis putative seven transmembrane proteins

A combination of multiple protein classification methods is described and used to identify a minimum set of 54 candidate seven transmembrane receptors in Arabidopsis thaliana.

Abstract

To identify divergent seven-transmembrane receptor (7TMR) candidates from the Arabidopsis thaliana genome, multiple protein classification methods were combined, including both alignment-based and alignment-free classifiers. This resolved problems in optimally training individual classifiers using limited and divergent samples, and increased stringency for candidate proteins. We identified 394 proteins as 7TMR candidates and highlighted 54 with corresponding expression patterns for further investigation.

Bioinformatics Plant biology

Background

Seven-transmembrane (7TM)-region containing proteins constitute the largest receptor superfamily in vertebrates and other metazoans. These cell-surface receptors are activated by a diverse array of ligands, and are involved in various signaling processes, such as cell proliferation, neurotransmission, metabolism, smell, taste, and vision. They are the central players in eukaryotic signal transduction. They are commonly referred to as G protein-coupled receptors (GPCRs) because most transduce extracellular signals into cellular physiological responses through the activation of heterotrimeric guanine nucleotide binding proteins (G proteins) 1. However, an increasing number of alternative 'G protein-independent' signaling mechanisms have been associated with groups of these 7TM proteins 2345. Thus, for precision and clarity, we refer to these proteins simply as 7TM receptors (7TMRs), and candidate proteins in organisms greatly divergent to humans are designated here as 7TM putative receptors (7TMpRs).

The human genome encodes approximately 800 or more 7TMRs, both with and without known cognate ligands (the latter are so-called orphan GPCRs); they thus constitute >1% of the gene complement 67. More than 1,000 genes or 5% of the Caenorhabditis elegans genome are predicted to encode 7TMRs; the majority of them appear to be chemoreceptors 8. Approximately 300 7TMR-encoding genes (about 1% to 2% of the genome) have been recognized in the Drosophila melanogaster genome 67. Compared to such large numbers of 7TMRs found in animal genomes, very few 7TMpRs have been reported in plants and fungi. Only 22 Arabidopsis 7TMpRs have been described so far. Fifteen of them constitute the 'mildew resistance locus O' (MLO) family, whose direct interaction with the G-protein α subunit (Gα) has not been shown 910. While another 7TMpR, GCR1 11, directly interacts with the plant Gα subunit GPA1 12, it has been shown that GCR1 can act independently of the heterotrimeric G-protein complex as well 2. Hsieh and Goodman 13 recently reported five expressed proteins predicted to have 7TM regions (heptahelical transmembrane proteins (HHPs) 1 to 5) but these, like the other 16, do not have candidate ligands. Finally, an unusual Regulator of G Signaling (RGS) protein (designated AtRGS1) has been predicted to have 7TM regions 14. RGS proteins function as a GTPase activating protein (GAP) to de-sensitize signaling by de-activating the Gα subunits of the heterotrimeric complex. Because Arabidopsis seedlings lacking AtRGS1 have reduced sensitivity to D-glucose 21415, the possibility exists that AtRGS1 is a novel D-glucose receptor having an agonist-regulated GAP function. Although we designate them 7TMpRs here, it should be noted that neither a ligand nor a full signaling cascade has been demonstrated yet for any of these plant proteins, and only for a barley MLO protein has the 7TM topology been experimentally confirmed 9.

None of the reported Arabidopsis 7TMpR proteins share substantial sequence similarity with known metazoan GPCRs constituting six different subfamilies. It appears that plant 7TMpRs dramatically diverged from known metazoan GPCRs over the 1.6 billion years since the plant and metazoan lineages bifurcated. It should be noted that Arabidopsis GCR1 shares weak but significant similarity with the cyclic AMP receptor, CAR1, found in the slime mold 21116. There is also very weak similarity to the Class B Secretin family GPCRs. However, other than GCR1, currently used search methods have not robustly identified plant 7TMpR proteins as candidate GPCRs. This great sequence divergence highlights the need for new approaches to identify divergent 7TMR candidates in non-metazoan genomes.

The human genome contains 16 Gα, 5 Gβ, and 12 Gγ genes. In stark contrast, both fungi and plants have much simpler G-protein coupled signaling systems. For example, the Arabidopsis genome contains one canonical Gα, one Gβ, and two Gγ genes 17. Similarly, a small number of G-proteins are found in fungi; there are two Gα, one Gβ, and one Gγ in Saccharomyces cerevisiae 181920 while Neurospora crassa and some fungi have more genes encoding each subunit 212223. Therefore, it may be reasonable to assume that plants and fungi have fewer GPCRs than human, and while approximately 200 Arabidopsis proteins were predicted to have 7TM regions, sequence divergence precludes unequivocal assignment of any as an orphan GPCR 2425. However, at least 61 7TMpRs have been recently predicted from the plant pathogenic fungus Magnaporthe grisea genome 26, raising the possibility that more divergent groups of 7TMpR proteins likely remain undiscovered in non-metazoan taxa.

In this report, we describe our comprehensive computational strategy for identifying 7TMpR candidates from the entire protein sequence set predicted from the A. thaliana genome, and compile their tissue-specific expression and co-expression patterns with G-proteins. To take advantage of different approaches, we combined multiple protein classification methods, including more specific (conservative) alignment-based classifiers and more sensitive alignment-free classifiers, to predict candidate 7TMpRs in divergent genomes more effectively.

Results and discussion

Identifying 7TMpR candidates using various protein classification methods

Among many protein classification methods commonly used, the current state-of-the-art and most used is the profile hidden Markov models (profile HMMs) 27. It is used to construct protein family databases such as Pfam 2829, SMART 3031, and Superfamily 32. However, profile HMMs and other currently used classification methods such as PROSITE 3334 and PRINTS 3536 share an important weakness. These methods rely on multiple alignments for generating their models (patterns, profile HMMs, and so on). Generating robust multiple alignments is difficult or impossible when extremely diverged sequences are included in the analysis; 7TMRs are one such protein family whose sequence similarities between subgroups can be lower than 25%. Furthermore, alignments are generated only from known related proteins (positive samples), and, therefore, no information from negative samples (unrelated protein sequences) is directly incorporated in the model building process. Identifiable 'hits' are, therefore, constrained by initial sampling bias, which becomes reinforced when models are iteratively rebuilt from accumulated sequences. Consequently, the predictive power, especially the sensitivity, of these classifiers decreases when they are applied against extremely diverged protein families.

To overcome this disadvantage and to increase sensitivities against such non-alignable similarities, several 'alignment-free' methods have been proposed recently. These methods quantify various properties of amino acid sequences and convert them into a descriptor array. Once multiple sequences with different lengths are transformed into a uniform matrix, various multivariate analysis methods can be applied. Kim et al. 37 and Moriyama and Kim 38 used parametric and non-parametric discriminant function analysis methods. Karchin et al. 39 incorporated profile HMMs with support vector machines (SVMs) using the Fisher kernel (SVM-Fisher) so that negative sample information can be taken into account when training the classifier. SVMs can be applied with completely 'alignment-free' sequence descriptors, for example, amino acid and dipeptide compositions. Such alignment-free classifiers are shown to outperform profile HMMs as well as Karchin et al.'s SVM-Fisher 4041 (PK Strope and EN Moriyama, submitted). Another multivariate method, partial least squares (PLS) regression, was used by Lapinsh et al. 42 with physico-chemical properties of amino acids. We recently re-evaluated the descriptors used with PLS and optimized them to discriminate 7TMRs from other proteins 43.

We applied these methods against the entire predicted protein sequence set derived from the A. thaliana genome. As shown in Table 1, among the 28,952 protein sequences, the Sequence Alignment and Modeling system (SAM), a profile HMM method, predicted only 16 (excluding one alternatively spliced gene sequence) as 7TMpR candidates. Fifteen of them are identified as MLO or similar to MLO and one as GCR1 in The Arabidopsis Information Resource (TAIR) 4445. It clearly shows that SAM is highly specific (discriminating) with no false positive, assuming that current annotations are correct. SAM failed to identify only one known MLO (MLO4: At1g11000). This protein, as well as AtRGS1 and five recently predicted 7TM proteins (HHP1-5), were among the 16 previously predicted Arabidopsis 7TMpRs not included in the randomly sampled 500 7TMR training sequences (see Materials and methods). Thus, we concluded that the predictive power of SAM alone is insufficient to identify highly diverged and potentially novel 7TMpR sequences.

Table 1

Numbers of 7TMpR candidates identified by various methods from the A. thaliana genome

Methods

Number of 7TMpR candidates*

HMMTOP

7TMs^†

236 (201)

6-8 TM^†

633 (545)

5-9 TMs^†

1,091 (957)

5-10 TMs^†

1,343 (1,179)

SAM

16 (15)

LDA

3,211 (2,935)

QDA

2,006 (1,820)

LOG

2,626 (2,394)

KNN (K = 5)

3,125 (2,839)

KNN (K = 10)

3,202 (2,906)

KNN (K = 15)

3,298 (3,004)

KNN (K = 20)

3,347 (3,043)

SVM-AA

2,263 (2,043)

SVM-di

2,004 (1,807)

PLS-ACC

2,671 (2,466)

*The numbers in parentheses show 7TMpR candidates after removing proteins derived from alternative splicing. ^†The numbers of TM regions predicted by HMMTOP.

The results obtained by SAM were compared with those obtained by alignment-free methods. As shown in Table 1, alignment-free methods (LDA, QDA, LOG, KNN, SVM with amino acid composition (SVM-AA), SVM with dipeptide composition (SVM-di), and PLS with amino acid properties (PLS-ACC)) predicted 2,000 to 3,400 proteins as 7TMpR candidates, which is about 10% of the entire predicted Arabidopsis proteome and about 30% to 50% of all possible transmembrane proteins (6,475 proteins) 2425. These alignment-free methods clearly call many false positives, and need further optimization to improve their discrimination power.

One advantage of alignment-free methods to be noted is their sensitivity against short or partial sequences 3738. Many of the 28,952 protein sequences used in this study are based only on ab initio gene prediction results, and hence are likely to contain various types of errors. If only a part of a 7TMR protein is predicted correctly, alignment-free methods could have a better chance to identify it.

Table 1 lists Arabidopsis proteins that were predicted to have five to ten transmembrane regions and bins them by the number of transmembrane regions. HMMTOP 2.0 4647 predicted 201 proteins as having 7TM regions. This number is close to a previous prediction (184 proteins) 2425. We should note, however, that no single method predicts 7TM regions from all known 7TMRs exactly (see Materials and methods). As mentioned above, it is also possible that some deduced Arabidopsis proteins we analyzed do not contain the entire correct coding region. There were 952 Arabidopsis proteins predicted to have five to nine TM regions. Based on the distribution of predicted TM numbers obtained from the entire GPCRDB entries, this range (5 to 9 TM regions) could cover almost all of the 7TMR candidates (99.1%; see Figure 1 and Materials and methods). The 22 previously predicted Arabidopsis 7TMpRs were predicted to have seven to ten TM regions (Figure 1). If we extend the range to 5 to 10 TM regions, the number of Arabidopsis 7TMpR candidates becomes 1,179 proteins.

Figure 1

Distribution of transmembrane numbers predicted by HMMTOP (black bars) and TMHMM (gray bars) from the 500 7TMR sample sequences

Distribution of transmembrane numbers predicted by HMMTOP (black bars) and TMHMM (gray bars) from the 500 7TMR sample sequences. Proportions (%) of the proteins predicted to have six to eight and five to nine TM regions by HMMTOP are shown at the top. The percentages shown in parentheses were obtained from the entire 7,674 7TMR dataset in GPCRDB. The numbers shown on the top of black bars are the number of previously predicted 22 Arabidopsis 7TMpR proteins.

Choosing 7TMpR candidates by combining prediction results

Among the ten alignment-free classifiers, LOG misclassified seven previously predicted Arabidopsis 7TMpRs. KNN with K set at 5, 10, and 15 missed one, while KNN with K set at 20 classified them all correctly (see Materials and methods on KNN). To reduce the number of false positives (non-7TMRs predicted as 7TMRs) as well as false negatives (7TMRs predicted as non-7TMRs) and to obtain a set of 7TMpR candidates with higher confidence, we examined combinations of the prediction results by the remaining six alignment-free methods (LDA, QDA, KNN with K = 20, SVM-AA, SVM-di, and PLS-ACC). There were 652 proteins predicted as 7TMpR candidates by all six methods (by choosing the strict intersection). Using the number of predicted TM regions to be 5 to 10, 394 (342 after removing duplicated entries due to alternative splicing) proteins were identified as 7TMR candidates. These Arabidopsis proteins are listed in Additional data file 1. Of the 22 previously predicted 7TMpRs, 20 were found in this list. Although HHP4 and HHP5 were not included in this list, both were identified by two of the alignment-free methods: KNN and SVM-AA. Note that RGS1 and five HHP (as well as nine MLO and GCR1) sequences were excluded from the training set, and these six were not identified as candidate 7TMpRs by SAM.

A further restriction to protein topology of exactly 7TM regions and an amino-terminus located extracellularly reduced the candidate number to 64 (54 excluding duplications due to alternative splicing). This set included nine of the 22 previously predicted 7TMpRs. These 54 7TMpR candidates are the first targets for our further analysis and are summarized in Table 2 (also listed in Additional data file 2). Eighteen are described as simply 'expressed proteins' in the TAIR database (except for AT3G26090, which encodes RGS1). Interestingly, one of them (AT5G27210) is known to have weak similarity to a mouse orphan 7TMR. While others are known to belong to certain protein families (for example, MtN3 family), in many cases, their molecular functions have not been identified, and further investigation on these 7TMpR candidates is warranted.

Table 2

Summary of the 54 7TMpR candidates identified in this study¹

Groups*

TAIR locus IDs

Multiple members from gene families

Nodulin MtN3 family proteins (8/17)

At1g21460, At3g16690, At3g28007, At3g48740, At4g25010, At5g13170, At5g23660, At5g50800

MLO proteins (7/15)

At1g11000 (MLO4), At1g26700 (MLO14), At1g42560 (MLO9), At2g33670 (MLO5), At2g44110 (MLO15), At4g24250 (MLO13), At5g53760 (MLO11)

Expressed protein family 1 (2/6)

At1g77220, At4g21570

GNS1/SUR4 membrane family proteins (3/4)

At1g75000, At3g06470, At4g36830

Perl1-like family protein (2/2)

At1g16560, At5g62130

TOM3 family proteins (3/3)

At1g14530, At2g02180, At4g21790

Expressed protein family 2 (3/5)

At1g10660, At2g47115, At5g62960

Expressed protein family 3 (2/4)

At3g09570, At5g42090

Expressed protein family 4 (2/5)

At1g49470, At5g19870

Expressed protein family 5 (2/5)

At3g63310, At4g02690

Single copy genes (8)

At1g48270 (GCR1), At1g57680, At2g41610, At2g31440, At3g04970, At3g26090 (RGS1), At3g59090, At4g20310

Single member from small gene families (8)

At2g01070, At3g19260, At2g35710, At2g16970, At1g15620, At1g63110, At4g36850, At5g27210

Single member from big gene families (4)

At1g71960, At3g01550, At5g23990, At5g37310

*The number of candidates identified in this study belonging to each group is shown in parentheses (the number of all proteins in each group is given after '/'). More detailed information is given in Additional data file 2.

The 54 proteins were grouped into families based on similarities to known protein sequences. Eight of the 54 7TMpR candidates, including GCR1 and RGS1, are encoded by single copy genes. In addition to the seven MLO proteins identified, there are eight MtN3 family members, two proteins of an unnamed family consisting of six expressed proteins, as well as multiple (two to three) members from smaller gene families (five or less). All members of the TOM3 family and the Perl1-like family, as well as the majority of the GNS/SUR4 family and an unnamed family consisting of five expressed proteins (expressed protein family 2) were included in the list. The identification of multiple members from these gene families using our alignment-free methods supported the consistency of this approach. However, for most of these families, not all members were found. Additionally, eight single representatives of small protein families consisting of two to five members and four single representatives of large protein families were found in the list. Some of these proteins, especially those from large protein families, may represent false positives as 7TMpR candidates. This 7TMR mining method can be refined, for example, by re-training models as well as using more flexible hierarchical classification.

The five predicted heptahelical proteins (HHP1-5) reported by Hsieh and Goodman 13 were identified by sequence similarity to human adiponectin receptors (AdipoRs) and membrane progestin receptors (mPRs) that share little sequence similarity to known GPCRs. HHP1-3 were identified in our initial list of 394 but were culled from the final list of 54 Arabidopsis 7TMpR candidates. This is because HMMTOP predicted HHP1, HHP2, HHP4, and HHP5 to have seven TM regions and intracellular amino termini, in contrast to known GPCRs. This unusual structural topology was also found in AdipoRs 1348. HHP3 had eight predicted TM regions. Of the 15 MLO proteins, 8 were also predicted to have 8 to 10 TM regions by HMMTOP (Figure 1). Recently, Benton et al. 49 experimentally showed that Drosophila odorant receptors, another extremely diverged 7TMR family, have intracellular amino termini. Among our 394 candidate list, 23 proteins were predicted to have seven TM regions and intracellular amino termini (Additional data file 1). Therefore, we consider these 54 as a minimum working set of 7TMpR candidates, and many of the other proteins included in the list of 394 should be examined in the second stage.

Expression patterns of genes encoding the 7TMpR candidates and G-protein subunits

We utilized the Meta-Analyzer server of the Genevestigator web site to study spatial expression patterns of Arabidopsis genes encoding the 7TMpR candidates and G-protein subunits. Note that the expression of MLO genes were not included in this analysis since we reported them recently 50. As is shown in Figure 2, expression patterns of analyzed 7TMpR candidates can be divided into two major groups; about half of them show distinct tissue specificity, whereas the other half either exhibit less distinct expression patterns or display ubiquitous expression. All genes encoding G-protein subunits fall into the latter major group. Ubiquitous expression of genes encoding G-protein subunits allows overlap with genes in both groups, and makes, in principle, co-functioning of G-proteins with these 7TMpR candidates spatially and temporally possible. All eight genes encoding the MtN3 family proteins appear to have distinct tissue specific expression. Among them, At3g48740 and At4g25010 have the highest sequence similarities to At5g23660 and At5g50800, respectively. Both pairs of genes share similar or overlapping expression patterns, suggesting relatedness/similarity of their functions. Confirming the actual functions of the 7TMpR candidates as GPCRs requires further extensive testing. A possible involvement of these candidate proteins in 'G protein-independent' signaling mechanisms also needs to be explored.

Figure 2

Expression patterns of Arabidopsis genes encoding 7TMpR candidates and G-protein subunits among tissues

Expression patterns of Arabidopsis genes encoding 7TMpR candidates and G-protein subunits among tissues. The figure was modified from an output of the Meta-Analyzer of Genevestigator (last updated in November 2005), which illustrates expression levels of each gene in different organs. Relative expression levels of a gene in different organs/tissues are given as heat maps in blue-scale coding that reflects absolute signal values, where darker colors represent stronger expression. All gene-level profiles are normalized for coloring such that, for each gene, the highest signal intensity obtains a value 100% (shown in the darkest blue and marked with an asterisk) and absence of signal obtains a value 0% (shown in white). All GeneChip data was processed using Affymetrix MAS5.0. Special precaution is required for gene expression in certain cell types (for example, pollen), since difference in normalization may achieve different results. Probe-sets of five 7TMpR candidates (At1g15620. At1g75000, At4g21570, At4g36850, and At5g23990) were not present in the 22K chip, and, therefore, their tissue-specific expression could not be assessed. For At2g35710, two probe-sets (265797_at^aand 265841_at^b) were designed on the chip. Gene names for those belonging to the MtN3 family are shown in boldface and marked with an asterisk. Genes encoding G-protein subunits (AGB1, GPA1, AGG1, and AGG2) as well as two reported 7TMpRs (RGS1 and GCR1) are labeled accordingly in boldface.

Conclusion

We show that the profile HMM protein classification method, currently one of the most used, is overly specific (conservative) when applied to extremely diverged 7TMpR proteins. Our premise is that there are more 7TMpRs yet to be identified in the A. thaliana and other genomes divergent to humans. The limitations were that the lack of available samples limits the effectiveness of profile HMM methods, and while alignment-free methods are more sensitive, they have high rates for false positives. The candidate 7TMpR proteins provided in this study, for example, can be included to expand the training set and re-iteration using refined training sets can be done to reduce false positive rates. However, this is possible only after these new candidates are confirmed as true positives experimentally.

The strategy we described here overcomes the 'chicken-or-egg' problem; predictions by multiple protein classification methods and the number of predicted transmembrane regions were used to identify a more likely reduced set of 7TMR candidates. By setting up various methods as hierarchical multiple filters, one can prioritize target protein sets for further experimental confirmation of their functions.

Materials and methods

Arabidopsis protein data

We downloaded 28,952 protein sequences from TIGR (Arabidopsis thaliana database release 5, dated 10 June 2004) 51. Among the 28,952 proteins, 2,760 are derived from alternative splicing.

Training data preparation for protein classification

Positive training samples (known 7TMR sequences) were obtained from GPCRDB (Information System for G Protein-Coupled Receptors, Release 9.0, last updated on 28 June 28 2005) 67. In the GPCRDB, 2,030 7TMRs (originally collected from the Swiss-Prot protein database) were grouped into six major classes (classes A to E plus the Frizzled/Smoothened family) and six putative families (ocular albinism proteins, insect odorant receptors, plant MLO receptors, nematode chemoreceptors, vomeronasal receptors, and taste receptors). Five hundred 7TMR sequences were randomly sampled and used as the positive samples. Note that 'putative/unclassified' (orphan) 7TMRs and bacteriorhodopsins were not included in this dataset. These 500 7TMRs included six of the15 known Arabidopsis MLO proteins. Among the 22 currently known Arabidopsis 7TMpRs, in addition to the nine MLO proteins, GCR1 as well as six recently identified Arabidopsis 7TMpRs (AtRGS1 and HHP1-5; GPCRDB does not list these proteins) were not included in the random 500 7TMR samples. Note that the 15 Arabidopsis 7TMpRs not included in the training set can be used to assess the classifier performance as test cases.

For negative samples, 500 non-7TMR sequences longer than 100 amino acids were randomly sampled from the Swiss-Prot section of the UniProt Knowledgebase 5253. The average length of the 500 non-7TMR sequences was 401 amino acids (with a maximum length of 2,512 amino acids). Positive and negative samples were combined to create a training dataset. Note that only positive samples were used to train the profile HMM classifier, SAM (see below).

Protein classification methods used

One alignment-based method (profile HMM) and four types of alignment-free multivariate methods were included in our analysis.

Profile hidden Markov models

Profile HMMs are full probabilistic representations of sequence profiles 27. Sample sequences need to be alignable, and thus only positive samples can be used for training. Two programs in SAM (version 3.4) 5455 were used: buildmodel to build profile HMMs with the nine-component Dirichlet mixture priors 56, and hmmscore to calculate scores and e-values. The 'calibration' option (for more accurate e-value calculation) and the fully local scoring option (-sw 2) were used. The e-value threshold was set at 0.01 for choosing 7TMR candidates.

Discriminant function analysis

Moriyama and Kim 38 described the three parametric (linear, quadratic, logistic) and nonparametric K-nearest neighbor methods that were shown to perform better than the profile HMM method. Therefore, we included these four alignment-free methods (LDA, QDA, LOG, and KNN) in our analysis. For KNN, K was set at 5, 10, 15, or 20, where K is the number of neighbors. The four variables used (amino acid index and three periodicity statistics) were described in Kim et al. 37. S-PLUS statistical package (Insightful Corporation, Seattle, WA, USA, version 6.1.2 for Linux) with the MASS module 57 was used for the classifier development.

Support vector machines with amino acid composition

SVMs are learning machines that make binary classifications based on a hyperplane separating a remapped instance space 58. A kernel function can be chosen so that the remapped instances on a multidimensional feature space are linearly separable. The radial basis kernel, exp(-γ||x - y||²), was used in this study. The parameter γ was set to 102 based on the median of Euclidean distances between positive examples and the nearest negative example as described in Jaakkola et al. 59. Simple 19 amino acid frequencies (the 20th amino acid frequency can be explained completely by the other 19) of each protein sequence were used as an input vector for SVMs. Programs svm_learn and svm_classify of the SVM^lightpackage version 5.0 60 were used for training and classification, respectively, by SVM. The default value of the regulatory parameter C (0.5006) was used with svm-learn. Our comparative analysis showed that SVM-AA performs better than profile HMMs when they are applied to remote similarity identification, the same problem we deal with in this study (PK Strope and EN Moriyama, submitted).

Support vector machines with dipeptide composition

We also included an SVM classifier with dipeptide composition 4041. The SVM^lightpackage version 5.0 60 was used for training and classification as before. The regulatory parameter C = 1 and the radial basis kernel function parameter γ = 90 were chosen by the grid analysis using 5-fold cross-validation.

Partial least squares with amino acid properties

PLS regression is a projection method that takes into account correlations between independent and dependent variables 61. We used the pls.pcr package, an R implementation developed by Wehrens and Mevik 6263, with the SIMPLS method, four latent variables, and cross-validation options. Each amino acid in the protein sequences was first converted to a set of 5 principal component scores developed from 12 physico-chemical properties. The auto/cross covariance (ACC) method developed by Wold et al. 64 was then applied to each of the converted sequences. ACC describes the average correlations between two residues a certain lag (amino acids) apart. The lag size of 30 was chosen for optimal classification performance. We found that the performance of PLS-ACC is robust even when only a small number of positive samples (5 or 10) are available for training. In contrast, the performance of profile HMMs suffered extremely when positive sample size was small. The 12 physico-chemical properties used and more details on the use of PLS in protein classification are described elsewhere 43. The cutoff value of 0.4999 was used for choosing 7TMR candidates in this study, which was determined as the average of the minimum error points 39 obtained from 500 replications of 10-fold cross-validation analysis using the training dataset.

Transmembrane region prediction

HMMTOP 2.0 4647 and TMHMM (originally as in 65 but implemented as S-TMHMM by 66) were used for predicting transmembrane regions. Figure 1 shows the numbers of TM regions predicted by the two methods for the 500 7TMR sequences used for classifier training. HMMTOP predicted 7TM regions from 433 7TMRs (86.6%), while only 165 7TMRs (33%) were predicted to have 7TM regions by TMHMM. HMMTOP predicted 97% or more of 7TMRs to have 6 to 8 TM regions, and with 5 to 9 TM regions more than 99% of 7TMRs were included. Using TMHMM, in order to include 97% of 7TMRs, the range of predicted TM numbers needs to be between 4 and 10. Therefore, we decided to use HMMTOP in our further analysis. With HMMTOP using the range of five to nine TM regions, we should be able to cover almost all possible 7TM proteins.

Grouping of the candidate proteins

The candidate proteins were grouped based on the e-values obtained by BLASTP protein similarity search 6768 against the Arabidopsis protein database using the default parameter set (for example, BLOSUM62) at the TAIR web site 45. The e-value threshold of 10^-20was used to identify protein families similar to the candidate proteins.

Expression patterns of genes encoding 7TMR candidates and G-protein subunits

Expression patterns of genes encoding 7TMpR candidates and G-protein subunits among tissues was studied by using the Meta-Analyzer server of the Genevestigator web site (last updated in November 2005) 6970. All data were generated using the 22K Affymetrix ATH1 Arabidopsis Genome array. Gene expression profiles based on microarray data were clustered according to similarity in expression patterns. Hierarchical clustering results were generated by default settings using pairwise Euclidean distances and the average linkage method.

Additional data files

The following additional data files are available with the online version of this paper. Additional data file 1 is the list of the 394 Arabidopsis thaliana 7TMpR candidates. Additional data file 2 lists the 54 7TMpR candidates identified in this study. These 7TMpR candidates were grouped based on their similarities with known protein families. HTML versions of the candidate lists with TAIR links and other supplementary data are available at 71.

Additional data file 1

The 394 A. thaliana 7TMpR candidates

The 394 A. thaliana 7TMpR candidates.

Click here for file

Additional data file 2

The 54 7TMpR candidates identified in this study

These 7TMpR candidates were grouped based on their similarities with known protein families. HTML versions of the candidate lists with TAIR links and other supplementary data are available at [72].

Click here for file

Acknowledgements

This work was partly funded by Nebraska EPSCoR Women in Science and NSF EPSCoR Type II grants (to ENM); Bioinformatics Interdisciplinary Research Scholars sponsored by NSF EPSCoR Infrastructure Improvement grant: Bioinformatics Research Laboratory (to PKS and SOO); and grants from the NIGMS (GM65989-01), the DOE (DE-FG02-05er15671), and the NSF (MCB-0209711) (to AMJ).

Seven-transmembrane receptors. Pierce KL Premont RT Lefkowitz RJ Nat Rev Mol Cell Biol 2002 3 639 650 10.1038/nrm908 12209124 GCR1 can act independently of heterotrimeric G-protein in response to brassinosteroids and gibberellins in <it>Arabidopsis</it> seed germination. Chen JG Pandey S Huang J Alonso JM Ecker JR Assmann SM Jones AM Plant Physiol 2004 135 907 915 514125 15181210 10.1104/pp.104.038992 The signal to move: <it>D. discoideum </it>go orienteering. Kimmel AR Parent CA Science 2003 300 1525 1527 10.1126/science.1085439 12791977 Transduction of receptor signals by beta-arrestins. Lefkowitz RJ Shenoy SK Science 2005 308 512 517 10.1126/science.1109237 15845844 Molecular mechanisms of ligand binding, signaling, and regulation within the superfamily of G-protein-coupled receptors: molecular modeling and mutagenesis approaches to receptor structure and function. Kristiansen K Pharmacol Ther 2004 103 21 80 10.1016/j.pharmthera.2004.05.002 15251227 GPCRDB information system for G protein-coupled receptors. Horn F Bettler E Oliveira L Campagne F Cohen FE Vriend G Nucleic Acids Res 2003 31 294 297 165550 12520006 10.1093/nar/gkg103 GPCRDB: Information System for G Protein-coupled Receptors http://www.gpcr.org/7tm/ Neurobiology of the <it>Caenorhabditis elegans </it>genome. Bargmann CI Science 1998 282 2028 2033 10.1126/science.282.5396.2028 9851919 Topology, subcellular localization, and sequence diversity of the Mlo family in plants. Devoto A Piffanelli P Nilsson I Wallin E Panstruga R von Heijne G Schulze-Lefert P J Biol Chem 1999 274 34993 35004 10.1074/jbc.274.49.34993 10574976 Molecular phylogeny and evolution of the plant-specific seven-transmembrane MLO family. Devoto A Hartmann HA Piffanelli P Elliott C Simmons C Taramino G Goh CS Cohen FE Emerson BC Schulze-Lefert P J Mol Evol 2003 56 77 88 10.1007/s00239-002-2382-5 12569425 Cloning of a putative G-protein-coupled receptor from <it>Arabidopsis thaliana</it>. Josefsson LG Rask L Eur J Biochem 1997 249 415 420 10.1111/j.1432-1033.1997.t01-1-00415.x 9370348 The <it>Arabidopsis </it>putative G protein-coupled receptor GCR1 interacts with the G protein alpha subunit GPA1 and regulates abscisic acid signaling. Pandey S Assmann SM Plant Cell 2004 16 1616 1632 490050 15155892 10.1105/tpc.020321 A novel gene family in <it>Arabidopsis </it>encoding putative heptahelical transmembrane proteins homologous to human adiponectin receptors and progestin receptors. Hsieh M-H Goodman HM J Exp Bot 2005 56 3137 3147 10.1093/jxb/eri311 16263907 A seven-transmembrane RGS protein that modulates plant cell proliferation. Chen J-G Willard FS Huang J Liang J Chasse SA Jones AM Siderovski DP Science 2003 301 1728 1731 10.1126/science.1087790 14500984 Role of a heterotrimeric G protein in regulation of <it>Arabidopsis </it>seed germination. Ullah H Chen JG Wang S Jones AM Plant Physiol 2002 129 897 907 161710 12068128 10.1104/pp.005017 Evidence for kinship between diverse G-protein coupled receptors. Josefsson LG Gene 1999 239 333 340 10.1016/S0378-1119(99)00392-3 10548735 Plants: the latest model system for G-protein research. Jones AM Assmann SM Embo Rep 2004 5 572 578 1299082 15170476 10.1038/sj.embor.7400174 Occurrence in <it>Saccharomyces cerevisiae </it>of a gene homologous to the cDNA coding for the alpha subunit of mammalian G proteins. Nakafuku M Itoh H Nakamura S Kaziro Y Proc Natl Acad Sci USA 1987 84 2140 2144 304604 3031665 10.1073/pnas.84.8.2140 Isolation of a second yeast <it>Saccharomyces cerevisiae </it>gene (GPA2) coding for guanine nucleotide-binding regulatory protein: studies on its structure and possible functions. Nakafuku M Obara T Kaibuchi K Miyajima I Miyajima A Itoh H Nakamura S Arai K Matsumoto K Kaziro Y Proc Natl Acad Sci USA 1988 85 1374 1378 279773 2830616 10.1073/pnas.85.5.1374 The STE4 and STE18 genes of yeast encode potential beta and gamma subunits of the mating factor receptor-coupled G protein. Whiteway M Hougan L Dignard D Thomas DY Bell L Saari GC Grant FJ O'Hara P MacKay VL Cell 1989 56 467 477 10.1016/0092-8674(89)90249-3 2536595 Overlapping functions for two G protein alpha subunits in <it>Neurospora crassa</it>. Baasiri RA Lu X Rowley PS Turner GE Borkovich KA Genetics 1997 147 137 145 1208097 9286674 Identification of a G protein alpha subunit from <it>Neurospora crassa </it>that is a member of the Gi family. Turner GE Borkovich KA J Biol Chem 1993 268 14805 14811 8325859 The genome sequence of the filamentous fungus <it>Neurospora crassa</it>. Galagan JE Calvo SE Borkovich KA Selker EU Read ND Jaffe D FitzHugh W Ma LJ Smirnov S Purcell S Nature 2003 422 859 868 10.1038/nature01554 12712197 ARAMEMNON, a novel database for <it>Arabidopsis</it> integral membrane proteins. Schwacke R Schneider A van der Graaff E Fischer K Catoni E Desimone M Frommer WB Flugge UI Kunze R Plant Physiol 2003 131 16 26 166783 12529511 10.1104/pp.011577 ARAMEMNON: Plant Membrane Protein Database http://aramemnon.botanik.uni-koeln.de Novel G-protein-coupled receptor-like proteins in the plant pathogenic fungus <it>Magnaporthe grisea</it>. Kulkarni R Thon M Pan H Dean R Genome Biol 2005 6 R24 1088943 15774025 10.1186/gb-2005-6-3-r24 Durbin R Eddy S Krogh A Mitchison G Biological Sequence Analysis: Probabilistic Models of Proteins and Nucleic Acids Cambridge: Cambridge University Press 1998 The Pfam protein families database. Bateman A Coin L Durbin R Finn RD Hollich V Griffiths-Jones S Khanna A Marshall M Moxon S Sonnhammer EL Nucleic Acids Res 2004 32 D138 141 308855 14681378 10.1093/nar/gkh121 Pfam: Database of Protein Families and HMMs http://pfam.janelia.org/ SMART 4.0: towards genomic data integration. Letunic I Copley RR Schmidt S Ciccarelli FD Doerks T Schultz J Ponting CP Bork P Nucleic Acids Res 2004 32 D142 144 308822 14681379 10.1093/nar/gkh088 SMART 4.0 http://smart.embl.de/ Assignment of homology to genome sequences using a library of hidden Markov models that represent all proteins of known structure. Gough J Karplus K Hughey R Chothia C J Mol Biol 2001 313 903 919 10.1006/jmbi.2001.5080 11697912 The PROSITE database. Hulo N Bairoch A Bulliard V Cerutti L De Castro E Langendijk-Genevaux PS Pagni M Sigrist CJ Nucleic Acids Res 2006 34 D227 230 1347426 16381852 10.1093/nar/gkj063 PROSITE: Database of Protein Families and Domains http://www.expasy.org/prosite/ PRINTS and its automatic supplement, prePRINTS. Attwood TK Bradley P Flower DR Gaulton A Maudling N Mitchell AL Moulton G Nordle A Paine K Taylor P Nucleic Acids Res 2003 31 400 402 165477 12520033 10.1093/nar/gkg030 PRINTS http://umber.sbs.man.ac.uk/dbbrowser/PRINTS/ Identification of novel multi-transmembrane proteins from genomic databases using quasi-periodic structural properties. Kim J Moriyama EN Warr CG Clyne PJ Carlson JR Bioinformatics 2000 16 767 775 10.1093/bioinformatics/16.9.767 11108699 Protein family classification with discriminant function analysis. Moriyama EN Kim J Genome Exploitation: Data Mining the Genome New York: Springer Gustafson JP, Shoemaker R, Snape JW 2005 121 132 Classifying G-protein coupled receptors with support vector machines. Karchin R Karplus K Haussler D Bioinformatics 2002 18 147 159 10.1093/bioinformatics/18.1.147 11836223 GPCRpred: an SVM-based method for prediction of families and subfamilies of G-protein coupled receptors. Bhasin M Raghava GP Nucleic Acids Res 2004 32 W383 389 441554 15215416 10.1093/nar/gkh001 GPCRpred http://www.imtech.res.in/raghava/gpcrpred/ Classification of G-protein coupled receptors by alignment-independent extraction of principal chemical properties of primary amino acid sequences. Lapinsh M Gutcaits A Prusis P Post C Lundstedt T Wikberg JES Protein Sci 2002 11 795 805 10.1110/ps.2500102 11910023 Protein family classification with partial least squares. Opiyo SO Moriyama EN J Proteome Res The <it>Arabidopsis </it>Information Resource (TAIR): a model organism database providing a centralized, curated gateway to Arabidopsis biology, research materials and community. Rhee SY Beavis W Berardini TZ Chen G Dixon D Doyle A Garcia-Hernandez M Huala E Lander G Montoya M Nucleic Acids Res 2003 31 224 228 165523 12519987 10.1093/nar/gkg076 The <it>Arabidopsis </it>Information Resource http://www.arabidopsis.org The HMMTOP transmembrane topology prediction server. Tusnady GE Simon I Bioinformatics 2001 17 849 850 10.1093/bioinformatics/17.9.849 11590105 HMMTOP http://www.enzim.hu/hmmtop Cloning of adiponectin receptors that mediate antidiabetic metabolic effects. Yamauchi T Kamon J Ito Y Tsuchida A Yokomizo T Kita S Sugiyama T Miyagishi M Hara K Tsunoda M Nature 2003 423 762 769 10.1038/nature01705 12802337 Atypical membrane topology and heteromeric function of <it>Drosophila</it> odorant receptors <it>in vivo</it>. Benton R Sachse S Michnick SW Vosshall LB PLoS Biol 2006 4 e20 1334387 16402857 10.1371/journal.pbio.0040020 Expression analysis of the AtMLO gene family encoding plant-specific seven-transmembrane domain proteins. Chen Z Hartmann HA Wu MJ Friedman EJ Chen JG Pulley M Schulze-Lefert P Panstruga R Jones AM Plant Mol Biol 2006 60 583 597 10.1007/s11103-005-5082-x 16525893 The Institute for Genomic Research (TIGR) <it>Arabidopsis thaliana</it> Database ftp site ftp://ftp.tigr.org/pub/data/a_thaliana/ath1/SEQUENCES/ATH1.pep.gz The Universal Protein Resource (UniProt). Bairoch A Apweiler R Wu CH Barker WC Boeckmann B Ferro S Gasteiger E Huang H Lopez R Magrane M Nucleic Acids Res 2005 33 D154 159 540024 15608167 10.1093/nar/gki070 UniProt: The Universal Protein Resource http://www.uniprot.org Hidden Markov models for sequence analysis: Extension and analysis of the basic method. Hughey R Krogh A Comput Appl Biosci 1996 12 95 107 8744772 SAM: Sequence Alignment and Modeling System http://www.cse.ucsc.edu/research/compbio/sam.html Dirichlet mixtures: a method for improving detection of weak but significant protein sequence homology. Sjolander K Karplus K Brown M Hughey R Krogh A Mian IS Haussler D Comput Appl Biosci 1996 12 327 345 8902360 S-plus MASS module http://www.stats.ox.ac.uk/pub/MASS4/ Vapnik VN The Nature of Statistical Learning Theory New York: Springer-Verlag 2 1999 A discriminative framework for detecting remote protein homologies. Jaakkola T Diekhans M Haussler D J Comput Biol 2000 7 95 114 10.1089/10665270050081405 10890390 Making large-Scale SVM learning practical. Joachims T Advances in Kernel Methods - Support Vector Learning Cambridge: MIT Press Schölkopf B, Burges C, Smola A 1999 169 184 Partial least squares regression: A tutorial. Geladi P Kowalski BR Anal Chim Acta 1986 185 1 17 10.1016/0003-2670(86)80028-9 R Development Core Team R: A Language and Environment for Statistical Computing Vienna, Austria: R Foundation for Statistical Computing 2005 pls: Partial Least Squares Regression (PLSR) and Principal Component Regression (PCR): R package version 1.2-1. http://mevik.net/work/software/pls.html DNA and peptide sequences and chemical processes multivariately modeled by principal component analysis and partial least-squares projections to latent structures. Wold S Jonsson J Sjostrom M Sandberg M Rannar S Anal Chim Acta 1993 277 239 253 10.1016/0003-2670(93)80437-P A hidden Markov model for predicting transmembrane helices in protein sequences. Sonnhammer EL von Heijne G Krogh A Proc Int Conf Intell Syst Mol Biol 1998 6 175 182 9783223 Best alpha-helical transmembrane protein topology predictions are achieved using hidden Markov models and evolutionary information. Viklund H Elofsson A Protein Sci 2004 13 1908 1917 10.1110/ps.04625404 15215532 Basic local alignment search tool. Altschul SF Gish W Miller W Myers EW Lipman DJ J Mol Biol 1990 215 403 410 10.1006/jmbi.1990.9999 2231712 BLAST http://www.ncbi.nlm.nih.gov/BLAST/ GENEVESTIGATOR. Arabidopsis microarray database and analysis toolbox. Zimmermann P Hirsch-Hoffmann M Hennig L Gruissem W Plant Physiol 2004 136 2621 2632 523327 15375207 10.1104/pp.104.046367 Genevestigator: Arabidopsis Microarray Database and Analysis Toolbox https://www.genevestigator.ethz.ch <it>Arabidopsis thaliana </it>7TMR Mining http://bioinfolab.unl.edu/emlab/at7tmr/index.html