
==== Front
Genome Biol Evol
Genome Biol Evol
gbe
Genome Biology and Evolution
1759-6653
Oxford University Press UK

39176990
10.1093/gbe/evae185
evae185
Letter
AcademicSubjects/SCI01130
AcademicSubjects/SCI01140
Evolution of Opsin Genes in Caddisflies (Insecta: Trichoptera)
https://orcid.org/0000-0001-6706-9190
Powell Ashlyn Department of Plant and Wildlife Sciences, Brigham Young University, Provo, UT, USA

https://orcid.org/0000-0001-8771-9154
Heckenhauer Jacqueline LOEWE Centre for Translational Biodiversity Genomics, Frankfurt, Germany
Senckenberg Research Institute and Natural History Museum Frankfurt, Frankfurt, Germany

https://orcid.org/0000-0002-6451-3425
Pauls Steffen U LOEWE Centre for Translational Biodiversity Genomics, Frankfurt, Germany
Senckenberg Research Institute and Natural History Museum Frankfurt, Frankfurt, Germany

https://orcid.org/0000-0002-3921-0908
Ríos-Touma Blanca Facultad de Ingenierías y Ciencias Aplicadas, Ingeniería Ambiental, Grupo de Investigación en Biodiversidad, Medio Ambiente y Salud, Universidad de Las Américas, Quito, Ecuador

https://orcid.org/0000-0002-6353-0450
Kuranishi Ryoichi B Graduate School of Science, Chiba University, Chiba, Japan
Kanagawa Institute of Technology, Kanagawa, Japan

https://orcid.org/0000-0003-1853-6340
Holzenthal Ralph W Department of Entomology, University of Minnesota, St Paul, MN, USA

https://orcid.org/0000-0002-7554-0816
Razuri-Gonzales Ernesto Senckenberg Research Institute and Natural History Museum Frankfurt, Frankfurt, Germany

https://orcid.org/0000-0001-6224-6303
Bybee Seth Department of Biology, Brigham Young University, Provo, UT, USA

https://orcid.org/0000-0002-4801-7579
Frandsen Paul B Department of Plant and Wildlife Sciences, Brigham Young University, Provo, UT, USA

Wheat Christopher Associate Editor
Corresponding authors: E-mails: ashlynpowell913@gmail.com; paul_frandsen@byu.edu.
9 2024
23 8 2024
23 8 2024
16 9 evae18519 8 2024
06 9 2024
© The Author(s) 2024. Published by Oxford University Press on behalf of Society for Molecular Biology and Evolution.
2024
https://creativecommons.org/licenses/by/4.0/ This is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0/), which permits unrestricted reuse, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Insects have evolved complex and diverse visual systems in which light-sensing protein molecules called “opsins” couple with a chromophore to form photopigments. Insect photopigments group into three major gene families based on wavelength sensitivity: long wavelength (LW), short wavelength (SW), and ultraviolet wavelength (UV). In this study, we identified 123 opsin sequences from whole-genome assemblies across 25 caddisfly species (Insecta: Trichoptera). We discovered the LW opsins have the most diversity across species and form two separate clades in the opsin gene tree. Conversely, we observed a loss of the SW opsin in half of the trichopteran species in this study, which might be associated with the fact that caddisflies are active during low-light conditions. Lastly, we found a single copy of the UV opsin in all the species in this study, with one exception: Athripsodes cinereus has two copies of the UV opsin and resides within a clade of caddisflies with colorful wing patterns.

caddisflies
evolution
opsins
visual systems
Trichoptera
insects
BYU College of Life Sciences Deutsche Forschungsgemeinschaft 10.13039/501100001659 502865717 Dirección General de Investigación Universidad de Las Américas 10.13039/100020906 National Science Foundation 10.13039/100000001 2217155
==== Body
pmcSignificance

While opsin evolution in some insect groups has been well-characterized, it has never been studied across caddisflies. Our findings provide insight into the diversity of opsins in caddisflies and form a basis for further research into the evolutionary drivers and complex visual systems in Trichoptera.

Introduction

Within the visual system, the ability to perceive light is critical and plays an essential role in the life histories of insects, including finding food, avoiding predators, and selecting a mate (van der Kooi et al. 2021). Light perception occurs primarily within three different types of visual organs in insects: the stemmata of larvae and the ocelli and compound eyes of adults (van der Kooi et al. 2021; Guignard et al. 2022). Upon light absorption, photoreceptors within the eyes—which contain opsin proteins and chromophores—change their configuration from a resting state to a signaling state, thereby indicating a physiological response (Shichida and Matsuyama 2009). Insect visual opsins form three major gene clades based on their peak wavelength sensitivity, namely, long wavelength (LW; 500 to 600 nm), short wavelength (SW; 400 to 500 nm), and ultraviolet wavelength (UV; 300 to 400 nm; Feuda et al. 2016; Lord et al. 2016; van der Kooi et al. 2021). Many insect groups possess an additional opsin type, Rhodopsin 7 (RH7), which does not have a known function in most insect groups but was found to be involved in circadian rhythms in Drosophila (Senthilan and Helfrich-Förster 2016; Ni et al. 2017).

