Open access peer-reviewed chapter

The Symbiome of Llaveia Cochineals (Hemiptera: Coccoidea: Monophlebidae) Includes a Gammaproteobacterial Cosymbiont Sodalis TME1 and the Known Candidatus Walczuchella monophlebidarum

Written By

Tania Rosas-Pérez, Arturo Vera-Ponce de León, Mónica Rosenblueth, Shamayim T. Ramírez-Puebla, Reiner Rincón-Rosales, Julio Martínez-Romero, Michael F. Dunn, Éva Kondorosi and Esperanza Martínez-Romero

Reviewed: 20 October 2016 Published: 12 April 2017

DOI: 10.5772/66442

From the Edited Volume

Insect Physiology and Ecology

Edited by Vonnie D.C. Shields

Chapter metrics overview

1,746 Chapter Downloads

View Full Metrics


The genome and transcriptome of the endosymbiotic flavobacterium Candidatus Walczuchella monophlebidarum revealed its role in the synthesis of essential amino acids for its host, the wax cochineal Llaveia axin axin. There were, however, missing genes in the endosymbiont for some biosynthetic pathways. Here, we characterized TME1, another cochineal symbiont that may metabolically complement Walczuchella. TME1 was ascribed to the gammaproteobacterial genus Sodalis on a phylogenomic basis using gene sequences from 143 proteins core genome sequences and the core average nucleotide identity (ANI) confirmed its position. Additionally, we describe Sodalis as a coherent genus. TME1 genome is around 3.4 Mb and has complete gene sequences for the biosynthesis of 10 essential amino acids, for polyamines, flagella, nitrate respiration, and detoxification among many others. Transcripts from ovaries and bacteriomes allowed the identification of differentially transcribed genes from the endosymbionts and host. Highly transcribed genes were identified in TME1 and transcripts involved in amino acid biosynthesis were found. We review here that cosymbionts that derived from different bacterial classes and genera seem to be advantageous for insects that have Flavobacteria as the primary endosymbionts.


  • endosymbionts
  • scale insect
  • Gammaproteobacteria
  • Sodalis-like
  • Alphaproteobacteria
  • fungi

1. Introduction

All organisms are inhabited by microbes that exert different effects on their hosts. In insects, there are many examples of beneficial associations with symbiotic microbes that have been linked to the insect ecological success. Symbionts that are vertically transmitted from mother to offspring and with an intrinsic interdependence with the insect host are considered as primary endosymbionts and they have reduced genomes [1, 2]; they do not grow on standard laboratory media. In theory, endosymbionts evolved from gut bacteria [3] that are largely more complex and may be determined by the diet and the environment. Primary endosymbionts may reside inside insect cells called bacteriocytes that may be found in specialized host structures called bacteriomes. Bacteriomes may be equivalent to plant-root nodules considering that they are host structures harboring particular bacterial species with specific roles [4]. But even in plants, cosymbionts have been encountered; for example, the slow-growing actinobacteria Micromonospora is found in nodules formed by Bradyrhizobium, Rhizobium, or Frankia in several legumes or actinorhizal roots, although Micromonospora is unable to form nodules [5]. Micromonospora has been reported to enhance nodulation and promote plant growth, may enhance plant defense responses, or inhibit pathogens [6].

In insects, cosymbiosis is not uncommon and there are cases in which two or more bacterial symbionts are found in the bacteriome [7, 8]. Additionally, other microbes including fungi may be found in the hemolymph or in different insect tissues [911]. Fungal symbionts may be found as well in specialized insect structures known as mycangia [12] or inside insect cells called mycetocytes [13].

In insects, primary bacterial endosymbionts synthesize essential amino acids or vitamins for their hosts and reside intracellularly in bacteriomes. In some cases, complementation of metabolic pathways seems to occur among different insect symbionts [1417]. Additionally, cosymbionts may have different roles, and some have been implicated in defense [1821], tolerance to stress [22], resistance to high temperatures [2325], to virus [2628], or may manipulate sex differentiation [29]. There is an example in which a secondary endosymbiont substituted a lost primary Buchnera symbiont in an aphid [30]. Among others, alpha, gamma and betaproteobacteria have been found as cosymbionts; for example, the primary endosymbiont Candidatus Sulcia muelleri (“Sulcia” from here on) (phylum Bacteroidetes, class Flavobacteria) with a highly reduced genome has betaproteobacteria as cosymbionts found in green rice leafhoppers [7], stinkbugs [31], and spittlebugs [32, 33]. In leafhoppers, the symbionts occupy different types of bacteriocytes that constitute the outer or inner regions of the bacteriome [7]. The Sulcia cosymbionts are Hodgkinia, Zinderia, Nasuia [34, 35] with very small genomes, and the gammaproteobacteria Baummania, Arsenophonus, or Sodalis, the latter considered as a new acquisition. Surprisingly a gammaproteobacterium may be found inside Sulcia cells and be transmitted to the next generation [36].

Scale insects (Hemiptera: Coccoidea) feed on plant sap, which is a nutritionally poor diet that lacks most of the essential amino acids. Therefore, these insects have built up symbiotic associations with bacteria that can synthesize them. Most of the scale insect families that have been analyzed, such as Monophlebidae, Coelostomidiidae, Orthezidae, Phenacoccinae from Pseudococcidae, Coccidae, Lecanodiaspididae, Diaspididae, and a clade of Eriococcidae, harbor flavobacteria as primary symbionts and enterobacteria as secondary symbionts. [3739]. It has been reported that the families from scale insects Dactylopiidae, some Eriococcidae, and Pseudococcinae from Pseudococcidae harbor different endosymbionts, which could indicate that they lost their flavobacteria and enterobacteria and acquired other endosymbionts [39]. Flavobacteria seem to be very ancient symbionts, perhaps starting symbiosis before the divergence of scale insects [39] (150–250 mya [40]). Although it has been suggested that Flavobacteria have cospeciated only within Monophlebidae, Coelostomidiidae, Ortheziidae, and Diaspididae [3841], and host switches seem to have occurred in the other families [39]. Otherwise, enterobacteria have undergone more evolutionary events (losses, duplications, and host switches). Some scale insects have enterobacteria closely related to Sodalis endosymbionts (Sodalis-like). But others may have symbionts closely related to Pantoea and Klebsiella [39].

Sodalis cosymbionts have been identified mainly by their 16S rRNA but also by other gene sequences. They have been found within various insect orders including Diptera, Coleoptera, Phthiraptera, and Hemiptera [4245]. The first described was S. glossinidius, the secondary symbiont of tsetse flies [46]. Later, bacteria with related gene sequences were referred as Sodalis-like [47] or Sodalis-affiliated but more recently several “Sodalis-like” bacteria and SOPE [48] are classified as Sodalis, others have been assigned to different genera. Still, scientists are in the process of making correct adscriptions for some of these bacteria [49].

