Blunt Snout Bream (Megalobrama amblycephala) MyD88 and TRAF6: Characterisation, Comparative Homology Modelling and Expression

MyD88 and TRAF6 play an essential role in the innate immune response in most animals. This study reports the full-length MaMyD88 and MaTRAF6 genes identified from the blunt snout bream (Megalobrama amblycephala) transcriptome profile. MaMyD88 is 2501 base pairs (bp) long, encoding a putative protein of 284 amino acids (aa), including the N-terminal DEATH domain of 78 aa and the C-terminal TIR domain of 138 aa. MaTRAF6 is 5474 bp long, encoding a putative protein of 542 aa, including the N-terminal low-complexity region, RING domain (40 aa), a coiled-coil region (64 aa) and C-terminal MATH domain (147 aa). Coding regions of MaMyD88 and MaTRAF6 genomic sequences consisted of five and six exons, respectively. Physicochemical and functional characteristics of the proteins were analysed. Alpha helices were dominant in the secondary structure of the proteins. Homology models of the MaMyD88 and MaTRAF6 domains were constructed applying the comparative modelling method. RT-qPCR was used to analyse the expression of MaMyD88 and MaTRAF6 mRNA transcripts in response to Aeromonas hydrophila challenge. Both genes were highly upregulated in the liver, spleen and kidney during the first 24 h after the challenge. While MyD88 and TRAF6 have been reported in various aquatic species, this is the first report and characterisation of these genes in blunt snout bream. This research also provides evidence of the important roles of these two genes in the blunt snout bream innate immune system.