Refine
Year of publication
Has Fulltext
- yes (48)
Is part of the Bibliography
- no (48)
Keywords
- Asian bush mosquito (2)
- Positive selection (2)
- genomics (2)
- land snails (2)
- 16S mtDNA (1)
- 16S rRNA (1)
- Adaptation (1)
- Adaptive sequence evolution (1)
- Alburnoides bipunctatus (1)
- Altitudinal (1)
- Art (Spezies) (1)
- B chromosome (1)
- Biological invasion (1)
- Carnivora (1)
- Charles (1)
- Chironomidae (1)
- Chironomus piger (1)
- Chironomus riparius (1)
- Climate change (1)
- Codon models (1)
- DNA-Barcoding (1)
- Darwin (1)
- Diferential delay equation (1)
- Diptera (1)
- Evolution (1)
- Evolutionary ecotoxicology (1)
- Fagaceae (1)
- Fungi (1)
- Generation time (1)
- Gradient (1)
- Helicellinae (1)
- Hi-C (1)
- Hybridisierung (1)
- Invasive species (1)
- Land invasion (1)
- McDonald-Kreitman test (1)
- MinION (1)
- Mitogenomics (1)
- Mollusca (1)
- Multigeneration life-cycle test (1)
- Oxidative stress (1)
- Panpulmonata (1)
- Pool-Seq (1)
- Population continuity (1)
- Population genomics (1)
- Rapoport effect (1)
- Reactive oxygen species (1)
- Rearing temperature (1)
- SARS-CoV-2 (1)
- SNP (1)
- SNPs (1)
- Stage-structured model (1)
- StagePop (1)
- Symbiosis (1)
- Vector mosquito (1)
- Water Framework Directive (1)
- acclimation (1)
- adaptation (1)
- admixture inference (1)
- allozymes (1)
- assembly (1)
- chromosomes (1)
- climate (1)
- climate change genomics (1)
- cuticular hydrocarbons (1)
- desaturase (1)
- developmental temperature (1)
- dispersal (1)
- ectotherm (1)
- elongase (1)
- environmental DNA (1)
- feeding (1)
- forest tree (1)
- formicine (1)
- fungi (1)
- gastropods (1)
- gene-flow (1)
- genome architecture (1)
- genome assembly and annotation (1)
- genomic base composition (1)
- hardwood (1)
- hidden speciation (1)
- hybrid assembly (1)
- insect (1)
- islands of divergence (1)
- landscape structure (1)
- latent factor mixed model (1)
- life cycle (1)
- macroecology (1)
- microsatellites (1)
- mitochondrial genome (1)
- monitoring (1)
- mutational spectrum bias (1)
- next generation sequencing (1)
- noninvasive species (1)
- population structure (1)
- quantitative shell variation (1)
- raccoon dog (Nyctereutes procyonoides) (1)
- random genetic drift (1)
- range expansion (1)
- reintroduction (1)
- repeat elements (1)
- repetitive DNA (1)
- reproductive isolation (1)
- reserve substance (1)
- short read assembly (1)
- speciation (1)
- species distribution models (1)
- temperature adaptation (1)
- transcriptomics (1)
- water parameter (1)
- whole genome pooled sequencing (1)
- yellow fever mosquito (1)
Institute
Background: As ectothermic animals, temperature influences insects in almost every aspect. The potential disease spreading Asian bush mosquito (Aedes japonicus japonicus) is native to temperate East Asia but invasive in several parts of the world. We report on the previously poorly understood temperature-dependence of its life history under laboratory conditions to understand invasion processes and to model temperature niches.
Results: To evaluate winter survival, eggs were exposed between 1 day and 14 days to low temperatures (5 °C, 0 °C, -5 °C and -9 °C). Hatching success was drastically decreased after exposure to 0 °C and -5 °C, and the minimal hatching success of 0% was reached at -9 °C after two days. We then exposed larvae to 14 temperatures and assessed their life trait parameters. Larval survival to adulthood was only possible between 10 °C and 31 °C. Based on this, we modelled the optimal (25 °C), minimal (7 °C) and maximal (31 °C) temperature for cumulative female survival. The time to adult emergence ranges from 12 days to 58 days depending on temperature. We used an age-at-emergence-temperature model to calculate the number of potential generations per year for the Asian bush mosquito in Germany with an average of 4.72 potential generations. At lower temperatures, individuals grew larger than at higher temperatures with female R1 length ranging from 3.04 ± 0.1 mm at 31 °C to 4.26 ± 0.2 mm at 15 °C.
Conclusions: Reduced egg hatch after exposure to sub-zero temperatures prohibits the establishment of the Asian bush mosquito in large parts of Germany. Larval overwintering is not possible at temperature ≤ 5 °C. The many potential generations displayed per year may contribute to the species’ invasion success. This study on the thermal ecology of the Asian bush mosquito adds to our knowledge on the temperature dependence of the species and data could be incorporated in epidemiological and population dynamic modelling.
Mutations are the ultimate basis of evolution, yet their occurrence rate is known only for few species. We directly estimated the spontaneous mutation rate and the mutational spectrum in the nonbiting midge C. riparius with a new approach. Individuals from ten mutation accumulation lines over five generations were deep genome sequenced to count de novo mutations that were not present in a pool of F1 individuals, representing parental genotypes. We identified 51 new single site mutations of which 25 were insertions or deletions and 26 single nucleotide mutations. This shift in the mutational spectrum compared to other organisms was explained by the high A/T content of the species. We estimated a haploid mutation rate of 2.1 × 10−9 (95% confidence interval: 1.4 × 10−9 – 3.1 × 10-9) that is in the range of recent estimates for other insects and supports the drift barrier hypothesis. We show that accurate mutation rate estimation from a high number of observed mutations is feasible with moderate effort even for nonmodel species.
Background: While research on the impact of global climate change (GCC) on ecosystems and species is flourishing, a fundamental component of biodiversity -- molecular variation -- has not yet received its due attention in such studies. Here we present a methodological framework for projecting the loss of intraspecific genetic diversity due to GCC.
Methods: The framework consists of multiple steps that and combines 1) hierarchical genetic clustering methods to define comparable units of inference, 2) species accumulation curves (SAC) to infer sampling completeness, and 3) species distribution modelling (SDM) to project the genetic diversity loss under GCC. We suggest procedures for existing data sets as well as specifically designed studies. We illustrate the approach with two worked examples from a land snail (Trochulus villosus) and a caddisfly (Smicridea (S.) mucronata).
Results: Sampling completeness was diagnosed on the third most coarse haplotype clade level for T. villosus and the second most coarse for S. mucronata. For both species, a substantial species range loss was projected under the chosen climate scenario. However, despite substantial differences in data set quality concerning spatial sampling and sampling depth, no loss of haplotype clades due to GCC was predicted for either species.
Conclusions: The suggested approach presents a feasible method to tap the rich resources of existing phylogeographic data sets and guide the design and analysis of studies explicitly designed to estimate the impact of GCC on a currently still neglected level of biodiversity.
Abstract
Divergence is mostly viewed as a progressive process often initiated by selection targeting individual loci, ultimately resulting in ever increasing genomic isolation due to linkage. However, recent studies show that this process may stall at intermediate stable equilibrium states without achieving complete genomic isolation. We tested the extent of genomic isolation between two recurrently hybridizing nonbiting midge sister taxa, Chironomus riparius and Chironomus piger, by analyzing the divergence landscape. Using a principal component‐based method, we estimated that only about 28.44% of the genomes were mutually isolated, whereas the rest was still exchanged. The divergence landscape was fragmented into isolated regions of on average 30 kb, distributed throughout the genome. Selection and divergence time strongly influenced lengths of isolated regions, whereas local recombination rate only had minor impact. Comparison of divergence time distributions obtained from several coalescence‐simulated divergence scenarios with the observed divergence time estimates in an approximate Bayesian computation framework favored a short and concluded divergence event in the past. Most divergence happened during a short time span about 4.5 million generations ago, followed by a stable equilibrium between mutual gene flow through ongoing hybridization for the larger part of the genome and isolation in some regions due to rapid purifying selection of introgression, supported by high effective population sizes and recombination rates.
Impact Summary
The process of speciation has fascinated biologists from early on. Prevailing theory suggested that gene flow among populations is the main obstacle for their divergence. Recently, it became clear that speciation with gene flow is possible under certain circumstances. However, it remains unclear how the divergence process proceeds in time, how widespread the phenomenon is, and whether it always and inevitably leads to complete isolation. Comparing the genomes of individuals of two regularly hybridizing sister taxa of nonbiting midges, we could show that they diverged during a short period millions of generations ago. Their divergence process apparently ceased before the entire genome was mutually isolated. The taxa remain distinct since, even though they share most of their genome. Our findings thus extend our view of the nature of species and the temporal dynamics of their divergence and describe novel approaches to analyze both current and past divergence processes.
In the course of global climate change, central Europe is experiencing more frequent and prolonged periods of drought. The drought years 2018 and 2019 affected European beeches (Fagus sylvatica L.) differently: even in the same stand, drought damaged trees neighboured healthy trees, suggesting that the genotype rather than the environment was responsible for this conspicuous pattern. We used this natural experiment to study the genomic basis of drought resistance with Pool-GWAS. Contrasting the extreme phenotypes identified 106 significantly associated SNPs throughout the genome. Most annotated genes with associated SNPs (>70%) were previously implicated in the drought reaction of plants. Non-synonymous substitutions led either to a functional amino acid exchange or premature termination. A SNP-assay with 70 loci allowed predicting drought phenotype in 98.6% of a validation sample of 92 trees. Drought resistance in European beech is a moderately polygenic trait that should respond well to natural selection, selective management, and breeding.
The gradual heterogeneity of climatic factors pose varying selection pressures across geographic distances that leave signatures of clinal variation in the genome. Separating signatures of clinal adaptation from signatures of other evolutionary forces, such as demographic processes, genetic drift, and adaptation to non-clinal conditions of the immediate local environment is a major challenge. Here, we examine climate adaptation in five natural populations of the harlequin fly Chironomus riparius sampled along a climatic gradient across Europe. Our study integrates experimental data, individual genome resequencing, Pool-Seq data, and population genetic modelling. Common-garden experiments revealed a positive correlation of population growth rates corresponding to the population origin along the climate gradient, suggesting thermal adaptation on the phenotypic level. Based on a population genomic analysis, we derived empirical estimates of historical demography and migration. We used an FST outlier approach to infer positive selection across the climate gradient, in combination with an environmental association analysis. In total we identified 162 candidate genes as genomic basis of climate adaptation. Enriched functions among these candidate genes involved the apoptotic process and molecular response to heat, as well as functions identified in other studies of climate adaptation in other insects. Our results show that local climate conditions impose strong selection pressures and lead to genomic adaptation despite strong gene flow. Moreover, these results imply that selection to different climatic conditions seems to converge on a functional level, at least between different insect species.
Genomic basis of ecological niche divergence among cryptic sister species of non-biting midges
(2013)
Background: There is a lack of understanding the evolutionary forces driving niche segregation of closely related organisms. In addition, pinpointing the genes driving ecological divergence is a key goal in molecular ecology. Here, larval transcriptome sequences obtained by next-generation-sequencing are used to address these issues in a morphologically cryptic sister species pair of non-biting midges (Chironomus riparius and C. piger).
Results: More than eight thousand orthologous open reading frames were screened for interspecific divergence and intraspecific polymorphisms. Despite a small mean sequence divergence of 1.53% between the sister species, 25.1% of 18,115 observed amino acid substitutions were inferred by α statistics to be driven by positive selection. Applying McDonald-Kreitman tests to 715 alignments of gene orthologues identified eleven (1.5%) genes driven by positive selection.
Conclusions: Three candidate genes were identified as potentially responsible for the observed niche segregation concerning nitrite concentration, habitat temperature and water conductivity. Additionally, signs of positive selection in the hydrogen sulfide detoxification pathway were detected, providing a new plausible hypothesis for the species’ ecological differentiation. Finally, a divergently selected, nuclear encoded mitochondrial ribosomal protein may contribute to reproductive isolation due to cytonuclear coevolution.
Background: Factors and processes shaping the population structure and spatial distribution of genetic diversity across a species' distribution range are important in determining the range limits. We comprehensively analysed the influence of recurrent and historic factors and processes on the population genetic structure, mating system and the distribution of genetic variability of the pulmonate freshwater snail Radix balthica. This analysis was based on microsatellite variation and mitochondrial haplotypes using Generalised Linear Statistical Modelling in a Model Selection framework. Results: Populations of R. balthica were found throughout North-Western Europe with range margins marked either by dispersal barriers or the presence of other Radix taxa. Overall, the population structure was characterised by distance independent passive dispersal mainly along a Southwest-Northeast axis, the absence of isolation-by-distance together with rather isolated and genetically depauperated populations compared to the variation present in the entire species due to strong local drift. A recent, climate driven range expansion explained most of the variance in genetic variation, reducing at least temporarily the genetic variability in this area. Other factors such as geographic marginality and dispersal barriers play only a minor role. Conclusions: To our knowledge, such a population structure has rarely been reported before. It might nevertheless be typical for passively dispersed, patchily distributed taxa (e.g. freshwater invertebrates). The strong local drift implied in such a structure is expected to erode genetic variation at both neutral and coding loci and thus probably diminish evolutionary potential. This study shows that the analysis of multiple factors is crucial for the inference of the processes shaping the distribution of genetic variation throughout species ranges. Additional files Additional file 1: Distribution of Radix taxa. Spatial distribution of the Radix MOTU as defined in Pfenninger et al. 2006 plus an additional, newly discovered taxon. This map is the basis for the inference of the species range of R. balthica. Additional file 2: Sampling site table and spatial distribution of diversity indices, selfing estimates and inferred population bottlenecks for R. balthica. Table of sampling site code, geographical position in decimal degrees latitude and longitude, number of individuals analysed with microsatellites (Nnuc), expected heterozygosity (HE) and standard deviation across loci, mean rarefied number of alleles per microsatellite locus (A) and their standard deviation, number of individuals analysed for mitochondrial variation (Nmt), rarefied number of mitochondrial COI haplotypes (Hmt), number of individuals measured for body size (Nsize). Figures A1 - A3 show a graphical representation of the spatial distribution of He, Hmt and, s, respectively. Additional file 3: Assessment of environmental marginality. PCA (principle component analysis) on 35 climatic parameters for the period from 1960 - 2000 from publicly availableWorldClim data. Additional file 4: Inference of a recent climate driven range expansion in R. balthica. Analysis of the freshwater benthos long term monitoring data of the Swedish national monitoring databases at the Swedish University of Agricultural Sciences SLU with canonical correspondence analysis.
Background: Until recently, read lengths on the Solexa/Illumina system were too short to reliably assemble transcriptomes without a reference sequence, especially for non-model organisms. However, with read lengths up to 100 nucleotides available in the current version, an assembly without reference genome should be possible. For this study we created an EST data set for the common pond snail Radix balthica by Illumina sequencing of a normalized transcriptome. Performance of three different short read assemblers was compared with respect to: the number of contigs, their length, depth of coverage, their quality in various BLAST searches and the alignment to mitochondrial genes. Results: A single sequencing run of a normalized RNA pool resulted in 16,923,850 paired end reads with median read length of 61 bases. The assemblies generated by VELVET, OASES, and SeqMan NGEN differed in the total number of contigs, contig length, the number and quality of gene hits obtained by BLAST searches against various databases, and contig performance in the mt genome comparison. While VELVET produced the highest overall number of contigs, a large fraction of these were of small size (< 200bp), and gave redundant hits in BLAST searches and the mt genome alignment. The best overall contig performance resulted from the NGEN assembly. It produced the second largest number of contigs, which on average were comparable to the OASES contigs but gave the highest number of gene hits in two out of four BLAST searches against different reference databases. A subsequent meta-assembly of the four contig sets resulted in larger contigs, less redundancy and a higher number of BLAST hits. Conclusion: Our results document the first de novo transcriptome assembly of a non-model species using Illumina sequencing data. We show that de novo transcriptome assembly using this approach yields results useful for downstream applications, in particular if a meta-assembly of contig sets is used to increase contig quality. These results highlight the ongoing need for improvements in assembly methodology. Keywords: next generation sequencing; short read assembly; Mollusca
Effects of seasonal or daily temperature variation on fitness and physiology of ectothermic organisms and their ways to cope with such variations have been widely studied. However, the way multivoltines organisms cope with temperature variations from one generation to the next is still not well understood. The aim of this study was to investigate whether the multivoltine midge Chironomus riparius Meigen (1803) responds mainly via acclimation as predicted by current theories or whether rapid genetic adaptation is involved. To investigate this issue, a common garden approach has been applied. A mix of larvae from five European populations was raised in the laboratory at three different pre‐exposure temperatures (PET): 14, 20, and 26°C. After three and five generations, respectively, larvae were exposed to three treatment temperatures (TT): 14, 20, and 26°C. Mortality was monitored for the first 48 hr and after emergence. After three generations, significant mortality rate differences depended on an interaction of PET and TT. This finding supports the hypothesis that chironomids respond rapidly to climatic variation via adaptive mechanisms and to a lesser extent via phenotypic plasticity. The result of the experiment indicates that three generations were sufficient to adapt to warm temperature, decreasing the mortality rate, highlighting the potential for chironomids to rapidly respond to seasonally changing conditions.
The success of social insects is largely intertwined with their highly advanced chemical communication system that facilitates recognition and discrimination of species and nest-mates, recruitment, and division of labor. Hydrocarbons, which cover the cuticle of insects, not only serve as waterproofing agents but also constitute a major component of this communication system. Two cryptic Crematogaster species, which share their nest with Camponotus ants, show striking diversity in their cuticular hydrocarbon (CHC) profile. This mutualistic system therefore offers a great opportunity to study the genetic basis of CHC divergence between sister species. As a basis for further genome-wide studies high-quality genomes are needed. Here, we present the annotated draft genome for Crematogaster levior A. By combining the three most commonly used sequencing techniques—Illumina, PacBio, and Oxford Nanopore—we constructed a high-quality de novo ant genome. We show that even low coverage of long reads can add significantly to overall genome contiguity. Annotation of desaturase and elongase genes, which play a role in CHC biosynthesis revealed one of the largest repertoires in ants and a higher number of desaturases in general than in other Hymenoptera. This may provide a mechanistic explanation for the high diversity observed in C. levior CHC profiles.
The European Beech is the dominant climax tree in most regions of Central Europe and valued for its ecological versatility and hardwood timber. Even though a draft genome has been published recently, higher resolution is required for studying aspects of genome architecture and recombination. Here we present a chromosome-level assembly of the more than 300 year-old reference individual, Bhaga, from the Kellerwald-Edersee National Park (Germany). Its nuclear genome of 541 Mb was resolved into 12 chromosomes varying in length between 28 Mb and 73 Mb. Multiple nuclear insertions of parts of the chloroplast genome were observed, with one region on chromosome 11 spanning more than 2 Mb of the genome in which fragments up to 54,784 bp long and covering the whole chloroplast genome were inserted randomly. Unlike in Arabidopsis thaliana, ribosomal cistrons are present in Fagus sylvatica only in four major regions, in line with FISH studies. On most assembled chromosomes, telomeric repeats were found at both ends, while centromeric repeats were found to be scattered throughout the genome apart from their main occurrence per chromosome. The genome- wide distribution of SNPs was evaluated using a second individual from Jamy Nature Reserve (Poland). SNPs, repeat elements and duplicated genes were unevenly distributed in the genomes, with one major anomaly on chromosome 4. The genome presented here adds to the available highly resolved plant genomes and we hope it will serve as a valuable basis for future research on genome architecture and for understanding the past and future of European Beech populations in a changing climate.
Background: The European beech is arguably the most important climax broad-leaved tree species in Central Europe, widely planted for its valuable wood. Here, we report the 542 Mb draft genome sequence of an up to 300-year-old individual (Bhaga) from an undisturbed stand in the Kellerwald-Edersee National Park in central Germany.
Findings: Using a hybrid assembly approach, Illumina reads with short- and long-insert libraries, coupled with long Pacific Biosciences reads, we obtained an assembled genome size of 542 Mb, in line with flow cytometric genome size estimation. The largest scaffold was of 1.15 Mb, the N50 length was 145 kb, and the L50 count was 983. The assembly contained 0.12% of Ns. A Benchmarking with Universal Single-Copy Orthologs (BUSCO) analysis retrieved 94% complete BUSCO genes, well in the range of other high-quality draft genomes of trees. A total of 62,012 protein-coding genes were predicted, assisted by transcriptome sequencing. In addition, we are reporting an efficient method for extracting high-molecular-weight DNA from dormant buds, by which contamination by environmental bacteria and fungi was kept at a minimum.
Conclusions: The assembled genome will be a valuable resource and reference for future population genomics studies on the evolution and past climate change adaptation of beech and will be helpful for identifying genes, e.g., involved in drought tolerance, in order to select and breed individuals to adapt forestry to climate change in Europe. A continuously updated genome browser and download page can be accessed from beechgenome.net, which will include future genome versions of the reference individual Bhaga, as new sequencing approaches develop.
Molluscs are the second most species-rich phylum in the animal kingdom, yet only 11 genomes of this group have been published so far. Here, we present the draft genome sequence of the pulmonate freshwater snail Radix auricularia. Six whole genome shotgun libraries with different layouts were sequenced. The resulting assembly comprises 4,823 scaffolds with a cumulative length of 910 Mb and an overall read coverage of 72×. The assembly contains 94.6% of a metazoan core gene collection, indicating an almost complete coverage of the coding fraction. The discrepancy of ∼690 Mb compared with the estimated genome size of R. auricularia (1.6 Gb) results from a high repeat content of 70% mainly comprising DNA transposons. The annotation of 17,338 protein coding genes was supported by the use of publicly available transcriptome data. This draft will serve as starting point for further genomic and population genetic research in this scientifically important phylum.
Precise estimates of genome sizes are important parameters for both theoretical and practical biodiversity genomics. We present here a fast, easy-to-implement and precise method to estimate genome size from the number of bases sequenced and the mean sequence coverage. To estimate the latter, we take advantage of the fact that a precise estimation of the Poisson distribution parameter lambda is possible from truncated data, restricted to the part of the coverage distribution representing the true underlying distribution. With simulations we could show that reasonable genome size estimates can be gained even from low-coverage (10X), highly discontinuous genome drafts. Comparison of estimates from a wide range of taxa and sequencing strategies with flow-cytometry estimates of the same individuals showed a very good fit and suggested that both methods yield comparable, interchangeable results.
Background: Driven by globalization, urbanization and climate change, the distribution range of invasive vector species has expanded to previously colder ecoregions. To reduce health-threatening impacts on humans, insect vectors are extensively studied. Population genomics can reveal the genomic basis of adaptation and help to identify emerging trends of vector expansion.
Results: By applying whole genome analyses and genotype-environment associations to populations of the main dengue vector Ae. aegypti, sampled along an altitudinal temperature gradient in Nepal (200- 1300m), we identify adaptive traits and describe the species’ genomic footprint of climate adaptation to colder ecoregions. We found two clusters of differentiation with significantly different allele frequencies in genes associated to climate adaptation between the highland population (1300m) and all other lowland populations (≤ 800 m). We revealed non-synonymous mutations in 13 of the candidate genes associated to either altitude, precipitation or cold tolerance and identified an isolation-by-environment differentiation pattern.
Conclusion: Other than the expected gradual differentiation along the altitudinal gradient, our results reveal a distinct genomic differentiation of the highland population. This finding either indicates a differential invasion history to Nepal or local high-altitude adaptation explaining the population’s phenotypic cold tolerance. In any case, this highland population can be assumed to carry pre-adapted alleles relevant for the species’ invasion into colder ecoregions worldwide that way expanding their climate niche.
The European Beech is the dominant climax tree in most regions of Central Europe and valued for its ecological versatility and hardwood timber. Even though a draft genome has been published recently, higher resolution is required for studying aspects of genome architecture and recombination. Here, we present a chromosome-level assembly of the more than 300 year-old reference individual, Bhaga, from the Kellerwald-Edersee National Park (Germany). Its nuclear genome of 541 Mb was resolved into 12 chromosomes varying in length between 28 and 73 Mb. Multiple nuclear insertions of parts of the chloroplast genome were observed, with one region on chromosome 11 spanning more than 2 Mb which fragments up to 54,784 bp long and covering the whole chloroplast genome were inserted randomly. Unlike in Arabidopsis thaliana, ribosomal cistrons are present in Fagus sylvatica only in four major regions, in line with FISH studies. On most assembled chromosomes, telomeric repeats were found at both ends, while centromeric repeats were found to be scattered throughout the genome apart from their main occurrence per chromosome. The genome-wide distribution of SNPs was evaluated using a second individual from Jamy Nature Reserve (Poland). SNPs, repeat elements and duplicated genes were unevenly distributed in the genomes, with one major anomaly on chromosome 4. The genome presented here adds to the available highly resolved plant genomes and we hope it will serve as a valuable basis for future research on genome architecture and for understanding the past and future of European Beech populations in a changing climate.
The Alpine Region, constituting the Alps and the Dinaric Alps, has played a major role in the formation of current patterns of biodiversity either as a contact zone of postglacial expanding lineages or as the origin of genetic diversity. In our study, we tested these hypotheses for two widespread, sympatric microgastropod taxa - Carychium minimum O.F. Müller, 1774 and Carychium tridentatum (Risso, 1826) (Gastropoda, Eupulmonata, Carychiidae) - by using COI sequence data and species potential distribution models analyzed in a statistical phylogeographical framework. Additionally, we examined disjunct transatlantic populations of those taxa from the Azores and North America. In general, both Carychium taxa demonstrate a genetic structure composed of several differentiated haplotype lineages most likely resulting from allopatric diversification in isolated refugial areas during the Pleistocene glacial periods. However, the genetic structure of Carychium minimum is more pronounced, which can be attributed to ecological constraints relating to habitat proximity to permanent bodies of water. For most of the Carychium lineages, the broader Alpine Region was identified as the likely origin of genetic diversity. Several lineages are endemic to the broader Alpine Region whereas a single lineage per species underwent a postglacial expansion to (re)colonize previously unsuitable habitats, e.g. in Northern Europe. The source populations of those expanding lineages can be traced back to the Eastern and Western Alps. Consequently, we identify the Alpine Region as a significant 'hot-spot' for the formation of genetic diversity within European Carychium lineages. Passive dispersal via anthropogenic means best explains the presence of transatlantic European Carychium populations on the Azores and in North America. We conclude that passive (anthropogenic) transport could mislead the interpretation of observed phylogeographical patterns in general.
One of the major problems in evolutionary biology is to elucidate the relationships between historical events and the tempo and mode of lineage divergence. The development of relaxed molecular clock models and the increasing availability of DNA sequences resulted in more accurate estimations of taxa divergence times. However, finding the link between competing historical events and divergence is still challenging. Here we investigate assigning constrained-age priors to nodes of interest in a time-calibrated phylogeny as a means of hypothesis comparison. These priors are equivalent to historic scenarios for lineage origin. The hypothesis that best explains the data can be selected by comparing the likelihood values of the competing hypotheses, modelled with different priors. A simulation approach was taken to evaluate the performance of the prior-based method and to compare it with an unconstrained approach. We explored the effect of DNA sequence length and the temporal placement and span of competing hypotheses (i.e. historic scenarios) on selection of the correct hypothesis and the strength of the inference. Competing hypotheses were compared applying a posterior simulation analogue of the Akaike Information Criterion and Bayes factors (obtained after calculation of the marginal likelihood with three estimators: Harmonic Mean, Stepping Stone and Path Sampling). We illustrate the potential application of the prior-based method on an empirical data set to compare competing geological hypotheses explaining the biogeographic patterns in Pleurodeles newts. The correct hypothesis was selected on average 89% times. The best performance was observed with DNA sequence length of 3500-10000 bp. The prior-based method is most reliable when the hypotheses compared are not temporally too close. The strongest inferences were obtained when using the Stepping Stone and Path Sampling estimators. The prior-based approach proved effective in discriminating between competing hypotheses when used on empirical data. The unconstrained analyses performed well but it probably requires additional computational effort. Researchers applying this approach should rely only on inferences with moderate to strong support. The prior-based approach could be applied on biogeographical and phylogeographical studies where robust methods for historical inferences are still lacking.