Automated identification of aquatic insects: A case study using deep learning and computer vision techniques

Deep learning techniques have recently found application in biodiversity research. Mayflies (Ephemeroptera), stoneflies (Plecoptera) and caddisflies (Trichoptera), often abbreviated as EPT, are frequently used for freshwater biomonitoring due to their large numbers and sensitivity to environmental changes. However, the morphological identification of EPT species is a challenging but fundamental task. Morphological identification of these freshwater insects is therefore not only extremely time-consuming and costly, but also often leads to misjudgments or generates datasets with low taxonomic resolution. Here, we investigated the application of deep learning to increase the efficiency and taxonomic resolution of biomonitoring programs. Our database contains 90 EPT taxa (genus or species level), with the number of images per category ranging from 21 to 300 (16,650 in total). Upon completion of training, a CNN (Convolutional Neural Network) model was created, capable of automatically classifying these taxa into their appropriate taxonomic categories with an accuracy of 98.7 %. Our model achieved a perfect classification rate of 100 % for 68 of the taxa in our dataset. We achieved noteworthy classification accuracy with morphologically closely related taxa within the training data (e.g., species of the genus Baetis, Hydropsyche, Perla). Gradient-weighted Class Activation Mapping (Grad-CAM) visualized the morphological features responsible for the classification of the treated species in the CNN models. Within Ephemeroptera, the head was the most important feature, while the thorax and abdomen were equally important for the classification of Plecoptera taxa. For the order Trichoptera, the head and thorax were almost equally important. Our database is recognized as the most extensive aquatic insect database, notably distinguished by its wealth of included categories (taxa). Our approach can help solve long-standing challenges in biodiversity research and address pressing issues in monitoring programs by saving time in sample identification.

A new species of Rhopalopsole (Plecoptera, Leuctridae) from Hainan Province, China

Background

Hainan Province is an island in the south of China and belongs to the Oriental Region. It has a unique geographical location and superior climatic conditions, providing a good living environment for Leuctridae insects. However, the species richness of the stonefly family Leuctridae in Hainan is low. Two species in total have been recorded, Rhopalopsolebawanglinga Li, Li & Yang, 2023 and Rhopalopsolehainana Li & Yang, 2010.

New information

A new species of Leuctridae (Plecoptera) from Wuzhi Mountains, Hainan Province of south China, Rhopalopsolewuzhishana sp. nov. is described and illustrated. We summarised the Leuctridae in Hainan Province and provide supplemental description and colour plates of Rhopalopsolehainana Li & Yang, 2010.

Utilising public sequence databases to investigate genetic diversity of stoneflies in Medvednica Nature Park

In Medvednica Nature Park, near Croatia's capital Zagreb, urbanisation significantly impacts the fauna. Comprehensive field research has never been conducted in this area, despite the presence of diverse microhabitats and the discovery of several rare species previously unknown in the Croatian fauna. This study provides the Park with first insight into the genetic and morphological diversity of stoneflies, one of the most endangered groups of organisms. Phylogenetic reconstructions and species delineation methods revealed intraspecific haplotype variation in most species (e.g. Brachypteraseticornis, Isoperlagrammatica and Leuctrabraueri), except for Leuctraprima. Additionally, our study has identified isolated populations that merit further in-depth investigation concerning morphology, genetics and ecology.

Sensory pathway in aquatic basal polyneoptera: Antennal sensilla and brain morphology in stoneflies

Aquatic insects represent a great portion of Arthropod diversity and the major fauna in inland waters. The sensory biology and neuroanatomy of these insects are, however, poorly investigated. This research aims to describe the antennal sensilla of nymphs of the stonefly Dinocras cephalotes using scanning electron microscopy and comparing them with the adult sensilla. Besides, central antennal pathways in nymphs and adults are investigated by neuron mass-tracing with tetramethylrhodamine, and their brain structures are visualized with an anti-synapsin antibody. No dramatic changes occur in the antennal sensilla during nymphal development, while antennal sensilla profoundly change from nymphs to adults when switching from an aquatic to an aerial lifestyle. However, similar brain structures are used in nymphs and adults to process diverging sensory information, perceived through different sensilla in water and air. These data provide valuable insights into the evolution of aquatic heterometabolous insects, maintaining a functional sensory system throughout development, including a distinct adaptation of the peripheral olfactory systems during the transition from detection of water-soluble chemicals to volatile compounds in the air. From a conservation biology perspective, the present data contribute to a better knowledge of the biology of stoneflies, which are very important bioindicators in rivers.

Mitochondrial genomes of Nemourinae species (Plecoptera: Nemouridae) and the phylogenetic implications

Currently, the classification system of 2 subfamilies within Nemouridae has been widely accepted. However, monophyly of 2 subfamilies has not been well supported by molecular evidence. To date, only mitogenomes from genus Nemoura of the subfamily Nemourinae were used in previous phylogenetic studies and produced conflicting results with morphological studies. Herein, we analyzed mitogenomes of 3 Nemourinae species to reveal their mitogenomic characteristics and to examine genus-level classification among Nemouridae. In this study, the genome organization of 3 mitogenomes is highly conserved in gene order, nucleotide composition, codon usage, and amino acid composition. In 3 Nemourinae species, there is a high variation in nucleotide diversity among the 13 protein-coding genes (PCGs). The Ka/Ks values for all PCGs were far lower than 1, indicating that these genes were evolving under purifying selection. The phylogenetic analyses highly support Nemurella as the sister group to Ostrocerca. Meanwhile, Nemoura is recovered as the sister group of Malenka; they are grouped with other Amphinemurinae and emerged from a paraphyletic Nemourinae. More molecular data from different taxonomic groups are needed to understand stoneflies phylogeny and evolution.

New insights into the chromosomes of stoneflies: I. Karyotype, C-banding and localization of ribosomal and telomeric DNA markers in Skwalacompacta (McLachlan, 1872) (Polyneoptera, Plecoptera, Perlodidae) from Siberia

This study provides data on chromosome number (2n♂♀=26), sex determination mechanism (XY♂/XX♀), C-banding pattern, distribution of clusters of telomeric TTAGG repeats and 18S ribosomal DNA in the karyotype of the stonefly Skwalacompacta (McLachlan, 1872). For the first time in the history of stoneflies cytogenetics, we provide photos of the chromosomes of the Plecoptera insects. The karyotype of males and females of S.compacta consists of 12 pairs of autosomes. Three pairs of large autosomes and four pairs of medium-sized autosomes are subacrocentric. The remaining pairs of autosomes are small, with unclear morphology. Pericentromeric C-bands were revealed in all autosomes. The sex chromosomes are also subacrocentric. The short arms of X and Y chromosomes are entirely heterochromatic and are rich in ribosomal DNA sequences. In the X chromosome this arm is larger than in the Y chromosome. It is likely that this arm associated with the nucleolar organizer (NOR). Telomeric DNA (TTAGG)n repeats were detected in the terminal regions of all chromosomes.

The First Complete Mitochondrial Genome of Genus Isocapnia (Plecoptera: Capniidae) and Phylogenetic Assignment of Superfamily Nemouroidea

Capniidae are a family of stoneflies, also known as snow flies, who emerge in winter. The phylogeny of Capniidae is widely accepted to be based on morphological analysis. Until now, only five Capniidae mitochondrial genomes have been sequenced so far. In addition, sampling is required to determine an accurate phylogenetic association because the generic classification of this family is still controversial and needs to be investigated further. In this study, the first mitogenome of genus Isocapnia was sequenced with a length of 16,200 bp and contained 37 genes, including a control region, two rRNAs, 22 tRNAs, and 13 PCGs, respectively. Twelve PCGs originated with the common start codon ATN (ATG, ATA, or ATT), while nad5 used GTG. Eleven PCGs had TAN (TAA or TAG) as their last codon; however, cox1 and nad5 had T as their final codon due to a shortened termination codon. All tRNA genes demonstrated the cloverleaf structure, which is distinctive for metazoans excluding the tRNASer1 (AGN) that missed the dihydrouridine arm. A Phylogenetic analysis of the superfamily Nemouroidea was constructed using thirteen PCGs from 32 formerly sequenced Plecoptera species. The Bayesian inference and maximum likelihood phylogeny tree structures derived similar results across the thirteen PCGs. Our findings strongly supported Leuctridae + ((Capniidae + Taeniopterygidae) + (Nemouridae + Notonemouridae)). Ultimately, the best well-supported generic phylogenetic relationship within Capniidae is as follows; (Isocapnia + (Capnia + Zwicknia) + (Apteroperla + Mesocapnia)). These findings will enable us to better understand the evolutionary relationships within the superfamily Nemouroidea and the generic classification and mitogenome structure of the family Capniidae.

