Bioinformatic identification of genes encoding Clq-domaincontaining proteins in zebrafish

C1q is the first subcomponent of classical pathway in the complement system and a major link between innate and acquired immunities. The globular （gC1q） domain similar with C1q was also found in many non-complement C1q-domain-containing （C1qDC） proteins which have similar crystal structure to that of the multifunctional tumor necrosis factor （TNF） ligand family, and also have diverse functions. In this study, we identified a total of 52 independent gene sequences encoding C1q-domain-containing proteins through comprehensive searches of zebrafish genome, cDNA and EST databases. In comparison to 31 orthologous genes in human and different numbers in other species, a significant selective pressure was suggested during vertebrate evolution. Domain organization of C1q-domain-containing （C1qDC） proteins mainly includes a leading signal peptide, a collagen-like region of variable length, and a C-terminal C1q domain. There are 11 highly conserved residues within the C1q domain, among which 2 are invariant within the zebrafish gene set. A more extensive database searches also revealed homologous C1qDC proteins in other vertebrates, invertebrates and even bacterium, but no homologous sequences for encoding C1qDC proteins were found in many species that have a more recent evolutionary history with zebrafish. Therefore, further studies on C1q-domain-containing genes among different species will help us understand evolutionary mechanism of innate and acquired immunities.     

The complete complement of C1q-domain-containing proteins in Homo sapiens

The C-terminal domains of the A, B, C chains of C1q subcomponent of C1 complex represent a common structural motif, the C1q domain, that is found in a diverse range of proteins. We analyzed the human genome for the complete complement of this family and have identified a total of 31 independent gene sequences. The predominant organization of C1q-domain-containing (C1qDC) proteins includes a leading signal peptide, a collagen-like region of variable length, and a C-terminal C1q domain. There are 15 highly conserved residues within the C1q domain, among which 8 are invariant within the human gene set and these are predicted to cluster within the hydrophobic core of the protein. We suggest a 3-subfamily classification based on sequence homology. For some C1qDC-encoding genes, strict orthology has been retained throughout vertebrate evolution and these examples suggest a highly specific functional role for C1qDC proteins that has been under significant selective pressure. Alternatively, individual species have co-opted C1qDC proteins for roles that are highly specific to their biology, suggesting an evolutionary strategy of gene duplication and functional diversification. A more extensive analysis of the evolutionary relationship of C1qDC proteins reveals an ancient rooting, with clear members found in eubacterial species. Curiously, we have been unable to identify C1qDC-encoding genes in many eukaryotic genomcs, such as Sacchromyces cerivisae and C. elegans, suggesting that the retention or loss of this gene family throughout evolution has been sporadic.     

C1q蛋白家族的结构、分布、分类和功能

C1q蛋白家族由众多含C1q结构域的蛋白组成, 从细菌到高等哺乳动物中都有分布。这类蛋白由一条信号肽、胶原样区(Collage-like region, CLR)和C1q球状结构域(Globular C1q domain, gC1q)组成。C1q蛋白家族根据其结构特点, 可分为三大类分子：C1q、C1q-like和ghC1q。C1q是补体经典途径的起始分子, 能够识别免疫复合物, 启动补体系统经典途径; 此外, 作为一种模式识别受体分子(Pattern recognition receptor, PRR), 它可以结合种类繁多的配体。C1q-like蛋白的结构类似于C1q分子, 含有CLR和gC1q结构域, 在水蛭中参与神经系统的修复, 在脊椎动物中实现从凝集素到免疫球蛋白结合分子的功能转变, 参与补体系统的激活。ghC1q蛋白只具有gC1q结构域和一段短的N末端序列, 包括分泌型蛋白(sghC1q)和非分泌型蛋白(cghC1q)。sghC1q在无脊椎动物固有免疫系统中发挥重要作用; 脊椎动物中的sghC1q可作为一类新型跨神经元调节因子, 在大脑的许多区域调节突触发育和突触可塑性。cghC1q基因最早可追溯至芽孢杆菌属的细菌中, 具有典型的gC1q果冻卷结构, 说明gC1q结构域有着非常悠久的进化历程且结构高度保守。文章对C1q蛋白家族的结构、分布、分类以及功能进行综述, 以期为从事该领域研究的科研人员提供有益参考。     

