A time-calibrated phylogeny of southern hemisphere stoneflies: Testing for Gondwanan origins

ebony underpins Batesian mimicry in melanic stoneflies

The evolution of Batesian mimicry - whereby harmless species avoid predation through their resemblance to harmful species - has long intrigued biologists. In rare cases, Batesian mimicry is linked to intraspecific colour variation, in which only some individuals within a population resemble a noxious 'model'. Here, we assess intraspecific colour variation within a widespread New Zealand stonefly, wherein highly melanized individuals of Zelandoperla closely resemble a chemically defended aposematic stonefly, Austroperla cyrene. We assess convergence in the colour pattern of these two species, compare their relative palatability to predators, and use genome-wide association mapping to assess the genetic basis of this resemblance. Our analysis reveals that melanized Zelandoperla overlap significantly with Austroperla in colour space but are significantly more palatable to predators, implying that they are indeed Batesian mimics. Analysis of 194,773 genome-wide SNPs reveals an outlier locus (ebony) strongly differentiating melanic versus non-melanic Zelandoperla. Genotyping of 338 specimens from a single Zelandoperla population indicates that ebony explains nearly 70% of the observed variance in melanism. As ebony has a well-documented role in insect melanin biosynthesis, our findings indicate this locus has a conserved function across deeply divergent hexapod lineages. Distributional records suggest a link between the occurrence of melanic Zelandoperla and the forested ecosystems where the model Austroperla is abundant, suggesting the potential for adaptive shifts in this system underpinned by environmental change.

South African Lagerstätte reveals middle Permian Gondwanan lakeshore ecosystem in exquisite detail

Continental ecosystems of the middle Permian Period (273–259 million years ago) are poorly understood. In South Africa, the vertebrate fossil record is well documented for this time interval, but the plants and insects are virtually unknown, and are rare globally. This scarcity of data has hampered studies of the evolution and diversification of life, and has precluded detailed reconstructions and analyses of ecosystems of this critical period in Earth’s history. Here we introduce a new locality in the southern Karoo Basin that is producing exceptionally well-preserved and abundant fossils of novel freshwater and terrestrial insects, arachnids, and plants. Within a robust regional geochronological, geological and biostratigraphic context, this Konservat- and Konzentrat-Lagerstätte offers a unique opportunity for the study and reconstruction of a southern Gondwanan deltaic ecosystem that thrived 266–268 million years ago, and will serve as a high-resolution ecological baseline towards a better understanding of Permian extinction events.

Biogeography and eye size evolution of the ogre-faced spiders

Net-casting spiders (Deinopidae) comprise a charismatic family with an enigmatic evolutionary history. There are 67 described species of deinopids, placed among three genera, Deinopis, Menneus, and Asianopis, that are distributed globally throughout the tropics and subtropics. Deinopis and Asianopis, the ogre-faced spiders, are best known for their giant light-capturing posterior median eyes (PME), whereas Menneus does not have enlarged PMEs. Molecular phylogenetic studies have revealed discordance between morphology and molecular data. We employed a character-rich ultra-conserved element (UCE) dataset and a taxon-rich cytochrome-oxidase I (COI) dataset to reconstruct a genus-level phylogeny of Deinopidae, aiming to investigate the group's historical biogeography, and examine PME size evolution. Although the phylogenetic results support the monophyly of Menneus and the single reduction of PME size in deinopids, these data also show that Deinopis is not monophyletic. Consequently, we formally transfer 24 Deinopis species to Asianopis; the transfers comprise all of the African, Australian, South Pacific, and a subset of Central American and Mexican species. Following the divergence of Eastern and Western deinopids in the Cretaceous, Deinopis/Asianopis dispersed from Africa, through Asia and into Australia with its biogeographic history reflecting separation of Western Gondwana as well as long-distance dispersal events.

Comparative analysis of mitochondrial genomes among the family Peltoperlidae (Plecoptera: Systellognatha) and phylogenetic implications

