Abstract
Motivation
Genome sequencing projects have revealed frequent gains and losses of genes between species. Previous versions of our software, Computational Analysis of gene Family Evolution ...(CAFE), have allowed researchers to estimate parameters of gene gain and loss across a phylogenetic tree. However, the underlying model assumed that all gene families had the same rate of evolution, despite evidence suggesting a large amount of variation in rates among families.
Results
Here, we present CAFE 5, a completely re-written software package with numerous performance and user-interface enhancements over previous versions. These include improved support for multithreading, the explicit modeling of rate variation among families using gamma-distributed rate categories, and command-line arguments that preclude the use of accessory scripts.
Availability and implementation
CAFE 5 source code, documentation, test data and a detailed manual with examples are freely available at https://github.com/hahnlab/CAFE5/releases.
Supplementary information
Supplementary data are available at Bioinformatics online.
The idea that species boundaries can be semipermeable to gene flow is now widely accepted but the evolutionary importance of introgressive hybridization remains unclear. Here we examine the genomic ...contribution of gene flow between two hybridizing chipmunk species, Tamias ruficaudus and T. amoenus. Previous studies have shown that ancient hybridization has resulted in complete fixation of introgressed T. ruficaudus mitochondrial DNA (mtDNA) in some populations of T. amoenus, but the extent of nuclear introgression is not known. We used targeted capture to sequence over 10,500 gene regions from multiple individuals of both species. We found that most of the nuclear genome is sorted between these species and that overall genealogical patterns do not show evidence for introgression. Our analysis rules out all but very minor levels of interspecific gene flow, indicating that introgressive hybridization has had little impact on the overall genetic composition of these species outside of the mitochondrial genome. Given that much of the evidence for introgression in animals has come from mtDNA, our results underscore that unraveling the importance introgressive hybridization during animal speciation will require a genome-wide perspective that is still absent for many species.
Full text
Available for:
BFBNIB, FZAB, GIS, IJS, KILJ, NLZOH, NMLJ, NUK, OILJ, PNG, SAZU, SBCE, SBMB, UL, UM, UPUK
Natural history museum collections provide unique resources for understanding how species respond to environmental change, including the abrupt, anthropogenic climate change of the past century. ...Ideally, researchers would conduct genome‐scale screening of museum specimens to explore the evolutionary consequences of environmental changes, but to date such analyses have been severely limited by the numerous challenges of working with the highly degraded DNA typical of historic samples. Here, we circumvent these challenges by using custom, multiplexed, exon capture to enrich and sequence ~11 000 exons (~4 Mb) from early 20th‐century museum skins. We used this approach to test for changes in genomic diversity accompanying a climate‐related range retraction in the alpine chipmunks (Tamias alpinus) in the high Sierra Nevada area of California, USA. We developed robust bioinformatic pipelines that rigorously detect and filter out base misincorporations in DNA derived from skins, most of which likely resulted from postmortem damage. Furthermore, to accommodate genotyping uncertainties associated with low‐medium coverage data, we applied a recently developed probabilistic method to call single‐nucleotide polymorphisms and estimate allele frequencies and the joint site frequency spectrum. Our results show increased genetic subdivision following range retraction, but no change in overall genetic diversity at either nonsynonymous or synonymous sites. This case study showcases the advantages of integrating emerging genomic and statistical tools in museum collection‐based population genomic applications. Such technical advances greatly enhance the value of museum collections, even where a pre‐existing reference is lacking and points to a broad range of potential applications in evolutionary and conservation biology.
Full text
Available for:
BFBNIB, FZAB, GIS, IJS, KILJ, NLZOH, NUK, OILJ, SAZU, SBCE, SBMB, UL, UM, UPUK
To date, exon capture has largely been restricted to species with fully sequenced genomes, which has precluded its application to lineages that lack high quality genomic resources. We developed a ...novel strategy for designing array-based exon capture in chipmunks (Tamias) based on de novo transcriptome assemblies. We evaluated the performance of our approach across specimens from four chipmunk species.
We selectively targeted 11,975 exons (~4 Mb) on custom capture arrays, and enriched over 99% of the targets in all libraries. The percentage of aligned reads was highly consistent (24.4-29.1%) across all specimens, including in multiplexing up to 20 barcoded individuals on a single array. Base coverage among specimens and within targets in each species library was uniform, and the performance of targets among independent exon captures was highly reproducible. There was no decrease in coverage among chipmunk species, which showed up to 1.5% sequence divergence in coding regions. We did observe a decline in capture performance of a subset of targets designed from a much more divergent ground squirrel genome (30 My), however, over 90% of the targets were also recovered. Final assemblies yielded over ten thousand orthologous loci (~3.6 Mb) with thousands of fixed and polymorphic SNPs among species identified.
Our study demonstrates the potential of a transcriptome-enabled, multiplexed, exon capture method to create thousands of informative markers for population genomic and phylogenetic studies in non-model species across the tree of life.
Full text
Available for:
DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK
For over 140 years, lichens have been regarded as a symbiosis between a single fungus, usually an ascomycete, and a photosynthesizing partner. Other fungi have long been known to occur as occasional ...parasites or endophytes, but the one lichen–one fungus paradigm has seldom been questioned. Here we show that many common lichens are composed of the known ascomycete, the photosynthesizing partner, and, unexpectedly, specific basidiomycete yeasts. These yeasts are embedded in the cortex, and their abundance correlates with previously unexplained variations in phenotype. Basidiomycete lineages maintain close associations with specific lichen species over large geographical distances and have been found on six continents. The structurally important lichen cortex, long treated as a zone of differentiated ascomycete cells, appears to consistently contain two unrelated fungi.
Full text
Available for:
BFBNIB, NMLJ, NUK, ODKLJ, PNG, SAZU, UL, UM, UPUK
Bark and ambrosia beetles are highly specialized weevils (Curculionidae) that have established diverse symbioses with fungi, most often from the order Ophiostomatales (Ascomycota, Sordariomycetes). ...The two types of beetles are distinguished by their feeding habits and intimacy of interactions with their symbiotic fungi. The tree tissue diet of bark beetles is facilitated by fungi, while ambrosia beetles feed solely on fungi that they farm. The farming life history strategy requires domestication of a fungus, which the beetles consume as their sole food source. Ambrosia beetles in the subfamily Platypodinae originated in the mid‐Cretaceous (119–88 Ma) and are the oldest known group of farming insects. However, attempts to resolve phylogenetic relationships and the timing of domestication events for fungal cultivars have been largely inconclusive. We sequenced the genomes of 12 ambrosia beetle fungal cultivars and bark beetle associates, including the devastating laurel wilt pathogen, Raffaelea lauricola, to estimate a robust phylogeny of the Ophiostomatales. We find evidence for contemporaneous diversification of the beetles and their associated fungi, followed by three independent domestication events of the ambrosia fungi genus Raffaelea. We estimate the first domestication of an Ophiostomatales fungus occurred ~86 Ma, 25 million years earlier than prior estimates and in close agreement with the estimated age of farming in the Platypodinae (96 Ma). Comparisons of the timing of fungal domestication events with the timing of beetle radiations support the hypothesis that the first large beetle radiations may have spread domesticated “ambrosia” fungi to other fungi‐associated beetle groups, perhaps facilitating the evolution of new farming lineages.
Full text
Available for:
BFBNIB, FZAB, GIS, IJS, KILJ, NLZOH, NUK, OILJ, SAZU, SBCE, SBMB, UL, UM, UPUK
Tree-killing bark beetles are the most economically important insects in conifer forests worldwide. However, despite >200 years of research, the drivers of population eruptions and crashes are still ...not fully understood and the existing knowledge is thus insufficient to face the challenges posed by the Anthropocene. We critically analyze potential biotic and abiotic drivers of population dynamics of an exemplary species, the European spruce bark beetle (ESBB) (Ips typographus) and present a multivariate approach that integrates the many drivers governing this bark beetle system. We call for hypothesis-driven, large-scale collaborative research efforts to improve our understanding of the population dynamics of this and other bark beetle pests. Our approach can serve as a blueprint for tackling other eruptive forest insects.
Bark beetles are currently causing unprecedented damage to European and North American forests.Their population dynamics rarely have been studied in a hypothesis-driven manner incorporating exogenous biotic variables.We propose a conceptual framework to reveal the drivers of bark beetle populations.This approach can be equally applied to other eruptive insect pests.
Full text
Available for:
GEOZS, IJS, IMTLJ, KILJ, KISLJ, NLZOH, NUK, OILJ, PNG, SAZU, SBCE, SBJE, UILJ, UL, UM, UPCLJ, UPUK, ZAGLJ, ZRSKP
Diverse insects are associated with ancient bacterial symbionts, whose genomes have often suffered drastic reduction and degeneration. In extreme cases, such symbiont genomes seem almost unable to ...sustain the basic cellular functioning, which comprises an open question in the evolution of symbiosis. Here,we report an insect group wherein an ancient symbiont lineage suffering massive genome erosion has experienced recurrent extinction and replacement by host-associated pathogenic microbes. Cicadas are associated with the ancient bacterial co-obligate symbionts Sulcia and Hodgkinia, whose streamlined genomes are specialized for synthesizing essential amino acids, thereby enabling the host to live on plant sap. However, our inspection of 24 Japanese cicada species revealed that while all species possessed Sulcia, only nine species retained Hodgkinia, and their genomes exhibited substantial structural instability. The remaining 15 species lacked Hodgkinia and instead harbored yeast-like fungal symbionts. Detailed phylogenetic analyses uncovered repeated Hodgkinia-fungus and fungus-fungus replacements in cicadas. The fungal symbionts were phylogenetically intermingled with cicada-parasitizing Ophiocordyceps fungi, identifying entomopathogenic origins of the fungal symbionts. Most fungal symbionts of cicadas were uncultivable, but the fungal symbiont of Meimuna opalifera was cultivable, possibly because it is at an early stage of fungal symbiont replacement. Genome sequencing of the fungal symbiont revealed its metabolic versatility, presumably capable of synthesizing almost all amino acids, vitamins, and other metabolites, which is more than sufficient to compensate for the Hodgkinia loss. These findings highlight a straightforward ecological and evolutionary connection between parasitism and symbiosis, which may provide an evolutionary trajectory to renovate deteriorated ancient symbiosis via pathogen domestication.
Full text
Available for:
BFBNIB, NMLJ, NUK, PNG, SAZU, UL, UM, UPUK
Maternally transmitted
infect about half of insect species, yet the predominant mode(s) of
acquisition remains uncertain. Species-specific associations could be old, with
and hosts codiversifying (
, ...cladogenic acquisition), or relatively young and acquired by horizontal transfer or introgression. The three
-clade hosts (
,
)
diverged ∼3 MYA and currently hybridize on the West African islands Bioko and São Tomé. Each species is polymorphic for nearly identical
that cause weak cytoplasmic incompatibility (CI)-reduced egg hatch when uninfected females mate with infected males.
-clade
are closely related to
Mel, globally polymorphic in
We use draft
and mitochondrial genomes to demonstrate that
-clade phylogenies for
and mitochondria tend to follow host nuclear phylogenies. However, roughly half of
individuals, sampled both inside and outside of the São Tomé hybrid zone, have introgressed
mitochondria. Both mitochondria and
possess far more recent common ancestors than the bulk of the host nuclear genomes, precluding cladogenic
acquisition. General concordance of
and mitochondrial phylogenies suggests that horizontal transmission is rare, but varying relative rates of molecular divergence complicate chronogram-based statistical tests. Loci that cause CI in
Mel are disrupted in
-clade
; but a second set of loci predicted to cause CI are located in the same WO prophage region. These alternative CI loci seem to have been acquired horizontally from distantly related
, with transfer mediated by flanking
-specific ISWpi1 transposons.
Many species have experienced dramatic changes in their abundance and distribution during recent climate change, but it is often unclear whether such ecological responses are accompanied by ...evolutionary change. We used targeted exon sequencing of 294 museum specimens (160 historic, 134 modern) to generate independent temporal genomic contrasts spanning a century of climate change (1911-2012) for two co-distributed chipmunk species: an endemic alpine specialist (Tamias alpinus) undergoing severe range contraction and a stable mid-elevation species (T. speciosus). Using a novel analytical approach, we reconstructed the demographic histories of these populations and tested for evidence of recent positive directional selection. Only the retracting species showed substantial population genetic fragmentation through time and this was coupled with positive selection and substantial shifts in allele frequencies at a gene, Alox15, involved in regulation of inflammation and response to hypoxia. However, these rapid population and gene-level responses were not detected in an analogous temporal contrast from another area where T. alpinus has also undergone severe range contraction. Collectively, these results highlight that evolutionary responses may be variable and context dependent across populations, even when they show seemingly synchronous ecological shifts. Our results demonstrate that temporal genomic contrasts can be used to detect very recent evolutionary responses within and among contemporary populations, even in the face of complex demographic changes. Given the wealth of specimens archived in natural history museums, comparative analyses of temporal population genomic data have the potential to improve our understanding of recent and ongoing evolutionary responses to rapidly changing environments.
Full text
Available for:
DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK