Reconstructing the phylogeny of aphids (Hemiptera: Aphididae) using DNA of the obligate symbiont Buchnera aphidicola

Taxonomic Revision of the Nearctic Genus Drepanaphis Del Guercio (Hemiptera, Aphididae: Drepanosiphinae)

The Nearctic aphid genus Drepanaphis Del Guercio, 1909, the largest within the subfamily Drepanosiphinae (Hemiptera: Aphididae), is characterised by distinctive dorsal abdominal tubercles. This study presents a comprehensive taxonomic revision of the genus, expanding the recognised species to 18, including the newly described Drepanaphis robinsoni Malik sp. nov. Detailed descriptions and figures for 44 morphs, encompassing alate viviparous females, oviparous females and males, are provided, with new identification keys for all known species and morphs. The sexual morphs of 15 species, particularly oviparous females, are documented for the first time. Morphometric and principal component analyses (PCA) are employed to distinguish the studied taxa. This study identifies and corrects numerous misidentifications in museum collections, previously labelled as D. acerifoliae, D. choanotricha, D. kanzensis, D. knowltoni, D. parva, D. sabrinae or D. tissoti. Furthermore, it revalidates the distinct status of D. nigricans and D. tissoti, which had been synonymised in earlier works. Current range maps for all species and images of key morphological features obtained through light and scanning electron microscopy are also presented, providing a more complete understanding of this understudied genus.

'Drifting' Buchnera genomes track the microevolutionary trajectories of their aphid hosts

Evolution of Buchnera-aphid host symbioses is often studied among species at macroevolutionary scales. Investigations within species offer a different perspective about how eco-evolutionary processes shape patterns of genetic variation at microevolutionary scales. Our study leverages new and publicly available whole-genome sequencing data to study Buchnera-aphid host evolution in Myzus persicae, the peach potato aphid, a globally invasive and polyphagous pest. Across 43 different asexual, clonally reproducing isofemale strains, we examined patterns of genomic covariation between Buchnera and their aphid host and considered the distribution of mutations in protein-coding regions of the Buchnera genome. We found Buchnera polymorphisms within aphid strains, suggesting the presence of genetically different Buchnera strains within the same clonal lineage. Genetic distance between pairs of Buchnera samples was positively correlated to genetic distance between their aphid hosts, indicating shared evolutionary histories. However, there was no segregation of genetic variation for both M. persicae and Buchnera with plant host (Brassicaceae and non-tobacco Solanaceae) and no associations between genetic and geographic distance at global or regional spatial scales. Abundance patterns of non-synonymous mutations were similar to synonymous mutations in the Buchnera genome, and both mutation classes had similar site frequency spectra. We hypothesize that a predominance of neutral processes results in the Buchnera of M. persicae to simply 'drift' with the evolutionary trajectory of their aphid hosts. Our study presents a unique microevolutionary characterization of Buchnera-aphid host genomic covariation across multiple aphid clones. This provides a new perspective on the eco-evolutionary processes generating and maintaining polymorphisms in a major pest aphid species and its obligate primary endosymbiont.

Repeated horizontal acquisition of lagriamide-producing symbionts in Lagriinae beetles

Microbial symbionts associate with multicellular organisms on a continuum from facultative associations to mutual codependency. In some of the oldest intracellular symbioses there is exclusive vertical symbiont transmission, and co-diversification of symbiotic partners over millions of years. Such symbionts often undergo genome reduction due to low effective population sizes, frequent population bottlenecks, and reduced purifying selection. Here, we describe multiple independent acquisition events of closely related defensive symbionts followed by genome erosion in a group of Lagriinae beetles. Previous work in Lagria villosa revealed the dominant genome-eroded symbiont of the genus Burkholderia produces the antifungal compound lagriamide and protects the beetle's eggs and larvae from antagonistic fungi. Here, we use metagenomics to assemble 11 additional genomes of lagriamide-producing symbionts from seven different host species within Lagriinae from five countries, to unravel the evolutionary history of this symbiotic relationship. In each host species, we detected one dominant genome-eroded Burkholderia symbiont encoding the lagriamide biosynthetic gene cluster (BGC). Surprisingly, however, we did not find evidence for host-symbiont co-diversification, or for a monophyly of the lagriamide-producing symbionts. Instead, our analyses support at least four independent acquisition events of lagriamide-encoding symbionts and subsequent genome erosion in each of these lineages. By contrast, a clade of plant-associated relatives retained large genomes but secondarily lost the lagriamide BGC. In conclusion, our results reveal a dynamic evolutionary history with multiple independent symbiont acquisitions characterized by high degree of specificity. They highlight the importance of the specialized metabolite lagriamide for the establishment and maintenance of this defensive symbiosis.

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.

Whole-genome sequence of the Cooley spruce gall adelgid, Adelges cooleyi (Hemiptera: Sternorrhyncha: Adelgidae)

The adelgids (Adelgidae) are a small family of sap-feeding insects, which, together with true aphids (Aphididae) and phylloxerans (Phylloxeridae), make up the infraorder Aphidomorpha. Some adelgid species are highly destructive to forest ecosystems such as Adelges tsugae, Adelges piceae, Adelges laricis, Pineus pini, and Pineus boerneri. Despite this, there are no high-quality genomic resources for adelgids, hindering advanced genomic analyses within Adelgidae and among Aphidomorpha. Here, we used PacBio continuous long-read and Illumina RNA-sequencing to construct a high-quality draft genome assembly for the Cooley spruce gall adelgid, Adelges cooleyi (Gillette), a gall-forming species endemic to North America. The assembled genome is 270.2 Mb in total size and has scaffold and contig N50 statistics of 14.87 and 7.18 Mb, respectively. There are 24,967 predicted coding sequences, and the assembly completeness is estimated at 98.1 and 99.6% with core BUSCO gene sets of Arthropoda and Hemiptera, respectively. Phylogenomic analysis using the A. cooleyi genome, 3 publicly available adelgid transcriptomes, 4 phylloxera transcriptomes, the Daktulosphaira vitifoliae (grape phylloxera) genome, 4 aphid genomes, and 2 outgroup coccoid genomes fully resolves adelgids and phylloxerans as sister taxa. The mitochondrial genome is 24 kb, among the largest in insects sampled to date, with 39.4% composed of noncoding regions. This genome assembly is currently the only genome-scale, annotated assembly for adelgids and will be a valuable resource for understanding the ecology and evolution of Aphidomorpha.

