Metabolic elevation in soft-tissue sarcomas (STS), as documented with 18F-Fluorodeoxyglucose positron emission tomography (18F-FDG-PET/CT) has been linked with cell proliferation, higher grade, and lower survivals. However, the recent diagnostic innovations (CINSARC gene-expression signature and tertiary lymphoid structure [TLS]) and therapeutic innovations (immune checkpoint inhibitors [ICIs]) for STS patients underscore the need to re-assess the role of 18F-FDG-PET/CT. Thus, in this correspondence, our objective was to investigate the correlations between STS metabolism as assessed by nuclear imaging, and the immune landscape as estimated by transcriptomics analysis, immunohistochemistry panels, and TLS assessment. Based on a prospective cohort of 85 adult patients with high-grade STS recruited in the NEOSARCOMICS trial (NCT02789384), we identified 3 metabolic groups according to 18F-FDG-PET/CT metrics (metabolic-low [60%], -intermediate [15.3%] and high [24.7%]). We found that T-cells CD8 pathway was significantly enriched in metabolic-high STS. Conversely, several pathways involved in antitumor immune response, cell differentiation and cell cycle, were downregulated in extreme metabolic-low STS. Next, multiplex immunofluorescence showed that densities of CD8+, CD14+, CD45+, CD68+, and c-MAF cells were significantly higher in the metabolic-high group compared to the metabolic-low group. Lastly, no association was found between metabolic group and TLS status. Overall, these results suggest that (i) rapidly proliferating and metabolically active STS can instigate a more robust immune response, thereby attracting immune cells such as T cells and macrophages, and (ii) metabolic activity and TLS could independently influence immune responses.

BACKGROUND: For decades, the basic treatment strategies of necrotizing soft tissue infections (NSTI) have remained unchanged, primarily relying on aggressive surgical removal of infected tissue, broad-spectrum antibiotics, and supportive intensive care. One treatment strategy that has been proposed as an adjunctive measure to improve patient outcomes is hyperbaric oxygen (HBO2) treatment. HBO2 treatment has been linked to several immune modulatory effects; however, investigating these effects is complicated due to the disease\'s acute life-threatening nature, metabolic and cell homeostasis dependent variability in treatment effects, and heterogeneity with respect to both patient characteristics and involved pathogens. To embrace this complexity, we aimed to explore the underlying biological mechanisms of HBO2 treatment in patients with NSTI on the gene expression level.
METHODS: We conducted an observational cohort study on prospective collected data, including 85 patients admitted to the intensive care unit (ICU) for NSTI. All patients were treated with one or two HBO2 treatments and had one blood sample taken before and after the intervention. Total RNAs from blood samples were extracted and mRNA purified with rRNA depletion, followed by whole-transcriptome RNA sequencing with a targeted sequencing depth of 20 million reads. A model for differentially expressed genes (DEGs) was fitted, and the functional aspects of the obtained set of genes was predicted with GO (Gene Ontology) and KEGG (Kyoto Encyclopedia of genes and Genomes) enrichment analyses. All analyses were corrected for multiple testing with FDR.
RESULTS: After sequential steps of quality control, a final of 160 biological replicates were included in the present study. We found 394 protein coding genes that were significantly DEGs between the two conditions with FDR < 0.01, of which 205 were upregulated and 189 were downregulated. The enrichment analysis of these DEGs revealed 20 GO terms in biological processes and 12 KEGG pathways that were significantly overrepresented in the upregulated DEGs, of which the term; \"adaptive immune response\" (GO:0002250) (FDR = 9.88E-13) and \"T cell receptor signaling pathway\" (hsa04660) (FDR = 1.20E-07) were the most significant. Among the downregulated DEGs two biological processes were significantly enriched, of which the GO term \"apoptotic process\" (GO:0006915) was the most significant (FDR = 0.001), followed by \"Positive regulation of T helper 1 cell cytokine production\" (GO:2000556), and \"NF-kappa B signaling pathway\" (hsa04064) was the only KEGG pathway that was significantly overrepresented (FDR = 0.001).
CONCLUSIONS: When one or two sessions of HBO2 treatment were administered to patients with a dysregulated immune response and systemic inflammation due to NSTI, the important genes that were regulated during the intervention were involved in activation of T helper cells and downregulation of the disease-induced highly inflammatory pathway NF-κB, which was associated with a decrease in the mRNA level of pro-inflammatory factors.
BACKGROUND: Biological material was collected during the INFECT study, registered at ClinicalTrials.gov (NCT01790698).