Insects typically possess one or more copies of each opsin type. Moreover, multiple cases of gene duplications and losses have been observed throughout insect opsin evolution (Spaethe and Briscoe 2004; Frentiu et al. 2007; Sison-Mangus et al. 2008; Briscoe et al. 2010; French et al. 2015; Futahashi et al. 2015; Feuda et al. 2016; Lord et al. 2016; Sharkey et al. 2017; Suvorov et al. 2017; Sondhi et al. 2021; Guignard et al. 2022; Friedrich 2023; Mulhair et al. 2023). These duplications, and subsequent diversification, of visual opsin genes are the primary mechanisms of evolution that lead to greater visual capacity and flexibility (Frentiu et al. 2007; Suvorov et al. 2017; Friedrich 2023) and are usually linked to a particular life-history strategy, living environment, or light condition (French et al. 2015; Futahashi et al. 2015; Lord et al. 2016; Sharkey et al. 2017, 2023; Sondhi et al. 2021; Guignard et al. 2022). While the evolution of opsins has been relatively well studied in some insect orders, opsin genes in caddisflies (Insecta: Trichoptera) have only been characterized in a single species as part of a broader comparative study across insects (Guignard et al. 2022).

As eggs, larvae, and pupae, caddisflies mainly inhabit the benthic zone of freshwater habitats, but as adults, they occupy terrestrial environments adjacent to freshwater (Morse et al. 2019). There are two monophyletic suborders within Trichoptera characterized by differences in habitat, morphology, and silk use: Annulipalpia (retreat making) and Integripalpia (cocoon and case making; Frandsen et al. 2024). Adult caddisflies resemble small moths; yet, while most species have wings and bodies covered in small hairs instead of scales, a few species have brightly colored wings with red, orange, green, or silver regions (Fig. 1) due to the development of hairs into scales (Holzenthal et al. 2007). Additionally, adult caddisflies possess varying eye sizes, some much larger than others (Fig. 1). Presumably, such varied environments—both aquatic and terrestrial—and wing colorations require a plastic and diverse visual system. To gain an understanding of opsin evolution in caddisflies, we analyzed the occurrence and phylogenetic relationships of opsin genes from whole-genome assemblies across 25 caddisfly species (Table 1), representing the major evolutionary lineages within the order.

Fig. 1. Adult caddisfly illustrations showing varying eye sizes and diverse wing colorations and patterns. a) Athripsodes cinereus (Leptoceridae); b) N. paramo (Leptoceridae); c) Nectopsyche nigricapilla (Leptoceridae); d) Nectopsyche ortizi (Leptoceridae); e) Banyallarga vicaria (Calamoceratidae); f) Phylloicus abdominalis (Calamoceratidae). *Species included in this study. Illustrations by Julie Martinez and Ralph Holzenthal.

Table 1 Locations and quality scores of genome assemblies for each species

