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

Characteristics and Comparative Analysis of Mitochondrial Genomes of the Aphid Genus Hyalopterus Koch (Hemiptera: Aphididae: Aphidinae)

Using Illumina sequencing technology, we generated complete mitochondrial genomes (mitogenomes) of three constituent species of the aphid genus Hyalopterus Koch, Hyalopterus amygdali (Blanchard), Hyalopterus arundiniformis Ghulamullah, and Hyalopterus pruni (Geoffroy). The sizes of the Hyalopterus mitogenomes range from 15,306 to 15,410 bp, primarily due to variations in the length of non-coding regions. The Hyalopterus mitogenomes consist of 37 coding genes arranged in the order of the ancestral insect mitogenome, a control region, and a repeat region between trnE and trnF. According to the COI-based analysis, one previously reported mitogenome of H. pruni should be assigned to H. arundiniformis. The gene order, nucleotide composition, and codon usage in the Hyalopterus mitogenomes are highly conserved and similar to those of other species of Aphidinae. The tandem repeat units differ in nucleotide composition, length, and copy number across three Hyalopterus species. Within the widespread Eurasian species H. arundiniformis, variation in repeat units among different geographic populations is observed, indicating that the repeat region may provide valuable insights for studying the intraspecific diversification of aphids. Phylogenetic analyses based on 28 complete mitogenomes of Aphidinae supported the monophyly of Aphidinae, Aphidini, Macrosiphini, and two subtribes of Aphidini. Hyalopterus was monophyletic. H. amygdali and H. pruni formed a sister group, while H. arundiniformis was placed basally. Characterization of the mitogenomes of Hyalopterus provides valuable resources for further comparative studies and for advancing our understanding of the aphid mitogenome architecture.

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.

Large-scale phylogenomics of aquatic bacteria reveal molecular mechanisms for adaptation to salinity

The crossing of environmental barriers poses major adaptive challenges. Rareness of freshwater-marine transitions separates the bacterial communities, but how these are related to brackish counterparts remains elusive, as do the molecular adaptations facilitating cross-biome transitions. We conducted large-scale phylogenomic analysis of freshwater, brackish, and marine quality-filtered metagenome-assembled genomes (11,248). Average nucleotide identity analyses showed that bacterial species rarely existed in multiple biomes. In contrast, distinct brackish basins cohosted numerous species, but their intraspecific population structures displayed clear signs of geographic separation. We further identified the most recent cross-biome transitions, which were rare, ancient, and most commonly directed toward the brackish biome. Transitions were accompanied by systematic changes in amino acid composition and isoelectric point distributions of inferred proteomes, which evolved over millions of years, as well as convergent gains or losses of specific gene functions. Therefore, adaptive challenges entailing proteome reorganization and specific changes in gene content constrains the cross-biome transitions, resulting in species-level separation between aquatic biomes.

The phylogeny and evolutionary ecology of hoverflies (Diptera: Syrphidae) inferred from mitochondrial genomes

Hoverflies (Diptera: Syrphidae) are a diverse group of pollinators and a major research focus in ecology, but their phylogenetic relationships remain incompletely known. Using a genome skimming approach we generated mitochondrial genomes for 91 species, capturing a wide taxonomic diversity of the family. To reduce the required amount of input DNA and overall cost of the library construction, sequencing and assembly was conducted on mixtures of specimens, which raises the problem of chimera formation of mitogenomes. We present a novel chimera detection test based on gene tree incongruence, but identified only a single mitogenome of chimeric origin. Together with existing data for a final set of 127 taxa, phylogenetic analysis on nucleotide and amino acid sequences using Maximum Likelihood and Bayesian Inference revealed a basal split of Microdontinae from all other syrphids. The remainder consists of several deep clades assigned to the subfamily Eristalinae in the current classification, including a clade comprising the subfamily Syrphinae (plus Pipizinae). These findings call for a re-definition of subfamilies, but basal nodes had insufficient support to allow such action. Molecular-clock dating placed the origin of the Syrphidae crown group in the mid-Cretaceous while the Eristalinae-Syrphinae clade likely originated near the K/Pg boundary. Transformation of larval life history characters on the tree suggests that Syrphidae initially had sap feeding larvae, which diversified greatly in diet and habitat association during the Eocene and Oligocene, coinciding with the diversification of angiosperms and the evolution of various insect groups used as larval host, prey, or mimicry models. Mitogenomes proved to be a powerful phylogenetic marker for studies of Syrphidae at subfamily and tribe levels, allowing dense taxon sampling that provided insight into the great ecological diversity and rapid evolution of larval life history traits of the hoverflies.

Analysis of the Genetic Diversity of Two Rhopalosiphum Species from China and Europe Based on Nuclear and Mitochondrial Genes

Population genetic studies can reveal clues about the evolution of adaptive strategies of aphid species in agroecosystems and demonstrate the influence of environmental factors on the genetic diversity and gene flow among aphid populations. To investigate the genetic diversity of two Rhopalosiphum aphid species from different geographical regions, 32 populations (n = 535) of the bird cherry-oat aphid (Rhopalosiphum padi Linnaeus) and 38 populations (n = 808) of the corn leaf aphid (Rhopalosiphum maidis Fitch) from China and Europe were analyzed using one nuclear (elongation factor-1 alpha) and two mitochondrial (cytochrome oxidase I and II) genes. Based on the COI-COII sequencing, two obvious clades between Chinese and European populations and a low level of gene flow (Nm = 0.15) were detected in R. padi, while no geographical-associated genetic variation was found for EF-1α in this species. All genes in R. maidis had low genetic variation, indicating a high level of gene flow (Nm = 5.31 of COI-COII and Nm = 2.89 of EF-1α). Based on the mitochondrial result of R. padi, we concluded that the long distance between China and Europe may be interrupting the gene flow. The discordant results of nuclear gene analyses in R. padi may be due to the slower evolution of nuclear genes compared to mitochondrial genes. The gene exchange may occur gradually with the potential for continuous migration of the aphid. This study facilitates the design of control strategies for these pests.