Spatial and phylogenetic structure of Alpine stonefly assemblages across seven habitats using DNA-species

Stream ecosystems are spatially heterogeneous, with many different habitat patches distributed within a small area. The influence of this heterogeneity on the biodiversity of benthic insect communities is well documented; however, studies of the role of habitat heterogeneity in species coexistence and assembly remain limited. Here, we investigated how habitat heterogeneity influences spatial structure (beta biodiversity) and phylogenetic structure (evolutionary processes) of benthic stonefly (Plecoptera, Insecta) communities. We sampled 20 sites along two Alpine rivers, including seven habitats in four different reaches (headwaters, meandering, bar-braided floodplain, and lowland spring-fed). We identified 21 morphological species and delineated 52 DNA-species based on sequences from mitochondrial cox1 and nuclear ITS markers. Using DNA-species, we first analysed the patterns of variation in richness, diversity, and assemblage composition by quantifing the contribution of each reach and habitat to the overall DNA-species diversity using an additive partition analysis and distance-based redundancy analysis. Using gene-tree phylogenies, we assessed whether environmental filtering could lead to the co-occurrence of DNA-species using a two-step analysis to detect a phylogenetic signal. All four reaches significantly contributed to DNA-species richness, with the meandering reach having the highest contribution. Habitats had an effect on DNA-species diversity, where glide, riffle and, pool influenced the spatial structure of stonefly assemblage possibly due to the high habitat heterogeneity. Among the habitats, the pool showed significant phylogenetic clustering, suggesting high levels of evolutionary adaptation and strong habitat filtering. This assemblage structure may be caused by long-term stability of the habitat and the similar requirements for co-occurring species. Our study shows the importance of different habitats for the spatial and phylogenetic structure of stonefly assemblage and sheds light on the habitat-specific diversity that may help improve conservation practices.

Additions to the fauna and biology of stoneflies (Plecoptera) in Taizi River Basin, Liaoning, with seven new species records to China

Background

An investigation report of stonefly fauna in Benxi Manchu Autonomous County, Liaoning Province, northeast China (used to be called Manchuria, now includes Liaoning, Jilin, Heilongjiang Provinces and parts of Inner Mongolia, which are adjacent to the Russian Far East and the Korean Peninsula). Materials were studied with field observation in 2018 and 2019.

New information

This paper records five families, nine genera and 14 species of stoneflies from Taizi River, Liaoning Province. Nine species have been recorded for the first time in China and the biology of several common species is also reported for the first time.

New Synonym and New Species Record of Filchneria (Plecoptera: Perlodidae) from China with a Morphological, Phylogenetic and Biogeographic Study on This Genus

The type species of Filchneria Klapálek, 1908, F. mongolica (Klapálek, 1901), is based on a single female collected from Mongolia, but it was considered the same as another species, F. songi from Qinling, China, when the genus Filchneria was proposed. This study narrates the story of these two species, which have been confused for a century. Until now, the distribution of F. mongolica has been confirmed only in Mongolia and Russia, and we recently recorded it for the first time in Inner Mongolia as a new species record in China. Additionally, the genus Sinoperlodes is a junior synonym of Filchneria, as demonstrated by both the morphological and molecular analysis. Phylogenetic analysis based on the subfamily Perlodinae is provided, along with morphological and biogeographic comparisons of Filchneria and its relatives.

New Synonym and New Species Record of Filchneria (Plecoptera: Perlodidae) from China with a Morphological, Phylogenetic and Biogeographic Study on This Genus

The type species of Filchneria Klapálek, 1908, F. mongolica (Klapálek, 1901), is based on a single female collected from Mongolia, but it was considered the same as another species, F. songi from Qinling, China, when the genus Filchneria was proposed. This study narrates the story of these two species, which have been confused for a century. Until now, the distribution of F. mongolica has been confirmed only in Mongolia and Russia, and we recently recorded it for the first time in Inner Mongolia as a new species record in China. Additionally, the genus Sinoperlodes is a junior synonym of Filchneria, as demonstrated by both the morphological and molecular analysis. Phylogenetic analysis based on the subfamily Perlodinae is provided, along with morphological and biogeographic comparisons of Filchneria and its relatives.

The mitochondrial genome analysis of a stonefly, Nemoura longicercia Okamoto, 1922 (Plecoptera: Nemouridae)

To better understand the diversity and phylogeny of Nemouridae, we sequenced and annotated the complete mitochondrial genome (mitogenome) of Nemoura longicercia Okamoto, 1922. The entire mitochondrial genome of N. longicercia was 15,728 bp. It contained 37 representative genes (22 transfer RNA genes (tRNAs), 13 protein-coding genes (PCGs), and two ribosomal RNA genes (rRNAs)) and a control region. In this genome, 23 genes were found on the heavy strand (H-strand) and 14 genes were on the light strand (L-strand). All PCGs began with the initiation codon ATN. Eleven PCGs stopped with the termination codon TAA or TAG, while COII and ND5 genes stopped with incomplete codon T. Bayesian’s inference (BI) and maximum-likelihood (ML) methods generated the identical tree topology across the PCGR dataset (13 PCGs plus two rRNAs). The monophyly of each subfamily was well-supported. This analysis supported the clade N. longicercia plus N. papilla as sister taxon to the clade N. meniscata plus N. nankinensis.

Perlodinellashennongjia sp. nov., a new species of Perlodinella Klapálek (Plecoptera, Perlodidae) from the central area of China

Background

Perlodinella Klapálek, 1912 is a small stonefly genus in the Palearctic areas of China and its biodiversity is underestimated.

New information

This paper reports a new species of Perlodinella, Perlodinellashennongjia sp. nov. in the Dajiuhu National Wetland Park of Shennongjia Forestry District, Hubei Province, central China. The description and illustrations of the new species are provided, based on male adults, female adults and eggs. The new species can be distinguished from its congeners by the characters of male and female genitalia and the egg structure.

A review of Leuctridae (Insecta, Plecoptera) in Wuyi Mountains, China

Background

Wuyi Mountains are located in the northern Oriental Region and the edge of the southern Palaearctic Region. They have a unique geographical location, complex landform and superior climatic conditions, providing a good ecological environment for Leuctridae species. However, due to the damage of some holotypes in the 20^th century, limited drawings and lack of colour figures, it is necessary to reorganise and supplement the preserved Leuctridae specimens from Wuyi Mountains.

New information

In this study, we found that there are twelve species of Leuctridae recorded in Wuyi Mountains, accounting for about 20% of the recorded species of Leuctridae in China. These records include two genera and five new distribution records species: one species of the genus Paraleuctra Hanson, 1941: Paraleuctraorientalis (Chu 1928) and eleven species of the genus Rhopalopsole Klapálek, 1912, including five new distribution records to Wuyi Mountains: Rhopalopsolefengyangshanensis Yang, Shi & Li, 2009; Rhopalopsolesinensis Yang & Yang, 1993; Rhopalopsoleyangdingi Sivec & Harper, 2008; Rhopalopsoleflata Yang & Yang, 1995; Rhopalopsolebasinigra Yang & Yang, 1995. Now a total of twelve species of Leuctridae have been recorded from Wuyi Mountains, Fujian Province of south-eastern China. In this paper, we also provide a key to the male, new images and some notes of these twelve species, except Rhopalopsolerecurvispina (Wu, 1949) and Rhopalopsolespiniplatta (Wu, 1949). We failed to collect these two species and we regard R.recurvispina as a nomen dubium, because there are no distinctive features that can be used to distinguish this species.

DNA barcoding for biodiversity assessment: Croatian stoneflies (Insecta: Plecoptera)

Background

