Molecular systematics of aphids and their primary endosymbionts

Comparative genomics of the primary endosymbiont Buchnera aphidicola in aphid hosts and their coevolutionary relationships

BACKGROUND

Coevolution between modern aphids and their primary obligate, bacterial endosymbiont, Buchnera aphidicola, has been previously reported at different classification levels based on molecular phylogenetic analyses. However, the Buchnera genome remains poorly understood within the Rhus gall aphids.

RESULTS

We assembled the complete genome of the endosymbiont Buchnera in 16 aphid samples, representing 13 species in all six genera of Rhus gall aphids by shotgun genome skimming method. We compared the newly assembled genomes with those from GenBank to comprehensively investigate patterns of coevolution between the bacteria Buchnera and their aphid hosts. Buchnera genomes were mostly collinear, and the pan-genome contained 684 genes, in which the core genome contained 256 genes with some lineages having large numbers of tandem gene duplications. There has been substantial gene-loss in each Buchnera lineage. We also reconstructed the phylogeny for Buchnera and their host aphids, respectively, using 72 complete genomes of Buchnera, along with the complete mitochondrial genomes and three nuclear genes of 31 corresponding host aphid accessions. The cophylogenetic test demonstrated significant coevolution between these two partner groups at individual, species, generic, and tribal levels.

CONCLUSIONS

Buchnera exhibits very high levels of genomic sequence divergence but relative stability in gene order. The relationship between the symbionts Buchnera and its aphid hosts shows a significant coevolutionary pattern and supports complexity of the obligate symbiotic relationship.

Ultraconserved elements from transcriptome and genome data provide insight into the phylogenomics of Sternorrhyncha (Insecta: Hemiptera)

Sternorrhyncha, one of the four major suborders of Hemiptera, is a phytophagous taxon inclusive of nearly 18 000 described species. The phylogenetic relationships within the taxon and the earliest-branching lineage of its infraorders remain incompletely understood. This study attempted to illuminate the phylogenetic relationships within Sternorrhyncha through the use of maximum likelihood, Bayesian inference and maximum parsimony analyses, employing ultraconserved element (UCE) data from 39 genomic and 62 transcriptomic datasets and thereby representing most families within the taxon. The probe set Hemiptera 2.7Kv1 was used to recover a total of 2731 UCE loci: from 547 to 1699 (with an average of 1084) across all genomic datasets and from 108 to 849 (with an average of 329) across all transcriptomic datasets. All three types of phylogenetic analyses employed in this study produced robust statistical support for Sternorrhyncha being a monophyletic group. The different methods of phylogenetic analysis produced inconsistent descriptions of topological structure at the infraorder level: while maximum likelihood and Bayesian inference analyses produced strong statistical evidence (100%) indicating the clade Psylloidea + Aleyrodoidea to be a sister of the clade Aphidoidea (Aphidomorpha) + Coccoidea (Coccomorpha), the maximum parsimony analysis failed to recover a similar result. Our results also provide detail on the phylogenetic relationships within each infraorder. This study presents the first use of UCE data to investigate the phylogeny of Sternorrhyncha. It also shows the viability of amalgamating genomic and transcriptomic data in studies of phylogenetic relationships, potentially highlighting a resource-efficient approach for future inquiries into diverse taxa through the integration of varied data sources.

Genome evolution of Buchnera aphidicola (Gammaproteobacteria): Insights into strand compositional asymmetry, codon usage bias, and phylogenetic implications

Taxa of Buchnera aphidicola (hereafter "Buchnera") are mutualistic intracellular symbionts of aphids, known for their remarkable biological traits such as genome reduction, strand compositional asymmetry, and symbiont-host coevolution. With the growing availability of genomic data, we performed a comprehensive analysis of 103 genomes of Buchnera strains from 12 host subfamilies, focusing on the genomic characterizations, codon usage patterns, and phylogenetic implications. Our findings revealed consistent features among all genomes, including small genome sizes, low GC contents, and gene losses. We also identified strong strand compositional asymmetries in all strains at the genome level. Further investigation suggested that mutation pressure may have played a crucial role in shaping codon usage of Buchnera. Moreover, the genomic asymmetries were reflected in asymmetric codon usage preferences within chromosomal genes. Notably, the levels of these asymmetries were varied among strains and were significantly influenced by the degrees of genome shrinkages. Lastly, our phylogenetic analyses presented an alternative topology of Aphididae, based on the Buchnera symbionts, providing robust confirmation of the paraphylies of Eriosomatinae, and Macrosiphini. Our objectives are to further understand the strand compositional asymmetry and codon usage bias of Buchnera taxa, and provide new perspectives for phylogenetic studies of Aphididae.

Portiera Gets Wild: Genome Instability Provides Insights into the Evolution of Both Whiteflies and Their Endosymbionts