European and Mediterranean Myzocallidini Aphid Species: DNA Barcoding and Remarks on Ecology with Taxonomic Modifications in An Integrated Framework

The genus Myzocallis Passerini (Hemiptera, Aphididae, Calaphidinae, Myzocallidini) is a rather primitive group of aphids currently comprising 45 species and 3 subspecies, subdivided into ten subgenera, three of them having a West Palaearctic distribution. The majority of the species inhabit Fagales plants and some of them are considered pests. Despite their ecological interest and the presence of some taxonomic controversies, there are only a few molecular studies on the group. Here, the main aims were to develop a DNA barcodes library for the molecular identification of West Palaearctic Myzocallis species, to evaluate the congruence among their morphological, ecological and DNA-based delimitation, and verify the congruence of the subgeneric subdivision presently adopted by comparing the results with those obtained for other Panaphidini species. These study findings indicate that Myzocallis (Agrioaphis) leclanti, originally described as a subspecies of M. (A.) castanicola and M. (M.) schreiberi, considered as a subspecies of M. (M.) boerneri, should be regarded at a rank of full species, and the subgenera Agrioaphis, Lineomyzocallis, Neomyzocallis, Pasekia were elevated to the rank of genus, while Myzocallis remain as such.

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.

Repeated dietary shifts in elapid snakes (Squamata: Elapidae) revealed by ancestral state reconstruction

Identifying the traits of ancestral organisms can reveal patterns and drivers of organismal diversification. Unfortunately, reconstructing complex multistate traits (such as diet) remains challenging. Adopting a ‘reconstruct, then aggregate’ approach in a maximum likelihood framework, we reconstructed ancestral diets for 298 species of elapid snakes. We tested whether different prey types were correlated with one another, tested for one-way contingency between prey type pairs, and examined the relationship between snake body size and dietary composition. We demonstrate that the evolution of diet was characterized by niche conservation punctuated by repeated dietary shifts. The ancestor of elapids most likely fed on reptiles and possibly amphibians, with deviations from this ancestral diet occurring repeatedly due to shifts into marine environments and changes in body size. Moreover, we demonstrate important patterns of prey use, including one-way dependency—most obviously the inclusion of eggs being dependent on a diet that already included the producers of those eggs. Despite imperfect dietary data, our approach produced a robust overview of dietary evolution. Given the paucity of natural history information for many organisms, our approach has the potential to increase the number of lineages to which ancestral state reconstructions of multistate traits can be robustly applied.

Characterization of Mariner transposons in seven species of Rhus gall aphids

Transposable elements (TEs), also known as jumping genes, are widely spread in the genomes of insects and play a considerable role in genomic evolution. Mariner/DD34D family belongs to class II transposable elements which is widely spread in the genomes of insects and have considerable role in genomic evolution. Mariner like elements (MLEs) were searched in the genomes of seven species of Rhus gall aphids belonging to six genera. In total, 121 MLEs were detected in the genomes of the seven investigated species of Rhus gall aphids, which showed a wide distribution in both close and distant related species. The sequences of MLEs ranged from 1 to 1.4 kb in length and the structural analysis of the MLEs showed that only five copies were potentially active with intact open reading frame (ORF) and terminal inverted repeats (TIRs). Phylogenetic analysis showed that all the 121 MLE sequences belonged to four subfamilies, i.e., Mauritiana, Drosophila, Vertumana and Irritans, among which Drosophila and Vertumana subfamilies were reported in aphids for the first time. Our present report revealed the diversity and distribution of MLEs in Rhus gall aphid genomes and expanded our understandings on the characterization of transposable elements in aphid genomes, which might be useful as genetic markers and tools and would play an important role in genomic evolution and adaptation of aphids.

A new species of aphid of the genus Nipponaphis (Hemiptera: Aphididae: Hormaphidinae) from China, inducing galls on the trunk of a witch-hazel (Hamamelidaceae)

Plants and insects have co-existed for millions of years. Although research has been conducted on various insect species that induce galls on various plant tissues, information is particularly scarce when it comes to insects that form galls on the tough trunk of their host plants. This contribution describes the gall-inducing aphid Nipponaphis hubeiensis sp. nov. from the Zhushan County, Shiyan City, Hubei Province of China. This aphid induces enclosed galls with woody external layer on the trunk of Sycopsis sinensis (Saxifragales: Hamamelidaceae), an uncommon ecological niche in the aphid-plant interaction system. Morphological features for the identification of new species are provided. In addition, a partial sequence of the nuclear gene EF1α was amplified and sequenced to construct a cluster graph. Based on the clustering graph combined with morphology traits, the gall-forming aphid was classified into Nipponaphis. The unique ecological habits of this new aphid will bring innovative perspectives to the study of the evolution and diversity in aphid-host interaction.

Complex evolution in Aphis gossypii group (Hemiptera: Aphididae), evidence of primary host shift and hybridization between sympatric species

Aphids provide a good model system to understand the ecological speciation concept, since the majority of the species are host-specific, and they spend their entire lifecycle on certain groups of host plants. Aphid species that apparently have wide host plant ranges have often turned out to be complexes of host-specialized biotypes. Here we investigated the various host-associated populations of the two recently diverged species, Aphis gossypii and A. rhamnicola, having multiple primary hosts, to understand the complex evolution with host-associated speciation. Using mitochondrial DNA marker and nine microsatellite loci, we reconstructed the haplotype network, and analyzed the genetic structure and relationships. Approximate Bayesian computation was also used to infer the ancestral primary host and host-associated divergence, which resulted in Rhamnus being the most ancestral host for A. gossypii and A. rhamnicola. As a result, Aphis gossypii and A. rhamnicola do not randomly use their primary and secondary host plants; rather, certain biotypes use only some secondary and specific primary hosts. Some biotypes are possibly in a diverging state through specialization to specific primary hosts. Our results also indicate that a new heteroecious race can commonly be derived from the heteroecious ancestor, showing strong evidence of ecological specialization through a primary host shift in both A. gossypii and A. rhamnicola. Interestingly, A. gossypii and A. rhamnicola shared COI haplotypes with each other, thus there is a possibility of introgression by hybridization between them by cross-sharing same primary hosts. Our results contribute to a new perspective in the study of aphid evolution by identifying complex evolutionary trends in the gossypii sensu lato complex.

