Enriching the gene set analysis of genome-wide data by incorporating directionality of gene expression and combining statistical hypotheses and methods

Decoding the epitranscriptional landscape from native RNA sequences

Traditional epitranscriptomics relies on capturing a single RNA modification by antibody or chemical treatment, combined with short-read sequencing to identify its transcriptomic location. This approach is labor-intensive and may introduce experimental artifacts. Direct sequencing of native RNA using Oxford Nanopore Technologies (ONT) can allow for directly detecting the RNA base modifications, although these modifications might appear as sequencing errors. The percent Error of Specific Bases (%ESB) was higher for native RNA than unmodified RNA, which enabled the detection of ribonucleotide modification sites. Based on the %ESB differences, we developed a bioinformatic tool, epitranscriptional landscape inferring from glitches of ONT signals (ELIGOS), that is based on various types of synthetic modified RNA and applied to rRNA and mRNA. ELIGOS is able to accurately predict known classes of RNA methylation sites (AUC > 0.93) in rRNAs from Escherichiacoli, yeast, and human cells, using either unmodified in vitro transcription RNA or a background error model, which mimics the systematic error of direct RNA sequencing as the reference. The well-known DRACH/RRACH motif was localized and identified, consistent with previous studies, using differential analysis of ELIGOS to study the impact of RNA m6A methyltransferase by comparing wild type and knockouts in yeast and mouse cells. Lastly, the DRACH motif could also be identified in the mRNA of three human cell lines. The mRNA modification identified by ELIGOS is at the level of individual base resolution. In summary, we have developed a bioinformatic software package to uncover native RNA modifications.

Integrative analyses of TEDDY Omics data reveal lipid metabolism abnormalities, increased intracellular ROS and heightened inflammation prior to autoimmunity for type 1 diabetes

BACKGROUND

The Environmental Determinants of Diabetes in the Young (TEDDY) is a prospective birth cohort designed to study type 1 diabetes (T1D) by following children with high genetic risk. An integrative multi-omics approach was used to evaluate islet autoimmunity etiology, identify disease biomarkers, and understand progression over time.

RESULTS

We identify a multi-omics signature that was predictive of islet autoimmunity (IA) as early as 1 year before seroconversion. At this time, abnormalities in lipid metabolism, decreased capacity for nutrient absorption, and intracellular ROS accumulation are detected in children progressing towards IA. Additionally, extracellular matrix remodeling, inflammation, cytotoxicity, angiogenesis, and increased activity of antigen-presenting cells are observed, which may contribute to beta cell destruction. Our results indicate that altered molecular homeostasis is present in IA-developing children months before the actual detection of islet autoantibodies, which opens an interesting window of opportunity for therapeutic intervention.

CONCLUSIONS

The approach employed herein for assessment of the TEDDY cohort showcases the utilization of multi-omics data for the modeling of complex, multifactorial diseases, like T1D.

Transcriptomic analyses of livers from mice exposed to 1,4-dioxane for up to 90 days to assess potential mode(s) of action underlying liver tumor development

1,4-Dioxane is a volatile organic compound with industrial and commercial applications as a solvent and in the manufacture of other chemicals. 1,4-Dioxane has been demonstrated to induce liver tumors in chronic rodent bioassays conducted at very high doses. The available evidence for 1,4-dioxane-induced liver tumors in rodents aligns with a threshold-dependent mode of action (MOA), with the underlying mechanism being less clear in the mouse than in rats. To gain a better understanding of the underlying molecular mechanisms related to liver tumor development in mice orally exposed to 1,4-dioxane, transcriptomics analysis was conducted on liver tissue collected from a 90-day drinking water study in female B6D2F₁/Crl mice (Lafranconi et al., 2020). Using tissue samples from female mice exposed to 1,4-dioxane in the drinking water at concentrations of 0, 40, 200, 600, 2,000 or 6,000 ppm for 7, 28, and 90 days, transcriptomic analyses demonstrate minimal treatment effects on global gene expression at concentrations below 600 ppm. At higher concentrations, genes involved in phase II metabolism and mitotic cell cycle checkpoints were significantly upregulated. There was an overall lack of enrichment of genes related to DNA damage response. The increase in mitotic signaling is most prevalent in the livers of mice exposed to 1,4-dioxane at the highest concentrations for 90 days. This finding aligns with phenotypic changes reported by Lafranconi et al. (2020) after 90-days of exposure to 6,000 ppm 1,4-dioxane in the same tissues. The transcriptomics analysis further supports overarching study findings demonstrating a non-mutagenic, threshold-based, mitogenic MOA for 1,4-dioxane-induced liver tumors.

Using Network-Based Machine Learning to Predict Transcription Factors Involved in Drought Resistance

Gene regulatory networks underpin stress response pathways in plants. However, parsing these networks to prioritize key genes underlying a particular trait is challenging. Here, we have built the Gene Regulation and Association Network (GRAiN) of rice (Oryza sativa). GRAiN is an interactive query-based web-platform that allows users to study functional relationships between transcription factors (TFs) and genetic modules underlying abiotic-stress responses. We built GRAiN by applying a combination of different network inference algorithms to publicly available gene expression data. We propose a supervised machine learning framework that complements GRAiN in prioritizing genes that regulate stress signal transduction and modulate gene expression under drought conditions. Our framework converts intricate network connectivity patterns of 2160 TFs into a single drought score. We observed that TFs with the highest drought scores define the functional, structural, and evolutionary characteristics of drought resistance in rice. Our approach accurately predicted the function of OsbHLH148 TF, which we validated using in vitro protein-DNA binding assays and mRNA sequencing loss-of-function mutants grown under control and drought stress conditions. Our network and the complementary machine learning strategy lends itself to predicting key regulatory genes underlying other agricultural traits and will assist in the genetic engineering of desirable rice varieties.

A pan-kidney cancer study identifies subtype specific perturbations on pathways with potential drivers in renal cell carcinoma

BACKGROUND

Renal cell carcinoma (RCC) is a complex disease and is comprised of several histological subtypes, the most frequent of which are clear cell renal cell carcinoma (ccRCC), papillary renal cell carcinoma (PRCC) and chromophobe renal cell carcinoma (ChRCC). While lots of studies have been performed to investigate the molecular characterizations of different subtypes of RCC, our knowledge regarding the underlying mechanisms are still incomplete. As molecular alterations are eventually reflected on the pathway level to execute certain biological functions, characterizing the pathway perturbations is crucial for understanding tumorigenesis and development of RCC.

METHODS

In this study, we investigated the pathway perturbations of various RCC subtype against normal tissue based on differential expressed genes within a certain pathway. We explored the potential upstream regulators of subtype-specific pathways with Ingenuity Pathway Analysis (IPA). We also evaluated the relationships between subtype-specific pathways and clinical outcome with survival analysis.

RESULTS

In this study, we carried out a pathway-based analysis to explore the mechanisms of various RCC subtypes with TCGA RNA-seq data. Both commonly altered pathways and subtype-specific pathways were detected. To identify the distinctive characteristics of each subtype, we focused on subtype-specific perturbed pathways. Specifically, we observed that some of the altered pathways were regulated by several recurrent upstream regulators which presenting different expression patterns among distinct RCC subtypes. We also noticed that a large number of perturbed pathways were controlled by the subtype-specific upstream regulators. Moreover, we also evaluated the relationships between perturbed pathways and clinical outcome. Prognostic pathways were identified and their roles in tumor development and progression were inferred.

CONCLUSIONS

In summary, we evaluated the relationships among pathway perturbations, upstream regulators and clinical outcome for differential subtypes in RCC. We hypothesized that the alterations of common upstream regulators as well as subtype-specific upstream regulators work together to affect the downstream pathway perturbations and drive cancer initialization and prognosis. Our findings not only increase our understanding of the mechanisms of various RCC subtypes, but also provide targets for personalized therapeutic intervention.

Traditional Human Populations and Nonhuman Primates Show Parallel Gut Microbiome Adaptations to Analogous Ecological Conditions

The results of this study highlight parallel gut microbiome traits in human and nonhuman primates, depending on subsistence strategy. Although these similarities have been reported before, the functional and ecological bases of this convergence are not fully understood.

Compared with urban-industrial populations, small-scale human communities worldwide share a significant number of gut microbiome traits with nonhuman primates. This overlap is thought to be driven by analogous dietary triggers; however, the ecological and functional bases of this similarity are not fully understood. To start addressing this issue, fecal metagenomes of BaAka hunter-gatherers and traditional Bantu agriculturalists from the Central African Republic were profiled and compared with those of a sympatric western lowland gorilla group (Gorillagorilla gorilla) across two seasons of variable dietary intake. Results show that gorilla gut microbiomes shared similar functional traits with each human group, depending on seasonal dietary behavior. Specifically, parallel microbiome traits were observed between hunter-gatherers and gorillas when the latter consumed more structural polysaccharides during dry seasons, while small-scale agriculturalist and gorilla microbiomes showed significant functional overlap when gorillas consumed more seasonal ripe fruit during wet seasons. Notably, dominance of microbial transporters, transduction systems, and gut xenobiotic metabolism was observed in association with traditional agriculture and energy-dense diets in gorillas at the expense of a functional microbiome repertoire capable of metabolizing more complex polysaccharides. Differential abundance of bacterial taxa that typically distinguish traditional from industrialized human populations (e.g., Prevotella spp.) was also recapitulated in the human and gorilla groups studied, possibly reflecting the degree of polysaccharide complexity included in each group’s dietary niche. These results show conserved functional gut microbiome adaptations to analogous diets in small-scale human populations and nonhuman primates, highlighting the role of plant dietary polysaccharides and diverse environmental exposures in this convergence.

