# The distribution of wolbachia in cubitermes (termitidae, termitinae) castes and colonies: a modelling approach

Virginie Roy, Marc Girondot, Myriam Harry

► To cite this version:

Virginie Roy, Marc Girondot, Myriam Harry. The distribution of wolbachia in cubitermes (termitidae, termitinae) castes and colonies: a modelling approach. PLoS ONE, 2015, 10 (2), e1371/journal.pone.0116070. This version is distributed under a Creative Commons Attribution 4.0 license.

HAL Id: hal-02636192

https://hal.inrae.fr/hal-02636192v1

Submitted on 27 May 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons CC BY 4.0 - Attribution - International License
<figure><img src="image_2.png" /></figure>

RESEARCH ARTICLE

# The Distribution of *Wolbachia* in *Cubitermes* (Termitidae, Termitinae) Castes and Colonies: A Modelling Approach

Virginie Roy¹, Marc Girondot², Myriam Harry³*

1 iEES—Institut d'écologie et des sciences de l'environnement de Paris, Département SOLéO, Université Paris-Est Créteil, Faculté des Sciences et Technologie, 61 avenue du Général de Gaulle, 94010 Créteil cedex, France, 2 Laboratoire d'Ecologie, Systématique et Evolution, Département d'Ecologie des Populations et des Communautés, Université Paris-Sud 11, Bâtiment 362, 91405 Orsay Cedex, France, 3 Laboratoire Evolution, Génomes et Spéciation, UPR 9034 CNRS, UR 072 IRD, Université Paris Sud-11, avenue de la Terrasse, Bâtiment 13, 91198 Gif sur Yvette, France/UFR de Sciences, Université Paris-Sud 11, 91400 Orsay, France

* myriam.harry@u-psud.fr

## Abstract

*Wolbachia* are endosymbiotic bacteria of arthropods and nematodes that are able to manipulate host reproduction. Although vertically transmitted via the cytoplasm in eggs, horizontal transmission of *Wolbachia* among and within arthropod species has been shown to be common. Eusocial insects represent interesting models for studying *Wolbachia* transmission due to colonial organization and close interaction between nestmates. Here we conducted a detailed screening of *Wolbachia* infection for 15 colonies of the very common soil-feeding termites *Cubitermes* spp. *affinis subarquatus* (Termitidae, Termitinae) that consist of four distinct phylogenetic species in the Lopé forest Reserve, Gabon. Infection tests showed that 50% of the individuals were *Wolbachia* positive (N = 555) with 90% of reproductives and 48% of offspring infected. White soldiers, which are transitional stages preceding mature soldiers, had a significantly higher mean infection rate (74%) than the other castes and stages (63%, 33% and 39% for larvae, workers and mature soldiers, respectively). We used a maximum likelihood method and Akaike's Information Criterion in order to explain the non-expected high rate of *Wolbachia* infection in white soldiers. The best model included a species effect for the stochastic loss of *Wolbachia* and a caste effect for the rate of gain. After fitting, the best model selected was for a species-specific rate of loss with a null rate of new gain for larvae, workers and soldiers and a probability of 0.72 whatever the species, that a white soldier becomes newly contaminated during that stage. The mean expected infection rate in white soldiers without a new gain was estimated to 17% instead of the 74% observed. Here we discuss the possible explanations to the high infection rate observed in white soldiers.

<figure><img src="image_3.png" /></figure>

OPEN ACCESS

**Citation:** Roy V, Girondot M, Harry M (2015) The Distribution of *Wolbachia* in *Cubitermes* (Termitidae, Termitinae) Castes and Colonies: A Modelling Approach. PLoS ONE 10(2): e0116070. doi:10.1371/journal.pone.0116070

**Academic Editor:** Kostas Bourtzis, International Atomic Energy Agency, AUSTRIA

**Received:** June 5, 2014

**Accepted:** December 5, 2014

**Published:** February 11, 2015

**Copyright:** © 2015 Roy et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Data Availability Statement:** All relevant data are within the paper.

**Funding:** MH and VR thank the University of Paris-Est Créteil for the financing of the ABI Prism automated sequencer (Bonus Qualité Recherche Grant and recurring budgets). The funder had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

**Competing Interests:** The authors have declared that no competing interests exist.
**Introduction**

*Wolbachia* is an endosymbiotic bacterium that infects a wide range of arthropod and nematode species [1–3]. This parasite has been shown to increase its transmission via manipulation of reproduction and sex ratio of arthropod hosts using different mechanisms: cytoplasmic incompatibility, parthenogenesis induction, male killing and feminization of genetic males [4–8]. *Wolbachia* infections were initially thought to be restricted to the germ-line tissue of their hosts, but results in a variety of insects and isopods demonstrated that they have a wide tissue tropism [9,10] and are detected in various somatic tissues comprising muscles, gut, head and salivary glands [5,11–13]. Although vertically transmitted via the cytoplasm in eggs (i.e. maternally inherited), the lack of congruence observed between the phylogenetic tree of *Wolbachia* and that of their hosts [3], and observations of natural interspecific and intraspecific transfer of the bacteria [14,15], suggested a degree of horizontal transmission both between and within arthropod host species.

Due to their colonial organization with permanent contact and trophic interactions between nestmates, eusocial insects are particularly interesting for studying *Wolbachia* transmission. *Wolbachia* prevalence in eusocial bees and wasps has received growing interest in recent years [16–19]. The bacterium appears to be widespread in ants with about one third of worldwide ant species being infected [2, 20–22]. Studies investigating infection rates in ant populations and colonies have provided interesting lines of research concerning dynamics of *Wolbachia* [20, 23–30]. In termite species, *Wolbachia* was described with a prevalence of 27% [31] and a great strain diversity was identified in this group [31–35]. However, to our knowledge no information is available concerning the infection status within termite colonies.

Termites (Blattodea: Termitoidae), albeit also eusocial insects, differ in many ways from eusocial Hymenoptera. In particular, most termite species possess a differentiated soldier caste and numerous molts throughout the lifecycle of individuals (i.e. termites are heterometabolous insects). Due to their non-cryptic colonial organization in epigeous nests, termites of the family Termitidae, subfamily Termitinae, are ideal models for studying *Wolbachia* distribution and exchanges. Colonies are headed by a royal pair, called the primary king and queen. Larval stages differentiate into two lines, a sexual line which produces nymphs i.e. future winged imagoes (during a short seasonal period), and a sterile line which produces true workers. A subset of the workers (males and/or females depending on the species, typically females for *Cubitermes* species) further differentiate into pre-soldiers (white soldiers), that constitute a transitional stage preceding that of mature soldier (Fig. 1a). In soil-feeding Termitinae, only the workers feed on organic matter and foster other members of the colony, which are unable to feed by themselves, by trophallaxis. Particularly, young larvae and white soldiers receive pure saliva, while mature soldiers receive regurgitated (stomodeal) food [36], as proctodeal trophallaxis was never reported in Termitidae [37,38].

Here we provide a detailed survey of *Wolbachia* infection in 15 termite colonies. Primary queens and kings, larvae, workers, white soldiers and mature soldiers were sampled in colonies of the very common soil-feeding termites *Cubitermes* spp. *affinis subarquatus* (Termitidae, Termitinae) that consist of four distinct phylogenetic species in the Lopé forest Reserve (Middle Ogooué, Gabon) [39]. DNA from individuals and pools of salivary glands were screened for the presence of *Wolbachia* using PCR assays targeting the *wsp* gene. *Wolbachia* sequences from *Cubitermes* spp. *affinis subarquatus* were previously obtained for *wsp* and four other genes: *ftsz*, *fbpA*, *coxA* and 16S rDNA [33]. Because experimental studies could not be conducted for soil-feeding termites owing to their complex ethology, ecology and life cycle, we modeled our data in order to explain the unexpectedly high rate of *Wolbachia* infection in some sterile adult stages such as white soldiers.
<figure><img src="image_4.png" /><figcaption>Fig 1. a) Development of non-sexual lines in the genus *Cubitermes* synthetized from various reports about *Cubitermes* species life cycle [36,79-81]. Larval stages comprise males and females, which both differentiate into a single stage of male and female workers (Worker stage I). A subset of the female workers further differentiate into pre-soldiers (i.e. white soldiers, transitional stage preceding mature soldiers), b) Left: fitted parameters for the selected model: probability of stochastic loss of *Wolbachia* between two stages according to the species and right: probability of gain of *Wolbachia* during stages. c) Observed (bold lines) and expected under the model (dotted lines) infection rates for the four species *Cubitermes* sp. A, B, C and D and stages/castes.</figcaption></figure>

