Bacterial community composition in blood-sucking arthropods can shift dramatically across time and space. We used 16S rRNA gene amplification and pyrosequencing to investigate the relative impact of vertebrate host-related, arthropod-related and environmental factors on bacterial community composition in fleas and ticks collected from rodents in southern Indiana (USA). Bacterial community composition was largely affected by arthropod identity, but not by the rodent host or environmental conditions. Specifically, the arthropod group (fleas vs ticks) determined the community composition of bacteria, where bacterial communities of ticks were less diverse and more dependent on arthropod traits—especially tick species and life stage—than bacterial communities of fleas. Our data suggest that both arthropod life histories and the presence of arthropod-specific endosymbionts may mask the effects of the vertebrate host and its environment.

Mutations frequently reoccur in the human mitochondrial DNA (mtDNA). However, it is unclear whether recurrent mtDNA nodal mutations (RNMs), that is, recurrent mutations in stems of unrelated phylogenetic nodes, are functional and hence selectively constrained. To answer this question, we performed comprehensive parsimony and maximum likelihood analyses of 9,868 publicly available whole human mtDNAs revealing 1,606 single nodal mutations (SNMs) and 679 RNMs. We then evaluated the potential functionality of synonymous, nonsynonymous and RNA SNMs and RNMs. For synonymous mutations, we have implemented the Codon Adaptation Index. For nonsynonymous mutations, we assessed evolutionary conservation, and employed previously described pathogenicity score assessment tools. For RNA genes’ mutations, we designed a bioinformatic tool which compiled evolutionary conservation and potential effect on RNA structure. While comparing the functionality scores of nonsynonymous and RNA SNMs and RNMs with those of disease-causing mtDNA mutations, we found significant difference (P < 0.001). However, 24 RNMs and 67 SNMs had comparable values with disease-causing mutations reflecting their potential function thus being the best candidates to participate in adaptive events of unrelated lineages. Strikingly, some functional RNMs occurred in unrelated mtDNA lineages that independently altered susceptibility to the same diseases, thus suggesting common functionality. To our knowledge, this is the most comprehensive analysis of selective signatures in the mtDNA not only within proteins but also within RNA genes. For the first time, we discover virtually all positively selected RNMs in our phylogeny while emphasizing their dual role in past evolutionary events and in disease today.