The hemi-metabolous aquatic order Plecoptera (stoneflies) constitutes an indispensable part of terrestrial and aquatic food webs due to their specific life cycle and habitat requirements. Stoneflies are considered one of the most sensitive groups to environmental changes in freshwater ecosystems and anthropogenic changes have caused range contraction of many species. Given the critical threat to stoneflies, the study of their distribution, morphological variability and genetic diversity should be one of the priorities in conservation biology. However, some aspects about stoneflies, especially a fully resolved phylogeny and their patterns of distribution are not well known. A study that includes comprehensive field research and combines morphological and molecular identification of stoneflies has not been conducted in Croatia so far. Thus, the major aim of this study was to regenerate a comprehensive and taxonomically well-curated DNA barcode database for Croatian stoneflies, to highlight the morphological variability obtained for several species and to elucidate results in light of recent taxonomy.

Methods

A morphological examination of adult specimens was made using basic characteristics for distinguishing species: terminalia in males and females, head and pronotum patterns, penial morphology, and egg structures. DNA barcoding was applied to many specimens to help circumscribe known species, identify cryptic or yet undescribed species, and to construct a preliminary phylogeny for Croatian stoneflies.

Results

Sequences (658 bp in length) of 74 morphospecies from all families present in Croatia were recovered from 87% of the analysed specimens (355 of 410), with one partial sequence of 605 bp in length for Capnopsis schilleri balcanica Zwick, 1984. A total of 84% morphological species could be unambiguously identified using COI sequences. Species delineation methods confirmed the existence of five deeply divergent genetic lineages, with monophyletic origin, which also differ morphologically from their congeners and represent distinct entities. BIN (Barcode Index Number) assignment and species delineation methods clustered COI sequences into different numbers of operational taxonomic units (OTUs). ASAP delimited 76 putative species and achieved a maximum match score with morphology (97%). ABGD resulted in 62 and mPTP in 61 OTUs, indicating a more conservative approach. Most BINs were congruent with traditionally recognized species. Deep intraspecific genetic divergences in some clades highlighted the need for taxonomic revision in several species-complexes and species-groups. Research has yielded the first molecular characterization of nine species, with most having restricted distributions and confirmed the existence of several species which had been declared extinct regionally.

Mitochondrial Genomes of the Genus Claassenia (Plecoptera: Perlidae) and Phylogenetic Assignment to Subfamily Perlinae

Mitochondrial genomes of three stoneflies, e.g., Claassenia magna Wu, 1948, Claassenia sp. 2 and Claassenia xucheni Chen, 2019 were sequenced in this study with 15,774, 15,777 and 15,746 bp in length, respectively. Each mitogenome contained 37 genes including 22 tRNAs, two ribosomal RNAs, 13 protein-coding genes (PCGs), and a noncoding control region (CR). In general, standard ATN start and TAN termination codons were evident in the PCGs. Although the dihydrouridine arm was absent in trnSer, the remaining 21 tRNAs displayed the characteristic cloverleaf secondary structure. Stem-loop structures were identified in the CRs of all three mitogenomes, but tandem repeats were only apparent in Claassenia xucheni. The mitogenomes of three Claassenia species were analyzed and compared with mitogenomes in 21 other stoneflies from the Perlidae and three Euholognatha species (Rhopalopsole bulbifera, Capnia zijinshana and Amphinemura longispina) as outgroups. Phylogenetic analyses using maximum likelihood and Bayesian inference. Phylogenetic analysis supported that Claassenia was recovered as the sister group of other Perlinae and Claassenia+Perlinae emerged from the paraphyletic Acroneuriinae. The final results supported that Claassenia was classified into subfamily Perlinae and proposed Claassenia represent a transitional group of the subfamilies Acroneuriinae and Perlinae. This study provided new molecular evidence for exploring the debatable taxonomic position of the genus Claassenia in Perlidae.

The fauna of aquatic invertebrates in the river impacted by wastewaters from the pulp and paper industry (Komi Republic)

BACKGROUND

Invertebrates are important elements of aquatic ecosystems and play a crucial role in the transformation of matter and energy in continental water bodies. Communities of aquatic invertebrates are characterised by high sensitivity to pollution by nutrients and toxic substances and acidification of water bodies; they serve as good bioindicators of the quality of the aquatic environment and impacts on hydroecosystems. All hydrobionts participate in the processes of self-purification of water bodies.The presented dataset provides information on the aquatic invertebrate community of a large northern river. During 2018-2020, we collected data on changes in the quantitative indicators of the development of benthic and planktonic communities, as well as the species diversity of their fauna. The dataset combines information about the occurrence and abundance of benthic and planktonic invertebrates and summarises data of aquatic invertebrate species found in the Vychegda River in the zone of influence from the pulp and paper mill.

NEW INFORMATION

The presented dataset is part of a monitoring programme of the river ecosystems in the production area of Mondi Syktyvkar JSC (the European North-East of Russia, Komi Republic). The dataset describes the structure of benthic invertebrate and plankton communities in the Northern Dvina River Basin. The data on the finding and abundance of large taxa of aquatic invertebrates and species of some groups: Oligochaeta, Cladocera, Copepoda, Rotifera, Ephemeroptera, Plecoptera and Trichoptera are presented. In total, the resource includes 8720 findings of invertebrates, of which 6041 are for zoobenthos organisms and 2679 for zooplankton organisms.

Genomic signatures of parallel alpine adaptation in recently evolved flightless insects

Natural selection along elevational gradients has potential to drive predictable adaptations across distinct lineages, but the extent of such repeated evolution remains poorly studied for many widespread alpine taxa. We present parallel genomic analyses of two recently evolved flightless alpine insect lineages to test for molecular signatures of repeated alpine adaptation. Specifically, we compare low-elevation vs. alpine stonefly ecotypes from parallel stream populations in which flightless upland ecotypes have been independently derived. We map 67,922 polymorphic genetic markers, generated across 176 Zelandoperla fenestrata specimens from two independent alpine stream populations in New Zealand's Rock and Pillar Range, to a newly developed plecopteran reference genome. Genome-wide scans revealed 31 regions with outlier single nucleotide polymorphisms (SNPs) differentiating lowland vs. alpine ecotypes in Lug Creek, and 37 regions with outliers differentiating ecotypes in Six Mile Creek. Of these regions, 13% (8/60) yielded outlier SNPs across both within-stream ecotype comparisons, implying comparable genomic shifts contribute to this repeated alpine adaptation. Candidate genes closely linked to repeated outlier regions include several with documented roles in insect wing-development (e.g., dishevelled), suggesting that they may contribute to repeated alpine wing reduction. Additional candidate genes have been shown to influence insect fecundity (e.g., ovo) and lifespan (e.g., Mrp4), implying that they might contribute to life history differentiation between upland and lowland ecotypes. Additional outlier genes have potential roles in the evolution of reproductive isolation among ecotypes (hedgehog and Desaturase 1). These results demonstrate how replicated outlier tests across independent lineages can potentially contribute to the discovery of genes underpinning repeated adaptation.

Vertical Migration of Adult Plecoptera and Trichoptera above Forested Headwater Streams

Stream insects are essential components of aquatic and terrestrial ecosystem structure and function. Terrestrial stages are important components of terrestrial food webs, and flight-capable individuals are responsible for long-distance dispersal. Horizontal migrations by flying or crawling adults away from stream channels that link insects to riparian food webs and movements across catchment boundaries are well established through empirical research, but studies examining vertical migration of adult stream insects into forest canopies are generally lacking. This study focused on differences in adult Plecoptera and Trichoptera abundance at ground level versus the riparian canopy and differences in abundances among summer and autumn sampling periods to empirically demonstrate use of canopy ecosystems by stream insects. Malaise traps at ground level and canopy traps placed 8 to 10 m above the stream at four sites in the Mosquito Creek watershed (Pennsylvania) were used to examine vertical migration. Larval assemblages were collected and compared to adult assemblage to investigate patterns of local migration in the catchment. We found significantly more stream insects at ground level than in the forest canopy for Trichoptera, Plecoptera, and all individual plecopteran families, but a meaningful number of individuals were found in the riparian canopy. Canopy abundances were similar to abundances captured in adjacent ground-level habitats in other studies. Comparisons of adult and larval abundances among sites, taxa, and stages indicated site- and taxon-specific patterns for vertical movement into riparian canopies. Demonstrating that adult stream insects utilize riparian forest canopies indicates that riparian forest conservation should be prioritized over reforestation and that several potential research questions exist to inform riparian management.

