Ecology has contrasting effects on genetic variation within species versus rates of molecular evolution across species in water beetles

Genetic diversity varies with species traits and latitude in predatory soil arthropods (Myriapoda: Chilopoda)

Aim

To investigate the drivers of intra-specific genetic diversity in centipedes, a group of ancient predatory soil arthropods.

Location

Asia, Australasia and Europe.

Time Period

Present.

Major Taxa Studied

Centipedes (Class: Chilopoda).

Methods

We assembled a database of 1245 mitochondrial cytochrome c oxidase subunit I sequences representing 128 centipede species from all five orders of Chilopoda. This sequence dataset was used to estimate genetic diversity for centipede species and compare its distribution with estimates from other arthropod groups. We studied the variation in centipede genetic diversity with species traits and biogeography using a beta regression framework, controlling for the effect of shared evolutionary history within a family.

Results

A wide variation in genetic diversity across centipede species (0-0.1713) falls towards the higher end of values among arthropods. Overall, 27.57% of the variation in mitochondrial COI genetic diversity in centipedes was explained by a combination of predictors related to life history and biogeography. Genetic diversity decreased with body size and latitudinal position of sampled localities, was greater in species showing maternal care and increased with geographic distance among conspecifics.

Main Conclusions

Centipedes fall towards the higher end of genetic diversity among arthropods, which may be related to their long evolutionary history and low dispersal ability. In centipedes, the negative association of body size with genetic diversity may be mediated by its influence on local abundance or the influence of ecological strategy on long-term population history. Species with maternal care had higher genetic diversity, which goes against expectations and needs further scrutiny. Hemispheric differences in genetic diversity can be due to historic climatic stability and lower seasonality in the southern hemisphere. Overall, we find that despite the differences in mean genetic diversity among animals, similar processes related to life-history strategy and biogeography are associated with the variation within them.

From DNA barcodes to ecology: Meta-analysis of central European beetles reveal link with species ecology but also to data pattern and gaps

DNA barcoding has been used worldwide to identify biological specimens and to delimit species. It represents a cost-effective, fast, and efficient way to assess biodiversity with help of the public Barcode of Life Database (BOLD) accounting for more than 236,000 animal species and more than 10 million barcode sequences. Here, we performed a meta-analysis of available barcode data of central European Coleoptera to detect intraspecific genetic patterns among ecological groups in relation to geographic distance with the aim to investigate a possible link between infraspecific variation and species ecology. We collected information regarding feeding style, body size, as well as habitat and biotope preferences. Mantel tests and two variants of Procrustes analysis, both involving the Principal Coordinates Neighborhood Matrices (PCNM) approach, were applied on genetic and geographic distance matrices. However, significance levels were too low to further use the outcome for further trait investigation: these were in mean for all ecological guilds only 7.5, 9.4, or 15.6% for PCNM + PCA, NMDS + PCA, and Mantel test, respectively, or at best 28% for a single guild. Our study confirmed that certain ecological traits were associated with higher species diversity and foster stronger genetic differentiation. Results suggest that increased numbers of species, sampling localities, and specimens for a chosen area of interest may give new insights to explore barcode data and species ecology for the scope of conservation on a larger scale. We performed a meta-analysis of available barcode data of central European beetles to detect intraspecific genetic patterns among ecological groups in relation to geographic distance, regarding feeding style, body size, as well as habitat and biotope preferences. Our study confirmed that certain ecological traits were associated with higher species diversity and foster stronger genetic differentiation. However, significance levels were too low to further use the outcome for further trait investigation.

The Quality of Sequence Data Affects Biodiversity and Conservation Perspectives in the Neotropical Damselfly Megaloprepus caerulatus

Ideally, the footprint of the evolutionary history of a species is drawn from integrative studies including quantitative and qualitative taxonomy, biogeography, ecology, and molecular genetics. In today’s research, species delimitations and identification of conservation units is often accompanied by a set of—at minimum—two sequence markers appropriate for the systematic level under investigation. Two such studies re-evaluated the species status in the world’s largest Odonata, the Neotropical damselfly Megaloprepus caerulatus. The species status of the genus Megaloprepus has long been debated. Despite applying a highly similar set of sequence markers, the two studies reached different conclusions concerning species status and population genetic relationships. In this study, we took the unique opportunity to compare the two datasets and analyzed the reasons for those incongruences. The two DNA sequence markers used (16S rDNA and CO1) were re-aligned using a strict conservative approach and the analyses used in both studies were repeated. Going step by step back to the first line of data handling, we show that a high number of unresolved characters in the sequence alignments as well as internal gaps are responsible for the different outcomes in terms of species delimitations and population genetic relationships. Overall, this study shows that high quality raw sequence data are an indispensable requirement, not only in odonate research.

