Coding and non-coding co-expression network analysis identifies key modules and driver genes associated with precursor lesions of gastric cancer

Background: Helicobacter pylori infection is the most important risk factor for gastric cancer (GC). Human gastric adenocarcinoma develops after long-term H. pylori infection via the Correa cascade. This carcinogenic pathway describes the progression from gastritis to atrophy, intestinal metaplasia (IM), dysplasia and GC. Patients with atrophy and intestinal metaplasia are considered to have precancerous lesions of GC (PLGC). H. pylori eradication and endoscopy surveillance are currently the only interventions for preventing GC. Better knowledge of the biology of human PLGC may help find stratification markers and contribute to better understanding of biological mechanisms. One way to achieve this is by using co-expression network analysis. Weighted gene co-expression network analysis (WGCNA) is often used to identify modules from co-expression networks and relate them to clinical traits. It also allows identification of driver genes that may be critical for PLGC. Aim: The purpose of this study was to identify co-expression modules and differential gene expression in dyspeptic patients at different stages of the Correa pathway. Methods: We studied 96 gastric biopsies from 78 patients that were clinically classified as: non-active (n = 10) and chronic-active gastritis (n = 20), atrophy (n = 12), and IM (n = 36). Gene expression of coding RNAs was determined by microarrays and non-coding RNAs by RNA-seq. The WGCNA package was used for network construction, module detection, module preservation and hub and driver gene selection. Results: WGCNA identified 20 modules for coding RNAs and 4 for each miRNA and small RNA class. Modules were associated with antrum and corpus gastric locations, chronic gastritis and IM. Notably, coding RNA modules correlated with the Correa cascade. One was associated with the presence of H. pylori. In three modules, the module eigengene (ME) gradually increased in the stages toward IM, while in three others the inverse relationship was found. One miRNA module was negatively correlated to IM and was used for a mRNA-miRNA integration analysis. WGCNA also uncovered driver genes. Driver genes show both high connectivity within a module and are significantly associated with clinical traits. Some of those genes have been previously involved in H. pylori carcinogenesis, but others are new. Lastly, using similar external transcriptomic data, we confirmed that the discovered mRNA modules were highly preserved. Conclusion: Our analysis captured co-expression modules that provide valuable information to understand the pathogenesis of the progression of PLGC.