Whiteflies (Hemiptera: Sternorrhyncha: Aleyrodidae) are a superfamily of small phloem-feeding insects. They rely on their primary endosymbionts "Candidatus Portiera aleyrodidarum" to produce essential amino acids not present in their diet. Portiera has been codiverging with whiteflies since their origin and therefore reflects its host's evolutionary history. Like in most primary endosymbionts, the genome of Portiera stays stable across the Aleyrodidae superfamily after millions of years of codivergence. However, Portiera of the whitefly Bemisia tabaci has lost the ancestral genome order, reflecting a rare event in the endosymbiont evolution: the appearance of genome instability. To gain a better understanding of Portiera genome evolution, identify the time point in which genome instability appeared and contribute to the reconstruction of whitefly phylogeny, we developed a new phylogenetic framework. It targeted five Portiera genes and determined the presence of the DNA polymerase proofreading subunit (dnaQ) gene, previously associated with genome instability, and two alternative gene rearrangements. Our results indicated that Portiera gene sequences provide a robust tool for studying intergenera phylogenetic relationships in whiteflies. Using these new framework, we found that whitefly species from the Singhiella, Aleurolobus, and Bemisia genera form a monophyletic tribe, the Aleurolobini, and that their Portiera exhibit genome instability. This instability likely arose once in the common ancestor of the Aleurolobini tribe (at least 70 Ma), drawing a link between the appearance of genome instability in Portiera and the switch from multibacteriocyte to a single-bacteriocyte mode of inheritance in this tribe.

Cophylogenetic analysis suggests cospeciation between the Scorpion Mycoplasma Clade symbionts and their hosts

Scorpions are predator arachnids of ancient origin and worldwide distribution. Two scorpion species, Vaejovis smithi and Centruroides limpidus, were found to harbor two different Mollicutes phylotypes: a Scorpion Mycoplasma Clade (SMC) and Scorpion Group 1 (SG1). Here we investigated, using a targeted gene sequencing strategy, whether these Mollicutes were present in 23 scorpion morphospecies belonging to the Vaejovidae, Carboctonidae, Euscorpiidae, Diplocentridae, and Buthidae families. Our results revealed that SMC is found in a species-specific association with Vaejovidae and Buthidae, whereas SG1 is uniquely found in Vaejovidae. SMC and SG1 co-occur only in Vaejovis smithi where 43% of the individuals host both phylotypes. A phylogenetic analysis of Mollicutes 16S rRNA showed that SMC and SG1 constitute well-delineated phylotypes. Additionally, we found that SMC and scorpion phylogenies are significantly congruent, supporting the observation that a cospeciation process may have occurred. This study highlights the phylogenetic diversity of the scorpion associated Mollicutes through different species revealing a possible cospeciation pattern.

Insect-bacteria parallel evolution in multiple-co-obligate-aphid association: a case in Lachninae (Hemiptera: Aphididae)

Parallel phylogenies between aphid and its obligate symbiont Buchnera are hot topics which always focused on aphid lower taxonomic levels. Symbionts in the subfamily Lachninae are special. Buchnera in many lachnine species has undergone functional and genome size reduction that was replaced by other co-obligate symbionts. In this study, we constructed the phylogenetic relationships of Lachninae with a combined dataset of five genes sequenced from Buchnera to estimate the effects of a dual symbiotic system in the aphid-Buchnera cospeciation association. The phylogeny of Buchnera in Lachninae was well-resolved in the combined dataset. Each of the genera formed strongly supported monophyletic groups, with the exception of the genus Cinara. The phylogeny based on sequences from Buchnera was divided into five tribes according to the clades of the Lachninae hosts tree, with the phylogenies of Buchnera and Lachninae being generally congruent. These results first provided evidence of parallel evolution at the aphid subfamily level comprehensively and supported the view that topological congruence between the phylogenies of Buchnera and Lachninae would not be interfered with the other co-obligate symbionts, such as Sarretia, in aphid-entosymbiont association. These results also provided new insight in understanding host-plant coevolution in lachnine lineages.

Polyphyly of the extinct family Oviparosiphidae and its implications for inferring aphid evolution (Hemiptera, Sternorrhyncha)

Aphidoidea, the so-called "true aphids" are one of the most challenging groups in terms of solving the phylogenetic relationships. Morphology-based analyses were strongly affected by widespread homoplasy, while the molecular-based attempts struggled with the lack of sufficient phylogenetic signal. Despite significant improvements, the higher classification still remains unresolved and rather controversial. However, the use of the fossil record, one of the most valuable sources of information, was mainly limited to calibration of a phylogenetic tree, without a direct inclusion into the analysis. The extinct family Oviparosiphidae has long been considered as the common ancestor of all recent Aphidoidea and it was used as a calibration point in several analyses, but it has been never analyzed in a phylogenetic context. The family has been treated as a monophyletic group purely based on the simultaneous presence of two abdominal structures, ovipositor and siphunculi. However, it has been shown recently that at least one more extinct lineage, present at the same time, was characterized by the same features. For these reasons, we performed a maximum parsimony analysis using morphological data for extinct aphid taxa to prove the monophyly of Oviparosiphidae. Our analysis shows that the presumed ancestor lineage of recent aphids is a polyphyletic group. Our results support the hypothesis of an early Mesozoic rapid radiation of aphids, which led to several different lineages characterized by both ovipositor and siphunculi. The results indicate the necessity of examining the other extinct families, and shows that the diversity of aphids before the Cretaceous Terrestrial Revolution (KTR) was higher than expected. Even though there is not enough data to perform a formal analysis, fossils seem to suggest a significant impact of the KTR on aphid diversification. Additionally, we have made a redescription of two genera and description of a new species, Vitimaphis subridens sp. nov.

The relationships within the Chaitophorinae and Drepanosiphinae (Hemiptera, Aphididae) inferred from molecular-based phylogeny and comprehensive morphological data

