Next Article in Journal
A Long-Term Demographic Analysis of Spotted Turtles (Clemmys guttata) in Illinois Using Matrix Models
Previous Article in Journal
Diversity of Species and Susceptibility Phenotypes toward Commercially Available Fungicides of Cultivable Fungi Colonizing Bones of Ursus spelaeus on Display in Niedźwiedzia Cave (Kletno, Poland)
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Aquatic Hemiptera in Southwest Cameroon: Biodiversity of Potential Reservoirs of Mycobacterium ulcerans and Multiple Wolbachia Sequence Types Revealed by Metagenomics

1
Laboratory for Emerging Infectious Diseases, University of Buea, PO Box 63, Buea, Cameroon
2
Department of Microbiology and Parasitology, University of Buea, PO Box 63, Buea, Cameroon
3
Institute of Infection and Global Health, University of Liverpool, Liverpool L3 5RF, UK
4
Institute of Integrative Biology, University of Liverpool, Liverpool L69 7ZB, UK
5
School of Natural Sciences and Psychology, Liverpool John Moores University, Liverpool L3 5UG, UK
6
Disease Control Department, London School of Hygiene and Tropical Medicine, Keppel Street, London WC1E 7HT, UK
*
Author to whom correspondence should be addressed.
Diversity 2019, 11(12), 225; https://doi.org/10.3390/d11120225
Submission received: 14 October 2019 / Revised: 14 November 2019 / Accepted: 22 November 2019 / Published: 25 November 2019
(This article belongs to the Special Issue Insect Symbiosis)

Abstract

:
Buruli ulcer (BU), caused by Mycobacterium ulcerans, is a neglected tropical disease associated with freshwater habitats. A variety of limnic organisms harbor this pathogen, including aquatic bugs (Hemiptera: Heteroptera), which have been hypothesized to be epidemiologically important reservoirs. Aquatic Hemiptera exhibit high levels of diversity in the tropics, but species identification remains challenging. In this study, we collected aquatic bugs from emerging foci of BU in the Southwest Region of Cameroon, which were identified using morphological and molecular methods. The bugs were screened for mycobacterial DNA and a selection of 20 mycobacteria-positive specimens from the families Gerridae and Veliidae were subjected to next-generation sequencing. Only one individual revealed putative M. ulcerans DNA, but all specimens contained sequences from the widespread alpha-proteobacterial symbiont, Wolbachia. Phylogenetic analysis placed the Wolbachia sequences into supergroups A, B, and F. Circularized mitogenomes were obtained for seven gerrids and two veliids, the first from these families for the African continent. This study suggests that aquatic Hemiptera may have a minor role (if any) in the spread of BU in Southwest Cameroon. Our metagenomic analysis provides new insights into the incursion of Wolbachia into aquatic environments and generated valuable resources to aid molecular taxonomic studies of aquatic Hemiptera.

Graphical Abstract

1. Introduction

Mycobacterium ulcerans is a slow-growing environmental pathogen that infects the skin and subcutaneous tissues, causing Buruli ulcer (BU; also known as Bairnsdale/Daintree ulcer in Australia) [1]. BU is an emerging, treatable but neglected skin infection that manifests as slowly developing, unspecific indolent nodules, papules, or induration by oedema, which can progress to necrotizing skin ulcers [2]. BU treatment outcomes can be good if the condition is diagnosed early, while late diagnosis can lead to extreme patient suffering and severe complications necessitating surgery or even amputation of limb(s) [3,4]. The incubation period of the disease varies from weeks to several months, with a median of four to five months [5,6]. Emergence of BU has been linked to the acquisition of a megaplasmid, pMUM, by Mycobacterium marinum, a far less virulent environmental pathogen associated with aquatic habitats [7]. The megaplasmid encodes for polyketide synthases and accessory enzymes for the biosynthesis of a macrolide toxin, mycolactone [8], which is responsible for tissue damage and local immune suppression in BU [9]. Phylogenomic evidence indicates that all mycolactone-producing mycobacteria (MPM) form a monophyletic group comprising three clonal lineages, which should be considered ecovars of M. ulcerans rather than distinct species (the so-called “Mycobacterium liflandii”, “Mycobacterium pseudoshottsii”, and “Mycobacterium shinshuense”) [7].
BU has been reported in over 35 countries [10], almost half of which are in Africa [2]. In Cameroon, the first cases of BU were diagnosed in 1969 in 47 patients residing in a confined area near the villages of Ayos and Akonlinga [11,12]. This was two decades after BU was first described in Australia [13] and some years after the first case report in Africa [14]. Although the number of BU cases in Cameroon has declined following active surveillance in line with recent global trends [15], it is still a major public health problem in this country. New foci of infection continue to emerge, while Bankim and Akonolinga remain as disease hot spots. From 2001 to 2014, the number of BU endemic health districts rose from two to 64 [15,16,17,18].
Mycobacterium ulcerans DNA has been detected in aquatic bugs (Hemiptera: Heteroptera) [19,20], other limnic organisms (e.g., aquatic plants, fish, tadpoles, and mosquitoes [21,22,23]), and environmental samples (freshwater and soil [24,25]) in several ecological studies, but its route(s) of transmission from these sources has not yet been fully elucidated [26,27,28]. It has been hypothesized that M. ulcerans can be transmitted by human-biting aquatic bugs, notably those in the families Naucoridae, Belostomatidae, Notonectidae, and Nepidae [29,30,31]. Although it has been clearly established that M. ulcerans specifically colonizes the salivary glands of biting bugs and they can transmit the pathogen to laboratory rodents [19,31], their role as disease vectors in nature remains controversial. Thus, while it is possible that bites from aquatic bugs could transmit M. ulcerans to humans [32], there is no evidence that this is epidemiologically relevant on a wide scale [33]. Nevertheless, as most families of aquatic bugs contain species (or morphs) with flight ability [34], the colonization of aquatic bugs by M. ulcerans may be important in the spread of the pathogen into new areas.
Previous studies in central Cameroon have implicated aquatic bugs as reservoirs of M. ulcerans, with colonization of the insects being detected in a BU-endemic area but not a proximate nonendemic area [19]. The current study investigated aquatic bugs from the previously undersampled Southwest Region of Cameroon for signatures of M. ulcerans in order to map environmental sources and potential reservoirs of the pathogen, as evaluating the presence of M. ulcerans in the environment is important in assessing the risks for human infection [26,28,34]. Precise species identification and delimitation of aquatic bugs in tropical regions, where species diversity is very high, can be challenging. Several authors have proposed the use of DNA barcoding targeting a standard region of the mitochondrial gene, cytochrome c oxidase subunit I (coi), in the taxonomy of true bugs [35,36]. Therefore, our study extended the previous use of DNA barcoding for identification of Cameroonian aquatic bugs [37] in a more focused geographical area, where BU appears to be emerging.
Recent reports of a high prevalence of the α-proteobacterial symbiont, Wolbachia, in aquatic insects [38] may also be of epidemiological relevance. Although Wolbachia is the most widespread animal symbiont on Earth, being found in > 50% of terrestrial arthropod species [39], the extent of its penetration into aquatic environments remains unclear. Wolbachia is known to affect the susceptibility of arthropod hosts to colonization by other microorganisms (either positively or negatively, depending on specific combinations of pathogen, vector, and Wolbachia strain [40,41,42]) and can profoundly impede vector competence if artificially introduced into naïve hosts [43]. Therefore, M. ulcerans and Wolbachia may interact within the microbiome of aquatic bugs, either in a facilitative or competitive manner as observed with Wolbachia in other systems [44,45,46], leading to potential impacts on dissemination of the pathogen in the host and perhaps the wider environment.
Here, we present the results of an intensive environmental and taxonomic survey of aquatic bugs in the Southwest Region of Cameroon. Through a combination of traditional morphology-based taxonomy and DNA barcoding coupled with metagenomics, we provide new insights into the diversity of aquatic Hemiptera in Cameroon and acquired the first complete mitochondrial genomes from bugs for the country. While our original goal to characterize environmental strains of M. ulcerans from aquatic bugs in this region was not achievable due to a virtual absence of the pathogen in our specimens, we demonstrate that serendipitously, they harbored surprisingly diverse Wolbachia sequence types.

2. Materials and Methods

2.1. Study Sites

A cross-sectional descriptive study, in which study sites were selected using a multistage sampling approach, was carried out in the Southwest Region of Cameroon. This region was chosen because it is BU-endemic, but very little information on the disease (especially environmental sources of M. ulcerans) exists for this part of the country. Localities endemic for BU have been identified in the Southwest Region, leading to the establishment of a BU detection and treatment center; however, research and control efforts have so far been focused only in the three main BU-endemic foci (located in other administrative regions): Akonlinga (Centre), Ayos (Centre), and Bankim (Adamawa) [16,47].
Six health districts (HDs) from a total of 18 in the Southwest Region were selected for the study based on the cumulative number of BU cases reported within a 14-year period (2001–2014) [15]. Two of the HDs (Mbonge and Ekondo-Titi) are referred to as mesoendemic foci in the present study because they reported cumulative BU cases of 75 and 37, respectively; values much lower than those reported by the three major endemic foci in Cameroon [Akonolinga (1081), Bankim (557), and Ayos (485)] within the same 14-year period [15]. The other four HDs (Limbe, Buea, Muyuka, and Kumba) are referred to as hypoendemic foci because each reported <10 cumulative BU cases within the stated period.
An initial scoping exercise was undertaken in the towns and accessible villages in each selected HD to i) introduce the study to the local population and ii) map the sampling points (also termed as water bodies or sampling sites), which was done with the assistance of the population. Most sites were shallow and slow flowing, but became flooded rapidly during heavy rain during the sampling period, increasing flow rates temporarily (Figure S1A–D)

2.2. Sample Collection and Transportation

Aquatic bugs were collected from each of the 25 water bodies (Figure 1) from 13–25th July 2017, with sampling taking place between the hours of 09:00 and 15:00. As the survey required the collection of representative samples of aquatic bugs in the brief time available between heavy rainfall, and the focus was on molecular screening rather than ecological research, the area or volume of water sampled at each location was not standardized. A team of four workers spent approximately 40 minutes sampling from the water surface to the substrate (a maximum depth of ~1.3 m) at each site using a professional three-piece collapsible hand net (250 mm wide frame, 1 mm in mesh size, and 300 mm deep) (EFE and GB Nets, Lostwithiel, UK). The samples were sorted in a white plastic tray (EFE and GB Nets) in order to remove extraneous vegetable matter and nonhemipteran animals. Putative hemipterans were transferred into a 50 mL polypropylene tube containing absolute acetone, kept in a cooler, and transported to the Laboratory for Emerging Infectious Diseases, University of Buea, for short-term storage. Subsequently, the acetone was decanted and all samples were shipped to the University of Liverpool, UK, under an import license issued by the UK Animal and Plant Health Agency. The acetone was replaced and the insects stored at 4 °C.

2.3. Aquatic Bugs Identification and Morphotyping

Each bug was numbered uniquely for taxa identification. As discussed in prior works, the identification of aquatic Hemiptera in West Africa remains problematic due to high species diversity and incomplete faunistic records [37,48]. Therefore, in case of specific relationships between bug populations and M. ulcerans, we proceeded to classify the material initially into morphotypes to allow retrospective matching between samples expended for molecular assays and those retained as voucher specimens. Using a stereomicroscope, taxonomic keys, written descriptions, specimen comparisons, image verification, and an extensive literature search [37,49,50,51,52,53,54,55,56,57,58,59,60,61,62,63,64,65,66,67,68,69], the aquatic bugs were placed into families, genera, and, wherever possible, species. Confirmation of the identity of selected specimens was done by morphological comparisons with material in the collections of the Natural History Museum, London.

2.4. DNA Extraction from Aquatic Bugs

The bugs were rinsed with distilled water to remove acetone and each specimen was cut into several pieces with a separate sterile disposable scalpel (Swann-Morton, Sheffield, England) before transfer into a 1.5 mL snap cap tube. Ammonium hydroxide (150 μL, 1.5 M) (Sigma-Aldrich) was added to the tube and boiled in a heating block at 100 °C for 15 min [70]. The tubes were centrifuged at 5,000× g for 1 min and reheated at 100 °C for 15 min with the lids open to evaporate ammonia and obtain a final volume of 70–100 µL. Insoluble material was removed after centrifugation at 10,000× g for 10 min and discarded. The concentration of DNA was determined by Quant-iT PicoGreen dsDNA quantification assay using a microplate fluorimeter (Infinite F200, Tecan) and Magellan Data Analysis Software (Tecan). The DNA samples were stored at 4 °C.

