Spatiotemporal genomic profiling of intestinal metaplasia reveals clonal dynamics of gastric cancer progression

Intestinal metaplasia (IM) is a pre-malignant condition of the gastric mucosa associated with increased gastric cancer (GC) risk. Analyzing 1,256 gastric samples (1,152 IMs) across 692 subjects from a prospective 10-year study, we identify 26 IM driver genes in diverse pathways including chromatin regulation (ARID1A) and intestinal homeostasis (SOX9). Single-cell and spatial profiles highlight changes in tissue ecology and IM lineage heterogeneity, including an intestinal stem-cell dominant cellular compartment linked to early malignancy. Expanded transcriptome profiling reveals expression-based molecular subtypes of IM associated with incomplete histology, antral/intestinal cell types, ARID1A mutations, inflammation, and microbial communities normally associated with the healthy oral tract. We demonstrate that combined clinical-genomic models outperform clinical-only models in predicting IMs likely to transform to GC. By highlighting strategies for accurately identifying IM patients at high GC risk and a role for microbial dysbiosis in IM progression, our results raise opportunities for GC precision prevention and interception.

Comprehensive molecular phenotyping of ARID1A-deficient gastric cancer reveals pervasive epigenomic reprogramming and therapeutic opportunities

OBJECTIVE

Gastric cancer (GC) is a leading cause of cancer mortality, with ARID1A being the second most frequently mutated driver gene in GC. We sought to decipher ARID1A-specific GC regulatory networks and examine therapeutic vulnerabilities arising from ARID1A loss.

DESIGN

Genomic profiling of GC patients including a Singapore cohort (>200 patients) was performed to derive mutational signatures of ARID1A inactivation across molecular subtypes. Single-cell transcriptomic profiles of ARID1A-mutated GCs were analysed to examine tumour microenvironmental changes arising from ARID1A loss. Genome-wide ARID1A binding and chromatin profiles (H3K27ac, H3K4me3, H3K4me1, ATAC-seq) were generated to identify gastric-specific epigenetic landscapes regulated by ARID1A. Distinct cancer hallmarks of ARID1A-mutated GCs were converged at the genomic, single-cell and epigenomic level, and targeted by pharmacological inhibition.

RESULTS

We observed prevalent ARID1A inactivation across GC molecular subtypes, with distinct mutational signatures and linked to a NFKB-driven proinflammatory tumour microenvironment. ARID1A-depletion caused loss of H3K27ac activation signals at ARID1A-occupied distal enhancers, but unexpectedly gain of H3K27ac at ARID1A-occupied promoters in genes such as NFKB1 and NFKB2. Promoter activation in ARID1A-mutated GCs was associated with enhanced gene expression, increased BRD4 binding, and reduced HDAC1 and CTCF occupancy. Combined targeting of promoter activation and tumour inflammation via bromodomain and NFKB inhibitors confirmed therapeutic synergy specific to ARID1A-genomic status.

CONCLUSION

Our results suggest a therapeutic strategy for ARID1A-mutated GCs targeting both tumour-intrinsic (BRD4-assocatiated promoter activation) and extrinsic (NFKB immunomodulation) cancer phenotypes.

Regulatory enhancer profiling of mesenchymal-type gastric cancer reveals subtype-specific epigenomic landscapes and targetable vulnerabilities

OBJECTIVE

Gastric cancer (GC) comprises multiple molecular subtypes. Recent studies have highlighted mesenchymal-subtype GC (Mes-GC) as a clinically aggressive subtype with few treatment options. Combining multiple studies, we derived and applied a consensus Mes-GC classifier to define the Mes-GC enhancer landscape revealing disease vulnerabilities.

DESIGN

Transcriptomic profiles of ~1000 primary GCs and cell lines were analysed to derive a consensus Mes-GC classifier. Clinical and genomic associations were performed across >1200 patients with GC. Genome-wide epigenomic profiles (H3K27ac, H3K4me1 and assay for transposase-accessible chromatin with sequencing (ATAC-seq)) of 49 primary GCs and GC cell lines were generated to identify Mes-GC-specific enhancer landscapes. Upstream regulators and downstream targets of Mes-GC enhancers were interrogated using chromatin immunoprecipitation followed by sequencing (ChIP-seq), RNA sequencing, CRISPR/Cas9 editing, functional assays and pharmacological inhibition.

RESULTS