IMPORTANCE The results of this study highlight parallel gut microbiome traits in human and nonhuman primates, depending on subsistence strategy. Although these similarities have been reported before, the functional and ecological bases of this convergence are not fully understood. Here, we show that this parallelism is, in part, likely modulated by the complexity of plant carbohydrates consumed and by exposures to diverse xenobiotics of natural and artificial origin. Furthermore, we discuss how divergence from these parallel microbiome traits is typically associated with adverse health outcomes in human populations living under culturally westernized subsistence patterns. This is important information as we trace the specific dietary and environmental triggers associated with the loss and gain of microbial functions as humans adapt to various dietary niches.

RNAflow: An Effective and Simple RNA-Seq Differential Gene Expression Pipeline Using Nextflow

RNA-Seq enables the identification and quantification of RNA molecules, often with the aim of detecting differentially expressed genes (DEGs). Although RNA-Seq evolved into a standard technique, there is no universal gold standard for these data's computational analysis. On top of that, previous studies proved the irreproducibility of RNA-Seq studies. Here, we present a portable, scalable, and parallelizable Nextflow RNA-Seq pipeline to detect DEGs, which assures a high level of reproducibility. The pipeline automatically takes care of common pitfalls, such as ribosomal RNA removal and low abundance gene filtering. Apart from various visualizations for the DEG results, we incorporated downstream pathway analysis for common species as Homo sapiens and Mus musculus. We evaluated the DEG detection functionality while using qRT-PCR data serving as a reference and observed a very high correlation of the logarithmized gene expression fold changes.

Unraveling Targetable Systemic and Cell-Type-Specific Molecular Phenotypes of Alzheimer's and Parkinson's Brains With Digital Cytometry

Alzheimer's disease (AD) and Parkinson's disease (PD) are the two most common neurodegenerative disorders worldwide, with age being their major risk factor. The increasing worldwide life expectancy, together with the scarcity of available treatment choices, makes it thus pressing to find the molecular basis of AD and PD so that the causing mechanisms can be targeted. To study these mechanisms, gene expression profiles have been compared between diseased and control brain tissues. However, this approach is limited by mRNA expression profiles derived for brain tissues highly reflecting their degeneration in cellular composition but not necessarily disease-related molecular states. We therefore propose to account for cell type composition when comparing transcriptomes of healthy and diseased brain samples, so that the loss of neurons can be decoupled from pathology-associated molecular effects. This approach allowed us to identify genes and pathways putatively altered systemically and in a cell-type-dependent manner in AD and PD brains. Moreover, using chemical perturbagen data, we computationally identified candidate small molecules for specifically targeting the profiled AD/PD-associated molecular alterations. Our approach therefore not only brings new insights into the disease-specific and common molecular etiologies of AD and PD but also, in these realms, foster the discovery of more specific targets for functional and therapeutic exploration.

BAHCC1 binds H3K27me3 via a conserved BAH module to mediate gene silencing and oncogenesis

Trimethylated histone H3 lysine 27 (H3K27me3) regulates gene repression, cell-fate determination and differentiation. We report that a conserved bromo-adjacent homology (BAH) module of BAHCC1 (BAHCC1^BAH) 'recognizes' H3K27me3 specifically and enforces silencing of H3K27me3-demarcated genes in mammalian cells. Biochemical, structural and integrated chromatin immunoprecipitation-sequencing-based analyses demonstrate that direct readout of H3K27me3 by BAHCC1 is achieved through a hydrophobic trimethyl-L-lysine-binding 'cage' formed by BAHCC1^BAH, mediating colocalization of BAHCC1 and H3K27me3-marked genes. BAHCC1 is highly expressed in human acute leukemia and interacts with transcriptional corepressors. In leukemia, depletion of BAHCC1, or disruption of the BAHCC1^BAH-H3K27me3 interaction, causes derepression of H3K27me3-targeted genes that are involved in tumor suppression and cell differentiation, leading to suppression of oncogenesis. In mice, introduction of a germline mutation at Bahcc1 to disrupt its H3K27me3 engagement causes partial postnatal lethality, supporting a role in development. This study identifies an H3K27me3-directed transduction pathway in mammals that relies on a conserved BAH 'reader'.

High temperature stress response is not sexually dimorphic at the whole-body level and is dependent on androgens to induce sex reversal

The understanding of the molecular and endocrine mechanisms behind environmentally-induced sex reversal in fish is of great importance in the context of predicting the potential effects of climate change, especially increasing temperature. Here, we demonstrate the global effects of high temperature on genome-wide transcription in medaka (Oryzias latipes) during early development. Interestingly, data analysis did not show sexual dimorphic changes, demonstrating that thermal stress is not dependent on genotypic sex. Additionally, our results revealed significant changes in several pathways under high temperature, such as stress response from brain, steroid biosynthesis, epigenetic mechanisms, and thyroid hormone biosynthesis, among others. These microarray data raised the question of what the exact molecular and hormonal mechanisms of action are for female-to-male sex reversal under high temperatures in fish. Complementary gene expression analysis revealed that androgen-related genes increase in females (XX) experiencing high water temperature. To test the involvement of androgens in thermal-induced sex reversal, an androgen antagonist was used to treat XX medaka under a high-temperature setup. Data clearly demonstrated failure of female-to-male sex reversal when androgen action is inhibited, corroborating the importance of androgens in environmentally-induced sex reversal.

Ocean acidification induces distinct transcriptomic responses across life history stages of the sea urchin Heliocidaris erythrogramma

Ocean acidification (OA) from seawater uptake of rising carbon dioxide emissions impairs development in marine invertebrates, particularly in calcifying species. Plasticity in gene expression is thought to mediate many of these physiological effects, but how these responses change across life history stages remains unclear. The abbreviated lecithotrophic development of the sea urchin Heliocidaris erythrogramma provides a valuable opportunity to analyse gene expression responses across a wide range of life history stages, including the benthic, post-metamorphic juvenile. We measured the transcriptional response to OA in H. erythrogramma at three stages of the life cycle (embryo, larva, and juvenile) in a controlled breeding design. The results reveal a broad range of strikingly stage-specific impacts of OA on transcription, including changes in the number and identity of affected genes; the magnitude, sign, and variance of their expression response; and the developmental trajectory of expression. The impact of OA on transcription was notably modest in relation to gene expression changes during unperturbed development and much smaller than genetic contributions from parentage. The latter result suggests that natural populations may provide an extensive genetic reservoir of resilience to OA. Taken together, these results highlight the complexity of the molecular response to OA, its substantial life history stage specificity, and the importance of contextualizing the transcriptional response to pH stress in light of normal development and standing genetic variation to better understand the capacity for marine invertebrates to adapt to OA.

ToxicoDB: an integrated database to mine and visualize large-scale toxicogenomic datasets

In the past few decades, major initiatives have been launched around the world to address chemical safety testing. These efforts aim to innovate and improve the efficacy of existing methods with the long-term goal of developing new risk assessment paradigms. The transcriptomic and toxicological profiling of mammalian cells has resulted in the creation of multiple toxicogenomic datasets and corresponding tools for analysis. To enable easy access and analysis of these valuable toxicogenomic data, we have developed ToxicoDB (toxicodb.ca), a free and open cloud-based platform integrating data from large in vitro toxicogenomic studies, including gene expression profiles of primary human and rat hepatocytes treated with 231 potential toxicants. To efficiently mine these complex toxicogenomic data, ToxicoDB provides users with harmonized chemical annotations, time- and dose-dependent plots of compounds across datasets, as well as the toxicity-related pathway analysis. The data in ToxicoDB have been generated using our open-source R package, ToxicoGx (github.com/bhklab/ToxicoGx). Altogether, ToxicoDB provides a streamlined process for mining highly organized, curated, and accessible toxicogenomic data that can be ultimately applied to preclinical toxicity studies and further our understanding of adverse outcomes.

Deciphering rice metabolic flux reprograming under salinity stress via in silico metabolic modeling

