Genomic characterization and expression pattern of Bf/C2 and C4 in miiuy croaker and molecular evolution analysis on mammals and fishes

The Establishment of Complement System Is from Gene Duplication and Domain Shuffling

The mammalian complement system constitutes a highly sophisticated body defense machinery. The evolutionary origin of the complement system can be traced to Coelenterata as the presence of the central component C3 and two activation proteases BF and MASP. In the present study, the main complement components were screened and analyzed from the genomes of different species in metazoan subphyla/phyla. C1q with classical domains can be traced to Annelida, and ficolin and MBL to Urochordata. C1r and C1s are only found in Chondrichthyes and even higher species, and MASP is traced to Coelenterata. In the evolutionary tree, C1r from Vertebrates is close to MASP1/2/3 from Deuterostomia and Coelenterata, and C1s from Vertebrates is close to MASP-like protease (MASPL) from Arthropoda, Mollusca, and Annelida. C2, BF, and DF can be traced to Mollusca, Coelenterata, and Porifera, respectively. There are no clear C2 and BF branches in the evolutionary tree. C3 can be traced to Coelenterata, and C4 and C5 are only in Chondrichthyes and even higher species. There are three clear C3, C4, and C5 branches in the evolutionary tree. C6-like (C6L) and C8 can be traced to Urochordata, and C7-like (C7L) can be traced to Cephalochordara. C6L, C7L, and C8 from Urochordata and Cephalochordara provide the structural conditions for the formation of Vertebrate MAC components. The findings unveil the evolutionary principles of the complement system and provide insight into its sophistication.

Structural comparison and expression function analysis of BF/C2 in Ctenopharyngodon idella and Squaliobarbus curriculus

Ctenopharyngodon idella and Squaliobarbus curriculus, members of the Cyprinidae family and Yaroideae subfamily, have shown different levels of resistance to grass carp reo virus (GCRV), with S. curriculus exhibiting higher resilience. In the pursuit to explore the distinctions in the structural and expression traits of BF/C2 (A,B) between the two species, we conducted an analysis involving the cloning and examination of various coding sequences (CDS). We successfully cloned the CDS of ci-BF/C2A and ci-BF/C2B from C. idella, which spanned 2259 bp and 2514 bp respectively, encoding 752 and 837 amino acids. Similarly, the CDS of sc-BF/C2A and sc-BF/C2B from S. curriculus were cloned, featuring lengths of 1353 bp and 2517 bp and encoding 450 and 838 amino acids, respectively. A chromosome collinearity assessment revealed that ci-BF/C2A demonstrated collinearity with sc-BF/C2A, a finding not replicated with ci-BF/C2B and sc-BF/C2B. Delving into gene structure, we discerned that ci-BF/C2A harbored a greater number of Tryp_SPc domains compared to sc-BF/C2A. Following this, we engineered and purified six prokaryotic recombinant proteins: CI-BF/C2A, CI-BF/C2A1 (a variant resulting from the deletion of the Tryp_SPc domain of CI-BF/C2A), CI-BF/C2A2 (representing the Tryp_SPc domain of CI-BF/C2A), CI-BF/C2B, SC-BF/C2A, and SC-BF/C2B. Through serum co-incubation tests with these recombinant proteins, we established the activation of the complement marker C3 in each case. Utilizing fluorescence quantitative expression analysis, we observed ubiquitous expression of ci-BF/C2A and ci-BF/C2B across all grass carp tissues, predominantly in the liver. This pattern mirrored in S. curriculus, where sc-BF/C2A was highly expressed in the gills, and sc-BF/C2B manifested notably in the liver. Kidney cell infection experiments on both species revealed enhanced resistance to GCRV post-incubation with the recombinant proteins. Notably, cells treated with SC-BF/C2 (A, B) exhibited pronounced resilience compared to those treated with CI-BF/C2 (A, B, A1, A2). However, cells incubated with CI-BF/C2A1 and CI-BF/C2A2 showed strengthen resistance relative to cells treated with CI-BF/C2A and CI-BF/C2B. In GCRV infection trials on grass carp, ci-BF/C2A and ci-BF/C2B expressions reached a zenith on the seventh day post-infection, highlighting a distinctive functional mode in immune defense against GCRV infection orchestrated by BF/C2. The empirical data underscores the pivotal role of the Tryp_SPc domain in immune responses to GCRV infection, pinpointing its influence on ci-BF/C2A expression. Conclusively, this investigation provides a foundational understanding of the unique immune function characteristics of BF/C2 in grass carp, paving the way for further scholarly exploration in this realm.