High resolution DNA barcode library for European butterflies reveals continental patterns of mitochondrial genetic diversity

The study of global biodiversity will greatly benefit from access to comprehensive DNA barcode libraries at continental scale, but such datasets are still very rare. Here, we assemble the first high-resolution reference library for European butterflies that provides 97% taxon coverage (459 species) and 22,306 COI sequences. We estimate that we captured 62% of the total haplotype diversity and show that most species possess a few very common haplotypes and many rare ones. Specimens in the dataset have an average 95.3% probability of being correctly identified. Mitochondrial diversity displayed elevated haplotype richness in southern European refugia, establishing the generality of this key biogeographic pattern for an entire taxonomic group. Fifteen percent of the species are involved in barcode sharing, but two thirds of these cases may reflect the need for further taxonomic research. This dataset provides a unique resource for conservation and for studying evolutionary processes, cryptic species, phylogeography, and ecology.

Two ways to be endemic. Alps and Apennines are different functional refugia during climatic cycles

Endemics co-occur because they evolved in situ and persist regionally or because they evolved ex situ and later dispersed to shared habitats, generating evolutionary or ecological endemicity centres, respectively. We investigate whether different endemicity centres can intertwine in the region ranging from Alps to Sicily, by studying their butterfly fauna. We gathered an extensive occurrence data set for butterflies of the study area (27,123 records, 269 species, in cells of 0.5 × 0.5 degrees of latitude-longitude). We applied molecular-based delimitation methods (GMYC model) to 26,557 cytochrome c oxidase subunit 1 (COI) sequences of Western Palearctic butterflies. We identified entities based on molecular delimitations and/or the checklist of European butterflies and objectively attributed occurrences to their most probable entity. We obtained a zoogeographic regionalisation based on the 69 endemics of the area. Using phylogenetic ANOVA we tested if endemics from different centres differ from each other and from nonendemics for key ecological traits and divergence time. Endemicity showed high incidence in the Alps and Southern Italy. The regionalisation separated the Alps from the Italian Peninsula and Sicily. The endemics of different centres showed a high turnover and differed in phylogenetic distances, phenology and distribution traits. Endemics are on average younger than nonendemics and the Peninsula-Sicily endemics also have lower variance in divergence than those from the Alps. The observed variation identifies Alpine endemics as paleoendemics, now occupying an ecological centre, and the Peninsula-Sicily ones as neoendemics, that diverged in the region since the Pleistocene. The results challenge the common view of the Alpine-Apennine area as a single "Italian refugium".

The Effects of Ecological Traits on the Rate of Molecular Evolution in Ray-Finned Fishes: A Multivariable Approach

Myriad environmental and biological traits have been investigated for their roles in influencing the rate of molecular evolution across various taxonomic groups. However, most studies have focused on a single trait, while controlling for additional factors in an informal way, generally by excluding taxa. This study utilized a dataset of cytochrome c oxidase subunit I (COI) barcode sequences from over 7000 ray-finned fish species to test the effects of 27 traits on molecular evolutionary rates. Environmental traits such as temperature were considered, as were traits associated with effective population size including body size and age at maturity. It was hypothesized that these traits would demonstrate significant correlations with substitution rate in a multivariable analysis due to their associations with mutation and fixation rates, respectively. A bioinformatics pipeline was developed to assemble and analyze sequence data retrieved from the Barcode of Life Data System (BOLD) and trait data obtained from FishBase. For use in phylogenetic regression analyses, a maximum likelihood tree was constructed from the COI sequence data using a multi-gene backbone constraint tree covering 71% of the species. A variable selection method that included both single- and multivariable analyses was used to identify traits that contribute to rate heterogeneity estimated from different codon positions. Our analyses revealed that molecular rates associated most significantly with latitude, body size, and habitat type. Overall, this study presents a novel and systematic approach for integrative data assembly and variable selection methodology in a phylogenetic framework.

