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

Submitted: June 3rd, 2016 Reviewed: October 20th, 2016 Published: April 12th, 2017

DOI: 10.5772/66442

Chapter metrics overview

1,479 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 Micromonosporais found in nodules formed by Bradyrhizobium, Rhizobium, or Frankiain several legumes or actinorhizal roots, although Micromonosporais unable to form nodules [5]. Micromonosporahas 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 Buchnerasymbiont in an aphid [30]. Among others, alpha, gamma and betaproteobacteria have been found as cosymbionts; for example, the primary endosymbiont CandidatusSulcia 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 Sulciacosymbionts 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 Sulciacells 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 Sodalisendosymbionts (Sodalis-like). But others may have symbionts closely related to Pantoeaand Klebsiella[39].

Sodaliscosymbionts 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 CandidatusWalczuchella 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 SodalisTME1. 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 axincollected in the state of Chiapas, Mexico, as described [50]. A photograph from L. axin axinfemale 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 Walczuchellagenome [50]. SodalisTME1 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 axinadult females on aJatropha curcasplant.

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 coliK-12 MG1655 was used as outgroup.

Comparative genomics was carried out with the following Sodalisgenomes: S. glossinidiusmorsitans from tsetse fly, Sodalis-like endosymbiont from the blood-feeding lice Proechinophthirus fluctus(an obligate ectoparasite of fur seals), S. pierantonusSOPE 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 Sodalisgenomes 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 axinfemale 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 melanogasterand 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.

DissectedL. axin axinadult 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 axinadult 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 Sodalisin 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 axinis closely related to other Sodalis-like endosymbionts, especially close to the free-living S. praecaptivus[63] (Figure 3). The small branches in the Sodalisgroup 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.Sodalisendosymbionts of plant feeding host: green; blood feeding host: red; free-living style: blue. * :SodalisTME1 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 Sodalisusing the same core genome used in the phylogenomic analysis. TME1 showed ANIs well over 95% that is used to delineate species with S. pierantoniusSOPE and S. praecaptivusHS1, but lower than 95% with S. glossinidiusmorsitans, 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. TME1Sodalis pierantonius SOPESodalis praecaptivus HS1Sodalis glossinidius str. morsitansSodalis-like str. PSPUSodalis-like str. SPI-1
Sodalisstr. TME1
Sodalis pierantoniusstr. SOPE98.45
Sodalis praecaptivusstr. HS198.5498.46
Sodalis glossinidiusstr. morsitans91.2491.0491.27
Sodalis-like str. PSPU89.8189.6789.8795.48
Sodalis-like str. SPI-192.2591.9492.1589.9185.86

Table 1.

Average nucleotide identity (ANI) percentage among Sodalisstrains. 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 SodalisTME1 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 SodalisTME1, S. pierantoniusstr. SOPE and S. praecaptivusstr. HS1 are genes for the synthesis of flagella and for nitrate reduction (narGHI) and nitrite reduction (nfrABCD). Maybe nitrate serves in Sodalisas an electron acceptor in anaerobiosis as occurs in bacterial symbionts of marine bivalves Lucinoma aequizonata[65]. SodalisTME1 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 Sodalisstrains show that allC and the alpha subunit for urease gene (ureC) orthologous were only present in SodalisTME1. Experimentally, uric acid and uricase activity were quantified in L. axin axinfemale 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 ofSodalisTME1.

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

Table 2.

Number of reads mapped to the reference genomes.

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

Table 3.

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

Walczuchellain 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 Walczuchellaand SodalisTME1 and the host L. axin axin.

From the ovary tissue, we found differential expression of Walczuchellagenes 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 Sodalisas well. They were obtained from different scale insects such as Insignortheziasp. and I. insignia, Icerya purchasi, Cripticeryasp., and Pseudococcus longispinusthat together with Llaveiawould 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 Sulciafound 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 Sodalisfound 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 Sodaliswith available genomes. ANI values and the phylogenomic analysis performed showed Sodalisas 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. glossinidiusfrom tsetse flies and Sodalisfrom the meadow spittlebug P. spumarius, group B is constituted by Sodalisfrom the fur seal P. fluctus, and group C contains the closely related TME1, the free-living S. praecaptivusand S. pierantoniusSOPE. The nucleotide sequence conservation among the group A symbiotic and free-living Sodalismay reflect that the former were recent acquisitions in insects without enough time for sequence divergence in their hosts. The presence of very similar Sodalisin 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 Walczuchellathat does not have complete pathways for the biosynthesis of all essential amino acids. Since SodalisTME1 has all enzymes for TCA it may complement this pathway in Walczuchella. It is worth noting that the flavobacterium CandidatusUzinura diaspidicola, an endosymbiont from the armored scale insect Aphytis melinusthat 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. glossinidiusthat 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. aeruginosastrain 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 axinfemales. 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]. SodalisTME1 has uricase and allantoinase-codifying genes, and the latter was highly expressed in bacteriomes suggesting that SodalisTME1 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 Sodalisin the bacteriome (our own unpublished results), thus we may suppose that Sodalisare localized in bacteriomes as Walczuchella. In Llaveia, in addition to Walczuchellaand Sodaliswe 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 Sodalisendosymbiont also harbor a Rickettsiathat is very abundant. The role of the very little abundant Rickettsia-like bacterium in Llaveiais unknown. Wolbachiais 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]. Sulciacosymbionts may belong to alpha, beta, or gammaproteobacteria, with alpha and betaproteobacteria looking like the oldest symbionts. It was reported that CandidatusZinderia 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 Sodalisseems 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 Sodaliscosymbiont 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 Sodalisthat 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 genusMicromonosporais widespread in legume root nodules: the example ofLupinus 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 ofMicromonosporaand 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 leafhopperNephotettix 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 aphidAcyrthosiphon 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 stinkbugCantao 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 melophagisp. nov.: phylogenetically independent comparative model to the tsetse fly symbiontSodalis 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 cicadellinicolaandSulcia 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 inPlanococcus citriin 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 whiteflyBemisia 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 aphidA. 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 aphidsAcyrthosiphon pisumunder 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 symbiontWolbachiainduces resistance to RNA viral infections inDrosophila 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 theGlossina pallidipeshytrosavirus 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.Wolbachiavariants induce differential protection to viruses inDrosophila 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 ofSodalissymbionts 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 "CandidatusSulcia muelleri" and "Ca. Nasuia deltocephalinicola" from the Pestiferous LeafhopperMacrosteles quadripunctulatus(Hemiptera: Cicadellidae). Genome Announcements. 2016;4pii: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,SodalisandRickettsia, associated with the seal louseProechinophthirus 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. “CandidatusCurculioniphilus buchneri” a novel clade of bacterial endocellular symbionts from weevels of the genusCurculio. 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. andSodalis glossinidius sp. nov., a microaerophilic secondary endosymbiont of the tsetse flyGlossina 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 genusColumbicola: 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 “CandidatusWalczuchella monophlebidarum” the flavobacterial endosymbiont ofLlaveia 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,Llaveiaspp. (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 withDactylopius(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 ofSodalis praecaptivus sp. nov., a close non-insect-associated member of theSodalis-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 clamLucinoma 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 ofSodalis glossinidiusin 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 ofUzinura, 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 endosymbiontSodalis glossinidiusutilizes 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 inSalmonella entericaserovar 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 ofDendroctonus rhizophagusandDendroctonus 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 symbioticafrom the aphidCinara 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

Submitted: June 3rd, 2016 Reviewed: October 20th, 2016 Published: April 12th, 2017