SNP detection of TLR8 gene, association study with susceptibility/resistance to GCRV and regulation on mRNA expression in grass carp, Ctenopharyngodon idella

Comparative Transcriptomic Analysis Revealed Potential Differential Mechanisms of Grass Carp Reovirus Pathogenicity

Grass carp reovirus (GCRV), one of the most serious pathogens threatening grass carp (Ctenopharyngodon idella), can lead to grass carp hemorrhagic disease (GCHD). Currently, GCRV can be divided into three genotypes, but the comparison of their pathogenic mechanisms and the host responses remain unclear. In this study, we utilized the Ctenopharyngodon idella kidney (CIK) model infected with GCRV to conduct comparative studies on the three genotypes. We observed a cytopathic effect (CPE) in the GCRV-I and GCRV-III groups, whereas the GCRV-II group did not show any CPE. Moreover, a consistent trend in the mRNA expression levels of antiviral-related genes across all experimental groups of CIK cells was detected via qPCR and further explored through RNA-seq analysis. Importantly, GO/KEGG enrichment analysis showed that GCRV-I, -II, and -III could all activate the immune response in CIK cells, but GCRV-II induced more intense immune responses. Intriguingly, transcriptomic analysis revealed a widespread down-regulation of metabolism processes such as steroid biosynthesis, butanoate metabolism, and N-Glycan biosynthesis in infected CIK cells. Overall, our results reveal the CIK cells showed unique responses in immunity and metabolism in the three genotypes of GCRV infection. These results provide a theoretical basis for understanding the pathogenesis and prevention and control methods of GCRV.

Novel single nucleotide polymorphisms of insulin-like growth factor-binding protein 7 (IGFBP7) gene significantly associated with growth traits in striped catfish (Pangasianodon hypophthalmus Sauvage, 1878)

Breeding program to improve economically important growth traits in striped catfish (Pangasianodon hypophthalmus) requires effective molecular markers. This study was conducted to identify single nucleotide polymorphisms (SNPs) of Insulin-like Growth Factor-Binding Protein 7 (IGFBP7) gene which plays multiple roles in regulating growth, energy metabolism and development. The association between SNPs in IGFBP7 gene and growth traits in striped catfish was analyzed in order to uncover the SNPs that have potential to be valuable markers for improving growth traits. Firstly, fragments of IGFBP7 gene from ten fast-growing fish and ten slow-growing fish were sequenced in order to discover SNPs. After filtering the detected SNPs, an intronic SNP (2060A > G) and two non-synonymous SNPs (344 T > C and 4559C > A) causing Leu78Pro and Leu189Met in protein, respectively, were subjected to further validated by individual genotyping in 70 fast-growing fish and 70 slow-growing fish using single base extension method. Our results showed that two SNPs (2060A > G and 4559 C > A (p. Leu189Met)) were significantly associated with the growth in P. hypophthalmus (p < 0.001), thus being candidate SNP markers for the growth traits of this fish. Moreover, linkage disequilibrium and association analysis with growth traits of haplotypes generated from the 3 filtered SNPs (344 T > C, 2060 A > G and 4559 C > A) were examined. These revealed that the non-coding SNP locus (2060A > G) had higher genetic diversity at which the G allele was predominant over the A allele in the fast-growing fish. Furthermore, the results of qPCR showed that expression of IGFBP7 gene with genotype GG (at locus 2060) in fast-growing group was significantly higher than that with genotype AA in slow-growing group (p < 0.05). Our study provides insights into the genetic variants of IGFBP7 gene and useful data source for development molecular marker for growth traits in breeding of the striped catfish.

Characterization of mandarin fish (Siniperca chuatsi) IL-6 and IL-6 signal transducer and the association between their SNPs and resistance to ISKNV disease