Crystal structure of zebrafish complement 1qA globular domain

C1q contains three globular domains (C1qgD) that are the key functional component of the classical complement system. C1qgD can interact with important immune molecules, including IgG and C‐reactive protein (CRP) to form defense systems to protect animals. Here, the first non‐mammalian structure, zebrafish C1qA globular domain (Dare‐C1qAgD) was solved. Although the overall architecture of Dare‐C1qAgD is similar to human C1qA, residues involved in C1qBgD, C1qCgD, and CRP binding are somewhat different while residues involved in IgG binding are not present in zebrafish. The structure gives insight into how human and fish C1qA evolved from an ancestral protein.     

Adiponectin binds C1q and activates the classical pathway of complement

The adipose-specific protein adiponectin binds to a number of target molecules, including damaged endothelium and the surface of apoptotic cells. However, the significance of this binding remains unclear. This study demonstrates the binding of purified C1q to recombinant adiponectin under physiological conditions, and the dependence of this upon Ca⁺⁺ and Mg⁺⁺. Binding was enhanced by metaperiodate-mediated destruction of glucosylgalactosyl sugars on adiponectin. Adiponectin was bound by the globular domain of the A chain of collagenase-digested C1q, and C1q binding induced deposition of C4 and C3 through activation of the classical complement pathway. After Western blotting, affinity-purified adiponectin from human serum bound C1q, whereas adiponectin in whole serum did not, unless pre-treated with metaperiodate. These results suggest adiponectin is member of the pattern-recognition family of defence collagens, able to bind target molecules and activate complement. It may therefore play an important role in innate immunity and autoimmune phenomena.     

Binding of complement subcomponent C1q to Streptococcus pyogenes: evidence for interactions with the M5 and FcRA76 proteins

Binding of C1q, the first component of the complement system, to some human pathogens has been earlier reported. In the present study, direct binding of C1q to group A streptococci (GAS) of various serotypes as well as some other Gram-positive and Gram-negative species was demonstrated. The interaction between C1q and GAS was investigated more in detail. In hot neutral extracts of a number of GAS strains two components of 64 and 52 kDa, respectively, bound C1q; alkaline and SDS extracts yielded the 52 kDa component as the main C1q-binding substance. Trypsin treatment of the SDS extracts of two GAS strains suggested the C1q-binding component(s) to be of protein nature. C1q-binding material purified from the SDS extract of an avirulent strain, type T27, was separated in 12% SDS-PAGE and probed in Western blot with human C1q and fibrinogen, conjugated to horse radish peroxidase (HRP) as well as rabbit IgG antibodies complexed to HRP (PAP system). The 52 kDa component was non-reactive with fibrinogen or rabbit IgG. However, C1q-binding components purified from the alkaline extracts of two M-positive strains revealed strong binding of either fibrinogen (type M5) or both fibrinogen and rabbit IgG (type M76); the molecular mass of these components, 55 kDa and 43–40 kDa, respectively, was in agreement with the reported molecular mass of the M5 and FcRA76 proteins. Our findings suggest that C1q may interact with GAS through certain M-family proteins as well as by a so far unidentified surface factor of protein nature occurring in most GAS strains. The involvement of M-family proteins, regarded as virulence factors of these organisms, may suggest the interaction of GAS with C1q as biologically important.     

Phylogenetic analyses of two "archaeal" genes in thermotoga maritima reveal multiple transfers between archaea and bacteria

