Deep Illumina sequencing reveals conserved and novel microRNAs in grass carp in response to grass carp reovirus infection

Deep sequencing identified miR-193b-3p as a positive regulator of autophagy targeting Akt3 in Ctenopharyngodon idella CIK cells during GCRV infection

Recent research has highlighted complex and close interaction between miRNAs, autophagy, and viral infection. In this study, we observed the autophagy status in CIK cells infected with GCRV at various time points. We found that GCRV consistently induced cellar autophagy from 0 h to 12 h post infection. Subsequently, we performed deep sequencing on CIK cells infected with GCRV at 0 h and 12 h respectively, identifying 38 DEMs and predicting 9581 target genes. With the functional enrichment analyses of GO and KEGG, we identified 35 autophagy-related target genes of these DEMs, among which akt3 was pinpointed as the most central hub gene using module assay of the PPI network. Then employing the miRanda and Targetscan programs for prediction, and verification through a double fluorescent enzyme system and qPCR method, we confirmed that miR-193 b-3p could target the 3'-UTR of grass carp akt3, reducing its gene expression. Ultimately, we illustrated that grass carp miR-193 b-3p could promote autophagy in CIK cells. Above results collectively indicated that miRNAs might play a critical role in autophagy of grass carp during GCRV infection and contributed significantly to antiviral immunity by targeting autophagy-related genes. This study may provide new insights into the intricate mechanisms involved in virus, autophagy, and miRNAs.

miR-2188-5p promotes GCRV replication by the targeted degradation of klf2a in Ctenopharyngodon idellus

Studies on host immunity evasion by aquatic viruses have largely focused on coding genes. There is accumulating evidence for the important biological functions of non-coding miRNAs in virus-host interactions. The regulatory functions of non-coding miRNAs in fish reovirus-host interactions remain unknown. Here, miR-2188-5p in grass carp (Ctenopharyngodon idellus), a miRNA specific to teleosts, was predicted to target the 3' UTR of the transcription factor klf2a. A correlation analysis and dual-luciferase reporter assay revealed that miR-2188-5p could induce the degradation of klf2a. The expression of miR-2188-5p induced the degradation of klf2a in a dose-dependent manner, suppressing the type I interferon response and promoting grass carp reovirus (GCRV) replication. As determined by a co-expression analysis, klf2a inhibited viral infection when miR-2188-5p was overexpressed. The targeted degradation of klf2a by miR-2188-5p could inhibit the type I interferon response and promote the replication of GCRV; however, this targeted degradation ability was insufficient to fully inhibit GCRV infection. These results provide novel insights into the regulatory effects and biological functions of non-coding miRNAs in fish-virus interactions.

Expression characteristics and interaction networks of microRNAs in spleen tissues of grass carp (Ctenopharyngodon idella)

The spleen is an important immune organ in fish. MicroRNAs (miRNAs) have been shown to play an important role in the regulation of immune function. However, miRNA expression profiles and their interaction networks associated with the postnatal late development of spleen tissue are still poorly understood in fish. The grass carp (Ctenopharyngodon idella) is an important economic aquaculture species in China. Here, two small RNA libraries were constructed from the spleen tissue of healthy grass carp at one-year-old and three-year-old. A total of 324 known conserved miRNAs and 9 novel miRNAs were identified by using bioinformatic analysis. Family analysis showed that 23 families such as let-7, mir-1, mir-10, mir-124, mir-8, mir-7, mir-9, and mir-153 were highly conserved between vertebrates and invertebrates. In addition, 14 families such as mir-459, mir-430, mir-462, mir-7147, mir-2187, and mir-722 were present only in fish. Expression analysis showed that the expression patterns of miRNAs in the spleen of one-year-old and three-year-old grass carp were highly consistent, and the percentage of miRNAs with TPM > 100 was above 39%. Twenty significant differentially expressed (SDE) miRNAs were identified. Gene ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis showed that these SDE miRNAs were primarily involved in erythrocyte differentiation, lymphoid organ development, immune response, lipid metabolic process, the B cell receptor signaling pathway, the T cell receptor signaling pathway, and the PPAR signaling pathway. In addition, the following miRNA-mRNA interaction networks were constructed: immune and hematopoietic, cell proliferation and differentiation, and lipid metabolism. This study determined the miRNA transcriptome as well as miRNA-mRNA interaction networks in normal spleen tissue during the late development stages of grass carp. The results expand the number of known miRNAs in grass carp and are a valuable resource for better understanding the molecular biology of the spleen development in grass carp.

An SNP at the target site of cid-miR-nov-1043 in the TOLLIP 3' UTR decreases mortality rate in grass carp subjected to ENU-induced mutagenesis following grass carp reovirus infection

N-ethyl-N-nitrosourea (ENU) selection is a useful technique to generate new mutations that may cause some functional changes in the gene. Through our previous genomic bulked segregant analysis (BSA), one single nucleotide polymorphism (SNP) at the 3' UTR of Toll interacting protein gene (TOLLIP^982T>C) was identified in grass carp (Ctenopharyngodon idella) subjected to ENU-induced mutagenesis. We found that the overexpression of cid-miR-nov-1043 mimics significantly suppressed the luciferase activity of the TOLLIP 3' UTR, but TOLLIP^982T>C mutation at the target site can decrease the binding affinity between the miRNA cid-miR-nov-1043 and TOLLIP 3' UTR, reducing the inhibition of TOLLIP mRNA transcription in grass carp subjected to ENU-induced mutagenesis. More importantly, we demonstrated that TOLLIP mRNA transcription levels in the gills, liver, kidney and the isolate white cells of the mutant grass carp were significantly (p < 0.01) higher than those in the corresponding tissues from the wild-type grass carp following infection with Grass Carp Reovirus (GCRV) for seven days, while the downstream gene of TOLLIP transforming growth factor β-activated kinase 1 (TAK1) and TAK1-binding protein 1 (TAB1), were higher expressed in wild-type grass carp. As a negative regulator in the pro-inflammatory pathway of NF-κB, TOLLIP inhibits the excessive inflammation in ENU grass carp after GCRV infection. Consistent with the TOLLIP expression, histopathological results demonstrated more severe inflammation in wild-type grass carp, compared to the TOLLIP^982T>C mutant grass carp on the seventh day. Severe inflammation will lead to thoroughly infiltration of chloride and inflammatory cells in the gill filaments. This seriously hindered the exchange of oxygen, which ultimately disrupted blood circulation. Meanwhile, the survival rate of the mutant grass carp was significantly (p < 0.01) higher than that of the wild-type grass carp, indicating that the TOLLIP^982T>C mutants showed strong anti-viral abilities. Our results revealed that an SNP in the TOLLIP 3' UTR may contribute to the suppression of serve inflammation subjected to ENU-induced mutagenesis following GCRV infection, which may be helpful for future resistant breeding development of grass carp.

Integrated mRNA and microRNA Transcriptome Sequencing Characterizes Sequence Variants and mRNA-microRNA Regulatory Networks in Grass Carp Fibroblasts Infected with Virulent and Attenuated GCRV

Grass carp hemorrhagic disease is a fatal disease caused by the grass carp reovirus (GCRV). The aberrant regulation of transcripts has been implicated in many types of diseases. In the present study, we characterized mRNA and miRNA transcriptomes of different virulent GCRVs using RNA sequencing (RNA-Seq). One hundred eighteen miRNAs were identified as being differentially expressed between different virulent viruses in grass carp fibroblasts. Eight miRNAs were selected to verify the RNA-Seq results using RT-PCR and mRNA methods. In total, 996 differentially expressed mRNA genes were identified in grass carp fibroblasts, while 901 miRNA-mRNA target pairs were observed to be inversely regulated in grass carp fibroblasts. Integrated mRNA/miRNA expression profiling analysis results showed that the most influenced processes were the immune response and cell death. Three miRNAs were shown to exhibit the same expression patterns when two different methods were used and had important functions during viral infection. These results provide insights into the miRNA-mediated regulation of mRNA and valuable resources on transcript variation and regulation during GCRV infection, which are potentially useful for mechanistic and drug studies.

Global profiling of lncRNAs-miRNAs-mRNAs reveals differential expression of coding genes and non-coding RNAs in the lung of beagle dogs at different stages of Toxocara canis infection

The roundworm Toxocara canis causes toxocariasis in dogs and larval migrans in humans. Better understanding of the lung response to T. canis infection could explain why T. canis must migrate to and undergoes part of its development inside the lung of the definitive host. In this study, we profiled the expression patterns of long non-coding RNAs (lncRNAs), microRNAs (miRNAs), and mRNAs in the lungs of Beagle dogs infected by T. canis, using high throughput RNA sequencing. At 24 h p.i., 1,012 lncRNAs, 393 mRNAs and 10 miRNAs were differentially expressed (DE). We also identified 883 DElncRNAs, 264 DEmRNAs and 20 DEmiRNAs at 96 h p.i., and 996 DElncRNAs, 342 DEmRNAs and eight DEmiRNAs at 36 days p.i., between infected and control dogs. Significant changes in the levels of expression of transcripts related to immune response and inflammation were associated with the antiparasitic response of the lung to T. canis. The remarkable increase in the expression of scgb1a1 at all time points after infection suggests the need for consistent moderation of the excessive inflammatory response. Also, upregulation of foxj1 at 24 h p.i., and downregulation of IL-1β and IL-21 at 96 h p.i., suggest an attenuation of the humoral immunity of infected dogs. These results indicate that T. canis pathogenesis in the lung is mediated through contributions from both pro-inflammatory and anti-inflammatory mechanisms. Competing endogenous RNA (ceRNA) network analysis revealed significant interactions between DElncRNAs, DEmiRNAs and DEmRNAs, and improved our understanding of the ceRNA regulatory mechanisms in the context of T. canis infection. These data provide comprehensive understanding of the regulatory networks that govern the lung response to T. canis infection and reveal new mechanistic insights into the interaction between the host and parasite during the course of T. canis infection in the canine.

