The effect of habitat type on speciation rates and range movements in aquatic beetles: inferences from species-level phylogenies

Spatio-temporal variation in water beetle assemblages across temperate freshwater ecosystems

Ecological communities are structured by several mechanisms, including temporal, spatial and environmental factors. However, the simultaneous effects of these factors have rarely been studied. Here, we investigated their role on water beetle assemblages sampled over a period of 18 years. Water beetles were sampled in the spring of each year in lotic and lentic water bodies from mainland region of Kalmar and Öland Island in southeastern Sweden. We assessed how past assemblage structure, environmental factors and spatial variables correlated with current assemblage structure using a variation partitioning approach. We also tested for correlates of temporal beta diversity of water beetle assemblages with multiple regressions. We found that past water beetle assemblage structure explained current water beetle assemblage structure better than the environmental and spatial correlates. We also observed that temporal beta diversity of water beetle assemblages was mainly due to species gain rather than to species loss. Finally, environmental variables (e.g., hydroperiod, habitat size and hydrology) and timespan between sampling events explained part of temporal beta diversity and contribution of species loss to total assemblage dissimilarity variation. Despite the fact that most variation remained unexplained, we found that ecological factors that have been thought to be important for water beetle richness and abundance in past studies (e.g. water body size, water permanence, shore slope, and whether the water body is lentic or lotic) were also correlated to temporal beta diversity. From a conservation point of view, our study suggest that temporal variability of assemblage structure should be included in biological monitoring because of its potential to predict current assemblage structure.

Habitat preference and diversification rates in a speciose lineage of diving beetles

The long-term geological stability of aquatic habitats has been demonstrated to be a determinant in the evolution of macroinvertebrate fauna, with species in running (lotic) waters having lower dispersal abilities, smaller ranges and higher gene flow between populations than species in standing (lentic) environments. Lotic species have been hypothesized to be more specialised, but the diversification dynamics of both habitat types have not been studied in detail. Using a speciose lineage of water beetles we test here whether diversification rates are related to the habitat preference of the species and its consequences on turnover, which we expect to be higher for lotic taxa. Moreover, we tested whether life in lotic environments is acting as an evolutionary dead-end as it is considered an ecological specialisation. We built a comprehensive molecular phylogeny with 473 terminals representing 421 of the 689 known species of the tribe Hydroporini (Coleoptera, Dytiscidae), using a combination of sequences from four mitochondrial and two nuclear genes plus 69 mitogenomes obtained with NGS. We found a general pattern of gradual acceleration of diversification rate with time, with 2-3 significant diversification shifts. However, habitat is not the main factor driving diversification in Hydroporini based on SecSSE analyses. The most recent common ancestor of Hydroporini was reconstructed as a lotic species, with multiple shifts to lentic environments. Most frequent transitions were estimated from lentic and lotic habitats to the category "both", followed by transitions from lotic to lentic and lentic to lotic respectively, although with very similar rates. Contrary to expectations, we found little evidence for differences in diversification dynamics between habitats, with lotic environments clearly not acting as evolutionary dead-ends in Hydroporini.

Going underwater: multiple origins and functional morphology of piercing-sucking feeding and tracheal system adaptations in water scavenger beetle larvae (Coleoptera: Hydrophiloidea)

Larvae of water scavenger beetles (Coleoptera: Hydrophiloidea) are adapted to a wide variety of aquatic habitats, but little is known about functional and evolutionary aspects of these adaptations. We review the functional morphology and evolution of feeding strategies of larvae of the families Hydrophilidae and Epimetopidae based on a detailed scanning electron microscope (SEM) analysis, analysis of video records of feeding behaviour and observations of living larvae. There are two main types of feeding mechanisms: chewing and piercing-sucking. The character mapping using the latest phylogenetic hypothesis for Hydrophiloidea infers the chewing system as the ancestral condition. The piercing-sucking mechanism evolved at least four times independently: once in Epimetopidae (Epimetopus) and three times in Hydrophilidae (Berosini: Berosus + Hemiosus; Laccobiini: Laccobius group; Hydrobiusini: Hybogralius). The piercing-sucking apparatus allows underwater extra-oral digestion and decreases the dependence of larvae on an aerial environment. A detailed study of the tracheal morphology of the piercing-sucking lineages reveals four independent origins of the apneustic respiratory system, all of them nested within lineages with piercing-sucking mouthparts. We conclude that piercing-sucking mouthparts represent a key innovation, which allows for the subsequent adaptation of the tracheal system, influences the diversification dynamics of the lineages and allows the shift to new adaptive zones.

Drivers of Benthic Macroinvertebrate Assemblages in Equatorial Alpine Rivers of the Rwenzoris (Uganda)

