In the competing endogenous RNA (ceRNA) hypothesis, different transcripts communicate through a competition for their common targeting microRNAs (miRNAs). Individual examples have clearly shown the functional importance of ceRNA in gene regulation and cancer biology. It remains unclear to what extent gene expression levels are regulated by ceRNA in general. One major hurdle to studying this problem is the intertwined connections in miRNA-target networks, which makes it difficult to isolate the effects of individual miRNAs.

Patient-derived tumor xenografts in mice are widely used in cancer research and have become important in developing personalized therapies. When these xenografts are subject to DNA sequencing, the samples could contain various amounts of mouse DNA. It has been unclear how the mouse reads would affect data analyses. We conducted comprehensive simulations to compare three alignment strategies at different mutation rates, read lengths, sequencing error rates, human-mouse mixing ratios and sequenced regions. We also sequenced a nasopharyngeal carcinoma xenograft and a cell line to test how the strategies work on real data.

Statistical models have been used to quantify the relationship between gene expression and transcription factor (TF) binding signals. Here we apply the models to the large-scale data generated by the ENCODE project to study transcriptional regulation by TFs. Our results reveal a notable difference in the prediction accuracy of expression levels of transcription start sites (TSSs) captured by different technologies and RNA extraction protocols. In general, the expression levels of TSSs with high CpG content are more predictable than those with low CpG content. For genes with alternative TSSs, the expression levels of downstream TSSs are more predictable than those of the upstream ones. Different TF categories and specific TFs vary substantially in their contributions to predicting expression. Between two cell lines, the differential expression of TSS can be precisely reflected by the difference of TF-binding signals in a quantitative manner, arguing against the conventional on-and-off model of TF binding. Finally, we explore the relationships between TF-binding signals and other chromatin features such as histone modifications and DNase hypersensitivity for determining expression. The models imply that these features regulate transcription in a highly coordinated manner.

We present a new method, OMSV, for accurately and comprehensively identifying structural variations (SVs) from optical maps. OMSV detects both homozygous and heterozygous SVs, SVs of various types and sizes, and SVs with or without creating or destroying restriction sites. We show that OMSV has high sensitivity and specificity, with clear performance gains over the latest method. Applying OMSV to a human cell line, we identified hundreds of SVs >2 kbp, with 68 % of them missed by sequencing-based callers. Independent experimental validation confirmed the high accuracy of these SVs. The OMSV software is available at http://yiplab.cse.cuhk.edu.hk/omsv/ .

Tumor subtype-specific metabolic reprogrammers could serve as targets of therapeutic intervention. Here we show that triple-negative breast cancer (TNBC) exhibits a hyper-activated cholesterol-biosynthesis program that is strongly linked to nuclear receptor RORγ, compared to estrogen receptor-positive breast cancer. Genetic and pharmacological inhibition of RORγ reduces tumor cholesterol content and synthesis rate while preserving host cholesterol homeostasis. We demonstrate that RORγ functions as an essential activator of the entire cholesterol-biosynthesis program, dominating SREBP2 via its binding to cholesterol-biosynthesis genes and its facilitation of the recruitment of SREBP2. RORγ inhibition disrupts its association with SREBP2 and reduces chromatin acetylation at cholesterol-biosynthesis gene loci. RORγ antagonists cause tumor regression in patient-derived xenografts and immune-intact models. Their combination with cholesterol-lowering statins elicits superior anti-tumor synergy selectively in TNBC. Together, our study uncovers a master regulator of the cholesterol-biosynthesis program and an attractive target for TNBC.

Nasopharyngeal cancer (NPC), endemic in Southeast Asia, lacks effective diagnostic and therapeutic strategies. Even in high-income countries the 5-year survival rate for stage IV NPC is less than 40%. Here we report high somatostatin receptor 2 (SSTR2) expression in multiple clinical cohorts comprising 402 primary, locally recurrent and metastatic NPCs. We show that SSTR2 expression is induced by the Epstein-Barr virus (EBV) latent membrane protein 1 (LMP1) via the NF-κB pathway. Using cell-based and preclinical rodent models, we demonstrate the therapeutic potential of SSTR2 targeting using a cytotoxic drug conjugate, PEN-221, which is found to be superior to FDA-approved SSTR2-binding cytostatic agents. Furthermore, we reveal significant correlation of SSTR expression with increased rates of survival and report in vivo uptake of the SSTR2-binding 68Ga-DOTA-peptide radioconjugate in PET-CT scanning in a clinical trial of NPC patients (NCT03670342). These findings reveal a key role in EBV-associated NPC for SSTR2 in infection, imaging, targeted therapy and survival.

