How accurately can we discriminate G-protein-coupled receptors as 7-tms TM protein sequences from other sequences?

The group of 2502 transmembrane (TM) protein sequences with seven TM segments (7-tms) registered in SWISS-PROT 46.0 contains 2200 G-protein-coupled receptors (GPCRs), indicating that GPCR candidates can be detected with a reliability of 87.9% in the eukaryotic genomes merely by correctly predicting the number of TM segments as 7-tms. The predictive accuracies of TM topology-prediction methods proposed so far are not as high as expected; even the best method, HMMTOP 2.0, can only achieve a capture rate of 7-tms sequences of 77.6%. It is necessary to improve this performance as much as possible, even if by only a few percentage points, in order to identify as many novel GPCR candidate genes as possible among the increasing number of newly sequenced genomes. In this study, we propose a simple but useful prediction method for detecting as many 7-tms TM protein sequences as GPCR candidates in eukaryotic genomes as possible. This is achieved by employing a two-step prediction procedure. The first step involves collecting 7-tms sequences by the best prediction method (HMMTOP 2.0), and the second involves picking up the remaining 7-tms sequences by the second-best method (TMHMM 2.0). By this procedure, the capture rate of 7-tms TM protein sequences in SWISS-PROT can be improved considerably from 77.6% to 84.5%, and the number of GPCR candidate sequences predicted as 7-tms in the human genome (Build 35) is increased from 790 (by HMMTOP 2.0) to 903. These 790 and 903 candidate sequences include, respectively, 587 and 636 of the known human GPCRs of the 717 registered in SWISS-PROT 46.0, demonstrating that the proposed combinatorial method is effective in detecting GPCR candidate genes in eukaryotic genomes.

Proteome-wide classification and identification of mammalian-type GPCRs by binary topology pattern.

G protein-coupled receptors (GPCRs), a large eukaryotic protein family, have proved difficult to comprehensively detect and functionally identify by homology searches and domain detection, because they are highly divergent and their sequences share strikingly little similarity. Transmembrane (TM) topology pattern analysis has been used to classify TM proteins, and such patterns are conserved within GPCRs of similar function. Here, we developed a stepwise binary topology pattern (BTP) method for GPCR classification and identification and used it to identify and classify mammalian-type GPCRs in the genomes of 10 different eukaryotic species. A binary topology pattern was obtained for each functional class or group by assigning binary loop threshold lengths of "0" (short loop) or "1" (long loop). The GPCR-classification ability of the BTP method had quite high accuracies for classifying GPCR functions at the class level (Classes A, B, C, Frizzled/Smoothened, Non-GPCR, based on the GPCRDB classification scheme), with many classes being classified with 100% accuracy. Sufficiently high accuracies were also maintained at the functional group level, 0.945 over 15 functional groups. Proteome-wide mammalian-type GPCR searches in 10 eukaryotic genomes (H. sapiens, M. musculus, F. rubripes, C. intestinalis, A. thaliana, D. melanogaster, A. gambiae, C. elegans, P. falciparum, S. cerevisiae) using the BTP method showed much higher classification/identification in non-mammalian genomes than typical BLAST searches, in which a higher number of sequences were classified as Non-GPCR. This stepwise BTP method should prove useful for the identification and functional classification of GPCRs from the genomes of a wide range of species.

Genomic analysis of immunity in a Urochordate and the emergence of the vertebrate immune system: "waiting for Godot".

Genome-wide sequence analysis in the invertebrate chordate, Ciona intestinalis, has provided a comprehensive picture of immune-related genes in an organism that occupies a key phylogenetic position in vertebrate evolution. The pivotal genes for adaptive immunity, such as the major histocompatibility complex (MHC) class I and II genes, T-cell receptors, or dimeric immunoglobulin molecules, have not been identified in the Ciona genome. Many genes involved in innate immunity have been identified, including complement components, Toll-like receptors, and the genes involved in intracellular signal transduction of immune responses, and show both expansion and unexpected diversity in comparison with the vertebrates. In addition, a number of genes were identified which predicted integral membrane proteins with extracellular C-type lectin or immunoglobulin domains and intracellular immunoreceptor tyrosine-based inhibitory motifs (ITIMs) and immunoreceptor tyrosine-based activation motifs (ITAMs) (plus their associated signal transduction molecules), suggesting that activating and inhibitory receptors have an MHC-independent function and an early evolutionary origin. A crucial component of vertebrate adaptive immunity is somatic diversification, and the recombination activating genes (RAG) and activation-induced cytidine deaminase (AID) genes responsible for the Generation of diversity are not present in Ciona. However, there are key V regions, the essential feature of an immunoglobulin superfamily VC1-like core, and possible proto-MHC regions scattered throughout the genome waiting for Godot.