The genome sequence of Thermotoga maritima revealed that 24% of its open reading frames (ORFs) showed the highest similarity scores to archaeal genes in BLAST analyses. Here we screened 16 strains from the genus Thermotoga and other related Thermotogales for the occurrence of two of these "archaeal" genes: the gene encoding the large subunit of glutamate synthase (gltB) and the myo-inositol 1P synthase gene (ino1). Both genes were restricted to the Thermotoga species within the Thermotogales. The distribution of the two genes, along with results from phylogenetic analyses, showed that they were acquired from Archaea during the divergence of the Thermotogales. Database searches revealed that three other bacteria-Dehalococcoides ethenogenes, Sinorhizobium meliloti, and Clostridium difficile-possess archaeal-type gltBs, and the phylogenetic analyses confirmed at least two lateral gene transfer (LGT) events between Bacteria and Archaea. These LGT events were also strongly supported by gene structure data, as the three domains in bacterial-type gltB are homologous to three independent ORFs in Archaea and Bacteria with archaeal-type gltBs. The ino1 gene has a scattered distribution among Bacteria, and apart from the Thermotoga strains it is found only in Aquifex aeolicus, D. ethenogenes, and some high-G+C Gram-positive bacteria. Phylogenetic analysis of the ino1 sequences revealed three highly supported prokaryotic clades, all containing a mixture of archaeal and bacterial sequences, and suggested that all bacterial ino1 genes had been recruited from archaeal donors. The Thermotoga strains and A. aeolicus acquired this gene independently from different archaeal species. Although transfer of genes from hyperthermophilic Archaea may have facilitated the evolution of bacterial hyperthermophily, between-domain transfers also affect mesophilic species. For hyperthermophiles, we hypothesize that LGT may be as much a consequence as the cause of adaptation to hyperthermophily.     

Mapping and conservation of the group-specific component gene in mouse

The group-specific component (GC), also known as the vitamin D-binding protein, transports vitamin D and its metabolites in plasma to target tissues throughout the body. The GC gene shares an evolutionary origin with genes encoding albumin (ALB) and alpha-fetoprotein (AFP). All three genes are descendants of an evolutionary ancestor that arose from an intragenic triplication. As a result, each gene is composed of three homologous domains. The study described here characterizes and compares mouse GC to the corresponding nucleotide and amino acid sequences of GC from human and rat. The deduced amino acid sequence of mouse GC was 78% identical to human and 91% identical to rat GC. The results suggest that, unlike the corresponding sequences in the ALB and AFP genes, chromosomal sequences encoding the first domain and the leader sequence of the GC gene have specifically been conserved throughout vertebrate evolution. Protection of domain I during evolution may correlate with an important functional aspect of its sequence. The mouse GC gene was mapped to chromosome 5, where the ALB and AFP genes are also located, demonstrating conservation of the three genes in vertebrate species.     

Sequences and domain structures of mammalian, avian, amphibian and teleost tropoelastins: Clues to the evolutionary history of elastins

Tropoelastin is the monomeric form of elastin, a polymeric extracellular matrix protein responsible for properties of extensibility and elastic recoil in connective tissues of most vertebrates. As an approach to investigate how sequence and structural characteristics of tropoelastin assist in polymeric assembly and account for the elastomeric properties of this polymer, and to better understand the evolutionary history of elastin, we have identified and characterized tropoelastins from frog (Xenopus tropicalis) and zebrafish (Danio rerio), comparing these to their mammalian and avian counterparts. Unlike other species, two tropoelastin genes were expressed in zebrafish. All tropoelastins shared a predominant and characteristic alternating domain arrangement, as well as the fundamental crosslinking sequence motifs. However, zebrafish and frog tropoelastins had several unusual characteristics, including increased exon numbers and protein molecular weights, and decreased hydropathies. For all tropoelastins there was evidence of evolutionary expansion of the proteins by extensive replication of a hydrophobic-crosslinking exon pair. This was particularly apparent for zebrafish and frog tropoelastin genes, where remnants of sequence similarity were also seen in introns flanking the replicated exon pair. While overall alignment of mammalian, avian, frog and zebrafish tropoelastin sequences was not possible because of sequence variability, the C-terminal exon was well-conserved in all species. In addition, good sequence alignment was possible for several exons just upstream of the putative region of replication, suggesting that these conserved domains may represent 'primordial' core sequences present in the ancestral sequence common to all tropoelastins and in some way essential to the structure/function of elastin.     

Origin and evolution of TRIM proteins: new insights from the complete TRIM repertoire of zebrafish and pufferfish

Tripartite motif proteins (TRIM) constitute a large family of proteins containing a RING-Bbox-Coiled Coil motif followed by different C-terminal domains. Involved in ubiquitination, TRIM proteins participate in many cellular processes including antiviral immunity. The TRIM family is ancient and has been greatly diversified in vertebrates and especially in fish. We analyzed the complete sets of trim genes of the large zebrafish genome and of the compact pufferfish genome. Both contain three large multigene subsets--adding the hsl5/trim35-like genes (hltr) to the ftr and the btr that we previously described--all containing a B30.2 domain that evolved under positive selection. These subsets are conserved among teleosts. By contrast, most human trim genes of the other classes have only one or two orthologues in fish. Loss or gain of C-terminal exons generated proteins with different domain organizations; either by the deletion of the ancestral domain or, remarkably, by the acquisition of a new C-terminal domain. Our survey of fish trim genes in fish identifies subsets with different evolutionary dynamics. trims encoding RBCC-B30.2 proteins show the same evolutionary trends in fish and tetrapods: they evolve fast, often under positive selection, and they duplicate to create multigenic families. We could identify new combinations of domains, which epitomize how new trim classes appear by domain insertion or exon shuffling. Notably, we found that a cyclophilin-A domain replaces the B30.2 domain of a zebrafish fintrim gene, as reported in the macaque and owl monkey antiretroviral TRIM5α. Finally, trim genes encoding RBCC-B30.2 proteins are preferentially located in the vicinity of MHC or MHC gene paralogues, which suggests that such trim genes may have been part of the ancestral MHC.     

Sm and Sm-like proteins assemble in two related complexes of deep evolutionary origin.

A group of seven Sm proteins forms a complex that binds to several RNAs in metazoans. All Sm proteins contain a sequence signature, the Sm domain, also found in two yeast Sm-like proteins associated with the U6 snRNA. We have performed database searches revealing the presence of 16 proteins carrying an Sm domain in the yeast genome. Analysis of this protein family confirmed that seven of its members, encoded by essential genes, are homologues of metazoan Sm proteins. Immunoprecipitation revealed that an evolutionarily related subgroup of seven Sm-like proteins is directly associated with the nuclear U6 and pre-RNase P RNAs. The corresponding genes are essential or required for normal vegetative growth. These proteins appear functionally important to stabilize U6 snRNA. The two last yeast Sm-like proteins were not found associated with RNA, and neither was essential for vegetative growth. To investigate whether U6-associated Sm-like protein function is widespread, we cloned several cDNAs encoding homologous human proteins. Two representative human proteins were shown to associate with U6 snRNA-containing complexes. We also identified archaeal proteins related to Sm and Sm-like proteins. Our results demonstrate that Sm and Sm-like proteins assemble in at least two functionally conserved complexes of deep evolutionary origin.     

The Antigen Receptor (NCCRP-1) on Catfish and Zebrafish Nonspecific Cytotoxic Cells Belongs to a New Gene Family Characterized by an F-Box-Associated Domain

