PAX9 regulates squamous cell differentiation and carcinogenesis in the oro-oesophageal epithelium

Epithelial and mesenchymal compartments of the developing bladder and urethra display spatially distinct gene expression patterns

The lower urinary tract is comprised of the bladder and urethra and develops from the cloaca, a transient endoderm-derived structure formed from the caudal hindgut. After cloacal septation to form the urogenital sinus and anorectal tract, the bladder gradually develops from the anterior portion of the urogenital sinus while the urethra elongates distally into the genital tubercle. The bladder is a target for regenerative and reconstructive therapies but engineering an impermeable bladder epithelial lining has proven challenging. Urethral epithelial function, including its role as an active immune barrier, is poorly studied and neglected in regenerative therapy. A deeper understanding of epithelial patterning of the urogenital sinus by the surrounding mesenchyme, also accounting for sex-specific differences, can inform regenerative therapies. In this study, we identified spatially distinct genes in the epithelial and mesenchymal compartments of the developing mouse bladder and urethra that could be potential drivers of patterning in the lower urinary tract. Our data revealed spatially restricted domains of transcription factor expression in the epithelium that corresponded with bladder or urethra-specific differentiation. Additionally, we identified the genes Wnt2, Klf4 and Pitx2 that localize to the mesenchyme of the developing bladder and could be potential drivers of bladder differentiation. Our data revealed an increase in the expression of several chemokine genes including Cx3cl1 and Cxcl14 in the developing urethral epithelium that correlated with an increase in epithelial-associated macrophages in the urethra. A survey of sex-specific differences in epithelial and mesenchymal compartments revealed several differentially expressed genes between the male and female urethra but few sex-specific differences in bladder. By comparing spatially distinct gene expression in the developing lower urinary tract, our study provides insights into the divergent differentiation trajectories of the fetal bladder and urethra that establish their adult functions.

Novel DNA methylation markers for early detection of gastric cardia adenocarcinoma and esophageal squamous cell carcinoma

Gastric cardia adenocarcinoma (GCA) and esophageal squamous cell carcinoma (ESCC) present significant health challenges in China, often diagnosed at advanced stages with poor prognoses. However, effective biomarkers for early detection remain elusive. This study aimed to integrate methylome and transcriptome data to identify DNA methylation markers for the early detection of GCA and ESCC. In the discovery stage, we conducted Infinium MethylationEPIC array analysis on 36 paired GCA and non-tumor adjacent tissues (NAT), identifying differentially methylated CpG sites (DMCs) between GCA/ESCC and NAT through combined analyses of in-house and publicly available data. In the validation stage, targeted pyrosequencing and quantitative real-time RT-PCR were performed on paired tumor and NAT samples from 50 GCA and 50 ESCC patients. In the application stage, an independent set of 438 samples, including GCA, ESCC, high- and low-grade dysplasia (HGD/LGD), and normal controls, was tested for selected DMCs using pyrosequencing. Our analysis validated three GCA-specific, two ESCC-specific, and one tumor-shared DMCs, exhibiting significant hypermethylation and decreased expression of target genes in tumor samples compared with NAT. Leveraging these DMCs, we developed a GCA-specific 4-marker panel (cg27284428, cg11798358, cg07880787, and cg00585116) with an area under the receiver operating characteristic curve (AUC) of 0.917, effectively distinguishing between cardia HGD/GCA patients and cardia LGD/normal controls. Similarly, an ESCC-specific 3-marker panel (cg14633892, cg04415798, and cg00585116) achieved an AUC of 0.865 in distinguishing esophageal HGD/ESCC cases. Furthermore, integrating cg00585116, age, and alcohol consumption yielded a tumor-shared logistic model with good discrimination for two cancer/HGD (AUC, 0.767; 95% confidence interval, 0.720-0.813). The mean AUC of the model after 5-fold cross-validation was 0.764. In summary, our study identifies novel DNA methylation markers capable of accurately distinguishing GCA/ESCC and HGD from LGD and normal controls. These findings offer promising prospects for targeted DNA methylation assays in future minimally invasive cancer screening methods.

Onco-Ontogeny of Squamous Cell Cancer of the First Pharyngeal Arch Derivatives

Head and neck squamous cell carcinoma (H&NSCC) is an anatomic, biological, and genetic complex disease. It involves more than 1000 genes implied in its oncogenesis; for this review, we limit our search and description to the genes implied in the onco-ontogeny of the derivates from the first pharyngeal arch during embryo development. They can be grouped as transcription factors and signaling molecules (that act as growth factors that bind to receptors). Finally, we propose the term embryo-oncogenesis to refer to the activation, reactivation, and use of the genes involved in the embryo's development during the oncogenesis or malignant tumor invasion and metastasis events as part of an onco-ontogenic inverse process.

Upregulation and functional roles of miR-450b in canine oral melanoma

