The Biological Hierarchy, Time, and Temporal 'Omics in Evolutionary Biology: A Perspective

Sequencing data-genomics, transcriptomics, epigenomics, proteomics, and metabolomics-have revolutionized biological research, enabling a more detailed study of processes, ranging from subcellular to evolutionary, that drive biological organization. These processes, collectively, are responsible for generating patterns of phenotypic variation and can operate over dramatically different timescales (milliseconds to billions of years). While researchers often study phenotypic variation at specific levels of biological organization to isolate processes operating at that particular scale, the varying types of sequence data, or 'omics, can also provide complementary inferences to link molecular and phenotypic variation to produce an integrated view of evolutionary biology, ranging from molecular pathways to speciation. We briefly describe how 'omics has been used across biological levels and then demonstrate the utility of integrating different types of sequencing data across multiple biological levels within the same study to better understand biological phenomena. However, single-time-point studies cannot evaluate the temporal dynamics of these biological processes. Therefore, we put forward temporal 'omics as a framework that can better enable researchers to study the temporal dynamics of target processes. Temporal 'omics is not infallible, as the temporal sampling regime directly impacts inferential ability. Thus, we also discuss the role the temporal sampling regime plays in deriving inferences about the environmental conditions driving biological processes and provide examples that demonstrate the impact of the sampling regime on biological inference. Finally, we forecast the future of temporal 'omics by highlighting current methodological advancements that will enable temporal 'omics to be extended across species and timescales. We extend this discussion to using temporal multi-omics to integrate across the biological hierarchy to evaluate and link the temporal dynamics of processes that generate phenotypic variation.

Proteome changes in an aquatic invertebrate consumer in response to different nutritional stressors

Nutrient imbalances in zooplankton are caused by the differences in elemental content of producers and the demand for elements in consumers, which alter the life-history traits in consumers. Changes in life-history traits are mediated through metabolic pathways that affect gene expression and the metabolome. However, less is known about proteomic changes to elemental-limitation in zooplankton. Here, we grew Daphnia pulex under high food quantity and quality (HF), low food quantity (LF), and phosphorus (P)-limited (PL) diets for six days and measured growth, elemental composition, and the proteome. Daphnids in both LF and PL diets grew less. Animals in LF diets had less carbon (C), while daphnids in PL diets had less P compared to HF fed animals. In total, we identified 1719 proteins that were used in a partial least squares regression discriminant analysis (PLS-DA). Focusing on a subset of the proteome, the PLS-DA resulted in a clear separation between animals fed HF diets and PL and LF diets. Many proteome changes in nutrient-limited diets are associated with growth, reproduction, lipid metabolism, and nutrient assimilation. Regardless of the limiting nutrient, there were less hemoglobin and small subunit processome component proteins compared to HF fed animals. Daphnids fed LF diets had less vitellogenin fused superoxide dismutase and more lipid-droplet hydrolase, whereas Daphnia fed PL diets had higher abundances of cytochrome P450 and serine protease. Our proteome results compliment other "omic" studies that could be used to study Daphnia physiology in lakes.

Metabolic Phenotype of Daphnia Under Hypoxia: Macroevolution, Microevolution, and Phenotypic Plasticity

Daphnia is a freshwater crustacean that is able to upregulate hemoglobin (Hb) in response to hypoxia, imparting a red color. We combine multiple field surveys across season with a lab experiment to evaluate changes in the metabolic phenotype of Daphnia in relation to environmental hypoxia. Looking at the zooplankton community, we found that D. pulicaria was restricted to lakes with a hypoxic hypolimnion. Comparing D. pulicaria with different amounts of Hb, red animals showed higher mRNA levels for several Hb genes, whereas most glycolytic genes showed red/pale differences of less than 50%. We also observed seasonal changes in the metabolic phenotype that differed between red and pale animals. Hb was upregulated early in the season in hypoxic lakes, and a relationship between Hb and lactate dehydrogenase only emerged later in the season in a temporal pattern that was lake specific. To evaluate whether these differences were due to specific lake environments or microevolutionary differences, we tested the induction of genes under controlled hypoxia in isofemale lines from each of four lakes. We found a strong response to 18 h hypoxia exposure in both Hb and lactate dehydrogenase mRNA, although the magnitude of the acute response was greater than the steady state differences in mRNA levels between pale and red Daphnia. The baseline expression of Hb and lactate dehydrogenase also varied between isofemale lines with different lake origins. These results, in combination with comparison of glycogen measurements, suggests that Hb functions primarily to facilitate oxygen delivery, mitigating systemic hypoxia, rather than an oxygen store. The combination of lab and field studies suggest that the metabolic phenotype of the animal is influenced by both microevolutionary differences (within and between lakes) as well as the spatial and temporal environmental heterogeneity of the lakes. The differences between Daphnia species, and the unexpected lack of hypoxia sensitivity of select glycolytic genes provide evidence of macroevolutionary differences in metabolic strategies to cope with hypoxia.

Role of sex-concordant gene expression in the coevolution of exaggerated male and female genitalia in a beetle group