How long is 3 km for a butterfly? Ecological constraints and functional traits explain high mitochondrial genetic diversity between Sicily and the Italian Peninsula

Populations inhabiting Mediterranean islands often show contrasting genetic lineages, even on islands that were connected to the mainland during glacial maxima. This pattern is generated by forces acting in historical and contemporary times. Understanding these phenomena requires comparative studies integrating genetic structure, functional traits and dispersal constraints. Using as a model the butterfly species living across the Messina strait (3 km wide) separating Sicily from the Italian Peninsula, we aimed to unravel the mechanisms limiting the dispersal of matrilines and generating genetic differentiation across a narrow sea strait. We analysed the mitochondrial COI gene of 84 butterfly species out of 90 documented in Sicily and compared them with populations from the neighbouring southern Italian Peninsula (1,398 sequences) and from the entire Palaearctic region (8,093 sequences). For each species, we regressed 13 functional traits and 2 ecological constraints to dispersal (winds experienced at the strait and climatic suitability) against genetic differentiation between Sicily and Italian Peninsula to understand the factors limiting dispersal. More than a third of the species showed different haplogroups across the strait and most of them also represented endemic haplogroups for this island. One fifth of Sicilian populations (and 32.3% of endemic lineages) had their closest relatives in distant areas, instead of the neighbouring Italian Peninsula, which suggests high relictuality. Haplotype diversity was significantly explained by the length of the flight period, an intrinsic phenology trait, while genetic differentiation was explained by both intrinsic traits (wingspan and degree of generalism) and contemporary local constraints (winds experienced at the strait and climatic suitability). A relatively narrow sea strait can produce considerable differentiation among butterfly matrilines and this phenomenon showed a largely deterministic fingerprint. Because of unfavourable winds, populations of the less dispersive Sicilian butterflies tended to differentiate into endemic variants or to maintain relict populations. Understanding these phenomena required the integration of DNA sequences, species traits and physical constraints for a large taxon at continental scale. Future studies may reveal if the patterns here shown for mitochondrial DNA are also reflected in the nuclear genome or, alternatively, are the product of limited female dispersal.

To remain or leave: Dispersal variation and its genetic consequences in benthic freshwater invertebrates

Variation in dispersal capacity may influence population genetic variation and relatedness of freshwater animals thus demonstrating how life-history traits influence patterns and processes that in turn influence biodiversity. The majority of studies have focused on the consequences of dispersal variation in taxa inhabiting riverine systems whose dendritic nature and upstream/downstream gradients facilitate characterizing populations along networks. We undertook extensive, large-scale investigations of the impacts of hydrological connectivity on population genetic variation in two freshwater bryozoan species whose dispersive propagules (statoblasts) are either attached to surfaces (Fredericella sultana) or are released as buoyant stages (Cristatella mucedo) and that live primarily in either lotic (F. sultana) or lentic environments (C. mucedo). Describing population genetic structure in multiple sites characterized by varying degrees of hydrological connectivity within each of three (or four) UK regions enabled us to test the following hypotheses: (1) genetic diversity and gene flow will be more influenced by hydrological connectivity in populations of C. mucedo (because F. sultana dispersal stages are retained); (2) populations of F. sultana will be characterized by greater genetic divergence than those of C. mucedo (reflecting their relative dispersal capacities); and (3) genetic variation will be greatest in F. sultana (reflecting a propensity for genetic divergence as a result of its low dispersal potential). We found that hydrological connectivity enhanced genetic diversity and gene flow among C. mucedo populations but not in F. sultana while higher overall measures of clonal diversity and greater genetic divergence characterized populations of F. sultana. We suggest that genetic divergence over time within F. sultana populations reflects a general constraint of releasing propagules that might eventually be swept to sea when taxa inhabit running waters. In contrast, taxa that primarily inhabit lakes and ponds may colonize across hydrologically connected regions, establishing genetically related populations. Our study contributes more nuanced views about drivers of population genetic structures in passively dispersing freshwater invertebrates as outlined by the Monopolization Hypothesis (Acta Oecologica, 23, 2002, 121) by highlighting how a range of demographic and evolutionary processes reflect life-history attributes of benthic colonial invertebrates (bryozoans) and cyclically parthenogenetic zooplankton. In addition, growing evidence that genetic divergence may commonly characterize populations of diverse groups of riverine taxa suggests that organisms inhabiting lotic systems may be particularly challenged by environmental change. Such change may predispose riverine populations to extinction as a result of genetic divergence combined with limited dispersal and gene flow.

