Exploring the gut DNA virome in fecal immunochemical test stool samples reveals associations with lifestyle in a large population-based study

Stool samples for fecal immunochemical tests (FIT) are collected in large numbers worldwide as part of colorectal cancer screening programs. Employing FIT samples from 1034 CRCbiome participants, recruited from a Norwegian colorectal cancer screening study, we identify, annotate and characterize more than 18000 DNA viruses, using shotgun metagenome sequencing. Only six percent of them are assigned to a known taxonomic family, with Microviridae being the most prevalent viral family. Linking individual profiles to comprehensive lifestyle and demographic data shows 17/25 of the variables to be associated with the gut virome. Physical activity, smoking, and dietary fiber consumption exhibit strong and consistent associations with both diversity and relative abundance of individual viruses, as well as with enrichment for auxiliary metabolic genes. We demonstrate the suitability of FIT samples for virome analysis, opening an opportunity for large-scale studies of this enigmatic part of the gut microbiome. The diverse viral populations and their connections to the individual lifestyle uncovered herein paves the way for further exploration of the role of the gut virome in health and disease.

Profiling small RNAs in fecal immunochemical tests: is it possible?

Fecal microRNAs represent promising molecules with potential clinical interest as non-invasive diagnostic and prognostic biomarkers. Colorectal cancer (CRC) screening based on the fecal immunochemical test (FIT) is an effective tool for prevention of cancer development. However, due to the poor sensitivity of FIT especially for premalignant lesions, there is a need for implementation of complementary tests. Improving the identification of individuals who would benefit from further investigation with colonoscopy using molecular analysis, such as miRNA profiling of FIT samples, would be ideal due to their widespread use. In the present study, we assessed the feasibility of applying small RNA sequencing to measure human miRNAs in FIT leftover buffer in samples from two European screening populations. We showed robust detection of miRNAs with profiles similar to those obtained from specimens sampled using the established protocol of RNA stabilizing buffers, or in long-term archived samples. Detected miRNAs exhibited differential abundances for CRC, advanced adenoma, and control samples that were consistent for FIT and RNA-stabilizing buffers. Interestingly, the sequencing data also allowed for concomitant evaluation of small RNA-based microbial profiles. We demonstrated that it is possible to explore the human miRNome in FIT leftover samples across populations and envision that the analysis of small RNA biomarkers can complement the FIT in large scale screening settings.

Associations of red and processed meat intake with screen-detected colorectal lesions

Limited data exist regarding the role of meat consumption in early-stage colorectal carcinogenesis. We examined associations of red and processed meat intake with screen-detected colorectal lesions in immunochemical fecal occult blood test (FIT)-positive participants, enrolled in the Norwegian CRCbiome study during 2017-2021, aged 55-77 years. Absolute and energy-adjusted intakes of red and processed meat (combined and individually) were assessed using a validated, semi-quantitative FFQ. Associations between meat intake and screen-detected colorectal lesions were examined using multinomial logistic regression analyses with adjustment for key covariates. Of 1162 participants, 319 presented with advanced colorectal lesions at colonoscopy. High v. low energy-adjusted intakes of red and processed meat combined, as well as red meat alone, were borderline to significantly positively associated with advanced colorectal lesions (OR of 1·24 (95 % CI 0·98, 1·57) and 1·34 (95 % CI 1·07, 1·69), respectively). A significant dose-response relationship was also observed for absolute intake levels (OR of 1·32 (95 % CI 1·09, 1·60) per 100 g/d increase in red and processed meat). For processed meat, no association was observed between energy-adjusted intakes and advanced colorectal lesions. A significant positive association was, however, observed for participants with absolute intake levels ≥ 100 v. < 50 g/d (OR of 1·19 (95 % CI 1·09, 1·31)). In summary, high intakes of red and processed meat were associated with presence of advanced colorectal lesions at colonoscopy in FIT-positive participants. The study demonstrates a potential role of dietary data to improve the performance of FIT-based screening.

The CRCbiome study: a large prospective cohort study examining the role of lifestyle and the gut microbiome in colorectal cancer screening participants

BACKGROUND

Colorectal cancer (CRC) screening reduces CRC incidence and mortality. However, current screening methods are either hampered by invasiveness or suboptimal performance, limiting their effectiveness as primary screening methods. To aid in the development of a non-invasive screening test with improved sensitivity and specificity, we have initiated a prospective biomarker study (CRCbiome), nested within a large randomized CRC screening trial in Norway. We aim to develop a microbiome-based classification algorithm to identify advanced colorectal lesions in screening participants testing positive for an immunochemical fecal occult blood test (FIT). We will also examine interactions with host factors, diet, lifestyle and prescription drugs. The prospective nature of the study also enables the analysis of changes in the gut microbiome following the removal of precancerous lesions.

METHODS

The CRCbiome study recruits participants enrolled in the Bowel Cancer Screening in Norway (BCSN) study, a randomized trial initiated in 2012 comparing once-only sigmoidoscopy to repeated biennial FIT, where women and men aged 50-74 years at study entry are invited to participate. Since 2017, participants randomized to FIT screening with a positive test result have been invited to join the CRCbiome study. Self-reported diet, lifestyle and demographic data are collected prior to colonoscopy after the positive FIT-test (baseline). Screening data, including colonoscopy findings are obtained from the BCSN database. Fecal samples for gut microbiome analyses are collected both before and 2 and 12 months after colonoscopy. Samples are analyzed using metagenome sequencing, with taxonomy profiles, and gene and pathway content as primary measures. CRCbiome data will also be linked to national registries to obtain information on prescription histories and cancer relevant outcomes occurring during the 10 year follow-up period.

DISCUSSION

The CRCbiome study will increase our understanding of how the gut microbiome, in combination with lifestyle and environmental factors, influences the early stages of colorectal carcinogenesis. This knowledge will be crucial to develop microbiome-based screening tools for CRC. By evaluating biomarker performance in a screening setting, using samples from the target population, the generalizability of the findings to future screening cohorts is likely to be high.

TRIAL REGISTRATION

ClinicalTrials.gov Identifier: NCT01538550 .

Patterns of genomic evolution in advanced melanoma

Genomic alterations occurring during melanoma progression and the resulting genomic heterogeneity between metastatic deposits remain incompletely understood. Analyzing 86 metastatic melanoma deposits from 53 patients with whole-exome sequencing (WES), we show a low branch to trunk mutation ratio and little intermetastatic heterogeneity, with driver mutations almost completely shared between lesions. Branch mutations consistent with UV damage indicate that metastases may arise from different subclones in the primary tumor. Selective gain of mutated BRAF alleles occurs as an early event, contrasting whole-genome duplication (WGD) occurring as a late truncal event in about 40% of cases. One patient revealed elevated mutational diversity, probably related to previous chemotherapy and DNA repair defects. In another patient having received radiotherapy toward a lymph node metastasis, we detected a radiotherapy-related mutational signature in two subsequent distant relapses, consistent with secondary metastatic seeding. Our findings add to the understanding of genomic evolution in metastatic melanomas.

As melanoma progresses, it evolves. Here, in advanced melanoma the authors study genomic evolution, highlighting trunk mutations dominated by the ultraviolet damage signature, common late truncal whole-genome duplication events, as well as selective copy number gain of mutant BRAF.

The MDM4 SNP34091 (rs4245739) C-allele is associated with increased risk of ovarian-but not endometrial cancer

The MDM4 protein (also known as MDMX or HDMX) is a negative regulator of p53, not only by direct interaction but also through its interaction with MDM2. Further, MDM4 overexpression and amplification have been observed in several cancer forms. Recently, a single nucleotide polymorphism (SNP) in the 3’ untranslated region of the MDM4 gene, SNP34091A > C (rs4245739) was reported to alter MDM4 messenger RNA (mRNA) stability by modulating a microRNA binding site, thereby leading to decreased MDM4 levels. In this case-control study, we aimed to evaluate the possible association between MDM4 SNP34091 status and cancer risk by comparing the genotype frequencies in large hospital-based cohorts of endometrial- (n = 1404) and ovarian (n = 1385) cancer patients with healthy female controls (n = 1870). Genotype frequencies were compared by odds ratio (OR) estimates and Fisher exact tests. We found that individuals harboring the MDM4 SNP34091AC/CC genotypes had a significantly elevated risk for serous ovarian cancer (SOC) in general and high-grade serous ovarian cancer (HGSOC) in particular (SOC: OR = 1.18., 95 % CI = 1.01–1.39; HGSOC: OR = 1.25, CI = 1.02–1.53). No association between SNP34091 genotypes and endometrial cancer risk was observed. Our data indicate the MDM4 SNP34091AC/CC genotypes to be associated with an elevated risk for SOC and in particular the HGSOC type.

MDM4 SNP34091 (rs4245739) and its effect on breast-, colon-, lung-, and prostate cancer risk

The MDM4 protein plays an important part in the negative regulation of the tumor suppressor p53 through its interaction with MDM2. In line with this, MDM4 amplification has been observed in several tumor forms. A polymorphism (rs4245739 A>C; SNP34091) in the MDM4 3′ untranslated region has been reported to create a target site for hsa‐miR‐191, resulting in decreased MDM4 mRNA levels. In this population‐based case–control study, we examined the potential association between MDM4 SNP34091, alone and in combination with the MDM2 SNP309T>G (rs2279744), and the risk of breast‐, colon‐, lung‐, and prostate cancer in Norway. SNP34091 was genotyped in 7,079 cancer patients as well as in 3,747 gender‐ and age‐matched healthy controls. MDM4 SNP34091C was not associated with risk for any of the tumor forms examined, except for a marginally significant association with reduced risk for breast cancer in a recessive model (OR = 0.77: 95% CI = 0.59–0.99). Stratifying according to MDM2 SNP309 status, we observed a reduced risk for breast cancer related to MDM4 SNP34091CC among individuals harboring the MDM2 SNP309GG genotype (OR = 0.41; 95% CI = 0.21–0.82). We conclude, MDM4 SNP34091 status to be associated with reduced risk of breast cancer, in particular in individuals carrying the MDM2 SNP309GG genotype, but not to be associated with either lung‐, colon‐ or prostate cancer.

Abstract 3082: miR-155-5p expression is associated with response to neoadjuvant paclitaxel treatment of locally advanced breast cancer

Introduction

Resistance to chemotherapy remains a prominent obstacle to cure in cancer. We have previously shown that genetic defects in key pathways involved in processes like apoptosis, cell-cycle arrest and senescence may be involved in resistance to chemotherapy in vivo.

MicroRNAs (miRNAs) are known to be crucial regulators of gene expression, and to potentially affect the response to drugs, including chemotherapy. Here, we aimed to explore the potential role of miRNAs in mechanisms of breast cancer chemoresistance in-vivo, by analyzing miRNA expression in breast cancer patients with recorded responses to chemotherapy.

Patient materials and Methods

Breast tumour biopsies were collected before chemotherapy in a prospective clinical study, where patients were randomized to receive epirubicin or paclitaxel monotherapy (n = 223). Information about the response to therapy, as defined by the UICC criteria, in addition to long-term outcome, was available for all patients. The global miRNA expression levels were analyzed in 50 breast tumour biopsies by deep sequencing. The patients were selected based on their response to first-line chemotherapy. To verify the results obtained by miRNA sequencing, selected miRNAs were quantified in all patients by qPCR.

Results

No global miRNA patterns were observed when performing unsupervised clustering of miRNA sequencing data. Among samples from paclitaxel-treated patients, we identified three differentially expressed miRNAs between response groups; the expression of miR-155 was significantly lower (p = 0.005) in patients with poor response to paclitaxel, and the expression of both miR-181a and miR-181b were significantly higher (p = 0.003 and 0.012, respectively) in patients with poor response to paclitaxel.

qPCR data validated the initial association identified between high miR-155-5p expression levels and response to paclitaxel treatment. Intriguingly though, high levels of miR-155 were also associated with a poor prognosis in the study population as a whole (p = 0.012), as well as with the presence of mutations in TP53 (p &lt; 0.001).

Conclusion

Here, we show that breast cancer patients responding to paclitaxel treatment have higher expression levels of miR-155-5p. Also, patients with elevated levels of miR-155-5p expression have a poor prognosis and a higher frequency of TP53 mutations.

Citation Format: Anne Hege Straume, Einar Birkeland, Stian Knappskog, Per Eystein Lønning. miR-155-5p expression is associated with response to neoadjuvant paclitaxel treatment of locally advanced breast cancer. [abstract]. In: Proceedings of the 106th Annual Meeting of the American Association for Cancer Research; 2015 Apr 18-22; Philadelphia, PA. Philadelphia (PA): AACR; Cancer Res 2015;75(15 Suppl):Abstract nr 3082. doi:10.1158/1538-7445.AM2015-3082

Abstract 2833: Mutation analysis of cancer drivers and DNA repair genes in chemosensitive versus resistant ovarian cancers

Background: The mechanisms underlying resistance to chemotherapy in ovarian cancer are incompletely understood. Identifying genetic alterations associated with treatment response is decisive in the determination of which patients may benefit from adjuvant chemotherapy.

Methods: Biopsies were collected from twenty patients diagnosed with ovarian cancer who were subjected to post-operative taxane- and platinum-containing chemotherapy. Patients were selected for genetic analyses based on response to chemotherapy, determined as time to relapse (10 sensitive and 10 resistant patients). A panel of 620 genes, including known cancer driver genes, as well as genes involved in DNA repair were analysed by massively parallel sequencing. Alignment and mutation calling was performed using MiSeq Reporter, with further manual filtering of variants to exclude common SNPs. Validation of low quality mutation calls was done by Sanger sequencing.

Results: A median of 6 genes (range: 3 - 45) per patient was found to harbour non-synonymous mutations. Among previously identified driver genes in ovarian cancer, we found mutations in TP53, BRCA1, CDK12, NF1 and CSMD3. These mutations were more common among patients with more advanced disease and higher grade. For example, TP53 mutations were found in 10 out of 12 patients with high grade, stage 3c or 4 disease, and in 2 out of 5 with lower stage and/or grade. One patient was found to have a tumor potentially of a hyper-mutator phenotype with 49 mutations in 45 genes identified within our gene panel. With respect to treatment efficacy, 73 and 40 genes were found to be mutated exclusively in patients with a good and poor response to treatment, respectively.

Conclusion: We describe the profile of mutations in cancer driver genes and DNA repair genes among patients suffering from ovarian cancer according to treatment response.

Citation Format: Einar Birkeland, Rakel Blaalid, Merete Bjørnslett, Anne Dørum, Per Eystein Lønning, Stian Knappskog. Mutation analysis of cancer drivers and DNA repair genes in chemosensitive versus resistant ovarian cancers. [abstract]. In: Proceedings of the 105th Annual Meeting of the American Association for Cancer Research; 2014 Apr 5-9; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2014;74(19 Suppl):Abstract nr 2833. doi:10.1158/1538-7445.AM2014-2833

Endometrial Carcinoma Recurrence Score (ECARS) validates to identify aggressive disease and associates with markers of epithelial-mesenchymal transition and PI3K alterations

PURPOSE

Our previously reported 29-gene expression signature identified an aggressive subgroup of endometrial cancer patients with PI3K activation. We here wanted to validate these findings by independent patient series.

PATIENTS AND METHODS

The 29-gene expression signature was assessed in fresh frozen tumor tissue from 280 primary endometrial carcinomas (three independent cohorts), 19 metastatic lesions and in 333 primary endometrial carcinomas using TCGA data, and expression was related to clinico-pathologic features and survival. The 29-gene signature was assessed by real-time quantitative PCR, DNA oligonucleotide microarrays, or RNA sequencing. PI3K alterations were assessed by immunohistochemistry, DNA microarrays, DNA sequencing, SNP arrays or fluorescence in situ hybridization. A panel of markers of epithelial-mesenchymal transition (EMT) was also correlated to the 29-gene signature score.

RESULTS

High 29-gene Endometrial Carcinoma Recurrence Score (ECARS) values consistently validated to identify patients with aggressive clinico-pathologic phenotype and reduced survival. Within the presumed favorable subgroups of low grade, endometrioid tumors confined to the uterus, high ECARS still predicted a poor prognosis. The score was higher in metastatic compared to primary lesions (P<0.001) and was significantly associated with potential measures of PI3K activation, markers of EMT and vascular invasion as an indicator of metastatic spread (all P<0.001).

CONCLUSIONS

ECARS validates to identify aggressive endometrial carcinomas in multiple, independent patients cohorts. The higher signature score in metastatic compared to primary lesions, and the potential link to PI3K activation and EMT, support further studies of ECARS in relation to response to PI3K and EMT inhibitors in clinical trials of metastatic endometrial carcinoma.

Low BRAF and NRAS expression levels are associated with clinical benefit from DTIC therapy and prognosis in metastatic melanoma

Metastatic melanoma is characterized by a poor response to chemotherapy. Furthermore, there is a lack of established predictive and prognostic markers. In this single institution study, we correlated mutation status and expression levels of BRAF and NRAS to dacarbazine (DTIC) treatment response as well as progression-free and overall survival in a cohort of 85 patients diagnosed with advanced melanoma. Neither BRAF nor NRAS mutation status correlated to treatment response. However, patients with tumors harboring NRAS mutations had a shorter overall survival (p < 0.001) compared to patients with tumors wild-type for NRAS. Patients having a clinical benefit (objective response or stable disease at 3 months) on DTIC therapy had lower BRAF and NRAS expression levels compared to patients progressing on therapy (p = 0.037 and 0.003, respectively). For BRAF expression, this association was stronger among patients with tumors wild-type for BRAF (p = 0.005). Further, low BRAF as well as NRAS expression levels were associated with a longer progression-free survival in the total population (p = 0.004 and <0.001, respectively). Contrasting low NRAS expression levels, which were associated with improved overall survival in the total population (p = 0.01), low BRAF levels were associated with improved overall survival only among patients with tumors wild-type for BRAF (p = 0.013). These findings indicate that BRAF and NRAS expression levels may influence responses to DTIC as well as prognosis in patients with advanced melanoma.

KRAS gene amplification and overexpression but not mutation associates with aggressive and metastatic endometrial cancer

BACKGROUND

Three quarter of endometrial carcinomas are treated at early stage. Still, 15 to 20% of these patients experience recurrence, with little effect from systemic therapies. Homo sapiens v-Ki-ras2 Kirsten rat sarcoma viral oncogenes homologue (KRAS) mutations have been reported to have an important role in tumorigenesis for human cancers, but there is limited knowledge regarding clinical relevance of KRAS status in endometrial carcinomas.

METHODS

We have performed a comprehensive and integrated characterisation of genome-wide expression related to KRAS mutations and copy-number alterations in primary- and metastatic endometrial carcinoma lesions in relation to clinical and histopathological data. A primary investigation set and clinical validation set was applied, consisting of 414 primary tumours and 61 metastatic lesions totally.

RESULTS

Amplification and gain of KRAS present in 3% of the primary lesions and 18% of metastatic lesions correlated significantly with poor outcome, high International Federation of Gynaecology and Obstetrics stage, non-endometrioid subtype, high grade, aneuploidy, receptor loss and high KRAS mRNA levels, also found to be associated with aggressive phenotype. In contrast, KRAS mutations were present in 14.7% of primary lesions with no increase in metastatic lesions, and did not influence outcome, but was significantly associated with endometrioid subtype, low grade and obesity.

CONCLUSION

These results support that KRAS amplification and KRAS mRNA expression, both increasing from primary to metastatic lesions, are relevant for endometrial carcinoma disease progression.

Loss of GPER identifies new targets for therapy among a subgroup of ERα-positive endometrial cancer patients with poor outcome

BACKGROUND

The G protein-coupled oestrogen receptor, GPER, has been suggested as an alternative oestrogen receptor. Our purpose was to investigate the potential of GPER as a prognostic and predictive marker in endometrial carcinoma and to search for new drug candidates to improve treatment of aggressive disease.

MATERIALS AND METHOD

A total of 767 primary endometrial carcinomas derived from three patient series, including an external dataset, were studied for protein and mRNA expression levels to investigate and validate if GPER loss identifies poor prognosis and new targets for therapy in endometrial carcinoma. Gene expression levels, according to ERα/GPER status, were used to search the connectivity map database for small molecular inhibitors with potential for treatment of metastatic disease for receptor status subgroups.

RESULTS

Loss of GPER protein is significantly correlated with low GPER mRNA, high FIGO stage, non-endometrioid histology, high grade, aneuploidy and ERα loss (all P-values ≤0.05). Loss of GPER among ERα-positive patients identifies a subgroup with poor prognosis that until now has been unrecognised, with reduced 5-year survival from 93% to 76% (P=0.003). Additional loss of GPER from primary to metastatic lesion counterparts further supports that loss of GPER is associated with disease progression.

CONCLUSION

These results support that GPER status adds clinically relevant information to ERα status in endometrial carcinoma and suggest a potential for new inhibitors in the treatment of metastatic endometrial cancers with ERα expression and GPER loss.