Dual RNA-Seq uncovers the function of an ABC transporter gene in the host-pathogen interaction between Epinephelus coioides and Pseudomonas plecoglossicida

As an important pathogen in aquaculture, Pseudomonas plecoglossicida has caused heavy losses. The expression of an ABC transporter gene-L321_23611 of P. plecoglossicida at 18 °C was found significant higher than those at 28 °C by RNA-seq and qRT-PCR. RNAi significantly reduced the content of L321_23611 mRNA in P. plecoglossicida with a maximal decrease of 89.2%. Compared with the wild type strain, the infection of L321_23611-RNAi strain resulted in the reduction in mortality and the onset time delay of a kind of marine teleosts, Epinephelus coioides. The results of dual RNA-seq showed that the RNAi of L321_23611 resulted in a significant change in both pathogen and host transcriptome in the spleens of infected E. coioides. The result of GO and KEGG analysis from dual RNA-seq data showed both host genes of chemokine signaling pathway, coagulation and complement system, hematopoietic cell lineage pathway as well as hemoglobin complex GO term and pathogenic genes of bacterial-type flagellum-dependent cell mortality GO term and flagellar assembly, biosynthesis of amino acids and lysine biosynthesis systems pathways were mainly affected by L321_23611 gene of P. plecoglossicida. The results indicated that: 1. ABC transporter gene-L321_23611 was a virulent gene of P. plecoglossicida. 2. Both the activation of the host immune pathways and depression of pathogenic virulence-related pathways facilitated E. coioides to remove L321_23611-RNAi strain than the wild type strain of P. plecoglossicida.

Investigating the Role of BATF3 in Grass Carp (Ctenopharyngodon idella) Immune Modulation: A Fundamental Functional Analysis

Basic leucine zipper transcription factor ATF-like (BATF)-3, belonging to activator protein 1 (AP-1) superfamily transcription factors, is essential for homeostatic development of CD8α⁺ classical dendritic cells activating CD8 T-cell responses to intracellular pathogens. In this study, the characteristics and cDNA cloning of the CiBATF3 molecule were described in grass carp (Ctenopharyngodon idella). CiBATF3 had abundant expression in immune-related organizations, including liver, spleen and gill, and grass carp reovirus (GCRV) infection had significantly changed its expression level. After Ctenopharyngodon idella kidney (CIK) cells were challenged with pathogen-associated molecular patterns (PAMPs), polyinosinic:polycytidylic acid (poly(I:C)) stimulation induced higher mRNA levels of CiBATF3 than that of lipopolysaccharide (LPS). Subcellular localization showed that CiBATF3-GFP was entirely distributed throughout cells and nuclear translocation of CiBATF3 was found after poly(I:C) treatment. Additionally, the interaction between CiBATF3 and interleukin 10 (IL-10) was proven by bimolecular fluorescence complementation (BiFC) system. The small interfering RNA (siRNA)-mediated CiBATF3 silencing showed that the mRNA of CiBATF3 and its downstream genes were down-regulated in vitro and in vivo. CiBATF3 played a negative regulatory role in the transcriptional activities of AP-1 and NF-κB reporter gene. In summary, the results may provide valuable information on fundamental functional mechanisms of CiBATF3.

Deep sequencing of small RNAs from 11 tissues of grass carp Ctenopharyngodon idella and discovery of sex-related microRNAs

This research identified 169 known microRNAs (miRNAs), 380 novel miRNAs, and 30,538 targets in 11 tissues (blood, brain, derma, gill, heart, intestine, kidney, liver, muscle, pronephros, and spleen) from grass carp Ctenopharyngodon idella with high-throughput sequencing (HTS). Transcripts per million (TPM) expression analysis detected 41 brain-enriched miRNAs (accounting for 61.19% of all tissue-enriched miRNAs). Real-time quantitative PCR (RT-qPCR) confirmed that 21 of 24 randomly selected tissue-enriched miRNAs from the TPM analysis were indeed tissue-enriched (P < 0.05), suggesting the HTS and TPM analyses were reliable. Nine of the 41 brain-enriched miRNAs are complementary to members of the double-sex and mab-3 related transcription factor family (dmrt) involved in sex differentiation. RT-qPCR revealed that cid-miR-138 was more highly expressed in testis than in ovary (P < 0.01), while the reverse was true for target gene dmrt4a (P < 0.01). This opposite expression pattern suggested the direct participation of cid-miR-138-dmrt4a in neuroendocrine mechanisms related to brain-pituitary networks during sex development. The discovery of miRNAs from 11 C. idella tissues expands the available fish miRNA database, and enhances our understanding of the role of sex-related miRNAs in tissue differentiation and maintenance of specific tissue functions in fishes.

