Molecular association and morphological characterisation of Himalopsyche larval types (Trichoptera, Rhyacophilidae)

Anna E. Hjalmarsson; Wolfram Graf; Sonja C. Jähnig; Simon Vitecek; Steffen U. Pauls

doi:10.3897/zookeys.773.24319

Abstract

Himalopsyche Banks, 1940 (Trichoptera, Rhyacophilidae) is a genus of caddisflies inhabiting mountain and alpine environments in Central and East Asia and the Nearctic. Of 53 known species, only five species have been described previously in the aquatic larval stage. We perform life stage association using three strategies (GMYC, PTP, and reciprocal monophyly) based on fragments of two molecular markers: the nuclear CAD, and the mitochondrial COI gene. A total of 525 individuals from across the range of Himalopsyche (Himalayas, Hengduan Shan, Tian Shan, South East Asia, Japan, and western North America) was analysed and 32 operational taxonomic units (OTUs) in our dataset delimited. Four distinct larval types of Himalopsyche are uncovered, and these are defined as the phryganea type, japonica type, tibetana type, and gigantea type and a comparative morphological characterisation of the larval types is presented. The larval types differ in a number of traits, most prominently in their gill configuration, as well as in other features such as setal configuration of the pronotum and presence/absence of accessory hooks of the anal prolegs.

Keywords

Caddisfly, GMYC, Hengduan Mountains, Himalaya, life stage association, PTP

Introduction

With ~15,000 described and around 50,000 presumed species, caddisflies are one of the larger insect orders and the largest primary aquatic insect order (de Moor and Ivanov 2008). Trichoptera have merolimnic life histories and the larvae have long been recognised as important ecological indicators (Resh and Unzicker 1975). They have several remarkable ecological traits, among others their diverse case-building behaviour or their ability to invade all types of aquatic habitats across the globe (except Antarctica; Holzenthal et al. 2007, Wiggins 1996). While renowned for their intricate cases, larvae of some caddis families roam freely, and only build pupal retreats. Among those, the family Rhyacophilidae Stephens, 1836 is particularly noteworthy for their diversity and ecological differentiation. This family entails the genera Himalopsyche Banks, 1940, Fansipangana Mey, 1996, Philocrena Lepneva, 1956, Phoupanpsyche Malicky, 2008, and Rhyacophila Pictet, 1834.

Species of the genus Himalopsyche are a particularly interesting group of caddisflies. They primarily inhabit mountain and alpine environments in Central and East Asia, although the genus also radiated into the Nearctic where it is represented by a single species, H. phryganea (Ross, 1941). Himalopsyche larvae mostly inhabit highly turbulent, fast-flowing streams, where they live as ferocious predators.

Including recent species descriptions, the genus Himalopsyche currently comprises 53 known species, adding to the status of the last major treatment of Himalopsyche by Schmid and Botosaneanu (1966). While our knowledge of the adult taxonomy of Himalopsyche is comparatively good, we know very little about the larval taxonomy and ecology of individual species. To date, five species have been described in the larval stage: H. japonica (Morton, 1900), H. phryganea, H. gigantea (Martynov, 1914), H. tibetana (Martynov, 1930), and H. acharai Malicky & Chantaramongkol, 1989 (Flint 1961, Graf and Sharma 1998, Lepneva 1945, 1970, Saito 1965, Tanida 1985, Thamsenanupap et al. 2005), as well as a larva corresponding to a hitherto unknown species, referred to as H. ‘larva hoplura’ (Lepneva 1945, 1970). Three distinct types of Himalopsyche larvae have been differentiated, most prominently based on their gill configuration, but also on differences in setal configuration of the pronotum and anal sclerites. In their comparative study, Graf and Sharma (1998) defined two of these larval types, Type A and Type B, which differ distinctly from the previously known larvae of H. phryganea and H. japonica. Type B could be assigned to H. tibetana. The species identity of larvae assigned to Type A could not be clarified before now, but this type shows similarities to larvae of H. gigantea. The larva of H. acharai was described by Thamsenanupap et al. (2005) and was compared to H. japonica and H. phryganea, but not with Type A and B sensuGraf and Sharma (1998). Thus, all Himalopsyche species known in the larval stage have never been compared and characterised simultaneously before.

Caddisflies are good biological indicators (e.g., Rainbow et al. 2012) and are essential elements in many standardised assessment systems, especially in North America, Europe, and Australia (e.g., Schmidt-Kloiber et al. 2008). However, the practical use of Trichoptera larvae as biological indicators is limited whenever taxonomic knowledge is poor. Higher-level taxonomic resolution (e.g., genus or family level) can often mask the variability of environment/species interactions that act at the species level and make these taxa valuable indicator species (Resh and Unzicker 1975, Verdonschot 2006, Ruiter et al. 2013).

Molecular data have proven successful in facilitating the association of larvae with adults in Trichoptera (e.g., Graf et al. 2005, Ruiter et al. 2013, Zhou et al. 2007). Since the earliest studies of this kind (Graf et al. 2005, Zhou et al. 2007) the mitochondrial cytochrome c oxidase (COI) and nuclear rDNA 28S genes have been the markers of choice for life stage associations. Zhou et al. (2007) suggested a protocol for life stage association based on COI (439 bp) and 28S (~430 bp) and emphasised that it is important to use more than a single genetic marker for life stage association. There are several reasons for this. First, data from different genes provide the possibility to cross-validate the results from each gene, and to identify potential contamination issues. Second, gene trees are likely to differ from the species phylogeny among close relatives because of incomplete lineage sorting (Funk and Omland 2003). Using more than one gene therefore increases the opportunities to cross-validate results from one gene with the other gene(s).

In this study, we associate larvae with adults based on molecular data from a nuclear and a mitochondrial gene, CAD and COI, respectively, and employ three different methods for life stage association to generate a consensus result. We then present the first comparative morphological characterisation of all known larval types of Himalopsyche.

Materials and methods

Sampling

Larval specimens used in this study were mostly collected in Nepal and China in 2011–2013; adult material was largely obtained from research collections (Suppl. material 1). In Nepal, samples were collected in 2012 and 2013 from the Langtang, Indrawati, and Arun catchments (Tachamo Shah et al. 2015) and in China from the Lancang and Jinsha catchments in 2011 and 2013 (Figure 1). Larvae were collected from streams using hand nets, kick nets and hand-picking. Adults were collected using light traps, and occasionally with hand nets and a Malaise trap. Light traps were either set up at dusk, using a white sheet and a UV lamp (active light trapping), or left over night with the lamp placed over a small pan filled with soap water (passive light trapping; Blahnik and Holzenthal 2004, Malicky 2004). Lamps used included 12V 15W Bioform blacklights, F8W/T5/BL350 Sylvania blacklights, and 22W circline BioQuip blacklights. All specimens were stored directly in 95% ethanol in the field, and refrigerated in the lab. Museum specimens of adults were borrowed for DNA sequencing from the private collection of Hans Malicky (Lunz am See, Austria), and from the Museum für Naturkunde (Berlin, Germany).

Figure 1.

Map showing sampling localities in A North America, and B East Asia. Colours/symbols indicate larval types of the OTUs. Elevation data from Jarvis et al. (2008), stream data from Andreadis et al. (2013).

Molecular data

