Joshua S. Rest

As species diverge, a wide range of evolutionary processes lead to changes in protein–protein interaction (PPI) networks and metabolic networks. The rate at which molecular networks evolve is an important question in evolutionary biology. Previous empirical work has focused on interactomes from model organisms to calculate rewiring rates, but this is limited by the relatively small number of species and sparse nature of network data across species. We present a proxy for variation in network topology: variation in drug–drug interactions (DDIs), obtained by studying drug combinations (DCs) across taxa. Here, we propose the rate at which DDIs change across species as an estimate of the rate at which the underlying molecular network changes as species diverge. We computed the evolutionary rates of DDIs using previously published data from a high-throughput study in gram-negative bacteria. Using phylogenetic comparative methods, we found that DDIs diverge rapidly over short evolutionary time periods, but that divergence saturates over longer time periods. In parallel, we mapped drugs with known targets in PPI and cofunctional networks. We found that the targets of synergistic DDIs are closer in these networks than other types of DCs and that synergistic interactions have a higher evolutionary rate, meaning that nodes that are closer evolve at a faster rate. Future studies of network evolution may use DC data to gain larger-scale perspectives on the details of network evolution within and between species.

We used long-read sequencing to produce a telomere-to-telomere genome assembly for the heterotrophic stramenopile protist Aurantiochytrium limacinum MYA-1381. Its ∼62 Mbp nuclear genome comprises 26 linear chromosomes with a novel configuration: subtelomeric rDNAs are interspersed with long repeated sequence elements denoted as LOng REpeated - TElomere And Rdna Spacers (LORE-TEARS). These repeats may play a role in chromosome end maintenance. A ∼300 Kbp circular herpesvirus-like genomic element is present at a high copy number. A 269 Kbp related virus-like element was found to reside between two complete sets of rRNA and LORE-TEAR sequences on one end of chromosome 15, indicating recent recombination between the viral and nuclear genome. Our data reveal new types of giant endogenous viral elements originating from herpes-like viruses and existing as either ‘stand-alone’ or integrated elements.

As species diverge, a wide range of evolutionary processes lead to changes in protein-protein interaction networks and metabolic networks. The rate at which biological networks evolve is an important question in evolutionary biology. Previous empirical work has focused on interactomes from model organisms to calculate rewiring rates, but this is limited by the relatively small number of species and sparse nature of network data across species. We present a proxy for variation in network topology: variation in drug-drug interactions (DDIs), obtained by studying drug combinations (DCs) across taxa. Here, we propose the rate at which DDIs change across species as an estimate of the rate at which the underlying biological network changes as species diverge. We computed the evolutionary rates of DDIs using previously published data from a high throughput study in gram-negative bacteria. Using phylogenetic comparative methods, we found that DDIs diverge rapidly over short evolutionary time periods, but that divergence saturates over longer time periods. In parallel, we mapped drugs with known targets in protein-protein interaction and co-functional networks. We found that the targets of synergistic DDIs are closer in these networks than other types of DCs and that synergistic interactions have a higher evolutionary rate, meaning that nodes that are closer evolve at a faster rate. Future studies of network evolution may use DC data to gain larger-scale perspectives on the details of network evolution within and between species.

Thraustochytrids (phylum: Labyrinthulomycota) are non-photosynthetic marine protists. Some thraustochytrids have crtIBY, a trifunctional fusion gene encoding a protein capable of β-carotene biosynthesis from geranylgeranyl pyrophosphate. Here we show that crtIBY is essential in, and encodes the sole pathway for, carotenoid biosynthesis in the thraustochytrid Aurantiochytrium limacinum ATCC MYA-1381. We explore the evolutionary origins of CrtIBY and discover that the closest related protein domains are present in a small but diverse group of other heterotrophic protists, including the apusomonad Thecamonas trahens and the dinoflagellates Oxyrrhis marina and Noctiluca scintillans. Each organism within this cluster also contains one or more β-carotene 15-15’ oxygenase genes (blh and rpe65), suggesting that acquisition of β-carotene biosynthesis genes may have been related to the production of retinal. Our findings support a novel origin of eukaryotic (apo)carotenoid biosynthesis by horizontal gene transfer from Actinobacteria, Bacteroidetes, and/or Archaea. This reveals a remarkable case of parallel evolution of eukaryotic (apo)carotenogenesis in divergent protistan lineages by repeated gene transfers.