Nowadays, the position of Peltoperlidae in Systellognatha has been resolved based on morphological analyses. However, there are different opinions based on molecular data. To date, only three peltoperlid mitogenomes are available, and more sampling is needed to obtain precise phylogenetic relationships. In this study, we obtained the complete mitogenomes of Cryptoperla kawasawai (15,832 bp) and Peltoperlopsis sagittata (15,756 bp). Our results show that gene content, gene order, DmTTF binding site, nucleotide composition, codon usage, ribonucleic acid (RNA) structure, and structural elements in the control region are highly conserved in peltoperlids. Heatmap analysis of codon usage shows that the AT-rich codons UUA, AUU, UUU, and AUA were commonly used codons in the Peltoperlidae. Evolutionary rate analyses of protein-coding genes reveal that different genes have been subject to different rates of molecular evolution correlated with the GC content. All tRNA genes in peltoperlid mitogenomes have a canonical cloverleaf secondary structure except for trnS1, whose dihydrouridine arm simply forms a loop. The control region of the family has several distinct structural characteristics and has the potential to serve as effective phylogenetic markers. Phylogenetic analyses support the monophyly of Perloidea, but the monophyly of Pteronarcyoidea is still not supported. The Peltoperlidae is placed as the earliest branch within the Systellognatha, and the estimated phylogenetic relationship is: Peltoperlidae + {(Styloperlidae + Pteronarcyidae) + [Perlidae + (Chloroperlidae + Perlodidae)]}. Our results provide new insight into the phylogeny of this group.

Systematic revision and phylogeny of Paragripopteryx Enderlein, 1909 (Plecoptera: Gripopterygidae)

Among the Neotropical gripopterygids, the genus Paragripopteryx occurs along the Brazilian Atlantic coast to Uruguay. Since its first recognition by Enderlein in 1909, the genus underwent a confusing taxonomic history with some combinations. In this study, we aim to revise Paragripopteryx and present the first morphology-based phylogeny for the genus. The analysis comprised 38 morphological characters and their respective states in 30 terminal taxa, including 13 Paragripopteryx species, two new species identified as belonging to Paragripopteryx, and 15 outgroup species, among which we can highlight 12 different South American genera and one Australian genus of Gripopterygidae. The cladistic analysis yielded a parsimonious tree for k = 3 (137 steps, consistency index = 0.445, and retention index = 0.591) where most Paragripopteryx are nested, except for Uruguayan Paragripopteryx munoai. We can then infer that in its current circumscription Paragripopteryx is polyphyletic. The following two species are described: Paragripopteryx dasalmas sp. nov. and Paragripopteryx ogum sp. nov. Paragripopteryx baratinii is designated as a nomen dubium. Additionally, we provide a key for species identification, updated geographical records, and illustrations for all species. As a corollary, our study gathers relevant morphological information that can help to better understand this genus and create foundations for the next steps.

Similar pattern, different paths: tracing the biogeographical history of Megaloptera (Insecta: Neuropterida) using mitochondrial phylogenomics

The sequential breakup of the supercontinent Pangaea since the Middle Jurassic is one of the crucial factors that has driven the biogeographical patterns of terrestrial biotas. Despite decades of effort searching for concordant patterns between diversification and continental fragmentation among taxonomic groups, increasing evidence has revealed more complex and idiosyncratic scenarios resulting from a mixture of vicariance, dispersal and extinction. Aquatic insects with discreet ecological requirements, low vagility and disjunct distributions represent a valuable model for testing biogeographical hypotheses by reconstructing their distribution patterns and temporal divergences. Insects of the order Megaloptera have exclusively aquatic larvae, their adults have low vagility, and the group has a highly disjunct geographical distribution. Here we present a comprehensive phylogeny of Megaloptera based on a large-scale mitochondrial genome sequencing of 99 species representing >90% of the world genera from all major biogeographical regions. Molecular dating suggests that the deep divergence within Megaloptera pre-dates the breakup of Pangaea. Subsequently, the intergeneric divergences within Corydalinae (dobsonflies), Chauliodinae (fishflies) and Sialidae (alderflies) might have been driven by both vicariance and dispersal correlated with the shifting continent during the Cretaceous, but with strikingly different and incongruent biogeographical signals. The austral distribution of many corydalids appears to be a result of colonization from Eurasia through southward dispersal across Europe and Africa during the Cretaceous, whereas a nearly contemporaneous dispersal via northward rafting of Gondwanan landmasses may account for the colonization of extant Eurasian alderflies from the south.

