Microbiota identified from preserved Anopheles

Mosquito species from the Anopheles gambiae complex and the Anopheles funestus group are dominant African malaria vectors. Mosquito microbiota play vital roles in physiology and vector competence. Recent research has focused on investigating the mosquito microbiota, especially in wild populations. Wild mosquitoes are preserved and transported to a laboratory for analyses. Thus far, microbial characterization post-preservation has been investigated in only Aedes vexans and Culex pipiens. Investigating the efficacy of cost-effective preservatives has also been limited to AllProtect reagent, ethanol and nucleic acid preservation buffer. This study characterized the microbiota of African Anopheles vectors: Anopheles arabiensis (member of the An. gambiae complex) and An. funestus (member of the An. funestus group), preserved on silica desiccant and RNAlater® solution. Microbial composition and diversity were characterized using culture-dependent (midgut dissections, culturomics, MALDI-TOF MS) and culture-independent techniques (abdominal dissections, DNA extraction, next-generation sequencing) from laboratory (colonized) and field-collected mosquitoes. Colonized mosquitoes were either fresh (non-preserved) or preserved for 4 and 12 weeks on silica or in RNAlater®. Microbiota were also characterized from field-collected An. arabiensis preserved on silica for 8, 12 and 16 weeks. Elizabethkingia anophelis and Serratia oryzae were common between both vector species, while Enterobacter cloacae and Staphylococcus epidermidis were specific to females and males, respectively. Microbial diversity was not influenced by sex, condition (fresh or preserved), preservative, or preservation time-period; however, the type of bacterial identification technique affected all microbial diversity indices. This study broadly characterized the microbiota of An. arabiensis and An. funestus. Silica- and RNAlater®-preservation were appropriate when paired with culture-dependent and culture-independent techniques, respectively. These results broaden the selection of cost-effective methods available for handling vector samples for downstream microbial analyses.

Culturomics, a culture-dependent technique that involves growing bacteria using nutrient media, is commonly used to characterize the mosquito's midgut bacteria as it is fast, cost-effective and provides reliable data [68][69][70][71]. Morphologically distinct colonies are isolated and subjected to matrix-assisted laser desorption ionization-time of flight (MALDI-TOF) mass spectrometry (MS). MALDI-TOF MS identifies bacteria based on their proteome, where proteins are cleaved into peptides and their molecular masses are used to create peptide mass fingerprints (PMFs) (reviewed by [71]). The PMFs of unknown bacteria are compared with the PMFs of known bacteria in a database for taxonomic identification [71].
However, certain bacteria cannot grow on selective media, and species identification is limited to the local database installed on the MALDI-TOF MS system [72,73]. Next-generation sequencing (NGS), a sensitive, culture-independent approach, addresses these downfalls (reviewed by [74]). NGS identifies bacteria based on their genome: conserved regions of the prokaryotic 16S ribosomal ribonucleic acid (rRNA) gene are used for amplification, and hypervariable regions of the gene are used to identify taxa [75][76][77]. Yet, NGS is vulnerable to bias and can be costly and time-consuming [74,[78][79][80][81].
As mosquito microbial studies are increasingly shifting to field-collected samples, field-caught mosquitoes are preserved and transported to a laboratory for analysis. This is because testing mosquitoes in field conditions is impractical due to the lack of a sterile environment and laboratory equipment. Field sites are also often far from suitably equipped laboratories. Although, the type of preservation method used is dependent on the type of downstream analysis being performed as certain preservatives are better suited for identifying specific entomological indicators [89].
The microbiota of Aedes vexans and Culex pipiens have been identified post-preservation from AllProtect reagent, ethanol, and NAP buffer [93]. However, the efficacy of other commonly used, cost-effective preservatives, such as silica and RNAlater ® , has not been investigated. Silica preserves large quantities of specimens and ensures long-term preservation at room temperature [89,92,95], while RNAlater ® preserves high-quality DNA and RNA and is most suitable for determining internal muscular anatomy [89]. Furthermore, the microbiota of preserved African Anopheles vectors has not been investigated.
Accordingly, this study assessed if the microbiota of African Anopheles vectors could be identified postpreservation from silica and RNAlater ® using culturedependent and culture-independent techniques. The microbiota of laboratory (colonized) Anopheles arabiensis (member of the An. gambiae complex) and An. funestus (member of the An. funestus group) were screened after preserving mosquitoes for up to 12 weeks with each preservative. Additionally, the microbiota of preserved field-collected An. arabiensis were characterized.

Biological material
Colonized mosquitoes were obtained from the Botha de Meillon Insectary, National Institute for Communicable Diseases (NICD), Johannesburg, South Africa. Two Anopheles species were used in this study: An. arabiensis (MBN colony) and An. funestus (FUMOZ colony). The MBN colony has mosquitoes from Mamfene, KwaZulu-Natal, South Africa, while the FUMOZ colony has mosquitoes from southern Mozambique [96,97]. Mosquitoes (20 female and 20 male per species per repeat; three biological repeats) were collected between 0-and 24-h postemergence (here forth called fresh samples). Additionally, 1 ml of each species' larval rearing water was collected (three biological repeats).
Field-caught An. arabiensis were collected between June and August 2019 from Mamfene. These samples had been preserved on silica in microcentrifuge tubes and were retrieved from the departmental archive. Species and Plasmodium-infection status were confirmed using multiplex polymerase chain reaction (PCR) [98][99][100] and enzyme-linked immunosorbent assay (ELISA) [101], respectively, as part of a departmental Sterile Insect Technique (SIT) project. Samples were retrieved after 8, 12 and 16 weeks of preservation.
As mosquito density varies per season, field-collected mosquitoes were low in numbers because June, July and August are the winter months in South Africa (collection numbers are highest in summer and lowest in winter [102]). Mosquito density may have also been exacerbated by the extensive drought in South Africa, which probably left mosquitoes without local breeding pools ([103]: filter between June and August 2019 to view the lack of rainfall). Therefore, only four females were retrieved per time-period, and bacterial identification was performed using culturomics due to the low cost for the limited number of samples. Additionally, as this sample size was low, it was not used to represent the wild mosquito population, and comparisons between field-collected and colonized mosquitoes could not be made.

Mosquito preservation
Colonized mosquitoes were preserved on silica or in RNAlater ® (20 females and 20 males per species, repeat, preservative, and preservation time-period; three biological repeats). Mosquitoes were preserved for 4 and 12 weeks per preservative. Prior to preservation, mosquitoes were immobilized at − 20 °C for 2 min. For silica preservation, mosquitoes were placed individually in 1.5-ml microcentrifuge tubes containing approximately five silica beads (silica gel blue self-indicator (copper sulphatebased), B&M Scientific, South Africa; cat no. CSGB0002) and were separated from silica using a piece of paper (Fig. 1a).
For RNAlater ® preservation, mosquitoes were surface sterilized in 70% ethanol (v/v), left to dry at room temperature on a sterile piece of Kimwipe ® (Kimberly-Clark, TX, USA; cat no. 34155), and individually submerged in 0.5 ml of RNAlater ® solution (Qiagen, Germany; cat no. 76106) in 1.5-ml microcentrifuge tubes (Fig. 1b). Samples preserved in RNAlater ® were stored according to the manufacturer's instructions (samples preserved for 4 weeks were stored at 4 °C, while samples preserved for 12 weeks were placed in RNAlater ® overnight at 4 °C and subsequently stored at − 80 °C). For each species, a set of fresh (non-preserved) samples was collected for comparison. All supplies (silica beads, microcentrifuge tubes, separating paper, etc.) were sterilized prior to use, and swabs of these supplies were cultured to ensure they were not contaminated.

Midgut dissection
Midgut dissections of fresh samples were performed aseptically per protocol by the WHO [104], where 4 µl of phosphate-buffered saline (PBS) was used per midgut; 20 midguts were pooled per species, sex, preservative, preservation time-period, and repeat.

Abdominal dissections of preserved samples
Due to desiccation, mosquito abdomens were shrivelled and brittle. Therefore, mosquitoes were surface sterilized twice in 70% ethanol (v/v), and abdominal segments I to V were dissected and placed in sterile microcentrifuge tubes containing 4 µl of PBS; 20 abdominal segments were pooled per sex, species, preservative, preservation time-period, and repeat. Abdominal segments were then homogenized using a TissueLyser II, followed by centrifugation. Negative controls were set up per group of samples using PBS were carried through during downstream analyses. Additionally, abdomens tore apart easily when mosquitoes were submerged in RNAlater ® and midguts could not be isolated. Thus, mosquitoes were removed from solution and blotted on tissue paper to remove excess RNAlater ® , followed by the adapted dissection method used for silica-preserved samples.

Culture-dependent bacterial identification
Each homogenate (10 µl) was plated on individual selective media agar plates (Table 1). Inoculates were plated and aerobically incubated for a minimum of 16 h at 37 °C. Plates without observable colonies were re-incubated for 16 h to account for slow-growing bacteria. Following incubation, isolates were distinguished morphologically, and distinct colonies were selected for repeated subculture by re-inoculation on fresh primary agar plates and incubation for a minimum of 16 h at 37 °C. Plates without observable colonies were re-incubated for 16 h to account for slow-growing bacteria. Negative controls were set up for each plate type, and during subsequent incubation periods, to ensure plates were not contaminated during incubation. Each plate type was also inoculated with PBS negative controls from dissections.

Mass spectrometry (MALDI-TOF MS)
Each ). An isolate was identified when spectra had a log score value ≥ 1.9 [81]. Every unidentified isolate was tested successively, where a portion of the same colony was placed on a new spot on the target plate and identified as described.

