Enhancing NSCLC recurrence prediction with PET/CT habitat imaging, ctDNA, and integrative radiogenomics-blood insights

While we recognize the prognostic importance of clinicopathological measures and circulating tumor DNA (ctDNA), the independent contribution of quantitative image markers to prognosis in non-small cell lung cancer (NSCLC) remains underexplored. In our multi-institutional study of 394 NSCLC patients, we utilize pre-treatment computed tomography (CT) and 18F-fluorodeoxyglucose positron emission tomography (FDG-PET) to establish a habitat imaging framework for assessing regional heterogeneity within individual tumors. This framework identifies three PET/CT subtypes, which maintain prognostic value after adjusting for clinicopathologic risk factors including tumor volume. Additionally, these subtypes complement ctDNA in predicting disease recurrence. Radiogenomics analysis unveil the molecular underpinnings of these imaging subtypes, highlighting downregulation in interferon alpha and gamma pathways in the high-risk subtype. In summary, our study demonstrates that these habitat imaging subtypes effectively stratify NSCLC patients based on their risk levels for disease recurrence after initial curative surgery or radiotherapy, providing valuable insights for personalized treatment approaches.

Synthetic PET from CT improves diagnosis and prognosis for lung cancer: Proof of concept

[18F]Fluorodeoxyglucose positron emission tomography (FDG-PET) and computed tomography (CT) are indispensable components in modern medicine. Although PET can provide additional diagnostic value, it is costly and not universally accessible, particularly in low-income countries. To bridge this gap, we have developed a conditional generative adversarial network pipeline that can produce FDG-PET from diagnostic CT scans based on multi-center multi-modal lung cancer datasets (n = 1,478). Synthetic PET images are validated across imaging, biological, and clinical aspects. Radiologists confirm comparable imaging quality and tumor contrast between synthetic and actual PET scans. Radiogenomics analysis further proves that the dysregulated cancer hallmark pathways of synthetic PET are consistent with actual PET. We also demonstrate the clinical values of synthetic PET in improving lung cancer diagnosis, staging, risk prediction, and prognosis. Taken together, this proof-of-concept study testifies to the feasibility of applying deep learning to obtain high-fidelity PET translated from CT.

Exploratory analysis of the cervix tumoral HPV antigen-specific T-cell repertoire during chemoradiation and after brachytherapy

BACKGROUND

Chemoradiation (CRT) may modulate the immune milieu as an in-situ vaccine. Rapid dose delivery of brachytherapy has unclear impact on T-cell repertoires. HPV-associated cancers express viral oncoproteins E6/E7, which enable tracking antigen/tumor-specific immunity during CRT.

METHODS

Thirteen cervical cancer patients on a multi-institutional prospective protocol from 1/2020-1/2023 underwent standard-of-care CRT with pulsed-dose-rate brachytherapy boost (2 fractions). Cervix swabs at various timepoints underwent multiplex DNA deep sequencing of the TCR-β/CDR3 region with immunoSEQ. Separately, HPV-responsive T-cell clones were also expanded ex vivo. Statistical analysis was via Mann-Whitney-U.

RESULTS

TCR productive clonality, templates, frequency, or rearrangements increased post-brachytherapy in 8 patients. Seven patients had E6/E7-responsive evolution over CRT with increased productive templates (ranges: 1.2-50.2 fold-increase from baseline), frequency (1.2-1.7), rearrangements (1.2-40.2), and clonality (1.2-15.4). Five patients had HPV-responsive clonal expansion post-brachytherapy, without changes in HPV non-responsive clones. Epitope mapping revealed VDJ rearrangements targeting cervical cancer-associated antigens in 5 patients. The only two patients with disease recurrence lacked response in all metrics. A lack of global TCR remodeling correlated with worse recurrence-free survival, p = 0.04.

CONCLUSION

CRT and brachytherapy alters the cervical cancer microenvironment to facilitate the expansion of specific T-cell populations, which may contribute to treatment efficacy.

HPV-related anal cancer is associated with changes in the anorectal microbiome during cancer development

Background

Squamous cell carcinoma of the anus (SCCA) is a rare gastrointestinal cancer. Factors associated with progression of HPV infection to anal dysplasia and cancer are unclear and screening guidelines and approaches for anal dysplasia are less clear than for cervical dysplasia. One potential contributing factor is the anorectal microbiome. In this study, we aimed to identify differences in anal microbiome composition in the settings of HPV infection, anal dysplasia, and anal cancer in this rare disease.

Methods

Patients were enrolled in two prospective studies. Patients with anal dysplasia were part of a cross-sectional cohort that enrolled women with high-grade lower genital tract dysplasia. Anorectal tumor swabs were prospectively collected from patients with biopsy-confirmed locally advanced SCCA prior to receiving standard-of-care chemoradiotherapy (CRT). Patients with high-grade lower genital tract dysplasia without anal dysplasia were considered high-risk (HR Normal). 16S V4 rRNA Microbiome sequencing was performed for anal swabs. Alpha and Beta Diversity and composition were compared for HR Normal, anal dysplasia, and anal cancer.

Results

60 patients with high-grade lower genital tract dysplasia were initially enrolled. Seven patients had concurrent anal dysplasia and 44 patients were considered HR Normal. Anorectal swabs from 21 patients with localized SCCA were included, sequenced, and analyzed in the study. Analysis of weighted and unweighted UniFrac distances demonstrated significant differences in microbial community composition between anal cancer and HR normal (p=0.018). LEfSe identified that all three groups exhibited differential enrichment of specific taxa. Peptoniphilus (p=0.028), Fusobacteria (p=0.0295), Porphyromonas (p=0.034), and Prevotella (p=0.029) were enriched in anal cancer specimens when compared to HR normal.

Conclusion

Although alpha diversity was similar between HR Normal, dysplasia and cancer patients, composition differed significantly between the three groups. Increased anorectal Peptoniphilus, Fusobacteria, Porphyromonas, and Prevotella abundance were associated with anal cancer. These organisms have been reported in various gastrointestinal cancers with roles in facilitating the proinflammatory microenvironment and neoplasia progression. Future work should investigate a potential role of microbiome analysis in screening for anal dysplasia and investigation into potential mechanisms of how these microbial imbalances influence the immune system and anal carcinogenesis.

Serial Genotyping of the Human Papillomavirus in Cervical Cancer: An Insight Into Virome Dynamics During Chemoradiation Therapy

PURPOSE

Human papillomavirus (HPV) is the primary driver of cervical cancer. Although studies in other malignancies correlated peripheral blood DNA clearance with favorable outcomes, research on the prognostic value of HPV clearance in gynecologic cancers using intratumoral HPV is scarce. We aimed to quantify the intratumoral HPV virome in patients undergoing chemoradiation therapy (CRT) and associate this with clinical characteristics and outcomes.

METHODS AND MATERIALS

This prospective study enrolled 79 patients with stage IB-IVB cervical cancer undergoing definitive CRT. Cervical tumor swabs collected at baseline and week 5 (end of intensity modulated radiation therapy) were sent for shotgun metagenome sequencing and processed via VirMAP, a viral genome sequencing and identification tool for all known HPV types. The data were categorized into HPV groups (16, 18, high risk [HR], and low risk [LR]). We used independent t tests and Wilcoxon signed-rank to compare continuous variables and χ² and Fisher exact tests to compare categorical variables. Kaplan-Meier survival modeling was performed with log-rank testing. HPV genotyping was verified using quantitative polymerase chain reaction to validate VirMAP results using receiver operating characteristic curve and Cohen's kappa.

RESULTS

