Tempo and mode of early gene loss in endosymbiotic bacteria from insects

Background  

Understanding evolutionary processes that drive genome reduction requires determining the tempo (rate) and the mode (size and types of deletions) of gene losses. In this study, we analysed five endosymbiotic genome sequences of the gamma-proteobacteria (three different Buchnera aphidicola strains, Wigglesworthia glossinidia, Blochmannia floridanus) to test if gene loss could be driven by the selective importance of genes. We used a parsimony method to reconstruct a minimal ancestral genome of insect endosymbionts and quantified gene loss along the branches of the phylogenetic tree. To evaluate the selective or functional importance of genes, we used a parameter that measures the level of adaptive codon bias in E. coli (i.e. codon adaptive index, or CAI), and also estimates of evolutionary rates (Ka) between pairs of orthologs either in free-living bacteria or in pairs of symbionts.     

Host-symbiont stability and fast evolutionary rates in an ant-bacterium association: cospeciation of camponotus species and their endosymbionts, candidatus blochmannia

Bacterial endosymbionts are widespread across several insect orders and are involved in interactions ranging from obligate mutualism to reproductive parasitism. Candidatus Blochmannia gen. nov. (Blochmannia) is an obligate bacterial associate of Camponotus and related ant genera (Hymenoptera: Formicidae). The occurrence of Blochmannia in all Camponotus species sampled from field populations and its maternal transmission to host offspring suggest that this bacterium is engaged in a long-term, stable association with its ant hosts. However, evidence for cospeciation in this system is equivocal because previous phylogenetic studies were based on limited gene sampling, lacked statistical analysis of congruence, and have even suggested host switching. We compared phylogenies of host genes (the nuclear EF-1alphaF2 and mitochondrial COI/II) and Blochmannia genes (16S ribosomal DNA [rDNA], groEL, gidA, and rpsB), totaling more than 7 kilobases for each of 16 Camponotus species. Each data set was analyzed using maximum likelihood and Bayesian phylogenetic reconstruction methods. We found minimal conflict among host and symbiont phylogenies, and the few areas of discordance occurred at deep nodes that were poorly supported by individual data sets. Concatenated protein-coding genes produced a very well-resolved tree that, based on the Shimodaira-Hasegawa test, did not conflict with any host or symbiont data set. Correlated rates of synonymous substitution (d(S)) along corresponding branches of host and symbiont phylogenies further supported the hypothesis of cospeciation. These findings indicate that Blochmannia-Camponotus symbiosis has been evolutionarily stable throughout tens of millions of years. Based on inferred divergence times among the ant hosts, we estimated rates of sequence evolution of Blochmannia to be approximately 0.0024 substitutions per site per million years (s/s/MY) for the 16S rDNA gene and approximately 0.1094 s/s/MY at synonymous positions of the genes sampled. These rates are several-fold higher than those for related bacteria Buchnera aphidicola and Escherichia coli. Phylogenetic congruence among Blochmannia genes indicates genome stability that typifies primary endosymbionts of insects.     

Endosymbiosis in statu nascendi: close phylogenetic relationship between obligately endosymbiotic and obligately free-living Polynucleobacter strains (Betaproteobacteria)

Bacterial strains affiliated to the phylogenetically shallow subcluster C (PnecC) of the Polynucleobacter cluster, which is characterized by a minimal 16S rRNA gene sequence similarity of approximately 98.5%, have been reported to occur as obligate endosymbionts of ciliates (Euplotes spp.), as well as to occur as free-living cells in the pelagic zone of freshwater habitats. We investigated if these two groups of closely related bacteria represent strains fundamentally differing in lifestyle, or if they simply represent different stages of a facultative endosymbiotic lifestyle. The phylogenetic analysis of 16S rRNA gene and 16S-23S ITS sequences of five endosymbiont strains from two different Euplotes species and 40 pure culture strains demonstrated host-species-specific clustering of the endosymbiont sequences within the PnecC subcluster. The sequences of the endosymbionts showed characteristics indicating an obligate endosymbiotic lifestyle. Cultivation experiments revealed fundamental differences in physiological adaptations, and determination of the genome sizes indicated a slight size reduction in endosymbiotic strains. We conclude that the two groups of PnecC bacteria represent obligately free-living and obligately endosymbiotic strains, respectively, and do not represent different stages of the same complex life cycle. These closely related strains occupy completely separated ecological niches. To our best knowledge, this is the closest phylogenetic relationship between obligate endosymbionts and obligately free-living bacteria ever revealed.     

A conservative test of genetic drift in the endosymbiotic bacterium Buchnera: slightly deleterious mutations in the chaperonin groEL

The obligate endosymbiotic bacterium Buchnera aphidicola shows elevated rates of sequence evolution compared to free-living relatives, particularly at nonsynonymous sites. Because Buchnera experiences population bottlenecks during transmission to the offspring of its aphid host, it is hypothesized that genetic drift and the accumulation of slightly deleterious mutations can explain this rate increase. Recent studies of intraspecific variation in Buchnera reveal patterns consistent with this hypothesis. In this study, we examine inter- and intraspecific nucleotide variation in groEL, a highly conserved chaperonin gene that is constitutively overexpressed in Buchnera. Maximum-likelihood estimates of nonsynonymous substitution rates across Buchnera species are strikingly low at groEL compared to other loci. Despite this evidence for strong purifying selection on groEL, our intraspecific analysis of this gene documents reduced synonymous polymorphism, elevated nonsynonymous polymorphism, and an excess of rare alleles relative to the neutral expectation, as found in recent studies of other Buchnera loci. Comparisons with Escherichia coli generally show patterns predicted by their differences in N(e). The sum of these observations is not expected under relaxed or balancing selection, selective sweeps, or increased mutation rate. Rather, they further support the hypothesis that drift is an important force driving accelerated protein evolution in this obligate mutualist.     

Genome rearrangement distances and gene order phylogeny in gamma-Proteobacteria

Genome rearrangements have been studied in 30 gamma-proteobacterial complete genomes by comparing the order of a reduced set of genes on the chromosome. This set included those genes fulfilling several characteristics, the main ones being that an ortholog was present in every genome and that none of them had been acquired by horizontal gene transfer. Genome rearrangement distances were estimated based on either the number of breakpoints or the minimal number of inversions separating two genomes. Breakpoint and inversion distances were highly correlated, indicating that inversions were the main type of rearrangement event in gamma-Proteobacteria. In general, the progressive increase in sequence-based distances between genome pairs was associated with the increase in their rearrangement-based distances but with several groups of distances not following this pattern. Compared with free-living enteric bacteria, the lineages of Pasteurellaceae were evolving, on average, to relatively higher rates of between 2.02 and 1.64, while the endosymbiotic bacterial lineages of Buchnera aphidicola and Wigglesworthia glossinidia were evolving at moderately higher rates of 1.38 and 1.35, respectively. Because we know that the rearrangement rate in the Bu. aphidicola lineage was close to zero during the last 100-150 Myr of evolution, we deduced that a much higher rate took place in the first period of lineage evolution after the divergence of the Escherichia coli lineage. On the other hand, the lineage of the endosymbiont Blochmannia floridanus did present an almost identical rate to free-living enteric bacteria, indicating that the increase in the genome rearrangement rate is not a general change associated with bacterial endosymbiosis. Phylogenetic reconstruction based on rearrangement distances showed a different topology from the one inferred by sequence information. This topology broke the proposed monophyly of the three endosymbiotic lineages and placed Bl. floridanus as a closer relative to E. coli than Yersinia pestis. These results indicate that the phylogeny of these insect endosymbionts is still an open question that will require the development of specific phylogenetic methods to confirm whether the sisterhood of the three endosymbiotic lineages is real or a consequence of a long-branch attraction phenomenon.     