doi:10.1371/journal.pone.0116070.g001

## Material and Methods

### Termite and tissue samples

Termites were sampled in 15 colonies of *Cubitermes* spp. A, B, C and D *affinis subarquatus* (Termitidae, Termitinae) in the Lopé forest Reserve (Middle Ogooué, Gabon) [39]. The authors are grateful to the authorities of Gabon represented by Mr E. Mamfoumbi of the ‘Ministère des Eaux et Forêts’ for permission to work in the Lopé Reserve. Termite samples were preserved in absolute ethanol at 4°C until DNA extraction. We tested 555 *Cubitermes* individuals for the presence of *Wolbachia*: 535 from neutral castes/stages (larvae, workers, white soldiers and mature soldiers) and all primary reproductives present at the time of the sampling i.e. 20 dealate queens and kings (Table 1). For DNA extraction of these samples, abdomens were removed to avoid contamination / PCR inhibition from gut content. Different tissues throughout the body were screened for two individuals per caste: head (without salivary glands), eggs,
<table><caption>Table 1. Prevalence of Wolbachia in 15 Cubitermes colonies and infection rates across different castes and stages.</caption><thead><tr><th>Termite species</th><th>Colony</th><th>Wolbachiasupergroup</th><th>AN wsp</th><th>infected queens</th><th>N<sub>O</sub></th><th>infected kings</th><th>N<sub>K</sub></th><th>infected larvae</th><th>N<sub>L</sub></th><th>infected workers</th><th>N<sub>W</sub></th><th>infected white soldiers</th><th>N<sub>WS</sub></th><th>infected mature soldiers</th><th>N<sub>S</sub></th><th>N<sub>tot</sub> infected N<sub>total</sub> (%)</th></tr></thead><tbody><tr><td>Cubitermes sp. A</td><td>T7</td><td>B</td><td>DQ124660</td><td>100%</td><td>1</td><td></td><td>83%</td><td>6</td><td>70%</td><td>10</td><td>50%</td><td>4</td><td>100%</td><td>6</td><td>21/27 (77.77%)</td></tr><tr><td>Cubitermes sp. A</td><td>T16</td><td>B</td><td>DQ124659</td><td>100%</td><td>1</td><td></td><td>87%</td><td>15</td><td>90%</td><td>10</td><td>100%</td><td>1</td><td>40%</td><td>10</td><td>28/37 (75.76%)</td></tr><tr><td>Cubitermes sp. A</td><td>T42</td><td>B</td><td>DQ124661</td><td>100%</td><td>1</td><td>100%</td><td>1</td><td>48%</td><td>21</td><td>40%</td><td>10</td><td>75%</td><td>4</td><td>80%</td><td>20</td><td>35/57 (61.40%)</td></tr><tr><td>Cubitermes sp. B</td><td>T5</td><td>A</td><td>DQ124665</td><td>100%</td><td>1</td><td>0%</td><td>1</td><td>0%</td><td>15</td><td>80%</td><td>5</td><td>0%</td><td>3</td><td>5/25 (20.00%)</td><td></td><td></td></tr><tr><td>Cubitermes sp. B</td><td>T26</td><td>A</td><td>DQ124667</td><td>100%</td><td>1</td><td></td><td>0%</td><td>9</td><td>75%</td><td>15</td><td>87%</td><td>4</td><td>10%</td><td>10</td><td>5/24 (20.83%)</td></tr><tr><td>Cubitermes sp. B</td><td>T34</td><td>A</td><td>DQ124666</td><td>100%</td><td>1</td><td></td><td>0%</td><td>19</td><td>80%</td><td>15</td><td>87%</td><td>5</td><td>7%</td><td>15</td><td>39/50 (78.00%)</td></tr><tr><td>Cubitermes sp. B</td><td>T37</td><td>A</td><td>DQ124663</td><td>100%</td><td>1</td><td>100%</td><td>1</td><td>0%</td><td>20</td><td>20%</td><td>5</td><td>7%</td><td>15</td><td>4/42 (9.52%)</td><td></td><td></td></tr><tr><td>Cubitermes sp. B</td><td>T38</td><td>A</td><td>DQ124664</td><td>100%</td><td>1</td><td>100%</td><td>1</td><td>21%</td><td>19</td><td>25%</td><td>12</td><td>100%</td><td>1</td><td>50%</td><td>14</td><td>17/48 (35.42%)</td></tr><tr><td>Cubitermes sp. B</td><td>T17</td><td>A</td><td>DQ124655</td><td>100%</td><td>1</td><td></td><td>88%</td><td>16</td><td>100%</td><td>1</td><td>0%</td><td>1</td><td>0%</td><td>1</td><td>16/19 (84.21%)</td></tr><tr><td>Cubitermes sp. C</td><td>T24</td><td>A</td><td>DQ124654</td><td>100%</td><td>1</td><td>100%</td><td>1</td><td>64%</td><td>25</td><td>30%</td><td>10</td><td>100%</td><td>2</td><td>30%</td><td>10</td><td>26/49 (53.06%)</td></tr><tr><td>Cubitermes sp. C</td><td>T31</td><td>A (wsp)</td><td>FJ851084</td><td></td><td></td><td></td><td>90%</td><td>10</td><td>11%</td><td>9</td><td>70%</td><td>30</td><td>17%</td><td>6</td><td>32/55 (58.18%)</td></tr><tr><td>Cubitermes sp. C</td><td>T45</td><td>A (wsp)</td><td>FJ851085</td><td></td><td></td><td></td><td>70%</td><td>10</td><td>0%</td><td>3</td><td>100%</td><td>2</td><td>0%</td><td>5</td><td>9/20 (45.00%)</td></tr><tr><td>Cubitermes sp. C</td><td>T46</td><td>A</td><td>DQ124658</td><td>100%</td><td>1</td><td></td><td>56%</td><td>9</td><td>0%</td><td>10</td><td>67%</td><td>3</td><td>0%</td><td>10</td><td>8/33 (24.24%)</td></tr><tr><td>Cubitermes sp. C</td><td>TX</td><td>A</td><td>DQ124656</td><td>100%</td><td>1</td><td>0%</td><td>1</td><td>80%</td><td>10</td><td>10%</td><td>10</td><td>90%</td><td>10</td><td>0%</td><td>10</td><td>19/42 (45.24%)</td></tr><tr><td>Cubitermes sp. D</td><td>TD1</td><td>B</td><td>DQ124662</td><td>100%</td><td>1</td><td>100%</td><td>1</td><td>10%</td><td>10</td><td>100%</td><td>5</td><td>40%</td><td>10</td><td>12/27 (44.44%)</td><td></td><td></td></tr><tr><td>N total</td><td></td><td></td><td></td><td>100%</td><td></td><td></td><td>71.43%</td><td>62.5% (54–71)</td><td>32.54% (25–40)</td><td>74.03% (64–84)</td><td>38.62% (31–47)</td><td></td><td></td><td></td><td>145 (49.73%)</td></tr></tbody></table>