Species	Accession number	Publication	Contig N50 (bp)	Compleasm (recovered BUSCOs)	
Leptonema lineaticorne	GCA_024500535.1	Heckenhauer et al. (2023)	14,931,587	C: 98.30% [S: 97.83%, D: 0.47%],
F: 0.66%, I: 0.00%, M: 1.04%	
Hydropsyche tenuis	GCA_009617725.1	Heckenhauer et al. (2019)	2,190,134	C: 97.93% [S: 97.65%, D: 0.28%],
F: 0.89%, I: 0.00%, M: 1.18%	
Parapsyche elsis	GCA_022651745.1	Frandsen et al. (2019), Heckenhauer et al. (2022)	5,591,679	C: 96.99% [S: 96.80%, D: 0.19%],
F: 0.80%, I: 0.00%, M: 2.21%	
Arctopsyche grandis	GCA_029955255.1	Frandsen et al. (2023)	6,470,670	C: 98.36% [S: 95.39%, D: 2.97%],
F: 0.71%, I: 0.00%, M: 0.94%	
Plectrocnemia conspersa	JAUTWC000000000	New	32,103,979	C: 97.74% [S: 97.32%, D: 0.42%],
F: 0.61%, I: 0.00%, M: 1.65%	
Philopotamus ludificatus	GCA_022495035.1	Heckenhauer et al. (2022)	35,449	C: 94.63% [S: 93.31%, D: 1.32%],
F: 3.15%, I: 0.05%, M: 2.17%	
Stenopsyche tienmushanensis	GCA_008973525.1	Luo et al. (2018)	1,296,863	C: 97.50% [S: 95.15%, D: 2.35%],
F: 1.04%, I: 0.00%, M: 1.46%	
Agraylea sexmaculata	GCA_022606485.1	Heckenhauer et al. (2022)	86,524	C: 95.67% [S: 90.96%, D: 4.71%],
F: 1.04%, I: 0.00%, M: 3.30%	
Glossosoma conforme	GCA_022606575.1	Heckenhauer et al. (2022)	2,212,131	C: 93.13% [S: 92.47%,D: 0.66%],
F: 0.66%, I: 0.00%, M: 6.21%	
Atopsyche davidsoni	GCA_022113835.1	Ríos-Touma et al. (2021)	14,095,054	C: 98.68% [S: 98.35%, D: 0.33%],
F: 0.61%, I: 0.00%, M: 0.71%	
Atopsyche callosa	Available on FigShare	New	25,586,909	C: 98.54% [S: 96.70%, D: 1.84%],
F: 0.66%, I: 0.00%, M: 0.80%	
Himalopsyche phryganea	GCA_022494535.1	Heckenhauer et al. (2019)	4,634,010	C: 98.21% [S: 97.83%,D: 0.38%],
F: 0.71%, I: 0.00%, M: 1.08%	
Himalopsyche tibetana	GCA_030503985.1	Heckenhauer et al. (2022)	28,889,006	C: 98.91% [S: 98.35%, D: 0.56%],
F: 0.56%, I: 0.00%, M: 0.52%	
Rhyacophila brunnea	Available on FigShare	New	1,306,779	C: 97.74% [S: 94.35%, D: 3.39%],
F: 1.22%, I: 0.00%, M: 1.04%	
Micrasema minimum	GCA_022494985.1	Heckenhauer et al. (2019)	69,526	C: 66.34% [S: 66.20%, D: 0.14%],
F: 7.25%, I: 0.00%, M: 26.41%	
Lepidostoma basale	GCA_022606425.2	Heckenhauer et al. (2019)	1,001,566	C: 98.16% [S: 97.55%, D: 0.61%],
F: 0.85%, I: 0.05%, M: 0.94%	
Drusus annulatus	GCA_022651775.1	Heckenhauer et al. (2019)	1,032,046	C: 97.92% [S: 97.50%, D: 0.42%],
F: 1.08%, I: 0.00%, M: 0.99%	
Halesus radiatus	GCA_022606495.2	Heckenhauer et al. (2019)	124,280	C: 91.20% [S: 89.69%, D: 1.51%],
F: 4.38%, I: 0.05%, M: 4.38%	
Hesperophylax magnus	GCA_026573805.1	Frandsen et al. (2023)	11,205,906	C: 98.68% [S: 96.56%, D: 2.12%],
F: 0.61%, I: 0.00%, M: 0.71%	
Agrypnia vestita	GCA_016648135.1	Olsen et al. (2021)	111,757	C: 91.25% [S: 83.10%, D: 8.15%],
F: 4.14%, I: 0.09%, M: 4.52%	
Eubasilissa regina	GCA_022840565.1	Kawahara et al. (2022)	29,378,647	C: 98.73% [S: 96.28%, D: 2.45%],
F: 0.80%, I: 0.00%, M: 0.47%	
Nectopsyche paramo	JAWQED000000000	New	1,090,281	C: 97.31% [S: 95.24%, D: 2.07%],
F: 0.80%, I: 0.00%, M: 1.88%	
Athripsodes cinereus	GCA_947579605.1	Darwin Tree of Life Project	948,465	C: 97.22% [S: 96.61%, D: 0.61%],
F: 1.08%, I: 0.00%, M: 1.69%	
Odontocerum albicorne	GCA_949825065.1	Darwin Tree of Life Project	387,033	C: 97.22% [S: 96.75%, D: 0.47%],
F: 1.18%, I: 0.00%, M: 1.60%	
Limnocentropus insolitus	Available on FigShare	New	33,230,923	C: 98.87% [S: 98.26%, D: 0.61%],
F: 0.56%, I: 0.00%, M: 0.56%	
Compleasm was run with the Endopterygota OrthoDB v10 BUSCO gene set and the results categories are as follows: C, complete; S, single; D, duplicated; F, fragmented Subclass 1 (only a portion of the gene is present in the assembly); I, fragmented Subclass 2 (different sections of the gene align to different locations in the assembly); M, missing.

Results

Opsin Distribution and Gene Tree

The number of opsin paralogs found within each species ranged from three to as many as nine (Fig. 3). Opsin diversity within Integripalpia (Fig. 3, Clade 2) was more variable than that of Annulipalpia (Fig. 3, Clade 1). For example, the number of LW opsins among species of Integripalpia ranged from one to five paralogs. While we recovered all three visual opsins in four of the seven species within “basal Integripalpia” (Fig. 3, Clade 2a), the three remaining species from the family Rhyacophilidae were found to have lost the SW opsin. Conversely, in the tube-case-making Integripalpia (Fig. 3, Clade 2b), the SW opsin was lost in all except Nectopsyche paramo and Athripsodes cinereus, both of which are from the family Leptoceridae.

The opsin sequences formed four distinct clades in the gene tree corresponding to the LW, SW, UV, and RH7 opsin groups (Fig. 2). Within each opsin clade, the arrangements of suborders were maintained, but occasionally, interspecific relationships of opsin sequences did not match those in the established species tree (Frandsen et al. 2024). We found only minor differences between the coding sequence (CDS) and the peptide gene trees, primarily within the LW and UV opsin groups (supplementary fig. S2, Supplementary Material online). However, none of these areas of incongruence occurred across opsin classes and they were in areas of the trees with lower bootstrap support values (supplementary fig. S2, Supplementary Material online). The LW opsins formed two distinct clades in the opsin gene tree, a phenomenon observed in other insect orders (Futahashi et al. 2015; Feuda et al. 2016; Lord et al. 2016; Sharkey et al. 2017, 2023; Suvorov et al. 2017; Sondhi et al. 2021; Guignard et al. 2022, Mulhair et al. 2023). All species had at least 1 LW opsin in the LW1 clade, while only 14 species—distributed across both suborders—had an additional LW2 opsin (supplementary fig. S1, Supplementary Material online).

