A vitamin-C-derived DNA modification catalysed by an algal TET homologue.

Methylation of cytosine to 5-methylcytosine (5mC) is a prevalent DNA modification found in many organisms. Sequential oxidation of 5mC by ten-eleven translocation (TET) dioxygenases results in a cascade of additional epigenetic marks and promotes demethylation of DNA in mammals1,2. However, the enzymatic activity and function of TET homologues in other eukaryotes remains largely unexplored. Here we show that the green alga Chlamydomonas reinhardtii contains a 5mC-modifying enzyme (CMD1) that is a TET homologue and catalyses the conjugation of a glyceryl moiety to the methyl group of 5mC through a carbon-carbon bond, resulting in two stereoisomeric nucleobase products. The catalytic activity of CMD1 requires Fe(II) and the integrity of its binding motif His-X-Asp, which is conserved in Fe-dependent dioxygenases3. However, unlike previously described TET enzymes, which use 2-oxoglutarate as a co-substrate4, CMD1 uses L-ascorbic acid (vitamin C) as an essential co-substrate. Vitamin C donates the glyceryl moiety to 5mC with concurrent formation of glyoxylic acid and CO2. The vitamin-C-derived DNA modification is present in the genome of wild-type C. reinhardtii but at a substantially lower level in a CMD1 mutant strain. The fitness of CMD1 mutant cells during exposure to high light levels is reduced. LHCSR3, a gene that is critical for the protection of C. reinhardtii from photo-oxidative damage under high light conditions, is hypermethylated and downregulated in CMD1 mutant cells compared to wild-type cells, causing a reduced capacity for photoprotective non-photochemical quenching. Our study thus identifies a eukaryotic DNA base modification that is catalysed by a divergent TET homologue and unexpectedly derived from vitamin C, and describes its role as a potential epigenetic mark that may counteract DNA methylation in the regulation of photosynthesis.

ZFP57 dictates allelic expression switch of target imprinted genes.

ZFP57 is a master regulator of genomic imprinting. It has both maternal and zygotic functions that are partially redundant in maintaining DNA methylation at some imprinting control regions (ICRs). In this study, we found that DNA methylation was lost at most known ICRs in Zfp57 mutant embryos. Furthermore, loss of ZFP57 caused loss of parent-of-origin-dependent monoallelic expression of the target imprinted genes. The allelic expression switch occurred in the ZFP57 target imprinted genes upon loss of differential DNA methylation at the ICRs in Zfp57 mutant embryos. Specifically, upon loss of ZFP57, the alleles of the imprinted genes located on the same chromosome with the originally methylated ICR switched their expression to mimic their counterparts on the other chromosome with unmethylated ICR. Consistent with our previous study, ZFP57 could regulate the NOTCH signaling pathway in mouse embryos by impacting allelic expression of a few regulators in the NOTCH pathway. In addition, the imprinted Dlk1 gene that has been implicated in the NOTCH pathway was significantly down-regulated in Zfp57 mutant embryos. Our allelic expression switch models apply to the examined target imprinted genes controlled by either maternally or paternally methylated ICRs. Our results support the view that ZFP57 controls imprinted expression of its target imprinted genes primarily through maintaining differential DNA methylation at the ICRs.

TransBorrow: genome-guided transcriptome assembly by borrowing assemblies from different assemblers.

RNA-seq technology is widely used in various transcriptomic studies and provides great opportunities to reveal the complex structures of transcriptomes. To effectively analyze RNA-seq data, we introduce a novel transcriptome assembler, TransBorrow, which borrows the assemblies from different assemblers to search for reliable subsequences by building a colored graph from those borrowed assemblies. Then, by seeding reliable subsequences, a newly designed path extension strategy accurately searches for a transcript-representing path cover over each splicing graph. TransBorrow was tested on both simulated and real data sets and showed great superiority over all the compared leading assemblers.

TransRef enables accurate transcriptome assembly by redefining accurate neo-splicing graphs.