Global dispersal and diversification in ground beetles of the subfamily Carabinae

The origin and diversification process of lineages of organisms that are currently widely distributed among continents is an interesting subject for exploring the evolutionary history of global species diversity. Ground beetles of the subfamily Carabinae are flightless except for one lineage, but nevertheless occur on all continents except Antarctica. Here, we used sequence data from ultraconserved elements to reconstruct the phylogeny, divergence time, biogeographical history, ancestral state of hind wings and changes in the speciation rate of Carabinae. Our results show that Carabinae originated in the Americas and diversified into four tribes during the period from the late Jurassic to the late Cretaceous, with two in South America (Celoglossini) and Australasia (Pamborini) and two in Laurasia (Cychrini and Carabini). The ancestral Carabinae were inferred to be winged; three of four tribes (Cychrini, Ceglossini and Pamborini) have completely lost their hind wings and flight capability. The remaining tribe, Carabini, diverged into the subtribes Carabina (wingless) and Calosomina (winged) in the Oligocene. Carabina originated in Europe, spread over Eurasia and diversified into approximately 1000 species, accounting for around 60% of all Carabinae species. Calosomina that were flight-capable dispersed from North America or Eurasia to South America, Australia, and Africa, and then flightless lineages evolved on oceanic islands and continental highlands. The speciation rate increased in the Cychrini and Carabini clades in Eurasia. Within Carabini, the speciation rate was higher for wingless than winged states. Our study showed that the global distribution of Carabinae resulted from ancient dispersal before the breakup of Gondwana and more recent dispersal through flight around the world. These patterns consequently illustrate the causal relationships of geographical history, evolution of flightlessness, and the global distribution and species diversity of Carabinae.

Rewinding the molecular clock in the genus Carabus (Coleoptera: Carabidae) in light of fossil evidence and the Gondwana split: A reanalysis

Molecular clocks have become powerful tools given increasing sequencing and fossil resources. However, calibration analyses outcomes depend on the choice of priors. Here, we revisited the seminal dating study published by Andújar and coworkers of the genus Carabus proposing that prior choices need re-evaluation. We hypothesized that reflecting fossil evidence and the Gondwanan split properly significantly rewinds the molecular clock. We re-used the dataset including five mitochondrial and four nuclear DNA fragments with a total length of 7888 nt. Fossil evidence for Oligocene occurrence of Calosoma was considered. Root age was set based on the fossil evidence of Harpalinae ground beetles in the Upper Cretaceous. Paleogene divergence of the outgroup taxa Ceroglossini and Pamborini is introduced as a new prior based on current paleontological and geological literature. The ultrametric time-calibrated tree of the extended nd5 dataset resulted in a median TMRCA Carabus of 53.92 Ma (HPD 95% 45.01–63.18 Ma), roughly 30 Ma older than in the Andújar study. The splits among C. rugosus and C. morbillosus (A), C. riffensis from the European Mesocarabus (B), and Eurycarabus and Nesaeocarabus (C) were dated to 17.58 (12.87–22.85), 24.14 (18.02–30.58), and 21.6 (16.44–27.43) Ma. They were decidedly older than those previously reported (7.48, 10.93, and 9.51 Ma). These changes were driven almost entirely by constraining the Carabidae time-tree root with a Harpalinae amber fossil at ~99 Ma. Utilizing the nd5 dating results of three well-supported Carabus clades as secondary calibration points for the complete MIT-NUC dataset led to a TMRCA of Carabus of 44.72 (37.54–52.22) Ma, compared with 25.16 Ma (18.41–33.04 Ma) in the previous study. Considering fossil evidence for Oligocene Calosoma and Late Cretaceous Harpalini together with the Gondwanan split as a new prior, our new approach supports the origin of genus Carabus in the Eocene. Our results are preliminary because of the heavy reliance on the nd5 gene, and thus will have to be tested with a sufficient set of nuclear markers. Additionally, uncertainties due to dating root age of the tree based on a single fossil and outgroup taxon affect the results. Improvement of the fossil database, particularly in the supertribe Carabitae, is needed to reduce these uncertainties in dating Carabus phylogeny.