Culture-independent bacterial identification
To supplement culturomics, midguts of fresh mosquitoes and abdomens of preserved mosquitoes (preserved on silica and in RNAlater ® for 4 and 12 weeks) were dissected and pooled (20 female and 20 male per species, preservative, preservation time-period, and repeat; three biological repeats) in sterile PBS as described.
Additionally, 1 ml of each species' larval rearing water was collected (three biological repeats). Bacterial DNA was extracted using the QIAamp ® DNA Microbiome Kit according to the manufacturer's instructions (Qiagen, Germany; cat no. 51704). Negative controls were set up during each stage of DNA extraction, and extraction was also performed on PBS negative controls from abdominal dissections. Prior to sequencing, DNA quality and purity were measured using the NanoDrop ™ 2000c spectrophotometer (Thermo Scientific, MA, USA; cat no. ND-2000C). Due to the high cost of NGS and the high number of negative controls, only experimental samples were sequenced. Nonetheless, negative controls were assessed using the NanoDrop ™ 2000c spectrophotometer and were cultured to ensure no contamination. Samples were sent to Macrogen Europe (Amsterdam, The Netherlands) for 16S rRNA gene sequencing targeting the V3-V4 regions with universal primers, Bakt 341F and Bakt 805R [110]. The Illumina MiSeq system was used to perform paired-end sequencing, and the Fast Length Adjustment of SHort reads (FLASH version 1.2.11) program was used to assemble reads [111]. Pre-processing (denoising) and clustering of sequences were performed with the CD-HIT-OTU and rDnaTools programs [112,113]. Diversity analyses and taxonomy assignments were performed with the Quantitative Insights Into Microbial Ecology (QIIME) program (see Additional file 1 for OTUs with taxonomy assignment) [114].

MALDI-TOF MS Library analysis
According to the MALDI-TOF MS MBT 7854 MSP Library, several bacteria are genetically indistinguishable from one another. Thus, the matching hints section of the library was used to compare MALDI-TOF MS and NGS results (see Additional file 2). Accordingly, indistinguishable bacteria were grouped, and a single bacterium was used to represent indistinguishable bacteria ( Table 2). Representative bacteria were chosen based on which indistinguishable bacteria were also present in NGS results.