Parallel Losses of Blue Opsin Correlate with Compensatory Neofunctionalization of UV-Opsin Gene Duplicates in Aphids and Planthoppers

Expanding on previous efforts to survey the visual opsin repertoires of the Hemiptera, this study confirms that homologs of the UV- and LW-opsin subfamilies are conserved in all Hemiptera, while the B-opsin subfamily is missing from the Heteroptera and subgroups of the Sternorrhyncha and Auchenorrhyncha, i.e., aphids (Aphidoidea) and planthoppers (Fulgoroidea), respectively. Unlike in the Heteroptera, which are characterized by multiple independent expansions of the LW-opsin subfamily, the lack of B-opsin correlates with the presence of tandem-duplicated UV-opsins in aphids and planthoppers. Available data on organismal wavelength sensitivities and retinal gene expression patterns lead to the conclusion that, in both groups, one UV-opsin paralog shifted from ancestral UV peak sensitivity to derived blue sensitivity, likely compensating for the lost B-opsin. Two parallel bona fide tuning site substitutions compare to 18 non-corresponding amino acid replacements in the blue-shifted UV-opsin paralogs of aphids and planthoppers. Most notably, while the aphid blue-shifted UV-opsin clade is characterized by a replacement substitution at one of the best-documented UV/blue tuning sites (Rhodopsin site 90), the planthopper blue-shifted UV-opsin paralogs retained the ancestral lysine at this position. Combined, the new findings identify aphid and planthopper UV-opsins as a new valuable data sample for studying adaptive opsin evolution.

Comparative transcriptomics of aphid species that diverged > 22 MYA reveals genes that are important for the maintenance of their symbiosis

Most plant-sap feeding insects have obligate relationships with maternally transmitted bacteria. Aphids require their nutritional endosymbiont, Buchnera aphidicola, for the production of essential amino acids. Such endosymbionts are harbored inside of specialized insect cells called bacteriocytes. Here, we use comparative transcriptomics of bacteriocytes between two recently diverged aphid species, Myzus persicae and Acyrthosiphon pisum, to identify key genes that are important for the maintenance of their nutritional mutualism. The majority of genes with conserved expression profiles in M. persicae and A. pisum are for orthologs previously identified in A. pisum to be important for the symbiosis. However, asparaginase which produces aspartate from asparagine was significantly up-regulated only in A. pisum bacteriocytes, potentially because Buchnera of M. persicae encodes its own asparaginase enzyme unlike Buchnera of A. pisum resulting in Buchnera of A. pisum to be dependent on its aphid host for aspartate. One-to-one orthologs that explained the most amount of variation for bacteriocyte specific mRNA expression for both species includes a collaborative gene for methionine biosynthesis, multiple transporters, a horizontally transmitted gene, and secreted proteins. Finally, we highlight species-specific gene clusters which may contribute to host adaptations and/or accommodations in gene regulation to changes in the symbiont or the symbiosis.

Microbiome of the Aphid Genus Dysaphis Börner (Hemiptera: Aphidinae) and Its Relation to Ant Attendance

Among mutualistic relationships of aphids with other organisms, there are two that seem to be of major importance: trophobiosis with ants and endosymbiosis of bacteria. While the former is well studied, the latter is the subject of an increasing amount of research constantly revealing new aspects of this symbiosis. Here, we studied the possible influence of ant attendance on the composition of aphid microbiota on primary and secondary hosts exploited by the aphid genus Dysaphis. The microbiome of 44 samples representing 12 aphid species was studied using an Illumina HiSeq 4000 with the V3-V4 region of 16S rRNA. The results showed a higher abundance of common facultative symbionts (Serratia, Regiella, Fukatsuia) in aphid species unattended by ants, but also on secondary hosts. However, in colonies attended by ants, the general species composition of bacterial symbionts was more rich in genera than in unattended colonies (Wolbachia, Gilliamella, Spiroplasma, Sphingomonas, Pelomonas). The results indicate a huge variability of facultative symbionts without clear correlation with ant attendance or aphid species. The possibility of multiple routes of bacterial infection mediated by ant-made environmental conditions is discussed.

Molecular and Developmental Signatures of Genital Size Macro-Evolution in Bugs

Our understanding of the genetic architecture of phenotypic traits has experienced drastic growth over the last years. Nevertheless, the majority of studies associating genotypes and phenotypes have been conducted at the ontogenetic level. Thus, we still have an elusive knowledge of how these genetic-developmental architectures evolve themselves and how their evolution is mirrored in the phenotypic change across evolutionary time. We tackle this gap by reconstructing the evolution of male genital size, one of the most complex traits in insects, together with its underlying genetic architecture. Using the order Hemiptera as a model, spanning over 350 million years of evolution, we estimate the correlation between genitalia and three features: development rate, body size, and rates of DNA substitution in 68 genes associated with genital development. We demonstrate that genital size macro-evolution has been largely dependent on body size and weakly influenced by development rate and phylogenetic history. We further revealed significant correlations between mutation rates and genital size for 19 genes. Interestingly, these genes have diverse functions and participate in distinct signaling pathways, suggesting that genital size is a complex trait whose fast evolution has been enabled by molecular changes associated with diverse morphogenetic processes. Our data further demonstrate that the majority of DNA evolution correlated with the genitalia has been shaped by negative selection or neutral evolution. Thus, in terms of sequence evolution, changes in genital size are predominantly facilitated by relaxation of constraints rather than positive selection, possibly due to the high pleiotropic nature of the morphogenetic genes.

Geographic isolation drives speciation in Nearctic aphids

Across herbivorous insect clades, species richness and host-use diversity tend to positively covary. This could be because host-use divergence drives speciation, or because it raises the ecological limits on species richness. To evaluate these hypotheses, we performed phylogenetic path model analyses of the species diversity of Nearctic aphids. Here, we show that variation in the species richness of aphid clades is caused mainly by host-use divergence, whereas variation in speciation rates is caused more by divergence in non-host-related niche variables. Aphid speciation is affected by both the evolution of host and non-host-related niche components, but the former is largely caused by the latter. Thus, our analyses suggest that host-use divergence can both raise the ecological limits on species richness and drive speciation, although in the latter case, host-use divergence tends to be a step along the causal path leading from non-host-related niche evolution to speciation.