Ancient relicts or recent immigrants? Different dating strategies alter diversification scenarios of New Zealand aquatic beetles (Coleoptera: Hydrophilidae: Berosus)

Dated species-level phylogenies are crucial for understanding the origin and evolutionary history of modern faunas, yet difficult to obtain due to the frequent absence of suitable age calibrations at species level. Substitution rates of related or more inclusive clades are often used to overcome this limitation but the accuracy of this approach remains untested. We compared tree dating based on substitution rates with analyses implementing fossil data by direct node-dating and indirect root-age constraints for the New Zealand endemic Berosus water beetles (Coleoptera: Hydrophilidae). The analysis based solely on substitution rates indicated a Miocene colonization of New Zealand and Pleistocene origin of species. By contrast, all analyses that implemented fossil data resulted in significantly older age estimates, indicating an ancient early Cenozoic origin of the New Zealand clade, diversification of species during or after the Oligocene transgression and Miocene-Pliocene origin of within-species population structure. Rate-calibrated time trees were incongruent with recently published Coleoptera time trees, the fossil record of Berosus and the distribution of outgroup species. Strong variation of substitution rates among Coleoptera lineages, as well as among lineages within the family Hydrophilidae, was identified as the principal reason for low accuracy of rate-calibrated analyses, resulting in underestimated node ages in Berosus. We provide evidence that Oligocene to Pliocene events, rather than the Pleistocene Glacial cycles, played an essential role in the formation of the modern New Zealand insect fauna.

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.

Perlidae (Plecoptera) from the Paranapiacaba Mountains, Atlantic Forest, Brazil: Diversity and implications of the integrative approach and teneral specimens on taxonomy

The study of complementary sources of biological variation (e.g. morphological, molecular) has allowed a better understanding of biodiversity through the construction of an integrative taxonomy. Using this approach, specimens from the Paranapiacaba Mountains, southeastern Brazil, were studied to update the knowledge on the stonefly family Perlidae from the region, characterize the species, and make associations between nymphs and adults using a fragment of the mitochondrial cytochrome c oxidase subunit I (COI) gene. The study also discusses the implications of integrative taxonomy and teneral specimens for the study of South American Perlidae. The molecular data were analyzed using Bayesian inference, Neighbor-joining, and delimiting species methods. Our results revealed that, in general, there was a morphological and molecular congruence between species. In the Paranapiacaba Mountains, three genera and 15 species were recorded: Anacroneuria boraceiensis Froehlich 2004, A. debilis (Pictet 1841) (new record), A. fiorentini De Ribeiro and Froehlich 2007 (new record), A. flintorum Froehlich 2002, A. iporanga Bispo and Froehlich 2004, A. itajaimirim Bispo and Froehlich 2004, A. polita (Burmeister 1913), A. subcostalis Klapálek 1921, A. tupi Bispo and Froehlich 2004 (with a description of the nymph), Kempnyia auberti Froehlich 1996, K. colossica (Navás 1934), K. flava Klapálek 1916, K. neotropica (Jacobson and Bianchi 1905) (including its new junior synonym K. petersorum Froehlich 1996), Kempnyia sp., and Macrogynoplax veneranda Froehlich 1984. COI sequences were obtained for 11 species, five of which had nymphs associated with adults. Among the five associated nymphs, the nymph of A. tupi is described here. The results of this study indicate that the color of adult teneral specimens differs from that of mature specimens. Given this, the synonym of K. neotropica and K. petersorum was proposed since these species have high morphological and molecular similarities and differ only in color patterns. In addition, the previous record of A. petersi Froehlich 2002 from the Paranapiacaba Mountains was invalidated since it was considered a teneral specimen of A. flintorum. These results suggest that the development of an integrative taxonomy is essential to continue advancing the study of Perlidae diversity in South America.