In fish, interleukin-6 (IL-6) is a very important immune-regulatory cytokine that plays a polyfunctional role in inflammation, metabolism, regeneration, and neural processes. IL-6 signal transducer (IL-6ST) is a specific receptor for IL-6 and expressed mainly in immune cells and hepatocytes. In this study, the complete cDNA and genomic DNA sequences of mandarin fish (Siniperca chuatsi) IL-6 and IL-6ST genes were identified and analyzed. Quantitative real-time PCR showed that IL-6 and IL-6ST were chiefly expressed in the immune organs. After challenge with infectious spleen and kidney necrosis virus (ISKNV), the expression levels of IL-6 were significantly up-regulated after 6 h and 24 h in the head kidney and spleen, respectively (p < 0.01), the peak value for both reached at 72 h, IL-6ST increased significantly after 120 h with a peak at 168 h in the head kidney (p < 0.01) and improved markedly at 168 h in the spleen (p < 0.01). Besides, IL-6 and IL-6ST have been identified 3 and 8 single nucleotide polymorphisms (SNPs), respectively. Statistical analysis showed that one SNP locus (1625C/T) in the coding region of IL-6 was significantly related to the resistance of mandarin fish against ISKNV. The 1625C→T locus in the coding region of IL-6 is a synonymous mutation; compared with the susceptible group, the frequency of allele T in the disease resistance group was significantly higher, which may be due to the rare codon produced by the mutation affecting translation. The involvement of IL-6 and IL-6ST in response to ISKNV infection in mandarin fish clearly indicate that the role of SNP markers in IL-6 was associated with the ISKNV resistance, which was demonstrated for the first time in our results. Thus, the current study may provide fundamental information for further breeding of mandarin fish with resistance to ISKNV infection.

R848 Is Involved in the Antibacterial Immune Response of Golden Pompano (Trachinotus ovatus) Through TLR7/8-MyD88-NF-κB-Signaling Pathway

R848 is an imidazoquinoline compound that is a specific activator of toll-like receptor (TLR) 7/8 and is often used in immunological research in mammals and teleosts. However, the immune responses initiated by R848 through the TLR7/8 pathway in response to bacterial infection remain largely unexplored in teleosts. In the current study, we investigated the antibacterial response and the participating signaling pathway initiated by R848 in golden pompano (Trachinotus ovatus). We found that R848 could stimulate the proliferation of head kidney lymphocytes (HKLs) in a dose-dependent manner, enhance the survival rate of HKLs, and inhibit the replication of bacteria in vivo. However, these effects induced by R848 were significantly reduced when chloroquine (CQ) was used to blocked endosomal acidification. Additionally, an in vivo study showed that R848 strengthened the antibacterial immunity of fish through a TLR7/8 and Myd88-dependent signaling pathway. A cellular experiment showed that Pepinh-MYD (a Myd88 inhibitor) significantly reduced the R848-mediated proliferation and survival of HKLs. Luciferase activity analysis showed that R848 enhanced the nuclear factor kappa B (NF-κB) activity, whereas this activity was reduced when CQ and Pepinh-MYD were present. Additionally, when an NF-κB inhibitor was present, the R848-mediated pro-proliferative and pro-survival effects on HKLs were significantly diminished. An in vivo study showed that knockdown of TLR7, TLR8, and Myd88 expression in golden pompano via siRNA following injection of R848 resulted in increased bacterial dissemination and colonization in fish tissues compared to that of fish injection of R848 alone, suggesting that R848-induced antibacterial immunity was significantly reduced. In conclusion, these results indicate that R848 plays an essential role in the antibacterial immunity of golden pompano via the TLR7/8-Myd88-NF-κB- signaling pathway.

Beta-adrenergic receptor stimulation influences the function of monocytes/macrophages in ayu (Plecoglossus altivelis)

Adrenergic receptors (ARs) are members of the G-protein-coupled receptor superfamily that can be categorized into αARs and βARs. The specific function of ARs in teleost monocytes/macrophages (MO/MФ) remains unknown. We determined the cDNA sequence of ARs from ayu (Plecoglossus altivelis; PaαAR and PaβAR). Sequence comparisons showed that PaαAR was most closely related to the αAR of the Japanese flounder and Nile tilapia, while PaβAR was most closely related to the βAR of Atlantic salmon. The AR transcripts were mainly expressed in the spleen, and their expression was altered in various tissues upon infection with Vibrio anguillarum. PaαAR and PaβAR proteins were upregulated in MO/MФ after infection, and PaβAR knockdown resulted in a pro-inflammatory status in ayu MO/MФ upon V. anguillarum infection and lowered the phagocytic activity of MO/MФ. Our results indicate that PaβAR plays the role of an anti-inflammatory mediator in the immune response of ayu against bacterial infection.