Anthropogenic evolution in an insect wing polymorphism following widespread deforestation

Anthropogenic environmental change can underpin major shifts in natural selective regimes, and can thus alter the evolutionary trajectories of wild populations. However, little is known about the evolutionary impacts of deforestation-one of the most pervasive human-driven changes to terrestrial ecosystems globally. Absence of forest cover (i.e. exposure) has been suggested to play a role in selecting for insect flightlessness in montane ecosystems. Here, we capitalize on human-driven variation in alpine treeline elevation in New Zealand to test whether anthropogenic deforestation has caused shifts in the distributions of flight-capable and flightless phenotypes in a wing-polymorphic lineage of stoneflies from the Zelandoperla fenestrata species complex. Transect sampling revealed sharp transitions from flight-capable to flightless populations with increasing elevation. However, these phenotypic transitions were consistently delineated by the elevation of local treelines, rather than by absolute elevation, providing a novel example of human-driven evolution in response to recent deforestation. The inferred rapid shifts to flightlessness in newly deforested regions have implications for the evolution and conservation of invertebrate biodiversity.

Estimates of resource transfer via winged adult insects from the hyporheic zone in a gravel-bed river

Hyporheic zone (HZ) locates below the riverbed providing habitat for macroinvertebrates from where the winged adult insects (i.e., hyporheic insects, HIs) emerge and bring out aquatic resources to the riparian zone. This study estimated mean daily flux as dry biomass (BM), carbon (C), and nitrogen (N) deriving from the dominant HI species Alloperla ishikariana (Plecoptera, Chloroperlidae) for a 4th-order gravel-bed river during the early-summer to summer periods. We hypothesized that HIs were an important contributor in total aquatic resources to the riparian zone. In 2017 and 2018, we set parallelly (May to August) and perpendicularly (June to October) oriented Malaise traps to catch the lateral and longitudinal directional dispersing winged adults of A. ishikariana, and other Ephemeroptera, Plecoptera, Trichoptera, and Diptera from the river and estimated the directional fluxes of them. We further split the directional fluxes as moving away or back to the channel (for lateral) and from down- to upstream or up- to downstream (for longitudinal). Alloperla ishikariana was similar to other Plecoptera species and differed clearly from Ephemeroptera and Trichoptera in directional characteristics of resources flux, suggesting that the extent and directions of HZ-derived resource transfer depend on taxon-specific flight behaviors of HIs. Contributions of A. ishikariana to the riparian zone in total aquatic C and N transfer seasonally varied and were lower in May (5%-6%) and August (2%-4%) and the highest in July (52%-70%). These conservative estimates largely increased (9% in May) after the supplementary inclusion of Diptera (Chironomidae and Tipulidae), part of which were considered HIs. We demonstrated that HZ could seasonally contribute a significant portion of aquatic resources to the riparian zone and highlighted the potential importance of HZ in nutrient balance in the river-riparian ecosystem.

Female answer specificity to male drumming calls in three closely related species of the stonefly genus Zwicknia (Plecoptera: Capniidae)

This study examines the rate of female answers to conspecific versus heterospecific male vibratory calls in three, closely related stonefly species: Zwicknia bifrons, Z. acuta, and Z. rupprechti. In a previous study those three species were recognized on the basis of their distinct male drumming calls along with differences in genital morphology and genetic divergence. During this study no-choice playback experiments using original male call samples from each species were performed, and the answer rate of females to conspecific and heterospecific signal variants was measured. Mixed effect logistic regression models were used to test if male call species identity had a statistically significant effect on female answer probability. Females answered conspecific male calls with significantly higher probability than heterospecific calls in all the three examined species, suggesting that the divergence of vibrational communication can be an important component of the prezygotic isolation between them. Low, but well detectable responsiveness to heterospecific calls was observable between Z. bifrons and Z. acuta, the two species closest to each other regarding mitochondrial genetic divergence and male call pattern similarity. Thus, our results are most congruent with a tight, gradual coevolution of male calls and female preferences.

Egg structure of five antarctoperlarian stoneflies (Insecta: Plecoptera, Antarctoperlaria)

The egg structures of five antarctoperlarian species - Stenoperla prasina of Eustheniidae; Austroperla cyrene of Austroperlidae; and Zelandobius truncus, Megaleptoperla grandis, and Acroperla trivacuata of Gripopterygidae, were examined in detail, and the groundplan of the egg structure was considered within the representative lineages of Antarctoperlaria and Plecoptera. The flattened egg shape and the circular arrangement of micropyles along the equator are regarded as potential autapomorphies for not only Eustheniidae but also for Eusthenioidea. Austroperlidae has eggs with thin, less-sclerotized chorion, a gelatinous layer on the surface, and micropyles roughly and randomly arranged along the equator. A significant ultrastructural difference between the attachment disc in Gripopterygidae and the anchor plate of arctoperlarian Systellognatha suggests that these structures were independently derived. The thin less-sclerotized chorion represents a groundplan feature in Plecoptera, along with micropyles arranged in a circle, including those circularly arranged along the equator of the egg. On the other hand, in contrast to previous understanding, the sclerotized hard chorion is regarded as a derived feature, having been independently acquired in each of Eustheniidae and Gripopterygidae of Antarctoperlaria and Systellognatha of Arctoperlaria.

Hemocyanin and hexamerins expression in response to hypoxia in stoneflies (Plecoptera, Insecta)

Many freshwater ecosystems worldwide undergo hypoxia events that can trigger physiological, behavioral, and molecular responses in many organisms. Among such molecular responses, the regulation of the hemocyanin (Hc) protein expression which plays a major role in oxygen transportation within aquatic insects remains poorly understood. The stoneflies (Plecoptera) are aquatic insects that possess a functional Hc in the hemolymph similar to crustacean that co-occurs with a nonfunctional Hc protein, hexamerins (Hx). However, the role of both proteins during hypoxia remains undetermined. Here, we evaluated the effect of hypoxia on the expression of Hc and Hx proteins via a comparison between hypoxia and normoxia amino acid sequence variation and protein expression pattern within 23 stonefly species. We induced short-term hypoxia in wild-caught stoneflies species, sequenced the target region of Hc and Hx by complementary DNA synthesis, characterized the protein biochemistry using sodium dodecyl sulfate-polyacrylamide gel electrophoresis, ultrafiltration, and polarographic fluorometric method, and amplified the genome region of the hypoxia-inducible factor (HIF) transcriptional response element that regulated Hc using genome walking library approach. We found a lack of Hc expression in all examined species during hypoxia conditions, despite recognition of the HIF gene region as a possible regulatory factor of Hc, suggesting that compensatory responses as metabolic changes or behavioral tracheal movements to enhance respiratory efficiency could be possible mechanics to compensate for hypoxia. A short Hc-like novel isoform was detected instead in these 23 species, possibly due to either protein degradation or alternative splicing mechanisms, suggesting that the protein could be performing a different function other than oxygen transportation. Hx during hypoxia was expressed and exhibited species-level amino acid changes, highlighting a possible role during hypoxia. Our results demonstrate that hypoxia could enable a similar potential adaptive response of multiple species regarding specific physiological requirements, thereby shedding light on community behavior in stress environments that may help us to improve conservation practices and biomonitoring.

Mountain stoneflies may tolerate warming streams: Evidence from organismal physiology and gene expression

Rapid glacier recession is altering the physical conditions of headwater streams. Stream temperatures are predicted to rise and become increasingly variable, putting entire meltwater-associated biological communities at risk of extinction. Thus, there is a pressing need to understand how thermal stress affects mountain stream insects, particularly where glaciers are likely to vanish on contemporary timescales. In this study, we measured the critical thermal maximum (CT_MAX ) of stonefly nymphs representing multiple species and a range of thermal regimes in the high Rocky Mountains, USA. We then collected RNA-sequencing data to assess how organismal thermal stress translated to the cellular level. Our focal species included the meltwater stonefly, Lednia tumana, which was recently listed under the U.S. Endangered Species Act due to climate-induced habitat loss. For all study species, critical thermal maxima (CT_MAX  > 20°C) far exceeded the stream temperatures mountain stoneflies experience (<10°C). Moreover, while evidence for a cellular stress response was present, we also observed constitutive expression of genes encoding proteins known to underlie thermal stress (i.e., heat shock proteins) even at low temperatures that reflected natural conditions. We show that high-elevation aquatic insects may not be physiologically threatened by short-term exposure to warm temperatures and that longer-term physiological responses or biotic factors (e.g., competition) may better explain their extreme distributions.