Fig. 2. Opsin gene tree. Maximum likelihood tree from caddisfly and outgroup opsin DNA sequences. Nodes are labeled with bootstrap values and branches are colored by opsin type. The letters at the end of the node labels (e.g. LW1a, LW1b, etc.) denote multiple copies of that opsin type in a species. Clades marked with a dot (A to E) were tested for episodic diversifying selection with BUSTED (supplementary table S3, Supplementary Material online).

Lineage-Specific Gene Duplication Events

We observed a few instances of paralogs that were paired in the opsin gene tree (Fig. 2) and adjacent to each other on the same contig in the genome assembly (supplementary table S1, Supplementary Material online), suggesting independent tandem duplication events. This was true for copies of the UV opsin in A. cinereus and the LW1 opsin in Stenopsyche tienmushanensis, Leptonema lineaticorne, Plectrocnemia conspersa, Hesperophylax magnus, Eubasilissa regina, O. albicorne, and A. cinereus (Fig. 2). Upon testing the branches of these duplicate opsins for selection with the Branch-Site Unrestricted Statistical Test for Episodic Diversification (BUSTED; Murrell et al. 2015; Kosakovsky Pond et al. 2020), we found evidence of episodic diversifying selection for three opsin paralogs: L. lineaticorne LW1c, S. tienmushanensis LW1b, and E. regina LW1b (supplementary table S3, Supplementary Material online).

We also found instances of gene duplication in the common ancestor of closely related species. For example, Clades A and B in Fig. 2 each have one or more LW1 opsin from the species L. lineaticorne, Hydropsyche tenuis, Parapsyche elsis, Arctopsyche grandis, and P. conspersa, all but the latter of which belong to the family Hydropsychidae. We tested for episodic diversifying selection using BUSTED (Murrell et al. 2015; Kosakovsky Pond et al. 2020) and found evidence of selection in the sequences in Clade A but not Clade B (Fig. 2, supplementary table S3, Supplementary Material online). Similarly, multiple duplication events occurred in the LW1 opsins in the lineage, leading to N. paramo and A. cinereus, both from the family Leptoceridae (Fig. 2, Clades C to E). Among these duplications, we found evidence of selection in the opsin sequences belonging to only Clade E (supplementary table S3, Supplementary Material online).

Discussion

We searched across 25 caddisfly genome assemblies to determine the number and phylogenetic relationships of opsins in Trichoptera. Our results suggest that caddisfly opsin evolution is likely driven by life-history strategies and ambient light conditions as found in other insect orders (Briscoe et al. 2010; French et al. 2015; Futahashi et al. 2015, Feuda et al. 2016; Lord et al. 2016; Sharkey et al. 2017; Suvorov et al. 2017; Sondhi et al. 2021; Guignard et al. 2022; Mulhair et al. 2023).

We found some incongruencies between the gene tree and the species tree. Given the relatively small number of characters compared with the reference species tree, which was generated from genome-wide data (Frandsen et al. 2024), this is not unexpected and could be due to stochastic error from an undersampling of characters. This is also evidenced by lower bootstrap values in areas of the tree that were incongruent with the species tree (Fig. 2).

The distribution of opsins within Annulipalpia, the fixed-retreat makers, was relatively invariable. Interestingly, their ecological distributions are also less varied; most species inhabit fast-moving streams as larvae and are short lived in the riparian zone as adults. While there were a few LW1 duplications in this suborder, the only loss we observed was the SW opsin gene in P. conspersa (Fig. 3). When searching for opsin genes, we found a sequence highly similar to the SW opsin in this species; however, we excluded it from the dataset due to the presence of stop codons. Given that the genome assembly of P. conspersa was of high quality, we hypothesize that this is likely a true loss (Table 1).

Fig. 3. Opsin counts by species. A bar plot of the number of opsin copies in each caddisfly genome. The bars are colored by opsin type. The species are ordered by the species phylogeny indicated below the bar plot, which is based on a recent study that examined the caddisfly phylogeny in depth (Frandsen et al. 2024). Families are labeled above the bar plot. Suborders are labeled within the phylogeny as follows: (1) Annulipalpia—retreat makers; (2) Integripalpia—cocoon- and tube-case makers; (a) (and Agraylea sexmaculata) basal Integripalpia—cocoon makers; (b) tube-case-making Integripalpia.

The SW opsin gene was also lost in Integripalpia within the basal cocoon-making family, Rhyacophilidae, and within the majority of the clade of tube-case makers (Fig. 3, Clade 2b). Loss of the SW opsin in other groups such as the American cockroach (French et al. 2015) and Neuropteroidea (Lord et al. 2016; Sharkey et al. 2017, 2023) is hypothesized to be associated with the low-light environments in which the ancestors of these insects lived. Caddisflies are primarily crepuscular and thus are most active during low-light conditions, which may be related to the loss of the SW opsin in Rhyacophilidae and in many tube-case-making caddisflies.