The flavobacteria endosymbiont Candidatus Walczuchella monophlebidarum (“Walczuchella” from here on) was sequenced from the giant wax cochineal Llaveia axin axin (Llave) (Coccoidea: Monophlebidae) [50]. This insect has been used to obtain a lacquer to coat traditional art crafts by native people in Mexico and Guatemala since pre-Hispanic times [51]. The flavobacterial genome revealed that the endosymbiont’s major role is to synthesize and provide amino acids to the insect host [50]. The Flavobacteria genome was obtained from the analysis of a metagenome of L. axin axin. From this metagenome, we could also ensemble sequences from other microorganisms. Here, we present the draft genome of another cosymbiont of Walczuchella, a Sodalis-like bacteria that is designated here as Sodalis TME1. We also present a comparison to the genomes of five other Sodalis, as well as preliminary data of a metatranscriptome performed in the bacteriome of L. axin axin adults and in the ovaries of senescent adults.


2. Materials and methods

DNA, sequencing, and assembly were performed from bacteriomes (Illumina HiSeq 2000) and from the homogenized of female adults (pyrosequencing) of L. axin axin collected in the state of Chiapas, Mexico, as described [50]. A photograph from L. axin axin female adults is shown in Figure 1. RAST and GosthKOALA from KEGG [52] were used for genomic and metabolic pathway annotation of the metagenomic data that was previously reported when we obtained the Walczuchella genome [50]. Sodalis TME1 genome sequence has been deposited at DDBJ/ENA/GenBank under the accession MNBX00000000. The version described in this chapter is MNBX01000000.

Figure 1.

L. axin axin adult females on a Jatropha curcas plant.

Comparative phylogenomic analysis was performed with 20 genomes of gammaproteobacteria from GeneBank. Gene calling of all genomes was performed using GeneMark version 2.5 [53]. The pangenome and core genome from orthologous genes of all strains were obtained by GET_HOMOLOGUES version 2.0 software [54] with -A -c -t 0 -M -n 35 and -A -c -t 0 -G -n 35 parameters. We selected a set of 143 unique single-copy orthologous genes from the core genome. Translated coding sequences of each gene were concatenated using BioEdit Version 7.2.5 and aligned with Clustal Omega version 1.2.1 [55]. Prottest3 version 3.4.2 [56] was used to select the best amino acid substitution model using the AICc correction. The edited alignment contained 47,803 amino acid positions. Maximum likelihood phylogeny was performed by PhyML software version 3.1 [57] using the CpREV model with the Shimodaira–Hasegawa-like procedure for internal branch support [58]. The genome of Escherichia coli K-12 MG1655 was used as outgroup.

Comparative genomics was carried out with the following Sodalis genomes: S. glossinidius morsitans from tsetse fly, Sodalis-like endosymbiont from the blood-feeding lice Proechinophthirus fluctus (an obligate ectoparasite of fur seals), S. pierantonus SOPE from rice weevils Sitophilus oryzae, the free-living S. praecaptivus, and Sodalis-like symbiont of the meadow spittlebug Philaenus spumarius. Orthologous genes and the core genomes were obtained by GET_HOMOLOGUES as described above. Core genome matrix was parsed from GET_HOMOLOGUES result, using the script. Shared genes between Sodalis-like TME1 and all other strains were retrieved by parsing the core matrix using custom perl scripts. Annotation of each gene cluster was carried out by BLASTp 2.2.30+ [59] searches against Uniref100 database. Furthermore, average nucleotide identity (ANI) was determined for all Sodalis genomes described above using the ANIcalculator software described by Varghese et al. [60] with the default parameters.

RNA was extracted from the bacteriome of L. axin axin female adults and from the ovaries of senescent female adults that do not possess the structure of the bacteriomes (bacteriomes degrade in senescent adults) (Figure 2). Sequencing of cDNA was performed by SOLID technology. The sequences were mapped to the genomes of Walczuchella, Sodalis-like TME1, and two insect reference genomes, Drosophila melanogaster and to the aphid Acyrthosiphon pisum. Differentially expressed genes were identified by comparing expression values between samples and using Kal's Z-test of proportions [61]. Genes with a change in the expression more than twofold and a p-value of <0.01 in the Z-test were considered as differentially expressed genes.

Figure 2.

Dissected L. axin axin adult females used for the metatranscriptome analysis. (A) early stage and (B) late stage or senescent adults.

To determine the uric acid and uricase activity, L. axin axin adult females were individually dissected under sterile conditions. Guts including the Malpighian tubules were extracted and metabolic activities were detected as described [62].


3. Results

We found gene sequences of an enterobacterium (gammaproteobacterium) related to Sodalis in the metagenome of the wax cochineal L. axin axin [50]. The phylogeny with a set of 143 conserved genes shows that the enterobacterium of L. axin axin is closely related to other Sodalis-like endosymbionts, especially close to the free-living S. praecaptivus [63] (Figure 3). The small branches in the Sodalis group may indicate that they have recently diverged while the large differences found in genome sizes among these endosymbionts indicate that evolution may be occurring mainly by genome reduction when compared to the larger genome of the free-living Sodalis (Figure 3).

Figure 3.

Maximum likelihood phylogeny of sequenced enterobacterial endosymbionts performed with 143 conserved genes. Sodalis endosymbionts of plant feeding host: green; blood feeding host: red; free-living style: blue. * : Sodalis TME1 used in this study. Scale bar indicates 1 % estimated sequence divergence. SH-aLRT values > 50 are indicated.

TME1 was compared with the ANI (average nucleotide identity) metric to other Sodalis using the same core genome used in the phylogenomic analysis. TME1 showed ANIs well over 95% that is used to delineate species with S. pierantonius SOPE and S. praecaptivus HS1, but lower than 95% with S. glossinidius morsitans, Sodalis-like SPU, and Sodalis-like SPI-1 (Table 1). There was a good correlation of the ANI values obtained and phylogenetic positions that allowed the identification of three groups within Sodalis (Figure 3 and Table 1).

Sodalis str. TME1 Sodalis pierantonius SOPE Sodalis praecaptivus HS1 Sodalis glossinidius str. morsitans Sodalis-like str. PSPU Sodalis-like str. SPI-1
Sodalis str. TME1
Sodalis pierantonius str. SOPE 98.45
Sodalis praecaptivus str. HS1 98.54 98.46
Sodalis glossinidius str. morsitans 91.24 91.04 91.27
Sodalis-like str. PSPU 89.81 89.67 89.87 95.48
Sodalis-like str. SPI-1 92.25 91.94 92.15 89.91 85.86