2.5. Coi Amplification (DNA Barcoding)

For DNA barcoding, the mitochondrial coi region was amplified. Unless stated otherwise, each conventional PCR run in this study had a final volume of 20 µL, comprising 10 µL BioMix Red master mix (2×), 1 µL of each primer from a 10 µM working stock (final concentration, 0.5 µM), 1 µL DNA template, and 7 µL nuclease-free water. All PCR amplifications were conducted in a Biometra TRIO thermal cycler (Analytik, Jena, Germany) and electrophoretic separations of PCR products were on a 1% agarose gel stained with SYBR Safe DNA gel stain (Invitrogen). The PCR products were visualized and photographed with a Safe Imager transilluminator (Invitrogen). A negative control, in which DNA template was replaced with nuclease-free water, was included in each PCR run.
The primers used for the amplification of the coi fragment (710 bp) were LCO1490: 5′-GGTCAACAAATCATAAAGATATTGG-3′ (sense) and HC02198: 5′-TAAACTTCAGGGTGACCAAAAAATCA-3′ (antisense) [71]. PCR amplifications comprised initial denaturation at 94 °C for 5 min, followed by 35 cycles of denaturation at 94 °C for 1 min, annealing at 50 °C for 1.5 min, and extension at 72 °C for 1 min. The final extension was at 72 °C for 5 min and the reactions were then held at 4 °C. Selected amplicons were directly sequenced in both directions by Macrogen (Seoul, South Korea) using the same primers used for PCR amplification.

2.6. Amplification of a Mycobacterial rpoB Gene Fragment by Quantitative PCR

Quantitative real-time PCR (qPCR) was used for rapid screening to detect and quantify mycobacterial DNA in a subset of aquatic bug specimens (see Results). The qPCR was designed to target a fragment (88 bp) of the RNA polymerase β-subunit (rpoB) gene that is 100% conserved in most nontuberculous mycobacteria, including MPM. The following primers and probes were synthesized by Eurofins Genomics (Ebersberg, Germany): forward primer (Myco-F: 5’-GATCTCCGACGGTGACAAGC-3’), reverse primer (Myco-R: 5’-CAGGAACGGCATGTCCTCG-3’) and Myco-probe (5’-6FAM-ACGGCAACAAGGGCGTCATCGGCAAGATCCT-BHQ-1-3’). A 20 μL reaction mix was prepared containing 1 μL DNA template, 10 μL SYBR Green master mix (2 ×) (Bioline), and final concentrations of 0.5 μM for each primer and 0.25 μM for the probe in PCR-grade water. Amplification proceeded on a CFB-3220 DNA Engine Opticon 2 System (Bio-Rad) at 50 °C for 3 min (holding) and 95 °C for 5 min (initial denaturation), followed by 40 cycles of 20 sec at 95 °C and 40 sec at 60 °C. To quantify the mycobacterial load of each sample, each run included a standard curve with DNA concentrations corresponding to 1,000,000 to 0.1 copies per reaction of a synthetic mycobacterial rpoB oligonucleotide standard (Eurofins Genomics), diluted in 100 ng/μL yeast tRNA (Invitrogen, Carlsbad, CA, USA) to prevent aggregation. Each DNA standard dilution and sample was assayed in duplicate. Opticon Monitor software (v. 3.1) was used for linear regression analysis based on tenfold dilutions of the standard. Each PCR run included a negative control (PCR-grade water) and a positive control (genomic DNA from Mycobacterium ulcerans NCTC 10417, Public Health England Culture Collections, Salisbury, UK).
The frequency of rpoB-positive bugs was compared between meso- and hypoendemic sites using a χ2 analysis at www.socscistatistics.com.

2.7. Confirmatory Conventional PCRs for Bacteria

For specimens positive with the rpoB qPCR, a 400 bp fragment of the IS2404 insertion sequence from MPM was amplified in a first-round PCR, while the nested round amplified a 150–200 bp fragment. Both rounds of PCR were performed with primers as previously described [72]. The PCR reactions were prepared as described above for coi, except for the nested round, in which 1 µL PCR product from the first-round PCR was used as a DNA template. The PCR amplifications were carried out under the following conditions: initial denaturation at 95 °C for 15 min followed by 25 cycles of denaturation at 94 °C for 30 sec, annealing at 64 °C for 1 min, and extension at 72 °C for 1 min. The final extension was at 72 °C for 10 min. The cycling parameters for the nested round were identical to those of the first round, except the number of cycles was raised to 30.
To amplify Wolbachia surface protein (wsp), primers wsp-81F and wsp-691R [73] were used at a final concentration of 0.5 µM in 20 µl reactions with BioMix Red master mix (see coi PCR above). The thermocycling conditions comprised 35 cycles at 94 °C for 1 min, 55 °C for 1 min, and 72 °C for 1 min.

2.8. Metagenomic Sequencing

Twenty gerrid and veliid specimens with the highest copy numbers of mycobacterial rpoB were subjected to next-generation sequencing on an Illumina HiSeq 4,000 platform (Illumina, Inc., San Diego, CA, USA). Illumina fragment libraries were prepared from genomic DNA using the Nextera XT kit as previously described [74], and sequenced as paired-ends (2 × 150 bp) over two lanes, with a minimum of 280 million clusters per lane. Base-calling, demultiplexing of indexed reads, and trimming of adapter sequences was conducted as reported previously [74]. Raw data were submitted to the Sequence Read Archive at the National Center for Biotechnology Information (NCBI) under BioProject PRJNA542604 and accession numbers SRX5884614–SRX5884633.

2.9. Taxonomic Assignment

Taxonomic read assignments were performed using k-mer-based classifiers Kraken v. 2.0.7 [75] and CLARK-S v. 1.2.5 [76] with default settings. The output from Kraken2 or CLARK-S was visualized using the Krona v. 2.7 [77] tool, creating hierarchical interactive pie charts. MetaPhlAn2 v. 2.6.0 [78] was also used to identify bacterial species and estimate their relative abundance across all 20 samples according to the database ‘mpa_v20_m200’. Velvet v. 1.2.10 [79] wrapped by VelvetOptimiser v. 2.2.6 was used to assemble the 20 bug genomes to contig level. According to these assemblies, GC-coverage plots (proportion of GC bases and Velvet node coverage) were generated using BlobTools v. 1.0.1 [80]. Taxonomic assignment of these contigs was achieved using megablastn to search against the nr database.

2.10. Mitogenome Assembly and Annotation

The mitogenome of each sample was assembled de novo using NOVOPlasty v. 2.7.2 [81], with the D. melanogaster mitogenome as a seed sequence (NC_024511.2). Each circularized mitogenome assembly obtained in the first run was used iteratively as a seed sequence to run NOVOPlasty again. In all the assembly processes, the k-mer length was set to 39 and other operation parameters in the configuration file were left on the default setting or set to fit the features of the reads. Protein-coding, rRNA, and tRNA genes were annotated with Mitos2 [82] by searching the Refseq 81 Metazoa database using the invertebrate mitochondrial genetic code. The protein-coding genes were manually added or edited, if necessary, after comparing with other hemipteran mitogenomes via tblastn. The tRNAs were identified with the tRNAscan-SE search server v. 2.0 [83] by setting sequence source to ‘other mitochondrial’ and genetic code for tRNA isotype prediction to ‘invertebrate mito’. Mitogenome sequences were compared using BLAST Ring Image Generator v. 0.95 [84] to Perittopus sp. (JQ910988.1) [85], Aquarius paludum (NC_012841.1) [86], and Gigantometra gigas (NC_041084.1) [87] reference sequences. Gene rearrangement analysis of the assembled mitogenomes were visualized using the Mauve genome aligner v. 25 February 2015 [88]. MrBayes v. 3.2.6 [89] was run with the generation parameter set to ‘1 in 100 samples’ and the substitution model set to ‘GTR + G + I’ (the first 25% of samples were discarded as burn-in) to build a mitogenome tree based on protein superalignments, as previously described [90]. The assembled mitogenomes were submitted to NCBI under accession numbers MN027271–MN027279.

2.11. Barcode Analysis of coi Sequences and Phylogenetics

All Sanger-sequenced coi reads were trimmed of ~50 poor-quality bases at the 5’ and 3’ end before assembly using the cap3 [91] program. The assemblies were then submitted to the Barcode of Life Data (BOLD) System [92] v. 4 under the project name “CBUG” and analyzed using the “Cluster Sequences” feature, which calculates pairwise distance following alignment with an amino acid-based hidden Markov model. Gaps were handled by pairwise deletion and the sequences were binned into operational taxonomic units (OTUs) by nearest-neighbor distance. The coi sequences were compared with reference gerrid, veliid, and mesoveliid sequences from NCBI and from a previous Cameroonian barcoding study of aquatic bugs [37], in which the data were available in their Supplementary Information but not BOLD Systems or NCBI. Alignments were performed using Mafft v. 7.402 [93] with the ‘--auto’ option and automatically trimmed with Gblocks v. 0.91b [94]. A phylogenetic tree for coi sequences was reconstructed using IQ-TREE v. 1.6.7 [95] with the parameters ‘-m MFP -alrt 1000 -bnni -bb 1000’ to determine the best-fit model for ultrafast bootstrap, and both Shimodaira–Hasegawa-like (SH)-like approximate likelihood ratio tests and ultrafast bootstrapping were conducted to access branch supports within one single run.

2.12. Analysis of Wolbachia Genes and Phylogenetics

Wolbachia reads assigned by Kraken2 were further assembled and scaffolded using SPAdes-3.7.1 [96] with the ‘--careful’ function. Wolbachia protein-coding genes in these assemblies were predicted using Prokka v1.13 [97] with default settings. One-to-one orthologs of Wolbachia protein-coding genes in five reference Wolbachia strains (GCF_000008025.1 GCF_000073005.1, GCF_000306885.1, GCF_000008385.1, GCF_001931755.2, and GCF_000829315.1), each from a different taxonomic “supergroup” (higher-level clade), were identified by using OrthoFinder 2.2.3 [98]. To investigate the Wolbachia supergroups in our samples, concatenated phylogenetic trees were constructed for multiple orthologs (depending on the recovery of Wolbachia protein-coding genes in each assembly) at the amino acid level using IQ-TREE as described above. A phylogenetic analysis for the Wolbachia wsp gene (trimmed PCR amplicons) was also performed using MrBayes with the same parameters as for the whole mitogenome tree. Velvet-assembled contigs classified as Wolbachia bacteriophage (phage WO) were aligned onto a complete phage genome (KX522565.1) from wVitA, a symbiont of the parasitoid wasp Nasonia vitripennis [99], to identify the presence of phage in sequenced bug samples.

3. Results

3.1. Distribution of Aquatic Bugs in the Health Districts

Of the 1,113 aquatic insects collected in the Southwest Region, 1,102 were identified as true bugs (Hemiptera) belonging to eight families and 29 putative genera. Specimens in each genus were further separated into morphotypes, generating a total of 110 (of which 34 were represented by only one specimen each). Subsequently, 331 individual specimens were retained for species confirmation and voucher specimens, while the remaining 771 were subjected to molecular analysis. Of the 25 sampling points, 17 (68%) were in hypoendemic zones, while the remaing eight (32%) were in mesoendemic areas. Similarly, of the 1,102 aquatic bugs captured, 714 (64.8%) were from the hypoendemic sites and 388 (35.2%) from the mesoendemic sites. The distribution of the aquatic bugs by HD is shown in Figure 2.

3.2. Composition of Aquatic Bug Taxa

Most bugs collected belonged to just two families, with 498 (45.2%) from the Gerridae (water striders) and 426 (38.7%) from the Veliidae (riffle bugs). The Notonectidae, Naucoridae, Belostomatidae, Mesoveliidae, Hydrometridae, and Nepidae were represented by 54 (4.9%), 52 (4.7%), 36 (3.3%), 18 (1.6%), 11 (1.0%), and 7 (0.6%) specimens, respectively (Figure 3). Overall, the population of aquatic bugs originating from families known to bite humans was very low (149/1,102, 15.5%) and distributed in the HDs according to Table S1.
Visual inspection of the morphotypes indicated that Gerridae were dominated by Limnogonus spp. and the Veliidae by Rhagovelia spp. (Table 1). The marked morphological diversity of these two genera, highlighting the complexity of the taxonomic assignment, is illustrated in Figure 4. Overall, the specimens were allocated to 42 putative species, some containing multiple morphotypes (Table 1). The family Gerridae was the most diverse, with 18 putative species, followed by Naucoridae (8 putative species), Notonectidae (6 putative species), Veliidae (4 putative species), Belostomatidae (3 putative species), and one putative species each for the families Mesoveliidae, Hydrometridae, and Nepidae. The Gerridae, Veliidae and Hydrometridae were present in all HDs (Table S2); indeed, six sampling points (three in Buea HD and one each in Muyuka, Kumba, and Mbonge) had only gerrids and veliids.
In the first round of M. ulcerans screening by rpoB qPCR, the human-biting bugs (Table S1) were not analyzed for the detection of M. ulcerans due to the high intraspecific morphological polymorphism observed and their relative scarcity. For example, the 54 notonectids, placed in six putative species, comprised 16 morphotypes with the most abundant morphotype having only five specimens.
Similarly, bugs in the families Mesoveliidae (18, 1.6%) and Hydrometridae (11, 1.0%) were also not analyzed for M. ulcerans carriage because of their small numbers (Table 1), although the former were not always easily differentiated from veliids (see Discussion).