Molecular phylogeny inferred from the mitochondrial genomes of Plecoptera with Oyamia nigribasis (Plecoptera: Perlidae)

In this study, the mitochondrial genome of the stonefly, Oyamia nigribasis Banks, 1920 (Plecoptera: Perlidae), was sequenced and compared with the mtDNA genomes of 38 other stoneflies and two Ephemerae. The O. nigribasis mitogenome is a circular 15,923 bp molecule that encodes a large, noncoding control region (CR) and 37 typical mtDNA genes; these include 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), and two ribosomal RNA genes (rRNAs), respectively. Most of the PCGs initiated with ATN and terminated with TAN. The dihydrouridine (DHU) arm of tRNASer (AGN) was missing, whereas the other 21 tRNAs all exhibited the typical cloverleaf secondary structure. Stem-loop (SL) structures and tandem repeats were identified in the CR. Phylogenetic analyses using Bayesian inference and maximum likelihood were undertaken to determine relationships between stoneflies. Results indicated that the Antarctoperlaria, which contains Gripopterygidae, was absolutely separated from Arctoperlaria; this finding agrees with morphology. Finally, the overall relationships could be summarized as follows ((((Notonemouridae + Nemouridae) + Leuctridae) + (Scopuridae + (Capniidae + Taeniopterygidae))) + (((Perlodidae + Chloroperlidae) + Perlidae) + (Pteronarcyidae + (Peltoperlidae + Styloperlidae))) + ((Diamphipnoidae + Eustheniidae) + Gripopterygidae)).

A first cladistic analysis of Antarctoperlinae (Plecoptera: Gripopterygidae) and a new micropterous species from Patagonia

Gripopterygidae is a diverse family of stoneflies, Plecoptera, distributed in the Southern Hemisphere. It has been traditionally divided into five subfamilies, but the monophyly of most of these is not supported by molecular the more comprehensive molecular analysis of the order. To test the monophyly of Antarctoperlinae, and to establish the phylogenetic position of a new Gripopterygidae species, we performed a morphological cladistic analysis including 38 morphological characters and 27 terminal taxa, with representatives of the four subfamilies present in South America and three Austroperlidae. Based on published information, we rooted the tree with Penturoperla barbata, Austroperlidae. As a result, Antarctoperlinae was recovered as polyphyletic, with Vesicaperla kuscheli and Plegoperla punctata, two members of the subfamily, placed outside the clade that includes the nine remaining genera of Antarctoperlinae. Vesicaperla also falls outside the family in previous molecular analysis. Based on this evidence, it should not be placed in Antarctoperlinae. Plegoperla punctata, known only from the type series, possesses many missing entries in our data matrix. Based on this, it seems convenient to maintain its subfamilial placement. In the tree obtained, the potential new species nests together with Pehuenioperla llaima. We thus accept it as a member of Pehuenioperla and describe it as P. microptera sp. nov.

Phylogenetic divergence of island biotas: Molecular dates, extinction, and "relict" lineages

