Cardinium symbiosis as a potential confounder of mtDNA based phylogeographic inference in Culicoides imicola (Diptera: Ceratopogonidae), a vector of veterinary viruses
Abstract Background: Culicoides imicola (Diptera: Ceratopogonidae) is an important Afrotropical and Palearctic vector of disease, transmitting viruses of animal health and economic significance. The apparent incursions of C. imicola into mainland Europe via wind-movement events has made it important to trace this species to better predict new areas of arbovirus outbreaks. A widely used method for tracking dispersal patterns of C. imicola employs a phylogeographic approach anchored on the mtDNA marker COI (cytochrome c oxidase subunit I). However, a problem with this approach is that maternally-inherited symbiotic bacteria can alter the frequency of COI mitochondrial haplotypes (mitotypes), masking the true patterns of movement and gene flow.Methods: In this study, the mtDNA structure of C. imicola in relation to Cardinium infection status was investigated through haplotype network analysis. COI Sanger sequences from infected and uninfected individuals were first compared, before extending the haplotype network to include mitotypes from a geographic range where Cardinium infection has previously been investigated.Results: The mitotype network of a South African population, containing both Cardinium-infected and uninfected individuals, demonstrated the presence of two broad mitotype groups. All Cardinium-infected specimens fell into one group (Fisher’s exact test, P<0.001) demonstrating a linkage disequilibrium between symbiont and mitochondria. Furthermore, by extending this haplotype network to include other C. imicola populations from the Mediterranean basin, we revealed mitotype variation between the Eastern and Western Mediterranean basins (EMB and WMB) mirrored Cardinium-infection heterogeneity. Conclusions: These observations suggest that the linkage disequilibrium of Cardinium and mitochondria reflects symbiont gene flow within the Mediterranean basin but may not assist in elucidating host gene flow. Subsequently, we urge caution on the single usage of the COI marker to determine population structure and movement in C. imicola, and instead suggest the complementary utilisation of additional molecular markers.