Comparative transcriptomics of anal fin pigmentation patterns in cichlid fishes.

BACKGROUND: Understanding the genetic basis of novel traits is a central topic in evolutionary biology. Two novel pigmentation phenotypes, egg-spots and blotches, emerged during the rapid diversification of East African cichlid fishes. Egg-spots are circular pigmentation markings on the anal fins of hundreds of derived haplochromine cichlids species, whereas blotches are patches of conspicuous anal fin pigmentation with ill-defined boundaries that occur in few species that belong to basal cichlid lineages. Both traits play an important role in the breeding behavior of this group of fishes. Knowledge about the origin, homology and underlying genetics of these pigmentation traits is sparse. RESULTS: Here, we present a comparative transcriptomic and differential gene expression analysis of egg-spots and blotches. We first conducted an RNA sequencing experiment where we compared egg-spot tissue with the remaining portion of egg-spot-free fin tissue using six individuals of Astatotilapia burtoni. We identified 1229 differentially expressed genes between the two tissue types. We then showed that rates of evolution of these genes are higher than average estimated on whole transcriptome data. Using quantitative real-time PCR, we found that 29 out of a subset of 46 differentially expressed genes showed an analogous expression pattern in another haplochromine species' egg-spots, Cynotilapia pulpican, strongly suggesting that these genes are involved in the egg-spot phenotype. Among these are the previously identified egg-spot gene fhl2a, two known patterning genes (hoxC12a and bmp3) as well as other pigmentation related genes such as asip. Finally, we analyzed the expression patterns of the same gene subset in two species that feature blotches instead of egg-spots, one haplochromine species (Pseudocrenilabrus philander) and one ectodine species (Callochromis macrops), revealing that the expression patterns in blotches and egg-spots are rather distinct. CONCLUSIONS: We identified several candidate genes that will serve as an important and useful resource for future research on the emergence and diversification of cichlid fishes' egg-spots. Only a limited degree of conservation of gene expression patterns was detected between the egg-spots of the derived haplochromines and blotches from ancestral haplochromines, as well as between the two types of blotches, suggesting an independent origin of these traits.

Population genetics and phylogeographic history of the insular lizard Podarcis lilfordi (Gunther, 1874) from the Balearic Islands based on genome-wide polymorphic data.

Islands provide a great system to explore the processes that maintain genetic diversity and promote local adaptation. We explored the genomic diversity of the Balearic lizard Podarcis lilfordi, an endemic species characterized by numerous small insular populations with large phenotypic diversity. Using the newly available genome for this species, we characterized more than 300,000 SNPs, merging genotyping-by-sequencing (GBS) data with previously published restriction site-associated DNA sequencing (RAD-Seq) data, providing a dataset of 16 island populations (191 individuals) across the range of species distribution (Menorca, Mallorca, and Cabrera). Results indicate that each islet hosts a well-differentiated population (F ST = 0.247 ± 0.09), with no recent immigration/translocation events. Contrary to expectations, most populations harbor a considerable genetic diversity (mean nucleotide diversity, P i = 0.144 ± 0.021), characterized by overall low inbreeding values (F IS < 0.1). While the genetic diversity significantly decreased with decreasing islet surface, maintenance of substantial genetic diversity even in tiny islets suggests variable selection or other mechanisms that buffer genetic drift. Maximum-likelihood tree based on concatenated SNP data confirmed the existence of the two major independent lineages of Menorca and Mallorca/Cabrera. Multiple lines of evidence, including admixture and root testing, robustly placed the origin of the species in the Mallorca Island, rather than in Menorca. Outlier analysis mainly retrieved a strong signature of genome differentiation between the two major archipelagos, especially in the sexual chromosome Z. A set of proteins were target of multiple outliers and primarily associated with binding and catalytic activity, providing interesting candidates for future selection studies. This study provides the framework to explore crucial aspects of the genetic basis of phenotypic divergence and insular adaptation.

Insular holobionts: persistence and seasonal plasticity of the Balearic wall lizard (Podarcis lilfordi) gut microbiota.

Background: Integrative studies of animals and associated microbial assemblages (i.e., the holobiont) are rapidly changing our perspectives on organismal ecology and evolution. Insular vertebrates provide ideal natural systems to understand patterns of host-gut microbiota coevolution, the resilience and plasticity these microbial communities over temporal and spatial scales, and ultimately their role in the host ecological adaptation. Methods: Here we used the endemic Balearic wall lizard Podarcis lilfordi to dissect the drivers of the microbial diversity within and across host allopatric populations/islets. By focusing on three extensively studied populations/islets of Mallorca (Spain) and fecal sampling from individually identified lizards along two years (both in spring and autumn), we sorted out the effect of islet, sex, life stage, year and season on the microbiota composition. We further related microbiota diversity to host genetics, trophic ecology and expected annual metabolic changes. Results: All the three populations showed a remarkable conservation of the major microbial taxonomic profile, while carrying their unique microbial signature at finer level of taxonomic resolution (Amplicon Sequence Variants (ASVs)). Microbiota distances across populations were compatible with both host genetics (based on microsatellites) and trophic niche distances (based on stable isotopes and fecal content). Within populations, a large proportion of ASVs (30-50%) were recurrently found along the four sampling dates. The microbial diversity was strongly marked by seasonality, with no sex effect and a marginal life stage and annual effect. The microbiota showed seasonal fluctuations along the two sampled years, primarily due to changes in the relative abundances of fermentative bacteria (mostly families Lachnospiraceae and Ruminococcaceae), without any major compositional turnover. Conclusions: These results support a large resilience of the major compositional aspects of the P. lilfordi gut microbiota over the short-term evolutionary divergence of their host allopatric populations (<10,000 years), but also indicate an undergoing process of parallel diversification of the both host and associated gut microbes. Predictable seasonal dynamics in microbiota diversity suggests a role of microbiota plasticity in the lizards' metabolic adaptation to their resource-constrained insular environments. Overall, our study supports the need for longitudinal and integrative studies of host and associated microbes in natural systems.

Chromosome-level genome assembly of Lilford's wall lizard, Podarcis lilfordi (Günther, 1874) from the Balearic Islands (Spain).

The Mediterranean lizard Podarcis lilfordi is an emblematic species of the Balearic Islands. The extensive phenotypic diversity among extant isolated populations makes the species a great insular model system for eco-evolutionary studies, as well as a challenging target for conservation management plans. Here we report the first high-quality chromosome-level assembly and annotation of the P. lilfordi genome, along with its mitogenome, based on a mixed sequencing strategy (10X Genomics linked reads, Oxford Nanopore Technologies long reads and Hi-C scaffolding) coupled with extensive transcriptomic data (Illumina and PacBio). The genome assembly (1.5 Gb) is highly contiguous (N50 = 90 Mb) and complete, with 99% of the sequence assigned to candidate chromosomal sequences and >97% gene completeness. We annotated a total of 25,663 protein-coding genes translating into 38,615 proteins. Comparison to the genome of the related species Podarcis muralis revealed substantial similarity in genome size, annotation metrics, repeat content, and a strong collinearity, despite their evolutionary distance (~18-20 MYA). This genome expands the repertoire of available reptilian genomes and will facilitate the exploration of the molecular and evolutionary processes underlying the extraordinary phenotypic diversity of this insular species, while providing a critical resource for conservation genomics.

Nuclear-mitochondrial sequences as witnesses of past interbreeding and population diversity in the jumping bristletail Mesomachilis.

Nuclear-mitochondrial sequences (NUMTs) can be especially powerful tools for evolutionary studies by providing a window into the nature of ancestral mitochondrial (mt) lineages. Here, we illustrate this property of NUMTs through a survey of such sequences in two species of the jumping bristletail genus Mesomachilis. An ∼1,800 basepair fragment encompassing three mt genes--COI, tRNA-Leu, and COII--was cloned twice for a sample of specimens, and colonies were screened. NUMTs in this genus are abundant and genetically diverse (up to 23% different from the associated mt sequence). Numerous independent nuclear integration events are scattered at different time depths, suggesting a continuous process of NUMT generation. By combining phylogenetic reconstruction with analyses of the pattern of nucleotide substitution on mt and NUMT branches, we inferred the dynamics of NUMT origin and evolution. The majority of NUMTs in the data set have apparently evolved from source mt lineages that were substantially different from the currently associated mt sequence, indicating NUMT transfer via matings among divergent lineages. The original source mt lineages either were unsampled or are extinct. Translocations of NUMTs to new nuclear genomes can preserve NUMTs following extinction of their ancestral mt lineages and can be used to detect mixing between divergent mt lineages and past levels of mt diversity.