Tracing and exploring the evolutionary origin and systematic function of fish complement C9

Complement C9, as a member of terminal complement component (TCC) protein, plays important roles in innate immunity. However, some complement components appear to show difference and evolutionary complexity between higher and lower vertebrates. Hence, it is essential to carry on a study of evolutionary origin and systematic function of C9 in fish and non-fish vertebrates. This study aims to explore the complement gene evolution and potential function in fish based on molecular and structural biology. Herein, we found complete divergence of C9 throughout the gene evolution. The optimal codons of C9 sequences tended to be closer to the genomes of lower vertebrates compared to higher vertebrates. Further, conserved amino acids in the C9 TMH1 region were identified, implying their potential functional association with MAC growth and pore formation. Transposons and simple repeats, as gene elements, exhibited a differential distribution in the genomic regions in different animal groups but were sparsely scattered around the sixth exon (TMH1 region). Notably, this demonstrated the regulatory complexity of the C9 gene in higher vertebrates. The negative selection pressures on fish and non-fish groups improved both the sequence conservation and similarity. Through gene/protein regulatory network and pathway analyses, the systematic function of C9 protein was showcased; thus, we could reveal the divergence of the systematic function of C9 across species from different evolutionary positions. In addition, more complicated functions of C9 in higher vertebrates could established by the altered spatial conformation of the protein. Collectively, the present study illustrates the C9 gene evolutionary process and the difference in its systematic function across multiple species. Such advances provide new insights for understanding the evolutionary and potential functions of complement C9.

Fish complement C4 gene evolution and gene/protein regulatory network analyses and simulated stereo conformation of C4-MASP-2 protein complex

Complement C4 is a central protein by acting as pivotal molecule in the activation of the complement system. More than a decade ago, C4 gene duplication had been found in several species including fish, revealing the evolutionary origin of C4 gene. However, the evolutionary pattern and systematic function of C4 are still limited. In this study, C4 D and H types in different species groups were completely diverged. The codon usage of C4 H type in higher vertebrates were much closer to their own genome environment, in contrast to lower vertebrates, suggesting that the evolution may provide the dynamic for homogeneous codon usage between specific gene and genome. Multiple C4 sequence alignment showed that the sequences were conserved among different species. However, sequence similarity was obviously different between species C4 D and H type. Negative selection pressure was found on C4 gene evolution and it may be one of the possible reasons for the sequence broad similarity and conservation among interspecies. Proteins from C4 protein-protein interaction (PPI) network were enriched in more hematopoiesis, infections, diseases and immune-related pathways in human than zebrafish. The result suggested that the functional complexities of C4 isotypes are distinct in species from different evolutionary positions. The simulated C4 protein structures between human and grass carp shared structural similarity and the stereo structures of grass carp C4-MASP-2 protein complexes were further simulated according to a study of human. These results suggested that the interaction between C4 and MASP-2 proteins may also exist in grass carp. Our results can provide an insight for the evolutionary process of C4 and better understanding to the potential mechanism of interaction between C4 and MASP-2 in fish species.

Complement proteins detected through iTRAQ-based proteomics analysis of serum from black carp Mylopharyngodon piceus in response to experimentally induced Aeromonas hydrophila infection

The black carp Mylopharyngodon piceus is one of the culturally important '4 famous domestic fishes' in China. Recently, infectious diseases caused by Aeromonas hydrophila have drastically altered the operation of the black carp farming industry. In the present study, isobaric tags for relative and absolute quantitation (iTRAQ) were combined with mass spectrometry analysis to screen for differentially abundant black carp serum proteins in response to experimentally induced A. hydrophila infection. A total of 86 differentially abundant proteins were quantified at 24 h post-infection, including 78 down-regulated proteins and 8 up-regulated proteins. The down-regulated proteins included complement C1q subcomponent subunit C, complement factor B/C2A, complement factor B/C2B, complement C3-Q1, complement C3, and complement C4-2. Bioinformatic analysis indicated that the differentially abundant proteins were mainly associated with complement and coagulation cascades (27.9%). Moreover, real-time PCR (qPCR) analysis revealed changes in the gene expression of both C3 and B/C2A in blood cells, liver, kidney, gills, and intestines of the black carp infected with A. hydrophila. However, mRNA expression levels did not consistently correlate with the corresponding protein levels. A polyclonal antibody was prepared using a synthetic C3 peptide. Immunofluorescence analysis showed that the expression of C3 in the kidney was increased with A. hydrophila infection. This work provides a useful characterization of the impact of A. hydrophila infection on the complement system of the black carp.