Canine oral melanoma (COM) is a common and highly aggressive disease with the potential to model human melanomas. Dysregulated microRNAs represent an interesting line of research for COM because they are implicated in tumor progression. One example is miR-450b, which has been investigated for its molecular mechanisms and biological functions in multiple human cancers, but not human or canine melanoma. Here, we aimed to investigate miR-450b as a potential diagnostic biomarker of COM and its functional roles in metastatic and non-metastatic forms of the disease. We investigated the expression of miR-450b and its target mRNA genes in clinical (tumor tissue and plasma) samples and metastatic and primary-tumor cell lines. Knockdown and overexpression experiments were performed to determine the influence of miR-450b on cell proliferation, migration, colony formation, and apoptosis. miR-450b was significantly upregulated in COM and differentiated between metastatic and non-metastatic tumors, and its potential as a biomarker of metastatic and non-metastatic COM was further confirmed in ROC analysis. miR-450b knockdown promoted cell proliferation, migration, and clonogenicity and inhibited apoptosis, whereas its overexpression yielded the reverse pattern. miR-450b directly binds 3' UTR of PAX9 mRNA and modulates its function leading to BMP4 downregulation and MMP9 upregulation at the transcript level. Furthermore, we surmised that miR-450b activates the Wnt signaling pathway based on gene ontology and enrichment analyses. We concluded that miR-450b has the potential as a diagnostic biomarker and could be a target candidate for COM treatment.

Stage II oesophageal carcinoma: peril in disguise associated with cellular reprogramming and oncogenesis regulated by pseudogenes

INTRODUCTION

Pseudogenes have been implicated for their role in regulating cellular differentiation and organismal development. However, their role in promoting cancer-associated differentiation has not been well-studied. This study explores the tumour landscape of oesophageal carcinoma to identify pseudogenes that may regulate events of differentiation to promote oncogenic transformation.

MATERIALS AND METHOD

De-regulated differentiation-associated pseudogenes were identified using DeSeq2 followed by 'InteractiVenn' analysis to identify their expression pattern. Gene expression dependent and independent enrichment analyses were performed with GSEA and ShinyGO, respectively, followed by quantification of cellular reprogramming, extent of differentiation and pleiotropy using three unique metrics. Stage-specific gene regulatory networks using Bayesian Network Splitting Average were generated, followed by network topology analysis. MEME, STREME and Tomtom were employed to identify transcription factors and miRNAs that play a regulatory role downstream of pseudogenes to initiate cellular reprogramming and further promote oncogenic transformation. The patient samples were stratified based on the expression pattern of pseudogenes, followed by GSEA, mutation analysis and survival analysis using GSEA, MAF and 'survminer', respectively.

RESULTS

Pseudogenes display a unique stage-wise expression pattern that characterizes stage II (SII) ESCA with a high rate of cellular reprogramming, degree of differentiation and pleiotropy. Gene regulatory network and associated topology indicate high robustness, thus validating high pleiotropy observed for SII. Pseudogene-regulated expression of SOX2, FEV, PRRX1 and TFAP2A in SII may modulate cellular reprogramming and promote oncogenesis. Additionally, patient stratification-based mutational analysis in SII signifies APOBEC3A (A3A) as a potential hallmark of homeostatic mutational events of reprogrammed cells which in addition to de-regulated APOBEC3G leads to distinct events of hypermutations. Further enrichment analysis for both cohorts revealed the critical role of combinatorial expression of pseudogenes in cellular reprogramming. Finally, survival analysis reveals distinct genes that promote poor prognosis in SII ESCA and patient-stratified cohorts, thus providing valuable prognostic bio-markers along with markers of differentiation and oncogenesis for distinct landscapes of pseudogene expression.

CONCLUSION

Pseudogenes associated with the events of differentiation potentially aid in the initiation of cellular reprogramming to facilitate oncogenic transformation, especially during SII ESCA. Despite a better overall survival of SII, patient stratification reveals combinatorial de-regulation of pseudogenes as a notable marker for a high degree of cellular differentiation with a unique mutational landscape.

Expression of paired box 9 defines an aggressive subset of lung adenocarcinoma preferentially occurring in smokers

AIMS

A distinct subset of lung adenocarcinomas (LADs), arising from a series of peripheral lung cells defined as the terminal respiratory unit (TRU), is characterised by thyroid transcription factor 1 (TTF-1) expression. The clinical relevance of transcription factors (TFs) other than TTF-1 remains unknown in LAD and was explored in the present study.

METHODS AND RESULTS