We identified and validated a 993-gene cancer-cell intrinsic Mes-GC classifier applicable to retrospective cohorts or prospective single samples. Multicohort analysis of Mes-GCs confirmed associations with poor patient survival, therapy resistance and few targetable genomic alterations. Analysis of enhancer profiles revealed a distinctive Mes-GC epigenomic landscape, with TEAD1 as a master regulator of Mes-GC enhancers and Mes-GCs exhibiting preferential sensitivity to TEAD1 pharmacological inhibition. Analysis of Mes-GC super-enhancers also highlighted NUAK1 kinase as a downstream target, with synergistic effects observed between NUAK1 inhibition and cisplatin treatment.

CONCLUSION

Our results establish a consensus Mes-GC classifier applicable to multiple transcriptomic scenarios. Mes-GCs exhibit a distinct epigenomic landscape, and TEAD1 inhibition and combinatorial NUAK1 inhibition/cisplatin may represent potential targetable options.

Chromatin Rewiring by Mismatch Repair Protein MSH2 Alters Cell Adhesion Pathways and Sensitivity to BET Inhibition in Gastric Cancer

Mutations in the DNA mismatch repair gene MSH2 are causative of microsatellite instability (MSI) in multiple cancers. Here, we discovered that besides its well-established role in DNA repair, MSH2 exerts a novel epigenomic function in gastric cancer. Unbiased CRISPR-based mass spectrometry combined with genome-wide CRISPR functional screening revealed that in early-stage gastric cancer MSH2 genomic binding is not randomly distributed but rather is associated specifically with tumor-associated super-enhancers controlling the expression of cell adhesion genes. At these loci, MSH2 genomic binding was required for chromatin rewiring, de novo enhancer-promoter interactions, maintenance of histone acetylation levels, and regulation of cell adhesion pathway expression. The chromatin function of MSH2 was independent of its DNA repair catalytic activity but required MSH6, another DNA repair gene, and recruitment to gene loci by the SWI/SNF chromatin remodeler SMARCA4/BRG1. Loss of MSH2 in advanced gastric cancers was accompanied by deficient cell adhesion pathway expression, epithelial-mesenchymal transition, and enhanced tumorigenesis in vitro and in vivo. However, MSH2-deficient gastric cancers also displayed addiction to BAZ1B, a bromodomain-containing family member, and consequent synthetic lethality to bromodomain and extraterminal motif (BET) inhibition. Our results reveal a role for MSH2 in gastric cancer epigenomic regulation and identify BET inhibition as a potential therapy in MSH2-deficient gastric malignancies.

SIGNIFICANCE

DNA repair protein MSH2 binds and regulates cell adhesion genes by enabling enhancer-promoter interactions, and loss of MSH2 causes deficient cell adhesion and bromodomain and extraterminal motif inhibitor synthetic lethality in gastric cancer.

Direct and indirect costs for hospitalized patients with dengue in Southern Sri Lanka

Background

The Southern Province of Sri Lanka is endemic with dengue, with frequent outbreaks and occurrence of severe disease. However, the economic burden of dengue is poorly quantified. Therefore, we conducted a cost analysis to assess the direct and indirect costs associated with hospitalized patients with dengue to households and to the public healthcare system.

Methods

From June 2017–December 2018, we prospectively enrolled children and adults with acute dengue hospitalized at the largest, public tertiary-care (1800 bed) hospital in the Southern Province, Sri Lanka. We administered a structured questionnaire to obtain information regarding direct costs spent by households on medical visits, medications, laboratory testing, and travel for seeking care for the illness. Indirect costs lost by households were estimated by identifying the days of work lost by patients and caregivers and school days lost by children. Direct hospital costs were estimated using gross costing approach and adjusted by multiplying by annual inflation rates in Sri Lankan rupees and converted to US dollars.

Results

A total of 1064 patients with laboratory-confirmed dengue were enrolled. The mean age (SD) was 35.9 years (15.6) with male predominance (66.2%). The mean durations of hospitalization for adults and paediatric patients were 3.86 (SD = 1.51) and 4 (SD = 1.32) days, respectively. The per-capita direct cost borne by the healthcare system was 233.76 USD, and was approximately 14 times greater than the per-capita direct cost borne by households (16.29 USD, SD = 14.02). The per-capita average number of loss of working days was 21.51 (SD = 41.71), with mean per-capita loss of income due to loss of work being 303.99 USD (SD = 569.77), accounting for over 70% of average monthly income. On average, 10.88 days (SD = 10.97) of school days were missed due to the dengue episode. School misses were expected to reduce future annual income of affected children by 0.44%.