For this study, we used partial sequence data of the single copy nuclear marker CAD and the mitochondrial COI. CAD has proved useful for insect phylogenetics (Bukontaite et al. 2014, Ekrem et al. 2010, Johanson and Malm 2010, Klopfstein et al. 2013, Moulton and Wiegmann 2004, Sharanowski et al. 2011, Wild and Maddison 2008) and in species delimitation (Foster et al. 2013, Song and Ahn 2014, Vitecek et al. 2017). For COI, we targeted the 658 bp “standard barcode region” fragment of COI that has been used extensively for species identification, but also for life stage association, population genetics, and phylogeny in insects (e.g., Hebert et al. 2004, Hjalmarsson et al. 2015, Ruiter et al. 2013, Zhou et al. 2007, 2016).

DNA was extracted from legs using one of the following methods: HotShot protocol (Montero-Pau et al. 2008; mainly used for larvae), Qiagen Dneasy Blood & Tissue Kit (Qiagen, mainly used for fresh adult material), or QIAamp DNA Micro Kit (Qiagen, mainly used for museum material). We used a combination of previously published and newly developed primers, specific to Himalopsyche. The primers for CAD were: 743nF-ino & 1028r-ino (850 bp; Johanson and Malm 2010) and C1Fb & C7Ra (758 bp; this work, Table 1). The COI primers were: HCO1490 &LCO2198 (658 bp; Folmer et al. 1994), and B1Fa & B3Ra (367 bp; this work, Table 1). Polymerase Chain Reaction (PCR) was performed using PeqGOLD Hot Start Taq Polymerase kits (PeqLab VWR) in standard reactions, for some protocols with the addition of BSA (Table 2). Sanger sequencing of PCR products was performed on a 3730XL DNA Analyzer (Applied Biosystems) at the Senckenberg Biodiversity Climate Research Centre Laboratory Centre. Sequences were assembled and edited in Geneious 7.0.6 (Kearse et al. 2012). Ambiguities were coded using IUPAC codes. Multiple sequence alignments of CAD and COI were made using the ClustalW algorithm (Thompson et al. 1994) as implemented in Geneious, and were checked for stop codons.

Table 1.

Download as

CSV

XLSX

Primers used for PCR and sequencing. Fragment lengths refer to the primer pairs 743nF-ino & 1028r-ino, C1Fb & C7Ra, HCO1490 & LCO2198, and B1Fa & B3Ra.

Gene	Primer	Sequence	Tm (°C)	Fragment length (bp)	Reference
CAD	743nF-ino	5’-GGIGTIACIACIGCITGYTTYGARCC-3’	52.4	850	Johanson and Malm 2010
CAD	1028r-ino	5’–TTRTTIGGIARYTGICCICCCAT–3’	42.1	850	Johanson and Malm 2010
CAD_internal	C1Fb	5’–TGYGTTGTRAAGATTCCGAG-3’	51.8	736	this work
CAD_internal	C7Ra	5’–TGTCCATTACAACCTCGAATG-3’	62.3	736	this work
COI	HCO1490	5’-GGTCAACAAATCATAAAGATATTGG-3’	53.2	658	Folmer et al. 1994
COI	LCO2198	5’-TAAACTTCAGGGTGACCAAAAAATCA-3’	51.0	658	Folmer et al. 1994
COI _internal	B1Fa	5’-ATTGCDACWGATCAWACAAA-3’	54.9	367	this work
COI _internal	B3Ra	5’-AAYGTARTWGTWACWGCTCA-3’	47.2	367	this work

Table 2.

Download as

CSV

XLSX

Protocols for 10 μL PCR reactions, using VWR peqGOLD Hot Taq DNA Polymerase kits. BSA = Bovine serum albumin. All numbers are given in μL.

Gene	Primer pair	Buffer S	dNTPS 2mM	BSA 20 mg/ml	Forward primer 10μM	Reverse primer 10μM	Taq	DNA	H₂O	Cycler
CAD	1028r-ino, 743nF-ino	1	1	0	0.25	0.25	0.1	1	6.4	5’ 95°C, 35x (45’’ 95°C, 45” 55°C, 60” 72°C) 5’ 72°C
CAD	C1Fb, C7Ra	1	1	0	0.25	0.25	0.1	1	6.4	5’ 95°C, 35x (30’’ 95°C, 30” 50°C, 45” 72°C) 5’ 72°C
COI	HCO1490, LCO2198	1	1	0	0.25	0.25	0.1	1	6.4	5’ 95°C, 5 x (30’’ 95°C, 1’ 44°C, 1’ 72°C), 15x (30’’ 95°C, 30’’ 48°C, 1’ 72°C), 20 x (30’’ 95°C, 30’’ 50°C, 1’ + (10’’ * n) 72°C), 5’ 72°C
COI	B1Fa, B3Ra	1	1	0.4	0.25	0.25	0.1	1	6.0	5’ 95°C, 35x (30’’ 95°C, 30” 45°C, 45” 72°C) 5’ 72°C

Figure 2.

Summary of GMYC results from the CAD and COI, indicating the overlap in results from the two genes. Colours indicate larval types, inferred from this study and the literature.

Life stage association

Adult males were identified to species based on morphology. The dataset included 38 adult species based on males, including four putative new species (Hjalmarsson submitted, Kuranishi et al. unpublished data). We used and compared the results of three phylogenetic association criteria: Poisson Tree Process (PTP; Zhang et al. 2013), General Mixed Yule Coalescent method (GMYC; Fujisawa and Barraclough 2013, Pons et al. 2006), and reciprocal monophyly. In total, we had five life stage association criteria (PTP and GYMC for each gene, and reciprocal monophyly). Life stage association was considered successful if at least three criteria were fulfilled. The life stage association criteria were defined as follows. PTP: larvae are conspecific with an adult if they form a PTP cluster containing only one and the same adult species. GMYC: larvae are conspecific with adults if they form a GMYC cluster containing only one and the same adult species. Reciprocal monophyly: species are considered reciprocally monophyletic if both genes return a monophylum containing one and the same adult species, and the same larval specimens. Reciprocal monophyly could only be tested for species with data from both genes and with >1 specimen per gene. We considered nodes as supported if they had posterior probability of at least 95%. We refer to groups of unresolved and paraphyletic species as ‘species complexes’.

Gene trees of all available specimens were reconstructed in MrBayes v3.2.6 (Ronquist et al. 2012). Alignments were partitioned per codon position with independent rates among partitions. Nucleotide substitution models were determined in PartitionFinder v2.1.1. (Lanfear et al. 2016; Table 3). Runs were generated for 10 to 50 million generations and were checked for convergence in Tracer v.1.6 (Rambaut et al. 2014). Fully resolved 50% majority rule consensus trees were generated using the ‘sumt’ command, with 25% burn-in. PTP and test of reciprocal monophyly were performed on the gene trees from MrBayes.

Table 3.

Download as

CSV

XLSX

Specifications of alignments used for gene tree reconstruction with BEAST and MrBayes.

Alignment	Number of sequences	Length (bp)	Variable sites	Parsimony informative sites	Missing data	Analysis	Substitution model
CAD	353	736	37,6%	28,0%	1.5%	MrBayes	1: GTR+I2: F813: HKY+G
CAD haplotypes	136	736	31.8%	26.4%	1.6%	BEAST	1: BMod2: BMod3: BMod
COI	451	658	40,0%	37,7%	18.5%	MrBayes	1: SYM+I+G2: F81+I3: GTR+I+G
COI haplotypes	183	658	39.2%	35.0%	15.9%	BEAST	1: BMod2: BMod3: BMod

