Brief Communication: The Predictable Network Topology of Evolutionary Genomic Constraint.

Large-scale comparative genomics studies offer valuable resources for understanding both functional and evolutionary rate constraints. It is suggested that constraint aligns with the topology of genomic networks, increasing toward the center, with intermediate nodes combining relaxed constraint with higher contributions to the phenotype due to pleiotropy. However, this pattern has yet to be demonstrated in vertebrates. This study shows that constraint intensifies toward the network's center in placental mammals. Genes with rate changes associated with emergence of hibernation cluster mostly toward intermediate positions, with higher constraint in faster-evolving genes, which is indicative of a "sweet spot" for adaptation. If this trend holds universally, network node metrics could predict high-constraint regions even in clades lacking empirical constraint data.

Functional genomics of abiotic environmental adaptation in lacertid lizards and other vertebrates.

Understanding the genomic basis of adaptation to different abiotic environments is important in the context of climate change and resulting short-term environmental fluctuations. Using functional and comparative genomics approaches, we here investigated whether signatures of genomic adaptation to a set of environmental parameters are concentrated in specific subsets of genes and functions in lacertid lizards and other vertebrates. We first identify 200 genes with signatures of positive diversifying selection from transcriptomes of 24 species of lacertid lizards and demonstrate their involvement in physiological and morphological adaptations to climate. To understand how functionally similar these genes are to previously predicted candidate functions for climate adaptation and to compare them with other vertebrate species, we then performed a meta-analysis of 1,100 genes under selection obtained from -omics studies in vertebrate species adapted to different abiotic factors. We found that the vertebrate gene set formed a tightly connected interactome, which was to 23% enriched in previously predicted functions of adaptation to climate, and to a large part (18%) involved in organismal stress response. We found a much higher degree of identical genes being repeatedly selected among different animal groups (43.6%), and of functional similarity and post-translational modifications than expected by chance, and no clear functional division between genes used for ectotherm and endotherm physiological strategies. In total, 171 out of 200 genes of Lacertidae were part of this network. These results highlight an important role of a comparatively small set of genes and their functions in environmental adaptation and narrow the set of candidate pathways and markers to be used in future research on adaptation and stress response related to climate change.

Thermal stress induces a positive phenotypic and molecular feedback loop in zebrafish embryos.

Aquatic organisms must cope with both rising and rapidly changing temperatures. These thermal changes can affect numerous traits, from molecular to ecological scales. Biotic stressors are already known to induce the release of chemical cues which trigger behavioural responses in other individuals. In this study, we infer whether fluctuating temperature, as an abiotic stressor, may similarly induce stress-like responses in individuals not directly exposed to the stressor. To test this hypothesis, zebrafish (Danio rerio) embryos were exposed for 24 h to fluctuating thermal stress, to medium in which another embryo was thermally stressed before ("stress medium"), and to a combination of these. Growth, behaviour, expression of molecular markers, and of whole-embryo cortisol were used to characterise the thermal stress response and its propagation between embryos. Both fluctuating high temperature and stress medium significantly accelerated development, by shifting stressed embryos from segmentation to pharyngula stages, and altered embryonic activity. Importantly, we found that the expression of sulfide:quinone oxidoreductase (SQOR), the antioxidant gene SOD1, and of interleukin-1ß (IL-1ß) were significantly altered by stress medium. This study illustrates the existence of positive thermal stress feedback loops in zebrafish embryos where heat stress can induce stress-like responses in conspecifics, but which might operate via different molecular pathways. If similar effects also occur under less severe heat stress regimes, this mechanism may be relevant in natural settings as well.

Aligning functional network constraint to evolutionary outcomes.

BACKGROUND: Functional constraint through genomic architecture is suggested to be an important dimension of genome evolution, but quantitative evidence for this idea is rare. In this contribution, existing evidence and discussions on genomic architecture as constraint for convergent evolution, rapid adaptation, and genic adaptation are summarized into alternative, testable hypotheses. Network architecture statistics from protein-protein interaction networks are then used to calculate differences in evolutionary outcomes on the example of genomic evolution in yeast, and the results are used to evaluate statistical support for these longstanding hypotheses. RESULTS: A discriminant function analysis lent statistical support to classifying the yeast interactome into hub, intermediate and peripheral nodes based on network neighborhood connectivity, betweenness centrality, and average shortest path length. Quantitative support for the existence of genomic architecture as a mechanistic basis for evolutionary constraint is then revealed through utilizing these statistical parameters of the protein-protein interaction network in combination with estimators of protein evolution. CONCLUSIONS: As functional genetic networks are becoming increasingly available, it will now be possible to evaluate functional genetic network constraint against variables describing complex phenotypes and environments, for better understanding of commonly observed deterministic patterns of evolution in non-model organisms. The hypothesis framework and methodological approach outlined herein may help to quantify the extrinsic versus intrinsic dimensions of evolutionary constraint, and result in a better understanding of how fast, effectively, or deterministically organisms adapt.