OPEN RESEARCH BADGES

This article has earned an Open Data Badge for making publicly available the digitally-shareable data necessary to reproduce the reported results. The data is available at https://doi.org/10.5061/dryad.1tm8705.

Is molecular evolution faster in the tropics?

The evolutionary speed hypothesis (ESH) suggests that molecular evolutionary rates are higher among species inhabiting warmer environments. Previously, the ESH has been investigated using small numbers of latitudinally-separated sister lineages; in animals, these studies typically focused on subsets of Chordata and yielded mixed support for the ESH. This study analyzed public DNA barcode sequences from the cytochrome c oxidase subunit I (COI) gene for six of the largest animal phyla (Arthropoda, Chordata, Mollusca, Annelida, Echinodermata, and Cnidaria) and paired latitudinally-separated taxa together informatically. Of 8037 lineage pairs, just over half (51.6%) displayed a higher molecular rate in the lineage inhabiting latitudes closer to the equator, while the remainder (48.4%) displayed a higher rate in the higher-latitude lineage. To date, this study represents the most comprehensive analysis of latitude-related molecular rate differences across animals. While a statistically-significant pattern was detected from our large sample size, our findings suggest that the EHS may not serve as a strong universal mechanism underlying the latitudinal diversity gradient and that COI molecular clocks may generally be applied across latitudes. This study also highlights the merits of using automation to analyze large DNA barcode datasets.

Barcoding of Chrysomelidae of Euro-Mediterranean area: efficiency and problematic species

Leaf beetles (Coleoptera: Chrysomelidae), with more than 37,000 species worldwide and about 2,300 in the Euro-Mediterranean region, are an ecological and economical relevant family, making their molecular identification of interest also in agriculture. This study, part of the Mediterranean Chrysomelidae Barcoding project (www.c-bar.org), aims to: (i) develop a reference Cytochrome c oxidase I (COI) library for the molecular identification of the Euro-Mediterranean Chrysomelidae; (ii) test the efficiency of DNA barcoding for leaf beetles identification; (iii) develop and compare optimal thresholds for distance-based identifications estimated at family and subfamily level, minimizing false positives and false negatives. Within this study, 889 COI nucleotide sequences of 261 species were provided; after the inclusion of information from other sources, a dataset of 7,237 sequences (542 species) was analysed. The average intra-interspecific distances were in the range of those recorded for Coleoptera: 1.6–24%. The estimated barcoding efficiency (~94%) confirmed the usefulness of this tool for Chrysomelidae identification. The few cases of failure were recorded for closely related species (e.g., Cryptocephalus marginellus superspecies, Cryptocephalus violaceus - Cryptocephalus duplicatus and some Altica species), even with morphologically different species sharing the same COI haplotype. Different optimal thresholds were achieved for the tested taxonomic levels, confirming that group-specific thresholds significantly improve molecular identifications.

Factors affecting interspecific differences in genetic divergence among populations of Anolis lizards in Cuba

BACKGROUND

Geographical patterns and degrees of genetic divergence among populations differ between species, reflecting relative potentials for speciation or cladogenesis and differing capacities for environmental adaptation. Identification of factors that contribute to genetic divergence among populations is important to the understanding of why some species exhibit greater interpopulation genetic divergence. In this study, we calculated the mean pairwise genetic distances among populations as species' average genetic divergence by a phylogeny using nuclear and mitochondrial genes of 303 individuals from 33 Cuban Anolis species and estimated species ages by another phylogeny using nuclear and mitochondrial genes of 51 Cuban and 47 non-Cuban Anolis species. We identified factors that influence species' differences in genetic divergence among 26 species of Anolis lizards from Cuba. Species ages, environmental heterogeneity within species ranges, and ecomorph types were considered as factors affecting average genetic divergences among populations.

RESULTS

The phylogenies presented in this study provide the most comprehensive sampling of Cuban Anolis species to date. The phylogeny showed more conservative evolution of Anolis ecomorphs within Cuba and identified twig anoles as a monophyletic group. Subsequent Phylogenetic Generalized Least Squares (PGLS) analyses showed that species age was positively correlated with species' average genetic divergence among populations.