3.3. Prevalence of Mycobacterial DNA in the Aquatic Bugs

The mycobacterial rpoB gene was detected and quantified in 8.9% of the 771 specimens, comprising 13.2% (53/405) from Gerridae and 4.4% (16/366) from Veliidae. The distribution of the positive specimens among the putative species is detailed in Table 2. The positive specimens comprised 7.9% (21/266) from mesoendemic sites and 9.5% (48/505) from hypoendemic sites (p = 0.457). None of the gerrid and veliid specimens positive for the mycobacterial rpoB gene were positive for the M. ulcerans-specific IS2404 sequence, although this region was successfully amplified from the positive control. To ensure potentially positive bug specimens had not been missed, DNA was also extracted from the 149 biting bugs (Table S1) and then assayed by the IS2404-nested PCR. Mycobacterium ulcerans DNA was not detected in any of these samples.

3.4. Bacterial Sequences in Gerrid and Veliid Metagenomic Datasets

As both epineustonic and nektonic bugs have been reported to have high infection rates with M. ulcerans in other BU-endemic regions of Cameroon (see Discussion), we applied next-generation genomic sequencing to the 20 bug specimens with the highest rpoB copy numbers to confirm that they were infected only with non-MPM species. Approximately 11–27 million reads were obtained per specimen (Table 3) and taxon-annotated GC-coverage plots (“blob analysis”) revealed that several of the samples contained apparent bacterial sequences at a relatively high coverage rate (Figure S2). In most samples, the CLARK-S tool identified ~10-fold fewer bacterial k-mers than did Kraken2 (data not shown), so was not used for downstream analyses of bacterial sequences. K-mers identified as originating from Mycobacterium spp. by Kraken2 exhibited low abundance (<0.5% of all bacterial read pairs) in all specimens except for a Metrocoris sp. (Gerridae) from Muyuka, in which mycobacterial sequences comprised almost 5% of all bacterial read pairs (Table 3). However, 98% of these mycobacterial k-mers were classified as deriving from Mycobacterium sp. WY10, a non-MPM soil-associated species [100], with no sequences classified as M. ulcerans in this specimen. Only one specimen, a Limnogonus hypoleucus from Limbe, contained M. ulcerans sequences at a proportion of >0.1% of all bacterial sequences (84 read pairs), accounting for ~50% of mycobacterial sequences (Table 3).
Strikingly, read pairs classified as Wolbachia were more abundant than mycobacterial sequences in 16/20 specimens, and in many cases, the differential in read pair count between the two bacterial genera was 2–3 logs (Table 3). However, Wolbachia read pairs accounted for a very wide range of bacterial sequences (0.05%–63.08%) between specimens. Using a conventional PCR targeting wsp, we were able to amplify this Wolbachia-specific gene from all 12 specimens with >6,000 Wolbachia read pairs, but from none of the samples with a lower Wolbachia DNA content (Table 3).
In order to classify the Wolbachia strains in the bug specimens, we applied three approaches. First, Sanger sequencing of the wsp amplicons and a Bayesian phylogenetic analysis suggested that most Wolbachia infections belonged to supergroup B, represented in the tree by strain wPip (Figure 5). Two specimens apparently deviated from this pattern, with a supergroup A sequence related to wMel being amplified from Rhagovelia sp. 2, while a more divergent wsp sequence was recovered from Microvelia sp. 1, which was nested between representative sequences from supergroups F and D (Figure 5). Second, as wsp is known to recombine [101], especially between supergroups A and B that coinfect numerous arthropod hosts, we also built trees from protein-coding Wolbachia orthologues recovered from the metagenomic datasets. These supported the supergroup assignments suggested by the wsp analysis in all cases, with the Wolbachia sequences from Microvelia sp. 1 indicating a placement in supergroup F (Figure S3).
Third, the metagenomic phylogenetic analysis tool MetaPhlAn2 identified matches with supergroup B Wolbachia genomes in almost all samples with amplifiable wsp, although a strong signal assigned to a supergroup A genome (wMelPop) was apparent in the Rhagovelia sp. 2 specimen (Figure 6). For the Microvelia sp. 1 sample, Wolbachia reads were detected but not classified as being associated with a specific sequenced strain. The classification of the Wolbachia reads by MetaPhlAn2 suggested that more than one strain might be present in several samples, especially the Rhagovelia spp (Figure 6). Finally, it was noteworthy that sequences matching the bacteriophage of Wolbachia (phage WO), which is known to carry the genes responsible for Wolbachia-mediated reproductive manipulations in arthropods [102], were identified on numerous contigs from the Rhagovelia sp. 2 specimen (Figure S4).

3.5. Assembly of Complete Mitogenomes from Gerrids and Veliids and DNA Barcoding

Of the 20 bug specimens subjected to next-generation sequencing, it was possible to assemble circularized mitogenomes in nine cases, of which only two were from veliids. Annotation of these genomes showed a gene order in both bug families that was typical of insects (Figure S5). The veliid mitogenomes appeared to be approximately 100 bp shorter than those from gerrids (Figure S5, Figure 7). However, due to a well-recognized limitation of Illumina sequencing technology [103], it was apparent that we achieved only low coverage of the control region (AT-rich region; see breaks in Figure 7 at ~15 kbp) compared with three available references from gerrids and veliids, especially for Metrocoris and Trepobates. Although putative differences in this region between species must be confirmed using alternative sequencing technologies, it appears to be potentially useful for discriminating between Limnogonus spp. (Figure 7A,B).
To evaluate the concordance between morphological identifications and DNA barcodes, 50 gerrid, veliid, and mesoveliid samples that produced bright coi amplicons by PCR were sequenced, and 43 samples generated unambiguous alignments of forward and reverse reads. These were supplemented by an additional nine coi sequences obtained from the Illumina sequencing effort, producing a total of 52 sequences across 20 morphotypes. Twenty-one Sanger-sequenced samples were flagged as containing stop codons by BOLD Systems but were retained for downstream analyses. Sequence clustering by BOLD Systems produced 19 putative OTUs. A phylogenetic analysis of the coi sequences from gerrids, veliids, and mesoveliids, incorporating sequences from the study of Ebong et al. and top BLASTn hits from NCBI, revealed incomplete resolution of bug families [e.g., note the position of “V5” veliid sequences from Ebong et al. and references for Stridulivelia spp., Microvelia americana, and Perittopus horvathi (all Veliidae) on the tree; Figure 8]. However, clustering of gerrid genera from our study (Trepobates, Eurymetropis, Metrocoris, Metrobates, and Limnogonus) and other species (L. hypoleucus, L. cereiventris, and to a large extent, L. intermedius) was clearly evident. The inclusion of barcodes containing stop codons frequently, but not always, generated spurious OTUs containing a single specimen in the BOLD cluster analysis as highlighted previously (Figure 8) [104]. This had a comparatively greater impact on the analysis of the smaller veliid dataset, although a more likely explanation for the aberrant placement of our Mesovelia specimens is specimen misidentification, since support for their positioning within a clade of Rhagovelia spp. was strong (Figure 8). In summary, the 20 morphotypes were collapsed into 12 high-confidence OTUs, once those that contained only sequences with stop codons (seven OTUs) were disqualified.
With the exception of the “V5” sequences from the Ebong et al. (2016) analysis [37] and our “Mesovelia” specimens, clustering at the family level in both studies was consistent (Figure 8). A phylogenetic analysis at the whole mitogenome level also clearly separated gerrid and veliid specimens (including references from NCBI) into family-specific groups, while Limnogonus spp. were resolved as a monophyletic clade with unambiguous segregation by species (Figure 9). The “V4” sequences from Ebong et al. [37] appear to belong to our Rhagovelia spp. OTU 5, and the “Ge6” sequences identified as L. hypoleucus in Ebong et al. [37] were positioned proximal to our own L. hypoleucus (OTU 10) (Figure 8). However, we did not identify any Angilia spp. (“V1”) in our study (Figure 8). Additional comparisons between the two works were not possible, as Ebong et al. [37] did not provide putative identifications for several sequence clusters and their data are not held in BOLD Systems. Nevertheless, it is striking that despite the much greater geographic coverage of the Ebong et al. study [37], extending across most of Cameroon’s administrative regions, they reported the initial identification of five gerrid genera, of which only three (all observed by us) were corroborated by DNA barcoding in their analysis. In our study located within one region, eight gerrid genera were recorded, of which five were selected for DNA barcoding and received strong phylogenetic support. Four of these genera (Trepobates, Eurymetropis, Metrocoris, and Metrobates) were not observed in the Ebong et al. study [37].

4. Discussion