An effect of 16S rRNA intercistronic variability on coevolutionary analysis in symbiotic bacteria: molecular phylogeny of Arsenophonus triatominarum

The genes of ribosomal RNA are the most popular and frequently used markers for bacterial phylogeny and reconstruction of insect-symbiont coevolution. In primary symbionts, such as Buchnera and Wigglesworthia, genome economization leads to the establishment of a single copy of these sequences. In phylogenetic studies, they provide sufficient information and yield phylogenetic trees congruent with host evolution. In contrast, other symbiotic lineages (e.g., the genus Arsenophonus) carry a higher number of rRNA copies in their genomes, which may have serious consequences for phylogenetic inference. In this study, we show that in Arsenophonus triatominarum the degree of heterogeneity can affect reconstruction of phylogenetic relationships and mask possible coevolution between the symbiont and its host. Phylogenetic arrangement of individual rRNA copies was used, together with a calculation of their divergence time, to demonstrate that the incongruent 16S rDNA trees and low nucleotide diversity in the secondary symbiont could be reconciled with the coevolutionary scenario.     

Multiple acquisition of methanogenic archaeal symbionts by anaerobic ciliates

Anaerobic heterotrichous ciliates (Armophoridae and Clevelandellidae) possess hydrogenosomes that generate molecular hydrogen and ATP. This intracellular source of hydrogen provides the basis for a stable endosymbiotic association with methanogenic archaea. We analyzed the SSU rRNA genes of 18 heterotrichous anaerobic ciliates and their methanogenic endosymbionts in order to unravel the evolution of this mutualistic association. Here, we show that the anaerobic heterotrichous ciliates constitute at least three evolutionary lines. One group consists predominantly of gut-dwelling ciliates, and two to three, potentially four, additional clades comprise ciliates that thrive in freshwater sediments. Their methanogenic endosymbionts belong to only two different taxa that are closely related to free-living methanogenic archaea from the particular ecological niches. The close phylogenetic relationships between the endosymbionts and free-living methanogenic archaea argue for multiple acquisitions from environmental sources, notwithstanding the strictly vertical transmission of the endosymbionts. Since phylogenetic analysis of the small-subunit (SSU) rRNA genes of the hydrogenosomes of these ciliates indicates a descent from the mitochondria of aerobic ciliates, it is likely that anaerobic heterotrichous ciliates hosted endosymbiotic methanogens prior to their radiation. Therefore, our data strongly suggest multiple acquisitions and replacements of endosymbiotic methanogenic archaea during their host's adaptation to the various ecological niches.     

Increased rates of sequence evolution in endosymbiotic bacteria and fungi with small effective population sizes

Mutualistic, maternally transmitted endosymbiotic microorganisms undergo severe population bottlenecks at each host generation, resulting in a reduction in effective population size (Ne). Previous studies of Buchnera, the primary endosymbiont of aphids, and of several other species of endosymbiotic bacteria have shown that these species exhibit an increase in the rate of substitution of slightly deleterious mutations, among other predicted effects of increased drift due to small Ne, such as reduced codon bias. However, these studies have been limited in taxonomic scope, and it was therefore not clear whether the increase in rate is a general feature of endosymbiont lineages. Here, we test the prediction that a long-term reduction in Ne causes an increase in substitution rate using DNA sequences of the 16S rRNA gene from 13 phylogenetically independent comparisons between taxonomically diverse endosymbiotic microorganisms and their free-living relatives. Maximum likelihood and distance-based methods both indicate a significant increase in substitution rate in a wide range of bacterial and fungal endosymbionts compared to closely related free-living lineages. We use the same data set to test whether 16S genes from endosymbionts display increased A + T content, another indicator of increased genetic drift, and find that there is no significant difference in base composition between endosymbiont and nonendosymbiont 16S genes. However, analysis of an additional data set of whole bacterial genomes demonstrates that, while host-dependent bacteria have significantly increased genomic A + T content, the base content of the 16S gene tends to vary less than that of the whole genome. It is possible that selection for stability of rRNA is strong enough to overcome the effects of drift toward increased A + T content in endosymbiont 16S genes, despite the reduced effective population sizes of these organisms.     

