Genomic sequence comparison, promoter activity, SNP detection of RIG-I gene and association with resistance/susceptibility to grass carp reovirus in grass carp (Ctenopharyngodon idella)

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.

Teleosts Genomics: Progress and Prospects in Disease Prevention and Control

Genome wide studies based on conventional molecular tools and upcoming omics technologies are beginning to gain functional applications in the control and prevention of diseases in teleosts fish. Herein, we provide insights into current progress and prospects in the use genomics studies for the control and prevention of fish diseases. Metagenomics has emerged to be an important tool used to identify emerging infectious diseases for the timely design of rational disease control strategies, determining microbial compositions in different aquatic environments used for fish farming and the use of host microbiota to monitor the health status of fish. Expounding the use of antimicrobial peptides (AMPs) as therapeutic agents against different pathogens as well as elucidating their role in tissue regeneration is another vital aspect of genomics studies that had taken precedent in recent years. In vaccine development, prospects made include the identification of highly immunogenic proteins for use in recombinant vaccine designs as well as identifying gene signatures that correlate with protective immunity for use as benchmarks in optimizing vaccine efficacy. Progress in quantitative trait loci (QTL) mapping is beginning to yield considerable success in identifying resistant traits against some of the highly infectious diseases that have previously ravaged the aquaculture industry. Altogether, the synopsis put forth shows that genomics studies are beginning to yield positive contribution in the prevention and control of fish diseases in aquaculture.

Grass Carp Laboratory of Genetics and Physiology 2 Serves As a Negative Regulator in Retinoic Acid-Inducible Gene I- and Melanoma Differentiation-Associated Gene 5-Mediated Antiviral Signaling in Resting State and Early Stage of Grass Carp Reovirus Infection

Laboratory of genetics and physiology 2 (LGP2) is a key component of RIG-I-like receptors (RLRs). However, the lack of the caspase recruitment domains (CARDs) results in its controversial functional performance as a negative or positive regulator in antiviral responses. Especially, no sufficient evidence uncovers the functional mechanisms of LGP2 in RLR signaling pathways in teleost. Here, negative regulation mechanism of LGP2 in certain situations in retinoic acid-inducible gene I (RIG-I) and melanoma differentiation-associated gene 5 (MDA5)-mediated antiviral responses was identified in Ctenopharyngodon idella kidney cells. LGP2 overexpression inhibits synthesis and phosphorylation of interferon regulatory factor 3/7 (IRF3/7), and mRNA levels and promoter activities of IFNs and NF-κBs in resting state and early phase of grass carp reovirus (GCRV) infection. Knockdown of LGP2 obtains opposite effects. Luciferase report assay indicates that LGP2 works at the upstream of RIG-I and MDA5. LGP2 binds to RIG-I and MDA5 with diverse domain preference and which is independent of GCRV infection. Furthermore, LGP2 restrains K63-linked ubiquitination of RIG-I and MDA5 in various degrees. These differences result in disparate repressive mechanisms of LGP2 to RIG-I- and MDA5-mediated signal activations of IFN-β promoter stimulator 1 and mediator of IRF3 activation. Interestingly, LGP2 also inhibits K48-linked RIG-I and MDA5 ubiquitination to suppress proteins degradation, which guarantees the basal protein levels for subsequently rapid signal activation. All these results reveal a mechanism that LGP2 functions as a suppressor in RLR signaling pathways to maintain cellular homeostasis in resting state and early phase during GCRV infection.

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 of the MmeHairy gene and expression analysis affected by two SNPs in the 3'-untranslated region in the clam Meretrix meretrix

As a bHLH transcriptional repressor, Hairy-related proteins can bind to DNA sites in target gene promoters and negatively regulate gene transcription. In the present study, the full-length cDNA of Hairy was obtained from the clam Meretrix meretrix (MmeHairy), and two SNPs in the 3'-untranslated region (UTR) of this gene, SNP1066 and 1067, were identified and characterized. Multiple sequence alignment and phylogenetic analysis revealed that MmeHairy belongs to the Hairy protein subfamily. Analysis of tissue expression patterns showed that the mRNA of MmeHairy had the highest expression level in the hepatopancreas. The expression levels of MmeHairy were up-regulated in the hepatopancreas after Vibrio challenge. Genotyping and quantitative analysis showed that the mRNA levels of MmeHairy were significantly different among individual clams with different genotypes at SNP1066 and 1067 (P < 0.05), which indicated that these two SNP loci may affect the expression of MmeHairy and could be used as candidate markers for future selection in M. meretrix breeding programs.

