Ecological Diversification into Broad Thermal Niches Revealed by Protein Resurrection and Proteome-Wide Ancestral Reconstruction During the Rapid Radiation of Alvinellid Worms

In the deep Pacific and Indian oceans, Alvinellid worms diversified about 100 million years ago to colonize a variety of hydrothermal vent environments. It has been suggested that the last common ancestor of this family was a thermophilic species. However, the evolutionary history of these worms is complex, with putative gene flow among ancestors. In this study, we investigated the early evolution of the family in relation to the diversification of thermal niches. We demonstrated that phylogenetic histories of alvinellid species greatly vary along chromosomes, possibly due to allele introgression between nascent ecotypes. We performed sequence reconstructions of ancestral cytosolic MDH and Cu/Zn SOD under the two most frequent phylogenies encountered along the genome. In silico simulations of folding stability for these two enzymes were highly correlated with their biophysical characterization (micro-calorimetry and differential scanning fluorimetry), and were used to predict the thermostability of thousands of alternatively reconstructed proteins. We further generalized the reconstruction at the proteome scale, taking the amino-acid usage bias as a proxy for folding stability in a new phylogenetic model accounting for amino-acid variations over time. Both approaches agree that the last common ancestors of Alvinellidae gained highly stable proteins, comparable to proteins of present-day thermotolerant species. Our computational protocol allowed us to validate this scenario under multiple phylogenetic hypotheses and different sequence reconstruction models. Considered together, Alvinellidae shed light on how metazoan species diversify in order to colonize different thermal niches, combining adaptive mutation and selection on genetic variants.