pathfindR: An R Package for Comprehensive Identification of Enriched Pathways in Omics Data Through Active Subnetworks

Plasma proteomic signature predicts who will get persistent symptoms following SARS-CoV-2 infection

BACKGROUND

The majority of those infected by ancestral Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) during the UK first wave (starting March 2020) did not require hospitalisation. Most had a short-lived mild or asymptomatic infection, while others had symptoms that persisted for weeks or months. We hypothesized that the plasma proteome at the time of first infection would reflect differences in the inflammatory response that linked to symptom severity and duration.

METHODS

We performed a nested longitudinal case-control study and targeted analysis of the plasma proteome of 156 healthcare workers (HCW) with and without lab confirmed SARS-CoV-2 infection. Targeted proteomic multiple-reaction monitoring analysis of 91 pre-selected proteins was undertaken in uninfected healthcare workers at baseline, and in infected healthcare workers serially, from 1 week prior to 6 weeks after their first confirmed SARS-CoV-2 infection. Symptom severity and antibody responses were also tracked. Questionnaires at 6 and 12 months collected data on persistent symptoms.

FINDINGS

Within this cohort (median age 39 years, interquartile range 30-47 years), 54 healthcare workers (44% male) had PCR or antibody confirmed infection, with the remaining 102 (38% male) serving as uninfected controls. Following the first confirmed SARS-CoV-2 infection, perturbation of the plasma proteome persisted for up to 6 weeks, tracking symptom severity and antibody responses. Differentially abundant proteins were mostly coordinated around lipid, atherosclerosis and cholesterol metabolism pathways, complement and coagulation cascades, autophagy, and lysosomal function. The proteomic profile at the time of seroconversion associated with persistent symptoms out to 12 months. Data are available via ProteomeXchange with identifier PXD036590.

INTERPRETATION

Our findings show that non-severe SARS-CoV-2 infection perturbs the plasma proteome for at least 6 weeks. The plasma proteomic signature at the time of seroconversion has the potential to identify which individuals are more likely to suffer from persistent symptoms related to SARS-CoV-2 infection.

FUNDING INFORMATION

The COVIDsortium is supported by funding donated by individuals, charitable Trusts, and corporations including Goldman Sachs, Citadel and Citadel Securities, The Guy Foundation, GW Pharmaceuticals, Kusuma Trust, and Jagclif Charitable Trust, and enabled by Barts Charity with support from University College London Hospitals (UCLH) Charity. This work was additionally supported by the Translational Mass Spectrometry Research Group and the Biomedical Research Center (BRC) at Great Ormond Street Hospital.

Proteomics endotyping of infants with severe bronchiolitis and risk of childhood asthma

BACKGROUND

Bronchiolitis is the leading cause of hospitalization in U.S. infants and a major risk factor for childhood asthma. Growing evidence supports clinical heterogeneity within bronchiolitis. We aimed to identify endotypes of infant bronchiolitis by integrating clinical, virus, and serum proteome data, and examine their relationships with asthma development.

METHODS

This is a multicenter prospective cohort study of infants hospitalized for physician-diagnosis of bronchiolitis. We identified bronchiolitis endotypes by applying unsupervised machine learning (clustering) approaches to integrated clinical, virus (respiratory syncytial virus [RSV], rhinovirus [RV]), and serum proteome data measured at hospitalization. We then examined their longitudinal association with the risk for developing asthma by age 6 years.

RESULTS

In 140 infants hospitalized with bronchiolitis, we identified three endotypes: (1) clinicalatopic virusRV proteomeNFκB-dysregulated , (2) clinicalnon-atopic virusRSV/RV proteomeTNF-dysregulated , and (3) clinicalclassic virusRSV proteomeNFκB/TNF-regulated endotypes. Endotype 1 infants were characterized by high proportion of IgE sensitization and RV infection. These endotype 1 infants also had dysregulated NFκB pathways (FDR < 0.001) and significantly higher risks for developing asthma (53% vs. 22%; adjOR 4.04; 95% CI, 1.49-11.0; p = 0.006), compared with endotype 3 (clinically resembling "classic" bronchiolitis). Likewise, endotype 2 infants were characterized by low proportion of IgE sensitization and high proportion of RSV or RV infection. These endotype 2 infants had dysregulated tumor necrosis factor (TNF)-mediated signaling pathway (FDR <0.001) and significantly higher risks for developing asthma (44% vs. 22%; adjOR 2.71; 95% CI, 1.03-7.11, p = 0.04).

CONCLUSION

In this multicenter cohort, integrated clustering of clinical, virus, and proteome data identified biologically distinct endotypes of bronchiolitis that have differential risks of asthma development.

Chemical inhibition of oxygen-sensing prolyl hydroxylases impairs angiogenic competence of human vascular endothelium through metabolic reprogramming

Endothelial cell (EC) metabolism has emerged as a driver of angiogenesis. While hypoxia inactivates the oxygen sensors prolyl-4 hydroxylase domain-containing proteins 1-3 (PHD1-3) and stimulates angiogenesis, the effects of PHDs on EC functions remain poorly defined. Here, we investigated the impact of chemical PHD inhibition by dimethyloxalylglycine (DMOG) on angiogenic competence and metabolism of human vascular ECs. DMOG reduced EC proliferation, migration, and tube formation capacities, responses that were associated with an unfavorable metabolic reprogramming. While glycolytic genes were induced, multiple genes encoding sub-units of mitochondrial complex I were suppressed with concurrent decline in nicotinamide adenine dinucleotide (NAD+) levels. Importantly, the DMOG-induced defects in EC migration could be partially rescued by augmenting NAD+ levels through nicotinamide riboside or citrate supplementation. In summary, by integrating functional assays, transcriptomics, and metabolomics, we provide insights into the effects of PHD inhibition on angiogenic competence and metabolism of human vascular ECs.

Survey for Computer-Aided Tools and Databases in Metabolomics

Metabolomics has advanced from innovation and functional genomics tools and is currently a basis in the big data-led precision medicine era. Metabolomics is promising in the pharmaceutical field and clinical research. However, due to the complexity and high throughput data generated from such experiments, data mining and analysis are significant challenges for researchers in the field. Therefore, several efforts were made to develop a complete workflow that helps researchers analyze data. This paper introduces a review of the state-of-the-art computer-aided tools and databases in metabolomics established in recent years. The paper provides computational tools and resources based on functionality and accessibility and provides hyperlinks to web pages to download or use. This review aims to present the latest computer-aided tools, databases, and resources to the metabolomics community in one place.

Cellular and Molecular Variations in Male and Female Murine Skeletal Muscle after Long-Term Feeding with a High-Fat Diet

Current information regarding the effects of a high-fat diet (HFD) on skeletal muscle is contradictory. This study aimed to investigate the effects of a long-term HFD on skeletal muscle in male and female mice at the morphological, cellular, and molecular levels. Adult mice of the C57BL/6 strain were fed standard chow or an HFD for 20 weeks. The tibialis anterior muscles were dissected, weighed, and processed for cellular and molecular analyses. Immunocytochemical and morphometric techniques were applied to quantify fiber size, satellite cells (SCs), and myonuclei. Additionally, PCR array and RT-qPCR tests were performed to determine the expression levels of key muscle genes. Muscles from HFD mice showed decreases in weight, SCs, and myonuclei, consistent with the atrophic phenotype. This atrophy was associated with a decrease in the percentage of oxidative fibers within the muscle. These findings were further confirmed by molecular analyses that showed significant reductions in the expression of Pax7, Myh1, and Myh2 genes and increased Mstn gene expression. Male and female mice showed similar trends in response to HFD-induced obesity. These findings indicate that the long-term effects of obesity on skeletal muscle resemble those of age-related sarcopenia.

Transcriptome analysis provides novel insights into the soil amendments induced response in continuously cropped Codonopsis tangshen

Codonopsis tangshen Oliv (C. tangshen) is an important Chinese traditional medicinal plant with various health benefits. However, the growth of C. tangshen are seriously affected by continuous cropping, which led to the decrease of the yield and quality. A field experiment was conducted to learn the effects of soil amendments on the growth of C. tangshen under continuous cropping condition, and the biological events which occurred at molecular level were investigated. The results indicated that the content of chlorophyll a (Chl a), chlorophyll b (Chl b), and carotenoid (Car) was significantly higher in SCPM (silicon-calcium-potassium-magnesium fertilizer), SCPMA (SCPM combined with azoxystrobin) and SCPMAOM (SCPM combined with azoxystrobin and organic manure) treatments. Moreover, the yield and the levels of alkaloid, polysaccharide, flavone and total protein in the treatments of SCPM, SCPMA and SCPMAOM were significantly higher than those in the control, and these indexes were all highest in the SCPMAOM treatment. RNA-sequencing (RNA-Seq) is an economical and efficient method to obtain genetic information for species with or without available genome data. In this study, RNA-Seq was performed to understand how continuously cropped C. tangshen responded to the soil amendments at the transcriptome level. The number of differentially expressed genes (DEGs) were as follows: CK vs. SCPM (719 up- and 1456 down-), CK vs. SCPMA (1302 up- and 1748 down-), CK vs. SCPMAOM (1274 up- and 1678 down-). The soil amendments affected the growth of C. tangshen mainly by regulating the genes involved in pathways of 'photosynthesis,' 'plant hormone signal transduction,' 'biosynthesis of unsaturated fatty acids,' 'phenylpropanoid biosynthesis,' and 'starch and sucrose metabolism,' etc. qRT-PCR was performed to validate the expressions of 10 target genes such as CP26, PsaF, and POX, etc., which verified the reliability of RNA-Seq results. Overall, this study revealed the roles and underlying mechanisms of the soil amendments in regulating the growth of continuously cropped C. tangshen at transcriptome level. These findings are beneficial for improving the continuous cropping tolerance and may be valuable for future genetic improvement of C. tangshen.