Rice is one of the most economically important commodities globally. However, rice plants are salt susceptible species in which high salinity can significantly constrain its productivity. Several physiological parameters in adaptation to salt stress have been observed, though changes in metabolic aspects remain to be elucidated. In this study, rice metabolic activities of salt-stressed flag leaf were systematically characterized. Transcriptomics and metabolomics data were combined to identify disturbed pathways, altered metabolites and metabolic hotspots within the rice metabolic network under salt stress condition. Besides, the feasible flux solutions in different context-specific metabolic networks were estimated and compared. Our findings highlighted metabolic reprogramming in primary metabolic pathways, cellular respiration, antioxidant biosynthetic pathways, and phytohormone biosynthetic pathways. Photosynthesis and hexose utilization were among the major disturbed pathways in the stressed flag leaf. Notably, the increased flux distribution of the photorespiratory pathway could contribute to cellular redox control. Predicted flux statuses in several pathways were consistent with the results from transcriptomics, end-point metabolomics, and physiological studies. Our study illustrated that the contextualized genome-scale model together with multi-omics analysis is a powerful approach to unravel the metabolic responses of rice to salinity stress.

Functional compensation precedes recovery of tissue mass following acute liver injury

The liver plays a central role in metabolism, protein synthesis and detoxification. It possesses unique regenerative capacity upon injury. While many factors regulating cellular proliferation during liver repair have been identified, the mechanisms by which the injured liver maintains vital functions prior to tissue recovery are unknown. Here, we identify a new phase of functional compensation following acute liver injury that occurs prior to cellular proliferation. By coupling single-cell RNA-seq with in situ transcriptional analyses in two independent murine liver injury models, we discover adaptive reprogramming to ensure expression of both injury response and core liver function genes dependent on macrophage-derived WNT/β-catenin signaling. Interestingly, transcriptional compensation is most prominent in non-proliferating cells, clearly delineating two temporally distinct phases of liver recovery. Overall, our work describes a mechanism by which the liver maintains essential physiological functions prior to cellular reconstitution and characterizes macrophage-derived WNT signals required for this compensation.

The liver possesses the ability to regenerate following sudden injury. Here, the authors use single-cell RNA-sequencing and in situ transcriptional analyses to identify a new phase of liver regeneration in mice aimed at maintaining essential functions throughout the regenerative process.

Evolution from adherent to suspension: systems biology of HEK293 cell line development

The need for new safe and efficacious therapies has led to an increased focus on biologics produced in mammalian cells. The human cell line HEK293 has bio-synthetic potential for human-like production attributes and is currently used for manufacturing of several therapeutic proteins and viral vectors. Despite the increased popularity of this strain we still have limited knowledge on the genetic composition of its derivatives. Here we present a genomic, transcriptomic and metabolic gene analysis of six of the most widely used HEK293 cell lines. Changes in gene copy and expression between industrial progeny cell lines and the original HEK293 were associated with cellular component organization, cell motility and cell adhesion. Changes in gene expression between adherent and suspension derivatives highlighted switching in cholesterol biosynthesis and expression of five key genes (RARG, ID1, ZIC1, LOX and DHRS3), a pattern validated in 63 human adherent or suspension cell lines of other origin.

Effects of Glycyrrhizin on Multi-Drug Resistant Pseudomonas aeruginosa

The effects of glycyrrhizin (GLY) on multi-drug resistant (MDR) systemic (MDR9) vs. ocular (B1045) Pseudomonas aeruginosa clinical isolates were determined. Proteomes of each isolate with/without GLY treatment were profiled using liquid chromatography mass spectrometry (LC-MS/MS). The effect of GLY on adherence of MDR isolates to immortalized human (HCET) and mouse (MCEC) corneal epithelial cells, and biofilm and dispersal was tested. Both isolates were treated with GLY (0.25 minimum inhibitory concentration (MIC), 10 mg/mL for MDR9 and 3.75 mg/mL for B1045) and subjected to proteomic analysis. MDR9 had a greater response to GLY (51% of identified proteins affected vs. <1% in B1045). In MDR9 vs. controls, GLY decreased the abundance of proteins for: antibiotic resistance, biofilm formation, and type III secretion. Further, antibiotic resistance and type III secretion proteins had higher control abundances in MDR9 vs. B1045. GLY (5 and 10 mg/mL) significantly reduced binding of both isolates to MCEC, and B1045 to HCET. MDR9 binding to HCET was only reduced at 10 mg/mL GLY. GLY (5 and 10 mg/mL) enhanced dispersal for both isolates, at early (6.5 h) but not later times (24–72 h). This study provides evidence that GLY has a greater effect on the proteome of MDR9 vs. B1045, yet it was equally effective at disrupting adherence and early biofilm dispersal.

Association between breast cancer risk and disease aggressiveness: Characterizing underlying gene expression patterns

The association between breast cancer risk defined by the Tyrer-Cuzick score (TC) and disease prognosis is not well established. Here, we investigated the relationship between 5-year TC and disease aggressiveness and then characterized underlying molecular processes. In a case-only study (n = 2474), we studied the association of TC with molecular subtypes and tumor characteristics. In a subset of patients (n = 672), we correlated gene expression to TC and computed a low-risk TC gene expression (TC-Gx) profile, that is, a profile expected to be negatively associated with risk, which we used to test for association with disease aggressiveness. We performed enrichment analysis to pinpoint molecular processes likely to be altered in low-risk tumors. A higher TC was found to be inversely associated with more aggressive surrogate molecular subtypes and tumor characteristics (P < .05) including Ki-67 proliferation status (P < 5 × 10^-07 ). Our low-risk TC-Gx, based on the weighted sum of 37 expression values of genes strongly correlated with TC, was associated with basal-like (P < 5 × 10^-13 ), HER2-enriched subtype (P < 5 × 10^-07 ) and worse 10-year breast cancer-specific survival (log-rank P < 5 × 10^-04 ). Associations between low-risk TC-Gx and more aggressive molecular subtypes were replicated in an independent cohort from The Cancer Genome Atlas database (n = 975). Gene expression that correlated with low TC was enriched in proliferation and oncogenic signaling pathways (FDR < 0.05). Moreover, higher proliferation was a key factor explaining the association with worse survival. Women who developed breast cancer despite having a low risk were diagnosed with more aggressive tumors and had a worse prognosis, most likely driven by increased proliferation. Our findings imply the need to establish risk factors associated with more aggressive breast cancer subtypes.

A Novel Open Access Web Portal for Integrating Mechanistic and Toxicogenomic Study Results

Applying toxicogenomics to improving the safety profile of drug candidates and crop protection molecules is most useful when it identifies relevant biological and mechanistic information that highlights risks and informs risk mitigation strategies. Pathway-based approaches, such as gene set enrichment analysis, integrate toxicogenomic data with known biological process and pathways. Network methods help define unknown biological processes and offer data reduction advantages. Integrating the 2 approaches would improve interpretation of toxicogenomic information. Barriers to the routine application of these methods in genome-wide transcriptomic studies include a need for “hands-on” computer programming experience, the selection of 1 or more analysis methods (eg pathway analysis methods), the sensitivity of results to algorithm parameters, and challenges in linking differential gene expression to variation in safety outcomes. To facilitate adoption and reproducibility of gene expression analysis in safety studies, we have developed Collaborative Toxicogeomics, an open-access integrated web portal using the Django web framework. The software, developed with the Python programming language, is modular, extensible and implements “best-practice” methods in computational biology. New study results are compared with over 4000 rodent liver experiments from Drug Matrix and open TG-GATEs. A unique feature of the software is the ability to integrate clinical chemistry and histopathology-derived outcomes with results from gene expression studies, leading to relevant mechanistic conclusions. We describe its application by analyzing the effects of several toxicants on liver gene expression and exemplify application to predicting toxicity study outcomes upon chronic treatment from expression changes in acute-duration studies.

Metabolic Responses of Carotenoid and Cordycepin Biosynthetic Pathways in Cordyceps militaris under Light-Programming Exposure through Genome-Wide Transcriptional Analysis

Cordyceps militaris is currently exploited for commercial production of specialty products as its biomass constituents are enriched in bioactive compounds, such as cordycepin. The rational process development is important for economically feasible production of high quality bioproducts. Light is an abiotic factor affecting the cultivation process of this entomopathogenic fungus, particularly in its carotenoid formation. To uncover the cell response to light exposure, this study aimed to systematically investigate the metabolic responses of C. militaris strain TBRC6039 using integrative genome-wide transcriptome and genome-scale metabolic network (GSMN)-driven analysis. The genome-wide transcriptome analysis showed 8747 expressed genes in the glucose and sucrose cultures grown under light-programming and dark conditions. Of them, 689 differentially expressed genes were significant in response to the light-programming exposure. Through integration with the GSMN-driven analysis using the improved network (iRT1467), the reporter metabolites, e.g., adenosine-5'-monophosphate (AMP) and 2-oxoglutarate, were identified when cultivated under the carotenoid-producing condition controlled by light-programming exposure, linking to up-regulations of the metabolic genes involved in glyoxalase system, as well as cordycepin and carotenoid biosynthesis. These results indicated that C. militaris had a metabolic control in acclimatization to light exposure through transcriptional co-regulation, which supported the cell growth and cordycepin production in addition to the accumulation of carotenoid as a photo-protective bio-pigment. This study provides a perspective in manipulating the metabolic fluxes towards the target metabolites through either genetic or physiological approaches.