Evolutionary Genetics of Microbial Symbioses.

Symbiosis is the living together of dissimilar organisms [...].

Accelerated microevolution in an outer membrane protein (OMP) of the intracellular bacteria Wolbachia.

BACKGROUND: Outer membrane proteins (OMPs) of Gram-negative bacteria are key players in the biology of bacterial-host interactions. However, while considerable attention has been given to OMPs of vertebrate pathogens, relatively little is known about the role of these proteins in bacteria that primarily infect invertebrates. One such OMP is found in the intracellular bacteria Wolbachia, which are widespread symbionts of arthropods and filarial nematodes. Recent experimental studies have shown that the Wolbachia surface protein (WSP) can trigger host immune responses and control cell death programming in humans, suggesting a key role of WSP for establishment and persistence of the symbiosis in arthropods. RESULTS: Here we performed an analysis of 515 unique alleles found in 831 Wolbachia isolates, to investigate WSP structure, microevolution and population genetics. WSP shows an eight-strand transmembrane beta-barrel structure with four extracellular loops containing hypervariable regions (HVRs). A clustering approach based upon patterns of HVR haplotype diversity was used to group similar WSP sequences and to estimate the relative contribution of mutation and recombination during early stages of protein divergence. Results indicate that although point mutations generate most of the new protein haplotypes, recombination is a predominant force triggering diversity since the very first steps of protein evolution, causing at least 50% of the total amino acid variation observed in recently diverged proteins. Analysis of synonymous variants indicates that individual WSP protein types are subject to a very rapid turnover and that HVRs can accommodate a virtually unlimited repertoire of peptides. Overall distribution of WSP across hosts supports a non-random association of WSP with the host genus, although extensive horizontal transfer has occurred also in recent times. CONCLUSIONS: In OMPs of vertebrate pathogens, large recombination impact, positive selection, reduced structural and compositional constraints, and extensive lateral gene transfer are considered hallmarks of evolution in response to the adaptive immune system. However, Wolbachia do not infect vertebrates. Here we predict that the rapid turnover of WSP loop motifs could aid in evading or inhibiting the invertebrate innate immune response. Overall, these features identify WSP as a strong candidate for future studies of host-Wolbachia interactions that affect establishment and persistence of this widespread endosymbiosis.

The mushroom habitat as an ecological arena for global exchange of Wolbachia.

Wolbachia infect a variety of arthropod and nematode hosts, but in arthropods, host phylogenetic relationships are usually poor predictors of strain similarity. This suggests that new infections are often established by horizontal transmission. To gain insight into the factors affecting the probability of horizontal transmission among host species, we ask how host phylogeny, geographical distribution and ecology affect patterns of Wolbachia strain similarity. We used multilocus sequence typing (MLST) to characterize Wolbachia strain similarity among dipteran hosts associated with fleshy mushrooms. Wolbachia Supergroup A was more common than Supergroup B in Diptera, and also more common in mycophagous than non-mycophagous Diptera. Within Supergroup A, host family within Diptera had no effect on strain similarity, and there was no tendency for Wolbachia strains from sympatric host species to be more similar to one another than to strains from hosts in different biogeographical realms. Supergroup A strains differed between mycophagous and non-mycophagous Diptera more than expected by chance, suggesting that ecological associations can facilitate horizontal transmission of Wolbachia within mycophagous fly communities. For Supergroup B, there were no significant associations between strain similarity and host phylogeny, biogeography, or ecology. We identified only two cases in which closely related hosts carried closely related Wolbachia strains, evidence that Wolbachia-host co-speciation or early introgression can occur but may not be a major contributor to overall strain diversity. Our results suggest that horizontal transmission of Wolbachia can be influenced by host ecology, thus leading to partial restriction of Wolbachia strains or strain groups to particular guilds of insects.

Microbial co-occurrence networks of gut microbiota reveal community conservation and diet-associated shifts in cichlid fishes.

