Next Article in Journal
Reviewing the WHO Tube Bioassay Methodology: Accurate Method Reporting and Numbers of Mosquitoes Are Key to Producing Robust Results
Previous Article in Journal
Lucilia sericata (Diptera: Calliphoridae) as Agent of Myiasis in a Goose in Italy and a Review of Myiasis by This Species in Birds
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Population Genetic Structure of the Bean Leaf Beetle Ootheca mutabilis (Coleoptera: Chrysomelidae) in Uganda

1
National Agricultural Research Organization (NARO), National Crops Resources Research Institute (NaCRRI), Kampala P.O. Box 7084, Uganda
2
College of Veterinary Medicine, Animal Resources and Biosecurity (CoVAB), Makerere University, Kampala P.O. Box 7062, Uganda
3
College of Natural Sciences (CoNAS), Makerere University, Kampala P.O. Box 7062, Uganda
4
Commonwealth Scientific and Industrial Research Organisation, Canberra, ACT 2601, Australia
5
International Institute of Tropical Agriculture, PMB 5320, Oyo Rd., Ibadan 20001, Nigeria
6
Natural Resources Institute, University of Greenwich, Central Avenue, Chatham Maritime, Kent ME4 4TB, UK
*
Author to whom correspondence should be addressed.
Submission received: 5 March 2022 / Revised: 26 May 2022 / Accepted: 4 June 2022 / Published: 14 June 2022

Abstract

:

Simple Summary

The bean leaf beetle (Ootheca mutabilis) has lately emerged as a major bean pest in Uganda, causing devastating crop losses. Despite its importance, little is known about its population genetic structure. We developed microsatellite DNA markers and combined them with partial mito-chondrial cytochrome oxidase subunit I gene sequences as a marker to examine the spatial pop-ulation genetic structure of 86 O. mutabilis samples from 16 populations. We developed a set of five highly polymorphic microsatellite DNA markers. From both types of markers, nearly all the genetic variation occurred within populations and there was no evidence of genetic differentiation in both markers. There was no isolation by distance between geographical and genetic distances for both markers except in one of the agro-ecological zones for mtCOI data. This information will assist in the design of O. mutabilis control strategies.

Abstract

Bean leaf beetle (BLB) (Ootheca mutabilis) has emerged as an important bean pest in Uganda, leading to devastating crop losses. There is limited information on the population genetic structure of BLB despite its importance. In this study, novel microsatellite DNA markers and the partial mitochondrial cytochrome oxidase subunit I (mtCOI) gene sequences were used to analyze the spatial population genetic structure, genetic differentiation and haplotype diversity of 86 O. mutabilis samples from 16 (districts) populations. We identified 19,356 simple sequence repeats (SSRs) (mono, di-, tri-, tetra-, penta-, and hexa-nucleotides) of which 81 di, tri and tetra-nucleotides were selected for primer synthesis. Five highly polymorphic SSR markers (4–21 alleles, heterozygosity 0.59–0.84, polymorphic information content (PIC) 50.13–83.14%) were used for this study. Analyses of the 16 O. mutabilis populations with these five novel SSRs found nearly all the genetic variation occurring within populations and there was no evidence of genetic differentiation detected for both types of markers. Also, there was no evidence of isolation by distance between geographical and genetic distances for SSR data and mtCOI data except in one agro-ecological zone for mtCOI data. Bayesian clustering identified a signature of admixture that suggests genetic contributions from two hypothetical ancestral genetic lineages for both types of markers, and the minimum-spanning haplotype network showed low differentiation in minor haplotypes from the most common haplotype with the most common haplotype occurring in all the 16 districts. A lack of genetic differentiation indicates unrestricted migrations between populations. This information will contribute to the design of BLB control strategies.

1. Introduction

