Prerpints.org logo
Preprint
Article

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

Altmetrics

Downloads

108

Views

65

Comments

0

A peer-reviewed article of this preprint also exists.

supplementary.zip (40.33KB )

This version is not peer-reviewed

Submitted:

24 June 2023

Posted:

26 June 2023

You are already at the latest version

Alerts
Abstract
Expanding 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 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, thereby 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. The combined findings identify aphid and planthopper UV-opsins as a new valuable data sample for studying adaptive opsin evolution.
Keywords: 
Subject: Biology and Life Sciences  -   Insect Science

1. Introduction

Public awareness of hemipteran insects is largely due to the detrimental impact of pest species such as aphids and stink bugs or the horrifying experience of exposure to bed bugs and alike [1,2,3]. Hemiptera, i.e. true bugs, however, diversified into far more than human nuisances. Its members conquered terrestrial and aquatic biotas by extracting nutrients from a diverse panel of animal and plant species with the help of sucking mouthparts, the morphological hallmark of this clade [4]. Select groups evolved traits as remarkable as acoustic long-distance communication (cicadas), water sliding (water striders), and diving (backswimmers and many other families). Our familiarity with hemipteran insects is also due to their expansive number of species. With over 95,000 documented species, the Hemiptera represent the largest order of direct-developing insects. In fact, their abundance is only challenged by the four megadiverse orders of holometabolous insects, ie. the Coleoptera (beetles, ~350,000 species), Diptera (true flies, ~160,000 species), Hymenoptera (bees and wasps, ~150,000 species), and Lepidoptera (moths and butterflies, ~150,000 species) [5].
Phylogenomic studies have begun to consolidate our understanding of hemipteran diversification [4,6,7,8], confirming the existence of three large suborders: the Sternorrhyncha (~12,500 species of aphids, whiteflies, and scale insects), Auchenorrhyncha (~ 42,000species of cicadas, leafhoppers, treehoppers, planthoppers, spittlebugs), and the Heteroptera (~45,000 species). The latter include families of well-known terrestrial species like the kissing bugs (Reduviidae), bedbugs (Cimicidae), plant bugs (Miridae), seed bugs (Lygaeidae), and stink bugs (Pentatomidae and related families) but also aquatic taxa such as water striders (Gerridae) and backswimmers (Notonectidae).
First insights into the genomic dimensions of hemipteran diversification have been gained through the genome sequence assemblies for a number of key species [4], starting with the pea aphid Acyrtosiphon pisum [9], the kissing bug Rhodnius prolixus [2] and the bed bug Cimex lectularius [1]. Subsequent efforts added genome assemblies from specialist species like the water strider Gerris buenoi [10], the established laboratory model Oncopeltus fasciatus [11], and the closely related pentatomid stink bug Halyomorpha halys [12]. Today, over 80 hempiteran genome assemblies are available through the National Center for Biotechnology Information (NCBI) [13] for deeper studies of gene families of interest.
The members of the opsin gene family encode light-sensitive G-protein coupled transmembrane receptors (GPCRs). Their capacity to absorb photon energy functions at the beginning of the phototransduction cascade in animal photoreceptors [14]. Paralogs of this ancient sensory gene family differ by wavelength-dictated photon absorption maxima, thereby providing optimized sensitivity to important spectra of the environment and, in many cases, facilitating the differentiation of wavelength inputs, i.e. color vision. Insects are characterized by three opsin subfamilies whose members are expressed in the photoreceptors of the peripheral visual system, i.e. visual opsins expressed in the compound eyes and the ocelli [15,16,17,18]. This includes the subfamily of green or long wavelength (LW) sensitive opsins with sensitivity peaks centered around 520 nm, the blue (B) sensitive opsins with sensitivity maxima around 450 nm, and UV-sensitive opsins with sensitivity maxima between 360 nm to 380 nm [18].
Color vision generally correlates with the conservation of all of three opsin subfamilies. Crepuscular or nocturnal species, however, are often characterized by the possession of UV- and LW-opsins only, i.e. the evolutionary loss of the B-opsin subfamily. Hemipteran examples of this paradigm include the kissing bug and the bed bug [1]. Less expected, however, was the genomic evidence of the lack of B-opsins in the pea aphid A. pisum [2,3,19]. Equally surprising, also the genome-wide analyses of the opsin repertoires of the water strider, G. buenoi, and the milkweed bug, O. fasciatus, failed to detect B-opsins [10,11]. Conistent with a general absence of B-opsins in the Heteroptera, [20] reported the conservation of LW- and UV-opsins but absence of B-opsin in plant bugs (Heteroptera: Miridae). Moreover, their finding of two LW-opsin paralogs in the plant bugs based on transcriptome evidence aligned well with the earlier reported expansion of the LW-opsin subfamily in the Miridae based on a Heteroptera-wide compilation genomic and transcriptomic opsin sequences [12].
While the lack of B-opsins from hemipteran taxa as distantly related as aphids and Heteroptera raised the possibility of an early loss of this subfamily during hemipteran diversification, the documentation of B-opsin in the green rice leafhopper Nephotettix cincticeps (Auchenorrhyncha: Cicadomorpha) pointed to a more complex scenario [21]. Moreover, transcriptome searches recovered a B-opsin homolog in the hackberry petiole gall psyllid Pachypsylla venusta (Sternorrhyncha) [10]. Corroborating the conservation of B-opsin in the Sternorrhyncha, a B-opsin homolog was recently also reported from the Asian citrus psyllid Diaphorina citri [22]. Combined, these data indicate that the B-opsin subfamily was conserved in early Hemiptera and lost multiple times during later diversification in select subgroups.
Here I present the results from an updated analysis of B-opsin conservation in the Hemiptera. This effort corroborates the existence of B-opsins in select lineages of the Sternorrhyncha and Auchenorrhyncha. In addition, B-opsin was also recovered from a representative of the moss bugs, i.e. the smaller hemipteran suborder Coleorrhyncha. The survey further reveals that the lack of B-opsin in subclades of the Sternorrhyncha and Auchenorrhyncha, i.e. in aphids (Aphidoidea) and planthoppers (Fulgoroidea), respectively, is correlated with the presence of tandem-duplicated UV-opsins. Protein sequence evolution evidence, previously published wavelength sensitivity data, and available data on the expression of opsin paralogs in the compound eye retina conclusively indicate that, in both cases, one clade of UV-opsin sister paralogs neofunctionalized by shifting from ancestral UV- to derived maximum blue sensitivity while the other clade retained ancestral maximum UV-sensitivity. These remarkably similar trajectories contrast with the Heteroptera where the lack of B-opsin is a trait of the entire taxon and accompanied by several independent expansions of the LW-opsin subfamily. Preliminary probing for tuning substitutions in the neofunctionalized B-opsins of aphids and planthoppers suggests a higher number of lineage-specific than parallel amino acid replacement substitutions. Taken together, these findings highlight Hemiptera as a pivotal group for studying the adaptive evolution of opsin GPCRs in addition to the ongoing studies in Coleoptera, Lepidoptera, and Diptera [23,24,25,26,27,28,29].

2. Results

Singleton B-opsin homologs in Sternorrhyncha, Auchenorrhyncha, and Coleorrhyncha

To explore the pervasiveness of B-opsin conservation in the Sternorrhyncha and Auchenorrhyncha, and to further scrutinize their absence from the Heteroptera, I used the B-opsin homologs of N. cincticeps [21] and D. citri [22] as queries in BLAST and tBLASTn searches in current gene, transcript, and genome assembly databases (Table 1). B-opsin family membership was probed by BLAST searches of candidate homologs against the Drosophila melanogaster genome assembly (NCBI Release 6.46) and analysis of relationships in a maximum likelihood estimated gene tree (Figure 1). These efforts identified singleton B-opsin homologs in additional Sternorrhyncha and Auchenorrhyncha species. In the Sternorrhyncha, B-opsins were detected in the hackberry nipplegall maker Pachypsylla celtidismamma and the silverleaf whitefly Bemisia tabaci (Table 1). In the Auchenorrhyncha, BLAST searches uncovered B-opsins in the green leafhopper Empoasca vitis and the tea green leafhopper Matsumurasca onukii, the spittlebug Clastoptera arizonana, and the cicada species Tamasa doddi and Burbunga queenslandica (Table 1). Moreover, B-opsin was also recovered from Xenophysella greensladeae, a member of the relatively small but likewise ancient hemipteran suborder Coleorrhyncha (moss bugs) (Figure 1).
Replicating the outcome of previous efforts, no B-opsins were recovered from heteropteran species [10,12]. Taken together, these findings confirmed the ancestral conservation of the B-opsin gene family in the Hemiptera and documented its continued conservation in the Coleorrhyncha besides select lineages of the Sternorrhyncha and Auchenorrhyncha.