The Chaitophorinae is a bionomically diverse Holarctic subfamily of Aphididae. The current classification includes two tribes: the Chaitophorini associated with deciduous trees and shrubs, and Siphini that feed on monocotyledonous plants. We present the first phylogenetic hypothesis for the subfamily, based on molecular and morphological datasets. Molecular analyses were based on the mitochondrial gene cytochrome oxidase subunit I (COI) and the nuclear gene elongation factor-1α (EF-1α). Phylogenetic inferences were obtained individually on each of genes and joined alignments using Bayesian inference (BI) and Maximum likelihood (ML). In phylogenetic trees reconstructed on the basis of nuclear and mitochondrial genes as well as a morphological dataset, the monophyly of Siphini and the genus Chaitophorus was supported. Periphyllus forms independent lineages from Chaitophorus and Siphini. Within this genus two clades comprising European and Asiatic species, respectively, were indicated. Concerning relationships within the subfamily, EF-1α and joined COI and EF-1α genes analysis strongly supports the hypothesis that Chaitophorini do not form a monophyletic clade. Periphyllus is a sister group to a clade containing Chaitophorus and Siphini. The Asiatic unit of Periphyllus also includes Trichaitophorus koyaensis. The analysis of morphological dataset under equally weighted parsimony also supports the view that Chaitophorini is an artificial taxon, as Lambersaphis pruinosae and Pseudopterocomma hughi, both traditionally included in the Chaitophorini, formed independent lineages. COI analyses support consistent groups within the subfamily, but relationships between groups are poorly resolved. These analyses were extended to include the species of closely related and phylogenetically unstudied subfamily Drepanosiphinae, which produced congruent results. Genera Drepanosiphum and Depanaphis are monophyletic and sister. The position of Yamatocallis tokyoensis differs in the molecular and morphological analyses, i.e. it is either an independent lineage (EF-1α, COI, joined COI and EF-1α genes) or is nested inside this unit (morphology). Our data also support separation of Chaitophorinae from Drepanosiphinae.

Phylogenetic Relationships of the Symbiotic Bacteria in the Aphid Sitobion avenae (Hemiptera: Aphididae)

Aphids have developed symbiotic associations with different bacterial species, and some morphological and molecular analyses have provided evidence of the host relationship between the primary symbiotic bacteria (Buchnera aphidicola) and the aphid while the contrary with the secondary symbiotic bacteria. In this study, we investigated the phylogenetic relationships of the bacterial endosymbionts in the aphid Sitobion avenae (F.). We characterized all bacterial endosymbionts in 10 genetically defined S. avenae clones by denaturing gradient gel electrophoresis and, from these clones, sequenced the 16S rRNA genes of both the primary endosymbiont, B. aphidicola (for the first time), and the secondary endosymbionts, Regiella insecticola and Hamiltonella defensa (for the first time). The phylogenetic analysis indicated that Buchnera from Sitobion related to those in Macrosiphoni. The analysis of the secondary endosymbionts indicated that there is no host relationship between H. defensa and R. insecticola from Sitobion and those from other aphid species. In this study, therefore, we identified further evidence for the relationship between Buchnera and its host and reported a relationship within the secondary endosymbionts of S. avenae from the same country, even though there were no relationships between the secondary bacteria and their host. We also discussed the diversity within the symbiotic bacteria in S. avenae clones.

Avoidance and Potential Remedy Solutions of Chimeras in Reconstructing the Phylogeny of Aphids Using the 16S rRNA Gene of Buchnera: A Case in Lachninae (Hemiptera)

It is known that PCR amplification of highly homologous genes from complex DNA mixtures can generate a significant proportion of chimeric sequences. The 16S rRNA gene is not only widely used in estimating the species diversity of endosymbionts in aphids but also used to explore the co-diversification of aphids and their endosymbionts. Thus, chimeric sequences may lead to the discovery of non-existent endosymbiont species and mislead Buchnera-based phylogenetic analysis that lead to false conclusions. In this study, a high probability (6.49%) of chimeric sequence occurrence was found in the amplified 16S rRNA gene sequences of endosymbionts from aphid species in the subfamily Lachninae. These chimeras are hybrid products of multiple parent sequences from the dominant species of endosymbionts in each corresponding host. It is difficult to identify the chimeric sequences of a new or unidentified species due to the high variability of their main parent, Buchnera aphidicola, and because the chimeric sequences can confuse the phylogenetic analysis of 16S rRNA gene sequences. These chimeras present a challenge to Buchnera-based phylogenetic research in aphids. Thus, our study strongly suggests that using appropriate methods to detect chimeric 16S rRNA sequences may avoid some false conclusions in endosymbiont-based aphid research.

Macroevolutionary patterns in the Aphidini aphids (Hemiptera: Aphididae): diversification, host association, and biogeographic origins

Background

Due to its biogeographic origins and rapid diversification, understanding the tribe Aphidini is key to understanding aphid evolution. Major questions about aphid evolution include origins of host alternation as well as age and patterns of diversification in relation to host plants. To address these questions, we reconstructed the phylogeny of the Aphidini which contains Aphis, the most diverse genus in the family. We used a combined dataset of one nuclear and four mitochondrial DNA regions. A molecular dating approach, calibrated with fossil records, was used to estimate divergence times of these taxa.

Principal Findings