Island formation is a key driver of biological evolution, and several studies have used geological ages of islands to calibrate rates of DNA change. However, many islands are home to "relict" lineages whose divergence apparently pre-dates island age. The geologically dynamic New Zealand (NZ) archipelago sits upon the ancient, largely submerged continent Zealandia, and the origin and age of its distinctive biota have long been contentious. While some researchers have interpreted NZ's biota as equivalent to that of a post-Oligocene island, a recent review of genetic studies identified a sizeable proportion of pre-Oligocene "relict" lineages, concluding that much of the biota survived an incomplete drowning event. Here, we assemble comparable genetic divergence data sets for two recently formed South Pacific archipelagos (Lord Howe; Chatham Islands) and demonstrate similarly substantial proportions of relict lineages. Similar to the NZ biota, our island reviews provide surprisingly little evidence for major genetic divergence "pulses" associated with island emergence. The dominance of Quaternary divergence estimates in all three biotas may highlight the importance of rapid biological turnover and new arrivals in response to recent climatic and/or geological disturbance and change. We provide a schematic model to help account for discrepancies between expected versus observed divergence-date distributions for island biotas, incorporating the effects of both molecular dating error and lineage extinction. We conclude that oceanic islands can represent both evolutionary "cradles" and "museums" and that the presence of apparently archaic island lineages does not preclude dispersal origins.

Insights into the molecular phylogeny of Rhaphidophoridae, an ancient, worldwide lineage of Orthoptera

We investigated the molecular phylogenetic divergence and historical biogeography of cave crickets belonging to the family Rhaphidophoridae (Orthoptera, Ensifera). We used taxa representative of most of the regions embraced by the family, considering samples of Macropathinae from Gondwana land (i.e., Tasmania, Australia, New Zealand, South Africa, and South America); Aemodogryllinae and Rhaphidophorinae from Southern-eastern Asia (i.e., India, Bhutan, China, Philippines and the Sulawesi islands); Dolichopodainae and Troglophilinae from the Mediterranean region and Ceuthophilinae from North America. Based on previous papers, we carried out an analysis of both mitochondrial and nuclear DNA sequences considering the ribosomal RNA units 12S, 16S, 18S, and 28S. To reconstruct phylogeny, we use cladistics, Maximum Likelihood (ML), and Bayesian analyses. All phylogenetic analyses showed the same highly supported topology generally congruent with the classical systematic arrangement at the level of each sub-family but strongly disagree with previous affinity hypotheses between sub-families based on morphological characters. Our results reveal a close affinity between Asiatic and Gondwanian taxa from one hand and between North American and Mediterranean ones from the other hand. Dating estimates indicated that Rhaphidophoridae originated in the Cretaceous period during the Mesozoic era with the ancestral area located both in the northern and southern hemisphere. A possible biogeographic scenario, reconstructed using S-DEC with RASP software, suggested that the current distribution of Rhaphidophoridae might be explained by a combination of both dispersal and vicariance events occurred especially in the ancestral populations. The radiation of Rhaphidophoridae started within the Pangaea, where the ancestor of Rhaphidophoridae occurred throughout an ancestral area including Australia, North America, and the Mediterranean region. The opening of the Atlantic Ocean promoted the divergence of North American and Mediterranean lineages while the differentiation of the southern lineages, spread from Australia, appears to be related to the fragmentation of Gondwana land.

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).

Molecular phylogeny and diversification timing of the Nemouridae family (Insecta, Plecoptera) in the Japanese Archipelago

The generation of the high species diversity of insects in Japan was profoundly influenced by the formation of the Japanese Archipelago. We explored the species diversification and biogeographical history of the Nemouridae Billberg, 1820 family in the Japanese Archipelago using mitochondrial DNA and nuclear DNA markers. We collected 49 species among four genera: Indonemoura Baumann, 1975; Protonemura Kempny, 1898; Amphinemura, Ris 1902 and Nemoura Latreille, 1796 in Japan, China, South Korea and North America. We estimated their divergence times-based on three molecular clock node calibrations-using Bayesian phylogeography approaches. Our results suggested that Japanese Archipelago formation events resulted in diversification events in the middle of the Cretaceous (<120 Ma), speciation in the Paleogene (<50 Ma) and intra-species diversification segregated into eastern and western Japan of the Fossa Magna region at late Neogene (20 Ma). The Indonemoura samples were genetically separated into two clades-that of Mainland China and that of Japan. The Japanese clade clustered with the Nemouridae species from North America, suggesting the possibility of a colonisation event prior to the formation of the Japanese Archipelago. We believe that our results enhanced the understanding both of the origin of the species and of local species distribution in the Japanese Archipelago.