Wolbachia supergroups were determined in Roy and Harry [33] for all colonies except T31 and T45. Colonies T31 and T45 were not submitted to a MLST analysis and only the wsp gene was sequenced. AN: Genbank Accession Numbers. N: number of individuals tested for each caste or stage. 95% confidence intervals are indicated in brackets. Tests for the significance of infection were performed via the Wald statistics using a generalized linear model (GLZ). Results obtained between castes are summarized using letters; infection rates that do not significantly differ share the same letter.

doi:10.1371/journal.pone.0116070.t001
ovarian tissues and Malpighian tubules were dissected from queens, head and fat tissues from kings, head, thorax and legs from workers, white soldiers and mature soldiers. Two pools of salivary glands from ten workers were also dissected (Table 2).

## DNA extraction

Total genomic DNA was individually isolated using Wilson buffer (Tris-hydrochloride 1M, EDTA 0.5M, sodium chloride 4.5M, sodium dodecylsulfate 20%, dithiothreitol, proteinase K), followed by a salting-out procedure. DNA from tissue samples and salivary glands was extracted using a 10% Chelex solution, which has been described as a means of increasing the signal from the PCR amplification of small amounts of DNA [40]. In order to discard *Wolbachia* false negatives, PCRs targeted on termite DNA were performed using forward primer SR-J-14233 (5'-AAG AGC GAC GGG CGA TGT GT-3') and reverse primer SR-N-14588 (5'-AAA CTA GGA TTA GAT ACC CTA TTA T-3') amplifying a fragment of 12S ribosomal RNA gene [41]. Low-quality DNA samples that did not amplify termite rDNA were rejected.

## PCR detection of *Wolbachia*

DNA from termites and tissues was screened for the presence of *Wolbachia* via Polymerase Chain Reaction (PCR) using wsp primers: wsp81F (5'-TGG TCC AAT AAG TGA TGA AGA AAC-3') and wsp691R (5'-AAA AAT TAA ACG CTA CTC CA-3') [42]. Each PCR reaction was performed in a total volume of 12.5µL, composed of 2X Taq PCR Master Mix (Qiagen, France), 0.5µM of each primer, 2.35mM MgCl2, 0.016mg/mL Bovin Serum Albumine and 2µL of template DNA. The conditions for amplification were initial denaturation at 95°C for 3 minutes then 35 cycles of (1) denaturation at 95°C for 30 seconds, (2) annealing at 45°C for 30 seconds, (3) extension at 72°C for 1.30 minutes and a final extension at 72°C for 10 minutes. A

<table><caption>Table 2. Tissue distribution of *Wolbachia* in Cubitermes castes and stages.</caption><thead><tr><th>Castes</th><th>Tissues</th><th>12S rRNA gene</th><th>Wolbachia</th></tr></thead><tbody><tr><td rowspan="4">Queen (n = 2)</td><td>Head</td><td>+</td><td>+</td></tr><tr><td>Eggs</td><td>+</td><td>+</td></tr><tr><td>Ovarian tissues</td><td>+</td><td>+</td></tr><tr><td>Malpighian tubes</td><td>+</td><td>+</td></tr><tr><td rowspan="2">King (n = 2)</td><td>Head</td><td>+</td><td>+</td></tr><tr><td>Fat tissues</td><td>+</td><td>+</td></tr><tr><td rowspan="2">Workers (n = 2)</td><td>Head</td><td>+</td><td>+</td></tr><tr><td>Thorax + legs</td><td>+</td><td>+</td></tr><tr><td rowspan="2">White soldiers (n = 2)</td><td>Head</td><td>+</td><td>+</td></tr><tr><td>Thorax + legs</td><td>+</td><td>+</td></tr><tr><td rowspan="2">Mature soldiers (n = 2)</td><td>Head</td><td>+</td><td>-</td></tr><tr><td>Thorax + legs</td><td>+</td><td>+</td></tr><tr><td>Pool of 10 workers (n = 2)</td><td>Salivary glands</td><td>+</td><td>+</td></tr></tbody></table>

Queens were from T5 (C. spA) and T42 (C. spB), kings were from T38 (C. spA) and T42 (C. spB), workers, white soldiers and mature soldiers were from T42 (C. spB) and workers' salivary glands were from T34 (C. spA) and T42 (C. spB). The occurrence of *Wolbachia* was tested using wsp PCR and recorded as positive (+) or negative (-). The DNA quality was tested by using 12S ribosomal RNA gene.

doi:10.1371/journal.pone.0116070.t002
negative (sterile water) and a positive control (known infected queen DNA) were included in each PCR run. Each sample was tested twice and was considered to be infected when at least one of the PCRs was positive. Sequencing of wsp gene was realized for up to three sequences per colony and for salivary glands using BigDye Terminator Cycle Sequencing kit version 1.1 (Applied Biosystems). Sequences were obtained using an automatic DNA sequencer (Applied-Biosystems, ABI PRISM 310). Accession numbers are listed in Table 1.

## Statistical analyses

Infection data were analyzed using a generalized linear model (GLZ). Because the dependent variable (infection status) could be either 0 (absent) or 1 (present), a binomial error structure and a logistic link function were used. Colony and caste were considered as independent variables (categorical factors) and the interaction of both variables was added to build the logit model. The Wald test was used to determine the significance of the effects in the model. Binomial confidence intervals for the average infection rates of the different castes were calculated using a normal approximation.

A model has been built as follows, assuming that caste/stage changes are independent of Wolbachia status. Let $I(s, c)$ being the rate of infection for the species $s$ and caste $c$ with species numbered from 1 to 4 and castes/stages from 1 (queen), 2 (larva), 3 (worker), 4 (white soldier) to 5 (mature soldier). Let $a(s, c)$ the probability that an individual for the species $s$ and caste $c$ had lost infection between caste $c-1$ and $c$. Let $b(s, c)$ the probability that an uninfected individual at caste $c-1$ becomes infected as being in caste $c$. Then the model of transmission rate between two castes can be expressed based on two different equations depending on the exact process for the gain of Wolbachia.

In the first case, the loss $a(s,c)$ is applied to the frequency of infected $I(s,c-1)$ individuals and the gain $b(s,c)$ is applied to the non-infected $(1-I(s,c-1))$ individuals at the previous stage $(c-1)$. The equation describing infection rate is therefore:

$$ 
\begin{aligned}
I(s, c) &= 1 - a(s, c)) \cdot I(s, c - 1) + b(s, c) \cdot (1 - I(s, c - 1)) \\
I(s, c) &= I(s, c - 1) \cdot (1 - a(s, c) - b(s, c)) + b(s, c)
\end{aligned}
\qquad (1) $$

Alternatively, the gain could occur later during the current stage. The difference with the previous model is that the proportion of uncontaminated individuals to which the gain is applied is the proportion of the current stage $(1-(I(s, c-1) (1-a(s,c))))$ and not the one at the previous stage $(1-(I(s,c-1)))$ as in equation (1). Then the equation describing infection rate becomes:

$$ 
\begin{aligned}
I(s, c) &= (1 - a(s, c)) \cdot I(s, c - 1) + b(s, c) \cdot (1 - ((1 - a(s, c)) \cdot I(s, c - 1))) \\
I(s, c) &= I(s, c - 1) \cdot (1 - a(s, c)) \cdot (1 - b(s, c)) + b(s, c)
\end{aligned}
\qquad (2) $$

Briefly, in the first model, gain is applied to individuals not infected at the previous stage. In the second model, individuals can lose infection during the current stage and can also gain infection.

The parameters $a(s, c)$ and $b(s, c)$ have been fitted using maximum likelihood with binomial link. This general model can be simplified letting $a(s, c)$ or $b(s, c)$ independent of the species for $a(., c)$ or $b(., c)$, independent of the caste for $a(s,.)$ or $b(s,.)$, or independent of both the caste and the species for $a(.,.)$ or $b(.,.)$. Furthermore, if $a(s, c) = 0$, no stochastic loss is allowed and if $b(s, c) = 0$, no gain is allowed.