Development of Vibrio spp. infection resistance related SNP markers using multiplex SNaPshot genotyping method in the clam Meretrix meretrix

The clam Meretrix meretrix is a commercially important mollusc species in the coastal areas of South and Southeast Asia. In the present study, large-scale SNPs were genotyped by the Multiplex SNaPshot genotyping method among the stocks of M. meretrix with different Vibrio spp. infection resistance profile. Firstly, the AUTOSNP software was applied to mine SNPs from M. meretrix transcriptome, and 323 SNP loci (including 120 indels) located on 64 contigs were selected based on Uniprot-GO associations. Then, 38 polymorphic SNP loci located on 15 contigs were genotyped successfully in the clam stocks with different resistance to Vibrio parahaemolyticus infection (11-R and 11-S groups). Pearson's Chi-square test was applied to compare the allele and genotype frequency distributions of the SNPs between the different stocks, and seven SNP markers located on three contigs were found to be associated with V. parahaemolyticus infection resistance trait. Haplotype-association analysis showed that six haplotypes had significantly different frequency distributions in 11-S and 11-R (P < 0.05). With selective genotyping between 09-R and 09-C populations, which had different resistance to Vibrio harveyi infection, four out of the seven selected SNPs had significantly different distributions (P < 0.05) and therefore they were considered to be associated with Vibrio spp. infection resistance. Sequence alignments and annotations indicated that the contigs containing the associated SNPs had high similarity to the immune related genes. All these results would be useful for the future marker-assisted selection of M. meretrix strains with high Vibrio spp. infection resistance.

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

Toll-like receptor 8 (TLR8), a prototypical intracellular member of TLR family, is generally linked closely to antiviral innate immune through recognizing viral nucleic acid. In this study, 5'-flanking region of Ctenopharyngodon idella TLR8 (CiTLR8), 671bp in length, was amplified and eight SNPs containing one SNP in the intron, three SNPs in the coding region (CDS) and four SNPs in the 3'-untranslated region (UTR) were identified and characterized. Of which 4062 A/T was significantly associated with the susceptibility/resistance to GCRV both in genotype and allele (P < 0.05), while 4168 C/T was extremely significantly associated with that (P < 0.01) according to the case (susceptibility)-control (resistance) analysis. Following the verification experiment, further analyses of mRNA expression, linkage disequilibrium (LD), haplotype and microRNA (miRNA) target site indicated that 4062 A/T and 4168 C/T in 3'-UTR might affect the miRNA regulation, while the exertion of antiviral effects of 4062 A/T might rely on its interaction with other SNPs. Additionally, the high-density of SNPs in 3'-UTR might reflect the specific biological functions of 3'-UTR. And also, the mutation of 747 A/G in intron changing the potential transcriptional factor-binding sites (TFBS) nearby might affect the expression of CiTLR8 transcriptionally or post-transcriptionally. Moreover, as predicted, the A/G transition of the only non-synonymous SNP (3846 A/G) in CDS causing threonine/alanine variation, could shorten the length of the α-helix and ultimately affect the integrity of the Toll-IL-1 receptor (TIR) domain. The functional mechanism of 3846 A/G might also involve a threonine phosphorylation signaling. This study may broaden the knowledge of TLR polymorphisms, lay the foundation for further functional research of CiTLR8 and provide potential markers as well as theoretical basis for resistance molecular breeding of grass carp against GCRV.

Insights into the antiviral immunity against grass carp (Ctenopharyngodon idella) reovirus (GCRV) in grass carp

Global fish production from aquaculture has rapidly grown over the past decades, and grass carp shares the largest portion. However, hemorrhagic disease caused by grass carp reovirus (GCRV) results in tremendous loss of grass carp (Ctenopharyngodon idella) industry. During the past years, development of molecular biology and cellular biology technologies has promoted significant advances in the understanding of the pathogen and the immune system. Immunoprophylaxis based on stimulation of the immune system of fish has also got some achievements. In this review, authors summarize the recent progresses in basic researches on GCRV; viral nucleic acid sensors, high-mobility group box proteins (HMGBs); pattern recognition receptors (PRRs), Toll-like receptors (TLRs) and retinoic acid inducible gene I- (RIG-I-) like receptors (RLRs); antiviral immune responses induced by PRRs-mediated signaling cascades of type I interferon (IFN-I) and IFN-stimulated genes (ISGs) activation. The present review also notices the potential applications of molecule genetic markers. Additionally, authors discuss the current preventive and therapeutic strategies (vaccines, RNAi, and prevention medicine) and highlight the importance of innate immunity in long term control for grass carp hemorrhagic disease.