Distinct stromal and immune cell interactions shape the pathogenesis of rheumatoid and psoriatic arthritis

OBJECTIVES

Immune and stromal cell communication is central in the pathogenesis of rheumatoid arthritis (RA) and psoriatic arthritis (PsA), however, the nature of these interactions in the synovial pathology of the two pathotypes can differ. Identifying immune-stromal cell crosstalk at the site of inflammation in RA and PsA is challenging. This study creates the first global transcriptomic analysis of the RA and PsA inflamed joint and investigates immune-stromal cell interactions in the pathogenesis of synovial inflammation.

METHODS

Single cell transcriptomic profiling of 178 000 synovial tissue cells from five patients with PsA and four patients with RA, importantly, without prior sorting of immune and stromal cells. This approach enabled the transcriptomic analysis of the intact synovial tissue and identification of immune and stromal cell interactions. State of the art data integration and annotation techniques identified and characterised 18 stromal and 14 immune cell clusters.

RESULTS

Global transcriptomic analysis of synovial cell subsets identifies actively proliferating synovial T cells and indicates that due to differential λ and κ immunoglobulin light chain usage, synovial plasma cells are potentially not derived from the local memory B cell pool. Importantly, we report distinct fibroblast and endothelial cell transcriptomes indicating abundant subpopulations in RA and PsA characterised by differential transcription factor usage. Using receptor-ligand interactions and downstream target characterisation, we identify RA-specific synovial T cell-derived transforming growth factor (TGF)-β and macrophage interleukin (IL)-1β synergy in driving the transcriptional profile of FAPα+THY1+ invasive synovial fibroblasts, expanded in RA compared with PsA. In vitro characterisation of patient with RA synovial fibroblasts showed metabolic switch to glycolysis, increased adhesion intercellular adhesion molecules 1 expression and IL-6 secretion in response to combined TGF-β and IL-1β treatment. Disrupting specific immune and stromal cell interactions offers novel opportunities for targeted therapeutic intervention in RA and PsA.

Signaling Pathways Associated with Chronic Wound Progression: A Systems Biology Approach

Previously we have shown that several oxidative stress-driven pathways in cutaneous chronic wounds are dysregulated in the first 48 h post-wounding. Here, we performed an RNASeq analysis of tissues collected up to day 20 after wounding, when we have determined full chronicity is established. Weighted Gene Correlation Network Analysis was performed in R segregating the genes into 14 modules. Genes in the modules significantly correlated (p < 0.05) to early and full chronicity were used for pathway analysis using pathfindR. In early chronicity, we observed enrichment of several pathways. Dysregulation of Ephrin/Eph signaling leads to growth cone collapse and impairs neuronal regeneration. Adra2b and Adra2a overexpression in early and full chronicity, respectively, decreased cAMP production and impaired re-epithelialization and granulation tissue formation. Several pathways involving a Smooth-muscle-actin (Acta1) were also enriched with Acta1 overexpression contributing to impaired angiogenesis. During full chronicity, the ‘JAK-STAT’ pathway was suppressed undermining host defenses against infection. Wnt signaling was also suppressed, impairing re-epithelialization and granulation tissue formation. Biomarkers of cancer such as overexpression of SDC1 and constitutive activation of ErbB2/HER2 were also identified. In conclusion, we show that during progression to full chronicity, numerous signaling pathways are dysregulated, including some related to carcinogenesis, suggesting that chronic wounds behave much like cancer. Experimental verification in vivo could identify candidates for treatment of chronic wounds.

The DCMU Herbicide Shapes T-cell Functions By Modulating Micro-RNA Expression Profiles

DCMU [N-(3,4-dichlorophenyl)-N-dimethylurea] or diuron is a widely used herbicide, which can cause adverse effects on human, especially on immune cells, due to their intrinsic properties and wide distribution. These cells are important for fighting not only against virus or bacteria but also against neoplastic cell development. We developed an approach that combines functional studies and miRNA and RNA sequencing data to evaluate the effects of DCMU on the human immune response against cancer, particularly the one carried out by CD8+ T cells. We found that DCMU modulates the expression of miRNA in a dose-dependent manner, leading to a specific pattern of gene expression and consequently to a diminished cytokine and granzyme B secretions. Using mimics or anti-miRs, we identified several miRNA, such as hsa-miR-3135b and hsa-miR-21-5p, that regulate these secretions. All these changes reduce the CD8+ T cells’ cytotoxic activity directed against cancer cells, in vitro and in vivo in a zebrafish model. To conclude, our study suggests that DCMU reduces T-cell abilities, participating thus to the establishment of an environment conducive to cancer development.

Effect of Diet Enriched With Hemp Seeds on Goat Milk Fatty Acids, Transcriptome, and miRNAs

In dairy ruminants, a diet supplemented with feed rich in unsaturated fatty acids can be an effective medium to increase the health-promoting properties of milk, although their effect on the pathways/genes involved in these processes has not been properly and completely defined to date. To improve our knowledge of the cell’s activity in specific conditions, next-generation RNA-sequencing technology was used to allow whole transcriptome characterization under given conditions. In addition to this, microRNAs (miRNAs) have recently been known as post-transcriptional regulators in fatty acid and cholesterol metabolism by targeting lipid metabolism genes. In this study, to analyze the transcriptome and miRNAs in goat milk after a supplemental diet enriched with linoleic acid (hemp seeds), next-generation RNA-sequencing was used in order to point out the general biological mechanisms underlying the effects related to milk fat metabolism. Ten pluriparous Alpine goats were fed with the same pretreatment diet for 40 days; then, they were arranged to two dietary treatments consisting of control (C) and hemp seed (H)-supplemented diets. Milk samples were collected at 40 (time point = T0) and 140 days of lactation (time point = T1). Milk fatty acid (FA) profiles revealed a significant effect of hemp seeds that determined a strong increment in the preformed FA, causing a reduction in the concentration of de-novo FA. Monounsaturated and polyunsaturated n−3 FAs were increased by hemp treatment, determining a reduction in the n−6/n−3 ratio. After removing milk fats and proteins, RNA was extracted from the milk cells and transcriptomic analysis was conducted using Illumina RNA-sequencing. A total of 3,835 genes were highly differentially expressed (p-value &lt; 0.05, fold change &gt; 1.5, and FDR &lt; 0.05) in the H group. Functional analyses evidenced changes in metabolism, immune, and inflammatory responses. Furthermore, modifications in feeding strategies affected also key transcription factors regulating the expression of several genes involved in milk fat metabolism, such as peroxisome proliferator-activated receptors (PPARs). Moreover, 38 (15 known and 23 novel) differentially expressed miRNAs were uncovered in the H group and their potential functions were also predicted. This study gives the possibility to improve our knowledge of the molecular changes occurring after a hemp seed supplementation in the goat diet and increase our understanding of the relationship between nutrient variation and phenotypic effects.

Hepatic Expression of the Na+-Taurocholate Cotransporting Polypeptide Is Independent from Genetic Variation

The hepatic Na+-taurocholate cotransporting polypeptide NTCP/SLC10A1 is important for the uptake of bile salts and selected drugs. Its inhibition results in increased systemic bile salt concentrations. NTCP is also the entry receptor for the hepatitis B/D virus. We investigated interindividual hepatic SLC10A1/NTCP expression using various omics technologies. SLC10A1/NTCP mRNA expression/protein abundance was quantified in well-characterized 143 human livers by real-time PCR and LC-MS/MS-based targeted proteomics. Genome-wide SNP arrays and SLC10A1 next-generation sequencing were used for genomic analyses. SLC10A1 DNA methylation was assessed through MALDI-TOF MS. Transcriptomics and untargeted metabolomics (UHPLC-Q-TOF-MS) were correlated to identify NTCP-related metabolic pathways. SLC10A1 mRNA and NTCP protein levels varied 44-fold and 10.4-fold, respectively. Non-genetic factors (e.g., smoking, alcohol consumption) influenced significantly NTCP expression. Genetic variants in SLC10A1 or other genes do not explain expression variability which was validated in livers (n = 50) from The Cancer Genome Atlas. The identified two missense SLC10A1 variants did not impair transport function in transfectants. Specific CpG sites in SLC10A1 as well as single metabolic alterations and pathways (e.g., peroxisomal and bile acid synthesis) were significantly associated with expression. Inter-individual variability of NTCP expression is multifactorial with the contribution of clinical factors, DNA methylation, transcriptional regulation as well as hepatic metabolism, but not genetic variation.

