Range-wide genetic structure and demographic history in the bat ectoparasite Cimex adjunctus

Host traits and environment interact to determine persistence of bat populations impacted by white-nose syndrome

Emerging infectious diseases have resulted in severe population declines across diverse taxa. In some instances, despite attributes associated with high extinction risk, disease emergence and host declines are followed by host stabilisation for unknown reasons. While host, pathogen, and the environment are recognised as important factors that interact to determine host-pathogen coexistence, they are often considered independently. Here, we use a translocation experiment to disentangle the role of host traits and environmental conditions in driving the persistence of remnant bat populations a decade after they declined 70-99% due to white-nose syndrome and subsequently stabilised. While survival was significantly higher than during the initial epidemic within all sites, protection from severe disease only existed within a narrow environmental space, suggesting host traits conducive to surviving disease are highly environmentally dependent. Ultimately, population persistence following pathogen invasion is the product of host-pathogen interactions that vary across a patchwork of environments.

Latrocimicinae completes the phylogeny of Cimicidae: meeting old morphologic data rather than modern host phylogeny

The family Cimicidae includes obligate hematophagous ectoparasites (bed bugs and their relatives) with high veterinary/medical importance. The evolutionary relationships of Cimicidae and their hosts have recently been reported in a phylogenetic context, but in the relevant study, one of the six subfamilies, the bat-specific Latrocimicinae, was not represented. In this study the only known species of Latrocimicinae, i.e., Latrocimex spectans, was analyzed with molecular and phylogenetic methods based on four (two nuclear and two mitochondrial) genetic markers. The completed subfamily-level phylogeny of Cimicidae showed that Latrocimicinae is most closely related to Haematosiphoninae (ectoparasites of birds and humans), with which it shares systematically important morphologic characters, but not hosts. Moreover, in the phylogenetic analyses, cimicid bugs that are known to infest phylogenetically distant bat hosts clustered together (e.g., Leptocimex and Stricticimex within Cacodminae), while cimicid subfamilies (Latrocimicinae, Primicimicinae) that are known to infest bat hosts from closely related superfamilies clustered distantly. In conclusion, adding Latrocimicinae significantly contributed to the resolution of the phylogeny of Cimicidae. The close phylogenetic relationship between Latrocimicinae and Haematosiphoninae is consistent with long-known morphologic data. At the same time, phylogenetic relationships of genera within subfamilies are inconsistent with the phylogeny of relevant hosts.

Sequencing whole mitochondrial genomes to assess genetic divergence between proposed silver-haired bat (Lasionycteris noctivagans) populations

The geographic distributions of eastern and western Lasionycteris noctivagans populations suggest they could be genetically isolated, but this has rarely been assessed using genetic data. Here, we evaluate this possibility by sequencing the complete mitochondrial genome of four silver-haired bats from eastern and western populations. The three usable mitogenomes were closely associated with other Vespertilionid bats and the phylogenetic tree revealed the two western individuals grouping together to form their own clade. Our results support the idea of small but significant genetic differences between eastern and western populations of these bats, but this should be tested further.

Bed Bugs (Hemiptera, Cimicidae): Overview of Classification, Evolution and Dispersion

The bed bugs (Cimex lectularius and C. hemipterus) have undergone a significant resurgence worldwide since the 1990s. A compilation of findings from a database, including 2650 scientific publications from seven major medical databases, allowed us to document main evolutionary events, from fossil evidence, dating from 11,000 years ago, until the present that has led to the current worldwide expansion of Cimicid species. We present the hypotheses on the possible dispersion pathways of bed bugs in light of the major historical and evolutionary events. A detailed classification of the Cimicidae family and finally, an illustrative map displaying the current distribution of known Cimex species in each geographical ecozone of Asia, Europe, Africa, the Americas, and Australia are presented.

Bed bugs: The move to humans as hosts

Cimicid insects, bed bugs and their allies, include about 100 species of blood-feeding ectoparasites. Among them, a few have become widespread and abundant pests of humans. Cimicids vary in their degree of specialization to hosts. Whereas most species specialize on insectivorous birds or bats, the common bed bug can feed on a range of distantly related host species, such as bats, humans, and chickens. We suggest that association with humans and generalism in bed bugs led to fundamentally different living conditions that fostered rapid growth and expansion of their populations. We propose that the evolutionary and ecological success of common bed bugs reflected exploitation of large homeothermic hosts (humans) that sheltered in buildings. This was a departure from congeners whose hosts are much smaller and often heterothermic. We argue that interesting insights into the biology of pest species may be obtained using an integrated view of their ecology and evolution.

The Klingon batbugs: Morphological adaptations in the primitive bat bugs, Bucimex chilensis and Primicimex cavernis, including updated phylogeny of Cimicidae