Do the toll-like receptors and complement systems play equally important roles in freshwater adapted Dolly Varden char (Salvelinus malma)?

Unlike the normal anadromous lifestyle, Chinese native Dolly Varden char (Salvelinus malma) is locked in land and lives in fresh water lifetime. To explore the effect of freshwater adaption on its immune system, we constructed a pooled cDNA library of hepatopancreas and spleen of Chinese freshwater Dolly Varden char (S. malma). A total of 27,829 unigenes were generated from 31,233 high-quality transcripts and 17,670 complete open reading frames (ORF) were identified. Totally 25,809 unigenes were successfully annotated and it classified more native than adaptive immunity-associated genes, and more genes involved in toll-like receptor signal pathway than those in complement and coagulation cascades (51 vs 3), implying the relative more important role of toll-like receptors than the complement system under bacterial injection for the freshwater Dolly Varden char. These huge different numbers of TLR and complement system identified in freshwater Dolly Varden char probably caused by distinct evolution pressure patterns between fish TLR and complement system, representative by TLR3 and TLR5 as well as C4 and C6, respectively, which were under purifying and positively selecting pressure, respectively. Further seawater adaptation experiment and the comparison study with our library will no doubt be helpful to elucidate the effect of freshwater adaption of Chinese native Dolly Varden char on its immune system.

Characterization, expression, and evolutionary analysis of new TLR3 and TLR5M genes cloned from the spiny eel Mastacembelus armatus

Toll-like receptors (TLRs) play an important role in innate and adaptive immunity. Here, we identify two new TLRs from the spiny eel Mastacembelus armatus (TLR3 and membrane TLR5M). Both MaTLR3 and MaTLR5M were expressed in all tested tissues; expression was highest in liver and spleen, respectively. After infection with Vibrio parahaemolyticus, expression of both TLRs fluctuated and differed significantly from controls at several time points. The predicted three-dimensional model of the MaTLR3 and MaTLR5M proteins indicates that most sites under positive selection were located in the extracellular domains of TLRs. Evolutionary analysis detected positively selected sites in the ancestral lineages of vertebrates, amphibians and reptiles. Multiple ML methods recovered 10 positively selected sites in teleost TLR3 and 24 in TLR5M, and most sites were located in leucine-rich repeat domain, possibly related to an "arms-race" co-evolution with pathogens.

Characterization and evolutionary analysis of duplicated C7 in miiuy croaker

The complement system, as one of the most sophisticated innate immune system, plays an important role in defense against invading microorganisms. The complement component C7 participates in the cytolytic phase of complement activation through a series of polymerization reactions with other terminal complement components. In this study, we derived two C7 genes from the whole genome of miiuy croaker which were the consequence of the fish-specific genome duplication. Our data showed that miiuy croaker C7-1 and C7-2 genes shared same structure domains. The analysis of gene synteny showed that high degree conserved of synteny was retained between miiuy croaker and other teleosts, and miiuy croaker had a relatively closer relationship with fugu. The expression of C7-1 and C7-2 in miiuy croaker healthy tissues revealed that they were ubiquitously expressed in all ten tested tissues. Besides, the immune response of C7-1 and C7-2 were different in spleen with Vibrio anguillarum, Staphylococcus aureus, poly I:C and LPS at 24 h post-injection, respectively. Furthermore, the expression patterns of C7-1 and C7-2 were different in liver, spleen and kidney after infected with V. anguillarum at different time-point. Evolutionary analysis showed that all the ancestral lineages underwent positive selection except for the ancestral lineages of fish C7-2, indicated that the ancestral lineages of fish C7-1 genes undertook more pressures than C7-2 in defense against the invading microorganisms. Meanwhile, a series of maximum likelihood methods were used to explore the evolutionary patterns on extant vertebrates' C7 genes. Three and one positive selection sites were found in extant mammalian C7 genes and fish C7-2 genes, but no positive selection site was found in extant fishes C7-1 genes. The result showed that extant fish C7-2 genes undertook more pressures compared with C7-1. In conclusion, fish C7-1 and C7-2 gene underwent different evolutionary patterns.