Optimal CD8+ T cell effector function requires costimulation-induced RNA-binding proteins that reprogram the transcript isoform landscape

Boosting T cell activation through costimulation directs defense against cancer and viral infections. Despite multiple studies targeting costimulation in clinical trials, the increased potency and reprogramming of T cells endowed by costimulation is poorly understood. Canonical dogma states that transcription mediates T cell activation. Here, we show that the spliceosome, controlling post-transcriptional alternative splicing and alternative polyadenylation, is the most enriched pathway in T cells after CD134/CD137 costimulation. Costimulation of CD8+ T cells significantly increases expression of 29 RNA-binding proteins while RNA-seq uncovers over 1000 differential alternative splicing and polyadenylation events. Using in vivo mouse and in vitro human models, we demonstrate that RNA-binding protein Tardbp is required for effector cytokine production, CD8+ T cell clonal expansion, and isoform regulation after costimulation. The prospect of immune response optimization through reprogramming of mRNA isoform production offered herein opens new avenues for experimentally and therapeutically tuning the activities of T cells.

Identification of HOX signatures contributing to oral cancer phenotype

The role of evolutionarily conserved homeobox-containing HOX genes as transcriptional regulators in the developmental specification of organisms is well known. The contribution of HOX genes involvement in oral cancer phenotype has yet to be fully ascertained. TCGA-HNSC HTSeq-counts and clinical data were retrieved from the GDC portal for oral cavity neoplasms. GEO datasets (GSE72627, GSE30784, GSE37991) were accessed and analyzed using GEO2R. Differential HOX gene expression was profiled using the DESeq2 R package with a log2 fold change cut-off (- 1 and + 1) and Benjamini-Hochberg p-adjusted value at ≤ 0.01. Gene set over-representation analysis and semantic analysis associated with the disease ontology was performed using the ClusterProfiler R package, and pathway over-representation analysis was performed using IMPaLa. HOX protein interaction network was constructed using the Pathfind R package. HOX phenotype associations were performed using Mammalian Phenotype Ontology, Human Phenotype Ontology, PhenGenI associations, Jensen tissues, and OMIM entries. Drug connectivity mapping was carried out with Dr. Insight R package. HOXA2 was upregulated in oral dysplasia but silenced during tumor progression. Loss of HOXB2 expression was consistent in the potentially malignant oral lesions as well as in the primary tumor. HOXA7, HOXA10, HOXB7, HOXC6, HOXC10, HOXD10, and HOXD11 were consistently upregulated from premalignancy to malignancy and were notably associated with risk factors. Overrepresentation analysis suggested HOXA10 was involved in the transcriptional misregulation contributing to the oral cancer phenotype. HOX genes subnetwork analysis showed crucial interactions with cell cycle regulators, growth responsive elements, and proto-oncogenes. Phenotype associations specific to the oral region involving HOX genes provide intrinsic cues to tumor development. The 5' HOX genes were aberrantly upregulated during oral carcinogenesis reflecting their posterior prevalence.

Brain aging is faithfully modelled in organotypic brain slices and accelerated by prions

Mammalian models are essential for brain aging research. However, the long lifespan and poor amenability to genetic and pharmacological perturbations have hindered the use of mammals for dissecting aging-regulatory molecular networks and discovering new anti-aging interventions. To circumvent these limitations, we developed an ex vivo model system that faithfully mimics the aging process of the mammalian brain using cultured mouse brain slices. Genome-wide gene expression analyses showed that cultured brain slices spontaneously upregulated senescence-associated genes over time and reproduced many of the transcriptional characteristics of aged brains. Treatment with rapamycin, a classical anti-aging compound, largely abolished the time-dependent transcriptional changes in naturally aged brain slice cultures. Using this model system, we discovered that prions drastically accelerated the development of age-related molecular signatures and the pace of brain aging. We confirmed this finding in mouse models and human victims of Creutzfeldt-Jakob disease. These data establish an innovative, eminently tractable mammalian model of brain aging, and uncover a surprising acceleration of brain aging in prion diseases.

Integrated bioinformatics analysis of validated and circulating miRNAs in ovarian cancer

Recent studies have focused on the early detection of ovarian cancer (OC) using tumor materials by liquid biopsy. The mechanisms of microRNAs (miRNAs) to impact OC and signaling pathways are still unknown. This study aims to reliably perform functional analysis of previously validated circulating miRNAs' target genes by using pathfindR. Also, overall survival and pathological stage analyses were evaluated with miRNAs' target genes which are common in the The Cancer Genome Atlas and GTEx datasets. Our previous studies have validated three downregulated miRNAs (hsa-miR-885-5p, hsa-miR-1909-5p, and hsalet7d-3p) having a diagnostic value in OC patients' sera, with high-throughput techniques. The predicted target genes of these miRNAs were retrieved from the miRDB database (v6.0). Active-subnetwork-oriented Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis was conducted by pathfindR using the target genes. Enrichment of KEGG pathways assessed by the analysis of pathfindR indicated that 24 pathways were related to the target genes. Ubiquitin-mediated proteolysis, spliceosome and Notch signaling pathway were the top three pathways with the lowest p-values (p < 0.001). Ninety-three common genes were found to be differentially expressed (p < 0.05) in the datasets. No significant genes were found to be significant in the analysis of overall survival analyses, but 24 genes were found to be significant with pathological stages analysis (p < 0.05). The findings of our study provide in-silico evidence that validated circulating miRNAs' target genes and enriched pathways are related to OC and have potential roles in theranostics applications. Further experimental investigations are required to validate our results which will ultimately provide a new perspective for translational applications in OC management.

Constitutive gene expression differs in three brain regions important for cognition in neophobic and non-neophobic house sparrows (Passer domesticus)

Neophobia (aversion to new objects, food, and environments) is a personality trait that affects the ability of wildlife to adapt to new challenges and opportunities. Despite the ubiquity and importance of this trait, the molecular mechanisms underlying repeatable individual differences in neophobia in wild animals are poorly understood. We evaluated wild-caught house sparrows (Passer domesticus) for neophobia in the lab using novel object tests. We then selected a subset of neophobic and non-neophobic individuals (n = 3 of each, all females) and extracted RNA from four brain regions involved in learning, memory, threat perception, and executive function: striatum, caudal dorsomedial hippocampus, medial ventral arcopallium, and caudolateral nidopallium (NCL). Our analysis of differentially expressed genes (DEGs) used 11,889 gene regions annotated in the house sparrow reference genome for which we had an average of 25.7 million mapped reads/sample. PERMANOVA identified significant effects of brain region, phenotype (neophobic vs. non-neophobic), and a brain region by phenotype interaction. Comparing neophobic and non-neophobic birds revealed constitutive differences in DEGs in three of the four brain regions examined: hippocampus (12% of the transcriptome significantly differentially expressed), striatum (4%) and NCL (3%). DEGs included important known neuroendocrine mediators of learning, memory, executive function, and anxiety behavior, including serotonin receptor 5A, dopamine receptors 1, 2 and 5 (downregulated in neophobic birds), and estrogen receptor beta (upregulated in neophobic birds). These results suggest that some of the behavioral differences between phenotypes may be due to underlying gene expression differences in the brain. The large number of DEGs in neophobic and non-neophobic birds also implies that there are major differences in neural function between the two phenotypes that could affect a wide variety of behavioral traits beyond neophobia.

Benefits and Challenges of Pre-clustered Network-Based Pathway Analysis

Functional analysis of gene sets derived from experiments is typically done by pathway annotation. Although many algorithms exist for analyzing the association between a gene set and a pathway, an issue which is generally ignored is that gene sets often represent multiple pathways. In such cases an association to a pathway is weakened by the presence of genes associated with other pathways. A way to counteract this is to cluster the gene set into more homogenous parts before performing pathway analysis on each module. We explored whether network-based pre-clustering of a query gene set can improve pathway analysis. The methods MCL, Infomap, and MGclus were used to cluster the gene set projected onto the FunCoup network. We characterized how well these methods are able to detect individual pathways in multi-pathway gene sets, and applied each of the clustering methods in combination with four pathway analysis methods: Gene Enrichment Analysis, BinoX, NEAT, and ANUBIX. Using benchmarks constructed from the KEGG pathway database we found that clustering can be beneficial by increasing the sensitivity of pathway analysis methods and by providing deeper insights of biological mechanisms related to the phenotype under study. However, keeping a high specificity is a challenge. For ANUBIX, clustering caused a minor loss of specificity, while for BinoX and NEAT it caused an unacceptable loss of specificity. GEA had very low sensitivity both before and after clustering. The choice of clustering method only had a minor effect on the results. We show examples of this approach and conclude that clustering can improve overall pathway annotation performance, but should only be used if the used enrichment method has a low false positive rate.

Unspecific CTL Killing Is Enhanced by High Glucose via TNF-Related Apoptosis-Inducing Ligand