Most generic divergences in Aphidini occurred in the Middle Tertiary, and species-level divergences occurred between the Middle and Late Tertiary. The ancestral state of host use for Aphidini was equivocal with respect to three states: monoecy on trees, heteroecy, and monoecy on grasses. The ancestral state of Rhopalosiphina likely included both heteroecy and monoecy, whereas that of Aphidina was most likely monoecy. The divergence times of aphid lineages at the generic or subgeneric levels are close to those of their primary hosts. The species-level divergences in aphids are consistent with the diversification of the secondary hosts, as a few examples suggest. The biogeographic origin of Aphidini as a whole was equivocal, but the major lineages within Aphidina likely separated into Nearctic, Western Palearctic, and Eastern Palearctic regions.

Conclusions

Most generic divergences in Aphidini occurred in the Middle Tertiary when primary hosts, mainly in the Rosaceae, were diverging, whereas species-level divergences were contemporaneous with diversification of the secondary hosts such as Poaceae in the Middle to Late Tertiary. Our results suggest that evolution of host alternation within Aphidini may have occurred during the Middle Tertiary (Oligocene) when the secondary hosts emerged.

Sequence conservation and functional constraint on intergenic spacers in reduced genomes of the obligate symbiont Buchnera

Analyses of genome reduction in obligate bacterial symbionts typically focus on the removal and retention of protein-coding regions, which are subject to ongoing inactivation and deletion. However, these same forces operate on intergenic spacers (IGSs) and affect their contents, maintenance, and rates of evolution. IGSs comprise both non-coding, non-functional regions, including decaying pseudogenes at varying stages of recognizability, as well as functional elements, such as genes for sRNAs and regulatory control elements. The genomes of Buchnera and other small genome symbionts display biased nucleotide compositions and high rates of sequence evolution and contain few recognizable regulatory elements. However, IGS lengths are highly correlated across divergent Buchnera genomes, suggesting the presence of functional elements. To identify functional regions within the IGSs, we sequenced two Buchnera genomes (from aphid species Uroleucon ambrosiae and Acyrthosiphon kondoi) and applied a phylogenetic footprinting approach to alignments of orthologous IGSs from a total of eight Buchnera genomes corresponding to six aphid species. Inclusion of these new genomes allowed comparative analyses at intermediate levels of divergence, enabling the detection of both conserved elements and previously unrecognized pseudogenes. Analyses of these genomes revealed that 232 of 336 IGS alignments over 50 nucleotides in length displayed substantial sequence conservation. Conserved alignment blocks within these IGSs encompassed 88 Shine-Dalgarno sequences, 55 transcriptional terminators, 5 Sigma-32 binding sites, and 12 novel small RNAs. Although pseudogene formation, and thus IGS formation, are ongoing processes in these genomes, a large proportion of intergenic spacers contain functional sequences.

Diversification of genes for carotenoid biosynthesis in aphids following an ancient transfer from a fungus

The pea aphid genome was recently found to harbor genes for carotenoid biosynthesis, reflecting an ancestral transfer from a fungus. To explore the evolution of the carotene desaturase gene family within aphids, sequences were retrieved from a set of 34 aphid species representing numerous deeply diverging lineages of aphids and analyzed together with fungal sequences retrieved from databases. All aphids have at least one copy of this gene and some aphid species have up to seven, whereas fungal genomes consistently have a single copy. The closest relatives of aphids, adelgids, also have carotene desaturase; these sequences are most closely related to those from aphids, supporting a shared origin from a fungal to insect transfer predating the divergence of adelgids and aphids. Likewise, all aphids, and adelgids, have carotenoid profiles that are consistent with their biosynthesis using the acquired genes of fungal origin rather than derivation from food plants. The carotene desaturase was acquired from a fungal species outside of Ascomycota or Basidiomycota and closest to Mucoromycotina among sequences available in databases. In aphids, an ongoing pattern of gene duplication is indicated by the presence of both anciently and recently diverged paralogs within genomes and by the presence of a high frequency of pseudogenes that appear to be recently inactivated. Recombination among paralogs is evident, making analyses of patterns of selection difficult, but tests of selection for a nonrecombining region indicates that duplications tend to be followed by bouts of positive selection. Species of Macrosiphini, which often show color polymorphisms, typically have a larger number of desaturase copies relative to other species sampled in the study. These results indicate that aphid evolution has been accompanied by ongoing evolution of carotenogenic genes, which have undergone duplication, recombination, and occasional positive selection to yield a wide variety of carotenoid profiles in different aphid species.

New Insights on the Evolutionary History of Aphids and Their Primary Endosymbiont Buchnera aphidicola

Since the establishment of the symbiosis between the ancestor of modern aphids and their primary endosymbiont, Buchnera aphidicola, insects and bacteria have coevolved. Due to this parallel evolution, the analysis of bacterial genomic features constitutes a useful tool to understand their evolutionary history. Here we report, based on data from B. aphidicola, the molecular evolutionary analysis, the phylogenetic relationships among lineages and a comparison of sequence evolutionary rates of symbionts of four aphid species from three subfamilies. Our results support previous hypotheses of divergence of B. aphidicola and their host lineages during the early Cretaceous and indicate a closer relationship between subfamilies Eriosomatinae and Lachninae than with the Aphidinae. They also reveal a general evolutionary pattern among strains at the functional level. We also point out the effect of lifecycle and generation time as a possible explanation for the accelerated rate in B. aphidicola from the Lachninae.

Evolutionary relationships among primary endosymbionts of the mealybug subfamily phenacoccinae (hemiptera: Coccoidea: Pseudococcidae)

