LRRCE: a leucine-rich repeat cysteine capping motif unique to the chordate lineage
© Park et al; licensee BioMed Central Ltd. 2008
Received: 09 July 2008
Accepted: 12 December 2008
Published: 12 December 2008
The small leucine-rich repeat proteins and proteoglycans (SLRPs) form an important family of regulatory molecules that participate in many essential functions. They typically control the correct assembly of collagen fibrils, regulate mineral deposition in bone, and modulate the activity of potent cellular growth factors through many signalling cascades. SLRPs belong to the group of extracellular leucine-rich repeat proteins that are flanked at both ends by disulphide-bonded caps that protect the hydrophobic core of the terminal repeats. A capping motif specific to SLRPs has been recently described in the crystal structures of the core proteins of decorin and biglycan. This motif, designated as LRRCE, differs in both sequence and structure from other, more widespread leucine-rich capping motifs. To investigate if the LRRCE motif is a common structural feature found in other leucine-rich repeat proteins, we have defined characteristic sequence patterns and used them in genome-wide searches.
The LRRCE motif is a structural element exclusive to the main group of SLRPs. It appears to have evolved during early chordate evolution and is not found in protein sequences from non-chordate genomes. Our search has expanded the family of SLRPs to include new predicted protein sequences, mainly in fishes but with intriguing putative orthologs in mammals. The chromosomal locations of the newly predicted SLRP genes would support the large-scale genome or gene duplications that are thought to have occurred during vertebrate evolution. From this expanded list we describe a new class of SLRP sequences that could be representative of an ancestral SLRP gene.
Given its exclusivity the LRRCE motif is a useful annotation tool for the identification and classification of new SLRP sequences in genome databases. The expanded list of members of the SLRP family offers interesting insights into early vertebrate evolution and suggests an early chordate evolutionary origin for the LRRCE capping motif.
The leucine-rich repeat (LRR) is a widespread structural motif of 20–30 amino acids easily identifiable at the primary structure level by the characteristic 11-residue hallmark sequence Lxx Lx Lxx Nx L, where x means any amino acid and the consensus Leu and Asn positions are often substituted by other hydrophobic residues such as Ile, Val, Phe, Cys, etc [1–5]. Proteins with LRR-architecture typically contain two or more LRRs in tandem and have been identified in all life forms, from viruses to eukaryotes . The continuously expanding LRR superfamily includes intracellular, extracellular and membrane-attached proteins characterized by a common modular architecture specially suited to favour protein-protein interactions [1–3, 5, 7–9]. These proteins participate in a variety of important biological functions, including among others cell adhesion and signalling, platelet aggregation, neural development, extracellular matrix assembly, bacterial pathogenicity, disease resistance and immune response [10–20]. LRR-containing proteins and domains form curved solenoid structures where each repeat is a turn of the solenoid. The concave side of the solenoid is defined by a parallel β-sheet interwoven with a variety of structures in the convex side which include α helices, 310 helices, polyproline II helices, tandems of β turns and short β strands [1–5, 9, 21]. The biological roles of LRR proteins and domains typically relate to their ability to engage in protein-protein interactions. However, some family members recognize other ligand types such as nucleic acids, lipopolysaccharides, lipopeptides, and even small compounds such as auxins [19, 22–27]. The sites for ligand recognition map preferentially but not exclusively to the concave sites of the LRR arched structures, as demonstrated by several crystal structures of LRR proteins in complex with their ligands (see  for a recent review). Recently, some LRR proteins have been shown to form highly stable dimers through their concave side [28–31] raising the possibility of alternative scenarios where LRR dimers are either the functional units or latent forms that require dissociation prior to ligand binding .
A distinct group of LRR proteins from the extracellular matrix forms the family known as small leucine-rich repeat proteins and proteoglycans (SLRPs) [10, 32–34]. These molecules are emerging as an important family or regulatory proteins with still undiscovered functions. They typically control the correct assembly of collagen fibrils, regulate mineral deposition in bone, and modulate the activity of potent cellular growth factors through signal transduction [10, 33–35]. SLRPs have in common clusters of cysteine residues flanking their LRR domains at both N- and C-termini. The crystal structures of the two most studied SLRPs, decorin and biglycan, have been recently determined [29, 31].
SLRPs have been traditionally classified into three classes (I, II and III) depending on their gene organisation, number of LRRs and spacing of cysteine residues at the amino-terminal cluster [10, 32, 33]. Other LRR molecules have been subsequently added to the family and two additional, non-canonical classes IV and V have been defined . Class IV and V SLRPs show clear differences with those of the three first classes in number of repeats and internal repeat structure [32, 34, 36–39]; their classification as SLRPs is due to functional similarity with canonical SLRPs, extracellular location, and presence of cysteine clusters flanking the LRR domain.
N-terminal capping (plant specific)
The cysteines form a disulphide knot that connects the antiparallel β strand to the first LRR. This motif is characteristic of all proteins from the SLRP family as well as secreted or membrane-bound LRR proteins. Its main structural elements appear to be maintained irrespective of the number and spacing of cysteines . A variation of the LRRNT capping motif specific to plants has also been described (LRRNP, Table 1). Its architecture and disulphide-bonding topology differs from that of LRRNT, as observed in the crystal structure of a polygalacturonase inhibitor from Phaseolus vulgaris .
The LRR C-terminal capping motif (LRRCT) contains normally four cysteines that stabilize the local structure with two disulphide bonds [44–47]. Characteristic of this motif is an α helix that covers the hydrophobic core of the last LRR (Figure 1b). This capping motif seems to occur slightly less often than LRRNT (Table 1) and appears to be exclusive to animal proteins. Many cysteine-capped LRR proteins have been automatically annotated as having either N-terminal or C-terminal capping structures, although close inspection of their sequences shows that in many cases both capping motifs are actually present.
In early sequence analyses, Kajava and Kobe classified the disulphide-bonded C-terminal capping motifs in four different subfamilies, named CF1 to CF4 [2, 3]. Sequences from the CF1 subfamily contain four cysteine residues and are the ones typically detected by current LRRCT database descriptors in automatic sequence annotation (Table 1). The CF2 subfamily is characterized by only two cysteines and was defined for some members of the SLRP family. The CF3 and CF4 subfamilies are specific to G-protein coupled receptors and plants, respectively. Protein sequence databases do not include separate descriptors for these four subfamilies, and sequences from the CF2 subfamily are not recognized by current LRRCT descriptors. The structure of the C-terminal capping motif from the CF2 subfamily was elucidated in the crystal structures of the protein cores of decorin and biglycan [29, 31], both being representative members of the SLRP family. This capping motif is structurally quite different from LRRCT . It extends to the last two LRRs, which are connected by a single disulphide bond (Figure 1c). The second-to-last LRR appears to be longer than all the other ones, and in the crystal structures of decorin and biglycan is extended laterally from the main LRR fold [29, 31]. We have previously named this longer, extended repeat as the "ear" repeat, and we will use the term LRRCE throughout this paper to designate the ear-containing LRR C-terminal capping motif. All extracellular LRR proteins currently classified as members of the SLRP family have N-terminal capping motifs of the LRRNT type, but their C-terminal motifs show variability. While chondroadherin and nyctalopin C-terminal sequences correspond to that of a typical LRRCT motif, the SLRPs from the canonical group (which here refers to classes I, II and III plus extracellular matrix protein 2, ECM2) show the LRRCE capping motif [29, 32]. Podocan, a recent addition to the SLRP family and class V representative, does not have any C-terminal capping [32, 36].
To investigate if the LRRCE motif is a common structural feature in other LRR proteins, we have defined characteristic sequence patterns based on a set of sequences from class I, II and III SLRPs and we have used those patterns in genome-wide searches. We present in this paper the results of this analysis.
Results and discussion
A comprehensive list of LRRCE-containing sequences
Sequences of LRRCE motifs from the UniProt non-redundant set were used in similarity searches with BLAST to retrieve additional LRRCE-containing sequences (see Methods). For most SLRPs, the LRRCE motifs are largely coded by the last exon from their corresponding genes. Thus, probe sequences just encompassing different LRRCE motifs were useful for quickly locating SLRP genes in genomes at early stages of annotation or for searching new SLRPs on the genomes of invertebrates and early chordates (see below). The extended list of 280 hits with their LRRCE sequences and accession codes, including those obtained from the NCBI and ENSEMBL databases, is provided in Additional File 1. The LRRCE regular expression pattern discussed earlier is consistent with all but seven of the sequences in this extended list. Three of these exceptions are incomplete sequences due to missing genomic data, and two more are predicted sequences with only one change with respect to the LRRCE pattern. Probably the only significant exceptions were the LRRCE sequences for the chicken and lizard homologues of ECM2, which contain an additional insertion (see Additional File 1). Many of the sequences found in similarity searches are from predicted model assemblies from genomes in early stages of annotation, and therefore some of the assignments should be considered preliminary.
Structure of the LRRCE capping motif
The LRRCE motif encompasses the ear repeat, which is extended laterally, the LRR following it, and the final β-strand closing the domain (Figure 2). The regular expression pattern runs from the beginning of the ear repeat to the second cysteine residue. Sequence conservation in LRRCE motifs across the different SLRPs follows largely structural dictates, with the highly conserved positions mainly corresponding to the core hydrophobic or asparagine residues characteristic of the LRR architecture, plus the two cysteine residues that are connected by a disulphide bond (Figure 2). Several additional positions show distinct preferences for polar or charged amino acids. The corresponding residues in the bovine decorin structure participate in a network of stabilizing charge-charge and hydrogen-bonding interactions between repeats. Thus, it is likely that similar interactions will be conserved in the other LRRCE-containing proteins to impart stability to the capping motif. Residue conservation in the ear itself is comparatively poor between closely related proteins (for example, between decorin and biglycan), but higher for the same protein across species (see examples in Figures 3 and Additional File 1). This pattern of conservation suggests that the ear extension contributes to the functional specialisation of the canonical SLRPs. For most sequences, the ear extension is 11–13 residues long (from the first cysteine to the first residue of the second LRR), and only keratocan and PRELP sequences show consistently long extensions. A buried lysine residue shown to stabilize the ear conformation in the crystal structure of decorin  is conserved as lysine or arginine in all class I and III sequences, whereas in class II sequences the same position is occupied by an aromatic or leucine residue. Fibromodulin, osteoglycin and epiphycan sequences show conserved N-linked glycosylation sites in their ear extensions. Additional potential N-linked glycosylation sites appear with varying degrees of conservation on different regions of the two repeats forming the LRRCE motif.
LRRCE motifs are always C-terminal
The LRRCE motif is a genuine C-terminal motif, with no instance of an equivalent architecture in the middle of an LRR domain or protein. The C-terminus is typically 9–15 residues away from the conserved second cysteine. The exception is the group of OMD sequences, which contain an extended C-terminal tail of about 60 residues after the LRRCE motif. This tail contains a stretch of negatively charged residues that presumably shares some functional role with the negatively charged glycosaminoglycan chains attached to the N-termini of decorin or biglycan, or the polyanionic stretches seen in the N-terminal region of asporin or preceding the LRR domain in ECM2.
The LRRCE motif can be related structurally to internal, disulphide-bonded LRR pairs that occasionally occur in LRR proteins. These LRRs do not show the lateral extensions characteristic of the LRRCE motifs. These intradomain, disulphide-bonded LRR pairs are much more widely distributed than LRRCE, as they occur in different LRR protein families from bacteria to humans. One such linkage can be seen in the three-dimensional structure of the ectodomain of Toll-like receptor 3 [49, 50]. This structure is formed by a tandem of 25 LRRs capped by LRRNT and LRRCT motifs, and contains one internal disulphide-bonded LRR pair.
The LRRCT and LRRCE structural motifs appear to be unrelated
There is no obvious relation between the LRRCT and LRRCE motifs, at least from sequences available to date. A conserved feature in the LRRCT capping motif is the presence of two or more Trp residues that contribute to maintain the hydrophobic core at the C-terminal end. A common sequence at the beginning of LRRCT motifs is NPWx Cx Cx3Wx3W, where the second and third Trp residues are located on the inner side of the α-helix characteristic of the LRRCT structure (Figure 1b). Conversely, Trp residues are not particularly conserved in the LRRCE sequences (Figure 3 and Additional File 1), and the only Trp residue in the structure of bovine decorin  is exposed to the solvent and does not participate in the hydrophobic core of the C-terminal end of the LRR structure.
A new class of canonical SLRPs, with sequences similar to ECM2
The LRRCE motif is unique to chordates
Given the exclusive vertebrate lineage of all sequences obtained in the UniProt set, searches with LRRCE sequences were carried out against several invertebrate genomes including Drosophila melanogaster (fruit fly), Apis mellifera (honeybee) , Anopheles gambiae (malaria mosquito), Aedes aegypti (yellowfever mosquito) , Caenorhabditis elegans (worm) and Strongylocentrotus purpuratus (sea urchin) . Of particular interest was the genome of the sea urchin, as this organism appears to have retained some of the genes later observed only in vertebrates . All the genomes investigated have large numbers of LRR proteins, many of those with LRRNT capping motifs. However, the searches failed to produce any true match of LRR proteins containing LRRCE motifs in these invertebrate genomes. Two LRR protein sequences from mosquito (Q16VM2_AEDAE and Q17NB1_AEDAE in TrEMBL) are currently annotated (incorrectly) as putative lumicans, although they lack recognizable LRRNT caps and neither their repeat lengths nor amino acid sequences correspond to those of lumican proteins. Nonchordate LRR protein sequences are often homologous to non-SLRP proteins such as toll-like receptors or slit proteins, and C-terminal disulphide-capping often occurs through the more common LRRCT type.
LRRCE motifs in early chordates: implications for SLRP evolution
An interesting possibility is that shorter, 12-LRR SLRPs (such as these from classes I and II), originated by exon skipping from an ancestral SLRP gene with 15 LRRs similar to SLRP1 from Ciona, in a manner illustrated by the short model assembly shown in Figure 5. Such exon skipping would have occurred after duplication of this ancestral SLRP gene. The two genes could have later evolved after additional tandem and large-scale duplication events into two independent lineages, one for ECM2 and ECM2-like proteins (discussed below), and another for the class I and class II SLRPs. Class III SLRPs would have originated by further exon skipping on the class I and II ancestor, after its divergence from the ECM2 lineage ancestor.
No gene containing an LRRCE motif could be detected in the currently available release of the amphioxus genome (Branchiostoma floridae, JGI version 1.0) . Sequence similarity searches using the protein sequences of the three SLRP proteins from Ciona (Figure 5 and ) produced putative orthologues of podocan and chondroadherin (data not shown), and partial hits to many non-SLRP protein sequences with LRR architecture (several cell-surface receptors, slit-like proteins, etc). However, the searches failed to return any clear orthologues of the canonical SLRPs or any sequence containing a LRRCE motif. Queries using sequences conforming to LRRCE motifs from Ciona yielded no hits either. The apparent absence of LRRCE-containing genes in Brachiostoma is intriguing. Cephalochordates (amphioxus) have been traditionally considered to be the primitive chordates that most resemble vertebrates, but this view has been very recently contradicted by the genome sequence data that suggests that tunicates (which include the ascidians), are more closely related to vertebrates than cephalochordates . This finding raises the possibility that the first LRRCE-containing SLRP genes appeared after the divergence of cephalochordates from the rest of the tunicate-vertebrate lineage. Alternatively, the gene equivalent to SLRP1 from Ciona was already present in a common ancestor of all chordates but may have been lost in the cephalochordate lineage.
Agnathans show already an expanded set of SLRP sequences
Representative examples of class I, II and III SLRP sequences in sea lamprey.
BGL2_PETMA (class I)
KERAL_PETMA (class II)
EPYL_PETMA (class III)
The LRRCE motif is a useful annotation tool for extending the SLRP family
The LRRCE capping motif is useful in sequence annotation as it appears to be exclusive to the canonical group of SLRPs. The presence of an LRRCE motif in a newly predicted protein sequence is sufficient for its quick identification and classification as a member of the canonical group of SLRPs (and also into one of its classes). The amino acid sequences of LRRCE motifs can be used as probes in similarity searches against translated genomic databases (TBLASTN). These searches can identify the locations and partial sequences of new putative SLRP genes in genomes at different levels of completion (as shown with the lamprey example). Using these probes, it is also possible to detect exons or domains that are missing in current assembly models of SLRP genes but nevertheless present in the genome (data not shown).
Newly predicted SLRP sequences related to ECM2 in fishes and mammals
Several ECM2-like genes have also been predicted both in fishes and mammalian genomes (Figures 4 and 10; see also Additional File 1 for sequences and accession codes). A particularly interesting example is that of a protein sequence similar to ECM2 that was first predicted in the genome of zebrafish (Danio rerio, accession codes Q1LYN3, XP_690561).
Two representative examples of ECMX predicted sequences (accession codes in parentheses).
Most interestingly, reciprocal interrogation of the human genome using the LRRCE motif from the zebrafish ECMX sequence gave a hit in chromosome X (accession code XP_001714654, currently known as hypothetical protein LOC389904), upstream of the human biglycan gene (Xq28). The predicted protein sequence of this putative novel SLRP is highly similar to the ECMX sequence from zebrafish, although different alternate model assemblies have slightly different number of repeats. This hypothetical protein, which we have named here ECMX owing to its similarity in repeat structure to ECM2 and its location in the X chromosome, has orthologues predicted in the genomes of orangutan, macaque, bovine, horse, dog, opossum (Monodelphis domestica), platypus (Ornithorhynchus anatinus), frog (Xenopus tropicalis), and several fishes (zebrafish, Gasterosteus aculeatus, Oryzias latipes, Takifugu rubripes and Tetraodon nigroviridis). The predicted sequence of the horse orthologue (Table 3) maintains the same repeat structure as the zebrafish one. Currently, it is unclear if the differences in repeat structure with the predicted human sequence are significant or if some of the predictions are partially incorrect. Partial transcription evidence for ECMX in humans (EST data) has been obtained from osteoarthritic cartilage and chondrosarcoma (accession codes BQ181183, BQ447619, BQ448435, BQ772123). There is also EST evidence for ECMX in zebrafish (accession codes EB980280, CR929461), Gasterosteus aculeatus (DW649744, DW649745), and Xenopus tropicalis (CN119819, CX371080, CX409086).
A second ECM2-like predicted sequence occurs upstream of the decorin gene in the genomes of zebrafish, Gasterosteus aculeatus, Oryzias latipes and Takifugu rubripes (Figure 11). This hypothetical protein, which we name here ECMZ, has a predicted sequence that appears to be related to the SLRP1 sequences from both Ciona species (Figures 4 and 10). These ECMZ sequences would therefore be the most similar to the ancient SLRP genes (or what is left of them) in vertebrates, and have not been retained in mammals, birds or reptiles. Given the number of ECM2-like sequences and their own clustering with the SLRP1 Ciona genes, away from classes I, II and III (Figure 4), we propose a new SLRP class that includes ECM2, ECMX, ECMZ and the Ciona SLRP1 sequences. Since these sequences appear to be more closely related to the ancestral gene from which all the canonical SLRPs derived, we name this new class "class A", for ancestral SLRPs. The similarity between the ECM2-ASPN gene tandem and the predicted ECMX-BGN and ECMZ-DCN tandems (Figure 11) strongly supports the notion of SLRP evolution by tandem and large-scale gene duplications as well as tandem gene migration [57, 64]. The duplicate copies of SLRP genes exclusive to ray-finned fishes (LUM2, FMOD2, OGN2, BGN2) appear in the same chromosomes as other SLRP sequences (Figure 11). These duplicates may be survivors of the proposed, fish-specific large-scale gene or genome duplication that would have occurred after the divergence between the actinopterygian (ray-finned fishes) and sarcopterygian (coelacanth, lungfish and all land vertebrates) lineages [65–68]. Not all the zebrafish SLRP genes shown in Figure 11 have been identified in the other fish genomes, reflecting the fact that different teleost fishes have retained different sets of duplicate genes [69, 70]. Comparison of the zebrafish and fugu genomes have revealed that despite a high degree of synteny and retention of a similar number of duplicates, in a significant number of cases, different paralogues have been preserved .
Finally, the chromosomal organization of the canonical SLRP genes in mammals follows the order shown in Figure 11 for the human and bovine genomes, where the class A paralogue is followed downstream by those from classes I, IIa, IIb and III. This organization would suggest that the ECMX-BGN pair might have been initially part of the FMOD-PRELP-OPTC gene cluster that currently is located in the chromosome 1 in humans, and later migrated to the X chromosome (as postulated for the BGN gene in ). A hypothetical class IIa gene downstream of ASPN would have disappeared completely, potentially by pseudogenisation. The fish SLRP genes share in part this chromosomal organization (Figure 11), and the class A SLRP upstream of the DCN gene can still be recognized in the predicted ECMZ sequences. An exception is presented by the LUM2 (class IIa) sequence in fishes, which appears intercalated between the FMOD (IIa) and PRELP (IIb) genes but is missing in non-fish genomes. The LUM2 sequence could have originated from additional local gene duplication in fishes or could have migrated from a different SLRP gene cluster. The presence of the additional fish SLRP duplicates (BGN2, FMOD2 and OGN2) points towards a complex history of local and large-scale duplications as well as gene migration for the surviving set of fish SLRPs.
Possible biological roles for LRRCE motifs
Probably the main role of the LRRCE motifs is to stabilize the LRR structures of the SLRPs by providing a disulphide-bonded C-terminal capping structure, in a similar manner to the LRRCT capping motifs present in many sequences of extracellular or membrane-associated LRR proteins. The requirement of disulphide-bonding integrity for SLRP biological actitvity has been demonstrated for decorin  and fibromodulin . Furthermore, thermal denaturation of decorin appears to be completely reversible as long as the disulphide bonds are not reduced .
Several diseased states have been associated to mutations occurring in the LRRCE regions of some SLRPs. Two different frame shift mutations on the decorin gene due to single base pair deletions in the LRRCE coding region have been linked to congenital stromal dystrophy of the cornea [73, 74]. Both mutations are predicted to result in a truncated decorin protein missing the last 33 C-terminal amino acids. Another mutation resulting in a premature stop codon in the ear extension of keratocan (R313X) has been linked to autosomal recessive cornea plana [75, 76]. These two truncations would eliminate most of the LRRCE structure in decorin and keratocan. Three amino acid substitutions in the LRRCE motif of opticin have been linked to high myopia , probably through disruption of the local tertiary structure. In all these examples, the predicted truncation or alteration of the LRRCE structure is likely to have a detrimental effect in the stability of the entire LRR domains of these SLRPs.
A direct role of the ear extensions in ligand binding remains an attractive yet still hypothetical scenario. The expanded set of LRRCE sequences presented here shows clear conservation trends across species in the ear extensions of a given SLRP, whereas these extensions are poorly conserved between closely related SLRPs. Thus, the ear extensions could help to differentiate the roles of SLRPs belonging to the same class. In known structures of LRR domains or proteins it is not uncommon to find extended repeats where the polypeptide chain loops out from the expected path of a regular LRR to rejoin it some residues later . These extensions may have functional significance, as shown for example by the so-called β-switch and β-finger in the structure of glycoprotein Ibα [44, 45]. Future biochemical and mutagenesis analyses on the SLRP ear extensions will be necessary to elucidate any functional role of these structures in ligand recognition and binding.
The LRRCE capping motif identified in the structures of the representative SLRPs decorin and biglycan appears to be a structural element exclusive to the main group of SLRPs, which includes the previously described classes I, II and III, plus a new class of ancestral genes that includes ECM2, the SLRP1 sequences from Ciona, and other ECM2-like sequences present mainly in fishes but with intriguing orthologues in mammals, including a yet uncharacterized new SLRP in the human X chromosome. The LRRCE motif appears to have evolved during early chordate evolution and is not found in non-chordate LRR protein sequences. Such evolutionary history is probably related to the known interactions of SLRPs with fibrilllar collagens and their regulation of collagen fibrillogenesis [10, 33, 34]. Given its exclusivity to the SLRP family, the LRRCE motif is a useful annotation tool for the identification and classification of new SLRP sequences in genome sequencing efforts. Analysis of LRRCE-containing sequences of organisms located phylogenetically between critical evolutionary events will provide useful clues for understanding the history of large-scale gene and genome duplications that appear to have occurred during vertebrate evolution. The expanded list of SLRP sequences, provided here for the first time, will facilitate the analysis of residue conservation trends in functionally significant sequence motifs, and ultimately will be useful for the elucidation of the full range of biological functions of this important family of extracellular matrix molecules.
Regular expression pattern and sequence retrieval
An initial set of 58 protein sequences annotated as SLRPs from classes I (21 sequences), II (24 sequences) and III (12 sequences) plus ECM2 (1 sequence) were selected from the Swiss-Prot database (52.1 release, http://www.ebi.ac.uk/swissprot) and aligned using CINEMA . This alignment was used together with the crystal structures of decorin and biglycan (PDB codes 1XKU and 2FT3 respectively), to define a characteristic regular expression pattern using PROSITE syntax . This pattern, designated as LRRCE hereupon, was used to retrieve additional sequences from the UniProt database (Swiss-Prot release 55.1, TrEMBL release 38.1)  using the ScanProsite tool . The LRRCE regular expression pattern was refined in iterative cycles until no further sequences were obtained. The final LRRCE expression pattern, in PROSITE syntax , was: [LIV]-X(2)-[LVIYFMA]-X-[LIFM]-X(2)-[NH]-X-[ILVF]-X(2)-[VIMFLY]-X(4)-[FIMLV]-C-X(7,20)-[LYIMV]-X(2)-[ILVTMF]-X-[LVMI]-X(2)-N-X-[IVLMAFT]-X(8,9)-[FYMPVAIS]-X-C.
Similarity searches and sequence alignment
Additional sequences were obtained through sequence similarity searches (BLAST and TBLASTN) on the NCBI  and ENSEMBL  databases. Sequences of LRRCE motifs from different organisms were used to query the different databases. Given the early stages of annotation of some of the genomes, some predicted sequences were manually corrected using supporting genome and EST data. The sequences were aligned using CLUSTALW  as implemented in the Kyoto University Bioinformatics Center http://align.genome.jp.
The phylogenetic analysis shown in Figure 4 was inferred from a multiple sequence alignment of the LRRCE motifs of the 280 sequences retrieved by ScanProsite and the similarity searches described above. The complete list of sequences used in this study and their accession numbers are provided in Additional File 1. Separate phylogenetic analyses were performed for each SLRP class on a reduced set of sequences from two mammals (human and bovine or horse), two birds (chicken and quail), lizard (anole), frog (Xenopus), four teleost fishes (zebrafish, stickleback, fugu and Tetraodon), lamprey, and the two Ciona species (Figures 6, 7, 8, 9, 10). For each class, phylogenetic analyses were inferred from a gap stripped multiple alignment of the selected sequences generated using CLUSTALW  and analyzed by three different independent phylogenetic methods. Neighbour-joining (NJ) trees and bootstrap replicates were generated using SEQBOOT, PROTDIST, NEIGHBOR and CONSENSE from the PHYLIP package  using default settings. Maximum Likelihood trees were inferred using PROML from the PHYLIP package using default settings. Bayesian tree inference values were produced using the MrBayes programme , where Markov Chain Monte Carlo analysis was performed for 100,000 generations using 6 chains. Clade-credibility values indicating statistically probable clades (>60%) are indicated from the three methods on the NJ trees (Figures 6, 7, 8, 9, 10), and UniProt organism abbreviations are used for sequence identification: XENTR, Xenopus tropicalis (western clawed frog); ANOCA, Anolis carolinensis (anole lizard); CHICK, Gallus gallus (chicken); COTJA, Coturnix japonica (Japanese quail); DANRE, Danio rerio (zebrafish); GASAC, Gasterosteus aculeatus (stickleback); TETNG, Tetraodon nigroviridis (green pufferfish); FUGRU, Fugu rubripes (Japanese pufferfish); PETMA, Petromyzon marinus (sea lamprey); CIOIN, Ciona intestinalis (transparent sea squirt); CIOSA, Ciona savignyi (Pacific transparent sea squirt).
HP acknowledges a Chevening Scholarship, funded by the Foreign and Commonwealth Office of the British Government, to support her postgraduate studies. JHJ acknowledges funding by a BBSRC PhD studentship. RBH, PNB and JB are members of the Wellcome Trust Centre for Cell-Matrix Research, supported by the Wellcome Trust.
- Kobe B, Deisenhofer J: Proteins with leucine-rich repeats. Curr Opin Struct Biol. 1995, 5 (3): 409-416. 10.1016/0959-440X(95)80105-7.PubMedView ArticleGoogle Scholar
- Kajava AV: Structural diversity of leucine-rich repeat proteins. J Mol Biol. 1998, 277 (3): 519-527. 10.1006/jmbi.1998.1643.PubMedView ArticleGoogle Scholar
- Kobe B, Kajava AV: The leucine-rich repeat as a protein recognition motif. Curr Opin Struct Biol. 2001, 11 (6): 725-732. 10.1016/S0959-440X(01)00266-4.PubMedView ArticleGoogle Scholar
- Enkhbayar P, Kamiya M, Osaki M, Matsumoto T, Matsushima N: Structural principles of leucine-rich repeat (LRR) proteins. Proteins. 2004, 54 (3): 394-403. 10.1002/prot.10605.PubMedView ArticleGoogle Scholar
- Bella J, Hindle KL, McEwan PA, Lovell SC: The leucine-rich repeat structure. Cell Mol Life Sci. 2008, 65 (15): 2307-2333. 10.1007/s00018-008-8019-0.PubMedView ArticleGoogle Scholar
- Buchanan SG, Gay NJ: Structural and functional diversity in the leucine-rich repeat family of proteins. Prog Biophys Mol Biol. 1996, 65 (1–2): 1-44. 10.1016/S0079-6107(96)00003-X.PubMedView ArticleGoogle Scholar
- Kobe B, Deisenhofer J: The leucine-rich repeat: a versatile binding motif. Trends Biochem Sci. 1994, 19 (10): 415-421. 10.1016/0968-0004(94)90090-6.PubMedView ArticleGoogle Scholar
- Kobe B, Deisenhofer J: A structural basis of the interactions between leucine-rich repeats and protein ligands. Nature. 1995, 374 (6518): 183-186. 10.1038/374183a0.PubMedView ArticleGoogle Scholar
- Kajava AV, Vassart G, Wodak SJ: Modeling of the three-dimensional structure of proteins with the typical leucine-rich repeats. Structure. 1995, 3 (9): 867-877. 10.1016/S0969-2126(01)00222-2.PubMedView ArticleGoogle Scholar
- Hocking AM, Shinomura T, McQuillan DJ: Leucine-rich repeat glycoproteins of the extracellular matrix. Matrix Biol. 1998, 17 (1): 1-19. 10.1016/S0945-053X(98)90121-4.PubMedView ArticleGoogle Scholar
- Kresse H, Schönherr E: Proteoglycans of the extracellular matrix and growth control. J Cell Physiol. 2001, 189 (3): 266-274. 10.1002/jcp.10030.PubMedView ArticleGoogle Scholar
- Martinon F, Tschopp J: NLRs join TLRs as innate sensors of pathogens. Trends Immunol. 2005, 26 (8): 447-454. 10.1016/j.it.2005.06.004.PubMedView ArticleGoogle Scholar
- Matilla A, Radrizzani M: The Anp32 family of proteins containing leucine-rich repeats. Cerebellum. 2005, 4 (1): 7-18. 10.1080/14734220410019020.PubMedView ArticleGoogle Scholar
- Chen Y, Aulia S, Li L, Tang BL: AMIGO and friends: an emerging family of brain-enriched, neuronal growth modulating, type I transmembrane proteins with leucine-rich repeats (LRR) and cell adhesion molecule motifs. Brain Res Rev. 2006, 51 (2): 265-274. 10.1016/j.brainresrev.2005.11.005.PubMedView ArticleGoogle Scholar
- Hohenester E, Hussain S, Howitt JA: Interaction of the guidance molecule Slit with cellular receptors. Biochem Soc Trans. 2006, 34 (Pt 3): 418-421.PubMedView ArticleGoogle Scholar
- McHale L, Tan X, Koehl P, Michelmore RW: Plant NBS-LRR proteins: adaptable guards. Genome Biol. 2006, 7 (4): 212-10.1186/gb-2006-7-4-212.PubMedView ArticleGoogle Scholar
- Pancer Z, Cooper MD: The evolution of adaptive immunity. Annu Rev Immunol. 2006, 24: 497-518. 10.1146/annurev.immunol.24.021605.090542.PubMedView ArticleGoogle Scholar
- Bierne H, Sabet C, Personnic N, Cossart P: Internalins: a complex family of leucine-rich repeat-containing proteins in Listeria monocytogenes. Microbes Infect. 2007, 9 (10): 1156-1166. 10.1016/j.micinf.2007.05.003.PubMedView ArticleGoogle Scholar
- Gay NJ, Gangloff M: Structure and function of Toll receptors and their ligands. Annu Rev Biochem. 2007, 76: 141-165. 10.1146/annurev.biochem.76.060305.151318.PubMedView ArticleGoogle Scholar
- Vanhoorelbeke K, Ulrichts H, Walle Van de G, Fontayne A, Deckmyn H: Inhibition of platelet glycoprotein Ib and its antithrombotic potential. Curr Pharm Des. 2007, 13 (26): 2684-2697. 10.2174/138161207781662867.PubMedView ArticleGoogle Scholar
- Kobe B, Deisenhofer J: Crystal structure of porcine ribonuclease inhibitor, a protein with leucine-rich repeats. Nature. 1993, 366 (6457): 751-756. 10.1038/366751a0.PubMedView ArticleGoogle Scholar
- Price SR, Evans PR, Nagai K: Crystal structure of the spliceosomal U2B"-U2A' protein complex bound to a fragment of U2 small nuclear RNA. Nature. 1998, 394 (6694): 645-650. 10.1038/29234.PubMedView ArticleGoogle Scholar
- Kim JI, Lee CJ, Jin MS, Lee CH, Paik SG, Lee H, Lee JO: Crystal structure of CD14 and its implications for lipopolysaccharide signaling. J Biol Chem. 2005, 280 (12): 11347-11351. 10.1074/jbc.M414607200.PubMedView ArticleGoogle Scholar
- Bell JK, Askins J, Hall PR, Davies DR, Segal DM: The dsRNA binding site of human Toll-like receptor 3. Proc Natl Acad Sci USA. 2006, 103 (23): 8792-8797. 10.1073/pnas.0603245103.PubMedView ArticleGoogle Scholar
- Jin MS, Kim SE, Heo JY, Lee ME, Kim HM, Paik SG, Lee H, Lee JO: Crystal structure of the TLR1-TLR2 heterodimer induced by binding of a tri-acylated lipopeptide. Cell. 2007, 130 (6): 1071-1082. 10.1016/j.cell.2007.09.008.PubMedView ArticleGoogle Scholar
- Kim HM, Park BS, Kim JI, Kim SE, Lee J, Oh SC, Enkhbayar P, Matsushima N, Lee H, Yoo OJ, Lee JO: Crystal structure of the TLR4-MD-2 complex with bound endotoxin antagonist Eritoran. Cell. 2007, 130 (5): 906-917. 10.1016/j.cell.2007.08.002.PubMedView ArticleGoogle Scholar
- Tan X, Calderon-Villalobos LI, Sharon M, Zheng C, Robinson CV, Estelle M, Zheng N: Mechanism of auxin perception by the TIR1 ubiquitin ligase. Nature. 2007, 446 (7136): 640-645. 10.1038/nature05731.PubMedView ArticleGoogle Scholar
- Scott PG, Grossmann JG, Dodd CM, Sheehan JK, Bishop PN: Light and X-ray scattering show decorin to be a dimer in solution. J Biol Chem. 2003, 278 (20): 18353-18359. 10.1074/jbc.M211936200.PubMedView ArticleGoogle Scholar
- Scott PG, McEwan PA, Dodd CM, Bergmann EM, Bishop PN, Bella J: Crystal structure of the dimeric protein core of decorin, the archetypal small leucine-rich repeat proteoglycan. Proc Natl Acad Sci USA. 2004, 101 (44): 15633-15638. 10.1073/pnas.0402976101.PubMedView ArticleGoogle Scholar
- Howitt JA, Clout NJ, Hohenester E: Binding site for Robo receptors revealed by dissection of the leucine-rich repeat region of Slit. EMBO J. 2004, 23 (22): 4406-4412. 10.1038/sj.emboj.7600446.PubMedView ArticleGoogle Scholar
- Scott PG, Dodd CM, Bergmann EM, Sheehan JK, Bishop PN: Crystal structure of the biglycan dimer and evidence that dimerization is essential for folding and stability of class I small leucine-rich repeat proteoglycans. J Biol Chem. 2006, 281 (19): 13324-13332. 10.1074/jbc.M513470200.PubMedView ArticleGoogle Scholar
- McEwan PA, Scott PG, Bishop PN, Bella J: Structural correlations in the family of small leucine-rich repeat proteins and proteoglycans. J Struct Biol. 2006, 155 (2): 294-305. 10.1016/j.jsb.2006.01.016.PubMedView ArticleGoogle Scholar
- Ameye L, Young MF: Mice deficient in small leucine-rich proteoglycans: novel in vivo models for osteoporosis, osteoarthritis, Ehlers-Danlos syndrome, muscular dystrophy, and corneal diseases. Glycobiology. 2002, 12 (9): 107R-116R. 10.1093/glycob/cwf065.PubMedView ArticleGoogle Scholar
- Schaefer L, Iozzo RV: Biological functions of the small leucine-rich proteoglycans: from genetics to signal transduction. J Biol Chem. 2008, 283 (31): 21305-21309. 10.1074/jbc.R800020200.PubMedView ArticleGoogle Scholar
- Waddington RJ, Roberts HC, Sugars RV, Schönherr E: Differential roles for small leucine-rich proteoglycans in bone formation. Eur Cell Mater. 2003, 6: 12-21.PubMedGoogle Scholar
- Ross MD, Bruggeman LA, Hanss B, Sunamoto M, Marras D, Klotman ME, Klotman PE: Podocan, a novel small leucine-rich repeat protein expressed in the sclerotic glomerular lesion of experimental HIV-associated nephropathy. J Biol Chem. 2003, 278 (35): 33248-33255. 10.1074/jbc.M301299200.PubMedView ArticleGoogle Scholar
- Ohta K, Lupo G, Kuriyama S, Keynes R, Holt CE, Harris WA, Tanaka H, Ohnuma S: Tsukushi functions as an organizer inducer by inhibition of BMP activity in cooperation with chordin. Dev Cell. 2004, 7 (3): 347-358. 10.1016/j.devcel.2004.08.014.PubMedView ArticleGoogle Scholar
- O'Connor E, Eisenhaber B, Dalley J, Wang T, Missen C, Bulleid N, Bishop PN, Trump D: Species specific membrane anchoring of nyctalopin, a small leucine-rich repeat protein. Hum Mol Genet. 2005, 14 (13): 1877-1887. 10.1093/hmg/ddi194.PubMedView ArticleGoogle Scholar
- Mochida Y, Parisuthiman D, Kaku M, Hanai J, Sukhatme VP, Yamauchi M: Nephrocan, a novel member of the small leucine-rich repeat protein family, is an inhibitor of transforming growth factor-β signaling. J Biol Chem. 2006, 281 (47): 36044-36051. 10.1074/jbc.M604787200.PubMedView ArticleGoogle Scholar
- Letunic I, Copley RR, Pils B, Pinkert S, Schultz J, Bork P: SMART 5: domains in the context of genomes and networks. Nucleic Acids Res. 2006, D257-D260. 10.1093/nar/gkj079. 34 Database
- Finn RD, Mistry J, Schuster-Böckler B, Griffiths-Jones S, Hollich V, Lassmann T, Moxon S, Marshall M, Khanna A, Durbin R, Eddy SR, Sonnhammer EL, Bateman A: Pfam: clans, web tools and services. Nucleic Acids Res. 2006, D247-D251. 10.1093/nar/gkj149. 34 Database
- Mulder NJ, Apweiler R, Attwood TK, Bairoch A, Bateman A, Binns D, Bork P, Buillard V, Cerutti L, Copley R, Courcelle E, Das U, Daugherty L, Dibley M, Finn R, Fleischmann W, Gough J, Haft D, Hulo N, Hunter S, Kahn D, Kanapin A, Kejariwal A, Labarga A, Langendijk-Genevaux PS, Lonsdale D, Lopez R, Letunic I, Madera M, Maslen J, McAnulla C, McDowall J, Mistry J, Mitchell A, Nikolskaya AN, Orchard S, Orengo C, Petryszak R, Selengut JD, Sigrist CJ, Thomas PD, Valentin F, Wilson D, Wu CH, Yeats C: New developments in the InterPro database. Nucleic Acids Res. 2007, D224-D228. 10.1093/nar/gkl841. 35 Database
- Di Matteo A, Federici L, Mattei B, Salvi G, Johnson KA, Savino C, De Lorenzo G, Tsernoglou D, Cervone F: The crystal structure of polygalacturonase-inhibiting protein (PGIP), a leucine-rich repeat protein involved in plant defense. Proc Natl Acad Sci USA. 2003, 100 (17): 10124-10128. 10.1073/pnas.1733690100.PubMedView ArticleGoogle Scholar
- Huizinga EG, Tsuji S, Romijn RA, Schiphorst ME, de Groot PG, Sixma JJ, Gros P: Structures of glycoprotein Iba and its complex with von Willebrand factor A1 domain. Science. 2002, 297 (5584): 1176-1179. 10.1126/science.107355.PubMedView ArticleGoogle Scholar
- Uff S, Clemetson JM, Harrison T, Clemetson KJ, Emsley J: Crystal structure of the platelet glycoprotein Iba N-terminal domain reveals an unmasking mechanism for receptor activation. J Biol Chem. 2002, 277 (38): 35657-35663. 10.1074/jbc.M205271200.PubMedView ArticleGoogle Scholar
- He XL, Bazan JF, McDermott G, Park JB, Wang K, Tessier-Lavigne M, He Z, Garcia KC: Structure of the Nogo receptor ectodomain: a recognition module implicated in myelin inhibition. Neuron. 2003, 38 (2): 177-185. 10.1016/S0896-6273(03)00232-0.PubMedView ArticleGoogle Scholar
- Barton WA, Liu BP, Tzvetkova D, Jeffrey PD, Fournier AE, Sah D, Cate R, Strittmatter SM, Nikolov DB: Structure and axon outgrowth inhibitor binding of the Nogo-66 receptor and related proteins. EMBO J. 2003, 22 (13): 3291-3302. 10.1093/emboj/cdg325.PubMedView ArticleGoogle Scholar
- UniProt Consortium: The universal protein resource (UniProt). Nucleic Acids Res. 2008, D190-D195. 36 Database
- Bell JK, Botos I, Hall PR, Askins J, Shiloach J, Segal DM, Davies DR: The molecular structure of the Toll-like receptor 3 ligand-binding domain. Proc Natl Acad Sci USA. 2005, 102 (31): 10976-10980. 10.1073/pnas.0505077102.PubMedView ArticleGoogle Scholar
- Choe J, Kelker MS, Wilson IA: Crystal structure of human toll-like receptor 3 (TLR3) ectodomain. Science. 2005, 309 (5734): 581-585. 10.1126/science.1115253.PubMedView ArticleGoogle Scholar
- Honeybee Genome Sequence Consortium: Insights into social insects from the genome of the honeybee Apis mellifera. Nature. 2006, 443 (7114): 931-949. 10.1038/nature05260.View ArticleGoogle Scholar
- Nene V, Wortman JR, Lawson D, Haas B, Kodira C, Tu ZJ, Loftus B, Xi Z, Megy K, Grabherr M, Ren Q, Zdobnov EM, Lobo NF, Campbell KS, Brown SE, Bonaldo MF, Zhu J, Sinkins SP, Hogenkamp DG, Amedeo P, Arensburger P, Atkinson PW, Bidwell S, Biedler J, Birney E, Bruggner RV, Costas J, Coy MR, Crabtree J, Crawford M, Debruyn B, Decaprio D, Eiglmeier K, Eisenstadt E, El-Dorry H, Gelbart WM, Gomes SL, Hammond M, Hannick LI, Hogan JR, Holmes MH, Jaffe D, Johnston JS, Kennedy RC, Koo H, Kravitz S, Kriventseva EV, Kulp D, Labutti K, Lee E, Li S, Lovin DD, Mao C, Mauceli E, Menck CF, Miller JR, Montgomery P, Mori A, Nascimento AL, Naveira HF, Nusbaum C, O'leary S, Orvis J, Pertea M, Quesneville H, Reidenbach KR, Rogers YH, Roth CW, Schneider JR, Schatz M, Shumway M, Stanke M, Stinson EO, Tubio JM, Vanzee JP, Verjovski-Almeida S, Werner D, White O, Wyder S, Zeng Q, Zhao Q, Zhao Y, Hill CA, Raikhel AS, Soares MB, Knudson DL, Lee NH, Galagan J, Salzberg SL, Paulsen IT, Dimopoulos G, Collins FH, Birren B, Fraser-Liggett CM, Severson DW: Genome sequence of Aedes aegypti, a major arbovirus vector. Science. 2007, 316 (5832): 1718-1723. 10.1126/science.1138878.PubMedView ArticleGoogle Scholar
- Sea Urchin Genome Sequencing Consortium: The genome of the sea urchin Strongylocentrotus purpuratus. Science. 2006, 314 (5801): 941-952. 10.1126/science.1133609.View ArticleGoogle Scholar
- Putnam NH, Butts T, Ferrier DE, Furlong RF, Hellsten U, Kawashima T, Robinson-Rechavi M, Shoguchi E, Terry A, Yu JK, Benito-Gutiérrez EL, Dubchak I, Garcia-Fernàndez J, Gibson-Brown JJ, Grigoriev IV, Horton AC, de Jong PJ, Jurka J, Kapitonov VV, Kohara Y, Kuroki Y, Lindquist E, Lucas S, Osoegawa K, Pennacchio LA, Salamov AA, Satou Y, Sauka-Spengler T, Schmutz J, Shin-I T, Toyoda A, Bronner-Fraser M, Fujiyama A, Holland LZ, Holland PW, Satoh N, Rokhsar DS: The amphioxus genome and the evolution of the chordate karyotype. Nature. 2008, 453 (7198): 1064-1071. 10.1038/nature06967.PubMedView ArticleGoogle Scholar
- Dehal P, Satou Y, Campbell RK, Chapman J, Degnan B, De Tomaso A, Davidson B, Di Gregorio A, Gelpke M, Goodstein DM, Harafuji N, Hastings KE, Ho I, Hotta K, Huang W, Kawashima T, Lemaire P, Martinez D, Meinertzhagen IA, Necula S, Nonaka M, Putnam N, Rash S, Saiga H, Satake M, Terry A, Yamada L, Wang HG, Awazu S, Azumi K, Boore J, Branno M, Chin-Bow S, DeSantis R, Doyle S, Francino P, Keys DN, Haga S, Hayashi H, Hino K, Imai KS, Inaba K, Kano S, Kobayashi K, Kobayashi M, Lee BI, Makabe KW, Manohar C, Matassi G, Medina M, Mochizuki Y, Mount S, Morishita T, Miura S, Nakayama A, Nishizaka S, Nomoto H, Ohta F, Oishi K, Rigoutsos I, Sano M, Sasaki A, Sasakura Y, Shoguchi E, Shin-i T, Spagnuolo A, Stainier D, Suzuki MM, Tassy O, Takatori N, Tokuoka M, Yagi K, Yoshizaki F, Wada S, Zhang C, Hyatt PD, Larimer F, Detter C, Doggett N, Glavina T, Hawkins T, Richardson P, Lucas S, Kohara Y, Levine M, Satoh N, Rokhsar DS: The draft genome of Ciona intestinalis: insights into chordate and vertebrate origins. Science. 2002, 298 (5601): 2157-2167. 10.1126/science.1080049.PubMedView ArticleGoogle Scholar
- Small KS, Brudno M, Hill MM, Sidow A: A haplome alignment and reference sequence of the highly polymorphic Ciona savignyi genome. Genome Biol. 2007, 8 (3): R41-10.1186/gb-2007-8-3-r41.PubMedView ArticleGoogle Scholar
- Huxley-Jones J, Robertson DL, Boot-Handford RP: On the origins of the extracellular matrix in vertebrates. Matrix Biol. 2007, 26 (1): 2-11. 10.1016/j.matbio.2006.09.008.PubMedView ArticleGoogle Scholar
- Matsushima N, Ohyanagi T, Tanaka T, Kretsinger RH: Super-motifs and evolution of tandem leucine-rich repeats within the small proteoglycans – biglycan, decorin, lumican, fibromodulin, PRELP, keratocan, osteoadherin, epiphycan, and osteoglycin. Proteins. 2000, 38 (2): 210-225. 10.1002/(SICI)1097-0134(20000201)38:2<210::AID-PROT9>3.0.CO;2-1.PubMedView ArticleGoogle Scholar
- Delsuc F, Brinkmann H, Chourrout D, Philippe H: Tunicates and not cephalochordates are the closest living relatives of vertebrates. Nature. 2006, 439 (7079): 965-968. 10.1038/nature04336.PubMedView ArticleGoogle Scholar
- Shintani S, Sato A, Toyosawa S, O'hUigin C, Klein J: Biglycan-like extracellular matrix genes of agnathans and teleosts. J Mol Evol. 2000, 51 (4): 363-373.PubMedGoogle Scholar
- Escriva H, Manzon L, Youson J, Laudet V: Analysis of lamprey and hagfish genes reveals a complex history of gene duplications during early vertebrate evolution. Mol Biol Evol. 2002, 19 (9): 1440-1450.PubMedView ArticleGoogle Scholar
- Kuratani S, Kuraku S, Murakami Y: Lamprey as an evo-devo model: lessons from comparative embryology and molecular phylogenetics. Genesis. 2002, 34 (3): 175-183. 10.1002/gene.10142.PubMedView ArticleGoogle Scholar
- Blaschke UK, Hedbom E, Bruckner P: Distinct isoforms of chicken decorin contain either one or two dermatan sulfate chains. J Biol Chem. 1996, 271 (48): 30347-30353. 10.1074/jbc.271.48.30347.PubMedView ArticleGoogle Scholar
- Henry SP, Takanosu M, Boyd TC, Mayne PM, Eberspaecher H, Zhou W, de Crombrugghe B, Hook M, Mayne R: Expression pattern and gene characterization of asporin. a newly discovered member of the leucine-rich repeat protein family. J Biol Chem. 2001, 276 (15): 12212-12221. 10.1074/jbc.M011290200.PubMedView ArticleGoogle Scholar
- Jaillon O, Aury JM, Brunet F, Petit JL, Stange-Thomann N, Mauceli E, Bouneau L, Fischer C, Ozouf-Costaz C, Bernot A, Nicaud S, Jaffe D, Fisher S, Lutfalla G, Dossat C, Segurens B, Dasilva C, Salanoubat M, Levy M, Boudet N, Castellano S, Anthouard V, Jubin C, Castelli V, Katinka M, Vacherie B, Biémont C, Skalli Z, Cattolico L, Poulain J, De Berardinis V, Cruaud C, Duprat S, Brottier P, Coutanceau JP, Gouzy J, Parra G, Lardier G, Chapple C, McKernan KJ, McEwan P, Bosak S, Kellis M, Volff JN, Guigó R, Zody MC, Mesirov J, Lindblad-Toh K, Birren B, Nusbaum C, Kahn D, Robinson-Rechavi M, Laudet V, Schachter V, Quétier F, Saurin W, Scarpelli C, Wincker P, Lander ES, Weissenbach J, Roest Crollius H: Genome duplication in the teleost fish Tetraodon nigroviridis reveals the early vertebrate proto-karyotype. Nature. 2004, 431 (7011): 946-957. 10.1038/nature03025.PubMedView ArticleGoogle Scholar
- Christoffels A, Koh EG, Chia JM, Brenner S, Aparicio S, Venkatesh B: Fugu genome analysis provides evidence for a whole-genome duplication early during the evolution of ray-finned fishes. Mol Biol Evol. 2004, 21 (6): 1146-1151. 10.1093/molbev/msh114.PubMedView ArticleGoogle Scholar
- Peer Van de Y: Tetraodon genome confirms Takifugu findings: most fish are ancient polyploids. Genome Biol. 2004, 5 (12): 250-10.1186/gb-2004-5-12-250.PubMedView ArticleGoogle Scholar
- Meyer A, Peer Van de Y: From 2R to 3R: evidence for a fish-specific genome duplication (FSGD). Bioessays. 2005, 27 (9): 937-945. 10.1002/bies.20293.PubMedView ArticleGoogle Scholar
- Taylor JS, Braasch I, Frickey T, Meyer A, Peer Van de Y: Genome duplication, a trait shared by 22000 species of ray-finned fish. Genome Res. 2003, 13 (3): 382-390. 10.1101/gr.640303.PubMedView ArticleGoogle Scholar
- Woods IG, Wilson C, Friedlander B, Chang P, Reyes DK, Nix R, Kelly PD, Chu F, Postlethwait JH, Talbot WS: The zebrafish gene map defines ancestral vertebrate chromosomes. Genome Res. 2005, 15 (9): 1307-1314. 10.1101/gr.4134305.PubMedView ArticleGoogle Scholar
- Scott PG, Winterbottom N, Dodd CM, Edwards E, Pearson CH: A role for disulphide bridges in the protein core in the interaction of proteodermatan sulphate and collagen. Biochem Biophys Res Commun. 1986, 138 (3): 1348-1354. 10.1016/S0006-291X(86)80431-4.PubMedView ArticleGoogle Scholar
- Scott PG, Nakano T, Dodd CM: Isolation and characterization of small proteoglycans from different zones of the porcine knee meniscus. Biochim Biophys Acta. 1997, 1336 (2): 254-262.PubMedView ArticleGoogle Scholar
- Bredrup C, Knappskog PM, Majewski J, Rødahl E, Boman H: Congenital stromal dystrophy of the cornea caused by a mutation in the decorin gene. Invest Ophthalmol Vis Sci. 2005, 46 (2): 420-426. 10.1167/iovs.04-0804.PubMedView ArticleGoogle Scholar
- Rødahl E, Van Ginderdeuren R, Knappskog PM, Bredrup C, Boman H: A second decorin frame shift mutation in a family with congenital stromal corneal dystrophy. Am J Ophthalmol. 2006, 142 (3): 520-521. 10.1016/j.ajo.2006.03.064.PubMedView ArticleGoogle Scholar
- Khan A, Al-Saif A, Kambouris M: A novel KERA mutation associated with autosomal recessive cornea plana. Ophthalmic Genet. 2004, 25 (2): 147-152. 10.1080/13816810490514397.PubMedView ArticleGoogle Scholar
- Matsushima N, Tachi N, Kuroki Y, Enkhbayar P, Osaki M, Kamiya M, Kretsinger RH: Structural analysis of leucine-rich-repeat variants in proteins associated with human diseases. Cell Mol Life Sci. 2005, 62 (23): 2771-2791. 10.1007/s00018-005-5187-z.PubMedView ArticleGoogle Scholar
- Majava M, Bishop PN, Hägg P, Scott PG, Rice A, Inglehearn C, Hammond CJ, Spector TD, Ala-Kokko L, Männikkö M: Novel mutations in the small leucine-rich repeat protein/proteoglycan (SLRP) genes in high myopia. Hum Mutat. 2007, 28 (4): 336-344. 10.1002/humu.20444.PubMedView ArticleGoogle Scholar
- Lord PW, Selley JN, Attwood TK: CINEMA-MX: a modular multiple alignment editor. Bioinformatics. 2002, 18 (10): 1402-1403. 10.1093/bioinformatics/18.10.1402.PubMedView ArticleGoogle Scholar
- Sigrist CJ, Cerutti L, Hulo N, Gattiker A, Falquet L, Pagni M, Bairoch A, Bucher P: PROSITE: a documented database using patterns and profiles as motif descriptors. Brief Bioinform. 2002, 3 (3): 265-274. 10.1093/bib/3.3.265.PubMedView ArticleGoogle Scholar
- de Castro E, Sigrist CJ, Gattiker A, Bulliard V, Langendijk-Genevaux PS, Gasteiger E, Bairoch A, Hulo N: ScanProsite: detection of PROSITE signature matches and ProRule-associated functional and structural residues in proteins. Nucleic Acids Res. 2006, W362-W365. 10.1093/nar/gkl124. 34 Web Server
- Ye J, McGinnis S, Madden TL: BLAST: improvements for better sequence analysis. Nucleic Acids Res. 2006, W6-W9. 10.1093/nar/gkl164. 34 Web Server
- Birney E, Andrews D, Caccamo M, Chen Y, Clarke L, Coates G, Cox T, Cunningham F, Curwen V, Cutts T, Down T, Durbin R, Fernandez-Suarez XM, Flicek P, Gräf S, Hammond M, Herrero J, Howe K, Iyer V, Jekosch K, Kähäri A, Kasprzyk A, Keefe D, Kokocinski F, Kulesha E, London D, Longden I, Melsopp C, Meidl P, Overduin B, Parker A, Proctor G, Prlic A, Rae M, Rios D, Redmond S, Schuster M, Sealy I, Searle S, Severin J, Slater G, Smedley D, Smith J, Stabenau A, Stalker J, Trevanion S, Ureta-Vidal A, Vogel J, White S, Woodwark C, Hubbard TJ: Ensembl 2006. Nucleic Acids Res. 2006, D556-D561. 10.1093/nar/gkj133. 34 Database
- Thompson JD, Higgins DG, Gibson TJ: CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 1994, 22 (22): 4673-4680. 10.1093/nar/22.22.4673.PubMedView ArticleGoogle Scholar
- Felsenstein J: PHYLIP (Phylogeny Inference Package). Version 3.6. 2006, Department of Genome Sciences, University of Washington, SeattleGoogle Scholar
- Huelsenbeck JP: MrBayes: Bayesian inference of phylogeny. 2000, Department of Biology, University of Rochester. Distributed by the authorGoogle Scholar
- Evans SV: SETOR: hardware-lighted three-dimensional solid model representations of macromolecules. J Mol Graph. 1993, 11 (2): 134-138. 10.1016/0263-7855(93)87009-T.PubMedView ArticleGoogle Scholar
- Crooks GE, Hon G, Chandonia JM, Brenner SE: WebLogo: a sequence logo generator. Genome Res. 2004, 14 (6): 1188-1190. 10.1101/gr.849004.PubMedView ArticleGoogle Scholar
This article is published under license to 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.