The Sub-Saharan alpine freshwater biodiversity is currently impacted by human settlements, climate change, agriculture, and mining activities. Because of the limited biodiversity studies in the region, a better understanding is needed of the important environmental variables affecting macroinvertebrate assemblages. In this paper, macroinvertebrate diversity responses to 18 environmental variables were studied at 30 sites along unique Rwenzori rivers at the equator in Uganda. We hypothesized that anthropogenic disturbance and local environmental variables affect macroinvertebrate diversity, irrespective of altitudinal gradients. Based on altitude and climate, the sites were subdivided into three altitude groups consisting of 10 sites each: upstream (US) 1400–1600 m.a.s.l.; midstream (MS) 1091–1399 m.a.s.l., and downstream (DS) 900–1090 m.a.s.l. A total of 44 macroinvertebrate families and 1623 individuals were identified. The macroinvertebrate diversity patterns were influenced by temperature, altitude, and latitude. Regression analysis revealed that temperature and nickel, were negative predictors of taxa richness. Nickel, which is released by mining activity, is detrimental to aquatic communities in Sub-Saharan alpine ecosystems. Significant longitudinal variation in macroinvertebrate diversity was observed between the sites, which were also affected by mineral and temperature gradients. Our study highlights the need for long-term monitoring in this region to detect and reduce the threats to river biodiversity from anthropogenic activity.

The tremendous diversity of Labiobaetis Novikova & Kluge in Indonesia (Ephemeroptera, Baetidae)

Material collected between 2010 and 2014 on the Indonesian islands of Sumatra, Bali, Sumba, Sumbawa, Sulawesi, and Seram unveiled the enormous diversity of Labiobaetis Novikova & Kluge in this country. Five species were reported from Indonesia previously (L. fulmeki (Ulmer), L. obscurum (Ulmer), L. necopinatum (Müller-Liebenau), L. ulmeri (Müller-Liebenau), and L. boettgeri (Ulmer)); all were described from adults only and no species were previously known at larval stage. We identified 18 new species by integrative taxonomy using genetic distance (COI, Kimura-2-parameter) and morphology, and they are described and illustrated based on their larvae. Another species, L. multus (Müller-Liebenau) from Malaysia, was also found in Indonesia, increasing the total number of species in Indonesia to 24. Seven morpho-groups of species are proposed based on morphological characters and a key to the larvae of all species from Indonesia and adjacent countries is provided. The total number of Labiobaetis species worldwide is augmented to 123. The examination of the new species allowed us to slightly modify the generic attributes of the larvae. The interspecific K2P distances are usually between 11% and 24%, the intraspecific distances are usually between 0% and 3%. The remarkable richness of the genus in Indonesia is discussed.

Dispersal limitation: Evolutionary origins and consequences in arthropods

Niche and dispersal ability are key traits for explaining the geographical structuring of species into discrete populations, and its evolutionary significance. Beyond their individual effects, the interplay between species niche and its geographic limits, together with the evolutionary lability of dispersal ability, can underpin trait diversification and speciation when exposed to gradients of selection. In this issue of Molecular Ecology, two complementary papers demonstrate how evolutionary lability for dispersal ability linked to niche shift can drive such a model in a context that includes selection. Both papers investigate the evolution of dispersal limitation in arthropods across altitudinal gradients, but using taxa with contrasting ecologies. McCulloch et al. (2019) investigate the evolution of wing loss at higher altitudes in stoneflies, a taxon inhabiting freshwater systems. Suzuki et al. (2019) report a similar phenomenon, but involving wing reduction at higher altitudes in scorpionflies, a taxon associated with moist terrestrial habitats. Here, we compare and contrast the results of both studies to explore their broader implications for understanding diversification and speciation within arthropods.

Water Beetles as Models in Ecology and Evolution

Beetles have colonized water many times during their history, with some of these events involving extensive evolutionary radiations and multiple transitions between land and water. With over 13,000 described species, they are one of the most diverse macroinvertebrate groups in most nonmarine aquatic habitats and occur on all continents except Antarctica. A combination of wide geographical and ecological range and relatively accessible taxonomy makes these insects an excellent model system for addressing a variety of questions in ecology and evolution. Work on water beetles has recently made important contributions to fields as diverse as DNA taxonomy, macroecology, historical biogeography, sexual selection, and conservation biology, as well as predicting organismal responses to global change. Aquatic beetles have some of the best resolved phylogenies of any comparably diverse insect group, and this, coupled with recent advances in taxonomic and ecological knowledge, is likely to drive an expansion of studies in the future.

The incredible diversity of Labiobaetis Novikova & Kluge in New Guinea revealed by integrative taxonomy (Ephemeroptera, Baetidae)

Material collected between 1999 and 2011 in Papua New Guinea and the Papua Province of Indonesia unveiled the enormous diversity of Labiobaetis on this island. Twenty-six new species were identified and delimited by integrative taxonomy using genetic distance (COI, Kimura-2-parameter) and morphology. These new species are described and illustrated based on larvae, augmenting the total number of Labiobaetis species on the island of New Guinea to 32. Seven morpho-groups of species are proposed based on morphological characters and a key to all New Guinea species is provided. The generic attributes of the larvae are summarised and slightly modified based on the examinations of the new species. Results on the genetics of most species (COI) are also provided. The interspecific K2P distances are between 13% and 32%, the intraspecific distances usually between 0% and 2%. Possible reasons for the remarkable richness of this genus in New Guinea are discussed.

Two new species of the megadiverse lentic diving-beetle genus Hydrovatus (Coleoptera, Dytiscidae) described from NE Thailand

Here we describe two new Hydrovatus species (Coleoptera: Dytiscidae: Hydroporinae: Hydrovatini) from the province of Khon Kaen, Isan region in NE Thailand. Hydrovatus is the third most species rich genus of diving beetles (Dytiscidae). It occurs on all continents except Antarctica and now numbers 210 currently recognized species. Both new species, Hydrovatusdiversipunctatussp. n. and Hydrovatusglobosussp. n., were collected at lights and are only known from the type locality “Khon Kaen” (a city and province). Diagnoses based on morphology for the separation from closely related species are given together with illustrations of male genitalia and habitus photos. We provide a determination key to Old World species of the pustulatus species group and to Oriental species of the oblongipennis species group.

Rates and patterns of molecular evolution in freshwater versus terrestrial insects

Insect lineages have crossed between terrestrial and aquatic habitats many times, for both immature and adult life stages. We explore patterns in molecular evolutionary rates between 42 sister pairs of related terrestrial and freshwater insect clades using publicly available protein-coding DNA sequence data from the orders Coleoptera, Diptera, Lepidoptera, Hemiptera, Mecoptera, Trichoptera, and Neuroptera. We furthermore test for habitat-associated convergent molecular evolution in the cytochrome c oxidase subunit I (COI) gene in general and at a particular amino acid site previously reported to exhibit habitat-linked convergence within an aquatic beetle group. While ratios of nonsynonymous-to-synonymous substitutions across available loci were higher in terrestrial than freshwater-associated taxa in 26 of 42 lineage pairs, a stronger trend was observed (20 of 31, p_binomial = 0.15, p_Wilcoxon = 0.017) when examining only terrestrial-aquatic pairs including fully aquatic taxa. We did not observe any widespread changes at particular amino acid sites in COI associated with habitat shifts, although there may be general differences in selection regime linked to habitat.

The origin of widespread species in a poor dispersing lineage (diving beetle genus Deronectes)

In most lineages, most species have restricted geographic ranges, with only few reaching widespread distributions. How these widespread species reached their current ranges is an intriguing biogeographic and evolutionary question, especially in groups known to be poor dispersers. We reconstructed the biogeographic and temporal origin of the widespread species in a lineage with particularly poor dispersal capabilities, the diving beetle genus Deronectes (Dytiscidae). Most of the ca. 60 described species of Deronectes have narrow ranges in the Mediterranean area, with only four species with widespread European distributions. We sequenced four mitochondrial and two nuclear genes of 297 specimens of 109 different populations covering the entire distribution of the four lineages of Deronectes, including widespread species. Using Bayesian probabilities with an a priori evolutionary rate, we performed (1) a global phylogeny/phylogeography to estimate the relationships of the main lineages within each group and root them, and (2) demographic analyses of the best population coalescent model for each species group, including a reconstruction of the geographical history estimated from the distribution of the sampled localities. We also selected 56 specimens to test for the presence of Wolbachia, a maternally transmitted parasite that can alter the patterns of mtDNA variability. All species of the four studied groups originated in the southern Mediterranean peninsulas and were estimated to be of Pleistocene origin. In three of the four widespread species, the central and northern European populations were nested within those in the northern areas of the Anatolian, Balkan and Iberian peninsulas respectively, suggesting a range expansion at the edge of the southern refugia. In the Mediterranean peninsulas the widespread European species were replaced by vicariant taxa of recent origin. The fourth species (D. moestus) was proven to be a composite of unrecognised lineages with more restricted distributions around the Western and Central Mediterranean. The analysis of Wolbachia showed a high prevalence of infection among Deronectes, especially in the D. aubei group, where all sequenced populations were infected with the only exception of the Cantabrian Mountains, the westernmost area of distribution of the lineage. In this group there was a phylogenetic incongruence between the mitochondrial and the nuclear sequence, although no clear pattern links this discordance to the Wolbachia infection. Our results suggest that, in different glacial cycles, populations that happened to be at the edge of the newly deglaciated areas took advantage of the optimal ecological conditions to expand their ranges to central and northern Europe. Once this favourable ecological window ended populations become isolated, resulting in the presence of closely related but distinct species in the Mediterranean peninsulas.

Phylogeography of Cephalotaxus oliveri (Cephalotaxaceae) in relation to habitat heterogeneity, physical barriers and the uplift of the Yungui Plateau

Habitat heterogeneity, physical barriers, and the uplift of the Yungui Plateau were found to deeply affect the phylogeographic pattern and evolutionary history of Cephalotaxus oliveri, a perennial conifer endemic to China. In this study, we explored the phylogeography using three chloroplast sequences (trnL-trnF, trnT-trnD and atpB-rbcL) in 22 natural populations of C. oliveri distributed throughout its range. The Yungui Plateau populations of C. oliveri were revealed to origin ca. 9.15Ma by molecular clock estimation, which is consistent with rapid uplift of the Qinghai-Tibetan Plateau (QTP) ca. 8-10Ma. Additionally, geological effects of the Yungui Plateau were suggested to promote the rapid intra-specific differentiation of C. oliveri in the Pliocene and Early Pleistocene. The relatively low level of genetic diversity (h=0.719, θ=1.17×10(-3)) and high population differentiation (NST=0.771 and GST=0.642) implied restricted gene flow among populations, which was confirmed by the Nested Clade Analysis (NCA). Mismatch distribution and haplotypes network provided evidences of recent demographic population expansion. Furthermore, the statistical dispersal-vicariance analysis indicated that the center of origin was in Central China. The comparison of haplotype distribution patterns indicated that the regions of HNHPS and HBLD were the potential refugia during the Pleistocene ice ages. Our results highlighted that habitat heterogeneity and physical barriers presenting in a species range can predict genetic patterns.