Selenium deficiency impaired immune function of the immune organs in young grass carp (Ctenopharyngodon idella)

This study aimed to investigate the effects of dietary selenium on resistance to skin haemorrhages and lesions and on immune function as well as the underlying mechanisms of those effects in the head kidney, spleen and skin of young grass carp (Ctenopharyngodon idella). A total of 540 healthy grass carp with initial body weight (226.48 ± 0.68 g) were randomly divided into six groups and fed six separate diets with graded dietary levels of selenium (0.025, 0.216, 0.387, 0.579, 0.795 and 1.049 mg/kg diet) for 80 days. After the feeding period, an immunization trial was performed by infection with Aeromonas hydrophila for 14 days. The results showed that, compared with the optimal selenium level, (1) selenium deficiency impaired the production of antibacterial compounds and immunoglobulins and down-regulated the transcript abundances of antimicrobial peptides and selenoproteins; (2) selenium deficiency aggravated inflammatory responses in part by up-regulating pro-inflammatory cytokines and down-regulating anti-inflammatory cytokines mRNA levels, which were partially related to [IKKα, β, γ/IκBα/NF-κB] signalling and [TOR/(S6K1, 4E-BP1)] signalling, respectively. Interestingly, selenium deficiency had no effect on the expression of TGF-β2, IL-4/13B, IL-10, IL-12p35, IL-15 (skin only) or 4E-BP2 in the head kidney, spleen and skin of young grass carp. Finally, based on the percent weight gain (PWG), the morbidity of skin haemorrhages and lesions, the ACP activity in the head kidney and the lysozyme activity in spleen, the optimal dietary selenium requirements for young grass carp were estimated to be 0.546-0.604 mg/kg diet. In summary, selenium deficiency decreased the growth performance and impaired the immune function in the head kidney, spleen and skin of young grass carp.

Global and Complement Gene-Specific DNA Methylation in Grass Carp after Grass Carp Reovirus (GCRV) Infection

Grass carp reovirus (GCRV) causes huge economic loss to the grass carp cultivation industry but the mechanism remains largely unknown. In this study, we investigated the global and complement gene-specific DNA methylation in grass carp after GCRV infection aimed to uncover the mechanism underlying GCRV infection. The global DNA methylation level was increased after GCRV infection. Expression levels of enzymes involved in DNA methylation including DNA methyltransferase (DNMT), ten-eleven translocation proteins (TETs), and glycine N-methyltransferase (GNMT) were significantly altered after GCRV infection. In order to investigate the relationship between the gene expression level and DNA methylation level, two representative complement genes, complement component 3 (C3) and kininogen-1 (KNG1), were selected for further analysis. mRNA expression levels of the two genes were significantly increased at 5 and 7 days after GCRV infection, whereas the DNA methylation level at the 5′ flanking regions of the two genes were down-regulated at the same time-points. Moreover, a negative correlation was detected between gene expression levels and DNA methylation levels of the two genes. Therefore, the current data revealed a global and complement gene-specific DNA methylation profile after GCRV infection. Our study would provide new insights into understanding the mechanism underlying GCRV infection.

Deep Circular RNA Sequencing Provides Insights into the Mechanism Underlying Grass Carp Reovirus Infection

Grass carp hemorrhagic disease, caused by the grass carp reovirus (GCRV), is a major disease that hampers the development of grass carp aquaculture in China. The mechanism underlying GCRV infection is still largely unknown. Circular RNAs (circRNAs) are important regulators involved in various biological processes. In the present study, grass carp were infected with GCRV, and spleen samples were collected at 0 (control), 1, 3, 5, and 7 days post-infection (dpi). Samples were used to construct and sequence circRNA libraries, and a total of 5052 circRNAs were identified before and after GCRV infection, of which 41 exhibited differential expression compared with controls. Many parental genes of the differentially expressed circRNAs are involved in metal ion binding, protein ubiquitination, enzyme activity, and nucleotide binding. Moreover, 72 binding miRNAs were predicted from the differentially expressed circRNAs, of which eight targeted genes were predicted to be involved in immune responses, blood coagulation, hemostasis, and complement and coagulation cascades. Upregulation of these genes may lead to endothelial and blood cell damage and hemorrhagic symptoms. Our results indicate that an mRNA–miRNA–circRNA network may be present in grass carp infected with GCRV, providing new insight into the mechanism underlying grass carp reovirus infection.