Differential interferon system gene expression profiles in susceptible and resistant gynogenetic clones of gibel carp challenged with herpesvirus CaHV

Interferon (IFN) system plays a vital role in the first line of defense against viruses. In this study, we first identified multiple transcripts of 15 IFN system genes, including PRRs (TLR2, TLR3, RIG-I, and LGP2), PRR-mediated IFN signal pathway (MyD88, MITA, and MAVS), IFN regulatory factors (IRF1, IRF3, IRF7, and IRF9), IFNs (IFNφ1 and IFNφ3), and ISGs (Mx and viperin), and one transcript of TLR9 in de novo transcriptome assembly data of gibel carp head-kidney. Multiple nucleotide alignments and phylogenetic analysis of common region showed that the transcripts of every of the 15 IFN system genes were classified into two homologs with distinctly divergent sequences, indicating that hexaploid gibel carp may be an allopolyploid. During Carassius auratus herpesvirus (CaHV) infection, gibel carp resistant clone H significantly suppressed CaHV replication with markedly less viral loads than those in highly susceptible clone A+ and moderately resistant clone F. Then, qPCR analyses were performed to reveal their differential and dynamic expression changes during CaHV infection in head kidney, spleen and liver among three gibel carp gynogenetic clones. Through qPCR and hierarchical clustering analysis, 8 genes, such as RIG-Is, LGP2s, IRF1-B, IRF3s, IRF7s, IRF9-B, Mxs, and viperins, were identified as candidate resistant-related genes. They remarkably increased their expression in immune tissues of three clones after CaHV infection. Significantly, the up-regulation folds of these genes in clone A+, F and H were related to their resistance ability to CaHV, progressively increasing from susceptible clone to resistant clone at 1 dpi. The positive correlation to the resistance ability suggested that resistant clone H immediately triggered stronger IFN response. IFNφ3 showed a different dynamic change and was sharply induced in moderately resistant clone F at 3 dpi. The other 5 IFN system genes (TLR2, TLR3, TLR9, MyD88, and MITA) maintained a low expression level after CaHV challenge. Interestingly, the A or B copies/homologs of almost these IFN system genes exhibited differential transcript abundance in immune tissue after CaHV challenge, suggesting A or B homologs might occur dominant or biased expression of homeologs during gibel carp evolution. These data provide candidate resistant-related genes for disease-resistance breeding of gibel carp.