There is now abundant evidence of rapid evolution in natural populations, but the genetic mechanisms of these changes remain unclear. One possible route to rapid evolution is through changes in the expression of genes that influence traits under selection. We examined contemporary evolutionary gene expression changes in plant populations responding to environmental fluctuations. We compared genome‐wide gene expression, using RNA‐seq, in two populations of Brassica rapa collected over four time points between 1997 and 2014, during which precipitation in southern California fluctuated dramatically and phenotypic and genotypic changes occurred. By combining transcriptome profiling with the resurrection approach, we directly examined evolutionary changes in gene expression over time. For both populations, we found a substantial number of differentially expressed genes between generations, indicating rapid evolution in the expression of many genes. Using existing gene annotations, we found that many changes occurred in genes involved in regulating stress responses and flowering time. These appeared related to the fluctuations in precipitation and were potentially adaptive. However, the evolutionary changes in gene expression differed across generations within and between populations, indicating largely independent evolutionary trajectories across populations and over time. Our study provides strong evidence for rapid evolution in gene expression, and indicates that changes in gene expression can be one mechanism of rapid evolutionary responses to selection episodes. This study also illustrates that combining resurrection studies with transcriptomics is a powerful approach for investigating evolutionary changes at the gene regulatory level, and will provide new insights into the genetic basis of contemporary evolution.

Marine protists are a polyphyletic group of organisms playing major roles in the ecology and biogeochemistry of the oceans, including performing much of Earth’s photosynthesis and driving the carbon, nitrogen, and silicon cycles. In addition, marine protists occupy key positions in the tree of life, including as the closest relatives of metazoans. Despite all the reasons to better understand them, knowledge of the cell biology of most marine protist lineages is sparse. This is beginning to change thanks to vibrant growth in the development of new model organisms. Here, we survey some recent advances in studying the cell biology of marine protists toward understanding the functional basis of their unique features, gaining new perspectives on universal eukaryotic biology, and for understanding homologous biology within metazoans and the evolution of metazoan traits.

The purpose of this protocol is to evaluate the significance of delta G values, as a strategy to determine the significance of secondary structure. Sequences are shuffled randomly 100x, maintaining base composition but not base order. The percentile of the delta G vlaue of the non-scrambled(original) sequence is obtained from the distribution created by the deltaG values of the 100X scrambled sequences.

Summary: The large amounts of transcriptome data available for Arabidopsis thaliana make a compelling case for the need to generalize results across studies and extract the most robust and meaningful information possible from them. The results of various studies seeking to identify water stress‐responsive genes only partially overlap. The aim of this work was to combine transcriptomic studies in a systematic way that identifies commonalities in response, taking into account variation among studies due to batch effects as well as sampling variation, while also identifying the effect of study‐specific variables, such as the method of applying water stress, and the part of the plant the mRNA was extracted from. We used meta‐analysis, the quantitative synthesis of independent research results, to summarize expression responses to water stress across studies, and meta‐regression to model the contribution of covariates that may affect gene expression. We found that some genes with small but consistent differential responses become evident only when results are synthesized across experiments, and are missed in individual studies. We also identified genes with expression responses that are attributable to use of different plant parts and alternative methods for inducing water stress. Our results indicate that meta‐analysis and meta‐regression provide a powerful approach for identifying a robust gene set that is less sensitive to idiosyncratic results and for quantifying study characteristics that result in contrasting gene expression responses across studies. Combining meta‐analysis with individual analyses may contribute to a richer understanding of the biology of water stress responses, and may prove valuable in other gene expression studies.

Recent studies have demonstrated adaptive evolutionary responses to climate change, but little is known about how these responses may influence ecological interactions with other organisms, including natural enemies. We used a resurrection experiment in the greenhouse to examine the effect of evolutionary responses to drought on the susceptibility of Brassica rapa plants to a fungal pathogen, Alternaria brassicae . In agreement with previous studies in this population, we found an evolutionary shift to earlier flowering postdrought, which was previously shown to be adaptive. Here, we report the novel finding that postdrought descendant plants were also more susceptible to disease, indicating a rapid evolutionary shift to increased susceptibility. This was accompanied by an evolutionary shift to increased specific leaf area (thinner leaves) following drought. We found that flowering time and disease susceptibility displayed plastic responses to experimental drought treatments, but that this plasticity did not match the direction of evolution, indicating that plastic and evolutionary responses to changes in climate can be opposed. The observed evolutionary shift to increased disease susceptibility accompanying adaptation to drought provides evidence that even if populations can rapidly adapt in response to climate change, evolution in other traits may have ecological effects that could make species more vulnerable.