BACKGROUND: The extent to which deterministic rather than stochastic processes guide gut bacteria co-existence and ultimately their assembling into a community remains largely unknown. Co-occurrence networks of bacterial associations offer a powerful approach to begin exploring gut microbial community structure, maintenance and dynamics, beyond compositional aspects alone. Here we used an iconic model system, the cichlid fishes, with their multiple lake assemblages and extraordinary ecological diversity, to investigate a) patterns of microbial associations that were robust to major phylogeographical variables, and b) changes in microbial network structure along dietary shifts. We tackled these objectives using the large gut microbiota sequencing dataset available (nine lakes from Africa and America), building geographical and diet-specific networks and performing comparative network analyses. RESULTS: Major findings indicated that lake and continental microbial networks were highly resembling in global topology and node taxonomic composition, despite the heterogeneity of the samples. A small fraction of the observed co-occurrences among operational taxonomic units (OTUs) was conserved across all lake assemblages. These were all positive associations and involved OTUs within the genera Cetobacterium and Turicibacter and several OTUs belonging to the families of Peptostreptococcaceae and Clostridiaceae (order Clostridiales). Mapping of diet contribution on the African Lake Tanganyika network (therefore excluding the geographic variable) revealed a clear community change from carnivores (C) to omnivores (O) to herbivores (H). Node abundances and effect size for pairwise comparisons between diets supported a strong contrasting pattern between C and H. Moreover, diet-associated nodes in H formed complex modules of positive interactions among taxonomically diverse bacteria (mostly Verrucomicrobia and Proteobacteria). CONCLUSIONS: Conservation of microbial network topologies and specific bacterial associations across distinct lake assemblages point to a major host-associated effect and potential deterministic processes shaping the cichlid gut microbiota. While the origin and biological relevance of these common associations remain unclear, their persistence suggests an important functional role in the cichlid gut. Among the very diverse cichlids of L. Tanganyika, diet nonetheless represents a major driver of microbial community changes. By intersecting results from predictive network inferences and experimental trials, future studies will be directed to explore the strength of these associations, predict the outcome of community alterations driven by diet and ultimately help understanding the role of gut microbiota in cichlid trophic diversification.

Phylogeography and Ecological Niche Shape the Cichlid Fish Gut Microbiota in Central American and African Lakes.

Cichlid fishes, with their repeated colonization of lakes and subsequent radiations at different scales of phylogenetic and ecological diversification, offer an excellent model system to understand the factors shaping the host-gut microbiota association in nature. Here, we characterized the gut microbiota of the Amphilophus species complex from Central America (known as the Midas cichlid complex), encompassing 158 wild specimens (13 species) collected from seven Nicaraguan lakes, and combined these data with previously published data from two African lakes (spanning 29 species). Our aim was to comprehensively explore trends in microbiota variation and persistence along the large spatial and temporal scales of cichlid diversification (from the oldest radiation in L. Tanganyika, 9-12 My old, to young ones in Nicaraguan crater lakes, <0.5 My old), in allopatry and sympatry (within and across lakes), and across the range of dietary niches (from highly specialized to generalist feeders). Despite their extraordinary diversity, cichlids shared a remarkably conserved microbial taxonomic profile, which argues for a primary role of the host genetics in the assembly and maintenance of these microbial communities. Within this partly constrained microbiota profile, geographic isolation (continent and lake) represented the first level of discrimination. For the Midas cichlid, a partial congruency was found between host microbiota and genetic distances, suggesting that microbial communities have partly diversified along their cichlid phylogeographic history of crater lake colonization. In sympatry (within lakes), the young and poorly ecologically diversified cichlid assemblages of Central American lakes display largely unresolved gut microbiotas (in terms of both alpha and beta diversities), whereas the phylogenetically and ecologically diverse species found in African lakes showed greater microbial interspecific diversity. This pattern largely points to the level of habitat segregation, trophic niche overlap, and reproductive barriers as major modulators of the gut microbiota connectivity among sympatric species.

Processes shaping gut microbiota diversity in allopatric populations of the endemic lizard Podarcis lilfordi from Menorcan islets (Balearic Islands).