Parallel losses of B-opsin in the Heteroptera, Auchenorrhyncha, and Sternorrhyncha

Parsimony-guided manual gene-species tree reconciliation [30] mapped the loss of B-opsin in the Heteroptera to the stem lineage of this group, as previously suggested [12] (Figure 2). Similarly consistent with previous conclusions, the global opsin gene tree provided also reproduced for the parallel expansions of the LW-opsin gene family in the Gerridae and Pentatomidae (Figure 2) [12].
In the Auchenorrhyncha, the exclusive conservation of B-opsins in members of the Cicadomorpha led to the conclusion that B-opsin was selectively lost in the lineage to the Fulgoroidea, i.e. planthoppers, somewhere between 200-300 million years ago (Figure 2). In the Sternorrhyncha, the conservation of B-opsin in the Psylloidea and Aleyrodidae suggested that the consistent lack of B-opsin in the five sampled aphid species was due its loss in the lineage to Aphidoidea, at least ~100 million years ago (Figure 2). Overall, these findings located a minimum of three parallel B-opsin losses in the hemipteran tree of life.

B-opsin losses correlate with UV-opsin duplications in the aphids and planthoppers

The second insight emerging from the opsin compilation survey and gene-species tree reconciliation was the correlation between B-opsin conservation and the number of UV- or LW-opsin homologs (Table 1 and Figure 2). As noted above, many heteropteran clades are characterized by the possession of tandem duplicated LW-opsin genes (Figure 2) [10,12]. Contrasting this pattern, Auchenorrhyncha and Sternorrhyncha species that lack B-opsins, i.e. planthoppers and aphids, are characterized by duplicated UV-opsin genes (Table 1 and Figure 2).
Previous studies reported the tandem linkage of the duplicated UV-opsins in the pea aphid [19]. The same constellation was encountered in other aphid species for which genome assemblies were available as well as for the B-opsin lacking member of the Auchenorrhyncha, the brown planthopper Nilaparvata lugens (Table 2). In the aphid species, interlocus distances ranged from 21,000bp to 30,000bp (Table 2). The UV-opsin paralogs of N. lugens were separated by less than 6,000 bp.

Accelerated protein sequence evolution rates in the putatively blue-shifted UV-opsin paralogs of aphids and planthoppers

A strong correlation between B-opsin losses and expansions in the LW- or UV-opsin gene families has also been found in beetles [23]. In this case, protein sequence analyses and experimental studies support the model that blue sensitivity is restored by the duplication and compensatory modification of LW- or UV-opsin homologs [31]. Intriguingly, the green peach aphid M. persicae was found to be green, blue, and UV-sensitive [32], despite its lack of B-opsin (Table 1). As a first step to explore the possibly compensatory nature of the UV-opsin duplications in the B-opsin lacking aphids and planthoppers, I probed for unequal protein sequence diversification rates. As one UV-opsin could be hypothesized to maintain ancestral UV peak sensitivity while the other paralog shifted to blue sensitivity as a derived trait, ancestrally functioning UV-opsin paralogs were expected to differ by fewer amino acid differences from the conserved singleton UV-opsins of B-opsin possessing outgroup species.
In support of this scenario, for both the B-opsin lacking aphid and planthopper species, the duplicated UV-opsin homologs clustered into two paralog clades in the global and suborder-specific opsin gene trees (Figure 1 and Figure 3). Moreover, in each case, one of UV-opsin paralog clade was characterized by a longer basal branch than the other cluster. Consistent with these gene tree topologies, analysis of pairwise protein sequence differences corroborated the distinction between presumptive ancestral (UV-anc) and derived (UV-der) sister paralog clusters. In the planthoppers, the UV-anc opsin paralogs differed on average by 0.502 (+/-0.005) substitutions per site from the singleton UV-opsin of M. onukii compared to 0.55 (+/-0.006) per site for the UV-der opsin paralogs. In the aphid species, the UV-anc opsin paralogs differed on average by 0.67 (+/-0.03) substitutions per site from the singleton UV-opsin of D. citri compared to 0.75 (+/-0.01) for the UV-der opsin paralogs (Supplementary data file S1). Relative rate tests, finally, rejected substitution rate constancy for UV-anc opsin and UV-der opsins combined but only marginally for the aphid UV-anc paralogs when the UV-anc and UV-der paralog clades were tested individually (Figure 3).

Protein sequence change at documented tuning sites

At the molecular level, the wavelength sensitivity peaks of specific opsin GPCRs are primarily determined by amino acid residues in the chromophore binding region, also called tuning sites [29]. Consistent with this mechanism, many studies have shown that amino acid residue changes at tuning sites diversified opsin wavelength specificities [31]. To probe the aphid and planthopper UV-der paralogs for tuning site substitutions, I surveyed 24 protein sequence sites that have been previously proposed or shown to represent tuning sites in UV-opsins of other insect orders [31,33].
In the aphid UV-opsin paralogs, consistent non-conservative amino acid residue differences between the UV-anc and UV-der opsin paralogs matched only two previously identified or proposed tuning sites, i.e. sites 90 and 105 as defined by homology to sites in the protein sequence of bovine Rhodopsin (Table 3). Similarly, in the planthopper UV-opsins, only two non-conservative amino acid residue differences between the UV-anc and UV-der opsin paralogs matched previously reported tuning sites. In this case, however, both the UV-anc paralog and the UV-der opsin paralog were each characterized by one derived amino acid residue state, i.e. at positions 293 and 294, respectively (Table 3).
While the overall low number of amino acid residue changes in previously reported tuning sites was notably low, the replacement of positively charged lysine (K) with charge neutral valine (V) at position 90 in the UV-der opsins of aphids was of exceptional interest as previously noted [19]. Comparative data from both vertebrates and invertebrates as well as wavelength sensitivity studies of recombinant opsins in Drosophila have shown that positively charged lysine at this position of the second GPCR transmembrane domain facilitates UV-sensitivity while negatively charged or neutral amino acid residues such as asparagine (N) or glutamate (E) promote blue sensitivity [34]. Consistent with this paradigm, the singleton outgroup UV-opsins and also the UV-anc opsin paralogs of both aphids and planthoppers were occupied by lysine (K) at site 90 (Table 3). The lysine (K) vs valine (V) replacement at position 90 in the UV-der opsins of aphids thus represented the highest confidence UV- to blue sensitivity shift-related amino acid change in the surveyed UV-opsin homolog dataset. At the same time, the contrasting conservation of the ancestral UV-sensitivity associated lysine state at position 90 in both the UV-anc and UV-der opsin paralogs of planthoppers represented evidence of convergent trajectories towards blue sensitivity in the planthopper and aphid UV-opsins at the protein sequence level.

Parallel vs diversified protein sequence changes

Given the evidence for B-opsin replacement through parallel UV-opsin duplicate neofunctionalization in the planthopper and aphids, scanning the UV-opsin paralogs of the two groups globally for parallel protein sequence changes offered the opportunity to identify new high-confidence tuning substitutions. In addition, this analysis was conducted to gauge the relative impacts of parallel vs convergent protein sequence changes during the presumed adaptive transition from UV- to blue sensitivity in the UV-der paralogs of aphids and planthoppers.
In the aphids, 11 sites were characterized by a derived non-conservative amino acid residue change in the UV-der paralogs vs a single derived non-conservative amino acid residue change in the UV-anc paralogs (Table 3). Similarly, in the planthopper UV-opsins, 11 sites were characterized by derived non-conservative amino acid residue states in the UV-der paralogs. The number of four derived non-conservative amino acid residue states in the UV-anc paralogs in the planthoppers, however, was notably higher compared to the aphid UV-anc paralogs. Most importantly, only two of the sites with derived non-conservative amino acid residue states matched between the UV-der paralogs of aphids and planthoppers (Table 3).
One of them, site 47, mapped to the first GPCR transmembrane domain based on multiple sequence alignment correspondence [19]. Phenylalanine (F) was the ancestral amino acid residue state in both aphids and planthoppers. The UV-der opsins of aphids harbored a leucine (L) at this position, while the UV-der opsins of planthoppers carried a tyrosine (Y) (Table 3). The second site was homologous to bovine Rhodopsin site 183, which resides in the second extracellular domain. Tyrosine (Y) represented the ancestral state in the UV-anc paralogs of both aphids and planthoppers which was, in both groups, replaced by phenylalanine (F) in the UV-der paralogs (Table 3).
Taken together, two protein sequence sites with bona fine parallel adaptive amino acid residue changes compared to 18 sites with non-shared putatively adaptive amino acid residue changes in the UV-der opsin paralogs of aphids and planthoppers. These findings corroborated the evidence that the parallel blue sensitivity shifts of the planthopper and aphid UV-der opsins were largely the outcome of convergent evolution at the protein sequence level.