Clade age and diversification rate variation explain disparity in species richness among water scavenger beetle (Hydrophilidae) lineages

Explaining the disparity of species richness across the tree of life is one of the great challenges in evolutionary biology. Some lineages are exceptionally species rich, while others are relatively species poor. One explanation for heterogeneity among clade richness is that older clades are more species rich because they have had more time to accrue diversity than younger clades. Alternatively, disparity in species richness may be due to among-lineage diversification rate variation. Here we investigate diversification in water scavenger beetles (Hydrophilidae), which vary in species richness among major lineages by as much as 20 fold. Using a time-calibrated phylogeny and comparative methods, we test for a relationship between clade age and species richness and for shifts in diversification rate in hydrophilids. We detected a single diversification rate increase in Megasternini, a relatively young and species rich clade whose diversity might be explained by the stunning diversity of ecological niches occupied by this clade. We find that Amphiopini, an old clade, is significantly more species poor than expected, possibly due to its restricted geographic range. The remaining lineages show a correlation between species richness and clade age, suggesting that both clade age and variation in diversification rates explain the disparity in species richness in hydrophilids. We find little evidence that transitions between aquatic, semiaquatic, and terrestrial habitats are linked to shifts in diversification rates.

The secondary contact zone of phylogenetic lineages of the Philaenus spumarius (Hemiptera: Aphrophoridae): an example of incomplete allopatric speciation

Previous studies on the phylogeography of the meadow spittlebug Philaenus spumarius (L.) (Hemiptera: Aphrophoridae) suggest the existence of a contact zone of its main phylogenetic lineages along mountain chains in Europe and western Asia. This study presents a detailed examination of the population genetics of P. spumarius within the Carpathian Mountains. The main objective was to determine whether the populations inhabiting that area consist of individuals belonging to different genetic units and whether the observed pattern could be an example of secondary contact zone which formed after incomplete allopatric speciation. Specimens from six transects across the Carpathian arc were examined. The mitochondrial phylogeography of the meadow spittlebug in the examined area clearly shows that individuals from both main clades meet and mix there. Representatives of all three main EF1-α clades were also found. The present distribution of the main clades with a zone of overlap along the mountain ranges may suggest that these phylogenetic lineages form a young hybrid zone. Moreover, a limited number of individuals were shown to possess heteroplasmic mitochondrial DNA, which gives additional support to intraspecific hybridization. P. spumarius could be used in future work as an excellent model species in investigating population genetics, intraspecific hybridization, and speciation in progress.

Freshwater biodiversity and aquatic insect diversification

Inland waters cover less than 1% of Earth's surface but harbor more than 6% of all insect species: Nearly 100,000 species from 12 orders spend one or more life stages in freshwater. Little is known about how this remarkable diversity arose, although allopatric speciation and ecological adaptation are thought to be primary mechanisms. Freshwater habitats are highly susceptible to environmental change and exhibit marked ecological gradients. Standing waters appear to harbor more dispersive species than running waters, but there is little understanding of how this fundamental ecological difference has affected diversification. In contrast to the lack of evolutionary studies, the ecology and habitat preferences of aquatic insects have been intensively studied, in part because of their widespread use as bioindicators. The combination of phylogenetics with the extensive ecological data provides a promising avenue for future research, making aquatic insects highly suitable models for the study of ecological diversification.

Loss of flight promotes beetle diversification

The evolution of flight is a key innovation that may enable the extreme diversification of insects. Nonetheless, many species-rich, winged insect groups contain flightless lineages. The loss of flight may promote allopatric differentiation due to limited dispersal power and may result in a high speciation rate in the flightless lineage. Here we show that loss of flight accelerates allopatric speciation using carrion beetles (Coleoptera: Silphidae). We demonstrate that flightless species retain higher genetic differentiation among populations and comprise a higher number of genetically distinct lineages than flight-capable species, and that the speciation rate with the flightless state is twice that with the flight-capable state. Moreover, a meta-analysis of 51 beetle species from 15 families reveals higher genetic differentiation among populations in flightless compared with flight-capable species. In beetles, which represent almost one-fourth of all described species, repeated evolution of flightlessness may have contributed to their steady diversification since the Mesozoic era.

The loss of flight in some insect lineages may promote allopatric differentiation and result in a high speciation rate. Here, using the carrion beetle, loss of flight is shown to accelerate allopatric speciation with higher genetic differentiation than for flight-capable species.

Phylogenetic inference of reciprocal effects between geographic range evolution and diversification

Geographic characters--traits describing the spatial distribution of a species-may both affect and be affected by processes associated with lineage birth and death. This is potentially confounding to comparative analyses of species distributions because current models do not allow reciprocal interactions between the evolution of ranges and the growth of phylogenetic trees. Here, we introduce a likelihood-based approach to estimating region-dependent rates of speciation, extinction, and range evolution from a phylogeny, using a new model in which these processes are interdependent. We demonstrate the method with simulation tests that accurately recover parameters relating to the mode of speciation and source-sink dynamics. We then apply it to the evolution of habitat occupancy in Californian plant communities, where we find higher rates of speciation in chaparral than in forests and evidence for expanding habitat tolerances.

Red drifters and dark residents: the phylogeny and ecology of a Plio-Pleistocene dragonfly radiation reflects Africa's changing environment (Odonata, Libellulidae, Trithemis)

In the last few million years, tropical Africa has experienced pronounced climatic shifts with progressive aridification. Such changes must have had a great impact on freshwater biota, such as Odonata. With about forty species, Trithemis dominates dragonfly communities across Africa, from rain-pools to streams, deserts to rainforests, and lowlands to highlands. Red-bodied species tend to favor exposed, standing and often temporary waters, have strong dispersal capacities, and some of the largest geographic ranges in the genus. Those in cooler habitats, like forest streams, are generally dark-bodied and more sedentary. We combined molecular analyses of ND1, 16S, and ITS (ITSI, 5.8S, and ITSII) with morphological, ecological, and geographical data for 81% of known Trithemis species, including three Asian and two Madagascan endemics. Using molecular clock analyses, the genus's origin was estimated 6-9Mya, with multiple lineages arising suddenly around 4Mya. Open stagnant habitats were inferred to be ancestral and the rise of Trithemis may have coincided with savannah-expansion in the late Miocene. The adaptation of red species to more ephemeral conditions leads to large ranges and limited radiation within those lineages. By contrast, three clades of dark species radiated in the Plio-Pleistocene, each within distinct ecological confines: (1) lowland streams, (2) highland streams, and (3) swampy habitats on alternating sides of the Congo-Zambezi watershed divide; together giving rise to the majority of species diversity in the genus. During Trithemis evolution, multiple shifts from open to more forested habitats and from standing to running waters occurred. Allopatry by habitat fragmentation may be the dominant force in speciation, but possibly genetic divergence across habitat gradients was also involved. The study demonstrates the importance of combining ecological and phylogenetic data to understand the origin of biological diversity under great environmental change.

Sympatric speciation in birds is rare: insights from range data and simulations

Sympatric speciation is now accepted as theoretically plausible and a likely explanation for divergence in a handful of taxa, but its contribution to large-scale patterns of speciation remains contentious. A major problem is that it is difficult to differentiate between alternate scenarios of geographic speciation when species ranges have shifted substantially in the past. Previous studies have searched for a signal of the geographic mode of speciation by testing for a correlation between time since speciation and range overlap. Here we use simulations to show that the proportion of species showing zero or complete range overlap are more reliable indicators of the geography of speciation than is the correlation between time since speciation and overlap. We then apply these findings to the distributions of 291 pairs of avian sister species. Although 49% of pairs show some overlap in their ranges, our simulations show that this is not surprising under allopatric models of speciation. More revealingly, less than 2% show complete range overlap. Our simulations demonstrate that the observed patterns are most consistent with a model in which allopatric speciation is dominant but in which sympatric speciation is also present and contributes 5% of speciation events.

Phylogenetic relationships of glassfrogs (Centrolenidae) based on mitochondrial and nuclear genes

Glassfrogs (family Centrolenidae) represent an exceptionally diverse group among Neotropical anurans, but their evolutionary relationships never have been assessed from a molecular perspective. Mitochondrial and nuclear markers were used to develop a novel hypothesis of centrolenid phylogeny. Ingroup sampling included 100 terminals, with 78 (53%) of the named species in the family, representing most of the phenotypic diversity described for the group. Thirty-five species representing taxa traditionally associated with glassfrogs were used as outgroups. Gene sampling consisted of complete or partial sequences of three mitochondrial (12S, 16S, ND1) and three nuclear markers (c-myc exon 2, RAG1, POMC) for a total of approximately 4362bp. Phylogenies were estimated using maximum parsimony, maximum likelihood, and Bayesian analyses for individual genes and combined datasets. The separate analysis of mitochondrial and nuclear datasets allowed us to clarify the relationships within glassfrogs; also, we corroborate the sister-group relationship between Allophryne ruthveni and glassfrogs. The new phylogeny differs significantly from all previous morphology-based hypotheses of relationships, and shows that hypotheses based on few traits are likely to misrepresent evolutionary history. Traits previously hypothesized as unambiguous synapomorphies are shown to be homoplastic, and all genera in the current taxonomy (Centrolene, Cochranella, Hyalinobatrachium, Nymphargus) are found to be poly- or paraphyletic. The new topology implies a South American origin of glassfrogs and reveals allopatric speciation as the most important speciation mechanism. The phylogeny profoundly affects the traditional interpretations of glassfrog taxonomy, character evolution, and biogeography-topics that now require more extensive evaluation in future studies.

Thermal tolerance, acclimatory capacity and vulnerability to global climate change