The Cimicidae is a family of blood-dependent ectoparasites in which dispersion capacity is greatly associated with host movements. Bats are the ancestral and most prevalent hosts for cimicids. Cimicids have a worldwide distribution matching that of their hosts, but the global classification is incomplete, especially for species outside the most common Cimicidae taxa. In this study, we place a little-studied cimicid species, Bucimex chilensis, within a comprehensive molecular phylogeny of Cimicidae by sequencing the genomic regions of this and other closely related species. For this study, we collected B. chilensis females from Myotis chiloensis in Tierra del Fuego, 1,300 km further south than previously known southernmost distribution boundary. We also sequenced COI regions from Primicimex cavernis, a species which together with B. chilensis comprise the entire subfamily Primiciminae. Using Bayesian posterior probability and maximum-likelihood approaches, we found that B. chilensis and P. cavernis clustered close to each other in the molecular analyses, receiving support from similar morphological features, agreeing with the morphology-based taxonomic placement of the two species within the subfamily Primiciminae. We also describe a previously unrecognized morphological adaptation of the tarsal structure, which allows the austral bat ectoparasite, B. chilensis, to cling on to the pelage of its known host, the Chilean myotis (Myotis chiloensis). Through a morphological study and behavioral observation, we elucidate how this tarsal structure operates, and we hypothesize that by clinging in the host pelage, B. chilensis is able to disperse effectively to new areas despite low host density. This is a unique feature shared by P. cavernis, the only other species in Primiciminae.

The Gambian epauletted fruit bat shows increased genetic divergence in the Ethiopian highlands and in an area of rapid urbanization

The Gambian epauletted fruit bat (Epomophorus gambianus) is an abundant species that roosts in both urban and rural settings. The possible role of E. gambianus as a reservoir host of zoonotic diseases underlines the need to better understand the species movement patterns. So far, neither observational nor phylogenetic studies have identified the dispersal range or behavior of this species. Comparative analyses of mitochondrial and nuclear markers from 20 localities across the known distribution of E. gambianus showed population panmixia, except for the populations in Ethiopia and southern Ghana (Accra and Ve-Golokwati). The Ethiopian population may be ancestral and is highly divergent to the species across the rest of its range, possibly reflecting isolation of an ancient colonization along an east-west axis. Mitochondrial haplotypes in the Accra population display a strong signature of a past bottleneck event; evidence of either an ancient or recent bottleneck using microsatellite data, however, was not detected. Demographic analyses identified population expansion in most of the colonies, except in the female line of descent in the Accra population. The molecular analyses of the colonies from Ethiopia and southern Ghana show gender dispersal bias, with the mitochondrial DNA fixation values over ten times those of the nuclear markers. These findings indicate free mixing of the species across great distances, which should inform future epidemiological studies.

Comparative phylogeography of parasitic Laelaps mites contribute new insights into the specialist-generalist variation hypothesis (SGVH)

Background

The specialist-generalist variation hypothesis (SGVH) in parasites suggests that, due to patchiness in habitat (host availability), specialist species will show more subdivided population structure when compared to generalist species. In addition, since specialist species are more prone to local stochastic extinction events with their hosts, they will show lower levels of intraspecific genetic diversity when compared to more generalist.

Results

To test the wider applicability of the SGVH we compared 337 cytochrome oxidase I mitochondrial DNA and 268 nuclear tropomyosin DNA sequenced fragments derived from two co-distributed Laelaps mite species and compared the data to 294 COI mtDNA sequences derived from the respective hosts Rhabdomys dilectus, R. bechuanae, Mastomys coucha and M. natalensis. In support of the SGVH, the generalist L. muricola was characterized by a high mtDNA haplotypic diversity of 0.97 (±0.00) and a low level of population differentiation (mtDNA F_st = 0.56, p < 0.05; nuDNA F_st = 0.33, P < 0.05) while the specialist L. giganteus was overall characterized by a lower haplotypic diversity of 0.77 (±0.03) and comparatively higher levels of population differentiation (mtDNA F_st = 0.87, P < 0.05; nuDNA F_st = 0.48, P < 0.05). When the two specialist L. giganteus lineages, which occur on two different Rhabdomys species, are respectively compared to the generalist parasite, L. muricola, the SGVH is not fully supported. One of the specialist L. giganteus species occurring on R. dilectus shows similar low levels of population differentiation (mtDNA F_st = 0.53, P < 0.05; nuDNA F_st = 0.12, P < 0.05) than that found for the generalist L. muricola. This finding can be correlated to differences in host dispersal: R. bechuanae populations are characterized by a differentiated mtDNA F_st of 0.79 (P < 0.05) while R. dilectus populations are less structured with a mtDNA F_st = 0.18 (P < 0.05).

Conclusions