3. Discussion

Elucidating and explaining the diversification of opsin genes is at the center of current endeavors to gain deeper insights into the visual evolution of insects [18,23,24,25,26,27,28,29]. The present study adds two new examples of maintained trichromacy via the neofunctionalization of UV-opsin paralogs that compensate for the evolutionary loss of the B-opsin subfamily in distantly related clades of the Hemiptera.

At least three B-opsin losses during hemipteran diversification

Previous studies reported the lack of B-opsin throughout the Heteroptera [12], in the aphid species A. pisum and M. persicae [19], and in the planthopper N. lugens [21]. In isolation, these findings were compatible with a very early loss of the B-opsin gene family in the stem lineage of the Hemiptera. This straightforward scenario, however, was contradicted by the existence of B-opsin orthologs in the green rice leafhopper N. cincticeps (Auchenorrhyncha: Cicadomorpha) [21], the hackberry petiole gall psyllid P. venusta (Sternorrhyncha) (Armisen et al. 2018), and the Asian citrus psyllid D. citri (Li et al. 2020). The present study documents the existence of B-opsins in further taxa of the Sternorrhyncha (Aleyrodoidea, Psylloidea), Auchenorrhyncha (Cercopoidea, Cicadoidea, Membracoidea), and the Coleorrhyncha (moss bugs). Combined, these findings reveal that the B-opsin gene family had remained conserved in early Hemiptera but was subsequently lost at least three independent times.
Given the current limitation of taxon sampling density, it is possible that future studies will detect additional opsin gene family changes in younger hemipteran taxa. At this point, however, it seems reasonable to assume that the results presented here conclusively capture the number of B-opsin gene losses in the largest hemipteran suborders, i.e. the Heteroptera, Auchenorrhyncha, and the Sternorrhyncha.
Some caution may be warranted concerning the actual number of B-opsin loss events in the Heteroptera. While the most parsimonious explanation for the absence of B-opsins in all Heteroptera is a single loss in the stem lineage of this group, less parsimonious scenarios of multiple early losses cannot be excluded with absolute certainty. In addition to the evolutionary ancientness of these possible events, this is because of the considerable number of independent, and in some cases very likely compensatory, LW-opsin duplications in this clade [12]. While resolving this ambiguity warrants further study, it is clear that the consequences of B-opsin loss in the Heteroptera differ in major interesting ways from aphids and planthoppers.

Timing B-opsin losses and UV-opsin duplications in planthoppers and aphids

In the aphids, the present study sampled five species from three subfamilies, i.e. the Aphidinae, Chaitophorinae, and Lachninae. Most molecular phylogenetic studies place these groups in distantly related branches of the aphid subfamily tree [35,36,37,38,39]. It, therefore, seems reasonable to postulate that the loss of B-opsin and correlated duplication of UV-opsin preceded the last common ancestor of aphids that has been estimated to have diverged about 100 million years ago [6]. At the same time, it remains possible that more comprehensive species sampling may produce evidence of a later loss of B-opsin in this group. Alternatively, the lack of B-opsins in aphids may even date back to a more ancient gene loss event, preceding the last common ancestor of aphids and scale insects (Coccoidea). Preliminary searches for B-opsin in this group failed to recover B-opsin homologs, despite the availability of high-coverage genome assembly data [40]. However, sequence information from additional scale insect taxa will need to become available before firm conclusions can be drawn on this issue.
In the planthopper case, the conservation of B-opsin in the Coleorrhyncha and all other clades of Auchenorrhyncha besides planthoppers (Fulgoroidea), i.e. froghoppers (Cercopoidea), cicadas (Cicadoidea), and treehoppers (Membracoidea) makes it straightforward to map the loss of B-opsin to the early evolution of planthoppers [6]. However, only one of the 13 families of planthoppers was sampled in the present study, i.e. the Delphacidae. Thus, also, in this case, broader species sampling needs to be awaited before this B-opsin loss event can be dated with ultimate confidence.

Combined evidence of B-opsin loss compensation through parallel UV-opsin neofunctionalization in aphids and planthoppers

Gene duplication events can lead to a variety of functionalization outcomes [41]. This includes the reduction of ancestral pleiotropy via differential inheritance of ancestral functions by the descendant paralogs, leading to subfunctionalization, or the acquisition of novel functions by one paralog, i.e. neofunctionalization. The latter outcome is frequently associated with dramatic protein sequence divergence between the emerging paralogs as one paralog changes more slowly under the continued constraints of ancestral functionality while the second paralog experiences a burst of adaptive protein sequence change while optimizing to a new function. In the evolutionary diversification of opsin genes, the neofunctionalization scenario applies most clearly when a newly emerged paralog shifts its ancestral peak sensitivity. In support of this notion, peak sensitivity-shifted neofunctionalized opsin paralogs have been found to differ more dramatically from the singleton homologs in outgroup species than the paralog that preserves ancestral peak sensitivity [42]. Based on this criterion, one UV-opsin paralog, labeled UV-der, of each of the aphid and planthopper UV-opsin sister paralog pairs emerged as putatively wavelength sensitivity shifted given its significantly higher overall amount of protein sequence change compared to the sister paralog. The latter, labeled UV-anc by the same logic, represents the candidate for maintaining ancestral UV-sensitivity.
In the long term, it will be desirable to obtain definitive evidence of compensatory blue-shifted UV-opsin gene family paralogs in the planthoppers and aphids by in vitro spectral sensitivity analyses [31,43]. This will be particularly attractive for exploring the two new candidate tuning sites 47 and 183, which are characterized by parallel amino acid substitutions in the aphid and planthopper UV-der paralogs. At this point, further support for the parallel compensatory UV-opsin neofunctionalization scenarios can be drawn from the findings of behavioral wavelength discrimination and electroretinogram (ERG) studies.
In the Sternorrhyncha, critical data are available from studies in the greenhouse whitefly Trialeurodes vaporariorum [44], the pea aphid A. pisum [45], and the green peach aphid M. persicae [32]. Following up on earlier behavioral and ERG studies that produced evidence of UV- and LW peak sensitivities in the greenhouse whitefly T. vaporariorum compound eyes [46], Stukenberg & Poehling (2019) demonstrated high-resolution behavioral wavelength discrimination assays that T. vaporariorum is trichromatic with peak sensitivities in the UV (340–370nm), blue (480–490nm), and green (510–520nm) [44]. Moreover, their data led to the model that blue-sensitivity boosts the accuracy of settling onto green food plant targets. A key observation leading to this conclusion was the observation of an activity-inhibiting effect of exposure to blue light [44].
While sequence data are not available for T. vaporariorum at this time, its trichromacy aligns well with the ancestral possession of the B-opsin in the Aleyrodidae as documented by the silverleaf whitefly B. tabaci (Table 1). The available ERG data for the B-opsin lacking A. pisum are considered ambiguous regarding the possibility of distinct spectral sensitivities in the blue and UV ranges [45]. ERG studies on the green peach aphid M. persicae, however, have shown that this species is characterized by a trichromatic profile of peak sensitivities, which match those of T. vaporariorum despite the lack of B-opsin [32]. It thus seems safe to assume that the detected blue sensitivity in aphids is supplied by the UV-der opsin paralogs as previously proposed [19].
In the Auchenorhyncha, behavioral and ERG data have been generated for the B-opsin possessing green rice leafhopper N. cincticeps [47] and the B-opsin-lacking brown planthopper N. lugens [21]. Phototaxis assays revealed highest attraction of N. cincticeps to green, blue, and yellow wavelengths [47]. ERG data revealed pronounced sensitivity peaks in the UV- and LW-ranges but the best modeling fit was obtained with a trichromatic composite of wavelength sensitivities centered on UV- (354nm), blue (449nm), and LW (527nm) absorption maxima [47]. The same primary UV- and LW-sensitivity peaks were reported for N. lugens as well as specific and pronounced attraction to blue light (470nm) besides UV (365nm, 385nm), and green (525nm) light [21]. Mirroring the situation in the Sternorrhyncha, these data are compatible with preserved trichromacy in N. lugens through compensatory neofunctionalization of the planthopper UV-der opsin paralog.
Another key requirement for color vision is the differential expression of opsin paralogs with different peak sensitivities in defined photoreceptor subtypes [26]. To date, retinal opsin gene expression data have been reported for N. lugens [21], N. cincticeps [47], and A. pisum [19]. Of these, N. cincticeps represents the ancestral state, possessing singleton homologs of the UV-, B-, and LW-opsin subfamilies [47]. In this species, the singleton LW-opsin homolog is expressed in a defined set of seven photoreceptors per ommatidium. The photoreceptor subtype not expression LW-opsin is characterized by the seemingly stochastic but mutually exclusive expression of the UV- or B-opsin homolog [47]. Taken together, the retinal mosaic defined by the differential expression of LW-, B- and UV-opsin is consistent with the behavioral and physiological evidence of trichromacy in N. cincticeps [47].
In the B-opsin lacking N. lugens and A. pisum, strong expression of the singleton LW-opsin was detected uniformly throughout the retina, suggesting expression in most, if not all, photoreceptors. The UV-opsin paralogs, in contrast, were detected in more selective patterns. In N. lugens, the candidate blue-shifted UV-der paralog (BPH_UVop1: XP_039285232) was detected in the ventral periphery of the compound eye while the ancestral UV-anc paralog (BPH_UVop2: XP_039285228) was detected in a subset of photoreceptors that were more evenly distributed across the retina [21]. In A. pisum, both the UV-der opsin (Ap-SWO2: XP_001951588) and the UV-anc opsin paralog (Ap-SWO3: XP_001951613) were detected in similar punctuated low density of photoreceptors across the retina [19], highly reminiscent of the retinal mosaic of N. cincticeps [47]. While the cellular resolution of the expression studies in N. lugens and A. pisum is too limited for conclusive inferences, the broader LW-opsin and the selective UV-opsin expression patterns in these species are consistent with the existence of LW-, B-, and UV-specific photoreceptors as in the ancestrally organized N. cincticeps [47].