Seventy-one LAD samples were subjected to high-throughput transcriptome screening of LAD using cap analysis gene expression (CAGE) sequencing data; CAGE provides genome-wide expression levels of the transcription start sites (TSSs). In total, 1083 invasive LAD samples were subjected to immunohistochemical examination for paired box 9 (PAX9) and TTF-1 expression levels. PAX9 is an endoderm development-associated TF that most strongly and inversely correlates with the expression of TTF-1 TSS subsets. Immunohistochemically, PAX9 expression was restricted to the nuclei of ciliated epithelial and basal cells in the bronchi and bronchioles and the nuclei of epithelial cells of the bronchial glands; moreover, PAX9 expression was observed in 304 LADs (28%). PAX9-positive LADs were significantly associated with heavy smoking, non-lepidic subtype, EGFR wild-type tumours and PD-L1 expression (all P < 0.0001). All these characteristics were opposite to those of TRU-type LADs with TTF-1 expression. PAX9 expression was an independent prognostic factor for decreased overall survival (P = 0.022).

CONCLUSIONS

Our results revealed that PAX9 expression defines an aggressive subset of LADs preferentially occurring in smokers that may arise from bronchial or bronchiolar cells.

Homeobox Genes in Odontogenic Lesions: A Scoping Review

BACKGROUND

Homeobox genes play crucial roles in tooth morphogenesis and development and thus mutations in homeobox genes cause developmental disorders such as odontogenic lesions. The aim of this scoping review is to identify and compile available data from the literatures on the topic of homeobox gene expression in odontogenic lesions.

METHOD

An electronic search to collate all the information on studies on homeobox gene expression in odontogenic lesions was carried out in four databases (PubMed, EBSCO host, Web of Science and Cochrane Library) with selected keywords. All papers which reported expression of homeobox genes in odontogenic lesions were considered.

RESULTS

A total of eleven (11) papers describing expression of homeobox genes in odontogenic lesions were identified. Methods of studies included next generation sequencing, microarray analysis, RT-PCR, Western blotting, in situ hybridization, and immunohistochemistry. The homeobox reported in odontogenic lesions includes LHX8 and DLX3 in odontoma; PITX2, MSX1, MSX2, DLX, DLX2, DLX3, DLX4, DLX5, DLX6, ISL1, OCT4 and HOX C in ameloblastoma; OCT4 in adenomatoid odontogenic tumour; PITX2 and MSX2 in primordial odontogenic tumour; PAX9 and BARX1 in odontogenic keratocyst; PITX2, ZEB1 and MEIS2 in ameloblastic carcinoma while there is absence of DLX2, DLX3 and MSX2 in clear cell odontogenic carcinoma.

CONCLUSIONS

This paper summarized and reviews the possible link between homeobox gene expression in odontogenic lesions. Based on the current available data, there are insufficient evidence to support any definite role of homeobox gene in odontogenic lesions.

Intrinsic Networks Regulating Tissue Repair: Comparative Studies of Oral and Skin Wound Healing

Wound repair is a systematic biological program characterized by four overlapping phases: hemostasis, inflammation, proliferation, and remodeling. Notwithstanding differences between species and distinct anatomical sites, the fundamental phases in the wound healing process are conserved among mammalian species. Oral wound healing is defined as an ideal wound healing model because it resolves rapidly and without scar formation. Understanding the regulation and contribution of the different molecular events that drive rapid wound healing in oral mucosa compared with those of the skin will help us define how these lesions heal more efficiently and may provide new therapeutic strategies that can be translated to the clinical settings for patients with chronic nonhealing wounds. Although all epithelial tissues have remarkable ability for tissue repair, the efficiency of such repair differs between epithelia (oral mucosa vs. cutaneous). This prompts the long-standing, fundamental biological and clinically relevant questions as to why and how does the oral mucosa achieve its enhanced wound healing capacity. In this review, we focus on (1) distinct innate wound healing capabilities of the oral mucosa, (2) lessons learned from comparative transcriptomic studies of oral mucosa versus skin, and (3) translation of findings to therapeutics for enhanced wound healing.

PAX9 in Cancer Development

Paired box 9 (PAX9) is a transcription factor of the PAX family functioning as both a transcriptional activator and repressor. Its functional roles in the embryonic development of various tissues and organs have been well studied. However, its roles and molecular mechanisms in cancer development are largely unknown. Here, we review the current understanding of PAX9 expression, upstream regulation of PAX9, and PAX9 downstream events in cancer development. Promoter hypermethylation, promoter SNP, microRNA, and inhibition of upstream pathways (e.g., NOTCH) result in PAX9 silencing or downregulation, whereas gene amplification and an epigenetic axis upregulate PAX9 expression. PAX9 may contribute to carcinogenesis through dysregulation of its transcriptional targets and related molecular pathways. In summary, extensive studies on PAX9 in its cellular and tissue contexts are warranted in various cancers, in particular, HNSCC, ESCC, lung cancer, and cervical SCC.

Computational identification of preneoplastic cells displaying high stemness and risk of cancer progression