The first invertebrate RIG-I-like receptor (RLR) homolog gene in the pacific oyster Crassostrea gigas

Retinoic acid-inducible gene I (RIG-I)-like receptor (RLR) is a pivotal receptor that detects numerous RNA and DNA viruses and mediates the innate induction of interferons and pro-inflammatory cytokines upon viral infection. In the present study, we cloned and characterized the first RIG-I gene in a marine mollusk, Crassostrea gigas, and designated it as CgRIG-I. The full-length CgRIG-I cDNA is 3436 bp, including 5'- and 3'-untranslated regions (UTRs) of 93 bp and 286 bp, respectively, and an open reading frame (ORF) of 3057 bp. The gene encodes a 1018 amino acid polypeptide with an estimated molecular mass of 116.5 kDa. SMART analysis showed that the CgRIG-I protein had the typical conserved domains, including the caspase activation and recruitment domains (CARDs), the RNA helicase domain and the C-terminal regulatory domain (RD). Phylogenetic analysis revealed that CgRIG-I was grouped into the clade of its vertebrate homologs. Moreover, CgRIG-I expression could be specifically increased after stimulation by poly(I:C) rather than by other PAMPs such as lipopolysaccharide (LPS), peptidoglycan (PGN), heat-killed Listeria monocytogenes (HKLM) and heat-killed Vibrio alginolyticus (HKVA). Meanwhile, six IRF, three STAT and one NF-κB predicted sites were identified in the CgRIG-I promoter, which was consistent with its high responsiveness to poly(I:C). In summary, this report provides the first CgRIG-I sequence of a mollusk, but its function in the antiviral immune response requires further investigation.

Identification of a novel N4BP1-like gene from grass carp (Ctenopharyngodon idella) in response to GCRV infection

Nedd4 binding protein 1 (N4BP1) has been identified as an interacting protein and a substrate of Nedd4 E3 ligase. However, the report about N4BP1's function is limit. In this study, a novel N4BP1 gene (CiN4BP1) was cloned from grass carp (Ctenopharyngodon idella). The full-length cDNA sequence of CiN4BP1 (3022 bp) included an open reading frame (ORF) of 2565 bp, which encoded a putative peptides of 854 amino acids containing one KH domain and one NYN domain. It was close homology (47% identify) to Oryzias latipes N4BP1. And mRNA expression of CiN4BP1 gene showed relatively high level in skin, gill, head kidney and spleen. After grass carp reovirus (GCRV) infection, CiN4BP1 was up-regulated in vivo and in vitro. Furthermore, overexpression of CiN4BP1 in CIK cells inhibited viral gene transcription. These data indicated that CiN4BP1 might play an important role in immune response to viral invasion.

Gene-based polymorphisms, genomic organization of interferon-β promoter stimulator 1 (IPS-1) gene and association study with the natural resistance to grass carp reovirus in grass carp Ctenopharyngodon idella

Interferon-β promoter stimulator 1 (IPS-1) plays a pivotal role in the production of type I interferon (IFN) and pro-inflammatory cytokines. Though its function in innate immunity was widely studied, its genetic polymorphisms and evolution information were rarely known. Herein, the present study firstly identified and characterized the genomic organization of CiIPS-1 (IPS-1 gene of grass carp, Ctenopharyngodon idella) and its nucleotide polymorphisms. CiIPS-1 consists of 8789 bp, with a 314 bp 5'-flanking region, 6 exons and 5 introns. After a challenge experiment, statistical analysis was performed to assess the association between the 17 polymorphisms of CiIPS-1 and the natural resistance of C. idella to grass carp reovirus (GCRV), which revealed that -3741 C/T in the intron 1, 933 C/G in the coding sequence and 2299 G/T in the last intron were significant (P < 0.05). The subsequent secondary challenge test further confirmed that -3741TT group were more resistant than -3741CC group to GCRV infection, while 933GG and 2299GG stocks were respectively more susceptible than 933CC and 2299TT stocks (P < 0.05). In addition, haplotype and polymorphism-polymorphism interaction analyses uncovered four significant haplotypes (P < 0.05) and two notable polymorphisms. These findings may provide basic data for the further functional research of CiIPS-1, and genetic markers for molecular and transgenic breeding of C. idella.