Parallel vs convergent events leading to the compensatory replacements of B-opsin with peak sensitivity shifted UV-opsins in aphids and planthoppers

Evolutionary change can lead to similar outcomes in different lineages by either equally congruent changes in homologous units or qualitatively unrelated changes in non-homologous units. The first type of scenario is commonly referred to as parallel evolution while the second one represents convergent evolution. In the case of molecular units like genes, similar outcomes can emerge in different lineages through the combination of parallel and convergent changes. Applied to the similar deployment of peak sensitivity modified UV-opsins for the replacement of B-opsins, the matching modification of homologous units constitutes parallel changes in aphids and planthopper evolution. At the protein sequence level, however, the amount of changes in non-homologous sites is considerably higher than the two protein sequence sites that changed in both lineages in parallel. Moreover, despite the very likely much younger age, the aphid UV-opsin sister paralogs diverged more dramatically based on the average pairwise evolutionary protein sequence divergence of 0.34 (+/-0.013) than the planthopper UV-opsin sister paralogs with an average pairwise evolutionary protein sequence divergence of 0.19 (+/-0.013) (Supplementary data Table S2).
While these differences may strike surprising given the similarity of the outcome, it is certainly partly explained by the at least 300 million years of evolutionary separation between the aphid and planthopper lineages [6]. The second explanation must reside in the large space of protein sequence variation that can modify the endogenous UV peak sensitivity of the chromophore to blue. The parallel cases of B-opsin loss compensation by UV-opsin gene duplicates thus constitute a valuable resource for elucidating the protein sequence combinatorics underlying adaptive opsin evolution.

4. Materials & methods

Homolog compilation

Previously identified visual opsin genes were extracted from the genome assembly reports of the stink bug H. halys [12], the mikweend bug O. fasciatus [11], the water strider G. buenoi [10], the bed louse C. lectularius [1], the kissing bug R. prolixus [2], and the pea aphid Acyrthosiphon pisum [9,48]. Additional genome assemblies interrogated for hemipteran visual opsins included that of the brown planthopper species N. lugens [49,50] and Laodelphax striatellus [51], the green leafhopper N. cincticeps [52], the Cedar bark aphid Cinara cedri [53], the Asian citrus psyllid D. citri [54], and the silverleaf whitefly B. tabaci [55].
Further homologs were identified through BLASTp searches of the NCBI non-redundant (nr) protein sequence dabase and tBLASTn searches of the NCBI Transcriptome Shotgun Assembly (TSA) database of the NCBI [13,56,57].

Sequence analyses

Protein sequence translations of identified opsin transcripts were produced with “Translate tool” of the Expasy suite hosted by the SIB: Swiss Institute of Bioinformatics [58]. Evolutionary protein sequence divergences were estimated with the MEGA software package [59] using the Jones-Taylor-Thornton (JTT) model and applying a gamma parameter value of 1. Pairwise relative rate tests [60] and maximum likelihood tests of molecular clock constancy were likewise conducted in the MEGA suite [59] on multiple sequence alignments generated with T-Coffee [61] followed by removal of ambiguous sites with TrimAL [62] using the “gappyout” stringency options.

Gene tree reconstruction

The global opsin protein sequence gene tree was estimated with RAxML [63] hosted in the CIPRES Science Gateway environment [64] from a multiple sequence alignment produced by T-Coffee [61] and cleared of ambiguous regions applying Gblocks default settings [65].

Supplementary Materials

The following supporting information can be downloaded at the website of this paper posted on Preprints.org.

Acknowledgments

I thank the anonymous reviewers for their helpful comments on the manuscript.