Mealybugs (Coccoidea: Pseudococcidae) are sap-sucking plant parasites that harbor bacterial endosymbionts within specialized organs. Previous studies have identified two subfamilies, Pseudococcinae and Phenacoccinae, within mealybugs and determined the primary endosymbionts (P-endosymbionts) of the Pseudococcinae to be Betaproteobacteria ("Candidatus Tremblaya princeps") containing Gammaproteobacteria secondary symbionts. Here, the P-endosymbionts of phenacoccine mealybugs are characterized based on 16S rRNA from the bacteria of 20 species of phenacoccine mealybugs and four outgroup Puto species (Coccoidea: Putoidae) and aligned to more than 100 published 16S rRNA sequences from symbiotic and free-living bacteria. Phylogenetic analyses recovered three separate lineages of bacteria from the Phenacoccinae, and these are considered to be the P-endosymbionts of their respective mealybug hosts, with those from (i) the mealybug genus Rastrococcus belonging to the Bacteroidetes, (ii) the subterranean mealybugs, tribe Rhizoecini, also within Bacteroidetes, in a clade sister to cockroach endosymbionts (Blattabacterium), and (iii) the remaining Phenacoccinae within the Betaproteobacteria, forming a well-supported sister group to "Candidatus Tremblaya princeps." Names are proposed for two strongly supported lineages: "Candidatus Brownia rhizoecola" for P-endosymbionts of Rhizoecini and "Candidatus Tremblaya phenacola" for P-endosymbionts of Phenacoccinae excluding Rastrococcus and Rhizoecini. Rates of nucleotide substitution among lineages of Tremblaya were inferred to be significantly faster than those of free-living Betaproteobacteria. Analyses also recovered a clade of Gammaproteobacteria, sister to the P-endosymbiont lineage of aphids ("Candidatus Buchnera aphidicola"), containing the endosymbionts of Putoidae, the secondary endosymbionts of pseudococcine mealybugs, and the endosymbionts of several other insect groups.

A phylogenetic analysis of armored scale insects (Hemiptera: Diaspididae), based upon nuclear, mitochondrial, and endosymbiont gene sequences

Armored scale insects (Hemiptera: Diaspididae) are among the most invasive insects in the world. They have unusual genetic systems, including diverse types of paternal genome elimination (PGE) and parthenogenesis. Intimate relationships with their host plants and bacterial endosymbionts make them potentially important subjects for the study of co-evolution. Here, we expand upon recent phylogenetic work (Morse and Normark, 2006) by analyzing armored scale and endosymbiont DNA sequences from 125 species of armored scale insect, represented by 253 samples and eight outgroup species. We used fragments of four different gene regions: the nuclear protein-coding gene Elongation Factor 1α (EF1α), the large ribosomal subunit (28S) rDNA, a mitochondrial region spanning parts of cytochrome oxidase I (COI) and cytochrome oxidase II (COII), and the small ribosomal subunit (16S) rDNA from the primary bacterial endosymbiont Uzinura diaspidicola. Maximum likelihood, and Bayesian analyses were performed producing highly congruent topological results. A comparison of two datasets, one with and one without missing data, found that missing data had little effect on topology. Our results broadly corroborate several major features of the existing classification, although we do not find any of the subfamilies, tribes or subtribes to be monophyletic as currently constituted. Using ancestral state reconstruction we estimate that the ancestral armored scale had the late PGE sex system, and it may as well have been pupillarial, though results differed between reconstruction methods. These results highlight the need for a complete revision of this family, and provide the groundwork for future taxonomic work in armored scale insects.

Aphid wing dimorphisms: linking environmental and genetic control of trait variation

Both genetic and environmental factors underlie phenotypic variation. While research at the interface of evolutionary and developmental biology has made excellent advances in understanding the contribution of genes to morphology, less well understood is the manner in which environmental cues are incorporated during development to influence the phenotype. Also virtually unexplored is how evolutionary transitions between environmental and genetic control of trait variation are achieved. Here, I review investigations into molecular mechanisms underlying phenotypic plasticity in the aphid wing dimorphism system. Among aphids, some species alternate between environmentally sensitive (polyphenic) and genetic (polymorphic) control of wing morph determination in their life cycle. Therefore, a traditional molecular genetic approach into understanding the genetically controlled polymorphism may provide a unique avenue into not only understanding the molecular basis of polyphenic variation in this group, but also the opportunity to compare and contrast the mechanistic basis of environmental and genetic control of similar dimorphisms.

Shifting the genomic gold standard for the prokaryotic species definition

DNA-DNA hybridization (DDH) has been used for nearly 50 years as the gold standard for prokaryotic species circumscriptions at the genomic level. It has been the only taxonomic method that offered a numerical and relatively stable species boundary, and its use has had a paramount influence on how the current classification has been constructed. However, now, in the era of genomics, DDH appears to be an outdated method for classification that needs to be substituted. The average nucleotide identity (ANI) between two genomes seems the most promising method since it mirrors DDH closely. Here we examine the work package JSpecies as a user-friendly, biologist-oriented interface to calculate ANI and the correlation of the tetranucleotide signatures between pairwise genomic comparisons. The results agreed with the use of ANI to substitute DDH, with a narrowed boundary that could be set at approximately 95-96%. In addition, the JSpecies package implemented the tetranucleotide signature correlation index, an alignment-free parameter that generally correlates with ANI and that can be of help in deciding when a given pair of organisms should be classified in the same species. Moreover, for taxonomic purposes, the analyses can be produced by simply randomly sequencing at least 20% of the genome of the query strains rather than obtaining their full sequence.