For each of the 28 tested models (14 with equation 1 and 14 with equation 2), the quasi-like- lihood of the observed rate at each stage 2 to 5 has been established and transformed to Akaike
Information Content (AIC) [43]. AIC has been corrected for overdispersion (ĉ = 2.9) and small effective size using QAICc [44]. This is a ranking measure that takes into account the quality of the fit of the model to the data while penalizing for the number of parameters used.

The models with the lowest values of QAICc were retained as good candidate models and ΔQAICc was calculated as the difference in value of QAICc between a particular model and the one with the lowest QAICc. Akaike weights ($w_i = \exp(-\Delta_{QAICc}/2)$) normalized to 1) were used to evaluate the relative support of various tested models [44]. Akaike weights can be directly interpreted as conditional probabilities for each model.

After a first round of selection, the selected model has been refined taking into account the fitted values of parameters. Parameters with overlapping confidence intervals were forced to be the same or null if they were not significantly different from 0. Standard errors on the parameter were estimated by inverting the Hessian matrix of second order derivatives of likelihood for each couple of parameters.

# Results

## Infection incidence in different colonies and castes

Wolbachia was detected in all 15 colonies. Infection tests showed that 50% (276/555) of individuals were Wolbachia positive with 90% (18/20) of the sampled queens and kings infected (Table 1). Global offspring infection rate was 48% (258/535) and per caste, 63% for larvae (90/144), 33% for workers (55/169), 74% for white soldiers (57/77) and 39% for mature soldiers (56/145). Mature soldiers and workers had a significantly lower mean infection rate than larvae and white soldiers (GLZ, p < 0.001). White soldiers had a significantly higher mean infection rate than larvae (GLZ, p = 0.0146).

Using a maximum likelihood method, two probabilities were estimated: the probability of an infected individual to become uninfected at the next stage and the probability of an uninfected individual to become infected at the next stage, using two different equations.

The best model has been selected using a two-step procedure. First the model describing the modality of gain has been selected (equation 1 or 2) and second the model has been simplified. The best model out of the 28 tested using Akaike weight ($p = 0.78$) includes a species effect for the stochastic loss of Wolbachia and a caste effect for the rate of gain described with the second equation (Table 3). The second selected model used also equation 2 ($p = 0.11$) and then the probability that the best model is described by equation 2 is around 0.95.

Then, for the selected model (equation 2), the rate of gain for larva, worker and soldier castes was fitted or fixed to 0. Among the new models, the best selected model ($p = 0.68$) was for a species-specific rate of loss (ranging from 0.63 to 0.23 for species C and A respectively) with a null rate of gain for larvae, workers and soldiers and a fitted probability of 0.72, whatever the species, that a white soldier becomes newly contaminated during that stage (Fig. 1b). Although a constant loss of Wolbachia from stage to stage was detected (48% loss in average between two stages for the four species), a higher infection than expected (obtained when the probability of new gain was fixed to zero) was observed in white soldiers (mean %obs = 74 ranging from 73% to 84% for species C and A, respectively, vs. mean %exp = 17 ranging from 4% to 44% for species C and A, respectively if no gain was allowed) and was then responsible for a higher infection than expected in mature soldiers (mean %obs = 38 ranging from 27% to 65% for species C and A vs. mean %exp = 11 ranging from 2% to 34% for species C and A, respectively) (Fig. 1c).

## Infection incidence in tissues and salivary glands

Wolbachia was detected in eggs, ovarian tissues, Malpighian tubules, head and fat tissues of queens and kings, and in thorax + legs of workers and soldiers (Table 2). The heads of white
<table><caption>**Table 3. Characteristics of the tested models.** Selected model based on QAICc is shown in bold.</caption><thead><tr><th></th><th>Model</th><th>Loss</th><th>Gain</th><th>-Ln L</th><th>Parameters</th><th>QAICc</th><th>Akaike weight</th></tr></thead><tbody><tr><td>1</td><td>1</td><td>a(..)</td><td>b(..)</td><td>390.21</td><td>1</td><td>209.12</td><td>0.0037</td></tr><tr><td>2</td><td>1</td><td>a(s,..)</td><td>b(..)</td><td>282.73</td><td>5</td><td>207.43</td><td>0.0086</td></tr><tr><td>3</td><td>1</td><td>a(s,..)</td><td>b(s,..)</td><td>275.11</td><td>8</td><td>212.66</td><td>0.0006</td></tr><tr><td>4</td><td>1</td><td>a(..)</td><td>b(s,..)</td><td>285.91</td><td>5</td><td>209.63</td><td>0.0029</td></tr><tr><td>5</td><td>1</td><td>a(., c)</td><td>b(..)</td><td>284.63</td><td>5</td><td>208.74</td><td>0.0045</td></tr><tr><td>6</td><td>1</td><td>a(., c)</td><td>b(., c)</td><td>275.37</td><td>8</td><td>212.84</td><td>0.0006</td></tr><tr><td>7</td><td>1</td><td>a(..)</td><td>b(., c)</td><td>279.66</td><td>5</td><td>205.32</td><td>0.0248</td></tr><tr><td>8</td><td>1</td><td>a(s, c)</td><td>b(s, c)</td><td>227.37</td><td>28</td><td>389.96</td><td>0.0000</td></tr><tr><td>9</td><td>1</td><td>a(., c)</td><td>b(s, c)</td><td>251.52</td><td>18</td><td>256.82</td><td>0.0000</td></tr><tr><td>10</td><td>1</td><td>a(s,..)</td><td>b(s, c)</td><td>249.77</td><td>18</td><td>255.62</td><td>0.0000</td></tr><tr><td>11</td><td>1</td><td>a(s, c)</td><td>b(., c)</td><td>236.95</td><td>18</td><td>246.79</td><td>0.0000</td></tr><tr><td>12</td><td>1</td><td>a(s, c)</td><td>b(s,..)</td><td>242.07</td><td>18</td><td>250.32</td><td>0.0000</td></tr><tr><td>13</td><td>1</td><td>a(s,..)</td><td>b(., c)</td><td>262.84</td><td>8</td><td>204.21</td><td>0.0433</td></tr><tr><td>14</td><td>1</td><td>a(., c)</td><td>b(s,..)</td><td>273.94</td><td>8</td><td>211.86</td><td>0.0009</td></tr><tr><td>15</td><td>2</td><td>a(..)</td><td>b(..)</td><td>298.05</td><td>2</td><td>209.76</td><td>0.0027</td></tr><tr><td>16</td><td>2</td><td>a(s,..)</td><td>b(..)</td><td>287.67</td><td>5</td><td>210.84</td><td>0.0016</td></tr><tr><td>17</td><td>2</td><td>a(s,..)</td><td>b(s,..)</td><td>285.17</td><td>8</td><td>219.60</td><td>0.0000</td></tr><tr><td>18</td><td>2</td><td>a(..)</td><td>b(s,..)</td><td>285.74</td><td>5</td><td>209.51</td><td>0.0031</td></tr><tr><td>19</td><td>2</td><td>a(., c)</td><td>b(..)</td><td>289.03</td><td>5</td><td>211.78</td><td>0.0010</td></tr><tr><td>20</td><td>2</td><td>a(., c)</td><td>b(., c)</td><td>274.78</td><td>8</td><td>212.44</td><td>0.0007</td></tr><tr><td>21</td><td>2</td><td>a(..)</td><td>b(., c)</td><td>275.17</td><td>5</td><td>202.23</td><td>0.1167</td></tr><tr><td>22</td><td>2</td><td>a(s, c)</td><td>b(s, c)</td><td>229.06</td><td>28</td><td>391.12</td><td>0.0000</td></tr><tr><td>23</td><td>2</td><td>a(., c)</td><td>b(s, c)</td><td>236.52</td><td>18</td><td>246.49</td><td>0.0000</td></tr><tr><td>24</td><td>2</td><td>a(s,..)</td><td>b(s, c)</td><td>234.59</td><td>18</td><td>245.17</td><td>0.0000</td></tr><tr><td>25</td><td>2</td><td>a(s, c)</td><td>b(., c)</td><td>237.00</td><td>18</td><td>246.83</td><td>0.0000</td></tr><tr><td>26</td><td>2</td><td>a(s, c)</td><td>b(s,..)</td><td>264.59</td><td>18</td><td>265.83</td><td>0.0000</td></tr><tr><td><b>27</b></td><td><b>2</b></td><td><b>a(s,..)</b></td><td><b>b(., c)</b></td><td><b>254.43</b></td><td><b>8</b></td><td><b>198.42</b></td><td><b>0.7838</b></td></tr><tr><td>28</td><td>2</td><td>a(., c)</td><td>b(s,..)</td><td>275.88</td><td>8</td><td>213.19</td><td>0.0005</td></tr><tr><td>29</td><td>Full</td><td></td><td></td><td>199.79</td><td>47</td><td></td><td></td></tr></tbody></table>