A phylogenomic study of endosymbiotic bacteria

Endosymbiotic bacteria of aphids, Buchnera aphidicola, and tsetse flies, Wigglesworthia glossinidia, are descendents of free-living gamma-Proteobacteria. The acceleration of sequence evolution in the endosymbiont genomes is here estimated from a phylogenomic analysis of the gamma-Proteobacteria. The tree topologies associated with the most highly conserved genes suggest that the endosymbionts form a sister group with Escherichia coli, Salmonella sp., and Yersinia pestis. Our results indicate that deviant tree topologies result from high substitution rates and biased nucleotide patterns, rather than from lateral gene transfer, as previously suggested. A reinvestigation of the relative rate increase in the endosymbiont genomes reveals variability among genes that correlate with host-associated metabolic dependencies. The conclusion is that host-level selection has retarded both the loss of genes and the acceleration of sequence evolution in endocellular symbionts.     

蚜虫与其初级内共生菌进化关系: 假说及演化机理

蚜虫是半翅目(Hemiptera)中一类取食植物韧皮部汁液的昆虫, 由于具有一些独特的生物学特征, 是研究重要适应进化问题的理想模型。蚜虫体内存在一类专性的胞内共生菌Buchnera, 对于蚜虫营养代谢和正常发育有重要贡献, 被称为蚜虫的初级内共生菌。蚜虫-Buchnera是研究共生关系的理想模型, 两者系统发育格局的研究有助于理解生物间专性共生关系的演化。本文系统综述了两者在不同分类水平（高级阶元、 低级阶元）上的系统发育关系研究。现有证据暗示： 两者在低级阶元水平上具有系统发育一致性, 而在高级阶元水平上可能没有平行演化关系, 这些对早期研究提出的平行演化假说提出了质疑。在现有研究的基础上, 本文建议从增加取样类群、 增加基因数目和数据量、 系统发育一致性检验等几个方面开展更深入的系统发育研究, 并开展Buchnera的转移实验, 从而检验Buchnera的横向转移及其基因在不同蚜虫支系中的进化速率一致性, 以便更客观地揭示蚜虫-Buchnera的进化关系。     

Assessment of bacterial endosymbiont diversity in <Emphasis Type="Italic">Otiorhynchus</Emphasis>spp. (Coleoptera: Curculionidae) larvae using a multitag 454 pyrosequencing approach

BACKGROUND: Weevils of the genus Otiorhynchus are regarded as devastating pests in a wide variety of horticultural crops worldwide. So far, little is known on the presence of endosymbionts in Otiorhynchus spp.. Investigation of endosymbiosis in this genus may help to understand the evolution of different reproductive strategies in these weevils (parthenogenesis or sexual reproduction), host-symbiont interactions, and may provide a future basis for novel pest management strategy development. Here, we used a multitag 454 pyrosequencing approach to assess the bacterial endosymbiont diversity in larvae of four economically important Otiorhynchus species. RESULTS: High-throughput tag-encoded FLX amplicon pyrosequencing of a bacterial 16S rDNA fragment was used to characterise bacterial communities associated with different Otiorhynchus spp. larvae. By sequencing a total of ~48,000 PCR amplicons, we identified 49 different operational taxonomic units (OTUs) as bacterial endosymbionts in the four studied Otiorhynchus species. More than 90% of all sequence reads belonged either to the genus Rickettsia or showed homology to the phylogenetic group of "Candidatus Blochmannia" and to endosymbionts of the lice Pedicinus obtusus and P. badii. By using specific primers for the genera Rickettsia and "Candidatus Blochmannia", we identified a new phylogenetic clade of Rickettsia as well as "Candidatus Nardonella" endosymbionts in Otiorhynchus spp. which are closely related to "Candidatus Blochmannia" bacteria. CONCLUSIONS: Here, we used multitag 454 pyrosequencing for assessment of insect endosymbiotic communities in weevils. As 454 pyrosequencing generates only quite short sequences, results of such studies can be regarded as a first step towards identifying respective endosymbiotic species in insects. In the second step of our study, we analysed sequences of specific gene regions for a more detailed phylogeny of selected endosymbiont genera. As a result we identified the presence of Rickettsia and "Candidatus Nardonella" endosymbionts in Otiorhynchus spp.. This knowledge is an important step in exploring bacteria-insect associations for potential use in insect pest control.     