Fine-scale cospeciation between Brachycaudus and Buchnera aphidicola: bacterial genome helps define species and evolutionary relationships in aphids

Aphids harbour an obligatory symbiont, Buchnera aphidicola, providing essential amino acids not supplied by their diet. These bacteria are transmitted vertically and phylogenic analyses suggest that they have 'cospeciated' with their hosts. We investigated this cospeciation phenomenon at a fine taxonomic level, within the aphid genus Brachycaudus. We used DNA-based methods of species delimitation in both organisms, to avoid biases in the definition of aphid and Buchnera species and to infer association patterns without the presumption of a specific interaction. Our results call into question certain 'taxonomic' species of Brachycaudus and suggest that B. aphidicola has diversified into independently evolving entities, each specific to a 'phylogenetic' Brachycaudus species. We also found that Buchnera and their hosts simultaneously diversified, in parallel. These results validate the use of Buchnera DNA data for inferring the evolutionary history of their host. The Buchnera genome evolves rapidly, making it the perfect tool for resolving ambiguities in aphid taxonomy. This study also highlights the usefulness of species delimitation methods in cospeciation studies involving species difficult to conceptualize--as is the case for bacteria--and in cases in which the taxonomy of the interacting organisms has not been determined independently and species definition depends on host association.

Evolution of the secondary symbiont "Candidatus serratia symbiotica" in aphid species of the subfamily lachninae

Buchnera aphidicola BCc, the primary endosymbiont of the aphid Cinara cedri (subfamily Lachninae), is losing its symbiotic capacity and might be replaced by the coresident "Candidatus Serratia symbiotica." Phylogenetic and morphological analyses within the subfamily Lachninae indicate two different "Ca. Serratia symbiotica" lineages and support the longtime coevolution of both symbionts in C. cedri.

Genetic evidence from mitochondrial, nuclear, and endosymbiont markers for the evolution of host plant associated species in the aphid genus Hyalopterus (Hemiptera: Aphididae)

Over the past several decades biologists' fascination with plant-herbivore interactions has generated intensive research into the implications of these interactions for insect diversification. The study of closely related phytophagous insect species or populations from an evolutionary perspective can help illuminate ecological and selective forces that drive these interactions. Here we present such an analysis for aphids in the genus Hyalopterus (Hemiptera: Aphididae), a cosmopolitan group that feeds on plants in the genus Prunus (Rosaceae). Hyalopterus currently contains two recognized species associated with different Prunus species, although the taxonomy and evolutionary history of the group is poorly understood. Using mitochondrial COI sequences, 16S rDNA sequences from the aphid endosymbiont Buchnera aphidicola, and nine microsatellite loci we investigated population structure in Hyalopterus from the most commonly used Prunus host species throughout the Mediterranean as well as in California, where the species H. pruni is an invasive pest. We found three deeply divergent lineages structured in large part by specific associations with plum, almond, and peach trees. There was no evidence that geographic or temporal barriers could explain the overall diversity in the genus. Levels of genetic differentiation are consistent with that typically attributed to aphid species and indicate divergence times older than the domestication of Prunus for agriculture. Interestingly, in addition to their typical hosts, aphids from each of the three lineages were frequently found on apricot trees. Apricot also appears to act as a resource mediated hybrid zone for plum and almond associated lineages. Together, results suggest that host plants have played a role in maintaining host-associated differentiation in Hyalopterus for as long as several million years, despite worldwide movement of host plants and the potential for ongoing hybridization.

The bacterial species dilemma and the genomic-phylogenetic species concept

The number of species of Bacteria and Archaea (ca 5000) is surprisingly small considering their early evolution, genetic diversity and residence in all ecosystems. The bacterial species definition accounts in part for the small number of named species. The primary procedures required to identify new species of Bacteria and Archaea are DNA-DNA hybridization and phenotypic characterization. Recently, 16S rRNA gene sequencing and phylogenetic analysis have been applied to bacterial taxonomy. Although 16S phylogeny is arguably excellent for classification of Bacteria and Archaea from the Domain level down to the family or genus, it lacks resolution below that level. Newer approaches, including multilocus sequence analysis, and genome sequence and microarray analyses, promise to provide necessary information to better understand bacterial speciation. Indeed, recent data using these approaches, while meagre, support the view that speciation processes may occur at the subspecies level within ecological niches (ecovars) and owing to biogeography (geovars). A major dilemma for bacterial taxonomists is how to incorporate this new information into the present hierarchical system for classification of Bacteria and Archaea without causing undesirable confusion and contention. This author proposes the genomic-phylogenetic species concept (GPSC) for the taxonomy of prokaryotes. The aim is twofold. First, the GPSC would provide a conceptual and testable framework for bacterial taxonomy. Second, the GPSC would replace the burdensome requirement for DNA hybridization presently needed to describe new species. Furthermore, the GPSC is consistent with the present treatment at higher taxonomic levels.

Comparative and functional genomics of closteroviruses