CONCLUSION

Although previous studies have focused on factors affecting genetic divergence within species, the present study showed for the first time that species differences in genetic divergence could be largely affected by species age.

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.

Slowly dispersing neotenic beetles can speciate on a penny coin and generate space-limited diversity in the tropical mountains

We demonstrate the controversial origin of a biological species within an area of a few kilometres in the absence of physical barriers. We employed nuclear rRNA/mitochondrial and genome-wide SNP approaches to infer relationships of four species of net-winged beetles characterised by female neoteny. Three species are distributed at low elevations and a single population colonised a 40 km(2) highland plateau and established distinct biological species despite incomplete genetic isolation. The speciation process is extreme in the highly localised spatial scale, due to the low dispersal power of neotenics, and provides clear support for a microallopatric model based on ecological conditions. In contrast with neutral evolution in a homogenous environment, as demonstrated by the genetic divergence and morphological similarity of two widely distributed low-mountain species, the environmental characteristics of the high-mountain plateau led to the origin of a species adapted to the local mimetic pattern and characterised by morphologically distinct genitalia. We conclude that the low dispersal propensity promotes neutral genetic differentiation in the first stage, but environmental characteristics play an important role the final phase of the speciation process. The unexpected speciation at such an extreme geographic scale points to the in situ origin and uniqueness of the mountain fauna.

Species mtDNA genetic diversity explained by infrapopulation size in a host-symbiont system

Understanding what shapes variation in genetic diversity among species remains a major challenge in evolutionary ecology, and it has been seldom studied in parasites and other host‐symbiont systems. Here, we studied mtDNA variation in a host‐symbiont non‐model system: 418 individual feather mites from 17 feather mite species living on 17 different passerine bird species. We explored how a surrogate of census size, the median infrapopulation size (i.e., the median number of individual parasites per infected host individual), explains mtDNA genetic diversity. Feather mite species genetic diversity was positively correlated with mean infrapopulation size, explaining 34% of the variation. As expected from the biology of feather mites, we found bottleneck signatures for most of the species studied but, in particular, three species presented extremely low mtDNA diversity values given their infrapopulation size. Their star‐like haplotype networks (in contrast with more reticulated networks for the other species) suggested that their low genetic diversity was the consequence of severe bottlenecks or selective sweeps. Our study shows for the first time that mtDNA diversity can be explained by infrapopulation sizes, and suggests that departures from this relationship could be informative of underlying ecological and evolutionary processes.

Rarity and Incomplete Sampling in DNA-Based Species Delimitation

DNA-based species delimitation may be compromised by limited sampling effort and species rarity, including "singleton" representatives of species, which hampers estimates of intra- versus interspecies evolutionary processes. In a case study of southern African chafers (beetles in the family Scarabaeidae), many species and subclades were poorly represented and 48.5% of species were singletons. Using cox1 sequences from >500 specimens and ∼100 species, the Generalized Mixed Yule Coalescent (GMYC) analysis as well as various other approaches for DNA-based species delimitation (Automatic Barcode Gap Discovery (ABGD), Poisson tree processes (PTP), Species Identifier, Statistical Parsimony), frequently produced poor results if analyzing a narrow target group only, but the performance improved when several subclades were combined. Hence, low sampling may be compensated for by "clade addition" of lineages outside of the focal group. Similar findings were obtained in reanalysis of published data sets of taxonomically poorly known species assemblages of insects from Madagascar. The low performance of undersampled trees is not due to high proportions of singletons per se, as shown in simulations (with 13%, 40% and 52% singletons). However, the GMYC method was highly sensitive to variable effective population size ([Formula: see text]), which was exacerbated by variable species abundances in the simulations. Hence, low sampling success and rarity of species affect the power of the GMYC method only if they reflect great differences in [Formula: see text] among species. Potential negative effects of skewed species abundances and prevalence of singletons are ultimately an issue about the variation in [Formula: see text] and the degree to which this is correlated with the census population size and sampling success. Clade addition beyond a limited study group can overcome poor sampling for the GMYC method in particular under variable [Formula: see text] This effect was less pronounced for methods of species delimitation not based on coalescent models.