Different responses in one-year-old and three-year-old grass carp reveal the mechanism of age restriction of GCRV infection

Grass carp is an important fish species in Chinese aquaculture, and can be afflicted by a hemorrhagic disease caused by the grass carp reovirus (GCRV). Interestingly, the affects of GCRV infection of grass carp are age-restricted, meaning that one-year-old grass carp can be infected and can suffer hemorrhagic disease, but three-year-old carp are not so afflicted. In this study, we investigated the mechanism responsible for this age-restricted pathology. We evaluated the relative copy number of GCRV RNA, the expression levels of proteins in blood, and changes in DNA methylation in carp from the two age groups after infection with GCRV. After GCRV infection, the relative copy number of GCRV RNA in three-year-old grass carp was significantly lower than in one-year-old carp. The differences in circulating protein levels mainly occurred in concentrated in complement and coagulation proteins, and the expression levels of these proteins were significantly higher in three-year-old grass carp than in one-year-old carp. Moreover, the expression levels of DNA methylation-related genes in the liver and spleen of one-year-old grass carp were significantly higher than those of three-year-old carp. These results suggested that as age of grass carp increases, faster and more efficient response of the immune system after viral infection, especially the complement system, and differences in DNA methylation may be important factors that affect the age restriction observed in GCRV infection. Our study provides new insights into the mechanisms underlying age restriction of GCRV infection.

Identification and immune-related analysis of SNPs in Litopenaeus vannamei Toll3 receptor

Tolls and Toll-like receptors (TLRs), as innate immune-recognition receptors that recognize molecular patterns associated with microbial pathogens, play a critical role in antimicrobial immune responses. Here, we report on single nucleotide polymorphisms (SNPs) of Litopenaeus vannamei Toll3 (LvToll3). Multiple sequence alignment of the L. vannamei Toll3 Leucine rich repeat C-terminal domain (LvToll3-LRR-CT) with other L. vannamei Tolls LRR-CT domains showed 39.23% - 43.96% homology at the nucleic acid level and 20.31% - 30.00% identity at the amino acid level. Analysis of different shrimp tissues by polymerase chain reaction denaturing gradient gel electrophoresis (PCR-DGGE) revealed that LvToll3-LRR-CT had genetic polymorphisms at both the genomic deoxyribonucleic acid (gDNA) and complementary deoxyribonucleic acid (cDNA) levels. Further, high-throughput sequencing analysis confirmed the presence of 8 non-synonymous SNP (nsSNP) and 1 nsSNPs with frequency greater than 1% at the gDNA level, while 13 nsSNPs and 2 nsSNPs with frequency greater than 1% at the cDNA level. In silico analysis revealed that the α-helix secondary structure and tertiary structure of LvToll3 changed when 3 SNPs (C2039T, T2041C, T2228C) were mutated. Interestingly, 2 novel bands on PCR-DGGE, which were identified as 2 nsSNPs (C2140A, T2186A) were observed following challenge with Streptococcus iniae but not with Vibrio parahaemolyticus or White spot syndrome virus (WSSV). Moreover, the secondary and tertiary structures of LvToll3 changed when the nsSNP T2186A was mutated. The present findings therefore provide novel insight into the molecular basis of shrimp innate immune response to pathogens through the generation of specific SNPs.

Molecular characterization and expression of TLR7 and TLR8 in barbel chub (Squaliobarbus curriculus): Responses to stimulation of grass carp reovirus and lipopolysaccharide

The barbel chub (Squaliobarbus curriculus) is a kind of small size commercial fish species that is widely spread in Asia and has shown significant resistance to disease. In this study, the full-length cDNA sequences of Toll-like receptor (TLR) 7 and 8 from S. curriculus, designated as ScTLR7 and ScTLR8, were cloned, and their differences in the structure and the responses to the grass carp (GCRV) infection and lipopolysaccharide stimulation were investigated. The full-length 3715 base pair (bp) cDNA of ScTLR7 contained a complete open reading frame of 3162 bp and encoded a putative polypeptide of 1053 amino acid residues. The full-length 4624 base pair (bp) cDNA of ScTLR8 contained a complete open reading frame of 3072 bp and encoded a putative polypeptide of 1023 amino acid residues. ScTLR7 and ScTLR8 consisted of N-terminal signal peptide, leucine-rich repeats (LRRs), and Toll/IL-1 Receptors domain. LRR motifs of ScTLR7 and ScTLR8 bend into horseshoe-like solenoid structure, while the number of LRRs between the two genes is different. Phylogenetic analysis showed that both the ScTLR7 and ScTLR8 were closely clustered with Ctenopharyngodon idellus and Megalobrama amblycephala. Quantitative real-time polymerase chain reaction analysis showed that the expression levels of ScTLR7 in S. curriculus were most abundant in the brain followed by the spleen and heart, and the lowest in the intestine. The highest expression level of ScTLR8 was observed in spleen and the lowest in the liver. After LPS stimulation, the relative expression levels of both ScTLR7 and ScTLR8 exhibited an overall trend of up-regulation. The expression levels of type I-IFN showed an overall trend of down-regulation at time points of 12, 24, 72 and 168 h compared to that of 6 h after LPS stimulation. Compared to 6 h post GCRV infection, the transcription level of ScTLR7 was up-regulated from 12 to 168 h, and transcription levels of ScTLR8, MyD88, and type I-IFN were firstly up-regulated from 12 to 72 h, and then down-regulated from 72 to 168 h. Correlation analysis showed that expression level of ScTLR7 in the spleen was significantly positively correlated with that of MyD88 (Pearson correlation coefficient: 0.909, P: 0.033), and a significantly positive correlation was also observed between expression levels of MyD88 and type I IFN (Pearson correlation coefficient: 0.962, P: 0.009), after GCRV stimulation. These results indicate that ScTLR7 and ScTLR8 may play important roles in the responses to the grass carp (GCRV) infection and lipopolysaccharide stimulation and trigger different downstream immune signal pathways.

SNP-based susceptibility-resistance association and mRNA expression regulation analyses of tlr7 to grass carp Ctenopharyngodon idella reovirus

Eleven single nucleotide polymorphisms (SNP) in Ctenopharyngodon idella toll-like receptor 7 (citlr7) gene, containing two in the 5'-flanking region, three within the single intron and six distributed in the coding sequence (CDS), were identified. A case-control study of 73 susceptible individuals and 67 resistant individuals was conducted to test the SNPs-based susceptibility-resistance association and mRNA expression of citlr7 to grass carp reovirus (GCRV), showing that both 820 A/G and 1726 A/G were significantly correlative sites in genotype (P < 0·05). Multifactor dimensionality reduction (MDR) analysis suggested the exertion of antiviral effects of 820 A/G might rely on SNPs interactions of citlr7 and C. idella toll-like receptor 8 (citlr8). Combining the mortality rate and citlr7 mRNA expression, it was suggested that 1726 GG-genotyped individuals might be more resistant than 1726 A/G genotyped individuals, indicating the selection on synonymous mutations in 1726 A/G might be susceptibility-resistance-type specific. In addition, haplotype analysis uncovered no significantly correlative haplotypes in citlr7. These findings may provide an in-depth insight for the further functional research of citlr7. The potential genetic markers identified may contribute to the molecular and transgenic breeding of C. idella.

MicroRNA regulation of Toll-like receptor signaling pathways in teleost fish

The innate immune system is the first line defense mechanism that recognizes, responds to, controls or eliminates invading pathogens. Toll-like receptors (TLRs) are a critical family of pattern recognition receptors (PRRs) tightly regulated by complex mechanisms involving many molecules to ensure a beneficial outcome in response to foreign invaders. MicroRNAs (miRNAs), a transcriptional and posttranscriptional regulator family in a wide range of biological processes, have been identified as new molecules related to the regulation of TLR-signaling pathways in immune responses. To date, at least 22 TLR types have been identified in more than a dozen different fish species. However, the functions and underlying mechanisms of miRNAs in the regulation of inflammatory responses related to the TLR-signaling pathway in fish is lacking. In this review, we summarize the regulation of miRNA expression profiles in the presence of TLR ligands or pathogen infections in teleost fish. We focus on the effects of miRNAs in regulating TLR-signaling pathways by targeting multiple molecules, including TLRs themselves, TLR-associated signaling proteins, and TLR-induced cytokines. An understanding of the relationship between the TLR-signaling pathways and miRNAs may provide new insights for drug intervention to manipulate immune responses in fish.

Basal polarization of the immune responses to Streptococcus agalactiae susceptible and resistant tilapia (Oreochromis niloticus)