Comparative analysis of amphibian genomes: An emerging resource for basic and applied research.

Amphibians are the most threatened group of vertebrates and are in dire need of conservation intervention to ensure their continued survival. They exhibit unique features including a high diversity of reproductive strategies, permeable and specialized skin capable of producing toxins and antimicrobial compounds, multiple genetic mechanisms of sex determination and in some lineages, the ability to regenerate limbs and organs. Although genomic approaches would shed light on these unique traits and aid conservation, sequencing and assembly of amphibian genomes has lagged behind other taxa due to their comparatively large genome sizes. Fortunately, the development of long-read sequencing technologies and initiatives has led to a recent burst of new amphibian genome assemblies. Although growing, the field of amphibian genomics suffers from the lack of annotation resources, tools for working with challenging genomes and lack of high-quality assemblies in multiple clades of amphibians. Here, we analyse 51 publicly available amphibian genomes to evaluate their usefulness for functional genomics research. We report considerable variation in genome assembly quality and completeness and report some of the highest transposable element and repeat contents of any vertebrate. Additionally, we detected an association between transposable element content and climatic variables. Our analysis provides evidence of conserved genome synteny despite the long divergence times of this group, but we also highlight inconsistencies in chromosome naming and orientation across genome assemblies. We discuss sequencing gaps in the phylogeny and suggest key targets for future sequencing endeavours. Finally, we propose increased investment in amphibian genomics research to promote their conservation.

Sex and gametogenesis stage are strong drivers of gene expression in Mytilus edulis exposed to environmentally relevant plasticiser levels and pH 7.7.

Plastic pollution and changes in oceanic pH are both pressing environmental issues. Little emphasis, however, has been placed on the influence of sex and gametogenesis stage when investigating the effects of such stressors. Here, we examined histology and molecular biomarkers of blue mussels Mytilus edulis exposed for 7 days to a pH 7.7 scenario (- 0.4 units) in combination with environmentally relevant concentrations (0, 0.5 and 50 µg/L) of the endocrine disrupting plasticiser di-2-ethylhexyl phthalate (DEHP). Through a factorial design, we investigated the gametogenesis cycle and sex-related expression of genes involved in pH homeostasis, stress response and oestrogen receptor-like pathways after the exposure to the two environmental stressors. As expected, we found sex-related differences in the proportion of developing, mature and spawning gonads in histological sections. Male gonads also showed higher levels of the acid-base regulator CA2, but females had a higher expression of stress response-related genes (i.e. sod, cat, hsp70). We found a significant effect of DEHP on stress response-related gene expression that was dependent on the gametogenesis stage, but there was only a trend towards downregulation of CA2 in response to pH 7.7. In addition, differences in gene expression between males and females were most pronounced in experimental conditions containing DEHP and/or acidified pH but never the control, indicating that it is important to consider sex and gametogenesis stage when studying the response of mussels to diverse stressors.

Heat induces multiomic and phenotypic stress propagation in zebrafish embryos.

Heat alters biology from molecular to ecological levels, but may also have unknown indirect effects. This includes the concept that animals exposed to abiotic stress can induce stress in naive receivers. Here, we provide a comprehensive picture of the molecular signatures of this process, by integrating multiomic and phenotypic data. In individual zebrafish embryos, repeated heat peaks elicited both a molecular response and a burst of accelerated growth followed by a growth slowdown in concert with reduced responses to novel stimuli. Metabolomes of the media of heat treated vs. untreated embryos revealed candidate stress metabolites including sulfur-containing compounds and lipids. These stress metabolites elicited transcriptomic changes in naive receivers related to immune response, extracellular signaling, glycosaminoglycan/keratan sulfate, and lipid metabolism. Consequently, non-heat-exposed receivers (exposed to stress metabolites only) experienced accelerated catch-up growth in concert with reduced swimming performance. The combination of heat and stress metabolites accelerated development the most, mediated by apelin signaling. Our results prove the concept of indirect heat-induced stress propagation toward naive receivers, inducing phenotypes comparable with those resulting from direct heat exposure, but utilizing distinct molecular pathways. Group-exposing a nonlaboratory zebrafish line, we independently confirm that the glycosaminoglycan biosynthesis-related gene chs1 and the mucus glycoprotein gene prg4a, functionally connected to the candidate stress metabolite classes sugars and phosphocholine, are differentially expressed in receivers. This hints at the production of Schreckstoff-like cues in receivers, leading to further stress propagation within groups, which may have ecological and animal welfare implications for aquatic populations in a changing climate.

Investigating the effects of arginine methylation inhibitors on microdissected brain tumour biopsies maintained in a miniaturised perfusion system.

Arginine methylation is a post-translational modification that consists of the transfer of one or two methyl (CH3) groups to arginine residues in proteins. Several types of arginine methylation occur, namely monomethylation, symmetric dimethylation and asymmetric dimethylation, which are catalysed by different protein arginine methyltransferases (PRMTs). Inhibitors of PRMTs have recently entered clinical trials to target several types of cancer, including gliomas (NCT04089449). People with glioblastoma (GBM), the most aggressive form of brain tumour, are among those with the poorest quality of life and likelihood of survival of anyone diagnosed with cancer. There is currently a lack of (pre)clinical research on the possible application of PRMT inhibitors to target brain tumours. Here, we set out to investigate the effects of clinically-relevant PRMT inhibitors on GBM biopsies. We present a new, low-cost, easy to fabricate perfusion device that can maintain GBM tissue in a viable condition for at least eight days post-surgical resection. The miniaturised perfusion device enables the treatment of GBM tissue with PRMT inhibitors ex vivo, and we observed a two-fold increase in apoptosis in treated samples compared to parallel control experiments. Mechanistically, we show thousands of differentially expressed genes after treatment, and changes in the type of arginine methylation of the RNA binding protein FUS that are consistent with hundreds of differential gene splicing events. This is the first time that cross-talk between different types of arginine methylation has been observed in clinical samples after treatment with PRMT inhibitors.

Consequences of combined exposure to thermal stress and the plasticiser DEHP in Mytilus spp. differ by sex.

Little is known about the combined effect of environmental factors and contaminants on commercially important marine species, and whether this effect differs by sex. In this study, blue mussels were exposed for seven days to both single and combined stressors (i.e., +3 °C elevated temperature and two environmentally relevant concentrations of the plastic softener DEHP, 0.5 and 50 µg/l) in a factorial design. Males were observed to be more sensitive to high temperature, demonstrated by the significant increase in out-of-season spawning gonads and higher gene expression of the antioxidant catalase and the estrogen receptor genes. On the other hand, while the gametogenesis cycle in females was more resilient than in males, DEHP exposure altered the estrogen-related receptor gene expression. We show that the combined stressors DEHP and increased temperature, in environmentally relevant magnitudes, have different consequences in male and female mussels, with the potential to impact the timing and breeding season success in Mytilus spp.

Transcriptomic Signatures of Experimental Alkaloid Consumption in a Poison Frog.

In the anuran family Dendrobatidae, aposematic species obtain their toxic or unpalatable alkaloids from dietary sources, a process known as sequestering. To understand how toxicity evolved in this family, it is paramount to elucidate the pathways of alkaloid processing (absorption, metabolism, and sequestering). Here, we used an exploratory skin gene expression experiment in which captive-bred dendrobatids were fed alkaloids. Most of these experiments were performed with Dendrobates tinctorius, but some trials were performed with D. auratus, D. leucomelas and Allobates femoralis to explore whether other dendrobatids would show similar patterns of gene expression. We found a consistent pattern of up-regulation of genes related to muscle and mitochondrial processes, probably due to the lack of mutations related to alkaloid resistance in these species. Considering conserved pathways of drug metabolism in vertebrates, we hypothesize alkaloid degradation is a physiological mechanism of resistance, which was evidenced by a strong upregulation of the immune system in D. tinctorius, and of complement C2 across the four species sampled. Probably related to this strong immune response, we found several skin keratins downregulated, which might be linked to a reduction of the cornified layer of the epidermis. Although not conclusive, our results offer candidate genes and testable hypotheses to elucidate alkaloid processing in poison frogs.