Genome Sequence of the Banana Aphid, Pentalonia nigronervosa Coquerel (Hemiptera: Aphididae) and Its Symbionts

The banana aphid, Pentalonia nigronervosa Coquerel (Hemiptera: Aphididae), is a major pest of cultivated bananas (Musa spp., order Zingiberales), primarily due to its role as a vector of Banana bunchy top virus (BBTV), the most severe viral disease of banana worldwide. Here, we generated a highly complete genome assembly of P. nigronervosa using a single PCR-free Illumina sequencing library. Using the same sequence data, we also generated complete genome assemblies of the P. nigronervosa symbiotic bacteria Buchnera aphidicola and Wolbachia. To improve our initial assembly of P. nigronervosa we developed a k-mer based deduplication pipeline to remove genomic scaffolds derived from the assembly of haplotigs (allelic variants assembled as separate scaffolds). To demonstrate the usefulness of this pipeline, we applied it to the recently generated assembly of the aphid Myzus cerasi, reducing the duplication of conserved BUSCO genes by 25%. Phylogenomic analysis of P. nigronervosa, our improved M. cerasi assembly, and seven previously published aphid genomes, spanning three aphid tribes and two subfamilies, reveals that P. nigronervosa falls within the tribe Macrosiphini, but is an outgroup to other Macrosiphini sequenced so far. As such, the genomic resources reported here will be useful for understanding both the evolution of Macrosphini and for the study of P. nigronervosa. Furthermore, our approach using low cost, high-quality, Illumina short-reads to generate complete genome assemblies of understudied aphid species will help to fill in genomic black spots in the diverse aphid tree of life.

New data on the aphid (Hemiptera, Aphididae) fauna of New Caledonia: some new biosecurity threats in a biodiversity hotspot

Thirty-three species of aphids are now established in New Caledonia. All species appear to have been introduced accidentally by human activity in the last century. Here, 17 aphid species are recorded for the first time: Aphiseugeniae, Aphisglycines, Aphisodinae, Aulacorthumsolani, Brachycaudushelichrysi, Cerataphisorchidearum, Greenideapsidii, Hyperomyzuscarduellinus, Hysteroneurasetariae, Lipaphispseudobrassicae, Micromyzuskatoi, Myzusornatus, Pentaloniacaladii, Rhopalosiphumnymphaeae, Rhopalosiphumrufiabdominale, Schizaphisrotundiventris, and Tetraneurafusiformis. Thirteen more species are also more or less regularly intercepted at the borders through biosecurity surveys, without further establishment. This demonstrates that aphids represent a major biosecurity threat, including a threat as potential plant virus vectors. The reinforcement of biosecurity is a priority for such biodiversity hotspots, from the perspectives of both agriculture and the native environment. Prioritisation and promotion of local development of vegetable and fruit production, rather than their risky importation from abroad, is desirable. Such an approach also should be promoted and extended to other Pacific islands, which all share the lack of native aphid fauna and their associated plant disease vector risks.

A large genomic insertion containing a duplicated follistatin gene is linked to the pea aphid male wing dimorphism

Wing dimorphisms have long served as models for examining the ecological and evolutionary tradeoffs associated with alternative phenotypes. Here, we investigated the genetic cause of the pea aphid (Acyrthosiphon pisum) male wing dimorphism, wherein males exhibit one of two morphologies that differ in correlated traits that include the presence or absence of wings. We mapped this trait difference to a single genomic region and, using third generation, long-read sequencing, we identified a 120 kb insertion in the wingless allele. This insertion includes a duplicated follistatin gene, which is a strong candidate gene in the minimal mapped interval to cause the dimorphism. We found that both alleles were present prior to pea aphid biotype lineage diversification, we estimated that the insertion occurred millions of years ago, and we propose that both alleles have been maintained in the species, likely due to balancing selection.

The Potyviruses: An Evolutionary Synthesis Is Emerging

In this review, encouraged by the dictum of Theodosius Dobzhansky that "Nothing in biology makes sense except in the light of evolution", we outline the likely evolutionary pathways that have resulted in the observed similarities and differences of the extant molecules, biology, distribution, etc. of the potyvirids and, especially, its largest genus, the potyviruses. The potyvirids are a family of plant-infecting RNA-genome viruses. They had a single polyphyletic origin, and all share at least three of their genes (i.e., the helicase region of their CI protein, the RdRp region of their NIb protein and their coat protein) with other viruses which are otherwise unrelated. Potyvirids fall into 11 genera of which the potyviruses, the largest, include more than 150 distinct viruses found worldwide. The first potyvirus probably originated 15,000-30,000 years ago, in a Eurasian grass host, by acquiring crucial changes to its coat protein and HC-Pro protein, which enabled it to be transmitted by migrating host-seeking aphids. All potyviruses are aphid-borne and, in nature, infect discreet sets of monocotyledonous or eudicotyledonous angiosperms. All potyvirus genomes are under negative selection; the HC-Pro, CP, Nia, and NIb genes are most strongly selected, and the PIPO gene least, but there are overriding virus specific differences; for example, all turnip mosaic virus genes are more strongly conserved than those of potato virus Y. Estimates of dN/dS (ω) indicate whether potyvirus populations have been evolving as one or more subpopulations and could be used to help define species boundaries. Recombinants are common in many potyvirus populations (20%-64% in five examined), but recombination seems to be an uncommon speciation mechanism as, of 149 distinct potyviruses, only two were clear recombinants. Human activities, especially trade and farming, have fostered and spread both potyviruses and their aphid vectors throughout the world, especially over the past five centuries. The world distribution of potyviruses, especially those found on islands, indicates that potyviruses may be more frequently or effectively transmitted by seed than experimental tests suggest. Only two meta-genomic potyviruses have been recorded from animal samples, and both are probably contaminants.

Horizontal Transfer of Bacterial Cytolethal Distending Toxin B Genes to Insects

Horizontal gene transfer events have played a major role in the evolution of microbial species, but their importance in animals is less clear. Here, we report horizontal gene transfer of cytolethal distending toxin B (cdtB), prokaryotic genes encoding eukaryote-targeting DNase I toxins, into the genomes of vinegar flies (Diptera: Drosophilidae) and aphids (Hemiptera: Aphididae). We found insect-encoded cdtB genes are most closely related to orthologs from bacteriophage that infect Candidatus Hamiltonella defensa, a bacterial mutualistic symbiont of aphids that confers resistance to parasitoid wasps. In drosophilids, cdtB orthologs are highly expressed during the parasitoid-prone larval stage and encode a protein with ancestral DNase activity. We show that cdtB has been domesticated by diverse insects and hypothesize that it functions in defense against their natural enemies.

Description of soybean aphid (Aphis glycines Matsumura) mitochondrial genome and comparative mitogenomics of Aphididae (Hemiptera: Sternorrhyncha)

The complete mitochondrial genome of the soybean aphid (Aphis glycines Matsumura), a major agricultural pest in the world, is described for the first time, which consists of 13 protein-coding genes, 22 tRNA genes, 2 rRNA genes, as well as a large repeat region between tRNA-Glu and tRNA-Phe, and an AT-rich control region. The 17,954 bp mtgenome is the largest reported from the family Aphididae, and its gene order follows the ancestral insect mtgenome except for the repeat region, which contains a 195 bp unit repeated 11.9 times, representing the highest reported repeats among the known aphid mtgenomes to date. A new molecular phylogeny of Aphidae is reconstructed based on all available aphid mtgenomes, and it is shown that the mtgenome data can robustly resolve relationships at the subfamily level, but do not have sufficient phylogenetic information to resolve deep relationships. A phylogeny-based comparative analysis of mtgenomes has been performed to investigate the evolution of the repeat region between tRNA-Glu and tRNA-Phe. So far, among aphids, 13 species are known to have this repeat region of variable lengths, and a phylogenetic analysis of the repeat region shows that a large proportion of the sequences are conserved across the phylogeny, suggesting that the repeat region evolved in the most recent common ancestor of Aphidinae and Eriosomatinae, and that it has gone through numerous episodes of lineage-specific losses and expansions. Combined together, this study provides novel insights into how the repeat regions have evolved within aphids.

Cross-Kingdom Analysis of Diversity, Evolutionary History, and Site Selection within the Eukaryotic Macrophage Migration Inhibitory Factor Superfamily

Macrophage migration inhibitory factors (MIF) are multifunctional proteins regulating major processes in mammals, including activation of innate immune responses. MIF proteins also play a role in innate immunity of invertebrate organisms or serve as virulence factors in parasitic organisms, raising the question of their evolutionary history. We performed a broad survey of MIF presence or absence and evolutionary relationships across 803 species of plants, fungi, protists, and animals, and explored a potential relation with the taxonomic status, the ecology, and the lifestyle of individual species. We show that MIF evolutionary history in eukaryotes is complex, involving probable ancestral duplications, multiple gene losses and recent clade-specific re-duplications. Intriguingly, MIFs seem to be essential and highly conserved with many sites under purifying selection in some kingdoms (e.g., plants), while in other kingdoms they appear more dispensable (e.g., in fungi) or present in several diverged variants (e.g., insects, nematodes), suggesting potential neofunctionalizations within the protein superfamily.

A Laterally Transferred Viral Gene Modifies Aphid Wing Plasticity

Organisms often respond to changing environments by altering development of particular traits. These plastic traits exhibit genetic variation; i.e., genotypes respond differently to the same environmental cues. Theoretical studies have demonstrated the importance of this variation, which is targeted by natural selection, in adapting plastic responses to maximize fitness [1, 2]. However, little is known about the underlying genetic mechanisms. We identify two laterally transferred genes that contribute to variation in a classic example of phenotypic plasticity: the pea aphid's ability to produce winged offspring in response to crowding. We discovered that aphid genotypes vary extensively for this trait and that aphid genes of viral origin are upregulated in response to crowding solely in highly inducible genotypes. We knocked down expression of these genes to demonstrate their functional role in wing plasticity. Through phylogenetic analysis, we found that these genes likely originated from a virus that infects rosy apple aphids and causes their hosts to produce winged offspring [3]. The function of these genes has therefore been retained following transfer to pea aphids. Our results uncover a novel role for co-opted viral genes, demonstrating that they are used to modulate ecologically relevant, plastic phenotypes. Our findings also address a critical question about the evolution of environmentally sensitive traits: whether the genes that control the expression of plastic traits also underlie variation in plasticity. The genes we identify originated from outside aphids themselves, and thus, our work shows that genes formerly unrelated to plasticity can fine-tune the strength of plastic responses to the environment.

microRNA regulation in an ancient obligate endosymbiosis