Despite evidence that organismal distributions are shifting in response to recent climatic warming, we have little information on direct links between species' physiology and vulnerability to climate change. We demonstrate a positive relationship between upper thermal tolerance and its acclimatory ability in a well-defined clade of closely related European diving beetles. We predict that species with the lowest tolerance to high temperatures will be most at risk from the adverse effects of future warming, since they have both low absolute thermal tolerance and poor acclimatory ability. Upper thermal tolerance is also positively related to species' geographical range size, meaning that species most at risk are already the most geographically restricted ones, being endemic to Mediterranean mountain systems. Our findings on the relationship between tolerance and acclimatory ability contrast with results from marine animals, suggesting that generalizations regarding thermal tolerance and responses to future rapid climate change may be premature.

Estimating diversification rates from phylogenetic information

Patterns of species richness reflect the balance between speciation and extinction over the evolutionary history of life. These processes are influenced by the size and geographical complexity of regions, conditions of the environment, and attributes of individuals and species. Diversity within clades also depends on age and thus the time available for accumulating species. Estimating rates of diversification is key to understanding how these factors have shaped patterns of species richness. Several approaches to calculating both relative and absolute rates of speciation and extinction within clades are based on phylogenetic reconstructions of evolutionary relationships. As the size and quality of phylogenies increases, these approaches will find broader application. However, phylogeny reconstruction fosters a perceptual bias of continual increase in species richness, and the analysis of primarily large clades produces a data selection bias. Recognizing these biases will encourage the development of more realistic models of diversification and the regulation of species richness.

Why are there so many insect species? Perspectives from fossils and phylogenies

Over half of all described species are insects, but until recently our understanding of the reasons for this diversity was based on very little macroevolutionary evidence. Here I summarize the hypotheses that have been posed, tests of these hypotheses and their results, and hence identify gaps in knowledge for future researchers to pursue. I focus on inferences from the following sources: (i) the fossil record, normally at family level, and (ii) insect phylogenies, sometimes combined with: (iii) the species richness of insect higher taxa, and (iv) current extinction risks. There is evidence that the species richness of insects has been enhanced by: (i) their relative age, giving time for diversification to take place; (ii) low extinction rates. There is little evidence that rates of origination have generally been high or that there are limits on numbers of species. However, the evidence on macroevolutionary rates is not yet so extensive or coherent as to present unequivocal messages. As regards morphological, ecological, or behavioural hypotheses, there is evidence that diversity has been enhanced by (iii) flight or properties resulting from it like enhanced dispersal, (iv) wing folding, and (v) complete metamorphosis. However, in all these cases the evidence is somewhat equivocal, either because of statistical issues or because evidence from different sources is conflicting. There is extensive evidence that diversity is affected by (vi) the ecological niche. Comparative studies indicate that phytophagy generally increases net diversification rates, and reduces extinction risk. However, niche specialization is also associated with an increase in extinction risk. Small body size (vii) is often associated with low extinction risk in comparative studies, but as yet there is no solid evidence that it consistently enhances net rates of diversification. Mouthpart diversity (viii) has generally increased over time in the insects, but cannot explain the apparent great increase in diversity seen in the Cretaceous and Tertiary. Sexual selection and sexual conflict (ix) are two processes that are widespread in insects, and there is comparative evidence linking both to increased diversification. Although some comparative evidence links tropical distributions (x) to increased rates of diversification, the extent to which latitudinal richness gradients are unusual in insects is equivocal. There is little to no direct evidence from fossils and phylogenies that insect diversity has generally been affected by (i) sensory- or neuro-sophistication, (ii) population size or density, (iii) generation time or fecundity, (iv) the presence of an exoskeleton or cuticle, (v) segmentation or appendage diversity, (vi) adaptability or genetic versatility, though all of these remain plausible hypotheses awaiting further tests. The data suggest that the insect body ground plan itself had no direct effect on insect diversity. Thus, whilst studies to date have given substantial understanding, substantial gaps still remain. Future challenges include: (i) interpreting conflicting messages from different sources of data; (ii) rating the importance of different hypotheses that are statistically supported; (iii) linking specific proximate to specific ultimate explanations and vice versa; and (iv) understanding how different ultimate hypotheses might be dependent on each other.

Birth-Death Models in Macroevolution

Birth-death models, and their subsets—the pure birth and pure death models—have a long history of use for informing thinking about macroevolutionary patterns. Here we illustrate with examples the wide range of questions they have been used to address, including estimating and comparing rates of diversification of clades, investigating the “shapes” of clades, and some rather surprising uses such as estimating speciation rates from data that are not resolved below the level of the genus. The raw data for inference can be the fossil record or the molecular phylogeny of a clade, and we explore the similarites and differences in the behavior of the birth-death models when applied to these different forms of data.

Habitat type predicts genetic population differentiation in freshwater invertebrates