Using phylogenetic path modelling, geographic isolation is revealed to be the main driver of Nearctic aphid speciation, with subsequent host-use divergence raising the limits on clade species richness.

Endosymbiotic Bacterial Diversity of Corn Leaf Aphid, Rhopalosiphum maidis Fitch (Hemiptera: Aphididae) Associated with Maize Management Systems

In this study, different maize fields cultivated under different management systems were sampled to test corn leaf aphid, Rhopalosiphum maidis, populations in terms of total and endosymbiotic bacterial diversity. Corn leaf aphid natural populations were collected from traditionally managed maize fields grown under high agricultural and natural landscape diversity as well as conventionally treated high-input agricultural fields grown in monoculture and with fertilizers use, hence with low natural landscape diversity. Total bacterial community assessment by DNA sequencing was performed using the Illumina MiSeq platform. In total, 365 bacterial genera were identified and 6 endosymbiont taxa. A high abundance of the primary endosymbiont Buchnera and secondary symbionts Serratia and Wolbachia were detected in all maize crops. Their frequency was found to be correlated with the maize management system used, probably with fertilizer input. Three other facultative endosymbionts (“Candidatus Hamiltonella”, an uncultured Rickettsiales genus, and Spiroplasma) were also recorded at different frequencies under the two management regimes. Principal components analyses revealed that the relative contribution of the obligate and dominant symbiont Buchnera to the aphid endosymbiotic bacterial community was 72%, whereas for the managed system this was only 16.3%. When facultative symbionts alone were considered, the effect of management system revealed a DNA diversity of 23.3%.

Evolution of DNA Methylation Across Ecdysozoa

DNA methylation is a crucial, abundant mechanism of gene regulation in vertebrates. It is less prevalent in many other metazoan organisms and completely absent in some key model species, such as Drosophila melanogaster and Caenorhabditis elegans. We report here a comprehensive study of the presence and absence of DNA methyltransferases (DNMTs) in 138 Ecdysozoa, covering Arthropoda, Nematoda, Priapulida, Onychophora, and Tardigrada. Three of these phyla have not been investigated for the presence of DNA methylation before. We observe that the loss of individual DNMTs independently occurred multiple times across ecdysozoan phyla. We computationally predict the presence of DNA methylation based on CpG rates in coding sequences using an implementation of Gaussian Mixture Modeling, MethMod. Integrating both analysis we predict two previously unknown losses of DNA methylation in Ecdysozoa, one within Chelicerata (Mesostigmata) and one in Tardigrada. In the early-branching Ecdysozoa Priapulus caudatus, we predict the presence of a full set of DNMTs and the presence of DNA methylation. We are therefore showing a very diverse and independent evolution of DNA methylation in different ecdysozoan phyla spanning a phylogenetic range of more than 700 million years.

OUP accepted manuscript

Homology of highly divergent genes often cannot be determined from sequence similarity alone. For example, we recently identified in the aphid Hormaphis cornu a family of rapidly evolving bicycle genes, which encode novel proteins implicated as plant gall effectors, and sequence similarity search methods yielded few putative bicycle homologs in other species. Coding sequence-independent features of genes, such as intron-exon boundaries, often evolve more slowly than coding sequences, however, and can provide complementary evidence for homology. We found that a linear logistic regression classifier using only structural features of bicycle genes identified many putative bicycle homologs in other species. Independent evidence from sequence features and intron locations supported homology assignments. To test the potential roles of bicycle genes in other aphids, we sequenced the genome of a second gall-forming aphid, Tetraneura nigriabdominalis and found that many bicycle genes are strongly expressed in the salivary glands of the gall forming foundress. In addition, bicycle genes are strongly overexpressed in the salivary glands of a non-gall forming aphid, Acyrthosiphon pisum, and in the non-gall forming generations of H. cornu. These observations suggest that Bicycle proteins may be used by multiple aphid species to manipulate plants in diverse ways. Incorporation of gene structural features into sequence search algorithms may aid identification of deeply divergent homologs, especially of rapidly evolving genes involved in host-parasite interactions.

Phylogeny of drepanosiphine aphids sensu lato (Hemiptera, Aphidoidea) inferred from molecular and morphological data

As the second largest and most diverse group in the superfamily Aphidoidea, the phylogeny of drepanosiphine aphids sensu lato (s.l.) is critical for discussing the evolution of aphids. However, the taxa composition and phylogenetic relationships of drepanosiphine aphids s.l. have not been fully elucidated to date. In this study, based on total-evidence analyses combining 4 molecular genes (3 mitochondrial, COI, tRNA-Leu/COII, and CytB; 1 nuclear, EF-1ɑ) and 64 morphological and biological characteristics, the phylogeny of this group was reconstructed for the first time at the subfamily level using different datasets, parsimonies and model-based methods. All of our phylogenetic inferences clearly indicated that the drepanosiphine aphids s.l. was not a monophyletic group and seemed to support the division of the drepanosiphine aphids s.l. into different groups classified at the subfamily level. Calaphidinae was also not a monophyletic group, and Saltusaphidinae was nested within this subfamily. Drepanosiphinae was not clustered with Chaitophorinae, which was inconsistent with the previous hypothesis of a close relationship between them, illustrating that their phylogeny remains controversial. Overall, some groups of drepanosiphine aphids s.l., including Phyllaphidinae, Macropodaphidinae, Pterastheniinae, Lizeriinae, Drepanosiphinae, Spicaphidinae, Saltusaphidinae, and Calaphidinae, clustered together and might constitute the actual drepanosiphine aphids s.l. To a certain extent, our results clarified the phylogenetic relationships among drepanosiphine aphids s.l. and confirmed their taxonomic status as subfamilies.

SLI1 confers broad-spectrum resistance to phloem-feeding insects