G9a is a lysine methyltransferase that regulates epigenetic modifications, transcription, and genome organization. However, whether these properties are dependent on one another or represent distinct functions of G9a remains unclear. In this study, we observe widespread DNA methylation loss in G9a depleted and catalytic mutant embryonic stem cells. Furthermore, we define how G9a regulates chromatin accessibility, epigenetic modifications, and transcriptional silencing in both catalytic-dependent and -independent manners. Reactivated retrotransposons provide alternative promoters and splice sites leading to the upregulation of neighboring genes and the production of chimeric transcripts. Moreover, while topologically associated domains and compartment A/B definitions are largely unaffected, the loss of G9a leads to altered chromatin states, aberrant CTCF and cohesin binding, and differential chromatin looping, especially at retrotransposons. Taken together, our findings reveal how G9a regulates the epigenome, transcriptome, and higher-order chromatin structures in distinct mechanisms.

Pulmonary lymphoepithelioma-like carcinoma (LELC) is a subtype of non-small cell lung cancer (NSCLC) characterized by marked lymphocytic infiltration and association with Epstein-Barr virus (EBV). The molecular basis underlying the disease remains unclear. We sought to study the molecular landscape by multiple approaches including whole genomic sequencing, capture-based targeted sequencing, fluorescent in situ hybridization and immunohistochemistry. Tumor cells from 57 EBV-positive pulmonary LELCs were isolated by careful microdissection prior to genomic sequencing. Integrated analysis revealed a distinct genomic landscape of low TP53 mutation rate (11%), low incidence of known drivers in the RTK/RAS/RAF (11%) and PI3K/AKT/mTOR pathways (7%), but enriched for loss-of-function mutations in multiple negative regulators of the NF-κB pathway. High level programmed cell death ligand-1 (PD-L1) expression was shown with 47% and 79% of the cases showing positive PD-L1 immunoreactivity at ≥50% and ≥1% tumor proportion score, respectively. Subsets of the patients with actionable fibroblast growth factor receptor 3 (FGFR3) aberrations (4%) and mismatch repair deficiency (4%) were potentially eligible for precision medicine. Pulmonary LELC showed a distinct genomic landscape, different from major NSCLC subtypes but resembled that of EBV-associated nasopharyngeal carcinoma. Our work facilitated the understanding of molecular basis underlying pulmonary LELC to explore potential therapeutic options.

Enhancers are important non-coding elements, but they have traditionally been hard to characterize experimentally. The development of massively parallel assays allows the characterization of large numbers of enhancers for the first time. Here, we developed a framework using Drosophila STARR-seq to create shape-matching filters based on meta-profiles of epigenetic features. We integrated these features with supervised machine-learning algorithms to predict enhancers. We further demonstrated that our model could be transferred to predict enhancers in mammals. We comprehensively validated the predictions using a combination of in vivo and in vitro approaches, involving transgenic assays in mice and transduction-based reporter assays in human cell lines (153 enhancers in total). The results confirmed that our model can accurately predict enhancers in different species without re-parameterization. Finally, we examined the transcription factor binding patterns at predicted enhancers versus promoters. We demonstrated that these patterns enable the construction of a secondary model that effectively distinguishes enhancers and promoters.

Long non-coding RNA HOX Transcript Antisense RNA (HOTAIR) is overexpressed in multiple cancers with diverse genetic profiles. Importantly, since HOTAIR heavily contributes to cancer progression by promoting tumor growth and metastasis, HOTAIR becomes a potential target for cancer therapy. However, the underlying mechanism leading to HOTAIR deregulation is largely unexplored. Here, we performed a pan-cancer analysis using more than 4,200 samples and found that intragenic exon CpG island (Ex-CGI) was hypermethylated and was positively correlated to HOTAIR expression. Also, we revealed that Ex-CGI methylation promotes HOTAIR expression through enhancing the transcription elongation process. Furthermore, we linked up the aberrant intragenic tri-methylation on H3 at lysine 4 (H3K4me3) and Ex-CGI DNA methylation in promoting transcription elongation of HOTAIR. Targeting the oncogenic CDK7-CDK9-H3K4me3 axis downregulated HOTAIR expression and inhibited cell growth in many cancers. To our knowledge, this is the first time that a positive feedback loop that involved CDK9-mediated phosphorylation of RNA Polymerase II Serine 2 (RNA PolII Ser2), H3K4me3, and intragenic DNA methylation, which induced robust transcriptional elongation and heavily contributed to the upregulation of oncogenic lncRNA in cancer has been demonstrated. Targeting the oncogenic CDK7-CDK9-H3K4me3 axis could be a novel therapy in many cancers through inhibiting the HOTAIR expression.