Data analyses, diversity indices and statistical analyses
For fresh and field-collected samples, results are presented as accumulative data across replicates. This is because sample sizes for these groups were lower than the sample sizes of preserved mosquitoes due to the inclusion of preservatives and preservation time-periods. There was also a contrast in microbial composition between fresh and preserved samples identified by culture-independent techniques, which may be attributed to potential contamination. Thus, when comparing group, for example, sex (female or male) irrespective of species (An. arabiensis or An. funestus), condition (fresh or preserved), preservative (silica or RNAlater ® ), preservation time-period (4 weeks or 12 weeks) or technique (culture-dependent or culture-independent), only commonly recurring bacteria (bacteria that appeared in at least 50% of replicates in the groups being compared) were reported. Species richness, bacterial diversity, relative abundance, and evenness were calculated. Species richness was measured as the number of operational taxonomic units (OTUs) per sample. Indices were reported per replicate, and the mean index per group was calculated (Additional file 3). Diversity was measured using the Shannon-Wiener diversity index (H): the higher the value of H, the more diverse the community [118]. To estimate relative abundance as a measure of species dominance, Simpson's reciprocal index (1/D) was calculated. 1/D measures the probability that a randomly selected species is the dominant species, where a score of one indicates the community is dominated by a single species [119]. E was used to estimate evenness, which ranges from zero to one with zero signifying no evenness and one signifying complete evenness [120].
Statistical analyses were performed at a 95% confidence interval assuming a 5% level of significance using STATA/IC version 16.1. As the data were not normally distributed (as per Shapiro-Wilk tests), non-parametric statistical analyses were performed. Two-sample Wilcoxon rank-sum (Mann-Whitney) tests were used to determine if diversity indices differed between (i) fresh females and males; (ii) fresh An. arabiensis and An. funestus; (iii) fresh mosquitoes and mosquitoes preserved for 4 weeks; (iv) fresh mosquitoes and mosquitoes preserved for 12 weeks; (v) mosquitoes preserved for 4 weeks and mosquitoes preserved for 12 weeks; (vi) silica preservation and RNAlater ® preservation; and, (vii) culturedependent and culture-independent techniques.