Table 1.

Average nucleotide identity (ANI) percentage among Sodalis strains. Values in bold are >95%. Colors correspond to green, plant-feeding host; red, blood-feeding host; blue, free-living style.

The draft assembly of the enterobacterial endosymbiont Sodalis TME1 genome consisted of 679 scaffolds with an N50 of 7713 and an average G + C content of 55.6%. The scaffolds sum 3.4 Mb [50]. A total of 3067 genes were identified to which a functional annotation was assigned. The functional categories more represented by the annotated genes were catabolic and cellular process as well as carbohydrate, amino acid and transcription DNA dependent metabolism (Figure 4). Interestingly, many phage-related sequences were found as well as genes for different multidrug efflux pumps and type III and IV secretion systems. TME1 has genes for polyamine biosynthesis and excretion as well as Ankyrin repeat domains and for a lactoyl-glutathione lyase that is a detoxifying enzyme [64]. Among the conserved genes in the core genome of Sodalis TME1, S. pierantonius str. SOPE and S. praecaptivus str. HS1 are genes for the synthesis of flagella and for nitrate reduction (narGHI) and nitrite reduction (nfrABCD). Maybe nitrate serves in Sodalis as an electron acceptor in anaerobiosis as occurs in bacterial symbionts of marine bivalves Lucinoma aequizonata [65]. Sodalis TME1 genome has genes for uric acid utilization such as uricase (uaZ), allantoinase (allB), allantoate deiminase (allC), and urease (ureC and ureD). Comparative genomics with all Sodalis strains show that allC and the alpha subunit for urease gene (ureC) orthologous were only present in Sodalis TME1. Experimentally, uric acid and uricase activity were quantified in L. axin axin female adults. We detected 5.86 ± 0.77 ng of uric acid per tissue µg−1 and 32.87 ± 5. 25 mU of uricase per tissue µg−1 in female cochineals.

Figure 4.

Gene functional categories of Sodalis TME1.

We obtained 11,042,037 and 11,042,428 reads from the cDNA sequence of the bacteriome and the ovaries, respectively. These two organs were selected for studying the differentially expressed genes because endosymbionts are transferred from bacteriomes to the ovaries for vertical transmission to their offspring. It was expected to find genes related to the migration of the endosymbionts from the bacteriome and the colonization of the ovaries. Reads mapped to the reference genomes are shown in Table 2. The number of genes that were statistically differentially expressed is shown in Table 3.

Reference genomeBacteriomeOvaries
Drosophila melanogaster (exons)2,019,5852,008,381
Acyrthosiphon pisum (mRNA refseq)3,082,3192,912,207
Sodalis TME1409,128483,601

Table 2.

Number of reads mapped to the reference genomes.

Reference genomeBacteriomeOvaries
Drosophila melanogaster (exons)494680
Acyrthosiphon pisum (mRNA refseq)244280
Sodalis TME16650

Table 3.

Number of genes differentially expressed according to Z-test (p < 0.01).

Walczuchella in the bacteriome tissue showed only two genes that exhibited differential expression, a putative hydrolase and the chaperone GroEL. Other genes showed a change in expression less than twofold compared to their expression in the ovary. The chaperonin GroES is almost at the limit for differential expression with 1.86-fold (Table 4).

WalczuchellaSodalis TME1Insect
Bacteriome (differential expression) (high RPKM)Putative hydrolaseT3SS-secreted effectorChaperon Hsp70
Chaperones GroEL, GroESAllantoinaseABC transporters
Hypothetical proteinsHypothetical proteinsAntiparasitic-like peptide
ATP synthase B subunitNAD biosynthesisAsparaginase
Amino acids biosynthesis genesFtsE cell division geneUnknown genes
Transcriptional regulationExtracellular glutamate receptor channel
Flagellum synthesisPhospholipids synthesis
Transcriptional regulation
Ovary (differential expression) (high RPKM)ATP synthase B and A subunits (pseudogenes)Hypothetical proteinsATPase subunit
AhpC oxidative stress geneNAD biosynthesisTransmembrane transporters of sugars and amino acids
GlycoproteaseFlagellum synthesisPeptidoglycan-binding protein
Amino acids biosynthesis genesFtsE cell division geneLysozyme
Cytochrome c oxydaseGlycolysisUnknown genes
SecY translocasePhage lysozymeTranscriptional regulation
Hypothetical proteinsTranscriptional regulationPhospholipids synthesis

Table 4.

Highly expressed and differentially expressed genes in the bacteriome and the ovaries in the endosymbionts Walczuchella and Sodalis TME1 and the host L. axin axin.

From the ovary tissue, we found differential expression of Walczuchella genes that code for some ATP synthase subunits (some of them annotated previously as pseudogenes), cytochrome c oxidase, also some genes of protein translocation systems, tryptophan, histidine and chorismate biosynthesis, one gene related to oxidative stress, and a gene that encodes a possible component of an ABC transporter (Table 4).

In the bacteriome, the enterobacterium TME1 showed very strong overexpression of a gene that codes an effector protein possibly secreted by the type III secretion system (TTSS), expressed 66.8-fold compared to its expression in the ovaries. Also, a gene that codes an allantoinase that participates in uric acid metabolism is highly overexpressed in the bacteriome, showing a 50-fold change. Other genes with overexpression in the bacteriome are four ABC transporters, a peroxidase, the heme synthase, two genes related to nucleotides biosynthesis, two genes related to lipid A biosynthesis, and two genes of the type III secretion system (Table 4).

In the ovary, TME1-overexpressed genes were related to NAD synthesis, carbohydrate metabolism, stress response, and some transporters and transcriptional regulators (Table 4).

Among the insect differentially expressed genes in the bacteriome there were 19 putative transporters (for amino acids, carbohydrates, vitamins, drugs, or unknown substrates), five genes related to defense systems including an antiparasitic peptide with identity to Drosomycin, three from D. melanogaster, two genes related to heat-shock response, an oxidative stress response gene, seven genes related to amino acid metabolism, and some genes related to lipid, carbohydrate, and vitamin metabolism (Table 4).

On the other hand, we found that in the insect, in the ovaries there was overexpression of 15 transporters, 17 immune response genes, some genes related to heat shock, desiccation, oxidative stress, and hypoxia response, and genes related to lipids, vitamins, carbohydrates, nucleotides, amino acids, and chitin synthesis and metabolism (Table 4).


4. Discussion