Interestingly, in contrast to most species within the tube-case makers, we were able to find the SW opsin in both N. paramo and A. cinereus. The latter was also the only species with two UV opsins and the highest number of LW opsins. Both N. paramo and A. cinereus belong to the family Leptoceridae, a family with known sexual eye dimorphism (Gullefors and Petersson 1993), long adult antennae, and many genera with intricate wing patterns and brightly colored or iridescent wing hairs and scales (Fig. 1; Holzenthal et al. 2007). The high number of opsins and the presence of the SW opsin could be related to the variety in wing coloration and patterns in some species of Leptoceridae. Future work should combine more sampling within this interesting area of the caddisfly phylogeny with gene expression and physiological data to better model the visual system of Trichoptera and to test hypotheses related to color vision and opsin diversity.

To further investigate the role of the duplication events that we observed in the LW and UV opsin groups, we tested paralogs in these areas of the tree for positive selection. We found evidence of episodic diversifying selection in some LW1 opsin paralogs but did not detect evidence of selection in the UV opsin paralogs (supplementary table S3, Supplementary Material online). In each instance when a paralog was found to be under selection, the duplicate paralog was not found to be under selection, possibly suggesting a route to neofunctionalization in those copies undergoing diversifying selection. Recent work in Lepidoptera and Hemiptera has identified instances of family- and species-specific duplications of visual opsins leading to adaptations that extend visual capacity (Spaethe and Briscoe 2004; Frentiu et al. 2007; Sison-Mangus et al. 2008; Briscoe et al. 2010; Feuda et al. 2016; Finkbeiner and Briscoe 2021; Friedrich 2023; Mulhair et al. 2023). Denser taxon sampling in future work can help clarify the evolutionary timing of duplication events and the mechanisms and role of selection that we uncovered.

Here, we conducted the first comprehensive study of visual opsins in Trichoptera. We found that the species with the highest diversity of opsins were derived from the group with known sexual eye dimorphism (Gullefors and Petersson 1993) and, which also contains species of the most colorful and intricately patterned wings, the Leptoceridae. Opsin evolution in caddisflies may also have been driven by life-history strategies and the low-light conditions during which caddisflies are active. The findings of this study provide a basis for future research on the diverse and complex visual systems in Trichoptera.

Materials and Methods

We assessed 25 species of Trichoptera using five newly assembled genomes (details in supplemental note S1, Supplementary Material online) and 20 publicly available genome assemblies (Table 1). Using the 1,000 opsin sequences found by Guignard et al. (2022), we performed a tBLASTn search of opsin sequences against each caddisfly genome, keeping hits with an e-value <10−40. The resulting opsin sequences were extracted from their corresponding genomes using Pyfaidx (Shirley et al. 2015). We filtered redundant hits from multiple queries by extracting the widest window from each contig and then classifying the gene phylogenetically downstream.