Resistance (R) genes usually compete in a coevolutionary arms race with reciprocal effectors to confer strain-specific resistance to pathogens or herbivorous insects. Here, we investigate the specificity of SLI1, a recently identified R gene in Arabidopsis that encodes a small heat shock-like protein involved in resistance to Myzus persicae aphids. In a panel with several aphid and whitefly species, SLI1 compromised reproductive rates of three species: the tobacco aphid M. persicae nicotianae, the cabbage aphid Brevicoryne brassicae and the cabbage whitefly Aleyrodes proletella. Electrical penetration graph recording of aphid behaviour, revealed shorter salivations and a 3-to-5-fold increase in phloem feeding on sli1 loss-of-function plants. The mustard aphid Lipaphis erysimi and Bemisia tabaci whitefly were not affected by SLI1. Unlike the other two aphid species, L. erysimi exhibited repetitive salivations preceding successful phloem feeding, indicating a role of salivary effectors in overcoming SLI1-mediated resistance. Microscopic characterization showed that SLI1 proteins localize in the sieve tubes of virtually all above- and below-ground tissues and co-localize with the aphid stylet tip after penetration of the sieve element plasma membrane. These observations reveal an unconventional R gene that escapes the paradigm of strain specificity and confers broad-spectrum quantitative resistance to phloem-feeding insects.

Natural selection shapes maintenance of orthologous sRNAs in divergent host-restricted bacterial genomes

Historically it has been difficult to study the evolution of bacterial small RNAs (sRNAs) across distantly related species. For example, identifying homologs of sRNAs is often difficult in genomes that have undergone multiple structural rearrangements. Also, some types of regulatory sRNAs evolve at rapid rates. The high degree of genomic synteny among divergent host-restricted bacterial lineages, including intracellular symbionts, is conducive to sRNA maintenance and homolog identification. In turn, symbiont genomes can provide us with novel insights into sRNA evolution. Here, we examine the sRNA expression profile of the obligate symbiont of psyllids, Carsonella ruddii, which has one of the smallest cellular genomes described. Using RNA-seq, we identified 36 and 32 antisense sRNAs (asRNAs) expressed by Carsonella from the psyllids Bactericera cockerelli (Carsonella-BC) and Diaphorina citri (Carsonella-DC), respectively. The majority of these asRNAs were associated with genes that are involved in essential amino acid biosynthetic pathways. Eleven of the asRNAs were conserved in both Carsonella lineages and the majority were maintained by selection. Notably, five of the corresponding coding sequences are also the targets of conserved asRNAs in a distantly related insect symbiont, Buchnera. We detected differential expression of two asRNAs for genes involved in arginine and leucine biosynthesis occurring between two distinct Carsonella-BC life stages. Using asRNAs identified in Carsonella, Buchnera, and Profftella which are all endosymbionts, and Escherichia coli, we determined that regions upstream of these asRNAs encode unique conserved patterns of AT/GC richness, GC skew, and sequence motifs which may be involved in asRNA regulation.

Microbiomes of willow-galling sawflies: effects of host plant, gall type, and phylogeny on community structure and function

While free-living herbivorous insects are thought to harbor microbial communities composed of transient bacteria derived from their diet, recent studies indicate that insects that induce galls on plants may be involved in more intimate host-microbe relationships. We used 16S rDNA metabarcoding to survey larval microbiomes of 20 nematine sawfly species that induce bud or leaf galls on 13 Salix species. The 391 amplicon sequence variants (ASVs) detected represented 69 bacterial genera in six phyla. Multi-variate statistical analyses showed that the structure of larval microbiomes is influenced by willow host species as well as by gall type. Nevertheless, a "core" microbiome composed of 58 ASVs is shared widely across the focal galler species. Within the core community, the presence of many abundant, related ASVs representing multiple distantly related bacterial taxa is reflected as a statistically significant effect of bacterial phylogeny on galler-microbe associations. Members of the core community have a variety of inferred functions, including degradation of phenolic compounds, nutrient supplementation, and production of plant hormones. Hence, our results support suggestions of intimate and diverse interactions between galling insects and microbes and add to a growing body of evidence that microbes may play a role in the induction of insect galls on plants.

Evolutionary Rates are Correlated Between Buchnera Endosymbionts and the Mitochondrial Genomes of Their Aphid Hosts

The evolution of bacterial endosymbiont genomes is strongly influenced by host-driven selection. Factors affecting host genome evolution will potentially affect endosymbiont genomes in similar ways. One potential outcome is correlations in molecular rates between the genomes of the symbiotic partners. Recently, we presented the first evidence of such correlations between the mitochondrial genomes of cockroaches and the genomes of their endosymbiont (Blattabacterium cuenoti). Here we investigate whether similar patterns are found in additional host-symbiont partners. We use partial genome data from multiple strains of the bacterial endosymbionts Buchnera aphidicola and Sulcia muelleri, and the mitochondrial genomes of their sap-feeding insect hosts. Both endosymbionts show phylogenetic congruence with the mitochondria of their hosts, a result that is expected due to their identical mode of inheritance. We compared root-to-tip distances and branch lengths of phylogenetically independent species pairs. Both analyses showed a highly significant correlation of molecular rates between the genomes of Buchnera and the mitochondrial genomes of their hosts. A similar correlation was detected between Sulcia and their hosts, but was not statistically significant. Our results indicate that evolutionary rate correlations between hosts and long-term symbionts may be a widespread phenomenon.

First Report on the Acrobat Ant Crematogaster scutellaris Storing Live Aphids in Its Oak-Gall Nests

Simple Summary

Galls represent an amazing microcosm which contains a variety of multiple interactions among different actors, and therefore, offers the opportunity to observe and investigate phenomena belonging to different areas of biology: from the development process, connected to the interaction between the galligenous agent and the host plant, to the moment of their colonization by different species, since some ants may provide defense against pathogens, certain phytophagous insects or favor mutualists. In the present work we describe some aspects of oak-gall colonization by different ant species, highlighting how the gall’s height on the plant influences ant colonization and how different ant species produce different nest architectures. The most relevant aspect, however, is the discovery of a novel ant-aphid relationship: the transport of living aphids into oak-gall nests. We found no evidence of immediate predation of these aphids inside the galls, so they are likely stored to overwinter due to a mutualistic relationship and/or serve as food storage. This is not only an interesting report on the mutualisms involving ants and their insect partners, but it may also have important consequences on the aphids’ phenology with the host plants. Once more, ants show their relevant impact on multitrophic interactions and ecosystem dynamics.

Abstract