Due to the annual cycle of the wax cochineal, we are only able to collect insects once a year during the rainy season. It is worth mentioning that in 2015 and 2016, we did not find cochineals in many of the places where we had collected previously. Considering the menace of mosquitoes transmitting Zika, or Chikungunya, extensive fumigations with chemical insecticides have been carried out in many places in Mexico, especially in Chiapas. The relation to the diminished populations of cochineals remains to be established.

A previous survey of symbiotic bacteria from scale insects in Mexico revealed the prevalence of Flavobacteria and Gammaproteobacteria [39]. Some of the Gammaproteobacteria had 16S ribosomal gene sequences closely related to those of TME1, and thus they may be considered as Sodalis as well. They were obtained from different scale insects such as Insignorthezia sp. and I. insignia, Icerya purchasi, Cripticerya sp., and Pseudococcus longispinus that together with Llaveia would be hosts for Sodalis.

While Flavobacteria and insects showed a co-divergent pattern of evolution, the phylogenetic relationships of the Gammaproteobacteria and insects were not parallel, indicating multiple enterobacterial transfers among the different hosts, and a more recent and less dependent symbiosis. In agreement, the genome size of the gammaproteobacterium TME1 is much larger than that from the primary endosymbiont from wax cochineals, the Flavobacteria Walczuchella, and also larger than those from other cosymbionts as the Betaproteobacteria that accompany the bacteroidete Sulcia found in some insects.

The genome from the gammaproteobacterium TME1 (3.4 Mb) is within the range of those from other Sodalis (1.4–4.7 Mb, Figure 3). There are very few genomes available from Sodalis, namely those from Sodalis found in blood-sucking insects as in lice [42] and tsetse flies [66], in plant-feeding insects as the rice weevils [44], in spittlebugs [45], and from a free-living bacterium [67]. The average nucleotide identity (ANI [68] being used for global genomic comparisons and considered now as a gold standard in prokaryote taxonomy [69]) was estimated for the Sodalis with available genomes. ANI values and the phylogenomic analysis performed showed Sodalis as a defined and coherent genus with three groups A–C. These groups could represent at least three different species according to the global standards [69]. Two of these groups were identified as different lineages by Lo et al. [49]. The phylogenetic groups that we described here have a 100 SH-like value support, group A contains S. glossinidius from tsetse flies and Sodalis from the meadow spittlebug P. spumarius, group B is constituted by Sodalis from the fur seal P. fluctus, and group C contains the closely related TME1, the free-living S. praecaptivus and S. pierantonius SOPE. The nucleotide sequence conservation among the group A symbiotic and free-living Sodalis may reflect that the former were recent acquisitions in insects without enough time for sequence divergence in their hosts. The presence of very similar Sodalis in distinct insect isolates reinforces the reports that indicate that they may frequently be transferred among hosts [39, 47].

TME1 has biosynthetic pathways for all essential amino acids and may supply the needs of the wax cochineal and of Walczuchella that does not have complete pathways for the biosynthesis of all essential amino acids. Since Sodalis TME1 has all enzymes for TCA it may complement this pathway in Walczuchella. It is worth noting that the flavobacterium Candidatus Uzinura diaspidicola, an endosymbiont from the armored scale insect Aphytis melinus that feed on parenchyma which may provide more nutrients than sap, supplies its host with all nutrients without the need of a cosymbiont [70]. Other armored scale insects have been reported to have a Sodalis-like endosymbiont [39].

In S. glossinidius that is a secondary symbiont of tsetse flies, a type III secretion system was found implicated in cell invasion and maybe required for colonizing the insect bacteriocytes [71]. Genes encoding for a similar system were found in TME1. Notably, genes that code for the type III secretion system (TTSS) as well as a gene coding for an effector protein that may be secreted by this system were among the most highly induced in the bacteriome of TME1. In Salmonella enterica, polyamines are required for full expression of TTSS and for some effector coding genes. Mutants in polyamine biosynthesis are affected in intracellular colonization and survival and may be complemented by adding polyamines to the medium [72]. Furthermore, the modulation of a TTSS by a spermidine transporter has been reported in Pseudomonas aeruginosa. Exogenous addition of spermidine to the wild P. aeruginosa strain increased the expression of genes that produce effector proteins [73]. TME1 has all genes for spermidine and putrescine biosynthesis as well as for the excretion of spermidine. Polyamines may regulate host defense responses as do some effectors secreted by TTSS. This remains to be tested.

Uric acid and uricase activity were detected in L. axin axin females. Uric acid is the final product of purine metabolism. Only few insects are capable of degrading uric acid into other products. In plant-feeding insects, bacterial and fungal symbionts are capable of recycling uric acid into other nitrogen sources [7476]. Sodalis TME1 has uricase and allantoinase-codifying genes, and the latter was highly expressed in bacteriomes suggesting that Sodalis TME1 could participate in providing nitrogen to the host by uric acid recycling.

By reverse transcriptase-polymerase chain reaction (RT-PCR) using primers targeted to Sodalis, we found sequences from Sodalis in the bacteriome (our own unpublished results), thus we may suppose that Sodalis are localized in bacteriomes as Walczuchella. In Llaveia, in addition to Walczuchella and Sodalis we found sequences of alphaproteobacteria that are related to Rickettsiales and several fungi that are reported elsewhere (Vera Ponce de León, submitted). Coincidently, the seal lice with a Sodalis endosymbiont also harbor a Rickettsia that is very abundant. The role of the very little abundant Rickettsia-like bacterium in Llaveia is unknown. Wolbachia is found in members of the Coelostomidiidae family [37] that is closely related to Monophlebidae insect family that contains the Mexican wax cochineals.

Here, we used the term symbiome [27] to refer to the group of primary and secondary (cosymbionts) endosymbionts (and/or their genomes), residing in a host. We consider that the term symbiome is more adequate than the terms endosymbiotic community or consortium that are sometimes used instead.