The GMYC method uses a haplotype-based ultrametric gene tree to determine the transition from inter- to intraspecific branching patterns (Fujisawa and Barraclough 2013; Pons et al. 2006). For this, we reconstructed chronograms in BEAST (see below). GMYC was performed separately on chronograms from each gene in R 3.2.3 (R Core Team 2015), using the splits package (Ezard et al. 2014). The single-threshold option of GMYC was used. It sets a single limit between inter- and intraspecific divergence patterns and has been shown to outperform the multiple-threshold option (Fujisawa and Barraclough 2013).

Haplotype-based chronograms require that identical sequences be removed from the alignment, leaving an alignment consisting only of unique sequences. Identical haplotypes were removed from the original alignments using collapsetypes_v4.6 (Chesters 2013), which outputs a reduced fasta-alignment and haplotype assignations of each sequence. Ultrametric gene trees based on haplotype alignments were reconstructed in BEAST2 v. 2.3.1 (Bouckaert et al. 2014). Model selection was done using bModeltest (Bouckaert and Drummond 2015), which estimates the best fitting model of sequence evolution simultaneously with the Bayesian tree search. The transition/transversion split option was chosen, which searches among 31 models of sequence evolution. The ‘empirical’ option was used for base frequencies. Alignments were partitioned per codon position with independent rates among partitions; trees and clocks were linked among partitions. Trees were reconstructed under a relaxed lognormal clock, and with a coalescent constant population tree prior. Priors were set to default settings, with infinity values replaced with hard bounds at 1000, to avoid improper priors. Independent analyses were executed with a run time of 0.5 to 1 billion generations each and runs were checked in Tracer. Maximum clade credibility trees were generated in TreeAnnotator, with 10–50% burn-in. PTP and GMYC analyses were performed on two independent trees, to check stability of the results. Rhyacophila polonica McLachlan, 1879 was included as outgroup for MrBayes analyses but was not included in BEAST analyses.

For this paper, we define the use of the terms ‘species’, ‘putative new species’, ‘cluster’ and ‘OTU’ (operational taxonomic unit) as follows: ‘Species’ refers to formally described morphological taxa, following established taxonomy. With ‘putative new species’ we mean morphologically distinct taxa that are still unknown in the literature. The term ‘cluster’ refers to specific results from one of the two analyses, outputting delimited GMYC and PTP ‘clusters’, respectively. For our consensus result from morphology, PTP, and GMYC based on CAD and COI we use the term OTU, which can represent single species or groups of species referred to as species complexes.

Comparative morphological studies

Comparative morphological analysis of larvae followed a standard procedure. We screened all larvae of each OTU for consistent morphological characters. Instar differentiation and thus assignment of most larvae to different instars is not possible with the currently available material. Therefore, general features commonly represented by all OTUs within one phylogenetic clade were considered as synapomorphies. Some characters were found present across all size classes of single OTUs and phylogenetic clades, e.g., distolateral accessory hooks on lateral plates of anal prolegs were consistently present in even the smallest instars.

Results

Datasets

The dataset comprised 525 Himalopsyche individuals (205 adults, 313 larvae and 8 pupae), and R. polonica as outgroup for MrBayes. We generated 352 Himalopsyche sequences of CAD (736 bp) and 450 Himalopsyche sequences of COI (658 bp). After haplotype reduction of the alignments, the CAD alignment had 136 unique haplotypes and the COI had 183 unique haplotypes. The total CAD alignment had 37.6% variable sites; the total COI haplotype alignment had slightly more with 40% (Table 3). For morphological treatment, we also included larvae collected in Japan and Thailand identified as H. japonica and H. acharai, respectively, although no molecular data for these specimens were available.

Life stage association

PTP delimited 29 clusters with CAD, and 62 with COI for the ingroup. GMYC delimited 27 clusters with CAD and 46–48 clusters with COI. The results from separate runs were stable except for GMYC with COI. This instability did not affect larval association and we hereafter only refer to COI run 1 which delimited 48 GMYC clusters (Table 4, Suppl. material 2, 3). We defined OTUs based on adult male morphospecies to the extent that it was possible. Clades containing several paraphyletic species were grouped into ‘species complexes’. In all cases but two, COI clusters were nested within CAD clusters, yielding an overall compatible result. The exceptions were the OTUsjaponica-complex (H. japonica and a putative new Japanese species [H. sp. n. 1529]), and triloba-complex (H. triloba (Hwang, 1958), H. efiel Malicky, 2012, H. hageni Banks, 1940, H. malenanda Schmid, 1963, H. maculipennis (Ulmer, 1905a), and H. yatrawalla Schmid & Botosaneanu, 1966), which both were paraphyletic in COI. The remaining species complexes were: martynovi-complex (H. martynovi Banks, 1940, H. epikur Malicky, 2011), and excisa-complex (H. excisa Ulmer, 1905b, H. placida Banks, 1947 and H. maitreya Schmid, 1963). Three OTUs were identified for which only larval material was available: H. sp. 1196 (L), H. sp. 1338 (L), and H. sp. 1254 (L). A H. platon Malicky, 2011 male formed a clade together with samples of larvae and females in COI but we lacked CAD data from the adult male so we cannot at this stage conclude whether this monophylum constitutes one or several species, and refer to this monophylum as platon-complex.

Table 4.

Download as

CSV

XLSX

Number of clusters delimited by PTP and GMYC. Number of PTP clusters refer to the ingroup only.

Gene	Method	Run	Chain length	Burn-in	ESS	Clusters
CAD	MrBayes & PTP	1	1* 10⁷	25%	All >200	29
CAD	MrBayes & PTP	2	1* 10⁷	25%	All >200	29
CAD	BEAST & GMYC	1	1* 10⁹	50%	All >200	27
CAD	BEAST & GMYC	2	5 * 10⁸	10%	Most >200	27
COI	MrBayes & PTP	1	5* 10⁷	25%	All >200	62
COI	MrBayes & PTP	2	1* 10⁷	25%	Most >200	62
COI	BEAST & GMYC	1	1* 10⁹	10%	All >200	48
COI	BEAST & GMYC	2	1 * 10⁹	10%	All >200	46

We could unambiguously associate 239 larvae and eight pupae to the following nine species: H. acharai, H. anomala Banks, 1940, H. digitata (Martynov, 1935), H. gregoryi (Ulmer, 1932), H. phryganea, H. sylvicola, H. tibetana, H. 677, and H. 685 (Table 5, Suppl. material 1). We could additionally associate 65 larvae to OTU-level for seven OTUS: excisa-complex, martynovi-complex, platon-complex, H. sp. 1196 (L), H. sp. 1254 (L), H. sp. 1338 (L), and triloba-complex.

Table 5.

Download as

CSV

XLSX

Overview of life stage association results for each OTU. Asterisks denote. The OTUsmartynovi-complex, H. 683, and H. 685, formed a single GMYC cluster in CAD, here referred to as ‘martynovi-clade’. Key: ^† indicates larvae without DNA data, * indicates conflicts in COI and CAD gene tree topologies.