Xylose Metabolism and the Effect of Oxidative Stress on Lipid and Carotenoid Production in Rhodotorula toruloides: Insights for Future Biorefinery

The use of cell factories to convert sugars from lignocellulosic biomass into chemicals in which oleochemicals and food additives, such as carotenoids, is essential for the shift toward sustainable processes. Rhodotorula toruloides is a yeast that naturally metabolises a wide range of substrates, including lignocellulosic hydrolysates, and converts them into lipids and carotenoids. In this study, xylose, the main component of hemicellulose, was used as the sole substrate for R. toruloides, and a detailed physiology characterisation combined with absolute proteomics and genome-scale metabolic models was carried out to understand the regulation of lipid and carotenoid production. To improve these productions, oxidative stress was induced by hydrogen peroxide and light irradiation and further enhanced by adaptive laboratory evolution. Based on the online measurements of growth and CO2 excretion, three distinct growth phases were identified during batch cultivations. Majority of the intracellular flux estimations showed similar trends with the measured protein levels and demonstrated improved NADPH regeneration, phosphoketolase activity and reduced β-oxidation, correlating with increasing lipid yields. Light irradiation resulted in 70% higher carotenoid and 40% higher lipid content compared to the optimal growth conditions. The presence of hydrogen peroxide did not affect the carotenoid production but culminated in the highest lipid content of 0.65 g/gDCW. The adapted strain showed improved fitness and 2.3-fold higher carotenoid content than the parental strain. This work presents a holistic view of xylose conversion into microbial oil and carotenoids by R. toruloides, in a process toward renewable and cost-effective production of these molecules.

CRISPR interference of nucleotide biosynthesis improves production of a single-domain antibody in Escherichia coli

Growth decoupling can be used to optimize the production of biochemicals and proteins in cell factories. Inhibition of excess biomass formation allows for carbon to be utilized efficiently for product formation instead of growth, resulting in increased product yields and titers. Here, we used CRISPR interference to increase the production of a single‐domain antibody (sdAb) by inhibiting growth during production. First, we screened 21 sgRNA targets in the purine and pyrimidine biosynthesis pathways and found that the repression of 11 pathway genes led to the increased green fluorescent protein production and decreased growth. The sgRNA targets pyrF, pyrG, and cmk were selected and further used to improve the production of two versions of an expression‐optimized sdAb. Proteomics analysis of the sdAb‐producing pyrF, pyrG, and cmk growth decoupling strains showed significantly decreased RpoS levels and an increase of ribosome‐associated proteins, indicating that the growth decoupling strains do not enter stationary phase and maintain their capacity for protein synthesis upon growth inhibition. Finally, sdAb production was scaled up to shake‐flask fermentation where the product yield was improved 2.6‐fold compared to the control strain with no sgRNA target sequence. An sdAb content of 14.6% was reached in the best‐performing pyrG growth decoupling strain.

Radiomics nomogram for prediction disease-free survival and adjuvant chemotherapy benefits in patients with resected stage I lung adenocarcinoma

Background

Robust imaging biomarkers are needed for risk stratification in stage I lung adenocarcinoma patients in order to select optimal treatment regimen. We aimed to construct and validate a radiomics nomogram for predicting the disease-free survival (DFS) of patients with resected stage I lung adenocarcinoma, and further identifying candidates benefit from adjuvant chemotherapy (ACT).

Methods

Using radiomics approach, we analyzed 554 patients' computed tomography (CT) images from three multicenter cohorts. Prognostic radiomics features were extracted from computed tomography (CT) images and selected using least absolute shrinkage and selection operator (LASSO) Cox regression model to build a radiomics signature for DFS stratification. The biological basis of radiomics was explored in the Radiogenomics dataset (n=79) by gene set enrichment analysis (GSEA). Then a nomogram that integrated the signature with these significant clinicopathologic factors in the multivariate analysis were constructed in the training cohort (n=238), and its prognostic accuracy was evaluated in the validation cohort (n=237). Finally, the predictive value of nomogram for ACT benefits was assessed.

Results

The radiomics signature with higher score was significantly associated with worse DFS in both the training and validation cohorts (P<0.001). The GSEA presented that the signature was highly correlated to characteristic metabolic process and immune system during cancer progression. Multivariable analysis revealed that age (P=0.031), pathologic TNM stage (P=0.043), histologic subtype (P=0.010) and the signature (P<0.001) were independently associated with patients' DFS. The integrated radiomics nomogram showed good discrimination performance, as well as good calibration and clinical utility, for DFS prediction in the validation cohort. We further found that the patients with high points (point ≥8.788) defined by the radiomics nomogram obtained a significant favorable response to ACT (P=0.04) while patients with low points (point <8.788) showed no survival difference (P=0.7).

Conclusions

The radiomics nomogram could be used for prognostic prediction and ACT benefits identification for patient with resected stage I lung adenocarcinoma.

Development and Validation of a Gene Signature for Prediction of Relapse in Stage I Testicular Germ Cell Tumors

Background: Testicular germ cell tumors (TGCTs) are commonly diagnosed tumors in young men. However, a satisfactory approach to predict relapse of stage I TGCTs is still lacking. Therefore, this study aimed to develop a robust risk score model for stage I TGCTs. Method: RNA-sequence data of stage I TGCTs and normal testis samples were downloaded and analyzed to identify different expression genes. Gene-based prognostic model was constructed in The Cancer Genome Atlas (TCGA) using least absolute shrinkage and selection operator (LASSO) regression analysis and validated in GSE99420 dataset. Potential biological functions of the genes in prognostic model were determined via Gene Set Enrichment Analysis (GSEA) between high-risk and low-risk patients. Results: A total of 9,391 differentially expressed genes and 84 prognosis-related genes were identified. An eight-gene-based risk score model was constructed to divide patients into high or low risk of relapse. The low-risk patients had a significantly better relapse-free survival (RFS) than high-risk patients in both training and validation cohorts (HR = 0.129, 95% CI = 0.059-0.284, P < 0.001; HR = 0.277, 95% CI = 0.116-0.661, P = 0.004, respectively). The area under the receiver operating characteristic curve (AUC) values at 5 years was 0.805 and 0.724 in the training and validation cohorts, respectively. Functional enrichment analyses showed that DNA replication, ribosome, cell cycle, and TGF-beta signaling pathway may contribute to the relapse process. Conclusion: In summary, our analysis provided a novel eight-gene signature that could predict RFS in stage I TGCT patients.

Estrogens decrease osteoclast number by attenuating mitochondria oxidative phosphorylation and ATP production in early osteoclast precursors

Loss of estrogens at menopause is a major cause of osteoporosis and increased fracture risk. Estrogens protect against bone loss by decreasing osteoclast number through direct actions on cells of the myeloid lineage. Here, we investigated the molecular mechanism of this effect. We report that 17β-estradiol (E₂) decreased osteoclast number by promoting the apoptosis of early osteoclast progenitors, but not mature osteoclasts. This effect was abrogated in cells lacking Bak/Bax-two pro-apoptotic members of the Bcl-2 family of proteins required for mitochondrial apoptotic death. FasL has been previously implicated in the pro-apoptotic actions of E₂. However, we show herein that FasL-deficient mice lose bone mass following ovariectomy indistinguishably from FasL-intact controls, indicating that FasL is not a major contributor to the anti-osteoclastogenic actions of estrogens. Instead, using microarray analysis we have elucidated that ERα-mediated estrogen signaling in osteoclast progenitors decreases "oxidative phosphorylation" and the expression of mitochondria complex I genes. Additionally, E₂ decreased the activity of complex I and oxygen consumption rate. Similar to E₂, the complex I inhibitor Rotenone decreased osteoclastogenesis by promoting osteoclast progenitor apoptosis via Bak/Bax. These findings demonstrate that estrogens decrease osteoclast number by attenuating respiration, and thereby, promoting mitochondrial apoptotic death of early osteoclast progenitors.

Histopathological and Molecular Signatures of a Mouse Model of Acute-on-Chronic Alcoholic Liver Injury Demonstrate Concordance With Human Alcoholic Hepatitis

Human alcoholic hepatitis (AH) carries a high mortality rate. AH is an acute-on-chronic form of liver injury characterized by hepatic steatosis, ballooned hepatocytes, neutrophil infiltration, and pericellular fibrosis. We aimed to study the pathogenesis of AH in an animal model which combines chronic hepatic fibrosis with intragastric alcohol administration. Adult male C57BL6/J mice were treated with CCl4 (0.2 ml/kg, 2×weekly by intraperitoneal injections for 6 weeks) to induce chronic liver fibrosis. Then, ethyl alcohol (up to 25 g/kg/day for 3 weeks) was administered continuously to mice via a gastric feeding tube, with or without one-half dose of CCl4. Liver and serum markers and liver transcriptome were evaluated to characterize acute-on-chronic-alcoholic liver disease in our model. CCl4 or alcohol treatment alone induced liver fibrosis or steatohepatitis, respectively, findings that were consistent with expected pathology. Combined treatment resulted in a marked exacerbation of liver injury, as evident by the development of inflammation, steatosis, and pericellular fibrosis, pathological features of human AH. E. coli and Candida were also detected in livers of mice cotreated with CCl4 and alcohol, indicating pathogen translocation from gut to liver, similar to human AH. Importantly, liver transcriptomic changes specific to combined treatment group demonstrated close concordance with pathways perturbed in patients with severe AH. Overall, mice treated with CCl4 and alcohol displayed key molecular and pathological characteristics of human AH-pericellular fibrosis, increased hepatic bacterial load, and dysregulation of the same molecular pathways. This model may be useful for developing therapeutics for AH.