The cosymbionts of different Flavobacteria in scale insects are diverse lineages of related Gammaproteobacteria [39]. Similarly, the cosymbionts of Sulcia (a flavobacterium as Walczuchella) are varied and may be different even in related hosts [7, 36]. Sulcia cosymbionts may belong to alpha, beta, or gammaproteobacteria, with alpha and betaproteobacteria looking like the oldest symbionts. It was reported that Candidatus Zinderia insecticola, the Betaproteobacteria of spittelbugs was probably substituted by a Sodalis-like symbiont in members of the Philaenini tribe of the spittelbugs [33, 45]. The displacement of betaproteobacterial cosymbionts by the gammaproteobacterium Sodalis seems recent and was described as an event “in statu nascendi” (in the stage of being born) in Cicadella viridis [77]. There are other examples where one endosymbiont may substitute another one or is on the way toward displacement of a highly reduced-genome endosymbiont [33, 7779]. Distinct (apparently replaceable) cosymbionts may fulfill the different needs of insects that may change overtime and conditions specially if the insect changes habit [22], otherwise there may be cosymbiont redundancy, with different bacteria performing the same or very similar role (e.g., the synthesis of essential amino acids). The Sodalis cosymbiont in the wax cochineals seems to be recently acquired as in C. viridis. The insect symbiome seems plastic or dynamic with cosymbionts playing a key role in this plasticity. Here, we enlarged the list of putative functions of Sodalis that may include uric acid recycling, polyamine biosynthesis, or detoxification.