Compositional variation of the gut microbiota across host allopatric populations can reflect both adaptation and stochasticity since the time of separation. Major factors shaping this variation include the host phylogeographic and demographic history, the microbiota inheritance, environmental inputs and dispersal of bacteria. Here we explored the impact of these factors in driving gut community diversity in seven allopatric populations of the omnivorous lizard Podarcis lilfordi from the Menorcan coastal islets, all descending from an ancestral mainland population. Using 16S rRNA Illumina sequencing, we showed that 'islet' and 'age' (time since islet separation from mainland) were the only significant variables in microbial community clustering, suggesting a partial islet-restricted diversification following these lizards phylogeography. Despite a significant variation, islets/populations were characterized by a remarkably low bacterial uniqueness (2.4% of total OTUs) and a minor differential enrichment of taxa, indicating a negligible impact of local inputs and important host common constraints. Overall, the extant pattern of similarity/dissimilarity among islets is compatible with partial retention of the ancestral mainland microbial pool, with differences among islets potentially explained by a differential loss of bacteria following population fragmentation and bottlenecks (i.e. ecological drift). While more quantitative data are needed to validate this hypothesis, this study unveils the importance of considering both neutral and niche-driven processes in driving contemporary patterns of gut metacommunity diversity.

New criteria for selecting the origin of DNA replication in Wolbachia and closely related bacteria.

BACKGROUND: The annotated genomes of two closely related strains of the intracellular bacterium Wolbachia pipientis have been reported without the identifications of the putative origin of replication (ori). Identifying the ori of these bacteria and related alpha-Proteobacteria as well as their patterns of sequence evolution will aid studies of cell replication and cell density, as well as the potential genetic manipulation of these widespread intracellular bacteria. RESULTS: Using features that have been previously experimentally verified in the alpha-Proteobacterium Caulobacter crescentus, the origin of DNA replication (ori) regions were identified in silico for Wolbachia strains and eleven other related bacteria belonging to Ehrlichia, Anaplasma, and Rickettsia genera. These features include DnaA-, CtrA- and IHF-binding sites as well as the flanking genes in C. crescentus. The Wolbachia ori boundary genes were found to be hemE and COG1253 protein (CBS domain protein). Comparisons of the putative ori region among related Wolbachia strains showed higher conservation of bases within binding sites. CONCLUSION: The sequences of the ori regions described here are only similar among closely related bacteria while fundamental characteristics like presence of DnaA and IHF binding sites as well as the boundary genes are more widely conserved. The relative paucity of CtrA binding sites in the ori regions, as well as the absence of key enzymes associated with DNA replication in the respective genomes, suggest that several of these obligate intracellular bacteria may have altered replication mechanisms. Based on these analyses, criteria are set forth for identifying the ori region in genome sequencing projects.

Convergence of gut microbiotas in the adaptive radiations of African cichlid fishes.

Ecoevolutionary dynamics of the gut microbiota at the macroscale level, that is, in across-species comparisons, are largely driven by ecological variables and host genotype. The repeated explosive radiations of African cichlid fishes in distinct lakes, following a dietary diversification in a context of reduced genetic diversity, provide a natural setup to explore convergence, divergence and repeatability in patterns of microbiota dynamics as a function of the host diet, phylogeny and environment. Here we characterized by 16S rRNA amplicon sequencing the gut microbiota of 29 cichlid species from two distinct lakes/radiations (Tanganyika and Barombi Mbo) and across a broad dietary and phylogenetic range. Within each lake, a significant deviation between a carnivorous and herbivorous lifestyle was found. Herbivore species were characterized by an increased bacterial taxonomic and functional diversity and converged in key compositional and functional community aspects. Despite a significant lake effect on the microbiota structure, this process has occurred with remarkable parallels in the two lakes. A metabolic signature most likely explains this trend, as indicated by a significant enrichment in herbivores/omnivores of bacterial taxa and functions associated with fiber degradation and detoxification of plant chemical compounds. Overall, compositional and functional aspects of the gut microbiota individually and altogether validate and predict main cichlid dietary habits, suggesting a fundamental role of gut bacteria in cichlid niche expansion and adaptation.

Gut Microbiota Dynamics during Dietary Shift in Eastern African Cichlid Fishes.