Complex and nonlinear climate-driven changes in freshwater insect communities over 42 years

The ongoing biodiversity crisis becomes evident in the widely observed decline in abundance and diversity of species, profound changes in community structure, and shifts in species' phenology. Insects are among the most affected groups, with documented decreases in abundance up to 76% in the last 25-30 years in some terrestrial ecosystems. Identifying the underlying drivers is a major obstacle as most ecosystems are affected by multiple stressors simultaneously and in situ measurements of environmental variables are often missing. In our study, we investigated a headwater stream belonging to the most common stream type in Germany located in a nature reserve with no major anthropogenic impacts except climate change. We used the most comprehensive quantitative long-term data set on aquatic insects available, which includes weekly measurements of species-level insect abundance, daily water temperature and stream discharge as well as measurements of additional physicochemical variables for a 42-year period (1969-2010). Overall, water temperature increased by 1.88 °C and discharge patterns changed significantly. These changes were accompanied by an 81.6% decline in insect abundance, but an increase in richness (+8.5%), Shannon diversity (+22.7%), evenness (+22.4%), and interannual turnover (+34%). Moreover, the community's trophic structure and phenology changed: the duration of emergence increased by 15.2 days, whereas the peak of emergence moved 13.4 days earlier. Additionally, we observed short-term fluctuations (<5 years) in almost all metrics as well as complex and nonlinear responses of the community toward climate change that would have been missed by simply using snapshot data or shorter time series. Our results indicate that climate change has already altered biotic communities severely even in protected areas, where no other interacting stressors (pollution, habitat fragmentation, etc.) are present. This is a striking example of the scientific value of comprehensive long-term data in capturing the complex responses of communities toward climate change.

Aquatic Insects Are Dramatically Underrepresented in Genomic Research

Aquatic insects comprise 10% of all insect diversity, can be found on every continent except Antarctica, and are key components of freshwater ecosystems. However, aquatic insect genome biology lags dramatically behind that of terrestrial insects. If genomic effort was spread evenly, one aquatic insect genome would be sequenced for every ~9 terrestrial insect genomes. Instead, ~24 terrestrial insect genomes have been sequenced for every aquatic insect genome. This discrepancy is even more dramatic if the quality of genomic resources is considered; for instance, while no aquatic insect genome has been assembled to the chromosome level, 29 terrestrial insect genomes spanning four orders have. We argue that a lack of aquatic insect genomes is not due to any underlying difficulty (e.g., small body sizes or unusually large genomes), yet it is severely hampering aquatic insect research at both fundamental and applied scales. By expanding the availability of aquatic insect genomes, we will gain key insight into insect diversification and empower future research for a globally important taxonomic group.

Tolerance of aquifer stoneflies to repeated hypoxia exposure and oxygen dynamics in an alluvial aquifer

Aquatic insects cope with hypoxia and anoxia using a variety of behavioral and physiological responses. Most stoneflies (Plecoptera) occur in highly oxygenated surface waters, but some species live underground in alluvial aquifers containing heterogeneous oxygen concentrations. Aquifer stoneflies appear to be supported by methane-derived food resources, which they may exploit using anoxia-resistant behaviors. We documented dissolved oxygen dynamics and collected stoneflies over 5 years in floodplain wells of the Flathead River, Montana. Hypoxia regularly occurred in two wells, and nymphs of Paraperla frontalis were collected during hypoxic periods. We measured mass-specific metabolic rates (MSMRs) at different oxygen concentrations (12, 8, 6, 4, 2, 0.5 mg l-1, and during recovery) for 111 stonefly nymphs to determine whether aquifer and benthic taxa differed in hypoxia tolerance. Metabolic rates of aquifer taxa were similar across oxygen concentrations spanning 2 to 12 mg l-1 (P>0.437), but the MSMRs of benthic taxa dropped significantly with declining oxygen (P<0.0001; 2.9-times lower at 2 vs. 12 mg l-1). Aquifer taxa tolerated short-term repeated exposure to extreme hypoxia surprisingly well (100% survival), but repeated longer-term (>12 h) exposures resulted in lower survival (38-91%) and lower MSMRs during recovery. Our work suggests that aquifer stoneflies have evolved a remarkable set of behavioral and physiological adaptations that allow them to exploit the unique food resources available in hypoxic zones. These adaptations help to explain how large-bodied consumers might thrive in the underground aquifers of diverse and productive river floodplains.

Remarkable anoxia tolerance by stoneflies from a floodplain aquifer

Alluvial aquifers are key components of river floodplains and biodiversity worldwide, but they contain extreme environmental conditions and have limited sources of carbon for sustaining food webs. Despite this, they support abundant populations of aquifer stoneflies that have large proportions of their biomass carbon derived from methane. Methane is typically produced in freshwater ecosystems in anoxic conditions, while stoneflies (Order: Plecoptera) are thought to require highly oxygenated water. The potential importance of methane-derived food resources raises the possibility that stonefly consumers have evolved anoxia-resistant behaviors and physiologies. Here we tested the anoxic and hypoxic responses of 2,445 stonefly individuals in three aquifer species and nine benthic species. We conducted experimental trials in which we reduced oxygen levels, documented locomotor activity, and measured survival rates. Compared to surface-dwelling benthic relatives, stoneflies from the alluvial aquifer on the Flathead River (Montana) performed better in hypoxic and anoxic conditions. Aquifer species sustained the ability to walk after 4-76 h of anoxia vs. 1 h for benthic species and survived on average three times longer than their benthic counterparts. Aquifer stoneflies also sustained aerobic respiration down to much lower levels of ambient oxygen. We show that aquifer taxa have gene sequences for hemocyanin, an oxygen transport respiratory protein, representing a possible mechanism for surviving low oxygen. This remarkable ability to perform well in low-oxygen conditions is unique within the entire order of stoneflies (Plecoptera) and uncommon in other freshwater invertebrates. These results show that aquifer stoneflies can exploit rich carbon resources available in anoxic zones, which may explain their extraordinarily high abundance in gravel-bed floodplain aquifers. These stoneflies are part of a novel food web contributing biodiversity to river floodplains.

The complete mitochondrial genome of Flavoperla biocellata Chu, 1929 (Plecoptera: Perlidae) and the phylogenetic analyses of Plecoptera

Of the roughly 400 species of Perlidae in the world, most species are widely distributed in the northern hemisphere, but a few can be found in South Africa and South America. There are only five species in the genus Flavoperla of the family Perlidae in China. To gain a better understanding of the architecture and evolution of mitochondrial genome in Flavoperla, the entire mitochondrial genome (mitogenome) of a Chinese Flavoperla biocellata Chu, 1929 from family Perlidae (Insecta: Plecoptera) was sequenced. The 15,805-bp long mitochondrial genome of F. biocellata contained 37 genes, including 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), two ribosomal RNA genes (rRNAs) and a putative control region (CR). The gene arrangement of F. biocellata was identical with that of other stoneflies and with the fly Drosophila yakuba. Most PCGs of F. biocellata used the standard ATN start codons and complete TAN termination codons. Twenty-one of the 22 tRNA genes exhibited cloverleaf secondary structures, but the dihydrouridine (DHU) arm of trnSer (AGN) was completely reduced. Phylogenetic analyses with both Bayesian inference (BI) and maximum likelihood methods (ML) generated similar topology, both supporting the monophyly of all stonefly families and the infraorder Systellognatha. The phylogenetic analysis based on mitochondrial genomic data from 30 stonefly species recovered a well-supported tree resolving higher-level relationships within Plecoptera. The northern hemisphere suborder Arctoperlaria divided into two groups, Euholognatha and Systellognatha. The southern hemisphere suborder Antarctoperlaria formed two clades: Eustheniidae+Diamphipnoidae and Austroperlidae+ Gripopterygidae; consistent with relationships proposed based on morphology. The final relationships within Plecoptera were recovered as (((Perlidae+(Perlodidae+Chloroperlidae))+(Pteronarcyidae+(Peltoperlidae+Styloperlidae))) +(Taeniopterygidae+(Capniidae+(Nemouridae+Notonemouridae))))+ (Gripopterygoidae+Eusthenioidae).