Evidence points towards the differentiation state of cells as a marker of cancer risk and progression. Measuring the differentiation state of single cells in a preneoplastic population could thus enable novel strategies for early detection and risk prediction. Recent maps of somatic mutagenesis in normal tissues from young healthy individuals have revealed cancer driver mutations, indicating that these do not correlate well with differentiation state and that other molecular events also contribute to cancer development. We hypothesized that the differentiation state of single cells can be measured by estimating the regulatory activity of the transcription factors (TFs) that control differentiation within that cell lineage. To this end, we present a novel computational method called CancerStemID that estimates a stemness index of cells from single-cell RNA-Seq data. CancerStemID is validated in two human esophageal squamous cell carcinoma (ESCC) cohorts, demonstrating how it can identify undifferentiated preneoplastic cells whose transcriptomic state is overrepresented in invasive cancer. Spatial transcriptomics and whole-genome bisulfite sequencing demonstrated that differentiation activity of tissue-specific TFs was decreased in cancer cells compared to the basal cell-of-origin layer and established that differentiation state correlated with differential DNA methylation at the promoters of these TFs, independently of underlying NOTCH1 and TP53 mutations. The findings were replicated in a mouse model of ESCC development, and the broad applicability of CancerStemID to other cancer-types was demonstrated. In summary, these data support an epigenetic stem-cell model of oncogenesis and highlight a novel computational strategy to identify stem-like preneoplastic cells that undergo positive selection.

PAX9 reactivation by inhibiting DNA methyltransferase triggers antitumor effect in oral squamous cell carcinoma

Aberrant DNA hypermethylation is associated with oral carcinogenesis. Procaine, a local anesthetic, is a DNA methyltransferase (DNMT) inhibitor that activates anticancer mechanisms. However, its effect on silenced tumor suppressor gene (TSG) activation and its biological role in oral squamous cell carcinoma (OSCC) remain unknown. Here, we report procaine inhibited DNA methylation by suppressing DNMT activity and increased the expression of PAX9, a differentiation gene in OSCC cells. Interestingly, the reactivation of PAX9 by procaine found to inhibit cell growth and trigger apoptosis in OSCC in vitro and in vivo. Likely, the enhanced PAX9 expression after exposure to procaine controls stemness and differentiation through the autophagy-dependent pathway in OSCC cells. PAX9 inhibition abrogated procaine-induced apoptosis, autophagy, and inhibition of stemness. In OSCC cells, procaine improved anticancer drug sensitivity through PAX9, and its deficiency significantly blunted the anticancer drug sensitivity mediated by procaine. Additionally, NRF2 activation by procaine facilitated the antitumor response of PAX9, and pharmacological inhibition of NRF2 by ML385 reduced death and prevented the decrease in the orosphere-forming potential of OSCC cells. Furthermore, procaine promoted antitumor activity in FaDu xenografts in athymic nude mice, and immunohistochemistry data showed that PAX9 expression was significantly enhanced in the procaine group compared to the vehicle control. In conclusion, PAX9 reactivation in response to DNMT inhibition could trigger a potent antitumor mechanism to provide a new therapeutic strategy for OSCC.

PAX9 functions as a tumor suppressor gene for cervical cancer via modulating cell proliferation and apoptosis

To investigate the effect of PAX9 on the progression of cervical cancer (CC). PAX9 expression was quantified in CC tissues and adjacent normal tissues, as well as human CC cell lines and human cervical epithelial cells (HCerEpiC). PAX9-overexpression lentiviral vectors were transfected into CC cell lines, followed by the measurement of proliferation and apoptosis and the quantification of apoptosis-related proteins. In vivo, mice were subcutaneously injected with CaSki cells transfected with PAX9-overexpression lentiviral vectors and control vectors. Then, the volume and weight of tumors were measured followed by hematoxylin and eosin (HE) staining, terminal deoxynucleotidyl transferase dUTP nick end labeling (TUNEL) staining, and immunohistochemistry. PAX9 expression in the CC tissues was lower than that in the adjacent normal tissues, which was correlated with the FIGO stage, tumor size, infiltration depth, parametrium invasion, lympho-vascular space invasion tumor-positive lymph nodes, and prognosis. Furthermore, PAX9 in CC cell lines was also lower than in HCerEpiC. PAX9 inhibits the CC cell proliferation and promotes the apoptosis, with the up-regulations of caspase-3, poly(ADP-ribose) polymerase (PARP), and Bax and the down-regulation of Bcl-2. In vivo experiments demonstrated that in the PAX9 group, the tumor weight and volume were lower than those in the vector group accompanying the decreased Ki-67, cleaved-caspase-3, and Bax expressions and the increased TUNEL and Bcl-2 expression. PAX9 was lowly expressed in the CC tissues and associated with the clinicopathological characteristics and prognosis. PAX9 could inhibit proliferation of CC cell lines and promote the apoptosis, thus suppressing the tumor growth in vivo, indicating its potential therapeutic role for CC treatment.

