Title: Phylogenomics reveals the slow-burning fuse of diatom evolution
Evolution is often uneven in its pace and outcomes, with long periods of stasis interrupted by abrupt increases in morphological and ecological disparity. With thousands of gene histories, phylogenomics can uncover the genomic signatures of these broad macroevolutionary trends. Diatoms are a species-rich lineage of microeukaryotes that contribute greatly to the global cycling of carbon, oxygen, and silica, which they use to build elaborately structured cell walls. We combined fossil information with newly sequenced transcriptomes from 181 diverse diatom species to reconstruct the pattern, timing, and genomic context of major evolutionary transitions. Diatoms originated 270 Mya, and after >100 My of relative stasis in morphology and ecology, a radiation near the Jurassic–Cretaceous boundary led to the diversity of habitats and cell wall architectures characteristic of modern diatoms. This transition was marked by a genome duplication and high levels of gene tree discordance. However, short generation times increase the probability of coalescence between speciation events, minimizing the impacts of incomplete lineage sorting and implicating sequence saturation and gene tree error as the main sources of discordance. Nevertheless, a rigorous tree-based approach to ortholog selection resulted in strongly supported relationships, including some that were uncertain previously. Three pulses of accelerated speciation were detected, two of which were associated with the evolution of novel traits and ecological transitions. The first 100 My of diatom evolution was a slow-burning fuse that led to a burst of innovations in ecology, morphology, and life history that are hallmarks of contemporary diatom assemblages.  more » « less
Award ID(s):
1353131 1651087
PAR ID:
10686866
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ; ; ; ; ; ;
Publisher / Repository:
NAS
Date Published:
Journal Name:
Proceedings of the National Academy of Sciences
Volume:
122
Issue:
22
ISSN:
0027-8424
Page Range / eLocation ID:
e2500153122
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Smith, Stacey (Ed.)
    Despite the obstacles facing marine colonists, most lineages of aquatic organisms have colonized and diversified in freshwaters repeatedly. These transitions can trigger rapid morphological or physiological change and, on longer timescales, lead to increased rates of speciation and extinction. Diatoms are a lineage of ancestrally marine microalgae that have diversified throughout freshwater habitats worldwide. We generated a phylogenomic data set of genomes and transcriptomes for 59 diatom taxa to resolve freshwater transitions in one lineage, the Thalassiosirales. Although most parts of the species tree were consistently resolved with strong support, we had difficulties resolving a Paleocene radiation, which affected the placement of one freshwater lineage. This and other parts of the tree were characterized by high levels of gene tree discordance caused by incomplete lineage sorting and low phylogenetic signal. Despite differences in species trees inferred from concatenation versus summary methods and codons versus amino acids, traditional methods of ancestral state reconstruction supported six transitions into freshwaters, two of which led to subsequent species diversification. Evidence from gene trees, protein alignments, and diatom life history together suggest that habitat transitions were largely the product of homoplasy rather than hemiplasy, a condition where transitions occur on branches in gene trees not shared with the species tree. Nevertheless, we identified a set of putatively hemiplasious genes, many of which have been associated with shifts to low salinity, indicating that hemiplasy played a small but potentially important role in freshwater adaptation. Accounting for differences in evolutionary outcomes, in which some taxa became locked into freshwaters while others were able to return to the ocean or become salinity generalists, might help further distinguish different sources of adaptive mutation in freshwater diatoms. 
    more » « less
  2. Abstract Diatoms, a major group of microalgae, play a critical role in global carbon cycling and primary production. Despite their ecological significance, comprehensive genomic resources for diatoms are limited. To address this, we have annotated previously unannotated genome assemblies of 49 diatom species. Genome assemblies were obtained from NCBI Datasets and processed for repeat elements using RepeatModeler2 and RepeatMasker. For gene prediction, BRAKER2 was employed in the absence of transcriptomic data, while BRAKER3 was utilised when transcriptome short read data were available from the Sequence Read Archive. The quality of genome assemblies and predicted protein sets was evaluated using BUSCO, ensuring high-quality genomic resources. Functional annotation was performed using EnTAP, providing insights into the biological roles of the predicted proteins. Our study enhances the genomic toolkit available for diatoms, facilitating future research in diatom biology, ecology, and evolution. 
    more » « less
  3. Piganeau, Gwenael (Ed.)
    Abstract Numerous factors shape the evolution of protein-coding genes, including shifts in the strength or type of selection following gene duplications or changes in the environment. Diatoms and other silicifying organisms use a family of silicon transporters (SITs) to import dissolved silicon from the environment. Freshwaters contain higher silicon levels than oceans, and marine diatoms have more efficient uptake kinetics and less silicon in their cell walls, making them better competitors for a scarce resource. We compiled SITs from 37 diatom genomes to characterize shifts in selection following gene duplications and marine–freshwater transitions. A deep gene duplication, which coincided with a whole-genome duplication, gave rise to two gene lineages. One of them (SIT1–2) is present in multiple copies in most species and is known to actively import silicon. These SITs have evolved under strong purifying selection that was relaxed in freshwater taxa. Episodic diversifying selection was detected but not associated with gene duplications or habitat shifts. In contrast, genes in the second SIT lineage (SIT3) were present in just half the species, the result of multiple losses. Despite conservation of SIT3 in some lineages for the past 90–100 million years, repeated losses, relaxed selection, and low expression highlighted the dispensability of SIT3, consistent with a model of deterioration and eventual loss due to relaxed selection on SIT3 expression. The extensive but relatively balanced history of duplications and losses, together with paralog-specific expression patterns, suggest diatoms continuously balance gene dosage and expression dynamics to optimize silicon transport across major environmental gradients. 
    more » « less
  4. ABSTRACT The application of high‐throughput sequencing to phylogenetic analyses is allowing authors to reconstruct the true evolutionary history of species. This work can illuminate specific mechanisms underlying divergence when combined with analyses of gene flow, recombination and selection. We conducted a phylogenomic analysis ofCatharus, a songbird genus with considerable potential for gene flow, variation in migratory behaviour and genomic resources. We documented discordance among trees constructed for mitochondrial, autosomal and sex (Z) chromosome partitions. Two trees were recovered on the Z. Both trees differed from the autosomes, one matched the mitochondria, and the other was unique to the Z. Gene flow with one species likely generated much of this discordance; substantial admixture betweenustulatusand the remaining species was documented and linked to at least two historic events. The tree unique to the Z likely reflects the true history ofCatharus; local genomic analyses recovered the same tree in autosomal regions with reduced admixture and recombination. Genes previously connected to migration were enriched in these regions suggesting transitions between migratory and non‐migratory states helped generate divergence. Migratory (vs. nonmigratory)Catharusformed a monophyletic clade in a subset of genomic regions. Gene flow was elevated in some of these regions suggesting adaptive introgression may have occurred, but the dominant pattern was of balancing selection maintaining ancestral polymorphisms important for olfaction and perhaps, by extension, adaptation to temperate climates. This work illuminates the evolutionary history of an important model in speciation and demonstrates how differential resistance to gene flow can affect local genomic patterns. 
    more » « less
  5. Nielsen, Rasmus (Ed.)
    Abstract An important goal of evolutionary genomics is to identify genomic regions whose substitution rates differ among lineages. For example, genomic regions experiencing accelerated molecular evolution in some lineages may provide insight into links between genotype and phenotype. Several comparative genomics methods have been developed to identify genomic accelerations between species, including a Bayesian method called PhyloAcc, which models shifts in substitution rate in multiple target lineages on a phylogeny. However, few methods consider the possibility of discordance between the trees of individual loci and the species tree due to incomplete lineage sorting, which might cause false positives. Here, we present PhyloAcc-GT, which extends PhyloAcc by modeling gene tree heterogeneity. Given a species tree, we adopt the multispecies coalescent model as the prior distribution of gene trees, use Markov chain Monte Carlo (MCMC) for inference, and design novel MCMC moves to sample gene trees efficiently. Through extensive simulations, we show that PhyloAcc-GT outperforms PhyloAcc and other methods in identifying target lineage-specific accelerations and detecting complex patterns of rate shifts, and is robust to specification of population size parameters. PhyloAcc-GT is usually more conservative than PhyloAcc in calling convergent rate shifts because it identifies more accelerations on ancestral than on terminal branches. We apply PhyloAcc-GT to two examples of convergent evolution: flightlessness in ratites and marine mammal adaptations, and show that PhyloAcc-GT is a robust tool to identify shifts in substitution rate associated with specific target lineages while accounting for incomplete lineage sorting. 
    more » « less