Comparative transcriptional analysis between virulent isolate HN1307 and avirulent isolate GD1108 of grass carp reovirus genotype II

As a widespread epidemic virus, genotype II of the grass carp reovirus poses a significant threat to the grass carp farming industry in China. Different genotype II isolates cause different degrees of virulence, although the underlying pathogenic mechanisms remain largely unknown. In this work, infections of grass carp with the virulent isolate grass carp reovirus (GCRV)-HN1307 and the avirulent isolate GCRV-GD1108 were performed to reveal a possible mutual transcriptional discrepancy. More differentially expressed genes (DEGs) were identified in the HN1307-infected group, which defined a grossly similar gene ontology (GO) pattern and different pathway landscape as the GD1108-infected group. Gene set enrichment analysis revealed that pathways related to innate immunity and metabolism were reciprocally activated and suppressed, respectively, following infection withHN1307, compared with GD1108. The trend analysis further indicated that immune-related pathways were involved in one of the four statistically significant profiles. Network analysis of transcription factor-gene interactions and protein-protein interactions on the immune-related profile suggested that among the core transcriptional factors (TFs) (UBTF, HCFC1, MAZ, MAX, and NRF1) and the hub proteins (Tlr3, Tlr7, Tlr9, Irf3, and Irf7), the latter were highly enriched in the toll-like receptor signaling pathway. Real-time quantitative PCR performed on the selected mRNAs validated the relative expression. This work will provide insights into the distinct transcriptional signatures from avirulent and virulent isolates of GCRV, which may contribute to the development of products for prevention.