A new winter stonefly (Plecoptera: Capniidae) from Shaanxi Province of China

A new winter stonefly species of the genus Capnia Pictet, 1841, Capnia s.l. bilobata sp. nov. is described and illustrated from Shaanxi Province of northwestern China. The new species is distinguished from all currently known Capnia s.l. by the male and female genitalic characters. The falciform projection on male tergum 9 and the unique-shaped female subgenital plate separate the new taxon from other congeners.

Mitochondrial genomes of the stoneflies Mesonemourametafiligera and Mesonemouratritaenia (Plecoptera, Nemouridae), with a phylogenetic analysis of Nemouroidea

In this study, two new mitochondrial genomes (mitogenomes) of Mesonemourametafiligera and Mesonemouratritaenia from the family Nemouridae (Insecta: Plecoptera) were sequenced. The Mesonemourametafiligera mitogenome was a 15,739 bp circular DNA molecule, which was smaller than that of M.tritaenia (15,778 bp) due to differences in the size of the A+T-rich region. Results show that gene content, gene arrangement, base composition, and codon usage were highly conserved in two species. Ka/Ks ratios analyses of protein-coding genes revealed that the highest and lowest rates were found in ND6 and COI and that all these genes were evolving under purifying selection. All tRNA genes in nemourid mitogenomes had a typical cloverleaf secondary structure, except for tRNA^Ser(AGN) which appeared to lack the dihydrouridine arm. The multiple alignments of nemourid lrRNA and srRNA genes showed that sequences of three species were highly conserved. All the A+T-rich region included tandem repeats regions and stem-loop structures. The phylogenetic analyses using Bayesian inference (BI) and maximum likelihood methods (ML) generated identical results. Amphinemurinae and Nemourinae were sister-groups and the family Nemouridae was placed as sister to Capniidae and Taeniopterygidae.

Ecosystem Services, Global Diversity, and Rate of Stonefly Species Descriptions (Insecta: Plecoptera)

Stoneflies (Insecta: Plecoptera) provide ecosystem services as indicators of water quality, as food for predators, as mediators of energy flow and nutrient cycling, and through cultural services related to recreation and artistic creativity. The Plecoptera Species File (PSF) aggregates stonefly nomenclature, distribution, and literature to help society and scientists understand the value of services stoneflies provide. Using PSF data, we examined global and regional diversity, compared species description rates, and predicted future species description numbers through the year 2100. Through 2018, extant species totaled 3718 with Temperate Asia having the greatest regional diversity at 1178 species. The Perlidae was the most species-rich of the 16 families at 1120 species. The recent global rate of species description was 43.6 species/yr, with Temperate Asia having the highest regional rate at 13.7 species/yr, followed by China and South America adding approximately 9.0 species/yr. We predicted that 1140 ± 130 new species would be described globally by 2050, and 2130 ± 330 by the year 2100, most of the increase occurring in China and South America. We discuss the possibility of reaching these predicted values.

The phylogeny and evolutionary timescale of stoneflies (Insecta: Plecoptera) inferred from mitochondrial genomes

Phylogenetic analysis based on mitochondrial genomic data from 25 stonefly species recovered a well-supported tree resolving higher-level relationships within Plecoptera (stoneflies). The monophyly of both currently recognized suborders was strongly supported, concordant with previous molecular analyses of Plecoptera. The southern hemisphere suborder Antarctoperlaria formed two clades: Eustheniidae + Diamphipnoidae and Austroperlidae + Gripopterygidae; consistent with relationships proposed based on morphology. The largely northern hemisphere suborder Arctoperlaria also divided into two groups, Euholognatha and Systellognatha, each composed of the five families traditionally assigned to each infraorder (the placement Scopuridae by mt genome data remains untested at this time). Within Euholognatha, strong support for the clade Nemouridae + Notonemouridae confirmed the northern origin of the currently southern hemisphere restricted Notonemouridae. Other family level relationships within the Arctoperlaria differ from those recovered by previous morphology and molecular based analyses. A fossil-calibrated divergence estimation suggests the formation of two suborders dates back to the Jurassic (181 Ma), with subsequent diversification of most stonefly families during the Cretaceous. This result confirms the hypothesis that initial divergence between the suborders was driven by the breakup of the supercontinent Pangaea into Laurasia and Gondwanaland (commencing 200 Ma and complete by 150 Ma).

Two new species and one new regional record of Indonemoura from Guangxi, China, with additions to larval characters (Plecoptera, Nemouridae)

Two new species of Nemouridae of the genus Indonemoura Baumann, 1975, Indonemouraquadratasp. n. and Indonemouraquadrispinasp. n., are described from Guangxi Zhuang Autonomous Region of southern China, on the basis of both sexes and larval stage. The affinities towards related species are discussed, together with generic characters of the larvae. Indonemourascalprata (Li & Yang, 2007) is recorded from Guangxi for the first time, and its hitherto unknown female is described.

Potential Effects of Future Climate Changes on Brazilian Cool-Adapted Stoneflies (Insecta: Plecoptera)

The continuous pursuit of welfare and economic development through the exploitation of natural resources by human societies consequently resulted in the ongoing process of climate change. Changes in the distribution of species towards the planet's poles and mountain tops are some of the expected to biological consequences of this process. Here, we assessed the potential effects of future climate change on four cool-adapted Gripopterygidae (Insecta: Plecoptera) species [Gripopteryx garbei Navás 1936, G. cancellata (Pictet 1841), Tupiperla gracilis (Burmeister 1839), and T. tessellata (Brauer 1866)] from Southeastern Brazilian Atlantic forest. As species adapted to cold conditions, in the future scenarios of climate change, we expected these organisms to shrink/change their distributions ranges towards areas with suitable climatic conditions in Southern Brazilian regions, when compared with their predicted distributions in present climatic conditions. We used seven principal components derived from 19 environmental variables from Worldclim database for the present scenario and also seven principal components obtained from 17 different Atmosphere-Ocean Global Circulation Models (AOGCMs), considering the most severe emission scenario for green-house gases to predict the species' distributions. Depending on the climatic scenario considered, there were polewards distribution range changes of the species. Additionally, we also observed an important decrease in the amount of protected modeled range for the species in the future scenarios. Considering that this Brazilian region may become hotter in the future and have its precipitation regime changed, as observed in the severe 2013-2014 drought, we believe these species adapted to high altitudes will be severely threatened in the future.

New species and emendations of Orphella: taxonomic and phylogenetic reassessment of the genus to establish the Orphellales, for stonefly gut fungi with a twist

We consolidate and present data for the sexual stages of five North American species of Orphella, fungal members of trichomycetes previously classified within Harpellales. Three species emendations accommodate the newly recognized characters, including not only the coiled zygospores and accompanying cells but also other morphological traits not provided in the original descriptions for O. avalonensis, O. haysii, and O. hiemalis. We describe three new species, Orphella cataloochensis from both the Smoky Mountains in USA and two provinces in Canada as well as O. pseudoavalonensis and O. pseudohiemalis, both from the Cascade Range, in Oregon, USA. Key morphological features for all known species are summarized and reviewed, with illustrations of some of the North American taxa to update and supplement the literature. The entire suite of morphological characters is discussed, with emphasis on species relationships and hypotheses on possible vicariant origins. We also present a molecular phylogeny based on nuc rDNA 18S and 28S, which supports Orphella as a lineage distinct from Harpellales, and we establish a new order, Orphellales, for it. With the combination of sexual features, now known for 12 of the 14 species of Orphella, and new molecular data, the group is now better characterized, facilitating and hopefully also promoting future studies toward a better understanding of their relationships, origins, and evolutionary history as stonefly gut-dwelling fungi.

Projection of invertebrate populations in the headwater streams of a temperate catchment under a changing climate

Climate change places considerable stress on riverine ecosystems by altering flow regimes and increasing water temperature. This study evaluated how water temperature increases under climate change scenarios will affect stream invertebrates in pristine headwater streams. The studied headwater-stream sites were distributed within a temperate catchment of Japan and had similar hydraulic-geographical conditions, but were subject to varying temperature conditions due to altitudinal differences (100 to 850 m). We adopted eight general circulation models (GCMs) to project air temperature under conservative (RCP2.6), intermediate (RCP4.5), and extreme climate scenarios (RCP8.5) during the near (2031-2050) and far (2081-2100) future. Using the water temperature of headwater streams computed by a distributed hydrological-thermal model as a predictor variable, we projected the population density of stream invertebrates in the future scenarios based on generalized linear models. The mean decrease in the temporally averaged population density of Plecoptera was 61.3% among the GCMs, even under RCP2.6 in the near future, whereas density deteriorated even further (90.7%) under RCP8.5 in the far future. Trichoptera density was also projected to greatly deteriorate under RCP8.5 in the far future. We defined taxa that exhibited temperature-sensitive declines under climate change as cold stenotherms and found that most Plecoptera taxa were cold stenotherms in comparison to other orders. Specifically, the taxonomic families that only distribute in Palearctic realm (e.g., Megarcys ochracea and Scopura longa) were selectively assigned, suggesting that Plecoptera family with its restricted distribution in the Palearctic might be a sensitive indicator of climate change. Plecoptera and Trichoptera populations in the headwaters are expected/anticipated to decrease over the considerable geographical range of the catchment, even under the RCP2.6 in the near future. Given headwater invertebrates play important roles in streams, such as contributing to watershed productivity, our results provide useful information for managing streams at the catchment-level.

Development and ultrastructure of the thickened serosa and serosal cuticle formed beneath the embryo in the stonefly Scopura montana Maruyama, 1987 (Insecta, Plecoptera, Scopuridae)

We aimed to describe the development and ultrastructure of the thickened serosa and serosal cuticle formed beneath the embryo of Plecoptera, using Scopura montana of Scopuridae as a euholognathan representative. Using transmission electron microscopy, we found that the egg membranes were composed of a thick exochorion, a thicker endochorion consisting of two sublayers, and an extremely thin vitelline membrane. The egg membrane construction represents a groundplan feature of the euholognathan egg membranes. The serosa converges beneath the embryo to form a thickened serosa, comprising cells in a radial arrangement, in association with the formation of the amnioserosal fold. The thickened serosa then deposits the thickened serosal cuticle, consisting of four layers differing in fine structure and electron density. After achieving its secretory function, the thickened serosa then disintegrates, and the liberated serosal cells float for a short period in the peripheral region of the egg inside. Collectively, our findings should provide the basis for further characterization of the serosal structures concerned, but we were unable to corroborate previous studies assigning the thickened serosa and serosal cuticle in Plecoptera to the water absorption function.

'Rolling' stoneflies (Insecta: Plecoptera) from mid-Cretaceous Burmese amber

This contribution describes seven new species of fossil stoneflies from Cretaceous Burmese amber, all of which are dedicated to present and past members of the Rolling Stones. Two species-Petroperla mickjaggeri gen. nov. sp. nov. and Lapisperla keithrichardsi gen. nov. sp. nov.-are placed in a new family Petroperlidae within the stemline of Systellognatha. The first Cretaceous larval specimen of Acroneuriinae, Electroneuria ronwoodi gen. nov. sp. nov., is also described along with another four new species that are placed within the Acroneuriinae genus Largusoperla Chen et al., 2018: Largusoperla charliewattsi sp. nov., Largusoperla brianjonesi sp. nov., Largusoperla micktaylori sp. nov., and Largusoperla billwymani sp. nov. Additional specimens of Acroneuriinae are described without formal assignment to new species due to insufficient preservation. Implications for stonefly phylogeny and palaeobiogeography are discussed.

Complete Mitochondrial Genome of Suwallia teleckojensis (Plecoptera: Chloroperlidae) and Implications for the Higher Phylogeny of Stoneflies

Stoneflies comprise an ancient group of insects, but the phylogenetic position of Plecoptera and phylogenetic relations within Plecoptera have long been controversial, and more molecular data is required to reconstruct precise phylogeny. Herein, we present the complete mitogenome of a stonefly, Suwallia teleckojensis, which is 16146 bp in length and consists of 13 protein-coding genes (PCGs), 2 ribosomal RNAs (rRNAs), 22 transfer RNAs (tRNAs) and a control region (CR). Most PCGs initiate with the standard start codon ATN. However, ND5 and ND1 started with GTG and TTG. Typical termination codons TAA and TAG were found in eleven PCGs, and the remaining two PCGs (COII and ND5) have incomplete termination codons. All transfer RNA genes (tRNAs) have the classic cloverleaf secondary structures, with the exception of tRNASer(AGN), which lacks the dihydrouridine (DHU) arm. Secondary structures of the two ribosomal RNAs were shown referring to previous models. A large tandem repeat region, two potential stem-loop (SL) structures, Poly N structure (2 poly-A, 1 poly-T and 1 poly-C), and four conserved sequence blocks (CSBs) were detected in the control region. Finally, both maximum likelihood (ML) and Bayesian inference (BI) analyses suggested that the Capniidae was monophyletic, and the other five stonefly families form a monophyletic group. In this study, S. teleckojensis was closely related to Sweltsa longistyla, and Chloroperlidae and Perlidae were herein supported to be a sister group.

Molecular phylogeny of Systellognatha (Plecoptera: Arctoperlaria) inferred from mitochondrial genome sequences

The infraorder Systellognatha is the most species-rich clade in the insect order Plecoptera and includes six families in two superfamilies: Pteronarcyoidea (Pteronarcyidae, Peltoperlidae, and Styloperlidae) and Perloidea (Perlidae, Perlodidae, and Chloroperlidae). To resolve the debatable phylogeny of Systellognatha, we carried out the first mitochondrial phylogenetic analysis covering all the six families, including three newly sequenced mitogenomes from two families (Perlodidae and Peltoperlidae) and 15 published mitogenomes. The three newly reported mitogenomes share conserved mitogenomic features with other sequenced stoneflies. For phylogenetic analyses, we assembled five datasets with two inference methods to assess their influence on topology and nodal support within Systellognatha. The results indicated that inclusion of the third codon positions of PCGs, exclusion of rRNA genes, the use of nucleotide datasets and Bayesian inference could improve the phylogenetic reconstruction of Systellognatha. The monophyly of Perloidea was supported in the mitochondrial phylogeny, but Pteronarcyoidea was recovered as paraphyletic and remained controversial. In this mitochondrial phylogenetic study, the relationships within Systellognatha were recovered as (((Perlidae + (Perlodidae + Chloroperlidae)) + (Pteronarcyidae + Styloperlidae)) + Peltoperlidae).

The first two mitochondrial genomes from Taeniopterygidae (Insecta: Plecoptera): Structural features and phylogenetic implications

The complete mitochondrial genomes (mitogenomes) of Taeniopteryx ugola and Doddsia occidentalis (Plecoptera: Taeniopterygidae) were firstly sequenced from the family Taeniopterygidae. The 15,353-bp long mitogenome of T. ugola and the 16,020-bp long mitogenome of D. occidentalis each contained 37 genes including 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), two ribosomal RNA genes (rRNAs) and a control region (CR). The mitochondrial gene arrangement of the two taeniopterygids and other stoneflies was identical with the putative ancestral mitogenome of Drosophila yakuba. Most PCGs used standard ATN start codons and TAN termination codons. Twenty-one of the 22 tRNAs in each mitogenome could fold into the cloverleaf secondary structures, while the dihydrouridine (DHU) arm of trnSer (AGN) was reduced or absent. Stem-loop (SL) structures, poly-T stretch, poly-[AT]_n stretch and tandem repeats were found in the CRs of the two mitogenomes. The phylogenetic analyses using Bayesian inference (BI) and maximum likelihood methods (ML) generated identical results, both supporting the monophyly of all stonefly families and the two infraorders, Systellognatha and Euholognatha. Taeniopterygidae was grouped with another two families from Euholognatha. The relationships within Plecoptera were recovered as (((Perlidae+Peltoperlidae)+((Pteronarcyidae+Chloroperlidae)+Styloperlidae))+((Capniidae+Taeniopterygidae)+Nemouridae))+Gripopterygidae.

Does wing reduction influence the relationship between altitude and insect body size? A case study using New Zealand's diverse stonefly fauna