doi:10.1371/journal.pone.0116070.t003

soldiers and workers, and the two pools of salivary glands were also infected. Although very short, the *wsp* sequences from the salivary glands and those from the entire individuals showed 100% identity (318/318 bp).

Discussion

Wolbachia detection

*Wolbachia* detection was assessed by the way of PCR assays, a method which has been considered as the most efficient and easily applicable to large samples, such as in this study (N = 555) [45]. Here termite colonies were screened for the presence of *Wolbachia* using PCR assays targeting the *wsp* (*Wolbachia* surface protein) gene, which is 10 times more variable in its DNA sequence than 16S rRNA [46]. Despite the fact that the *wsp* gene was shown to be unreliable as a phylogenetic tool due to frequent recombination [47–49], it is still widely used for the detection of *Wolbachia* infection [50,51]. Particularly, the validity of *wsp* to diagnose *Wolbachia* A and B supergroups (i.e. those found in *Cubitermes* spp.) is well documented [45,52,53]. A negative point concerning the sensitivity
of *wsp* gene was recently evoked by Simoes et al. [45] as *wsp* primers were shown to aspecifically detect DNA from bacterial genera that are phylogenetically close to *Wolbachia* such as *Rickettsia*, *Ana-plasma* or *Ehrlichia*. In order to counter this problem, 1–3 *wsp* sequences per colony were obtained from different individuals and none of them showed aspecific amplification. Because nestmate sequences were identical to that of the queen, and because the gain is null in workers that are the sole stage to take food directly from the environment, transfer from infected humus was discarded and the current assumption is that the *Wolbachia* are inherited from the queen.

## Species pattern

All tested colonies were infected in a geographic range of about 0.1ha. As shown by Roy and Harry [33], *Cubitermes* spA and spD harbored B-supergroup *Wolbachia* while *Cubitermes* spB and spC harbored B-supergroup *Wolbachia*. The strict pattern “one *Wolbachia* strain/one *Cubitermes* species” suggests that interspecific horizontal transfers are not common. Compared to *Cubitermes* spB and spC hosting A *Wolbachia* supergroup (Fig. 1c), *Cubitermes* spA hosting B *Wolbachia* supergroup showed a particularly high level of infection for larvae (73%) (a stage which could give either future swarming females able to transmit the bacterium or sterile workers), and a lower rate of loss according to the model. Differences in transmission efficiency linked to the bacterial B-strain harbored by this species could be invoked. *Cubitermes* spD *affinis subarquatus*, hosting also B-*Wolbachia*, showed a moderate level of infection for larvae (40%) but is represented by a unique colony and rates are hardly interpretable in this way. To date, only one study conducted on *Aedes albopictus* demonstrated such a higher transmission efficiency of group B *Wolbachia* relative to supergroup A [54].

## *Wolbachia* loss in adult individuals

Results of the model indicate that rate of loss of *Wolbachia* between two successive stages is constant and independent on the species and caste/stages. Workers and mature soldiers showed a lower infection rate than larvae. This result is particularly clear in *Cubitermes* spC, where five on six colonies had sampled larvae and workers stages and a sharp decrease in observed infection rates. As shown by previous reports on flies, mosquitoes, planthoppers and ants, host age may influence *Wolbachia* infection levels and tissue tropism [12, 26, 30, 55, 56]. Particularly, infection data in ant colonies showed lower infection rates in adult workers compared to worker pupae [26,30]. It was suggested that selection did not favor heritable symbionts conservation in sterile castes (workers and soldiers) because infection could induce a decrease in colony productivity, and that symbiont curing from adult workers may be an adaptation that ensures efficient production of sexual females [2]. Hard environmental conditions were also proposed to explain the elimination of *Wolbachia* in foraging ant workers [26,57]. It is known that intense heat for short periods reduce *Wolbachia* densities and even cure some insects of their bacterial infections in laboratory environments [57–59] and these factors were proposed to have similar effects in natural insect populations [60].

Termite soldiers represent a true caste which is unique among social insects in function and development [61]. Dynamics of *Wolbachia* infection in the soldier caste were explored for the first time in this study and we showed that infection rate is higher in white soldiers than in mature soldiers. Here the white soldier-mature soldier difference is principally driven by *Cubitermes* spB, C and D since one colony in *Cubitermes* spA (T7) showed a reverse pattern.

## The case of white soldiers

In order to explain the higher infection rate in white soldiers, various hypotheses could be examined. Some of them can be discarded such as transfer from infected humus in white soldiers,
because they do not feed by themselves, or the detection of *Wobachia* in alien white soldiers coming from colonies with a high rate of infection. However, colony fusion was never reported in *Cubitermes* species and if, would impact all stages and not only white soldiers. Although inquilinism is important in this genus [62] and particularly common in *Cubitermes* spp. *affinis subarquatus* in the Lopé Reserve (MH, personal observations), soldiers of inquiline genera other than *Cubitermes* are easily morphologically distinguished and no alien genotype was detected in other castes [39].

Another hypothesis could be a *Wolbachia* detection bias due to an increase of bacteria density in white soldiers compared to workers and mature soldiers. In the mosquito *Aedes albopictus* it was shown that the bacterium linked its own replication to that of its host cell, *Wolbachia* densities in eggs being greatest during embryonation and declining throughout diapause [63]. White soldiers are the site of multiple cell divisions: hard mandibles of the future soldier are constructed, salivary glands change their shape and many other morphological modifications occur [36]. White soldiers could constitute a preferential stage for *Wolbachia* replication, thus facilitating the detection of infections. In order to explore this hypothesis, sensitivity tests were carried out using serial dilutions of a number of infected samples from 1 down to 1:100 and these tests showed positive amplification of *wsp* gene, indicating that molecular tools used to detect infection are very powerful. Experiments using real-time quantitative PCR to compare *Wolbachia* density in white soldiers and workers could complement these tests.

The 28 tested models (14 for equation 1 and 14 for equation 2) indicated that *Wolbachia* transfer via horizontal routes between members of the colony (equation 2) has 95% of chance to be the best models among those tested to explain the observed pattern of infection. Among the 14 models, the ones with only stochastic loss of *Wolbachia* and no gain are not supported at all (p<0.001). It was found that *Wolbachia* occurred in reproductive organs of queens and kings but also in somatic tissues of all individuals (head, thorax and legs). The broad tissue trophism of *Wolbachia* is known in many arthropod species [11–13,27,64]. The high rates of infection in white soldiers, associated with the detection of *Wolbachia* in salivary glands of workers, could be the result of a recontamination by infected saliva via trophallaxis between feeding workers and white soldiers. An experimental study concerning trophallaxis and feeding in *Cubitermes fungifaber* demonstrated that white soldiers received salivary secretions, while mature soldiers received regurgitated food (stomodeal) from workers [65]. Field observations also showed that Termitidae workers fed all the colony members but that only white soldiers and young larvae were heavily loaded with pure saliva [36,66].