A basic challenge in evolutionary biology is to establish links between ecology and evolution of species. One important link is the habitat template. It has been hypothesized, that the spatial and temporal settings of a habitat strongly influence the evolution of species dispersal propensity. Here, we evaluate the importance of the habitat type on genetic population differentiation of species using freshwater habitats as a model system. Freshwater habitats are either lentic (standing) or lotic (running). On average, lotic habitats are more stable and predictable over space and time than lentic habitats. Therefore, lentic habitats should favour the evolution of higher dispersal propensity which ensures population survival of lentic species. To test for such a relationship, we used extensive data on species' genetic population differentiation of lentic and lotic freshwater invertebrates retrieved from published allozyme studies. Overall, we analysed more than 150 species from all over the world. Controlling for several experimental, biological and geographical confounding effects, we always found that lentic invertebrates exhibit, on average, lower genetic population differentiation than lotic species. This pattern was consistent across insects, crustaceans and molluscs. Our results imply fundamental differences in genetic population differentiation among species adapted to either lentic or lotic habitats. We propose that such differences should occur in a number of other habitat types that differ in spatio-temporal stability. Furthermore, our results highlight the important role of lotic habitats as reservoirs for evolutionary processes and the potential for rapid speciation.

Trans-oceanic and endemic origins of the small minnow mayflies (Ephemeroptera, Baetidae) of Madagascar

We investigated the relative importance of dispersal and vicariance in forming the Madagascar insect fauna, sequencing approximately 2300bp from three rRNA gene regions to investigate the phylogeny of Afrotropical small minnow mayflies (Ephemeroptera: Baetidae). Six lineages contained trans-oceanic sister taxa, and variation in genetic divergence between sister taxa revealed relationships that range from very recent dispersal to ancient vicariance. Dispersal was most recent and frequent in species that spend the larval stage in standing water, adding to evidence that these evolutionarily unstable habitats may select for ecological traits that increase dispersal in insects. Ancestral state likelihood analysis suggested at least one Afrotropical lineage had its origin in Madagascar, demonstrating that unidirectional dispersal from a continental source may be too simplistic. We conclude that the Malagasy mayfly fauna should be considered in a biogeographical context that extends beyond Madagascar itself, encompassing trans-oceanic dispersal within multiple lineages.

Speciation and diversity on tropical rocky shores: a global phylogeny of snails of the genus Echinolittorina

A phylogenetic approach to the origin and maintenance of species diversity ideally requires the sampling of all species within a clade, confirmation that they are evolutionarily distinct entities, and knowledge of their geographical distributions. In the marine tropics such studies have mostly been of fish and reef-associated organisms, usually with high dispersal. In contrast, snails of the genus Echinolittorina (Littorinidae) are restricted to rocky shores, have a four-week pelagic development (and recorded dispersal up to 1400 km), and show different evolutionary patterns. We present a complete molecular phylogeny of Echinolittorina, derived from Bayesian analysis of sequences from nuclear 28S rRNA and mitochondrial 12S rRNA and COI genes (nodal support indicated by posterior probabilities, maximum likelihood, and neighbor-joining bootstrap). This consists of 59 evolutionarily significant units (ESUs), including all 50 known taxonomic species. The 26 ESUs found in the Indo-West Pacific region form a single clade, whereas the eastern Pacific and Atlantic species are basal. The earliest fossil occurred in the Tethys during the middle Eocene and we suggest that the Indo-West Pacific clade has been isolated since closure of the Tethyan seaway in the early Miocene. The geographical distributions of all species (based on more than 3700 locality records) appear to be circumscribed by barriers of low temperature, unsuitable sedimentary habitat, stretches of open water exceeding about 1400 km, and differences in oceanographic conditions on the continuum between oceanic and continental. The geographical ranges of sister species show little or no overlap, indicating that the speciation mode is predominantly allopatric. Furthermore, range expansion following speciation appears to have been limited, because a high degree of allopatry is maintained through three to five branching points of the phylogeny. This may be explained by infrequent long-distance colonization, habitat specialization on the oceanic/continental gradient, and perhaps by interspecific competition. In the eastern Pacific plus Atlantic we identify five cases of divergence on either side of the Isthmus of Panama, but our estimates of their ages pre-date the emergence of the Isthmus. There are three examples of sister relationships between species in the western Atlantic and eastern Atlantic, all resulting from dispersal to the east. Within the Indo-West Pacific, we find no geographical pattern of speciation events; narrowly endemic species of recent origin are present in both peripheral and central parts of the region. Evidence from estimated divergence times of sister species, and from a plot of the number of lineages over time, suggest that there has been no acceleration of diversification during the glacio-eustatic cycles of the Plio-Pleistocene. In comparison with reefal organisms, species of Echinolittorina on rocky shores may be less susceptible to extinction or isolation during sea-level fluctuations. The species richness of Echinolittorina in the classical biogeographic provinces conforms to the common pattern of highest diversity (11 species) in the central "East Indies Triangle" of the Indo-West Pacific, with a subsidiary focus in the eastern Pacific and western Atlantic, and lowest diversity in the eastern Atlantic. The diversity focus in the East Indies Triangle is produced by a mosaic of restricted allopatric species and overlap of a few widespread ones, and is the result of habitat specialization rather than historical vicariance. This study emphasizes the plurality of biogeographic histories and speciation patterns in the marine tropics.

Phylogeny and historical biogeography of Agabinae diving beetles (Coleoptera) inferred from mitochondrial DNA sequences

