Refine
Year of publication
Document Type
- Preprint (30) (remove)
Language
- English (30)
Has Fulltext
- yes (30)
Is part of the Bibliography
- no (30) (remove)
Keywords
- Chiton (1)
- Chocó rainforest (1)
- East Africa (1)
- Ecuador (1)
- Giraffa (1)
- Mollusca (1)
- Phylogenomics (1)
- Polyplacophora (1)
- aboveground biomass (1)
- biodiversity (1)
Institute
- Institut für Ökologie, Evolution und Diversität (30) (remove)
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.
Seed harvesting from wild plant populations is key for ecological restoration, but may threaten the persistence of source populations. Consequently, several countries have set guidelines limiting the proportions of harvestable seeds. Here, we use high-resolution data from 298 plant species to model the demographic consequences of seed harvesting. We find that the current guidelines only protect some species, but are insufficient or overly restrictive for others. We show that the maximum possible fraction of seed harvesting is strongly associated with harvesting frequency and generation time of the target species, ranging from 100% in long-lived species to <1% in the most annuals. Our results provide quantitative basis to guide seed harvesting legislation based on species’ generation time and harvesting regime.
Background: In the speciation continuum the strength of reproductive isolation varies, and species boundaries are blurred by gene flow. Interbreeding among giraffe (Giraffa spp.) in captivity is known and anecdotal reports of natural hybrids exist. In Kenya, Nubian (G. camelopardalis camelopardalis), reticulated (G. reticulata), and Masai giraffe sensu stricto (G. tippelskirchi tippelskirchi) are parapatric, and thus the country might be a melting pot for these taxa. We analyzed 128 genomes of wild giraffe, 113 newly sequenced, representing these three taxa.
Results: We found varying levels of Nubian ancestry in 13 reticulated giraffe sampled across the Laikipia Plateau most likely reflecting historical gene flow between these two lineages. Although comparatively weaker signs of ancestral gene flow and potential mitochondrial introgression from reticulated into Masai giraffe were also detected, estimated admixture levels between these two lineages are minimal. Importantly, contemporary gene flow between East African giraffe lineages was not statistically significant. Effective population sizes have declined since the Late Pleistocene, more severely for Nubian and reticulated giraffe.
Conclusions: Despite historically hybridizing, these three giraffe lineages have maintained their overall genomic integrity suggesting effective reproductive isolation, consistent with the previous classification of giraffe into four species.
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.
Tree bark constitutes ideal habitat for microbial communities, because it is a stable substrate, rich in micro-niches. Bacteria, fungi, and terrestrial microalgae together form microbial communities, which in turn support more bark-associated organisms, such as mosses, lichens, and invertebrates, thus contributing to forest biodiversity. We have a limited understanding of the diversity and biotic interactions of the bark-associated microbiome, as investigations have mainly focussed on agriculturally relevant systems and on single taxonomic groups. Here we implemented a multi-kingdom metabarcoding approach to analyse diversity and community structure of the green algal, bacterial, and fungal components of the bark-associated microbial communities of beech, the most common broadleaved tree of Central European forests. We identified the most abundant taxa, hub taxa, and co-occurring taxa. We found that tree size (as a proxy for age) is an important driver of community assembly, suggesting that environmental filtering leads to less diverse fungal and algal communities over time. Conversely, forest management intensity had negligible effects on microbial communities on bark. Our study suggests the presence of undescribed, yet ecologically meaningful taxa, especially in the fungi, and highlights the importance of bark surfaces as a reservoir of microbial diversity. Our results constitute a first, essential step towards an integrated framework for understanding microbial community assembly processes on bark surfaces, an understudied habitat and neglected component of terrestrial biodiversity. Finally, we propose a cost-effective sampling strategy to study bark-associated microbial communities across large spatial or environmental scales.
From hunting and foraging to clearing land for agriculture, humans modify forest biodiversity, landscapes, and climate. Forests constantly undergo disturbance–recovery dynamics and understanding them is a major objective of ecologists and conservationists. Chronosequences are a useful tool for understanding global restoration efforts. They represent a space-for-time substitution approach suited for the quantification of the resistance of ecosystem properties to withstand disturbance and the resilience of these properties until reaching pre-disturbance levels. Here we introduce a newly established chronosequence with 62 plots (50 ⍰ 50 m) in active cacao plantations and pastures, early and late regeneration, and mature old-growth forests, across a 200 km2 area in the extremely wet Chocó rainforest. Our chronosequence covers by far the largest total area of plots compared to others in the Neotropics. Plots ranged from 159–615 masl in a forested landscape with 74 ± 2.8 % forest cover within a 1-km radius including substantial old-growth forest cover. Land-use legacy and regeneration time were not confounded by elevation. We tested how six forest structure variables (maximum tree height and DBH, basal area, number of stems, vertical vegetation heterogeneity, and light availability), aboveground biomass (AGB), and rarefied tree species richness change along our chronosequence. Forest structure variables, AGB, and tree species richness increased with regeneration time and are predicted to reach similar levels to those in old-growth forests after ca. 30–116, 202, and 108 yrs, respectively. Compared to previous work in the Neotropics, old-growth forests in Canandé accumulate high AGB that takes one of the largest time spans reported until total recovery. Our chronosequence comprises one of the largest tree species pools, covers the largest total area of regenerating and old-growth forests, and has higher forest cover than other Neotropical chronosequences. Hence, our chronosequence can be used to determine the time for recovery and stability (resistance and resilience) of different taxa and ecosystem functions, including species interaction networks. This integrative effort will ultimately help to understand how one of the most diverse forests on the planet recovers from large-scale disturbances.
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.
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.
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.
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.
Mollusca is the second-largest animal phylum with over 100,000 species among eight distinct taxonomic classes. Across 1000 living species in the class Polyplacophora, chitons have a relatively constrained morphology but with some notable deviations. Several genera possess “shell eyes”, true eyes with a lens and retina that are embedded within the dorsal shells, which represent the most recent evolution of animal eyes. The phylogeny of major chiton clades is mostly well established, in a set of superfamily and higher-level taxa supported by various approaches including multiple gene markers, mitogenome-phylogeny and phylotranscritomic approaches as well as morphological studies. However, one critical lineage has remained unclear: Schizochiton was controversially suggested as a potential independent origin of chiton shell eyes. Here, with the draft genome sequencing of Schizochiton incisus (superfamily Schizochitonoidea) plus assembly of transcriptome data from other polyplacophorans, we present phylogenetic reconstructions using both mitochondrial genomes and phylogenomic approaches with multiple methods. Phylogenetic trees from mitogenomic data are inconsistent, reflecting larger scale confounding factors in molluscan mitogenomes. A consistent robust topology was generated with protein coding genes using different models and methods. Our results support Schizochitonoidea is a sister group to other Chitonoidea in Chitonina, in agreement with established classification. This suggests that the earliest origin of shell eyes is in Schizochitonoidea, which were also gained secondarily in other genera in Chitonoidea. Our results have generated a holistic review of the internal relationship within Polyplacophora, and a better understanding on the evolution of Polyplacophora.
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.
Intraspecific genomic variability affects a species’ adaptive potential towards climatic conditions. Variation in gene content across populations and environments may point at genomic adaptations to specific environments. The lichen symbiosis, a stable association of fungal and photobiont partners, offers an excellent system to study environmentally driven gene content variation. Many species have remarkable environmental tolerances, and often form populations in different climate zones. Here we combine comparative and population genomics to assess the presence and absence of genes in high elevation and low elevation genomes of two lichenized fungi of the genus Umbilicaria. The two species have non-overlapping ranges, but occupy similar climatic niches in North America (U. phaea) and Europe (U. pustulata): high elevation populations are located in the cold temperate zone and low elevation populations in the Mediterranean zone. We assessed gene content variation along replicated elevation gradients in each of the two species, based on a total of 2050 individuals across 26 populations. Specifically, we assessed shared orthologs across species within the same climate zone, and tracked which genes increase or decrease in abundance within populations along elevation. In total, we found 16 orthogroups with shared orthologous genes in genomes at low elevation and 13 at high elevation. Coverage analysis revealed one ortholog that is exclusive to genomes at low elevation. Conserved domain search revealed domains common to the protein kinases (PKs) superfamily. We traced the discovered ortholog in populations along five replicated elevation gradients on both continents. The protein kinase gene linearly declined in abundance with increasing elevation, and was absent in the highest populations. We consider the parallel loss of an ortholog in two species and in two geographic settings a rare find, and a step forward in understanding the genomic underpinnings of climatic tolerances in lichenized fungi. In addition, the tracking of gene content variation provides a widely applicable framework for retrieving biogeographical determinants of gene presence/absence patterns. Our work provides insights into gene content variation of lichenized fungi in relation to climatic gradients, suggesting a new research direction with implications for understanding evolutionary trajectories of complex symbioses in relation to climatic change.
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.
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.
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.
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.
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.
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.
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.