The power and limitations of gene expression pathway analyses toward predicting population response to environmental stressors

FLEXIBILITY IN GENE COEXPRESSION AT DEVELOPMENTAL AND EVOLUTIONARY TIMESCALES

The explosion of next-generation sequencing technologies has allowed researchers to move from studying single genes, to thousands of genes, and thereby to also consider the relationships within gene networks. Like others, we are interested in understanding how developmental and evolutionary forces shape the expression of individual genes, as well as the interactions among genes. To this end, we characterized the effects of genetic background and developmental environment on brain gene coexpression in two parallel, independent evolutionary lineages of Trinidadian guppies (Poecilia reticulata). We asked whether connectivity patterns among genes differed based on genetic background and rearing environment, and whether a gene's connectivity predicted its propensity for expression divergence. In pursuing these questions, we confronted the central challenge that standard approaches fail to control the Type I error and/or have low power in the presence of high dimensionality (i.e., large number of genes) and small sample size, as in many gene expression studies. Using our data as a case study, we detail central challenges, discuss sample size guidelines, and provide rigorous statistical approaches for exploring coexpression differences with small sample sizes. Using these approaches, we find evidence that coexpression relationships differ based on both genetic background and rearing environment. We report greater expression divergence in less connected genes and suggest this pattern may arise and be reinforced by selection.

Genomic insights into variation in thermotolerance between hybridizing swordtail fishes

Understanding how organisms adapt to changing environments is a core focus of research in evolutionary biology. One common mechanism is adaptive introgression, which has received increasing attention as a potential route to rapid adaptation in populations struggling in the face of ecological change, particularly global climate change. However, hybridization can also result in deleterious genetic interactions that may limit the benefits of adaptive introgression. Here, we used a combination of genome-wide quantitative trait locus mapping and differential gene expression analyses between the swordtail fish species Xiphophorus malinche and X. birchmanni to study the consequences of hybridization on thermotolerance. While these two species are adapted to different thermal environments, we document a complicated architecture of thermotolerance in hybrids. We identify a region of the genome that contributes to reduced thermotolerance in individuals heterozygous for X. malinche and X. birchmanni ancestry, as well as widespread misexpression in hybrids of genes that respond to thermal stress in the parental species, particularly in the circadian clock pathway. We also show that a previously mapped hybrid incompatibility between X. malinche and X. birchmanni contributes to reduced thermotolerance in hybrids. Together, our results highlight the challenges of understanding the impact of hybridization on complex ecological traits and its potential impact on adaptive introgression.

The master regulator OxyR orchestrates bacterial oxidative stress response genes in space and time

Bacteria employ diverse gene regulatory networks to survive stress, but deciphering the underlying logic of these complex networks has proved challenging. Here, we use time-resolved single-cell imaging to explore the functioning of the E. coli regulatory response to oxidative stress. We observe diverse gene expression dynamics within the network. However, by controlling for stress-induced growth-rate changes, we show that these patterns involve just three classes of regulation: downregulated genes, upregulated pulsatile genes, and gradually upregulated genes. The two upregulated classes are distinguished by differences in the binding of the transcription factor, OxyR, and appear to play distinct roles during stress protection. Pulsatile genes activate transiently in a few cells for initial protection of a group of cells, whereas gradually upregulated genes induce evenly, generating a lasting protection involving many cells. Our study shows how bacterial populations use simple regulatory principles to coordinate stress responses in space and time. A record of this paper's transparent peer review process is included in the supplemental information.

E.PathDash, pathway activation analysis of publicly available pathogen gene expression data

E.PathDash facilitates re-analysis of gene expression data from pathogens clinically relevant to chronic respiratory diseases, including a total of 48 studies, 548 samples, and 404 unique treatment comparisons. The application enables users to assess broad biological stress responses at the KEGG pathway or gene ontology level and also provides data for individual genes. E.PathDash reduces the time required to gain access to data from multiple hours per data set to seconds. Users can download high-quality images such as volcano plots and boxplots, differential gene expression results, and raw count data, making it fully interoperable with other tools. Importantly, users can rapidly toggle between experimental comparisons and different studies of the same phenomenon, enabling them to judge the extent to which observed responses are reproducible. As a proof of principle, we invited two cystic fibrosis scientists to use the application to explore scientific questions relevant to their specific research areas. Reassuringly, pathway activation analysis recapitulated results reported in original publications, but it also yielded new insights into pathogen responses to changes in their environments, validating the utility of the application. All software and data are freely accessible, and the application is available at scangeo.dartmouth.edu/EPathDash.

IMPORTANCE

Chronic respiratory illnesses impose a high disease burden on our communities and people with respiratory diseases are susceptible to robust bacterial infections from pathogens, including Pseudomonas aeruginosa and Staphylococcus aureus, that contribute to morbidity and mortality. Public gene expression datasets generated from these and other pathogens are abundantly available and an important resource for synthesizing existing pathogenic research, leading to interventions that improve patient outcomes. However, it can take many hours or weeks to render publicly available datasets usable; significant time and skills are needed to clean, standardize, and apply reproducible and robust bioinformatic pipelines to the data. Through collaboration with two microbiologists, we have shown that E.PathDash addresses this problem, enabling them to elucidate pathogen responses to a variety of over 400 experimental conditions and generate mechanistic hypotheses for cell-level behavior in response to disease-relevant exposures, all in a fraction of the time.

Metabolic differentiation of brushtail possum populations resistant and susceptible to plant toxins revealed via differential gene expression

The Australian brushtail possum (Trichosurus vulpecula) is adapted to a wide range of food plants across its range and is exposed to numerous physiological challenges. Populations that are resistant to the plant toxin sodium fluoroacetate are of particular interest as this compound has been used since the 1940s for vertebrate pest management around the world. Candidate gene identification is an important first step in understanding how spatial populations have responded to local selection resulting in local physiological divergence. We employ differential gene expression of liver samples from wild-caught brushtail possums from toxin-resistant and toxin-susceptible populations to identify candidate genes that might be involved in metabolic pathways associated with toxin-resistance. This allowed us to identify genetic pathways involved in resistance to the plant toxin sodium fluoroacetate in Western Australian possums but not those originally from south eastern Australia. We identified differentially expressed genes in the liver that are associated with cell signalling, encapsulating structure, cell mobility, and tricarboxylic acid cycle. The gene expression differences detected indicate which metabolic pathways are most likely to be associated with sodium fluoroacetate resistance in these marsupials and we provide a comprehensive list of candidate genes and pathways to focus on for future studies.

Comparative Transcriptomic Reveals Greater Similarities in Response to Temperature Than to Invasive Alien Predator in the Damselfly Ischnura elegans Across Different Geographic Scales

The impact of global changes on populations may not be necessarily uniform across a species' range. Here, we aim at comparing the phenotypic and transcriptomic response to warming and an invasive predator cue in populations across different geographic scales in the damselfly Ischnura elegans. We collected adult females in two ponds in southern Poland (central latitude) and two ponds in southern Sweden (high latitude). We raised their larvae in growth chambers and exposed them to combination of temperature and a predator cue released by the crayfish Orconectes limosus. When larvae reached the prefinal larval stage, they were phenotyped for traits related to growth and size and collected for a gene expression analysis. High-latitude populations exhibited greater phenotypic and transcriptomic variation than central-latitude populations. Across latitudes and ponds, temperature generally increased growth rate and the predator cue decreased mass, but the effects of temperature were also pond-specific. Comparison of the transcriptomic profiles revealed a greater overlap in the response to temperature across latitudes and ponds, especially for pathway-related oxidative stress and sugar and lipid metabolism. The transcriptomic response to a predator cue and to the interaction temperature × predator cue was more pond-specific and overlapped only for few genes and pathways related to cuticle, development and signal transduction. We demonstrated that central- and high-latitude populations may partially respond through similar mechanisms to warming and, to a lower extent to a predator cue and to the interaction temperature × predator cue. For the predator cue and the interaction, the large fraction of ponds-specific genes suggests local adaptation. We show that high-latitude populations were generally more plastic at the phenotypic and transcriptomic level and may be more capable to cope with environmental changes than their central-latitude counterparts.

MICROPHERRET: MICRObial PHEnotypic tRait ClassifieR using Machine lEarning Techniques

BACKGROUND

In recent years, there has been a rapid increase in the number of microbial genomes reconstructed through shotgun sequencing, and obtained by newly developed approaches including metagenomic binning and single-cell sequencing. However, our ability to functionally characterize these genomes by experimental assays is orders of magnitude less efficient. Consequently, there is a pressing need for the development of swift and automated strategies for the functional classification of microbial genomes.

RESULTS

The present work leverages a suite of supervised machine learning algorithms to establish a range of 86 metabolic and other ecological functions, such as methanotrophy and plastic degradation, starting from widely obtainable microbial genome annotations. Tests performed on independent datasets demonstrated robust performance across complete, fragmented, and incomplete genomes above a 70% completeness level for most of the considered functions. Application of the algorithms to the Biogas Microbiome database yielded predictions broadly consistent with current biological knowledge and correctly detecting functionally-related nuances of archaeal genomes. Finally, a case study focused on acetoclastic methanogenesis demonstrated how the developed machine learning models can be refined or expanded with models describing novel functions of interest.

CONCLUSIONS

The resulting tool, MICROPHERRET, incorporates a total of 86 models, one for each tested functional class, and can be applied to high-quality microbial genomes as well as to low-quality genomes derived from metagenomics and single-cell sequencing. MICROPHERRET can thus aid in understanding the functional role of newly generated genomes within their micro-ecological context.

From whole bodies to single cells: A guide to transcriptomic approaches for ecology and evolutionary biology

RNA sequencing (RNAseq) methodology has experienced a burst of technological developments in the last decade, which has opened up opportunities for studying the mechanisms of adaptation to environmental factors at both the organismal and cellular level. Selecting the most suitable experimental approach for specific research questions and model systems can, however, be a challenge and researchers in ecology and evolution are commonly faced with the choice of whether to study gene expression variation in whole bodies, specific tissues, and/or single cells. A wide range of sometimes polarised opinions exists over which approach is best. Here, we highlight the advantages and disadvantages of each of these approaches to provide a guide to help researchers make informed decisions and maximise the power of their study. Using illustrative examples of various ecological and evolutionary research questions, we guide the readers through the different RNAseq approaches and help them identify the most suitable design for their own projects.

Identification and validation of a novel risk model based on cuproptosis‑associated m6A for head and neck squamous cell carcinoma

BACKGROUND

Head and neck squamous cell carcinoma (HNSCC) is a prevalent cancer with a poor survival rate due to anatomical limitations of the head and a lack of reliable biomarkers. Cuproptosis represents a novel cellular regulated death pathway, and N6-methyladenosine (m6A) is the most common internal RNA modification in mRNA. They are intricately connected to tumor formation, progression, and prognosis. This study aimed to construct a risk model for HNSCC using a set of mRNAs associated with m6A regulators and cuproptosis genes (mcrmRNA).

METHODS

RNA-seq and clinical data of HNSCC patients from The Cancer Genome Atlas (TCGA) database were analyzed to develop a risk model through the least absolute shrinkage and selection operator (LASSO) analysis. Survival analysis and receiver operating characteristic (ROC) analysis were performed for the high- and low-risk groups. Additionally, the model was validated using the GSE41613 dataset from the Gene Expression Omnibus (GEO) database. GSEA and CIBERSORT were applied to investigate the immune microenvironment of HNSCC.

RESULTS

A risk model consisting of 32 mcrmRNA was developed using the LASSO analysis. The risk score of patients was confirmed to be an independent prognostic indicator by multivariate Cox analysis. The high-risk group exhibited a higher tumor mutation burden. Additionally, CIBERSORT analysis indicated varying levels of immune cell infiltration between the two groups. Significant disparities in drug sensitivity to common medications were also observed. Enrichment analysis further unveiled significant differences in metabolic pathways and RNA processing between the two groups.

CONCLUSIONS

Our risk model can predict outcomes for HNSCC patients and offers valuable insights for personalized therapeutic approaches.

Investigation of the Alternations in the Muscle Quality of Swimming Crab (Ovalipes punctatus) during Cold-Chain Transportation Using Physicochemical and TMT-Based Quantitative Proteomic Analysis

Physicochemical properties and protein alterations in Ovalipes punctatus during cold-chain transportation were examined via sensory scores, water-holding capacity (WHC), glucose (GLU) content, catalase (CAT) activity, urea nitrogen (UN) content, and tandem mass tag (TMT)-based proteomic analysis. The results revealed that sensory characteristics and texture of crab muscle deteriorated during transportation. Proteomic analysis revealed 442 and 470 different expressed proteins (DEPs) in crabs after 18 h (FC) and 36 h (DC) of transportation compared with live crabs (LC). Proteins related to muscle structure and amino acid metabolism significantly changed, as evidenced by the decreased WHC and sensory scores of crab muscle. Glycolysis, calcium signaling, and peroxisome pathways were upregulated in the FC/LC comparison, aligning with the changes in GLU content and CAT activity, revealing the stress response of energy metabolism and immune response in crabs during 0-18 h of transportation. The downregulated tricarboxylic acid (TCA) cycle and carcinogenesis-reactive oxygen species pathways were correlated with the decreasing trend in CAT activity, suggesting a gradual retardation in both energy and antioxidant metabolism in crabs during 18-36 h of transportation. Furthermore, the regulated purine nucleoside metabolic and nucleoside diphosphate-related processes, with the increasing changes in UN content, revealed the accumulation of metabolites in crabs.

PRANA: an R package for differential co-expression network analysis with the presence of additional covariates

BACKGROUND

Advances in sequencing technology and cost reduction have enabled an emergence of various statistical methods used in RNA-sequencing data, including the differential co-expression network analysis (or differential network analysis). A key benefit of this method is that it takes into consideration the interactions between or among genes and do not require an established knowledge in biological pathways. As of now, none of existing softwares can incorporate covariates that should be adjusted if they are confounding factors while performing the differential network analysis.

RESULTS

We develop an R package PRANA which a user can easily include multiple covariates. The main R function in this package leverages a novel pseudo-value regression approach for a differential network analysis in RNA-sequencing data. This software is also enclosed with complementary R functions for extracting adjusted p-values and coefficient estimates of all or specific variable for each gene, as well as for identifying the names of genes that are differentially connected (DC, hereafter) between subjects under biologically different conditions from the output.

CONCLUSION

Herewith, we demonstrate the application of this package in a real data on chronic obstructive pulmonary disease. PRANA is available through the CRAN repositories under the GPL-3 license: https://cran.r-project.org/web/packages/PRANA/index.html .

Identification of biomarkers and sex differences in the placenta of fetal growth restriction

AIM

Fetal growth restriction (FGR) can lead to short-term and long-term impairments in the fetus. The placenta functions as an exchanger for substance transport, playing a critical role in fetal growth. However, the mechanism from the placental standpoint is still not fully understood. In this study, we aimed to investigate the pathophysiological mechanisms in the placenta that mediated the development of FGR and sex differences.

METHODS

We analyzed the gene expression profiles of GSE100415 containing specific normotensive FGR placental samples and GSE114691 with canonical samples using three different methods, differentially expressed gene analysis, weighted gene co-expression network analysis, and gene set enrichment analysis. Gene enrichment was performed, including the gene ontology and pathway from the Kyoto Encyclopedia of Genes and Genomes. The important process was then validated in pregnant Wistar rats subcutaneously administered dexamethasone (0.2 mg/kg/d) or saline from gestation Day 9 to 21.

RESULTS

Our results revealed little difference between the comparison of normal and normotensive FGR placental samples but confirmed the sex difference. Further analyses of the canonical samples identified the occurrence of vascular dysfunction, which was validated by the calculation of the vascular lumen area, showing that the vascular lumen in the FGR group was more than in the control. We also discovered 17 significantly expressed genes from the involved eigengenes.

CONCLUSION

Our study provides an important theoretical and experimental basis to reevaluate the development of FGR from the placental standpoint and suggests a series of biomarkers for future clinical use.

Ex vivo exposures to arsenite and its methylated trivalent metabolites alter gene transcription in mouse sperm cells

We have previously reported that preconception exposure to iAs may contribute to the development of diabetes in mouse offspring by altering gene expressions in paternal sperm. However, the individual contributions of iAs and its methylated metabolites, monomethylated arsenic (MAs) and dimethylated arsenic (DMAs), to changes in the sperm transcriptome could not be determined because all three As species are present in sperm after in vivo iAs exposure. The goal of the present study was to assess As species-specific effects using an ex vivo model. We exposed freshly isolated mouse sperm to either 0.1 or 1 μM arsenite (iAsIII) or the methylated trivalent arsenicals, MAsIII and DMAsIII, and used RNA-sequencing to identify differentially expressed genes, enriched pathways, and associated protein networks. For all arsenicals tested, the exposures to 0.1 μM concentrations had greater effects on gene expression than 1 μM exposures. Transcription factor AP-1 and B cell receptor complexes were the most significantly enriched pathways in sperm exposed to 0.1 μM iAsIII. The Mre11 complex and Antigen processing were top pathways targeted by exposure to 0.1 μM MAsIII and DMAsIII, respectively. While there was no overlap between gene transcripts altered by ex vivo exposures in the present study and those altered by in vivo exposure in our prior work, several pathways were shared, including PI3K-Akt signaling, Focal adhesion, and Extracellular matrix receptor interaction pathways. Notably, the protein networks associated with these pathways included those with known roles in diabetes. This study is the first to assess the As species-specific effects on sperm transcriptome, linking these effects to the diabetogenic effects of iAs exposure.

Gene expression studies of plastic and evolutionary responses to global warming

Phenotypic plasticity can be a rapid response for coping with global warming, yet may be insufficient to protect species from extinction. Evolutionary adaptation may reinforce adaptive or oppose maladaptive plastic responses. With advances in technology whole transcriptomes can provide us with an unprecedented overview of genes and functional processes underlying the interplay between plasticity and evolution. We advocate that insects provide ideal opportunities to study plasticity in non-adapted and thermally adapted populations to infer reaction norms across the whole transcriptome ('reactionomes'). This can advance our understanding of how the interplay between plasticity and evolution shapes responses to warming. So far, a limited number of studies suggest predominantly maladaptive plastic responses to novel environments that are reduced with time, but much more research is needed to infer general patterns.

Cross-feeding interactions in short chain gaseous alkane-driven perchlorate and selenate reduction

Short chain gaseous alkanes (SCGAs) mainly consist of methane (CH₄), ethane (C₂H₆), propane (C₃H₈) and butane (C₄H₁₀). The first three SCGAs have been shown to remove perchlorate (ClO₄^-) and selenate (SeO₄^2-), yet it is unknown whether C₄H₁₀ is available to reduce these contaminants. This study demonstrated that C₄H₁₀ fed biofilms were capable of reducing ClO₄^- and SeO₄^2- to chloride (Cl^-) and elemental selenium (Se⁰), respectively, by employing two independent membrane biofilms reactors (MBfRs). Batch tests showed that C₄H₁₀ and oxygen fed biofilms had much higher ClO₄^- and SeO₄^2- reduction rates and enhanced expression levels of bmoX and pcrA than that without C₄H₁₀ or O₂. Polyhydroxyalkanoates (PHA) accumulated in the biofilms when C₄H₁₀ was supplied, and they decomposed for driving ClO₄^- and SeO₄^2- reduction when C₄H₁₀ was absent. Moreover, we revisited the literature and found that a cross-feeding pathway seems to be universal in microaerobic SCGA-driven perchlorate and selenate reduction processes. In the ClO₄^--reducing MBfRs, Mycobacterium primarily conducts C₂H₆ and C₃H₈ oxidation in synergy with Dechloromonas who performs perchlorate reduction, while both Mycobacterium and Rhodococcus carried out C₄H₁₀ oxidation with perchlorate-respiring Azospira as the partner. In the SeO₄^2--reducing MBfRs, Mycobacterium oxidized C₂H₆ solely or oxidized C₃H₈ jointly with Rhodococcus, while Burkholderiaceae likely acted as the selenate-reducing bacterium. When C₄H₁₀ was supplied as the electron donor, both Mycobacterium and Rhodococcus conducted C₄H₁₀ oxidation in synergy with unknow selenate-reducing bacterium. Collectively, we confirm that from CH₄ to C₄H₁₀, all SCGAs could be utilized as electron donors for bio-reduction process. These findings offer insights into SCGA-driven bio-reduction processes, and are helpful in establishing SCGA-based technologies for groundwater remediation.

Mitochondrial volume density and evidence for its role in adaptive divergence in response to thermal tolerance in threespine stickleback

Phenotypic plasticity is predicted to permit persistence in new environments, and may subsequently evolve to enhance fitness. Colonizing environments with lower winter temperatures can lead to the evolution of lower critical thermal minima; the corresponding physiological traits associated with temperature tolerance are predicted to involve mitochondrial function. Threespine stickleback (Gasterosteus aculeatus) have colonized freshwater lakes along the Pacific Northwest. These freshwater populations are known to exhibit cold-induced increases in mitochondrial volume density in pectoral muscle, but whether such plasticity evolved before or after colonization is uncertain. Here, we measure critical thermal minima (CTmin) in one marine and one freshwater population of threespine stickleback, and mitochondrial volume density in pectoral and cardiac tissue of both populations acclimated to different temperature treatments (6.2, 14.5 and 20.6 ℃). Mitochondrial volume density increased with cold acclimation in pectoral muscle; cardiac muscle was non-plastic but had elevated mitochondrial volume densities compared to pectoral muscle across all temperature treatments. There were no differences in the levels of plasticity between marine and freshwater stickleback, but neither were there differences in CTmin. Importantly, marine stickleback exhibited plasticity under low-salinity conditions, suggesting that marine stickleback had at least one necessary phenotype for persistence in freshwater environments before colonization occurred.

Genetic compensation rather than genetic assimilation drives the evolution of plasticity in response to mild warming across latitudes in a damselfly

Global warming is causing plastic and evolutionary changes in the phenotypes of ectotherms. Yet, we have limited knowledge on how the interplay between plasticity and evolution shapes thermal responses and underlying gene expression patterns. We assessed thermal reaction norm patterns across the transcriptome and identified associated molecular pathways in northern and southern populations of the damselfly Ischnura elegans. Larvae were reared in a common garden experiment at the mean summer water temperatures experienced at the northern (20°C) and southern (24°C) latitudes. This allowed a space-for-time substitution where the current gene expression levels at 24°C in southern larvae are a proxy for the expected responses of northern larvae under gradual thermal evolution to the predicted 4°C warming. Most differentially expressed genes showed fixed differences across temperatures between latitudes, suggesting that thermal genetic adaptation will mainly evolve through changes in constitutive gene expression. Northern populations also frequently showed plastic responses in gene expression to mild warming, while southern populations were much less responsive to temperature. Thermal responsive genes in northern populations showed to a large extent a pattern of genetic compensation, namely gene expression that was induced at 24°C in northern populations remained at a lower constant level in southern populations, and were associated with metabolic and translation pathways. There was instead little evidence for genetic assimilation of an initial plastic response to mild warming. Our data therefore suggest that genetic compensation rather than genetic assimilation may drive the evolution of plasticity in response to mild warming in this damselfly species.