Urinary bladder cancer is a common cancer worldwide. Currently, the modality of treating and monitoring bladder cancer is wide. Nonetheless, the high recurrence rate of non-muscle-invasive bladder cancer after surgical resection is still unsatisfactory. Hereby, our study demonstrated whether the intra-operative and post-operative environments will affect bladder cancer recurrence utilizing in vitro cell line model. Bladder cancer cell lines were submerged in four different irrigating fluids for assessing their tumorigenic properties. Our results showed that sterile water performed the best in terms of the magnitude of cytotoxicity to cell lines. Besides, we also investigated cytotoxic effects of the four irrigating agents as well as mitomycin C (MMC) in normothermic and hyperthermic conditions. We observed that sterile water and MMC had an increased cytotoxic effect to bladder cancer cell lines in hyperthermic conditions. Altogether, our results could be translated into clinical practice in the future by manipulating the intra-operative and post-operative conditions in order to lower the chance of residual cancer cells reimplant onto the bladder, which in turns, reducing the recurrence rate of bladder cancers.

Tumor innervation is a common phenomenon with unknown mechanism. Here, we discovered a direct mechanism of tumor-associated macrophage (TAM) for promoting de novo neurogenesis via a subset showing neuronal phenotypes and pain receptor expression associated with cancer-driven nocifensive behaviors. This subset is rich in lung adenocarcinoma associated with poorer prognosis. By elucidating the transcriptome dynamics of TAM with single-cell resolution, we discovered a phenomenon "macrophage to neuron-like cell transition" (MNT) for directly promoting tumoral neurogenesis, evidenced by macrophage depletion and fate-mapping study in lung carcinoma models. Encouragingly, we detected neuronal phenotypes and activities of the bone marrow-derived MNT cells (MNTs) in vitro. Adoptive transfer of MNTs into NOD/SCID mice markedly enhanced their cancer-associated nocifensive behaviors. We identified macrophage-specific Smad3 as a pivotal regulator for promoting MNT at the genomic level; its disruption effectively blocked the tumor innervation and cancer-dependent nocifensive behaviors in vivo. Thus, MNT may represent a precision therapeutic target for cancer pain.

Therapy-induced tumour microenvironment (TME) remodelling poses a major hurdle for cancer cure. As the majority of patients with hepatocellular carcinoma (HCC) exhibits primary or acquired resistance to antiprogrammed cell death (ligand)-1 (anti-PD-[L]1) therapies, we aimed to investigate the mechanisms underlying tumour adaptation to immune-checkpoint targeting.

Immunodynamics in the tumor microenvironment can be precisely examined by using multiple antigen identification approaches. Here, we present a protocol for capturing expression levels of multiple target proteins in the same specimen at single-cell resolution using a tyramide signal amplification-based immunofluorescent multiplexing system. We describe steps for tumor tissue microarray preparation, multiplex immunohistochemistry staining, image acquisition, and quantification. This protocol can quantify immune cells in tissues from patients or experimental disease models at a protein level. For complete details on the use and execution of this protocol, please refer to Chung et al. (2023),1 Tang et al. (2022),2 and Tang et al. (2022).3.

The precision-recall curve (PRC) and the area under it (AUPRC) are useful for quantifying classification performance. They are commonly used in situations with imbalanced classes, such as cancer diagnosis and cell type annotation. We evaluated 10 popular tools for plotting PRC and computing AUPRC, which were collectively used in >3,000 published studies. We found the AUPRC values computed by the tools rank classifiers differently and some tools produce overly-optimistic results.

Natural small compounds comprise most cellular molecules and bind proteins as substrates, products, cofactors, and ligands. However, a large-scale investigation of in vivo protein-small metabolite interactions has not been performed. We developed a mass spectrometry assay for the large-scale identification of in vivo protein-hydrophobic small metabolite interactions in yeast and analyzed compounds that bind ergosterol biosynthetic proteins and protein kinases. Many of these proteins bind small metabolites; a few interactions were previously known, but the vast majority are new. Importantly, many key regulatory proteins such as protein kinases bind metabolites. Ergosterol was found to bind many proteins and may function as a general regulator. It is required for the activity of Ypk1, a mammalian AKT/SGK kinase homolog. Our study defines potential key regulatory steps in lipid biosynthetic pathways and suggests that small metabolites may play a more general role as regulators of protein activity and function than previously appreciated.