The largest extant RNA genomes are found in two diverse families of positive-strand RNA viruses, the animal Coronaviridae and the plant Closteroviridae. Comparative analysis of the viruses from the latter family reveals three levels of gene conservation. The most conserved gene module defines RNA replication and is shared with plant and animal viruses in the alphavirus-like superfamily. A module of five genes that function in particle assembly and transport is a hallmark of the family Closteroviridae and was likely present in the ancestor of all three closterovirus genera. This module includes a homologue of Hsp70 molecular chaperones and three diverged copies of the capsid protein gene. The remaining genes show dramatic variation in their numbers, functions, and origins among closteroviruses within and between the genera. Proteins encoded by these genes include suppressors of RNA silencing, RNAse III, papain-like proteases, the AlkB domain implicated in RNA repair, Zn-ribbon-containing protein, and a variety of proteins with no detectable homologues in the current databases. The evolutionary processes that have shaped the complex and fluid genomes of the large RNA viruses might be similar to those that have been involved in evolution of genomic complexity in other divisions of life.

Biology bacteriocyte-associated endosymbionts of plant sap-sucking insects

Psyllids, whiteflies, aphids, and mealybugs are members of the suborder Sternorrhyncha and share a common property, namely the utilization of plant sap as their food source. Each of these insect groups has an obligatory association with a different prokaryotic endosymbiont, and the association is the result of a single infection followed by maternal, vertical transmission of the endosymbionts. The result of this association is the domestication of the free-living bacterium to serve the purposes of the host, namely the synthesis of essential amino acids. This domestication is probably in all cases accompanied by a major reduction in genome size. The different properties of the genomes and fragments of the genomes of these endosymbionts suggest that there are different constraints on the permissible evolutionary changes that are probably a function of the gene repertoire of the endosymbiont ancestor and the gene losses that occurred during the reduction of genome size.

Rarity associated with specific ecological niches in the bacterial world: the 'Synergistes' example

The 'Synergistes' group, which apparently represents an as yet unnamed division of the bacteria, was explored in 93 anaerobic environments (guts, soils, digestors, etc.). From 16S rDNA gene-targeted polymerase chain reaction (PCR) assays, this group appeared to be present in 90% of the anaerobic microbial ecosystems analysed. The phylogeny of 103 16S rDNA sequences from 30 ecosystems showed a strong link between 16S rDNA sequences and given ecosystems. 'Synergistes' 16S rDNA sequences from animal sources (termites, guinea pigs, pigs, birds, etc.) formed clustered phylogenetical groups. 'Synergistes' groups were also associated either with anaerobic digestors and soils or with thermophilic conditions. Sequences available from the DNA database were consistent with the results. These results show the wide diversity of the 'Synergistes' division as well as the specific ecological niche of each 16S rDNA sequences.

Cospeciation between the primary endosymbionts of mealybugs and their hosts

Mealybugs have an association with prokaryotic endosymbionts that are located in specialized cells called bacteriocytes. In order to compare the phylogeny of the host with that of the previously published phylogeny of the endosymbionts, 3.1 to 3.2 kilobase DNA fragments containing mitochondrial cytB (part), nd1,16S ribosomal DNA(rDNA), and 12S rDNA (part) were amplified and sequenced. A phylogenetic analysis of the data and a comparison with the trees obtained from endosymbiont genes and host 18S and 28S rDNA indicated that all the trees were similar. This result is consistent with an infection of a mealybug ancestor with a precursor of the endosymbiont followed by the vertical transmission of the endosymbiont to progeny. Comparison of the guanine + cytosine (G + C) contents of the mealybug mitochondrial genes with the same genes from other members of Sternorrhyncha and Arthropoda indicated that the mealybug genes had unusually low G + C contents in their DNAs (10.2 to 11.1 mol%).

Evolution of the leucine gene cluster in Buchnera aphidicola: insights from chromosomal versions of the cluster

In Buchnera aphidicola strains associated with the aphid subfamilies Thelaxinae, Lachninae, Pterocommatinae, and Aphidinae, the four leucine genes (leuA, -B, -C, and -D) are located on a plasmid. However, these genes are located on the main chromosome in B. aphidicola strains associated with the subfamilies Pemphiginae and Chaitophorinae. The sequence of the chromosomal fragment containing the leucine cluster and flanking genes has different positions in the chromosome in B. aphidicola strains associated with three tribes of the subfamily Pemphiginae and one tribe of the subfamily Chaitophorinae. Due to the extreme gene order conservation of the B. aphidicola genomes, the variability in the position of the leucine cluster in the chromosome may be interpreted as resulting from independent insertions from an ancestral plasmid-borne leucine gene. These findings do not support a chromosomal origin for the leucine genes in the ancestral B. aphidicola and do support a back transfer evolutionary scenario from a plasmid to the main chromosome.

Molecular systematics of aphids (Homoptera: Aphididae): new insights from the long-wavelength opsin gene

Viviparous aphids (Aphididae) constitute a monophyletic group within the Homoptera with more than 4000 extant species worldwide but higher diversity in temperate regions. Several aspects of their biology account for attention paid to this group of insects. Their plant-sap-sucking way of feeding with many species transmitting viruses to crop plants has important implications on crop management strategies. Cyclical parthenogenesis associated in many groups to host alternation and elaborate polyphenisms is of special interests for evolutionists. Finally, the ancient association of most aphid species with intracellular endosymbiotic bacteria (Buchnera sp.) has also received much attention from evolutionists interested in mechanisms involved in the symbiotic process. Knowing the phylogenetic relationships among major aphid taxa is of special interest to evolutionists interested in the above issues. However, until recently, molecular approaches to aphid phylogeny were absent and discussions on the evolution of aphid life-cycles and on evolutionary aspects of their symbiotic association with Buchnera were framed by morphology-based phylogenies. Recently, two reports using molecular approaches attempted to address the yet unresolved phylogeny of Aphididae with limited although somehow different conclusions. In the present report we study the utility of the long-wave opsin gene in resolving phylogenetic relationships among seven subfamilies within the Aphididae. Our results corroborate some previously proposed relationships and suggest a revision of some others. In particular, our data support grouping the analysed aphid species into three main clades, being the subfamily Lachninae one of them, which contradicts its generally accepted sistership relationship with the subfamily Aphidinae. Moreover, our data also suggest a basal position of Lachninae which has implications on current discussions about the ancestrality of conifer-feeding in modern aphids.

