Approaches for integrating heterogeneous RNA-seq data reveal cross-talk between microbes and genes in asthmatic patients.

Sputum induction is a non-invasive method to evaluate the airway environment, particularly for asthma. RNA sequencing (RNA-seq) of sputum samples can be challenging to interpret due to the complex and heterogeneous mixtures of human cells and exogenous (microbial) material. In this study, we develop a pipeline that integrates dimensionality reduction and statistical modeling to grapple with the heterogeneity. LDA(Latent Dirichlet allocation)-link connects microbes to genes using reduced-dimensionality LDA topics. We validate our method with single-cell RNA-seq and microscopy and then apply it to the sputum of asthmatic patients to find known and novel relationships between microbes and genes

Association between Tumor Microbiome and Hypoxia across Anatomic Subsites of Head and Neck Cancers

Purpose/Objective(s): Microbiome has been shown to affect tumorigenesis by promoting inflammation. However, the association between the upper aerodigestive microbiome and head and neck squamous cell carcinoma (HNSCC) is not well established. Hypoxia is a modifiable factor associated with poor radiation response. Our study analyzed the HNSCC tumor samples from The Cancer Genome Atlas (TCGA) to investigate the relationship between different HNSCC tumor subsites, hypoxia, and local tumor microbiome composition. Results: A total of 357 patients were included [Oral cavity (OC) = 226, Oropharynx (OPx) = 53, and Larynx/Hypopharynx (LHPx) = 78], of which 12.8%, 71.7%, and 10.3%, respectively, were HPV positive. The mean (SD) hypoxia scores were 30.18 (11.10), 24.31 (14.13), and 29.53 (12.61) in OC, OPx, and LHPx tumors, respectively, with higher values indicating greater hypoxia. The hypoxia score was significantly higher for OC tumors compared to OPx (p = 0.044) and LHPx (p = 0.002). There was no significant correlation between hypoxia and HPV status. Pseudomonas sp. in OC, Actinomyces sp. and Sulfurimonas sp. in OPx, and Filifactor, Pseudomonas and Actinomyces sp. in LHPx had the strongest association with the hypoxia score. Materials/Methods: Tumor RNAseq samples from TCGA were processed, and the R package “tmesig” was used to calculate gene expression signature, including the Buffa hypoxia (BH) score, a validated hypoxia signature using 52 hypoxia-regulated genes. Microbe relative abundances were modeled with primary tumor location and a high vs. low tertile BH score applying a gamma-distributed generalized linear regression using the “stats” package in R, with adjusted p-value < 0.05 considered significant. Conclusions: In our study, oral cavity tumors were found to be more hypoxic compared to other head and neck subsites, which could potentially contribute to their radiation resistance. For each subsite, distinct microbial populations were over-represented in hypoxic tumors in a subsite-specific manner. Further studies focusing on an association between microbiome, hypoxia, and patient outcomes are warranted