The catfish nonspecific cytotoxic cell receptor protein (NCCRP-1) provides an important function in target cell recognition and activation of cytotoxicity. This report identifies and characterizes a zebrafish orthologue of the catfish NCCRP-1. The zebrafish NCCRP-1 cDNA contains an open reading frame that encodes a predicted protein of 237 amino acids with a MW of 27 kDa and a pI of 5.5. Sequence similarities comparisons show that the NCCRP-1 receptors from these two phylogenetically distant species share a high degree of identity. These results suggested that NCCRP-1 performs a crucial function in innate immunity in teleosts. Further, a zebrafish 17-mer peptide corresponding to the catfish NCCRP-1 antigen-binding domain inhibited (catfish) cytotoxicity toward conventional tumor target cells (HL-60). These data appeared to indicate that the zebrafish NCCRP-1 protein may function as an antigen recognition molecule and, as such, may participate in innate immunity in teleosts. A homology search of the zebrafish NCCRP-1 protein revealed that it shares a significant level of identity with another group of proteins belonging to an F-box subfamily. These proteins share an F-box domain in the N terminus (not present in NCCRP-1) and an extremely conserved C-terminal region that has been termed the F-box-associated domain (FBA). The FBA is currently of unknown function. A new gene family is proposed in this work, based on similarities in the FBA sequences with the catfish and zebrafish NCCRP-1 peptides. This new gene family includes several F-box domain-containing proteins and a predicted C. elegans protein. Received: 20 June 2001 / Accepted: 31 August 2001     

Confirmation of a second natural preQ1 aptamer class in Streptococcaceae bacteria

Bioinformatics searches of eubacterial genomes have yielded many riboswitch candidates where the identity of the ligand is not immediately obvious on examination of associated genes. One of these motifs is found exclusively in the family Streptococcaceae within the 5' untranslated regions (UTRs) of genes encoding the hypothetical membrane protein classified as COG4708 or DUF988. While the function of this protein class is unproven, a riboswitch binding the queuosine biosynthetic intermediate pre-queuosine(1) (preQ(1)) has been identified in the 5' UTR of homologous genes in many Firmicute species of bacteria outside of Streptococcaceae. Here we show that a representative of the COG4708 RNA motif from Streptococcus pneumoniae R6 also binds preQ(1). Furthermore, representatives of this RNA have structural and molecular recognition characteristics that are distinct from those of the previously described preQ(1) riboswitch class. PreQ(1) is the second metabolite for which two or more distinct classes of natural aptamers exist, indicating that natural aptamers utilizing different structures to bind the same metabolite may be more common than is currently known. Additionally, the association of preQ(1) binding RNAs with most genes encoding proteins classified as COG4708 strongly suggests that these proteins function as transporters for preQ(1) or another queuosine biosynthetic intermediate.     

Engineering by homologous recombination: exploring sequence and function within a conserved fold

In nature similar protein folds accommodate distant sequences and support diverse functions. This observation coupled with the recognition that proteins can tolerate many homologous substitutions inspires protein engineers to use recombination to search for new functions within sequences encoding structurally related molecules. These searches have led to proteins with novel activities, diversified specificities and greater stabilities. Computational methods that exploit structural and evolutionary information are being used to design highly mutated yet still natively folded chimeric proteins and protein libraries.     

Signatures of selection acting on the innate immunity gene Toll-like receptor 2 (TLR2) during the evolutionary history of rodents

Patterns of selection acting on immune defence genes have recently been the focus of considerable interest. Yet, when it comes to vertebrates, studies have mainly focused on the acquired branch of the immune system. Consequently, the direction and strength of selection acting on genes of the vertebrate innate immune defence remain poorly understood. Here, we present a molecular analysis of selection on an important receptor of the innate immune system of vertebrates, the Toll-like receptor 2 (TLR2), across 17 rodent species. Although purifying selection was the prevalent evolutionary force acting on most parts of the rodent TLR2, we found that codons in close proximity to pathogen-binding and TLR2-TLR1 heterodimerization sites have been subject to positive selection. This indicates that parasite-mediated selection is not restricted to acquired immune system genes like the major histocompatibility complex, but also affects innate defence genes. To obtain a comprehensive understanding of evolutionary processes in host-parasite systems, both innate and acquired immunity thus need to be considered.