Researchers have long been intrigued by evolutionary processes that explain biological diversity. Numerous studies have reported strong associations between animal body size and altitude, but insect analyses have often yielded equivocal results. Here, we analyze a collection database of New Zealand's diverse endemic stonefly fauna (106 species across 21 genera) to test for relationships between altitude and plecopteran body size. This insect assemblage includes a variety of wing-reduced (26 spp) and fully winged (80 spp) taxa and covers a broad range of altitudes (0-2,000 m). We detected significant relationships between altitude and body size for wing-reduced, but not fully winged, stonefly taxa. These results suggest that, while the maintenance of flight apparatus might place a constraint on body size in some fully winged species, the loss of flight may free insects from this evolutionary constraint. We suggest that rapid switches in insect dispersal ability may facilitate rapid evolutionary shifts across a number of biological attributes and may explain the inconsistent results from previous macroecological analyses of insect assemblages.

Differences in protein expression among five species of stream stonefly (Plecoptera) along a latitudinal gradient in Japan

Proteome variation among natural populations along an environmental gradient may provide insights into how the biological functions of species are related to their local adaptation. We investigated protein expression in five stream stonefly species from four geographic regions along a latitudinal gradient in Japan with varying climatic conditions. The extracted proteins were separated by two-dimensional gel electrophoresis and identified by matrix-assisted laser desorption/ionization of time-of-flight (MALDI TOF/TOF), yielding 446 proteins. Low interspecies variation in the proteome profiles was observed among five species within geographical regions, presumably due to the co-occurring species sharing the environments. However, large spatial variations in protein expression were found among four geographic regions, suggesting strong regulation of protein expression in heterogeneous environments, where the spatial variations were positively correlated with water temperature. We identified 21 unique proteins expressed specifically in a geographical region and six common proteins expressed throughout all regions. In warmer regions, metabolic proteins were upregulated, whereas proteins related to cold stress, the photoperiod, and mating were downregulated. Oxygen-related and energy-production proteins were upregulated in colder regions with higher altitudes. Thus, our proteomic approach is useful for identifying and understanding important biological functions related to local adaptations by populations of stoneflies.

Application of Aquatic Insects (Ephemeroptera, Plecoptera And Trichoptera) In Water Quality Assessment of Malaysian Headwater

The Ephemeroptera, Plecoptera and Trichoptera (EPT) community structure and the specific sensitivity of certain EPT genera were found to be influenced by water parameters in the rivers of Gunung Jerai Forest Reserve (GJFR) in the north of peninsular Malaysia. The scores of EPT taxa richness of >10 in all rivers indicated all rivers’ habitats were non-impacted, having good water quality coinciding with Class I and Class II of Malaysian water quality index (WQI) classification of potable water. The abundance of EPT was very high in Teroi River (9,661 individuals) but diversity was lower (22 genera) than Tupah River which was highly diverse (28 genera) but lower in abundance (4,263 individuals). The lowest abundance and moderate diversity was recorded from Batu Hampar River (25 genera). Baetis spp. and Thalerosphyrus spp., Neoperla spp. and Cheumatopsyche spp. were the most common genera found. Classification for all rivers using EPT taxa Richness Index and WQI gave different category of water quality, respectively. The WQI classified Tupah and Batu Hampar rivers into Class II and Teroi River (Class I) was two classes above the classification of the EPT taxa Richness Index.

Climate-induced glacier and snow loss imperils alpine stream insects

Climate warming is causing rapid loss of glaciers and snowpack in mountainous regions worldwide. These changes are predicted to negatively impact the habitats of many range-restricted species, particularly endemic, mountaintop species dependent on the unique thermal and hydrologic conditions found only in glacier-fed and snow melt-driven alpine streams. Although progress has been made, existing understanding of the status, distribution, and ecology of alpine aquatic species, particularly in North America, is lacking, thereby hindering conservation and management programs. Two aquatic insects - the meltwater stonefly (Lednia tumana) and the glacier stonefly (Zapada glacier) - were recently proposed for listing under the U.S. Endangered Species Act due to climate-change-induced habitat loss. Using a large dataset (272 streams, 482 total sites) with high-resolution climate and habitat information, we describe the distribution, status, and key environmental features that limit L. tumana and Z. glacier across the northern Rocky Mountains. Lednia tumana was detected in 113 streams (175 sites) within Glacier National Park (GNP) and surrounding areas. The probability of L. tumana occurrence increased with cold stream temperatures and close proximity to glaciers and permanent snowfields. Similarly, densities of L. tumana declined with increasing distance from stream source. Zapada glacier was only detected in 10 streams (24 sites), six in GNP and four in mountain ranges up to ~600 km southwest. Our results show that both L. tumana and Z. glacier inhabit an extremely narrow distribution, restricted to short sections of cold, alpine streams often below glaciers predicted to disappear over the next two decades. Climate warming-induced glacier and snow loss clearly imperils the persistence of L. tumana and Z. glacier throughout their ranges, highlighting the role of mountaintop aquatic invertebrates as sentinels of climate change in mid-latitude regions.

First Mitochondrial Genome from Nemouridae (Plecoptera) Reveals Novel Features of the Elongated Control Region and Phylogenetic Implications

The complete mitochondrial genome (mitogenome) of Nemoura nankinensis (Plecoptera: Nemouridae) was sequenced as the first reported mitogenome from the family Nemouridae. The N. nankinensis mitogenome was the longest (16,602 bp) among reported plecopteran mitogenomes, and it contains 37 genes including 13 protein-coding genes (PCGs), 22 transfer RNA (tRNA) genes and two ribosomal RNA (rRNA) genes. Most PCGs used standard ATN as start codons, and TAN as termination codons. All tRNA genes of N. nankinensis could fold into the cloverleaf secondary structures except for trnSer (AGN), whose dihydrouridine (DHU) arm was reduced to a small loop. There was also a large non-coding region (control region, CR) in the N. nankinensis mitogenome. The 1751 bp CR was the longest and had the highest A+T content (81.8%) among stoneflies. A large tandem repeat region, five potential stem-loop (SL) structures, four tRNA-like structures and four conserved sequence blocks (CSBs) were detected in the elongated CR. The presence of these tRNA-like structures in the CR has never been reported in other plecopteran mitogenomes. These novel features of the elongated CR in N. nankinensis may have functions associated with the process of replication and transcription. Finally, phylogenetic reconstruction suggested that Nemouridae was the sister-group of Capniidae.

Antennal sensilla of the stonefly Dinocras cephalotes (Plecoptera: Perlidae)

Plecoptera, one of the most primitive groups of Neoptera, are important aquatic insects usually employed as bioindicators of high water quality. Notwithstanding the well-developed antennae of the adult, its sensory abilities are so far not well known. The present paper describes at ultrastructural level under scanning and transmission electron microscopy the antennal sensilla of the adult stonefly Dinocras cephalotes (Plecoptera, Perlidae). Adult males and females show a filiform antenna constituted of a scape, a pedicel and a flagellum composed of very numerous segments with no clear sexual dimorphism in the number and distribution of the antennal sensilla. The most represented sensilla are sensilla trichodea, with different length, whose internal structure reveal their mechanosensory function, sensilla chaetica, with an apical pore, with an internal structure revealing a typical gustatory function, porous pegs representing single-walled olfactory sensilla, digitated pegs with hollow cuticular spoke channels representing double-walled olfactory sensilla, pegs in pits for which we hypothesize a thermo-hygrosensory function. The diversity of described sensilla is discussed in relation to known biological aspects of the studied species. This opens new perspectives in the study of the behavior of these aquatic insects during their adult stage.

A new species of Neoperla from China, with a redescription of the female of N. mnong Stark, 1987 (Plecoptera, Perlidae)

A new species of the Neoperla clymene group (Plecoptera, Perlidae), Neoperla chebalinga sp. n. from Guangdong Province of southern China is described, illustrated, and compared with related taxa. The new species is characterized by the slender aedeagal tube, strongly sclerotized dorsally, and weakly sclerotized ventrally with an upcurved, medial, finger-like membranous lobe. Additionally the aedeagal sac gradually tapers to a blunt apex with a dorsoapical patch of spines. A supplementary description of the female of Neoperla mnong Stark, 1987 from Guangdong Province, China is also given.