This study provides new data about the role of ants in mutualistic interactions with aphids mediated by galls. We focused our investigation on galls induced by the cynipid Andricus kollari by conducting a survey and a subsequent experiment in an Italian oak forest. The ants Crematogaster scutellaris, Colobopsis truncata and Temnothorax italicus frequently used the galls as nests: Crematogaster scutellaris occupied galls which were located higher on the oak trees, while C. truncata and T. italicus were located at lower positions. In addition, galls occupied by C. scutellaris showed varied internal architecture in relation to the colony composition. Importantly, field surveys revealed for the first time that C. scutellaris nest galls also contained live individuals of the non-galligenous aphid Panaphis juglandis. Field experiments suggested that the ants actively seek, collect and stock live aphids. No signs of predation and injuries were detected on the stored aphids, which were probably kept for safe overwintering, though we cannot exclude a possible occasional use as food. This report reveals a possible novel relationship which could have important consequences on the phenology and presence of aphids on the host plant.

The First Complete Mitochondrial Genome of Lachninae Species and Comparative Genomics Provide New Insights into the Evolution of Gene Rearrangement and the Repeat Region

Complete mitochondrial genomes are valuable resources for different research fields such as genomics, molecular evolution and phylogenetics. The subfamily Lachninae represents one of the most ancient evolutionary lineages of aphids. To date, however, no complete Lachninae mitogenome is available in public databases. Here we report the Stomaphis sinisalicis mitogenome, representing the first complete mitogenome of Lachninae. The S. sinisalicis mitogenome is consist of 13 protein-coding genes (PCGs), two rRNA genes (rRNAs), 22 tRNA genes (tRNAs), a control region and a large tandem repeat region. Strikingly, the mitogenome exhibits a novel, highly rearranged gene order between trnE and nad1 compared with that of other aphids. The presence of repeat region in the basal Lachninae may further indicate it is probably an ancestral feature of aphid mitogenomes. Collectively, this study provides new insights on mitogenome evolution and valuable data for future comparative studies across different insect lineages.

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.

A chromosome-level genome assembly of the woolly apple aphid, Eriosoma lanigerum Hausmann (Hemiptera: Aphididae)

Woolly apple aphid (WAA, Eriosoma lanigerum Hausmann) (Hemiptera: Aphididae) is a major pest of apple trees (Malus domestica, order Rosales) and is critical to the economics of the apple industry in most parts of the world. Here, we generated a chromosome-level genome assembly of WAA-representing the first genome sequence from the aphid subfamily Eriosomatinae-using a combination of 10X Genomics linked-reads and in vivo Hi-C data. The final genome assembly is 327 Mb, with 91% of the assembled sequences anchored into six chromosomes. The contig and scaffold N50 values are 158 kb and 71 Mb, respectively, and we predicted a total of 28,186 protein-coding genes. The assembly is highly complete, including 97% of conserved arthropod single-copy orthologues based on Benchmarking Universal Single-Copy Orthologs (busco) analysis. Phylogenomic analysis of WAA and nine previously published aphid genomes, spanning four aphid tribes and three subfamilies, reveals that the tribe Eriosomatini (represented by WAA) is recovered as a sister group to Aphidini + Macrosiphini (subfamily Aphidinae). We identified syntenic blocks of genes between our WAA assembly and the genomes of other aphid species and find that two WAA chromosomes (El5 and El6) map to the conserved Macrosiphini and Aphidini X chromosome. Our high-quality WAA genome assembly and annotation provides a valuable resource for research in a broad range of areas such as comparative and population genomics, insect-plant interactions and pest resistance management.

Whole genome sequence of the soybean aphid, Aphis glycines