Deep Proteomic Deconvolution of Interferons and HBV Transfection Effects on a Hepatoblastoma Cell Line

Interferons are commonly utilized in the treatment of chronic hepatitis B virus (HBV) infection but are not effective for all patients. A deep understanding of the limitations of interferon treatment requires delineation of its activity at multiple "omic" levels. While myriad studies have characterized the transcriptomic effects of interferon treatment, surprisingly, few have examined interferon-induced effects at the proteomic level. To remedy this paucity, we stimulated HepG2 cells with both IFN-α and IFN-λ and performed proteomic analysis versus unstimulated cells. Alongside, we examined the effects of HBV transfection in the same cell line, reasoning that parallel IFN and HBV analysis might allow determination of cases where HBV transfection counters the effects of interferons. More than 6000 proteins were identified, with multiple replicates allowing for differential expression analysis at high confidence. Drawing on a compendium of transcriptomic data, as well as proteomic half-life data, we suggest means by which transcriptomic results diverge from our proteomic results. We also invoke a recent multiomic study of HBV-related hepatocarcinoma (HCC), showing that despite HBV's role in initiating HCC, the regulated proteomic landscapes of HBV transfection and HCC do not strongly align. Special focus is applied to the proteasome, with numerous components divergently altered under IFN and HBV-transfection conditions. We also examine alterations of other protein groups relevant to HLA complex peptide display, unveiling intriguing alterations in a number of ubiquitin ligases. Finally, we invoke genome-scale metabolic modeling to predict relevant alterations to the metabolic landscape under experimental conditions. Our data should be useful as a resource for interferon and HBV researchers.

The Presence of Active Brown Adipose Tissue Determines Cold-Induced Energy Expenditure and Oxylipin Profiles in Humans

BACKGROUND

Accumulating evidence links brown adipose tissue (BAT) to increased cold-induced energy expenditure (CIEE) and regulation of lipid metabolism in humans. BAT has also been proposed as a novel source for biologically active lipid mediators including polyunsaturated fatty acids (PUFAs) and oxylipins. However, little is known about cold-mediated differences in energy expenditure and various lipid species between individuals with detectable BAT positive (BATpos) and those without BAT negative (BATneg).

METHODS

Here we investigated a unique cohort of matched BATpos and BATneg individuals identified by 18F-fluorodeoxyglucose positron emission tomography combined with computed tomography ([18F]-FDG PET/CT). BAT function, CIEE, and circulating oxylipins, were analyzed before and after short-term cold exposure using [18F]-FDG PET/CT, indirect calorimetry, and high-resolution mass spectrometry, respectively.

RESULTS

We found that active BAT is the major determinant of CIEE since only BATpos individuals experienced significantly increased energy expenditure in response to cold. A single bout of moderate cold exposure resulted in the dissipation of an additional 20 kcal excess energy in BATpos but not in BATneg individuals. The presence of BAT was associated with a unique systemic PUFA and oxylipin profile characterized by increased levels of anti-inflammatory omega-3 fatty acids as well as cytochrome P450 products but decreased concentrations of some proinflammatory hydroxyeicosatetraenoic acids when compared with BATneg individuals. Notably, cold exposure raised circulating levels of various lipids, including the recently identified BAT-derived circulating factors (BATokines) DiHOME and 12-HEPE, only in BATpos individuals.

CONCLUSIONS

In summary, our data emphasize that BAT in humans is a major contributor toward cold-mediated energy dissipation and a critical organ in the regulation of the systemic lipid pool.

Revealing links between gut microbiome and its fungal community in Type 2 Diabetes Mellitus among Emirati subjects: A pilot study

Type 2 diabetes mellitus (T2DM) drastically affects the population of Middle East countries with an ever-increasing number of overweight and obese individuals. The precise links between T2DM and gut microbiome composition remain elusive in these populations. Here, we performed 16 S rRNA and ITS2- gene based microbial profiling of 50 stool samples from Emirati adults with or without T2DM. The four major enterotypes initially described in westernized cohorts were retrieved in this Emirati population. T2DM and non-T2DM healthy controls had different microbiome compositions, with an enrichment in Prevotella enterotype in non-T2DM controls whereas T2DM individuals had a higher proportion of the dysbiotic Bacteroides 2 enterotype. No significant differences in microbial diversity were observed in T2DM individuals after controlling for cofounding factors, contrasting with reports from westernized cohorts. Interestingly, fungal diversity was significantly decreased in Bacteroides 2 enterotype. Functional profiling from 16 S rRNA gene data showed marked differences between T2DM and non-T2DM controls, with an enrichment in amino acid degradation and LPS-related modules in T2DM individuals, whereas non-T2DM controls had increased abundance of carbohydrate degradation modules in concordance with enterotype composition. These differences provide an insight into gut microbiome composition in Emirati population and its potential role in the development of diabetes mellitus.

MYO5B mutations in pheochromocytoma/paraganglioma promote cancer progression

Identification of additional cancer-associated genes and secondary mutations driving the metastatic progression in pheochromocytoma and paraganglioma (PPGL) is important for subtyping, and may provide optimization of therapeutic regimens. We recently reported novel recurrent nonsynonymous mutations in the MYO5B gene in metastatic PPGL. Here, we explored the functional impact of these MYO5B mutations, and analyzed MYO5B expression in primary PPGL tumor cases in relation to mutation status. Immunohistochemistry and mRNA expression analysis in 30 PPGL tumors revealed an increased MYO5B expression in metastatic compared to non-metastatic cases. In addition, subcellular localization of MYO5B protein was altered from cytoplasmic to membranous in some metastatic tumors, and the strongest and most abnormal expression pattern was observed in a paraganglioma harboring a somatic MYO5B:p.G1611S mutation. In addition to five previously discovered MYO5B mutations, the present study of 30 PPGL (8 previous and 22 new samples) also revealed two, and hence recurrent, mutations in the gene paralog MYO5A. The three MYO5B missense mutations with the highest prediction scores (p.L587P, p.G1611S and p.R1641C) were selected and functionally validated using site directed mutagenesis and stable transfection into human neuroblastoma cells (SK-N-AS) and embryonic kidney cells (HEK293). In vitro analysis showed a significant increased proliferation rate in all three MYO5B mutated clones. The two somatically derived mutations, p.L587P and p.G1611S, were also found to increase the migration rate. Expression analysis of MYO5B mutants compared to wild type clones, demonstrated a significant enrichment of genes involved in migration, proliferation, cell adhesion, glucose metabolism, and cellular homeostasis. Our study validates the functional role of novel MYO5B mutations in proliferation and migration, and suggest the MYO5-pathway to be involved in the malignant progression in some PPGL tumors.

Up to 25% of pheochromocytoma/paraganglioma (PPGL) cases develop metastatic disease with poor outcome and few treatment options. The disease mechanism is not fully understood, and to date there are no reliable markers to predict malignancy. We have recently discovered novel missense mutations in the non-conventional myosin 5 gene (MYO5B), an endosomal transport protein, which we now show enhances progression and migration in PPGLs. MYO5B mutations were preferentially found in patients with metastatic disease and SDH deficiency (germline SDHB-mutations). Abolished SDH activity result in a metabolic switch to aerobic glycolysis requiring increased glucose consumption. Since the MYO5B mutations were found to drive progression through downstream up-regulation of glucose metabolism genes, e.g. glucagon, we hypothesize that these mutations may fuel the pseudohypoxic state by altering glucose uptake in cancer cells. Our result is the first to link the myosin 5 genes to PPGL tumorigenesis. Further, it shows that the tumor progression route in PPGL is complex, with contribution from several genetic factors. An increasing number of studies show dysregulation and importance of the MYO5-proteins in cancer, but little is still known about the precise role and mechanism of mutations, hence more research in this area is needed.

High-throughput, microscope-based sorting to dissect cellular heterogeneity