References

  1. Benoit, J.B.; Adelman, Z.N.; Reinhardt, K.; Dolan, A.; Poelchau, M.; Jennings, E.C.; Szuter, E.M.; Hagan, R.W.; Gujar, H.; Shukla, J.N.; et al. Unique Features of a Global Human Ectoparasite Identified through Sequencing of the Bed Bug Genome. Nat. Commun. 2016, 7, 10165. [Google Scholar] [CrossRef]
  2. Mesquita, R.D.; Vionette-Amaral, R.J.; Lowenberger, C.; Rivera-Pomar, R.; Monteiro, F.A.; Minx, P.; Spieth, J.; Carvalho, A.B.; Panzera, F.; Lawson, D.; et al. Genome of Rhodnius Prolixus, an Insect Vector of Chagas Disease, Reveals Unique Adaptations to Hematophagy and Parasite Infection. Proc. Natl. Acad. Sci. U. S. A. 2015, 112, 14936–14941. [Google Scholar] [CrossRef]
  3. Liu, Q.; Guo, Y.; Zhang, Y.; Hu, W.; Li, Y.; Zhu, D.; Zhou, Z.; Wu, J.; Chen, N.; Zhou, X.-N. A Chromosomal-Level Genome Assembly for the Insect Vector for Chagas Disease, Triatoma Rubrofasciata. Gigascience 2019, 8. [Google Scholar] [CrossRef]
  4. Panfilio, K.A.; Angelini, D.R. By Land, Air, and Sea: Hemipteran Diversity through the Genomic Lens. Curr Opin Insect Sci 2018, 25, 106–115. [Google Scholar] [CrossRef] [PubMed]
  5. Condamine, F.L.; Clapham, M.E.; Kergoat, G.J. Global Patterns of Insect Diversification: Towards a Reconciliation of Fossil and Molecular Evidence? Sci. Rep. 2016, 6, 19208. [Google Scholar] [CrossRef] [PubMed]
  6. Johnson, K.P.; Dietrich, C.H.; Friedrich, F.; Beutel, R.G.; Wipfler, B.; Peters, R.S.; Allen, J.M.; Petersen, M.; Donath, A.; Walden, K.K.O.; et al. Phylogenomics and the Evolution of Hemipteroid Insects. Proc. Natl. Acad. Sci. U. S. A. 2018, 115, 12775–12780. [Google Scholar] [CrossRef] [PubMed]
  7. Kieran, T.J.; Gordon, E.R.L.; Forthman, M.; Hoey-Chamberlain, R.; Kimball, R.T.; Faircloth, B.C.; Weirauch, C.; Glenn, T.C. Insight from an Ultraconserved Element Bait Set Designed for Hemipteran Phylogenetics Integrated with Genomic Resources. Mol. Phylogenet. Evol. 2019, 130, 297–303. [Google Scholar] [CrossRef]
  8. Li, H.; Leavengood, J.M., Jr.; Chapman, E.G.; Burkhardt, D.; Song, F.; Jiang, P.; Liu, J.; Zhou, X.; Cai, W. Mitochondrial Phylogenomics of Hemiptera Reveals Adaptive Innovations Driving the Diversification of True Bugs. Proc. Biol. Sci. 2017, 284, 20171223. [Google Scholar] [CrossRef]
  9. International Aphid Genomics Consortium. Genome Sequence of the Pea Aphid Acyrthosiphon Pisum. PLoS Biol. 2010, 8, e1000313. [Google Scholar] [CrossRef]
  10. Armisen, D.; Rajakumar, R.; Friedrich, M.; Benoit, J.B.; Robertson, H.M.; Panfilio, K.A.; Ahn, S.-J.; Poelchau, M.F.; Chao, H.; Dinh, H.; et al. The Genome of the Water Strider Gerris Buenoi Reveals Expansions of Gene Repertoires Associated with Adaptations to Life on the Water. BMC Genomics 2018, 19, 832. [Google Scholar] [CrossRef]
  11. Panfilio, K.A.; Vargas Jentzsch, I.M.; Benoit, J.B.; Erezyilmaz, D.; Suzuki, Y.; Colella, S.; Robertson, H.M.; Poelchau, M.F.; Waterhouse, R.M.; Ioannidis, P.; et al. Molecular Evolutionary Trends and Feeding Ecology Diversification in the Hemiptera, Anchored by the Milkweed Bug Genome. Genome Biol. 2019, 20, 64. [Google Scholar] [CrossRef]
  12. Sparks, M.E.; Bansal, R.; Benoit, J.B.; Blackburn, M.B.; Chao, H.; Chen, M.; Cheng, S.; Childers, C.; Dinh, H.; Doddapaneni, H.V.; et al. Brown Marmorated Stink Bug, Halyomorpha Halys (Stål), Genome: Putative Underpinnings of Polyphagy, Insecticide Resistance Potential and Biology of a Top Worldwide Pest. BMC Genomics 2020, 21, 227. [Google Scholar] [CrossRef] [PubMed]
  13. Pruitt, K.D.; Tatusova, T.; Maglott, D.R. NCBI Reference Sequence (RefSeq): A Curated Non-Redundant Sequence Database of Genomes, Transcripts and Proteins. Nucleic Acids Res. 2005, 33, D501–D504. [Google Scholar] [CrossRef] [PubMed]
  14. McCulloch, K.J.; Macias-Muñoz, A.; Briscoe, A.D. Insect Opsins and Evo-Devo: What Have We Learned in 25 Years? Philos. Trans. R. Soc. Lond. B Biol. Sci. 2022, 377, 20210288. [Google Scholar] [CrossRef]
  15. Henze, M.J.; Oakley, T.H. The Dynamic Evolutionary History of Pancrustacean Eyes and Opsins. Integr. Comp. Biol. 2015, 55, 830–842. [Google Scholar] [CrossRef] [PubMed]
  16. Roberts, N.S.; Hagen, J.F.D.; Johnston, R.J., Jr. The Diversity of Invertebrate Visual Opsins Spanning Protostomia, Deuterostomia, and Cnidaria. Dev. Biol. 2022, 492, 187–199. [Google Scholar] [CrossRef]
  17. Porter, M.L.; Cronin, T.W.; McClellan, D.A.; Crandall, K.A. Molecular Characterization of Crustacean Visual Pigments and the Evolution of Pancrustacean Opsins. Mol. Biol. Evol. 2007, 24, 253–268. [Google Scholar] [CrossRef]
  18. Briscoe, A.D.; Chittka, L. The Evolution of Color Vision in Insects. Annu. Rev. Entomol. 2001, 46, 471–510. [Google Scholar] [CrossRef]
  19. Collantes-Alegre, J.M.; Mattenberger, F.; Barberà, M.; Martínez-Torres, D. Characterisation, Analysis of Expression and Localisation of the Opsin Gene Repertoire from the Perspective of Photoperiodism in the Aphid Acyrthosiphon Pisum. J. Insect Physiol. 2018, 104, 48–59. [Google Scholar] [CrossRef]
  20. Xu, P.; Lu, B.; Chao, J.; Holdbrook, R.; Liang, G.; Lu, Y. The Evolution of Opsin Genes in Five Species of Mirid Bugs: Duplication of Long-Wavelength Opsins and Loss of Blue-Sensitive Opsins. BMC Ecol Evol 2021, 21, 66. [Google Scholar] [CrossRef]
  21. Matsumoto, Y.; Wakakuwa, M.; Yukuhiro, F.; Arikawa, K.; Hiroaki, N. Attraction to Different Wavelength Light Emitting Diodes (LEDs), the Compound Eye Structure, and Opsin Genes in Nilaparvata Lugens. Jap. J. Appl. Entomol. Zool. 2014, 58, 111–118. [Google Scholar] [CrossRef]
  22. Li, C.; Tian, F.; Lin, T.; Wang, Z.; Liu, J.; Zeng, X. The Expression and Function of Opsin Genes Related to the Phototactic Behavior of Asian Citrus Psyllid. Pest Manag. Sci. 2020, 76, 1578–1587. [Google Scholar] [CrossRef] [PubMed]
  23. Sharkey, C.R.; Fujimoto, M.S.; Lord, N.P.; Shin, S.; McKenna, D.D.; Suvorov, A.; Martin, G.J.; Bybee, S.M. Overcoming the Loss of Blue Sensitivity through Opsin Duplication in the Largest Animal Group, Beetles. Sci. Rep. 2017, 7, 8. [Google Scholar] [CrossRef] [PubMed]
  24. Liénard, M.A.; Bernard, G.D.; Allen, A.; Lassance, J.-M.; Song, S.; Childers, R.R.; Yu, N.; Ye, D.; Stephenson, A.; Valencia-Montoya, W.A.; et al. The Evolution of Red Color Vision Is Linked to Coordinated Rhodopsin Tuning in Lycaenid Butterflies. Proc. Natl. Acad. Sci. U. S. A. 2021, 118, e2008986118. [Google Scholar] [CrossRef]
  25. McCulloch, K.J.; Macias-Muñoz, A.; Mortazavi, A.; Briscoe, A.D. Multiple Mechanisms of Photoreceptor Spectral Tuning in Heliconius Butterflies. Mol. Biol. Evol. 2022, 39, msac067. [Google Scholar] [CrossRef]
  26. Wernet, M.F.; Perry, M.W.; Desplan, C. The Evolutionary Diversity of Insect Retinal Mosaics: Common Design Principles and Emerging Molecular Logic. Trends Genet. 2015, 31, 316–328. [Google Scholar] [CrossRef]
  27. Feuda, R.; Goulty, M.; Zadra, N.; Gasparetti, T.; Rosato, E.; Pisani, D.; Rizzoli, A.; Segata, N.; Ometto, L.; Stabelli, O.R. Phylogenomics of Opsin Genes in Diptera Reveals Lineage-Specific Events and Contrasting Evolutionary Dynamics in Anopheles and Drosophila. Genome Biol. Evol. 2021, 13, evab170. [Google Scholar] [CrossRef]
  28. Wakakuwa, M.; Terakita, A.; Koyanagi, M.; Stavenga, D.G.; Shichida, Y.; Arikawa, K. Evolution and Mechanism of Spectral Tuning of Blue-Absorbing Visual Pigments in Butterflies. PLoS One 2010, 5, e15015. [Google Scholar] [CrossRef]
  29. van der Kooi, C.J.; Stavenga, D.G.; Arikawa, K.; Belušič, G.; Kelber, A. Evolution of Insect Color Vision: From Spectral Sensitivity to Visual Ecology. Annu. Rev. Entomol. 2021, 66, 435–461. [Google Scholar] [CrossRef]
  30. Goodman, M.; Czelusniak, J.; Moore, G.W.; Romero-Herrera, A.E.; Matsuda, G. Fitting the Gene Lineage into Its Species Lineage, a Parsimony Strategy Illustrated by Cladograms Constructed from Clobin Sequences. Syst. Biol. 1979, 28, 132–163. [Google Scholar] [CrossRef]
  31. Sharkey, C.R.; Blanco, J.; Lord, N.P.; Wardill, T.J. Jewel Beetle Opsin Duplication and Divergence Is the Mechanism for Diverse Spectral Sensitivities. Mol. Biol. Evol. 2023, 40, msad023. [Google Scholar] [CrossRef] [PubMed]
  32. Kirchner, S.M.; Döring, T.F.; Saucke, H. Evidence for Trichromacy in the Green Peach Aphid, Myzus Persicae (Sulz.) (Hemiptera: Aphididae). J. Insect Physiol. 2005, 51, 1255–1260. [Google Scholar] [CrossRef] [PubMed]
  33. Lord, N.P.; Plimpton, R.L.; Sharkey, C.R.; Suvorov, A.; Lelito, J.P.; Willardson, B.M.; Bybee, S.M. A Cure for the Blues: Opsin Duplication and Subfunctionalization for Short-Wavelength Sensitivity in Jewel Beetles (Coleoptera: Buprestidae). BMC Evol. Biol. 2016, 16, 107. [Google Scholar] [CrossRef]
  34. Salcedo, E.; Zheng, L.; Phistry, M.; Bagg, E.E.; Britt, S.G. Molecular Basis for Ultraviolet Vision in Invertebrates. J. Neurosci. 2003, 23, 10873–10878. [Google Scholar] [CrossRef]
  35. Biello, R.; Singh, A.; Godfrey, C.J.; Fernández, F.F.; Mugford, S.T.; Powell, G.; Hogenhout, S.A.; Mathers, T.C. A Chromosome-Level Genome Assembly of the Woolly Apple Aphid, Eriosoma Lanigerum Hausmann (Hemiptera: Aphididae). Mol. Ecol. Resour. 2021, 21, 316–326. [Google Scholar] [CrossRef] [PubMed]
  36. Nováková, E.; Hypša, V.; Klein, J.; Foottit, R.G.; von Dohlen, C.D.; Moran, N.A. Reconstructing the Phylogeny of Aphids (Hemiptera: Aphididae) Using DNA of the Obligate Symbiont Buchnera Aphidicola. Mol. Phylogenet. Evol. 2013, 68, 42–54. [Google Scholar] [CrossRef] [PubMed]
  37. Owen, C.L.; Miller, G.L. Phylogenomics of the Aphididae: Deep Relationships between Subfamilies Clouded by Gene Tree Discordance, Introgression and the Gene Tree Anomaly Zone. Syst. Entomol. 2022, 47, 470–486. [Google Scholar] [CrossRef]
  38. Ortiz-Rivas, B.; Martínez-Torres, D. Combination of Molecular Data Support the Existence of Three Main Lineages in the Phylogeny of Aphids (Hemiptera: Aphididae) and the Basal Position of the Subfamily Lachninae. Mol. Phylogenet. Evol. 2010, 55, 305–317. [Google Scholar] [CrossRef]
  39. Chen, J.; Wang, Y.; Jiang, L.; Qiao, G. Mitochondrial Genome Sequences Effectively Reveal Deep Branching Events in Aphids (Insecta: Hemiptera: Aphididae). Zool. Scr. 2017, 46, 706–717. [Google Scholar] [CrossRef]
  40. Wang, W.; Chen, X.; Bashir, N.H.; Lu, Q.; Zhang, J.; Ling, X.; Ding, W.; Chen, H. High-Resolution Lac Insect Genome Assembly Provides Genetic Insights into Lac Synthesis and Evolution of Scale Insects. bioRxiv 2023, 2023.02.05.526168. [Google Scholar] [CrossRef]
  41. Kuzmin, E.; Taylor, J.S.; Boone, C. Retention of Duplicated Genes in Evolution. Trends Genet. 2022, 38, 59–72. [Google Scholar] [CrossRef] [PubMed]
  42. Weadick, C.J.; Chang, B.S.W. Complex Patterns of Divergence among Green-Sensitive (RH2a) African Cichlid Opsins Revealed by Clade Model Analyses. BMC Evol. Biol. 2012, 12, 206. [Google Scholar] [CrossRef] [PubMed]
  43. Liénard, M.A.; Valencia-Montoya, W.A.; Pierce, N.E. Molecular Advances to Study the Function, Evolution and Spectral Tuning of Arthropod Visual Opsins. Philos. Trans. R. Soc. Lond. B Biol. Sci. 2022, 377, 20210279. [Google Scholar] [CrossRef]
  44. Stukenberg, N.; Poehling, H.-M. Blue–green Opponency and Trichromatic Vision in the Greenhouse Whitefly ( Trialeurodes Vaporariorum ) Explored Using Light Emitting Diodes. Ann. Appl. Biol. 2019, 175, 146–163. [Google Scholar] [CrossRef]
  45. Döring, T.F.; Kirchner, S.M.; Skorupski, P.; Hardie, J. Spectral Sensitivity of the Green Photoreceptor of Winged Pea Aphids. Physiol. Entomol. 2011, 36, 392–396. [Google Scholar] [CrossRef]
  46. Mellor, H.E.; Bellingham, J.; Anderson, M. Spectral Efficiency of the Glasshouse Whitefly Trialeurodes Vaporariorum and Encarsia Formosa Its Hymenopteran Parasitoid. Entomol. Exp. Appl. 1997, 83, 11–20. [Google Scholar] [CrossRef]
  47. Wakakuwa, M.; Stewart, F.; Matsumoto, Y.; Matsunaga, S.; Arikawa, K. Physiological Basis of Phototaxis to near-Infrared Light in Nephotettix Cincticeps. J. Comp. Physiol. A Neuroethol. Sens. Neural Behav. Physiol. 2014, 200, 527–536. [Google Scholar] [CrossRef]
  48. Li, Y.; Park, H.; Smith, T.E.; Moran, N.A. Gene Family Evolution in the Pea Aphid Based on Chromosome-Level Genome Assembly. Mol. Biol. Evol. 2019, 36, 2143–2156. [Google Scholar] [CrossRef]
  49. Yin, C.; Shen, G.; Guo, D.; Wang, S.; Ma, X.; Xiao, H.; Liu, J.; Zhang, Z.; Liu, Y.; Zhang, Y.; et al. InsectBase: A Resource for Insect Genomes and Transcriptomes. Nucleic Acids Res. 2016, 44, D801–D807. [Google Scholar] [CrossRef]
  50. Ma, W.; Xu, L.; Hua, H.; Chen, M.; Guo, M.; He, K.; Zhao, J.; Li, F. Chromosomal-Level Genomes of Three Rice Planthoppers Provide New Insights into Sex Chromosome Evolution. Mol. Ecol. Resour. 2021, 21, 226–237. [Google Scholar] [CrossRef]
  51. Zhu, J.; Jiang, F.; Wang, X.; Yang, P.; Bao, Y.; Zhao, W.; Wang, W.; Lu, H.; Wang, Q.; Cui, N.; et al. Genome Sequence of the Small Brown Planthopper, Laodelphax Striatellus. Gigascience 2017, 6, 1–12. [Google Scholar] [CrossRef] [PubMed]
  52. Yan, B.; Yu, X.; Dai, R.; Li, Z.; Yang, M. Chromosome-Level Genome Assembly of Nephotettix Cincticeps (Uhler, 1896) (Hemiptera: Cicadellidae: Deltocephalinae). Genome Biol. Evol. 2021, 13, evab236. [Google Scholar] [CrossRef]
  53. Julca, I.; Marcet-Houben, M.; Cruz, F.; Vargas-Chavez, C.; Johnston, J.S.; Gómez-Garrido, J.; Frias, L.; Corvelo, A.; Loska, D.; Cámara, F.; et al. Phylogenomics Identifies an Ancestral Burst of Gene Duplications Predating the Diversification of Aphidomorpha. Mol. Biol. Evol. 2020, 37, 730–756. [Google Scholar] [CrossRef] [PubMed]
  54. Shippy, T.D.; Miller, S.; Massimino, C.; Vosburg, C.; Hosmani, P.S.; Flores-Gonzalez, M.; Mueller, L.A.; Hunter, W.B.; Benoit, J.B.; Brown, S.J.; et al. Annotating Genes in Diaphorina Citri Genome Version 3. protocols. io 2020. [Google Scholar]
  55. Chen, W.; Hasegawa, D.K.; Kaur, N.; Kliot, A.; Pinheiro, P.V.; Luan, J.; Stensmyr, M.C.; Zheng, Y.; Liu, W.; Sun, H.; et al. The Draft Genome of Whitefly Bemisia Tabaci MEAM1, a Global Crop Pest, Provides Novel Insights into Virus Transmission, Host Adaptation, and Insecticide Resistance. BMC Biol. 2016, 14, 110. [Google Scholar] [CrossRef] [PubMed]
  56. Altschul, S.F.; Koonin, E.V. Iterated Profile Searches with PSI-BLAST—a Tool for Discovery in Protein Databases. Trends Biochem. Sci. 1998, 23, 444–447. [Google Scholar] [CrossRef] [PubMed]
  57. McGinnis, S.; Madden, T.L. BLAST: At the Core of a Powerful and Diverse Set of Sequence Analysis Tools. Nucleic Acids Res. 2004, 32, W20–W25. [Google Scholar] [CrossRef]
  58. Artimo, P.; Jonnalagedda, M.; Arnold, K.; Baratin, D.; Csardi, G.; de Castro, E.; Duvaud, S.; Flegel, V.; Fortier, A.; Gasteiger, E.; et al. ExPASy: SIB Bioinformatics Resource Portal. Nucleic Acids Res. 2012, 40, 597–603. [Google Scholar] [CrossRef]
  59. Tamura, K.; Stecher, G.; Peterson, D.; Filipski, A.; Kumar, S. MEGA6: Molecular Evolutionary Genetics Analysis Version 6.0. Mol. Biol. Evol. 2013, 30, 2725–2729. [Google Scholar] [CrossRef]
  60. Tajima, F. Simple Methods for Testing the Molecular Evolutionary Clock Hypothesis. Genetics 1993, 135, 599–607. [Google Scholar] [CrossRef]
  61. Notredame, C.; Higgins, D.G.; Heringa, J. T-Coffee: A Novel Method for Fast and Accurate Multiple Sequence Alignment. J. Mol. Biol. 2000, 302, 205–217. [Google Scholar] [CrossRef] [PubMed]
  62. Capella-Gutiérrez, S.; Silla-Martínez, J.M.; Gabaldón, T. trimAL: A Tool for Automated Alignment Trimming in Large-Scale Phylogenetic Analyses. Bioinformatics 2009, 25, 1972–1973. [Google Scholar] [CrossRef] [PubMed]
  63. Stamatakis, A. RAxML Version 8: A Tool for Phylogenetic Analysis and Post-Analysis of Large Phylogenies. Bioinformatics 2014, 30, 1312–1313. [Google Scholar] [CrossRef]
  64. Miller, M.A.; Pfeiffer, W.; Schwartz, T. Creating the CIPRES Science Gateway for Inference of Large Phylogenetic Trees. In Proceedings of the 2010 Gateway Computing Environments Workshop (GCE); IEEE, 2010; pp. 1–8. [Google Scholar] [CrossRef]
  65. Talavera, G.; Castresana, J. Improvement of Phylogenies after Removing Divergent and Ambiguously Aligned Blocks from Protein Sequence Alignments. Syst. Biol. 2007, 56, 564–577. [Google Scholar] [CrossRef] [PubMed]
