METHODS: Using a mouse model, the miRNA expressions in liver during DENV-1 infection was investigated using high throughput miRNA sequencing. The differential expressions of miRNAs were then validated by qPCR, followed by target genes prediction. The identified miRNA targets were subjected to gene ontology (GO) annotation and pathway enrichment analysis to elucidate the potential biological pathways and molecular mechanisms associated with DENV-1 infection.
RESULTS: A total of 224 and 372 miRNAs out of 433 known mouse miRNAs were detected in the livers of DENV-1-infected and uninfected mice, respectively; of these, 207 miRNAs were present in both libraries. The miR-148a-3p and miR-122-5p were the two most abundant miRNAs in both groups. Thirty-one miRNAs were found to have at least 2-fold change in upregulation or downregulation, in which seven miRNAs were upregulated and 24 miRNAs were downregulated in the DENV-1-infected mouse livers. The miR-1a-3p was found to be the most downregulated miRNA in the DENV-1-infected mouse livers, with a significant fold change of 0.10. To validate the miRNA sequencing result, the expression pattern of 12 miRNAs, which were highly differentially expressed or most abundant, were assessed by qPCR and nine of them correlated positively with the one observed in deep sequencing. In silico functional analysis revealed that the adaptive immune responses involving TGF-beta, MAPK, PI3K-Akt, Rap1, Wnt and Ras signalling pathways were modulated collectively by 23 highly differentially expressed miRNAs during DENV-1 infection.
CONCLUSION: This study provides the first insight into the global miRNA expressions of mouse livers in response to DENV-1 infection in vivo and the possible roles of miRNAs in modulating the adaptive immune responses during DENV-1 infection.
Results: We generated 43 Gb of short Illumina reads and 9 Gb of long Nanopore reads, representing approximate genome coverage of 54× and 11×, respectively, based on the range of estimated k-mer-predicted genome sizes of between 791 and 967 Mbp. The final assembled genome is contained in 6404 scaffolds with an accumulated length of 880 Mb (96.3% BUSCO-calculated genome completeness). Compared with the Illumina-only assembly, the hybrid approach generated 94% fewer scaffolds with an 18-fold increase in N50 length (401 kb) and increased the genome completeness by an additional 16%. A total of 27 240 high-quality protein-coding genes were predicted from the clown anemonefish, 26 211 (96%) of which were annotated functionally with information from either sequence homology or protein signature searches.
Conclusions: We present the first genome of any anemonefish and demonstrate the value of low coverage (∼11×) long Nanopore read sequencing in improving both genome assembly contiguity and completeness. The near-complete assembly of the A. ocellaris genome will be an invaluable molecular resource for supporting a range of genetic, genomic, and phylogenetic studies specifically for clownfish and more generally for other related fish species of the family Pomacentridae.