Mass rearing the black soldier fly, Hermetia illucens, for waste bioremediation and valorisation is gaining traction on a global scale. While the health and productivity of this species are underpinned by associations with microbial taxa, little is known about the factors that govern gut microbiome assembly, function, and contributions towards host phenotypic development in actively feeding larvae. In the present study, a 16S rDNA gene sequencing approach applied to a study system incorporating both feed substrate and genetic variation is used to address this knowledge gap. It is determined that the alpha diversity of larval gut bacterial communities is driven primarily by features of the larval feed substrate, including the diversity of exogenous bacterial populations. Microbiome beta diversity, however, demonstrated patterns of differentiation consistent with an influence of diet, larval genetic background, and a potential interaction between these factors. Moreover, evidence for an association between microbiome structure and the rate of larval fat accumulation was uncovered. Taxonomic enrichment analysis and clustering of putative functional gut profiles further suggested that feed-dependent turnover in microbiome communities is most likely to impact larval characteristics. Taken together, these findings indicate that host–microbiome interactions in this species are complex yet relevant to larval trait emergence.

No abstract

The micro-evolutionary forces that shape genetic diversity during domestication have been assessed in many plant and animal systems. However, the impact of these processes on gene expression, and consequent functional adaptation to artificial environments, remains under-investigated. In this study, whole-transcriptome dynamics associated with the early stages of domestication of the black soldier fly (BSF), Hermetia illucens, were assessed. Differential gene expression (DGE) was evaluated in relation to (i) generational time within the cultured environment (F2 vs. F3), and (ii) two selection strategies [no artificial selective pressure (NS); and selection for greater larval mass (SEL)]. RNA-seq was conducted on 5th instar BSF larvae (n = 36), representing equal proportions of the NS (F2 = 9; F3 = 9) and SEL (F2 = 9; F3 = 9) groups. A multidimensional scaling plot revealed greater gene expression variability within the NS and F2 subgroups, while the SEL group clustered separately with lower levels of variation.Comparisons between generations revealed 898 differentially expressed genes (DEGs; FDR-corrected p < 0.05), while between selection strategies, 213 DEGs were observed (FDR-corrected p < 0.05). Enrichment analyses revealed that metabolic, developmental, and defence response processes were over-expressed in the comparison between F2 and F3 larvae, while metabolic processes were the main differentiating factor between NS and SEL lines. This illustrates the functional adaptations that occur in BSF colonies across generations due to mass rearing; as well as highlighting genic dynamics associated with artificial selection for production traits that might inform future selective breeding strategies.