Cancer CRC: A Comprehensive Cancer Core Transcriptional Regulatory Circuit Resource and Analysis Platform

A core transcriptional regulatory circuit (CRC) is a group of interconnected auto-regulating transcription factors (TFs) that form loops and can be identified by super-enhancers (SEs). Studies have indicated that CRCs play an important role in defining cellular identity and determining cellular fate. Additionally, core TFs in CRCs are regulators of cell-type-specific transcriptional regulation. However, a global view of CRC properties across various cancer types has not been generated. Thus, we integrated paired cancer ATAC-seq and H3K27ac ChIP-seq data for specific cell lines to develop the Cancer CRC (http://bio.liclab.net/Cancer_crc/index.html). This platform documented 94,108 cancer CRCs, including 325 core TFs. The cancer CRC also provided the “SE active core TFs analysis” and “TF enrichment analysis” tools to identify potentially key TFs in cancer. In addition, we performed a comprehensive analysis of core TFs in various cancer types to reveal conserved and cancer-specific TFs.

PAX9 Determines Epigenetic State Transition and Cell Fate in Cancer

Abnormalities in genetic and epigenetic modifications can lead to drastic changes in gene expression profiles that are associated with various cancer types. Small cell lung cancer (SCLC) is an aggressive and deadly form of lung cancer with limited effective therapies currently available. By utilizing a genome-wide CRISPR-Cas9 dropout screen in SCLC cells, we identified paired box protein 9 (PAX9) as an essential factor that is overexpressed in human malignant SCLC tumor samples and is transcriptionally driven by the BAP1/ASXL3/BRD4 epigenetic axis. Genome-wide studies revealed that PAX9 occupies distal enhancer elements and represses gene expression by restricting enhancer activity. In multiple SCLC cell lines, genetic depletion of PAX9 led to significant induction of a primed-active enhancer transition, resulting in increased expression of a large number of neural differentiation and tumor-suppressive genes. Mechanistically, PAX9 interacted and cofunctioned with the nucleosome remodeling and deacetylase (NuRD) complex at enhancers to repress nearby gene expression, which was reversed by pharmacologic HDAC inhibition. Overall, this study provides mechanistic insight into the oncogenic function of the PAX9/NuRD complex epigenetic axis in human SCLC and suggests that reactivation of primed enhancers may have potential therapeutic efficacy in treating SCLC expressing high levels of PAX9. SIGNIFICANCE: A genome-wide screen in small cell lung cancer reveals PAX9/NuRD-mediated epigenetic enhancer silencing and tumor progression, supporting the development of novel personalized therapeutic approaches targeting the PAX9-regulated network.

ANANSE: an enhancer network-based computational approach for predicting key transcription factors in cell fate determination

Proper cell fate determination is largely orchestrated by complex gene regulatory networks centered around transcription factors. However, experimental elucidation of key transcription factors that drive cellular identity is currently often intractable. Here, we present ANANSE (ANalysis Algorithm for Networks Specified by Enhancers), a network-based method that exploits enhancer-encoded regulatory information to identify the key transcription factors in cell fate determination. As cell type-specific transcription factors predominantly bind to enhancers, we use regulatory networks based on enhancer properties to prioritize transcription factors. First, we predict genome-wide binding profiles of transcription factors in various cell types using enhancer activity and transcription factor binding motifs. Subsequently, applying these inferred binding profiles, we construct cell type-specific gene regulatory networks, and then predict key transcription factors controlling cell fate transitions using differential networks between cell types. This method outperforms existing approaches in correctly predicting major transcription factors previously identified to be sufficient for trans-differentiation. Finally, we apply ANANSE to define an atlas of key transcription factors in 18 normal human tissues. In conclusion, we present a ready-to-implement computational tool for efficient prediction of transcription factors in cell fate determination and to study transcription factor-mediated regulatory mechanisms. ANANSE is freely available at https://github.com/vanheeringen-lab/ANANSE.

FGFR2 maintains cancer cell differentiation via AKT signaling in esophageal squamous cell carcinoma

Fibroblast growth factors (FGFs) and their receptors (FGFRs) are important for signaling to maintain cancer stem-like cells (CSCs) in esophageal squamous cell carcinoma (ESCC). However, which FGF receptor, 1, 2, 3, 4, and L1, is essential or whether FGFRs have distinct different roles in ESCC-CSCs is still in question. This study shows that FGFR2, particularly the IIIb isoform, is highly expressed in non-CSCs. Non-CSCs have an epithelial phenotype, and such cells are more differentiated in ESCC. Further, FGFR2 induces keratinocyte differentiation through AKT but not MAPK signaling and diminishes CSC populations. Conversely, knockdown of FGFR2 induces epithelial-mesenchymal transition (EMT) and enriches CSC populations in ESCC. Finally, data analysis using The Cancer Genome Atlas (TCGA) dataset shows that expression of FGFR2 significantly correlated with cancer cell differentiation in clinical ESCC samples. The present study shows that each FGFR has a distinct role and FGFR2-AKT signaling is a key driver of keratinocyte differentiation in ESCC. Activation of FGFR2-AKT signaling could be a future therapeutic option targeting CSC in ESCC.

Genome-Wide DNA Methylation Profiling of Esophageal Squamous Cell Carcinoma from Global High-Incidence Regions Identifies Crucial Genes and Potential Cancer Markers

Epigenetic mechanisms such as aberrant DNA methylation (DNAme) are known to drive esophageal squamous cell carcinoma (ESCC), yet they remain poorly understood. Here, we studied tumor-specific DNAme in ESCC cases from nine high-incidence countries of Africa, Asia, and South America. Infinium MethylationEPIC array was performed on 108 tumors and 51 normal tissues adjacent to the tumors (NAT) in the discovery phase, and targeted pyrosequencing was performed on 132 tumors and 36 NAT in the replication phase. Top genes for replication were prioritized by weighting methylation results using RNA-sequencing data from The Cancer Genome Atlas and GTEx and validated by qPCR. Methylome analysis comparing tumor and NAT identified 6,796 differentially methylated positions (DMP) and 866 differential methylated regions (DMR), with a 30% methylation (Δβ) difference. The majority of identified DMPs and DMRs were hypermethylated in tumors, particularly in promoters and gene-body regions of genes involved in transcription activation. The top three prioritized genes for replication, PAX9, SIM2, and THSD4, had similar methylation differences in the discovery and replication sets. These genes were exclusively expressed in normal esophageal tissues in GTEx and downregulated in tumors. The specificity and sensitivity of these DNAme events in discriminating tumors from NAT were assessed. Our study identified novel, robust, and crucial tumor-specific DNAme events in ESCC tumors across several high-incidence populations of the world. Methylome changes identified in this study may serve as potential targets for biomarker discovery and warrant further functional characterization. SIGNIFICANCE: This largest genome-wide DNA methylation study on ESCC from high-incidence populations of the world identifies functionally relevant and robust DNAme events that could serve as potential tumor-specific markers. GRAPHICAL ABSTRACT: http://cancerres.aacrjournals.org/content/canres/81/10/2612/F1.large.jpg.

The clinical significance and correlative signaling pathways of paired box gene 9 in development and carcinogenesis

Paired box 9 (PAX9) gene belongs to the PAX family, which encodes a family of metazoan transcription factors documented by a conserved DNA binding paired domain 128-amino-acids, critically essential for physiology and development. It is primarily expressed in embryonic tissues, such as the pharyngeal pouch endoderm, somites, neural crest-derived mesenchyme, and distal limb buds. PAX9 plays a vital role in craniofacial development by maintaining the odontogenic potential, mutations, and polymorphisms associated with the risk of tooth agenesis, hypodontia, and crown size in dentition. The loss-of-function of PAX9 in the murine model resulted in a short life span due to the arrest of cleft palate formation and skeletal abnormalities. According to recent studies, the PAX9 gene has a significant role in maintaining squamous cell differentiation, odontoblast differentiation of pluripotent stem cells, deregulation of which is associated with tumor initiation, and malignant transformation. Moreover, PAX9 contributes to promoter hypermethylation and alcohol- induced oro-esophageal squamous cell carcinoma mediated by downregulation of differentiation and apoptosis. Likewise, PAX9 activation is also reported to be associated with drug sensitivity. In summary, this current review aims to understand PAX9 function in the regulation of development, differentiation, and carcinogenesis, along with the underlying signaling pathways for possible cancer therapeutics.

Alcohol drinking inhibits NOTCH-PAX9 signaling in esophageal squamous epithelial cells

Alcohol drinking has been established as a major risk factor for esophageal diseases. Our previous study showed that ethanol exposure inhibited PAX9 expression in human esophageal squamous epithelial cells in vitro and in vivo. In this study, we aimed to investigate the molecular pathways through which alcohol drinking suppresses PAX9 in esophageal squamous epithelial cells. We first demonstrated the inhibition of NOTCH by ethanol exposure in vitro. NOTCH regulated PAX9 expression in KYSE510 and KYSE410 cells in vitro and in vivo. RBPJ and NOTCH intracellular domain (NIC) D1 ChIP-PCR confirmed Pax9 as a direct downstream target of NOTCH signaling in mouse esophagus. NOTCH inhibition by alcohol drinking was further validated in mouse esophagus and human tissue samples. In conclusion, ethanol exposure inhibited NOTCH signaling and thus suppressed PAX9 expression in esophageal squamous epithelial cells in vitro and in vivo. Our data support a novel mechanism of alcohol-induced esophageal injury through the inhibition of NOTCH-PAX9 signaling. © 2020 The Pathological Society of Great Britain and Ireland. Published by John Wiley & Sons, Ltd.

NOTCH and Esophageal Squamous Cell Carcinoma

Esophageal squamous cell carcinoma (ESCC) is a deadly disease that requires extensive research on its mechanisms, prevention, and therapy. Recent studies have shown that NOTCH mutations are commonly seen in human ESCC. This chapter summarizes our current understanding of the NOTCH pathway in normal esophagus and in ESCC. In normal esophagus, NOTCH pathway regulates the development of esophageal squamous epithelium, in particular, squamous differentiation. Exposure to extrinsic and intrinsic factors, such as gastroesophageal reflux, alcohol drinking, and inflammation, downregulates the NOTCH pathway and thus inhibits squamous differentiation of esophageal squamous epithelial cells. In ESCC, NOTCH plays a dual role as both a tumor suppressor pathway and an oncogenic pathway. In summary, further studies are warranted to develop NOTCH activators for the prevention of ESCC and NOTCH inhibitors for targeted therapy of a subset of ESCC with activated NOTCH pathway.

Identification of key regulators in parathyroid adenoma using an integrative gene network analysis

Parathyroid adenoma (PA) is marked by a certain benign outgrowth in the surface of parathyroid glands. The transcriptome analysis of parathyroid adenomas can provide a deep insight into actively expressed genes and transcripts. Hence, we analyzed and compared the gene expression profiles of parathyroid adenomas and healthy parathyroid gland tissues from Gene Expression Omnibus (GEO) database. We identified a total of 280 differentially expressed genes (196 up-regulated, 84 down-regulated), which are involved in a wide array of biological processes. We further constructed a gene interaction network and analyzed its topological properties to know the network structure and its hidden mechanism. This will help to understand the molecular mechanisms underlying parathyroid adenoma development. We thus identified 13 key regulators (PRPF19, SMC3, POSTN, SNIP1, EBF1, MEIS2, PAX9, SCUBE2, WNT4, ARHGAP10, DOCK5, CAV1 and VSIR), which are deep-rooted from top to bottom in the gene interaction network forming a backbone for the network. The structural features of the network are probably maintained by crosstalk between important genes within the network along with associated functional modules.Thus, gene-expression profiling and network approach could be used to provide an independent platform to glen insights from available clinical data.

Modes of genetic adaptations underlying functional innovations in the rumen

The rumen is the hallmark organ of ruminants and hosts a diverse ecosystem of microorganisms that facilitates efficient digestion of plant fibers. We analyzed 897 transcriptomes from three Cetartiodactyla lineages: ruminants, camels and cetaceans, as well as data from ruminant comparative genomics and functional assays to explore the genetic basis of rumen functional innovations. We identified genes with relatively high expression in the rumen, of which many appeared to be recruited from other tissues. These genes show functional enrichment in ketone body metabolism, regulation of microbial community, and epithelium absorption, which are the most prominent biological processes involved in rumen innovations. Several modes of genetic change underlying rumen functional innovations were uncovered, including coding mutations, genes newly evolved, and changes of regulatory elements. We validated that the key ketogenesis rate-limiting gene (HMGCS2) with five ruminant-specific mutations was under positive selection and exhibits higher synthesis activity than those of other mammals. Two newly evolved genes (LYZ1 and DEFB1) are resistant to Gram-positive bacteria and thereby may regulate microbial community equilibrium. Furthermore, we confirmed that the changes of regulatory elements accounted for the majority of rumen gene recruitment. These results greatly improve our understanding of rumen evolution and organ evo-devo in general.

NRF2/ACSS2 axis mediates the metabolic effect of alcohol drinking on esophageal squamous cell carcinoma

Alcohol drinking is a leading risk factor for the development of esophageal squamous cell carcinoma (ESCC). However, the molecular mechanisms of alcohol-associated ESCC remain poorly understood. One of the most commonly mutated genes in ESCC is nuclear factor erythroid 2 like 2 (NFE2L2 or NRF2), which is a critical transcription factor regulating oxidative stress response and drug detoxification. When NRF2 is hyperactive in cancer cells, however, it leads to metabolic reprogramming, cell proliferation, chemoradioresistance, and poor prognosis. In this study, hyperactive NRF2 was found to up-regulate acetyl-CoA synthetase short-chain family members 2 (ACSS2), an enzyme that converts acetate to acetyl-CoA, in ESCC cells and mouse esophagus. We also showed that knockdown of NRF2 or ACSS2 led to decreased ACSS2 expression, which in turn reduced the levels of acetyl-CoA and ATP with or without ethanol exposure. In addition, ethanol exposure enhanced lipid synthesis in ESCC cells. Moreover, we observed a change in the metabolic profile of ESCC cells exposed to ethanol as a result of their NRF2 or ACSS2 status. We further showed that ACSS2 contributed to the invasive capability of NRF2high ESCC cells exposed to ethanol. In conclusion, the NRF2/ACSS2 axis mediates the metabolic effect of alcohol drinking on ESCC.

Paired box 9 regulates VSMC phenotypic transformation, proliferation, and migration via sonic hedgehog

AIMS

Vascular smooth muscle cells (VSMCs) play a crucial role in the progression of atherosclerosis. Paired box 9 (Pax9) is a member of the Pax gene family which participates in the development of various tissues and organs. However, the effect of Pax9 on atherosclerosis and VSMCs and the underlying mechanisms remain unclear.

MAIN METHODS

Western blotting was performed to assess Pax9 expression in atherosclerosis and VSMCs. Pax9 siRNA and overexpression plasmid were constructed to explore the biological function. Cell proliferation assay, phalloidin staining, and Transwell assay, accompanied by the sonic hedgehog (Shh) signaling pathway antagonist, cyclopamine (5 μM) and agonist, SAG (100 nM), were used to evaluate the VSMC phenotype, proliferation, and migration, as well as explore the associated mechanisms.

KEY FINDINGS

We first discovered Pax9 to be significantly increased in atherosclerotic mice and platelet-derived growth factor-BB (PDGF-BB)-induced VSMCs. Pax9 knockdown inhibited the phenotypic transformation, proliferation, and migration of VSMCs, whereas the opposite effect was observed when Pax9 was overexpressed. Next, we established that Shh was activated in PDGF-BB-induced VSMCs. Moreover, Pax9 overexpression further activated Shh and exacerbated the phenotypic transformation, proliferation, and migration of PDGF-BB-induced VSMCs. These changes were effectively inhibited by treatment with the Shh signaling pathway antagonist. Consistently, Pax9 knockdown down-regulated Shh expression and inhibited the phenotypic transformation, proliferation, and migration of PDGF-BB-induced VSMCs. Treatment with the Shh signaling pathway agonist prevented these changes.

SIGNIFICANCE

Pax9 regulated VSMC phenotypic transformation, proliferation, and migration via Shh, which may represent a novel target for the treatment of atherosclerosis.

Transcriptomic analysis to affirm the regulatory role of long non-coding RNA in horn cancer of Indian zebu cattle breed Kankrej (Bos indicus)

Long non-coding RNA (lncRNA) was previously considered as a non-functional transcript, which now established as part of regulatory elements of biological events such as chromosome structure, remodeling, and regulation of gene expression. The study presented here showed the role of lncRNA through differential expression analysis on cancer-related coding genes in horn squamous cell carcinoma of Indian zebu cattle. A total of 10,360 candidate lncRNAs were identified and further analyzed for its coding potential ability using three tools (CPC, CPAT, and PLEK) that provide 8862 common lncRNAs. Pfam analysis of these common lncRNAs gave 8612 potential candidates for lncRNA differential expression analysis. Differential expression analysis showed a total of 59 significantly differentially expressed genes and 19 lncRNAs. Pearson's correlation analysis was used to identify co-expressed mRNA-lncRNAs to established relation of the regulatory role of lncRNAs in horn cancer. We established a positive relation of seven upregulated (XLOC_000016, XLOC_002198, XLOC_002851, XLOC_ 007383, XLOC_010701, XLOC_010272, and XLOC_011517) and one downregulated (XLOC_011302) lncRNAs with eleven genes that are related to keratin family protein, keratin-associated protein family, cornifelin, corneodesmosin, serpin family protein, and metallothionein that have well-established role in squamous cell carcinoma through cellular communication, cell growth, cell invasion, and cell migration. These biological events were found to be related to the MAPK pathway of cell cycle regulation indicating the role of lncRNAs in manipulating cell cycle regulation during horn squamous cell carcinomas that will be useful in identifying molecular portraits related to the development of horn cancer.

Large-scale ruminant genome sequencing provides insights into their evolution and distinct traits

The ruminants are one of the most successful mammalian lineages, exhibiting morphological and habitat diversity and containing several key livestock species. To better understand their evolution, we generated and analyzed de novo assembled genomes of 44 ruminant species, representing all six Ruminantia families. We used these genomes to create a time-calibrated phylogeny to resolve topological controversies, overcoming the challenges of incomplete lineage sorting. Population dynamic analyses show that population declines commenced between 100,000 and 50,000 years ago, which is concomitant with expansion in human populations. We also reveal genes and regulatory elements that possibly contribute to the evolution of the digestive system, cranial appendages, immune system, metabolism, body size, cursorial locomotion, and dentition of the ruminants.