Although many insects are associated with obligate bacterial endosymbionts, the mechanisms by which these host/endosymbiont associations are regulated remain mysterious. While microRNAs (miRNAs) have been recently identified as regulators of host/microbe interactions, including host/pathogen and host/facultative endosymbiont interactions, the role miRNAs may play in mediating host/obligate endosymbiont interactions is virtually unknown. Here, we identified conserved miRNAs that potentially mediate symbiotic interactions between aphids and their obligate endosymbiont, Buchnera aphidicola. Using small RNA sequence data from Myzus persicae and Acyrthosiphon pisum, we annotated 93 M. persicae and 89 A. pisum miRNAs, among which 69 were shared. We found 14 miRNAs that were either highly expressed in aphid bacteriome, the Buchnera-housing tissue, or differentially expressed in bacteriome vs. gut, a non-Buchnera-housing tissue. Strikingly, 10 of these 14 miRNAs have been implicated previously in other host/microbe interaction studies. Investigating the interaction networks of these miRNAs using a custom computational pipeline, we identified 103 miRNA::mRNA interactions shared between M. persicae and A. pisum. Functional annotation of the shared mRNA targets revealed only two over-represented cluster of orthologous group categories: amino acid transport and metabolism, and signal transduction mechanisms. Our work supports a role for miRNAs in mediating host/symbiont interactions between aphids and their obligate endosymbiont Buchnera. In addition, our results highlight the probable importance of signal transduction mechanisms to host/endosymbiont coevolution.

Diversity and evolution of mariner-like elements in aphid genomes

Background

Although transposons have been identified in almost all organisms, genome-wide information on mariner elements in Aphididae remains unknown. Genomes of Acyrthosiphon pisum, Diuraphis noxia and Myzus persicae belonging to the Macrosiphini tribe, actually available in databases, have been investigated.

Results

A total of 22 lineages were identified. Classification and phylogenetic analysis indicated that they were subdivided into three monophyletic groups, each of them containing at least one putative complete sequence, and several non-autonomous sublineages corresponding to Miniature Inverted-Repeat Transposable Elements (MITE), probably generated by internal deletions. A high proportion of truncated and dead copies was also detected. The three clusters can be defined from their catalytic site: (i) mariner DD34D, including three subgroups of the irritans subfamily (Macrosiphinimar, Batmar-like elements and Dnomar-like elements); (ii) rosa DD41D, found in A. pisum and D. noxia; (iii) a new clade which differs from rosa through long TIRs and thus designated LTIR-like elements. Based on its catalytic domain, this new clade is subdivided into DD40D and DD41D subgroups. Compared to other Tc1/mariner superfamily sequences, rosa DD41D and LTIR DD40-41D seem more related to maT DD37D family.

Conclusion

Overall, our results reveal three clades belonging to the irritans subfamily, rosa and new LTIR-like elements. Data on structure and specific distribution of these transposable elements in the Macrosiphini tribe contribute to the understanding of their evolutionary history and to that of their hosts.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3856-6) contains supplementary material, which is available to authorized users.

Vat, an Amazing Gene Conferring Resistance to Aphids and Viruses They Carry: From Molecular Structure to Field Effects

We review half a century of research on Cucumis melo resistance to Aphis gossypii from molecular to field levels. The Vat gene is unique in conferring resistance to both A. gossypii and the viruses it transmits. This double phenotype is aphid clone-dependent and has been observed in 25 melon accessions, mostly from Asia. It is controlled by a cluster of genes including CC-NLR, which has been characterized in detail. Copy-number polymorphisms (for the whole gene and for a domain that stands out in the LLR region) and single-nucleotide polymorphisms have been identified in the Vat cluster. The role of these polymorphisms in plant/aphid interactions remains unclear. The Vat gene structure suggests a functioning with separate recognition and response phases. During the recognition phase, the VAT protein is thought to interact (likely indirectly) with an aphid effector introduced during cell puncture by the aphid. A few hours later, several miRNAs are upregulated in Vat plants. Peroxidase activity increases, and callose and lignin are deposited in the walls of the cells adjacent to the stylet path, disturbing aphid behavior. In aphids feeding on Vat plants, Piwi-interacting RNA-like sequences are abundant and the levels of other miRNAs are modified. At the plant level, resistance to aphids is quantitative (aphids escape the plant and display low rates of reproduction). Resistance to viruses is qualitative and local. Durability of NLR genes is highly variable. A. gossypii clones are adapted to Vat resistance, either by introducing a new effector that interferes with the deployment of plant defenses, or by adapting to the defenses it triggered. Viruses transmitted in a non-persistent manner cannot adapt to Vat resistance. At population level, Vat reduces aphid density and genetic diversity. The durability of Vat resistance to A. gossypii populations depends strongly on the agro-ecosystem, including, in particular, the presence of other cucurbit crops serving as alternative hosts for adapted clones in fall and winter. At the crop level, Vat resistance decreases the intensity of virus epidemics when A. gossypii is the main aphid vector in the crop environment.

The Complete Mitochondrial Genome of Mindarus keteleerifoliae (Insecta: Hemiptera: Aphididae) and Comparison with Other Aphididae Insects

The mitogenome of Mindarus keteleerifoliae Zhang (Hemiptera: Aphididae) is a 15,199 bp circular molecule. The gene order and orientation of M. keteleerifoliae is similarly arranged to that of the ancestral insect of other aphid mitogenomes, and, a tRNA isomerism event maybe identified in the mitogenome of M. keteleerifoliae. The tRNA-Trp gene is coded in the J-strand and the same sequence in the N-strand codes for the tRNA-Ser gene. A similar phenomenon was also found in the mitogenome of Eriosoma lanigerum. However, whether tRNA isomers in aphids exist requires further study. Phylogenetic analyses, using all available protein-coding genes, support Mindarinae as the basal position of Aphididae. Two tribes of Aphidinae were recovered with high statistical significance. Characteristics of the M. keteleerifoliae mitogenome revealed distinct mitogenome structures and provided abundant phylogenetic signals, thus advancing our understanding of insect mitogenomic architecture and evolution. But, because only eight complete aphid mitogenomes, including M. keteleerifoliae, were published, future studies with larger taxon sampling sizes are necessary.

Amino acid transporter expansions associated with the evolution of obligate endosymbiosis in sap-feeding insects (Hemiptera: sternorrhyncha)

Background

