- 1Key Laboratory of Bio-Resources and Eco-Environment of Ministry of Education, College of Life Sciences, Sichuan University, Chengdu, China
- 2Key Laboratory of Sichuan for Fishes Conservation and Utilization in the Upper Reaches of the Yangtze River, College of Life Science, Neijiang Normal University, Neijiang, China
- 3Shenzhen Key Lab of Marine Genomics, Guangdong Provincial Key Lab of Molecular Breeding in Marine Economic Animals, BGI Academy of Marine Sciences, BGI Marine, BGI, Shenzhen, China
Toll‐like receptor 1 (TLR1) mediates the innate immune response to a variety of microbes through recognizing cell wall components (such as bacterial lipoproteins) in mammals. However, the detailed molecular mechanism of TLR1 involved in pathogen immunity in the representative hybrid yellow catfish (Pelteobagrus fulvidraco ♀ × P. vachelli ♂) has not been well studied. In the present study, we identified the TLR1 gene from the hybrid yellow catfish, and further comparative synteny data from multiple species confirmed that the TLR1 gene is highly conserved in teleosts. Phylogenetic analysis revealed distinguishable TLR1s in diverse taxa, suggesting consistence in evolution of the TLR1 proteins with various species. Structural prediction indicated that the three-dimensional structures of TLR1 proteins are relatively conserved among different taxa. Positive selection analysis showed that purifying selection dominated the evolutionary process of TLR1s and TLR1-TIR domain in both vertebrates and invertebrates. Expression pattern analysis based on the tissue distribution showed that TLR1 mainly transcribed in the gonad, gallbladder and kidney, and the mRNA levels of TLR1 in kidney were remarkably up-regulated after Aeromonas hydrophila stimulation, indicating that TLR1 participates in the inflammatory responses to exogenous pathogen infection in hybrid yellow catfish. Homologous sequence alignment and chromosomal location indicated that the TLR signaling pathway is very conserved in the hybrid yellow catfish. The expression patterns of TLR signaling pathway related genes (TLR1- TLR2 - MyD88 - FADD - Caspase 8) were consistent after pathogen stimulation, revealing that the TLR signaling pathway is triggered and activated after A. hydrophila infection. Our findings will lay a solid foundation for better understanding the immune roles of TLR1 in teleosts, as well as provide basic data for developing strategies to control disease outbreak in hybrid yellow catfish.
1 Introduction
Yellow catfish (Pelteobagrus fulvidraco) has become an important economic cultured fish species in China and East Asia due to its rapid growth speed, good taste, and high nutritional values (1). However, intensive aquaculture leads to crowding and excessive stress, which usually result in outbreaks of bacterial infections with high mortality of cultured fishes (2). Aeromonas hydrophila is a typical rod-shaped gram-negative bacterium with worldwide distribution in aquatic environments (3, 4). Hydrophilic bacteria are well-studied as fish pathogens, which may cause a series of diseases, including active Aeromonas sepsis in carp, tilapia, perch, catfish and salmon (5).
A. hydrophila is the pathogenic bacteria of bacterial hemorrhagic sepsis caused by multidrug-resistant (MDR), which has increased exponentially in the past decade and has reached an alarming rate, leading to a major problem in the aquaculture industry of China (6). At the same time, the outbreak of active motile Aeromonad septicemia (MAS) in various fishes caused by A. hydrophila has attracted worldwide attentions (7). In previous studies, researchers observed that tilapia and channel catfish infected with A. hydrophila would be damaged in the gills, liver, and intestines, leading to histopathological changes in the infected organs (8, 9).
In the infection process, the innate immune system detect exogenous bacteria via sensing cell wall-related components, including lipopolysaccharides, lipoproteins, peptidoglycans, and flagellin (10). Each of these molecules is defined as a pathogen-associated molecular pattern (PAMP) and induces an inflammatory response to realize the host defense (11, 12). The first characterized and most widely studied pattern recognition receptor (PRR) in vertebrates and invertebrates is Toll-like receptors (TLRs) (13, 14). TLRs are important parts of PRR, participating in innate immune defense against exogenous pathogen invasion (15). TLRs represent a large superfamily of type I transmembrane glycoproteins, some of which are shared in multiple species, while others are more restricted in their limited distributions (16).
After pathogen recognition, the Toll/Interleukin-1 receptor (TIR) domain of the TLR triggers downstream signal transductions. Previous findings showed that almost all the TLRs activate the common signaling pathway via the TIR-containing linker MyD88, which leads to the activation of NF-κB to promote the transcription of cytokine genes for involvement in the inflammation (17). The interleukin-1 receptor (IL-1R)/TLR superfamily was first identified in 1998 as a protein family containing the Toll-IL-1 receptor domain (18). In that time, many orphan receptors were found in the IL-1R branch, and TLRs had not been proven to be the key innate immune receptor for sensing microbial products (19). Various members of the TLR1 family are usually identified in the genomes of birds and mammals, and they seem to be the consequence of successive rounds of tandem gene replication from the ancestral gene (20). Previous studies on TLR1 have mainly focused on mammals, and it has been reported that TLR2 dimerizes with TLR1 to recognize microbial triacyl lipoproteins, and with TLR6 to recognize diacyl lipopeptides present in mycoplasma, lipoteichoic acid in gram-positive bacteria, or zymosan in yeast (21, 22). The interaction of TLR1 and TLR2 enhances activation of NF-κB in response to synthetic lipopeptides (23). In addition, acylated lipoprotein analogs are recognized by TLR1 at the first time (23). TLR1 can also exert its anti-fungal host defense effect through β-defensin 3, which subsequently activates immune cells through TLR1/TLR2 (24). Certain mutations in the TLR1 affect the susceptibility and immune response to exogenous bacterial pathogens (16).
It has been reported that proteins in living cells interact specifically or non-specifically with a large number of biomolecules (25). To understand the behavior of proteins under macromolecular crowding conditions in cells, it is essential to observe their spatial structure at the atomic level in the physiological environment (26). Areas with essential functions tend to evolve more slowly, and a comparative sequence analysis would be able to identify and characterize the regulatory regions of the genome that have functioned well (27, 28). The evolution rate dN/dS, representing the ratio of non-synonymous substitution rate to synonymous substitution rate, is calculated from the selection pressure of protein coding genes for helping to uncover the mechanism of molecular evolution (29). These above-mentioned methods were employed in our present research in accordance with our previous reports (30, 31).
TLR1 gene has been cloned and characterized in various fish species, such as Chinese perch (Siniperca chuatsi), orange-spotted grouper (Epinephelus coioides), Dabry’s sturgeon (Acipenser dabryanus) and Atlantic salmon (Salmo salar) (32–35). However, the specific molecular mechanism of TLR1 involved in pathogen immune responses remains unknown in hybrid yellow catfish. The published genome at the chromosome-level of the yellow catfish (36) provides a good reference for our present research. In this study, we attempt to characterize the TLR1 gene (including prediction of protein structures, construction of a phylogenetic tree, and performance of natural selection analysis) for quantification of tissue distribution and its response to bacterial pathogen stimulation, which are helpful for revealing the molecular mechanism of TLR1 in the immune response in hybrid yellow catfish.
2 Materials and methods
2.1 Identification of TLR1 from yellow catfish genome
The genome data of yellow catfish were downloaded from NCBI (GenBank assembly accession number: PRJNA494039), and the BLAST search was performed based on protein sequence similarity to identify the TLR1 gene. As previously reported (36), the TLR1 protein was annotated. The detailed sequence was available for characterization.
2.2 Sequence processing and gene synteny analysis
Genomes and corresponding annotation files of TLR1 in nine representative species, including banded blenny (Salarias fasciatus), Crimson tide cichlid (Pundamilia nyererei), yellow perch (Perca flavescens), Japanese medaka (Oryzias latipes), black rockcod (Notothenia coriiceps), yellow catfish (P. fulvidraco), Tibetan plateau loach (Triplophysa tibetana), and bicolor damselfish (Stegastes partitus) were downloaded from the NCBI databases. Amino acid sequences of representative TLR1 in different fish species were identified by annotation and confirmed with conserved functional motifs and domains. A comparative genomic survey was conducted to recognize the genetic loci of TLR1 genes in yellow catfish as well as other examined fishes.
2.3 Phylogenetic analysis and functional domain identification
To explore the origin and evolution of vertebrate TLR1 genes, we downloaded TLR1 protein sequences of 71 representative species from the NCBI database (Supplementary Table 1). A phylogenetic tree of these TLR1 genes was constructed using the Maximum likelihood method, and the final topology was evaluated with the Poisson correction method (37). Gaps were removed from the entire alignment. Bootstrap values were calculated for evaluation of the stability and reliability of the tree with 1,000 duplications (38). Domain architectures of these TLR1 proteins were predicted using the Simple Modular Architecture Research Tool (39, 40).
2.4 Predicted spatial structures of the TLR1 proteins
In this study, the latest version of DeepMind AlphaFold2 (https://github.com/deepmind/alphafold, accessed on 11 August 2022) was used to predict the spatial structures of TLR1 proteins in multiple species, and the corresponding data were visualized using PyMol 2.5 (41, 42).
2.5 Modes and strength of natural selection on TLR1 genes in invertebrates and vertebrates
Synonymous substitutions in protein coding sequences do not result in amino acid variation, which occurs in the third or sometimes the first base of the examined codon (43). Therefore, the gene-level method based on the ratio of non-synonymous (Ka) and synonymous (Ks) substitution rates was used here to detect the potential positive selection signals that are closely related to protein-coding genes in various species (44). The Ka/Ks ratio of 1, <1 or >1 in the coding sequence of the studied protein can be interpreted as neutral mutation, negative (purified) selection, or positive (diversified) selection, respectively. The genomic data of five vertebrate classes (including Mammalia, Amphibia, Aves, Reptilia, and Osteichthyes; 103 species in total) and two invertebrate classes (Insecta and Bivalvia; 13 species in total) from NCBI were examined to determine the origin history of vertebrate-like TLR1s.
To detect the TIR domain architectures of all Toll-like candidates, HMMER was used to blast against the Pfam database (45). Full-length TLR proteins and the TIR domains of each Toll-like candidate were collected for subsequent analyses. For the evolutionary landscape of each TLR gene, the coding sequence alignments were generated by MUSCLE (align codons) (46). Meanwhile, the codeml method implemented in the PAML v4.7 package (47) was employed in this study to estimate dN/dS, dN, number of nonsynonymous substitutions per nonsynonymous site, dS, and number of synonymous substitutions were calculated using Datamonkey webserver (http://classic.datamonkey.org.php, accessed on 11 August 2022) (48, 49).
2.6 Identification of chromosome localization of genes involved in the TLR signaling pathway
Genes involved in the TLR signaling pathway were identified by using the KEGG (https://www.genome.jp/kegg/) database. Subsequently, the genome (GCA_022655615.1) of yellow catfish were downloaded and used as the reference genome. Then, TLR1, TLR2, Caspase 8, MyD88, FADD, TOLLIP isoform 1 and TOLLIP isoform 2 genes were aligned to the reference genome through tblastn to determine the corresponding position and copy number on the chromosome. Furthermore, the mg2c (version 2.0) online mapping website (http://mg2c.iask.in/mg2c_v2.0/, accessed on 11 August 2022) were employed to visualize the location information of the TLR1, TLR2, Caspase 8, MyD88, FADD, TOLLIP isoform 1 and TOLLIP isoform 2 gene.
2.7 Sample processing and pathogen challenge
The original group of hybrid yellow catfish were provided by a local fish farm in the suburb of Chengdu City, Sichuan Province, China. We randomly selected 6 fishes from thirty fishes (100 ± 2.0g) for studying the tissue distribution of TLR1 mRNA. Another thirty fishes (100 ± 2.0g) were selected randomly for studying the effect of bacterial infection on TLR1 expression levels. In order to carry out the bacterial infection experiment, we purchased A. hydrophila from Shanghai Luwei Technology Co. Ltd. (Shanghai, China). Before the practical experiment, fishes were temporarily fed indoors in a circulating water tank at 28°C for 7 days with commercial feed (protein content at about 40%; Lianyungang Tongwei Feed Co. Ltd., Jiangsu, China). The strain of A. hydrophila was inoculated on Luria-Bertani (LB) solid medium and incubated at 37°C for 24 h. A single clone was picked for inoculation in 3 mL LB liquid medium. After incubation for 24 h on a shaker, 1 × 106 CFU/mL bacteria were placed into a turbidimeter tube for turbidity measurement. Number of bacteria was measured using the routine plate counting method.
Six fishes were randomly selected and intraperitoneally injected with 50 μL of phosphate buffered saline (PBS, pH 7.2) as the control group. Another group of 24 fishes were intraperitoneally injected with 50 μL of formalin-treated hydrophile with a concentration of 1.5×107 CFU/Ml in PBS (as the experimental group). Six fishes were randomly selected from the experimental group 24 h after the injection (hpi). The sampled fishes were anesthetized with 300 mg/L MS-222 (Sigma-Aldrich, St. Louis, MO, USA), and then the kidney, skin, liver and gills were collected for subsequent experiments. All animal experiments in this study were approved by the Ethics Committee of the College of Life Sciences, Sichuan University (Permit No. SCU221208001).
2.8 RNA extraction and reverse transcription PCR (RT-PCR)
According to the manufacturer’s protocol, TRIZOL reagent (Invitrogen, Carlsbad, CA, USA) was used to extract total RNA from various tissues at different stages after the bacterial infection. The RNA samples were processed with DNase to remove potential genomic DNA contamination. Subsequently, agarose gel electrophoresis was performed and a NanoDrop 2000 spectrophotometer (Thermo Scientific, Waltham, MA, USA) was used to evaluate the quality of the separated RNAs. Reverse transcription of RNA samples was conducted using Quantscript Reverse Transcriptase Kit (Tiangen Biotech, Beijing, China). The obtained cDNA solution was used as the template for PCRs.
2.9 Quantitative real-time PCR (qRT-PCR)
A CFX Connect Real-Time PCR Detection system (Bio-Rad, Hercules, CA, USA) was used for the qRT-PCR in this study. Three fishes were used for qRT-PCR analysis. The primer sequence for TLR1, TLR2, Caspase 8, MyD88, FADD, TOLLIP isoform 1 and TOLLIP isoform 2 were listed in Table 1. β-actin was selected and used as reference gene and its amplification effect was realized with specific primers (see Table 1 for sequences). In the PCRs, SuperReal PreMix Plus (SYBR Green) (Tiangen Biotech) was used according to the manufacturer’s protocol. The PCR conducted using the following system, SYBR Premix Ex Taq (from Tiangen Biotech) was 10 μL, the cDNA used here was 40 ng, the final concentration of the primers (referred to anti and sense primer) was 0.2 μM and the making up the final volume of the system to 20 μL with double distilled water. Then the PCR were performed on a Bio-Rad T100 Thermal Cycler amplifier (Bio-Rad, Hercules, CA, USA). with the procedure, initiated 1 min at 95°C, then initiated 5 s at 95°C, annealed 20 s at 60°C and extended 20 s at 72°C. The PCR ends after 40 cycles in the last three steps. At the end of each PCR run, the melting curve analysis was performed in the range of 55°C to 99°C. For each sample, quantitative RT-PCR was performed in triplicate.
2.10 Statistical analysis
Statistical analysis was performed using the popular GraphPad Prim 9 software (GraphPad Software, San Diego, CA, USA) and corresponding calculation was conducted using the routine 2 -ΔΔCT method (50). All data are represented as mean ± SEM (n=3), and the normality and homogeneity of the variance were checked before statistical analysis. A one-way analysis of variance (ANOVA) was performed to compare the differences in expression levels between different samples.
3 Results
3.1 Protein sequence and multiple sequence alignment
Our sequence analysis results show that cDNA of the 9 TLR1s range from 2,385 to 3,557 bp, among which 2,385-2,514 bp open reading frames were predicted to encode proteins with 794-837 amino acids (Table 2). Multiple protein sequences alignment of TLR1 from ten fishes (including yellow catfish) are summarized in Figure 1. Our results showed that a transmembrane domain and the TIR domain were conserved in various fishes, and there is a leucine-rich repeat (LRR) domain in the C-terminus of TLR1s (see more details in Figure 1). Our analysis of amino acid sequence homology in different species indicated that TLR1 in yellow catfish exhibited somewhat higher levels of sequence similarity compared to most of the other species TLR1s (Table 3). According to similarity analyses in this study, TLR1 in yellow catfish is closely related to TLR1 in striped catfish (Pangasianodon hypophthalmus) with similarity of 80.83% (Table 3).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1163781/fimmu-14-1163781-HTML/image_m/fimmu-14-1163781-g001.jpg)
Figure 1 Multiple protein sequences alignment of the TLR1 proteins in nine bony fishes. The black box represents a leucine-rich repeat (LRR) domain, the solid purple box represents the transmembrane (TM) domain, and the green box represents the TIR domain.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1163781/fimmu-14-1163781-HTML/image_m/fimmu-14-1163781-t003.jpg)
Table 3 Amino Acid similarity of yellow catfish TLR1 to other species TLR1s: To determine the percentage similarity, the TLR1 sequence of yellow catfish was aligned with other species orthologues using CLUSTAL W multiple alignment.
3.2 Comparative synteny and spatial structures of the TLR1 proteins
The genetic synteny analysis showed that five genes (including Scf25, Ics10, TLR1, β-klotho, and Ubi-E2), in great majority of bony fish showed that they formed clusters (Figure 2A), although chromosomal rearrangements occurred approximately 465 million years ago (51) due to the differentiation of vertebrates (52, 53). Previous findings revealed that β-Klotho (encoded by Klb) is an obligate co-receptor that mediates FGF21 and FGF15/19 signaling for regulation of glucose and lipid metabolism (54). Bindings of ubiquitin to the substrate is generally considered to occur through the formation of iso-peptide bonds between the C-terminal glycine residue of ubiquitin and the lysine residue of the substrate (55). Individual member of the conserved family of ubiquitin-conjugating enzymes (E2s) mediates the ubiquitination and turnover of the specific substrates in the ubiquitin-dependent degradation pathway (56, 57). However, the current understanding of Scf25 and Ics10 in bony fishes is limited.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1163781/fimmu-14-1163781-HTML/image_m/fimmu-14-1163781-g002.jpg)
Figure 2 (A) Comparative synteny of TLR1 genes in nine bony fish genomes. Genes and intergenic regions are highlighted with colorful blocks and solid lines, respectively. (B) Predicted spatial structures of TLR1 proteins from eight animal taxa. These taxa include Mammalia, Amphibia, Aves, Reptilia, Chondrichthyes, Osteichthyes, Insecta, and Bivalvia. One representative species of each taxon was selected for detailed prediction.
AlphaFold is a protein structure prediction tool based on deep learning. It has achieved successes in its highly accurate structure predictions (58). Previous studies showed that the AlphaFold architecture performed high-precision training using only supervised learning of PDB data (59). Multiple sequence alignment (MSA) was added to AlphaFold2 to integrate protein structure and biological information into deep learning algorithms. Our prediction results of the TLR1 spatial structures from eight species are provided in Figure 2B. It seems that these spatial structures of TLR1 proteins were mostly horseshoe-shaped, and the corresponding TIR domains are also conserved (in green color on the upper right of Figure 2B). However, the reptile TLR1was significantly different from others. The 3D structure of TLR1 in Pinta Island tortoise (Chelonoidis abingdonii) looks like a more opened while shorter semicircle (Figure 2B, VI), although the overall shape of TLRs is similar with the TIR domain on the top of the semicircle (see Figure 2B).
3.3 Phylogenetic relationship and domains of TLR1s
In order to understand the relationship among various TLR1 proteins in both vertebrates and invertebrates, a phylogenetic tree of 71 TLR1 amino acid sequences was constructed using the maximum likelihood method (60, 61). The TLR1 protein sequence from stony coral (Stylophora pistillata) was used as the outgroup. According to the well-supported phylogenetic topology (Figure 3), our results showed that all the examined TLR1 genes are divided into invertebrate and vertebrate subgroups.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1163781/fimmu-14-1163781-HTML/image_m/fimmu-14-1163781-g003.jpg)
Figure 3 Phylogenetic relationships of TLR1 proteins in various fishes. The phylogenetic tree was constructed with the maximum likelihood method based on a dataset of protein sequences. The stony coral (Stylophora pistillata) was used as the outgroup. Genbank accession numbers of these protein sequences are provided in Supplemental Table 1. Low complexity, low complexity region; LRR, leucine-rich repeat; LRR_CT, leucine-rich repeat C-terminal; LRR_NT, leucine-rich repeat N-terminal; LRR_TYP, leucine-rich repeats, typical (most populated) subfamily; TIR, Toll/interleukin-1 receptor domain; TR, transmembrane region.
The vertebrate subgroup is further divided into six main branches (including amphibians, mammals, reptiles, birds, bony fishes and cartilaginous fishes), and the invertebrate subgroup is further divided into corals, bivalves and insects. We determined that the topologies generated by the maximum likelihood (ML) (61) and neighbor-joining (NJ) (62, 63) methods (data not shown) are similar, indicating relative stability of the phylogenetic tree. Obviously, the evolution of TLR1s is consistent with the evolution of species.
A TLR1 protein typically contains TIR, LRR, LRR_NT, LRR_CT, LRR_TPY, transmembrane domain and low repeat domain (Figure 3), but not all TLR1s have the seven domains. However, the TIR domain is present in all the 71 studied species (Figure 3), which is the signature signal transduction domain of any TLR and its linker as a scaffold for assembling protein complexes of innate immune signals (64). The TLR1 is missing in some bony fishes (the blue column in Figure 3). In fact, this polymorphism of the transmembrane domain was reported to affect the innate immune response to bacterial lipopeptides and the susceptibility to multiple pathogens (65).
3.4 Natural selection for TLR1 and TIR domain in both vertebrates and invertebrates and codon-based positive selection
To identify the strength of natural selection of TLR1 in both vertebrates and invertebrates, a selection analysis was performed. The average dN/dS values of TLR1 genes in each taxon were 0.278 (Mammalia), 0.329 (Aves), 0.294 (Reptilia), 0.240 (Osteichthyes), 0.177 (Amphibian), 0.282 (Bivalvia), and 0.406 (Insecta) respectively (Figure 4A); the average dN/dS values of TLR1-TIR in each taxon were 0.108 (Amphibia), 0.223 (Mammalia), 0.116 (Aves), 0.335 (Reptilia), 0.100 (Insecta), 0.304 (Osteichthyes), and 0.862 (Bivalvia) (Figure 4A) respectively. The dN/dS values of TLR1 and TLR1-TIR in invertebrate lineages are significantly greater than those in vertebrates (Figures 4A, B). The results of PSSs/non-PSSs are summarized (see Figures 5A, B). It seems that TLR1 genes in the invertebrate lineages exhibited a faster rate of evolution than vertebrates.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1163781/fimmu-14-1163781-HTML/image_m/fimmu-14-1163781-g004.jpg)
Figure 4 Natural selection analysis of seven animal taxa. (A) The dN/dS values of TLR1 and TLR1-TIR domain in representative Mammalia (n=326), Amphibia (n=29), Aves (n=326), Reptilia (n=135), Osteichthyes (n= 301), Insecta (n=16) and Bivalvia (n= 22), “n” is the number of species analyzed for each taxon. (B) Comparison of the average dN/dS values of TLR1 and TLR1-TIR domain in the examined taxa. *p<0.05; **p<0.01; ***p<0.001; ****p<0.0001. ns, no significant.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1163781/fimmu-14-1163781-HTML/image_m/fimmu-14-1163781-g005.jpg)
Figure 5 Analysis of positive selection site ratio in seven animal taxa. (A) The PSS/non-PSS value of TLR1 and TLR1-TIR domain in representative Mammalia, Amphibia, Aves, Reptilia, Osteichthyes, Insecta and Bivalvia. (B) Comparison of average PSS/non-PSS values of TLR1 and TLR1-TIR domain in the examined taxa.
3.5 Distribution of TLR1 gene in different tissues
The tissue distribution of sampling organs of hybrid yellow catfish used in this study were exhibited in Figure 6A. The tissues-specific expression levels (Figure 6B) demonstrated that the TLR1 gene was widely transcribed in all the 11 tested tissues (including gallbladder, intestine, liver, stomach, heart, gills, spleen, kidney, skin, gonad and brain), but exhibited somehow tissue preference. The TLR1 gene has high mRNA levels in the gonad, gallbladder and kidney, which are important components of the immune system in fishes.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1163781/fimmu-14-1163781-HTML/image_m/fimmu-14-1163781-g006.jpg)
Figure 6 (A) Tissue distribution of sampling organs used in this study. The fish model here referred to the hybrid yellow catfish. (B) Relative transcription patterns of TLR1 gene in 11 examined tissues of the hybrid yellow catfish. β-actin was selected and used as the reference gene. Each error bar represents a standard error of the mean values (n = 3).
3.6 Responses of TLR1 gene to exogenous A. hydrophila infection
After A. hydrophila infection, expression levels of TLR1 gene at different time points in hybrid yellow catfish were measured. Our results showed that TLR1 gene responded to bacterial attack, but the response pattern is tissue-specific (Figure 7). TLR1 mRNA levels in kidney were remarkably up-regulated after the bacterial infection.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1163781/fimmu-14-1163781-HTML/image_m/fimmu-14-1163781-g007.jpg)
Figure 7 Relative transcription levels of TLR1 gene in the kidney of hybrid yellow catfish infection of exogenous A.hydrophila. β-actin was selected and used as the reference gene. Different letters above the bars represent significant differences among the examined groups. Each error bar represents a standard error of the mean values (n = 3).
3.7 Chromosome location of genes involved in the TLR signaling pathway
The related genes involved in TLR signaling pathway are shown in Figure 8A, in which (TLR1-TLR2) - MyD88 - FADD- Caspase 8 composed of TLR signaling pathway finally induces apoptosis (Figure 8). Our results showed that there are 26 chromosomes in the genome of the yellow catfish, among which the longest chromosome is Chr 1 with the length of 52.4 Mb, and the shortest chromosome is Chr 26 with the length of 16.9 Mb (Figure 9). And we found that TLR1 was located at 5815026-5818519 of chromosome 15 with a single copy in the yellow catfish (GCA_022655615.1) genome (Figure 9 and Table 4). TOLLIP and FADD located on chromosome 2. Caspase 8, TLR2 and MyD88 located on chromosome 5, 18 and 25, respectively. Refer to Table 4 for the exact location of TLR signaling pathway genes on different chromosomes.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1163781/fimmu-14-1163781-HTML/image_m/fimmu-14-1163781-g008.jpg)
Figure 8 The diagram of TLR signaling pathways in hybrid yellow catfish, from top to bottom are extracellular, cytoplasmic and nuclear.
3.8 Expression profiling of genes involved in the TLR signaling pathway in response to A. hydrophila challenge in the in hybrid yellow catfish
Then the expression levels of TLR2, MyD88, FADD, TOLLIP isoform 1 and TOLLIP isoform 2 involved in the TLR signaling pathway were tracked at six time points (0 h, 3 h, 6 h, 12 h, 24 h, 48 h). Our results showed that transcript levels of all these genes were upregulated after the A. hydrophila stimulation (Figure 10). Among the six time points, the expression pattern of TLR2 was wavy, and reached the maximum at 6 h (Figure 10A). MyD88 and FADD exhibited a similar expression pattern, which gradually increased before the 12 h, reached the maximum at 12 h, and showed a process of first down-regulation and then up-regulation in the subsequent 24 h and 48 h (Figures 10B, C). The expression of Caspase 8 gene reached the maximum at 3 H, and showed a downward trend at the subsequent time points (Figure 10D). The expression levels of TOLLIP isoform 1 and TOLLIP isoform 2, two isoforms of TOLLIP, also exhibited a similar expression pattern, which reached the highest values at 12 H and 6 H, respectively (Figures 10E, F).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1163781/fimmu-14-1163781-HTML/image_m/fimmu-14-1163781-g010.jpg)
Figure 10 Relative expression levels of TLR signal pathway related genes in kidney after the A. hydrophila challenge. Transcriptional changes of hybrid yellow catfish TLR2 (A), MyD88 (B), FADD (C), caspase 8 (D), TOLLIP isoform 1 (E) and TOLLIP isoform 2 (F) in kidney were analyzed at 0, 3, 6, 12, 24 and 48 H after A. hydrophila challenge.
4 Discussion
In this study, the hybrid yellow catfish TLR1 gene was identified for further characterization and investigation. Protein sequence analysis showed that a transmembrane domain and the TIR domain were conserved in various fishes (Figure 1). Our comparative synteny data from multiple species confirmed that the identified TLR1 gene do exist in the yellow catfish genome (Figure 2A). Predicted spatial structures of each TLR1 is typical horseshoe shaped with well-positioned pockets (Figure 2B) to accommodate the triacylated lipopeptide ligand (59). It seems that the TLR1s from different taxa are relatively conserved, although reptilia present many variances (Figure 2B, VI). The phylogenetic tree of TLR1 proteins from 71 representative species is divided into two branches of vertebrates and invertebrates (Figure 3). Among these species, there are three classes of invertebrates (Anthozoa, Bivalvia, and Insecta) and six classes of vertebrates (Mammalia, Amphibia, Aves, Reptilia, Chondrichthyes, and Osteichthyes). Our genomics survey showed that the main functional domains of TLR1 proteins are LRR, LRR-C-terminal (LRR_CT) motifs and LRR-N-terminal (LRRNT) motif in the extracellular domain, the transmembrane domain and a TIR in the cytoplasmic domain, low complexity region and LRR_TYP motif (LRR-TYP represents typical LRRs, whereas LRR represents LRR outliers) (Figure 3). TLR1 proteins of the analyzed 71 species share the TIR -domain, while most TLR1 proteins in invertebrates lost the LRR_CT (Figure 3). Positive selection analysis showed that purifying selection dominated the evolutionary process of TLR1s and TLR1-TIR domain in both vertebrates and invertebrates (Figure 4). The results of PSSs/non-PSSs show that TLR1 genes in the invertebrate lineages exhibited a faster rate of evolution than vertebrates (Figure 5). qRT-PCR was applied to detect tissue distribution of the TLR1 transcript in yellow catfish, and our results demonstrated detection in all the 11 selected tissues while high levels in the immune system (Figure 6). Chromosomal mapping results showed that TLR1, TLR2, caspase 8, MyD88, FADD, TOLLIP isoform 1 and TOLLIP isoform 2 involved in the TLR signaling pathway were single-copy genes, implying that the genes of the TLR signaling pathway was relatively conserved in the yellow catfish (Figure 9 and Table 4). After infection with A. hydrophila, the transcriptional levels of TLR1, TLR2, caspase 8, MyD88, FADD, TOLLIP isoform 1 and TOLLIP isoform 2 were all significantly up-regulated in the kidney, indicating that bacterial infection induced the activation of TLR signaling pathway in the yellow catfish (Figure 10).
TIR domain is the representative signal domain of TLRs and their adaptors, serving as a scaffold for assembling protein complexes in the innate immune signal transductions (64, 66). Proteins with TIR domain have also been identified in plants, where they mediate disease resistance (67); however, TIR domain containing proteins in bacteria are involved in virulence (61). In this study, TLR1s of all the 71 examined species shared the TIR domain, implying that TIR may be a functional domain necessary for TLR1 to trigger the immune response in yellow catfish. Effect of the loss of LRR_CT in vertebrate TLR1 is unknown, presumably because genes are accompanied by a certain extent in the process of evolution. Interestingly, previous studies reported that TLRs form homodimers and heterodimers within the membrane, and the single transmembrane domain of these receptors has been involved in the dimerization and corresponding functions (68, 69). The loss of transmembrane domain in TLR1 proteins of some bony fishes may inactivate the formation of dimers in the membrane and further affect its function. We speculate that this evolution may be related to the aquatic environment for fishes’ growth. The small variances in protein spatial structure of TLR1s from different taxa indicated that the function of this protein in various species is relatively conserved (70, 71).
In general, the TLR proteins in prototype metazoans are structurally characterized by three typical domains, including hydrophobic tandem LRR extracellular domain (ECD, mediating PAMP recognition), a short transmembrane (TM) domain, and the intracellular TIR (for signal transmission to downstream pathway components) (72). In our present study, the pattern and strength of natural selection that affects the TLR1s and TLR1-TIRs in both vertebrates and invertebrates were well analyzed. The values of dN/dS (ω) indicate that all the animal TLR1 genes are under purification selection (the value ranges from 0.177 in Amphibian to 0.406 in Insecta for TLR1 genes; Figure 4). A number of functional residues of TLR1 proteins in different species are highly conserved, although they are located directly at the host-environment interface, providing a rigid structural framework for identifying the conserved pathogen-associated molecular patterns (73). The cytoplasmic TIR signal domain induces a signal cascade when the TLR recognizes a specific ligand (74). With average dN/dS (ω) values of TLR1-TIR ranging from 0.100 for Insecta to 0.862 for Bivalvia, it seems that TLR1-TIR exhibited better evolutionary flexibility than TLR1 (Figure 4); that is to say, the differences in the activity levels among different taxa TLR1-TIR domains and downstream molecular components are particularly obvious. The mean values of PSS/non-PSS ranged from 17/1258 in Insecta to 42/780 in Aves (Figure 5), suggesting that TLR1 in Aves is potentially more tolerant for non-synonymous mutations, which may be subjected to positive selection and fixation in Aves (75).
Previous studies showed that tissue distribution analysis helps us to understand the relative expression levels of target genes in various tissues (76, 77), and to determine the detailed biological processes with functional genes. Our results in this study revealed widespread presence of TLR1 transcript in all the 11 collected tissues (Figure 6), indicating that TLR1 is ubiquitous in hybrid yellow catfish. Meanwhile, TLR1 exhibited a similar transcription pattern in Acipenser dabryanus (76). Interestingly, the mRNA levels of TLR1 in immune system were much higher (Figure 6), indicating that TLR1 may mainly function for immunity response in hybrid yellow catfish. Our further data demonstrated that the transcription levels of TLR1 increased significantly after the infection of exogenous A. hydrophila when compared with the control group (Figure 7), revealing that TLR1 participated in bacterial infection in hybrid yellow catfish.
In mammals, apoptosis of the infected cells prevented the spread of microbes throughout the whole body, a strategy that has been maintained throughout the evolution. Detection of specific glycopeptides in the cell wall peptidoglycan triggers the activation of caspase-8-dependent apoptosis, thereby enhancing the clearance of bacteria in the infected cells (78). TLR2, in conjunction with TLR1, plays an important role in the innate immune response by recognizing microbial lipoproteins and lipopeptides in the process of infection (21). Apoptosis signals triggered by TLR2 were functions with MyD88, and this pathway involves FAS-associated death domain protein (FADD) and caspase 8. Furthermore, the binding of MyD88 to FADD is sufficient to induce the cell apoptosis (79). However, TLR signaling pathway consisting of (TLR1-2) - MyD88 - FADD - Caspase 8 has not been systematically studied in teleosts. In this study, TLR1, TLR2, caspase 8, MyD88, FADD and TOLLIP were successfully identified in the yellow catfish, and we proved that there is an only single copy of these genes through the homologous sequence alignment and chromosome location, indicating that TLR signaling pathways are very conserved in the yellow catfish. The expression profiling of TLR signaling pathway genes after pathogen stimulation showed that the infection of pathogens induced the activation of (TLR1-2) - MyD88 - FADD - Caspase 8 involved signaling pathway, and inhibits the spread of bacteria in the body through cell apoptosis.
5 Conclusions
In summary, the TLR1 gene of hybrid yellow catfish was characterized using the comparative genomic survey. The phylogenetic analysis of TLR1 revealed the evolutionary relationship of the selected animals. Our study showed that TLR signaling pathway is very conserved and indispensable in bacterial infection in the hybrid yellow catfish. Meanwhile, our research provides a reference dataset for in-depth studies on the molecular mechanisms of (TLR1-2) - MyD88 - FADD - Caspase 8 mediating apoptosis pathway in bony fish after the bacterial pathogen infection.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding authors.
Ethics statement
The animal study was reviewed and approved by the Ethics Committee of the College of Life Sciences, Sichuan University.
Author contributions
SG designed the experiments. SG, WG and MZ performed the experiments. SG analyzed the data. SG, WG, MZ, QY, FL, LC, Z-YW, YJ and PX contributed reagents/materials tools. SG contributed to the discussion. SG wrote the draft manuscript. ZS, QS and Z-YW supervised the project and revised the manuscript. Z-YW acquired the fund. All authors contributed to the article and approved the submitted version.
Funding
This work was financially supported by the Natural Science Fund of Sichuan Province of China (No. 2023NSFSC1221).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fimmu.2023.1163781/full#supplementary-material
References
1. Geng Y, Liu D, Han S, Zhou Y, Wang KY, Huang XL, et al. Outbreaks of vibriosis associated with Vibrio mimicus in freshwater catfish in China. Aquaculture (2014) 433:82–4. doi: 10.1016/j.aquaculture.2014.05.053
2. Liu F, Shi H-Z, Guo Q-S, Yu Y-B, Wang A-M, Lv F, et al. Effects of astaxanthin and emodin on the growth, stress resistance and disease resistance of yellow catfish (Pelteobagrus fulvidraco). Fish Shellfish Immunol (2016) 51:125–35. doi: 10.1016/j.fsi.2016.02.020
3. Fernandes DC, Eto SF, Funnicelli MIG, Fernandes CC, Charlie-Silva I, Belo MAA, et al. Immunoglobulin y in the diagnosis of Aeromonas hydrophila infection in Nile tilapia (Oreochromis niloticus). Aquaculture (2019) 500:576–85. doi: 10.1016/j.aquaculture.2018.10.045
4. Guo S, Zeng M, Gao W, Li F, Wei X, Shi Q, et al. Toll-like receptor 3 in the hybrid yellow catfish (Pelteobagrus fulvidraco♀× p. vachelli♂): Protein structure, evolution and immune response to exogenous Aeromonas hydrophila and poly (I: C) stimuli. Animals (2023) 13:288. doi: 10.3390/ani13020288
5. Citterio B, Biavasco F. Aeromonas hydrophila virulence. Virulence (2015) 6:417–8. doi: 10.1080/21505594.2015.1058479
6. Zhu W, Zhou S, Chu W. Comparative proteomic analysis of sensitive and multi-drug resistant Aeromonas hydrophila isolated from diseased fish. Microb Pathogen (2020) 139:103930. doi: 10.1016/j.micpath.2019.103930
7. Pang M, Jiang J, Xie X, Wu Y, Dong Y, Kwok AHY, et al. Novel insights into the pathogenicity of epidemic Aeromonas hydrophila ST251 clones from comparative genomics. Sci Rep (2015) 5:9833. doi: 10.1038/srep09833
8. Abdelhamed H, Ibrahim I, Baumgartner W, Lawrence ML, Karsi A. Characterization of histopathological and ultrastructural changes in channel catfish experimentally infected with virulent. Aeromonas hydrophila. Front Microbiol (2017) 8:1519. doi: 10.3389/fmicb.2017.01519
9. AlYahya SA, Ameen F, Al-Niaeem KS, Al-Sa'adi BA, Hadi S, Mostafa AA. Histopathological studies of experimental Aeromonas hydrophila infection in blue tilapia, oreochromis aureus. Saudi J Biol Sci (2018) 25:182–5. doi: 10.1016/j.sjbs.2017.10.019
10. Dickson K, Lehmann C. Inflammatory response to different toxins in experimental sepsis models. Int J Mol Sci (2019) 20:4341. doi: 10.3390/ijms20184341
11. Cabău G, Crișan TO, Klück V, Popp RA, Joosten LAB. Urate-induced immune programming: Consequences for gouty arthritis and hyperuricemia. Immunol Rev (2020) 294:92–105. doi: 10.1111/imr.12833
12. Kieser KJ, Kagan JC. Multi-receptor detection of individual bacterial products by the innate immune system. Nat Rev Immunol (2017) 17:376–90. doi: 10.1038/nri.2017.25
13. Takeda K, Akira S. Toll-like receptors in innate immunity. Int Immunol (2005) 17:1–14. doi: 10.1093/intimm/dxh186
14. Liu H, Guo S, He Y, Shi Q, Yang M, You X. Toll protein family structure, evolution and response of the whiteleg shrimp (Litopenaeus vannamei) to exogenous iridescent virus. J Fish Dis (2021) 44:1131–45. doi: 10.1111/jfd.13374
15. Faure-Dupuy S, Vegna S, Aillot L, Dimier L, Esser K, Broxtermann M, et al. Characterization of pattern recognition receptor expression and functionality in liver primary cells and derived cell lines. J Innate Immun (2018) 10:339–48. doi: 10.1159/000489966
16. Huang Y, Temperley ND, Ren L, Smith J, Li N, Burt DW. Molecular evolution of the vertebrate TLR1 gene family - a complex history of gene duplication, gene conversion, positive selection and co-evolution. BMC Evol Biol (2011) 11:149. doi: 10.1186/1471-2148-11-149
17. Kawai T, Akira S. Toll-like receptor downstream signaling. Arthritis Res Ther (2004) 7:12. doi: 10.1186/ar1469
18. Rock FL, Hardiman G, Timans JC, Kastelein RA, Bazan JF. A family of human receptors structurally related to Drosophila toll. Proc Nat Acad Sci (1998) 95:588–93. doi: 10.1073/pnas.95.2.588
19. O'Neill LAJ. The interleukin-1 receptor/Toll-like receptor superfamily: 10 years of progress. Immunol Rev (2008) 226:10–8. doi: 10.1111/j.1600-065X.2008.00701.x
20. Xing Q, Liao H, Xun X, Wang J, Zhang Z, Yang Z, et al. Genome-wide identification, characterization and expression analyses of TLRs in yesso scallop (Patinopecten yessoensis) provide insight into the disparity of responses to acidifying exposure in bivalves. Fish Shellfish Immunol (2017) 68:280–8. doi: 10.1016/j.fsi.2017.07.020
21. Jin MS, Kim SE, Heo JY, Lee ME, Kim HM, Paik S-G, et al. Crystal structure of the TLR1-TLR2 heterodimer induced by binding of a tri-acylated lipopeptide. Cell (2007) 130:1071–82. doi: 10.1016/j.cell.2007.09.008
22. Kang JY, Nan X, Jin MS, Youn S-J, Ryu YH, Mah S, et al. Recognition of lipopeptide patterns by toll-like receptor 2-toll-like receptor 6 heterodimer. Immunity (2009) 31:873–84. doi: 10.1016/j.immuni.2009.09.018
23. Takeuchi O, Sato S, Horiuchi T, Hoshino K, Takeda K, Dong Z, et al. Cutting edge: role of toll-like receptor 1 in mediating immune response to microbial lipoproteins. J Immunol (2002) 169:10. doi: 10.4049/jimmunol.169.1.10
24. Plantinga TS, Johnson MD, Scott WK, van de Vosse E, Velez Edwards DR, Smith PB, et al. Toll-like receptor 1 polymorphisms increase susceptibility to candidemia. J Infect Dis (2012) 205:934–43. doi: 10.1093/infdis/jir867
25. Wirth AJ, Gruebele M. Quinary protein structure and the consequences of crowding in living cells: Leaving the test-tube behind. BioEssays (2013) 35:984–93. doi: 10.1002/bies.201300080
26. Ikeya T, Güntert P, Ito Y. Protein structure determination in living cells. Int J Mol Sci (2019) 20:2442. doi: 10.3390/ijms20102442
27. Cooper GM, Sidow A. Genomic regulatory regions: Insights from comparative sequence analysis. Curr Opin Genet Dev (2003) 13:604–10. doi: 10.1016/j.gde.2003.10.001
28. Faksri K, Xia E, Ong RT-H, Tan JH, Nonghanphithak D, Makhao N, et al. Comparative whole-genome sequence analysis of v isolated from tuberculous meningitis and pulmonary tuberculosis patients. Sci Rep (2018) 8:4910. doi: 10.1038/s41598-018-23337-y
29. Spielman SJ, Wilke CO. The relationship between dN/dS and scaled selection coefficients. Mol Biol Evol (2015) 32:1097–108. doi: 10.1093/molbev/msv003
30. Guo S, Wen Z, Zhang X, Li F, Cui H, Lin X, et al. Characterization of five caspase genes and their transcriptional changes in response to exogenous iridescent virus challenge in the whiteleg shrimp (Litopenaeus vannamei). Aquaculture (2021) 534:736192. doi: 10.1016/j.aquaculture.2020.736192
31. Wen Z, Li Y, Bian C, Shi Q, Li Y. Genome-wide identification of a novel elovl4 gene and its transcription in response to nutritional and osmotic regulations in rabbitfish (Siganus canaliculatus). Aquaculture (2020) 529:735666. doi: 10.1016/j.aquaculture.2020.735666
32. Han P, Wang S, Zhang Q, Zhang S, Shao R, Xu W, et al. Molecular characterization and expression analysis of TLR1 and TLR4 from the endangered fish dabry's sturgeon (Acipenser dabryanus). Dev Comp Immunol (2018) 86:180–8. doi: 10.1016/j.dci.2018.05.009
33. Salazar C, Haussmann D, Kausel G, Figueroa J. Molecular cloning of Salmo salar toll-like receptors (TLR1, TLR22, TLR5M and TLR5S) and expression analysis in SHK-1 cells during Piscirickettsia salmonis infection. J Fish Dis (2016) 39:239–48. doi: 10.1111/jfd.12354
34. Wang KL, Chen SN, Huo HJ, Nie P. Identification and expression analysis of sixteen toll-like receptor genes, TLR1, TLR2a, TLR2b, TLR3, TLR5M, TLR5S, TLR7–9, TLR13a–c, TLR14, TLR21–23 in mandarin fish siniperca chuatsi. Dev Comp Immunol (2021) 121:104100. doi: 10.1016/j.dci.2021.104100
35. Wei YC, Pan TS, Chang MX, Huang B, Xu Z, Luo TR, et al. Cloning and expression of toll-like receptors 1 and 2 from a teleost fish, the orange-spotted grouper. Epinephelus coioides. Vet Immunol Immunop (2011) 141:173–82. doi: 10.1016/j.vetimm.2011.02.016
36. Gong G, Dan C, Xiao S, Guo W, Huang P, Xiong Y, et al. Chromosomal-level assembly of yellow catfish genome using third-generation DNA sequencing and Hi-c analysis. GigaScience (2018) 7:giy120. doi: 10.1093/gigascience/giy120
37. Telles GP, Araújo GS, Walter MEMT, Brigido MM, Almeida NF. Live neighbor-joining. BMC Bioinf (2018) 19:172. doi: 10.1186/s12859-018-2162-x
38. Hoang DT, Chernomor O, Von Haeseler A, Minh BQ, Vinh LS. UFBoot2: improving the ultrafast bootstrap approximation. Mol Biol Evol (2018) 35:518–22. doi: 10.1093/molbev/msx281
39. El-Gebali S, Mistry J, Bateman A, Eddy SR, Luciani A, Potter SC, et al. The pfam protein families database in 2019. Nucleic Acids Res (2019) 47:D427–D32. doi: 10.1093/nar/gky995
40. Letunic I, Khedkar S, Bork P. SMART: recent updates, new developments and status in 2020. Nucleic Acids Res (2021) 49:D458–D60. doi: 10.1093/nar/gkaa937
41. Yuan S, Chan HCS, Hu Z. Using PyMOL as a platform for computational drug design. WIREs Comput Mol Sci (2017) 7:e1298. doi: 10.1002/wcms.1298
42. Yang Q, Wei C, Guo S, Liu J, Tao Y. Cloning and characterization of al-lactate dehydrogenase gene from ruminococcaceae bacterium CPB6. World J Microbiol Biotechnol (2020) 36:1–10. doi: 10.1007/s11274-020-02958-4
43. Fay JC, Wu C-I. Sequence divergence, functional constraint, and selection in protein evolution. Annu Rev Genomics Hum Genet (2003) 4:213–35. doi: 10.1146/annurev.genom.4.020303.162528
44. Zhang Z, Li J, Zhao X-Q, Wang J, Wong GK-S, Yu J. KaKs_Calculator: calculating ka and ks through model selection and model averaging. Genom Proteom Bioinf (2006) 4:259–63. doi: 10.1016/S1672-0229(07)60007-2
45. Johnson LS, Eddy SR, Portugaly E. Hidden Markov model speed heuristic and iterative HMM search procedure. BMC Bioinf (2010) 11:431. doi: 10.1186/1471-2105-11-431
46. Edgar RC. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res (2004) 32:1792–7. doi: 10.1093/nar/gkh340
47. Yang Z. PAML 4: phylogenetic analysis by maximum likelihood. Mol Biol Evol (2007) 24:1586–91. doi: 10.1093/molbev/msm088
48. Delport W, Poon AFY, Frost SDW, Kosakovsky Pond SL. Datamonkey 2010: a suite of phylogenetic analysis tools for evolutionary biology. Bioinformatics (2010) 26:2455–7. doi: 10.1093/bioinformatics/btq429
49. Pond SLK, Frost SDW. Datamonkey: rapid detection of selective pressure on individual sites of codon alignments. Bioinformatics (2005) 21:2531–3. doi: 10.1093/bioinformatics/bti320
50. Berkman SJ, Roscoe EM, Bourret JC. Comparing self-directed methods for training staff to create graphs using graphpad prism. J Appl Behav Anal (2019) 52:188–204. doi: 10.1002/jaba.522
51. Li X, Hamyat M, Liu C, Ahmad S, Gao X, Guo C, et al. Identification and characterization of the WOX family genes in five solanaceae species reveal their conserved roles in peptide signaling. Genes (2018) 9:260. doi: 10.3390/genes9050260
52. dos Reis M, Thawornwattana Y, Angelis K, Telford Maximilian J, Donoghue Philip CJ, Yang Z. Uncertainty in the timing of origin of animals and the limits of precision in molecular timescales. Curr Biol (2015) 25:2939–50. doi: 10.1016/j.cub.2015.09.066
53. Inoue JG, Miya M, Lam K, Tay B-H, Danks JA, Bell J, et al. Evolutionary origin and phylogeny of the modern holocephalans (Chondrichthyes: Chimaeriformes): A mitogenomic perspective. Mol Biol Evol (2010) 27:2576–86. doi: 10.1093/molbev/msq147
54. Shi SY, Lu Y-W, Richardson J, Min X, Weiszmann J, Richards WG, et al. A systematic dissection of sequence elements determining β-klotho and FGF interaction and signaling. Sci Rep (2018) 8:11045. doi: 10.1038/s41598-018-29396-5
55. Scaglione KM, Basrur V, Ashraf NS, Konen JR, Elenitoba-Johnson KSJ, Todi SV, et al. The ubiquitin-conjugating enzyme (E2) Ube2w ubiquitinates the n terminus of substrates. J Biol Chem (2013) 288:18784–8. doi: 10.1074/jbc.C113.477596
56. Arora H, Wilcox SM, Johnson LA, Munro L, Eyford BA, Pfeifer CG, et al. The ATP-binding cassette gene ABCF1 functions as an E2 ubiquitin-conjugating enzyme controlling macrophage polarization to dampen lethal septic shock. Immunity (2019) 50:418–31.e6. doi: 10.1016/j.immuni.2019.01.014
57. Ling J, Cheloha RW, McCaul N, Sun Z-YJ, Wagner G, Ploegh HL. A nanobody that recognizes a 14-residue peptide epitope in the E2 ubiquitin-conjugating enzyme UBC6e modulates its activity. Mol Immunol (2019) 114:513–23. doi: 10.1016/j.molimm.2019.08.008
58. Ruff KM, Pappu RV. AlphaFold and implications for intrinsically disordered proteins. J Mol Biol (2021) 433:167208. doi: 10.1016/j.jmb.2021.167208
59. Jumper J, Evans R, Pritzel A, Green T, Figurnov M, Ronneberger O, et al. Highly accurate protein structure prediction with AlphaFold. Nature (2021) 596:583–9. doi: 10.1038/s41586-021-03819-2
60. Streit RL, Luginbuhl TE. Maximum likelihood method for probabilistic multihypothesis tracking. Orlando FL, United States: SPIE's International Symposium on Optical Engineering and Photonics in Aerospace Sensing (1994) 2235:394–405. doi: 10.1117/12.179066
61. Waldhuber A, Snyder GA, Römmler F, Cirl C, Müller T, Xiao TS, et al. A comparative analysis of the mechanism of toll-like receptor-disruption by tir-containing protein c from uropathogenic escherichia coli. Pathogens (2016) 5:25. doi: 10.3390/pathogens5010025
62. Gascuel O, Steel M. Neighbor-joining revealed. Mol Biol Evol (2006) 23:1997–2000. doi: 10.1093/molbev/msl072
63. Sheneman L, Evans J, Foster JA. Clearcut: a fast implementation of relaxed neighbor joining. Bioinformatics (2006) 22:2823–4. doi: 10.1093/bioinformatics/btl478
64. Essuman K, Summers DW, Sasaki Y, Mao X, Yim AKY, DiAntonio A, et al. TIR domain proteins are an ancient family of NAD+-consuming enzymes. Curr Biol (2018) 28:421–30.e4. doi: 10.1016/j.cub.2017.12.024
65. Hawn TR, Misch EA, Dunstan SJ, Thwaites GE, Lan NTN, Quy HT, et al. A common human TLR1 polymorphism regulates the innate immune response to lipopeptides. Eur J Immunol (2007) 37:2280–9. doi: 10.1002/eji.200737034
66. O'Neill LAJ, Golenbock D, Bowie AG. The history of toll-like receptors — redefining innate immunity. Nat Rev Immunol (2013) 13:453–60. doi: 10.1038/nri3446
67. Askarian F, van Sorge NM, Sangvik M, Beasley FC, Henriksen JR, Sollid JUE, et al. A Staphylococcus aureus TIR domain protein virulence factor blocks TLR2-mediated NF-κB signaling. J Innate Immun (2014) 6:485–98. doi: 10.1159/000357618
68. Godfroy JI III, Roostan M, Moroz YS, Korendovych IV, Yin H. Isolated toll-like receptor transmembrane domains are capable of oligomerization. PLoS One (2012) 7:e48875. doi: 10.1371/journal.pone.0048875
69. Reuven EM, Fink A, Shai Y. Regulation of innate immune responses by transmembrane interactions: Lessons from the TLR family. Biochim Biophys Acta (2014) 1838:1586–93. doi: 10.1016/j.bbamem.2014.01.020
70. Mahley RW, Weisgraber KH, Huang Y. Apolipoprotein e: structure determines function, from atherosclerosis to alzheimer's disease to AIDS. J Lipid Res (2009) 50:S183–S8. doi: 10.1194/jlr.R800069-JLR200
71. Senior AW, Evans R, Jumper J, Kirkpatrick J, Sifre L, Green T, et al. Improved protein structure prediction using potentials from deep learning. Nature (2020) 577:706–10. doi: 10.1038/s41586-019-1923-7
72. Kawai T, Akira S. The role of pattern-recognition receptors in innate immunity: Update on toll-like receptors. Nat Commun (2010) 11:373–84. doi: 10.1038/ni.1863
73. Werling D, Jann OC, Offord V, Glass EJ, Coffey TJ. Variation matters: TLR structure and species-specific pathogen recognition. Trends Immunol (2009) 30:124–30. doi: 10.1016/j.it.2008.12.001
74. Fekonja O, Benčina M, Jerala R. Toll/interleukin-1 receptor domain dimers as the platform for activation and enhanced inhibition of toll-like receptor signaling. J Biol Chem (2012) 287:30993–1002. doi: 10.1074/jbc.M112.376186
75. Barreiro LB, Ben-Ali M, Quach H, Laval G, Patin E, Pickrell JK, et al. Evolutionary dynamics of human toll-like receptors and their different contributions to host defense. PLoS Genet (2009) 5:e1000562. doi: 10.1371/journal.pgen.1000562
76. Ramsköld D, Wang ET, Burge CB, Sandberg R. An abundance of ubiquitously expressed genes revealed by tissue transcriptome sequence data. PLoS Comput Biol (2009) 5:e1000598. doi: 10.1371/journal.pcbi.1000598
77. Zhou Q, Liu M, Xia X, Gong T, Feng J, Liu W, et al. A mouse tissue transcription factor atlas. Nat Commun (2017) 8:15089. doi: 10.1038/ncomms15089
78. Salaun B, Romero P, Lebecque S. Toll-like receptors’ two-edged sword: When immunity meets apoptosis. Eur J Immunol (2007) 37:3311–8. doi: 10.1002/eji.200737744
Keywords: hybrid yellow catfish, innate immunity, Aeromonas hydrophila, TLR signaling pathway, expression profile
Citation: Guo S, Gao W, Zeng M, Liu F, Yang Q, Chen L, Wang Z, Jin Y, Xiang P, Chen H, Wen Z, Shi Q and Song Z (2023) Characterization of TLR1 and expression profiling of TLR signaling pathway related genes in response to Aeromonas hydrophila challenge in hybrid yellow catfish (Pelteobagrus fulvidraco ♀ × P. vachelli ♂). Front. Immunol. 14:1163781. doi: 10.3389/fimmu.2023.1163781
Received: 11 February 2023; Accepted: 17 March 2023;
Published: 28 March 2023.
Edited by:
Akira Goto, (INSERM), FranceReviewed by:
Cristobal Espinosa Ruiz, University of Murcia, SpainSilvia Leonor Mundo, University of Buenos Aires, Argentina
Copyright © 2023 Guo, Gao, Zeng, Liu, Yang, Chen, Wang, Jin, Xiang, Chen, Wen, Shi and Song. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Zhengyong Wen, zhengyong_wen@126.com; Qiong Shi, shiqiong@genomics.cn; Zhaobin Song, zbsong@scu.edu.cn