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)
Institute
Background: The faunal and floral relationship of northward-drifting India with its neighboring continents is of general biogeographic interest as an important driver of regional biodiversity. However, direct biogeographic connectivity of India and Southeast Asia during the Cenozoic remains largely unexplored. We investigate timing, direction and mechanisms of faunal exchange between India and Southeast Asia, based on a molecular phylogeny, molecular clock-derived time estimates and biogeographic reconstructions of the Asian freshwater crab family Gecarcinucidae. Results: Although the Gecarcinucidae are not an element of an ancient Gondwana fauna, their subfamily Gecarcinucinae, and probably also the Liotelphusinae, evolved on the Indian Subcontinent and subsequently dispersed to Southeast Asia. Estimated by a model testing approach, this dispersal event took place during the Middle Eocene, and thus before the final collision of India and the Tibet-part of Eurasia. Conclusions: We postulate that the India and Southeast Asia were close enough for exchange of freshwater organisms during the Middle Eocene, before the final Indian--Eurasian collision. Our data support geological models that assume the Indian plate having tracked along Southeast Asia during its move northwards.
The Asian bush mosquito (Aedes japonicus japonicus, Theobald 1901) is an invasive culicid species which originates in Asia but is nowadays present in northern America and Europe. It is a competent vector for several human disease pathogens. In addition to the public health threat, this invasive species may also be an ecological threat for native container-breeding mosquitoes which share a similar larval habitat. Therefore, it is of importance to gain knowledge on ecological and eco-toxicological features of the Asian bush mosquito. However, optimal laboratory feeding conditions have not yet been established. Standardized feeding methods will be needed in assessing the impact of insecticides or competitional strength of this species. To fill this gap, we performed experiments on food quality and quantity for Ae. j. japonicus larvae. We found out that the commercial fish food TetraMin (Tetra, Melle, Germany) in a dose of 10 mg per larva is the most suitable food tested. We also suggest a protocol with a feeding sequence of seven portions for all larval stages of this species.
The data provided is related to the article "Phylogenetic analyses of gazelles reveal repeated transitions of key ecological traits and provide novel insights into the origin of the genus Gazella". The data is based on 48 tissue samples of all nine extant species of the genus Gazella, namely Gazella gazella, Gazella arabica, Gazella bennettii, Gazella cuvieri, Gazella dorcas, Gazella leptoceros, Gazella marica, Gazella spekei, and Gazella subgutturosa and four related taxa (Saiga tatarica, Antidorcas marsupialis, Antilope cervicapra and Eudorcas rufifrons). It comprises alignments of sequences of a cytochrome b data set and of six nuclear intron markers. For the latter new primers were designed based on cattle and sheep genomes. Based on these alignments phylogenetic trees were inferred using Bayesian Inference and Maximum Likelihood methods. Furthermore, ancestral character states (inferred with BayesTraits 1.0) and ancestral ranges based on a Dispersal-Extinction-Cladogenesis model were estimated and results׳ files were stored within this article.
Active species reintroduction is an important conservation tool when aiming for the restoration of biological communities and ecosystems. The effective monitoring of reintroduction success is a crucial factor in this process. Here, we used a combination of environmental DNA (eDNA) techniques and species distribution models (SDMs) to evaluate the success of recent reintroductions of the freshwater fish Alburnoides bipunctatus in central Germany. We built SDMs without and with eDNA presence data to locate further suitable reintroduction sites and potentially overlooked populations of the species. We successfully detected eDNA of A. bipunctatus at all reintroduction sites, as well as several adjacent sites mostly in downstream direction, which supports the success of reintroduction efforts. eDNA‐based species detection considerably improved SDMs for A. bipunctatus, which allowed to identify species presence in previously unknown localities. Our results confirm the usefulness of eDNA techniques as standard tool to monitor reintroduced fish populations. We propose that combining eDNA with SDMs is a highly effective approach for long‐term monitoring of reintroduction success in aquatic species.
Background One of the central issues in ecology is the question what allows sympatric occurrence of closely related species in the same general area? The non-biting midges Chironomus riparius and C. piger, interbreeding in the laboratory, have been shown to coexist frequently despite of their close relatedness, similar ecology and high morphological similarity. Methodology/Principal Findings In order to investigate factors shaping niche partitioning of these cryptic sister species, we explored the actual degree of reproductive isolation in the field. Congruent results from nuclear microsatellite and mitochondrial haplotype analyses indicated complete absence of interspecific gene-flow. Autocorrelation analysis showed a non-random spatial distribution of the two species. Though not dispersal limited at the scale of the study area, the sister species occurred less often than expected at the same site, indicating past or present competition. Correlation and multiple regression analyses suggested the repartition of the available habitat along water chemistry gradients (nitrite, conductivity, CaCO3), ultimately governed by differences in summer precipitation regime. Conclusions We show that these morphologically cryptic sister species partition their niches due to a certain degree of ecological distinctness and total reproductive isolation in the field. The coexistence of these species provides a suitable model system for the investigation of factors shaping the distribution of closely related, cryptic species.
Bears are iconic mammals with a complex evolutionary history. Natural bear hybrids and studies of few nuclear genes indicate that gene flow among bears may be more common than expected and not limited to polar and brown bears. Here we present a genome analysis of the bear family with representatives of all living species. Phylogenomic analyses of 869 mega base pairs divided into 18,621 genome fragments yielded a well-resolved coalescent species tree despite signals for extensive gene flow across species. However, genome analyses using different statistical methods show that gene flow is not limited to closely related species pairs. Strong ancestral gene flow between the Asiatic black bear and the ancestor to polar, brown and American black bear explains uncertainties in reconstructing the bear phylogeny. Gene flow across the bear clade may be mediated by intermediate species such as the geographically wide-spread brown bears leading to large amounts of phylogenetic conflict. Genome-scale analyses lead to a more complete understanding of complex evolutionary processes. Evidence for extensive inter-specific gene flow, found also in other animal species, necessitates shifting the attention from speciation processes achieving genome-wide reproductive isolation to the selective processes that maintain species divergence in the face of gene flow.
Background: The invasive temperate mosquito Aedes japonicus japonicus is a potential vector for various infectious diseases and therefore a target of vector control measures. Even though established in Germany, it is unclear whether the species has already reached its full distribution potential. The possible range of the species, its annual population dynamics, the success of vector control measures and future expansions due to climate change still remain poorly understood. While numerous studies on occurrence have been conducted, they used mainly presence data from relatively few locations. In contrast, we used experimental life history data to model the dynamics of a continuous stage-structured population to infer potential seasonal densities and ask whether stable populations are likely to establish over a period of more than one year. In addition, we used climate change models to infer future ranges. Finally, we evaluated the effectiveness of various stage-specific vector control measures.
Results: Aedes j. japonicus has already established stable populations in the southwest and west of Germany. Our models predict a spread of Ae. j. japonicus beyond the currently observed range, but likely not much further eastwards under current climatic conditions. Climate change models, however, will expand this range substantially and higher annual densities can be expected. Applying vector control measures to oviposition, survival of eggs, larvae or adults showed that application of adulticides for 30 days between late spring and early autumn, while ambient temperatures are above 9 °C, can reduce population density by 75%. Continuous application of larvicide showed similar results in population reduction. Most importantly, we showed that with the consequent application of a mixed strategy, it should be possible to significantly reduce or even extinguish existing populations with reasonable effort.
Conclusion: Our study provides valuable insights into the mechanisms concerning the establishment of stable populations in invasive species. In order to minimise the hazard to public health, we recommend vector control measures to be applied in ‘high risk areas’ which are predicted to allow establishment of stable populations to establish.
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.
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.
There is increasing evidence that rapid phenotypic adaptation of quantitative traits is not uncommon in nature. However, the circumstances under which rapid adaptation of polygenic traits occurs are not yet understood. Building on previous concepts of soft selection, i.e. frequency and density dependent selection, I developed and tested the hypothesis that adaptation speed of a polygenic trait depends on the number of offspring per breeding pair in a randomly mating diploid population.
Using individual based modelling on a range of offspring per parent (2–200) in populations of various size (100–10000 individuals), I could show that the by far largest proportion of variance (42%) was explained by the offspring number, regardless of genetic trait architecture (10–50 loci, different locus contribution distributions). In addition, it was possible to identify the majority of the responsible loci and account for even more of the observed phenotypic change with a moderate population size.
The simulation results suggest that offspring numbers may a crucial factor for the adaptation speed of quantitative loci. Moreover, as large offspring numbers translates to a large phenotypic variance in the offspring of each parental pair, this genetic bet hedging strategy increases the chance to contribute to the next generation in unpredictable environments.
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 non-biting 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 (DNMs) 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 point 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 x 10−9 (95% confidence interval: 1.4 x 10−9 – 3.1 x 10−9) which 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 non-model species.
Positive selection on panpulmonate mitogenomes provide new clues on adaptations to terrestrial life
(2016)
Background: Transitions from marine to intertidal and terrestrial habitats resulted in a significant adaptive radiation within the Panpulmonata (Gastropoda: Heterobranchia). This clade comprises several groups that invaded the land realm independently and in different time periods, e.g., Ellobioidea, Systellomatophora, and Stylommatophora. Thus, mitochondrial genomes of panpulmonate gastropods are promising to screen for adaptive molecular signatures related to land invasions.
Results: We obtained three complete mitochondrial genomes of terrestrial panpulmonates, i.e., the ellobiid Carychium tridentatum, and the stylommatophorans Arion rufus and Helicella itala. Our dataset consisted of 50 mitogenomes comprising almost all major panpulmonate lineages. The phylogenetic tree based on mitochondrial genes supports the monophyly of the clade Panpulmonata. Terrestrial lineages were sampled from Ellobioidea (1 sp.) and Stylommatophora (9 spp.). The branch-site test of positive selection detected significant non-synonymous changes in the terrestrial branches leading to Carychium (Ellobiodea) and Stylommatophora. These convergent changes occurred in the cob and nad5 genes (OXPHOS complex III and I, respectively).
Conclusions: The convergence of the non-synonymous changes in cob and nad5 suggest possible ancient episodes of positive selection related to adaptations to non-marine habitats. The positively selected sites in our data are in agreement with previous results in vertebrates suggesting a general pattern of adaptation to the new metabolic requirements. The demand for energy due to the colonization of land (for example, to move and sustain the body mass in the new habitat) and the necessity to tolerate new conditions of abiotic stress may have changed the physiological constraints in the early terrestrial panpulmonates and triggered adaptations at the mitochondrial level.
Background Cryptic species are two or more distinct but morphologically similar species that were classified as a single species. During the past two decades we observed an exponential growth of publications on cryptic species. Recently published reviews have demonstrated cryptic species have profound consequences on many biological disciplines. It has been proposed that their distribution is non-random across taxa and biomes. Results We analysed a literature database for the taxonomic and biogeographical distribution of cryptic animal species reports. Results from regression analysis indicate that cryptic species are almost evenly distributed among major metazoan taxa and biogeographical regions when corrected for species richness and study intensity. Conclusion This indicates that morphological stasis represents an evolutionary constant and that cryptic metazoan diversity does predictably affect estimates of earth´s animal diversity. Our findings have direct theoretical and practical consequences for a number of prevailing biological questions with regard to global biodiversity estimates, conservation efforts and global taxonomic initiatives.
Molluscs are the second most species-rich phylum in the animal kingdom, yet only eleven 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 72x. 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 to 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.
Active transposable elements (TEs) may result in divergent genomic insertion and abundance patterns among conspecific populations. Upon secondary contact, such divergent genetic backgrounds can theoretically give rise to classical Dobzhansky-Muller incompatibilities (DMI), a way how TEs can contribute to the evolution of endogenous genetic barriers and eventually population divergence. We investigated whether differential TE activity created endogenous selection pressures among conspecific populations of the non-biting midge Chironomus riparius, focussing on a Chironomus-specific TE, the minisatellite-like Cla-element, whose activity is associated with speciation in the genus. Using an improved and annotated draft genome for a genomic study with five natural C. riparius populations, we found highly population-specific TE insertion patterns with many private insertions. A highly significant correlation of pairwise population FST from genome-wide SNPs with the FST estimated from TEs suggests drift as the major force driving TE population differentiation. However, the significantly higher Cla-element FST level due to a high proportion of differentially fixed Cla-element insertions indicates that segregating, i.e. heterozygous insertions are selected against. With reciprocal crossing experiments and fluorescent in-situ hybridisation of Cla-elements to polytene chromosomes, we documented phenotypic effects on female fertility and chromosomal mispairings that might be linked to DMI in hybrids. We propose that the inferred negative selection on heterozygous Cla-element insertions causes endogenous genetic barriers and therefore acts as DMI among C. riparius populations. The intrinsic genomic turnover exerted by TEs, thus, may have a direct impact on population divergence that is operationally different from drift and local adaptation.
Bears are iconic mammals with a complex evolutionary history. Natural bear hybrids and studies of few nuclear genes indicate that gene flow among bears may be more common than expected and not limited to the closely related polar and brown bears. Here we present a genome analysis of the bear family with representatives of all living species. Phylogenomic analyses of 869 mega base pairs divided into 18,621 genome fragments yielded a well-resolved coalescent species tree despite signals for extensive gene flow across species. However, genome analyses using three different statistical methods show that gene flow is not limited to closely related species pairs. Strong ancestral gene flow between the Asiatic black bear and the ancestor to polar, brown and American black bear explains numerous uncertainties in reconstructing the bear phylogeny. Gene flow across the bear clade may be mediated by intermediate species such as the geographically wide-spread brown bears leading to massive amounts of phylogenetic conflict. Genome-scale analyses lead to a more complete understanding of complex evolutionary processes. The increasing evidence for extensive inter-specific gene flow, found also in other animal species, necessitates shifting the attention from speciation processes achieving genome-wide reproductive isolation to the selective processes that maintain species divergence in the face of gene flow.
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.
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. By applying whole genome analyses and genotype-environment associations to populations of the main dengue vector Aedes aegypti, sampled along an altitudinal gradient in Nepal (200–1300 m), we identify putatively adaptive traits and describe the species' genomic footprint of climate adaptation to colder ecoregions. We found two differentiated clusters with significantly different allele frequencies in genes associated to climate adaptation between the highland population (1300 m) and all other lowland populations (≤800 m). We revealed nonsynonymous mutations in 13 of the candidate genes associated to either altitude, precipitation or cold tolerance and identified an isolation-by-environment differentiation pattern. Other than the expected gradual differentiation along the altitudinal gradient, our results reveal a distinct genomic differentiation of the highland population. Local high-altitude adaptation could be one explanation of the population's phenotypic cold tolerance. Carrying alleles relevant for survival under colder climate increases the likelihood of this highland population to a worldwide expansion into other colder ecoregions.
Highlights
• BaP exposure increases the mutation rate of C. riparius.
• BaP exposure is detrimental for the fitness and the population dynamics of C. riparius.
• Multi-generational studies are essential to assess evolutionary implications of anthropogenic substances on biodiversity.
Abstract
The release of polycyclic aromatic hydrocarbons (PAHs) into the environment is posing a threat to ecosystems and human health. Benzo(a)pyrene (BaP) is considered a biomarker of PAH exposure and is classified as a Group 1 carcinogen. However, it was not known whether BaP is mutagenic, i.e. induces inherited germline mutations. In this study, we used a recently established method, which combines short-term mutation accumulation lines (MAL) with whole genome sequencing (WGS) to assess mutagenicity in the non-biting midge Chironomus riparius. The mutagenicity analysis was supplemented by an evaluation of the development of population fitness in three successive generations in the case of chronic exposure to BaP at a high concentration (100 μg/L). In addition, the level of ROS-induced oxidative stress was examined in vivo. Exposure to the higher BaP concentration led to an increase in germline mutations relative to the control, while the lower concentration showed no mentionable effect. Against expectations, BaP exposure decreased ROS-level compared to the control and is thus probably not responsible for the increased mutation rate. Likewise, the higher BaP concentration decreased fitness measured as population growth rate per day (PGR) significantly over all generations, without signs of rapid evolutionary adaptations. Our results thus highlighted that high BaP exposure may influence the evolutionary trajectory of organisms.
Background: Many fungal species occur across a variety of habitats. Particularly lichens, fungi forming symbioses with photosynthetic partners, have evolved remarkable tolerances for environmental extremes. Despite their ecological importance and ubiquity, little is known about the genetic basis of adaption in lichen populations. Here we studied patterns of genome-wide differentiation in the lichen-forming fungus Lasallia pustulata along an altitudinal gradient in the Mediterranean region. We resequenced six populations as pools and identified highly differentiated genomic regions. We then detected gene-environment correlations while controlling for shared population history and pooled sequencing bias, and performed ecophysiological experiments to assess fitness differences of individuals from different environments.
Results: We detected two strongly differentiated genetic clusters linked to Mediterranean and temperate-oceanic climate, and an admixture zone, which coincided with the transition between the two bioclimates. High altitude individuals showed ecophysiological adaptations to wetter and more shaded conditions. Highly differentiated genome regions contained a number of genes associated with stress response, local environmental adaptation, and sexual reproduction.
Conclusions: Taken together our results provide evidence for a complex interplay between demographic history and spatially varying selection acting on a number of key biological processes, suggesting a scenario of ecological speciation.