Genetic diversity and phylogenetic relationships of tsetse flies of the palpalis group in Congo Brazzaville based on mitochondrial cox1 gene sequences

Despite the morphological characterization established in the 1950s and 1960s, the identity of extant taxa that make up Glossina fuscipes (s.l.) in the Congo remains questionable. Previous claims of overlap between G. fuscipes (believed to be G. f. quanzensis) and G. palpalis palpalis around Brazzaville city further complicate the taxonomic status and population dynamics of the two taxa. This study aimed to determine the phylogenetic relationships between G. fuscipes (s.l.) and G. p. palpalis and to assess genetic variation among G. fuscipes (s.l.) populations in Congo Brazzaville.

Phylogenetic analysis revealed minor differences between G. fuscipes (s.l.) and G. p. palpalis. Genetic diversity of G. fuscipes (s.l.) was high in the populations sampled except one. Genetic differentiation ranged from moderate to high among subpopulations. There was a restricted gene flow between G. fuscipes (s.l.) populations in the north and central part of the country. Genetic signatures based on cox1 showed recent expansion and recovery of G. fuscipes (s.l.) populations from previous bottlenecks. To fully understand the species distribution limits, we recommend further studies involving a wider sampling scheme including the swampy Mossaka focus for G. fuscipes (s.l.) and the entire range of G. p. palpalis in South Congo.

We collected 263 G. fuscipes (s.l.) from northern and central regions, and 65 G. p. palpalis from southern part of the country. The mitochondrial cytochrome c oxidase subunit 1 (cox1) gene was amplified using taxa-specific primer pairs. Sequence data were analyzed in DnaSP and Arlequin to assess the genetic diversity, differentiation and demographic history of G. fuscipes (s.l.) populations.

The general BLAST analysis yielded a similarity of 99% for G. fuscipes (s.l.) and G. p. palpalis. BLASTn analysis for G. fuscipes (s.l.) showed > 98% identity with GenBank sequences for G. fuscipes (s.l.), with BEMB population showing 100% similarity with G. f. fuscipes. Glossina fuscipes (s.l.) populations showed high haplotype diversity (H = 46, Hd = 0.884), moderate nucleotide diversity ( = 0.012) and moderate (FST = 0.072) to high (FST = 0.152) genetic differentiation. Most of the genetic variation (89.73%) was maintained within populations. The mismatch analysis and neutrality tests indicated recent tsetse population expansions. Conclusions: Phylogenetic analysis revealed minor differences between G. fuscipes (s.l.) and G. p. palpalis. Genetic diversity of G. fuscipes (s.l.) was high in the populations sampled except one. Genetic differentiation ranged from moderate to high among subpopulations. There was a restricted gene flow between G. fuscipes (s.l.) populations in the north and central part of the country. Genetic signatures based on cox1 showed recent expansion and recovery of G. fuscipes (s.l.) populations from previous bottlenecks. To fully understand the species distribution limits, we recommend further studies involving a wider sampling scheme including the swampy Mossaka focus for G. fuscipes (s.l.) and the entire range of G. p. palpalis in South Congo.