Genotyping-by-sequencing supports a genetic basis for wing reduction in an alpine New Zealand stonefly

Wing polymorphism is a prominent feature of numerous insect groups, but the genomic basis for this diversity remains poorly understood. Wing reduction is a commonly observed trait in many species of stoneflies, particularly in cold or alpine environments. The widespread New Zealand stonefly Zelandoperla fenestrata species group (Z. fenestrata, Z. tillyardi, Z. pennulata) contains populations ranging from fully winged (macropterous) to vestigial-winged (micropterous), with the latter phenotype typically associated with high altitudes. The presence of flightless forms on numerous mountain ranges, separated by lowland fully winged populations, suggests wing reduction has occurred multiple times. We use Genotyping by Sequencing (GBS) to test for genetic differentiation between fully winged (n = 62) and vestigial-winged (n = 34) individuals, sampled from a sympatric population of distinct wing morphotypes, to test for a genetic basis for wing morphology. While we found no population genetic differentiation between these two morphotypes across 6,843 SNP loci, we did detect several outlier loci that strongly differentiated morphotypes across independent tests. These findings indicate that small regions of the genome are likely to be highly differentiated between morphotypes, suggesting a genetic basis for wing reduction. Our results provide a clear basis for ongoing genomic analysis to elucidate critical regulatory pathways for wing development in Pterygota.

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.

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).

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.

Diversification shifts in leafroller moths linked to continental colonization and the rise of angiosperms

Tectonic dynamics and niche availability play intertwined roles in determining patterns of diversification. Such drivers explain the current distribution of many clades, whereas events such as the rise of angiosperms can have more specific impacts, such as on the diversification rates of herbivores. The Tortricidae, a diverse group of phytophagous moths, are ideal for testing the effects of these determinants on the diversification of herbivorous clades. To estimate ancestral areas and diversification patterns in Tortricidae, a complete tribal-level dated tree was inferred using molecular markers (one mitochondrial and five nuclear) and calibrated using fossil constraints. We found that Tortricidae diverged from their sister group c. 120 Myr ago (Ma) and diversified c. 97 Ma, a timeframe synchronous with the rise of angiosperms in the Early-mid Cretaceous. Ancestral areas analysis, based on updated Wallace's biogeographical regions, supports the hypothesis of a Gondwanan origin of Tortricidae in the South American plate. We also detected an increase in speciation rate that coincided with the peak of angiosperm diversification in the Cretaceous. This in turn probably was further heightened by continental colonization of the Palaeotropics when angiosperms became dominant by the end of the Late Cretaceous.

Phylogenomic analyses of Crassiclitellata support major Northern and Southern Hemisphere clades and a Pangaean origin for earthworms

BACKGROUND

Earthworms (Crassiclitellata) are a diverse group of annelids of substantial ecological and economic importance. Earthworms are primarily terrestrial infaunal animals, and as such are probably rather poor natural dispersers. Therefore, the near global distribution of earthworms reflects an old and likely complex evolutionary history. Despite a long-standing interest in Crassiclitellata, relationships among and within major clades remain unresolved.

METHODS

In this study, we evaluate crassiclitellate phylogenetic relationships using 38 new transcriptomes in combination with publicly available transcriptome data. Our data include representatives of nearly all extant earthworm families and a representative of Moniligastridae, another terrestrial annelid group thought to be closely related to Crassiclitellata. We use a series of differentially filtered data matrices and analyses to examine the effects of data partitioning, missing data, compositional and branch-length heterogeneity, and outgroup inclusion.