These findings suggest that in ectoparasites, host specificity and the vagility of the host are both important drivers for parasite dispersal. It is proposed that the SGHV hypothesis should also incorporate reference to host dispersal since in our case only the specialist species who occur on less mobile hosts showed more subdivided population structure when compared to generalist species.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1245-7) contains supplementary material, which is available to authorized users.

Host association influences variation at salivary protein genes in the bat ectoparasite Cimex adjunctus

Parasite-host relationships create strong selection pressures that can lead to adaptation and increasing specialization of parasites to their hosts. Even in relatively loose host-parasite relationships, such as between generalist ectoparasites and their hosts, we may observe some degree of specialization of parasite populations to one of the multiple potential hosts. Salivary proteins are used by blood-feeding ectoparasites to prevent hemostasis in the host and maximize energy intake. We investigated the influence of association with specific host species on allele frequencies of salivary protein genes in Cimex adjunctus, a generalist blood-feeding ectoparasite of bats in North America. We analysed two salivary protein genes: an apyrase, which hydrolyses ATP at the feeding site and thus inhibits platelet aggregation, and a nitrophorin, which brings nitrous oxide to the feeding site, inhibiting platelet aggregation and vasoconstriction. We observed more variation at both salivary protein genes among parasite populations associated with different host species than among populations from different spatial locations associated with the same host species. The variation in salivary protein genes among populations on different host species was also greater than expected under a neutral scenario of genetic drift and gene flow. Finally, host species was an important predictor of allelic divergence in genotypes of individual C. adjunctus at both salivary protein genes. Our results suggest differing selection pressures on these two salivary protein genes in C. adjunctus depending on the host species.

Population structure in two geographically sympatric and congeneric ectoparasites (Cimex adjunctus and Cimex lectularius) in the North American Great Lakes region

Subdivided populations can be described by different models of population structure that reflect population organization, dynamics, and connectivity. We used genetic data to investigate population structure in two geographically sympatric, congeneric species of generalist ectoparasites of warm-blooded animals. We characterized the spatial genetic structure of the eastern bat bug (Cimex adjunctus Barber, 1939), an understudied and fairly abundant species, using microsatellite markers at a spatial scale representing contemporary dispersal of the species. We found seven genetic clusters, global [Formula: see text] of 0.2, 33% of genetic variation among sites, and nonsignificant isolation-by-distance. We compared these results with the common bed bug (Cimex lectularius L., 1758), a closely related but conversely well-known species, in the same geographic area. We found stronger genetic structuring in C. lectularius than in C. adjunctus, with 11 genetic clusters, [Formula: see text] of 0.7, 57% of genetic variation among sites, and significant but weak isolation-by-distance (R2 = 0.09). These results suggest that while both species can be described as having classic metapopulation structure, C. adjunctus leans more towards a patchy population and C. lectularius leans more towards a nonequilibrium metapopulation. The difference in population structure between these species may be attributable to differences in movement potential and extinction–colonization dynamics.

Comparative analysis of landscape effects on spatial genetic structure of the big brown bat and one of its cimicid ectoparasites

Identification of landscape features that correlate with genetic structure permits understanding of factors that may influence gene flow in a species. Comparing effects of the landscape on a parasite and host provides potential insights into parasite‐host ecology. We compared fine‐scale spatial genetic structure between big brown bats (Eptesicus fuscus) and their cimicid ectoparasite (Cimex adjunctus; class Insecta) in the lower Great Lakes region of the United States, in an area of about 160,000 km². We genotyped 142 big brown bat and 55 C. adjunctus samples at eight and seven microsatellite loci, respectively, and inferred effects of various types of land cover on the genetic structure of each species. We found significant associations between several land cover types and genetic distance in both species, although different land cover types were influential in each. Our results suggest that even in a parasite that is almost entirely reliant on its hosts for dispersal, land cover can affect gene flow differently than in the hosts, depending on key ecological aspects of both species.

Host association and selection on salivary protein genes in bed bugs and related blood-feeding ectoparasites

Reciprocal selective pressures can drive coevolutionary changes in parasites and hosts, and result in parasites that are highly specialized to their hosts. Selection and host co-adaptation are better understood in endoparasites than in ectoparasites, whose life cycles may be more loosely linked to that of their hosts. Blood-feeding ectoparasites use salivary proteins to prevent haemostasis in the host, and maximize energy intake. Here we looked for signals of selection in salivary protein genes of ectoparasite species from a single genus (Cimex) that associate with a range of hosts including mammals (bats and humans) and birds (swallows). We analysed two genes that code for salivary proteins that inhibit platelet aggregation and vasoconstriction and may directly affect the efficiency of blood feeding in these species. Significant positive selection was detected at five codons in one gene in all bat-associated species groups. Our results suggest association with bats, versus humans or swallows, has posed a selective pressure on the salivary apyrase gene in species of Cimex.