From phylogenetics to phylogenomics: the evolutionary relationships of insect endosymbiotic gamma-Proteobacteria as a test case

The increasing availability of complete genome sequences and the development of new, faster methods for phylogenetic reconstruction allow the exploration of the set of evolutionary trees for each gene in the genome of any species. This has led to the development of new phylogenomic methods. Here, we have compared different phylogenetic and phylogenomic methods in the analysis of the monophyletic origin of insect endosymbionts from the gamma-Proteobacteria, a hotly debated issue with several recent, conflicting reports. We have obtained the phylogenetic tree for each of the 579 identified protein-coding genes in the genome of the primary endosymbiont of carpenter ants, Blochmannia floridanus, after determining their presumed orthologs in 20 additional Proteobacteria genomes. A reference phylogeny reflecting the monophyletic origin of insect endosymbionts was further confirmed with different approaches, which led us to consider it as the presumed species tree. Remarkably, only 43 individual genes produced exactly the same topology as this presumed species tree. Most discrepancies between this tree and those obtained from individual genes or by concatenation of different genes were due to the grouping of Xanthomonadales with beta-Proteobacteria and not to uncertainties over the monophyly of insect endosymbionts. As previously noted, operational genes were more prone to reject the presumed species tree than those included in information-processing categories, but caution should be exerted when selecting genes for phylogenetic inference on the basis of their functional category assignment. We have obtained strong evidence in support of the monophyletic origin of gamma-Proteobacteria insect endosymbionts by a combination of phylogenetic and phylogenomic methods. In our analysis, the use of concatenated genes has shown to be a valuable tool for analyzing primary phylogenetic signals coded in the genomes. Nevertheless, other phylogenomic methods such as supertree approaches were useful in revealing alternative phylogenetic signals and should be included in comprehensive phylogenomic studies.     

毛蚜属Chaitophorus蚜虫(半翅目:蚜科:毛蚜亚科)与初级内共生菌Buchnera的演化关系

【目的】蚜虫体内共生菌种类丰富,二者关系十分密切。几乎所有蚜虫都具有一类专性的初级内共生菌Buchnera aphidicola,二者的专性共生关系使蚜虫-Buchnera成为研究共生关系演化的理想模型。本研究对蚜虫-Buchnera在低级阶元水平上的"平行演化假说"进行了验证。【方法】选取在杨属Populus或柳属Salix植物上营同寄主全周期生活的毛蚜属Chaitophorus蚜虫作为研究对象,基于不同来源的分子标记(蚜虫线粒体基因、核基因和内共生菌基因),运用最大似然法和贝叶斯法重建蚜虫和Buchnera的系统树,并利用Tree Map、Jane和Para Fit检验二者是否具有协同系统发生关系。【结果】Tree Map和Jane分析检测到毛蚜属蚜虫与Buchnera具有显著的共成种信号,Para Fit分析结果表明二者的总体关联极为显著。【结论】毛蚜属蚜虫与其初级内共生菌Buchnera在种级及以下水平上符合"平行演化假说",并且二者的演化关系不会受到寄主植物差异的影响。     

Phylogenetic analysis of a highly specific association between ectosymbiotic, sulfur-oxidizing bacteria and a marine nematode.