Kaposi sarcoma (KS) is the most common cancer in people living with HIV (PLWH) in many countries where KS-associated herpesvirus is endemic. Treatment has changed little in 20 years, but the disease presentation has. This prospective cohort study enrolled 122 human immunodeficiency virus (HIV) positive KS patients between 2017 and 2019 in Malawi. Participants were treated with bleomycin, vincristine and combination antiretroviral therapy, the local standard of care. One-year overall survival was 61%, and progression-free survival was 58%. The 48-week complete response rate was 35%. RNAseq (n = 78) differentiated two types of KS lesions, those with marked endothelial characteristics and those enriched in inflammatory transcripts. This suggests that different KS lesions are in different disease states consistent with the known heterogeneous clinical response to treatment. In contrast to earlier cohorts, the plasma HIV viral load of KS patients in our study was highly variable. A total of 25% of participants had no detectable HIV; all had detectable KSHV viral load. Our study affirms that many KS cases today develop in PLWH with well-controlled HIV infection and that different KS lesions have differing molecular compositions. Further studies are needed to develop predictive biomarkers for this disease.

Fasciolosis is a globally widespread trematodiasis with a major economic and veterinary impact. Therefore, this disease is responsible for millions of dollars in losses to the livestock industry, and also constitutes an emerging human health problem in endemic areas. The ubiquitous nature of Fasciola hepatica, the main causative agent, is one of the key factors for the success of fasciolosis. Accordingly, this parasite is able to subsist in a wide variety of ecosystems and hosts, thanks to the development of a plethora of strategies for adaption and immune evasion. Fasciolosis comprises a growing concern due to its high prevalence rates, together with the emergence of strains of the parasite resistant to the treatment of choice (triclabendazole). These facts highlight the importance of developing novel control measures which allow for an effective protection against the disease before F. hepatica settles in a niche inaccessible to the immune system. However, knowledge about the initial phases of the infection, including the migration mechanisms of the parasite and the early innate host response, is still scarce. Recently, our group developed an in vitro host-parasite interaction model that allowed the early events to be unveiled after the first contact between the both actors. This occurs shortly upon ingestion of F. hepatica metacercariae and the emergence of the newly excysted juveniles (FhNEJ) in the host duodenum. Here, we present a transcriptomic analysis of such model using an approach based on RNA sequencing (RNA-Seq), which reveals changes in gene expression related to proteolysis and uptake of metabolites in FhNEJ. Additionally, contact with the parasite triggered changes in host intestinal cells related to pseudogenes expression and host defence mechanisms, including immune response, among others. In sum, these results provide a better understanding of the early stages of fasciolosis at molecular level, and a pool of targets that could be used in future therapeutic strategies against the disease.

In differential gene expression data analysis, one objective is to identify groups of co-expressed genes from a large dataset in order to detect the association between such a group of genes and an experimental condition. This is often done through a clustering approach, such as k-means or bipartition hierarchical clustering, based on particular similarity measures in the grouping process. In such a dataset, the gene differential expression itself is an innate attribute that can be used in the feature extraction process. For example, in a dataset consisting of multiple treatments versus their controls, the expression of a gene in each treatment would have three possible behaviors, upregulated, downregulated, or unchanged. We present in this chapter, a differential expression feature extraction (DEFE) method by using a string consisting of three numerical values at each character to denote such behavior, i.e., 1 = up, 2 = down, and 0 = unchanged, which results in up to 3B differential expression patterns across all B comparisons. This approach has been successfully applied in many research projects, and among these, we demonstrate the strength of DEFE in a case study on RNA-sequencing (RNA-seq) data analysis of wheat challenged with the phytopathogenic fungus, Fusarium graminearum. Combinations of multiple schemes of DEFE patterns revealed groups of genes putatively associated with resistance or susceptibility to FHB.

Host responses to vaccines are complex but important to investigate. To facilitate the study, we have developed a tool called Vaccine Induced Gene Expression Analysis Tool (VIGET), with the aim to provide an interactive online tool for users to efficiently and robustly analyze the host immune response gene expression data collected in the ImmPort/GEO databases. VIGET allows users to select vaccines, choose ImmPort studies, set up analysis models by choosing confounding variables and two groups of samples having different vaccination times, and then perform differential expression analysis to select genes for pathway enrichment analysis and functional interaction network construction using the Reactome\'s web services. VIGET provides features for users to compare results from two analyses, facilitating comparative response analysis across different demographic groups. VIGET uses the Vaccine Ontology (VO) to classify various types of vaccines such as live or inactivated flu vaccines, yellow fever vaccines, etc. To showcase the utilities of VIGET, we conducted a longitudinal analysis of immune responses to yellow fever vaccines and found an intriguing complex activity response pattern of pathways in the immune system annotated in Reactome, demonstrating that VIGET is a valuable web portal that supports effective vaccine response studies using Reactome pathways and ImmPort data.

The etiology and most risk factors for a sporadic first primary neoplasm in childhood or subsequent second primary neoplasms are still unknown. One established causal factor for therapy-associated second primary neoplasms is the exposure to ionizing radiation during radiation therapy as a mainstay of cancer treatment. Second primary neoplasms occur in 8% of all cancer survivors within 30 years after the first diagnosis in Germany, but the underlying factors for intrinsic susceptibilities have not yet been clarified. Thus, the purpose of this nested case-control study was the investigation and comparison of gene expression and affected pathways in primary fibroblasts of childhood cancer survivors with a first primary neoplasm only or with at least one subsequent second primary neoplasm, and controls without neoplasms after exposure to a low and a high dose of ionizing radiation.
Primary fibroblasts were obtained from skin biopsies from 52 adult donors with a first primary neoplasm in childhood (N1), 52 with at least one additional primary neoplasm (N2+), as well as 52 without cancer (N0) from the KiKme study. Cultured fibroblasts were exposed to a high [2 Gray (Gy)] and a low dose (0.05 Gy) of X-rays. Messenger ribonucleic acid was extracted 4 h after exposure and Illumina-sequenced. Differentially expressed genes (DEGs) were computed using limma for R, selected at a false discovery rate level of 0.05, and further analyzed for pathway enrichment (right-tailed Fisher\'s Exact Test) and (in-) activation (z ≥|2|) using Ingenuity Pathway Analysis.
After 0.05 Gy, least DEGs were found in N0 (n = 236), compared to N1 (n = 653) and N2+ (n = 694). The top DEGs with regard to the adjusted p-value were upregulated in fibroblasts across all donor groups (SESN1, MDM2, CDKN1A, TIGAR, BTG2, BLOC1S2, PPM1D, PHLDB3, FBXO22, AEN, TRIAP1, and POLH). Here, we observed activation of p53 Signaling in N0 and to a lesser extent in N1, but not in N2+. Only in N0, DNA (excision-) repair (involved genes: CDKN1A, PPM1D, and DDB2) was predicted to be a downstream function, while molecular networks in N2+ were associated with cancer, as well as injury and abnormalities (among others, downregulation of MSH6, CCNE2, and CHUK). After 2 Gy, the number of DEGs was similar in fibroblasts of all donor groups and genes with the highest absolute log2 fold-change were upregulated throughout (CDKN1A, TIGAR, HSPA4L, MDM2, BLOC1SD2, PPM1D, SESN1, BTG2, FBXO22, PCNA, and TRIAP1). Here, the p53 Signaling-Pathway was activated in fibroblasts of all donor groups. The Mitotic Roles of Polo Like Kinase-Pathway was inactivated in N1 and N2+. Molecular Mechanisms of Cancer were affected in fibroblasts of all donor groups. P53 was predicted to be an upstream regulator in fibroblasts of all donor groups and E2F1 in N1 and N2+. Results of the downstream analysis were senescence in N0 and N2+, transformation of cells in N0, and no significant effects in N1. Seven genes were differentially expressed in reaction to 2 Gy dependent on the donor group (LINC00601, COBLL1, SESN2, BIN3, TNFRSF10A, EEF1AKNMT, and BTG2).
Our results show dose-dependent differences in the radiation response between N1/N2+ and N0. While mechanisms against genotoxic stress were activated to the same extent after a high dose in all groups, the radiation response was impaired after a low dose in N1/N2+, suggesting an increased risk for adverse effects including carcinogenesis, particularly in N2+.

Direct sequencing of single molecules through nanopores allows for accurate quantification and full-length characterization of native RNA or complementary DNA (cDNA) without amplification. Both nanopore-based native RNA and cDNA approaches involve complex transcriptome procedures at a lower cost. However, there are several differences between the two approaches. In this study, we perform matched native RNA sequencing and cDNA sequencing to enable relevant comparisons and evaluation. Using Saccharomyces cerevisiae, a eukaryotic model organism widely used in industrial biotechnology, two different growing conditions are considered for comparison, including the poly-A messenger RNA isolated from yeast cells grown in minimum media under respirofermentative conditions supplemented with glucose (glucose growth conditions) and from cells that had shifted to ethanol as a carbon source (ethanol growth conditions). Library preparation for direct RNA sequencing is shorter than that for direct cDNA sequencing. The sequence characteristics of the two methods were different, such as sequence yields, quality score of reads, read length distribution, and mapped on reference ability of reads. However, differential gene expression analyses derived from the two approaches are comparable. The unique feature of direct RNA sequencing is RNA modification; we found that the RNA modification at the 5\' end of a transcript was underestimated due to the 3\' bias behavior of the direct RNA sequencing. Our comprehensive evaluation from this work could help researchers make informed choices when selecting an appropriate long-read sequencing method for understanding gene functions, pathways, and detailed functional characterization.

High NUE (nitrogen use efficiency) has great practical significance for sustainable crop production. Wheat is one of the main cultivated crops worldwide for human food and nutrition. However, wheat grain productivity is dependent upon cultivars with high NUE in addition to the application of nitrogen fertilizers. In order to understand the molecular mechanisms exhibiting a high NUE response, a comparative transcriptomics study was carried out through RNA-seq analysis to investigate the gene expression that regulates NUE, in root and shoot tissue of N-efficient (PBW677) and N-inefficient (703) cultivars under optimum and nitrogen (N) stress. Differentially expressed gene analysis revealed a total of 2,406 differentially expressed genes (DEGs) present in both the contrasting cultivars under N stress. The efficient genotype PBW677 had considerably more abundant DEGs with 1,653 (903 roots +750 shoots) compared to inefficient cultivar PBW703 with 753 (96 roots +657 shoots). Gene ontology enrichment and pathway analysis of these DEGs suggested that the two cultivars differed in terms of adaptive mechanism. Gene enrichment analysis revealed that among the upregulated and downregulated genes the overrepresented and underrepresented gene categories belonged to biological processes like DNA binding, response to abiotic stimulus, photosynthesis, carbon fixation, carbohydrate metabolic process, nitrogen compound metabolic process, nitrate transport, and translation in cultivar PBW677, while the enriched biological processes were nucleosome assembly, chromatin remodeling, DNA packaging, lipid transport, sulfur compound metabolic process, protein modifications, and protein folding and refolding in N inefficient cultivar PBW703. We found several transcription factors (MYB, WRKY, RING finger protein, zinc finger protein, transporters, NRT1, amino acid transporters, sugar), protein kinases, and genes involved in N absorption, transportation, and assimilation to be highly expressed in high NUE cultivar PBW677. In our study, we report 13 potential candidate genes which showed alternate gene expression in the two contrasting cultivars under study. These genes could serve as potential targets for future breeding programs.

OBJECTIVE: Transcriptional profiling of immune cells is an indispensable tool in biomedical research; however, heterogenous sample types routinely used in transcriptomic studies may mask important cell type-specific transcriptional differences. Techniques to isolate desired cell types are used to overcome this limitation. We sought to evaluate the use of immunomagnetic B cell isolation on RNA quality and transcriptional output. Additionally, we aimed to develop a B cell gene signature representative of a freshly isolated B cell population to be used as a tool to verify isolation efficacy and to provide a transcriptional standard for evaluating maintenance or deviation from traditional B cell identity.
RESULTS: We found RNA quality and RNA-sequencing output to be comparable between donor-matched PBMC, whole blood, and B cells following negative selection by immunomagnetic B cell isolation. Transcriptional analysis enabled the development of an 85 gene B cell signature. This signature effectively clustered isolated B cells from heterogeneous sample types in our study and naïve and memory B cells when applied to transcriptional data from a published source. Additionally, by identifying B cell signature genes whose functional role in B cells is currently unknown, our gene signature has uncovered areas for future investigation.