TNF-related apoptosis inducing ligand (TRAIL) is expressed on cytotoxic T lymphocytes (CTLs) and TRAIL is linked to progression of diabetes. However, the impact of high glucose on TRAIL expression and its related killing function in CTLs still remains largely elusive. Here, we report that TRAIL is substantially up-regulated in CTLs in environments with high glucose (HG) both in vitro and in vivo. Non-mitochondrial reactive oxygen species, NFκB and PI3K/Akt are essential in HG-induced TRAIL upregulation in CTLs. TRAIL^high CTLs induce apoptosis of pancreatic beta cell line 1.4E7. Treatment with metformin and vitamin D reduces HG-enhanced expression of TRAIL in CTLs and coherently protects 1.4E7 cells from TRAIL-mediated apoptosis. Our work suggests that HG-induced TRAIL^high CTLs might contribute to the destruction of pancreatic beta cells in a hyperglycemia condition.

Antibody-mediated blockade of the IL23 receptor destabilizes intratumoral regulatory T cells and enhances immunotherapy

Regulatory T cells (Treg) can impede antitumor immunity and currently represent a major obstacle to effective cancer immunotherapy. Targeting tumor-infiltrating regulatory Treg while sparing systemic Treg represents an optimal approach to this problem. Here, we provide evidence that the interleukin 23 receptor (IL23R) expressed by tumor-infiltrating Treg promotes suppressive activity. Disruption of the IL23R results in increased responsiveness of destabilized Treg to the IL12 cytokine, the production of γ-interferon, and the recruitment of CD8 T cells that inhibit tumor growth. Since the Treg destabilization pathway that is initiated by IL23R blockade is distinct and independent from the destabilization pathway coupled to glucocorticoid-induced TNFR-related protein (GITR) activation, we examined the impact of the coordinate induction of the two destabilization pathways on antitumor immune responses. Combined GITR and IL23R antibody treatment of mice inoculated with MC38 tumors resulted in robust and synergistic antitumor responses. These findings indicate that the delineation of independent Treg destabilization pathways may allow improved approaches to the development of combination immunotherapy for cancers.

Proteomic profiles in cerebrospinal fluid predicted death and disability in term infants with perinatal asphyxia: A pilot study

Aim

Perinatal asphyxia, resulting in hypoxic‐ischaemic encephalopathy (HIE), has been associated with high mortality rates and severe lifelong neurodevelopmental disabilities. Our aim was to study the association between the proteomic profile in cerebrospinal fluid (CSF) and the degree of HIE and long‐term outcomes.

Methods

We prospectively enrolled 18‐term born infants with HIE and 10‐term born controls between 2000 and 2004 from the Karolinska University Hospital. An antibody suspension bead array and FlexMap3D analysis was used to characterise 178 unique brain‐derived and inflammation associated proteins in their CSF.

Results

Increased CSF concentrations of several brain‐specific proteins were observed in the proteome of HIE patients compared with the controls. An upregulation of neuroinflammatory pathways was also noted and this was confirmed by pathway analysis. Principal component analysis revealed a gradient from favourable to unfavourable HIE grades and outcomes. The proteins that provided strong predictors were structural proteins, including myelin basic protein and alpha‐II spectrin. The functional proteins included energy‐related proteins like neuron‐specific enolase and synaptic regulatory proteins. Increased CSF levels of 51 proteins correlated with adverse outcomes in infants with HIE.

Conclusion

Brain‐specific proteins and neuroinflammatory mediators in CSF may predict HIE degrees and outcomes after perinatal asphyxia.

An RNA-sequencing transcriptome of the rodent Schwann cell response to peripheral nerve injury

BACKGROUND

The important contribution of glia to mechanisms of injury and repair of the nervous system is increasingly recognized. In stark contrast to the central nervous system (CNS), the peripheral nervous system (PNS) has a remarkable capacity for regeneration after injury. Schwann cells are recognized as key contributors to PNS regeneration, but the molecular underpinnings of the Schwann cell response to injury and how they interact with the inflammatory response remain incompletely understood.

METHODS

We completed bulk RNA-sequencing of Schwann cells purified acutely using immunopanning from the naïve and injured rodent sciatic nerve at 3, 5, and 7 days post-injury. We used qRT-PCR and in situ hybridization to assess cell purity and probe dataset integrity. Finally, we used bioinformatic analysis to probe Schwann cell-specific injury-induced modulation of cellular pathways.

RESULTS

Our data confirm Schwann cell purity and validate RNAseq dataset integrity. Bioinformatic analysis identifies discrete modules of genes that follow distinct patterns of regulation in the 1st days after injury and their corresponding molecular pathways. These findings enable improved differentiation of myeloid and glial components of neuroinflammation after peripheral nerve injury and highlight novel molecular aspects of the Schwann cell injury response such as acute downregulation of the AGE/RAGE pathway and of secreted molecules Sparcl1 and Sema5a.

CONCLUSIONS

We provide a helpful resource for further deciphering the Schwann cell injury response and a depth of transcriptional data that can complement the findings of recent single cell sequencing approaches. As more data become available on the response of CNS glia to injury, we anticipate that this dataset will provide a valuable platform for understanding key differences in the PNS and CNS glial responses to injury and for designing approaches to ameliorate CNS regeneration.

In silico identification of novel biomarkers for key players in transition from normal colon tissue to adenomatous polyps

Adenomatous polyps of the colon are the most common neoplastic polyps. Although most of adenomatous polyps do not show malign transformation, majority of colorectal carcinomas originate from neoplastic polyps. Therefore, understanding of this transformation process would help in both preventive therapies and evaluation of malignancy risks. This study uncovers alterations in gene expressions as potential biomarkers that are revealed by integration of several network-based approaches. In silico analysis performed on a unified microarray cohort, which is covering 150 normal colon and adenomatous polyp samples. Significant gene modules were obtained by a weighted gene co-expression network analysis. Gene modules with similar profiles were mapped to a colon tissue specific functional interaction network. Several clustering algorithms run on the colon-specific network and the most significant sub-modules between the clusters were identified. The biomarkers were selected by filtering differentially expressed genes which also involve in significant biological processes and pathways. Biomarkers were also validated on two independent datasets based on their differential gene expressions. To the best of our knowledge, such a cascaded network analysis pipeline was implemented for the first time on a large collection of normal colon and polyp samples. We identified significant increases in TLR4 and MSX1 expressions as well as decrease in chemokine profiles with mostly pro-tumoral activities. These biomarkers might appear as both preventive targets and biomarkers for risk evaluation. As a result, this research proposes novel molecular markers that might be alternative to endoscopic approaches for diagnosis of adenomatous polyps.

Placental Abnormalities are Associated With Specific Windows of Embryo Culture in a Mouse Model

Assisted Reproductive Technologies (ART) employ gamete/embryo handling and culture in vitro to produce offspring. ART pregnancies have an increased risk of low birth weight, abnormal placentation, pregnancy complications, and imprinting disorders. Embryo culture induces low birth weight, abnormal placental morphology, and lower levels of DNA methylation in placentas in a mouse model of ART. Whether preimplantation embryos at specific stages of development are more susceptible to these perturbations remains unresolved. Accordingly, we performed embryo culture for several discrete periods of preimplantation development and following embryo transfer, assessed fetal and placental outcomes at term. We observed a reduction in fetal:placental ratio associated with two distinct windows of preimplantation embryo development, one prior to the morula stage and the other from the morula to blastocyst stage, whereas placental morphological abnormalities and reduced imprinting control region methylation were only associated with culture prior to the morula stage. Extended culture to the blastocyst stage also induces additional placental DNA methylation changes compared to embryos transferred at the morula stage, and female concepti exhibited a higher loss of DNA methylation than males. By identifying specific developmental windows of susceptibility, this study provides a framework to optimize further culture conditions to minimize risks associated with ART pregnancies.

Avian auditory hair cell regeneration is accompanied by JAK/STAT-dependent expression of immune-related genes in supporting cells

The avian hearing organ is the basilar papilla that, in sharp contrast to the mammalian cochlea, can regenerate sensory hair cells and thereby recover from deafness within weeks. The mechanisms that trigger, sustain and terminate the regenerative response in vivo are largely unknown. Here, we profile the changes in gene expression in the chicken basilar papilla after aminoglycoside antibiotic-induced hair cell loss using RNA-sequencing. We identified changes in gene expression of a group of immune-related genes and confirmed with single-cell RNA-sequencing that these changes occur in supporting cells. In situ hybridization was used to further validate these findings. We determined that the JAK/STAT signaling pathway is essential for upregulation of the damage-response genes in supporting cells during the second day after induction of hair cell loss. Four days after ototoxic damage, we identified newly regenerated, nascent auditory hair cells that express genes linked to termination of the JAK/STAT signaling response. The robust, transient expression of immune-related genes in supporting cells suggests a potential functional involvement of JAK/STAT signaling in sensory hair cell regeneration.

Biometal Dyshomeostasis in Olfactory Mucosa of Alzheimer's Disease Patients