We performed gene prediction using AUGUSTUS v3.4.0 (Stanke, et al. 2006), followed by a BLASTp search against the online NCBI databases, maintaining only hits similar to other insect opsins. We then manually checked the annotations in Geneious Prime v2023.0.4 (https://www.geneious.com), using outgroup sequences as a guide, to ensure that the entire gene was correctly annotated (see supplementary note S1, Supplementary Material online, for more details). We included opsin sequences from a variety of insect orders for outgroup comparison: Lepidoptera (Danaus plexippus, Heliconius melpomene, Manduca sexta, Papilio xuthus), Odonata (Sympetrum frequens), Diptera (Drosophila melanogaster), and Coleoptera (Tribolium castaneum), all accessed through GenBank (Benson et al. 2013). We performed additional searches for opsin genes to further verify the absence of the SW and LW2 opsins in many species (see supplementary note S1, Supplementary Material online, for more details). We also provide supplementary tables (supplementary table S2a to c, Supplementary Material online) with information on the completeness of each visual opsin sequence.

Opsin Gene Tree Reconstruction

To determine phylogenetic relationships among opsin sequences, we first aligned the opsin peptide sequences using MAFFT v7.487 (Katoh and Standley 2013) and created a codon alignment with PAL2NAL v14.1 (Suyama et al. 2009). We performed phylogenetic reconstruction on both the CDS and the peptide alignments by first selecting the best substitution model using ModelFinder (Kalyaanamoorthy et al. 2017; Minh et al. 2020) and then performing a maximum likelihood tree search with 1,000 UltraFast bootstrap replicates corrected with the bootstrap nearest neighbor interchange option enabled to guard against the risk of overestimating bootstrap support (-bb 1000 -bnni). We viewed the resulting trees in FigTree v1.4.4 (Rambaut 2018). To highlight differences between the CDS and the peptide trees, we created a face-to-face comparison in R with ggtree v3.2.1 (Yu et al. 2016) and ggplot2 v3.3.5 (Wickham 2016; supplementary fig. S2, Supplementary Material online).

Selection Analysis

To assess the duplications in the LW and UV opsin groups, we created separate codon alignments for both the LW and the UV opsin groups with MAFFT v7.487 (Katoh and Standley 2013) and PAL2NAL v14.1 (Suyama et al. 2009) and then tested for episodic diversifying selection using BUSTED as implemented in HyPhy (Murrell et al. 2015; Kosakovsky Pond et al. 2020). We tested branches of species-specific duplicated opsins individually as well as five deeper duplication events (Fig. 2, Clades A to E) and reported the resulting P-values in supplementary table S3, Supplementary Material online.

Supplementary Material

evae185_Supplementary_Data

Acknowledgments

The authors thank the BYU Office of Research Computing for the compute resources that we used in this study. They thank Ed Wilcox from the BYU DNA Sequencing Center for sequencing the new genomes presented here.

Supplementary Material

Supplementary material is available at Genome Biology and Evolution online.

Funding

The authors thank the BYU College of Life Sciences for providing funding through the College Undergraduate Research Award Program. J.H. acknowledges funding from the Deutsche Forschungsgemeinschaft, Project number 502865717, B.R.-T., S.U.P., R.W.H., and P.B.F. acknowledge funding from the Dirección General de Investigación, Universidad de Las Américas (Ecuador): “Montane freshwater diversity, from taxonomy to functional genomics, an approximation from Trichoptera-Part II” (AMB.BRT.23.02). Funding was also provided by the National Science Foundation, MCB 2217155 to P.B.F.

Authors Contributions

Conceptualization: A.P. and P.B.F. Formal analysis: A.P., J.H., E.R.-G., and P.B.F. Resources: J.H., B.R.-T., R.B.K., and P.B.F. Data curation: A.P., J.H., and E.R.-G. Writing—original draft prep: A.P. Writing—review and editing: all authors. Visualization: A.P. and R.W.H. Supervision: P.B.F. Funding acquisition: A.P., J.H., S.U.P., B.R.-T., and P.B.F.

Data Availability

Genome assemblies are available on GenBank at their respective accession numbers. Gene alignments for the opsins, gff files, tree files, and new genome assemblies are available on FigShare at the following DOI: https://doi.org/10.6084/m9.figshare.24164217.v1.
==== Refs
Literature Cited

Benson  DA, Cavanaugh  M, Clark  K, Karsch-Mizrachi  I, Lipman  DJ, Ostell  J, Sayers  EW. GenBank. Nucleic Acids Res.  2013:41 (D1 ):D36–D42. 10.1093/nar/gks1195.23193287
Briscoe  AD, Bybee  SM, Bernard  GD, Yuan  F, Sison-Mangus  MP, Reed  RD, Warren  AD, Llorente-Bousquets  J, Chiao  C-C. Positive selection of a duplicated UV-sensitive visual pigment coincides with wing pigment evolution in Heliconius butterflies. Proc Natl Acad Sci U S A.  2010:107 (8 ):3628–3633. 10.1073/pnas.0910085107.20133601
Feuda  R, Marletaz  F, Bentley  MA, Holland  PW. Conservation, duplication, and divergence of five opsin genes in insect evolution. Genome Biol Evol.  2016:8 (3 ):579–587. 10.1093/gbe/evw015.26865071
Finkbeiner  SD, Briscoe  AD. True UV color vision in a female butterfly with two UV opsins. J Exp Biol. 2021:224 (18 ):jeb242802. 10.1242/jeb.242802.34587624
Frandsen  PB, Bursell  MG, Taylor  AM, Wilson  SB, Steeneck  A, Stewart  RJ. Exploring the underwater silken architectures of caddisworms: comparative silkomics across two caddisfly suborders. Philos Trans R Soc Lond B Biol Sci. 2019:374 (1784 ):20190206. 10.1098/rstb.2019.0206.31495307
Frandsen  PB, Holzenthal  RW, Espeland  M, Breinholt  J, Thomas Thorpe  JA, Simon  S, Kawahara  AY, Plotkin  D, Hotaling  S, Li  Y, et al  Phylogenomics recovers multiple origins of portable case making in caddisflies (Insecta: Trichoptera), nature's underwater architects. Proc Biol Sci.  2024:291 (2026 ):20240514. 10.1098/rspb.2024.0514.38955232
Frandsen  PB, Hotaling  S, Powell  A, Heckenhauer  J, Kawahara  AY, Baker  RH, Hayashi  CY, Ríos-Touma  B, Holzenthal  R, Pauls  SU, et al  Allelic resolution of insect and spider silk genes reveals hidden genetic diversity. Proc Natl Acad Sci U S A.  2023:120 (18 ):e2221528120. 10.1073/pnas.2221528120.37094147
French  AS, Meisner  S, Liu  H, Weckström  M, Torkkeli  PH. Transcriptome analysis and RNA interference of cockroach phototransduction indicate three opsins and suggest a major role for TRPL channels. Front Physiol.  2015:6 :207. 10.3389/fphys.2015.00207.26257659
Frentiu  FD, Bernard  GD, Sison-Mangus  MP, Brower  AVZ, Briscoe  AD. Gene duplication is an evolutionary mechanism for expanding spectral diversity in the long-wavelength photopigments of butterflies. Mol Biol Evol.  2007:24 (9 ):2016–2028. 10.1093/molbev/msm132.17609538
Friedrich  M . Parallel losses of blue opsin correlate with compensatory neofunctionalization of UV-opsin gene duplicates in aphids and planthoppers. Insects. 2023:14 (9 ):774. 10.3390/insects14090774.37754742
Futahashi  R, Kawahara-Miki  R, Kinoshita  M, Yoshitake  K, Yajima  S, Arikawa  K, Fukatsu  T. Extraordinary diversity of visual opsin genes in dragonflies. Proc Natl Acad Sci U S A.  2015:112 (11 ):E1247–E1256. 10.1073/pnas.1424670112.25713365
Guignard  Q, Allison  JD, Slippers  B. The evolution of insect visual opsin genes with specific consideration of the influence of ocelli and life history traits. BMC Ecol Evol.  2022:22 (1 ):2. 10.1186/s12862-022-01960-8.34996358
Gullefors  B, Petersson  E. Sexual dimorphism in relation to swarming and pair formation patterns in Leptocerid caddisflies (Trichoptera: Leptoceridae). J Insect Behav.  1993:6 (5 ):563–577. 10.1007/BF01048123.
Heckenhauer  J, Frandsen  PB, Gupta  DK, Paule  J, Prost  S, Schell  T, Schneider  JV, Stewart  RJ, Pauls  SU. Annotated draft genomes of two caddisfly species Plectrocnemia conspersa CURTIS and Hydropsyche tenuis NAVAS (Insecta: Trichoptera). Genome Biol Evol.  2019:11 (12 ):3445–3451. 10.1093/gbe/evz264.31774498
Heckenhauer  J, Frandsen  PB, Sproul  JS, Li  Z, Paule  J, Larracuente  AM, Maughan  PJ, Barker  MS, Schneider  JV, Stewart  RJ, et al  Genome size evolution in the diverse insect order Trichoptera. GigaScience. 2022:11 :giac011. 10.1093/gigascience/giac011.35217860
Heckenhauer  J, Stewart  RJ, Ríos-Touma  B, Powell  A, Dorji  T, Frandsen  PB, Pauls  SU. Characterization of the primary structure of the major silk gene, h-fibroin, across caddisfly (Trichoptera) suborders. iScience. 2023:26 (8 ):107253. 10.1016/j.isci.2023.107253.37529107
Holzenthal  R, Blahnik  R, Prather  A, Kjer  K. Order Trichoptera kirby, 1813 (Insecta), caddisflies. Zootaxa. 2007:1668 (1 ):639–698. 10.11646/zootaxa.1668.1.29.
Kalyaanamoorthy  S, Minh  BQ, Wong  TKF, von Haeseler  A, Jermiin  LS. ModelFinder: fast model selection for accurate phylogenetic estimates. Nat Methods.  2017:14 (6 ):587–589. 10.1038/nmeth.4285.28481363
Katoh  K, Standley  DM. MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol Biol Evol.  2013:30 (4 ):772–780. 10.1093/molbev/mst010.23329690
Kawahara  AY, Storer  CG, Markee  A, Heckenhauer  J, Powell  A, Plotkin  D, Hotaling  S, Cleland  TP, Dikow  RB, Dikow  T, et al  Long-read HiFi sequencing correctly assembles repetitive heavy fibroin silk genes in new moth and caddisfly genomes. Gigabyte. 2022:2022 :gigabyte64. 10.46471/gigabyte.64.36824508
Kosakovsky Pond  SL, Poon  AFY, Velazquez  R, Weaver  S, Hepler  NL, Murrell  B, Shank  SD, Magalis  BR, Bouvier  D, Nekrutenko  A, et al  Hyphy 2.5—a customizable platform for evolutionary hypothesis testing using phylogenies. Mol Biol Evol.  2020:37 (1 ):295–299. 10.1093/molbev/msz197.31504749
Lord  NP, Plimpton  RL, Sharkey  CR, Suvorov  A, Lelito  JP, Willardson  BM, Bybee  SM. A cure for the blues: opsin duplication and subfunctionalization for short-wavelength sensitivity in jewel beetles (Coleoptera: Buprestidae). BMC Evol Biol.  2016:16 (1 ):107. 10.1186/s12862-016-0674-4.27193495
Luo  S, Tang  M, Frandsen  PB, Stewart  RJ, Zhou  X. The genome of an underwater architect, the caddisfly Stenopsyche tienmushanensis Hwang (Insecta: Trichoptera). GigaScience. 2018:7 (12 ):giy143. 10.1093/gigascience/giy143.30476205
Minh  BQ, Schmidt  HA, Chernomor  O, Schrempf  D, Woodhams  MD, von Haeseler  A, Lanfear  R. IQ-TREE 2: new models and efficient methods for phylogenetic inference in the genomic era. Mol Biol Evol.  2020:37 (5 ):1530–1534. 10.1093/molbev/msaa015.32011700
Morse  JC, Frandsen  PB, Graf  W, Thomas  JA. Diversity and ecosystem services of Trichoptera. Insects. 2019:10 (5 ):125. 10.3390/insects10050125.31052441
Mulhair  PO, Crowley  L, Boyes  DH, Lewis  OT, Holland  PWH. Opsin gene duplication in Lepidoptera: retrotransposition, sex linkage, and gene expression. Mol Biol Evol.  2023:40 (11 ):msad241. 10.1093/molbev/msad241.37935057
Murrell  B, Weaver  S, Smith  MD, Wertheim  JO, Murrell  S, Aylward  A, Eren  K, Pollner  T, Martin  DP, Smith  DM, et al  Gene-wide identification of episodic selection. Mol Biol Evol.  2015:32 (5 ):1365–1371. 10.1093/molbev/msv035.25701167
Ni  JD, Baik  LS, Holmes  TC, Montell  C. A rhodopsin in the brain functions in circadian photoentrainment in Drosophila. Nature. 2017:545 (7654 ):340–344. 10.1038/nature22325.28489826
Olsen  LK, Heckenhauer  J, Sproul  JS, Dikow  RB, Gonzalez  VL, Kweskin  MP, Taylor  AM, Wilson  SB, Stewart  RJ, Zhou  X, et al  Draft genome assemblies and annotations of Agrypnia vestita Walker, and Hesperophylax magnus Banks reveal substantial repetitive element expansion in tube case-making caddisflies (Insecta: Trichoptera). Genome Biol Evol.  2021:13 (3 ):evab013. 10.1093/gbe/evab013.33501983
Rambaut  A . 2018. FigTree. Version 1.4.4. [Computer software]. http://tree.bio.ed.ac.uk/software/figtree.
Ríos-Touma  B, Holzenthal  RW, Rázuri-Gonzales  E, Heckenhauer  J, Pauls  SU, Storer  CG, Frandsen  PB. De novo genome assembly and annotation of an Andean caddisfly, Atopsyche davidsoni Sykora, 1991, a model for genome research of high-elevation adaptations. Genome Biol Evol.  2021:14 (1 ):evab286. 10.1093/gbe/evab286.
Senthilan  PR, Helfrich-Förster  C. Rhodopsin 7—the unusual rhodopsin in Drosophila. PeerJ. 2016:4 :e2427. 10.7717/peerj.2427.27651995
Sharkey  CR, Blanco  J, Lord  NP, Wardill  TJ. Jewel beetle opsin duplication and divergence is the mechanism for diverse spectral sensitivities. Mol Biol Evol.  2023:40 (2 ):msad023. 10.1093/molbev/msad023.36721951
Sharkey  CR, Fujimoto  MS, Lord  NP, Shin  S, McKenna  DD, Suvorov  A, Martin  GJ, Bybee  SM. Overcoming the loss of blue sensitivity through opsin duplication in the largest animal group, beetles. Sci Rep.  2017:7 (1 ):8. 10.1038/s41598-017-00061-7.28127058
Shichida  Y, Matsuyama  T. Evolution of opsins and phototransduction. Philos Trans R Soc Lond B Biol Sci. 2009:364 (1531 ):2881–2895. 10.1098/rstb.2009.0051.19720651
Shirley  M, Ma  Z, Pedersen  B, Wheelan  S. Efficient “pythonic” access to FASTA files using pyfaidx. PeerJ Prepr. 2015:3 :e970v1. 10.7287/peerj.preprints.970v1.
Sison-Mangus  MP, Briscoe  AD, Zaccardi  G, Knüttel  H, Kelber  A. The lycaenid butterfly Polyommatus icarus uses a duplicated blue opsin to see green. J Exp Biol. 2008:211 (3 ):361–369. 10.1242/jeb.012617.18203991
Sondhi  Y, Ellis  EA, Bybee  SM, Theobald  JC, Kawahara  AY. Light environment drives evolution of color vision genes in butterflies and moths. Commun Biol.  2021:4 (1 ):177. 10.1038/s42003-021-01688-z.33564115
Spaethe  J, Briscoe  AD. Early duplication and functional diversification of the opsin gene family in insects. Mol Biol Evol.  2004:21 (8 ):1583–1594. 10.1093/molbev/msh162.15155799
Stanke  M, Keller  O, Gunduz  I, Hayes  A, Waack  S, Morgenstern  B. AUGUSTUS: ab initio prediction of alternative transcripts. Nucleic Acids Res.  2006:34 (Web Server ):W435–W439. 10.1093/nar/gkl200.16845043
Suvorov  A, Jensen  NO, Sharkey  CR, Fujimoto  MS, Bodily  P, Wightman  HMC, Ogden  TH, Clement  MJ, Bybee  SM. Opsins have evolved under the permanent heterozygote model: insights from phylotranscriptomics of Odonata. Mol Ecol.  2017:26 (5 ):1306–1322. 10.1111/mec.13884.27758014
Suyama  M, Torrents  D, Bork  P. PAL2NAL: robust conversion of protein sequence alignments into the corresponding codon alignments. Nucleic Acids Res.  2009:34 (Web Server ):W609–W612. 10.1093/nar/gkl315.
The Darwin Tree of Life Project Consortium . Sequence locally, think globally: the Darwin Tree of Life Project. Proc Natl Acad Sci U S A.  2022:119 (4 ):e2115642118. 10.1073/pnas.2115642118.35042805
van der Kooi  CJ, Stavenga  DG, Arikawa  K, Belusic  G, Kelber  A. Evolution of insect color vision: from spectral sensitivity to visual ecology. Annu Rev Entomol.  2021:66 (1 ):435–461. 10.1146/annurev-ento-061720-071644.32966103
Wickham  H . Ggplot2: elegant graphics for data analysis. New York: Springer-Verlag; 2016. ISBN 978-3-319-24277-4 . https://ggplot2.tidyverse.org.
Yu  G, Smith  DK, Zhu  H, Guan  Y, Lam  TT-Y. Ggtree: an r package for visualization and annotation of phylogenetic trees with their covariates and other associated data. Methods Ecol Evol.  2016:8 (1 ):28–36. 10.1111/2041-210X.12628.
