Institut für Ökologie, Evolution und Diversität
Refine
Year of publication
Document Type
- Preprint (30) (remove)
Language
- English (30) (remove)
Has Fulltext
- yes (30)
Is part of the Bibliography
- no (30)
Keywords
- Chiton (1)
- Chocó rainforest (1)
- East Africa (1)
- Ecuador (1)
- Giraffa (1)
- Mollusca (1)
- Phylogenomics (1)
- Polyplacophora (1)
- aboveground biomass (1)
- biodiversity (1)
Plant communities provide floral resource-landscapes for pollinators. Yet, it is insufficiently understood how these landscapes shape pollinator-mediated interactions among multiple plant species. Here, we study how pollinators and the seed set of plants respond to the distribution of a floral resource (nectar sugar) in space and across plant species, inflorescences and flowering phenologies. In a global biodiversity hotspot, we quantified floral resource-landscapes on 27 sites of 4 ha comprising 127,993 shrubs of 19 species. Visitation rates of key bird pollinators strongly depended on the phenology of site-scale resource amounts. Seed set of focal plants increased with resources of conspecific neighbours and with site-scale resources, notably with heterospecific resources of lower quality (less sugar per inflorescence). Floral resources are thus a common currency determining how multiple plant species interact via pollinators. These interactions may alter conditions for species coexistence in plant communities and cause community-level Allee effects that promote extinction cascades.
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.
Phylogenetic analyses of nuclear and mitochondrial genomes have shown that polar bears captured the mitochondrial genome of brown bears some 160,00 years ago. This hybridization event likely led to an extinction of the original polar bear mitochondrial genome. However, parts of the mitochondrial DNA occasionally integrates into the nuclear genome, forming pseudogenes called numts (nuclear mitochondrial integrations). Screening the polar bear genome for numts, we identified only 13 such integrations. Analyses of whole-genome sequences from additional polar bears, brown and American black bears as well as the giant panda indicates that the discovered numts entered the bear lineage before the initial ursid radiation some 14 million years ago. Our findings suggests a low integration rate of numts in the bear lineage and a complete loss of the original polar bear mitochondrial genome.
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.
Compared to sequence analyses, phylogenetic reconstruction from transposable elements (TEs) offers an additional perspective to study evolutionary processes. However, detecting phylogenetically informative TE insertions requires tedious experimental work, limiting the power of phylogenetic inference. Here, we analyzed the genomes of seven bear species using high throughput sequencing data to detect thousands of TE insertions. The newly developed pipeline for TE detection called TeddyPi (TE detection and discovery for Phylogenetic Inference) obtained 150,513 high-quality TE insertions in the genomes of ursine and tremarctine bears. By integrating different TE insertion callers and using a stringent filtering approach, the TeddyPi pipeline produced highly reliable TE insertion calls, which were confirmed by extensive in vitro validation experiments. Screening for single nucleotide substitutions in the flanking regions of the TEs show that these substitutions correlate with the phylogenetic signal from the TE insertions. Our phylogenomic analyses show that TEs are a major driver of genomic variation in bears and enabled phylogenetic reconstruction of a well-resolved species tree, even with strong signals for incomplete lineage sorting and introgression. The analyses show that the Asiatic black, sun and sloth bear form a monophyletic clade. TeddyPi is open source and can be adapted to various TE and structural variation callers. The pipeline makes it easy to confidently extract thousands of TE insertions even from low coverage genomes of non-model organisms, opening new possibilities for biologists to study phylogenies, evolutionary processes as well as rates and patterns of (retro-)transposition and structural variation.
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.
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.
The basidiomycete smut fungi are predominantly plant parasitic, causing severe losses in some crops. Most species feature a saprotrophic haploid yeast stage, and several smut fungi are only known from this stage, with some isolated from habitats without suitable hosts, e.g. from Antarctica. Thus, these species are generally believed to be apathogenic, but recent findings that some of these might have a plant pathogenic sexual counterpart, casts doubts on the validity of this hypothesis. Here, four Pseudozyma genomes were re-annotated and compared to published smut pathogens and the well-characterised effector gene Pep1 from these species was checked for its ability to complement a Pep1 deletion strain of Ustilago maydis. It was found that 113 high-confidence putative effector proteins were conserved among smut and Pseudozyma genomes. Among these were several validated effector proteins, including Pep1. By genetic complementation we show that Pep1 homologs from the supposedly apathogenic yeasts restore virulence in Pep1-deficient mutants Ustilago maydis. Thus, it is concluded that Pseudozyma species have retained a suite of effectors. This hints at the possibility that Pseudozyma species have kept an unknown plant pathogenic stage for sexual recombination or that these effectors have positive effects when colonising plant surfaces.
All giraffe (Giraffa) were previously assigned to a single species (G. Camelopardalis) and nine subspecies. However, multi-locus analyses of all subspecies have shown that there are four genetically distinct clades and suggest four giraffe species. This conclusion might not be fully accepted due to limited data and lack of explicit gene flow analyses. Here we present an extended study based on 21 independent nuclear loci from 137 individuals. Explicit gene flow analyses identify less than one migrant per generation, including between the closely related northern and reticulated giraffe. Thus, gene flow analyses and population genetics of the extended dataset confirm four genetically distinct giraffe clades and support four independent giraffe species. The new findings call for a revision of the IUCN classification of giraffe taxonomy. Three of the four species are threatened with extinction, mostly occurring in politically unstable regions, and as such, require the highest conservation support possible.
Consistent individual differences in behaviour (animal personality) are widespread throughout the Animal Kingdom. This includes variation in risk-taking versus risk-averse behavioural tendencies. Variation in several personality dimensions is associated with distinct fitness consequences and thus, may become a target of natural and/or sexual selection. However, the link between animal personality and mate choice—as a major component of sexual selection—remains understudied. We asked (1) whether females and males of the livebearing fish Poecilia mexicana prefer risk-taking mating partners (directional mating preference), (2) or if their preferences are dependent on the choosing individual’s own personality type (assortative mating). We characterized each test subject for its risk-taking behaviour, assessed as the time to emerge from shelter and enter an unknown area. In dichotomous association preference tests, we offered two potential mating partners that differed in risk-taking behaviour but were matched for other phenotypic traits (body size, shape, and colouration). Females, but not males, exhibited a strong directional preference for risk-taking over risk-averse mating partners. At the same time, the strength of females’ preferences correlated positively with their own risk-taking scores. Our study is the first to demonstrate that a strong overall preference for risk-taking mating partners does not preclude effects of choosing individuals’ own personality type on (subtle) individual variation in mating preferences. More generally, two different preferences functions appear to interact to determine the outcome of individual mate choice decisions.
Ongoing climate change is a major threat to biodiversity and impacts on species distributions and abundances are already evident. Heterogenous responses of species due to varying abiotic tolerances and dispersal abilities have the potential to further amplify or ameliorate these impacts through changes in species assemblages. Here we investigate the impacts of climate change on terrestrial bird distributions and, subsequently, on species richness as well as on different aspects of phylogenetic diversity of species assemblages across the globe. We go beyond previous work by disentangling the potential impacts on assemblage phylogenetic diversity of species gains vs. losses under climate change and compare the projected impacts to randomized assemblage changes.
We show that climate change might not only affect species numbers and composition of global species assemblages but could also have profound impacts on assemblage phylogenetic diversity, which, across extensive areas, differ significantly from random changes. Both the projected impacts on phylogenetic diversity and on phylogenetic structure vary greatly across the globe. Projected increases in the evolutionary history contained within species assemblages, associated with either increasing phylogenetic diversification or clustering, are most frequent at high northern latitudes. By contrast, projected declines in evolutionary history, associated with increasing phylogenetic over-dispersion or homogenisation, are projected across all continents.
The projected widespread changes in the phylogenetic structure of species assemblages show that changes in species richness do not fully reflect the potential threat from climate change to ecosystems. Our results indicate that the most severe changes to the phylogenetic diversity and structure of species assemblages are likely to be caused by species range shifts rather than range reductions and extinctions. Our findings highlight the importance of considering diverse measures in climate impact assessments and the value of integrating species-specific responses into assessments of entire community changes.
Wildfires are relatively rare in subarctic tundra ecosystems, but they can strongly change ecosystem properties. Short-term fire effects on subarctic tundra vegetation are well documented, but long-term vegetation recovery has been studied less. The frequency of tundra fires will increase with climate warming. Understanding the long-term effects of fire is necessary to predict future ecosystem changes.
We used a space-for-time approach to assess vegetation recovery after fire over more than four decades. We studied soil and vegetation patterns on three large fire scars (>44, 28 and 12 years old) in dry, lichen-dominated forest tundra in Western Siberia. On 60 plots, we determined soil temperature and permafrost thaw depth, sampled vegetation and measured plant functional traits. We assessed trends in NDVI to support the field-based results on vegetation recovery.
Soil temperature, permafrost thaw depth and total vegetation cover had recovered to pre-fire levels after >44 years, as well as total vegetation cover. In contrast, after >44 years, functional groups had not recovered to the pre-fire state. Burnt areas had lower lichen and higher bryophyte and shrub cover. The dominating shrub species, Betula nana, exhibited a higher vitality (higher specific leaf area and plant height) on burnt compared with control plots, suggesting a fire legacy effect in shrub growth. Our results confirm patterns of shrub encroachment after fire that were detected before in other parts of the Arctic and Subarctic. In the so far poorly studied Western Siberian forest tundra we demonstrate for the first time, long-term fire-legacies on the functional composition of relatively dry shrub- and lichen-dominated vegetation.
The plant family Brassicaceae includes some of the most studied hosts of plant microbiomes, targeting microbial diversity, community assembly rules, and effects on host performance. Compared to bacteria, eukaryotes in the brassicaceous microbiome remain understudied, especially under natural settings. Here, we assessed the impact of host identity and age on the assembly of fungal and oomycete root communities, using DNA metabarcoding of roots and associated soil of three annual co-habiting Brassicaceae collected at two time points. Our results showed that fungal communities are more diverse and structured than those of oomycetes. In both cases, plant identity and sampling time had little influence on community variation, whereas root/soil compartment had a strong effect by exerting control on the entry of soil microorganisms into the roots. The enrichment in roots of specific fungi suggests a specialization towards the asymptomatic colonization of plant tissues, which could be relevant to host’s fitness and health.
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 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.
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.
Feeding exclusively on blood, vampire bats represent the only obligate sanguivorous lineage among mammals. To uncover genomic changes associated with adaptations to this unique dietary specialization, we generated a new haplotype-resolved reference-quality genome of the common vampire bat (Desmodus rotundus) and screened 26 bat species for genes that were specifically lost in the vampire bat lineage. We discovered previously-unknown gene losses that relate to metabolic and physiological changes, such as reduced insulin secretion (FFAR1, SLC30A8), limited glycogen stores (PPP1R3E), and a distinct gastric physiology (CTSE). Other gene losses likely reflect the biased nutrient composition (ERN2, CTRL) and distinct pathogen diversity of blood (RNASE7). Interestingly, the loss of REP15 likely helped vampire bats to adapt to high dietary iron levels by enhancing iron excretion and the loss of the 24S-hydroxycholesterol metabolizing enzyme CYP39A1 could contribute to their exceptional cognitive abilities. Finally, losses of key cone phototransduction genes (PDE6H, PDE6C) suggest that these strictly-nocturnal bats completely lack cone-based vision. These findings enhance our understanding of vampire bat biology and the genomic underpinnings of adaptations to sanguivory.
Background Transposable elements (TEs) are an important source of genome plasticity across the tree of life. Accumulating evidence suggests that TEs may not be randomly distributed in the genome. Drift and natural selection are important forces shaping TE distribution and accumulation, acting directly on the TE element or indirectly on the host species. Fungi, with their multifaceted phenotypic diversity and relatively small genome size, are ideal models to study the role of TEs in genome evolution and their impact on the host’s ecological and life history traits. Here we present an account of all TEs found in a high-quality reference genome of the lichen-forming fungus Umbilicaria pustulata, a macrolichen species comprising two climatic ecotypes: Mediterranean and cold-temperate. We trace the occurrence of the newly identified TEs in populations along three replicated elevation gradients using a Pool-Seq approach, to identify TE insertions of potential adaptive significance.
Results We found that TEs cover 21.26 % of the 32.9 Mbp genome, with LTR Gypsy and Copia clades being the most common TEs. Out of a total of 182 TE copies we identified 28 insertions displaying consistent insertion frequency differences between the two host ecotypes across the elevation gradients. Most of the highly differentiated insertions were located near genes, indicating a putative function.
Conclusions This pioneering study into the content and climate niche-specific distribution of TEs in a lichen-forming fungus contributes to understanding the roles of TEs in fungal evolution. Particularly, it may serve as a foundation for assessing the impact of TE dynamics on fungal adaptation to the abiotic environment, and the impact of TE activity on the evolution and maintenance of a symbiotic lifestyle.