Mutualistic obligate endosymbioses shape the evolution of endosymbiont genomes, but their impact on host genomes remains unclear. Insects of the sub-order Sternorrhyncha (Hemiptera) depend on bacterial endosymbionts for essential amino acids present at low abundances in their phloem-based diet. This obligate dependency has been proposed to explain why multiple amino acid transporter genes are maintained in the genomes of the insect hosts. We implemented phylogenetic comparative methods to test whether amino acid transporters have proliferated in sternorrhynchan genomes at rates grater than expected by chance.

Results

By applying a series of methods to reconcile gene and species trees, inferring the size of gene families in ancestral lineages, and simulating the null process of birth and death in multi-gene families, we uncovered a 10-fold increase in duplication rate in the AAAP family of amino acid transporters within Sternorrhyncha. This gene family expansion was unmatched in other closely related clades lacking endosymbionts that provide essential amino acids.

Conclusions

Our findings support the influence of obligate endosymbioses on host genome evolution by both inferring significant expansions of gene families involved in symbiotic interactions, and discovering increases in the rate of duplication associated with multiple emergences of obligate symbiosis in Sternorrhyncha.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-015-0315-3) contains supplementary material, which is available to authorized users.

Characterization of mariner-like transposons of the mauritiana Subfamily in seven tree aphid species

Mariner-like elements (MLEs) are Class II transposons present in all eukaryotic genomes in which MLEs have been searched for. This article reports the detection of MLEs in seven of the main fruit tree aphid species out of eight species studied. Deleted MLE sequences of 916-919 bp were characterized, using the terminal-inverted repeats (TIRs) of mariner elements belonging to the mauritiana Subfamily as primers. All the sequences detected were deleted copies of full-length elements that included the 3'- and 5'-TIRs but displayed internal deletions affecting Mos1 activity. Networks based on the mtDNA cytochrome oxidase subunit-I (CO-I) and MLE sequences were incongruent, suggesting that mutations in transposon sequences had accumulated before speciation of tree aphid species occurred, and that they have been maintained in this species via vertical transmissions. This is the first evidence of the widespread occurrence of MLEs in aphids.

Molecular and morphological differentiation between Aphis gossypii Glover (Hemiptera, Aphididae) and related species, with particular reference to the North American Midwest

The cotton aphid, Aphis gossypii, is one of the most biologically diverse species of aphids; a polyphagous species in a family where most are host specialists. It is economically important and belongs to a group of closely related species that has challenged aphid taxonomy. The research presented here seeks to clarify the taxonomic relationships and status of species within the Aphid gossypii group in the North American Midwest. Sequences of the mitochondrial cytochrome oxidase 1 (COI), nuclear elongation factor 1-α (EF1-α), and nuclear sodium channel para-type (SCP) genes were used to differentiate between Aphid gossypii and related species. Aphis monardae, previously synonymised with Aphid gossypii, is re-established as a valid species. Phylogenetic analyses support the close relationship of members of the Aphid gossypii group native to North America (Aphid forbesi, Aphid monardae, Aphid oestlundi, Aphid rubifolii, and Aphid rubicola), Europe (Aphid nasturtii, Aphid urticata and Aphid sedi), and Asia (Aphid agrimoniae, Aphid clerodendri, Aphid glycines, Aphid gossypii, Aphid hypericiphaga, Aphid ichigicola, Aphid ichigo, Aphid sanguisorbicola, Aphid sumire and Aphid taraxicicola). The North American species most closely related to Aphid gossypii are Aphid monardae and Aphid oestlundi. The cosmopolitan Aphid gossypii and Aphid sedi identified in the USA are genetically very similar using COI and EF1-α sequences, but the SCP gene shows greater genetic distance between them. We present a discussion of the biological and morphological differentiation of these species.

Molecular characterization of the Aphis gossypii olfactory receptor gene families

The cotton aphid, Aphis gossypii Glover, is a polyphagous pest that inflicts great damage to cotton yields worldwide. Antennal olfaction, which is extremely important for insect survival, mediates key behaviors such as host preference, mate choice, and oviposition site selection. In insects, odor detection is mediated by odorant receptors (ORs) and ionotropic receptors (IRs), which ensure the specificity of the olfactory sensory neuron responses. In this study, our aim is to identify chemosensory receptors in the cotton aphid genome, as a means to uncover olfactory encoding of the polyphagous feeding habits as well as to aid the discovery of new targets for behavioral interference. We identified a total of 45 candidate ORs and 14 IRs in the cotton aphid genome. Among the candidate AgoORs, 9 are apparent pseudogenes, while 19 can be clustered with ORs from the pea aphid, forming 16 AgoOR/ApOR orthologous subgroups. Among the candidate IRs, we identified homologs of the two highly conserved co-receptors IR8a and IR25a; no AgoIR retain the complete glutamic acid binding domain, suggesting that putative AgoIRs bind different ligands. Our results provide the necessary information for functional characterization of the chemosensory receptors of A. gossypii, with potential for new or refined applications of semiochemicals-based control of this pest insect.

Aphis (Hemiptera: Aphididae) species groups found in the Midwestern United States and their contribution to the phylogenetic knowledge of the genus

A phylogeny of the genus Aphis Linnaeus, 1 758 was built primarily from specimens collected in the Midwest of the United States. A data matrix was constructed with 68 species and 41 morphological characters with respective character states of alate and apterous viviparous females. Dendrogram topologies of analyses performed using UPGMA (Unweighted Pair Group Method with Arithmetic Mean), Maximum Parsimony and Bayesian analysis of Cytochrome Oxidase I, Elongation Factor 1-α and primary endosymbiont Buchnera aphidicola 16S sequences were not congruent. Bayesian analysis strongly supported most terminal nodes of the phylogenetic trees. The phylogeny was strongly supported by EF1-α, and analysis of COI and EF1-α molecular data combined with morphological characters. It was not supported by single analysis of COI or Buchnera aphidicola 16S. Results from the Bayesian phylogeny show 4 main species groups: asclepiadis, fabae, gossypii, and middletonii. Results place Aphis and species of the genera Protaphis Börner, 1952, Toxoptera Koch, and Xerobion Nevsky, 1928 in a monophyletic clade. Morphological characters support this monophyly as well. The phylogeny shows that the monophyletic clade of the North American middletonii species group belong to the genus Protaphis: P. debilicornis (Gillette & Palmer, 1929), comb. nov., P. echinaceae (Lagos and Voegtlin, 2009), comb. nov., and P. middletonii (Thomas, 1879). The genus Toxoptera should be considered a subgenus of Aphis (stat. nov.). The analysis also indicates that the current genus Iowana Frison, 1954 should be considered a subgenus of Aphis (stat. nov.).