Conclusions

Dengue requiring hospitalization had a substantial economic burden on the public healthcare system in Sri Lanka and the affected households. These findings emphasize the importance of strengthening dengue control activities and improved use of hospital-based resources for care to reduce the economic impact of dengue in Sri Lanka.

Integrative epigenomic and high-throughput functional enhancer profiling reveals determinants of enhancer heterogeneity in gastric cancer

BACKGROUND

Enhancers are distal cis-regulatory elements required for cell-specific gene expression and cell fate determination. In cancer, enhancer variation has been proposed as a major cause of inter-patient heterogeneity-however, most predicted enhancer regions remain to be functionally tested.

METHODS

We analyzed 132 epigenomic histone modification profiles of 18 primary gastric cancer (GC) samples, 18 normal gastric tissues, and 28 GC cell lines using Nano-ChIP-seq technology. We applied Capture-based Self-Transcribing Active Regulatory Region sequencing (CapSTARR-seq) to assess functional enhancer activity. An Activity-by-contact (ABC) model was employed to explore the effects of histone acetylation and CapSTARR-seq levels on enhancer-promoter interactions.

RESULTS

We report a comprehensive catalog of 75,730 recurrent predicted enhancers, the majority of which are GC-associated in vivo (> 50,000) and associated with lower somatic mutation rates inferred by whole-genome sequencing. Applying CapSTARR-seq to the enhancer catalog, we observed significant correlations between CapSTARR-seq functional activity and H3K27ac/H3K4me1 levels. Super-enhancer regions exhibited increased CapSTARR-seq signals compared to regular enhancers, even when decoupled from native chromatin contexture. We show that combining histone modification and CapSTARR-seq functional enhancer data improves the prediction of enhancer-promoter interactions and pinpointing of germline single nucleotide polymorphisms (SNPs), somatic copy number alterations (SCNAs), and trans-acting TFs involved in GC expression. We identified cancer-relevant genes (ING1, ARL4C) whose expression between patients is influenced by enhancer differences in genomic copy number and germline SNPs, and HNF4α as a master trans-acting factor associated with GC enhancer heterogeneity.

CONCLUSIONS

Our results indicate that combining histone modification and functional assay data may provide a more accurate metric to assess enhancer activity than either platform individually, providing insights into the relative contribution of genetic (cis) and regulatory (trans) mechanisms to GC enhancer functional heterogeneity.

Evaluation of the probiotic and functional potential of Lactobacillus agilis 32 isolated from pig manure

Escherichia coli is a symbiotic bacterium in humans and animals and an important pathogen of humans and animals. Prevention and suppression of E. coli infection is of great concern. In this study, we isolated a strain of Lactobacillus agilis 32 from pig manure and evaluated its biological characteristics, and found that its bacterial survival rate was 25% after 4 h of treatment at pH 2, and under the condition of 0·5% bile concentration, its survival rate exceeds 30%. In addition, L. agilis 32 has a cell surface hydrophobicity of 77·8%, and exhibits 67·1% auto-aggregation and 63·2% aggregation with Enterotoxigenic E. coli 10 (ETEC 10). FITC fluorescence labelling showed that the fluorescence intensity of cecum was significantly higher than that of duodenum, jejunum or colon (P < 0·05), but no significant difference from ileum. Lactobacillus agilis 32 bacterial culture and CFS showed average inhibition zone diameters of 14·2 and 15·4 mm respectively. Lactobacillus agilis 32 CFS treatment can significantly reduce the pathogenicity of ETEC 10. These results show that L. agilis 32 is an active and potential probiotic, and it has a good antibacterial effect on ETEC10, which provides basic research for probiotics to prevent and treat intestinal diarrhoea pathogen infection.

Genomic and epigenomic EBF1 alterations modulate TERT expression in gastric cancer