Many DNA methylome profiling methods cannot distinguish between 5-methylcytosine (5mC) and 5-hydroxymethylcytosine (5hmC). Because 5mC typically acts as a repressive mark whereas 5hmC is an intermediate form during active demethylation, the inability to separate their signals could lead to incorrect interpretation of the data. Is the extra information contained in 5hmC signals worth the additional experimental and computational costs? Here we combine whole-genome bisulfite sequencing (WGBS) and oxidative WGBS (oxWGBS) data in various human tissues to investigate the quantitative relationships between gene expression and the two forms of DNA methylation at promoters, transcript bodies, and immediate downstream regions. We find that 5mC and 5hmC signals correlate with gene expression in the same direction in most samples. Considering both types of signals increases the accuracy of expression levels inferred from methylation data by a median of 18.2% as compared to having only WGBS data, showing that the two forms of methylation provide complementary information about gene expression. Differential analysis between matched tumor and normal pairs is particularly affected by the superposition of 5mC and 5hmC signals in WGBS data, with at least 25%-40% of the differentially methylated regions (DMRs) identified from 5mC signals not detected from WGBS data. Our results also confirm a previous finding that methylation signals at transcript bodies are more indicative of gene expression levels than promoter methylation signals. Overall, our study provides data for evaluating the cost-effectiveness of some experimental and analysis options in the study of DNA methylation in normal and cancer samples.

The lack of representative nasopharyngeal carcinoma (NPC) models has seriously hampered research on EBV carcinogenesis and preclinical studies in NPC. Here we report the successful growth of five NPC patient-derived xenografts (PDXs) from fifty-eight attempts of transplantation of NPC specimens into NOD/SCID mice. The take rates for primary and recurrent NPC are 4.9% and 17.6%, respectively. Successful establishment of a new EBV-positive NPC cell line, NPC43, is achieved directly from patient NPC tissues by including Rho-associated coiled-coil containing kinases inhibitor (Y-27632) in culture medium. Spontaneous lytic reactivation of EBV can be observed in NPC43 upon withdrawal of Y-27632. Whole-exome sequencing (WES) reveals a close similarity in mutational profiles of these NPC PDXs with their corresponding patient NPC. Whole-genome sequencing (WGS) further delineates the genomic landscape and sequences of EBV genomes in these newly established NPC models, which supports their potential use in future studies of NPC.

Diabetes and obesity are complex diseases associated with insulin resistance and fatty liver. The latter is characterized by dysregulation of the Akt, AMP-activated protein kinase (AMPK), and IGF-I pathways and expression of microRNAs (miRNAs). In China, multicomponent traditional Chinese medicine (TCM) has been used to treat diabetes for centuries. In this study, we used a three-herb, berberine-containing TCM to treat male Zucker diabetic fatty rats. TCM showed sustained glucose-lowering effects for 1 week after a single-dose treatment. Two-week treatment attenuated insulin resistance and fatty degeneration, with hepatocyte regeneration lasting for 1 month posttreatment. These beneficial effects persisted for 1 year after 1-month treatment. Two-week treatment with TCM was associated with activation of AMPK, Akt, and insulin-like growth factor-binding protein (IGFBP)1 pathways, with downregulation of miR29-b and expression of a gene network implicated in cell cycle, intermediary, and NADPH metabolism with normalization of CYP7a1 and IGFBP1 expression. These concerted changes in mRNA, miRNA, and proteins may explain the sustained effects of TCM in favor of cell survival, increased glucose uptake, and lipid oxidation/catabolism with improved insulin sensitivity and liver regeneration. These novel findings suggest that multicomponent TCM may be a useful tool to unravel genome regulation and expression in complex diseases.

Transcription factors function by binding different classes of regulatory elements. The Encyclopedia of DNA Elements (ENCODE) project has recently produced binding data for more than 100 transcription factors from about 500 ChIP-seq experiments in multiple cell types. While this large amount of data creates a valuable resource, it is nonetheless overwhelmingly complex and simultaneously incomplete since it covers only a small fraction of all human transcription factors.