Dynamic recruitment of amino acid transporters to the insect/symbiont interface

Symbiosis is well known to influence bacterial symbiont genome evolution and has recently been shown to shape eukaryotic host genomes. Intriguing patterns of host genome evolution, including remarkable numbers of gene duplications, have been observed in the pea aphid, a sap-feeding insect that relies on a bacterial endosymbiont for amino acid provisioning. Previously, we proposed that gene duplication has been important for the evolution of symbiosis based on aphid-specific gene duplication in amino acid transporters (AATs), with some paralogs highly expressed in the cells housing symbionts (bacteriocytes). Here, we use a comparative approach to test the role of gene duplication in enabling recruitment of AATs to bacteriocytes. Using genomic and transcriptomic data, we annotate AATs from sap-feeding and non sap-feeding insects and find that, like aphids, AAT gene families have undergone independent large-scale gene duplications in three of four additional sap-feeding insects. RNA-seq differential expression data indicate that, like aphids, the sap-feeding citrus mealybug possesses several lineage-specific bacteriocyte-enriched paralogs. Further, differential expression data combined with quantitative PCR support independent evolution of bacteriocyte enrichment in sap-feeding insect AATs. Although these data indicate that gene duplication is not necessary to initiate host/symbiont amino acid exchange, they support a role for gene duplication in enabling AATs to mediate novel host/symbiont interactions broadly in the sap-feeding suborder Sternorrhyncha. In combination with recent studies on other symbiotic systems, gene duplication is emerging as a general pattern in host genome evolution.

Cannibalism in invasive, native and biocontrol populations of the harlequin ladybird

Background

Cannibalism is widespread in both vertebrates and invertebrates but its extent is variable between and within species. Cannibalism depends on population density and nutritional conditions, and could be beneficial during colonisation of new environments. Empirical studies are needed to determine whether this trait might facilitate invasion of a new area in natural systems. We investigated whether the propensity for cannibalism in H. axyridis differs both between native and invasive populations and between invasive populations from the core and from the front of the invasive area in Western Europe. We also compared the propensity for cannibalism of these natural populations with that of laboratory-reared biocontrol populations. We measured the cannibalism rates of eggs by first instar larvae and adult females at two different individual densities of ladybirds from three types of population (invasive, native and biocontrol), in laboratory-controlled conditions.

Results

Cannibalism was significantly greater in larvae from invasive populations compared to native or biocontrol populations, but there was no difference in cannibalism rates between populations from the core or front of the invaded range. Cannibalism was significantly lower in larvae from biocontrol populations compared to wild (invasive and native) populations. No differences in cannibalism rates of adult females were found between any populations. While high population density significantly increased cannibalism in both larvae and adults, the norm of reaction of cannibalism to individual density did not change significantly during the invasion and/or laboratory rearing processes.

Conclusion

This study is the first to provide evidence for a higher propensity for cannibalism in invasive populations compared to native ones. Our experiments also shed light on the difference in cannibalism evolution with respect to life stages. However, we are still at an early stage in understanding the underlying mechanisms and several different research perspectives are needed to determine whether the higher propensity for cannibalism is a general feature of the invasion process.

A total-evidence phylogenetic analysis of Hormaphidinae (Hemiptera: Aphididae), with comments on the evolution of galls

A phylogenetic analysis of Hormaphidinae is presented based on a total-evidence approach. Four genes (two mitochondrial, COI and CytB, and two nuclear, EF-1α and LWO) are combined with 65 morphological and seven biological characters. Sixty-three hormaphidine species representing three tribes and 36 genera as well as nine outgroups are included. Parsimony and model-based approaches are used, and several support values and implied weighting schemes are explored to assess clade stability. The monophyly of Hormaphidinae and Nipponaphidini is supported, but Cerataphidini and Hormaphidini are not recovered as monophyletic. Based on the parsimony hypothesis from the total-evidence analysis, the phylogenetic relationships within Hormaphidinae are discussed. Cerataphidini is re-delimited to exclude Doraphis and Tsugaphis, and Hormaphidini is redefined to include Doraphis. Ceratocallis Qiao & Zhang is established as a junior synonym of Ceratoglyphina van der Goot, syn. nov. Lithoaphis quercisucta Qiao, Guo & Zhang is transferred to the genus Neohormaphis Noordam as Neohormaphis quercisucta (Qiao, Guo & Zhang) comb. nov. Galls have evolved independently within three tribes of Hormaphidinae. In Cerataphidini, pseudogalls are ancestral, both single-cavity and multiple-cavity galls have evolved once, and galls appear to have evolved towards greater complexity. Galling on secondary hosts has evolved twice in hormaphidines.

Comparative analysis of mitochondrial genomes of five aphid species (Hemiptera: Aphididae) and phylogenetic implications

Insect mitochondrial genomes (mitogenomes) are of great interest in exploring molecular evolution, phylogenetics and population genetics. Only two mitogenomes have been previously released in the insect group Aphididae, which consists of about 5,000 known species including some agricultural, forestry and horticultural pests. Here we report the complete 16,317 bp mitogenome of Cavariella salicicola and two nearly complete mitogenomes of Aphis glycines and Pterocomma pilosum. We also present a first comparative analysis of mitochondrial genomes of aphids. Results showed that aphid mitogenomes share conserved genomic organization, nucleotide and amino acid composition, and codon usage features. All 37 genes usually present in animal mitogenomes were sequenced and annotated. The analysis of gene evolutionary rate revealed the lowest and highest rates for COI and ATP8, respectively. A unique repeat region exclusively in aphid mitogenomes, which included variable numbers of tandem repeats in a lineage-specific manner, was highlighted for the first time. This region may have a function as another origin of replication. Phylogenetic reconstructions based on protein-coding genes and the stem-loop structures of control regions confirmed a sister relationship between Cavariella and pterocommatines. Current evidence suggest that pterocommatines could be formally transferred into Macrosiphini. Our paper also offers methodological instructions for obtaining other Aphididae mitochondrial genomes.