AVPL is a doctoral student from Programa de Doctorado en Ciencias Biomédicas, Universidad Nacional Autónoma de México (UNAM) and received the 331625 fellowship from Consejo Nacional de Ciencia y Tecnología (CONACyT).


  1. 1. McCutcheon JP. The bacterial essence of tiny symbiont genomes. Current Opinion in Microbiology. 2010;13:73–78. DOI: 10.1016/j.mib.2009.12.002.
  2. 2. Martínez-Cano DJ, Reyes-Prieto M, Martínez-Romero E, Partida-Martínez LP, Latorre A, Moya A, Delaye L. Evolution of small prokaryotic genomes. Frontiers in Microbiology. 2015; 5:742. DOI: 10.3389/fmicb.2014.00742.
  3. 3. Husník F, Chrudimský T, Hypša V. Multiple origins of endosymbiosis within the Enterobacteriaceae (γ-Proteobacteria): convergence of complex phylogenetic approaches. BMC Biology. 2011;9:87. DOI: 10.1186/1741-7007-9-87.
  4. 4. Ramírez-Puebla ST, Servín-Garcidueñas LE, Jiménez-Marín B, Bolaños LM, Rosenblueth M, Martínez J, Rogel MA, Ormeño-Orrillo E, Martínez-Romero E. Gut and root microbiota commonalities. Applied and Environmental Microbiology. 2013; 79:2–9. DOI: 10.1128/AEM.02553-12.
  5. 5. Trujillo ME, Alonso-Vega P, Rodríguez R, Carro L, Cerda E, Alonso P, Martínez-Molina E. The genus Micromonospora is widespread in legume root nodules: the example of Lupinus angustifolius. ISME Journal. 2010; 4:1265–1281. DOI: 10.1038/ismej.2010.55.
  6. 6. Trujillo ME, Riesco R, Benito P, Carro L. Endophytic Actinobacteria and the interaction of Micromonospora and nitrogen fixing plants. Frontiers in Microbiology. 2015; 6:1341. DOI: 10.3389/fmicb.2015.01341.
  7. 7. Noda H, Watanabe K, Kawai S, Yukuhiro F, Miyoshi T, Tomizawa M, Koizumi Y, Nikoh N, Fukatsu T. Bacteriome-associated endosymbionts of the green rice leafhopper Nephotettix cincticeps (Hemiptera: Cicadellidae). Applied Entomology and Zoology. 2012; 47:217–225. DOI:10.1007/s13355-012-0110-1
  8. 8. Bennett GM, Moran NA. Small, smaller, smallest: the origins and evolution of ancient dual symbioses in a Phloem-feeding insect. Genome Biology and Evolution. 2013; 5:1675–1688. DOI: 10.1093/gbe/evt118.
  9. 9. Tsuchida T, Koga R, Meng XY, Matsumoto T, Fukatsu T. Characterization of a facultative endosymbiotic bacterium of the pea aphid Acyrthosiphon pisum. Microbial Ecology. 2005; 49:126–133. DOI: 10.1007/s00248-004-0216-2.
  10. 10. Kaiwa N, Hosokawa T, Kikuchi Y, Nikoh N, Meng XY, Kimura N, Ito M, Fukatsu T. Primary gut symbiont and secondary, Sodalis-allied symbiont of the Scutellerid stinkbug Cantao ocellatus. Applied and Environmental Microbiology. 2010; 76:3486–3494. DOI: 10.1128/AEM.00421-10.
  11. 11. Chrudimský T, Husník F, Nováková E, Hypša V. Candidatus Sodalis melophagi sp. nov.: phylogenetically independent comparative model to the tsetse fly symbiont Sodalis glossinidius. PLoS One. 2012; 7:e40354. DOI: 10.1371/journal.pone.0040354.
  12. 12. Happ GM, Happ CM, Barras SJ. Fine structure of the prothoracic mycangium, a chamber for the culture of symbiotic fungi, in the southern pine beetle, Dendroctonus frontalis. Tissue and Cell. 1971; 3:295–308. DOI: 10.1016/S0040-8166(71)80024-1
  13. 13. Kuechler SM, Gibbs G, Burckhardt D, Dettner K, Hartung V. Diversity of bacterial endosymbionts and bacteria-host co-evolution in Gondwanan relict moss bugs (Hemiptera: Coleorrhyncha: Peloridiidae). Environmental Microbiology. 2013; 15:2031–2042. DOI: 10.1111/1462-2920.12101.
  14. 14. Wu D, Daugherty SC, Van Aken SE, Pai GH, Watkins KL, Khouri H, Tallon LJ, Zaborsky JM, Dunbar HE, Tran PL, Moran NA, Eisen JA. Metabolic complementarity and genomics of the dual bacterial symbiosis of sharpshooters. PLoS Biology. 2006; 4:e188. DOI: 10.1371/journal.pbio.0040188.
  15. 15. Cottret L, Milreu PV, Acuña V, Marchetti-Spaccamela A, Stougie L, Charles H, Sagot MF. Graph-based analysis of the metabolic exchanges between two co-resident intracellular symbionts, Baumannia cicadellinicola and Sulcia muelleri, with their insect host, Homalodisca coagulata. PLoS Computational Biology. 2010; 6:e1000904.
  16. 16. López-Madrigal S, Latorre A, Porcar M, Moya A, Gil R. Mealybugs nested endosymbiosis: going into the ‘matryoshka’ system in Planococcus citri in depth. BMC Microbiology 2013; 13:74. DOI: 10.1186/1471-2180-13-74
  17. 17. Rao Q, Rollat-Farnier PA, Zhu DT, Santos-Garcia D, Silva FJ, Moya A, Latorre A, Klein CC, Vavre F, Sagot MF, Liu SS, Mouton L, Wang XW. Genome reduction and potential metabolic complementation of the dual endosymbionts in the whitefly Bemisia tabaci. BMC Genomics. 2015; 16:226. DOI: 10.1186/s12864-015-1379-6.
  18. 18. Scarborough CL, Ferrari J, Godfray HC. Aphid protected from pathogen by endosymbiont. Science. 2005; 310:1781–1781. DOI: 10.1126/science.1120180.
  19. 19. Moran NA, Tran P, Gerardo NM. Symbiosis and insect diversification: an ancient symbiont of sap-feeding insects from the bacterial phylum Bacteroidetes. Applied Environmental Microbiology. 2005; 71:8802–8810. DOI:10.1128/AEM.71.12.8802-8810.2005.
  20. 20. Nakabachi A, Ueoka R, Oshima K, Teta R, Mangoni A, Gurgui M, Oldham NJ, van Echten-Deckert G, Okamura K, Yamamoto K, Inoue H, Ohkuma M, Hongoh Y, Miyagishima SY, Hattori M, Piel J, Fukatsu T. Defensive bacteriome symbiont with a drastically reduced genome. Current Biology. 2013; 23:1478–1484. DOI: 10.1016/j.cub.2013.06.027.
  21. 21. Leclair M, Pons I, Mahéo F, Morlière S, Simon J-C, Outreman Y. Diversity in symbiont consortia in the pea aphid complex is associated with large phenotypic variation in the insect host. Evolutionary Ecology. 2016; DOI:10.1007/s10682-016-9856-1 (published on line).
  22. 22. Oliver KM, Degnan PH, Burke GR, Moran NA. Facultative symbionts in aphids and the horizontal transfer of ecologically important traits. Annual Review of Entomology. 2010; 55:247–266. DOI 10.1146/annurev-ento-112408-085305.
  23. 23. Chen DQ, Montllor CB, Purcell A H. Fitness effects of two facultative endosymbiotic bacteria on the pea aphid, Acyrthosiphon pisum, and the blue alfalfa aphid A. kondoi. Entomologia Experimentalis et Applicata. 2000; 95:315–323. DOI: 10.1046/j.1570-7458.2000.00670.x
  24. 24. Russell JA, Moran NA. Costs and benefits of symbiont infection in aphids: variation among symbionts and across temperatures. Proceedings of the Royal Society B: Biological Sciences. 2006; 273:603–610. DOI:10.1098/rspb.2005.3348.
  25. 25. Montllor C B, Maxmen A, Purcell AH. Facultative bacterial endosymbionts benefit pea aphids Acyrthosiphon pisum under heat stress. Ecological Entomology. 2002; 27:189–195. DOI:10.1046/j.1365-2311.2002.00393.x
  26. 26. Teixeira L, Ferreira A, Ashburner M. The bacterial symbiont Wolbachia induces resistance to RNA viral infections in Drosophila melanogaster. PLoS Biology. 2008; 6:e2. DOI: 10.1371/journal.pbio.1000002
  27. 27. Bouclas DG, Kariithi HM, Bourtzis K, Schneider DI, Kelley K, Miller WJ, Parker AG, Abd-Alla AM. Transgenerational transmission of the Glossina pallidipes hytrosavirus depends on the presence of a functional symbiome. PLoS One. 2013; 8:e61150. DOI: 10.1371/journal.pone.0061150.
  28. 28. Chrostek E, Marialva MS, Esteves SS, Weinert LA, Martinez J, Jiggins FM, Teixeira L. Wolbachia variants induce differential protection to viruses in Drosophila melanogaster: a phenotypic and phylogenomic analysis. PLoS Genetics. 2013; 9:e1003896. DOI: 10.1371/journal.pgen.1003896
  29. 29. Werren JH, Baldo L, Clark ME. Wolbachia: master manipulators of invertebrate biology. Nature Reviews Microbiology. 2008; 6:741–751. DOI: 10.1038/nrmicro1969.
  30. 30. Koga R, Tsuchida T, Fukatsu T. Changing partners in an obligate symbiosis: a facultative endosymbiont can compensate for loss of the essential endosymbiont Buchnera in an aphid. Proceedings of the Royal Society B: Biological Sciences. 2003; 270:2543–2550. DOI: 10.1098/rspb.2003.2537
  31. 31. Hosokawa T, Kaiwa N, Matsuura Y, Kikuchi Y, Fukatsu T. Infection prevalence of Sodalis symbionts among stinkbugs. Zoological Letters. 2015; 1:5. DOI: 10.1186/s40851-014-0009-5.
  32. 32. McCutcheon JP, Moran NA. Functional convergence in reduced genomes of bacterial symbionts spanning 200 My of evolution. Genome Biology and Evolution. 2010: 2:708–718. DOI:10.1093/gbe/evq055
  33. 33. Koga R, Bennett GM, Cryan JR, Moran NA. Evolutionary replacement of obligate symbionts in an ancient and diverse insect lineage. Environmental Microbiology. 2013; 15:2073–2081. DOI: 10.1111/1462-2920.12121.
  34. 34. Van Leuven JT, Meister RC, Simon C, McCutcheon JP. Sympatric speciation in a bacterial endosymbiont results in two genomes with the functionality of one. Cell. 2014; 158:1270–1280. DOI: 10.1016/j.cell.2014.07.047.
  35. 35. Bennett GM, Abbà S, Kube M, Marzachì C. Complete genome sequences of the obligate symbionts "Candidatus Sulcia muelleri" and "Ca. Nasuia deltocephalinicola" from the Pestiferous Leafhopper Macrosteles quadripunctulatus (Hemiptera: Cicadellidae). Genome Announcements. 2016; 4 pii:e01604–15. DOI: 10.1128/genomeA.01604-15.
  36. 36. Kobiałka M, Michalik A, Walczak M, Junkiert Ł, Szklarzewicz T. Sulcia symbiont of the leafhopper Macrosteles laevis (Ribaut, 1927) (Insecta, Hemiptera, Cicadellidae: Deltocephalinae) harbors Arsenophonus bacteria. Protoplasma. 2016; 253:903–912. DOI: 10.1007/s00709-015-0854-x.
  37. 37. Dhami MK, Turner AP, Deines P, Beggs JR, Taylor MW. Ultrastructural and molecular characterization of a bacterial symbiosis in the ecologically important scale insect family Coelostomidiidae. FEMS Microbiology Ecology. 2012; 81:537–546. DOI: 10.1111/j.1574-6941.2012.01378.x
  38. 38. Gruwell ME, Morse GE, Normark BB. Phylogenetic congruence of armored scale insects (Hemiptera: Diaspididae) and their primary endosymbionts from the phylum Bacteroidetes. Molecular Phylogenetics and Evolution. 2007; 44:267–280. DOI: org/10.1016/j.ympev.2007.01.014
  39. 39. Rosenblueth M, Sayavedra L, Sámano‐Sánchez H, Roth A., Martínez‐Romero E. Evolutionary relationships of flavobacterial and enterobacterial endosymbionts with their scale insect hosts (Hemiptera: Coccoidea). Journal of Evolutionary Biology. 2012; 25:2357–2368. DOI: 10.1111/j.1420-9101.2012.02611.x
  40. 40. Koteja J. Essay on the prehistory of the scale insects (Homoptera, Coccinea). In Annales Zoologici. Państwowe Wydawnictwo Naukowe: Warzawa, Wroclaw; 1985; 38. p. 461–503.
  41. 41. Dhami MK, Buckley TR, Beggs JR, Taylor MW. Primary symbiont of the ancient scale insect family Coelostomidiidae exhibits strict cophylogenetic patterns. Symbiosis. 2013; 61:77–91. DOI: 10.1007/s13199-013-0257-8
  42. 42. Boyd BM, Allen JM, Koga R, Fukatsu T, Sweet AD, Johnson KP, Reed DL. Two bacterial genera, Sodalis and Rickettsia, associated with the seal louse Proechinophthirus fluctus (Phthiraptera: Anoplura). Applied and Environmental Microbiology. 2016; 82:3185–3197. DOI: 10.1128/AEM.00282-16.
  43. 43. Toju H, Koga R, Nikoh N, Meng XY, Kimura N, Fukatsu T. “Candidatus Curculioniphilus buchneri” a novel clade of bacterial endocellular symbionts from weevels of the genus Curculio. Applied and Environmental Microbiology. 2010; 76:275–282. DOI: 10.1128/AEM.02154-09
  44. 44. Oakeson KF, Gil R, Clayton AL, Dunn DM, von Niederhausern AC, Hamil C, Aoyagi A, Duval B, Baca A, Silva FJ, Vallier A, Jackson DG, Latorre A, Weiss RB, Heddi A, Moya A, Dale C. Genome degeneration and adaptation in a nascent stage of symbiosis. Genome Biology and Evolution. 2014; 6:76–93. DOI: 10.1093/gbe/evt210.
  45. 45. Koga R, Moran NA. Swapping symbionts in spittlebugs: evolutionary replacement of a reduced genome symbiont. ISME Journal. 2014; 8:1237–1246. DOI: 10.1038/ismej.2013.235.
  46. 46. Dale C, Maudlin I. Sodalis gen. nov. and Sodalis glossinidius sp. nov., a microaerophilic secondary endosymbiont of the tsetse fly Glossina morsitans morsitans. International Journal of Systematic Bacteriology. 1999; 49:267–275. DOI: 10.1099/00207713-49-1-267.
  47. 47. Smith WA, Oakeson KF, Johnson KP, Reed DL, Carter T, Smith KL, Koga R, Fukatsu T, Clayton DH, Dale C. Phylogenetic analysis of symbionts in feather-feeding lice of the genus Columbicola: evidence for repeated symbiont replacements. BMC Evolutionary Biology. 2013; 13:109. DOI: 10.1186/1471-2148-13-109.
  48. 48. Gil García R, Belda E, Gosalbes MJ, Delaye L, Vallier A, Vincent Monégat C, Heddi A, Silva GJ, Moya A, Latorre A. Massive presence of insertion sequences in the genome of SOPE, the primary endosymbiont of the rice weevil "Sitophilus oryzae". International Microbiology. 2008; 11:4–48. DOI: 10.2436/20.1501.01.43.
  49. 49. Lo WS, Huang YY, Kuo CH. Winding paths to simplicity: genome evolution in facultative insect symbionts. FEMS Microbiology Reviews. 2016; pii:fuw028. DOI: org/10.1093/femsre/fuw028
  50. 50. Rosas-Pérez T, Rosenblueth M, Rincón-Rosales R, Mora J, Martínez-Romero E. Genome sequence of “Candidatus Walczuchella monophlebidarum” the flavobacterial endosymbiont of Llaveia axin axin (Hemiptera: Coccoidea: Monophlebidae). Genome Biology and Evolution. 2014; 6:714–726. DOI: 10.1093/gbe/evu049
  51. 51. Williams M, MacVean CM. Ethnococcidology: use of the giant margarodids, Llaveia spp. (Homoptera: Coccoidea: Margarodidae), by indigenous peoples of Mesoamerica in their culture, medicine and arts. Israel Journal of Entomology. 1995; 29:147–148.
  52. 52. Kanehisa M, Sato Y, Morishima K. BlastKOALA and GhostKOALA: KEGG tools for functional characterization of genome and metagenome sequences. Journal of Molecular Biology. 2016; 428:726–731. DOI: 10.1016/j.jmb.2015.11.006
  53. 53. Besemer J, Borodovsky M. GeneMark: Web software for gene finding in prokaryotes, eukaryotes and viruses. Nucleic Acids Research. 2005; 33:451–454. DOI:10.1093/nar/gki487
  54. 54. Contreras-Moreira B, Vinuesa P. GET_HOMOLOGUES, a versatile software package for scalable and robust microbial pangenome analysis. Applied and Environmental Microbiology. 2013; 79:7696–7701. DOI:10.1128/AEM.02411-13
  55. 55. Sievers F, Wilm A, Dineen D, Gibson TJ, Karplus K, Li W, Higgins DG. Fast, scalable generation of high-quality protein multiple sequence alignments using Clustal Omega. Molecular Systems Biology. 2011; 7:539. DOI:10.1038/msb.2011.75
  56. 56. Darriba D, Taboada GL, Doallo R, Posada D. ProtTest 3: fast selection of best-fit models of protein evolution. Bioinformatics. 2011; 27:1164–1165. DOI:10.1093/bioinformatics/btr088
  57. 57. Guindon S, Dufayard J-F, Lefort V, Anisimova M, Hordijk W, Gascuel O. New algorithms and methods to estimate maximum-likelihood phylogenies: assessing the performance of PhyML 3.0. Systematic Biology. 2010; 59:307–321. DOI:10.1093/sysbio/syq010
  58. 58. Shimodaira H, Hasegawa M. Letter to the editor multiple comparisons of log-likelihoods with applications to phylogenetic inference. Molecular Biology and Evolution. 1999; 16:1114–1116. DOI:10.1177/0148607109348061
  59. 59. Camacho C, Coulouris G, Avagyan V, Ma N, Papadopoulos J, Bealer K, Madden TL. BLAST+: architecture and applications. BMC Bioinformatics. 2009; 10:421. DOI:10.1186/1471-2105-10-421
  60. 60. Varghese NJ, Mukherjee S, Ivanova N, Konstantinidis KT, Mavrommatis K, Kyrpides NC, Pati A. Microbial species delineation using whole genome sequences. Nucleic Acids Research. 2015; 43:6761–6771. DOI: 10.1093/nar/gkv657.
  61. 61. Kal AJ, van Zonneveld AJ, Benes V, van den Berg M, Koerkamp MG, Albermann K, et al. Dynamics of gene expression revealed by comparison of serial analysis of gene expression transcript profiles from yeast grown on two different carbon sources. Molecular Biology of the Cell. 1999; 10:1859–1872. DOI:10.1091/mbc.10.6.1859
  62. 62. Vera-Ponce de León A, Sanchez-Flores A, Rosenblueth M, Martínez-Romero E. Fungal community associated with Dactylopius (Hemiptera: Coccoidea: Dactylopiidae) and its role in uric acid metabolism. Frontiers in Microbiology. 2016; 7:954. DOI: 10.3389/fmicb.2016.00954
  63. 63. Chari A, Oakeson KF, Enomoto S, Jackson DG, Fisher MA, Dale C. Phenotypic characterization of Sodalis praecaptivus sp. nov., a close non-insect-associated member of the Sodalis-allied lineage of insect endosymbionts. International Journal of Systematic and Evolutionary Microbiology. 2015; 65:1400–1405. DOI: 10.1099/ijs.0.000091.
  64. 64. Thornalley PJ. The glyoxalase system: new developments towards functional characterization of a metabolic pathway fundamental to biological life. The Biochemical Journal. 1990; 269:1–11. DOI: 10.1042/bj2690001
  65. 65. Hentschel U, Hand SC, Felbeck H. The contribution of nitrate respiration to the energy budget of the symbiont-containing clam Lucinoma aequizonata: A calorimetric study. The Journal of Experimental Biology. 1996; 199:427–433.
  66. 66. Toh H, Weiss BL, Perkin SA, Yamashita A, Oshima K, Hattori M, Aksoy S. Massive genome erosion and functional adaptations provide insights into the symbiotic lifestyle of Sodalis glossinidius in the tsetse host. Genome Research. 2006; 16:149–156. DOI:10.1101/gr.4106106
  67. 67. Clayton AL, Oakeson KF, Gutin M, Pontes A, Dunn DM, von Niederhausern AC, Weiss RB, Fisher M, Dale C. A novel human-infection-derived bacterium provides insights into the evolutionary origins of mutualistic insect–bacterial symbioses. PLoS Genetics 2012; 8:e1002990. DOI: 10.1371/journal.pgen.1002990
  68. 68. Konstantinidis KT, Tiedje JM. Genomic insights that advance the species definition for prokaryotes. Proceedings of the National Academy of Sciences of the United States of America. 2005; 102:2567–2572. DOI:10.1073/pnas.0409727102.
  69. 69. Richter M, Rosselló-Móra R. Shifting the genomic gold standard for the prokaryotic species definition. Proceedings of the National Academy of Sciences of the United States of America. 2009; 106:19126–19131. DOI: 10.1073/pnas.0906412106.
  70. 70. Sabree ZL, Huang CY, Okusu A, Moran NA, Normark BB. The nutrient supplying capabilities of Uzinura, an endosymbiont of armoured scale insects. Environmental Microbiology. 2013; 15:1988–1999. DOI: 10.1111/1462-2920.12058.
  71. 71. Dale C, Young SA, Haydon DT, Welburn SC. The insect endosymbiont Sodalis glossinidius utilizes a type III secretion system for cell invasion. Proceedings of the National Academy of Sciences of the United States of America. 2001; 98:1883–1888. DOI:10.1073/pnas.98.4.1883
  72. 72. Jelsbak L, Thomsen LE, Wallrodt I, Jensen PR, Olsen JE. Polyamines are required for virulence in Salmonella enterica serovar Typhimurium. PLoS One. 2012; 7:e36149. DOI: 10.1371/journal.pone.0036149.
  73. 73. Zhou L, Wang J, Zhang LH. Modulation of bacterial Type III secretion system by a spermidine transporter dependent signaling pathway. PLoS One. 2007;2:e1291. DOI: 10.1371/journal.pone.0001291.
  74. 74. Potrikus CJ, Breznak JA. Anaerobic degradation of uric Acid by gut bacteria of termites. Applied and Environmental Microbiology. 1980; 40(1):125–132.
  75. 75. Morales-Jiménez J, Vera-Ponce de León A, García-Domínguez A, Martínez-Romero E, Zúñiga G, Hernández-Rodríguez C. Nitrogen-fixing and uricolytic bacteria associated with the gut of Dendroctonus rhizophagus and Dendroctonus valens (Curculionidae: Scolytinae). Microbial Ecology. 2013; 66:200–210. DOI:10.1007/s00248-013-0206-3
  76. 76. Patiño-Navarrete R, Piulachs M-D, Belles X, Moya A, Latorre A, Peretó J. The cockroach Blattella germanica obtains nitrogen from uric acid through a metabolic pathway shared with its bacterial endosymbiont. Biology Letters. 2014; 10:7–10. DOI:10.1098/rsbl.2014.0407.
  77. 77. Michalik A, Jankowska W, Kot M, Gołas A, Szklarzewicz T. Symbiosis in the green leafhopper, Cicadella viridis (Hemiptera, Cicadellidae). Association in statu nascendi? Arthropod Structure and Development. 2014; 43:579–587. DOI: 10.1016/j.asd.2014.07.005.
  78. 78. Lamelas A, Gosalbes MJ, Manzano-Marín A, Peretó J, Moya A, Latorre A. Serratia symbiotica from the aphid Cinara cedri: a missing link from facultative to obligate insect endosymbiont. PLoS Genetics. 2011; 7:e1002357. DOI: 10.1371/journal.pgen.1002357.
  79. 79. Lefèvre C, Charles H, Vallier A, Delobel B, Farrell B, Heddi A. Endosymbiont phylogenesis in the dryophthoridae weevils: evidence for bacterial replacement. Molecular Biology and Evolution. 2004; 21:965–973. DOI:10.1093/molbev/msh063

Written By

Tania Rosas-Pérez, Arturo Vera-Ponce de León, Mónica Rosenblueth, Shamayim T. Ramírez-Puebla, Reiner Rincón-Rosales, Julio Martínez-Romero, Michael F. Dunn, Éva Kondorosi and Esperanza Martínez-Romero

Reviewed: 20 October 2016 Published: 12 April 2017