Microbial composition of colonized and field-collected Anopheles arabiensis
Anopheles arabiensis were predominantly colonized by Proteobacteria, irrespective of technique, followed by Bacteroidetes (Fig. 2a). While culture-dependent results estimated Elizabethkingia as the dominant genus in both An. arabiensis sexes, culture-independent results identified Serratia as the dominant genus (Fig. 2b-e).
Culture-dependent results identified Serratia fonticola as a commonly recurring bacterium, but S. fonticola was not identified by NGS. Thus, it is likely that S. fonticola identified by culturomics was actually Serratia oryzae because S. oryzae was identified as a commonly recurring bacterium by NGS and S. oryzae cannot be detected by the MBT 7854 MSP Library. Results of a pairwise alignment (performed using National Center for Table 2 Genetically indistinguishable bacteria grouped in this study A single bacterium was used to represent bacteria that were indistinguishable by the MBT 7854 MSP Library * E. anophelis is misidentified as E. meningoseptica by MALDI-TOF MS [115] ** R. ornithinolytica is misidentified as Klebsiella pneumonia or K. oxytoca [116,117]   An. arabiensis (Table 3; see Additional file 6 for replicate  details). Irrespective of condition, E. anophelis and S. oryzae were recurring bacteria that were common between females and males. Enterobacter cloacae was specific to  [121]: http:// www. ncbi. nlm. nih. gov/ BLAST/) of 16S rRNA sequences revealed 96.75% sequence similarity between S. fonticola (accession number: CP011254.1) and S. oryzae (accession number: NR_157762.1). Meanwhile, a pairwise alignment between S. fonticola and Serratia liquefaciens (accession number: UGYL01000001.1) and S. fonticola and Serratia marcescens (accession number: CP063354.1; the other Serratia species identified by NGS in this study), revealed 94.03% and 85.53% sequence similarity, respectively. Therefore, S. oryzae was used to represent S. fonticola throughout this study.
Aeromonas hydrophila, E. anophelis and S. oryzae were identified as common bacteria between female and male An. arabiensis and their larval rearing water, irrespective of identification technique (Fig. 3). Staphylococcus epidermidis was common between both sexes and the larval rearing water in culture-dependent results, but only in males and the larval rearing water in cultureindependent results. Raoultella ornithinolytica was common between both sexes in culture-dependent results, but only in females in culture-independent results. Some bacteria were identified in a sex-and/or techniquedependent manner (see Additional file 4 for an overview of replicates).
Proteobacteria remained the dominant phylum in preserved An. arabiensis (Additional file 5). Elizabethkingia anophelis and S. oryzae were recurring bacteria that were common in silica-preserved and RNAlater ® -preserved females and S. epidermidis was specific to males (Table 3; see Additional file 7 for replicate details). These sex-specific bacteria were also recurring bacteria in RNAlater ®preserved (and not silica-preserved) samples. Overall, E. anophelis and S. oryzae were common between fresh and preserved female and male An. arabiensis.
Field-collected An. arabiensis were predominantly colonized by Firmicutes (Fig. 4a). The dominant genus in field-collected An. arabiensis was Staphylococcus (Fig. 4b). Although Firmicutes was the dominant phylum in samples preserved for 8 and 12 weeks, samples preserved for 16 weeks were predominated by Proteobacteria (Fig. 4c). Staphylococcus epidermidis and Staphylococcus hominis were common between samples preserved for 8 and 12 weeks while samples preserved for 16 weeks did not have bacteria in common with samples preserved for 8 and 12 weeks (Fig. 4d).

Microbial composition of colonized Anopheles funestus
Anopheles funestus were predominantly colonised by Proteobacteria, irrespective of technique, followed by Bacteroidetes (Fig. 5a). While both identification techniques agree that Serratia is the dominant genus in females, culture-dependent results represent Aeromonas as the dominant genus in males and culture-independent results represent Elizabethkingia as the dominant genus in males (Fig. 5b-e).
Aeromonas hydrophila was identified as a common bacterium between An. funestus females and males and their larval rearing water, irrespective of identification technique (Fig. 6). Serratia oryzae was common between both sexes in culture-dependent results but common in both sexes and the larval rearing water in culture-independent results. Staphylococcus epidermidis was common between both sexes and the larval rearing water in culture-dependent results but only in males and the larval rearing water in culture-independent results. Some bacteria were identified in a sex-and/or techniquedependent manner (see Additional file 8 for an overview of replicates).
Proteobacteria remained the dominant phylum in preserved An. funestus (Additional file 5). Serratia oryzae was a recurring bacterium that was common in silica-preserved and RNAlater ® -preserved An. funestus (Table 4; see Additional file 9 for replicate details). Irrespective of condition, S. oryzae was a recurring bacterium common between females and males. Enterobacter cloacae was specific to females and S. epidermidis was specific to males (Table 4; see Additional file 10 for replicate details). These sex-specific bacteria were also recurring bacteria in RNAlater ® -preserved (and not silica-preserved) samples. Overall, S. oryzae was a common bacterium between fresh and preserved An. funestus.

Comparison of microbial composition between Anopheles arabiensis and Anopheles funestus
Collectively, An. arabiensis and An. funestus were predominantly colonized by Proteobacteria (Additional file 11), irrespective of sex, condition or identification technique (see Additional file 12 for replicate details). Both species were also predominantly colonized by bacteria belonging to the Elizabethkingia and Serratia genera. Elizabethkingia anophelis and S. oryzae were recurring bacteria common between both species irrespective of sex, and E. cloacae was dominant in An. funestus. Enterobacter cloacae was specific to females, while S. epidermidis was specific to males, irrespective of species.

Microbial diversity per sex, species, condition, and technique
Males of both species had higher bacterial species richness than females irrespective of species, condition or identification technique (Table 5). Nonetheless, An. arabiensis and An. funestus had overall comparable species richness.
There were no significant differences amongst diversity indices between female and male mosquitoes (Additional file 13, 6A-C). There were also no significant differences amongst H and E indices between An. arabiensis and An. funestus. There was a significant difference in 1/D between species (P = 0.0209), where An. arabiensis had a higher 1/D than An. funestus (Additional file 13, 6D-F).
Diversity indices were comparable between fresh mosquitoes and mosquitoes preserved for 4 weeks (Additional file 14, 7A-C), fresh mosquitoes and mosquitoes preserved for 12 weeks (Additional file 14, 7D-F), and mosquitoes preserved for 4 and 12 weeks (Additional file 14, 7G-I). There were no significant differences amongst diversity indices between preservatives (Additional file 14, 7J-L). There were, however, significant differences in diversity indices between identification techniques. Culture-independent techniques estimated higher H (P = 0.0200) and 1/D (P = 0.0053) indices than culture-dependent techniques (Fig. 7A-B). Culture-dependent techniques estimated a higher E index (P = 0.0053) than culture-independent techniques (Fig. 7C).

Discussion
Using culture-dependent and culture-independent techniques, this study characterized the microbiota of colonized adult Anopheles preserved either by desiccation in silica or fixation in RNAlater ® . This study echoes notions in the existing literature on the Anopheles microbiota: (i) culture-dependent and culture-independent techniques are complementary; (ii) microbiota are influenced by the mosquito's environment (laboratory vs field); and, (iii) some microbiota are species-and/or sex-specific. Colonized Anopheles in this study had abdomens predominated by Proteobacteria, which is the most common phylum reported in Anopheles studies [2,5,10,18,49,53,63,[122][123][124]. Previous studies have reported Acinetobacter, Bacillus, Corynebacterium, Prevotella, Pseudomonas, Thorsellia, and Veillonella as common genera in colonized Anopheles, while this study reported Elizabethkingia and Serratia as common genera in colonized An. arabiensis and An. funestus [8,18,39,49,88,125].   Firmicutes, which corresponds with studies on the midguts of field-collected Anopheles albimanus from Colombia and Anopheles from Vietnam [53,85]. These studies also identified Staphylococcus and Bacillus as dominant genera colonizing the midguts of field-caught Anopheles, which is consistent with the results of field-collected An. arabiensis in this study.
The dominant bacteria common between colonized An. arabiensis and An. funestus were E. anophelis and S. oryzae. Elizabethkingia anophelis was first isolated from An. gambiae midguts and has been well documented in Anopheles [126][127][128]. Elizabethkingia anophelis protects mosquitoes from infection and positively impacts fecundity in Anopheles [129,130]. As E. anophelis was isolated from An. arabiensis and An. funestus in this study, this bacterium could have similar roles in both vector species. Meanwhile, S. oryzae has not been documented in Anopheles midguts and, thus, its role remains unknown.
Staphylococcus epidermidis was dominant in male Anopheles and has previously been isolated from the midguts of field-collected Anopheles pharoensis and from the salivary glands of colonized An. arabiensis [131]. However, its role is yet to be elucidated. Enterobacter cloacae was dominant in female Anopheles, as well as An. funestus, and is a known Anopheles midgut symbiont [132]. Enterobacter cloacae has been shown to influence vector immunity: E. cloacae affects the development of Plasmodium berghei and P. falciparum in Anopheles stephensi, as well as Plasmodium vivax in An. albimanus [48,133]. Additionally, this bacterium has been tested for paratransgenesis, the genetic modification of symbiotic bacteria to express anti-Plasmodium effector molecules [134,135], in An. stephensi [122]. Thus, E. cloacae may play an immunological role in female Anopheles and in An. funestus. The aforementioned bacteria could also have sex-specific roles in these anopheline species, possibly for the digestion of different food sources (i.e., digestion of blood by females and sugar by males) [16][17][18].
Furthermore, this study shows that neither condition (fresh or preserved), preservative (silica or RNAlater ® ), nor preservation time-period (4 weeks or 12 weeks) influenced microbial composition or diversity. Thus, silica or RNAlater ® are efficient, cost-effective alternatives to previously investigated preservatives; that is, AllProtect reagent, ethanol, and NAP buffer [93].
It was hypothesized that silica desiccation compacts the midgut to form a rigid, secure biofilm around the bacteria; previous studies report that bacterial species encapsulated in silica gel are preserved, along with their biological activity [136,137]. During preservation, the bacteria enter a state of dormancy and when released into saline solution and plated on selective media, the bacteria exit dormancy and acquire nutrients for active growth. Therefore, using silica preservation in combination with culture-dependent techniques is useful because it distinguishes bacteria capable of entering and exiting dormancy, whereas culture-independent techniques cannot distinguish between live and dead bacteria [138]. However, silica preservation may only be suitable for specific bacteria.
It was also hypothesized that preserving mosquitoes in RNAlater ® causes midguts to become pulpous, and RNAlater ® solution may come into contact with the bacteria and inhibit growth when homogenates are placed on nutrient agar because RNAlater ® is bacteriostatic [138,139]. This may account for the few bacteria identified from RNAlater ® -preserved mosquitoes by culture-dependent techniques in comparison to culture-independent techniques. However, since RNAlater ® preserves DNA in high quality, this preservation method is suitable when paired with culture-independent techniques.
Culture-independent techniques identified a richer and more diverse composition of bacteria than culturedependent techniques, which is expected since cultureindependent methods are highly sensitive [138]. As demonstrated, the type of bacterial identification technique affects microbial composition and diversity: culture-independent techniques estimated higher species richness, diversity, relative abundance, and microbial community evenness than culture-dependent technique. Yet, contamination cannot be ruled out.
Nevertheless, the combination of preservatives and identification methods is useful for identifying Anopheles midgut bacteria as it provides a large dataset of overlapping bacteria and can be used for future studies investigating fresh and preserved Anopheles. This may broaden the knowledge on the Anopheles microbial community and could aid future investigations elucidating the role that specific midgut bacteria play in vector species. It could also be used to compare the microbiota of preserved P. falciparum-infected and -uninfected vector species, thereby providing insight into the vector-microbiota-pathogen relationship.
A limit of this study is that the sample size of colonized mosquitoes was larger than the sample size of field-collected mosquitoes as the latter was scarce. Additionally, as there is no way of telling mosquito age upon collection, the age of field-collected Anopheles was unclear. Therefore, since microbiota change during development, the age of field-collected mosquitoes used in this study most likely influenced the bacteria identified [34,140,141]. The conditions that mosquitoes were handled during collection may have also affected the types of bacteria that were preserved (i.e., depending on the time it took for mosquitoes to be immobilised and placed on silica after field collection, this may have affected the bacterial community).
Although culture-independent procedures (DNA extraction and NGS) were performed at the same time for all samples, different generations of samples were collected. This may account for the difference in microbial composition between fresh and preserved samples. As culture-dependent techniques were performed at different times, but with the same generations of mosquitoes, this may have introduced batch effects. Thus, either generational effects, batch effects, contamination, or a combination of these, may account for the lack of uniformity amongst replicates.
Further, pooling does not provide a true representation of the mosquito microbiota because there is high variability between individual mosquitoes [93]. This may have also limited this study. The aforementioned limitations also limit the conclusions, and further investigation (investigating mosquito microbiota individually and increasing the overall sample sizes) is recommended.
As many midgut bacteria are acquired from the environment, identifying preserved microbiota from Diptera in an area can be used to study ecological changes in an environment over time (i.e., if there is a change in an environment, it would be worthwhile investigating if there is also a change in the microbiota of the Diptera inhabiting that area). This could aid in understanding changes in environmental bacteria and the effect that these changes have on the midguts of local Diptera, and on an ecological system as a whole. In addition, extending preservation studies using silica and RNAlater ® to other Diptera may provide insight into Diptera-pathogen relationships and aid studies investigating symbiotic control to reduce disease transmission.