Olfactory function, orchestrated by the cells of the olfactory mucosa at the rooftop of the nasal cavity, is disturbed early in the pathogenesis of Alzheimer's disease (AD). Biometals including zinc and calcium are known to be important for sense of smell and to be altered in the brains of AD patients. Little is known about elemental homeostasis in the AD patient olfactory mucosa. Here we aimed to assess whether the disease-related alterations to biometal homeostasis observed in the brain are also reflected in the olfactory mucosa. We applied RNA sequencing to discover gene expression changes related to metals in olfactory mucosal cells of cognitively healthy controls, individuals with mild cognitive impairment and AD patients, and performed analysis of the elemental content to determine metal levels. Results demonstrate that the levels of zinc, calcium and sodium are increased in the AD olfactory mucosa concomitantly with alterations to 17 genes related to metal-ion binding or metal-related function of the protein product. A significant elevation in alpha-2-macroglobulin, a known metal-binding biomarker correlated with brain disease burden, was observed on the gene and protein levels in the olfactory mucosa cells of AD patients. These data demonstrate that the olfactory mucosa cells derived from AD patients recapitulate certain impairments of biometal homeostasis observed in the brains of patients.

Synovial Fibroblast Sialylation Regulates Cell Migration and Activation of Inflammatory Pathways in Arthritogenesis

Synovial fibroblasts have emerged as critical underlying factors to perpetuate chronic joint inflammation in Rheumatoid Arthritis. Like any other cell, synovial fibroblasts are covered with a complex layer of glycans that can change in response to extracellular signals, such as inflammation. We have previously shown that inflammatory synovial fibroblasts show decreased levels of sialic acid, but our understanding of sialic acid-dependent pathophysiological pathways in these stromal cells is still very limited. In this report, we used in vivo and in vitro studies with exogenous sialidases and RNA sequencing to investigate the responses of murine synovial fibroblasts upon desialylation. Our results show that hyposialylated fibroblasts present a dysregulated migratory ability and an activated phenotype characterized by the expression of inflammatory mediators, such as cytokines and chemokines, and anti-viral related mechanisms. Removal of surface sialic acid also affected the expression of sialyltransferases, revealing the existence of a positive feedback to sustain reduced sialylation. Moreover, we demonstrate that synovial fibroblasts subsets have distinct sialyltransferase expression profiles, both in healthy and arthritic mice. These findings underline the ability of sialic acid to modulate homeostatic and inflammatory responses in non-immune synovial fibroblasts, suggesting that sialylation plays a key role in perpetuating local inflammation in the arthritic joint.

Inside the Joint of Inflammatory Arthritis Patients: Handling and Processing of Synovial Tissue Biopsies for High Throughput Analysis

Inflammatory arthritis is a chronic systemic autoimmune disease of unknown etiology, which affects the joints. If untreated, these diseases can have a detrimental effect on the patient's quality of life, leading to disabilities, and therefore, exhibit a significant socioeconomic impact and burden. While studies of immune cell populations in arthritis patient's peripheral blood have been informative regarding potential immune cell dysfunction and possible patient stratification, there are considerable limitations in identifying the early events that lead to synovial inflammation. The joint, as the site of inflammation and the local microenvironment, exhibit unique characteristics that contribute to disease pathogenesis. Understanding the contribution of immune and stromal cell interactions within the inflamed joint has been met with several technical challenges. Additionally, the limited availability of synovial tissue biopsies is a key incentive for the utilization of high-throughput techniques in order to maximize information gain. This review aims to provide an overview of key methods and novel techniques that are used in the handling, processing and analysis of synovial tissue biopsies and the potential synergy between these techniques. Herein, we describe the utilization of high dimensionality flow cytometric analysis, single cell RNA sequencing, ex vivo functional assays and non-intrusive metabolic characterization of synovial cells on a single cell level based on fluorescent lifetime imaging microscopy. Additionally, we recommend important points of consideration regarding the effect of different storage and handling techniques on downstream analysis of synovial tissue samples. The introduction of new powerful techniques in the study of synovial tissue inflammation, brings new challenges but importantly, significant opportunities. Implementation of novel approaches will accelerate our path toward understanding of the mechanisms involved in the pathogenesis of inflammatory arthritis and lead to the identification of new avenues of therapeutic intervention.

Computing microRNA-gene interaction networks in pan-cancer using miRDriver

DNA copy number aberrated regions in cancer are known to harbor cancer driver genes and the short non-coding RNA molecules, i.e., microRNAs. In this study, we integrated the multi-omics datasets such as copy number aberration, DNA methylation, gene and microRNA expression to identify the signature microRNA-gene associations from frequently aberrated DNA regions across pan-cancer utilizing a LASSO-based regression approach. We studied 7294 patient samples associated with eighteen different cancer types from The Cancer Genome Atlas (TCGA) database and identified several cancer-specific and common microRNA-gene interactions enriched in experimentally validated microRNA-target interactions. We highlighted several oncogenic and tumor suppressor microRNAs that were cancer-specific and common in several cancer types. Our method substantially outperformed the five state-of-art methods in selecting significantly known microRNA-gene interactions in multiple cancer types. Several microRNAs and genes were found to be associated with tumor survival and progression. Selected target genes were found to be significantly enriched in cancer-related pathways, cancer hallmark and Gene Ontology (GO) terms. Furthermore, subtype-specific potential gene signatures were discovered in multiple cancer types.

EGCG Prevents the Onset of an Inflammatory and Cancer-Associated Adipocyte-like Phenotype in Adipose-Derived Mesenchymal Stem/Stromal Cells in Response to the Triple-Negative Breast Cancer Secretome

Background: Triple-negative breast cancer (TNBC) cells secretome induces a pro-inflammatory microenvironment within the adipose tissue, which hosts both mature adipocytes and adipose-derived mesenchymal stem/stromal cells (ADMSC). The subsequent acquisition of a cancer-associated adipocyte (CAA)-like phenotype is, however, unknown in ADMSC. While epidemiological studies suggest that consuming a polyphenol-rich diet reduces the incidence of some obesity-related cancers, the chemopreventive impact of green tea-derived epigallocatechin-3-gallate (EGCG) against the cues that trigger the CAA phenotype remain undocumented in ADMSC. Methods: Human ADMSC were exposed to human TNBC-derived MDA-MB-231 conditioned media (TNBC cells secretome) supplemented or not with EGCG. Differential gene expression was assessed through RNA-Seq analysis and confirmed by RT-qPCR. Protein expression levels and the activation status of signal transduction pathways mediators were determined by Western blotting. ADMSC chemotaxis was assessed by a real-time cell migration assay. Results: The TNBC cells secretome induced in ADMSC the expression of the CAA cytokines CCL2, CCL5, IL-1β, and IL-6, and of immunomodulators COX2, HIF-1α, VEGFα, and PD-L1. The epithelial-to-mesenchymal biomarker Snail was found to control the CAA phenotype. EGCG inhibited the induction of CAA genes and the activation status of Smad2 and NF-κB. The induced chemotactic response was also inhibited by EGCG. Conclusion: The induction of an inflammatory and CAA-like phenotype in ADMSC can be triggered by the TNBC cells secretome, while still efficiently prevented by diet-derived polyphenols.

Disruption of STIM1-mediated Ca2+ sensing and energy metabolism in adult skeletal muscle compromises exercise tolerance, proteostasis, and lean mass

OBJECTIVE

Stromal interaction molecule 1 (STIM1) is a single-pass transmembrane endoplasmic/sarcoplasmic reticulum (E/SR) protein recognized for its role in a store operated Ca2+ entry (SOCE), an ancient and ubiquitous signaling pathway. Whereas STIM1 is known to be indispensable during development, its biological and metabolic functions in mature muscles remain unclear.

METHODS

Conditional and tamoxifen inducible muscle STIM1 knock-out mouse models were coupled with multi-omics tools and comprehensive physiology to understand the role of STIM1 in regulating SOCE, mitochondrial quality and bioenergetics, and whole-body energy homeostasis.

RESULTS

This study shows that STIM1 is abundant in adult skeletal muscle, upregulated by exercise, and is present at SR-mitochondria interfaces. Inducible tissue-specific deletion of STIM1 (iSTIM1 KO) in adult muscle led to diminished lean mass, reduced exercise capacity, and perturbed fuel selection in the settings of energetic stress, without affecting whole-body glucose tolerance. Proteomics and phospho-proteomics analyses of iSTIM1 KO muscles revealed molecular signatures of low-grade E/SR stress and broad activation of processes and signaling networks involved in proteostasis.

CONCLUSION

These results show that STIM1 regulates cellular and mitochondrial Ca2+ dynamics, energy metabolism and proteostasis in adult skeletal muscles. Furthermore, these findings provide insight into the pathophysiology of muscle diseases linked to disturbances in STIM1-dependent Ca2+ handling.

Sexual dimorphism in neurological function after SCI is associated with disrupted neuroinflammation in both injured spinal cord and brain

Whereas human spinal cord injury (SCI) is more common in men, the prevalence is growing in women. However, little is known about the effect of biological sex on brain dysfunction and injury mechanisms. To model the highest per capita rate of injury (ages between 16 and 30 years old) in humans, in the present study, young adult or a young/middle-aged male and female C57BL/6 mice were subjected to moderate contusion SCI. When mice were injured at 10-12-week-old, transcriptomic analysis of inflammation-related genes and flow cytometry revealed a more aggressive neuroinflammatory profile in male than females following 3 d SCI, ostensibly driven by sex-specific changes myeloid cell function rather than cell number. Female mice were generally more active at baseline, as evidenced by greater distance traveled in the open field. After SCI, female mice had more favorable locomotor function than male animals. At 13 weeks post-injury, male mice showed poor performance in cognitive and depressive-like behavioral tests, while injured female mice showed fewer deficits in these tasks. However, when injured at 6 months old followed by 8 months post-injury, male mice had considerably less inflammatory activation compared with female animals despite having similar or worse outcomes in affective, cognitive, and motor tasks. Collectively, these findings indicate that sex differences in functional outcome after SCI are associated with the age at onset of injury, as well as disrupted neuroinflammation not only at the site of injury but also in remote brain regions. Thus, biological sex should be considered when designing new therapeutic agents.