RNA-seq technology is widely employed in various research areas related to transcriptome analyses, and the identification of all the expressed transcripts from short sequencing reads presents a considerable computational challenge. In this study, we introduce TransRef, a new computational algorithm for accurate transcriptome assembly by redefining a novel graph model, the neo-splicing graph, and then iteratively applying a constrained dynamic programming to reconstruct all the expressed transcripts for each graph. When TransRef is utilized to analyze both real and simulated datasets, its performance is notably better than those of several state-of-the-art assemblers, including StringTie2, Cufflinks and Scallop. In particular, the performance of TransRef is notably strong in identifying novel transcripts and transcripts with low-expression levels, while the other assemblers are less effective.

TransCirc: an interactive database for translatable circular RNAs based on multi-omics evidence.

TransCirc (https://www.biosino.org/transcirc/) is a specialized database that provide comprehensive evidences supporting the translation potential of circular RNAs (circRNAs). This database was generated by integrating various direct and indirect evidences to predict coding potential of each human circRNA and the putative translation products. Seven types of evidences for circRNA translation were included: (i) ribosome/polysome binding evidences supporting the occupancy of ribosomes onto circRNAs; (ii) experimentally mapped translation initiation sites on circRNAs; (iii) internal ribosome entry site on circRNAs; (iv) published N-6-methyladenosine modification data in circRNA that promote translation initiation; (v) lengths of the circRNA specific open reading frames; (vi) sequence composition scores from a machine learning prediction of all potential open reading frames; (vii) mass spectrometry data that directly support the circRNA encoded peptides across back-splice junctions. TransCirc provides a user-friendly searching/browsing interface and independent lines of evidences to predicte how likely a circRNA can be translated. In addition, several flexible tools have been developed to aid retrieval and analysis of the data. TransCirc can serve as an important resource for investigating the translation capacity of circRNAs and the potential circRNA-encoded peptides, and can be expanded to include new evidences or additional species in the future.

In vivo CRISPR screening unveils histone demethylase UTX as an important epigenetic regulator in lung tumorigenesis.

Lung cancer is the leading cause of cancer-related death worldwide. Inactivation of tumor suppressor genes (TSGs) promotes lung cancer malignant progression. Here, we take advantage of the clustered regularly interspaced short palindromic repeats (CRISPR)/Cas9-mediated somatic gene knockout in a KrasG12D/+ mouse model to identify bona fide TSGs. From individual knockout of 55 potential TSGs, we identify five genes, including Utx, Ptip, Acp5, Acacb, and Clu, whose knockout significantly promotes lung tumorigenesis. These candidate genes are frequently down-regulated in human lung cancer specimens and significantly associated with survival in patients with lung cancer. Through crossing the conditional Utx knockout allele to the KrasG12D/+ mouse model, we further find that Utx deletion dramatically promotes lung cancer progression. The tumor-promotive effect of Utx knockout in vivo is mainly mediated through an increase of the EZH2 level, which up-regulates the H3K27me3 level. Moreover, the Utx-knockout lung tumors are preferentially sensitive to EZH2 inhibitor treatment. Collectively, our study provides a systematic screening of TSGs in vivo and identifies UTX as an important epigenetic regulator in lung tumorigenesis.

Engineering the Berreman mode in mid-infrared polar materials.

We demonstrate coupling to and control over the broadening and dispersion of a mid-infrared leaky mode, known as the Berreman mode, in samples with different dielectric environments. We fabricate subwavelength films of AlN, a mid-infrared epsilon-near-zero material that supports the Berreman mode, on materials with a weakly negative permittivity, strongly negative permittivity, and positive permittivity. Additionally, we incorporate ultra-thin AlN layers into a GaN/AlN heterostructure, engineering the dielectric environment above and below the AlN. In each of the samples, coupling to the Berreman mode is observed in angle-dependent reflection measurements at wavelengths near the longitudinal optical phonon energy. The measured dispersion of the Berreman mode agrees well with numerical modes. Differences in the dispersion and broadening for the different materials is quantified, including a 13 cm-1 red-shift in the energy of the Berreman mode for the heterostructure sample.

Landscape of transcriptional deregulation in lung cancer.

BACKGROUND: Lung cancer is a very heterogeneous disease that can be pathologically classified into different subtypes including small-cell lung carcinoma (SCLC), lung adenocarcinoma (LUAD), lung squamous cell carcinoma (LUSC) and large-cell carcinoma (LCC). Although much progress has been made towards the oncogenic mechanism of each subtype, transcriptional circuits mediating the upstream signaling pathways and downstream functional consequences remain to be systematically studied. RESULTS: Here we trained a one-class support vector machine (OC-SVM) model to establish a general transcription factor (TF) regulatory network containing 325 TFs and 18724 target genes. We then applied this network to lung cancer subtypes and identified those deregulated TFs and downstream targets. We found that the TP63/SOX2/DMRT3 module was specific to LUSC, corresponding to squamous epithelial differentiation and/or survival. Moreover, the LEF1/MSC module was specifically activated in LUAD and likely to confer epithelial-to-mesenchymal transition, known important for cancer malignant progression and metastasis. The proneural factor, ASCL1, was specifically up-regulated in SCLC which is known to have a neuroendocrine phenotype. Also, ID2 was differentially regulated between SCLC and LUSC, with its up-regulation in SCLC linking to energy supply for fast mitosis and its down-regulation in LUSC linking to the attenuation of immune response. We further described the landscape of TF regulation among the three major subtypes of lung cancer, highlighting their functional commonalities and specificities. CONCLUSIONS: Our approach uncovered the landscape of transcriptional deregulation in lung cancer, and provided a useful resource of TF regulatory network for future studies.

The RNA-binding protein QKI suppresses cancer-associated aberrant splicing.

Lung cancer is the leading cause of cancer-related death worldwide. Aberrant splicing has been implicated in lung tumorigenesis. However, the functional links between splicing regulation and lung cancer are not well understood. Here we identify the RNA-binding protein QKI as a key regulator of alternative splicing in lung cancer. We show that QKI is frequently down-regulated in lung cancer, and its down-regulation is significantly associated with a poorer prognosis. QKI-5 inhibits the proliferation and transformation of lung cancer cells both in vitro and in vivo. Our results demonstrate that QKI-5 regulates the alternative splicing of NUMB via binding to two RNA elements in its pre-mRNA, which in turn suppresses cell proliferation and prevents the activation of the Notch signaling pathway. We further show that QKI-5 inhibits splicing by selectively competing with a core splicing factor SF1 for binding to the branchpoint sequence. Taken together, our data reveal QKI as a critical regulator of splicing in lung cancer and suggest a novel tumor suppression mechanism involving QKI-mediated regulation of the Notch signaling pathway.

NEDD9 promotes lung cancer metastasis through epithelial-mesenchymal transition.

Metastasis is the major cause for high mortality of lung cancer with the underlying mechanisms poorly understood. The scaffolding protein neural precursor cell expressed, developmentally down-regulated 9 (NEDD9) has been identified as a pro-metastasis gene in several types of cancers including melanoma and breast cancer. However, the exact role and related mechanism of NEDD9 in regulating lung cancer metastasis still remain largely unknown. Here, we demonstrate that NEDD9 knockdown significantly inhibits migration, invasion and metastasis of lung cancer cells in vitro and in vivo. The pro-metastasis role of Nedd9 in lung cancer is further supported by studies in mice models of spontaneous cancer metastasis. Moreover, we find that NEDD9 promotes lung cancer cell migration and invasion through the induction of epithelial-mesenchymal transition (EMT) potentially via focal adhesion kinase activation. More importantly, NEDD9 expression inversely correlates with E-cadherin expression in human lung cancer specimens, consistent with the findings from in vitro studies. Taken together, this study highlights that NEDD9 is an important mediator promotes lung cancer metastasis via EMT.

EML4-ALK fusions drive lung adeno-to-squamous transition through JAK-STAT activation.

Human lung adenosquamous cell carcinoma (LUAS), containing both adenomatous and squamous pathologies, exhibits strong cancer plasticity. We find that ALK rearrangement is detectable in 5.1-7.5% of human LUAS, and transgenic expression of EML4-ALK drives lung adenocarcinoma (LUAD) formation initially and squamous transition at late stage. We identify club cells as the main cell-of-origin for squamous transition. Through recapitulating lineage transition in organoid system, we identify JAK-STAT signaling, activated by EML4-ALK phase separation, significantly promotes squamous transition. Integrative study with scRNA-seq and immunostaining identify a plastic cell subpopulation in ALK-rearranged human LUAD showing squamous biomarker expression. Moreover, those relapsed ALK-rearranged LUAD show notable upregulation of squamous biomarkers. Consistently, mouse squamous tumors or LUAD with squamous signature display certain resistance to ALK inhibitor, which can be overcome by combined JAK1/2 inhibitor treatment. This study uncovers strong plasticity of ALK-rearranged tumors in orchestrating phenotypic transition and drug resistance and proposes a potentially effective therapeutic strategy.

Adeno-to-squamous transition drives resistance to KRAS inhibition in LKB1 mutant lung cancer.

KRASG12C inhibitors (adagrasib and sotorasib) have shown clinical promise in targeting KRASG12C-mutated lung cancers; however, most patients eventually develop resistance. In lung patients with adenocarcinoma with KRASG12C and STK11/LKB1 co-mutations, we find an enrichment of the squamous cell carcinoma gene signature in pre-treatment biopsies correlates with a poor response to adagrasib. Studies of Lkb1-deficient KRASG12C and KrasG12D lung cancer mouse models and organoids treated with KRAS inhibitors reveal tumors invoke a lineage plasticity program, adeno-to-squamous transition (AST), that enables resistance to KRAS inhibition. Transcriptomic and epigenomic analyses reveal ΔNp63 drives AST and modulates response to KRAS inhibition. We identify an intermediate high-plastic cell state marked by expression of an AST plasticity signature and Krt6a. Notably, expression of the AST plasticity signature and KRT6A at baseline correlates with poor adagrasib responses. These data indicate the role of AST in KRAS inhibitor resistance and provide predictive biomarkers for KRAS-targeted therapies in lung cancer.

MicroRNA-143 (miR-143) regulates cancer glycolysis via targeting hexokinase 2 gene.

High glycolysis, well known as "Warburg effect," is frequently observed in a variety of cancers. Whether the deregulation of miRNAs contributes to the Warburg effect remains largely unknown. Because miRNA regulates gene expression at both mRNA and protein levels, we constructed a gene functional association network, which allows us to detect the gene activity instead of gene expression, to integratively analyze the microarray data for gene expression and miRNA expression profiling and identify glycolysis-related gene-miRNA pairs deregulated in cancer. Hexokinase 2 (HK2), coding for the first rate-limiting enzyme of glycolysis, is among the top list of genes predicted and potentially regulated by multiple miRNAs including miR-143. Interestingly, miR-143 expression was inversely associated with HK2 protein level but not mRNA level in human lung cancer samples. miR-143, down-regulated by mammalian target of rapamycin activation, reduces glucose metabolism and inhibits cancer cell proliferation and tumor formation through targeting HK2. Collectively, we have not only established a novel methodology for gene-miRNA pair prediction but also identified miR-143 as an essential regulator of cancer glycolysis via targeting HK2.

ANCCA protein expression is a novel independent poor prognostic marker in surgically resected lung adenocarcinoma.

BACKGROUND: AAA+ nuclear coregulator cancer associated (ANCCA) is found to be overexpressed in various cancer types and could play a role in common and fundamental cellular processes. A recent study suggested that ANCCA was a likely driver whose expression explained the behavior of differentially expressed proliferation-related genes in lung adenocarcinoma. However, protein expression of ANCCA in lung adenocarcinoma and its association with clinicopathologic parameters and commonly reported driver mutations remains unexplored. METHODS: ANCCA expression was evaluated by immunohistochemistry in 143 surgically resected lung adenocarcinomas and was correlated with clinicopathologic and molecular variables including adenocarcinoma histologic subtypes, tumor, node, metastasis status, relapse-free survival, overall survival, EGFR mutations, KRAS mutations, HER2 mutations and ALK fusions. RESULTS: Positive ANCCA expression was significantly associated with male sex, smokers, poorly differentiated tumors, nonlepidic predominant subtype, more advanced T stage, lymph nodal metastasis and late disease stage. Cox multivariate analysis revealed that ANCCA-positive expression was an independent predictor of worse relapse-free survival [hazard ratio (HR) 1.736, 95 % confidence interval (CI) 1.075-2.804; P = .024) and overall survival (HR 7.758, 95 % CI 2.955-20.370; P < .001). The addition of ANCCA protein expression to the prognostic model using pathologic stage markedly improved the prognostic accuracy; the concordance index increased from .692 to .788, and the Akaike information criterion decreased from 354.20 to 336.11. CONCLUSIONS: We have identified ANCCA protein expression as a novel independent poor prognostic indicator in lung adenocarcinoma. Prospective studies are warranted to validate its potential prognostic value in combination with the current staging system.

LKB1 inhibits lung cancer progression through lysyl oxidase and extracellular matrix remodeling.

LKB1 loss-of-function mutations, observed in ∼30% of human lung adenocarcinomas, contribute significantly to lung cancer malignancy progression. We show that lysyl oxidase (LOX), negatively regulated by LKB1 through mTOR-HIF-1α signaling axis, mediates lung cancer progression. Inhibition of LOX activity dramatically alleviates lung cancer malignancy progression. Up-regulated LOX expression triggers excess collagen deposition in Lkb1-deficient lung tumors, and thereafter results in enhanced cancer cell proliferation and invasiveness through activation of ß1 integrin signaling. High LOX level and activity correlate with poor prognosis and metastasis. Our findings provide evidence of how LKB1 loss of function promotes lung cancer malignancy through remodeling of extracellular matrix microenvironment, and identify LOX as a potential target for disease treatment in lung cancer patients.

Oxidative stress-triggered Wnt signaling perturbation characterizes the tipping point of lung adeno-to-squamous transdifferentiation.

Lkb1 deficiency confers the Kras-mutant lung cancer with strong plasticity and the potential for adeno-to-squamous transdifferentiation (AST). However, it remains largely unknown how Lkb1 deficiency dynamically regulates AST. Using the classical AST mouse model (Kras LSL-G12D/+;Lkb1flox/flox, KL), we here comprehensively analyze the temporal transcriptomic dynamics of lung tumors at different stages by dynamic network biomarker (DNB) and identify the tipping point at which the Wnt signaling is abruptly suppressed by the excessive accumulation of reactive oxygen species (ROS) through its downstream effector FOXO3A. Bidirectional genetic perturbation of the Wnt pathway using two different Ctnnb1 conditional knockout mouse strains confirms its essential role in the negative regulation of AST. Importantly, pharmacological activation of the Wnt pathway before but not after the tipping point inhibits squamous transdifferentiation, highlighting the irreversibility of AST after crossing the tipping point. Through comparative transcriptomic analyses of mouse and human tumors, we find that the lineage-specific transcription factors (TFs) of adenocarcinoma and squamous cell carcinoma form a "Yin-Yang" counteracting network. Interestingly, inactivation of the Wnt pathway preferentially suppresses the adenomatous lineage TF network and thus disrupts the "Yin-Yang" homeostasis to lean towards the squamous lineage, whereas ectopic expression of NKX2-1, an adenomatous lineage TF, significantly dampens such phenotypic transition accelerated by the Wnt pathway inactivation. The negative correlation between the Wnt pathway and AST is further observed in a large cohort of human lung adenosquamous carcinoma. Collectively, our study identifies the tipping point of AST and highlights an essential role of the ROS-Wnt axis in dynamically orchestrating the homeostasis between adeno- and squamous-specific TF networks at the AST tipping point.

Counteracting lineage-specific transcription factor network finely tunes lung adeno-to-squamous transdifferentiation through remodeling tumor immune microenvironment.

Human lung adenosquamous cell carcinoma (LUAS), containing both adenomatous and squamous pathologies, harbors strong plasticity and is significantly associated with poor prognosis. We established an up-to-date comprehensive genomic and transcriptomic landscape of LUAS in 109 Chinese specimens and demonstrated LUAS development via adeno-to-squamous transdifferentiation. Unsupervised transcriptomic clustering and dynamic network biomarker analysis identified an inflammatory subtype as the critical transition stage during LUAS development. Dynamic dysregulation of the counteracting lineage-specific transcription factors (TFs), containing adenomatous TFs NKX2-1 and FOXA2, and squamous TFs TP63 and SOX2, finely tuned the lineage transition via promoting CXCL3/5-mediated neutrophil infiltration. Genomic clustering identified the most malignant subtype featured with STK11-inactivation, and targeting LSD1 through genetic deletion or pharmacological inhibition almost eradicated STK11-deficient lung tumors. These data collectively uncover the comprehensive molecular landscape, oncogenic driver spectrum and therapeutic vulnerability of Chinese LUAS.

HLungDB: an integrated database of human lung cancer research.

The human lung cancer database (HLungDB) is a database with the integration of the lung cancer-related genes, proteins and miRNAs together with the corresponding clinical information. The main purpose of this platform is to establish a network of lung cancer-related molecules and to facilitate the mechanistic study of lung carcinogenesis. The entries describing the relationships between molecules and human lung cancer in the current release were extracted manually from literatures. Currently, we have collected 2585 genes and 212 miRNA with the experimental evidences involved in the different stages of lung carcinogenesis through text mining. Furthermore, we have incorporated the results from analysis of transcription factor-binding motifs, the promoters and the SNP sites for each gene. Since epigenetic alterations also play an important role in lung carcinogenesis, genes with epigenetic regulation were also included. We hope HLungDB will enrich our knowledge about lung cancer biology and eventually lead to the development of novel therapeutic strategies. HLungDB can be freely accessed at http://www.megabionet.org/bio/hlung.

Loss of TSC1/TSC2 sensitizes immune checkpoint blockade in non-small cell lung cancer.

Tuberous sclerosis complex subunit 1 (TSC1) and 2 (TSC2) are frequently mutated in non-small cell lung cancer (NSCLC), however, their effects on antitumor immunity remained unexplored. A CRISPR screening in murine KrasG12D/Trp53-/- (KP) model identified Tsc1 and Tsc2 as potent regulators of programmed cell death ligand 1 (Pd-l1) expression in vitro and sensitivity to anti-programmed cell death receptor 1 (PD-1) treatment in vivo. TSC1 or TSC2 knockout (KO) promoted the transcriptional and membrane expression of PD-L1 in cell lines. TSC2-deficient tumors manifested an inflamed microenvironment in patient samples and The Cancer Genome Atlas dataset. In syngeneic murine models, KP-Tsc2-KO tumors showed notable response to anti-PD-1 antibody treatment, but Tsc2-wild-type tumors did not. Patients with TSC1/TSC2-mutant NSCLC receiving immune checkpoint blockade (ICB) had increased durable clinical benefit and survival. Collectively, TSC1/TSC2 loss defines a distinct subtype of NSCLC characterized as inflamed tumor microenvironment and superior sensitivity to ICB.

c-CSN: Single-cell RNA Sequencing Data Analysis by Conditional Cell-specific Network.

The rapid advancement of single-cell technologies has shed new light on the complex mechanisms of cellular heterogeneity. However, compared to bulk RNA sequencing (RNA-seq), single-cell RNA-seq (scRNA-seq) suffers from higher noise and lower coverage, which brings new computational difficulties. Based on statistical independence, cell-specific network (CSN) is able to quantify the overall associations between genes for each cell, yet suffering from a problem of overestimation related to indirect effects. To overcome this problem, we propose the c-CSN method, which can construct the conditional cell-specific network (CCSN) for each cell. c-CSN method can measure the direct associations between genes by eliminating the indirect associations. c-CSN can be used for cell clustering and dimension reduction on a network basis of single cells. Intuitively, each CCSN can be viewed as the transformation from less "reliable" gene expression to more "reliable" gene-gene associations in a cell. Based on CCSN, we further design network flow entropy (NFE) to estimate the differentiation potency of a single cell. A number of scRNA-seq datasets were used to demonstrate the advantages of our approach. 1) One direct association network is generated for one cell. 2) Most existing scRNA-seq methods designed for gene expression matrices are also applicable to c-CSN-transformed degree matrices. 3) CCSN-based NFE helps resolving the direction of differentiation trajectories by quantifying the potency of each cell. c-CSN is publicly available at https://github.com/LinLi-0909/c-CSN.