OTU	Samples	Species	Larvae	Larval type	CAD units (PTP/GMYC)	COI units (PTP/GMYC)	Association with CAD (PTP/GMYC)	Association with COI (PTP/GMYC)	Reciprocal monophyly (both genes)	Number of criteria fulfilled
H. acharai	4	H. acharai	yes	japonica type	1/1	1/1	yes/yes	yes/yes	yes	5
H. anomala	6	H. anomala	yes	tibetana type	1/1	4/3	yes/yes	no/yes	yes	4
H. auricularis	2	H. auricularis	no	unknown	1/1	2/1	–	–	not applicable	–
H. biansata	2	H. biansata	no	unknown	no data	1/1	–	–	not applicable	–
H. diehli	1	H. diehli	no	unknown	1/1	1/1	–	–	not applicable	–
H. digitata	34	H. digitata	yes	tibetana type	2/1	2/2	yes/yes	no/no	yes	3
H. eos	2	H. eos	no	unknown	no data	1/1	–	–	not applicable	–
excisa-complex	31	H. excisa, H. placida, H. maitreya	yes	tibetana type	1/1	5/2	no/no	no/no	no	0
H. gigantea	2	H. gigantea	no	gigantea type	no data	1/1	–	–	not applicable	–
H. gregoryi	77	H. gregoryi	yes	tibetana type	1/1	1/1	yes/yes	yes/yes	yes	5
H. gyamo	1	H. gyamo	no	unknown	1/1	1/1	–	–	not applicable	–
H. horai	5	H. horai	no	unknown	1/1	1/1	–	–	yes	–
japonica-complex*	8	H. japonica, H. sp. n. 1529	yes^†	japonica type	2/1	4/3	–	–	no	–
H. kuldschensis	2	H. kuldschensis	no	unknown	1/1	1/1	–	–	yes	–
H. lanceolata	4	H. lanceolata	no	unknown	1/1	1/1	–	–	yes	–
H. lepcha	5	H. lepcha	no	unknown	1/1	2/2	–	–	yes	–
H. lua	5	H. lua	no	unknown	no data	2/1	–	–	not applicable	–
martynovi-complex	30	H. martynovi, H. epikur	yes	tibetana type	martynovi–clade	1/1	no/no	no/no	no	0
H. navasi	4	H. navasi	no	unknown	1/1	4/4	–	–	not applicable	–
H. phryganea	4	H. phryganea	yes	phryganea type	1/1	2/2	yes/yes	yes/yes	yes	5
platon-complex	27	H. platon, H.. sp (F), H.. sp (L)	yes	tibetana type	1/1	4/3	no/no	no/no	yes	1
H. sp. 1196 (L)	17	NA	yes	tibetana type	1/1	2/1	–	–	yes	–
H. sp. 1254 (L)	4	NA	yes	tibetana type	1/1	no data	–	–	not applicable	–
H. sp. 1338 (L)	1	NA	yes	gigantea type	1/1	1/1	–	–	not applicable	–
H. sylvicola	33	H. sylvicola	yes	gigantea type	1/1	6/3	yes/yes	no/yes	yes	5
H. tibetana	58	H. tibetana	yes	tibetana type	1/1	2/2	yes/yes	yes/yes	yes	5
H. todma	4	H. todma	no	unknown	1/1	no data	–	–	not applicable	–
triloba-complex*	67	H. triloba, H. hageni, H. macilipennis, H. malenanda, H. efiel, H. yatrawalla	yes	gigantea type	2/2	5/3	no/no	no/no	no	0
H. yongma	1	H. yongma	no	unknown	1/1	1/1	–	–	not applicable	–
H. 677	41	H. 677	yes	tibetana type	1/1	1/1	yes/yes	yes/yes	yes	5
H. 683	8	H. 683	no	tibetana type	martynovi–clade	1/1	–	–	yes	–
H. 685	34	H. 685	yes	tibetana type	martynovi–clade	1/1	no/no	yes/yes	yes	3

Morphology of Himalopsyche larvae

Synapomorphic larval characters of Himalopsyche according to Lepneva (1970), Ross (1956), Flint (1961), Ulmer (1957), and the present study are:

• Mandibles with prominent lateral protuberances (Figure 14, arrow)

• 2^nd and 3^rd leg with anterodorsal single coxal gills (e.g., Figure 18, arrow a, Figure 20, arrow a)

• Abdomen with prominent and complex gills consisting of a multitude of gill filaments positioned on several bases or a single base which can be slightly to distinctly protuberant (e.g. Figure 3).

• Anal proleg with two proximal accessory hooks fused with lateral sclerites (e.g. Figure 38, arrow d, Figure 41, arrow d).

Figures 3–6.

Habitus of Himalopsyche larvae. 3 H. phryganea 4 H. gregoryi 5 H. japonica 6 H. sylvicola. Arrow A points to anterior.

Based on these characters, larvae of Himalopsyche can easily be differentiated from the closely related genera Philocrena and Rhyacophila. The larvae of the monotypic genera Fansipangana and Phoupanpsyche are unknown. Within Himalopsyche, the four different larval types can easily be differentiated based on distinct character states.

phryganea type

Himalopsyche phryganea is the only species known with this larval type. Larvae of H. phryganea were described and illustrated by Flint (1961) and Wiggins (1996). Larvae of the phryganea type are characterised by the following set of characters:

Thorax. Pronotum with a single row of long dark setae along the entire anterior margin; short light recumbent setae concentrated at anterolateral pronotal edges; Sa1 present as a transversal band of 4–5 setae, Sa2 absent (Figure 7); legs without dorsal fringe of setae (Figure 11). Gills. Ventral gills at meso- and metathorax absent (Figure 23); thoracic and abdominal gills arranged on a single suboval, slightly protuberant base extending obliquely from anterodorsal to mediolateral position (Figs 3, 15, 16, 25, 29, 30). Abdomen. Abdomen without ventral protuberances (Figure 30); single ventral medial sclerite on abdomen III-VII, oval, transversally elongated (Figure 30, arrow). Anal prolegs. Stout, distolateral accessory hook absent, dorsal plate with rounded central protuberance (Figs 37–38, arrow b); dorsal spine on basal anal claw dark (Figure 37, arrow c).

Figures 7–10.

Pronotum of Himalopsyche larvae. 7 H. phryganea 8 H. anomala 9 H. japonica 10 H. sylvicola. Arrow A points to anterior.

Figures 11–14.

Mesolegs and head of Himalopsyche larvae. 11 Mesoleg of H. gregoryi 12 Mesoleg of H. sylvicola, arrows indicate dorsal fringe of setae on legs 13 Mesoleg of H. sylvicola. Arrows indicate pennate setae on coxa and femora 14 Head of H. phryganea, arrow points to lateral protuberances on mandible.

Figures 15–22.

Thorax of Himalopsyche larvae in dorsal and lateral view. 15, 16 H. phryganea 17, 18 H. gregoryi, arrow a points to anterodorsal single coxal gill 19, 20 H. japonica, arrow a indicates anterodorsal single coxal gill; arrow b indicates the ventral gills 21, 22 H. sylvicola. Arrows A points to anterior.

tibetana type

One species with larvae of the tibetana type has been described in the larval stage: H. tibetana (Graf and Sharma 1998, Type B). The following OTUs could be assigned to the tibetana type: H. anomala, H. digitata, excisa-complex, H. gregoryi, martynovi-complex, platon-complex, H.. sp. 1196 (L), and H.. sp. 1254 (L), H. tibetana, H. 677, and H. 685. Larvae of the tibetana type are characterised by the following set of characters:

Thorax. Pronotum with two rows of setae along the anterior edge, anteriormost row of setae short, light, recumbent and posterior row setae longer, black; Sa 1 present as a transversal band of 4–5 setae, Sa2 absent (Figure 8); legs without dorsal fringe of setae (Figure 11). Gills. Ventral gills at meso- and metathorax absent (Figure 23); thoracic gills arranged on anterodorsal and anterolateral bases, abdominal gills arranged on anterodorsal, anterolateral and posterolateral bases (Figs 4, 17, 18, 26, 31, 32); anterodorsal and anterolateral gill bases on abdomen I close, seemingly fused. Abdomen. Without ventral protuberances; single ventral medial sclerite on abdomen III-VII oval, transversally elongated (Figure 32, arrow). Anal prolegs. More elongated, with distolateral accessory hook (Figs 39–40, arrow a), dorsal plate with central hook (Figs 39–40, arrow b), dorsal spine on basal anal claw dark (Figs 39–40, arrow c).

Figures 23–28.

Thorax ventral, abdomen dorsal. 23 H. gregoryi 24 H. japonica, arrows indicate ventral gills 25 H. phryganea 26 H. gregoryi, 27 H. acharai 28 H. sylvicola. Arrows A points to anterior.

Figures 29–36.

Abdomen lateral, ventral. 29, 30 H. phryganea, arrow points to ventral medial sclerite 31, 32 H. gregoryi, arrow points to ventral medial sclerite 33, 34 H. japonica, arrow points to ventral protuberance 35, 36 H. sylvicola, arrow a indicates small rounded protuberance on lateral process, while arrow b indicates the ventral medial sclerite. Roman numbers indicate abdominal segments. Arrow A points to anterior.

japonica type

Two species of the japonica type have been described in the larval stage: H. acharai by Thamsenanupap et al. (2005), and H. japonica by Saito (1965) and Tanida (1985). Larvae of the japonica type are characterised by the following set of characters:

Thorax. Pronotum with a single row of setae along the anterior edge; Sa1 present as a transversal row of 2–4 dark setae medially, Sa2 present as a group (sometimes arranged as sagittal band) of 2–4 setae (Figure 9); legs without dorsal fringe of setae (Figure 11). Gills. Ventral gills at meso- and metathorax present (Figure 24, arrows, Figure 20, arrow b); thoracic and abdominal gills arranged on two joint bases, one anterodorsal and one on a small lateral protuberance, extending posterolaterally (Figs 5, 19, 20, 27, 33, 34); abdomen I with gills on a single base only. Abdomen. With ventral protuberances (Fig. 34, arrow); ventral medial sclerites absent (Figure 34). Anal prolegs. Stout, distolateral accessory hook present (Figs 41–42, arrow a), dorsal plate flat without rounded central protuberance; dorsal spine on basal anal claw yellowish (Figure 41, arrow c).

Figures 37–44.

Anal prolegs, lateral, caudal. 37, 38 H. phryganea 39 H. gregoryi 40 H. tibetana 41, 42 H. japonica 43, 44 H. sylvicola. Key: arrows a distolateral accessory hook. arrows b protuberance/hook on dorsal plate. arrows c dorsal spine on basal anal claw. arrows d proximal accessory hooks fused with lateral sclerites. Arrow A points to anterior.

gigantea type

Larvae of this type have been described by Lepneva (1945, 1970), Schmid and Botosaneanu (1966), and Graf and Sharma (1998, Type A). Both Graf and Sharma (1998) and Schmid and Botosaneanu (1966) described larvae of unknown species identity. We assigned the following OTUs to the gigantea type: H. sylvicola, triloba-complex, and H.. sp. 1338 (L). Larvae of the gigantea type are characterised by the following set of characters:

Thorax. Pronotum with a single row of setae along the anterior edge; Sa1 present as a transversal band of 2–4 setae, Sa2 present as a sagittal band of 7–9 dark setae, prominent (10); legs with dorsal fringe of setae (Figure 12), with pennate setae on coxa and femora (Figure 13). Gills. Ventral gills at meso- and metathorax absent (Figure 23); thoracic and abdominal gills arranged on large lateral processes (conical processes sensuGraf and Sharma 1998), lateral processes slightly smaller on meso-and metathorax and distinctly smaller on abdomen I (Figs 6, 21, 22, 28, 35, 36); lateral processes with small rounded protuberances proximoventrally (Figure 36, arrow a). Abdomen. Without ventral protuberances; single ventral medial sclerite on abdomen II-VII suboval, thin, transversally elongated (Figure 36, arrow b). Anal prolegs. Stout, without distolateral accessory hook, dorsal plate flat without protuberance or hook (Figs 43–44); dorsal spine on basal anal claw yellowish (Figs 43–44, arrow c).

Discussion

We used several life stage association strategies based on two genes in a comparative setting: PTP, GMYC and reciprocal monophyly. All strategies acknowledge a successful association in cases of sequence identity, but differ in that PTP and GMYC use branch lengths to estimate which clades represent distinct units, and that reciprocal monophyly requires congruence of gene trees. The importance of using more than one gene (Dupuis et al. 2012) was apparent in this study, as the two genes yielded somewhat different, albeit overall compatible, results. Delimitation results based on COI tended to split species into several units, especially with PTP, where CAD did not. This could be explained by the quicker coalescence time in COI, with its matrilineal inheritance leading to an effective population size 1/4 of that of nuclear genes. The PTP and GMYC methods search for the transition between inter- and intraspecific branching patterns, and if coalescence of COI lineages occurs within isolated populations, then the PTP and GMYC methods cannot distinguish between population and species signals. CAD was generally congruent with established taxonomy except in the species complexes, which were unresolvable by either gene. The major assumption, and limitation, of all association methods employed here is that they assume monophyletic gene trees. Closely related species may be difficult to separate genetically due to incomplete lineage sorting, and recent or ongoing gene flow (e.g., Kutschera et al. 2014). Therefore, recent divergences will always be problematic for PTP and GMYC analyses. GMYC has, for example, been shown to be an accurate and conservative method under conditions where effective sample sizes are low, and divergence times between species are high, since such conditions yield monophyletic species (Esselstyn et al. 2012, Fujisawa and Barraclough 2013). To resolve the species complexes more genes should be used, ideally under a multi-species coalescence approach, such as BPP (Yang 2015), or STACEY (Jones 2017). These methods combine data from several genes to estimate the species borders, accounting for incomplete lineage sorting and conflicting information among genes.

We observed clear morphological differences in morphology between the larval types. Within larval types, however, the examined material did not show any stable and reliable morphological characters to delimitate larvae at species-level. More material, particularly of last instar larvae, and ideally from numerous sites is required to better assess interspecific, intraspecific and ontological variation in the here described as well as other morphological characters. Only then will it be possible to assess if the observed morphological variation is useful for delimiting species in the larval stage.

Conclusion

In this study, we present the characteristic morphological differences of four larval types of Himalopsyche. Life-stage association based on molecular data enabled us to do this, as it provided an OTU assignation for over 300 larvae. We found little or no morphological differences among species within the same type. Once we are able to better discern organisms at lower taxonomic rank (by morphology or molecular association, e.g., in high-throughput barcoding studies), aquatic insects and other benthic invertebrates can become much more valuable for biological monitoring in poorly studied regions. More generally, it is essential to be able to distinguish taxa at the lowest taxonomic resolution (i.e., to species) to understand their ecology and evolution. This is also of great relevance when developing tools to assess ecological status to ensure sustainable use and management of natural resources.

Author contributions