The findings of our study contrast markedly with others from Cameroon in unearthing little evidence for the presence of M. ulcerans in aquatic bugs. In order to determine if this is a genuine insight into distinctive ecological characteristics of the pathogen in Southwest Cameroon, the possibility of artefacts or biases arising from our methods used in the field and/or laboratory need to be considered. We sought to screen the bug specimens rapidly with a general mycobacterial rpoB qPCR assay that allowed us to assess bacterial load, not only qualitative presence, for downstream next-generation sequencing. This was originally intended to provide data on the environmental strain(s) of M. ulcerans in the Southwest Region. However, the fact that the rpoB-positive specimens were negative with the nested IS2404 assay, which targets a highly repeated insertion sequence [7], strongly suggested that the bugs contained non-MPM DNA. This lack of M. ulcerans was confirmed in most cases following whole genome sequencing of 20 bug specimens with relatively high rpoB copy numbers. Thus, it seems unlikely that M. ulcerans infections in our collection were missed.
In common with other surveys of M. ulcerans infection in aquatic invertebrates [19,33,105], we did not attempt to remove external contamination from the insects’ surface. There were two reasons for this: (a) in systematic studies, “surface sterilization” of insects has been shown to have no significant impact on the microbiome profiles determined by next-generation sequencing [106]; and (b) in predatory species, the impact of ingested organisms within gut contents is likely to have a much larger influence on microbiome data than surface contamination. The latter would have required painstaking dissection of hundreds of fresh specimens to mitigate even in part. Thus, although ~80 read pairs were classified as originating from M. ulcerans in one gerrid from a hypoendemic site (Limbe), this extremely low signal (not confirmed by the IS2404 assay) could originate from degraded M. ulcerans DNA from gut contents or nucleic acids bound to the insect cuticle. Nevertheless, it is suggestive of the presence of M. ulcerans at that site.
Since the laboratory procedures between our study and those reporting substantial infection rates of M. ulcerans in aquatic bugs were similar, then potential differences in the field sampling strategy should be considered. The most obvious is that our study was much smaller than those of Marion et al. [19] and Garchitorena et al. [105] in terms of the total number of bugs collected, and we assayed insects individually rather than in small pools. The previous studies were also undertaken in highly endemic areas for BU, with >10-fold more recorded human cases than even our “mesoendemic” sites. The influence of seasonality is also potentially important, although in areas of high endemicity, the positivity rate for M. ulcerans has been reported to be greatest in July to August [19,105], which was a similar period to our field survey in the Southwest.
A further important difference between our study and prior screens of aquatic bugs for M. ulcerans was evident in the dominant hemipteran families found at our sampling sites. Following the publication of the hypothesis that biting bugs are vectors of M. ulcerans [29], subsequent studies have tended to focus on biting taxa at the expense of nonbiting bugs, whereas we paid greater attention to gerrids and veliids, as they were much more abundant in our survey. A recent ecological analysis linked infection of aquatic Hemiptera with M. ulcerans to the nektonic lifestyle of predatory species that feed on animals their own size or even larger, occupy aquatic vegetation, and live towards the bottom of the water column [34]. However, it is possible that these conclusions are a result of confirmation bias resulting from an emphasis in the literature on biting species in the epidemiology of BU. In accordance with this interpretation, when gerrids, and even phytophagous taxa such as the Corixidae, have been screened previously for M. ulcerans in Cameroon (albeit with smaller sample sizes than biting bugs), they have been found to be positive at rates equivalent or higher than those of other bug taxa [19]. Moreover, the first isolation and characterization of the pathogen from a nonclinical source was from a gerrid in Benin [20], indicating that this aquatic bug family is at least as likely to have a symbiotic relationship with M. ulcerans as any other. Data from veliids in the BU literature are too scant to draw conclusions [33], but as they have a similar ecology to gerrids, there is no a priori reason to believe that they cannot be colonized by M. ulcerans. These considerations all emphasize that the lack of evidence for M. ulcerans in aquatic Hemiptera in Southwest Cameroon is not the result of methodological flaws; indeed, several studies on a much larger scale have fundamentally challenged the hypothesis that these insects have any substantive role in BU epidemiology [26,33].
The sequencing of the bug specimens in our study provided an opportunity to investigate other components of the microbiome in these relatively little-studied insects. In this small sample, Wolbachia was found not only to be the dominant symbiont but also to be represented by multiple sequence types. While Wolbachia is ubiquitous in arthropods from terrestrial environments, very few data on its distribution in aquatic species existed until recently, apart from analyses of holometabolous, medically important vectors with aquatic larval stages (principally mosquitoes [107], but also blackflies [108,109], and some limited surveys in aquatic crustaceans [110,111]). However, a 2017 meta-analysis of Wolbachia prevalence in aquatic insects (n = 228 species) estimated that >50% of species are infected, with Hemiptera exhibiting a particularly high prevalence (69%) [38]. Importantly, assessing Wolbachia infection by molecular methods alone (including metagenomics) leads to potential pitfalls in interpretation, primarily due to high rates of lateral gene transfer from Wolbachia into the host genome, which can lead to Wolbachia relics or “genomic fossils” in species that may have lost the infection [112,113,114]. Host genomes containing sequences from more than one supergroup have even been reported [115]. Therefore, our putative Wolbachia infections in Cameroonian gerrids and veliids require confirmation by microscopy for whole bacteria or deeper sequencing to recover closed bacterial genomes.
Whether all of the Wolbachia infections in the aquatic bugs are extant or not, in our very small sample of fully sequenced specimens, evidence for three Wolbachia supergroups was apparent. Supergroups A and B are the most prevalent clades in arthropods worldwide, forming the basis of the so-called Wolbachia “pandemic” [116,117]. Coinfection by both clades in a single host has been described repeatedly [118,119], with some evidence for this in our own specimens. The presence of supergroup F-like sequences in Microvelia sp. 1 was a more unusual finding, as this clade is less widespread than the A and B supergroups, although it is also the only one known from both arthropod and nematode hosts [120,121]. As the Wolbachia signal in this specimen was relatively low (~6,000 read pairs, or 0.04% of the total), the possibility that the clade F sequences originate from gut contents or an arthropod or nematode parasite or parasitoid cannot be discounted at this stage [122,123]. A priority for future research on the association of aquatic bugs with M. ulcerans will be to determine if Wolbachia influences this relationship; e.g., whether its presence might explain the limited evidence of M. ulcerans infection in our specimens compared with other locations in Cameroon. Notably, while most research on pathogen protection by Wolbachia in arthropods has focused on artificial infections for vector-borne disease control [43], Wolbachia has also been reported to influence susceptibility to pathogens in some natural systems [41,42].
A final contribution of our study was to extend the integrative taxonomy of aquatic Hemiptera in Cameroon following a previous published work [37]. The Ebong et al. study [37] included a much greater number of specimens used for DNA barcoding (188) than our own from across the whole country. However, our intensive field survey in the Southwest produced a similar number of barcodes to Ebong et al. [37] for the Gerridae and Veliidae and they were selected from a larger sample size for these taxa. We obtained barcodes from four gerrid genera not observed in the Ebong et al. [37] study (Trepobates, Eurymetropis, Metrocoris, and Metrobates) and complete mitogenomes for Trepobates, Metrocoris, and Limnogonus; the first to be published worldwide. The barcoding effort showed some limitations, particularly with respect to the amplification of putative nuclear mitochondrial transfers (which accumulate mutations) and/or the impact of high rates of heteroplasmy [104]. However, with stringent interpretation of barcode segregation, biases caused by these artifacts can be minimized, and the effects of mixing Sanger and next-generation sequencing technologies should also be considered carefully [87].
We noted that separation of aquatic bug families by coi sequences was not perfect in all cases, especially when considering reference sequences deposited in NCBI. Misidentifications are an obvious explanation for this discrepancy, and probably also explain the aberrant placement of Mesovelia spp. in our study. Unfortunately, although previous studies using a larger number of characters (e.g., whole mitogenomes) have been published for Hemiptera, representation of aquatic taxa has been patchy, with no inclusion of mesoveliid specimens for instance [85,86]. While our whole mitogenome analysis for the small number of gerrid and veliid assemblies available separated the two families unambiguously, future efforts should focus on obtaining complete mitogenomes from mesoveliids too.
At the genus and species level, the coi barcoding proved its usefulness for gerrids in particular, but sampling of a more diverse range of genera will be required before its utility with veliid specimens can be firmly established. A very large barcoding effort for Hemiptera found that the minimum interspecific distance was >3% for 77% of congeneric species pairs, although a case of identical barcodes between different genera of mirid bugs was uncovered [36]. An important consideration in species infected with symbionts that can cause cytoplasmic incompatibility, such as Wolbachia, is introgression of mitochondrial haplotypes following symbiont-driven selective sweeps [124]. In some cases, this introgression can penetrate across hybrid zones during breakdown of reproductive barriers [125,126,127]. In this respect, it is noteworthy that a phylogenetic analysis of Limnogonus spp. using two mitochondrial and one nuclear marker found the two subgenera (Limnogonus sensu stricto and the exclusively Afrotropical Limnogonoides) to be paraphyletic, as well as weak bootstrap support for clades and species within Limnogonoides [48]. Although L. guttatus was not included in that analysis, it is thought to be closely related to L. hypoleucus and L. poissoni [48], and the extent to which these species might be able to hybridize in nature, potentially sharing Wolbachia strains, should be explored in future.

5. Conclusions

In this survey of M. ulcerans infection of aquatic Hemiptera in Southwest Cameroon, an emerging focus of BU, almost no evidence for a role for these insects as an environmental reservoir of the pathogen was uncovered. Therefore, we recommend that future epidemiological studies in the region sample more broadly from a variety of animal, plant, and inanimate sources. As aquatic Hemiptera from Cameroon are likely to remain of interest for both basic biodiversity studies and comparative ecological analyses on the distribution of M. ulcerans in the environment across West Africa, the role of Wolbachia in these species, and its potential impacts on hemipteran population structure and microbiomes, should be evaluated further.

Supplementary Materials

The following are available online at https://www.mdpi.com/1424-2818/11/12/225/s1, Figure S1: Representative sampling points, Figure S2: Representative blob plot from sample 15 (Rhagovelia sp. 2 from Mbonge), Figure S3: Concatenated tree based on 26 protein-coding Wolbachia orthologs for sample 5, Figure S4: Contigs from sample 15 (Rhagovelia sp. 2 from Mbonge) mapping to a reference bacteriophage of Wolbachia (WO phage genome), Figure S5: Alignment of two veliid (Rhagovelia spp.) and seven gerrid (Limnogonus, Metrocoris and Trepobates spp.) mitogenomes, Table S1: The distribution of human-biting bugs captured in the health districts. Table S2: Distribution of aquatic Hemiptera in the Southwest Region of Cameroon.

Author Contributions

Conceptualization, S.N.E, R.J.P., L.M.N., A.C.D., and B.L.M.; methodology, S.N.E., X.D., R.J.P., A.C.D., and B.L.M.; validation, S.N.E., R.J.P., X.D., and B.L.M.; formal analysis, S.N.E., X.D., and B.L.M.; investigation, S.N.E., A.J.K., C.S.H., R.J.P., and B.L.M.; resources, S.N.E., A.C.D., L.M.N., R.N.N., and B.L.M.; data curation, S.N.E., X.D., and B.L.M.; writing—original draft preparation, S.N.E., X.D., and B.L.M.; writing—review and editing, S.N.E., A.C.D., R.J.P., X.D., and B.L.M.; visualization, S.N.E., X.D., and B.L.M.; supervision, S.N.E., L.M.N., R.N.N., R.J.P., and B.L.M.; project administration, S.N.E., L.M.N., R.M.N., and B.L.M.; funding acquisition, S.N.E., L.M.N., and B.L.M.

Funding