The Agabinae, with more than 350 species, is one of the most diverse lineages of diving beetles (Dytiscidae). Using the mitochondrial genes 16S rRNA and cytochrome oxidase I we present a phylogenetic analysis based on 107 species drawn mostly from the four main Holarctic genera. Two of these genera (Ilybius and Ilybiosoma) are consistently recovered as monophyletic with strong support, Platambus is never recovered as monophyletic, and Agabus is found paraphyletic with respect to several of the species groups of Platambus. Basal relationships among the main lineages are poorly defined, although within each of them relationships are in general robust and very consistent across the parameter space, and in agreement with previous morphological analyses. In each of the two most diverse lineages (Ilybius and Agabus including part of Platambus) there is a basal split between Palearctic and Nearctic clades, estimated to have occurred in the late Eocene. The Palearctic clade in turn splits into a Western Palearctic clade and a clade containing mostly Eastern Palearctic species, and assumed to be ancestrally Eastern Palearctic but with numerous transitions to a Holarctic or Nearctic distribution. These results suggest an asymmetry in the colonization routes, as there are very few cases of transcontinental range expansions originating from the Nearctic or the Western Palearctic. According to standard clock estimates, we do not find any transcontinental shift during the Pliocene, but numerous speciation events within each of the continental or subcontinental regions.

Speciation of Iberian diving beetles in Pleistocene refugia (Coleoptera, Dytiscidae)

The Mediterranean basin is an area of high diversity and endemicity, but the age and origin of its fauna are still largely unknown. Here we use species-level phylogenies based on approximately 1300 base pairs of the genes 16S rRNA and cytochrome oxidase I to establish the relationships of 27 of the 34 endemic Iberian species of diving beetles in the family Dytiscidae, and to investigate their level of divergence. Using a molecular clock approach, 18-19 of these species were estimated to be of Pleistocene origin, with four to six of them from the Late Pleistocene ( approximately 100 000 years). A second, lower speciation frequency peak was assigned to Late Miocene or Early Pliocene. Analysis of the distributional ranges showed that endemic species placed in the tip nodes of the trees are significantly more likely to be allopatric with their sisters than endemic species at lower node levels. Allopatric sister species are also significantly younger than sympatric clades, in agreement with an allopatric mode of speciation and limited subsequent range movement. These results strongly suggest that for some taxa Iberian populations were isolated during the Pleistocene long enough to speciate, and apparently did not expand their ranges to recolonize areas north of the Pyrenees. This is in contradiction to observations from fossil beetles in areas further north, which document large range movements associated with the Pleistocene glacial cycles hypothesized to suppress population isolation and allopatric speciation.

Mitochondrial DNA phylogeography and population history of Meladema diving beetles on the Atlantic Islands and in the Mediterranean basin (Coleoptera, Dytiscidae)

The phylogeny and population history of Meladema diving beetles (Coleoptera, Dytiscidae) were examined using mitochondrial DNA sequence from 16S ribosomal RNA and cytochrome oxidase I genes in 51 individuals from 22 populations of the three extant species (M. imbricata endemic to the western Canary Islands, M. lanio endemic to Madeira and M. coriacea widespread in the Western Mediterranean and on the western Canaries), using a combination of phylogenetic and nested clade analyses. Four main lineages are observed within Meladema, representing the three recognized species plus Corsican populations of M. coriacea. Phylogenetic analyses demonstrate the sister relationship of the two Atlantic Island taxa, and suggest the possible paraphyly of M. coriacea. A molecular clock approach reveals that speciation within the genus occurred in the Early Pleistocene, indicating that the Atlantic Island endemics are not Tertiary relict taxa as had been proposed previously. Our results point to past population bottlenecks in all four lineages, with recent (Late-Middle Pleistocene) range expansion in non-Corsican M. coriacea and M. imbricata. Within the Canary Islands, M. imbricata seems to have independently colonized La Gomera and La Palma from Tenerife (although a colonization of La Palma from La Gomera cannot be discarded), and M. coriacea has independently colonized Tenerife and Gran Canaria from separate mainland lineages. In the Mediterranean basin this species apparently colonized Corsica on a single occasion, relatively early in its evolutionary history (Early Pleistocene), and has colonized Mallorca recently on multiple occasions. On the only island where M. coriacea and M. imbricata are broadly sympatric (Tenerife), we report evidence of bidirectional hybridization between the two species.

Recent Diversification Rates in North American Tiger Beetles Estimated from a Dated mtDNA Phylogenetic Tree

Species-level phylogenies derived from DNA sequence data provide a tool for estimating diversification rates and how these rates change over time, but to date there have been few empirical studies, particularly on insect groups. We use a densely sampled phylogenetic tree based on mitochondrial DNA to investigate diversification rates in the North American tiger beetles (genus Cicindela). Using node ages estimated from sequence data and calibrated by biogeographical evidence, we estimate an average per-lineage diversification rate of at least 0.22 +/- 0.08 species/Myr over the time interval since the most recent colonization that led to a radiation within the continent. In addition, we find evidence for a weak, recent increase in the net diversification rate. This is more consistent with a late Pleistocene increase in the speciation rate than with a constant rate of background extinction, but the results are sensitive to the dating method and taxon sampling. We discuss practical limitations to phylogenetic studies of diversification rates.