Transcriptional reactivation of telomerase catalytic subunit (TERT) is a frequent hallmark of cancer, occurring in 90% of human malignancies. However, specific mechanisms driving TERT reactivation remain obscure for many tumor types and in particular gastric cancer (GC), a leading cause of global cancer mortality. Here, through comprehensive genomic and epigenomic analysis of primary GCs and GC cell lines, we identified the transcription factor early B cell factor 1 (EBF1) as a TERT transcriptional repressor and inactivation of EBF1 function as a major cause of TERT upregulation. Abolishment of EBF1 function occurs through 3 distinct (epi)genomic mechanisms. First, EBF1 is epigenetically silenced via DNA methyltransferase, polycomb-repressive complex 2 (PRC2), and histone deacetylase activity in GCs. Second, recurrent, somatic, and heterozygous EBF1 DNA-binding domain mutations result in the production of dominant-negative EBF1 isoforms. Third, more rarely, genomic deletions and rearrangements proximal to the TERT promoter remobilize or abolish EBF1-binding sites, derepressing TERT and leading to high TERT expression. EBF1 is also functionally required for various malignant phenotypes in vitro and in vivo, highlighting its importance for GC development. These results indicate that multimodal genomic and epigenomic alterations underpin TERT reactivation in GC, converging on transcriptional repressors such as EBF1.

DNA epigenetic signature predictive of benefit from neoadjuvant chemotherapy in oesophageal adenocarcinoma: results from the MRC OE02 trial

BACKGROUND

DNA methylation signatures describing distinct histological subtypes of oesophageal cancer have been reported. We studied DNA methylation in samples from the MRC OE02 phase III trial, which randomised patients with resectable oesophageal cancer to surgery alone (S) or neoadjuvant chemotherapy followed by surgery (CS).

AIM

The aim of the study was to identify epigenetic signatures predictive of chemotherapy benefit in patients with oesophageal adenocarcinoma (OAC) from the OE02 trial and validate the findings in an independent cohort.

METHODS

DNA methylation was analysed using the Illumina GoldenGate platform on surgically resected OAC specimens from patients in the OE02 trial. Cox proportional hazard analysis was performed to select probes predictive of survival in the CS arm. Non-negative matrix factorisation was used to perform clustering and delineate DNA methylation signatures. The findings were validated in an independent cohort of patients with gastroesophageal adenocarcinoma treated with neoadjuvant chemotherapy.

RESULTS

A total of 229 patients with OAC were analysed from the OE02 trial (118 in the CS arm and 111 in the S arm). There was no difference in DNA methylation status between the CS and S arms. A metagene signature was created by dichotomising samples into two clusters. In cluster 1, patients in the CS arm had significant overall survival (OS) benefit (median OS CS: 931 days vs. S: 536 days [HR: 1.54, P = 0.031]). In cluster 2, patients in the CS arm had similar (or worse) OS compared with patients in the S arm (CS: 348 days vs. S: 472 days [HR: 0.70, P = 0.1], and test of interaction was significant (p = 0.005). In the validation cohort (n = 13), there was no difference in DNA methylation status in paired pre- and post-treatment samples. When the epigenetic signature was applied, cluster 1 samples had better OS (median OS, cluster 1: 1174 days vs. cluster 2: 392 days, HR: 3.47, p = 0.059) CONCLUSIONS: This is the first and largest study of DNA methylation in patients with OAC uniformly treated in a randomised phase III trial. We identified an epigenetic signature that may serve as a predictive biomarker for chemotherapy benefit in OAC.

DNA methylation signature predictive of benefit from neoadjuvant chemotherapy in esophageal adenocarcinoma: Results from the MRC OEO2 phase III trial

43

Background: Platinum and 5-Fluorouracil (5FU) neoadjuvant chemotherapy followed by surgery is one of the standard approaches for patients with resectable EAC. To date, there are no predictive biomarkers of chemotherapy benefit. We hypothesize that DNA methylation of genes in key biologic and oncogenic pathways predict for chemotherapy benefit in EAC. Methods: In the OE02 trial, 802 patients with resectable esophageal carcinoma were randomised to surgery alone (S) versus two cycles of cisplatin and 5FU chemotherapy followed by surgery (CS). DNA was extracted from 213 EAC resection specimens (110 from the (CS) arm, 103 from the (S) arm). DNA methylation was analyzed at 1505 CpG sites within 807 genes using the Illumina GoldenGate platform. Cox proportional hazard analysis was performed to identify predictive markers of survival in (CS) arm; non-negative matrix factorization (NMF) was used to delineate methylation signatures. Results: Methylation status of 1505 CpG sites had no statistical difference between the (CS) and (S) arms. In the (CS) arm, 87 (5.7%) CpG sites were initially identified as promising candidates in univariate analysis (p < 0.05 cutoff). NMF generated a 4 CpG site signature which divided patients into poor risk and good risk. Genes involved in the signature include RUNX1T1, CCND2, MST1R and MMP14. Survival was significantly different between poor risk and good risk in (CS) arm (HR 0.32, 95% CI: 0.21 to 0.52, p < 0.0001). No difference in survival was detected in the surgery arm (HR 1.12, 95% CI: 0.76 to 1.80, p = 0.48), suggesting the signature served as a predictive and not prognostic biomarker. Methylation signature remained an independent predictor of survival in multivariate analysis with clinicopathologic factors (along with age and vascular invasion). Conclusions: Chemotherapy does not appear to change methylation status of EAC. Hypermethylation of RUNX1T1, CCND2 and hypomethylation of MST1R and MMP14 leads to significantly decreased benefit from chemotherapy in EA. We describe an epigenetic signature which may serve as a predictive biomarker for chemotherapy benefit using data form the largest bank of DNA methylation in EA reported to date.

DEF: an automated dead-end filling approach based on quasi-endosymbiosis

Motivation

Gap filling for the reconstruction of metabolic networks is to restore the connectivity of metabolites via finding high-confidence reactions that could be missed in target organism. Current methods for gap filling either fall into the network topology or have limited capability in finding missing reactions that are indirectly related to dead-end metabolites but of biological importance to the target model.

Results

We present an automated dead-end filling (DEF) approach, which is derived from the wisdom of endosymbiosis theory, to fill gaps by finding the most efficient dead-end utilization paths in a constructed quasi-endosymbiosis model. The recalls of reactions and dead ends of DEF reach around 73% and 86%, respectively. This method is capable of finding indirectly dead-end-related reactions with biological importance for the target organism and is applicable to any given metabolic model. In the E. coli iJR904 model, for instance, about 42% of the dead-end metabolites were fixed by our proposed method.

Availabilty and Implementaion

DEF is publicly available at http://bis.zju.edu.cn/DEF/.

Contact

mchen@zju.edu.cn

Supplimentary Information

Supplementary data are available at Bioinformatics online.

Determination of an arylether antiarrhythmic and its N-dealkyl metabolite in rat plasma and hepatic microsomal incubates using liquid chromatography-tandem mass spectrometry

A method was developed and validated for the quantification of (+/-)-trans-[2-morpholino-1-(1-naphthalene-ethyloxy]cyclohexane monohydrochloride (RSD1070) and its N-dealkyl metabolite in rat plasma and hepatic microsomal incubates. Chromatographic separations were achieved using reversed-phase high-performance liquid chromatography coupled with positive ion electrospray ionization and detection by tandem mass spectrometry. The assay was linear from 2.5 to 100 ng/ml and this range was used for validation. Inter- and intra-assay variability (n=6), extraction recovery, and stability in plasma were assessed. The estimated limit of quantitation was in the range 2.5-3 ng/ml for both analytes in rat plasma. The analytical method was used in a pharmacokinetic study of RSD1070 in rats after a single i.v. bolus of 12 mg/kg.

Does pheromone biology of Lambdina athasaria and L. pellucidaria contribute to their reproductive isolation?

Recently, 7-methylheptadecane and 7,11-dimethylheptadecane have been reported as sex pheromone components of both spring hemlock looper (SHL), Lambdina athasaria, and pitch pine looper (PPL), Lambdina pellucidaria. Our objective was to test the hypothesis that SHL and PPL are reproductively isolated, in part, through species specificity in: (1) absolute configuration of pheromone components, (2) diel periodicity of pheromonal communication, and/or (3) seasonal flight period. In coupled gas chromatographic-electroantennographic detection (GC-EAD) analyses of stereoselectively synthesized (7S)- and (7R)-7-methylheptadecane [7S; 7R] as well as (7S,11S)-, (7R,11R)-, and (meso-7,11)-7,11-dimethylheptadecane [7S,11S; 7R,11R; meso-7,11], only 7S and meso-7,11 elicited responses by male SHL and PPL antennae. In field experiments, male SHL and PPL were attracted only to lures containing 7S plus meso-7,11. In hourly recordings of trap-captured males, SHL and PPL in their respective habitats were trapped between 24:00 and 03:00 hr. Capture of both SHL and PPL in pheromone-baited traps throughout June indicated overlapping seasonal flight periods. These findings of identical absolute configuration of pheromoal components, diel periodicity of pheromonal communication, and overlap of seasonal flight periods support synonymy of SHL and PPL. Finite taxonomic classification of PPL and SHL must await careful assessment of further criteria, such as morphometrics, molecular comparisons and ecological analyses.