Preprints 77547 g001
Figure 1. Global maximum likelihood gene tree of visual opsins from select hemipteran and outgroup species. Numbers at branches represent nonparametric bootstrap support from 100 replications. tree. Scale bar corresponds to 0.25 substitutions per amino acid site. Opsin homologs from the outgroup species Drosophila melanogaster (Dmel) and Frankliniella occidentalis (Focc) indicated by dark grey bolt font. Abbreviations of included hemipteran species: Apis = Acyrthosiphon pisum, Btab = Bemisia tabaci, Bque = Burbunga queenslandica, Cari = Clastoptera arizonana, Dcit = Diaphorina citri, Cced = Cinara cedri, Evit = Empoasca vitis, Gbue = Gerris buenoi, Hal = Halyomorpha halys, Lstri = Laodelphax striatellus, Ncin = Nephotettix cincticeps, Monu = Matsumurasca onukii, Mper = Myzus persicae, Nlug = Nilaparvata lugens, Pcel = Pachypsylla celtidismamma, Rmai = Rhopalosiphum maidis, Sfla = Sipha flava, Sfur = Sogatella furcifera, Tdod = Tamasa doddi, Xgre = Xenophysella greensladeae.
Figure 1. Global maximum likelihood gene tree of visual opsins from select hemipteran and outgroup species. Numbers at branches represent nonparametric bootstrap support from 100 replications. tree. Scale bar corresponds to 0.25 substitutions per amino acid site. Opsin homologs from the outgroup species Drosophila melanogaster (Dmel) and Frankliniella occidentalis (Focc) indicated by dark grey bolt font. Abbreviations of included hemipteran species: Apis = Acyrthosiphon pisum, Btab = Bemisia tabaci, Bque = Burbunga queenslandica, Cari = Clastoptera arizonana, Dcit = Diaphorina citri, Cced = Cinara cedri, Evit = Empoasca vitis, Gbue = Gerris buenoi, Hal = Halyomorpha halys, Lstri = Laodelphax striatellus, Ncin = Nephotettix cincticeps, Monu = Matsumurasca onukii, Mper = Myzus persicae, Nlug = Nilaparvata lugens, Pcel = Pachypsylla celtidismamma, Rmai = Rhopalosiphum maidis, Sfla = Sipha flava, Sfur = Sogatella furcifera, Tdod = Tamasa doddi, Xgre = Xenophysella greensladeae.
Preprints 77547 g001
Preprints 77547 g002
Figure 2. Gene species reconciliation tree of visual opsin gene family evolution in the Hemiptera. Phylogenetic relationships and clade ages based on Johnson et al. (2018). Purple, blue, and green phylograms represent the UV-, B-, and LW-opsin subfamily gene trees. Stop signs indicate hypothesized time points of inferred B-opsin loss events. Burst signs indicate hypothesized time points of inferred UV-opsin duplication events. Mya = Million years ago.
Figure 2. Gene species reconciliation tree of visual opsin gene family evolution in the Hemiptera. Phylogenetic relationships and clade ages based on Johnson et al. (2018). Purple, blue, and green phylograms represent the UV-, B-, and LW-opsin subfamily gene trees. Stop signs indicate hypothesized time points of inferred B-opsin loss events. Burst signs indicate hypothesized time points of inferred UV-opsin duplication events. Mya = Million years ago.
Preprints 77547 g002
Preprints 77547 g003
Figure 3. Protein sequence diversification rate differences between ancestral and putatively blue-shifted UV-opsin paralogs of planthoppers and aphids. Maximum likelihood gene trees of the (a) aphid and (b) planthopper UV-opsin clades. Purple and blue branches highlight the presumed ancestrally UV-sensitive UV-anc and blue-sensitivity shifted UV-der paralog lineages, respectively. Scale bars represent 0.01 substitutions per site. See Figure 1 for species abbreviations. (a) Tested individually for molecular clock compatible divergence rates, the UV-opsin homologs in the aphid UV-anc and UV-der paralog groups have diverged at significantly (p=0.023) and non-significantly (p = 0.59) different rates (ns). Clock-like protein sequence divergence is also rejected when the members of both clades are tested together (p=0.026). The singleton UV-opsin homologs of Pachypsylla celtidismamma (Pcel), Diaphorina citri (Dcit), and Frankliniella occidentalis (Focc) (Thysanoptera) were included as outgroup sequences in the molecular clock tests. (b) Members of the planthopper UV-anc and UV-der paralog groups have not significantly diverged from other members within their paralog groups based pairwise relative rate test results (p>0.5). The singleton UV-opsin homolog of Empoasca vitis (Evit) was used as outgroup in the pairwise relative rate tests. However, clock-like protein sequence divergence is rejected when the members of both clades are tested together (p=0.0068). The singleton UV-opsin homologs of Nephotettix cincticeps (Ncin), Empoasca vitis, Matsumurasca onukii (Monu), and Frankliniella occidentalis were included as outgroup sequences in the molecular clock tests.
Figure 3. Protein sequence diversification rate differences between ancestral and putatively blue-shifted UV-opsin paralogs of planthoppers and aphids. Maximum likelihood gene trees of the (a) aphid and (b) planthopper UV-opsin clades. Purple and blue branches highlight the presumed ancestrally UV-sensitive UV-anc and blue-sensitivity shifted UV-der paralog lineages, respectively. Scale bars represent 0.01 substitutions per site. See Figure 1 for species abbreviations. (a) Tested individually for molecular clock compatible divergence rates, the UV-opsin homologs in the aphid UV-anc and UV-der paralog groups have diverged at significantly (p=0.023) and non-significantly (p = 0.59) different rates (ns). Clock-like protein sequence divergence is also rejected when the members of both clades are tested together (p=0.026). The singleton UV-opsin homologs of Pachypsylla celtidismamma (Pcel), Diaphorina citri (Dcit), and Frankliniella occidentalis (Focc) (Thysanoptera) were included as outgroup sequences in the molecular clock tests. (b) Members of the planthopper UV-anc and UV-der paralog groups have not significantly diverged from other members within their paralog groups based pairwise relative rate test results (p>0.5). The singleton UV-opsin homolog of Empoasca vitis (Evit) was used as outgroup in the pairwise relative rate tests. However, clock-like protein sequence divergence is rejected when the members of both clades are tested together (p=0.0068). The singleton UV-opsin homologs of Nephotettix cincticeps (Ncin), Empoasca vitis, Matsumurasca onukii (Monu), and Frankliniella occidentalis were included as outgroup sequences in the molecular clock tests.
Preprints 77547 g003
Table 1. Compilation of hemipteran visual opsin genes.
Table 1. Compilation of hemipteran visual opsin genes.
Species Diaphorina citri Superfamily Family LW Opsins B Opsins UV Opsins References
Nilaparvata lugens Auchenorrhyncha Fulgoroidea Delphacidae XP_039295943 - XP_039285232
XP_039285228		 Matsumoto et al.(2014)
Ma et al. (2021)
Laodelphax striatellus Auchenorrhyncha Fulgoroidea Delphacidae RZF36139 - BAO03863
BAO03862		 Ma et al. (2021)
Zhu et al. (2017)
Matsumoto et al.(2014)
Sogatella furcifera Auchenorrhyncha Fulgoroidea Delphacidae BAO03858 - BAO03860
BAO03859		 Matsumoto et al.(2014)
Ma et al. (2021)
Clastoptera arizonana Auchenorrhyncha Cercopoidea Clastopteridae GEDC01016786 GEDC01026123 GEDC01015247 BioProject: PRJNA303152
Tamasa doddi Auchenorrhyncha Cicadoidea Cicadidae GEIK01032296 GEIK01032035 GEIK01048237 Johnson et al. (2018)
Burbunga queenslandica Auchenorrhyncha Cicadoidea Cicadidae GEJZ01024619 GEJZ01025003 GEJZ01025175 Johnson et al. (2018)
Nephotettix cincticeps Auchenorrhyncha Membracoidea Cicadellidae BAO03864 BAO03865 BAO03866 Matsumoto et al.(2014)
Empoasca vitis Auchenorrhyncha Membracoidea Cicadellidae AHK12248 AIY53970 AHN82368 unpublished
Matsumurasca onukii Auchenorrhyncha Membracoidea Cicadellidae QCC70516 QCC70517 QCC70518 unpublished
Xenophysella greensladeae Coleorrhyncha - Peloridiidae GAYI02055701 GAYI02059430 GAYI02063107 Misof et al. (2014)
Gerris buenoi Heteroptera Gerroidea Gerridae LW1 GBUE012030
LW2 GBUE012031
LW3 GBUE012032
LW4 GBUE012033		 - GBUE021384 Armisen et al. (2018)
Halyomorpha halys Heteroptera Pentatomoidea Pentatomidae XP_014274078
XP_014274077
XP_014274075		 - XP_014273911 Sparks et al. (2020)
Acyrthosiphon pisum Sternorrhyncha Aphidoidea Aphididae XP_008185403 - XP_001951588
XP_001951613		 Li Y et al. (2019)
Myzus persicae Sternorrhyncha Aphidoidea Aphididae XP_022159972 - XP_022162222
XP_022162214		 Rispe et al. (2008)
Sipha flava Sternorrhyncha Aphidoidea Aphididae XP_025414175 - XP_025410308
XP_025410414		 BioProject PRJNA479456
Cinara cedri Sternorrhyncha Aphidoidea Aphididae LW_VVC43438 - VVC32069
VVC32067		 Julca et al. (2020)
Rhopalosiphum maidis Sternorrhyncha Aphidoidea Aphididae XP_026807241 - XP_026808077
XP_026807991		 BioProject PRJNA503266
Bemisia tabaci Sternorrhyncha Aleyrodoidea Aleyrodidae XP_018901683 XP_018913779 XP_018909731 Chen et al. (2016)
Pachypsylla celtidismamma Sternorrhyncha Psylloidea Aphalaridae GCYC01060233 GCYC01076696 GCYC01081548 Johnson et al. (2018)
Diaphorina Sternorrhyncha Psylloidea Liviidae GDFN01027368 XP_008487287 KAI5717139 Johnson et al. (2018)
Table 2. Genomic linkages of duplicated planthopper and aphid UV-opsin paralogs.
Table 2. Genomic linkages of duplicated planthopper and aphid UV-opsin paralogs.
Genus UV-Opsins Locus ID Linkage Group start (bp) end (bp) interlocus distance (bp) # intervening genes
Nilaparvata lugens XP_039285228 LOC111058697 5 75,381,302 75,410,087
XP_039285232 LOC111051439 5 75,416,027 75,437,887 5,940 0
Acyrthosiphon pisum XP_001951588 LOC100163348 A3 27,235,720 27,247,458
XP_001951613 LOC100161312 A3 27,217,337 27,226,019 30,121 0
Myzus persicae XP_022162222 LOC111028006 NW_019101158 450,218 461,548
XP_022162214 LOC111027999 NW_019101158 432,490 442,835 29,058 1
Rhopalosiphum maidis XP_026808077 LOC113550447 NC_040877 46,583,828 46,592,054
XP_026807991 LOC113550397 NC_040878 46,570,436 46,578,338 21,618 0
Table 3. Candidate tuning site amino acid changes in aphid and planthopper UV-opsin paralogs. Site numbers represent homologous sites in bovine Rhodopsin (NP_001014890) or, when denoted with an asterisk, UV-opsin 2 (WCQ76394) of the jewel beetle species Chrysochroa rajah [31]. Sites that correspond to sites that have been previously reported to affect spectral sensitivity are indicated in bolt print in column 1. Light grey background indicates derived amino acid residue changes in putatively ancestrally UV-sensitive opsin paralogs (UV-anc). Cells with dark grey backgrounds and white bolt font indicate derived amino acid residue changes in the putatively derived blue-sensitive UV-opsin paralogs (UV-der). Ancestral amino acid residue states in aphids defined by amino acid residue states in the singleton UV-opsin orthologs of P. celtidismamma and D. citri (see Figure 3a). Ancestral states in planthopper UV-opsins defined by amino acid residue states in the singleton UV-opsin orthologs of N. cincticeps, E. vitis, M. onukii (see Figure 3b). See Supplementary data file S2 for a complete table of investigated sites and Supplementary data file S3 for annotated multiple sequence alignment.
Table 3. Candidate tuning site amino acid changes in aphid and planthopper UV-opsin paralogs. Site numbers represent homologous sites in bovine Rhodopsin (NP_001014890) or, when denoted with an asterisk, UV-opsin 2 (WCQ76394) of the jewel beetle species Chrysochroa rajah [31]. Sites that correspond to sites that have been previously reported to affect spectral sensitivity are indicated in bolt print in column 1. Light grey background indicates derived amino acid residue changes in putatively ancestrally UV-sensitive opsin paralogs (UV-anc). Cells with dark grey backgrounds and white bolt font indicate derived amino acid residue changes in the putatively derived blue-sensitive UV-opsin paralogs (UV-der). Ancestral amino acid residue states in aphids defined by amino acid residue states in the singleton UV-opsin orthologs of P. celtidismamma and D. citri (see Figure 3a). Ancestral states in planthopper UV-opsins defined by amino acid residue states in the singleton UV-opsin orthologs of N. cincticeps, E. vitis, M. onukii (see Figure 3b). See Supplementary data file S2 for a complete table of investigated sites and Supplementary data file S3 for annotated multiple sequence alignment.
Aphids Planthoppers
Site ancestral UV-anc UV-der ancestral UV-anc UV-der
15 L L M L L L
29* P A P T T L
33 Y Y Y Y Y H
45* E E E A E A
47* D S R E D E
43 L M S I/L L V
47 F F L F F Y
52 V/I/L L V L L F
63 L/I C C C C S
84 F F F F F L
85 V/M V/L/C S M/L M/L L
87 M M V M M M
90 K K V K K K
94 F F L F F F
102 G G K(T) G G G
105 Q Q P L T/S A
125 G/S S G A S S
183 Y Y F Y Y F
201 R R,K K R Q R
204 V L/V V V V T
244 (D)S Q Q Q Q M
264 S S/A A S S A
273 L M M(L) L M L
286 G G I G G G/V
293 C I V L C L
294 C,T,A F F T T A
295 C C A C C C
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.
Prerpints.org logo

Preprints.org is a free preprint server supported by MDPI in Basel, Switzerland.

facebook logotwitter logolinkedin logo
MDPI Initiatives
Important Links
Subscribe

Disclaimer

Terms of Use

Privacy Policy

© 2024 MDPI (Basel, Switzerland) unless otherwise stated