DNA barcoding and elucidation of cryptic aphid species (Hemiptera: Aphididae) in India

Rapid, precise and timely identification of invasive pest insects such as aphids is important and a challenge worldwide due to their complex life cycles, parthenogenetic reproduction, sex and colour morphs. In this respect, DNA barcoding employing a 658 bp fragment of 5′ region of the mitochondrial cytochrome oxidase I (CO-I) gene is an effective tool in addressing the above. In the present study, we employed CO-I for discriminating 142 individuals representing 32 species of aphids from India. Sequence analyses revealed that the intraspecific and interspecific distances ranged from zero to 3.8% and 2.31 to 18.9%, respectively. In addition, the study also showed for the first time the prevalence of three cryptic species, namely Brevicoryne brassicae (Linnaeus), Hyperomyzus carduellinus (Theobald) and Brachycaudus helichrysi (Kaltenbach) from India. Our work has clearly demonstrated that DNA barcoding is an efficient and accurate method for identification of aphid species (including cryptic species), an approach that potentially could play an important role in formulating viable pest management strategies, more especially biocontrol.

Evolutionary history of assassin bugs (insecta: hemiptera: Reduviidae): insights from divergence dating and ancestral state reconstruction

Assassin bugs are one of the most successful clades of predatory animals based on their species numbers (∼6,800 spp.) and wide distribution in terrestrial ecosystems. Various novel prey capture strategies and remarkable prey specializations contribute to their appeal as a model to study evolutionary pathways involved in predation. Here, we reconstruct the most comprehensive reduviid phylogeny (178 taxa, 18 subfamilies) to date based on molecular data (5 markers). This phylogeny tests current hypotheses on reduviid relationships emphasizing the polyphyletic Reduviinae and the blood-feeding, disease-vectoring Triatominae, and allows us, for the first time in assassin bugs, to reconstruct ancestral states of prey associations and microhabitats. Using a fossil-calibrated molecular tree, we estimated divergence times for key events in the evolutionary history of Reduviidae. Our results indicate that the polyphyletic Reduviinae fall into 11-14 separate clades. Triatominae are paraphyletic with respect to the reduviine genus Opisthacidius in the maximum likelihood analyses; this result is in contrast to prior hypotheses that found Triatominae to be monophyletic or polyphyletic and may be due to the more comprehensive taxon and character sampling in this study. The evolution of blood-feeding may thus have occurred once or twice independently among predatory assassin bugs. All prey specialists evolved from generalist ancestors, with multiple evolutionary origins of termite and ant specializations. A bark-associated life style on tree trunks is ancestral for most of the lineages of Higher Reduviidae; living on foliage has evolved at least six times independently. Reduviidae originated in the Middle Jurassic (178 Ma), but significant lineage diversification only began in the Late Cretaceous (97 Ma). The integration of molecular phylogenetics with fossil and life history data as presented in this paper provides insights into the evolutionary history of reduviids and clears the way for in-depth evolutionary hypothesis testing in one of the most speciose clades of predators.

Parallel molecular evolution in an herbivore community

Numerous insects have independently evolved the ability to feed on plants that produce toxic secondary compounds called cardenolides and can sequester these compounds for use in their defense. We surveyed the protein target for cardenolides, the alpha subunit of the sodium pump, Na(+),K(+)-ATPase (ATPα), in 14 species that feed on cardenolide-producing plants and 15 outgroups spanning three insect orders. Despite the large number of potential targets for modulating cardenolide sensitivity, amino acid substitutions associated with host-plant specialization are highly clustered, with many parallel substitutions. Additionally, we document four independent duplications of ATPα with convergent tissue-specific expression patterns. We find that unique substitutions are disproportionately associated with recent duplications relative to parallel substitutions. Together, these findings support the hypothesis that adaptation tends to take evolutionary paths that minimize negative pleiotropy.

Host shifts from Lamiales to Brassicaceae in the sawfly genus Athalia

Plant chemistry can be a key driver of host shifts in herbivores. Several species in the sawfly genus Athalia are important economic pests on Brassicaceae, whereas other Athalia species are specialized on Lamiales. These host plants have glucosides in common, which are sequestered by larvae. To disentangle the possible direction of host shifts in this genus, we examined the sequestration specificity and feeding deterrence of iridoid glucosides (IGs) and glucosinolates (GSs) in larvae of five species which either naturally sequester IGs from their hosts within the Plantaginaceae (Lamiales) or GSs from Brassicaceae, respectively. Furthermore, adults were tested for feeding stimulation by a neo-clerodane diterpenoid which occurs in Lamiales. Larvae of the Plantaginaceae-feeders did not sequester artificially administered p-hydroxybenzylGS and were more deterred by GSs than Brassicaceae-feeders were by IGs. In contrast, larvae of Brassicaceae-feeders were able to sequester artificially administered catalpol (IG), which points to an ancestral association with Lamiales. In line with this finding, adults of all tested species were stimulated by the neo-clerodane diterpenoid. Finally, in a phylogenetic tree inferred from genetic marker sequences of 21 Athalia species, the sister species of all remaining 20 Athalia species also turned out to be a Lamiales-feeder. Fundamental physiological pre-adaptations, such as the establishment of a glucoside transporter, and mechanisms to circumvent activation of glucosides by glucosidases are therefore necessary prerequisites for successful host shifts between Lamiales and Brassicaceae.