The phylogenetic relationship of chemoautotrophic, sulfur-oxidizing, ectosymbiotic bacteria growing on a marine nematode, a Laxus sp. (formerly a Catanema sp.), to known endosymbionts and free-living bacteria was determined. Comparative 16S rRNA sequencing was used to investigate the unculturable nematode epibionts, and rRNA-targeted oligonucleotide hybridization probes were used to identify the ectosymbionts in situ. Both analyses revealed a remarkably specific and stable symbiosis. Unique hybridization of a specific probe to the ectosymbionts indicated that only one species of bacteria was present and growing on the cuticle of the nematode. Distance and parsimony methods used to infer phylogenetic trees both placed the nematode ectosymbionts at the base of a branch containing chemoautotrophic, sulfur-oxidizing endosymbionts of three bivalve families and of the tube worm Riftia pachyptila. The most closely related free-living bacteria were chemoautotrophic sulfur oxidizers belonging to the genus Thiomicrospira. Furthermore, our results suggested that a second, only distantly related group of thioautotrophic endosymbionts has as its deepest branch surface-colonizing bacteria belonging to the genus Thiothrix, some of which are capable of sulfur-oxidizing chemoautotrophic growth.     

Accurate Phylogenetic Tree Reconstruction from Quartets: A Heuristic Approach

Supertree methods construct trees on a set of taxa (species) combining many smaller trees on the overlapping subsets of the entire set of taxa. A ‘quartet’ is an unrooted tree over taxa, hence the quartet-based supertree methods combine many -taxon unrooted trees into a single and coherent tree over the complete set of taxa. Quartet-based phylogeny reconstruction methods have been receiving considerable attentions in the recent years. An accurate and efficient quartet-based method might be competitive with the current best phylogenetic tree reconstruction methods (such as maximum likelihood or Bayesian MCMC analyses), without being as computationally intensive. In this paper, we present a novel and highly accurate quartet-based phylogenetic tree reconstruction method. We performed an extensive experimental study to evaluate the accuracy and scalability of our approach on both simulated and biological datasets.     

On the use of nucleic acid sequences to infer early branchings in the tree of life

Simplifying assumptions made in various tree reconstruction methods-- notably rate constancy among nucleotide sites, homogeneity, and stationarity of the substitutional processes--are clearly violated when nucleotide sequences are used to infer distant relationships. Use of tree reconstruction methods based on such oversimplified assumptions can lead to misleading results, as pointed out by previous authors. In this paper, we made use of a (discretized) gamma distribution to account for variable rates of substitution among sites and built models that allowed for unequal base frequencies in different sequences. The models were nonhomogeneous Markov-process models, assuming different patterns of substitution in different parts of the tree. Data of the small-subunit rRNAs from four species were analyzed, where base frequencies were quite different among sequences and rates of substitution were highly variable at sites. Parameters in the models were estimated by maximum likelihood, and models were compared by the likelihood-ratio test. The nonhomogeneous models provided significantly better fit to the data than homogeneous models despite their involvement of many parameters. They also appeared to produce reasonable estimation of the phylogenetic tree; in particular, they seemed able to identify the root of the tree.      

玉米蚜体内参与传毒的共生菌groEL基因的克隆和原核表达

以玉米蚜杨凌生物型为材料,设计特异性引物采用PCR的方法在国内首先克隆了一种玉米蚜体内参与传毒的共生菌groEL基因,序列测定结果表明：玉米蚜杨凌生物型共生菌groEL奏长为1647bp,编码548个氨基酸,登录Genebank,序列号为AF387863。构建了该基因的原核表达载体,用pBV221表达出63KDa的非融合目的蛋白,用pET-3a表达出69KDa的融合蛋白,二者均有较高的表达量。     

GenNon-h: Generating multiple sequence alignments on nonhomogeneous phylogenetic trees