Microscopy is a powerful tool for characterizing complex cellular phenotypes, but linking these phenotypes to genotype or RNA expression at scale remains challenging. Here, we present Visual Cell Sorting, a method that physically separates hundreds of thousands of live cells based on their visual phenotype. Automated imaging and phenotypic analysis directs selective illumination of Dendra2, a photoconvertible fluorescent protein expressed in live cells; these photoactivated cells are then isolated using fluorescence-activated cell sorting. First, we use Visual Cell Sorting to assess hundreds of nuclear localization sequence variants in a pooled format, identifying variants that improve nuclear localization and enabling annotation of nuclear localization sequences in thousands of human proteins. Second, we recover cells that retain normal nuclear morphologies after paclitaxel treatment, and then derive their single-cell transcriptomes to identify pathways associated with paclitaxel resistance in cancers. Unlike alternative methods, Visual Cell Sorting depends on inexpensive reagents and commercially available hardware. As such, it can be readily deployed to uncover the relationships between visual cellular phenotypes and internal states, including genotypes and gene expression programs.

Glucagon Resistance at the Level of Amino Acid Turnover in Obese Subjects With Hepatic Steatosis

Glucagon secretion is regulated by circulating glucose, but it has turned out that amino acids also play an important role and that hepatic amino acid metabolism and glucagon are linked in a mutual feedback cycle, the liver-α-cell axis. On the basis of this knowledge, we hypothesized that hepatic steatosis might impair glucagon's action on hepatic amino acid metabolism and lead to hyperaminoacidemia and hyperglucagonemia. We subjected 15 healthy lean and 15 obese steatotic male participants to a pancreatic clamp with somatostatin and evaluated hepatic glucose and amino acid metabolism when glucagon was at basal levels and at high physiological levels. The degree of steatosis was evaluated from liver biopsy specimens. Total RNA sequencing of liver biopsy specimens from the obese steatotic individuals revealed perturbations in the expression of genes predominantly involved in amino acid metabolism. This group was characterized by fasting hyperglucagonemia, hyperaminoacidemia, and no lowering of amino acid levels in response to high levels of glucagon. Endogenous glucose production was similar between lean and obese individuals. Our results suggest that hepatic steatosis causes resistance to the effect of glucagon on amino acid metabolism. This results in increased amino acid concentrations and increased glucagon secretion, providing a likely explanation for fatty liver-associated hyperglucagonemia.

Multi-omics systems toxicology study of mouse lung assessing the effects of aerosols from two heat-not-burn tobacco products and cigarette smoke

•

Multi-omics systems toxicology study, comprising five omics data modalities.

•

Multi-Omics Factor Analysis and multi-modality functional network interpretation.

•

Cigarettes smoke (CS) induced complex immunoregulatory interactions across molecular layers.

•

Aerosols from two heat-not-burn tobacco products had less impact on lungs than CS.

Cigarette smoke (CS) causes adverse health effects and, for smoker who do not quit, modified risk tobacco products (MRTPs) can be an alternative to reduce the risk of developing smoking-related diseases. Standard toxicological endpoints can lack sensitivity, with systems toxicology approaches yielding broader insights into toxicological mechanisms. In a 6-month systems toxicology study on ApoE^−/− mice, we conducted an integrative multi-omics analysis to assess the effects of aerosols from the Carbon Heated Tobacco Product (CHTP) 1.2 and Tobacco Heating System (THS) 2.2—a potential and a candidate MRTP based on the heat-not-burn (HnB) principle—compared with CS at matched nicotine concentrations. Molecular exposure effects in the lungs were measured by mRNA/microRNA transcriptomics, proteomics, metabolomics, and lipidomics. Integrative data analysis included Multi-Omics Factor Analysis and multi-modality functional network interpretation. Across all five data modalities, CS exposure was associated with an increased inflammatory and oxidative stress response, and lipid/surfactant alterations. Upon HnB aerosol exposure these effects were much more limited or absent, with reversal of CS-induced effects upon cessation and switching to CHTP 1.2. Functional network analysis revealed CS-induced complex immunoregulatory interactions across the investigated molecular layers (e.g., itaconate, quinolinate, and miR-146) and highlighted the engagement of the heme–Hmox–bilirubin oxidative stress axis by CS. This work exemplifies how multi-omics approaches can be leveraged within systems toxicology studies and the generated multi-omics data set can facilitate the development of analysis methods and can yield further insights into the effects of toxicological exposures on the lung of mice.

Assessment of the Mode of Action Underlying the Effects of GenX in Mouse Liver and Implications for Assessing Human Health Risks

GenX is an alternative to environmentally persistent long-chain perfluoroalkyl and polyfluoroalkyl substances. Mice exposed to GenX exhibit liver hypertrophy, elevated peroxisomal enzyme activity, and other apical endpoints consistent with peroxisome proliferators. To investigate the potential role of peroxisome proliferator-activated receptor alpha (PPARα) activation in mice, and other molecular signals potentially related to observed liver changes, RNA sequencing was conducted on paraffin-embedded liver sections from a 90-day subchronic toxicity study of GenX conducted in mice. Differentially expressed genes were identified for each treatment group, and gene set enrichment analysis was conducted using gene sets that represent biological processes and known canonical pathways. Peroxisome signaling and fatty acid metabolism were among the most significantly enriched gene sets in both sexes at 0.5 and 5 mg/kg GenX; no pathways were enriched at 0.1 mg/kg. Gene sets specific to the PPARα subtype were significantly enriched. These findings were phenotypically anchored to histopathological changes in the same tissue blocks: hypertrophy, mitoses, and apoptosis. In vitro PPARα transactivation assays indicated that GenX activates mouse PPARα. These results indicate that the liver changes observed in GenX-treated mice occur via a mode of action (MOA) involving PPARα, an important finding for human health risk assessment as this MOA has limited relevance to humans.

Cofactor Engineering Redirects Secondary Metabolism and Enhances Erythromycin Production in Saccharopolyspora erythraea

Saccharopolyspora erythraea is used for industrial erythromycin production. To explore the physiological role of intracellular energy state in metabolic regulation by S. erythraea, we initially overexpressed the F₁ part of the endogenous F₁F₀-ATPase in the high yielding erythromycin producing strain E3. The F₁-ATPase expression resulted in lower [ATP]/[ADP] ratios, which was accompanied by a strong increase in the production of a reddish pigment and a decreased erythromycin production. Subsequent transcriptional analysis revealed that the lower intracellular [ATP]/[ADP] ratios exerted a pleotropic regulation on the metabolism of S. erythraea. The lower [ATP]/[ADP] ratios induced physiological changes to restore the energy balance, mainly via pathways that tend to produce ATP or regenerate NADH. The F₁-ATPase overexpression strain exhibited a state of redox stress, which was correlated to an alteration of electron transport at the branch of the terminal oxidases, and S. erythraea channeled the enhanced glycolytic flux toward a reddish pigment in order to reduce NADH formation. The production of erythromycin was decreased, which is in accordance with the net ATP requirement and the excess NADH formed through this pathway. Partial growth inhibition by apramycin increased the intracellular [ATP]/[ADP] ratios and demonstrated a positive correlation between [ATP]/[ADP] ratios and erythromycin synthesis. Finally, overexpression of the entire F₁F₀-ATPase complex resulted in 28% enhanced erythromycin production and markedly reduced pigment synthesis in E3. The work illustrates a feasible strategy to optimize the distribution of fluxes in secondary metabolism.

The reduction of DSS-induced colitis severity in mice exposed to cigarette smoke is linked to immune modulation and microbial shifts

Exposure to cigarette smoke (CS) causes detrimental health effects, increasing the risk of cardiovascular, pulmonary diseases and carcinogenesis in exposed individuals. The impact of CS on Inflammatory Bowel Disease (IBD) has been established by a number of epidemiological and clinical studies. In fact, CS is associated with a higher risk of developing Crohn's disease (CD) while inversely correlates with the development, disease risks, and relapse rate of ulcerative colitis (UC). To investigate the effect of CS exposure on experimental colitis, we performed a comprehensive and integrated comparative analysis of colon transcriptome and microbiome in mice exposed to dextran sodium sulfate (DSS) and CS. Colon transcriptome analysis revealed that CS downregulated specific pathways in a concentration-dependent manner, affecting both the inflammatory state and composition of the gut microbiome. Metagenomics analysis demonstrated that CS can modulate DSS-induced dysbiosis of specific bacterial genera, contributing to resolve the inflammation or accelerate recovery. The risks of smoking far outweigh any possible benefit, thus smoking cessation must always be encouraged because of its significant health benefits. However, the inverse association between active smoking and the development of UC cannot be ignored and the present study lays the foundation for investigating potential molecular mechanisms responsible for the attenuation of colitis by certain compounds of tobacco when decoupled from combustion.

Transcriptomic abnormalities in peripheral blood in bipolar disorder, and discrimination of the major psychoses