Various experimental studies described horizontal transmission of *Wolbachia* during close contact of individuals. Rigaud and Juchault [67] reported direct transfer of the symbiont in the isopod *Armadillum vulgare* via blood-to-blood contacts, probably during moulting or by the way of predator injuries. More recently, Le Clec'h et al. [68] proposed cannibalism and predation as paths for horizontal passage of *Wolbachia* between isopods. Direct *Wolbachia* transmission between *Trichogramma* wasp larvae sharing a same food source has also been described when superparasitism occurred, probably by consumption of infected larvae by uninfected ones or during blood-to-blood contact [14,15]. *Wolbachia* has been detected in the salivary glands of several insects, such as mosquitoes [69] and insects feeding on host plants [70,71]. In these last examples, saliva was suggested to be a medium for microorganism transmission, including *Wolbachia*, in an insect-host plant-insect pathway.

In termites, recurrent transfer via feeding workers may be of functional importance and produce other effects than those currently known for *Wolbachia*. In *Acromyrmex* ants, it has been suggested that the high intensity of *Wolbachia* found in the haemolymph and faeces could enhance transmission via parasitoid and the faecal-oral route, or have a role in modulating the immune response of the host [72]. Similarly, it has been proposed that *Wolbachia* found extracellularly in the gut may have a nutritional role [73]. Here the
high proportion of infected white soldiers could represent a side-effect of the specialization of *Wolbachia* in the alimentary canal, due to the same “trophically dependent” behavior of white soldiers and larvae (comprising larvae destined to develop into sexual individuals). The vertical transmission would be then supplemented with horizontal transfers via trophallactic exchanges between infected individuals and future imagoes. Examining the prevalence of infection in the reproductive lineage (nymphs and alates) would probably provide useful information concerning the transmission patterns and dynamics of *Wolbachia* infections. These results suggest new ways of transmission to explore in other social insects and more generally in arthropods since such close interactions between individuals are not limited to termites: trophallaxy and brood care that are common in social insects also occur in burying beetles and earwigs, where larvae or nymphs are fed with regurgitations.

A set of explanations implying patterns of caste differentiation has to be mentioned but cannot be tested with current data because of lack of degree of freedom in the models (more parameters to estimate than information) or lack of information (i.e. absolute age of individuals). In *Cubitermes* species, only a small subset of workers is destined to transition to white soldiers as soldiers represent a small proportion of individuals relative to workers (1.6% in *C. fungifaber*, [74], 0.8–5.1% in *C. sankurensis* [75], 1.4–8.4% in *C. subcrenulatus* [76]). Regarding the terminology of Oster and Wilson [77], these species with a minimal investment in defence are classified as tychophile species [78]. *Cubitermes* soldiers differentiate from the molt of workers that are not morphologically different from other workers [74] and there is no obvious way of knowing whether an individual is destined to become a soldier or not.

In all models, caste/stage changes are assumed to be independent of *Wolbachia* status. However, it remains possible that white soldiers have elevated levels of infection compared to workers because workers infected with *Wolbachia* preferentially develop into white soldiers. Similarly, lower infection rates detected in workers compared to larvae could be explained by the fact that the subset of the larvae transforming into nymphs (i.e. future swarming alates), and not in workers, would be preferentially infected by *Wolbachia*. This raises interesting questions about the benefits and costs of *Wolbachia* in different members of a colony.

White soldiers could also show elevated levels of infection compared to workers because old workers, which would be more likely to have lost *Wolbachia* if *Wolbachia* decline is constant over time, outnumber younger white soldiers (Fig. 1). Our model can appear quite simplified because *Wolbachia* loss and/or gain are treated stage per stage rather than a timely process (e.g. workers are treated as a stage preceding white soldiers). Although very interesting, it seems nevertheless difficult to test or model these hypotheses because individuals could not be traced for infection and life trajectory.

If caste/stage changes are independent of *Wolbachia* status, the more reliable explanation to the high infection rate observed in white soldiers is that *Wolbachia* transfer could be a general phenomenon involving direct transmission from individual to individual via saliva during trophallactic exchanges in termites.

## Acknowledgments

Thanks are due to K. Abernethy and L. White for field facilities and scientific collaboration. We would also like to thank three anonymous reviewers for their many helpful comments. We would like to dedicate this paper to the memory of our colleague D. Lachaise for the fruitful discussions that we have had on this model.
Author Contributions

Conceived and designed the experiments: MH VR MG. Performed the experiments: VR. Analyzed the data: MH VR MG. Contributed reagents/materials/analysis tools: MH VR MG. Wrote the paper: MH VR MG.

References

1. Fenn K, Blaxter M (2006) *Wolbachia* genomes: revealing the biology of parasitism and mutualism. *Trends Parasitol* 22: 60–65. PMID: 16406333

2. Russell JA, Funaro CF, Giraldo YM, Goldman-Huertas B, Suh D, et al. (2012) A veritable menagerie of heritable bacteria from ants, butterflies, and beyond: broad molecular surveys and a systematic review. *PLoS ONE* 7: e51027. doi: 10.1371/journal.pone.0051027 PMID: 23284655

3. Werren JH, Zhang W, Guo LR (1995) Evolution and phylogeny of *Wolbachia*: reproductive parasites of Arthropods. *Proc R Soc Lond, B, Biol Sci* 261: 55–63.

4. Charlat S, Bourtzis K, Merçot H (2002) *Wolbachia*-induced cytoplasmic incompatibility. In: Seckbach J, editor. Symbiosis. Dordrecht, the Netherlands: Kluwer Academic Publisher. pp. 624–644.

5. Fast EM, Toomey ME, Panaram K, Desjardins D, Kolaczyk ED, et al. (2011) *Wolbachia* enhance *Drosophila* stem cell proliferation and target the germline stem cell niche. *Science* 334: 990–992. doi: 10.1126/science.1209609 PMID: 22021671

6. Hoffmann AA, Turelli M (1997) Cytoplasmic incompatibility in insects. In: O'Neill SL, Hoffmann AA, Werren JH, editors. Influential passengers: inherited microorganisms and arthropod reproduction. New York: Oxford University Press. pp. 42–80.

7. Rousset F, Solignac M (1995) Evolution of single and double *Wolbachia* symbioses during speciation in the *Drosophila simulans* complex. Proc Natl Acad Sci USA 92: 6389–6393. PMID: 7604001

8. Werren JH, Baldo L, Clark ME (2008) *Wolbachia*: master manipulators of invertebrate biology. Nat Rev Microbiol 6: 741–751. doi: 10.1038/nrmicro1969 PMID: 18794912

9. Frydman HM, Li JM, Robson DN, Wieschaus E (2006) Somatic stem cell niche tropism in *Wolbachia*. *Nature* 441: 509–512. PMID: 16724067

10. Pfarr KM, Hoerauf A (2007) A niche for *Wolbachia*. *Trends Parasitol* 23: 5–7. PMID: 17113827

11. Cheng Q, Ruel TD, Zhou W, Moloo SK, Majiwa P, et al. (2000) Tissue distribution and prevalence of *Wolbachia* infections in tsetse flies, *Glossina* spp. Med Vet Entomol 14: 44–50. PMID: 10759311

12. Dobson SL, Bourtzis K, Braig HR, Jones BF, Zhou W, et al. (1999) *Wolbachia* infections are distributed throughout insect somatic and germ line tissues. Insect Biochem Mol Biol 29: 153–160. PMID: 10196738

13. Espino CI, Gomez T, Gonzalez G, do Santos MFB, Solano J, et al. (2009) Multi-organ and fecal detection of *Wolbachia* in triatomine *Rhodnius pallescens* (Hemiptera, Reduviidae). Appl Environ Microbiol 75: 547–550. doi: 10.1128/AEM.01665-08 PMID: 19028913

14. Huigens ME, de Almeida RP, Boons PA, Luck RF, Stouthamer R (2004) Natural interspecific and intra-specific horizontal transfer of parthenogenesis-inducing *Wolbachia* in *Trichogramma* wasps. Proc R Soc Lond, B, Biol Sci 271: 509–515.