AEH, SCJ, SP, and WG developed and planned the paper. AEH, SCJ, and SP conducted field work. AEH generated molecular, geographic and elevational data and performed molecular association analyses. AEH, SV and WG studied morphology of larvae; AEH drafted larval type descriptions with support from SV and WG. AEH and WG photographed larvae. AEH wrote the initial manuscript and all authors contributed substantially to further versions of the manuscript.

Acknowledgements

We thank Hans Malicky (Lunz am See), Wolfram Mey (Berlin), Bob Wisseman (Oregon), R. B. Kuranishi (Chiba), and Hiroyuki Nishimoto (Komaki) for providing valuable specimens for the study. Field work and molecular work were funded by DFG Grant PA1617/2-1. WG was financially supported by the EC-funded Synthesys Access programme (DE-TAF-5887: Morphological characterisation and description of molecularly associated and differentiated clades of larvae of genus Himalopsyche (Trichoptera, Rhyacophilidae)). SCJ acknowledges funding from the German Federal Ministry of Education and Research (BMBF) for the “GLANCE” project (Global change effects in river ecosystems; 01LN1320A). We thank Ram Devi Tachamo Shah, Deep Narayan Shah, Fengqing Li, Subodh Sharma, Qinghua Cai, Xiaoli Tong, and Felicitas Erzinger (née Hoppeler) for invaluable contributions in the field, as well as the Department of National Parks and Wildlife Conservation (DNPWC) Nepal for providing the research permits.

References

Andreadis KM, Schumann GJ-P, Tamlin Pavelsky T (2013) A simple global river bankfull width and depth database. Water Resources Research 49: 7164–7168. https://doi.org/10.1002/wrcr.20440

Banks N (1940) Report on certain groups of neuropteroid insects from Szechwan, China. Proceedings of the United Stated National Museum 88: 173–220. https://doi.org/10.5479/si.00963801.88-3079.173

Banks N (1947) Some neuropterous insects from Szechwan, China. Fieldiana Zoology 31(12): 97–107. https://doi.org/10.5962/bhl.title.2821

Blahnik RJ, Holzenthal RW (2004) Collection and curation of Trichoptera, with emphasis on pinned material. Nectopsyche, the Neotropical Trichoptera Newsletter 1: 8–20. http://hdl.handle.net/11299/190744

Bouckaert R, Heled J, Kühnert D, Vaughan T, Wu C-H, Xie D, Suchard MA, Rambaut A, Drummond AJ (2014) BEAST 2: A Software Platform for Bayesian Evolutionary Analysis. PLoS Computational Biology 10: e1003537. https://doi.org/10.1371/journal.pcbi.1003537

Bouckaert RR, Drummond AJ (2015) bModelTest: Bayesian phylogenetic site model averaging and model comparison. BMC Evolutionary Biology 17: 42. https://doi.org/10.1186/s12862-017-0890-6

Bukontaite R, Miller KB, Bergsten J (2014) The utility of CAD in recovering Gondwanan vicariance events and the evolutionary history of Aciliini (Coleoptera: Dytiscidae). BMC Evolutionary Biology 14: 5. https://doi.org/10.1186/1471-2148-14-5

Chesters D (2013) Collapsetypes. http://sourceforge.net/projects/collapsetypes [computer software]

Dupuis JR, Roe AD, Sperling FAH (2012) Multi-locus species delimitation in closely related animals and fungi: one marker is not enough. Molecular Ecology 21: 4422–4436. https://doi.org/10.1111/j.1365-294X.2012.05642.x

Ekrem T, Willassen E, Stur E (2010) Phylogenetic utility of five genes for dipteran phylogeny: A test case in the Chironomidae leads to generic synonymies. Molecular Phylogenetics and Evolution 57: 561–571. https://doi.org/10.1016/j.ympev.2010.06.006

Esselstyn JA, Evans BJ, Sedlock JL, Anwarali Khan FA, Heaney LR (2012) Single-locus species delimitation: a test of the mixed Yule-coalescent model, with an empirical application to Philippine round-leaf bats. Proceedings of the Royal Society B: Biological Sciences 279: 3678–3686. https://doi.org/10.1098/rspb.2012.0705

Ezard T, Fujisawa T, Barraclough TG (2014) splits: SPecies’ LImits by Threshold Statistics. R package version 1.0-19/r51. https://R-Forge.R-project.org/projects/splits/

Flint OS (1961) The presumed larva of Himalopsche phryganea (Ross) (Trichoptera: Rhyacophilidae). The Pan-Pacific Entomologist 37: 199–202.

Folmer O, Black M, Hoeh W, Lutz R, Vrijenhoek R (1994) DNA primers for amplification of mitochondrial cytochrome C oxidase subunit I from diverse metazoan invertebrates. Molecular Marine Biology and Biotechnology 3: 294–299. https://www.ncbi.nlm.nih.gov/pubmed/7881515

Foster P, Bergo E, Bourke B, Oliveira T, Nagaki S, Sant’Ana D, Sallum M (2013) Phylogenetic analysis and DNA-based species confirmation in Anopheles (Nyssorhynchus). PLoS ONE 8: e54063. https://doi.org/10.1371/journal.pone.0054063

Fujisawa T, Barraclough TG (2013) Delimiting species using single-locus data and the generalized mixed yule coalescent approach: A revised method and evaluation on simulated data sets. Systematic Biology 62: 707–724. https://doi.org/10.1093/sysbio/syt033

Funk DJ, Omland KE (2003) Species-level paraphyly and polyphyly: Frequency, causes, and consequences, with insights from animal mitochondrial DNA. Annual Review of Ecology, Evolution, and Systematics 34: 397–423. https://doi.org/10.1146/annurev.ecolsys.34.011802.132421

Graf W, Sharma S (1998) The larva of Himalopsyche tibetana (Insecta: Trichoptera, Rhyacophilidae) with some ecological notes on the Genus Himalopsyche Banks, 1940 from Nepal. In: Chalise SR, Khanal NR (Eds) International Conference on Ecohydrology of High Mountain Areas.International Centre for Integrated Mountain Development, Kathmandu, 561–571.

Graf W, Lubini V, Pauls SU (2005) Larval description of Drusus muelleri McLachlan, 1868 (Trichoptera: Limnephilidae) with some notes on its ecology and systematic position within the genus Drusus. Annales de Limnologie - International Journal of Limnology 41: 93–98. https://doi.org/10.1051/limn/2005012

Hebert PDN, Stoeckle MY, Zemlak TS, Francis CM (2004) Identification of birds through DNA barcodes. PLoS Biology 2(10): e312. https://doi.org/10.1371/journal.pbio.0020312

Hjalmarsson AE (submitted) Delimitation and description of three new species of Himalopsyche (Trichoptera: Rhyacophilidae) from the Hengduan Mountains, China.

Hjalmarsson AE, Bergsten J, Monaghan MT (2015) Dispersal is linked to habitat use in 59 species of water beetles (Coleoptera: Adephaga) on Madagascar. Ecography 38: 732–739. https://doi.org/10.1111/ecog.01138

Holzenthal RW, Blahnik RJ, Prather AL, Kjer KM (2007) Order Trichoptera Kirby, 1813 (Insecta), Caddisflies. In: Zhang Z-Q, Shear WA (Eds) Linnaeus Tercentenary: Progress in Invertebrate Taxonomy.Zootaxa 1668: 639–698. http://www.mapress.com/zootaxa/2007f/zt01668p698.pdf

Hwang C (1958) Descriptions of Chinese caddisflies (Trichoptera). Acta Zoologica Sinica 10: 279–284.