RESULTS AND DISCUSSION

We recover a consistent, strongly supported ingroup topology irrespective of differences in methodology. The topology supports two major earthworm clades, each of which consists of a Northern Hemisphere subclade and a Southern Hemisphere subclade. Divergence time analysis results are concordant with the hypothesis that these north-south splits are the result of the breakup of the supercontinent Pangaea.

CONCLUSIONS

These results support several recently proposed revisions to the classical understanding of earthworm phylogeny, reveal two major clades that seem to reflect Pangaean distributions, and raise new questions about earthworm evolutionary relationships.

Cretaceous origin of giant rhinoceros beetles (Dynastini; Coleoptera) and correlation of their evolution with the Pangean breakup

The giant rhinoceros beetles (Dynastini, Scarabaeidae, Coleoptera) are distributed in tropical and temperate regions in Asia, America and Africa. Recent molecular phylogenetic studies have revealed that the giant rhinoceros beetles can be divided into three clades representing Asia, America and Africa. Although a correlation between their evolution and the continental drift during the Pangean breakup was suggested, there is no accurate divergence time estimation among the three clades based on molecular data. Moreover, there is a long chronological gap between the timing of the Pangean breakup (Cretaceous: 110-148 Ma) and the emergence of the oldest fossil record (Oligocene: 33 Ma). In this study, we estimated their divergence times based on molecular data, using several combinations of fossil calibration sets, and obtained robust estimates. The inter-continental divergence events among the clades were estimated to have occurred about 99 Ma (Asian clade and others) and 78 Ma (American clade and African clade), both of which are after the Pangean breakup. These estimates suggest their inter-continental divergences occurred by overseas sweepstakes dispersal, rather than by vicariances of the population caused by the Pangean breakup.

Incorporating fossils into hypotheses of insect phylogeny

Fossils represent stem and crown lineages, and their inclusion in phylogenetic reconstruction influences branch lengths, topology, and divergence time estimation. In addition, paleontological data may inform trends in morphological evolution as well as biogeographic history. Here we review the incorporation of fossils in studies of insect evolution, from morphological analyses to combined 'total evidence' node dating analyses. We discuss challenges associated with fossil based phylogenetics, and suggest best practices for use in tree reconstruction.

The first fossil salmonfly (Insecta: Plecoptera: Pteronarcyidae), back to the Middle Jurassic

BACKGROUND

The fossil record of Plecoptera (stoneflies) is considered relatively complete, with stem-groups of each of the three major lineages, viz. Antarctoperlaria, Euholognatha and Systellognatha (and some of their families) represented in the Mesozoic. However, the family Pteronarcyidae (the salmonflies; including two genera, Pteronarcys and Pteronarcella) has no fossil record to date, and the family has been suggested to have diverged recently.

RESULTS

In this paper, we report on a set of specimens belonging to a new fossil species of stonefly, discovered from the Middle Jurassic Daohugou locality (China). Our comparative analysis of wing venation and body characters demonstrates that the new species belongs to the Pteronarcyidae, and is more closely related to Pteronarcys than to Pteronarcella. However, it differs from all known species of the former genus. It is therefore assigned to a new genus and named Pteroliriope sinitshenkovae gen. et sp. nov. under the traditional nomenclatural procedure. The cladotypic nomenclatural procedure is also employed, with the resulting combination Pteroliriope nec Pteronarcys sinitshenkovae sp. nov.

CONCLUSIONS

The first discovery of a fossil member of the Pteronarcyidae demonstrates that the corresponding lineage is not a very recent offshoot but was already present ca. 165 million years ago. This discovery concurs with the view that divergence of most stonefly families took place very early, probably in the Triassic, or even in the Permian. This contribution demonstrates the need for (re-)investigations of the systematics of fossil stoneflies to refine divergence date estimates for Plecoptera lineages.