The gut microbiota structure reflects both a host phylogenetic history and a signature of adaptation to the host ecological, mainly trophic niches. African cichlid fishes, with their array of closely related species that underwent a rapid dietary niche radiation, offer a particularly interesting system to explore the relative contribution of these two factors in nature. Here we surveyed the host intra- and interspecific natural variation of the gut microbiota of five cichlid species from the monophyletic tribe Perissodini of lake Tanganyika, whose members transitioned from being zooplanktivorous to feeding primarily on fish scales. The outgroup riverine species Astatotilapia burtoni, largely omnivorous, was also included in the study. Fusobacteria, Firmicutes and Proteobacteria represented the dominant components in the gut microbiota of all 30 specimens analysed according to two distinct 16S rRNA markers. All members of the Perissodini tribe showed a homogenous pattern of microbial alpha and beta diversities, with no significant qualitative differences, despite changes in diet. The recent diet shift between zooplantkon- and scale-eaters simply reflects on a significant enrichment of Clostridium taxa in scale-eaters where they might be involved in the scale metabolism. Comparison with the omnivorous species A. burtoni suggests that, with increased host phylogenetic distance and/or increasing herbivory, the gut microbiota begins differentiating also at qualitative level. The cichlids show presence of a large conserved core of taxa and a small set of core OTUs (average 13-15%), remarkably stable also in captivity, and putatively favoured by both restricted microbial transmission among related hosts (putatively enhanced by mouthbrooding behavior) and common host constraints. This study sets the basis for a future large-scale investigation of the gut microbiota of cichlids and its adaptation in the process of the host adaptive radiation.

Comparative transcriptomics of Eastern African cichlid fishes shows signs of positive selection and a large contribution of untranslated regions to genetic diversity.

The hundreds of endemic species of cichlid fishes in the East African Great Lakes Tanganyika, Malawi, and Victoria are a prime model system in evolutionary biology. With five genomes currently being sequenced, eastern African cichlids also represent a forthcoming genomic model for evolutionary studies of genotype-to-phenotype processes in adaptive radiations. Here we report the functional annotation and comparative analyses of transcriptome data sets for two eastern African cichlid species, Astatotilapia burtoni and Ophthalmotilapia ventralis, representatives of the modern haplochromines and ectodines, respectively. Nearly 647,000 expressed sequence tags were assembled in more than 46,000 contigs for each species using the 454 sequencing technology, largely expanding the current sequence data set publicly available for these cichlids. Total predicted coverage of their proteome diversity is approximately 50% for both species. Comparative qualitative and quantitative analyses show very similar transcriptome data for the two species in terms of both functional annotation and relative abundance of gene ontology terms expressed. Average genetic distance between species is 1.75% when all transcript types are considered including nonannotated sequences, 1.33% for annotated sequences only including untranslated regions, and decreases to nearly half, 0.95%, for coding sequences only, suggesting a large contribution of noncoding regions to their genetic diversity. Comparative analyses across the two species, tilapia and the outgroup medaka based on an overlapping data set of 1,216 genes (∼526 kb) demonstrate cichlid-specific signature of disruptive selection and provide a set of candidate genes that are putatively under positive selection. Overall, these data sets offer the genetic platform for future comparative analyses in light of the upcoming genomes for this taxonomic group.

Modes of acquisition of Wolbachia: horizontal transfer, hybrid introgression, and codivergence in the Nasonia species complex.

Wolbachia are maternally inherited bacteria that infect a large number of insects and are responsible for different reproductive alterations of their hosts. One of the key features of Wolbachia biology is its ability to move within and between host species, which contributes to the impressive diversity and range of infected hosts. Using multiple Wolbachia genes, including five developed for Multi-Locus Sequence Typing (MLST), the diversity and modes of movement of Wolbachia within the wasp genus Nasonia were investigated. Eleven different Wolbachia were found in the four species of Nasonia, including five newly identified infections. Five infections were acquired by horizontal transmission from other insect taxa, three have been acquired by hybridization between two Nasonia species, which resulted in a mitochondrial-Wolbachia sweep from one species to the other, and at least three have codiverged during speciation of their hosts. The results show that a variety of transfer mechanisms of Wolbachia are possible even within a single host genus. Codivergence of Wolbachia and their hosts is uncommon and provides a rare opportunity to investigate long-term Wolbachia evolution within a host lineage. Using synonymous divergence among codiverging infections and host nuclear genes, we estimate Wolbachia mutation rates to be approximately one-third that of the nuclear genome.

Specialization and geographic isolation among Wolbachia symbionts from ants and lycaenid butterflies.