Disulfiram inhibits neutrophil extracellular trap formation protecting rodents from acute lung injury and SARS-CoV-2 infection

Severe acute lung injury has few treatment options and a high mortality rate. Upon injury, neutrophils infiltrate the lungs and form neutrophil extracellular traps (NETs), damaging the lungs and driving an exacerbated immune response. Unfortunately, no drug preventing NET formation has completed clinical development. Here, we report that disulfiram — an FDA-approved drug for alcohol use disorder — dramatically reduced NETs, increased survival, improved blood oxygenation, and reduced lung edema in a transfusion-related acute lung injury (TRALI) mouse model. We then tested whether disulfiram could confer protection in the context of SARS-CoV-2 infection, as NETs are elevated in patients with severe COVID-19. In SARS-CoV-2–infected golden hamsters, disulfiram reduced NETs and perivascular fibrosis in the lungs, and it downregulated innate immune and complement/coagulation pathways, suggesting that it could be beneficial for patients with COVID-19. In conclusion, an existing FDA-approved drug can block NET formation and improve disease course in 2 rodent models of lung injury for which treatment options are limited.

Muscle transcriptomics shows overexpression of cadherin 1 in inclusion body myositis

Objective

This study aimed to elucidate the molecular features of inclusion body myositis (IBM).

Methods

We performed RNA sequencing analysis of muscle biopsy samples from 67 participants, consisting of 58 myositis patients with the pathological finding of CD8‐positive T cells invading non‐necrotic muscle fibers expressing major histocompatibility complex class I (43 IBM, 6 polymyositis, and 9 unclassifiable myositis), and 9 controls.

Results

Cluster analysis, principal component analysis, and pathway analysis showed that differentially expressed genes and pathways identified in IBM and polymyositis were mostly comparable. However, pathways related to cell adhesion molecules were upregulated in IBM as compared with polymyositis and controls (p < 0.01). Notably, CDH1, which encodes the epidermal cell junction protein cadherin 1, was overexpressed in the muscles of IBM, which was validated by another RNA sequencing dataset from previous publications. Western blotting confirmed the presence of mature cadherin 1 protein in the muscles of IBM. Immunohistochemical staining confirmed the positivity for anti‐cadherin 1 antibody in the muscles of IBM, whereas there was no muscle fiber positive for anti‐cadherin 1 antibody in immune‐mediated necrotizing myopathy, antisynthetase syndrome, and controls. The fibers stained with anti‐cadherin 1 antibody did not have rimmed vacuoles or abnormal protein accumulation. Experimental skeletal muscle regeneration and differentiation systems showed that CDH1 is expressed during skeletal muscle regeneration and differentiation.

Interpretation

CDH1 was detected as a differentially expressed gene, and immunohistochemistry showed that cadherin 1 exists in the muscles of IBM, whereas it was rarely seen in those of other idiopathic inflammatory myopathies. Cadherin 1 upregulation in muscle could provide a valuable clue to the pathological mechanisms of IBM. ANN NEUROL 2022;91:317–328

Analysis of the effects of aryl hydrocarbon receptor expression on cancer cell invasion via three-dimensional microfluidic invasion assays

The aryl hydrocarbon receptor (AHR) is a ligand-activated transcription factor that binds to xenobiotics and activates expression of response elements to metabolize these compounds. The AHR pathway has been associated with a long list of diseases including cancer; however, it is debated whether AHR is tumorigenic or tumour-inhibiting. In particular, there are contradictory reports in the literature regarding the effects of AHR expression level on metastatic breast cancer. Here we used a 3D invasion assay called cell invasion in digital microfluidic microgel systems (CIMMS) to study the effect of AHR expression on invasion. In this study, MDA-MB-231 cells with stable knockout of AHR (AHRko) showed enhanced invasive characteristics and reduced proliferation, and cells with transient overexpression of AHR showed reduced invasiveness. Overexpression of AHR with a mutation in the DNA binding domain showed no difference in invasiveness compared to control, which suggests that the changes in invasiveness are related to the expression of AHR. CIMMS also allowed for extraction of sub-populations of invaded cells for RNA sequencing experiments. A comparison of the transcriptomes of invaded subpopulations of wild-type and AHRko cells identified 1809 genes that were differentially expressed, with enriched pathways including cell cycle, proliferation, survival, immunoproteasome activation, and activation of matrix metalloproteases. In sum, the data reported here for MDA-MB-231 cells suggests some new interpretations of the discrepancy in the literature on the role of AHR in breast cancer. We propose that the unique combination of functional discrimination with transcriptome profiling provided by CIMMS will be valuable for a wide range of mechanistic invasion-biology studies in the future.

Whole Transcriptome Profiling of Adrenocortical Tumors Using Formalin-Fixed Paraffin-Embedded Samples

Whole transcriptome profiling is a promising technique in adrenal studies; however, whole transcriptome profiling of adrenal disease using formalin-fixed paraffin-embedded (FFPE) samples has to be further explored. The aim of this study was to evaluate the utility of transcriptome data from FFPE samples of adrenocortical tumors. We performed whole transcriptome profiling of FFPE and fresh frozen samples of adrenocortical carcinoma (ACC, n = 3), aldosterone-producing adenoma (APA, n = 3), and cortisol-producing adenoma (CPA, n = 3), and examined the similarity between the transcriptome data. We further examined whether the transcriptome data of FFPE samples could be used to distinguish tumor types and detect marker genes. The number of read counts was smaller in FFPE samples than in fresh frozen samples (P < 0.01), while the number of genes detected was similar (P = 0.39). The gene expression profiles of FFPE and fresh frozen samples were highly correlated (r = 0.93, P < 0.01). Tumor types could be distinguished by consensus clustering and principal component analysis using transcriptome data from FFPE samples. In the differential expression analysis between ACC and APA-CPA, known marker genes of ACC (e.g., CCNB2, TOP2A, and MAD2L1) were detected in FFPE samples of ACC. In the differential expression analysis between APA and CPA, known marker genes of APA (e.g., CYP11B2, VSNL1, and KCNJ5) were detected in the APA of FFPE samples. The results suggest that FFPE samples may be a reliable alternative to fresh frozen samples for whole transcriptome profiling of adrenocortical tumors.

An Integrative Multi-Omics Analysis of The Molecular Links between Aging and Aggressiveness in Thyroid Cancers

Aging modifies risk in all cancers, but age is used as a clinical staging criterion uniquely in thyroid cancer (TC). The molecular drivers of age-dependent TC onset and aggressiveness remain poorly understood. We applied an integrative, multi-omics data analysis approach to characterize these signatures. Our analysis reveals that aging, independent of BRAF^V600E mutational status, drives a significant accumulation of aggressiveness-related markers and poorer survival outcomes, most noticeably at age 55 and over. We identified that chromosomal alterations in loci 1p/1q as aging-associated drivers of aggressiveness, and that depleted infiltration with tumor surveillant CD8+T and follicular helper T cells, dysregulation of proteostasis- and senescence-related processes, and ERK1/2 signaling cascade are key features of the aging thyroid and TC onset/progression and aggressiveness in aging patients but not in young individuals. A panel of 23 genes, including those related to cell division such as CENPF, ERCC6L, and the kinases MELK and NEK2, were identified and rigorously characterized as aging-dependent and aggressiveness-specific markers. These genes effectively stratified patients into aggressive clusters with distinct phenotypic enrichment and genomic/transcriptomic profiles. This panel also showed excellent performance in predicting metastasis stage, BRAF^V600E, TERT promoter mutation, and survival outcomes and was superior to the American Thyroid Association (ATA) methodology in predicting aggressiveness risk. Our analysis established clinically relevant biomarkers for TC aggressiveness factoring in aging as an important component.

GeneTonic: an R/Bioconductor package for streamlining the interpretation of RNA-seq data

BACKGROUND

The interpretation of results from transcriptome profiling experiments via RNA sequencing (RNA-seq) can be a complex task, where the essential information is distributed among different tabular and list formats-normalized expression values, results from differential expression analysis, and results from functional enrichment analyses. A number of tools and databases are widely used for the purpose of identification of relevant functional patterns, yet often their contextualization within the data and results at hand is not straightforward, especially if these analytic components are not combined together efficiently.

RESULTS

We developed the GeneTonic software package, which serves as a comprehensive toolkit for streamlining the interpretation of functional enrichment analyses, by fully leveraging the information of expression values in a differential expression context. GeneTonic is implemented in R and Shiny, leveraging packages that enable HTML-based interactive visualizations for executing drilldown tasks seamlessly, viewing the data at a level of increased detail. GeneTonic is integrated with the core classes of existing Bioconductor workflows, and can accept the output of many widely used tools for pathway analysis, making this approach applicable to a wide range of use cases. Users can effectively navigate interlinked components (otherwise available as flat text or spreadsheet tables), bookmark features of interest during the exploration sessions, and obtain at the end a tailored HTML report, thus combining the benefits of both interactivity and reproducibility.