Jarvis A, Reuter HI, Nelson A, Guevara E (2008) Hole-filled SRTM for the globe Version 4. Available from the CGIAR-CSI SRTM 90m Database. http://srtm.csi.cgiar.org

Johanson KA, Malm T (2010) Testing the monophyly of Calocidae (Insecta: Trichoptera) based on multiple molecular data. Molecular Phylogenetics and Evolution 54: 535–541. https://doi.org/10.1016/j.ympev.2009.09.025

Jones G (2017) Algorithmic improvements to species delimitation and phylogeny estimation under the multispecies coalescent. Journal of Mathematical Biology 74: 447–467. https://doi.org/10.1007/s00285-016-1034-0

Kearse M, Moir R, Wilson A, Stones-Havas S, Cheung M, Sturrock S, Buxton S, Cooper A, Markowitz S, Duran C, Thierer T, Ashton B, Mentjies P, Drummond A (2012) Geneious Basic: an integrated and extendable desktop software platform for the organization and analysis of sequence data. Bioinformatics 28: 1647–1649. https://doi.org/10.1093/bioinformatics/bts199

Klopfstein S, Vilhelmsen L, Heraty JM, Sharkey M, Ronquist F (2013) The Hymenopteran Tree of Life: Evidence from protein-coding genes and objectively aligned ribosomal data. PLoS ONE 8: e69344. https://doi.org/10.1371/journal.pone.0069344

Kutschera VE, Bidon T, Hailer F, Rodi JL, Fain SR, Janke A (2014) Bears in a forest of gene trees: Phylogenetic inference is complicated by incomplete lineage sorting and gene flow. Molecular Biology and Evolution 31: 2004–2017. https://doi.org/10.1093/molbev/msu186

Lanfear R, Frandsen PB, Wright AM, Senfeld T, Calcott B (2016) PartitionFinder 2: new methods for selecting partitioned models of evolution formolecular and morphological phylogenetic analyses. Molecular biology and evolution. https://doi.org/10.1093/molbev/msw260

Lepneva SG (1945) Remarkable larvae of the genus Rhyacophila Pictet (Trichoptera) from streams of Central Asia. Revue Entomologique de l’URSS 28: 64–74.

Lepneva SG (1956) Novaia licinka ruceinika (Trichoptera, Rhyacophilidre) iz vysokogorii Zakav Kaziia. Revue Entomologique de l’URSS 35(4): 899–911.

Lepneva SG (1970) Fauna of the U.S.S.R, Trichoptera, Volume II No. 1, Larvae and Pupae of Annulipalpia. Smithsonian Institution and the National Science Foundation, Washington, DC, 638 pp.

Malicky H (2004) Atlas of European Trichoptera (2^nd edn). Springer Netherlands, Dordrecht.

Malicky H (2008) Beschreibungen von neuen Trichopteren aus Asien. Braueria 35: 45–57.

Malicky H (2011) Neue Trichopteren aus Europa und Asien. Braueria 38: 23–43.

Malicky H (2012) Neue asiatische Köcherfliegen aus neuen Ausbeuten (Insecta, Trichoptera) Linzer biologische Beiträge 44(2): 1263–1310.

Malicky H, Chantaramongkol P (1989) Einige Rhyacophilidae aus Thailand (Rhyacophilidae). Entomologische Zeitschrift 99(3): 17–24.

Martynov AV (1914) Contributions a la faune des Trichoptères des possessions Russes dans l’Asie Centrale. Annuaire du Musée Zoologique de L'Académie Impériale des sciences de Petrograd 19(3): 401–437.

Martynov AV (1930) On the Trichopterous fauna of China and Eastern Tibet. Proceedings of the Zoological Society of London 5: 65–112.

Martynov AV (1935) A collection of Trichoptera from the Indian Museum. Records of the Indian Museum 37: 93–209.

Martynov AV (1936) On a collection of Trichoptera from the Indian Museum. Records of the Indian Museum.38: 239–306.

McLachlan R (1879) A monographic revision and synopsis of the Trichoptera of the European fauna. Part VIII, John van Vorrt, London, 702 pp.

Mey W (1996) Die Köcherfliegenfauna des Fan Si Pan-Massivs in Nord-Vietnam. Beiträge zur Entomologie 46(1): 36–65. https://doi.org/10.21248/contrib.entomol.46.1.39-65

Montero-Pau J, Gómez A, Muñoz J (2008) Application of an inexpensive and high-throughput genomic DNA extraction method for the molecular ecology of zooplanktonic diapausing eggs. Limnology and Oceanography: Methods 6: 218–222. https://doi.org/10.4319/lom.2008.6.218

de Moor FC, Ivanov VD (2008) Global diversity of caddisflies (Trichoptera: Insecta) in freshwater. Hydrobiologia 595: 393–407. https://doi.org/10.1007/s10750-007-9113-2

Morton KJ (1900) Descriptions of new species of Oriental Rhyacophila. Transactions of the Entomological Society of London 1: 1–7.

Moulton JK, Wiegmann BM (2004) Evolution and phylogenetic utility of CAD (rudimentary) among Mesozoic-aged Eremoneuran Diptera (Insecta). Molecular Phylogenetics and Evolution 31: 363–378. https://doi.org/10.1016/S1055-7903(03)00284-7

Pictet FJ (1834) Recherches pour servir à l'histoire et à l'anatomie des Phryganides. Abraham Cherbuliez, Geneve, 300 pp. https://doi.org/10.5962/bhl.title.8547

Pons J, Barraclough TG, Gomez-Zurita J, Cardoso A, Duran DP, Hazell S, Kamoun S, Sumlin WD, Vogler AP (2006) Sequence-based species delimitation for the DNA Taxonomy of undescribed insects. Systematic Biology 55: 595–609. https://doi.org/10.1080/10635150600852011

R Core Team (2015) R: A language and environment for statistical computing. - R Foundation for Statistical Computing, Vienna, Austria. http://www.R-project.org/

Rainbow PS, Hildrew AG, Smith BD, Geatches T, Luoma SN (2012) Caddisflies as biomonitors identifying thresholds of toxic metal bioavailability that affect the stream benthos. Environmental Pollution 166: 196–207. https://doi.org/10.1016/j.envpol.2012.03.017

Rambaut A, Suchard MA, Xie D, Drummond AJ (2014) Tracer v1.6. http://tree.bio.ed.ac.uk/software/tracer/

Resh VH, Unzicker JD (1975) Water quality monitoring and aquatic organisms: the importance of species identification. Journal of the Water Pollution Control Federation 47: 9–19. http://www.jstor.org/stable/25038592

Ronquist F, Teslenko M, van der Mark P, Ayres DL, Darling A, Höhna S, Larget B, Liu L, Suchard MA, Huelsenbeck JP (2012) MrBayes 3.2: Efficient Bayesian Phylogenetic Inference and Model Choice Across a Large Model Space. Systematic Biology 61(3): 539–542. https://doi.org/10.1093/sysbio/sys029

Ross HH (1941) Descriptions and Records of North American Trichoptera. Transactions of the American Entomological Society 67(1): 35–126.

Ross HH (1956) Evolution and classification of the mountain caddisflies. The University of Illinois Press, Urbana, 213 pp.

Ruiter DE, Boyle EE, Zhou X (2013) DNA barcoding facilitates associations and diagnoses for Trichoptera larvae of the Churchill (Manitoba, Canada) area. BMC Ecology 13: 5. https://doi.org/10.1186/1472-6785-13-5