15. Huigens ME, Luck RF, Klaassen RHG, Maas MFPM, Timmermans MJTN, et al. (2000) Infectious parthenogenesis. *Nature* 405: 178–179. PMID: 10821272

16. Gerth M, GeißLer A, Bleidorn C (2011) *Wolbachia* infections in bees (Anthophila) and possible implications for DNA barcoding. Syst Biodivers 9: 319–327.

17. Gerth M, Röthe J, Bleidorn C (2013) Tracing horizontal *Wolbachia* movements among bees (Anthophila): a combined approach using MLST data and host phylogeny. Mol Ecol 22: 6149–6162. doi: 10.1111/mec.12549 PMID: 24118435

18. Guruprasad NM, Sumithra ML, Puttaraju HP (2012) Long PCR: a sensitive PCR protocol for amplification of *Wolbachia* endosymbiont in Indian honey bees. J Entomol Res 36: 119–122.

19. Pattabhiramaiah M, Brueckner D, Witzel KP, Junier P, Reddy MS (2011) Prevalence of *Wolbachia* in the European honeybee, *Apis mellifera carnaica*. World Appl Sci J 15: 1503–1506.

20. Frost CL, Fernandez-Marin H, Smith JE, Hughes WOH (2010) Multiple gains and losses of *Wolbachia* symbionts across a tribe of fungus-growing ants. Mol Ecol 19: 4077–4085. doi: 10.1111/j.1365-294X.2010.04764.x PMID: 20738784

21. Russell JA, Goldman-Huertas B, Moreau CS, Baldo L, Stahlhut JK, et al. (2009) Specialization and geographic isolation among *Wolbachia* symbionts from ants and lycaenid butterflies. Evolution 63: 624–640. doi: 10.1111/j.1558-5646.2008.00579.x PMID: 19054050
22. Wenseleers T, Ito F, Van Borm S, Huybrechts R, Volckaert F, et al. (1998) Widespread occurrence of the microorganism *Wolbachia* in ants. P Roy Soc B-Biol Sci 265: 1447-1452. PMID: 9721689

23. Ahrens M, Shoemaker DD (2005) Evolutionary history of *Wolbachia* infections in the fire ant *Solenopsis invicta*. BMC Evol Biol 5: 35. PMID: 15927071

24. Bouwma AM, Ahrens ME, DeHeer CJ, Shoemaker DD (2006) Distribution and prevalence of *Wolbachia* in introduced populations of the fire ant *Solenopsis invicta*. Insect Mol Biol 15: 89-93. PMID: 16469072

25. Kautz S, Rubin BER, Moreau CS (2013) Bacterial infections across the ants: frequency and prevalence of *Wolbachia*, *Spiroplasma*, and *Asaia*. *Psyche* 2013: 11.

26. Keller L, Liautard C, Reuter M, Brown WD, Sundstrom L, et al. (2001) Sex ratio and *Wolbachia* infection in the ant *Formica exsecta*. Heredity 87: 227-233. PMID: 11703514

27. Shoemaker DD, Ross KG, Keller L, Vargo EL, Werren JH (2000) *Wolbachia* infections in native and introduced populations of fire ants (Solenopsis spp.). Insect Mol Biol 9: 661-673. PMID: 11122476

28. Tsutsui ND, Kauppinen SN, Oyafuso AF, Grosberg RK (2003) The distribution and evolutionary history of *Wolbachia* infection in native and introduced populations of the invasive argentine ant (*Linepithema humile*). Mol Ecol 12: 3057-3068. PMID: 14629385

29. Viljakainen L, Reuter M, Pamilo P (2008) *Wolbachia* transmission dynamics in *Formica* wood ants. BMC Evol Biol 8: 55. doi: 10.1186/1471-2148-8-55 PMID: 18291041

30. Wenseleers T, Sundstrom L, Billen J (2002) Deleterious *Wolbachia* in the ant *Formica truncorum*. Proc R Soc Lond, B, Biol Sci 269: 623-629.

31. Lo N, Evans TA (2007) Phylogenetic diversity of the intracellular symbiont *Wolbachia* in termites. Mol Phylogenet Evol 44: 461-466. PMID: 17174110

32. Bordenstein S, Rosengaus R (2005) Discovery of a novel *Wolbachia* supergroup in Isoptera. Curr Microbiol 51: 393-398. PMID: 16252129

33. Roy V, Harry M (2007) Diversity of *Wolbachia* isolated from the *Cubitermes* sp. *affinis* *subarquatus* complex of species (Termitidae), revealed by multigene phylogenies. FEMS Microbiol Lett 274: 102-111. PMID: 17663704

34. Salunke B, Salunkhe R, Patole M, Shouche Y (2010) *Wolbachia* and termite association: present status and future implications. J Biosci 35: 171-175. PMID: 20689172

35. Salunke BK, Salunkhe RC, Dhotre DP, Khandagale AB, Walujkar SA, et al. (2010) Diversity of *Wolbachia* in *Odontotermes* spp. (Termitidae) and *Coptotermes heimi* (Rhinotermitidae) using the multigene approach. FEMS Microbiol Lett 307: 55-64. doi: 10.1111/j.1574-6968.2010.01960.x PMID: 20402779

36. Grassé PP (1982) Termitologia. Tome I. Anatomie-Physiologie-Reproduction des termites. Paris: Masson. 715 p.

37. Eggleton P (2011) An introduction to termites: biology, taxonomy and functional morphology. Biology of termites: a modern synthesis. Springer. pp. 1-26.

38. Nalepa C, Bignell D, Bandi C (2001) Detritivory, coprophagy, and the evolution of digestive mutualisms in Dictyoptera. Insectes Soc 48: 194-201.

39. Roy V, Demanche C, Livet A, Harry M (2006) Genetic differentiation in the soil-feeding termite *Cubitermes* sp. *affinis* *subarquatus*: occurrence of cryptic species revealed by nuclear and mitochondrial markers. BMC Evol Biol 6: 102. PMID: 17123444

40. Walsh PS, Metzger DA, Higuchi R (1991) Chelex 100 as a medium for simple extraction of DNA for PCR-based typing from forensic material. Biotechniques 10: 506-513. PMID: 1867860

41. Simon C, Frati F, Beckenbach A, Crespi B, Liu H, et al. (1994) Evolution, weighting, and phylogenetic utility of mitochondrial gene sequences and a compilation of conserved polymerase chain reaction primers. Ann Entomol Soc Am 87: 651-701.

42. Braig HR, Zhou W, Dobson SL, O'Neill SL (1998) Cloning and characterization of a gene encoding the major surface protein of the bacterial endosymbiont *Wolbachia pipientis*. J Bacteriol 180: 2373-2378. PMID: 9573188

43. Akaike H (1974) A new look at the statistical model identification. IEEE Trans Autom Control 19: 716-723.

44. Burnham KP, Anderson DR (2002) Model selection and multimodel inference: A practical information-theoretic approach. New York: Springer-Verlag.

45. Simoes PM, Mialdea G, Reiss D, Sagot MF, Charlat S (2011) *Wolbachia* detection: an assessment of standard PCR Protocols. Mol Ecol Resour 11: 567-572. doi: 10.1111/j.1755-0998.2010.02955.x PMID: 21481216

46. Zhou W, Rousset F, O'Neill S (1998) Phylogeny and PCR-based classification of *Wolbachia* strains using wsp gene sequences. Proc R Soc Lond, B, Biol Sci 265: 509-515.
47. Baldo L, Bordenstein S, Wernegreen JJ, Werren JH (2006) Widespread recombination throughout *Wolbachia* genomes. Mol Biol Evol 23: 437–449. PMID: 16267140

48. Baldo L, Dunning Hotopp JC, Jolley KA, Bordenstein SR, Biber SA, et al. (2006) Multilocus Sequence Typing System for the endosymbiont *Wolbachia pipientis*. Appl Environ Microbiol 72: 7098–7110. PMID: 16936055

49. Baldo L, Lo N, Werren JH (2005) Mosaic nature of the *Wolbachia* surface protein. J Bacteriol 187: 5406–5418. PMID: 16030235

50. Hughes GL, Allsopp PG, Brumbley SM, Woolfit M, McGraw EA, et al. (2011) Variable infection frequency and high diversity of multiple strains of *Wolbachia pipientis* in *Perkinsiella* planthoppers. Appl Environ Microbiol 77: 2165–2168. doi: 10.1128/AEM.02878-10 PMID: 21278277

51. Plantard O, Bouju-Albert A, Malard M-A, Hermouet A, Capron G, et al. (2012) Detection of *Wolbachia* in the tick *Ixodes ricinus* is due to the presence of the Hymenoptera endoparasitoid *Ixodiphagus hookeri*. PLoS ONE 7: e30692. doi: 10.1371/journal.pone.0030692 PMID: 22292021

52. Kikuchi Y, Fukatsu T (2003) Diversity of *Wolbachia* endosymbionts in heteropteran bugs. Appl Environ Microbiol 69: 6082–6090. PMID: 14532065

53. Mikac KM (2007) PCR confirms multiple *Wolbachia* strain infection in Australian and international populations of the invasive stored-product psocid *Liposcelis bostrychophila* Badonnel. J Stored Prod Res 43: 594–597.

54. Kittayapong P, Baisley KJ, Sharpe RG, Baimai V, O'Neill SL (2002) Maternal transmission efficiency of *Wolbachia* superinfections in *Aedes albopictus* populations in Thailand. Am J Trop Med Hyg 66: 103–107. PMID: 12135258

55. Kittayapong P, Mongkalangoon P, Baimai V, O'Neill SL (2002) Host age effect and expression of cytoplasmic incompatibility in field populations of *Wolbachia*-superinfected *Aedes albopictus*. Heredity 88: 270–274. PMID: 11920134

56. Noda H, Koizumi Y, Zhang Q, Deng K (2001) Infection density of *Wolbachia* and incompatibility level in two planthopper species, *Laodelphax striatellus* and *Sogatella furcifera*. Insect Biochem Molec 31: 727–737. PMID: 11267910

57. Hurst GD, Johnson AP, Schulenburg JH, Fuyama Y (2000) Male-killing *Wolbachia* in *Drosophila*: a temperature-sensitive trait with a threshold bacterial density. Genetics 156: 699–709. PMID: 11014817

58. Snook RR, Cleland SY, Wolfner MF, Karr TL (2000) Offsetting effects of *Wolbachia* infection and heat shock on sperm production in *Drosophila simulans*: analyses of fecundity, fertility and accessory gland proteins. Genetics 155: 167–178. PMID: 10790392

59. Stouthamer R, Breeuwer JAJ, Hurst GDD (1999) *Wolbachia pipientis*: microbial manipulator of Arthropod reproduction. Annu Rev Microbiol 53: 71–102. PMID: 10547686

60. Keller GP, Windsor DM, Saucedo JM, Werren JH (2004) Reproductive effects and geographical distributions of two *Wolbachia* strains infecting the Neotropical beetle, *Chelymorpha alternans* Boh. (Chrysomelidae, Cassidinae). Mol Ecol 13: 2405–2420. PMID: 15245413

61. Korb J, Hartfelder K (2008) Life history and development—a framework for understanding developmental plasticity in lower termites. Biological Reviews 83: 295–313. PMID: 18979593

62. Dejean A, Ruelle J (1995) Importance of *Cubitermes* termitaries as shelter for alien incipient termite societies. Insectes Soc 42: 129–136.

63. Ruang-Areerate T, Kittayapong P, McGraw EA, Baimai V, O'Neil S (2004) *Wolbachia* replication and host cell division in *Aedes albopictus*. Curr Microbiol 49: 10–12. PMID: 15297923

64. Min KT, Benzer S (1997) *Wolbachia*, normally a symbiont of *Drosophila*, can be virulent, causing degeneration and early death. Proc Natl Acad Sci USA 94: 10792–10796. PMID: 9380712

65. Alibert J (1963) Echanges trophallactiques chez un termite supérieur. Contamination par le phosphore radio-actif de la population d'un nid de *Cubitermes fungifaber*. Insectes Soc 10: 1–12.

66. Noirot C (1969) Formation of castes in higher termites. In: Krishna K, Weesner FM, editors. Biology of Termites. New York: Academic Press. pp. 311–350.

67. Rigaud T, Juchault P (1995) Success and failure of horizontal transfers of feminizing *Wolbachia* endosymbionts in woodlice. J Evol Biol 8: 249–255.

68. Le Clec'h W, Chevalier FD, Genty L, Bertaux J, Bouchon D, et al. (2013) Cannibalism and predation as paths for horizontal passage of *Wolbachia* between terrestrial Isopods. PLoS ONE 8: e60232. doi: 10.1371/journal.pone.0060232 PMID: 23593179

69. Tsai KH, Huang CG, Wu WJ, Chuang CK, Lin CC, et al. (2006) Parallel infection of Japanese encephalitis virus and *Wolbachia* within cells of mosquito glands. J Med Entomol 43: 752–756. PMID: 16892635
70. Mitsuhashi W, Saiki T, Wei W, Kawakita H, Sato M (2002) Two novel strains of *Wolbachia* coexisting in both species of mulberry leafhoppers. Insect Mol Biol 11: 577–584. PMID: 12421415

71. Sintupachee S, Milne JR, Poonchaisri S, Baimai V, Kittayapong P (2006) Closely related *Wolbachia* strains within the pumpkin arthropod community and the potential for horizontal transmission via the plant. Microb Ecol 51: 294–301. PMID: 16598632

72. Frost CL, Pollock SW, Smith JE, Hughes WO (2014) Wolbachia in the flesh: symbiont intensities in germ-line and somatic tissues challenge the conventional view of *Wolbachia* transmission routes. PloS ONE 9: e95122. doi: 10.1371/journal.pone.0095122 PMID: 24988478

73. Andersen SB, Boye M, Nash DR, Boomsma JJ (2012) Dynamic *Wolbachia* prevalence in *Acromyrmex* leaf-cutting ants: potential for a nutritional symbiosis. J Evol Biol 25: 1340–1350. doi: 10.1111/j.1420-9101.2012.02521.x PMID: 22530696

74. Han S, Noirot C (1983) Développement de la jeune colonie chez *Cubitermes fungifaber* (Sjöstedt) (Isoptera, Termitidae). Ann Soc Entomol Fr 19: 413–420.

75. Bouillon A, Lekie R (1964) Populations, rythme d'activité diurne et cycle de croissance du nid de *Cubitermes sankurensis* Wasmann (Isoptera, Termitinae). In: Bouillon A, editor. Études sur les termites africains, Chapter 19. Paris: Masson. pp. 197–213.

76. Bodot P (1970) Composition des colonies de *Cubitermes subcrenulatus* Silvestri (Isoptera Termitidae). C R Acad Sci Hebd Seances Acad Sci D 271: 327–330. PMID: 4317829

77. Oster GF, Wilson EO (1978) Caste and ecology in the social insects: Princeton University Press. PMID: 740003

78. Noirot C, Darlington J (2000) Termite nests: architecture, regulation and defence. In: Abe T, Bignell DE, Higashi M, editors. Termites: evolution, sociality, symbioses, ecology. Dordrecht, the Netherlands: Kluwer Academic Publishers. pp. 121–139.

79. Han SH (1987) Fondation et croissance des colonies de termites supérieurs: Doctoral dissertation, Dijon.

80. Noirot C (1955) Recherches sur le polymorphisme des termites supérieurs (Termitidae): Vol. 11, Masson et Cie, Paris.

81. Williams RM (1959) Colony development in *Cubitermes ugandensis* Fuller (Isoptera: Termitidae). Insectes Soc 6: 291–304.