This work was funded by the Medical Research Foundation, UK, under an Africa Research Excellence Fund (AREF) Research Development Fellowship grant (number ESEMUMRF-157-0010-F-ESEMU) awarded to S.N.E. We thank Alan Bannister (University of Liverpool) for assistance with the photography of the bugs (Figure 4).

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Tai, A.Y.C.; Athan, E.; Friedman, N.D.; Hughes, A.; Walton, A.; O’Brien, D.P. Increased severity and spread of Mycobacterium ulcerans, Southeastern Australia. Emerg. Infect. Dis. 2018, 24, 58–64. [Google Scholar] [CrossRef]
  2. Röltgen, K.; Pluschke, G. Epidemiology and disease burden of Buruli ulcer: A review. Res. Rep. Trop. Med. 2015, 6, 59–73. [Google Scholar]
  3. Toutous Trellu, L.; Nkemenang, P.; Comte, E.; Ehounou, G.; Atangana, P.; Mboua, D.J.; Rusch, B.; Njih Tabah, E.; Etard, J.F.; Mueller, Y.K. Differential diagnosis of skin ulcers in a Mycobacterium ulcerans endemic area: Data from a prospective study in Cameroon. PLoS Negl. Trop. Dis. 2016, 10, e0004385. [Google Scholar] [CrossRef] [PubMed]
  4. Jacobsen, K.H.; Padgett, J.J. Risk factors for Mycobacterium ulcerans infection. Int. J. Infect. Dis. 2010, 14, e677–e681. [Google Scholar] [CrossRef] [PubMed]
  5. Loftus, M.J.; Trubiano, J.A.; Tay, E.L.; Lavender, C.J.; Globan, M.; Fyfe, J.A.M.; Johnson, P.D.R. The incubation period of Buruli ulcer (Mycobacterium ulcerans infection) in Victoria, Australia remains similar despite changing geographic distribution of disease. PLoS Negl. Trop. Dis. 2018, 12, e0006323. [Google Scholar] [CrossRef] [PubMed]
  6. Trubiano, J.A.; Lavender, C.J.; Fyfe, J.A.; Bittmann, S.; Johnson, P.D. The incubation period of Buruli ulcer (Mycobacterium ulcerans infection). PLoS Negl. Trop. Dis. 2013, 7, e2463. [Google Scholar] [CrossRef] [PubMed]
  7. Doig, K.D.; Holt, K.E.; Fyfe, J.A.; Lavender, C.J.; Eddyani, M.; Portaels, F.; Yeboah-Manu, D.; Pluschke, G.; Seemann, T.; Stinear, T.P. On the origin of Mycobacterium ulcerans, the causative agent of Buruli ulcer. BMC Genom. 2012, 13, 258. [Google Scholar] [CrossRef] [PubMed]
  8. George, K.M.; Chatterjee, D.; Gunawardana, G.; Welty, D.; Hayman, J.; Lee, R.; Small, P.L. Mycolactone: A polyketide toxin from Mycobacterium ulcerans required for virulence. Science 1999, 283, 854–857. [Google Scholar] [CrossRef]
  9. Fraga, A.G.; Cruz, A.; Martins, T.G.; Torrado, E.; Saraiva, M.; Pereira, D.R.; Meyers, W.M.; Portaels, F.; Silva, M.T.; Castro, A.G.; et al. Mycobacterium ulcerans triggers T-cell immunity followed by local and regional but not systemic immunosuppression. Infect. Immun. 2011, 79, 421–430. [Google Scholar] [CrossRef]
  10. Simpson, H.; Deribe, K.; Tabah, E.N.; Peters, A.; Maman, I.; Frimpong, M.; Ampadu, E.; Phillips, R.; Saunderson, P.; Pullan, R.L.; et al. Mapping the global distribution of Buruli ulcer: A systematic review with evidence consensus. Lancet Glob. Health 2019, 7, e912–e922. [Google Scholar] [CrossRef]
  11. Pouillot, R.; Matias, G.; Wondje, C.M.; Portaels, F.; Valin, N.; Ngos, F.; Njikap, A.; Marsollier, L.; Fontanet, A.; Eyangoh, S. Risk factors for Buruli ulcer: A case control study in Cameroon. PLoS Negl. Trop. Dis. 2007, 1, e101. [Google Scholar] [CrossRef] [PubMed]
  12. Ravisse, P. Skin ulcer caused by Mycobacterium ulcerans in Cameroon. I. Clinical, epidemiological and histological study. Bull. Soc. Pathol. Exot. 1977, 70, 109–124. [Google Scholar]
  13. MacCallum, P.; Tolhurst, J.C.; Buckle, G.; Sissons, H.A. A new mycobacterial infection in man. J. Pathol. Bacteriol. 1948, 60, 93–122. [Google Scholar] [CrossRef] [PubMed]
  14. Dodge, O.G.; Lunn, H.F. Buruli ulcer: A mycobacterial skin ulcer in a Uganda child. J. Trop. Med. Hyg. 1962, 65, 139–142. [Google Scholar]
  15. Tabah, E.N.; Nsagha, D.S.; Bissek, A.C.; Njamnshi, A.K.; Bratschi, M.W.; Pluschke, G.; Um Boock, A. Buruli ulcer in Cameroon: The development and impact of the national control programme. PLoS Negl. Trop. Dis. 2016, 10, e0004224. [Google Scholar] [CrossRef]
  16. Akoachere, J.F.; Nsai, F.S.; Ndip, R.N. A community based study on the mode of transmission, prevention and treatment of Buruli ulcers in Southwest Cameroon: Knowledge, attitude and practices. PLoS ONE 2016, 11, e0156463. [Google Scholar] [CrossRef]
  17. Marion, E.; Landier, J.; Boisier, P.; Marsollier, L.; Fontanet, A.; Le Gall, P.; Aubry, J.; Djeunga, N.; Umboock, A.; Eyangoh, S. Geographic expansion of Buruli ulcer disease, Cameroon. Emerg. Infect. Dis. 2011, 17, 551–553. [Google Scholar] [CrossRef]
  18. Noeske, J.; Kuaban, C.; Rondini, S.; Sorlin, P.; Ciaffi, L.; Mbuagbaw, J.; Portaels, F.; Pluschke, G. Buruli ulcer disease in Cameroon rediscovered. Am. J. Trop. Med. Hyg. 2004, 70, 520–526. [Google Scholar] [CrossRef]
  19. Marion, E.; Eyangoh, S.; Yeramian, E.; Doannio, J.; Landier, J.; Aubry, J.; Fontanet, A.; Rogier, C.; Cassisa, V.; Cottin, J.; et al. Seasonal and regional dynamics of M. ulcerans transmission in environmental context: Deciphering the role of water bugs as hosts and vectors. PLoS Negl. Trop. Dis. 2010, 4, e731. [Google Scholar] [CrossRef]
  20. Portaels, F.; Meyers, W.M.; Ablordey, A.; Castro, A.G.; Chemlal, K.; de Rijk, P.; Elsen, P.; Fissette, K.; Fraga, A.G.; Lee, R.; et al. First cultivation and characterization of Mycobacterium ulcerans from the environment. PLoS Negl. Trop. Dis. 2008, 2, e178. [Google Scholar] [CrossRef]
  21. Marsollier, L.; Stinear, T.; Aubry, J.; Saint Andre, J.P.; Robert, R.; Legras, P.; Manceau, A.L.; Audrain, C.; Bourdon, S.; Kouakou, H.; et al. Aquatic plants stimulate the growth of and biofilm formation by Mycobacterium ulcerans in axenic culture and harbor these bacteria in the environment. Appl. Environ. Microbiol. 2004, 70, 1097–1103. [Google Scholar] [CrossRef] [PubMed]
  22. Willson, S.J.; Kaufman, M.G.; Merritt, R.W.; Williamson, H.R.; Malakauskas, D.M.; Benbow, M.E. Fish and amphibians as potential reservoirs of Mycobacterium ulcerans, the causative agent of Buruli ulcer disease. Infect. Ecol. Epidemiol. 2013, 3, 19946. [Google Scholar]
  23. Johnson, P.D.; Azuolas, J.; Lavender, C.J.; Wishart, E.; Stinear, T.P.; Hayman, J.A.; Brown, L.; Jenkin, G.A.; Fyfe, J.A. Mycobacterium ulcerans in mosquitoes captured during outbreak of Buruli ulcer, Southeastern Australia. Emerg. Infect. Dis. 2007, 13, 1653–1660. [Google Scholar] [CrossRef] [PubMed]
  24. Ross, B.C.; Johnson, P.D.; Oppedisano, F.; Marino, L.; Sievers, A.; Stinear, T.; Hayman, J.A.; Veitch, M.G.; Robins-Browne, R.M. Detection of Mycobacterium ulcerans in environmental samples during an outbreak of ulcerative disease. Appl. Environ. Microbiol. 1997, 63, 4135–4138. [Google Scholar]
  25. Tian, R.B.; Niamke, S.; Tissot-Dupont, H.; Drancourt, M. Detection of Mycobacterium ulcerans DNA in the environment, Ivory Coast. PLoS ONE 2016, 11, e0151567. [Google Scholar] [CrossRef]
  26. Garchitorena, A.; Ngonghala, C.N.; Texier, G.; Landier, J.; Eyangoh, S.; Bonds, M.H.; Guegan, J.F.; Roche, B. Environmental transmission of Mycobacterium ulcerans drives dynamics of Buruli ulcer in endemic regions of Cameroon. Sci. Rep. 2015, 5, 18055. [Google Scholar] [CrossRef]
  27. Ohtsuka, M.; Kikuchi, N.; Yamamoto, T.; Suzutani, T.; Nakanaga, K.; Suzuki, K.; Ishii, N. Buruli ulcer caused by Mycobacterium ulcerans subsp shinshuense: A rare case of familial concurrent occurrence and detection of insertion sequence 2404 in Japan. JAMA Dermatol. 2014, 150, 64–67. [Google Scholar] [CrossRef]
  28. Garchitorena, A.; Guegan, J.F.; Leger, L.; Eyangoh, S.; Marsollier, L.; Roche, B. Mycobacterium ulcerans dynamics in aquatic ecosystems are driven by a complex interplay of abiotic and biotic factors. Elife 2015, 4, e07616. [Google Scholar] [CrossRef]
  29. Portaels, F.; Elsen, P.; Guimaraes-Peres, A.; Fonteyne, P.A.; Meyers, W.M. Insects in the transmission of Mycobacterium ulcerans infection. Lancet 1999, 353, 986. [Google Scholar] [CrossRef]
  30. Mosi, L.; Williamson, H.; Wallace, J.R.; Merritt, R.W.; Small, P.L. Persistent association of Mycobacterium ulcerans with West African predaceous insects of the family Belostomatidae. Appl. Environ. Microbiol. 2008, 74, 7036–7042. [Google Scholar] [CrossRef]
  31. Marsollier, L.; Robert, R.; Aubry, J.; Saint Andre, J.P.; Kouakou, H.; Legras, P.; Manceau, A.L.; Mahaza, C.; Carbonnelle, B. Aquatic insects as a vector for Mycobacterium ulcerans. Appl. Environ. Microbiol. 2002, 68, 4623–4628. [Google Scholar] [CrossRef] [PubMed]
  32. Marion, E.; Chauty, A.; Yeramian, E.; Babonneau, J.; Kempf, M.; Marsollier, L. A case of guilt by association: Water bug bite incriminated in M. ulcerans infection. Int. J. Mycobacteriol. 2014, 3, 158–161. [Google Scholar] [CrossRef] [PubMed]
  33. Benbow, M.E.; Williamson, H.; Kimbirauskas, R.; McIntosh, M.D.; Kolar, R.; Quaye, C.; Akpabey, F.; Boakye, D.; Small, P.; Merritt, R.W. Aquatic invertebrates as unlikely vectors of Buruli ulcer disease. Emerg. Infect. Dis. 2008, 14, 1247–1254. [Google Scholar] [CrossRef]
  34. Ebong, S.M.A.; Garcia-Pena, G.E.; Pluot-Sigwalt, D.; Marsollier, L.; Le Gall, P.; Eyangoh, S.; Guegan, J.F. Ecology and feeding habits drive infection of water bugs with Mycobacterium ulcerans. Ecohealth 2017, 14, 329–341. [Google Scholar] [CrossRef] [PubMed]
  35. Jung, S.; Duwal, R.K.; Lee, S. COI barcoding of true bugs (Insecta, Heteroptera). Mol. Ecol. Resour. 2011, 11, 266–270. [Google Scholar] [CrossRef] [PubMed]
  36. Park, D.S.; Foottit, R.; Maw, E.; Hebert, P.D.N. Barcoding bugs: DNA-based identification of the true bugs (Insecta: Hemiptera: Heteroptera). PLoS ONE 2011, 6, e18749. [Google Scholar] [CrossRef]
  37. Ebong, S.M.A.; Petit, E.; Le Gall, P.; Chen, P.P.; Nieser, N.; Guilbert, E.; Njiokou, F.; Marsollier, L.; Guegan, J.F.; Pluot-Sigwalt, D.; et al. Molecular species delimitation and morphology of aquatic and sub-aquatic bugs (Heteroptera) in Cameroon. PLoS ONE 2016, 11, e0154905. [Google Scholar]
  38. Sazama, E.J.; Bosch, M.J.; Shouldis, C.S.; Ouellette, S.P.; Wesner, J.S. Incidence of Wolbachia in aquatic insects. Ecol. Evol. 2017, 7, 1165–1169. [Google Scholar] [CrossRef]
  39. Weinert, L.A.; Araujo-Jnr, E.V.; Ahmed, M.Z.; Welch, J.J. The incidence of bacterial endosymbionts in terrestrial arthropods. Proc. Biol. Sci. 2015, 282, 20150249. [Google Scholar] [CrossRef]
  40. Teixeira, L.; Ferreira, A.; Ashburner, M. The bacterial symbiont Wolbachia induces resistance to RNA viral infections in Drosophila melanogaster. PLoS Biol. 2008, 6, e2. [Google Scholar] [CrossRef] [Green Version]
  41. Zele, F.; Nicot, A.; Berthomieu, A.; Weill, M.; Duron, O.; Rivero, A. Wolbachia increases susceptibility to Plasmodium infection in a natural system. Proc. Biol. Sci. 2014, 281, 20132837. [Google Scholar] [CrossRef] [PubMed]
  42. Graham, R.I.; Grzywacz, D.; Mushobozi, W.L.; Wilson, K. Wolbachia in a major African crop pest increases susceptibility to viral disease rather than protects. Ecol. Lett. 2012, 15, 993–1000. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Kamtchum-Tatuene, J.; Makepeace, B.L.; Benjamin, L.; Baylis, M.; Solomon, T. The potential role of Wolbachia in controlling the transmission of emerging human arboviral infections. Curr. Opin. Infect. Dis. 2017, 30, 108–116. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Jeffries, C.L.; Lawrence, G.G.; Golovko, G.; Kristan, M.; Orsborne, J.; Spence, K.; Hurn, E.; Bandibabone, J.; Tantely, L.M.; Raharimalala, F.N.; et al. Novel Wolbachia strains in Anopheles malaria vectors from sub-Saharan Africa. Wellcome Open Res. 2018, 3, 113. [Google Scholar] [CrossRef]
  45. Fromont, C.; Adair, K.L.; Douglas, A.E. Correlation and causation between the microbiome, Wolbachia and host functional traits in natural populations of drosophilid flies. Mol. Ecol. 2019, 28, 1826–1841. [Google Scholar] [CrossRef]
  46. Cheng, D.; Chen, S.; Huang, Y.; Pierce, N.E.; Riegler, M.; Yang, F.; Zeng, L.; Lu, Y.; Liang, G.; Xu, Y. Symbiotic microbiota may reflect host adaptation by resident to invasive ant species. PLoS Pathog. 2019, 15, e1007942. [Google Scholar] [CrossRef] [Green Version]
  47. Bolz, M.; Bratschi, M.W.; Kerber, S.; Minyem, J.C.; Um Boock, A.; Vogel, M.; Bayi, P.F.; Junghanss, T.; Brites, D.; Harris, S.R.; et al. Locally confined clonal complexes of Mycobacterium ulcerans in two Buruli ulcer endemic regions of Cameroon. PLoS Negl. Trop. Dis. 2015, 9, e0003802. [Google Scholar] [CrossRef] [Green Version]
  48. Damgaard, J.; Buzzetti, F.M.; Mazzucconi, S.A.; Weir, T.A.; Zettel, H. A molecular phylogeny of the pan-tropical pond skater genus Limnogonus Stål 1868 (Hemiptera-Heteroptera: Gerromorpha-Gerridae). Mol. Phylogenet. Evol. 2010, 57, 669–677. [Google Scholar] [CrossRef]
  49. Tchakonté, S.; Ajeagah, G.A.; Tchatcho, N.L.N.; Camara, A.I.; Diomandé, D.; Ngassam, P. Stream’s water quality and description of some aquatic species of Coleoptera and Hemiptera (Insecta) in Littoral Region of Cameroon. Biodivers. J. 2015, 6, 27–40. [Google Scholar]
  50. Mbogho, A.Y.; Sites, R.W. Naucoridae Leach, 1815 (Hemiptera: Heteroptera) of Tanzania. Afr. Invertebr. 2013, 54, 513–542. [Google Scholar] [CrossRef]
  51. Ebong, S.M.; Eyangoh, S.; Marion, E.; Landier, J.; Marsollier, L.; Guegan, J.F.; Legall, P. Survey of water bugs in Bankim, a new Buruli ulcer endemic area in Cameroon. J. Trop. Med. 2012, 2012, 123843. [Google Scholar] [CrossRef] [PubMed]
  52. Chen, P.P.; Nieser, N.; Zettel, H. The Aquatic and Semiaquatic Bugs (Heteroptera—Nepomorpha & Gerromorpha) of Malesia; Brill: Leiden, The Netherlands, 2005. [Google Scholar]
  53. Andersen, N.M.; Weir, T.A. Australian Water Bugs. Their Biology and Identification (Hemiptera-Heteroptera, Gerromorpha & Nepomorpha); Apollo Books; CSIRO Publishing: Clayton, Australia, 2004; Volume 14. [Google Scholar]
  54. Andersen, N.M. The Semiaquatic Bugs (Hemiptera, Gerromorpha). Phylogeny, Adaptations, Biogeography, and Classification; Scandinavian Science Press Ltd.: Ganløse, Denmark, 1982; Volume 3. [Google Scholar]
  55. De Moor, I.J.; Day, J.A.; De Moor, F.C. Insecta II. Hemiptera, Megaloptera, Neuroptera, Trichoptera and Lepidoptera; Water Research Commission: Pretoria, South Africa, 2003; Volume 8. [Google Scholar]
  56. Chen, P.P.; Zettel, H. Five new species of the Halobatinae genus Metrocoris Mayr, 1865 (Insecta: Hemiptera: Gerridae) from continental Asia. Ann. Naturhist. Mus. Wien Serie B 1999, 101, 13–32. [Google Scholar]
  57. Hecher, C.; Zettel, H. Faunistical and morphological notes on Limnogonus subgenus Limnogonoides Poisson 1965 (Heteroptera: Gerridae). Linz. Biol. Beitr. 1996, 28, 325–333. [Google Scholar]
  58. Schuh, R.T.; Slater, J.A. True Bugs of the World (Hemiptera: Heteroptera) Classification and Natural History; Cornell University Press: Ithaca, NY, USA, 1995. [Google Scholar]
  59. Dethier, M. Introduction à la morphologie, la biologie et la classification des Hétéroptères. Bull. Romand Entomol. 1981, 1, 11–16. [Google Scholar]
  60. Poisson, R.A. Contribution à l’étude des Hydrocorises de Madagascar (mission J. Millot et R. Paulian 1949). Mém. Inst. Sci. Madag. 1952, 1, 24–70. [Google Scholar]
  61. Poisson, R.A. Contribution à la faune du Cameroun Hémiptères aquatiques. Faune Colon. Fr. 1929, 3, 135–164. [Google Scholar]
  62. Poisson, R.A. De voyage ML. Chopard en Côte d’Ivoire (1938–39)—Hémiptères aquatiques. Rev. Fr. Entomol. 1941, 8, 77–82. [Google Scholar]
  63. Poisson, R.A. Hydrocorises du Cameroun. Mission, J. Carayon 1947. Rev. Fr. Entomol. 1948, 3, 167–177. [Google Scholar]
  64. Poisson, R.A.; de Voyage, M.P.P. Grasse en Afrique Occidentale Française. Hémiptères aquatiques. Ann. Soc. Entomol. Fr. 1937, 106, 115–132. [Google Scholar]
  65. Poisson, R.A. Quelques Hydrocorises nouveaux de l’Afrique du Sud (mission Suédoise Brink et Rüdebeck). Bull. Soc. Sci. Bretagne 1955, 30, 129–140. [Google Scholar]
  66. Poisson, R.A. Contribution a l’etude des Hydrocorises de Madagascar (Heteroptera). 4e memoire. Mém. Inst. Sci. Madag. 1956, 7, 243–265. [Google Scholar]
  67. Poisson, R.A. Catalogue des Hétéroptères Hydrocorises Africano-Malgaches de la famille des Nepidae (Ltreille) 1802. Bull. IFAN Tome 1965, 27, 229–269. [Google Scholar]
  68. Poisson, R.A. Catalogue des insectes Hétéroptères Gerridae Leach, 1807 Africano-Malgaches. Bull. IFAN Tome 1965, 27, 1466–1503. [Google Scholar]
  69. Poisson, R.A. Catalogue des insectes Hétéroptères Notonectidae Leach, 1815 Africano-Malgaches. Bull. IFAN Tome 1966, 28, 729–768. [Google Scholar]
  70. Ammazzalorso, A.D.; Zolnik, C.P.; Daniels, T.J.; Kolokotronis, S.O. To beat or not to beat a tick: Comparison of DNA extraction methods for ticks (Ixodes scapularis). PeerJ 2015, 3, e1147. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  71. Folmer, O.; Black, M.; Hoeh, W.; Lutz, R.; Vrijenhoek, R. DNA primers for amplification of mitochondrial cytochrome c oxidase subunit I from diverse metazoan invertebrates. Mol. Mar. Biol. Biotechnol. 1994, 3, 294–299. [Google Scholar] [PubMed]
  72. Ablordey, A.; Amissah, D.A.; Aboagye, I.F.; Hatano, B.; Yamazaki, T.; Sata, T.; Ishikawa, K.; Katano, H. Detection of Mycobacterium ulcerans by the loop mediated isothermal amplification method. PLoS Negl. Trop. Dis. 2012, 6, e1590. [Google Scholar] [CrossRef] [PubMed]
  73. Zhou, W.; Rousset, F.; O’Neil, S. Phylogeny and PCR-based classification of Wolbachia strains using wsp gene sequences. Proc. Biol. Sci. 1998, 265, 509–515. [Google Scholar] [CrossRef]
  74. Dong, X.F.; Armstrong, S.D.; Xia, D.; Makepeace, B.L.; Darby, A.C.; Kadowaki, T. Draft genome of the honey bee ectoparasitic mite, Tropilaelaps mercedesae, is shaped by the parasitic life history. Gigascience 2017, 6, gix008. [Google Scholar] [CrossRef] [Green Version]
  75. Wood, D.E.; Salzberg, S.L. Kraken: Ultrafast metagenomic sequence classification using exact alignments. Genome Biol. 2014, 15, R46. [Google Scholar] [CrossRef] [Green Version]
  76. Ounit, R.; Lonardi, S. Higher classification sensitivity of short metagenomic reads with CLARK-S. Bioinformatics 2016, 32, 3823–3825. [Google Scholar] [CrossRef] [Green Version]
  77. Ondov, B.D.; Bergman, N.H.; Phillippy, A.M. Interactive metagenomic visualization in a web browser. BMC Bioinform. 2011, 12, 385. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  78. Truong, D.T.; Franzosa, E.A.; Tickle, T.L.; Scholz, M.; Weingart, G.; Pasolli, E.; Tett, A.; Huttenhower, C.; Segata, N. Metaphlan2 for enhanced metagenomic taxonomic profiling. Nat. Methods 2015, 12, 902–903. [Google Scholar] [CrossRef] [PubMed]
  79. Zerbino, D.R.; Birney, E. Velvet: Algorithms for de novo short read assembly using de Bruijn graphs. Genome Res. 2008, 18, 821–829. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  80. Laetsch, D.R.; Blaxter, M.L. Blobtools: Interrogation of genome assemblies. F1000Research 2017, 6, 1287. [Google Scholar] [CrossRef]
  81. Dierckxsens, N.; Mardulyn, P.; Smits, G. Novoplasty: De novo assembly of organelle genomes from whole genome data. Nucleic Acids Res. 2017, 45, e18. [Google Scholar]
  82. Bernt, M.; Donath, A.; Juhling, F.; Externbrink, F.; Florentz, C.; Fritzsch, G.; Putz, J.; Middendorf, M.; Stadler, P.F. Mitos: Improved de novo metazoan mitochondrial genome annotation. Mol. Phylogenet. Evol. 2013, 69, 313–319. [Google Scholar] [CrossRef]
  83. Lowe, T.M.; Chan, P.P. tRNAscan-SE on-line: Integrating search and context for analysis of transfer RNA genes. Nucleic Acids Res. 2016, 44, W54–W57. [Google Scholar] [CrossRef]
  84. Alikhan, N.F.; Petty, N.K.; Ben Zakour, N.L.; Beatson, S.A. Blast ring image generator (BRIG): Simple prokaryote genome comparisons. BMC Genom. 2011, 12, 402. [Google Scholar] [CrossRef] [Green Version]
  85. Li, H.; Leavengood, J.M., Jr.; Chapman, E.G.; Burkhardt, D.; Song, F.; Jiang, P.; Liu, J.; Zhou, X.; Cai, W. Mitochondrial phylogenomics of Hemiptera reveals adaptive innovations driving the diversification of true bugs. Proc. Biol. Sci. 2017, 284, 20171223. [Google Scholar] [CrossRef]
  86. Hua, J.; Li, M.; Dong, P.; Cui, Y.; Xie, Q.; Bu, W. Phylogenetic analysis of the true water bugs (Insecta: Hemiptera: Heteroptera: Nepomorpha): Evidence from mitochondrial genomes. BMC Evol. Biol. 2009, 9, 134. [Google Scholar] [CrossRef] [Green Version]
  87. Sun, X.; Wang, Y.; Chen, P.P.; Wang, H.; Lixiang, L.; Ye, Z.; Wu, Y.; Li, T.; Bu, W.; Xie, Q. Biased heteroplasmy within the mitogenomic sequences of Gigantometra gigas revealed by Sanger and high-throughput methods. Zool. Syst. 2018, 43, 356–386. [Google Scholar]
  88. Darling, A.C.; Mau, B.; Blattner, F.R.; Perna, N.T. Mauve: Multiple alignment of conserved genomic sequence with rearrangements. Genome Res. 2004, 14, 1394–1403. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  89. Ronquist, F.; Teslenko, M.; van der Mark, P.; Ayres, D.L.; Darling, A.; Hohna, S.; Larget, B.; Liu, L.; Suchard, M.A.; Huelsenbeck, J.P. Mrbayes 3.2: Efficient Bayesian phylogenetic inference and model choice across a large model space. Syst. Biol. 2012, 61, 539–542. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  90. Dong, X.; Chaisiri, K.; Xia, D.; Armstrong, S.D.; Fang, Y.; Donnelly, M.J.; Kadowaki, T.; McGarry, J.W.; Darby, A.C.; Makepeace, B.L. Genomes of trombidid mites reveal novel predicted allergens and laterally transferred genes associated with secondary metabolism. Gigascience 2018, 7, giy127. [Google Scholar] [CrossRef] [Green Version]
  91. Huang, X.; Madan, A. Cap3: A DNA sequence assembly program. Genome Res. 1999, 9, 868–877. [Google Scholar] [CrossRef] [Green Version]
  92. Ratnasingham, S.; Hebert, P.D. Bold: The Barcode of Life data system (http://www.Barcodinglife.Org). Mol. Ecol. Notes 2007, 7, 355–364. [Google Scholar] [CrossRef] [Green Version]
  93. Katoh, K.; Standley, D.M. MAFFT multiple sequence alignment software version 7: Improvements in performance and usability. Mol. Biol. Evol. 2013, 30, 772–780. [Google Scholar] [CrossRef] [Green Version]
  94. Talavera, G.; Castresana, J. Improvement of phylogenies after removing divergent and ambiguously aligned blocks from protein sequence alignments. Syst. Biol. 2007, 56, 564–577. [Google Scholar] [CrossRef] [Green Version]
  95. Nguyen, L.T.; Schmidt, H.A.; von Haeseler, A.; Minh, B.Q. Iq-tree: A fast and effective stochastic algorithm for estimating maximum-likelihood phylogenies. Mol. Biol. Evol. 2015, 32, 268–274. [Google Scholar] [CrossRef]
  96. Bankevich, A.; Nurk, S.; Antipov, D.; Gurevich, A.A.; Dvorkin, M.; Kulikov, A.S.; Lesin, V.M.; Nikolenko, S.I.; Pham, S.; Prjibelski, A.D.; et al. Spades: A new genome assembly algorithm and its applications to single-cell sequencing. J. Comput. Biol. 2012, 19, 455–477. [Google Scholar] [CrossRef] [Green Version]
  97. Seemann, T. Prokka: Rapid prokaryotic genome annotation. Bioinformatics 2014, 30, 2068–2069. [Google Scholar] [CrossRef] [PubMed]
  98. Emms, D.M.; Kelly, S. Orthofinder: Solving fundamental biases in whole genome comparisons dramatically improves orthogroup inference accuracy. Genome Biol. 2015, 16, 157. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  99. Bordenstein, S.R.; Bordenstein, S.R. Eukaryotic association module in phage WO genomes from Wolbachia. Nat. Commun. 2016, 7, 13155. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  100. Gu, H.; Chen, Y.; Liu, X.; Wang, H.; Shen-Tu, J.; Wu, L.; Zeng, L.; Xu, J. The effective migration of Massilia sp. Wf1 by Phanerochaete chrysosporium and its phenanthrene biodegradation in soil. Sci. Total Environ. 2017, 593, 695–703. [Google Scholar] [CrossRef] [PubMed]
  101. Baldo, L.; Werren, J.H. Revisiting Wolbachia supergroup typing based on WSP: Spurious lineages and discordance with MLST. Curr. Microbiol. 2007, 55, 81–87. [Google Scholar] [CrossRef]
  102. LePage, D.P.; Metcalf, J.A.; Bordenstein, S.R.; On, J.; Perlmutter, J.I.; Shropshire, J.D.; Layton, E.M.; Funkhouser-Jones, L.J.; Beckmann, J.F.; Bordenstein, S.R. Prophage WO genes recapitulate and enhance Wolbachia-induced cytoplasmic incompatibility. Nature 2017, 543, 243–247. [Google Scholar] [CrossRef] [Green Version]
  103. Bronstein, O.; Kroh, A.; Haring, E. Mind the gap! The mitochondrial control region and its power as a phylogenetic marker in echinoids. BMC Evol. Biol. 2018, 18, 80. [Google Scholar] [CrossRef]
  104. Song, H.; Buhay, J.E.; Whiting, M.F.; Crandall, K.A. Many species in one: DNA barcoding overestimates the number of species when nuclear mitochondrial pseudogenes are coamplified. Proc. Natl. Acad. Sci. USA 2008, 105, 13486–13491. [Google Scholar] [CrossRef] [Green Version]
  105. Garchitorena, A.; Roche, B.; Kamgang, R.; Ossomba, J.; Babonneau, J.; Landier, J.; Fontanet, A.; Flahault, A.; Eyangoh, S.; Guegan, J.F.; et al. Mycobacterium ulcerans ecological dynamics and its association with freshwater ecosystems and aquatic communities: Results from a 12-month environmental survey in Cameroon. PLoS Negl. Trop. Dis. 2014, 8, e2879. [Google Scholar] [CrossRef]
  106. Hammer, T.J.; Dickerson, J.C.; Fierer, N. Evidence-based recommendations on storing and handling specimens for analyses of insect microbiota. PeerJ 2015, 3, e1190. [Google Scholar] [CrossRef] [Green Version]
  107. Iturbe-Ormaetxe, I.; Walker, T.; O’Neill, S.L. Wolbachia and the biological control of mosquito-borne disease. EMBO Rep. 2011, 12, 508–518. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  108. Crainey, J.L.; Wilson, M.D.; Post, R.J. Phylogenetically distinct Wolbachia gene and pseudogene sequences obtained from the African onchocerciasis vector Simulium squamosum. Int. J. Parasitol. 2010, 40, 569–578. [Google Scholar] [CrossRef] [PubMed]
  109. Woodford, L.; Bianco, G.; Ivanova, Y.; Dale, M.; Elmer, K.; Rae, F.; Larcombe, S.D.; Helm, B.; Ferguson, H.M.; Baldini, F. Vector species-specific association between natural Wolbachia infections and avian malaria in black fly populations. Sci. Rep. 2018, 8, 4188. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  110. Cordaux, R.; Pichon, S.; Hatira, H.B.; Doublet, V.; Greve, P.; Marcade, I.; Braquart-Varnier, C.; Souty-Grosset, C.; Charfi-Cheilchrouha, F.; Bouchon, D. Widespread Wolbachia infection in terrestrial isopods and other crustaceans. Zookeys 2012, 176, 123–131. [Google Scholar] [CrossRef]
  111. Baltanas, A.; Zabal-Aguirre, M.; Pita, M.; Lopez-Fernandez, C. Wolbachia identified in a new crustacean host: An explanation of the prevalence of asexual reproduction in non-marine ostracods? Fund. Appl. Limnol. 2007, 169, 217–221. [Google Scholar] [CrossRef]
  112. Aikawa, T.; Anbutsu, H.; Nikoh, N.; Kikuchi, T.; Shibata, F.; Fukatsu, T. Longicorn beetle that vectors pinewood nematode carries many Wolbachia genes on an autosome. Proc. R. Soc. B Biol. Sci. 2009, 276, 3791–3798. [Google Scholar] [CrossRef] [Green Version]
  113. McNulty, S.N.; Foster, J.M.; Mitreva, M.; Hotopp, J.C.D.; Martin, J.; Fischer, K.; Wu, B.; Davis, P.J.; Kumar, S.; Brattig, N.W.; et al. Endosymbiont DNA in endobacteria-free filarial nematodes indicates ancient horizontal genetic transfer. PLoS ONE 2010, 5, e11029. [Google Scholar] [CrossRef]
  114. Koutsovoulos, G.; Makepeace, B.; Tanya, V.N.; Blaxter, M. Palaeosymbiosis revealed by genomic fossils of Wolbachia in a strongyloidean nematode. PLoS Genet. 2014, 10, e1004397. [Google Scholar] [CrossRef] [Green Version]
  115. Funkhouser-Jones, L.J.; Sehnert, S.R.; Martinez-Rodriguez, P.; Toribio-Fernandez, R.; Pita, M.; Bella, J.L.; Bordenstein, S.R. Wolbachia co-infection in a hybrid zone: Discovery of horizontal gene transfers from two Wolbachia supergroups into an animal genome. PeerJ 2015, 3, e1479. [Google Scholar] [CrossRef] [Green Version]
  116. LePage, D.; Bordenstein, S.R. Wolbachia: Can we save lives with a great pandemic? Trends Parasitol. 2013, 29, 385–393. [Google Scholar] [CrossRef] [Green Version]
  117. Gerth, M.; Gansauge, M.T.; Weigert, A.; Bleidorn, C. Phylogenomic analyses uncover origin and spread of the Wolbachia pandemic. Nat. Commun. 2014, 5, 5117. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  118. Armbruster, P.; Damsky, W.E.; Giordano, R.; Birungi, J.; Munstermann, L.E.; Conn, J.E. Infection of New- and Old-world Aedes albopictus (Diptera: Culicidae) by the intracellular parasite Wolbachia: Implications for host mitochondrial DNA evolution. J. Med. Entomol. 2003, 40, 356–360. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  119. Ellegaard, K.M.; Klasson, L.; Naslund, K.; Bourtzis, K.; Andersson, S.G.E. Comparative genomics of Wolbachia and the bacterial species concept. PLoS Genet. 2013, 9, e1003381. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  120. Casiraghi, M.; Bordenstein, S.R.; Baldo, L.; Lo, N.; Beninati, T.; Wernegreen, J.J.; Werren, J.H.; Bandi, C. Phylogeny of Wolbachia pipientis based on glta, groel and ftsz gene sequences: Clustering of arthropod and nematode symbionts in the F supergroup, and evidence for further diversity in the Wolbachia tree. Microbiol. SGM 2005, 151, 4015–4022. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  121. Hosokawa, T.; Koga, R.; Kikuchi, Y.; Meng, X.Y.; Fukatsu, T. Wolbachia as a bacteriocyte-associated nutritional mutualist. Proc. Natl. Acad. Sci. USA 2010, 107, 769–774. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  122. Tijsse-Klasen, E.; Braks, M.; Scholte, E.J.; Sprong, H. Parasites of vectors—Ixodiphagus hookeri and its Wolbachia symbionts in ticks in the Netherlands. Parasit. Vectors 2011, 4, 228. [Google Scholar] [CrossRef] [Green Version]
  123. Chrostek, E.; Gerth, M. Is Anopheles gambiae a natural host of Wolbachia? MBio 2019, 10, e00784-19. [Google Scholar] [CrossRef] [Green Version]
  124. Cariou, M.; Duret, L.; Charlat, S. The global impact of Wolbachia on mitochondrial diversity and evolution. J. Evol. Biol. 2017, 30, 2204–2210. [Google Scholar] [CrossRef] [Green Version]
  125. Zabal-Aguirre, M.; Arroyo, F.; Garcia-Hurtado, J.; de la Torre, J.; Hewitt, G.M.; Bella, J.L. Wolbachia effects in natural populations of Chorthippus parallelus from the Pyrenean hybrid zone. J. Evol. Biol. 2014, 27, 1136–1148. [Google Scholar] [CrossRef]
  126. Raychoudhury, R.; Baldo, L.; Oliveira, D.C.; Werren, J.H. Modes of acquisition of Wolbachia: Horizontal transfer, hybrid introgression, and codivergence in the Nasonia species complex. Evolution 2009, 63, 165–183. [Google Scholar] [CrossRef]
  127. Miyata, M.; Konagaya, T.; Yukuhiro, K.; Nomura, M.; Kageyama, D. Wolbachia-induced meiotic drive and feminization is associated with an independent occurrence of selective mitochondrial sweep in a butterfly. Biol. Lett. 2017, 13, 20170153. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. Sample collection sites in the Southwest Region of Cameroon.
Figure 1. Sample collection sites in the Southwest Region of Cameroon.
Diversity 11 00225 g001
Figure 2. The distribution of the aquatic bugs collected in the health districts.
Figure 2. The distribution of the aquatic bugs collected in the health districts.
Diversity 11 00225 g002
Figure 3. Aquatic bug families identified in this study and their distribution between six locations (health districts).
Figure 3. Aquatic bug families identified in this study and their distribution between six locations (health districts).
Diversity 11 00225 g003
Figure 4. Representative specimens of the most common genera of aquatic bugs identified in this study: Limnogonus (Gerridae) and Rhagovelia (Veliidae). Note the marked morphlogical diversity within each genus.
Figure 4. Representative specimens of the most common genera of aquatic bugs identified in this study: Limnogonus (Gerridae) and Rhagovelia (Veliidae). Note the marked morphlogical diversity within each genus.
Diversity 11 00225 g004
Figure 5. Bayesian tree of Wolbachia surface protein gene sequences from 12 aquatic bug specimens compared with those representing five major Wolbachia supergroups: wMel (from Drosophila melanogaster), wPip (from Culex pipiens), wOo (from Onchocerca ochengi), wBm (from Brugia malayi), wFol (from Folsomia candida), and wCle (from Cimex lectularius). Letters in parentheses after Wolbachia strain names refer to supergroup designations; suffixes after bug names in square brackets refer to location codes (see Table 3).
Figure 5. Bayesian tree of Wolbachia surface protein gene sequences from 12 aquatic bug specimens compared with those representing five major Wolbachia supergroups: wMel (from Drosophila melanogaster), wPip (from Culex pipiens), wOo (from Onchocerca ochengi), wBm (from Brugia malayi), wFol (from Folsomia candida), and wCle (from Cimex lectularius). Letters in parentheses after Wolbachia strain names refer to supergroup designations; suffixes after bug names in square brackets refer to location codes (see Table 3).
Diversity 11 00225 g005
Figure 6. Metagenomic phylogenetic analysis of Wolbachia sequences from 20 aquatic bug genomes as determined by MetaPhlAn2 [78]. Wolbachia reference genomes are from symbionts of N. vitripennis (wVitB), Culex pipiens molestus (wPip_Mol), Diaphorina citri (wDc), D. melanogaster (wMelPop), Hypolimnas bolina (wBol1), and Culex pipiens pipiens (wPip).
Figure 6. Metagenomic phylogenetic analysis of Wolbachia sequences from 20 aquatic bug genomes as determined by MetaPhlAn2 [78]. Wolbachia reference genomes are from symbionts of N. vitripennis (wVitB), Culex pipiens molestus (wPip_Mol), Diaphorina citri (wDc), D. melanogaster (wMelPop), Hypolimnas bolina (wBol1), and Culex pipiens pipiens (wPip).
Diversity 11 00225 g006
Figure 7. BLAST Ring Image Generator [84] plots for circularized mitogenomes of Gerridae (A,B) and Veliidae (C). Details of reference genomes from National Center for Biotechnology Information (NCBI) are provided in the center of each map.
Figure 7. BLAST Ring Image Generator [84] plots for circularized mitogenomes of Gerridae (A,B) and Veliidae (C). Details of reference genomes from National Center for Biotechnology Information (NCBI) are provided in the center of each map.
Diversity 11 00225 g007
Figure 8. Phylogenetic tree of coi sequences from gerrids (green), veliids (blue), and mesoveliids (dark red). Sequences from the current study are marked with an asterisk (*) and show location codes in square brackets, references from NCBI show accession numbers, and sequences from the prior study of Ebong et al. [37] are represented by short alphanumeric codes. The “Cluster Sequences” feature of BOLD Systems [92] was used to generate the operational taxonomic units (OTUs) and nearest-neighbor distances (in parentheses). Dubious OTUs (labelled in red text) are associated with stop codons in coi sequences, indicated by a caret symbol (^) after the putative species name. Numbers on the nodes of the tree are in the format “SH-like approximate likelihood ratio test support (%)/ultrafast bootstrap support (%)”.
Figure 8. Phylogenetic tree of coi sequences from gerrids (green), veliids (blue), and mesoveliids (dark red). Sequences from the current study are marked with an asterisk (*) and show location codes in square brackets, references from NCBI show accession numbers, and sequences from the prior study of Ebong et al. [37] are represented by short alphanumeric codes. The “Cluster Sequences” feature of BOLD Systems [92] was used to generate the operational taxonomic units (OTUs) and nearest-neighbor distances (in parentheses). Dubious OTUs (labelled in red text) are associated with stop codons in coi sequences, indicated by a caret symbol (^) after the putative species name. Numbers on the nodes of the tree are in the format “SH-like approximate likelihood ratio test support (%)/ultrafast bootstrap support (%)”.
Diversity 11 00225 g008
Figure 9. Bayesian tree of whole mitogenomes from nine Cameroonian aquatic bugs (Gerridae: Limnogonus spp., Metrocoris sp., Trepobates sp.; Veliidae; Rhagovelia spp.) plus references from NCBI (Gerridae: Aquarius paludum, Gigantometra gigas; Veliidae, Perittopus sp.). The tree was constructed using protein-coding sequences only.
Figure 9. Bayesian tree of whole mitogenomes from nine Cameroonian aquatic bugs (Gerridae: Limnogonus spp., Metrocoris sp., Trepobates sp.; Veliidae; Rhagovelia spp.) plus references from NCBI (Gerridae: Aquarius paludum, Gigantometra gigas; Veliidae, Perittopus sp.). The tree was constructed using protein-coding sequences only.
Diversity 11 00225 g009
Table 1. Composition of aquatic bugs based on morphological characteristics and confirmation at the Natural History Museum, London.
Table 1. Composition of aquatic bugs based on morphological characteristics and confirmation at the Natural History Museum, London.
Family (n)GenusPutative SpeciesTotalNo. of Morphotypes
Gerridae (498)LimnogonusLimnogonus (L.) curriei Bergroth 191621
Limnogonus (L.) guttatus Poisson 1948451
Limnogonus (L.) hypoleucus Gerstaecker 187341
Limnogonus (L.) intermedius Poisson 1941571
Limnogonus (L.) poissoni Andersen 1975831
Limnogonus (s. str.) cereinventris Signoret 1862301
Limnogonus (L.) spp.417
Limnogonus (s. str.) sp.11
AquariusAquarius remigis31
Aquarius spp.63
TenagogonusTenagogonus olbovitlotus61
Tenagogonus spp.84
TrepobatesTrepobates spp.8913
MetrobatesMetrobates spp.5510
MetrocorisMetrocoris spp.164
NeogerrisNeogerris severini11
RhagadotarsusRhagadotarsus sp.11
EurymetraEurymetra sp.11
EurymetropsisEurymetropsis sp.501
Veliidae (426)RhagoveliaRhagovelia reitteri Reuters, 18821841
Rhagovelia spp.1959
MicroveliaMicrovelia spp.476
Notonectidae (54)WalambianisopsWalambianisops spp.205
AnisopsAnisops spp.113
ParanisopsParanisops sp.11
EnitharesEnithares sp.52
NotonectaNotonecta sp.21
NychiaNychia sp.154
Naucoridae (52)AneurocorisAneurocoris sp.251
IllyocorisIllyocoris sp.62
LaccocorisLaccocoris sp.51
MaccrocolisMaccrocolis laticollis21
NaucorisNaucoris obscuratus21
Naucoris sp.93
NeomaccrocorisNeomaccrocoris parviceps ocellatus21
Neomaccrocoris parviceps parviceps11
Belostomatidae (36)BelostomaBelostoma sp.11
DiplonychusDiplonychus sp.346
LethocerusLethocerus sp.11
Mesoveliidae (18)MesoveliaMesovelia sp.183
Hydrometridae (11)HydrometraHydrometra sp.111
Nepidae (7)RanatraRanatra sp.71
Total1102110
Table 2. Detection and quantification of rpoB gene fragments in aquatic bugs from the study sites.
Table 2. Detection and quantification of rpoB gene fragments in aquatic bugs from the study sites.
Family (n).GenusPutative SpeciesTotal AnalyzedrpoB Positive (%)Mean rpoB Copy no.
Gerridae (405)LimnogonusLimnogonus (L.) guttatus Poisson 1948397 (17.9)63
Limnogonus (L.) poissoni Andersen 19757615 (19.7)143
Limnogonus (L.) intermedius Poisson 1941507 (14)82
Limnogonus (L.) curriei Bergroth 191620 (0.0)0
Limnogonus (s. str.) cereinventris Signoret 1862247 (29.2)80
Limnogonus (L.) hypoleucus Gerstaecker 187321 (50)58
Limnogonus spp.340 (0.0)0
AquariusAquarius spp.50 (0.0)0
TenagogonusTenagogonus spp.82 (25.0)13
MetrobatesMetrobates spp.425 (11.9)47
TrepobatesTrepobates spp.704 (5.7)50
EurymetropsisEurymetropsis sp.464 (8.7)112
MetrocorisMetrocoris spp.71 (14.3)91
Veliidae (366)RhagoveliaRhagovelia reitteri Reuters 18821713 (1.8)44
Rhagovelia spp.16411 (6.7)111
MicroveliaMicrovelia spp.312 (6.5)33
Total77169 (8.9)
Table 3. Distribution of sequences classified as originating from Mycobacterium spp. and Wolbachia in 20 aquatic bug genomes as determined by Kraken2.
Table 3. Distribution of sequences classified as originating from Mycobacterium spp. and Wolbachia in 20 aquatic bug genomes as determined by Kraken2.
SampleHealth DistrictLocation CodeSpeciesTotal ReadsAll Bacterial Read Pairs (% Total)Mycobacterium Read Pairs (% Bacterial)M. Ulcerans Read Pairs (% Bacterial)Wolbachia Read Pairs (% Bacterial)Wsp PCR
1BueaB05Tenagogonus sp. 314,767,85258,114 (0.39%)82 (0.14%)0 (0.00%)125 (0.22%)-
2BueaB01Limnogonus guttatus13,155,71175,982 (0.58%)62 (0.08%)1 (0.00%)17,278 (22.74%)+
3BueaB02Limnogonus cereiventris27,361,770161,821 (0.59%)53 (0.03%)1 (0.00%)30,423 (18.80%)+
4BueaB04Limnogonus cereiventris12,617,27557,802 (0.46%)18 (0.03%)0 (0.00%)13,120 (22.70%)+
5BueaB02Microvelia sp. 115,199,889156,452 (1.03%)35 (0.02%)0 (0.00%)6,017 (3.85%)+
6KumbaK02Metrobates sp. 711,913,06961,477 (0.52%)73 (0.12%)1 (0.00%)29 (0.05%)-
7KumbaK02Metrobates sp. 611,847,97062,737 (0.53%)57 (0.09%)0 (0.00%)70 (0.11%)-
8KumbaK02Metrobates sp. 812,711,11066,881 (0.53%)68 (0.10%)0 (0.00%)32 (0.05%)-
9KumbaK01Rhagovelia reitteri12,615,049168,912 (1.34%)148 (0.09%)0 (0.00%)106,556 (63.08%)+
10LimbeL02Limnogonus hypoleucus13,728,59869,491 (0.51%)172 (0.25%)89 (0.13%)7470 (10.75%)+
11MbongeMb03Limnogonus poissoni13,897,92787,297 (0.63%)67 (0.08%)0 (0.00%)24,629 (28.21%)+
12MbongeMb03Limnogonus poissoni13,479,48064,224 (0.48%)75 (0.12%)0 (0.00%)13,345 (20.78%)+
13MbongeMb02Limnogonus intermedius15,587,895113,606 (0.73%)61 (0.05%)0 (0.00%)9359 (8.24%)+
14MbongeMb04Limnogonus poissoni14,300,93782,433 (0.58%)57 (0.07%)0 (0.00%)13,412 (16.27%)+
15MbongeMb06Rhagovelia sp. 214,197,736129,978 (0.92%)60 (0.05%)0 (0.00%)66,515 (51.17%)+
16MuyukaMu02Metrocoris sp. 114,136,65372,422 (0.51%)3598 (4.97%)*0 (0.00%)142 (0.20%)-
17MuyukaMu05Eurymetropsis sp.14,455,52662,105 (0.43%)90 (0.14%)0 (0.00%)48 (0.08%)-
18MuyukaMu01Trepobates sp. 513,054,26857,197 (0.44%)37 (0.06%)0 (0.00%)5326 (9.31%)-
19MuyukaMu05Rhagovelia sp. 115,150,283138,234 (0.91%)59 (0.04%)0 (0.00%)5186 (3.75%)-
20MuyukaMu05Rhagovelia sp. 114,978,846131,153 (0.88%)100 (0.08%)0 (0.00%)65,354 (49.83%)+
* 98% of these reads were classified as deriving from Mycobacterium sp. WY10.

Share and Cite

MDPI and ACS Style

Esemu, S.N.; Dong, X.; Kfusi, A.J.; Hartley, C.S.; Ndip, R.N.; Ndip, L.M.; Darby, A.C.; Post, R.J.; Makepeace, B.L. Aquatic Hemiptera in Southwest Cameroon: Biodiversity of Potential Reservoirs of Mycobacterium ulcerans and Multiple Wolbachia Sequence Types Revealed by Metagenomics. Diversity 2019, 11, 225. https://doi.org/10.3390/d11120225

AMA Style

Esemu SN, Dong X, Kfusi AJ, Hartley CS, Ndip RN, Ndip LM, Darby AC, Post RJ, Makepeace BL. Aquatic Hemiptera in Southwest Cameroon: Biodiversity of Potential Reservoirs of Mycobacterium ulcerans and Multiple Wolbachia Sequence Types Revealed by Metagenomics. Diversity. 2019; 11(12):225. https://doi.org/10.3390/d11120225

Chicago/Turabian Style

Esemu, Seraphine N., Xiaofeng Dong, Achah J. Kfusi, Catherine S. Hartley, Roland N. Ndip, Lucy M. Ndip, Alistair C. Darby, Rory J. Post, and Benjamin L. Makepeace. 2019. "Aquatic Hemiptera in Southwest Cameroon: Biodiversity of Potential Reservoirs of Mycobacterium ulcerans and Multiple Wolbachia Sequence Types Revealed by Metagenomics" Diversity 11, no. 12: 225. https://doi.org/10.3390/d11120225

APA Style

Esemu, S. N., Dong, X., Kfusi, A. J., Hartley, C. S., Ndip, R. N., Ndip, L. M., Darby, A. C., Post, R. J., & Makepeace, B. L. (2019). Aquatic Hemiptera in Southwest Cameroon: Biodiversity of Potential Reservoirs of Mycobacterium ulcerans and Multiple Wolbachia Sequence Types Revealed by Metagenomics. Diversity, 11(12), 225. https://doi.org/10.3390/d11120225

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop