CRISPRi-mediated functional analysis of NKX2-1-binding sites in the lung

Identification of oxidative stress signatures of lung adenocarcinoma and prediction of patient prognosis or treatment response with single-cell RNA sequencing and bulk RNA sequencing data

Lung adenocarcinoma (LUAD), the most common subtype of lung cancer globally, has seen improved prognosis with advancements in diagnostic, surgical, radiotherapy, and molecular therapy techniques, while its 5-year survival rate remains low. Molecular biomarkers provide prognostic value. Oxidative stress factors, such as reactive nitrogen species and ROS, are crucial in various stages of tumor progression, influencing cell transformation, proliferation, angiogenesis, and metastasis. ROS demonstrate dual roles, affecting tumor cells, hypoxia sensitivity, and the microenvironment. Comprehensive analysis of oxidative stress in LUAD has not been conducted to date. Therefore, we systematically investigated the regulatory patterns of oxidative stress in LUAD based on oxidative stress-related genes and correlated these patterns with cellular infiltration characteristics of the tumor immune microenvironment. The model utilizes single-factor Cox analysis to screen key differential genes with prognostic value and employs least absolute shrinkage and selection operator (LASSO) penalized Cox regression analysis to construct a prognostic-related prediction model. Ten candidate genes were selected based on this model. The risk score was constructed using the coefficients and expression levels of these ten genes. Furthermore, the impact of this risk score on overall survival (OS) was determined. Two genes with the most significant differential expression, SFTPB and S100P, were selected through qRT-PCR. Cell experiments including CCK-8, Edu, transwell assays confirmed their effects on lung cancer cells growth, consistent with the results of bioinformatics analysis. These findings suggested that this model held potential clinical value for evaluating the prognosis of lung adenocarcinoma.

Dosage amplification dictates oncogenic regulation by the NKX2-1 lineage factor in lung adenocarcinoma

Amplified oncogene expression is a critical and widespread driver event in cancer, yet our understanding of how amplification-mediated elevated dosage mediates oncogenic regulation is limited. Here, we find that the most significant focal amplification event in lung adenocarcinoma (LUAD) targets a lineage super-enhancer near the NKX2-1 lineage transcription factor. The NKX2-1 super-enhancer is targeted by focal and co-amplification with NKX2-1, and activation or repression controls NKX2-1 expression. We find that NKX2-1 is a widespread dependency in LUAD cell lines, where NKX2-1 pioneers enhancer accessibility to drive a lineage addicted state in LUAD, and NKX2-1 confers persistence to EGFR inhibitors. Notably, we find that oncogenic NKX2-1 regulation requires expression above a minimum dosage threshold-NKX2-1 dosage below this threshold is insufficient for cell viability, enhancer remodeling, and TKI persistence. Our data suggest that copy-number amplification can be a gain-of-function alteration, wherein amplification elevates oncogene expression above a critical dosage required for oncogenic regulation and cancer cell survival.

FOXA2 Cooperates with Mutant KRAS to Drive Invasive Mucinous Adenocarcinoma of the Lung

The endoderm-lineage transcription factor FOXA2 has been shown to inhibit lung tumorigenesis in in vitro and xenograft studies using lung cancer cell lines. However, FOXA2 expression in primary lung tumors does not correlate with an improved patient survival rate, and the functional role of FOXA2 in primary lung tumors remains elusive. To understand the role of FOXA2 in primary lung tumors in vivo, here, we conditionally induced the expression of FOXA2 along with either of the two major lung cancer oncogenes, EGFRL858R or KRASG12D, in the lung epithelium of transgenic mice. Notably, FOXA2 suppressed autochthonous lung tumor development driven by EGFRL858R, whereas FOXA2 promoted tumor growth driven by KRASG12D. Importantly, FOXA2 expression along with KRASG12D produced invasive mucinous adenocarcinoma (IMA) of the lung, a fatal mucus-producing lung cancer comprising approximately 5% of human lung cancer cases. In the mouse model in vivo and human lung cancer cells in vitro, FOXA2 activated a gene regulatory network involved in the key mucous transcription factor SPDEF and upregulated MUC5AC, whose expression is critical for inducing IMA. Coexpression of FOXA2 with mutant KRAS synergistically induced MUC5AC expression compared with that induced by FOXA2 alone. ChIP-seq combined with CRISPR interference indicated that FOXA2 bound directly to the enhancer region of MUC5AC and induced the H3K27ac enhancer mark. Furthermore, FOXA2 was found to be highly expressed in primary tumors of human IMA. Collectively, this study reveals that FOXA2 is not only a biomarker but also a driver for IMA in the presence of a KRAS mutation.

SIGNIFICANCE

FOXA2 expression combined with mutant KRAS drives invasive mucinous adenocarcinoma of the lung by synergistically promoting a mucous transcriptional program, suggesting strategies for targeting this lung cancer type that lacks effective therapies.

Expression landscapes in non-small cell lung cancer shaped by the thyroid transcription factor 1