In East and Southern Africa, the bean leaf beetle (Ootheca mutabilis (Coleoptera: Chrysomelidae)) is an economically important pest of the common bean (Phaseolus vulgaris L.) [1]. In the recent past, Ootheca species have been classified as field pests of common bean (P. vulgaris L.) in Uganda’s Northern and Eastern regions [2,3]. The larvae consume and destroy root tissue, while the adults consume and skeletonize leaves, flowers, and immature pods [4]. During heavy infestations, complete crop losses may occur [1,5]. Bean leaf beetles have also been linked to virus transmission in cowpeas [1,6,7]. Ootheca mutabilis is the most common Ootheca species (70.3%) in all agro-ecological zones of Uganda [8], despite being reported to be primarily found in lowlands [1,4,8]. However, in Uganda, its abundance and damage to common bean vary depending on location and season [8]. Furthermore, O. mutabilis comes in a variety of colours and is frequently confused with other Ootheca species such as O. proteus, O. bennigseni, and O. orientalis [9]. Some previously described Ootheca species have been re-described and re-assigned to other species, resulting in revisions to the morphological identification records [9]. Presently, the only reliable morphological identification approach is dissection and analysis of male genitalia [9] which, excludes species-level identification of females and other growth stages.
Microsatellite genetic markers, commonly known as simple sequence repeat (SSR) markers, are comprised of tandemly repeated nucleotide motifs of 1–6 base pairs (bp). They have been frequently employed in evolutionary genetic studies [10,11] to infer life history, including mating behavior [10,12,13], mating frequencies [14], gene flow patterns [12], dispersal [15,16], range shifting [15], and host-shifting [17]. The advantages of employing microsatellites in population genetics related research are that they are relatively common, co-dominant, ubiquitous, and substantially polymorphic [11]. However, because of their interaction with mobile elements [18,19], it can be difficult to segregate them into certain insect orders [18,20,21]. Microsatellites have been used in a variety of insects, including the European stag beetle [22] and the red flour beetle [23].
High throughput sequencing (HTS) enables cost-effective and rapid identification of a large number of SSR loci even when only a fraction of the sequencing run is used [24]. The data generated from HTS can be applied in de novo whole-genome sequencing of complex genomes. While exploring this approach, small to large coverage can be performed depending on the researcher’s availability of resources [24].
Recent research on O. mutabilis has revealed information on its distribution [4,8], morphological appearance [9], abundance [4,8], and yield losses [8] which, are some of the important information in designing its control strategies. However, there is no available information on its population genetic structure despite being one of the most important considerations in designing long-lasting area-wide pest management strategies. This limits the development of effective area-wide management techniques for these damaging agricultural pests. Understanding BLB population structure, namely if discrete and genetically unique subpopulations exist, is thus necessary for the development of effective control approaches for these destructive pests.
This study aimed to determine the spatial population genetic structure of O. mutabilis in Uganda using microsatellite DNA markers and partial mitochondrial DNA cytochrome c oxidase subunit I (mtCOI). Since there were no existing genetic resources specific to O. mutabilis that could be utilized to determine its population genetic structure, we developed SSR markers which we deployed as a genetic tool. Given the desirable attributes of microsatellites, such as codominance, high polymorphism level, and data reproducibility can be instrumental in revealing small-scale resolution of the demographic events [25,26]. Mitochondrial markers are usually helpful in the analysis of earlier phylogeographic events and large-scale patterns of genetic diversity [27,28]. Combining the two markers that experience different modes of inheritance and degrees of polymorphism provides more information that is unattainable through studies where only one of the markers is used. Molecular diagnostic approaches are required to help delineate BLB species status, with the use of partial mitochondrial DNA cytochrome c oxidase subunit I (mtCOI) sequences being a preferred method globally and successfully employed in our Ugandan laboratory for other agricultural pests (for example, [29,30,31]. Additionally, the mtCOI partial gene was utilized to barcode the morphologically indistinguishable BLBs prior to population genetic structure analysis, given the unreliability of morphological species identification based on colour patterns alone [9] purposely to avoid misidentification.

2. Materials and Methods

2.1. Sample Collection

During the 2016–2017 and 2018 rain seasons, we collected BLBs from 17 districts (Figure 1, Table 1 and Table 2) in five bean production agro-ecological zones (Table 1). Adult BLBs were sampled in farmers’ common bean fields. We selected at least one district from each agro-ecological zone, from which two sub-counties were selected. We picked about 10 insects from ten gardens in each sub-county. Each beetle sample was immediately placed individually in a screw-capped 2 mL vial containing 95% ethanol, and the ethanol was replaced at least twice at day intervals to avoid DNA degradation caused by adult beetle secretions. Before further analysis, samples were transported to the National Crops Resources Research Institute (NaCRRI), Namulonge in Wakiso district, and kept in boxes at room temperature. During sample collection, a GPS coordinate was taken for each garden from which samples were picked, and each GPS point was regarded as one sample since only one sample (BLB) was considered for every GPS point where BLBs were recovered. From every farmer’s field where BLBs were found, at least one BLB was picked. In one of the districts (Lwengo) selected for sample collection in the Central Wooden Savannah agro-ecological zone, BLBs were not recovered at the time of collection (Table 2). This, therefore, reduced the number of districts as well as the number of samples in this agro-ecological zone.

2.2. Identification of Bean Leaf Beetles

Representatives of BLB samples were identified morphologically at the Universitätkoblenz-Landau Institut für Integrierte Naturwissenschaften Abteilung Biologie Universitätsstraße 156,070 Koblenz, Germany, purposely to establish the identity of Ootheca species which information would help us in mass collecting of the true Ootheca species for laboratory analysis. Before the collection of the samples for laboratory analysis, various leaf beetles collected from the bean plants were sorted to pick representatives of each morphotype which were considered for morphological analysis. As described by [9], the BLB samples studied for the population genetic structure were selected for DNA analysis based on colour patterns of the elytra, head, thorax, abdomen, and legs. In this regard, we selected 99 samples based on their appearance as follows: M1 (O. mutabilis with elytra upper half black and lower half brownish) (21 beetles), M2 (O. mutabilis with black elytra) (39 beetles), and M3 (O. mutabilis with brownish elytra) (39 beetles) (Figure 2). This total included the samples found to be O. proteus before DNA barcoding as they could not be at all distinguished from O. mutabilis. All the BLBs with different colour appearances as explained above and shown in Figure 2, as well as reported by [9], were included in the analysis. This colour distinction was made to examine potential genetic variations between O. mutabilis morphotypes. At least, every district from which BLB samples were recovered, was considered as well as some of the samples collected. The number of beetles per district considered for analysis was different because after DNA barcoding, non-O. mutabilis samples were removed (Hoima and Nakasongola) (Table 2). Also, because we wanted to consider different morphotypes, in some districts, samples were increased so that different morphotypes would be included since some districts had one morphotype and others had two or all three morphotypes. Also, during sample selection in the laboratory, one sample was considered per GPS point.

2.3. DNA Isolation and Quantification

For DNA isolation, we used the Qiagen DNeasy Blood & Tissue Kit (Hilden, Germany) in accordance with the manufacturer’s instructions. Following isolation, each DNA sample was quantified using an Agilent Technologies NanoDrop 2000 spectrophotometer (Waldbronn, Germany) and the quality was confirmed using 1% agarose gel electrophoresis.
To confirm their identity, all samples selected based on colour were barcoded using the mtCOI partial gene primers. The PCR primers used were BLB-LCO: 5′-GGTCAACAAATCATAAAGATATTGG-3′ and BLB-HCO: 5′-TAAACTTCAGGGTGACCAAAAAATCA-3′, which amplify a 710 bp fragment (M. Otim unpublished). Each reaction was carried out in a 25 µL volume comprising 1 µL of template DNA, 1 µL of 10 pmol/µL primer, 2.5 µL of 10X DreamTaq Green buffer, 0.5 µL of dNTP (10 mM), 0.25 µL (1.25 units) of Taq DNA polymerase (5 U/µL), 2.5 µL 5% tween20, and 16.25 µL of nuclease-free water. The PCR conditions were as follows: 1 cycle of denaturation at 95 °C for 2 min, 35 cycles of 20 s at 95 °C, 30 s at 52 °C annealing temperature, 1 min at 72 °C, and a final extension cycle at 72 °C for 10 min, after which reactions were held at 4 °C. All samples were tested for amplification success on 1.3% w/v agarose in TAE buffer and stained with ethidium bromide as stated above. Sequences generated from mtCOI PCR products were processed with Pregap4 and Gap4 [32] and compared to NCBI sequences of O. mutabilis (KY574530.1, KY574526.1, KY574527.1) and O. proteus (KY574525.1, KY574524.1, KY574523.1, KY574522.1).

2.4. Genome Sequencing, Quality Check and Raw Read Assembly

GENEWIZ (www.genewiz.com was contracted to carry out high throughput sequencing (HTS). DNA was isolated from three ethanol-preserved BLBs denoted by the letters M1, M2, and M3 (Figure 2). The Illumina HiSeq 2500 system was used to construct and sequence whole-genome DNA libraries, with an insert size of 300 to 400 bp and 2 × 150 paired-end reads. Fast QC v0.11.7 [33] was used to check the quality of raw reads before processing them for de novo assembly. Following QC, raw reads were processed, with adaptor sequences trimmed, duplicate sequences eliminated, and the sequences assembled de novo. All raw read processing and de novo assembly were performed in Geneious v10.2 [34] using the default settings (i.e., allow gaps—maximum gaps per read 20%, ignore words repeated more than 1000 times, do not merge variants with coverage over approximately 6, merge homopolymer variants).

2.5. Microsatellite Prediction, Primer Design and Blast Search of Microsatellite Sequences in GenBank

WebSat [35] was used to identify microsatellites in the assembled contigs, and Primer3 [36] was used to design primers. The primer design parameters included a primer size range of 18 nucleotides at the lowest, 22 at the optimum, and 27 at the maximum. Primer Tm was 57 °C at minimum, 60 °C at optimum, and 68 °C at maximum. The primer GC% ranged from 40 to 80. The product size range was 100–400 bp. After designing the primers, 81 desalted primer pairs were ordered from Macrogen Europe (dna.macrogen.com). Primers were designed for all di-nucleotide, tri-nucleotide, and tetra-nucleotide microsatellite loci identified. Microsatellite sequences were compared with microsatellite sequences in the NCBI GenBank using BLASTX and BLASTN [37] to find out if they shared similarities with sequences from other insects and putative retrotransposable elements (e.g., [18]).

2.6. Microsatellite DNA Marker PCR Optimization, Polymorphism Testing, Primer Labelling and Fragment Analysis

Each microsatellite primer pair was optimized for amplification prior to being evaluated for polymorphism on eight O. mutabilis samples (at least one sample from every population). Each locus was amplified in a 12.5 µL PCR reaction containing 0.5 µL of 50 ng DNA template, 0.5 µL of 10 pmol/L primer, 1.25 µL of 10X DreamTaq green buffer, 0.25 µL of DreamTaq dNTP (10 mM), 0.125 µL of DreamTaq DNA polymerase (5 U/µL), 1.25 µL% Tween20, and 8.125 µL of nuclease-free water. The PCR conditions were as follows: initial denaturation (4 min 94 °C), 35 amplification cycles (20 s 94 °C, 30 s annealing temperature, 45 s 72 °C), and a final extension (10 min 72 °C). Then, the reactions were held at 4 °C. The amplicon was run on a 3% w/v agarose in 1X TAE buffer at first, then on 6% acrylamide gels for 6 h at 120 V with a 100 bp DNA ladder in 1X TAE buffer. Polyacrylamide gels were stained for at least 20 min in an ethidium bromide (0.5 g/mL) solution before being de-stained in distilled water. The U-genius gel documentation system was used to visualize fragment sizes (www.syngene.com. Each locus was tested at least twice for reproducibility. Following primer optimization, five primers were selected based on polymorphism on polyacrylamide gels. These primers were ordered for synthesis from Macrogen Europe (dna.macrogen-europe.com). Each forward primer was labelled with 6-FAM or HEX fluorescent dyes (dna.macrogen-europe.com) (Table 3). The PCR reactions were carried out in single reactions, and then the 6-FAM and HEX PCR products with similar amplicon product size ranges were pooled and run as multiplex. The Applied Biosystems 3730XL DNA Analyzer was used to analyze fragments (outsourced to Macrogen Europe).

2.7. Genotyping and Data Scoring

Genotyping was performed using GeneMarker v2.6.3 [38], and alleles were scored based on their size. MICROCHECKER [39] was used to assess the accuracy of allele scoring. MolKin v3.0 [40] was used to determine polymorphic information content (PIC), observed and expected heterozygosities (Ho, He).

2.8. Population Genetic Structure and Differentiation Analysis

During data analysis, different districts were treated as different populations leading to a total of 16 populations as detailed (Table 2). The elimination of morphologically indistinguishable samples from the analysis reduced the originally selected samples from 99 to 86.
DnaSP v6 was used to compute mtCOI partial gene diversity components such as the number of haplotypes, haplotype diversity, nucleotide diversity, and polymorphic site estimates [41].
We investigated population genetic structure by AMOVA that partitions total variance into covariance components. It then verifies the hierarchical or non-hierarchical variation distribution (i.e., among populations (fixation index (FST), among populations within groups (FSC) and among groups (FCT)) [42]. Analysis of molecular variance (AMOVA) for both sets of data (each data set analysed independently) was performed in Arlequin version 3.5 using 1000 permutations, all at a 0.05 significant level. Analysis of molecular variance (AMOVA) was performed using both non-hierarchical (i.e., all populations in one group) and hierarchical designs (i.e., populations subdivided into five agro-ecological zones as detailed in Table 2). The spatial analysis of molecular variance for both types of markers was calculated using SAMOVA 2.0 software [43] to obtain inferences about the population groups. The SAMOVA for both types of markers were analyzed using 1000 simulated annealing processes by varying K (number of groups) from 2 to 10. The best K value was selected according to when FCT reached a plateau. The best K grouping was used to calculate AMOVA in Arlequin to estimate genetic differentiation.
Geographic structures of O. mutabilis populations for both types of markers were investigated using the Bayesian approach implemented in STRUCTURE v2.3.4 [44]. STRUCTURE uses a coalescent genetic technique to group similar multilocus genotypes into inferred ancestral genetic clusters (K), regardless of an individual’s geographical origin. Using the admixture model, we conducted ten independent runs for each value of K ranging from 1 to 5. Each run consisted of a burn-in of 50,000 steps followed by 100,000 Markov-Chain Monte Carlo (MCMC) repetitions. For each potential value of K, ten replicates were used. The mtCOI sequences were first processed so that only haplotypes are used in the inference as explained in the STRUCTURE user manual [44]. The LOCPRIOR command was used to perform the STRUCTURE runs, and the genotypes defined were based on the geographic location of the O. mutabilis samples. The true value of K was estimated using the program STRUCTURE HARVESTER [45], as described in [46]. CLUMPAK [47] was used to visualize the structural results.
The mtCOI haplotype network was inferred using the minimum spanning network approach described in POPART [48] using 2000 bootstraps based on sequence alignment exported from DnaSP v6 as a haplotype nexus file.

2.9. Isolation by Distance (IBD) Analysis

Evidence of IBD for both types of markers was tested using the Mantel test [49] with 9999 permutations in GenAIEx v6.503 [50]. Isolation by distance was calculated for different agro-ecological zones as detailed in Table 2.

3. Results

3.1. Quality Check of NGS Sequences, De Novo Assembly, SSR Prediction and Primer Design

A total of 272,853,156 raw reads were generated, with equal forward and reverse reads. The length of all sequences was 151 bp, with a GC content of 34%. QC results included: basic statistics passed; per base sequence quality passed; per tile sequence quality passed; per sequence quality scores passed; per base sequence content passed; and, sequence length distribution passed. There was no over-representation of the sequences. The assembly of all sequences in the selected part of the reads resulted in a total of 282,696 contigs. From the assembled contigs, a total of 19,356 SSR were identified, including (i) mononucleotides (14,629; 75.6%), (ii) di-nucleotides (2780; 14.4%), (iii) tri-nucleotides (1288; 6.7%), (iv) tetra-nucleotides (352; 1.8%), (v) penta-nucleotides (258; 1.3%), and (vi) hexa-nucleotides (49; 0.3%).

3.2. Microsatellite PCR Optimization and Polymorphism Testing

Sixty-five of the 81 loci analyzed showed multiple fragments on agarose gels and were therefore removed from the analysis. On agarose gels, a total of 16 loci were observed as a single band, and they also appeared polymorphic on 6% acrylamide gels. Six of the 16 polymorphic loci were eliminated due to inadequate PCR amplification (i.e., fuzzy bands in some samples, and failed PCR amplification in others). Two of the ten loci had primers that overlapped the (GT) SSR units and were thus eliminated. We eliminated three of the remaining eight loci due to low repeatability. As a result, five loci were chosen for genotyping (representing a 6% success rate) and labelled with fluorescent dyes (Table 3). The number of SSR alleles ranged from 4 to 21, with an average of 11.6. (Table 3).

3.3. Fragment Analysis and Allele Scoring

A locus was judged to have low polymorphism if the PIC was less than 25% and high polymorphism if the PIC was greater than 50% [51]. Locus BLB2_om66 had the lowest PIC of 50.1% while BLB2_om33 had the highest PIC of 83.1% (Table 4). The five loci employed in this investigation were, therefore, all were considered to be of high polymorphism, with an average PIC of 69.1%.
At the population level, the average PIC based on all loci across all populations was 53.59% (Table 5), while the population average heterozygosity for all loci was 0.73, with BLB2_om33 having the highest heterozygosity of 0.84 and BLB2_om66 having the lowest heterozygosity of 0.59. (Table 4). Observed heterozygosity was higher than expected heterozygosity in all populations, ranging from 0.75–0.84, with an average observed heterozygosity of 0.80 (Table 5). The loci BLB2_om17 and BLB2_om32 had an excess of homozygosity, which could indicate the potential presence of null alleles or allele drop-out.

3.4. Population Genetic Structure, Differentiation and Gene Flow

The 86 O. mutabilis mtCOI partial gene sequences were analyzed, and 21 segregating sites (S) with an average of 0.827 nucleotide differences were found (Kt). We found 20 haplotypes with an estimated haplotype diversity (H, also known as gene diversity, a measure of the probability that two random alleles are different [52]) of 0.51 and a low nucleotide diversity Pi (π, i.e., the average number of nucleotide differences per site in pairwise DNA sequence comparison [52]) of 0.00127. We found moderate haplotype diversity among the 20 haplotypes (GenBank: MW278873-MW278892) (Figure 3).
The minimum spanning haplotype network (Figure 3a) revealed one major haplotype (haplotype 1) that was detected in all five agro-ecologies (A, B, C, D, and E) and 69.77% of the individuals (n = 60). The second most prevalent haplotype was haplotype 14, which was present in 5.81% (n = 5). Three haplotypes (6, 7 and 18) had 2.30% with two individuals each. The remaining 15 haplotypes each had n = 1. (i.e., 1.16% each) (Figure 3a). Haplotypes 6, 7 and 18 with the third-highest frequency (2.30% each) varied in agro-ecologies although most of the individuals (n = 4) belonged to agro-ecological zone A (Figure 3a).
In the 16 populations, the number of segregating sites ranged from 0 to 6 with Lira, Amuru, and Bulisa having the highest number of haplotypes, haplotype diversity ranged from 0 to 0.9 among the districts with Lira district having the highest and nucleotide diversity was low in all the populations (districts) (Table 6).
Another minimum spanning haplotype network (Figure 3b) for the samples from different districts was inferred which resulted in all the individuals from all different districts sharing haplotype 1. In both haplotype networks, all other haplotypes originate from haplotype 1 (Figure 3b) indicating low differentiation from the most common haplotype.
The neighbour-joining phylogenetic tree constructed (with one sample (representative) from each haplotype) clearly clustered together all the O. mutabilis haplotypes with the O. mutabilis reference sequences (KY574530.1, KY574526.1, KY574527.1) with high node support (100%), but with distinct intra-species differences (Figure 4). The O. proteus samples that were morphologically indistinguishable from O. mutabilis were clearly clustered together, with O. proteus reference sequences (KY574525.1, KY574524.1, KY574523.1, KY574522.1) likewise with intra-species differences (Figure 4).
Hierarchical analysis of molecular variance (AMOVA) for both types of markers resulted in all the genetic variation (100%) occurring within the populations. Non-hierarchical AMOVA for microsatellites had the highest genetic variation occurring within the populations and the lower genetic variation occurred among the populations (Table 7(a1)) except for non-hierarchical AMOVA for mtCOI where all the genetic variation occurred within the populations (Table 7(b1)). For SSRs, the fixation indices among groups, among populations within groups, and within populations as calculated with hierarchical AMOVA were FCT = 0.02073 (p = 0.05), FSC = −0.01523 (p = 0.05), and FST = 0.00582 (p = 0.05) (Table 7(a2)) indicating that there is low differentiation and relatively high gene flow. The same result trend was also detected in the mtCOI marker for hierarchical AMOVA (Table 7(b2)).
From the analysis by SAMOVA, for both microsatellites and mtCOI data, the FCT values reached a plateau at K = 2, FCT = 0.4, p = 0.058 and 0.094, p = 0.05474 for SSRs and mtCOI respectively (Table 8(a1) and 8(a2) respectively) by SAMOVA. The suggested structure by SAMOVA for SSR data resulted in the highest genetic variation occurring within populations and for mtCOI data, all the genetic variation occurred within populations (Table 8(a3) and 8(a4) respectively).
Based on the results from AMOVA for both markers, FST (0.00069) (SSR markers) and FST (−0.01605) (mtCOI), genetic differentiation was low.
Based on both SSR and mtCOI partial gene markers, STRUCTURE analysis identified K = 2 as the likely optimal number of ancestral genetic clusters, with all individuals categorized as one dominant cluster (blue colour) and different degrees of genetic contributions (orange colour) from a second minor genetic cluster (Figure 5a,b).

3.5. Isolation by Distance (IBD)

Isolation by distance (IBD) was calculated for each agro-ecological zone independently. During this inference, agro-ecological zones were regarded as populations as detailed in Table 2. There was no evidence of IBD (SSR data) detected (Mantel test for agro-ecological zones A1 to E1 respectively: ((A`1) r = −0.061, p = 0.222; (B1) r = 0.107, p = 0.205; (C1) r = 0.204, p = 0.332; (D1) r = −0.038, p = 0.373; (E1) r = −0.073, p = 30.363 (Figure 6)). Isolation by distance was also calculated for mtCOI partial gene data for each agro-ecological zone independently. There was no evidence of IBD detected also ((A2) r = −0.044, p = 0.342; (B2) r = −0.128, p = 0.685; (C2) r = −0.050, p = 0.755; (E2) r = −0.153, p = 0.273) except in one of the populations ((D2) agro-ecological zones ((D). Agro-ecological zone D2) ((r = 0.161, p = 0.073).

4. Discussion

Prior to designing robust and effective pest management strategies, information about the level of individuals’ interaction among locations [53], in addition to pests’ life cycles and their predators [54], need to be availed. The former (pests’ interaction among locations) can be availed through the analysis of the population genetic structure of the pest in question [55].
Our study is the first to report the spatial population genetic structure of O. mutabilis in Uganda. Based on both types of markers, we found that the highest genetic variation occurred within the populations, low among groups and none among groups within populations. While sample sizes could account for the low genetic variation detected among groups and none were detected among groups within-population level, this was likely also an indication that there was intermixing or interbreeding between individuals from different populations. This aspect can be attributed to the capabilities for migration of our studied species as beetles are reported to fly short distances, although such trivial flights can also lead to the covering of longer distances [55,56]. Our results concur with a previous report which showed that highly dispersing phytophagous arthropods are characterized by homogenizing effects due to gene flow over distant localities [57]. Another leaf beetle, Cerotoma trifurcata, from the mid-western United States [55] exhibited similar high variability within samples, while the lowest degree of variation was among populations, indicating substantial levels of gene flow. Low levels of variation among populations can be attributed to the presence of few or no geographical barriers that impede gene flow. The absence of geographical barriers among C. trifucarta communities in the mid-western United States facilitated dispersal and reduced geographical fragmentation and genetic differentiation [55].
There was no evidence of genetic differentiation for O. mutabilis for both microsatellite and mtCOI partial gene molecular markers indicated by insignificant FST values. This suggests that migration occurs across O. mutabilis populations, particularly those separated by large geographic distances. Local and long-distance migrations of O. mutabilis may be aided by the presence of other host plants and staggered bean planting. Beans are grown at different times of the year in the study area as one of the cultural practices to manage BLBs [8], and this practice has the potential to drive BLB migration between gardens. However, no information is available to confirm whether the presence of beans in the gardens attracts BLBs feeding on beans from distant gardens. It is also possible that BLBs prefer certain bean varieties over others, causing them to migrate between gardens. As a result, it will also be necessary to investigate if volatiles in bean plants can attract BLBs and encourage their migration as it has been reported to occur in Spodoptera littoralis caterpillars while feeding on maize [58].
Bean leaf beetles oviposit in the soil, and the eggs hatch into larvae, pupae, and finally adults [5,8]. However, it is unknown whether these eggs can be carried from one garden to another via farm tools such as pangas, hoes, gum boots, and so on.
Our study areas were limited only to Uganda, however, BLBs have been reported to occur elsewhere including in countries from East Africa (e.g., Tanzania) and West Africa (e.g., Senegal) [9]. It would be a greater opportunity to understand their genetic status as the East African Rift Valleys have been shown to support population substructure and/or early speciation in both invertebrates (e.g., [59,60,61]) and vertebrates [1,47]. It remains to be investigated whether the geographical distribution of O. mutabilis, which, crossed the Rift Valley, may have resulted in comparable incipient speciation as observed in the cassava whitefly Bemisia ‘SSA1’ species using a whole-genome analysis technique [59]. Mitochondrial DNA markers, such as the partial COI gene, are ineffective for distinguishing between closely related subspecies [59,62,63].
STRUCTURE analyses from SSR markers revealed genetic mixing from two ancestral genetic lineages in our O. mutabilis samples. Based on limited nuclear (SSR) markers and the maternally inherited mtCOI marker, our findings showed that the genetic composition of the current populations could be explained as the outcome of genetic contributions from two hypothetical ancestral O. mutabilis genetic clusters. This study detected no evidence of genetic differentiation indicating that the Ugandan population likely represents a single panmictic population. To gain a better knowledge of the evolutionary genetics and landscape adaptability of this major agricultural coleopteran pest complex, an in-depth population structure study based on whole-genome sequencing would be desirable.
The inability to distinguish species based on physical characteristics is a hindrance to effective pest management. Both O. mutabilis and O. proteus are important bean pest species that are morphologically indistinguishable and can only be separated by dissection and study of the male genitalia [9], a process that does not distinguish between female species. The mtCOI gene can help identify morphologically challenging (e.g., [31,64]) and cryptic species (e.g., [65]) including Coleoptera (e.g., [66]). Molecular identification using the mtCOI partial gene, as described in this study, successfully distinguished between O. mutabilis and O. proteus that could not be distinguished by colour appearances alone during laboratory sample selection and should be used in future evolutionary genetic studies of Ootheca species.
Our IBD analysis results also showed that O. mutabilis has had little or no barriers to free migration in all populations. In research conducted by Krell et al. [56], a single C. trifurcata beetle travelled 4.9 km. These trivial flights in search of mates, oviposition sites, and food [56] may contribute to greater regional coverage over time, resulting in minimal genetic differentiation and substantial gene flow [55]. Behavioural research findings may aid in identifying features that contribute to analogous IBD and gene flow findings in our target species.

5. Conclusions

This study found that almost all the genetic variation in our target species, O. mutabilis, occurred within populations which, we attributed to dispersal that facilitated genetic mixing between populations. Ootheca mutabilis were divided into two population groups by SAMOVA at K = 2. When inferred, the structure suggested by SAMOVA resulted in almost all the genetic variation occurring within the populations. There was no evidence of genetic differentiation as seen from insignificant FST and we attributed it to gene flow between different O. mutabilis populations. On inferring the Bayesian clustering, STRUCTURE at K = 2 categorized all the samples as one dominant cluster (admixture) and different degrees of genetic contributions from a second minor genetic cluster. Therefore, the Ugandan populations of O. mutabilis likely represent a single panmictic population with genetic contributions from two ancestral lineages. Our study provides a baseline for future evolutionary and functional genomic studies to generate a better understanding of host-plant adaptation, insecticide resistance management, and the development of integrated pest management control measures for this important pest. Future research on O. mutabilis should increase coverage of samples from other African countries via a whole-genome sequencing approach.

Author Contributions

Conceptualization, M.H.O.; Data curation, D.K., W.T.T. and M.H.O.; Formal analysis, D.K., W.T.T. and I.S.K.; Funding acquisition, M.H.O. and S.E.S.; Investigation, D.K., V.P.A., W.T.T., P.N. and M.H.O.; Methodology, D.K., V.P.A., W.T.T., P.N., G.S., S.E.S. and M.H.O.; Project administration, S.T.N. and M.H.O.; Resources, D.K., S.O., S.T.N. and M.H.O.; Supervision, S.T.N. and M.H.O.; Visualization, D.K., M.H.O. and W.T.T.; Writing—original draft, D.K., V.P.A., W.T.T., P.N., P.P., G.S., S.E.S. and M.H.O.; Writing—review & editing, D.K., V.P.A., W.T.T., P.N., P.P., S.O., S.T.N., I.S.K., G.S., S.E.S. and M.H.O. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Bill & Melinda Gates Foundation’s Programme for Emerging Agricultural Research Leaders (OPP1131470_2015).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Acknowledgments

We thank the local communities for their support and advice during sample collection.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Minja, M. Promotion of integrated pest management strategies of major insect pests of Phaseolus beans in hillside systems in eastern and southern Africa. In Crop Protection Programme; Final Technical Report; DFID; CIAT: Kigali, Rwanda, 2005. [Google Scholar]
  2. Adipala, E.; Omongo, C.A.; Sabiti, A.; Obuo, J.E.; Edema, R.; Bua, B.; Atyang, A.; Nsubuga, E.N.; Ogenga-latigo, M.W. Pests and diseases on cowpea in Uganda: Experiences from a diagnostic survey. Afr. Crop Sci. J. 1999, 7, 465–478. [Google Scholar] [CrossRef]
  3. Kyamanywa, S.; Mukibi, J.; Otim, M. Use of trap crops for management of bean leaf beetles (Ootheca spp.) in Apac district of Uganda. In African Crop Science Conference Proceedings; African Crop Science Society: Kampala, Uganda, 2001; Volume 5, pp. 167–170. [Google Scholar]
  4. Halerimana, C.; Kyamanywa, S.; Olaboro, S.; Paparu, P.; Nkalubo, S.T.; Colvin, J.; Cheke, R.A.; Wagner, T.; Seal, S.E.; Kriticos, D.J.; et al. Distribution and relative abundance of bean leaf beetles (Ootheca spp.) (Insecta: Coleoptera: Chrysomelidae) in Uganda. Insects 2021, 12, 1048. [Google Scholar] [CrossRef] [PubMed]
  5. Ampofo, J.K.O.; Massomo, S.M. Host effects on bean foliage beetle (Coleoptera: Chrysomelidae) emergence pattern in northern Tanzania. Bean Improv. Cooperative. Annu. Rep. 1999, 109, 110. [Google Scholar]
  6. Allen David, J.; Ampofo, J.K.O.; Wortmann Charles, S. Pests, Diseases, and Nutritional Disorders of the Common Bean in Africa: A Field Guide; CIAT Publication No. 260; Centro Internacional de Agricultura Tropical (CIAT): Wageningen, NE, USA; Technical Centre for Agricultural and Rural Cooperation (TCA): Cali, CO, USA, 1996; 132p. [Google Scholar]
  7. Sastry, K.S. Seed-Borne Plant Virus Diseases; Springer Science & Business Media: Berlin/Heidelberg, Germany, 2013; ISBN 978-81-322-0813-6. (eBook). [Google Scholar]
  8. Halerimana, C. Distribution of bean leaf beetles and associated yield losses. Master’s Thesis, Makerere University, Kampala, Uganda, 2019. [Google Scholar]
  9. Kortenhaus, S.; Wagner, T. Revision of Ootheca Chevrolat, 1837 from tropical Africa–Redescriptions, descriptions of new species and identification key (Coleoptera: Chrysomelidae: Galerucinae). Zootaxa 2010, 2659, 1–52. [Google Scholar] [CrossRef]
  10. Crozier, R.H.; Oldroyd, B.P.; Tay, W.T.; Kaufmann, B.E.; Johnson, R.N.; Carew, M.E.; Jennings, K.M. Molecular advances in understanding social insect population structure. Electrophoresis 1997, 18, 1672–1675. [Google Scholar] [CrossRef] [PubMed]
  11. Moges, A.D.; Admassu, B.; Belew, D.; Yesuf, M.; Njuguna, J.; Kyalo, M.; Ghimire, S.R. Development of microsatellite markers and analysis of genetic diversity and population structure of Colletotrichum gloeosporioides from Ethiopia. PLoS ONE 2016, 11, e0151257. [Google Scholar] [CrossRef] [PubMed]
  12. Tay, W.T.; Crozier, R.H. Mating behaviour of Rhytidoponera sp. 12 ants inferred from microsatellite analysis. Mol. Ecol. 2001, 10, 167–173. [Google Scholar] [CrossRef]
  13. Tay, W.T.; Miettinen, M.; Kaitala, A. Do male golden egg bugs carry eggs they have fertilized? A microsatellite analysis. Behav. Ecol. 2003, 14, 481–485. [Google Scholar] [CrossRef] [Green Version]
  14. Hammond, R.L.; Bourke, A.F.; Bruford, M.W. Mating frequency and mating system of the polygynous ant, Leptothorax acervorum. Mol. Ecol. 2001, 10, 2719–2728. [Google Scholar] [CrossRef]
  15. Renault, D. A Review of the phenotypic traits associated with insect dispersal polymorphism, and experimental designs for sorting out resident and disperser phenotypes. Insects 2020, 11, 214. [Google Scholar] [CrossRef] [Green Version]
  16. Lombaert, E.; Boll, R.; Lapchin, L. Dispersal strategies of phytophagous insects at a local scale: Adaptive potential of aphids in an agricultural environment. BMC Evol. Biol. 2006, 6, 75. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Roberts, J.M.K.; Anderson, D.; Tay, W.T. Multiple host shifts by the emerging honeybee parasite. Varroa Jacobsoni. Mol. Ecol. 2015, 24, 2379–2391. [Google Scholar] [CrossRef] [PubMed]
  18. Tay, W.T.; Behere, G.T.; Batterham, P.; Heckel, D.G. Generation of microsatellite repeat families by RTE retrotransposons in lepidopteran genomes. BMC Evol. Biol. 2010, 10, 144. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Gordon, K.; Tay, W.T.; Collinge, D.; Williams, A.; Batterham, P. Genetics and molecular biology of the major crop pest genus Helicoverpa. In Molecular Biology and Genetics of the Lepidoptera; Goldsmith, M.R., Marec, F., Eds.; CRC Press: Boca Raton, FL, USA, 2009; pp. 219–238. [Google Scholar]
  20. Sinama, M.; Dubut, V.; Costedoat, C.; Gilles, A.; Junker, M.; Malausa, T.; Jean-françois, M.; Nève, G.; Pech, N.; Schmitt, T.; et al. Challenges of microsatellite development in Lepidoptera: Euphydryas aurinia (Nymphalidae) as a case study. EJE 2011, 108, 261–266. [Google Scholar] [CrossRef] [Green Version]
  21. Zhang, D.X. Lepidopteran microsatellite DNA: Redundant but promising. Trends Ecol. Evol. 2004, 19, 507–509. [Google Scholar] [CrossRef]
  22. Mckeown, N.; Harvey, D.; Healey, A.; Skujina, I.; Cox, K.; Gange, A.C.; Shaw, P. Isolation and characterisation of the first microsatellite markers for the European stag beetle, Lucanus cervus (Coleoptera: Lucanidae). Eur. J. Entomol. 2018, 115, 620–623. [Google Scholar] [CrossRef]
  23. Drury, D.W.; Siniard, A.L.; Wade, M.J. Genetic differentiation among wild populations of Tribolium castaneum estimated using microsatellite markers. J. Hered. 2009, 100, 732–741. [Google Scholar] [CrossRef]
  24. Ekblom, R.; Galindo, J. Applications of next generation sequencing in molecular ecology of non-model organisms. Heredity 2011, 107, 1. [Google Scholar] [CrossRef] [Green Version]
  25. Weber, J.L.; Wong, C. Mutation of human short tandem repeats. Hum. Mol. Genet. 1993, 2, 1123–1128. [Google Scholar] [CrossRef]
  26. Hancock, J.M. Microsatellites and other simple sequences: Genomic context and mutational mechanisms. In Microsatellites Evolution, Applications; Goldstein, D.B., Schlötterer, C., Eds.; New York Oxford University Press: New York, NY, USA, 1999; pp. 1–9. [Google Scholar]
  27. Randi, E.; Pierpaoli, M.; Danilkin, A. Mitochondrial DNA polymorphism in populations of Siberian and European roe deer (Capreolus pygargus and C. capreolus). Heredity 1998, 80, 429–437. [Google Scholar] [CrossRef] [PubMed]
  28. Papura, D.; Burban, C.; van Helden, M.; Giresse, X.; Nusillard, B.; Guillemaud, T.; Kerdelhué, C. Microsatellite and Mitochondrial Data Provide Evidence for a Single Major Introduction for the Neartic Leafhopper Scaphoideus titanus in Europe. PLoS ONE 2012, 7, e36882. [Google Scholar] [CrossRef] [Green Version]
  29. Mugerwa, H.; Colvin, J.; Alicai, T.; Omongo, C.A.; Kabaalu, R.; Visendi, P.; Seal, S.E. Genetic diversity of whitefly (Bemisia spp.) on crop and uncultivated plants in Uganda: Implications for the control of this devastating pest species complex in Africa. J. Pest Sci. 2021, 94, 1307–1330. [Google Scholar] [CrossRef] [PubMed]
  30. Otim, M.H.; Adumo Aropet, S.; Opio, M.; Kanyesigye, D.; Nakelet Opolot, H.; Tay, W.T. Parasitoid distribution and parasitism of the fall armyworm Spodoptera frugiperda (Lepidoptera: Noctuidae) in different maize producing regions of Uganda. Insects 2021, 12, 121. [Google Scholar] [CrossRef] [PubMed]
  31. Otim, M.H.; Tay, W.T.; Walsh, T.K.; Kanyesigye, D.; Adumo, S.; Abongosi, J.; Ochen, S.; Sserumaga, J.; Alibu, S.; Abalo, G. Detection of sister-species in invasive populations of the fall armyworm Spodoptera frugiperda (Lepidoptera: Noctuidae) from Uganda. PLoS ONE 2018, 13, e0194571. [Google Scholar]
  32. Staden, R.; Beal, K.F.; Bonfield, J.K. The staden package, 1998. In Bioinformatics Methods and Protocols; Springer: Berlin/Heidelberg, Germany, 2000; pp. 115–130. [Google Scholar]
  33. Andrews, S. FastQC: A Quality Control Tool for High Throughput Sequence Data. 2010. Available online: https://www.bioinformatics.babraham.ac.uk/projects/fastqc/ (accessed on 15 February 2018).
  34. Kearse, M.; Moir, R.; Wilson, A.; Stones-havas, S.; Cheung, M.; Sturrock, S.; Buxton, S.; Cooper, A.; Markowitz, S.; Duran, C. Geneious Basic: An integrated and extendable desktop software platform for the organization and analysis of sequence data. Bioinformatics 2012, 28, 1647–1649. [Google Scholar] [CrossRef] [PubMed]
  35. Martins, W.S.; Lucas, D.C.S.; De Souza Neves, K.F.; Bertioli, D.J. WebSat A web software for microsatellite marker development. Bioinformation 2009, 3, 282. [Google Scholar] [CrossRef] [Green Version]
  36. Rozen, S.; Skaletsky, H. Primer3 on the WWW for general users and for biologist programmers. Methods Mol. Biol. 2000, 132, 365–386. [Google Scholar]
  37. Taubs, G. Sense from sequences: Stephen F. Altschul on bettering BLAST. Sci. Watch. 2000, 11, 3–4. [Google Scholar]
  38. Hulce, D.; Li, X.; Snyder-Leiby, T.; Liu, C.J. GeneMarker® genotyping software: Tools to increase the statistical power of DNA fragment analysis. J. Biomol. Tech. 2011, 22, S35. [Google Scholar]
  39. Van Oosterhout, C.; Hutchinson, W.F.; Wills, D.P.; Shipley, P. MICRO-CHECKER: Software for identifying and correcting genotyping errors in microsatellite data. Mol. Ecol. Notes 2004, 4, 535–538. [Google Scholar] [CrossRef]
  40. Gutiérrez, J.P.; Royo, L.; Álvarez, I.; Goyache, F. MolKin v2. 0: A computer program for genetic analysis of populations using molecular coancestry information. J. Hered. 2005, 96, 718–721. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  41. Rozas, J.; Ferrer-Mata, A.; Sánchez-Delbarrio, J.C.; Guirao-Rico, S.; Librado, P.; Ramos-Onsins, S.E.; Sánchez-Gracia, A. DnaSP 6: DNA sequence polymorphism analysis of large data sets. Mol. Biol. Evol. 2017, 34, 3299–3302. [Google Scholar] [CrossRef] [PubMed]
  42. Excoffier, L.; Lischer, H.E.L. Arlequin suite ver 3.5: A new series of programs to perform population genetics analyses under Linux and Windows. Mol. Ecol. Res. 2010, 10, 564–567. [Google Scholar] [CrossRef] [PubMed]
  43. Dupanloup, I.; Schneidera, S.; Excoffier, N.L. A simulated annealing approach to define the genetic structure of populations. Mol. Ecol. 2002, 11, 2571–2581. [Google Scholar] [CrossRef] [PubMed]
  44. Pritchard, J.K.; Stephens, M.; Donnelly, P. Inference of population structure using multilocus genotype data. Genetics 2000, 155, 945–959. [Google Scholar] [CrossRef]
  45. Earl, D.A.; Vonholdt, B.M. Structure Harvester: A website and program for visualizing Structure output and implementing the Evanno method. Conserv. Genet. Resour. 2012, 4, 359–361. [Google Scholar] [CrossRef]
  46. Evanno, G.; Regnaut, S.; Goudet, J. Detecting the number of clusters of individuals using the software Structure: A simulation study. Mol. Ecol. 2005, 14, 2611–2620. [Google Scholar] [CrossRef] [Green Version]
  47. Kopelman, N.M.; Mayzel, J.; Jakobsson, M.; Rosenberg, N.A.; Mayrose, I. Clumpak: A program for identifying clustering modes and packaging population structure inferences across K. Mol. Ecol. Resour. 2015, 15, 1179–1191. [Google Scholar] [CrossRef] [Green Version]
  48. Leigh, J.W.; Bryant, D. POPART: Full-feature software for haplotype network construction. Methods Ecol. Evol. 2015, 6, 1110–1116. [Google Scholar] [CrossRef]
  49. Mantel, N.; Haenszel, W. Statistical aspects of the analysis of data from retrospective studies of disease. J. Natl. Cancer Inst. 1959, 22, 719–748. [Google Scholar] [PubMed]
  50. Peakall, R.; Smouse, P.E. GENALEX 6: Genetic analysis in Excel. Population genetic software for teaching and research. Mol. Ecol. Resour. 2006, 6, 288–295. [Google Scholar] [CrossRef]
  51. Shete, S.; Tiwari, H.; Elston, R.C. On estimating the heterozygosity and polymorphism information content value. Theor. Popul. Biol. 2000, 57, 265–271. [Google Scholar] [CrossRef] [PubMed]
  52. Nei, M. Estimation of average heterozygosity and genetic distance from a small number of individuals. Genetics 1978, 23, 341–369. [Google Scholar] [CrossRef] [PubMed]
  53. De León, J.H.; Jones, W.A.; Morgan, D.J. Population genetic structure of Homalodisca coagulata (Homoptera: Cicadellidae), the vector of the bacterium Xylella fastidiosa causing Pierce’s disease in grapevines. Ann. Entomol. Soc. Am. 2004, 97, 809–818. [Google Scholar]
  54. Alam, M.Z.; Crump, A.R.; Haque, M.M.; Islam, M.S.; Hossain, E.; Hasan, S.B.; Hasan, S.B.; Hossain, M.S. Effects of Integrated Pest Management on Pest Damage and Yield Components in a Rice Agro-Ecosystem in the Barisal Region of Bangladesh. Front. Environ. Sci. 2016, 4, 22. [Google Scholar] [CrossRef] [Green Version]
  55. Tiroesele, B.; Skoda, S.R.; Hunt, T.E.; Lee, D.J.; Molina-ochoa, J.; Foster, J.E. Population structure, genetic variability, and gene flow of the bean leaf beetle, Cerotoma trifurcata, in the Midwestern United States. J. Insect Sci. 2014, 14, 62. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  56. Krell, R.K.; Wilson, T.A.; Pedigo, L.P.; Rice, M.E. Characterization of bean leaf beetle (Coleoptera: Chrysomelidae) flight capacity. J. Kans. Entomol. Soc. 2003, 406–416. [Google Scholar]
  57. The influence of dispersal and diet breadth on patterns of genetic isolation by distance in phytophagous insects. Am. Nat. 1998, 152, 428–446. [CrossRef]
  58. von Mérey, G.E.; Veyrat, N.; D’Alessandro, M.; Turlings, T.C.J. Herbivore-induced maize leaf volatiles affect attraction and feeding behavior of Spodoptera littoralis caterpillars. Front. Plant Sci. 2013, 4, 209. [Google Scholar] [CrossRef] [Green Version]
  59. Elfekih, S.; Tay, W.T.; Polaszek, A.; Gordon, K.H.J.; Kunz, D.; Macfadyen, S.; De Barro, P.J. On species delimitation, hybridization and population structure of cassava whitefly in Africa. Sci. Rep. 2021, 11, 7923. [Google Scholar] [CrossRef]
  60. Lehmann, T.; Hawley, W.; Grebert, H.; Danga, M.; Atieli, F.; Collins, F. The rift valley complex as a barrier to gene flow for Anopheles gambiae in Kenya. J. Hered. 1999, 90, 613–621. [Google Scholar] [CrossRef] [Green Version]
  61. Sezonlin, M.; Dupas, S.; Le Rü, B.; Le Gall, P.; Moyal, P.; Calatayud, P.A.; Giffard, I.; Faure, N.; Silvain, J.F. Phylogeography and population genetics of the maize stalk borer Busseola fusca (Lepidoptera, Noctuidae) in sub-Saharan Africa. Mol. Ecol. 2006, 15, 407–420. [Google Scholar] [CrossRef] [PubMed]
  62. Anderson, C.; Tay, W.T.; Mcgaughran, A.; Gordon, K.; Walsh, T.K. Population structure and gene flow in the global pest, Helicoverpa armigera. Mol. Ecol. 2016, 25, 5296–5311. [Google Scholar] [CrossRef] [Green Version]
  63. Behere, G.T.; Tay, W.T.; Russell, D.A.; Heckel, D.G.; Appleton, B.R.; Kranthi, K.R.; Batterham, P. Mitochondrial DNA analysis of field populations of Helicoverpa armigera (Lepidoptera: Noctuidae) and of its relationship to H. zea. BMC Evol. Biol. 2007, 7, 117. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  64. Behere, G.; Tay, W.; Russell, D.; Batterham, P. Molecular markers to discriminate among four pest species of Helicoverpa (Lepidoptera: Noctuidae). Bull. Entomol. Res. 2008, 98, 599. [Google Scholar] [CrossRef] [PubMed]
  65. Elfekih, S.; Tay, W.T.; Gordon, K.; Court, L.N.; De Barro, P.J. Standardized molecular diagnostic tool for the identification of cryptic species within the Bemisia tabaci complex. Pest Manag. Sci. 2018, 74, 170–173. [Google Scholar] [CrossRef]
  66. Tay, W.T.; Beckett, S.J.; De Barro, P.J. Phosphine resistance in Australian Cryptolestes species (Coleoptera: Laemophloeidae): Perspectives from mitochondrial DNA cytochrome oxidase I analysis. Pest Manag. Sci. 2016, 72, 1250–1259. [Google Scholar] [CrossRef]
Figure 1. Districts of Uganda collection locations of bean leaf beetle samples used in the study.
Figure 1. Districts of Uganda collection locations of bean leaf beetle samples used in the study.
Insects 13 00543 g001
Figure 2. Different colourmorphs of Ootheca mutabilis. M1: O. mutabilis with elytra upper half black and lower half brownish; M2: O. mutabilis with black elytra; M3: O. mutabilis with brown elytra.
Figure 2. Different colourmorphs of Ootheca mutabilis. M1: O. mutabilis with elytra upper half black and lower half brownish; M2: O. mutabilis with black elytra; M3: O. mutabilis with brown elytra.
Insects 13 00543 g002
Figure 3. Minimum spanning haplotype network showing evolutionary relationships among haplotypes in different agro-ecological zones (a), and different districts (b). Each small black line along a connecting line represents a change of one base pair. Haplotypes are colour coded according to the population. In haplotype network (a) (red represents population (A) Northern moist farmlands, green represents population (B) Western mid-altitude farmlands, purple represents population (C) Central wooden savannah, yellow represents population (D) Southern and Eastern Lake Kyoga basin and pink represents population (E) North-western farmlands). Circle sizes correspond to the haplotype numbers.
Figure 3. Minimum spanning haplotype network showing evolutionary relationships among haplotypes in different agro-ecological zones (a), and different districts (b). Each small black line along a connecting line represents a change of one base pair. Haplotypes are colour coded according to the population. In haplotype network (a) (red represents population (A) Northern moist farmlands, green represents population (B) Western mid-altitude farmlands, purple represents population (C) Central wooden savannah, yellow represents population (D) Southern and Eastern Lake Kyoga basin and pink represents population (E) North-western farmlands). Circle sizes correspond to the haplotype numbers.
Insects 13 00543 g003
Figure 4. Neighbour Joining phylogenetic tree of the 20 O. mutabilis haplotypes (red) found in the study, and reference sequences downloaded from NCBI as follows: O. mutabilis (blue), O. proteus (pink), O. proteus (green, morphologically similar samples to O. mutabilis separated after DNA-barcoding) and the outgroup Longitarsus tabidus.
Figure 4. Neighbour Joining phylogenetic tree of the 20 O. mutabilis haplotypes (red) found in the study, and reference sequences downloaded from NCBI as follows: O. mutabilis (blue), O. proteus (pink), O. proteus (green, morphologically similar samples to O. mutabilis separated after DNA-barcoding) and the outgroup Longitarsus tabidus.
Insects 13 00543 g004
Figure 5. (a). Population structure across five and sixteen populations agro-ecological zones and districts respectively of Ootheca mutabilis obtained using structure (for SSR data and (b)) for mtCO1 partial gene sequences both at K = 2. Vertical bars represent individuals. Numbers on the horizontal axis represent agro-ecological zones; (1) Northern moist farmlands (A), (2) Western mid-altitude farmlands (B), (3) Central wooden savannah (C), (4) Southern and Eastern Lake Kyoga basin (D), and (5) North-western farmlands (E).
Figure 5. (a). Population structure across five and sixteen populations agro-ecological zones and districts respectively of Ootheca mutabilis obtained using structure (for SSR data and (b)) for mtCO1 partial gene sequences both at K = 2. Vertical bars represent individuals. Numbers on the horizontal axis represent agro-ecological zones; (1) Northern moist farmlands (A), (2) Western mid-altitude farmlands (B), (3) Central wooden savannah (C), (4) Southern and Eastern Lake Kyoga basin (D), and (5) North-western farmlands (E).
Insects 13 00543 g005
Figure 6. Correlation between geographical and genetic distances of Ootheca mutabilis samples for each independent agro-ecological zone. (A1,A2) Northern moist farmlands (Agro-ecological zone A), (B1,B2) Western mid-altitude farmlands (Agro-ecological zone B), (C1,C2) Central wooden savannah (Agro-ecological zone C), (D1,D2) Southern and Eastern Lake Kyoga basin (Agro-ecological zone D) and (E1,E2) North-western farmlands (Agro-ecological zone E)). (A1E1) represent SSR data and (A2E2) represent mtCOI partial gene data.
Figure 6. Correlation between geographical and genetic distances of Ootheca mutabilis samples for each independent agro-ecological zone. (A1,A2) Northern moist farmlands (Agro-ecological zone A), (B1,B2) Western mid-altitude farmlands (Agro-ecological zone B), (C1,C2) Central wooden savannah (Agro-ecological zone C), (D1,D2) Southern and Eastern Lake Kyoga basin (Agro-ecological zone D) and (E1,E2) North-western farmlands (Agro-ecological zone E)). (A1E1) represent SSR data and (A2E2) represent mtCOI partial gene data.
Insects 13 00543 g006aInsects 13 00543 g006bInsects 13 00543 g006cInsects 13 00543 g006d
Table 1. Ugandan districts from which the samples were picked, GPS coordinates for each district and the total number of samples picked per district.
Table 1. Ugandan districts from which the samples were picked, GPS coordinates for each district and the total number of samples picked per district.
DistrictLatitudeLongitudeCollectorsSamples Collected
Dokolo2.0114833.1367Charles. H, Dalton. K, Sam. O, Sekandi. W58
Lira2.5019632.91012Charles. H, Dalton. K, Sam. O, Sekandi. W41
Oyam2.355932.60652Charles. H, Dalton. K, Sam. O, Sekandi. W30
Apac1.8844632.37174Charles. H, Dalton. K, Sam. O, Sekandi. W35
Amuru2.8149231.98196Charles. H, Dalton. K, Sam. O, Sekandi. W45
Gulu2.9603232.41548Charles. H, Dalton. K, Sam. O, Sekandi. W36
Nwoya2.6247332.14631Charles. H, Dalton. K, Sam. O, Sekandi. W42
Bulisa1.7619731.43002Charles. H, Dalton. K, Sam. O, Sekandi. W45
Hoima1.5009731.33689Charles. H, Dalton. K, Sam. O, Sekandi. W20
Nakasongola1.4692732.2695Charles. H, Dalton. K, Sam. O, Sekandi. W25
Amuria2.057433.50213Charles. H, Dalton. K, Sam. O, Sekandi. W37
Soroti5.3712021.94900Charles. H, Dalton. K, Sam. O, Sekandi. W22
Adjumani3.2563333.7836Charles. H, Dalton. K, Sam. O, Sekandi. W30
Zombo2.5159231.00431Charles. H, Dalton. K, Sam. O, Sekandi. W20
Koboko3.3820631.06935Charles. H, Dalton. K, Sam. O, Sekandi. W21
Moyo3.7036131.67643Charles. H, Dalton. K, Sam. O, Sekandi. W24
Arua3.1530931.01043Charles. H, Dalton. K, Sam. O, Sekandi. W23
Table 2. Agro-ecological zones, number of samples of O. mutabilis and O. proteus found per agro-ecological zone after DNA barcoding with mtCOI partial gene and Ugandan districts belonging to particular agro-ecological zones where bean leaf beetle samples were collected.
Table 2. Agro-ecological zones, number of samples of O. mutabilis and O. proteus found per agro-ecological zone after DNA barcoding with mtCOI partial gene and Ugandan districts belonging to particular agro-ecological zones where bean leaf beetle samples were collected.
Population CodeAgro-Ecological ZoneNumber of O. mutabilis Samples AnalyzedDistricts
O. mutabilisO. proteus
ANorthern moist farmlands50Dokolo
50Lira
60Oyam
50Apac
80Amuru
40Gulu
70Nwoya
BWestern mid-altitude farmlands100Bulisa
19Hoima
CCentral wooden savannah43Nakasongola
00Lwengo
DSouthern and Eastern Lake Kyoga basin90Amuria
40Soroti
ENorth-western farmlands50Adjumani
30Zombo
40Koboko
30Moyo
40Arua
Table 3. Characteristics of the five microsatellite loci developed and used for the Ootheca mutabilis population genetic study. NA (Number of alleles).
Table 3. Characteristics of the five microsatellite loci developed and used for the Ootheca mutabilis population genetic study. NA (Number of alleles).
Locus NameMotifSize (bp)Primer Sequence 5′ Insects 13 00543 i001 3′GenBank Accession NumberNATm (°C)Fluorescent Label
BLB2_om1(GAT)2(CAA)11343–365F: TCAACTACCACCATCACAAACC
R: CAATGTGGAGCAACTACGTCAT
MT0740939585′6-FAM
BLB2_om17(CTT)10368–396F: CCAATCCGCTTCTCTATATCCA
R: GGAGCAATGTTATGCCTGATTT
MT07409416575′6-FAM
BLB2_om32(GACG)6160–195F: CATATAGCGAAAACCCGAAATC
R:AGAAGTACAAGTATGGCCCGAA
MT07409621585′6-FAM
BLB2_om33(ACA)5.(ACG).(ACA)16256–288F: ATTGAAAGTTGTATCGGTCGCT
R: CTTGACATGAAAACGAGATCCA
MT0740954585′HEX
BLB2_om66(AGT)2(AGC)7337–345F: CTATGGTCGTTTTCTCCGACAT
R: GACGTTTCTTCTCGGTTGTAGC
MT0740978605′HEX
Table 4. Diversity indices of the five microsatellite loci developed for Ootheca mutabilis. Heterozygosity and polymorphic information content are He and PIC, respectively.
Table 4. Diversity indices of the five microsatellite loci developed for Ootheca mutabilis. Heterozygosity and polymorphic information content are He and PIC, respectively.
Locus NameHePIC
BLB_om10.7575.0%
BLB_om170.6558.8%
BLB_om320.8078.5%
BLB_om330.8483.1%
BLB_om660.5950.1%
Average0.7369.1%
Table 5. Gene diversity in the five populations of Ootheca mutabilis. Alphabet letters represent populations; (A) Northern moist farmlands, (B) Western mid-altitude farmlands, (C) Central wooden savannah, (D) Southern and Eastern Lake Kyoga basin and € North-western farmlands. Observed heterozygosity (Ho); expected heterozygosity (He), and polymorphic information content (PIC) are shown.
Table 5. Gene diversity in the five populations of Ootheca mutabilis. Alphabet letters represent populations; (A) Northern moist farmlands, (B) Western mid-altitude farmlands, (C) Central wooden savannah, (D) Southern and Eastern Lake Kyoga basin and € North-western farmlands. Observed heterozygosity (Ho); expected heterozygosity (He), and polymorphic information content (PIC) are shown.
Population CodeHoHePIC
A0.820.7256.97%
B0.800.6851.48%
C0.750.6647.42%
D0.840.7357.72%
E0.780.7054.36%
Average0.800.7053.59%
Table 6. Summary of genetic diversity indices for 658 bp fragment of the mtCOI partial gene. Number of haplotypes (h), Haplotype diversity (Hd), Average number of differences (K), and Nucleotide diversity (π) for each population are provided.
Table 6. Summary of genetic diversity indices for 658 bp fragment of the mtCOI partial gene. Number of haplotypes (h), Haplotype diversity (Hd), Average number of differences (K), and Nucleotide diversity (π) for each population are provided.
PopulationNumber of SequencesNumber of Segregating SiteshHdKπ
Dokolo5120.40.40
Lira5340.91.20
Oyam6230.60.670
Apac5120.40.40
Amuru8640.641.680
Gulu4120.50.50
Nwoya7430.521.140
Bulisa10440.530.960
Nakasongola4220.510
Amuria9230.560.610
Soroti4120.50.50
Adjumani501000
Zombo3120.670.670
Koboko4230.8310
Moyo3220.671.330
Arua4220.510
Table 7. (a1): Non-hierarchical AMOVA for all the sixteen populations of O. mutabilis, and (a2) hierarchical AMOVA for all the sixteen districts grouped into five agro-ecological zones based on SSR data; (b1) non-hierarchical AMOVA for all sixteen populations of O. mutabilis, and (b2) hierarchical AMOVA for all the sixteen populations grouped into their relative agro-ecological zones based on mtCOI data.
Table 7. (a1): Non-hierarchical AMOVA for all the sixteen populations of O. mutabilis, and (a2) hierarchical AMOVA for all the sixteen districts grouped into five agro-ecological zones based on SSR data; (b1) non-hierarchical AMOVA for all sixteen populations of O. mutabilis, and (b2) hierarchical AMOVA for all the sixteen populations grouped into their relative agro-ecological zones based on mtCOI data.
(a1)
Source of VariationSum of SquaresVariance ComponentsPercentage Variation (%)
Among populations27.240.00122 Va0.07
Within populations264.051.780 Vb99.93
Total291.291.790100
(a2)
Source of VariationSum of SquaresVariance ComponentsFixation Indices
Among groups10.080.03722 VaFST = 0.00582, p = 0.00880
Among populations within groups17.16−0.027 VbFSC = −0.01523, p = 0.42326
Within populations264.051.785 VcFCT = 0.02073, p = 0.00587
Total291.291.795
(b1)
Source of Variationd.fSum of SquaresVariance Components
Among populations155.83−0.00671 Va
Within populations70264.050.42450 Vb
Total85291.290.418
(b2)
Source of Variationd.fSum of SquaresVariance ComponentsFixation Indices
Among groups41.710.00448 VaFST = −0.01333, p = 0.73900
Among populations within groups114.12−0.01007 VbFSC = −0.02429, p = 0.76051
Within populations7029.720.42450 VcFCT = 0.01070, p = 0.30108
Total8535.550.419
Table 8. (a1). Spatial analysis of molecular variance (SAMOVA) of SSR data and (a2) SAMOVA of mtCOI data showing the genetic diversity partitions in the O. mutabilis populations. (a3) and (a4) Analysis of molecular variance (AMOVA) of the genetic structure as suggested by SAMOVA for SSR and mtCOI data respectively.
Table 8. (a1). Spatial analysis of molecular variance (SAMOVA) of SSR data and (a2) SAMOVA of mtCOI data showing the genetic diversity partitions in the O. mutabilis populations. (a3) and (a4) Analysis of molecular variance (AMOVA) of the genetic structure as suggested by SAMOVA for SSR and mtCOI data respectively.
(a1)
Source of Variationd.fSum of SquaresVariance Components
Among groups157,895.354273.04972 Va
Among populations within groups14143,369.46387.00479 Vb
Within populations156935,627.35997.61091 Vc
Total1711,136,892.1110,657.67
FCT = 0.40094, p = 0.05767
(a2)
Source of Variationd.fSum of SquaresVariance Components
Among groups10.9284.9040.04245 Va
Among populations within groups144.9−0.01442 Vb
Within populations7029.720.42450 Vc
Total8535.550.45
FCT = 0.09381, p = 0.05474
(a3)
Source of VariationSum of SquaresVariance ComponentsPercentage Variation (%)
Among populations27.240.00122 Va0.07
Within populations264.051.78482 Vb99.93
Total291.291.79100
FST = 0.00069
(a4)
Source of Variationd.f.Sum of SquaresVariance ComponentsPercentage Variation (%)
Among populations155.83−0.00671 Va−1.61
Within populations7029.720.42450 Vb101.61
Total8535.550.42100
FST = −0.01605
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kanyesigye, D.; Alibu, V.P.; Tay, W.T.; Nalela, P.; Paparu, P.; Olaboro, S.; Nkalubo, S.T.; Kayondo, I.S.; Silva, G.; Seal, S.E.; et al. Population Genetic Structure of the Bean Leaf Beetle Ootheca mutabilis (Coleoptera: Chrysomelidae) in Uganda. Insects 2022, 13, 543. https://0-doi-org.brum.beds.ac.uk/10.3390/insects13060543

AMA Style

Kanyesigye D, Alibu VP, Tay WT, Nalela P, Paparu P, Olaboro S, Nkalubo ST, Kayondo IS, Silva G, Seal SE, et al. Population Genetic Structure of the Bean Leaf Beetle Ootheca mutabilis (Coleoptera: Chrysomelidae) in Uganda. Insects. 2022; 13(6):543. https://0-doi-org.brum.beds.ac.uk/10.3390/insects13060543

Chicago/Turabian Style

Kanyesigye, Dalton, Vincent Pius Alibu, Wee Tek Tay, Polycarp Nalela, Pamela Paparu, Samuel Olaboro, Stanley Tamusange Nkalubo, Ismail Siraj Kayondo, Gonçalo Silva, Susan E. Seal, and et al. 2022. "Population Genetic Structure of the Bean Leaf Beetle Ootheca mutabilis (Coleoptera: Chrysomelidae) in Uganda" Insects 13, no. 6: 543. https://0-doi-org.brum.beds.ac.uk/10.3390/insects13060543

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