We performed a transcriptome-wide meta-analysis and gene co-expression network analysis to identify genes and gene networks dysregulated in the peripheral blood of bipolar disorder (BD) cases relative to unaffected comparison subjects, and determined the specificity of the transcriptomic signatures of BD and schizophrenia (SZ). Nineteen genes and 4 gene modules were significantly differentially expressed in BD cases. Thirteen gene modules were shown to be differentially expressed in a combined case-group of BD and SZ subjects called "major psychosis", including genes biologically linked to apoptosis, reactive oxygen, chromatin remodeling, and immune signaling. No modules were differentially expressed between BD and SZ cases. Machine-learning classifiers trained to separate diagnostic classes based solely on gene expression profiles could distinguish BD cases from unaffected comparison subjects with an area under the curve (AUC) of 0.724, as well as BD cases from SZ cases with AUC = 0.677 in withheld test samples. We introduced a novel and straightforward method called "polytranscript risk scoring" that could distinguish BD cases from unaffected subjects (AUC = 0.672) and SZ cases (AUC = 0.607) significantly better than expected by chance. Taken together, our results highlighted gene expression alterations common to BD and SZ that involve biological processes of inflammation, oxidative stress, apoptosis, and chromatin regulation, and highlight disorder-specific changes in gene expression that discriminate the major psychoses.

FMN reduces Amyloid-β toxicity in yeast by regulating redox status and cellular metabolism

Alzheimer's disease (AD) is defined by progressive neurodegeneration, with oligomerization and aggregation of amyloid-β peptides (Aβ) playing a pivotal role in its pathogenesis. In recent years, the yeast Saccharomyces cerevisiae has been successfully used to clarify the roles of different human proteins involved in neurodegeneration. Here, we report a genome-wide synthetic genetic interaction array to identify toxicity modifiers of Aβ42, using yeast as the model organism. We find that FMN1, the gene encoding riboflavin kinase, and its metabolic product flavin mononucleotide (FMN) reduce Aβ42 toxicity. Classic experimental analyses combined with RNAseq show the effects of FMN supplementation to include reducing misfolded protein load, altering cellular metabolism, increasing NADH/(NADH + NAD+) and NADPH/(NADPH + NADP+) ratios and increasing resistance to oxidative stress. Additionally, FMN supplementation modifies Htt103QP toxicity and α-synuclein toxicity in the humanized yeast. Our findings offer insights for reducing cytotoxicity of Aβ42, and potentially other misfolded proteins, via FMN-dependent cellular pathways.

EFMviz: A COBRA Toolbox extension to visualize Elementary Flux Modes in Genome-Scale Metabolic Models

Elementary Flux Modes (EFMs) are a tool for constraint-based modeling and metabolic network analysis. However, systematic and automated visualization of EFMs, capable of integrating various data types is still a challenge. In this study, we developed an extension for the widely adopted COBRA Toolbox, EFMviz, for analysis and graphical visualization of EFMs as networks of reactions, metabolites and genes. The analysis workflow offers a platform for EFM visualization to improve EFM interpretability by connecting COBRA toolbox with the network analysis and visualization software Cytoscape. The biological applicability of EFMviz is demonstrated in two use cases on medium (Escherichia coli, iAF1260) and large (human, Recon 2.2) genome-scale metabolic models. EFMviz is open-source and integrated into COBRA Toolbox. The analysis workflows used for the two use cases are detailed in the two tutorials provided with EFMviz along with the data used in this study.

Probing Carbon Utilization of Cordyceps militaris by Sugar Transportome and Protein Structural Analysis

Beyond comparative genomics, we identified 85 sugar transporter genes in Cordyceps militaris, clustering into nine subfamilies as sequence- and phylogenetic-based functional classification, presuming the versatile capability of the fungal growths on a range of sugars. Further analysis of the global gene expression patterns of C. militaris showed 123 genes were significantly expressed across the sucrose, glucose, and xylose cultures. The sugar transporters specific for pentose were then identified by gene-set enrichment analysis. Of them, the putative pentose transporter, CCM_06358 gene, was highest expressed in the xylose culture, and its functional role in xylose transport was discovered by the analysis of conserved structural motifs. In addition, a battery of molecular modeling methods, including homology modeling, transport pathway analysis, residue interaction network combined with molecular mechanics Poisson–Boltzmann surface area simulation (MM-PBSA), was implemented for probing the structure and function of the selected pentose transporter (CCM_06358) as a representative of sugar transportome in C. militaris. Considering the network bottlenecks and structural organizations, we further identified key amino acids (Phe38 and Trp441) and their interactions with other residues, contributing the xylose transport function, as verified by binding free energy calculation. The strategy used herein generated remarkably valuable biological information, which is applicable for the study of sugar transportome and the structure engineering of targeted transporter proteins that might link to the production of bioactive compounds derived from xylose metabolism, such as cordycepin.

Multi-omics characterization of a diet-induced obese model of non-alcoholic steatohepatitis

To improve the understanding of the complex biological processes underlying the development of non-alcoholic steatohepatitis (NASH), a multi-omics approach combining bulk RNA-sequencing based transcriptomics, quantitative proteomics and single-cell RNA-sequencing was used to characterize tissue biopsies from histologically validated diet-induced obese (DIO) NASH mice compared to chow-fed controls. Bulk RNA-sequencing and proteomics showed a clear distinction between phenotypes and a good correspondence between mRNA and protein level regulations, apart from specific regulatory events discovered by each technology. Transcriptomics-based gene set enrichment analysis revealed changes associated with key clinical manifestations of NASH, including impaired lipid metabolism, increased extracellular matrix formation/remodeling and pro-inflammatory responses, whereas proteomics-based gene set enrichment analysis pinpointed metabolic pathway perturbations. Integration with single-cell RNA-sequencing data identified key regulated cell types involved in development of NASH demonstrating the cellular heterogeneity and complexity of NASH pathogenesis.

A 6-month systems toxicology inhalation study in ApoE-/- mice demonstrates reduced cardiovascular effects of E-vapor aerosols compared with cigarette smoke

Smoking cigarettes is harmful to the cardiovascular system. Considerable attention has been paid to the reduced harm potential of alternative nicotine-containing inhalable products such as e-cigarettes. We investigated the effects of E-vapor aerosols or cigarette smoke (CS) on atherosclerosis progression, cardiovascular function, and molecular changes in the heart and aorta of female apolipoprotein E-deficient (ApoE^-/-) mice. The mice were exposed to aerosols from three different E-vapor formulations: 1) carrier (propylene glycol and vegetable glycerol), 2) base (carrier and nicotine), or 3) test (base and flavor) or to CS from 3R4F reference cigarettes for up to 6 mo. Concentrations of CS and base or test aerosols were matched at 35 µg nicotine/L. Exposure to CS, compared with sham-exposed fresh air controls, accelerated atherosclerotic plaque formation, whereas no such effect was seen for any of the three E-vapor aerosols. Molecular changes indicated disease mechanisms related to oxidative stress and inflammation in general, plus changes in calcium regulation, and altered cytoskeletal organization and microtubule dynamics in the left ventricle. While ejection fraction, fractional shortening, cardiac output, and isovolumic contraction time remained unchanged following E-vapor aerosols exposure, the nicotine-containing base and test aerosols caused an increase in isovolumic relaxation time similar to CS. A nicotine-related increase in pulse wave velocity and arterial stiffness was also observed, but it was significantly lower for base and test aerosols than for CS. These results demonstrate that in comparison with CS, E-vapor aerosols induce substantially lower biological responses associated with smoking-related cardiovascular diseases.NEW & NOTEWORTHY Analysis of key urinary oxidative stress markers and proinflammatory cytokines showed an absence of oxidative stress and inflammation in the animals exposed to E-vapor aerosols. Conversely, animals exposed to conventional cigarette smoke had high urinary levels of these markers. When compared with conventional cigarette smoke, E-vapor aerosols induced smaller atherosclerotic plaque surface area and volume. Systolic and diastolic cardiac function, as well as endothelial function, were further significantly less affected by electronic cigarette aerosols than conventional cigarette smoke. Molecular analysis demonstrated that E-vapor aerosols induce significantly smaller transcriptomic dysregulation in the heart and aorta compared with conventional cigarette smoke.

Physiological responses of Saccharomyces cerevisiae to industrially relevant conditions: Slow growth, low pH, and high CO2 levels

Engineered strains of Saccharomyces cerevisiae are used for industrial production of succinic acid. Optimal process conditions for dicarboxylic‐acid yield and recovery include slow growth, low pH, and high CO₂. To quantify and understand how these process parameters affect yeast physiology, this study investigates individual and combined impacts of low pH (3.0) and high CO₂ (50%) on slow‐growing chemostat and retentostat cultures of the reference strain S. cerevisiae CEN.PK113‐7D. Combined exposure to low pH and high CO₂ led to increased maintenance‐energy requirements and death rates in aerobic, glucose‐limited cultures. Further experiments showed that these effects were predominantly caused by low pH. Growth under ammonium‐limited, energy‐excess conditions did not aggravate or ameliorate these adverse impacts. Despite the absence of a synergistic effect of low pH and high CO₂ on physiology, high CO₂ strongly affected genome‐wide transcriptional responses to low pH. Interference of high CO₂ with low‐pH signaling is consistent with low‐pH and high‐CO₂ signals being relayed via common (MAPK) signaling pathways, notably the cell wall integrity, high‐osmolarity glycerol, and calcineurin pathways. This study highlights the need to further increase robustness of cell factories to low pH for carboxylic‐acid production, even in organisms that are already applied at industrial scale.

A multi-omics investigation of the molecular characteristics and classification of six metabolic syndrome relevant diseases

Metabolic syndrome (MTS) is a cluster of concurrent metabolic abnormal conditions. MTS and its component metabolic diseases are heterogeneous and closely related, making their relationships complicated, thus hindering precision treatment. Methods: We collected seven groups of samples (group a: healthy individuals; group b: obesity; group c: MTS; group d: hyperglycemia, group e: hypertension, group f: hyperlipidemia; group g: type II diabetes, n=7 for each group). We examined the molecular characteristics of each sample by metabolomic, proteomic and peptidomic profiling analysis. The differential molecules (including metabolites, proteins and peptides) between each disease group and the healthy group were recognized by statistical analyses. Furthermore, a two-step clustering workflow which combines multi-omics and clinical information was used to redefine molecularly and clinically differential groups. Meanwhile, molecular, clinical, network and pathway based analyses were used to identify the group-specific biological features. Results: Both shared and disease-specific molecular profiles among the six types of diseases were identified. Meanwhile, the patients were stratified into three distinct groups which were different from original disease definitions but presented significant differences in glucose and lipid metabolism (Group 1: relatively favorable metabolic conditions; Group 2: severe dyslipidemia; Group 3: dysregulated insulin and glucose). Group specific biological signatures were also systematically described. The dyslipidemia group showed higher levels in multiple lipid metabolites like phosphatidylserine and phosphatidylcholine, and showed significant up-regulations in lipid and amino acid metabolism pathways. The glucose dysregulated group showed higher levels in many polypeptides from proteins contributing to immune response. The another group, with better glucose/lipid metabolism ability, showed higher levels in lipid regulating enzymes like the lecithin cholesterol acyltransferase and proteins involved in complement and coagulation cascades. Conclusions: This multi-omics based study provides a general view of the complex relationships and an alternative classification for various metabolic diseases where the cross-talk or compensatory mechanism between the immune and metabolism systems plays a critical role.

Structural, functional, and molecular impact on the cardiovascular system in ApoE-/- mice exposed to aerosol from candidate modified risk tobacco products, Carbon Heated Tobacco Product 1.2 and Tobacco Heating System 2.2, compared with cigarette smoke

AIM

To investigate the molecular, structural, and functional impact of aerosols from candidate modified risk tobacco products (cMRTP), the Carbon Heated Tobacco Product (CHTP) 1.2 and Tobacco Heating System (THS) 2.2, compared with that of mainstream cigarette smoke (CS) on the cardiovascular system of ApoE-/- mice.

METHODS

Female ApoE-/- mice were exposed to aerosols from THS 2.2 and CHTP 1.2 or to CS from the 3R4F reference cigarette for up to 6 months at matching nicotine concentrations. A Cessation and a Switching group (3 months exposure to 3R4F CS followed by filtered air or CHTP 1.2 for 3 months) were included. Cardiovascular effects were investigated by echocardiographic, histopathological, immunohistochemical, and transcriptomics analyses.

RESULTS

Continuous exposure to cMRTP aerosols did not affect atherosclerosis progression, heart function, left ventricular (LV) structure, or the cardiovascular transcriptome. Exposure to 3R4F CS triggered atherosclerosis progression, reduced systolic ejection fraction and fractional shortening, caused heart LV hypertrophy, and initiated significant dysregulation in the transcriptomes of the heart ventricle and thoracic aorta. Importantly, the structural, functional, and molecular changes caused by 3R4F CS were improved in the smoking cessation and switching groups.

CONCLUSION

Exposure to cMRTP aerosols lacked most of the CS exposure-related functional, structural, and molecular effects. Smoking cessation or switching to CHTP 1.2 aerosol caused similar recovery from the 3R4F CS effects in the ApoE-/- model, with no further acceleration of plaque progression beyond the aging-related rate.

The coding and non-coding transcriptional landscape of subependymal giant cell astrocytomas

Tuberous sclerosis complex (TSC) is an autosomal dominantly inherited neurocutaneous disorder caused by inactivating mutations in TSC1 or TSC2, key regulators of the mechanistic target of rapamycin complex 1 (mTORC1) pathway. In the CNS, TSC is characterized by cortical tubers, subependymal nodules and subependymal giant cell astrocytomas (SEGAs). SEGAs may lead to impaired circulation of CSF resulting in hydrocephalus and raised intracranial pressure in patients with TSC. Currently, surgical resection and mTORC1 inhibitors are the recommended treatment options for patients with SEGA. In the present study, high-throughput RNA-sequencing (SEGAs n = 19, periventricular control n = 8) was used in combination with computational approaches to unravel the complexity of SEGA development. We identified 9400 mRNAs and 94 microRNAs differentially expressed in SEGAs compared to control tissue. The SEGA transcriptome profile was enriched for the mitogen-activated protein kinase (MAPK) pathway, a major regulator of cell proliferation and survival. Analysis at the protein level confirmed that extracellular signal-regulated kinase (ERK) is activated in SEGAs. Subsequently, the inhibition of ERK independently of mTORC1 blockade decreased efficiently the proliferation of primary patient-derived SEGA cultures. Furthermore, we found that LAMTOR1, LAMTOR2, LAMTOR3, LAMTOR4 and LAMTOR5 were overexpressed at both gene and protein levels in SEGA compared to control tissue. Taken together LAMTOR1-5 can form a complex, known as the 'Ragulator' complex, which is known to activate both mTORC1 and MAPK/ERK pathways. Overall, this study shows that the MAPK/ERK pathway could be used as a target for treatment independent of, or in combination with mTORC1 inhibitors for TSC patients. Moreover, our study provides initial evidence of a possible link between the constitutive activated mTORC1 pathway and a secondary driver pathway of tumour growth.

Postmortem Cortical Transcriptomics of Lewy Body Dementia Reveal Mitochondrial Dysfunction and Lack of Neuroinflammation

OBJECTIVE

Prevalence of Lewy body dementias (LBD) is second only to Alzheimer's disease (AD) among people with neurodegenerative dementia. LBD cause earlier mortality, more intense neuropsychiatric symptoms, more caregivers' burden, and higher costs than AD. The molecular mechanisms underlying LBD are largely unknown. As advancing molecular level mechanistic understanding is essential for identifying reliable peripheral biomarkers and novel therapeutic targets for LBD, the authors aimed to identify differentially expressed genes (DEG), and dysfunctional molecular networks in postmortem LBD brains.

METHODS

The authors investigated the transcriptomics of postmortem anterior cingulate and dorsolateral prefrontal cortices of people with pathology-verified LBD using next-generation RNA-sequencing. The authors verified the identified DEG using high-throughput quantitative polymerase chain reactions. Functional implications of identified DEG and the consequent metabolic reprogramming were evaluated by Ingenuity pathway analyses, genome-scale metabolic modeling, reporter metabolite analyses, and in silico gene silencing.

RESULTS

The authors identified and verified 12 novel DEGs (MPO, SELE, CTSG, ALPI, ABCA13, GALNT6, SST, RBM3, CSF3, SLC4A1, OXTR, and RAB44) in LBD brains with genome-wide statistical significance. The authors documented statistically significant down-regulation of several cytokine genes. Identified dysfunctional molecular networks highlighted the contributions of mitochondrial dysfunction, oxidative stress, and immunosenescence toward neurodegeneration in LBD.

CONCLUSION

Our findings support that chronic microglial activation and neuroinflammation, well-documented in AD, are notably absent in LBD. The lack of neuroinflammation in LBD brains was corroborated by statistically significant down-regulation of several inflammatory markers. Identified DEGs, especially down-regulated inflammatory markers, may aid distinguishing LBD from AD, and their biomarker potential warrant further investigation.

Effects of overexpression of STB5 in Saccharomyces cerevisiae on fatty acid biosynthesis, physiology and transcriptome

Microbial conversion of biomass to fatty acids (FA) and products derived thereof is an attractive alternative to the traditional oleochemical production route from animal and plant lipids. This study examined if NADPH-costly FA biosynthesis could be enhanced by overexpressing the transcription factor Stb5 in Saccharomyces cerevisiae. Stb5 activates expression of multiple genes encoding enzymes within the pentose phosphate pathway (PPP) and other NADPH-producing reactions. Overexpression of STB5 led to a decreased growth rate and an increased free fatty acid (FFA) production during growth on glucose. The improved FFA synthetic ability in the glucose phase was shown to be independent of flux through the oxidative PPP. RNAseq analysis revealed that STB5 overexpression had wide-ranging effects on the transcriptome in the batch phase, and appeared to cause a counterintuitive phenotype with reduced flux through the oxidative PPP. During glucose limitation, when an increased NADPH supply is likely less harmful, an overall induction of the proposed target genes of Stb5 (eg. GND1/2, TAL1, ALD6, YEF1) was observed. Taken together, the strategy of utilizing STB5 overexpression to increase NADPH supply for reductive biosynthesis is suggested to have potential in strains engineered to have strong ability to consume excess NADPH, alleviating a potential redox imbalance.