TTF-1-expressing non-small cell lung cancer (NSCLC) is one of the most prevalent lung cancer types worldwide. However, theparadoxical activity of TTF-1 in both lung carcinogenesis and tumor suppression is believed to be context-dependentwhich calls for a deeper understanding about the pathological expression of TTF-1. In addition, the expression circuitry of TTF-1-target genes in NSCLC has not been well examined which necessitates to revisit the involvement of TTF-1- in a multitude of oncologic pathways. We used RNA-seq and clinical data of patients from The Cancer Genome Atlas (TCGA) and Genotype-Tissue Expression (GTEx), including ChIP-seq data from different NSCLC cell lines, and mapped the proteome of NSCLC tumor. Our analysis showed significant variability in TTF-1 expression among NSCLC,and further clarified that this variability is orchestrated at the transcriptional level. We also found that high TTF-1 expression could negatively influence the survival outcomes of stage 1 LUAD which may be attributed to growth factor receptor-driven activation of mitogenic and angiogenic pathways. Mechanistically, TTF-1 may also control the genes associated with pathways involved in acquired TKI drug resistance or response to immune checkpoint inhibitors. Lastly, proteome-based biomarker discovery in stage 1 LUAD showed that TTF-1 positivity is potentially associated with the upregulation of several oncogenes which includes interferon proteins, MUC1, STAT3, and EIF2AK2. Collectively, this study highlights the potential involvement of TTF-1 in cell proliferation, immune evasion, and angiogenesis in early-stage NSCLC.

Development of a Lung Vacancy Mouse Model through CRISPR/Cas9-Mediated Deletion of Thyroid Transcription Factor 1 Exon 2

A developmental niche vacancy in host embryos is necessary for stem cell complementation-based organ regeneration (SCOG). Thyroid transcription factor 1 (TTF-1) is a tissue-specific transcription factor that regulates the embryonic development and differentiation of the thyroid and, more importantly, lungs; thus, it has been considered as a master gene to knockout in order to develop a lung vacancy host. TTF-1 knockout mice were originally produced by inserting a stop codon in Exon 3 of the gene (E3stop) through embryonic stem cell-based homologous recombination. The main problems of utilizing E3stop host embryos for lung SCOG are that these animals all have a tracheoesophageal fistula (TEF), which cannot be corrected by donor stem cells, and most of them have monolateral sac-like lungs. To improve the mouse model towards achieving SCOG-based lung generation, in this project, we used the CRISPR/Cas9 tool to remove Exon 2 of the gene by zygote microinjection and successfully produced TTF-1 knockout (E2del) mice. Similar to E3stop, E2del mice are birth-lethal due to retarded lung development with sac-like lungs and only a rudimentary bronchial tree, increased basal cells but without alveolar type II cells and blood vessels, and abnormal thyroid development. Unlike E3stop, 57% of the E2del embryos presented type I tracheal agenesis (TA, a kind of human congenital malformation) with a shortened trachea and clear separations of the trachea and esophagus, while the remaining 43% had TEF. Furthermore, all the E2del mice had bilateral sac-like lungs. Both TA and bilateral sac-like lungs are preferred in SCOG. Our work presents a new strategy for producing SCOG host embryos that may be useful for lung regeneration.

A Novel CRISPR Interference Effector Enabling Functional Gene Characterization with Synthetic Guide RNAs

While CRISPR interference (CRISPRi) systems have been widely implemented in pooled lentiviral screening, there has been limited use with synthetic guide RNAs for the complex phenotypic readouts enabled by experiments in arrayed format. Here we describe a novel deactivated Cas9 fusion protein, dCas9-SALL1-SDS3, which produces greater target gene repression than first or second generation CRISPRi systems when used with chemically modified synthetic single guide RNAs (sgRNAs), while exhibiting high target specificity. We show that dCas9-SALL1-SDS3 interacts with key members of the histone deacetylase and Swi-independent three complexes, which are the endogenous functional effectors of SALL1 and SDS3. Synthetic sgRNAs can also be used with in vitro-transcribed dCas9-SALL1-SDS3 mRNA for short-term delivery into primary cells, including human induced pluripotent stem cells and primary T cells. Finally, we used dCas9-SALL1-SDS3 for functional gene characterization of DNA damage host factors, orthogonally to small interfering RNA, demonstrating the ability of the system to be used in arrayed-format screening.

Regulation of mi-RNAs Target Cancer Genes Between Exercise and Non-exercise in Rat Rheumatoid Arthritis Induction: Pilot Study

Introduction:

Rheumatoid arthritis is associated with various cancers. Many studies have investigated physical exercise interventions as health improvements to ameliorate the risk of cancer during rheumatoid arthritis diagnosis. Recently, microRNAs were used as biomarkers for health assessment and cancer prediction in rheumatoid arthritis patients.

Methods:

The effects of exercise interventions on serum microRNAs were investigated in pristane-induced arthritis (PIA) rat models. Twelve Sprague-Dawley male rats were divided into 4 groups including non-exercise without PIA (N-EX), non-exercise with PIA (N-EX + PIA), exercise without PIA (EX) and exercise with PIA (EX + PIA). Blood samples were collected at the end of the study period to analyze miRNA biomarkers and target cancer gene predictions.

Results:

Four significant Rattus norvegicus (rno-microRNAs) may purpose as tumor suppressors were identified as potential target cancer gene candidate expressions within the 4 comparative interventional exercise groups. One rno-microRNA and target cancer gene candidate was up-regulated and 3 rno-microRNAs and their target cancer genes were down-regulated.

Conclusions:

Exercise interventions affected rno-miRNAs regulated target cancer gene candidates ITPR3, SOCS6, ITGA6, and NKX2-1 as biomarkers for cancer prognosis in rheumatoid arthritis diagnosis.

CRISPRi links COVID-19 GWAS loci to LZTFL1 and RAVER1

Background

To identify host genetic variants (SNPs) associated with COVID-19 disease severity, a number of genome-wide association studies (GWAS) have been conducted. Since most of the identified variants are located at non-coding regions, such variants are presumed to affect the expression of neighbouring genes, thereby influencing COVID-19 disease severity. However, it remains largely unknown which genes are influenced by such COVID-19 GWAS loci.

Methods

CRISPRi (interference)-mediated gene expression analysis was performed to identify genes functionally regulated by COVID-19 GWAS loci by targeting regions near the loci (SNPs) in lung epithelial cell lines. The expression of CRISPRi-identified genes was investigated using COVID-19-contracted human and monkey lung single-nucleus/cell (sn/sc) RNA-seq datasets.

Findings

CRISPRi analysis indicated that a region near rs11385942 at chromosome 3p21.31 (locus of highest significance with COVID-19 disease severity at intron 5 of LZTFL1) significantly affected the expression of LZTFL1 (P<0.05), an airway cilia regulator. A region near rs74956615 at chromosome 19p13.2 (locus located at the 3’ untranslated exonic region of RAVER1), which is associated with critical illness in COVID-19, affected the expression of RAVER1 (P<0.05), a coactivator of MDA5 (IFIH1), which induces antiviral response genes, including ICAM1. The sn/scRNA-seq datasets indicated that the MDA5/RAVER1-ICAM1 pathway was activated in lung epithelial cells of COVID-19-resistant monkeys but not those of COVID-19-succumbed humans.

Interpretation

Patients with risk alleles of rs11385942 and rs74956615 may be susceptible to critical illness in COVID-19 in part through weakened airway viral clearance via LZTFL1-mediated ciliogenesis and diminished antiviral immune response via the MDA5/RAVER1 pathway, respectively.

Funding

NIH.

TGF-β signaling and Creb5 cooperatively regulate Fgf18 to control pharyngeal muscle development

The communication between myogenic cells and their surrounding connective tissues is indispensable for muscle morphogenesis. During late embryonic development in mice, myogenic progenitors migrate to discrete sites to form individual muscles. The detailed mechanism of this process remains unclear. Using mouse levator veli palatini (LVP) development as a model, we systematically investigated how a distinct connective tissue subpopulation, perimysial fibroblasts, communicates with myogenic cells to regulate mouse pharyngeal myogenesis. Using single-cell RNAseq data analysis, we identified that TGF-β signaling is a key regulator for the perimysial fibroblasts. Loss of TGF-β signaling in the neural crest-derived palatal mesenchyme leads to defects in perimysial fibroblasts and muscle malformation in the soft palate in Osr2^Cre;Tgfbr1^fl/fl mice. In particular, Creb5, a transcription factor expressed in the perimysial fibroblasts, cooperates with TGF-β signaling to activate expression of Fgf18. Moreover, Fgf18 supports pharyngeal muscle development in vivo and exogenous Fgf18 can partially rescue myogenic cell numbers in Osr2^Cre;Tgfbr1^fl/fl samples, illustrating that TGF-β-regulated Fgf18 signaling is required for LVP development. Collectively, our findings reveal the mechanism by which TGF-β signaling achieves its functional specificity in defining the perimysial-to-myogenic signals for pharyngeal myogenesis.

Differential Regulation of Human Surfactant Protein A Genes, SFTPA1 and SFTPA2, and Their Corresponding Variants

The human SFTPA1 and SFTPA2 genes encode the surfactant protein A1 (SP-A1) and SP-A2, respectively, and they have been identified with significant genetic and epigenetic variability including sequence, deletion/insertions, and splice variants. The surfactant proteins, SP-A1 and SP-A2, and their corresponding variants play important roles in several processes of innate immunity as well in surfactant-related functions as reviewed elsewhere [1]. The levels of SP-A have been shown to differ among individuals both under baseline conditions and in response to various agents or disease states. Moreover, a number of agents have been shown to differentially regulate SFTPA1 and SFTPA2 transcripts. The focus in this review is on the differential regulation of SFTPA1 and SFTPA2 with primary focus on the role of 5′ and 3′ untranslated regions (UTRs) and flanking sequences on this differential regulation as well molecules that may mediate the differential regulation.