One of the highest priority areas for improvement is the development of effective strategies for decreasing disease mortality levels in aquaculture production, a better understanding of the components of the fish immune system and their functions in the context of pathogen invasion is needed. Tilapia is the most common fish in South China, and Streptococcus agalactiae has become the most serious disease problem for tilapia industry in China. Here, we profiled gene expression differences between tilapia differing in their susceptibility to S. agalactiae both basally (before infection) and at three early timepoints post-infection (5 h, 50 h, and 7 d). Between group comparisons revealed 5756 unique genes differentially expressed greater than 2-fold at one or more timepoints. And the resistant fish showed much more strong ability in pathogen recognition, antigen presentation, immune activation, while the susceptible fish showed fast activation of apoptosis. Taken together, the immune profiles expand our knowledge for molecular mechanisms for disease resistance, as well as provide solid molecular resources for further identification of the candidate markers for disease-resistant selection and evaluation of disease prevention and treatment options for tilapia industry.

SNPs in the Toll1 receptor of Litopenaeus vannamei are associated with immune response

Tolls and Toll-like receptors (TLRs) are important regulators in the innate immune system and their genetic variations usually affect the host's susceptibility/resistance to pathogen infections. In this study, we report on the single nucleotide polymorphisms (SNPs) of Toll1 in Litopenaeus vannamei (LvToll1) and how this is associated with immune response. PCR-DGGE analysis revealed genetic polymorphisms in LvToll1 at both the genomic DNA (gDNA) and cDNA levels. Using high-throughput sequencing, 223 SNPs were identified at the gDNA level, of which 145 were non-synonymous SNP (nsSNP), with 3 nsSNPs having frequency over 1%. On the other hand, 60 SNPs were identified at the cDNA level including 38 nsSNPs and 4 nsSNPs with frequency over 1%. Upon challenging shrimps with Streptococcus iniae, Vibrio parahaemolyticus and white spot syndrome virus (WSSV), LvToll1 was shown to generate 6, 4 and 4 novel bands, respectively when analyzed with PCR-DGGE. Sequencing analysis of these bands showed that they contained 6, 4 and 2 nsSNPs, respectively. Moreover, the nsSNP C1526T was detected in S. iniae-resistant but not in susceptible shrimps. Most significantly, the C1526T mutation could shorten the α-helix of the LRR domain and was predicted to affect the function of LvToll1, indicating that SNP C1526T might be associated with shrimp's resistance to pathogen infections. In sum, our findings here reveal that the genetic polymorphisms of Toll receptor are linked with the immune response to pathogen infections in L. vannamei.

Multiple functional SNPs in differentially expressed genes modify risk and survival of non-small cell lung cancer in Chinese female non-smokers

DNA genotype can affect gene expression, and gene expression can influence the onset and progression of diseases. Here we conducted a comprehensive study, we integrated analysis of gene expression profile and single nucleotide polymorphism (SNP) microarray data in order to scan out the critical genetic changes that participate in the onset and development of non-small cell lung cancer (NSCLC). Gene expression profile datasets were downloaded from the GEO database. Firstly, differentially expressed genes (DEGs) between NSCLC samples and adjacent normal samples were identified. Next, by STRING database, protein-protein interaction (PPI) network was constructed. At the same time, hub genes in PPI network were identified. Then, some functional SNPs in hub genes that may affect gene expression have been annotated. Finally, we carried a study to explore the relationship between functional SNPs and NSCLC risk and overall survival in Chinese female non-smokers. A total of 488 DEGs were identified in our study. There are 29 proteins with a higher degree of connectivity in the PPI network, including FOS, IL6 and MMP9. By using database annotation, we got 8 candidate functional SNPs that may affect the expression level of hub proteins. In the case-control study, we found that rs4754-T allele, rs959173-C allele and rs2239144-G allele were the protective allele of NSCLC risk. In dominant model, rs4754-CT+TT genotype were associated with a shorter survival time. In general, our study provides a novel research direction in the field of multi-omic data integration, and helps us find some critical genetic changes in disease.

Transcriptome profiling of developing spleen tissue and discovery of immune-related genes in grass carp (Ctenopharyngodon idella)