Saito T (1965) Larvae of two fresh water insects, Himalopsyche japonica and Neohapalothrix kanii. Report of the Japanese Society of Systematic Zoology 1: 14–16.

Schmid F (1963) Quelques Himalopsyche indiennes (Trichopères, Rhyacophilidae). Bonner zoologische Beiträge 14: 206–223.

Schmid F, Botosaneanu L (1966) Le genre Himalopsyche Banks (Trichoptera, Rhyacophilidae). Annals of the Entomological Society of Quebec 11: 123–176.

Schmidt-Kloiber A, Graf W, Moog O, Lorenz A, Hering D (2008) The Indicator Database for European Freshwater Invertebrates. Proceedings of the first conference on faunistics and zoogeography of European Trichoptera 55: 85–88.

Sharanowski BJ, Dowling APG, Sharkey MJ (2011) Molecular phylogenetics of Braconidae (Hymenoptera: Ichneumonoidea), based on multiple nuclear genes, and implications for classification. Systematic Entomology 36: 549–572. https://doi.org/10.1111/j.1365-3113.2011.00580.x

Song JH, Ahn KJ (2014) Species delimitation in the Aleochara fucicola species complex (Coleoptera: Staphylinidae: Aleocharinae) and its phylogenetic relationships. Zoologica Scripta 43: 629–640. https://doi.org/10.1111/zsc.12077

Stephens JF (1836) Illustrations of British entomology; or, a synopsis of indigenous insects: containing their generic and specific distinctions; with an account of their metamorphoses, times of appearance, localities, food, and economy, as far as practicable. Vol. VI. Mandibulata. Baldwin & Cradock, London, 240 pp.

Tachamo Shah RD, Sharma S, Haase P, Jähnig SC, Pauls SU (2015) The climate sensitive zone along an altitudinal gradient in central Himalayan rivers: a useful concept to monitor climate change impacts in mountain regions. Climatic Change 132(2): 265–278. https://doi.org/10.1007/s10584-015-1417-z

Tanida K (1985) Trichoptera. In: Kawai T (Ed.) An illustrated book of aquatic insects of Japan.Tokai University Press, Tokyo, 167–215.

Thamsenanupap P, Chantaramongkol P, Malicky H (2005) Description of caddis larvae (Trichoptera) from northern Thailand of the genera Himalopsyche (Rhyacophilidae), Arctopsyche (Arctopsychidae), cf. Eoneureclipsis (Psychomiidae) and Inthanopsyche (Odontoceridae). Braueria 32: 7–11.

Thompson JD, Higgins DG, Gibson TJ (1994) CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Research 22: 4673–4680. https://doi.org/10.5962/bhl.title.8547

Ulmer G (1905a) Zur Kenntnis aussereuropäischer Trichopteren. Stettin Entomologische Zeitung 66: 107–110.

Ulmer G (1905b) Neue und wenig bekannte Trichopteren der Museen zu Brüssel und Paris. Annales de la Société entomologique de Belgique.49: 17–42.

Ulmer G (1932) Aquatic Insects of China Article III. Neue Chinesiche Trichopteren, nebst Übersicht über die bisher aus China bekannten Arten. Peking Natural History Bulletin 7: 39–70.

Ulmer G (1957) Köcherfliegen (Trichopteren) von den Sund a-Inseln. Teil III. Larven und puppen der Annulipalpia. Archiv für Hydrobiologie - Supplements 23: 109–470.

Verdonschot PFM (2006) Beyond Masses and Blooms: The indicative value of Oligochaetes. Hydrobiologia 564: 127–142. https://doi.org/10.1007/s10750-005-1714-z

Vitecek S, Kučinić M, Previšić A, Živić I, Stojanović K, Keresztes L, Bálint M, Hoppeler F, Waringer J, Graf W, Pauls SU (2017) Integrative taxonomy by molecular species delimitation: multi-locus data corroborate a new species of Balkan Drusinae micro-endemics. BMC Evolutionary Biology 17: 129. https://doi.org/10.1186/s12862-017-0972-5

Wiggins GB (1996) Larvae of the North American Caddisfly Genera (Trichoptera). 2nd ed. University of Toronto Press, Toronto, 457 pp.

Wild AL, Maddison DR (2008) Evaluating nuclear protein-coding genes for phylogenetic utility in beetles. Molecular Phylogenetics and Evolution 48: 877–891. https://doi.org/10.1016/j.ympev.2008.05.023

Yang Z (2015) The BPP program for species tree estimation and species delimitation. Current Zoology 61(5): 854–865. https://doi.org/10.1093/czoolo/61.5.854

Zhang J, Kapli P, Pavlidis P, Stamatakis A (2013) A general species delimitation method with applications to phylogenetic placements. Bioinformatics 29(22): 2869–2876https://doi.org/10.1093/bioinformatics/btt499

Zhou X, Kjer KM, Morse JC (2007) Associating larvae and adults of Chinese Hydropsychidae caddisflies (Insecta: Trichoptera) using DNA sequences. Journal of the North American Benthological Society 26: 719–742. https://doi.org/10.1899/06-089.1

Zhou X, Frandsen PB, Holzenthal RW, Beet CR, Bennett KR, Blahnik RJ, Bonada N, Cartwright D, Chuluunbat S, Cocks GV, Collins GE, deWaard J, Dean J, Flint OS, Hausmann A, Hendrich L, Hess M, Hogg ID, Kondratieff BC, Malicky H, Milton MA, Morinière J, Morse JC, Ngera Mwangi F, Pauls SU, Razo Gonzalez M, Rinne A, Robinson JL, Salokannel J, Shackleton M, Smith B, Stamatakis A, StClair R, Thomas JA, Zamora-Muñoz C, Ziesmann T, Kjer KM (2016) The Trichoptera barcode initiative: a strategy for generating a species-level Tree of Life. Philosophical Transactions B 371: 20160025. https://doi.org/10.1098/rstb.2016.0025

Supplementary materials

Supplementary material 1

List of specimens

Anna E. Hjalmarsson, Wolfram Graf, Sonja C. Jähnig, Simon Vitecek, Steffen U. Pauls

This dataset is made available under the Open Database License (http://opendatacommons.org/licenses/odbl/1.0/). The Open Database License (ODbL) is a license agreement intended to allow users to freely share, modify, and use this Dataset while maintaining this same freedom for others, provided that the original source and author(s) are credited.

Download file (61.60 kb)

Supplementary material 2

Figure S1. Life stage association results based on CAD

Anna E. Hjalmarsson, Wolfram Graf, Sonja C. Jähnig, Simon Vitecek, Steffen U. Pauls

This dataset is made available under the Open Database License (http://opendatacommons.org/licenses/odbl/1.0/). The Open Database License (ODbL) is a license agreement intended to allow users to freely share, modify, and use this Dataset while maintaining this same freedom for others, provided that the original source and author(s) are credited.

Download file (6.13 MB)

Supplementary material 3

Figure S2. Life stage association results based on COI

Anna E. Hjalmarsson, Wolfram Graf, Sonja C. Jähnig, Simon Vitecek, Steffen U. Pauls

This dataset is made available under the Open Database License (http://opendatacommons.org/licenses/odbl/1.0/). The Open Database License (ODbL) is a license agreement intended to allow users to freely share, modify, and use this Dataset while maintaining this same freedom for others, provided that the original source and author(s) are credited.

Download file (6.60 MB)