CONCLUSION

GeneTonic is distributed as an R package in the Bioconductor project ( https://bioconductor.org/packages/GeneTonic/ ) under the MIT license. Offering both bird's-eye views of the components of transcriptome data analysis and the detailed inspection of single genes, individual signatures, and their relationships, GeneTonic aims at simplifying the process of interpretation of complex and compelling RNA-seq datasets for many researchers with different expertise profiles.

Glutamine-fructose-6-phosphate transaminase 2 (GFPT2) is upregulated in breast epithelial-mesenchymal transition and responds to oxidative stress

Breast cancer cells that have undergone partial epithelial–mesenchymal transition (EMT) are believed to be more invasive than cells that have completed EMT. To study metabolic reprogramming in different mesenchymal states, we analyzed protein expression following EMT in the breast epithelial cell model D492 with single-shot LFQ supported by a SILAC proteomics approach. The D492 EMT cell model contains three cell lines: the epithelial D492 cells, the mesenchymal D492M cells, and a partial mesenchymal, tumorigenic variant of D492 that overexpresses the oncogene HER2. The analysis classified the D492 and D492M cells as basal-like and D492HER2 as claudin-low. Comparative analysis of D492 and D492M to tumorigenic D492HER2 differentiated metabolic markers of migration from those of invasion. Glutamine-fructose-6-phosphate transaminase 2 (GFPT2) was one of the top dysregulated enzymes in D492HER2. Gene expression analysis of the cancer genome atlas showed that GFPT2 expression was a characteristic of claudin-low breast cancer. siRNA-mediated knockdown of GFPT2 influenced the EMT marker vimentin and both cell growth and invasion in vitro and was accompanied by lowered metabolic flux through the hexosamine biosynthesis pathway (HBP). Knockdown of GFPT2 decreased cystathionine and sulfide:quinone oxidoreductase (SQOR) in the transsulfuration pathway that regulates H₂S production and mitochondrial homeostasis. Moreover, GFPT2 was within the regulation network of insulin and EGF, and its expression was regulated by reduced glutathione (GSH) and suppressed by the oxidative stress regulator GSK3-β. Our results demonstrate that GFPT2 controls growth and invasion in the D492 EMT model, is a marker for oxidative stress, and associated with poor prognosis in claudin-low breast cancer.

•

GFPT2 is upregulated following EMT.

•

GFPT2 is a marker for claudin-low breast cancer.

•

GFPT2 affects vimentin, cell proliferation, and cell invasion.

•

GFPT2 responds to oxidative stress.

•

GFPT2 is regulated by insulin and EGF.

Transcriptomics of single dose and repeated carbon black and ozone inhalation co-exposure highlight progressive pulmonary mitochondrial dysfunction

BACKGROUND

Air pollution is a complex mixture of particles and gases, yet current regulations are based on single toxicant levels failing to consider potential interactive outcomes of co-exposures. We examined transcriptomic changes after inhalation co-exposure to a particulate and a gaseous component of air pollution and hypothesized that co-exposure would induce significantly greater impairments to mitochondrial bioenergetics. A whole-body inhalation exposure to ultrafine carbon black (CB), and ozone (O3) was performed, and the impact of single and multiple exposures was studied at relevant deposition levels. C57BL/6 mice were exposed to CB (10 mg/m3) and/or O3 (2 ppm) for 3 h (either a single exposure or four independent exposures). RNA was isolated from lungs and mRNA sequencing performed using the Illumina HiSeq. Lung pathology was evaluated by histology and immunohistochemistry. Electron transport chain (ETC) activities, electron flow, hydrogen peroxide production, and ATP content were assessed.

RESULTS

Compared to individual exposure groups, co-exposure induced significantly greater neutrophils and protein levels in broncho-alveolar lavage fluid as well as a significant increase in mRNA expression of oxidative stress and inflammation related genes. Similarly, a significant increase in hydrogen peroxide production was observed after co-exposure. After single and four exposures, co-exposure revealed a greater number of differentially expressed genes (2251 and 4072, respectively). Of these genes, 1188 (single exposure) and 2061 (four exposures) were uniquely differentially expressed, with 35 mitochondrial ETC mRNA transcripts significantly impacted after four exposures. Both O3 and co-exposure treatment significantly reduced ETC maximal activity for complexes I (- 39.3% and - 36.2%, respectively) and IV (- 55.1% and - 57.1%, respectively). Only co-exposure reduced ATP Synthase activity (- 35.7%) and total ATP content (30%). Further, the ability for ATP Synthase to function is limited by reduced electron flow (- 25%) and translation of subunits, such as ATP5F1, following co-exposure.

CONCLUSIONS

CB and O3 co-exposure cause unique transcriptomic changes in the lungs that are characterized by functional deficits to mitochondrial bioenergetics. Alterations to ATP Synthase function and mitochondrial electron flow underly a pathological adaptation to lung injury induced by co-exposure.

XDeathDB: a visualization platform for cell death molecular interactions

Lots of cell death initiator and effector molecules, signalling pathways and subcellular sites have been identified as key mediators in both cell death processes in cancer. The XDeathDB visualization platform provides a comprehensive cell death and their crosstalk resource for deciphering the signaling network organization of interactions among different cell death modes associated with 1461 cancer types and COVID-19, with an aim to understand the molecular mechanisms of physiological cell death in disease and facilitate systems-oriented novel drug discovery in inducing cell deaths properly. Apoptosis, autosis, efferocytosis, ferroptosis, immunogenic cell death, intrinsic apoptosis, lysosomal cell death, mitotic cell death, mitochondrial permeability transition, necroptosis, parthanatos, and pyroptosis related to 12 cell deaths and their crosstalk can be observed systematically by the platform. Big data for cell death gene-disease associations, gene-cell death pathway associations, pathway-cell death mode associations, and cell death-cell death associations is collected by literature review articles and public database from iRefIndex, STRING, BioGRID, Reactom, Pathway's commons, DisGeNET, DrugBank, and Therapeutic Target Database (TTD). An interactive webtool, XDeathDB, is built by web applications with R-Shiny, JavaScript (JS) and Shiny Server Iso. With this platform, users can search specific interactions from vast interdependent networks that occur in the realm of cell death. A multilayer spectral graph clustering method that performs convex layer aggregation to identify crosstalk function among cell death modes for a specific cancer. 147 hallmark genes of cell death could be observed in detail in these networks. These potential druggable targets are displayed systematically and tailoring networks to visualize specified relations is available to fulfil user-specific needs. Users can access XDeathDB for free at https://pcm2019.shinyapps.io/XDeathDB/ .

Proteomic profiling of MIS-C patients indicates heterogeneity relating to interferon gamma dysregulation and vascular endothelial dysfunction

Multi-system Inflammatory Syndrome in Children (MIS-C) is a major complication of Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) infection in pediatric patients. Weeks after an often mild or asymptomatic initial infection with SARS-CoV-2 children may present with a severe shock-like picture and marked inflammation. Children with MIS-C present with varying degrees of cardiovascular and hyperinflammatory symptoms. Here we perform a comprehensive analysis of the plasma proteome of more than 1400 proteins in children with SARS-CoV-2. We hypothesize that the proteome would reflect heterogeneity in hyperinflammation and vascular injury, and further identify pathogenic mediators of disease. We show that protein signatures demonstrate overlap between MIS-C, and the inflammatory syndromes macrophage activation syndrome (MAS) and thrombotic microangiopathy (TMA). We demonstrate that PLA2G2A is an important marker of MIS-C that associates with TMA. We find that IFNγ responses are dysregulated in MIS-C patients, and that IFNγ levels delineate clinical heterogeneity.

Proteomic Analysis of Mesenchymal Stromal Cell-Derived Extracellular Vesicles and Reconstructed Membrane Particles

Extracellular vesicles (EV) derived from mesenchymal stromal cells (MSC) are a potential therapy for immunological and degenerative diseases. However, large-scale production of EV free from contamination by soluble proteins is a major challenge. The generation of particles from isolated membranes of MSC, membrane particles (MP), may be an alternative to EV. In the present study we generated MP from the membranes of lysed MSC after removal of the nuclei. The yield of MP per MSC was 1 × 10⁵ times higher than EV derived from the same number of MSC. To compare the proteome of MP and EV, proteomic analysis of MP and EV was performed. MP contained over 20 times more proteins than EV. The proteins present in MP evidenced a multi-organelle origin of MP. The projected function of the proteins in EV and MP was very different. Whilst proteins in EV mainly play a role in extracellular matrix organization, proteins in MP were interconnected in diverse molecular pathways, including protein synthesis and degradation pathways and demonstrated enzymatic activity. Treatment of MSC with IFNγ led to a profound effect on the protein make up of EV and MP, demonstrating the possibility to modify the phenotype of EV and MP through modification of parent MSC. These results demonstrate that MP are an attractive alternative to EV for the development of potential therapies. Functional studies will have to demonstrate therapeutic efficacy of MP in preclinical disease models.

