Mealybugs nested endosymbiosis: going into the ‘matryoshka’ system in Planococcus citri in depth
© López-Madrigal et al.; licensee BioMed Central Ltd. 2013
Received: 17 January 2013
Accepted: 25 March 2013
Published: 1 April 2013
In all branches of life there are plenty of symbiotic associations. Insects are particularly well suited to establishing intracellular symbiosis with bacteria, providing them with metabolic capabilities they lack. Essential primary endosymbionts can coexist with facultative secondary symbionts which can, eventually, establish metabolic complementation with the primary endosymbiont, becoming a co-primary. Usually, both endosymbionts maintain their cellular identity. An exception is the endosymbiosis found in mealybugs of the subfamily Pseudoccinae, such as Planococcus citri, with Moranella endobia located inside Tremblaya princeps.
We report the genome sequencing of M. endobia str. PCVAL and the comparative genomic analyses of the genomes of strains PCVAL and PCIT of both consortium partners. A comprehensive analysis of their functional capabilities and interactions reveals their functional coupling, with many cases of metabolic and informational complementation. Using comparative genomics, we confirm that both genomes have undergone a reductive evolution, although with some unusual genomic features as a consequence of coevolving in an exceptional compartmentalized organization.
M. endobia seems to be responsible for the biosynthesis of most cellular components and energy provision, and controls most informational processes for the consortium, while T. princeps appears to be a mere factory for amino acid synthesis, and translating proteins, using the precursors provided by M. endobia. In this scenario, we propose that both entities should be considered part of a composite organism whose compartmentalized scheme (somehow) resembles a eukaryotic cell.
KeywordsNested endosymbiosis Planococcus citri Moranella endobia Tremblaya princeps functional complementation
Symbiosis is a widespread natural phenomenon that has been postulated as one of the main sources of evolutionary innovation [1, 2], and it is an example of compositional evolution involving the combination of systems of independent genetic material . Many insects have established mutualistic symbiotic relationships, particularly with intracellular bacteria that inhabit specialized cells of the animal host (bacteriocytes). In most insect-bacteria endosymbioses described to date, host insects have unbalanced diets, poor in essential nutrients that are supplemented by their endosymbionts. Attending to their dispensability for host survival, we distinguish between primary (P) or obligate, and secondary (S) or facultative endosymbionts. P-endosymbionts are essential for host fitness and reproduction, and maternally transmitted through generations, while S-symbionts are not essential and can experience horizontal transfer. The genomes of P-endosymbionts usually exhibit an increase in their A + T content and undergo great size reduction, among other changes. The main evolutionary forces accounting for these features are relaxation of purifying selection on genes rendered unnecessary in the enriched intracellular environment, and random genetic drift due to a strong population bottlenecking throughout intergenerational transmission of the bacteria . P and S symbionts can coexist in the same host. When an S-symbiont is also present, the irreversible genomic degenerative process could lead to the loss of some P-endosymbiont metabolic capabilities needed by the host. In this situation, two outcomes are possible: the host insect can recruit those functions from the S-symbiont, which then becomes a co-primary endosymbiont, establishing metabolic complementation with the former P-endosymbiont to fulfill the host needs or [5–8]; alternatively, the S-symbiont may replace its neighbor .
Mealybugs (Hemiptera: Sternorrhyncha: Pseudoccidae) form one of the largest families of scale insects, including many agricultural pest species that cause direct crops damage or vector plant diseases while feeding on sap . All mealybug species analyzed so far possess P-endosymbionts. Two subfamilies have been identified, Phenacoccinae and Pseudococcinae , the latter having been studied in greater depth, all of which live in symbiosis with the β-proteobacterium “Candidatus Tremblaya princeps” (T. princeps from now on, for the sake of simplicity). Universal presence, along with the cocladogenesis of endosymbionts and host insects, led to T. princeps being considered the mealybug P-endosymbiont . However, recently, other P-endosymbionts from the β-proteobacteria and Bacteroidetes groups have been identified in the subfamily Phenacoccinae . Most genera of the subfamily Pseudococcinae also harbor additional γ-proteobacteria endosymbionts that, due to their discontinuous presence and polyphyletic origin, have been considered as S-symbionts . An unprecedented structural organization of the endosymbionts of the citrus mealybug Planococcus citri was revealed by von Dohlen and coworkers : each T. princeps cell harbors several S-endosymbiont cells, being the first known case of prokaryote-prokaryote endocelullar symbiosis. The S-endosymbiont has recently been named “Candidatus Moranella endobia” (M. endobia from now on) . The dynamics of both endosymbiont populations throughout the insect life-cycle and their differential behavior depending on host sex  suggest that both play an important role in their hosts’ nutritional and reproductive physiology, putting into question the secondary role of M. endobia.
The sequencing of two fragments of the genome of T. princeps from the pineapple mealybug, Dysmicoccus brevipes, showed a set of unexpected genomic features compared with that found in most P-endosymbiont reduced genomes. This species presents a rather high genomic G + C content – a rare condition among P-endosymbionts with the only known exception being “Candidatus Hodgkinia cicadicola” (P-endosymbiont of the cicada Diceroprocta semicincta) –, a partial genomic duplication including the ribosomal operon and neighbor genes, and low gene density. All other sequenced genomes from endosymbionts having a long relationship with their host maintain a single set of rRNA genes, therefore these data suggested an unprecedented complexity for this P-endosymbiont genome, an unexpected finding for a long co-evolutionary process, as already elucidated for this symbiotic system . However, the recent sequencing of two strains of T. princeps from P. citri (PCIT and PCVAL) has shown that it is, in fact, the smallest (139 kb) and most simplified bacterial genome described to date [16, 19]. Functional analysis reveals that the genetic repertoire of T. princeps is unable to sustain cellular life, according to Gil et al. (2004) , and that it entirely depends on M. endobia for many essential functions. Even though most of its genome is occupied by ribosomal genes and genes involved in the biosynthesis of essential amino acids, T. princeps likely depends on its symbiotic consortium partner to build its own ribosomes and for amino acid production [16, 19].
The work published by McCutcheon and von Dohlen  mainly focused on the analysis of the T. princeps genome and detangling the amino acid biosynthetic pathways in which all three partners (T. princeps, M. endobia and the host) appear to be involved. However, the characteristics and functionality of the M. endobia genome, as well as other possible modes of complementation between the two endosymbionts, have remained largely unexplored. In this work we present a comprehensive analysis of the predicted consortium functional capabilities and interactions, thus offering new insights into how this bacterial consortium may function internally. Additionally, we have performed a comparative analysis of both endosymbiont genomes in two P. citri strains, PCIT  and PCVAL ( and this work). Our analysis suggests that both genomes have undergone reductive evolution, albeit with some unusual genomic features, probably as a consequence of their unprecedented compartmentalized organization.
Results and discussion
Main features and genomic variability between two strains of P. citrinested endosymbionts
Main genomic features of the two strains of the P. citri endosymbiotic consortium already sequenced
GenBank accession number
Genome size (bp)
Total gene number
Small RNA genes
19 (CDS) 6 (tRNA)
19 (CDS) 4 (tRNA)
Overall gene density (%)
Average ORF length (bp)
Average IGRs (bp)
G + C content (%)
Both consortium partners lack a canonical oriC, which is consistent with the absence of dnaA, similarly to many other reduced endosymbiont genomes already sequenced (e.g., Blochmannia floridanus, Wigglesworthia glossinidia, Carsonella rudii, Hodgkinia cidadicola, Zinderia insecticola, and Sulcia muelleri). This has been considered an indication that the endosymbionts rely on their host for the control of their own replication . Another shared genomic characteristic of both endosymbionts is their low gene density (already noticed in  for T. princeps) and the large average length of the intergenic regions, in which no traces of homology with coding regions of other bacteria can be found. Although these traits are unusual in bacterial endosymbionts, they have also been described for Serratia symbiotica SCc, the co-primary endosymbiont of Buchnera aphidicola in the aphid Cinara cedri. This non-coding DNA is probably the remnant of ancient pseudogenes that are gradually being eroded .
The reductive process affecting both genomes has led to the loss of most regulatory functions. Thus, they lack most regulatory genes and some genes have lost regulatory domains. This is the case of metL and adk from T. princeps, which have lost the regulatory ‘ATC’ domain, or the loss of the ‘HTH’ domain of birA, the ‘PNPase C’ domain of rne and the ‘DEAD box A’ of dead in the case of M. endobia. Additionally, many other genes have been shortened due to frameshifts or the presence of premature stop codons, in comparison with their orthologs in free-living relatives (e.g. sspB, rplQ, rplO and aroC in T. princeps; thiC, ybgI, yacG, ygbQ, ftsL, ftsY and tilS in M. endobia). In some cases, the shortening removes some non-essential protein domains completely (e.g., engA, rpoA and rpoD in T. princeps; secA, aceF, yebA and metG in M. endobia). The loss of the ‘anticodon binding domain of tRNA’ and ‘putative tRNA binding domain’ of metG, encoding methionyl-tRNA synthetase is common to other endosymbionts with reduced genomes.
T. princeps genome comparison
The genome alignment of both T. princeps strains showed a high degree of identity at the sequence level (99.98%, being 138,903 bp identical), which is coherent with their evolutionary proximity and extreme genome reduction. Although we also detected the 7,032-bp region flanked by 71-bp inverted repeats described in the strain PCIT , we only found one orientation in the population used for genome sequencing.
The genome of strain PCVAL only differs in 4 nucleotides in length from strain PCIT , involving five short indel events of one (4 cases) or two nucleotides (1 case). Additionally, 23 nucleotide substitutions were detected. Transitions represent 43.5% (10/23) of the total substitutions. Although the number of mutations is too small to be representative and, therefore, it is difficult to draw clear conclusions, it is noteworthy that all indels plus 87% of the detected substitutions between both strains are located in the coding fraction of the genome, in spite of its low coding density. One of the detected indels affects the start codon of aroC, involved in the biosynthetic pathway of aromatic amino acids, which is then changed to a GTG start codon. Two other short deletions yield the loss (AT) and recovery (T) of the reading frame of ilvD, needed for the synthesis of isoleucine and valine. The non-inactivating character of these mutations on genes involved in biosynthetic pathways of essential amino acids without an ortholog in the genome of M. endobia, corroborates their importance for the bacterial partnership. The other two indels, as well as 20 out of 23 of the observed substitutions, were located at the 3′ end of rplQ, which suggests that this region could be a mutational hot-spot. To confirm this point, we analyzed the original P. citri DNA samples used in the genome sequencing experiments by PCR amplification of the rplQ and flanking ITS regions, as well as new DNA samples obtained from individual insects cultivated in Almassora (Spain) and from environmental colonies collected in Murcia (Spain). Although all three samples were obtained from different plant hosts and separated by more than 300 Km, they were identical. Since we have no direct availability of the PCIT strain, it is feasible that the Spanish and American populations differ.
M. endobia genomes comparison
The alignment of both genomes of M. endobia showed that the genome of strain PCVAL is 65 nucleotides shorter than that of PCIT, and allowed the identification of 262 substitutions. Among them, 90.1% were G/C↔A/T changes, with only 18 A↔T changes and 8 G↔C changes, which is additional indirect evidence of the mutational bias towards A/T already observed in the codon usage analysis (Additional file 2). As expected for a neutral process, the mutational bias affected both strains equally, being the changes G/C↔A/T evenly distributed (50.4% A/T in strain PCIT and 49.5% in PCVAL). Regarding the genome distribution of the polymorphisms, 47% of them (123) map onto IGRs, and 4.5% (12) onto 10 pseudogenes. The 139 substitutions detected in the coding fraction affect only 111 out of the 406 orthologous genes. Among these substitutions, 77 are synonymous (dS = 0.0011 ± 0,0001), and 62 non-synonymous (dN = 0.0005 ± 0,0000), with a ω = 0.44, suggesting the action of purifying selection. It is worth noticing that about 75% of them affect functional domains, suggesting that many putatively functional genes accumulate mutations, which also justifies the maintenance of a minimal set of molecular chaperones to help in the proper folding of the encoded proteins.
Additionally, 60 indels were detected between both M. endobia strains, with a mean size of 5.4 nucleotides, although there is a great variance, between 1 and 75 nucleotides. Results showed 58.3% (35/60) of the indels affect homopolymers of A (22/39), T (12/36) and, less frequently, G (5/37) and C (3/35), which is consistent with the higher proportion of A and T homopolymers. This fact may be related with the above-mentioned A/T mutational bias. Although artifacts due to sequencing errors cannot be ruled out, given that PCVAL genomes were assembled based on 454 sequencing data, there are several pieces of evidence that indicate that the observed indels may be real. First, although homopolymers can be found both in coding and non-coding regions, most indels affect the non-coding parts of the genome. Second, even when A/T homopolymers are quite abundant in the M. endobia genome (844 cases equal to or bigger than 6 nucleotides), only a small fraction of them are affected by indels (29 cases, representing 3.4%). Finally, the coverage of the affected regions was always higher than 27X, and the PCVAL reads polymorphism was almost null. The remaining indels affect microsatellites of 2 to 8 nucleotides with a small number of copies. Forty-seven indels (78.3%) map onto intergenic regions, pseudogenes (2 in ΨpdxB, 1 in ΨprfC) or the non-functional part of shortened genes (dnaX), and only 13 indels (21.7%) map onto coding regions. Most of these are located on the 3′ end of the affected gene, causing enlargement or shortening of the ORFs compared with the orthologous gene in other γ-proteobacteria. Thus, glyQ (involved in translation) and ptsI (participating in the incorporation of sugars to the intermediary metabolism) are enlarged in strain PCVAL, while rppH (involved in RNA catabolism) is shortened in this strain without affecting described functional domains. Conversely, the shortening of fis (encoding a bacterial regulatory protein) in PCVAL, and of yicC (unknown function) and panC (involved in the metabolism of cofactors and vitamins, a function that is incomplete in M. endobia) in PCIT, affect some functional domains, although their activity might not be compromised. Finally, amino acid losses without frameshift were observed in PCVAL (relative to PCIT) for the loci holC (encoding subunit chi of DNA polymerase III), rluB (involved in ribosome maturation), surA (encoding a chaperone involved in proper folding of external membrane proteins), and pitA (encoding an inorganic phosphate transporter). None of the corresponding functional domains were affected in the first two cases, while the indel polymorphisms mapped inside the ‘PPIC-type PPIASE’ domain in surA, which appears to be dispensable for the chaperone qualities of the protein . Therefore, it seems that most (if not all) changes that could affect the functions of the encoded proteins have been removed by the action of purifying selection.
Functional analysis of the nested consortium
Most endosymbiotic systems analyzed to date at the genomic level have a nutritional basis, and many of them involve the biosynthesis of essential amino acids that are in short supply in the host diet. The metabolic pathways leading to amino acid biosynthesis in the T. princeps-M. endobia consortium found in P. citri were recently analyzed in detail by McCutcheon and von Dohlen  and, therefore, they will not be dealt with in this study. These authors also stated that T. princeps is unable to perform DNA replication, recombination or repair by itself, and the same applies to translation. They speculate that a passive mechanism such as cell lysis could provide T. princeps with the needed gene products from M. endobia. Our present work provides a detailed analysis of the M. endobia functional capabilities, based on a functional analysis of its genome, regarding informational functions or other intermediate metabolism pathways beyond amino acids biosynthesis. In the following sections these functional capabilities will be analyzed in a comprehensive manner, considering both endosymbiotic partners, in order to identify putative additional levels of complementation between them.
DNA repair and recombination
Contrary to what is found in bacterial endosymbionts with similarly reduced genomes, M. endobia has quite a complete set of genes for DNA repair and recombination, while none were annotated in the T. princeps genome [16, 19]. Although it has lost the nucleotide excision repair genes (only uvrD is present), M. endobia retains a base excision repair system (the DNA glycosylases encoded by mutM and ung plus xth, the gene encoding exonuclease III, involved in the repair of sites where damaged bases have been removed). The mismatch repair system is also almost complete, since only mutH, encoding the endonuclease needed in this process to cleave the unmethylated strand, has been lost. Additionally, M. endobia also retains almost the entire molecular machinery for homologous recombination (recABCGJ, ruvABC, priAB), which could be responsible for the concerted evolution of the duplications in both genomes. In the absence of recD, the RecBC enzyme can still promote recombination, since it retains helicase and RecA loading activity. The missing exonuclease V activity can be replaced by other exonucleases with ssDNA degradation activity in the 5′ → 3′ sense, such of RecJ , which has been preserved. The final step in homologous recombination requires the reloading of origin-independent replication machinery. Two replisome reloading systems have been described in E. coli, one of which requires the participation of PriA, PriB and DnaT , and it appears that helicase DnaB loading and unwinding of a replication fork is dependent upon the activities of DnaT and DnaC, among other restart proteins. These last two proteins are the only two elements of the replisome that are not encoded in the M. endobia genome. However, mutations in dnaC which have the ability to bypass such requirements in the loading of DnaB have been described , and dnaC is also absent in other reduced genomes that have been characterized (e.g. Blochmannia floridanus, Wigglesworthia glossinidia or Mycoplasma genitalium). Additionally, the role of DnaT in primosome assembly has not been fully elucidated . Therefore, it cannot be ruled out that dnaT is not essential for the functioning of the homologous recombination system in this bacterial consortium.
Even though most genes present in the T. princeps genome are involved in RNA metabolism (78 out of 116 genes, occupying 35% of its genome length and 49% of its coding capacity) [16, 19], it still seems to depend on M. endobia for transcription and translation. Thus, T. princeps encodes every essential subunit of the core RNA polymerase (rpoBCA) and a single sigma factor (rpoD), but no other genes involved in the basic transcription machinery or in RNA processing and degradation are present in its genome. On the other hand, M. endobia possesses a minimal but yet complete transcription machinery  plus some additional genes, including the ones encoding the ω subunit of the RNA polymerase (rpoZ), the sigma-32 factor (rpoH), and the transcription factor Rho. It also presents several genes involved in the processing and degradation of functional RNAs, i.e. pnp, rnc (processing of rRNA and regulatory antisense RNAs), hfq (RNA chaperone), rne, orn, rnr (rRNA maturation and mRNA regulation in stationary phase), and rppH (mRNA degradation). It is surprising that the small genome of M. endobia has also retained several transcriptional regulators, the functions of which are not yet fully understood, and which are absent in other endosymbionts with reduced genomes. These include CspB and CspC (predicted DNA-binding transcriptional regulators under stress conditions), and NusB, which is required in E. coli for proper transcription of rRNA genes, avoiding premature termination . cpxR, encoding the cytoplasmic response regulator of the two-component signal transduction system Cpx, the stress response system that mediates adaptation to envelope protein misfolding , is also preserved, while the companion sensor kinase cpxA appears to be a pseudogene. This might be an indication of a constitutive activation of the regulatory protein.
Finally, as it was already stated, ribosomes are the best preserved molecular machinery in T. princeps[16, 19]. In addition to two copies of the ribosomal 23S-16S operon, it encodes 49 out of 56 ribosomal proteins needed to make a complete ribosome. On the other hand, M. endobia has also retained a full set of ribosomal proteins and also presents two copies of the 23S and 5S rRNA genes. The high redundancy of rRNA and ribosomal protein genes might indicate that ribosomes from both members of the consortium are not exchangeable, or that redundancy is needed to achieve proper levels of ribosomal components for cell functioning. Both genomes encode the tmRNA, a molecule needed to solve problems that arise during translation while only M. endobia encodes ribosome maturation proteins and translational factors.
Protein processing, folding and secretion
As compared with their orthologs in free-living relatives, both endosymbionts have retained at least a minimal set of chaperones  required for the proper folding of functional proteins in both members of the consortium. This is consistent with the presence of proteins accumulating non-synonymous substitutions. Some proteins can also be exported across the inner and outer membranes via typical gram-negative secretion systems (reviewed in ) encoded exclusively in the M. endobia genome. As other endosymbionts with similarly reduced genomes, M. endobia has retained a fully functional Sec translocation complex . It also encodes Ffh, which together with 4.5S RNA forms the signal recognition particle (SRP), needed to bind the signal sequence of the proteins targeted for secretion through this system and to drive them to FtsY, the SRP receptor. Although in other endosymbionts there is an alternative system to assist proteins in their secretion, in which the proteins are recognized by the SecB chaperone after translation, this system cannot be functional in this consortium, because secB appears to be a pseudogene .
T. princeps has almost null metabolic capacities, except for the production of essential amino acids, as described elsewhere . Only M. endobia encodes a phosphotransferase system (PTS) for the uptake of hexose as carbon source, and it is predicted to perform glycolysis, transform pyruvate into acetate, and use it to feed the pathway for fatty acids biosynthesis, similarly to that described for B. aphidicola BCc, with highly reduced metabolic capabilities . However, the pentose phosphate pathway appears to be incomplete, since only zwf, pgl and tkt have been preserved, while talA appears to be a pseudogene. Interestingly, T. princeps has retained a transaldolase TalB, which along with transketolase (Tkt) creates a reversible link between the pentose phosphate pathway and glycolysis, revealing another possible case of metabolic complementation between both bacteria.
Regarding the tricarboxylic acid (TCA) cycle, only mdh (encoding malate dehydrogenase) has been preserved in T. princeps, while M. endobia has retained only the genes that encode succinyl-CoA synthetase. This is the only step that has been maintained in S. symbiotica SCc , where the authors indicate that it must have been retained because it is necessary for lysine biosynthesis. Nevertheless, this cannot be the case in this consortium, since lysine biosynthesis cannot be accomplished.
As in other endosymbionts, NAD+ can be regenerated by the action of the NADH-quinone oxidoreductase encoded by the nuo operon. But, in the absence of ATP synthase coupled to the electron transport chain, the whole consortium relies on substrate-level phosphorylation as a source of ATP. Acetyl-CoA can also be a source of ATP thanks to the presence of the genes ackA and pta.
The consortium also shares with other endosymbiotic bacteria with reduced genomes the incapability to synthesize nucleotides de novo. T. princeps has completely lost all genes involved in this function, while M. endobia retains a metabolic capacity similar to B. aphidicola BCc . All triphosphate nucleotides could be obtained by phosphorylation from diphosphate nucleotides via pyruvate kinase A (pykA), while deoxynucleotides could be obtained via ribonucleoside diphosphate reductase 1 (whose subunits are encoded by nrdA and nrdB). The only preserved diphosphate kinase is adenylate kinase (adk), while cytidylate kinase appears to be a pseudogene. Although it has been described that at least one purine and one pyrimidine kinase are needed to phosphorylate all dinucleotides, the fact that Adk is the same kinase that has been preserved in B. aphidicola BCc might be an indication that, in endosymbiotic bacteria, this enzyme can act on both nucleotide types. The presence of dut guarantees that the thymidylate nucleotides can also be synthesized using dUTP as a primary source.
The endosymbiotic system has almost completely lost the ability to synthesize vitamins and cofactors. Yet, the importance of the [Fe-S] clusters in this consortium is revealed by the presence of complete machinery for the assembly of such components, a complex system that is not fully preserved in other reduced genomes of endosymbiotic bacteria. The [Fe-S] clusters are one of the most ubiquitous and functionally versatile prosthetic groups in nature . Although it is known that these clusters can spontaneously be assembled from the required components under the proper conditions, it is not an efficient procedure in vivo. In E. coli, their assembly requires a complex machinery and it is achieved by two sets of proteins, the Suf (sufABCDSE) and the Isc (iscSUA) systems. Both members of the consortium are involved in the maintenance of this machinery, revealing another possible case of metabolic complementation. The complete suf operon is present in the genome of M. endobia. Regarding the Isc system, both partners of the consortium retain iscS, and T. princeps also encodes iscU, while they both lack iscA. However, IscA belongs to the HesB family of proteins, and a hesB gene has been identified in T. princeps. Additionally, ErpA, an A-type iron-sulfur protein that can bind both [2Fe-2S] and [4Fe-4S] clusters, is present in M. endobia.
The cell envelope structure is usually highly simplified in Gram-negative endosymbiotic bacteria, which lack most (if not all) of the genes needed for the biosynthesis of murein and lipopolysaccharides, and these two bacteria are not an exception. In fact, T. princeps has lost all the genes involved in these functions, and M. endobia has also lost many pathways, although it still retains some peptidoglycan synthetases and hydrolases needed for septum formation during cell division. It is noteworthy that this is the first analyzed case of an endosymbiont with a highly reduced genome that retains the ability to synthesize lipid IVA, the biosynthetic precursor of lipopolysaccharydes.
Only M. endobia has preserved genes related to cellular transport, which must ensure proper exchange of metabolites with the host cell and between both endosymbionts. Many nutrients pass the outer membrane of Gram-negative bacteria via a family of integral outer-membrane proteins (OMPs). The only OMP encoded in the consortium genomes is OmpF, the protein that forms osmotically regulated pores for the passage of small solutes such as sugars, ions and amino acids, with a preference for cationic molecules. Its proper functioning might be essential for the system, since bamA (yaeT) and bamD (yfiO), coding for the essential components of the assembly machinery of beta-barrel OMPs, as well as bamB (yfgL), the gene encoding an additional lipoprotein of the system, have been preserved . Additionally, it also retained the two chaperones Skp and SurA, which prevent folding and aggregation of OMPs in the periplasm during passage through the Sec translocon, and assist in their folding once they reach the assembly machinery in the outer membrane, respectively. Although DegP, the protease and chaperone identified to be involved in the degradation of misfolded OMPs, is not present, M. endobia encodes DegQ, another periplasmic protease which exhibits functional overlap with its homolog DegP [43, 44].
Only a limited set of active transporters are encoded in the M. endobia genome. Those include a phosphotransferase system for the transport of hexoses, ABC transporters for zinc, glutathione, lipopolysaccharides and lipidA, as well as a low-affinity inorganic phosphate transporter. Additionally, the M. endobia genome also codes for two channels associated with osmotic stress response, MscL and YbaL, which are absent in all Sternorrhyncha endosymbiont genomes sequenced so far. It is worth mentioning that, in addition to low molecular weight molecules, such as ions, metabolites and osmoprotectants, MscL is reported to be involved in the excretion of some small cytoplasmic proteins [45–47]. Therefore, it cannot be ruled out that the preservation of this mechanosensitive channel is an essential part of this peculiar endosymbiont nested system. MscL might be involved in the exchange of molecules between the two bacteria.
The detailed analysis of the functional capabilities of the two components of the nested endosymbiosis in P. citri suggests the existence of an intricate case of complementation, involving not only metabolic but also informational functions. Thus, despite the fact that M. endobia resembles B. aphidicola BCc , another endosymbiont with a highly reduced genome, in many functions such as transport, biosynthesis of cellular envelope and nucleotides, and its incapability to synthesize ATP coupled to the electron transport chain, it possesses particular characteristics that might be related to its coevolution with T. princeps. While complementation for amino acid biosynthesis has been described in other endosymbiotic systems, this is the first case in which all energy sources appear to be provided only by one of the partners, similarly to what happens in the eukaryotic cell, where the mitochondria is in charge of this function. Additionally, two genes encoding channels associated with osmotic stress response (mscL and ybaL) have been preserved in its genome. The fact that this kind of molecule has not been identified in other P-endosymbionts with reduced genomes might indicate their connection with special requirements of nested endosymbiosis, and might be involved in the exchange of molecules between both partners.
On the other hand, T. princeps does not resemble any known organelle, but it would not be reasonable to consider it, in a strict sense, as a living organism, since it has lost many essential genes involved in informational functions, as well as most metabolic pathways except for the ability to synthesize most essential amino acids, some of which require the cooperation of M. endobia and the host . T. princeps retains most, but not all, of the translation machinery, for which it also seems to depend on M. endobia, even though almost half of its coding capacity is devoted to this function [16, 19]. Additionally, it is unable to replicate on its own, although one can hypothesize that composite DNA and RNA polymerases (made of subunits encoded in both genomes) perform this function. T. princeps appears to be completely dependent on M. endobia for the synthesis of ATP, nucleotides or its cellular envelope, but still retains a complete set of molecular chaperones and proteins needed for the synthesis of [Fe-S] clusters.
Another intriguing fact revealed by our analysis is the overrepresentation of tRNAs genes in the M. endobia genome. This fact, together with the duplication in the rRNA operon in both genomes, appears to indicate an important translational activity in which both endosymbionts seem to be engaged. However, it lacks tRNA-Lys(AAG) which, surprisingly, has two functional copies in the small genome of T. princeps. This might be an indication that there is a mutual exchange of molecules between both compartments, although further studies are required to demonstrate this.
Nature is prolific in instances of symbiotic cooperation to give rise to new organisms, and new discoveries are always possible. Taking into consideration the deduced exceptional complementation inferred for this endosymbiotic system, we propose that T. princeps and M. endobia should be considered part of a new composite organism rather than a bacterial consortium.
Insect sample collection and DNA extraction
A population of P. citri from an initial sample obtained from a Cactaceae at the Botanical Garden of the Universitat de Valencia (Valencia, Spain) was reared in the laboratory at room temperature, fed on fresh pumpkins and used for genome sequencing. Two other populations of P. citri were used for additional experiments. One of them was obtained from a melon field in Murcia (Spain), the second one from a cultured population reared on germinated potatoes at the “Centro de Sanidad Vegetal” (Generalitat Valenciana) in Almassora (Castelló, Spain).
Total DNA enriched in bacterial endosymbionts was extracted from viscera of 20–30 adult female insects in sterile conditions and mechanically homogenized. In order to reduce insect DNA contamination, the samples were subjected to consecutive centrifugations at 1150 g and 1300 g for 10 minutes, and genomic DNA was obtained from the supernatant following a CTAB (Cetyltrimethylammonium bromide) extraction method .
Genome sequencing and assembly
The purified genomic DNA was shotgun sequenced using 454/Roche GS-FLX Titanium technology at the Genomics and Health area of the Public Health Research Center (CSISP, Generalitat Valenciana). One half-plate single-ends, and one-fourth plate paired-ends (3 kb of fragment size) sequencing experiments were performed, yielding a total of 1.3 million reads. Sequences of eukaryotic origin were eliminated after a taxonomic assignation process by Galaxy . Filtered reads were automatically assembled by MIRA  and the resulting contigs were manually edited with the Gap4 program from the Staden package software . The remaining gaps in the genome of M. endobia str. PCVAL were closed by ABI sequencing of PCR products obtained with designed primers, at the sequencing facility of the Universitat de València. Potential oriC on both genomes were sought with the OriginX program .
Total DNA samples obtained from the P. citri populations from Murcia and Almassora were used to further analyze the rplQ region from the T. princeps genome. The region comprised between genes rpoA and aroK was amplified and sequenced using the primers rpoA-F (5′-TGCCAGGCCTAGTGCTAAACATCA-3′) and aroK-R (5′-TGTCGCCAGGACTGCTATCAATGT-3′).
Gene annotation and functional analysis
ARAGORN , tRNAscan , and Rfam  sowftware packages were used for RNA genes prediction. Coding genes were annotated by BASys (Bacterial Annotation System, , RAST  and refined by BLAST searches . Finally, functional domain studies in Pfam database  were performed when coding-genes functionality assessment was required. Artemis  and MEGA5  programs were used for genome statistics calculation and codon usage analysis. Metabolic capabilities were analyzed with Blast2Go  and KAAS  programs. Functional information from the BioCyc , KEEG  and BRENDA  databases were also used in this context. Genome alignments were performed using MAFFT .
Annotated ORFs were considered as functional genes following two non-exclusionary criteria: the conservation of at least 80% of the sequence length of the closest orthologs found by BLAST in non-redundant databases, and/or the maintenance of the essential functional domains detected by Pfam .
The genome sequence of M. endobia strain PCVAL has been deposited at the GenBank (accession number CP003881). The GenBank accession numbers of the other three genome sequences used in this study are as follows: “Ca. Tremblaya princeps” str. PCVAL, CP002918; “Ca. Tremblaya princeps” str. PCIT, CP002244; M. endobia strain PCIT, CP002243.
We thank Dr. Ferran Garcia (Universitat Politècnica de Valencia, Spain) and Alberto García (Centro de Sanidad Vegetal, Generalitat Valenciana, Almassora, Spain) for providing mealybug samples. Financial support was provided by grants BFU2009-12895-C02-01/BMC (Ministerio de Ciencia e Innovación, Spain) and BFU2012-39816-C02-01 (Ministerio de Economía y Competitividad, Spain) to A. Latorre and by grant Prometeo/2009/092 (Conselleria d’Educació, Generalitat Valenciana, Spain) to A. Moya. S. López-Madrigal is a recipient of a fellowship from the Ministerio de Educación (Spain).
- Moya A, Pereto J, Gil R, Latorre A: Learning how to live together: genomic insights into prokaryote-animal symbioses. Nat Rev Genet. 2008, 9: 218-229. 10.1038/nrg2319.PubMedView Article
- McCutcheon JP, Moran NA: Extreme genome reduction in symbiotic bacteria. Nat Rev Microbiol. 2011, 10: 13-26.PubMed
- Watson RA: The impact of sex, symbiosis and modularity on the gradualist framework of evolution. 2006, Cambridge (Massachusetts): The MIT Press
- Gil R, Latorre A, Moya A: Evolution of prokaryote-animal symbiosis from a genomics perspective. (Endo)symbiotic Methanogenic Archaea. Edited by: Hackstein JHP. 2010, Berlin Heidelberg: Springer, 207-233. [Steinbüchel A (Series Editor): Microbiology Monographs, vol. 19]View Article
- Lamelas A, Gosalbes MJ, Manzano-Marin A, Pereto J, Moya A, Latorre A: Serratia symbiotica from the aphid Cinara cedri: a missing link from facultative to obligate insect endosymbiont. PLoS Genet. 2011, 7: e1002357-10.1371/journal.pgen.1002357.PubMedPubMed CentralView Article
- Wu D, Daugherty SC, Van Aken SE, Pai GH, Watkins KL, Khouri H, Tallon LJ, Zaborsky JM, Dunbar HE, Tran PL: Metabolic complementarity and genomics of the ual bacterial symbiosis of sharpshooters. PLoS Biol. 2006, 4: e188-10.1371/journal.pbio.0040188.PubMedPubMed CentralView Article
- McCutcheon JP, McDonald BR, Moran NA: Origin of an alternative genetic code in the extremely small and GC-rich genome of a bacterial symbiont. PLoS Genet. 2009, 5: e1000565-10.1371/journal.pgen.1000565.PubMedPubMed CentralView Article
- McCutcheon JP, Moran NA: Functional convergence in reduced genomes of bacterial symbionts spanning 200 MY of evolution. Genome Biol Evol. 2010, 2: 708-718.PubMedPubMed Central
- Lefevre C, Charles H, Vallier A, Delobel B, Farrell B, Heddi A: Endosymbiont phylogenesis in the Dryophthoridae weevils: evidence for bacterial replacement. Mol Biol Evol. 2004, 21: 965-973. 10.1093/molbev/msh063.PubMedView Article
- ScaleNet. http://www.sel.barc.usda.gov/scalenet/scalenet.htm,
- Hardy NB, Gullan PJ, Hodgson CJ: A subfamily-level classification of mealybugs (Hemiptera: Pseudococcidae) based on integrated molecular and morphological data. Syst Entomol. 2008, 33: 51-71. 10.1111/j.1365-3113.2007.00408.x.View Article
- Munson MA, Baumann P, Moran NA: Phylogenetic relationships of the endosymbionts of mealybugs (Homoptera: Pseudococcidae) based on 16S rDNA sequences. Mol Phylogenet Evol. 1992, 1: 26-30. 10.1016/1055-7903(92)90032-C.PubMedView Article
- Gruwell ME, Hardy NB, Gullan PJ, Dittmar K: Evolutionary relationships among primary endosymbionts of the mealybug subfamily Phenacoccinae (Hemiptera: Coccoidea: Pseudococcidae). Appl Environ Microbiol. 2010, 76: 7521-7525. 10.1128/AEM.01354-10.PubMedPubMed CentralView Article
- Thao ML, Gullan PJ, Baumann P: Secondary (gamma-Proteobacteria) endosymbionts infect the primary (beta-Proteobacteria) endosymbionts of mealybugs multiple times and coevolve with their hosts. Appl Environ Microbiol. 2002, 68: 3190-3197. 10.1128/AEM.68.7.3190-3197.2002.PubMedPubMed CentralView Article
- Von Dohlen CD, Kohler S, Alsop ST, McManus WR: Mealybug betaproteobacterial endosymbionts contain gamma-proteobacterial symbionts. Nature. 2001, 412: 433-436. 10.1038/35086563.PubMedView Article
- McCutcheon JP, Von Dohlen CD: An interdependent metabolic patchwork in the nested symbiosis of mealybugs. Curr Biol. 2011, 21: 1366-1372. 10.1016/j.cub.2011.06.051.PubMedPubMed CentralView Article
- Kono M, Koga R, Shimada M, Fukatsu T: Infection dynamics of coexisting beta and gammaproteobacteria in the nested endosymbiotic system of mealybugs. Appl Environ Microbiol. 2008, 74: 4175-4184. 10.1128/AEM.00250-08.PubMedPubMed CentralView Article
- Baumann L, Thao ML, Hess JM, Johnson MW, Baumann P: The genetic properties of the primary endosymbionts of mealybugs differ from those of other endosymbionts of plant sap-sucking insects. Appl Environ Microbiol. 2002, 68: 3198-3205. 10.1128/AEM.68.7.3198-3205.2002.PubMedPubMed CentralView Article
- Lopez-Madrigal S, Latorre A, Porcar M, Moya A, Gil R: Complete genome sequence of “Candidatus Tremblaya princeps” strain PCVAL, an intriguing translational machine below the living-cell status. J Bacteriol. 2011, 193: 5587-5588. 10.1128/JB.05749-11.PubMedPubMed CentralView Article
- Gil R, Latorre A, Moya A: Bacterial endosymbionts of insects: insights from comparative genomics. Environ Microbiol. 2004, 6: 1109-1122. 10.1111/j.1462-2920.2004.00691.x.PubMedView Article
- Gil R, Silva FJ, Zientz E, Delmotte F, Gonzalez-Candelas F, Latorre A, Rausell C, Kamerbeek J, Gadau J, Holldobler B, Van Ham RCHJ, Gross R, Moya A: The genome sequence of Blochmannia floridanus: Comparative analysis of reduced genomes. Proc Nat Acad Sci USA. 2003, 100: 9388-9393. 10.1073/pnas.1533499100.PubMedPubMed CentralView Article
- Akman L, Yamashita A, Watanabe H, Oshima K, Shiba T, Hattori M, Aksoy S: Genome sequence of the endocellular obligate symbiont of tsetse flies, Wigglesworthia glossinidia. Nat Genet. 2002, 32: 402-407. 10.1038/ng986.PubMedView Article
- Nakabachi A, Yamashita A, Toh H, Ishikawa H, Dunbar HE, Moran NA, Hattori M: The 160-kilobase genome of the bacterial endosymbiont Carsonella. Science. 2006, 314: 267-10.1126/science.1134196.PubMedView Article
- McCutcheon JP, McDonald BR, Moran NA: Convergent evolution of metabolic roles in bacterial co-symbionts of insects. Proc Nat Acad Sci USA. 2009, 106: 15394-15399. 10.1073/pnas.0906424106.PubMedPubMed CentralView Article
- McCutcheon JP, Moran NA: Parallel genomic evolution and metabolic interdependence in an ancient symbiosis. Proc Nat Acad Sci USA. 2007, 104: 19392-19397. 10.1073/pnas.0708855104.PubMedPubMed CentralView Article
- Gil R, Latorre A: Factors behind junk DNA in bacteria. Genes. 2012, 3: 634-650. 10.3390/genes3040634.PubMedPubMed CentralView Article
- Hershberg R, Petrov DA: Evidence that mutation is universally biased towards AT in bacteria. PLoS Genet. 2010, 6: e1001115-10.1371/journal.pgen.1001115.PubMedPubMed CentralView Article
- Hildebrand F, Meyer A, Eyre-Walker A: Evidence of selection upon genomic GC-content in bacteria. PLoS Genet. 2010, 6: e1001107-10.1371/journal.pgen.1001107.PubMedPubMed CentralView Article
- Behrens S, Maier R, De Cock H, Schmid FX, Gross CA: The SurA periplasmic PPIase lacking its parvulin domains functions in vivo and has chaperone activity. EMBO J. 2001, 20: 285-294. 10.1093/emboj/20.1.285.PubMedPubMed CentralView Article
- Dermic D: Functions of multiple exonucleases are essential for cell viability, DNA repair and homologous recombination in recD mutants of Escherichia coli. Genetics. 2006, 172: 2057-2069.PubMed
- Heller RC, Marians KJ: The disposition of nascent strands at stalled replication forks dictates the pathway of replisome loading during restart. Mol Cell. 2005, 17: 733-743. 10.1016/j.molcel.2005.01.019.PubMedView Article
- Xu L, Marians KJ: Purification and characterization of DnaC810, a primosomal protein capable of bypassing PriA function. J Biol Chem. 2000, 275: 8196-8205. 10.1074/jbc.275.11.8196.PubMedView Article
- Fraser CM, Gocayne JD, White O, Adams MD, Clayton RA, Fleischmann RD, Bult CJ, Kerlavage AR, Sutton G, Kelley JM, Fritchman JL, Weidman JF, Small KV, Sandusky M, Fuhrmann J, Nguyen D, Utterback TR, Saudek DM, Phillips CA, Merrick JM, Tomb JF, Dougherty BA, Bott KF, Hu PC, Lucier TS, Peterson SN, Smith HO, Hutchison CA, Venter JC: The minimal gene complement of Mycoplasma genitalium. Science. 1995, 270: 397-403. 10.1126/science.270.5235.397.PubMedView Article
- Lopper M, Boonsombat R, Sandler SJ, Keck JL: A hand-off mechanism for primosome assembly in replication restart. Mol Cell. 2007, 26: 781-793. 10.1016/j.molcel.2007.05.012.PubMedPubMed CentralView Article
- Gil R, Silva FJ, Pereto J, Moya A: Determination of the core of a minimal bacterial gene set. Microbiol Mol Biol Rev. 2004, 68: 518-537. 10.1128/MMBR.68.3.518-537.2004.PubMedPubMed CentralView Article
- Quan S, Zhang N, French S, Squires CL: Transcriptional polarity in rRNA operons of Escherichia coli nusA and nusB mutant strains. J Bacteriol. 2005, 187: 1632-1638. 10.1128/JB.187.5.1632-1638.2005.PubMedPubMed CentralView Article
- Price NL, Raivio TL: Characterization of the Cpx regulon in Escherichia coli strain MC4100. J Bacteriol. 2009, 191: 1798-1815. 10.1128/JB.00798-08.PubMedPubMed CentralView Article
- Tseng TT, Tyler BM, Setubal JC: Protein secretion systems in bacterial-host associations, and their description in the Gene Ontology. BMC Microbiol. 2009, 9 (Suppl 1): S2-10.1186/1471-2180-9-S1-S2.PubMedPubMed CentralView Article
- Perez-Brocal V, Gil R, Ramos S, Lamelas A, Postigo M, Michelena JM, Silva FJ, Moya A, Latorre A: A small microbial genome: the end of a long symbiotic relationship?. Science. 2006, 314: 312-313. 10.1126/science.1130441.PubMedView Article
- Johnson DC, Dean DR, Smith AD, Johnson MK: Structure, function, and formation of biological iron-sulfur clusters. Annu Rev Biochem. 2005, 74: 247-281. 10.1146/annurev.biochem.74.082803.133518.PubMedView Article
- Beinert H: Iron-sulfur proteins: ancient structures, still full of surprises. J Biol Inorg Chem. 2000, 5: 2-15. 10.1007/s007750050002.PubMedView Article
- Malinverni JC, Werner J, Kim S, Sklar JG, Kahne D, Misra R, Silhavy TJ: YfiO stabilizes the YaeT complex and is essential for outer membrane protein assembly in Escherichia coli. Mol Microbiol. 2006, 61: 151-164. 10.1111/j.1365-2958.2006.05211.x.PubMedView Article
- Singh N, Kuppili RR, Bose K: The structural basis of mode of activation and functional diversity: a case study with HtrA family of serine proteases. Arch Biochem Biophys. 2011, 516: 85-96. 10.1016/j.abb.2011.10.007.PubMedView Article
- Sawa J, Malet H, Krojer T, Canellas F, Ehrmann M, Clausen T: Molecular adaptation of the DegQ protease to exert protein quality control in the bacterial cell envelope. J Biol Chem. 2011, 286: 30680-30690. 10.1074/jbc.M111.243832.PubMedPubMed CentralView Article
- Ajouz B, Berrier C, Garrigues A, Besnard M, Ghazi A: Release of thioredoxin via the mechanosensitive channel MscL during osmotic downshock of Escherichia coli cells. J Biol Chem. 1998, 273: 26670-26674. 10.1074/jbc.273.41.26670.PubMedView Article
- Berrier C, Garrigues A, Richarme G, Ghazi A: Elongation factor Tu and DnaK are transferred from the cytoplasm to the periplasm of Escherichia coli during osmotic downshock presumably via the mechanosensitive channel mscL. J Bacteriol. 2000, 182: 248-251. 10.1128/JB.182.1.248-251.2000.PubMedPubMed CentralView Article
- van den Bogaart G, Krasnikov V, Poolman B: Dual-color fluorescence-burst analysis to probe protein efflux through the mechanosensitive channel MscL. Biophys J. 2007, 92: 1233-1240. 10.1529/biophysj.106.088708.PubMedPubMed CentralView Article
- Ausubel F: Short Protocols in Molecular Biology: A Compendium of Methods from Current Protocols in Molecular Biology. 1999, New York: Wiley, 44
- Galaxy. http://galaxyproject.org,
- MIRA. http://sourceforge.net/projects/mira-assembler,
- Staden R, Beal KF, Bonfield JK: The Staden package, 1998. Methods Mol Biol. 2000, 132: 115-130.PubMed
- Worning P, Jensen LJ, Hallin PF, Staerfeldt HH, Ussery DW: Origin of replication in circular prokaryotic chromosomes. Environ Microbiol. 2006, 8: 353-361. 10.1111/j.1462-2920.2005.00917.x.PubMedView Article
- Laslett D, Canback B: ARAGORN, a program to detect tRNA genes and tmRNA genes in nucleotide sequences. Nucleic Acids Res. 2004, 32: 11-16. 10.1093/nar/gkh152.PubMedPubMed CentralView Article
- Lowe TM, Eddy SR: tRNAscan-SE: a program for improved detection of transfer RNA genes in genomic sequence. Nucleic Acids Res. 1997, 25: 955-964.PubMedPubMed CentralView Article
- Gardner PP, Daub J, Tate JG, Nawrocki EP, Kolbe DL, Lindgreen S, Wilkinson AC, Finn RD, Griffiths-Jones S, Eddy SR, Bateman A: Rfam: updates to the RNA families database. Nucleic Acids Res. 2009, 37 (Database issue): D136-D140.PubMedPubMed CentralView Article
- Van Domselaar GH, Stothard P, Shrivastava S, Cruz JA, Guo A, Dong X, Lu P, Szafron D, Greiner R, Wishart DS: BASys: a web server for automated bacterial genome annotation. Nucleic Acids Res. 2005, 33 (Web Server issue): W455-W459.PubMedPubMed CentralView Article
- Aziz RK, Bartels D, Best AA, DeJongh M, Disz T, Edwards RA, Formsma K, Gerdes S, Glass EM, Kubal M, Meyer F, Olsen GJ, Olson R, Osterman AL, Overbeek RA, McNeil LK, Paarmann D, Paczian T, Parrello B, Pusch GD, Reich C, Stevens R, Vassieva O, Vonstein V, Wilke A, Zagnitko O: The RAST Server: rapid annotations using subsystems technology. BMC Genomics. 2008, 9: 75-10.1186/1471-2164-9-75.PubMedPubMed CentralView Article
- Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997, 25: 3389-3402. 10.1093/nar/25.17.3389.PubMedPubMed CentralView Article
- Pfam. http://pfam.sanger.ac.uk,
- Carver T, Berriman M, Tivey A, Patel C, Bohme U, Barrell BG, Parkhill J, Rajandream MA: Artemis and ACT: viewing, annotating and comparing sequences stored in a relational database. Bioinformatics. 2008, 24: 2672-2676. 10.1093/bioinformatics/btn529.PubMedPubMed CentralView Article
- Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S: MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol. 2011, 28: 2731-2739. 10.1093/molbev/msr121.PubMedPubMed CentralView Article
- Blast2Go. http://www.blast2go.com,
- Moriya Y, Itoh M, Okuda S, Yoshizawa AC, Kanehisa M: KAAS: an automatic genome annotation and pathway reconstruction server. Web Server issue. 2007, 35: W182-W185.
- BioCyc. http://biocyc.org,
- KEEG. http://www.genome.jp/kegg,
- BRENDA. http://www.brenda-enzymes.info,
- Katoh K, Misawa K, Kuma K, Miyata T: MAFFT: a novel method for rapid multiple sequence alignment based on fast Fourier transform. Nucleic Acids Res. 2002, 30: 3059-3066. 10.1093/nar/gkf436.PubMedPubMed CentralView Article
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.