Detection and phylogenetic analysis of Wolbachia in Collembola

Wolbachia are obligatory, cytoplasmatically inherited alpha-Proteobacteria which are known for infecting the reproductive tissues of many arthropods. Their prevalence in the large group of Collembola, however, is not known, except for PCR detection in the parthenogenetically reproducing species Folsomia candida (Order: Entomobryomorpha; Family: Isotomidae). In this study, fluorescence in situ hybridization on microscopic sections of F. candida specimens indicated that Wolbachia-related bacteria were restricted to tissues of the ovary and brain. PCR with primers designed to detect 16S rRNA genes of Wolbachia were positive with specimens from all of five geographically independent F. candida breeding stocks and with three parthenogenetic species from another order (Poduromorpha; Family Tullbergiidae), i.e. Mesaphorura italica, M. macrochaeta and Paratullbergia callipygos. In contrast, negative results were obtained with the two sexually reproducing species, Isotoma viridis (Isotomidae) and Protaphorura fimata (Poduromorpha; Onychiuridae). The ftsZ gene of Wolbachia could be PCR-amplified from all Wolbachia-positive hosts with the exception of M. macrochaeta. The phylogenetic distances of the ftsZ and 16S rRNA gene sequences reflected the phylogenetic distances of the host organisms but the sequences of Wolbachia were relatively closely related, indicating that Wolbachia infections took place after the Collembola had diversified. Our study confirms a monophyletic branch (supergroup E) of Collembola colonizing Wolbachia and indicates that this group is a sister group of supergroup A, the latter harbouring a high diversity of host organisms within the group of insects.

A comparative analysis of morphological and ecological characters of European aphids and lycaenids in relation to ant attendance

Ants are a major environmental factor for many insect species. For example, aphids and lycaenids have evolved an array of associations with ants ranging from obligate myrmecophily to the avoidance of contact. Here we (1) analyze the predictive power of different ecological and morphological traits for explaining the strength of the association between ants and aphids/lycaenids and (2) contrast different taxonomic levels with respect to the variance explained by ant attendance. Data come from a literature survey including 112 species of aphids and 103 species of lycaenids from Europe. For aphids, feeding on woody plant parts is positively associated with ant attendance, while a high degree of mobility, feeding in isolation, and the possession of wings in the adult stage are negatively associated with ant attendance. In lycaenids, feeding on inflorescences and feeding on Fabaceae host plants is closely associated with ant attendance, while living in forests bears a smaller likelihood to establish mutualistic relationships. Body size always appeared to be a poor predictor for the degree of ant attendance. Overall, in both insect groups less than 10% of the variation in the ecological traits recorded is explained by the different modes of ant association. When decomposing the variance in traits explained by ant attendance at different taxonomic levels, aphids and lycaenids show contrasting results. In aphids, most variance in the degree of ant attendance is explained at the subfamily level and least at the species level. The opposite is true for lycaenids, where most variance is explained at the lowest taxonomic level. Possible mechanisms explaining these different patterns of associations with ants are suggested.

Side-stepping secondary symbionts: widespread horizontal transfer across and beyond the Aphidoidea

To elucidate the co-evolutionary relationships between phloem-feeding insects and their secondary, or facultative, bacterial symbionts, we explore the distributions of three such microbes--provisionally named the R-type (or PASS, or S-sym), T-type (or PABS), and U-type--across a number of aphid and psyllid hosts through the use of diagnostic molecular screening techniques and DNA sequencing. Although typically maternally transmitted, phylogenetic and pairwise divergence analyses reveal that these bacteria have been independently acquired by a variety of unrelated insect hosts, indicating that horizontal transfer has helped to shape their distributions. Based on the high genetic similarity between symbionts in different hosts, we argue that transfer events have occurred recently on an evolutionary timescale. In several instances, however, closely related symbionts associate with related hosts, suggesting that horizontal transfer between distant relatives may be rarer than transmission between close relatives. Our findings on the prevalence of these symbionts within many aphid taxa, along with published observations concerning their effects on host fitness, imply a significant role of facultative symbiosis in aphid ecology and evolution.

Molecular characterization of the leucine cluster in Buchnera sp. strain PSY, a primary endosymbiont of the aphid Pemphigus spyrothecae

Buchnera strains from most aphid subfamilies studied to date have been found to carry the leucine gene cluster (leuA, -B, -C, and -D) on a plasmid, an organization unique among bacteria. Here, however, we demonstrate a classical chromosomal location of the cluster in Buchnera sp. strain PSY from the aphid Pemphigus spyrothecae (subfamily Pemphiginae). The genes that flank leuABCD in Buchnera sp. strain PSY appear to be adjacent in the genome of Buchnera sp. strain APS, a strain carrying a leucine plasmid. We propose that the presence of a leucine plasmid predates the diversification of symbiotic Buchnera and that the chromosomal location observed in Buchnera sp. strain PSY arose by a transfer of the leucine genes from a plasmid to the chromosome.