- Research article
- Open Access
LRRCE: a leucine-rich repeat cysteine capping motif unique to the chordate lineage
BMC Genomics volume 9, Article number: 599 (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.
Many LRR proteins other than SLRPs are flanked at the N- and C-termini by disulphide-bonded caps that are thought to protect the hydrophobic core of the first and last LRRs [2, 3]. Both N-terminal and C-terminal capping motifs are described in databases for protein domain identification and analysis such as SMART , Pfam  and InterPro  (Table 1). In the LRR N-terminal capping motif (LRRNT), a single β-strand runs antiparallel to the main β-sheet and is followed by a short LRR of 20 or 21 residues (Figure 1a). The consensus sequence contains 4 cysteines in a Cx n Cx Cx m C pattern, with n and m being variable numbers.
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
The final regular expression pattern for the LRRCE capping motif and its mapping to the three-dimensional structure of bovine decorin are shown in Figure 2. Using the LRRCE regular expression pattern, a total of 175 sequences were retrieved by ScanProsite from the UniProt database (Swiss-Prot release 55.1, TrEMBL release 38.1) . Splice variants were excluded from the search and sequence duplicates were filtered, resulting in a non-redundant set of 110 UniProt sequences. Figure 3 shows a selection of aligned LRRCE sequences from this non-redundant set (an extended list of sequences including accession codes is available in Additional File 1), plus two sequences from urochordates (discussed later). The LRRCE regular expression pattern was accurate: all sequences in this UniProt set were of proteins with LRR architecture and had a repeat structure that identified them as canonical SLRPs, with the LRRNT capping motif at the N-terminus and the LRRCE capping motif at the C-terminus. The pattern was also comprehensive: all proteins known to belong to the canonical group of SLRPs or annotated as similar to them were retrieved in the search.
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
A cluster analysis of the extended set of LRRCE-containing sequences is shown in the form of a phylogenetic tree inferred from a multiple sequence alignment of the LRRCE motifs (Figure 4 and Additional File 2). The sequences group themselves in the canonical SLRP classes (I, II and III), plus an additional cluster containing ECM2 and related sequences. Class I includes asporin (ASPN), decorin (DCN) and biglycan (BGN) sequences. Class II contains two groups and could in fact be subdivided into IIa with lumican (LUM) and fibromodulin (FMOD) sequences, and IIb with osteomodulin (OMD), keratocan (KERA) and prolargin (PRELP) sequences. Class III includes sequences for opticin (OPTC), epiphycan (EPYC) and osteoglycin (OGN), also known as mimecan. The new "class A" is discussed later and includes the sequences of extracellular matrix protein 2 (ECM2), ECM2-like protein from the X chromosome (ECMX), ECM2-like predicted proteins upstream of the DCN gene in fish genomes (ECMZ), and the small leucine-rich repeat protein from Ciona intestinalis and Ciona savignyi (SLRP1) (see below). The same clustering was obtained using the complete SLRP sequences (not shown), indicating that the LRRCE sequences on their own are useful for canonical group SLRP classification.
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
Similarity searches using LRRCE sequences were then carried out against the genomes of the early chordates amphioxus (lancelet), Branchiostoma floridae , and the ascidians Ciona intestinalis  and Ciona savignyi . Only one gene containing the LRRCE motif was found for each Ciona species. This gene, referred to here as SLRP1, has been proposed as representative of the ancestor of all modern canonical SLRPs . Two model assemblies have been proposed for the Ciona intestinalis gene resulting in proteins with different numbers of LRRs (Figure 5). These two assemblies only differ in the prediction of one additional exon-intron boundary in the short model, which results in the skipping of one and a half exons from the long model. The long model assembly for the Ciona intestinalis SLRP1 gene has been confirmed by EST data (Figure 5). The resulting protein sequence has 15 LRRs in which repeats alternate their lenghts following a predominant 21-24-26 pattern; this repeat structure is very similar to that seen in the ECM2 sequences from vertebrates . The cluster analysis of LRRCE motifs (Figure 4) places the SLRP1 sequences of Ciona intestinalis and Ciona savigny in the same group as ECM2 and related sequences. The alternating pattern of repeat lengths (21-24-26) is common to SLRP1 as well as a Ciona gene representative of the podocan ancestor , suggesting that such an alternating repeat sequence was already present in a common precursor of these two lineages. This observation is consistent with the concept of tandem LRR supermotifs characteristic of the evolutionary history of the SLRP family .
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
Searches for LRRCE motifs against the current release of the sea lamprey genome Petromyzon marinus (Genome Sequence Center, Washington University) produced six LRRCE-containing sequences. Two of them correspond to previously reported biglycan-like proteins , whereas the four additional hits correspond to predicted partial sequences similar to decorin (two sequences), epiphycan (one sequence) and keratocan (one sequence). We have named these sequences biglycan-like proteins 1 and 2 (BGL1 and BGL2), decorin-like proteins 1 and 2 (DCL1 and DCL2), epiphycan-like protein (EPYL) and keratocan-like protein (KERAL) (Table 2). Although these data should still be considered preliminary, the sequences represent the earliest examples to date of class I, II and III SLRPs, suggesting that the divergence of the SLRP ancestor into three classes, following gene duplication, occurred before the lamprey-gnathostome split. Completion of the lamprey genome and possible identification of additional copies of EPYL and KERAL genes will clarify the relation between individual SLRP gene duplications and the large-scale gene or genome duplication that is thought to have occurred before the divergence of lampreys from gnathostomes [57, 61, 62].
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).
Figures 6, 7, 8, 9, 10 show partial phylogenetic trees for each SLRP class, inferred from multiple sequence alignments of LRRCE-containing sequences from a subset of genomes (see Methods). Not all SLRP sequences have been found in all genomes. This could be owing to incomplete coverage in some of these genomes or to genuine absence of particular genes due to gene loss. For example, it has been known for some time that chickens do not have a BGN gene , and unsurprisingly no such gene has been found in searches against the chicken genome. Interestingly, the BGN gene is present in fishes, reptiles and mammals. The genomes of ray-finned fishes (zebrafish, stickleback, medaka, and the pufferfishes) appear to have an extended set of SLRPs, with additional copies of fibromodulin (FMOD2), lumican (LUM2), osteoglycin (OGN2) or biglycan (BGN2) genes (Figures 6, 7, 8, 9), and yet seem to have mostly lost the genes for OMD and OPTC (with the exception of zebrafish, which appears to have retained a particularly complete set of SLRP genes).
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).
This sequence, referred to here as ECMX for reasons explained later, is predicted upstream of the biglycan gene (Figure 11), in a completely analogous manner to the ECM2-ASPN gene tandem seen in mammalian and fish genomes (Figure 11). The predicted ECMX sequence for zebrafish (Table 3) shows a LRR structure highly similar to ECM2, with 15 LRRs and a sequence of repeat lengths highly reminiscent of the 21-24-26 pattern mentioned earlier.
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).
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.
Kajava AV: Structural diversity of leucine-rich repeat proteins. J Mol Biol. 1998, 277 (3): 519-527. 10.1006/jmbi.1998.1643.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
Matilla A, Radrizzani M: The Anp32 family of proteins containing leucine-rich repeats. Cerebellum. 2005, 4 (1): 7-18. 10.1080/14734220410019020.
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.
Hohenester E, Hussain S, Howitt JA: Interaction of the guidance molecule Slit with cellular receptors. Biochem Soc Trans. 2006, 34 (Pt 3): 418-421.
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.
Pancer Z, Cooper MD: The evolution of adaptive immunity. Annu Rev Immunol. 2006, 24: 497-518. 10.1146/annurev.immunol.24.021605.090542.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
Sea Urchin Genome Sequencing Consortium: The genome of the sea urchin Strongylocentrotus purpuratus. Science. 2006, 314 (5801): 941-952. 10.1126/science.1133609.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
Felsenstein J: PHYLIP (Phylogeny Inference Package). Version 3.6. 2006, Department of Genome Sciences, University of Washington, Seattle
Huelsenbeck JP: MrBayes: Bayesian inference of phylogeny. 2000, Department of Biology, University of Rochester. Distributed by the author
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.
Crooks GE, Hon G, Chandonia JM, Brenner SE: WebLogo: a sequence logo generator. Genome Res. 2004, 14 (6): 1188-1190. 10.1101/gr.849004.
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.
HP is responsible for creating the sequence patterns, data collection and phylogenetic analysis. JHJ contributed to the production and analysis of phylogenetic trees and to the evolutionary and biological insights of the manuscript. RPBH and PNB contributed to the evolutionary and biological insights of the manuscript. TKA contributed to the data analysis and helped in the supervision of the bioinformatic analysis. JB conceived the study, analyzed the data and wrote the manuscript. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 2: high-resolution version of Figure 4. Larger version of Figure 4, with legible sequence names at the end of the phylogenetic tree branches. (TIFF 2 MB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Park, H., Huxley-Jones, J., Boot-Handford, R.P. et al. LRRCE: a leucine-rich repeat cysteine capping motif unique to the chordate lineage. BMC Genomics 9, 599 (2008). https://doi.org/10.1186/1471-2164-9-599
- Accession Code
- Branchiostoma Floridae
- Canonical Group
- Regular Expression Pattern