Joint-Tissue Integrative Analysis Identified Hundreds of Schizophrenia Risk Genes

Genome-wide association studies (GWAS) have identified a large number of schizophrenia risk variants, and most of them are mapped to noncoding regions. By leveraging multiple joint-tissue gene expression data and GWAS data, we herein performed a transcriptome-wide association study (TWAS) and Mendelian randomization (MR) analysis and identified 144 genes whose mRNA levels were related to genetic risk of schizophrenia. Most of these genes exhibited diametrically opposite trends of expression in prenatal and postnatal brain tissues, despite that their expression levels in dorsolateral prefrontal cortex (DLPFC) tissues did not significantly differ between schizophrenics and healthy controls. We then found significant enrichment of these genes in dopamine-related pathways that were repeatedly implicated in schizophrenia pathogenesis and in the action of antipsychotic drugs. Gene expression analysis using single cell RNA-sequencing (scRNA-seq) data of mid-gestation fetal brains further revealed enrichment of these genes in glutamatergic excitatory neurons and cycling progenitors. These lines of evidence, in consistency with previous findings, confirmed the polygenic nature of schizophrenia and highlighted involvement of early neurodevelopment aberrations in this disorder. Further investigations using advanced algorithms in both bulk brain tissues and in single cells and at different developmental stages are necessary to characterize transcriptomic features of schizophrenia pathogenesis along brain development.

Impaired HMG-CoA Reductase Activity Caused by Genetic Variants or Statin Exposure: Impact on Human Adipose Tissue, β-Cells and Metabolome

Inhibition of 3-hydroxy-3-methyl-glutaryl-CoA (HMG-CoA) reductase is associated with an increased risk of new-onset type 2 diabetes. We studied the association of genetic or pharmacological HMG-CoA reductase inhibition with plasma and adipose tissue (AT) metabolome and AT metabolic pathways. We also investigated the effects of statin-mediated pharmacological inhibition of HMG-CoA reductase on systemic insulin sensitivity by measuring the HOMA-IR index in subjects with or without statin therapy. The direct effects of simvastatin (20-250 nM) or its active metabolite simvastatin hydroxy acid (SA) (8-30 nM) were investigated on human adipocyte glucose uptake, lipolysis, and differentiation and pancreatic insulin secretion. We observed that the LDL-lowering HMGCR rs12916-T allele was negatively associated with plasma phosphatidylcholines and sphingomyelins, and HMGCR expression in AT was correlated with various metabolic and mitochondrial pathways. Clinical data showed that statin treatment was associated with HOMA-IR index after adjustment for age, sex, BMI, HbA1c, LDL-c levels, and diabetes status in the subjects. Supra-therapeutic concentrations of simvastatin reduced glucose uptake in adipocytes and normalized fatty acid-induced insulin hypersecretion from β-cells. Our data suggest that inhibition of HMG-CoA reductase is associated with insulin resistance. However, statins have a very mild direct effect on AT and pancreas, hence, other tissues as the liver or muscle appear to be of greater importance.

Whole-Genome Comparative Analysis Reveals Association Between Salmonella Genomic Variation and Egg Production Systems

Non-typhoidal Salmonella, particularly Salmonella enterica serovar Typhimurium (S. Typhimurium), is the predominant endemic serovar in the Australian egg production industry and is one of the most frequently reported serovars in foodborne infections in Australia. This study was conducted to investigate the genomic characteristics of Salmonella isolated from retail table eggs in Western Australia and to identify the impact of production systems on genomic characteristics of Salmonella such as virulence and antimicrobial resistance. A total of 40 non-typhoidal Salmonella isolates [S. Typhimurium isolates (n = 28) and Salmonella Infantis isolates (n = 12)] sourced from retail eggs produced by different production systems (barn-laid, cage, and free-range) in Western Australia were sequenced by whole-genome sequencing. The isolates were de novo assembled, annotated, and analyzed. The results indicated an association between Salmonella genomic variation and the system used to raise poultry for egg production (p-value < 0.05). All but one of the S. Infantis isolates were recovered from eggs collected from poultry raised under barn and cage production systems. A higher proportion (83.3%) of S. Typhimurium isolates were recovered from the eggs produced by free-range production system as compared with those produced under barn (76.9%) and cage production systems (53.3%). Our analysis indicated that Salmonella isolated from the eggs produced by barn and cage production systems had more virulence genes than the isolates of the free-range produced eggs. A low carriage of antimicrobial-resistant gene was detected in the isolates of this study. We have built a Salmonella genomics database and characteristics-linked gene pools to facilitate future study, characterization, and tracing of Salmonella outbreaks.

DNA Methylation Profiling for the Diagnosis and Prognosis of Patients with Nontuberculous Mycobacterium Lung Disease

The incidence of nontuberculous Mycobacterium (NTM) lung disease is rapidly increasing; however, its diagnosis and prognosis remain unclear while selecting patients who will respond to appropriate treatment. Differences in DNA methylation patterns between NTM patients with good or poor prognosis could provide important therapeutic targets. We used the Illumina MethylationEPIC (850k) DNA methylation microarray to determine the pattern between differentially methylated regions (DMRs) in NTM patients with good or poor prognosis (n = 4/group). Moreover, we merged and compared 20 healthy controls from previous Illumina Methylation450k DNA methylation microarray data. We selected and visualized the DMRs in the form of heatmaps, and enriched terms associated with these DMRs were identified by functional annotation with the “pathfinder” package. In total, 461 and 293 DMRs (|Log2 fold change| > 0.1 and P < 0.03) were more methylated in patients with four poor and four good prognoses, respectively. Furthermore, 337 and 771 DMRs (|Log2 fold change| > 0.08 and P < 0.001) were more methylated in eight NTM patients and 20 healthy controls, respectively. TGFBr1 was significantly less methylated, whereas HLA-DR1 and HLA-DR5 were more methylated in patients with poor prognosis (compared to those with good prognosis). LRP5, E2F1, and ADCY3 were the top three less-methylated genes in NTM patients (compared with the controls). The mTOR and Wnt signaling pathway-related genes were less methylated in patients with NTM. Collectively, genes related to Th1- cell differentiation, such as TGFBr1 and HLA-DR, may be used as biomarkers for predicting the treatment response in patients with NTM lung disease.

Bioinformatic Analysis Reveals Central Role for Tumor-Infiltrating Immune Cells in Uveal Melanoma Progression

Tumor-infiltrating immune cells are capable of effective cancer surveillance, and their abundance is linked to better prognosis in numerous tumor types. However, in uveal melanoma (UM), extensive immune infiltrate is associated with poor survival. This study aims to decipher the role of different tumor-infiltrating cell subsets in UM in order to identify potential targets for future immunotherapeutic treatment. We have chosen the TCGA-UVM cohort as a training dataset and GSE22138 as a testing dataset by mining publicly available databases. The abundance of 22 immune cell types was estimated using CIBERSORTx. Then, to determine the significance of tumor-infiltrating cell subsets in UM, we built a multicell type prognostic signature, which was validated in the testing cohort. The created signature was built upon the negative prognostic role of CD8+ T cells and M0 macrophages and the positive role of neutrophils. Based on the created signature score, we divided the patients into low- and high-risk groups. Kaplan-Meier, Cox, and ROC analyses demonstrated superior performance of our risk score compared to either clinical or pathologic characteristics of both cohorts. Further, we found the molecular pathways associated with cancer immunoevasion and metastasis to be enriched in the high-risk group, explaining both the lack of adequate immune surveillance despite increased infiltration of CD8+ T cells as well as the higher metastatic potential. Genes associated with tryptophan metabolism (IDO1 and KYNU) and metalloproteinases were among the most differentially expressed between the high- and low-risk groups. Our correlation analyses interpreted in context of published in vitro data strongly suggest the central role of CD8+ T cells in shifting the UM tumor microenvironment towards suppressive and metastasis-promoting. Therefore, we propose further investigations of IDO1 and metalloproteinases as novel targets for immunotherapy in lymphocyte-rich metastatic UM patients.

NEBULA is a fast negative binomial mixed model for differential or co-expression analysis of large-scale multi-subject single-cell data

The increasing availability of single-cell data revolutionizes the understanding of biological mechanisms at cellular resolution. For differential expression analysis in multi-subject single-cell data, negative binomial mixed models account for both subject-level and cell-level overdispersions, but are computationally demanding. Here, we propose an efficient NEgative Binomial mixed model Using a Large-sample Approximation (NEBULA). The speed gain is achieved by analytically solving high-dimensional integrals instead of using the Laplace approximation. We demonstrate that NEBULA is orders of magnitude faster than existing tools and controls false-positive errors in marker gene identification and co-expression analysis. Using NEBULA in Alzheimer's disease cohort data sets, we found that the cell-level expression of APOE correlated with that of other genetic risk factors (including CLU, CST3, TREM2, C1q, and ITM2B) in a cell-type-specific pattern and an isoform-dependent manner in microglia. NEBULA opens up a new avenue for the broad application of mixed models to large-scale multi-subject single-cell data.