Seasonal fungal biome of mosquitoes determined via ITS region sequencing
Comparison of Shannon diversity between samples groups
Changes between season/location (A), season (B), and field vs. lab (C) are shown. Significance determined via Kruskal-Wallis chi-square with Dunn’s multiple comparison adjustment where applicable.
Comparison of mean number of sequence variants (A) and genera (B) for each season/location
Significance determined via Kruskal-Wallis chi-square with Dunn’s multiple comparison adjustment where applicable.
High presence of sample singletons at the Genus (A) and sequence variant (B) level
Genus-level DESeq2 analysis of samples between seasons and locations
A negative log2 fold change represents a genus-agglomerated taxon that loses read prevalence from the Wet to Dry Season, and a positive log2 fold change represents a genus that increases during this transition. p-adj is the Benjamini-Hochberg adjusted p-value accounting for a false discovery rate of 0.1. Each panel compares transition from Wet (left) to Dry (right, panels A and B) or from Field (left) to Laboratory (right, panel C).
Random forest supervised learning to discriminate season and location for each sample group
Analysis is performed on all samples (pooled and individuals) for each location. The top ten variables (genera) important to the created model is shown in Panel B. Significantly different center points in the ordination via PERMANOVA are present in Panel C. Prevalence of the genera for each sample group is presented in Panel D.
Sequences variants present in 19 most abundant genera
Phylogenetic tree representing each 16S sequence variant (A) and OTU-level (B, 97% identity) for the 19 most prevalent taxa found in the Anopheles coluzzii samples.
Supplementary information on the 19 most abundant genera in the study
Presence in arthropods is from Minard et al. unless otherwise noted. Non-insect environment and presence in arthropods is not comprehensive.