Some sexual traits, including genitalia, have undergone coevolutionary diversification toward exaggerated states in both sexes among closely related species, but the underlying genetic mechanisms that allow correlated character evolution between the sexes are poorly understood. Here, we studied interspecific differences in gene expression timing profiles involved in the correlated evolution of corresponding male and female genital parts in three species of ground beetle in Carabus (Ohomopterus). The male and female genital parts maintain morphological matching, whereas large interspecific variation in genital part size has occurred in the genital coevolution between the sexes toward exaggeration. We analyzed differences in gene expression involved in the interspecific differences in genital morphology using whole transcriptome data from genital tissues during genital morphogenesis. We found that the gene expression variance attributed to sex was negligible for the majority of differentially expressed genes, thus exhibiting sex-concordant expression, although large variances were attributed to stage and species differences. For each sex, we obtained co-expression gene networks and hub genes from differentially expressed genes between species that might be involved in interspecific differences in genital morphology. These gene networks were common to both sexes, and both sex-discordant and sex-concordant gene expression were likely involved in species-specific genital morphology. In particular, the gene expression related to exaggerated genital size showed no significant intersexual differences, implying that the genital sizes in both sexes are controlled by the same gene network with sex-concordant expression patterns, thereby facilitating the coevolution of exaggerated genitalia between the sexes while maintaining intersexual matching.

Next-generation sequencing of DNA from resting eggs: signatures of eutrophication in a lake's sediment

Hatching resting stages of ecologically important organisms such as Daphnia from lake sediments, referred to as resurrection ecology, is a powerful approach to assess changes in alleles and traits over time. However, the utility of the approach is constrained by a few obstacles, including low and/or biased hatching among genotypes. Here, we eliminated such bottlenecks by investigating DNA sequences isolated directly (i.e. without hatching) from resting eggs found in the sediments of Lake Constance spanning pre-, peri-, and post-eutrophication. While we expected genome-wide changes, we specifically expected changes in alleles related to pathways involved in mitigating effects of cyanobacterial toxins. We used pairwise F_ST-analyses to identify transcripts that showed strongest divergence among the four different populations and a clustering analysis to identify correlations between allele frequency shifts and changes in abiotic and biotic lake parameters. In a cluster that correlated with the increased abundance of cyanobacteria in Lake Constance we find genes that have been reported earlier to be differentially expressed in response to the cyanobacterial toxin microcystin and to microcystin-free cyanobacteria. We further reveal the enrichment of gene ontology terms that have been shown to be involved in microcystin-related responses in other organisms but not yet in Daphnia and as such are candidate loci for adaptation of natural Daphnia populations to increased cyanobacterial abundances. In conclusion this approach of investigating DNA extracted from Daphnia resting stages allowed to determine frequency changes of loci in a natural population over time.

Insights into the genetic basis of predator-induced response in Daphnia galeata

Phenotypic plastic responses allow organisms to rapidly adjust when facing environmental challenges-these responses comprise morphological, behavioral but also life-history changes. Alteration of life-history traits when exposed to predation risk have been reported often in the ecological and genomic model organism Daphnia. However, the molecular basis of this response is not well understood, especially in the context of fish predation. Here, we characterized the transcriptional profiles of two Daphnia galeata clonal lines with opposed life histories when exposed to fish kairomones. First, we conducted a differential gene expression, identifying a total of 125 candidate transcripts involved in the predator-induced response, uncovering substantial intraspecific variation. Second, we applied a gene coexpression network analysis to find clusters of tightly linked transcripts revealing the functional relations of transcripts underlying the predator-induced response. Our results showed that transcripts involved in remodeling of the cuticle, growth, and digestion correlated with the response to environmental change in D. galeata. Furthermore, we used an orthology-based approach to gain functional information for transcripts lacking gene ontology (GO) information, as well as insights into the evolutionary conservation of transcripts. We could show that our candidate transcripts have orthologs in other Daphnia species but almost none in other arthropods. The unique combination of methods allowed us to identify candidate transcripts, their putative functions, and evolutionary history associated with predator-induced responses in Daphnia. Our study opens up to the question as to whether the same molecular signature is associated with fish kairomones-mediated life-history changes in other Daphnia species.

Centennial clonal stability of asexual Daphnia in Greenland lakes despite climate variability

Climate and environmental condition drive biodiversity at many levels of biological organization, from populations to ecosystems. Combined with paleoecological reconstructions, palaeogenetic information on resident populations provides novel insights into evolutionary trajectories and genetic diversity driven by environmental variability. While temporal observations of changing genetic structure are often made of sexual populations, little is known about how environmental change affects the long-term fate of asexual lineages. Here, we provide information on obligately asexual, triploid Daphnia populations from three Arctic lakes in West Greenland through the past 200-300 years to test the impact of environmental change on the temporal and spatial population genetic structure. The contrasting ecological state of the lakes, specifically regarding salinity and habitat structure may explain the observed lake-specific clonal composition over time. Palaeolimnological reconstructions show considerable regional environmental fluctuations since 1,700 (the end of the Little Ice Age), but the population genetic structure in two lakes was almost unchanged with at most two clones per time period. Their local populations were strongly dominated by a single clone that has persisted for 250-300 years. We discuss possible explanations for the apparent population genetic stability: (a) persistent clones are general-purpose genotypes that thrive under broad environmental conditions, (b) clonal lineages evolved subtle genotypic differences unresolved by microsatellite markers, or (c) epigenetic modifications allow for clonal adaptation to changing environmental conditions. Our results motivate research into the mechanisms of adaptation in these populations, as well as their evolutionary fate in the light of accelerating climate change in the polar regions.