ABSTRACT: BACKGROUND: A number of software packages are available to generate DNA multiple sequence alignments (MSAs) evolved under continuous-time Markov processes on phylogenetic trees. On the other hand, methods of simulating the DNA MSA directly from the transition matrices do not exist. Moreover, existing software restricts to the time-reversible models and it is not optimized to generate nonhomogeneous data (i.e. placing distinct substitution rates at different lineages). RESULTS: We present the first package designed to generate MSAs evolving under discrete-time Markov processes on phylogenetic trees, directly from probability substitution matrices. Based on the input model and a phylogenetic tree in the Newick format (with branch lengths measured as the expected number of substitutions per site), the algorithm produces DNA alignments of desired length. GenNon-h is publicly available for download. CONCLUSION: The software presented here is an efficient tool to generate DNA MSAs on a given phylogenetic tree. GenNon-h provides the user with the nonstationary or nonhomogeneous phylogenetic data that is well suited for testing complex biological hypotheses, exploring the limits of the reconstruction algorithms and their robustness to such models.     

Inferring pattern and process: maximum-likelihood implementation of a nonhomogeneous model of DNA sequence evolution for phylogenetic analysis

A nonhomogeneous, nonstationary stochastic model of DNA sequence evolution allowing varying equilibrium G + C contents among lineages is devised in order to deal with sequences of unequal base compositions. A maximum-likelihood implementation of this model for phylogenetic analyses allows handling of a reasonable number of sequences. The relevance of the model and the accuracy of parameter estimates are theoretically and empirically assessed, using real or simulated data sets. Overall, a significant amount of information about past evolutionary modes can be extracted from DNA sequences, suggesting that process (rates of distinct kinds of nucleotide substitutions) and pattern (the evolutionary tree) can be simultaneously inferred. G + C contents at ancestral nodes are quite accurately estimated. The new method appears to be useful for phylogenetic reconstruction when base composition varies among compared sequences. It may also be suitable for molecular evolution studies.      

The obligate mutualist Wigglesworthia glossinidia influences reproduction, digestion, and immunity processes of its host, the tsetse fly

Tsetse flies (Diptera: Glossinidae) are vectors for trypanosome parasites, the agents of the deadly sleeping sickness disease in Africa. Tsetse also harbor two maternally transmitted enteric mutualist endosymbionts: the primary intracellular obligate Wigglesworthia glossinidia and the secondary commensal Sodalis glossinidius. Both endosymbionts are transmitted to the intrauterine progeny through the milk gland secretions of the viviparous female. We administered various antibiotics either continuously by per os supplementation of the host blood meal diet or discretely by hemocoelic injections into fertile females in an effort to selectively eliminate the symbionts to study their individual functions. A symbiont-specific PCR amplification assay and fluorescence in situ hybridization analysis were used to evaluate symbiont infection outcomes. Tetracycline and rifampin treatments eliminated all tsetse symbionts but reduced the fecundity of the treated females. Ampicillin treatments did not affect the intracellular Wigglesworthia localized in the bacteriome organ and retained female fecundity. The resulting progeny of ampicillin-treated females, however, lacked Wigglesworthia but still harbored the commensal Sodalis. Our results confirm the presence of two physiologically distinct Wigglesworthia populations: the bacteriome-localized Wigglesworthia involved with nutritional symbiosis and free-living Wigglesworthia in the milk gland organ responsible for maternal transmission to the progeny. We evaluated the reproductive fitness, longevity, digestion, and vectorial competence of flies that were devoid of Wigglesworthia. The absence of Wigglesworthia completely abolished the fertility of females but not that of males. Both the male and female Wigglesworthia-free adult progeny displayed longevity costs and were significantly compromised in their blood meal digestion ability. Finally, while the vectorial competence of the young newly hatched adults without Wigglesworthia was comparable to that of their wild-type counterparts, older flies displayed higher susceptibility to trypanosome infections, indicating a role for the mutualistic symbiosis in host immunobiology. The ability to rear adult tsetse that lack the obligate Wigglesworthia endosymbionts will now enable functional investigations into this ancient symbiosis.