Wolbachia are the most prevalent and influential bacteria described among the insects to date. But despite their significance, we lack an understanding of their evolutionary histories. To describe the evolution of symbioses between Wolbachia and their hosts, we surveyed global collections of two diverse families of insects, the ants and lycaenid butterflies. In total, 54 Wolbachia isolates were typed using a Multi Locus Sequence Typing (MLST) approach, in which five unlinked loci were sequenced and analyzed to decipher evolutionary patterns. AMOVA and phylogenetic analyses demonstrated that related Wolbachia commonly infect related hosts, revealing a pattern of host association that was strongest among strains from the ants. A review of the literature indicated that horizontal transfer is most successful when Wolbachia move between related hosts, suggesting that patterns of host association are driven by specialization on a common physiological background. Aside from providing the broadest and strongest evidence to date for Wolbachia specialization, our findings also reveal that strains from New World ants differ markedly from those in ants from other locations. We, therefore, conclude that both geographic and phylogenetic barriers have promoted evolutionary divergence among these influential symbionts.

Wolbachia: master manipulators of invertebrate biology.

Wolbachia are common intracellular bacteria that are found in arthropods and nematodes. These alphaproteobacteria endosymbionts are transmitted vertically through host eggs and alter host biology in diverse ways, including the induction of reproductive manipulations, such as feminization, parthenogenesis, male killing and sperm-egg incompatibility. They can also move horizontally across species boundaries, resulting in a widespread and global distribution in diverse invertebrate hosts. Here, we review the basic biology of Wolbachia, with emphasis on recent advances in our understanding of these fascinating endosymbionts.

Insight into the routes of Wolbachia invasion: high levels of horizontal transfer in the spider genus Agelenopsis revealed by Wolbachia strain and mitochondrial DNA diversity.

The pandemic distribution of Wolbachia (alpha-proteobacteria) across arthropods is largely due to the ability of these maternally inherited endosymbionts to successfully shift hosts across species boundaries. Yet it remains unclear whether Wolbachia has preferential routes of transfer among species. Here, we examined populations of eight species of the North American funnel-web spider genus Agelenopsis to evaluate whether Wolbachia show evidence for host specificity and the relative contribution of horizontal vs. vertical transmission of strains within and among related host species. Wolbachia strains were characterized by multilocus sequence typing (MLST) and Wolbachia surface protein (WSP) sequences, and analysed in relation to host phylogeny, mitochondrial diversity and geographical range. Results indicate that at least three sets of divergent Wolbachia strains invaded the genus Agelenopsis. After each invasion, the Wolbachia strains preferentially shuffled across species of this host genus by horizontal transfer rather than cospeciation. Decoupling of Wolbachia and host mitochondrial haplotype (mitotypes) evolutionary histories within single species reveals an extensive contribution of horizontal transfer also in the rapid dispersal of Wolbachia among conspecific host populations. These findings provide some of the strongest evidence to support the association of related Wolbachia strains with related hosts by means of both vertical and horizontal strain transmission. Similar analyses across a broader range of invertebrate taxa are needed, using sensitive methods for strain typing such as MLST, to determine if this pattern of Wolbachia dispersal is peculiar to Agelenopsis (or spiders), or is in fact a general pattern in arthropods.

Revisiting Wolbachia supergroup typing based on WSP: spurious lineages and discordance with MLST.

The obligate intracellular bacteria Wolbachia are taxonomically subdivided into eight supergroups (named A-H). Supergroup typing of strains has been mostly based on phylogenetic inference of the Wolbachia surface protein (wsp), a gene that recently has been shown to experience high rates of recombination. This brings into question its suitability not only for microtaxonomy, but also for supergroup classification of the genus. A Multilocus Sequence Typing (MLST) scheme for Wolbachia has recently been developed that types strains based on five conserved genes, thus providing a rigorous supergroup annotation of strains. Here we report striking discrepancies in supergroup designation between MLST and wsp inferences, and propose a revision of current methods for Wolbachia supergroup typing. Transfer of whole wsp gene sequences between supergroups A and B has occurred. Furthermore, as a result of intragenic recombination, wsp phylogeny creates spurious basal lineages that are not supported by MLST. For example, the proposed supergroup G, based upon wsp alone, likely represents only a wsp recombinant clade. Removal of supergroup G is advised until and unless the existence of this lineage is substantiated by other sequence information (e.g., MLST). We recommend a full characterization MLST for a correct strain typing, while, based on the current data set, use of a single MLST gene can be effective for supergroup designation of A and B strains. Finally, we note that the sharing of wsp sequences between A and B strains indicates a strong genetic cohesiveness of Wolbachia strains, supporting designation of these bacteria within the same species, W. pipientis.