Aphids are emerging as model organisms for both basic and applied research. Of the 5,000 estimated species, only three aphids have published whole genome sequences: the pea aphid Acyrthosiphon pisum, the Russian wheat aphid, Diuraphis noxia, and the green peach aphid, Myzus persicae. We present the whole genome sequence of a fourth aphid, the soybean aphid (Aphis glycines), which is an extreme specialist and an important invasive pest of soybean (Glycine max). The availability of genomic resources is important to establish effective and sustainable pest control, as well as to expand our understanding of aphid evolution. We generated a 302.9 Mbp draft genome assembly for Ap. glycines using a hybrid sequencing approach. This assembly shows high completeness with 19,182 predicted genes, 92% of known Ap. glycines transcripts mapping to contigs, and substantial continuity with a scaffold N50 of 174,505 bp. The assembly represents 95.5% of the predicted genome size of 317.1 Mbp based on flow cytometry. Ap. glycines contains the smallest known aphid genome to date, based on updated genome sizes for 19 aphid species. The repetitive DNA content of the Ap. glycines genome assembly (81.6 Mbp or 26.94% of the 302.9 Mbp assembly) shows a reduction in the number of classified transposable elements compared to Ac. pisum, and likely contributes to the small estimated genome size. We include comparative analyses of gene families related to host-specificity (cytochrome P450's and effectors), which may be important in Ap. glycines evolution. This Ap. glycines draft genome sequence will provide a resource for the study of aphid genome evolution, their interaction with host plants, and candidate genes for novel insect control methods.

MinYS: mine your symbiont by targeted genome assembly in symbiotic communities

Most metazoans are associated with symbionts. Characterizing the effect of a particular symbiont often requires getting access to its genome, which is usually done by sequencing the whole community. We present MinYS, a targeted assembly approach to assemble a particular genome of interest from such metagenomic data. First, taking advantage of a reference genome, a subset of the reads is assembled into a set of backbone contigs. Then, this draft assembly is completed using the whole metagenomic readset in a de novo manner. The resulting assembly is output as a genome graph, enabling different strains with potential structural variants coexisting in the sample to be distinguished. MinYS was applied to 50 pea aphid resequencing samples, with variable diversity in symbiont communities, in order to recover the genome sequence of its obligatory bacterial symbiont, Buchnera aphidicola. It was able to return high-quality assemblies (one contig assembly in 90% of the samples), even when using increasingly distant reference genomes, and to retrieve large structural variations in the samples. Because of its targeted essence, it outperformed standard metagenomic assemblers in terms of both time and assembly quality.

Soybean aphid biotype 1 genome: Insights into the invasive biology and adaptive evolution of a major agricultural pest

The soybean aphid, Aphis glycines Matsumura (Hemiptera: Aphididae) is a serious pest of the soybean plant, Glycine max, a major world-wide agricultural crop. We assembled a de novo genome sequence of Ap. glycines Biotype 1, from a culture established shortly after this species invaded North America. 20.4% of the Ap. glycines proteome is duplicated. These in-paralogs are enriched with Gene Ontology (GO) categories mostly related to apoptosis, a possible adaptation to plant chemistry and other environmental stressors. Approximately one-third of these genes show parallel duplication in other aphids. But Ap. gossypii, its closest related species, has the lowest number of these duplicated genes. An Illumina GoldenGate assay of 2380 SNPs was used to determine the world-wide population structure of Ap. Glycines. China and South Korean aphids are the closest to those in North America. China is the likely origin of other Asian aphid populations. The most distantly related aphids to those in North America are from Australia. The diversity of Ap. glycines in North America has decreased over time since its arrival. The genetic diversity of Ap. glycines North American population sampled shortly after its first detection in 2001 up to 2012 does not appear to correlate with geography. However, aphids collected on soybean Rag experimental varieties in Minnesota (MN), Iowa (IA), and Wisconsin (WI), closer to high density Rhamnus cathartica stands, appear to have higher capacity to colonize resistant soybean plants than aphids sampled in Ohio (OH), North Dakota (ND), and South Dakota (SD). Samples from the former states have SNP alleles with high F_ST values and frequencies, that overlap with genes involved in iron metabolism, a crucial metabolic pathway that may be affected by the Rag-associated soybean plant response. The Ap. glycines Biotype 1 genome will provide needed information for future analyses of mechanisms of aphid virulence and pesticide resistance as well as facilitate comparative analyses between aphids with differing natural history and host plant range.

Phylogeny of the Subtribe Monoctonina (Hymenoptera, Braconidae, Aphidiinae)

Members of the Monoctonina subtribe have long been neglected in applied studies of the subfamily Aphidiinae, due to their low economic importance, as they do not parasitize pests of cultivated plants. Consequently, data about this group are scarce, including its taxonomy and phylogeny. In the present study, we explore inter- and intraspecific genetic variation of Monoctonina species, including genera Monoctonus Haliday 1833, Monoctonia Starý 1962, Falciconus Mackauer 1959 and Harkeria Cameron 1900. We employ two molecular markers, the barcode region of the mitochondrial cytochrome c oxidase subunit I (COI) and the D2 region of the 28S nuclear gene (28S rDNA), to analyze genetic structuring and phylogeny of all available Monoctonina species, and combine them with morphological data for an integrative approach. We report one new species, and three potentially new species which can be formally described when further specimens are available. Analysis of phylogenetic relationships within the subtribe shows a basal position for the genera Falciconus and Monoctonia, and the close relatedness of Harkeria and Monoctonus.

Phylogenomics Identifies an Ancestral Burst of Gene Duplications Predating the Diversification of Aphidomorpha

Aphids (Aphidoidea) are a diverse group of hemipteran insects that feed on plant phloem sap. A common finding in studies of aphid genomes is the presence of a large number of duplicated genes. However, when these duplications occurred remains unclear, partly due to the high relatedness of sequenced species. To better understand the origin of aphid duplications we sequenced and assembled the genome of Cinara cedri, an early branching lineage (Lachninae) of the Aphididae family. We performed a phylogenomic comparison of this genome with 20 other sequenced genomes, including the available genomes of five other aphids, along with the transcriptomes of two species belonging to Adelgidae (a closely related clade to the aphids) and Coccoidea. We found that gene duplication has been pervasive throughout the evolution of aphids, including many parallel waves of recent, species-specific duplications. Most notably, we identified a consistent set of very ancestral duplications, originating from a large-scale gene duplication predating the diversification of Aphidomorpha (comprising aphids, phylloxerids, and adelgids). Genes duplicated in this ancestral wave are enriched in functions related to traits shared by Aphidomorpha, such as association with endosymbionts, and adaptation to plant defenses and phloem-sap-based diet. The ancestral nature of this duplication wave (106-227 Ma) and the lack of sufficiently conserved synteny make it difficult to conclude whether it originated from a whole-genome duplication event or, alternatively, from a burst of large-scale segmental duplications. Genome sequencing of other aphid species belonging to different Aphidomorpha and related lineages may clarify these findings.

Genome Evolution of the Obligate Endosymbiont Buchnera aphidicola

An evolutionary consequence of uniparentally transmitted symbiosis is degradation of symbiont genomes. We use the system of aphids and their maternally inherited obligate endosymbiont, Buchnera aphidicola, to explore the evolutionary process of genome degradation. We compared complete genome sequences for 39 Buchnera strains, including 23 newly sequenced symbiont genomes from diverse aphid hosts. We reconstructed the genome of the most recent shared Buchnera ancestor, which contained 616 protein-coding genes, and 39 RNA genes. The extent of subsequent gene loss varied across lineages, resulting in modern genomes ranging from 412 to 646 kb and containing 354-587 protein-coding genes. Loss events were highly nonrandom across loci. Genes involved in replication, transcription, translation, and amino acid biosynthesis are largely retained, whereas genes underlying ornithine biosynthesis, stress responses, and transcriptional regulation were lost repeatedly. Aside from losses, gene order is almost completely stable. The main exceptions involve movement between plasmid and chromosome locations of genes underlying tryptophan and leucine biosynthesis and supporting nutrition of aphid hosts. This set of complete genomes enabled tests for signatures of positive diversifying selection. Of 371 Buchnera genes tested, 29 genes show strong support for ongoing positive selection. These include genes encoding outer membrane porins that are expected to be involved in direct interactions with hosts. Collectively, these results indicate that extensive genome reduction occurred in the ancestral Buchnera prior to aphid diversification and that reduction has continued since, with losses greater in some lineages and for some loci.

The complete mitochondrial genome of Greenidea ficicola (Hemiptera: Aphididae: Greenideinae), a pest of Ficus

In this study, the complete mitochondrial genome of the pest aphid Greenidea ficicola was determined. The mitogenome was 17,361 bp in length, containing 13 protein-coding genes, 22 tRNA genes, 2 rRNA genes, 1 long control region located between srRNA and tRNA^Ile , and a repeat region located between tRNA^Glu and tRNA^Phe . Thirteen protein-coding genes have typical ATN start codon and TAA termination codon. All tRNAs were predicted to contain typical clover-leaf secondary structures except tRNA^Ser(gct). The length of lrRNA and srRNA are 1270 bp and 773 bp, respectively. Phylogenetic analysis shows that Greenideinae speices form a highly supported clade.

The complete mitochondrial genome of Paracolopha morrisoni (Baker, 1919) (Hemiptera: Aphididae)

We have determined the mitochondrial genome of Paracolopha morrisoni (Baker, 1919), a gall-forming aphid collected from Korea. The circular mitogenome of Paracolopha morrisoni is 16,330 bp long including 13 protein-coding genes, 2 ribosomal RNA genes, 22 transfer RNAs, and a single large non-coding region of 932 bp. The base composition was AT-biased (84.9%). Gene order of P. morrisoni is identical to all other aphid mitochondrial genomes. Phylogenetic trees show that P. morrisoni is sister to Eriosoma lanigerum both belonging to tribe. The mitochondrial genome of P. morrisoni will be useful in understanding the genetic backgrounds of the species.

Congruent phylogenetic relationships of Melaphidina aphids (Aphididae: Eriosomatinae: Fordini) according to nuclear and mitochondrial DNA data with taxonomic implications on generic limits

Melaphidina aphids (Rhus-gall aphids; Eriosomatinae: Fordini) comprise five genera from eastern Asia and one monotypic genus from eastern North America. Melaphidina are unique in feeding on plant species of Rhus subgenus Rhus (Anacardiaceae), on which they form galls during the summer. The phylogenetic relationships among some species of Melaphidina aphids remain controversial. In this study, we sought to resolve the backbone phylogeny of Melaphidina aphids by sampling 15 accessions representing all six genera, all species, and all subspecies except Meitanaphis microgallis using 20 gene regions: five nuclear genes as well as 13 protein-coding genes and two rRNA genes of the mitochondrial genome. Phylogenetic analyses included Bayesian and maximum likelihood methods. Independent analyses of nuclear and mitochondrial genes returned congruent topologies, and analyses of all gene regions combined showed well-supported relationships among Melaphidina species. In particular, these were: (1) Nurudea (excluding N. ibofushi) is sister to a clade composed of the five remaining genera; (2) the monotypic North American genus Melaphis is sister to a clade comprising the four remaining genera; and (3) (Schlechtendalia + N. ibofushi) is sister to the clade (Floraphis (Meitanaphis + Kaburagia). Our results support the transfer of Meitanaphis flavogallis to Kaburagia as an additional subspecies or species, and the recognition of Floraphis as a distinct genus. This study provides important molecular resources for subsequent evolutionary studies using more nuclear genes on the Melaphidina aphids and their close allies.

A large-scale field study of bacterial communities in cereal aphid populations across Morocco

Insects are frequently associated with bacteria that can have significant ecological and evolutionary impacts on their hosts. To date, few studies have examined the influence of environmental factors to microbiome composition of aphids. The current work assessed the diversity of bacterial communities of five cereal aphid species (Sitobion avenae, Rhopalosiphum padi, R. maidis, Sipha maydis and Diuraphis noxia) collected across Morocco, covering a wide range of environmental conditions. We aimed to test whether symbiont combinations are host or environment specific. Deep 16S rRNA sequencing enabled us to identify 17 bacterial operational taxonomic units (OTUs). The obligate symbiont Buchnera aphidicola was represented by five OTUs with multiple haplotypes in many single samples. Facultative endosymbionts were presented by a high prevalence of Regiella insecticola and Serratia symbiotica in S. avenae and Si. maydis, respectively. In addition to these symbiotic partners, Pseudomonas, Acinetobacter, Pantoea, Erwinia and Staphyloccocus were also identified in aphids, suggesting that the aphid microbiome is not limited to the presence of endosymbiotic bacteria. Beside a significant association between host species and bacterial communities, an inverse correlation was also found between altitude and α-diversity. Overall, our results support that symbiont combinations are mainly host specific.

The mitochondrial genome of Greenidea psidii van der Goot (Hemiptera: Aphididae: Greenideinae) and comparisons with other Aphididae aphids

The complete mitochondrial genome of Greenidea psidii was sequenced and compared with the genomes of other aphid species. The G. psidii mitogenome is 16,202 bp long with an A + T content of 85.4%, comprising 37 genes arranged in the same order as the inferred insect ancestral arrangement, a control region and a special repeat region between trnE and trnF. G. psidii is the first representative possessing such a repeat region from the subfamily Greenideinae. These repeat motifs found in Aphidinae, Eriosomatinae and Greenideinae are species-specific, differing in nucleotide composition, length and copy number across different species. All reported complete aphid mitogenomes are A- and C-skewed in the whole genomes and show opposite skewness for C and G between the majority and minority strands. The ratios of nonsynonymous to synonymous substitution rates indicated that the evolution of the aphid mitogenomes has been dominated by purifying selection. Based on whole mitogenome sequences, the phylogenetic relationships among all aphid species with complete mitogenomes were investigated. G. psidii was robustly clustered with Cervaphis quercus, confirming its affiliation with the subfamily Greenideinae.

A Freeloader? The Highly Eroded Yet Large Genome of the Serratia symbiotica Symbiont of Cinara strobi

Genome reduction is pervasive among maternally inherited bacterial endosymbionts. This genome reduction can eventually lead to serious deterioration of essential metabolic pathways, thus rendering an obligate endosymbiont unable to provide essential nutrients to its host. This loss of essential pathways can lead to either symbiont complementation (sharing of the nutrient production with a novel co-obligate symbiont) or symbiont replacement (complete takeover of nutrient production by the novel symbiont). However, the process by which these two evolutionary events happen remains somewhat enigmatic by the lack of examples of intermediate stages of this process. Cinara aphids (Hemiptera: Aphididae) typically harbor two obligate bacterial symbionts: Buchnera and Serratia symbiotica. However, the latter has been replaced by different bacterial taxa in specific lineages, and thus species within this aphid lineage could provide important clues into the process of symbiont replacement. In the present study, using 16S rRNA high-throughput amplicon sequencing, we determined that the aphid Cinara strobi harbors not two, but three fixed bacterial symbionts: Buchnera aphidicola, a Sodalis sp., and S. symbiotica. Through genome assembly and genome-based metabolic inference, we have found that only the first two symbionts (Buchnera and Sodalis) actually contribute to the hosts' supply of essential nutrients while S. symbiotica has become unable to contribute towards this task. We found that S. symbiotica has a rather large and highly eroded genome which codes only for a few proteins and displays extensive pseudogenization. Thus, we propose an ongoing symbiont replacement within C. strobi, in which a once "competent" S. symbiotica does no longer contribute towards the beneficial association. These results suggest that in dual symbiotic systems, when a substitute cosymbiont is available, genome deterioration can precede genome reduction and a symbiont can be maintained despite the apparent lack of benefit to its host.

Is there a relationship between the morphology of the forewing axillary sclerites and the way the wing folds in aphids (Aphidomorpha, Sternorrhyncha, Hemiptera)?

The present study describes the relationship between the morphology of the forewing axillary sclerites and the way the wings fold among 24 aphid genera as compared to a representative of coccids. Architecture of the forewing base was imaged with scanning electron and optical (fluorescence) microscopy. Significant differences in morphology of axillary sclerites between aphid species were observed, despite their belonging to one infraorder. Detailed description of 41 features of axillary sclerites was made. There was no difference between axillaries of viviparous (Aphididae) and oviparous (Adelges sp., Phylloxera sp.) species. No clear relationship between morphology of the axillary sclerites and the wing folding could be confirmed. Instead, the thorax structure determines the way the wing folds in aphids. Phylogenetic analysis based on our results cannot be conducted at this stage of study. To show how three-dimensional the structures are and how difficult to describe, a short animation of Aphis fabae (Aphididae) wing base was added. This is a preliminary study about morphology of axillary sclerites among aphids.

Can social partnerships influence the microbiome? Insights from ant farmers and their trophobiont mutualists

Mutualistic interactions with microbes have played a crucial role in the evolution and ecology of animal hosts. However, it is unclear what factors are most important in influencing particular host–microbe associations. While closely related animal species may have more similar microbiota than distantly related species due to phylogenetic contingencies, social partnerships with other organisms, such as those in which one animal farms another, may also influence an organism's symbiotic microbiome. We studied a mutualistic network of Brachymyrmex and Lasius ants farming several honeydew‐producing Prociphilus aphids and Rhizoecus mealybugs to test whether the mutualistic microbiomes of these interacting insects are primarily correlated with their phylogeny or with their shared social partnerships. Our results confirm a phylogenetic signal in the microbiomes of aphid and mealybug trophobionts, with each species harbouring species‐specific endosymbiont strains of Buchnera (aphids), Tremblaya and Sodalis (mealybugs), and Serratia (both mealybugs and aphids) despite being farmed by the same ants. This is likely explained by strict vertical transmission of trophobiont endosymbionts between generations. In contrast, our results show the ants’ microbiome is possibly shaped by their social partnerships, with ants that farm the same trophobionts also sharing strains of sugar‐processing Acetobacteraceae bacteria, known from other honeydew‐feeding ants and which likely reside extracellularly in the ants’ guts. These ant–microbe associations are arguably more “open” and subject to horizontal transmission or social transmission within ant colonies. These findings suggest that the role of social partnerships in shaping a host's symbiotic microbiome can be variable and is likely dependent on how the microbes are transmitted across generations.

Evolutionary loss and replacement of Buchnera, the obligate endosymbiont of aphids

Symbiotic interactions between organisms create new ecological niches. For example, many insects survive on plant-sap with the aid of maternally transmitted bacterial symbionts that provision essential nutrients lacking in this diet. Symbiotic partners often enter a long-term relationship in which the co-evolutionary fate of lineages is interdependent. Obligate symbionts that are strictly maternally transmitted experience genetic drift and genome degradation, compromising symbiont function and reducing host fitness unless hosts can compensate for these deficits. One evolutionary solution is the acquisition of a novel symbiont with a functionally intact genome. Whereas almost all aphids host the anciently acquired bacterial endosymbiont Buchnera aphidicola (Gammaproteobacteria), Geopemphigus species have lost Buchnera and instead contain a maternally transmitted symbiont closely related to several known insect symbionts from the bacterial phylum Bacteroidetes. A complete genome sequence shows the symbiont has lost many ancestral genes, resulting in a genome size intermediate between that of free-living and symbiotic Bacteroidetes. The Geopemphigus symbiont retains biosynthetic pathways for amino acids and vitamins, as in Buchnera and other insect symbionts. This case of evolutionary replacement of Buchnera provides an opportunity to further understand the evolution and functional genomics of symbiosis.

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.

Transovarial Transmission of Symbionts in Insects

Many insects, on account of their unbalanced diet, live in obligate symbiotic associations with microorganisms (bacteria or yeast-like symbionts), which provide them with substances missing in the food they consume. In the body of host insect, symbiotic microorganisms may occur intracellularly (e.g., in specialized cells of mesodermal origin termed bacteriocytes, in fat body cells, in midgut epithelium) or extracellularly (e.g., in hemolymph, in midgut lumen). As a rule, symbionts are vertically transmitted to the next generation. In most insects, symbiotic microorganisms are transferred from mother to offspring transovarially within female germ cells. The results of numerous ultrastructural and molecular studies on symbiotic systems in different groups of insects have shown that they have a large diversity of symbiotic microorganisms and different strategies of their transmission from one generation to the next. This chapter reviews the modes of transovarial transmission of symbionts between generations in insects.

All 37 Mitochondrial Genes of Aphid Aphis craccivora Obtained from Transcriptome Sequencing: Implications for the Evolution of Aphids

The availability of mitochondrial genome data for Aphididae, one of the economically important insect pest families, in public databases is limited. The advent of next generation sequencing technology provides the potential to generate mitochondrial genome data for many species timely and cost-effectively. In this report, we used transcriptome sequencing technology to determine all the 37 mitochondrial genes of the cowpea aphid, Aphis craccivora. This method avoids the necessity of finding suitable primers for long PCRs or primer-walking amplicons, and is proved to be effective in obtaining the whole set of mitochondrial gene data for insects with difficulty in sequencing mitochondrial genome by PCR-based strategies. Phylogenetic analyses of aphid mitochondrial genome data show clustering based on tribe level, and strongly support the monophyly of the family Aphididae. Within the monophyletic Aphidini, three samples from Aphis grouped together. In another major clade of Aphididae, Pterocomma pilosum was recovered as a potential sister-group of Cavariella salicicola, as part of Macrosiphini.