Comparison of genomic and transcriptional microbiome analysis in gastric cancer patients and healthy individuals

Abstract

BACKGROUND Helicobacter pylori and the stomach microbiome play a crucial role in gastric carcinogenesis, and detailed characterization of the microbiome is necessary for a better understanding of the pathophysiology of the disease. There are two common modalities for microbiome analysis: DNA (16S rRNA gene) and RNA (16S rRNA transcript) sequencing. The implications from the use of one or another sequencing approach on the characterization and comparability of the mucosal microbiome in gastric cancer (GC) are poorly studied. AIM To characterize the microbiota of GC using 16S rRNA gene and its transcript and determine difference in the bacterial composition. METHODS In this study, 316 DNA and RNA samples extracted from 105 individual stomach biopsies were included. The study cohort consisted of 29 healthy control individuals and 76 patients with GC. Gastric tissue biopsy samples were collected from damaged mucosa and healthy mucosa at least 5 cm from the tumor tissue. From the controls, healthy stomach mucosa biopsies were collected. From all biopsies RNA and DNA were extracted. RNA was reverse transcribed into cDNA. V1-V2 region of bacterial 16S rRNA gene from all samples were amplified and sequenced on an Illumina MiSeq platform. Bray-Curtis algorithm was used to construct sample-similarity matrices abundances of taxonomic ranks in each sample type. For significant differences between groups permutational multivariate analysis of variance and Mann-Whitney test followed by false-discovery rate test were used. RESULTS Microbial analysis revealed that only a portion of phylotypes (18%-30%) overlapped between microbial profiles obtained from DNA and RNA samples. Detailed analysis revealed differences between GC and controls depending on the chosen modality, identifying 17 genera at the DNA level and 27 genera at the RNA level. Ten of those bacteria were found to be different from the control group at both levels. The key taxa showed congruent results in various tests used; however, differences in 7 bacteria taxa were found uniquely only at the DNA level, and 17 uniquely only at the RNA level. Furthermore, RNA sequencing was more sensitive for detecting differences in bacterial richness, as well as differences in the relative abundance of Reyranella and Sediminibacterium according to the type of GC. In each study group (control, tumor, and tumor adjacent) were found differences between DNA and RNA bacterial profiles. CONCLUSION Comprehensive microbial study provides evidence for the effect of choice of sequencing modality on the microbiota profile, as well as on the identified differences between case and control.