Grass carp Ctenopharyngodon idella is an important freshwater aquaculture species. However, studies regarding transcriptomic profiling of developing spleen tissue in the grass carp are lacking. Here, the transcriptome sequencing from the spleen tissue of one-year-old (cis1) and three-year-old (cis3) grass carp was performed using Illumina paired-end sequencing technology. The de novo assemblies yielded 48,970 unigenes with average lengths of 1264.51 bp from the two libraries. The assembled unigenes were evaluated and functionally annotated by comparing with sequences in major public databases including Nr, COG, Swiss-Prot, KEGG, Pfam and GO. Comparative analysis of expression levels revealed that a total of 38,254 unigenes were expressed in both the cis1 and cis3 libraries, while 4356 unigenes were expressed only in the cis1 library, and 3312 unigenes were expressed only in the cis3 library. Meanwhile, 1782 unigenes (including 903 down-regulated and 879 up-regulated unigenes) were differentially expressed between the two developmental stages of the grass carp spleen. Based on GO and KEGG enrichment analysis, these differentially expressed genes widely participated in the regulation of immunity and response in the grass carp. Moreover, the main components of six immune-related pathways were identified, including complement and coagulation cascades, Toll-like receptor signaling, B-cell receptor signaling, T-cell receptor signaling, antigen processing and presentation, and chemokine signaling. Finally, two identified transcripts including TLR 8 and complement component C8 were validated for reliability by RT-PCR. Collectively, the results obtained in this study will provide a basis for the study of molecular mechanisms in grass carp spleen development.

Identification and expression analysis of toll-like receptor genes (TLR8 and TLR9) in mucosal tissues of turbot (Scophthalmus maximus L.) following bacterial challenge

Mucosal immune system is one of the most important components in the innate immunity and constitutes the front line of host defense against infection, especially for teleost, which are living in the pathogen-rich aquatic environment. The pathogen recognition receptors (PRRs), which can recognize the conserved pathogen-associated molecular patterns (PAMPs) of bacteria, are considered as one of the most important component for pathogen recognition and immune signaling pathways activation in mucosal immunity. In this regard, we sought to identify TLR8 and TLR9 in turbot (Scophthalmus maximus), as well as their mucosal expression patterns following different bacterial infection in mucosal tissues for the first time. The full-length TLR8 transcript consists of an open reading frame (ORF) of 3108 bp encoding the putative peptide of 1035 amino acids. While the TLR9 was 6730 bp long, containing a 3168 bp ORF that encodes 1055 amino acids. The phylogenetic analysis revealed both TLR8 and TLR9 showed the closest relationship to large yellow croaker. Moreover, both TLR8 and TLR9 could be detected in all examined healthy turbot tissues, with the lowest expression level in liver and a relatively moderate expression pattern in healthy mucosal tissues. Distinct expression patterns of TLR8 and TLR9 were comparatively observed in the mucosal tissues (intestine, gill and skin) following Vibrio anguillarum and Streptococcus iniae infection, suggesting their different roles for mucosal immunity. Further functional studies are needed to better characterize TLR8 and TLR9 and their family members, to better understand the ligand specificity and to identify their roles in different mucosal tissues in protecting fish from the pathogenically hostile environment.

Sensors of Infection: Viral Nucleic Acid PRRs in Fish

Viruses produce nucleic acids during their replication, either during genomic replication or transcription. These nucleic acids are present in the cytoplasm or endosome of an infected cell, or in the extracellular space to be sensed by neighboring cells during lytic infections. Cells have mechanisms of sensing virus-generated nucleic acids; these nucleic acids act as flags to the cell, indicating an infection requiring defense mechanisms. The viral nucleic acids are called pathogen-associated molecular patterns (PAMPs) and the sensors that bind them are called pattern recognition receptors (PRRs). This review article focuses on the most recent findings regarding nucleic acids PRRs in fish, including: Toll-like receptors (TLRs), RIG-I-like receptors (RLRs), cytoplasmic DNA sensors (CDSs) and class A scavenger receptors (SR-As). It also discusses what is currently known of the downstream signaling molecules for each PRR family and the resulting antiviral response, either type I interferons (IFNs) or pro-inflammatory cytokine production. The review highlights what is known but also defines what still requires elucidation in this economically important animal. Understanding innate immune systems to virus infections will aid in the development of better antiviral therapies and vaccines for the future.