At baseline, 42%, 12%, 25%, and 16% of patients were positive for HPV 16, HPV 18, HPV HR, and HPV LR, respectively, and 8% were HPV negative. HPV type was associated with insurance status and CRT response. Patients with HPV 16+ and other HPV HR+ tumors were significantly more likely to have a complete response to CRT versus patients with HPV 18 and HPV LR/HPV-negative tumors. Overall HPV viral loads predominantly decreased throughout CRT, except for HPV LR viral load.

CONCLUSIONS

Rarer, less well-studied HPV types in cervical tumors are clinically significant. HPV 18 and HPV LR/negative tumors are associated with poor CRT response. This feasibility study provides a framework for a larger study of intratumoral HPV profiling to predict outcomes in patients with cervical cancer.

Feasibility of a novel non-invasive swab technique for serial whole-exome sequencing of cervical tumors during chemoradiation therapy

Background

Clinically relevant genetic predictors of radiation response for cervical cancer are understudied due to the morbidity of repeat invasive biopsies required to obtain genetic material. Thus, we aimed to demonstrate the feasibility of a novel noninvasive cervical swab technique to (1) collect tumor DNA with adequate throughput to (2) perform whole-exome sequencing (WES) at serial time points over the course of chemoradiation therapy (CRT).

Methods

Cervical cancer tumor samples from patients undergoing chemoradiation were collected at baseline, at week 1, week 3, and at the completion of CRT (week 5) using a noninvasive swab-based biopsy technique. Swab samples were analyzed with whole-exome sequencing (WES) with mutation calling using a custom pipeline optimized for shallow whole-exome sequencing with low tumor purity (TP). Tumor mutation changes over the course of treatment were profiled.

Results

216 samples were collected and successfully sequenced for 70 patients (94% of total number of tumor samples collected). A total of 33 patients had a complete set of samples at all four time points. The mean mapping rate was 98% for all samples, and the mean target coverage was 180. Estimated TP was greater than 5% for all samples. Overall mutation frequency decreased during CRT but mapping rate and mean target coverage remained at >98% and >180 reads at week 5.

Conclusion

This study demonstrates the feasibility and application of a noninvasive swab-based technique for WES analysis which may be applied to investigate dynamic tumor mutational changes during treatment to identify novel genes which confer radiation resistance.

Fusobacterium is enriched in oral cancer and promotes induction of programmed death-ligand 1 (PD-L1)

Recently, increased number of studies have demonstrated a relationship between the oral microbiome and development of head and neck cancer, however, there are few studies to investigate the role of oral bacteria in the context of the tumor microenvironment in a single head and neck subsite. Here, paired tumor and adjacent normal tissues from thirty-seven oral tongue squamous cell carcinoma (SCC) patients were subjected to 16S rRNA gene sequencing and whole exome sequencing (WES), in addition to RNA sequencing for tumor samples. We observed that Fusobacterium was significantly enriched in oral tongue cancer and that Rothia and Streptococcus were enriched in adjacent normal tissues. A decrease in alpha diversity was found in tumor when compared to adjacent normal tissues. While increased Fusobacterium in tumor samples was not associated with changes in immune cell infiltration, it was associated with increased PD-L1 mRNA expression. Therefore, we examined the effects of Fusobacterium on PD-L1 expression in head and neck SCC cell lines. We demonstrated that infection with Fusobacterium species can increase both PD-L1 mRNA and surface PD-L1 protein expression on head and neck cancer cell lines. The correlation between Fusobacterium and PD-L1 expression in oral tongue SCC, in conjunction with the ability of the bacterium to induce PD-L1 expression in vitro suggests a potential role for Fusobacterium on modulation of the tumor immune microenvironment in head and neck cancer.

Metagenomes of rectal swabs in larger, advanced stage cervical cancers have enhanced mucus degrading functionalities and distinct taxonomic structure

Background

Gut microbiome community composition differs between cervical cancer (CC) patients and healthy controls, and increased gut diversity is associated with improved outcomes after treatment. We proposed that functions of specific microbial species adjoining the mucus layer may directly impact the biology of CC.

Method

Metagenomes of rectal swabs in 41 CC patients were examined by whole-genome shotgun sequencing to link taxonomic structures, molecular functions, and metabolic pathway to patient’s clinical characteristics.

Results

Significant association of molecular functions encoded by the metagenomes was found with initial tumor size and stage. Profiling of the molecular function abundances and their distributions identified 2 microbial communities co-existing in each metagenome but having distinct metabolism and taxonomic structures. Community A (Clostridia and Proteobacteria predominant) was characterized by high activity of pathways involved in stress response, mucus glycan degradation and utilization of degradation byproducts. This community was prevalent in patients with larger, advanced stage tumors. Conversely, community B (Bacteroidia predominant) was characterized by fast growth, active oxidative phosphorylation, and production of vitamins. This community was prevalent in patients with smaller, early-stage tumors.

Conclusions

In this study, enrichment of mucus degrading microbial communities in rectal metagenomes of CC patients was associated with larger, more advanced stage tumors.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12885-022-09997-0.

Expansion of Candidate HPV-Specific T Cells in the Tumor Microenvironment during Chemoradiotherapy Is Prognostic in HPV16+ Cancers

Human papillomavirus (HPV) infection causes 600,000 new cancers worldwide each year. HPV-related cancers express the oncogenic proteins E6 and E7, which could serve as tumor-specific antigens. It is not known whether immunity to E6 and E7 evolves during chemoradiotherapy or affects survival. Using T cells from 2 HPV16⁺ patients, we conducted functional T-cell assays to identify candidate HPV-specific T cells and common T-cell receptor motifs, which we then analyzed across 86 patients with HPV-related cancers. The HPV-specific clones and E7-related T-cell receptor motifs expanded in the tumor microenvironment over the course of treatment, whereas non-HPV-specific T cells did not. In HPV16⁺ patients, improved recurrence-free survival was associated with HPV-responsive T-cell expansion during chemoradiotherapy.

Whole-Genome Sequencing of Common Salivary Gland Carcinomas: Subtype-Restricted and Shared Genetic Alterations

PURPOSE

Salivary gland carcinomas (SGCs) are pathologically classified into several widely diverse subtypes, of which adenoid cystic carcinoma (ACC), mucoepidermoid carcinoma (MEC), and salivary duct carcinoma (SDC) are the most commonly encountered. A comparative genetic analysis of these subtypes provides detailed information on the genetic alterations that are associated with their tumorigenesis and may lead to the identification of biomarkers to guide tumor-specific clinical trials.

EXPERIMENTAL DESIGN

Whole-genome sequencing of 58 common SGCs (20 ACCs, 20 SDCs, and 18 MECs) was performed to catalog structural variations, copy number, rearrangements, and driver mutations. Data were bioinformatically analyzed and correlated with clinicopathologic parameters, and selected targets were validated.

RESULTS

Novel and recurrent type-specific and shared genetic alterations were identified within and among 3 subtypes. Mutually exclusive canonical fusion and nonfusion genomic alterations were identified in both ACC and MEC. In ACCs, loss of chromosome 12q was dominant in MYB or MYBL1 fusion-positive tumors and mutations of NOTCH pathway were more common in these fusion negatives. In MECs, CRTC1-MAML2 fusion-positive tumors showed frequent BAP1 mutation, and tumors lacking this fusion were enriched with LRFN1 mutation. SDCs displayed considerable genetic instability, lacked recurrent chromosomal rearrangements, and demonstrated nonoverlapping TP53 mutation and ERBB2 amplification in a subset of tumors. Limited genetic alterations, including focal amplifications of 8q21-q23, were shared by all subtypes and were associated with poor survival.

CONCLUSIONS

This study delineates type-specific and shared genetic alterations that are associated with early phenotypic commitment and the biologic progression of common SGCs. These alterations, upon validation, could serve as biomarkers in tumor-specific clinical trials.

Gut microbiome diversity is an independent predictor of survival in cervical cancer patients receiving chemoradiation

Diversity of the gut microbiome is associated with higher response rates for cancer patients receiving immunotherapy but has not been investigated in patients receiving radiation therapy. Additionally, current studies investigating the gut microbiome and outcomes in cancer patients may not have adjusted for established risk factors. Here, we sought to determine if diversity and composition of the gut microbiome was independently associated with survival in cervical cancer patients receiving chemoradiation. Our study demonstrates that the diversity of gut microbiota is associated with a favorable response to chemoradiation. Additionally, compositional variation among patients correlated with short term and long-term survival. Short term survivor fecal samples were significantly enriched in Porphyromonas, Porphyromonadaceae, and Dialister, whereas long term survivor samples were significantly enriched in Escherichia Shigella, Enterobacteriaceae, and Enterobacteriales. Moreover, analysis of immune cells from cervical tumor brush samples by flow cytometry revealed that patients with a high microbiome diversity had increased tumor infiltration of CD4+ lymphocytes as well as activated subsets of CD4 cells expressing ki67+ and CD69+ over the course of radiation therapy. Modulation of the gut microbiota before chemoradiation might provide an alternative way to enhance treatment efficacy and improve treatment outcomes in cervical cancer patients.

Effect of Antibiotics on Gut and Vaginal Microbiomes Associated with Cervical Cancer Development in Mice

Antibiotics affect microbial diversity in the gut, leading to dysbiosis and impaired immunity. However, the impact of antibiotics on microbial communities at other sites, such as vagina is less understood. It is also not clear whether changes induced by antibiotics in both microbiomes affect the development of cervical cancer. In this study, we utilized the murine model to evaluate these questions. We show that oral application of broad-spectrum antibiotics in mice changed not only diversity, but composition and sharing of gut and vaginal microbiomes in mice and influenced cervical cancer development in an orthotopic tumor model. Antibiotics decreased richness and diversity indexes in the gut but increased them in the vagina. Some beneficial taxa, such as Bacteroides, Ruminococcaceae, and Lachnospiraceae increased their abundance in the vagina while other pathogenic species, such as Proteobacteria, were decreased. As a result of the changes, mice with greater richness and diversity of the vaginal microbiome after antibiotics exposure were less likely developed tumors. No association between richness and diversity of the gut microbiome and tumor development was identified.

Proteogenomic Analysis of Salivary Adenoid Cystic Carcinomas Defines Molecular Subtypes and Identifies Therapeutic Targets

PURPOSE

Salivary gland adenoid cystic carcinoma (ACC) has heterogeneous clinical behavior. Currently, all patients are treated uniformly, and no standard-of-care systemic therapy exists for metastatic ACC. We conducted an integrated proteogenomic analyses of ACC tumors to identify dysregulated pathways and propose a classification with therapeutic implications.

EXPERIMENTAL DESIGN

RNA/DNA sequencing of 54 flash-frozen salivary ACCs and reverse phase protein array (RPPA) in 38 specimens were performed, with validation by Western blotting and/or IHC. Three independent ACC cohorts were used for validation.

RESULTS

Both unbiased RNA sequencing (RNA-seq) and RPPA analysis revealed two molecular subtypes: ACC-I (37%) and ACC-II (63%). ACC-I had strong upregulation of MYC, MYC target genes, and mRNA splicing, enrichment of NOTCH-activating mutations, and dramatically worse prognosis. ACC-II exhibited upregulation of TP63 and receptor tyrosine kinases (AXL, MET, and EGFR) and less aggressive clinical course. TP63 and MYC were sufficient to assign tumors to ACC subtypes, which was validated in one independent cohort by IHC and two additional independent cohorts by RNA-seq. Furthermore, IHC staining for MYC and P63 protein levels can be used to identify ACC subtypes, enabling rapid clinical deployment to guide therapeutic decisions. Our data suggest a model in which ACC-I is driven by MYC signaling through either NOTCH mutations or direct amplification, which in turn suppress P63 signaling observed in ACC-II, producing unique therapeutic vulnerabilities for each subtype.

CONCLUSIONS

Cooccurrence of multiple actionable protein/pathways alterations in each subtype indicates unique therapeutic vulnerabilities and opportunities for optimal combination therapy for this understudied and heterogeneous disease.

Neoadjuvant Chemotherapy Increases Cytotoxic T Cell, Tissue Resident Memory T Cell, and B Cell Infiltration in Resectable NSCLC

INTRODUCTION

The combination of programmed cell death protein-1 or programmed death-ligand 1 immune checkpoint blockade and chemotherapy has revolutionized the treatment of advanced NSCLC, but the mechanisms underlying this synergy remain incompletely understood. In this study, we explored the relationships between neoadjuvant chemotherapy and the immune microenvironment (IME) of resectable NSCLC to identify novel mechanisms by which chemotherapy may enhance the effect of immune checkpoint blockade.

METHODS

Genomic, transcriptomic, and immune profiling data of 511 patients treated with neoadjuvant chemotherapy followed by surgery (NCT) versus upfront surgery (US) were compared with determined differential characteristics of the IMEs derived from whole-exome sequencing (NCT = 18; US = 73), RNA microarray (NCT = 45; US = 202), flow cytometry (NCT = 17; US = 39), multiplex immunofluorescence (NCT = 10; US = 72), T-cell receptor sequencing (NCT = 16 and US = 63), and circulating cytokines (NCT = 18; US = 73).

RESULTS

NCT was associated with increased infiltration of cytotoxic CD8⁺ T cells and CD20⁺ B cells. Moreover, NCT was associated with increases in CD8⁺CD103⁺ and CD4⁺CD103⁺PD-1⁺TIM3^- tissue resident memory T cells. Gene expression profiling supported memory function of CD8⁺ and CD4⁺ T cells. However, NCT did not affect T-cell receptor clonality, richness, or tumor mutational burden. Finally, NCT was associated with decreased plasma BDNF (TrkB) at baseline and week 4 after surgery.

CONCLUSIONS

Our study supports that, in the context of resectable NSCLC, neoadjuvant chemotherapy promotes antitumor immunity through T and B cell recruitment in the IME and through a phenotypic change toward cytotoxic and memory CD8⁺ and CD4⁺ memory helper T cells.

Abstract PR03: Nongenomic BAP1 aberrancy drives highly aggressive cutaneous melanoma phenotype

The purpose of this study was to determine the role of BAP1 levels in cutaneous melanoma (CM). BAP1 is a tumor suppressor in which loss of heterozygosity (LOH) from mutation and copy number alteration is well described in germline and somatic cancers. Although BAP1 genomic alterations in CM are extremely rare (2% of 665 samples from 5 datasets), marked variability in BAP1 expression is observed in CM. We show that low nuclear BAP1 levels portend a significantly worse clinical outcome in stage III CM (n=37, log rank p ≤0.01 for both overall survival and progression-free survival). Gene Set Enrichment Analysis (GSEA) revealed low BAP1 expression to be most highly ranked with an increased epithelial–mesenchymal transition (EMT) gene expression profile in CM tumors (n=379, FDR q = 1.34E-26) and cell lines (n=53, FDR q = 2.86E-116). We identify the expression of ZEB1, a master regulator of EMT, to be significantly associated with low BAP1 expression in CM tumors and cell lines (p= 1.5E-04 and 3.3E-05, respectively). Analysis of the BAP1 promoter indicates three canonical ZEB1 binding sites. Functional experiments show ZEB1 to bind to the BAP1 promoter, and luciferase activity assays indicate that ZEB1 acts as a transcriptional suppressor of BAP1 expression with differential utilization of the promoter binding sites. Targeted reduction of endogenous ZEB1 caused increased BAP1 levels, while targeted reduction of BAP1 did not modulate ZEB1 levels, consistent with ZEB1 having a suppressive effect on BAP1. Phenotypically, targeted reduction of BAP1 in CM cells resulted in a switch from a more differentiated, melanocytic state, to a less differentiated, more migratory and invasive phenotype. Extinguishing melanocyte-specific BAP1 in mice with a BRAF V600E mutant genetic background resulted in the emergence of primary melanoma tumors, with a marked EMT gene expression profile, and resultant metastases. Given the phenotypic changes associated with BAP1 levels in our mouse and human studies, we then tested the effect of modulating BAP1 on BRAF targeted therapy. Exogenous expression of BAP1 sensitized BRAF inhibitor (vemurafenib)-resistant melanoma cells, while the targeted reduction of BAP1 desensitized BRAF inhibitor-sensitive melanoma cells. BRAF mutant/BAP1 loss mice failed to exhibit a marked response to vemurafenib treatment compared to control mice. These data implicate regulation of BAP1 to be a major mechanism that characterizes a highly malignant and treatment-resistant subset of tumors. Our study indicates that nongenomic reduction in BAP1 through ZEB1 transcriptional modulation may be a key factor in aggressive CM.

This abstract is also being presented as Poster A30.

Citation Format: Haifeng Zhu, Junna Oba, Xiaoxing Yu, Caitlin A. Creasy, Marie-Andrée Forget, Fernando Carapeto, Cara L. Haymaker, Chang-Jiun Wu, Tatiana V. Karpinets, Wei-Lien Wang, Michael T. Tetzlaff, Alexander J. Lazar, Gordon B. Mills, Amanda R. Moore, Yu Chen, Jianhua Zhang, Jeffrey E. Gershenwald, Jennifer A. Wargo, Chantale Bernatchez, Patrick Hwu, P. Andrew Futreal, Scott E. Woodman. Nongenomic BAP1 aberrancy drives highly aggressive cutaneous melanoma phenotype [abstract]. In: Proceedings of the AACR Special Conference on Melanoma: From Biology to Target; 2019 Jan 15-18; Houston, TX. Philadelphia (PA): AACR; Cancer Res 2020;80(19 Suppl):Abstract nr PR03.

Abstract B11: Comparisons of microbial composition identified via 16s rRNA sequencing and whole-genome sequence-based analysis using gut samples for patients with cervical cancer

Introduction: 16S RNA sequencing and whole-genome sequencing (WGS) techniques have been previously shown to have a significant degree of correlation, particularly for higher-order-level taxa. However, the large majority of those studies utilize data derived from samples collected in similar but not identical contexts. This project provides a unique opportunity in that both 16S RNA and WGS sequencing datasets were derived from a single fecal sample originating from each patient while undergoing chemoradiation therapy. Using this high-quality information, we investigated the correlation of these two datasets in terms of microbial composition and taxa abundances.

Methods: Forty-one rectal swab samples were collected at the time of pelvic examination before the initiation of chemoradiation treatment and sequenced via Illumina platforms. Alpha and beta diversity was evaluated using PCoA ordination derived from several distance metrics (Bray-Curtis, Jaccard, Weighted and Unweighted UniFrac, and Euclidean). Diversity measures from different data analysis tools were also compared to highlight differences in data processing. Taxa with the highest relative abundances within each set of sequencing data were used to evaluate consensus between the datasets.

Results: Alpha diversity assessed by different measures from WGS analysis correlated (Spearman rho value) with the same measures derived from 16S rRNA sequencing. These measures include OTU abundances (rho = 0.763, p = 6.61e−09), Shannon diversity (rho = 0.809, p = 4.40e−09), and InvSimpson (rho = 0.797, p = 2.11e−08). Beta diversity analysis highlighted a greater level of similarity within WGS samples than 16sRNA samples and a high degree of dissimilarity between the two datasets. Among the top fourteen most abundant taxa identified in 16S rRNA sequencing, the four best correlated with WGS results were Faecalibacterium, rho = 0.475, p = 1.68e−03; Mollicutes, rho = 0.302, p = 5.48e−02; Clostridiales, rho = 0.021, p = 8.98e−01; and Bacteria, rho = 0.016, p = 9.23e−01. The remaining ten taxa were negatively correlated.

Conclusions: The use of whole-genome sequencing approaches in microbiome analysis is expanding in order to provide greater depth of information about microbial function. However, comparisons of results from 16S rRNA sequencing and WGS will need to be evaluated with caution due to differences in the interpretation of relative abundances of specific taxa. Future efforts may help to better harmonize these analyses to facilitate comparisons across platforms.

Citation Format: Greyson Biegert, Xiaogang Wu, Tatiana V. Karpinets, Melissa Mezzari, Jianhua Zhang, Lauren Colbert, Ann Klopp. Comparisons of microbial composition identified via 16s rRNA sequencing and whole-genome sequence-based analysis using gut samples for patients with cervical cancer [abstract]. In: Proceedings of the AACR Special Conference on the Microbiome, Viruses, and Cancer; 2020 Feb 21-24; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2020;80(8 Suppl):Abstract nr B11.

Microbial Diversity and Composition Is Associated with Patient-Reported Toxicity during Chemoradiation Therapy for Cervical Cancer

PURPOSE

Patients receiving pelvic radiation for cervical cancer experience high rates of acute gastrointestinal (GI) toxicity. The association of changes in the gut microbiome with bowel toxicity from radiation is not well characterized.

METHODS AND MATERIALS

Thirty-five patients undergoing definitive chemoradiation therapy (CRT) underwent longitudinal sampling (baseline and weeks 1, 3, and 5) of the gut microbiome and prospective assessment of patient-reported GI toxicity. DNA was isolated from stool obtained at rectal examination and analyzed with 16S rRNA sequencing. GI toxicity was assessed with the Expanded Prostate Cancer Index Composite instrument to evaluate frequency, urgency, and discomfort associated with bowel function. Shannon diversity index was used to characterize alpha (within sample) diversity. Weighted UniFrac principle coordinates analysis was used to compare beta (between sample) diversity between samples using permutational multivariate analysis of variance. Linear discriminant analysis effect size highlighted microbial features that best distinguish categorized patient samples.

RESULTS

Gut microbiome diversity continuously decreased over the course of CRT, with the largest decrease at week 5. Expanded Prostate Cancer Index Composite bowel function scores also declined over the course of treatment, reflecting increased symptom burden. At all individual time points, higher diversity of the gut microbiome was linearly correlated with better patient-reported GI function, but baseline diversity was not predictive of eventual outcome. Patients with high toxicity demonstrated different compositional changes during CRT in addition to compositional differences in Clostridia species.

CONCLUSIONS

Over time, increased radiation toxicity is associated with decreased gut microbiome diversity. Baseline diversity is not predictive of end-of-treatment bowel toxicity, but composition may identify patients at risk for developing high toxicity.

Linking Associations of Rare Low-Abundance Species to Their Environments by Association Networks

Studies of microbial communities by targeted sequencing of rRNA genes lead to recovering numerous rare low-abundance taxa with unknown biological roles. We propose to study associations of such rare organisms with their environments by a computational framework based on transformation of the data into qualitative variables. Namely, we analyze the sparse table of putative species or OTUs (operational taxonomic units) and samples generated in such studies, also known as an OTU table, by collecting statistics on co-occurrences of the species and on shared species richness across samples. Based on the statistics we built two association networks, of the rare putative species and of the samples respectively, using a known computational technique, Association networks (Anets) developed for analysis of qualitative data. Clusters of samples and clusters of OTUs are then integrated and combined with metadata of the study to produce a map of associated putative species in their environments. We tested and validated the framework on two types of microbiomes, of human body sites and that of the Populus tree root systems. We show that in both studies the associations of OTUs can separate samples according to environmental or physiological characteristics of the studied systems.

Gut microbiome modulates response to anti-PD-1 immunotherapy in melanoma patients

Preclinical mouse models suggest that the gut microbiome modulates tumor response to checkpoint blockade immunotherapy; however, this has not been well-characterized in human cancer patients. Here we examined the oral and gut microbiome of melanoma patients undergoing anti-programmed cell death 1 protein (PD-1) immunotherapy (n = 112). Significant differences were observed in the diversity and composition of the patient gut microbiome of responders versus nonresponders. Analysis of patient fecal microbiome samples (n = 43, 30 responders, 13 nonresponders) showed significantly higher alpha diversity (P < 0.01) and relative abundance of bacteria of the Ruminococcaceae family (P < 0.01) in responding patients. Metagenomic studies revealed functional differences in gut bacteria in responders, including enrichment of anabolic pathways. Immune profiling suggested enhanced systemic and antitumor immunity in responding patients with a favorable gut microbiome as well as in germ-free mice receiving fecal transplants from responding patients. Together, these data have important implications for the treatment of melanoma patients with immune checkpoint inhibitors.

Abstract 2672: Response to anti-PD-1 based therapy in metastatic melanoma patients is associated with the diversity and composition of the gut microbiome

Background: Melanoma therapy has benefitted greatly from immune checkpoint blockade, although responses are variable and not always durable. There is a growing appreciation of the role of the microbiome in cancer-related outcomes and recent evidence in murine models suggests that modulation of the gut microbiome may enhance responses to immune checkpoint blockade in melanoma. However this has not been investigated in patients. Here, we demonstrate that differential bacterial “signatures” exist in the gut microbiome of responders (R) and non-responders (NR) to anti-PD-1 therapy, and that insights gained could be used to derive actionable strategies to enhance responses.

Methods: We collected buccal (n=105) and stool (n=53) samples from a cohort of anti-PD-1 treated metastatic melanoma patients (n=110). Patients were classified as either R or NR based on RECIST criteria, and 16S rDNA, and whole-genome shotgun sequencing was performed to characterize the diversity, composition and functional capabilities of the microbiomes. Immune profiling (via 7-marker IHC panel of CD3, CD8, PD-1, PD-L1, Granzyme B, RORγT and FoxP3) and cytokine analyses were also performed on available tumors and serum samples at baseline.

Results: In these studies, we observed significant differences in the diversity and composition of the gut microbiome in R versus NR to PD-1 blockade at baseline, but no clear differences in buccal microbiomes. Specifically, R had a significantly higher alpha diversity compared to NR (p=0.017), and the Ruminococcaceae family of the Clostridiales order was enriched in R whereas Prevotellaceae family of the Bacteroidales order was enriched in NR. Immune profiling demonstrated significantly increased immune infiltrates in baseline tumor samples of R, with a positive correlation between CD8, CD3, PD-1 and FoxP3 T-cell density and abundance of specific bacteria enriched in R (e.g. Faecalibacterium). Low diversity was also associated with elevated levels of chronic inflammation markers in the serum at baseline. Lastly, we saw differentially abundant metabolic pathways in the gut microbiomes of R (pyrimidine nucleotide biosynthesis, fatty acid biosynthesis, shikimate pathway) vs NR (Tricarboxylic acid cycle, assimilatory sulphate and nitrate reduction, tryptophan biosynthesis).

Conclusion: Differences exist in the diversity and composition of the gut microbiome in R vs NR to anti-PD-1 therapy and these microbiota could bridge the gap between host metabolism and anti-tumor immunity. These results have far-reaching implications and suggest that modifications to the gut microbiome could potentially enhance therapeutic responses to immune checkpoint blockade.

Citation Format: Vancheswaran Gopalakrishnan, Christine Spencer, Alexandre Reuben, Peter Prieto, Diego Vicente, Tatiana V. Karpinets, Courtney W. Hudgens, Diane S. Hutchinson, Michael Tetzlaff, Alexander Lazar, Michael A. Davies, Jeffrey E. Gershenwald, Robert Jenq, Patrick Hwu, Padmanee Sharma, James Allison, Andrew Futreal, Nadim Ajami, Joseph Petrosino, Carrie Daniel-MacDougall, Jennifer A. Wargo. Response to anti-PD-1 based therapy in metastatic melanoma patients is associated with the diversity and composition of the gut microbiome [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2017; 2017 Apr 1-5; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2017;77(13 Suppl):Abstract nr 2672. doi:10.1158/1538-7445.AM2017-2672

Oncogenic Kras drives invasion and maintains metastases in colorectal cancer

Human colorectal cancer (CRC) is a major cause of cancer mortality and frequently harbors activating mutations in the KRAS gene. To understand the role of oncogenic KRAS in CRC, we engineered a mouse model of metastatic CRC that harbors an inducible oncogenic Kras allele (Kras^mut ) and conditional null alleles of Apc and Trp53 (iKAP). The iKAP model recapitulates tumor progression from adenoma through metastases. Whole-exome sequencing revealed that the Kras^mut allele was heterogenous in primary tumors yet homogenous in metastases, a pattern consistent with activated Kras^mut signaling being a driver of progression to metastasis. System-level and functional analyses revealed the TGF-β pathway as a key mediator of Kras^mut -driven invasiveness. Genetic extinction of Kras^mut resulted in specific elimination of the Kras^mut subpopulation in primary and metastatic tumors, leading to apoptotic elimination of advanced invasive and metastatic disease. This faithful CRC model provides genetic evidence that Kras^mut drives CRC invasion and maintenance of metastases.

Novel algorithmic approach predicts tumor mutation load and correlates with immunotherapy clinical outcomes using a defined gene mutation set

BACKGROUND

While clinical outcomes following immunotherapy have shown an association with tumor mutation load using whole exome sequencing (WES), its clinical applicability is currently limited by cost and bioinformatics requirements.

METHODS

We developed a method to accurately derive the predicted total mutation load (PTML) within individual tumors from a small set of genes that can be used in clinical next generation sequencing (NGS) panels. PTML was derived from the actual total mutation load (ATML) of 575 distinct melanoma and lung cancer samples and validated using independent melanoma (n = 312) and lung cancer (n = 217) cohorts. The correlation of PTML status with clinical outcome, following distinct immunotherapies, was assessed using the Kaplan-Meier method.

RESULTS

PTML (derived from 170 genes) was highly correlated with ATML in cutaneous melanoma and lung adenocarcinoma validation cohorts (R2 = 0.73 and R2 = 0.82, respectively). PTML was strongly associated with clinical outcome to ipilimumab (anti-CTLA-4, three cohorts) and adoptive T-cell therapy (1 cohort) clinical outcome in melanoma. Clinical benefit from pembrolizumab (anti-PD-1) in lung cancer was also shown to significantly correlate with PTML status (log rank P value < 0.05 in all cohorts).

CONCLUSIONS

The approach of using small NGS gene panels, already applied to guide employment of targeted therapies, may have utility in the personalized use of immunotherapy in cancer.

Metabolic environments and genomic features associated with pathogenic and mutualistic interactions between bacteria and plants

Genomic characteristics discriminating parasitic and mutualistic relationship of bacterial symbionts with plants are poorly understood. This study comparatively analyzed the genomes of 54 mutualists and pathogens to discover genomic markers associated with the different phenotypes. Using metabolic network models, we predict external environments associated with free-living and symbiotic lifestyles and quantify dependences of symbionts on the host in terms of the consumed metabolites. We show that specific differences between the phenotypes are pronounced at the levels of metabolic enzymes, especially carbohydrate active, and protein functions. Overall, biosynthetic functions are enriched and more diverse in plant mutualists whereas processes and functions involved in degradation and host invasion are enriched and more diverse in pathogens. A distinctive characteristic of plant pathogens is a putative novel secretion system with a circadian rhythm regulator. A specific marker of plant mutualists is the co-residence of genes encoding nitrogenase and ribulose bisphosphate carboxylase/oxygenase (RuBisCO). We predict that RuBisCO is likely used in a putative metabolic pathway to supplement carbon obtained heterotrophically with low-cost assimilation of carbon from CO2. We validate results of the comparative analysis by predicting correct phenotype, pathogenic or mutualistic, for 20 symbionts in an independent set of 30 pathogens, mutualists, and commensals.

Binding Motifs in Bacterial Gene Promoters Modulate Transcriptional Effects of Global Regulators CRP and ArcA

Bacterial gene regulation involves transcription factors (TF) that bind to DNA recognition sequences in operon promoters. These recognition sequences, many of which are palindromic, are known as regulatory elements or transcription factor binding sites (TFBS). Some TFs are global regulators that can modulate the expression of hundreds of genes. In this study we examine global regulator half-sites, where a half-site, which we shall call a binding motif (BM), is one half of a palindromic TFBS. We explore the hypothesis that the number of BMs plays an important role in transcriptional regulation, examining empirical data from transcriptional profiling of the CRP and ArcA regulons. We compare the power of BM counts and of full TFBS characteristics to predict induced transcriptional activity. We find that CRP BM counts have a nonlinear effect on CRP-dependent transcriptional activity and predict this activity better than full TFBS quality or location.

Analyzing large biological datasets with association networks

Due to advances in high-throughput biotechnologies biological information is being collected in databases at an amazing rate, requiring novel computational approaches that process collected data into new knowledge in a timely manner. In this study, we propose a computational framework for discovering modular structure, relationships and regularities in complex data. The framework utilizes a semantic-preserving vocabulary to convert records of biological annotations of an object, such as an organism, gene, chemical or sequence, into networks (Anets) of the associated annotations. An association between a pair of annotations in an Anet is determined by the similarity of their co-occurrence pattern with all other annotations in the data. This feature captures associations between annotations that do not necessarily co-occur with each other and facilitates discovery of the most significant relationships in the collected data through clustering and visualization of the Anet. To demonstrate this approach, we applied the framework to the analysis of metadata from the Genomes OnLine Database and produced a biological map of sequenced prokaryotic organisms with three major clusters of metadata that represent pathogens, environmental isolates and plant symbionts.

Enhancing a Pathway-Genome Database (PGDB) to capture subcellular localization of metabolites and enzymes: the nucleotide-sugar biosynthetic pathways of Populus trichocarpa

Understanding how cellular metabolism works and is regulated requires that the underlying biochemical pathways be adequately represented and integrated with large metabolomic data sets to establish a robust network model. Genetically engineering energy crops to be less recalcitrant to saccharification requires detailed knowledge of plant polysaccharide structures and a thorough understanding of the metabolic pathways involved in forming and regulating cell-wall synthesis. Nucleotide-sugars are building blocks for synthesis of cell wall polysaccharides. The biosynthesis of nucleotide-sugars is catalyzed by a multitude of enzymes that reside in different subcellular organelles, and precise representation of these pathways requires accurate capture of this biological compartmentalization. The lack of simple localization cues in genomic sequence data and annotations however leads to missing compartmentalization information for eukaryotes in automatically generated databases, such as the Pathway-Genome Databases (PGDBs) of the SRI Pathway Tools software that drives much biochemical knowledge representation on the internet. In this report, we provide an informal mechanism using the existing Pathway Tools framework to integrate protein and metabolite sub-cellular localization data with the existing representation of the nucleotide-sugar metabolic pathways in a prototype PGDB for Populus trichocarpa. The enhanced pathway representations have been successfully used to map SNP abundance data to individual nucleotide-sugar biosynthetic genes in the PGDB. The manually curated pathway representations are more conducive to the construction of a computational platform that will allow the simulation of natural and engineered nucleotide-sugar precursor fluxes into specific recalcitrant polysaccharide(s).

Database URL: The curated Populus PGDB is available in the BESC public portal at http://cricket.ornl.gov/cgi-bin/beocyc_home.cgi and the nucleotide-sugar biosynthetic pathways can be directly accessed at http://cricket.ornl.gov:1555/PTR/new-image?object=SUGAR-NUCLEOTIDES.

BESC knowledgebase public portal

The BioEnergy Science Center (BESC) is undertaking large experimental campaigns to understand the biosynthesis and biodegradation of biomass and to develop biofuel solutions. BESC is generating large volumes of diverse data, including genome sequences, omics data and assay results. The purpose of the BESC Knowledgebase is to serve as a centralized repository for experimentally generated data and to provide an integrated, interactive and user-friendly analysis framework. The Portal makes available tools for visualization, integration and analysis of data either produced by BESC or obtained from external resources.

AVAILABILITY

http://besckb.ornl.gov.

Mutant alcohol dehydrogenase leads to improved ethanol tolerance in Clostridium thermocellum

Clostridium thermocellum is a thermophilic, obligately anaerobic, gram-positive bacterium that is a candidate microorganism for converting cellulosic biomass into ethanol through consolidated bioprocessing. Ethanol intolerance is an important metric in terms of process economics, and tolerance has often been described as a complex and likely multigenic trait for which complex gene interactions come into play. Here, we resequence the genome of an ethanol-tolerant mutant, show that the tolerant phenotype is primarily due to a mutated bifunctional acetaldehyde-CoA/alcohol dehydrogenase gene (adhE), hypothesize based on structural analysis that cofactor specificity may be affected, and confirm this hypothesis using enzyme assays. Biochemical assays confirm a complete loss of NADH-dependent activity with concomitant acquisition of NADPH-dependent activity, which likely affects electron flow in the mutant. The simplicity of the genetic basis for the ethanol-tolerant phenotype observed here informs rational engineering of mutant microbial strains for cellulosic ethanol production.

CAZymes Analysis Toolkit (CAT): web service for searching and analyzing carbohydrate-active enzymes in a newly sequenced organism using CAZy database

The Carbohydrate-Active Enzyme (CAZy) database provides a rich set of manually annotated enzymes that degrade, modify, or create glycosidic bonds. Despite rich and invaluable information stored in the database, software tools utilizing this information for annotation of newly sequenced genomes by CAZy families are limited. We have employed two annotation approaches to fill the gap between manually curated high-quality protein sequences collected in the CAZy database and the growing number of other protein sequences produced by genome or metagenome sequencing projects. The first approach is based on a similarity search against the entire nonredundant sequences of the CAZy database. The second approach performs annotation using links or correspondences between the CAZy families and protein family domains. The links were discovered using the association rule learning algorithm applied to sequences from the CAZy database. The approaches complement each other and in combination achieved high specificity and sensitivity when cross-evaluated with the manually curated genomes of Clostridium thermocellum ATCC 27405 and Saccharophagus degradans 2-40. The capability of the proposed framework to predict the function of unknown protein domains and of hypothetical proteins in the genome of Neurospora crassa is demonstrated. The framework is implemented as a Web service, the CAZymes Analysis Toolkit, and is available at http://cricket.ornl.gov/cgi-bin/cat.cgi.

Shewanella knowledgebase: integration of the experimental data and computational predictions suggests a biological role for transcription of intergenic regions

Shewanellae are facultative γ-proteobacteria whose remarkable respiratory versatility has resulted in interest in their utility for bioremediation of heavy metals and radionuclides and for energy generation in microbial fuel cells. Extensive experimental efforts over the last several years and the availability of 21 sequenced Shewanella genomes made it possible to collect and integrate a wealth of information on the genus into one public resource providing new avenues for making biological discoveries and for developing a system level understanding of the cellular processes. The Shewanella knowledgebase was established in 2005 to provide a framework for integrated genome-based studies on Shewanella ecophysiology. The present version of the knowledgebase provides access to a diverse set of experimental and genomic data along with tools for curation of genome annotations and visualization and integration of genomic data with experimental data. As a demonstration of the utility of this resource, we examined a single microarray data set from Shewanella oneidensis MR-1 for new insights into regulatory processes. The integrated analysis of the data predicted a new type of bacterial transcriptional regulation involving co-transcription of the intergenic region with the downstream gene and suggested a biological role for co-transcription that likely prevents the binding of a regulator of the upstream gene to the regulator binding site located in the intergenic region.

Database URL:http://shewanella-knowledgebase.org:8080/Shewanella/ or http://spruce.ornl.gov:8080/Shewanella/

Shewregdb: Database and visualization environment for experimental and predicted regulatory information in Shewanella oneidensis mr-1

Shewanella oneidensis MR-1 is an important model organism for environmental research as it has an exceptional metabolic and respiratory versatility regulated by a complex regulatory network. We have developed a database to collect experimental and computational data relating to regulation of gene and protein expression, and, a visualization environment that enables integration of these data types. The regulatory information in the database includes predictions of DNA regulator binding sites, sigma factor binding sites, transcription units, operons, promoters, and RNA regulators including non-coding RNAs, riboswitches, and different types of terminators.

Conserved synteny at the protein family level reveals genes underlying Shewanella species' cold tolerance and predicts their novel phenotypes

Bacteria of the genus Shewanella can thrive in different environments and demonstrate significant variability in their metabolic and ecophysiological capabilities including cold and salt tolerance. Genomic characteristics underlying this variability across species are largely unknown. In this study, we address the problem by a comparison of the physiological, metabolic, and genomic characteristics of 19 sequenced Shewanella species. We have employed two novel approaches based on association of a phenotypic trait with the number of the trait-specific protein families (Pfam domains) and on the conservation of synteny (order in the genome) of the trait-related genes. Our first approach is top-down and involves experimental evaluation and quantification of the species’ cold tolerance followed by identification of the correlated Pfam domains and genes with a conserved synteny. The second, a bottom-up approach, predicts novel phenotypes of the species by calculating profiles of each Pfam domain among their genomes and following pair-wise correlation of the profiles and their network clustering. Using the first approach, we find a link between cold and salt tolerance of the species and the presence in the genome of a Na⁺/H⁺ antiporter gene cluster. Other cold-tolerance-related genes include peptidases, chemotaxis sensory transducer proteins, a cysteine exporter, and helicases. Using the bottom-up approach, we found several novel phenotypes in the newly sequenced Shewanella species, including degradation of aromatic compounds by an aerobic hybrid pathway in Shewanella woodyi, degradation of ethanolamine by Shewanella benthica, and propanediol degradation by Shewanella putrefaciens CN32 and Shewanella sp. W3-18-1.

Phenotype fingerprinting suggests the involvement of single-genotype consortia in degradation of aromatic compounds by Rhodopseudomonas palustris

Anaerobic degradation of complex organic compounds by microorganisms is crucial for development of innovative biotechnologies for bioethanol production and for efficient degradation of environmental pollutants. In natural environments, the degradation is usually accomplished by syntrophic consortia comprised of different bacterial species. This strategy allows consortium organisms to reduce efforts required for maintenance of the redox homeostasis at each syntrophic level. Cellular mechanisms that maintain the redox homeostasis during the degradation of aromatic compounds by one organism are not fully understood. Here we present a hypothesis that the metabolically versatile phototrophic bacterium Rhodopseudomonas palustris forms its own syntrophic consortia, when it grows anaerobically on p-coumarate or benzoate as a sole carbon source. We have revealed the consortia from large-scale measurements of mRNA and protein expressions under p-coumarate, benzoate and succinate degrading conditions using a novel computational approach referred as phenotype fingerprinting. In this approach, marker genes for known R. palustris phenotypes are employed to determine the relative expression levels of genes and proteins in aromatics versus non-aromatics degrading condition. Subpopulations of the consortia are inferred from the expression of phenotypes and known metabolic modes of the R. palustris growth. We find that p-coumarate degrading conditions may lead to at least three R. palustris subpopulations utilizing p-coumarate, benzoate, and CO₂ and H₂. Benzoate degrading conditions may also produce at least three subpopulations utilizing benzoate, CO₂ and H₂, and N₂ and formate. Communication among syntrophs and inter-syntrophic dynamics in each consortium are indicated by up-regulation of transporters and genes involved in the curli formation and chemotaxis. The N₂-fixing subpopulation in the benzoate degrading consortium has preferential activation of the vanadium nitrogenase over the molybdenum nitrogenase. This subpopulation in the consortium was confirmed in an independent experiment by consumption of dissolved nitrogen gas under the benzoate degrading conditions.

A unifying concept for the dependence of whole-crop N : P ratio on biomass: theory and experiment

BACKGROUND AND AIMS

Numerous estimates have been made of the concentrations of N and P required for good growth of crop species but they have not been defined by any unifying model. The aim of the present study was to develop such a model for the dependence of the N : P ratio on crop mass, to test its validity and to use it to identify elements of similarity between different crop species and wild plants.

METHODS

A model was derived between plant N : P ratio (Rw) and its dry biomass per unit area (W) during growth with near optimum nutrition by considering that plants consist of growth-related tissue and storage-related tissue with N : P ratios Rg and Rs, respectively. Testing and calibration against experimental data on different crop species led to a simple equation between Rw and W which was tested against independent experimental data.

KEY RESULTS

The validity of the model and equation was supported by 365 measurements of Rw in 38 field experiments on crops. Rg and Rs remained approximately constant throughout growth, with average values of 11.8 and 5.8 by mass. The model also approximately predicted the relationships between leaf N and P concentrations in 124 advisory estimates on immature tissues and in 385 wild species from published global surveys.

CONCLUSIONS

The N : P ratio of the biomass of very different crops, during growth with near optimum levels of nutrients, is defined entirely in terms of crop biomass, an average N : P ratio of the storage/structure-related tissue of the crop and an average N : P ratio of the growth-related tissue. The latter is similar to that found in leaves of many wild plant species, and even micro-organisms and terrestrial and freshwater autotrophs.

RNA:protein ratio of the unicellular organism as a characteristic of phosphorous and nitrogen stoichiometry and of the cellular requirement of ribosomes for protein synthesis

BACKGROUND

Mean phosphorous:nitrogen (P:N) ratios and relationships of P:N ratios with the growth rate of organisms indicate a surprising similarity among and within microbial species, plants, and insect herbivores. To reveal the cellular mechanisms underling this similarity, the macromolecular composition of seven microorganisms and the effect of specific growth rate (SGR) on RNA:protein ratio, the number of ribosomes, and peptide elongation rate (PER) were analyzed under different conditions of exponential growth.

RESULTS

It was found that P:N ratios calculated from RNA and protein contents in these particular organisms were in the same range as the mean ratios reported for diverse organisms and had similar positive relationships with growth rate, consistent with the growth-rate hypothesis. The efficiency of protein synthesis in microorganisms is estimated as the number of active ribosomes required for the incorporation of one amino acid into the synthesized protein. This parameter is calculated as the SGR:PER ratio. Experimental and theoretical evidence indicated that the requirement of ribosomes for protein synthesis is proportional to the RNA:protein ratio. The constant of proportionality had the same values for all organisms, and was derived mechanistically from the characteristics of the protein-synthesis machinery of the cell (the number of nucleotides per ribosome, the average masses of nucleotides and amino acids, the fraction of ribosomal RNA in the total RNA, and the fraction of active ribosomes). Impairment of the growth conditions decreased the RNA:protein ratio and increased the overall efficiency of protein synthesis in the microorganisms.

CONCLUSION

Our results suggest that the decrease in RNA:protein and estimated P:N ratios with decrease in the growth rate of the microorganism is a consequence of an increased overall efficiency of protein synthesis in the cell resulting from activation of the general stress response and increased transcription of cellular maintenance genes at the expense of growth related genes. The strong link between P:N stoichiometry, RNA:protein ratio, ribosomal requirement for protein synthesis, and growth rate of microorganisms indicated by the study could be used to characterize the N and P economy of complex ecosystems such as soils and the oceans.

Histone H3 lysine 9 and H4 lysine 20 trimethylation and the expression of Suv4-20h2 and Suv-39h1 histone methyltransferases in hepatocarcinogenesis induced by methyl deficiency in rats

The field of cancer epigenetics has received much attention in recent years. However, the relationship of cancer epigenetics with cancer etiology is not clear. Recent studies suggest the involvement of altered DNA methylation and histone modifications in the emergence of epigenetically reprogrammed cells with specific tumor-related phenotypes at premalignant stages of tumor development. In this study, we used a methyl-deficient model of rodent hepatocarcinogenesis to examine the roles of DNA, histone H3 lysine 9 and histone H4 lysine 20 methylation, and the level of the expression of Suv39h1 and Suv4-20h2 histone methyltransferases in the carcinogenic process. We demonstrated that the development of liver tumors was characterized by progressive demethylation of DNA repeats, decrease in histone H4 lysine 20 trimethylation, and a gradual decrease in the expression of Suv4-20h2 histone methyltransferase. A prominent increase in the trimethylation of histone H3 lysine 9 and in the expression of Suv39h1 histone methyltransferase was observed in preneoplastic nodules and liver tumors indicating the promotional role of these epigenetic alterations at later stages of carcinogenesis. The appearance of tumor-specific epigenetic alterations (demethylation of repetitive elements, loss of histone H4 lysine 20 trimethylation, altered expression of Suv4-20h2 and Suv39h1 histone methyltransferases) at preneoplastic stages of hepatocarcinogenesis provides experimental support for the epigenetic hypothesis of tumorigenesis that considers stress-induced epigenetic reprogramming of the cell as an important prerequisite to succeeding mutations.

Tumorigenesis: the adaptation of mammalian cells to sustained stress environment by epigenetic alterations and succeeding matched mutations

Recent studies indicate that during tumorigenic transformations, cells may generate mutations by themselves as a result of error-prone cell division with participation of error-prone polymerases and aberrant mitosis. These mechanisms may be activated in cells by continuing proliferative and survival signaling in a sustained stress environment (SSE). The paper hypothesizes that long-term exposure to this signaling epigenetically reprograms the genome of some cells and, in addition, leads to their senescence. The epigenetic reprogramming results in: (i) hypermethylation of tumor-suppressor genes involved in the onset of cell-cycle arrest, apoptosis and DNA repair; (ii) hypomethylation of proto-oncogenes associated with persistent proliferative activity; and (iii) the global demethylation of the genome and activation of DNA repeats. These epigenetic changes in the proliferating cells associate with their replicative senescence and allow the reprogrammed senescent cells to overcome the cell-cycle arrest and to activate error-prone replications. It is hypothesized that the generation of mutations in the error-prone replications of the epigenetically reprogrammed cells is not random. The mutations match epigenetic alterations in the cellular genome, namely gain of function mutations in the case of hypomethylation and loss of functions in the case of hypermethylation. In addition, continuing proliferation of the cells imposed by signaling in SSE speeds up the natural selection of the mutant cells favoring the survival of the cells with mutations that are beneficial in the environment. In this way, a stress-induced replication of the cells epigenetically reprograms their genome for quick adaptation to stressful environments providing an increased rate of mutations, epigenetic tags to beneficial mutations and quick selection process. In combination, these processes drive the origin of the transformed mammalian cells, cancer development and progression. Support from genomic, biochemical and medical studies for the proposed hypothesis, and its implementations are discussed.

Model of the developing tumorigenic phenotype in mammalian cells and the roles of sustained stress and replicative senescence

The molecular mechanisms that drive mammalian cells to the development of cancer are the subject of intense biochemical, genetic and medical studies. But for the present, there is no comprehensive model that might serve as a general framework for the interpretation of experimental data. This paper is an attempt to create a conceptual model of the mechanism of the developing tumorigenic phenotype in mammalian cells, defined as having high genomic instability and proliferative activity. The basic statement in the model is that mutations acquired by tumor cells are not caused directly by external DNA damaging agents, but instead are produced by the cell itself as an output of a Mutator Response similar to the bacterial "SOS response" and characterized by the initiation of error-prone cell cycle progression and an elevated rate of mutation. This response may be induced in arrested mammalian cells by intracellular and extracellular proliferative signals combined with blocked apoptosis. The mutant cells originated by this response are subjected to natural selection via apoptosis and turnover. This selection process favors the survival of cells with high proliferative activity and the suppression of apoptosis resulting in the long run in the appearance of immortalized cells with high proliferative activity. Either a sustained stressful environment accompanied by continuing apoptotic cell death, or replicative senescence, provides conditions suitable for activation of the Mutator Response, namely the emergence of arrested cells with blocked apoptosis and the induction of proliferative signal. It also accelerates the selection process by providing continuing cell turnover. The proposed mechanism is described at the level of involved metabolic pathways and proteins and substantiated by the related experimental data available in the literature.

Tailored gene array databases: applications in mechanistic toxicology

MOTIVATION

The development of an annotated global database suitable for a wide range of investigations is a challenging and labor-intensive task. Thus, the development of databases tailored for specific applications remains necessary. For example, in the field of toxicology, no annotated gene array databases are now available that may assist in the correlation of changes in gene activity to cellular functions and processes associated with the toxic response.

RESULTS

As an example of a tailored annotated database, an attempt was made to systematize available biological information on genes present on the Affymetrix Rat Toxicology U34 GeneChip, with a focus on how the gene products relate to liver cells and their response to chemical toxins. The information collected was imbedded in a local relational database to analyze data obtained in toxicological gene array experiments with hydrazine-exposed hepatocytes. The advantages and benefits of the tailored database in the biological interpretation of the results are demonstrated.