Linking chromatin acylation mark-defined proteome and genome in living cells.

A generalizable strategy with programmable site specificity for in situ profiling of histone modifications on unperturbed chromatin remains highly desirable but challenging. We herein developed a single-site-resolved multi-omics (SiTomics) strategy for systematic mapping of dynamic modifications and subsequent profiling of chromatinized proteome and genome defined by specific chromatin acylations in living cells. By leveraging the genetic code expansion strategy, our SiTomics toolkit revealed distinct crotonylation (e.g., H3K56cr) and ß-hydroxybutyrylation (e.g., H3K56bhb) upon short chain fatty acids stimulation and established linkages for chromatin acylation mark-defined proteome, genome, and functions. This led to the identification of GLYR1 as a distinct interacting protein in modulating H3K56cr's gene body localization as well as the discovery of an elevated super-enhancer repertoire underlying bhb-mediated chromatin modulations. SiTomics offers a platform technology for elucidating the "metabolites-modification-regulation" axis, which is widely applicable for multi-omics profiling and functional dissection of modifications beyond acylations and proteins beyond histones.

Understanding the Function of Mammalian Sirtuins and Protein Lysine Acylation.

Protein lysine acetylation is an important posttranslational modification that regulates numerous biological processes. Targeting lysine acetylation regulatory factors, such as acetyltransferases, deacetylases, and acetyl-lysine recognition domains, has been shown to have potential for treating human diseases, including cancer and neurological diseases. Over the past decade, many other acyl-lysine modifications, such as succinylation, crotonylation, and long-chain fatty acylation, have also been investigated and shown to have interesting biological functions. Here, we provide an overview of the functions of different acyl-lysine modifications in mammals. We focus on lysine acetylation as it is well characterized, and principles learned from acetylation are useful for understanding the functions of other lysine acylations. We pay special attention to the sirtuins, given that the study of sirtuins has provided a great deal of information about the functions of lysine acylation. We emphasize the regulation of sirtuins to illustrate that their regulation enables cells to respond to various signals and stresses.

Histone H3 lysine 27 crotonylation mediates gene transcriptional repression in chromatin.

Histone lysine acylation, including acetylation and crotonylation, plays a pivotal role in gene transcription in health and diseases. However, our understanding of histone lysine acylation has been limited to gene transcriptional activation. Here, we report that histone H3 lysine 27 crotonylation (H3K27cr) directs gene transcriptional repression rather than activation. Specifically, H3K27cr in chromatin is selectively recognized by the YEATS domain of GAS41 in complex with SIN3A-HDAC1 co-repressors. Proto-oncogenic transcription factor MYC recruits GAS41/SIN3A-HDAC1 complex to repress genes in chromatin, including cell-cycle inhibitor p21. GAS41 knockout or H3K27cr-binding depletion results in p21 de-repression, cell-cycle arrest, and tumor growth inhibition in mice, explaining a causal relationship between GAS41 and MYC gene amplification and p21 downregulation in colorectal cancer. Our study suggests that H3K27 crotonylation signifies a previously unrecognized, distinct chromatin state for gene transcriptional repression in contrast to H3K27 trimethylation for transcriptional silencing and H3K27 acetylation for transcriptional activation.

CDYL1-dependent decrease in lysine crotonylation at DNA double-strand break sites functionally uncouples transcriptional silencing and repair.

Previously, we showed that CDYL1 is recruited to DNA double-strand breaks (DSBs) to promote homologous recombination (HR) repair and foster transcriptional silencing. However, how CDYL1 elicits DSB-induced silencing is not fully understood. Here, we identify a CDYL1-dependent local decrease in the transcriptionally active marks histone lysine crotonylation (Kcr) and crotonylated lysine 9 of H3 (H3K9cr) at AsiSI-induced DSBs, which correlates with transcriptional silencing. Mechanistically, we reveal that CDYL1 crotonyl-CoA hydratase activity counteracts Kcr and H3K9cr at DSB sites, which triggers the eviction of the transcription elongation factor ENL and fosters transcriptional silencing. Furthermore, genetic inhibition of CDYL1 hydratase activity blocks the reduction in H3K9cr and alleviates DSB-induced silencing, whereas HR efficiency unexpectedly remains intact. Therefore, our results functionally uncouple the repair and silencing activity of CDYL1 at DSBs. In a broader context, we address a long-standing question concerning the functional relationship between HR repair and DSB-induced silencing, suggesting that they may occur independently.

Decrotonylation of cGAS K254 prompts homologous recombination repair by blocking its DNA binding and releasing PARP1.

Cyclic GMP-AMP synthase (cGAS), a cytosolic DNA sensor, also exhibits nuclear genomic localization and is involved in DNA damage signaling. In this study, we investigated the impact of cGAS crotonylation on the regulation of the DNA damage response, particularly homologous recombination repair, following exposure to ionizing radiation (IR). Lysine 254 of cGAS is constitutively crotonylated by the CREB-binding protein; however, IR-induced DNA damage triggers sirtuin 3 (SIRT3)-mediated decrotonylation. Lysine 254 decrotonylation decreased the DNA-binding affinity of cGAS and inhibited its interaction with PARP1, promoting homologous recombination repair. Moreover, SIRT3 suppression led to homologous recombination repair inhibition and markedly sensitized cancer cells to IR and DNA-damaging chemicals, highlighting SIRT3 as a potential target for cancer therapy. Overall, this study revealed the crucial role of cGAS crotonylation in the DNA damage response. Furthermore, we propose that modulating cGAS and SIRT3 activities could be potential strategies for cancer therapy.

Dynamic profiling and functional interpretation of histone lysine crotonylation and lactylation during neural development.

Metabolites such as crotonyl-CoA and lactyl-CoA influence gene expression by covalently modifying histones, known as histone lysine crotonylation (Kcr) and lysine lactylation (Kla). However, the existence patterns, dynamic changes, biological functions and associations of these modifications with histone lysine acetylation and gene expression during mammalian development remain largely unknown. Here, we find that histone Kcr and Kla are widely distributed in the brain and undergo global changes during neural development. By profiling the genome-wide dynamics of H3K9ac, H3K9cr and H3K18la in combination with ATAC and RNA sequencing, we reveal that these marks are tightly correlated with chromatin state and gene expression, and extensively involved in transcriptome remodeling to promote cell-fate transitions in the developing telencephalon. Importantly, we demonstrate that global Kcr and Kla levels are not the consequence of transcription and identify the histone deacetylases (HDACs) 1-3 as novel 'erasers' of H3K18la. Using P19 cells as an induced neural differentiation system, we find that HDAC1-3 inhibition by MS-275 pre-activates neuronal transcriptional programs by stimulating multiple histone lysine acylations simultaneously. These findings suggest that histone Kcr and Kla play crucial roles in the epigenetic regulation of neural development.

Auto-Kla: a novel web server to discriminate lysine lactylation sites using automated machine learning.

Recently, lysine lactylation (Kla), a novel post-translational modification (PTM), which can be stimulated by lactate, has been found to regulate gene expression and life activities. Therefore, it is imperative to accurately identify Kla sites. Currently, mass spectrometry is the fundamental method for identifying PTM sites. However, it is expensive and time-consuming to achieve this through experiments alone. Herein, we proposed a novel computational model, Auto-Kla, to quickly and accurately predict Kla sites in gastric cancer cells based on automated machine learning (AutoML). With stable and reliable performance, our model outperforms the recently published model in the 10-fold cross-validation. To investigate the generalizability and transferability of our approach, we evaluated the performance of our models trained on two other widely studied types of PTM, including phosphorylation sites in host cells infected with SARS-CoV-2 and lysine crotonylation sites in HeLa cells. The results show that our models achieve comparable or better performance than current outstanding models. We believe that this method will become a useful analytical tool for PTM prediction and provide a reference for the future development of related models. The web server and source code are available at http://tubic.org/Kla and https://github.com/tubic/Auto-Kla, respectively.

RMTLysPTM: recognizing multiple types of lysine PTM sites by deep analysis on sequences.

Post-translational modification (PTM) occurs after a protein is translated from ribonucleic acid. It is an important living creature life phenomenon because it is implicated in almost all cellular processes. Identification of PTM sites from a given protein sequence is a hot topic in bioinformatics. Lots of computational methods have been proposed, and they provide good performance. However, most previous methods can only tackle one PTM type. Few methods consider multiple PTM types. In this study, a multi-label classification model, named RMTLysPTM, was developed to recognize four types of lysine (K) PTM sites, including acetylation, crotonylation, methylation and succinylation. The surrounding sites of a lysine site were selected to constitute a peptide segment, representing the lysine at the center. Deep analysis was conducted to count the distribution of 2-residues with fixed location across the four types of lysine PTM sites. By aggregating the distribution information of 2-residues in one peptide segment, the peptide segment was encoded by informative features. Furthermore, a prediction engine that can precisely capture the traits of the above representations was designed to recognize the types of lysine PTM sites. The cross-validation results on two datasets (Qiu and CPLM training datasets) suggested that the model had extremely high performance and RMTLysPTM had strong generalization ability by testing it on protein Q16778 and CPLM testing datasets. The model was found to be generally superior to all previous models and those using popular methods and features. A web server was set up for RMTLysPTM, and it can be accessed at http://119.3.127.138/.

Protein Crotonylation Promotes Osteogenic Differentiation of Periodontal Ligament Stem Cells via the PI3K-AKT Pathway.

Posttranslational modifications (PTMs) are crucial regulatory mechanisms for cellular differentiation and organismal development. Acylation modification is one of the main PTMs that plays a pivotal role in regulating the osteogenic differentiation of mesenchymal stem cells and is a focal point of research in bone tissue regeneration. However, its mechanism remains incompletely understood. This article aims to investigate the impact of protein crotonylation on osteogenic differentiation in periodontal ligament stem cells (PDLSCs) and elucidate its underlying mechanisms. Western blot analysis identified that the modification level of acetylation, crotonylation, and succinylation were significantly upregulated after osteogenic induction of PDLSCs. Subsequently, sodium crotonate (NaCr) was added to the medium and acyl-CoA synthetase short-chain family member 2 (ACSS2) was knocked down by short hairpin RNA plasmids to regulate the total level of protein crotonylation. The results indicated that treatment with NaCr promoted the expression of osteogenic differentiation-related factors in PDLSCs, whereas silencing ACSS2 had the opposite effect. In addition, mass spectrometry analysis was used to investigate the comprehensive analysis of proteome-wide crotonylation in PDLSCs under osteogenic differentiation. The analysis revealed that the level of protein crotonylation related to the PI3K-AKT signaling pathway was significantly upregulated in PDLSCs after osteogenic induction. Treatment with NaCr and silencing ACSS2 affected the activation of the PI3K-AKT signaling pathway. Collectively, our study demonstrates that protein crotonylation promotes osteogenic differentiation of PDLSCs via the PI3K-AKT pathway, providing a novel targeting therapeutic approach for bone tissue regeneration.

Chromodomain Protein CDYL Acts as a Crotonyl-CoA Hydratase to Regulate Histone Crotonylation and Spermatogenesis.

Lysine crotonylation (Kcr) is a newly identified histone modification that is associated with active transcription in mammalian cells. Here we report that the chromodomain Y-like transcription corepressor CDYL negatively regulates histone Kcr by acting as a crotonyl-CoA hydratase to convert crotonyl-CoA to ß-hydroxybutyryl-CoA. We showed that the negative regulation of histone Kcr by CDYL is intrinsically linked to its transcription repression activity and functionally implemented in the reactivation of sex chromosome-linked genes in round spermatids and genome-wide histone replacement in elongating spermatids. Significantly, Cdyl transgenic mice manifest dysregulation of histone Kcr and reduction of male fertility with a decreased epididymal sperm count and sperm cell motility. Our study uncovers a biochemical pathway in the regulation of histone Kcr and implicates CDYL-regulated histone Kcr in spermatogenesis, adding to the understanding of the physiology of male reproduction and the mechanism of the spermatogenic failure in AZFc (Azoospermia Factor c)-deleted infertile men.

Protein lysine four-carbon acylations in health and disease.

Lysine acylation, a type of posttranslational protein modification sensitive to cellular metabolic states, influences the functions of target proteins involved in diverse cellular processes. Particularly, lysine butyrylation, crotonylation, ß-hydroxybutyrylation, and 2-hydroxyisobutyrylation, four types of four-carbon acylations, are modulated by intracellular concentrations of their respective acyl-CoAs and sensitive to alterations of nutrient metabolism induced by cellular and/or environmental signals. In this review, we discussed the metabolic pathways producing these four-carbon acyl-CoAs, the regulation of lysine acylation and deacylation, and the functions of individual lysine acylation.

Glutaryl-CoA dehydrogenase suppresses tumor progression and shapes an anti-tumor microenvironment in hepatocellular carcinoma.

BACKGROUND & AIMS: Crotonylation, a crotonyl-CoA-based non-enzymatic protein translational modification, affects diverse biological processes, such as spermatogenesis, tissue injury, inflammation, and neuropsychiatric diseases. Crotonylation is decreased in hepatocellular carcinomas (HCCs), but the mechanism remains unknown. In this study, we aim to describe the role of glutaryl-CoA dehydrogenase (GCDH) in tumor suppression. METHODS: Three cohorts containing 40, 248 and 17 pairs of samples were used to evaluate the link between GCDH expression levels and clinical characteristics of HCC, as well as responses to anti-programmed cell death protein 1 (PD-1) treatment. Subcutaneous xenograft, orthotopic xenograft, Trp53Δhep/Δhep; MYC- and Ctnnboe;METoe-driven mouse models were adopted to validate the effects of GCDH on HCC suppression. RESULTS: GCDH depletion promoted HCC growth and metastasis, whereas its overexpression reversed these processes. As GCDH converts glutaryl-CoA to crotonyl-CoA to increase crotonylation levels, we performed lysine crotonylome analysis and identified the pentose phosphate pathway (PPP) and glycolysis-related proteins PGD, TKT, and ALDOC as GCDH-induced crotonylation targets. Crotonyl-bound targets showed allosteric effects that controlled their enzymatic activities, leading to decreases in ribose 5-phosphate and lactate production, further limiting the Warburg effect. PPP blockade also stimulated peroxidation, synergizing with senescent modulators to induce senescence in GCDHhigh cells. These cells induced the infiltration of immune cells by the SASP (senescence-associated secretory cell phenotype) to shape an anti-tumor immune microenvironment. Meanwhile, the GCDHlow population was sensitized to anti-PD-1 therapy. CONCLUSION: GCDH inhibits HCC progression via crotonylation-induced suppression of the PPP and glycolysis, resulting in HCC cell senescence. The senescent cell further shapes an anti-tumor microenvironment via the SASP. The GCDHlow population is responsive to anti-PD-1 therapy because of the increased presence of PD-1+CD8+ T cells. IMPACT AND IMPLICATIONS: Glutaryl-CoA dehydrogenase (GCDH) is a favorable prognostic indicator in liver, lung, and renal cancers. In addition, most GCDH depletion-induced toxic metabolites originate from the liver, accumulate locally, and cannot cross the blood-brain barrier. Herein, we show that GCDH inhibits hepatocellular carcinoma (HCC) progression via crotonylation-induced suppression of the pentose phosphate pathway and glycolysis, resulting in HCC cell senescence. We also found that more PD-1+CD8+ T cells are present in the GCDHlow population, who are thus more responsive to anti-PD-1 therapy. Given that the GCDHlow and GCDHhigh HCC population can be distinguished based on serum glucose and ammonia levels, it will be worthwhile to evaluate the curative effects of pro-senescent and immune-therapeutic strategies based on the expression levels of GCDH.

DeepCap-Kcr: accurate identification and investigation of protein lysine crotonylation sites based on capsule network.

Lysine crotonylation (Kcr) is a posttranslational modification widely detected in histone and nonhistone proteins. It plays a vital role in human disease progression and various cellular processes, including cell cycle, cell organization, chromatin remodeling and a key mechanism to increase proteomic diversity. Thus, accurate information on such sites is beneficial for both drug development and basic research. Existing computational methods can be improved to more effectively identify Kcr sites in proteins. In this study, we proposed a deep learning model, DeepCap-Kcr, a capsule network (CapsNet) based on a convolutional neural network (CNN) and long short-term memory (LSTM) for robust prediction of Kcr sites on histone and nonhistone proteins (mammals). The proposed model outperformed the existing CNN architecture Deep-Kcr and other well-established tools in most cases and provided promising outcomes for practical use; in particular, the proposed model characterized the internal hierarchical representation as well as the important features from multiple levels of abstraction automatically learned from a small number of samples. The trained model was well generalized in other species (papaya). Moreover, we showed the features and properties generated by the internal capsule layer that can explore the internal data distribution related to biological significance (as a motif detector). The source code and data are freely available at https://github.com/Jhabindra-bioinfo/DeepCap-Kcr.

Adapt-Kcr: a novel deep learning framework for accurate prediction of lysine crotonylation sites based on learning embedding features and attention architecture.

Protein lysine crotonylation (Kcr) is an important type of posttranslational modification that is associated with a wide range of biological processes. The identification of Kcr sites is critical to better understanding their functional mechanisms. However, the existing experimental techniques for detecting Kcr sites are cost-ineffective, to a great need for new computational methods to address this problem. We here describe Adapt-Kcr, an advanced deep learning model that utilizes adaptive embedding and is based on a convolutional neural network together with a bidirectional long short-term memory network and attention architecture. On the independent testing set, Adapt-Kcr outperformed the current state-of-the-art Kcr prediction model, with an improvement of 3.2% in accuracy and 1.9% in the area under the receiver operating characteristic curve. Compared to other Kcr models, Adapt-Kcr additionally had a more robust ability to distinguish between crotonylation and other lysine modifications. Another model (Adapt-ST) was trained to predict phosphorylation sites in SARS-CoV-2, and outperformed the equivalent state-of-the-art phosphorylation site prediction model. These results indicate that self-adaptive embedding features perform better than handcrafted features in capturing discriminative information; when used in attention architecture, this could be an effective way of identifying protein Kcr sites. Together, our Adapt framework (including learning embedding features and attention architecture) has a strong potential for prediction of other protein posttranslational modification sites.

The Plant PTM Viewer 2.0: in-depth exploration of plant protein modification landscapes.

Post-translational modifications (PTMs) greatly increase protein diversity and functionality. To help the plant research community interpret the ever-increasing number of reported PTMs, the Plant PTM Viewer (https://www.psb.ugent.be/PlantPTMViewer) provides an intuitive overview of plant protein PTMs and the tools to assess it. This update includes 62 novel PTM profiling studies, adding a total of 112 000 modified peptides reporting plant PTMs, including 14 additional PTM types and three species (moss, tomato, and soybean). Furthermore, an open modification re-analysis of a large-scale Arabidopsis thaliana mass spectrometry tissue atlas identified previously uncharted landscapes of lysine acylations predominant in seed and flower tissues and 3-phosphoglycerylation on glycolytic enzymes in plants. An extra 'Protein list analysis' tool was developed for retrieval and assessing the enrichment of PTMs in a protein list of interest. We conducted a protein list analysis on nuclear proteins, revealing a substantial number of redox modifications in the nucleus, confirming previous assumptions regarding the redox regulation of transcription. We encourage the plant research community to use PTM Viewer 2.0 for hypothesis testing and new target discovery, and also to submit new data to expand the coverage of conditions, plant species, and PTM types, thereby enriching our understanding of plant biology.

Marine algae-derived oligosaccharide via protein crotonylation of key targeting for management of type 2 diabetes mellitus in the elderly.

Global aging is a tendency of the world, as is the increasing prevalence of diabetes, and the two are closely linked. In our early research, Enteromorpha prolifera oligosaccharide (EPO) possesses the excellent ability of anti-oxidative, anti-inflammatory, and anti-diabetic. We aim to further explore the deeper mechanism of how EPO delays aging and regulates glycometabolism. EPO effectively impacts crotonylation procession to enhance glucose metabolism and reduce cell senescence in aging diabetic rats. Crotonylation modification of XPO1 influences the expression of critical genes, including p53, CDK1, and CCNB1, which affect cell cycle regulation and aging. Additionally, EPO improves glucose metabolism by inhibiting the crotonylation modification of HSPA8-K126 and activating the AKT pathway. EPO promotes crotonylation of histones in intestinal cells, influencing the aging process by increasing the butyric acid-producing bacteria Ruminococcaceae. The observed enhancement in pyrimidine metabolism underscores EPO's potential role in regulating intestinal health, presenting a promising avenue for delaying aging. In summary, our findings affirm EPO as a naturally bioactive ingredient with significant potential for anti-aging and antidiabetic interventions.

Deep-Kcr: accurate detection of lysine crotonylation sites using deep learning method.

As a newly discovered protein posttranslational modification, histone lysine crotonylation (Kcr) involved in cellular regulation and human diseases. Various proteomics technologies have been developed to detect Kcr sites. However, experimental approaches for identifying Kcr sites are often time-consuming and labor-intensive, which is difficult to widely popularize in large-scale species. Computational approaches are cost-effective and can be used in a high-throughput manner to generate relatively precise identification. In this study, we develop a deep learning-based method termed as Deep-Kcr for Kcr sites prediction by combining sequence-based features, physicochemical property-based features and numerical space-derived information with information gain feature selection. We investigate the performances of convolutional neural network (CNN) and five commonly used classifiers (long short-term memory network, random forest, LogitBoost, naive Bayes and logistic regression) using 10-fold cross-validation and independent set test. Results show that CNN could always display the best performance with high computational efficiency on large dataset. We also compare the Deep-Kcr with other existing tools to demonstrate the excellent predictive power and robustness of our method. Based on the proposed model, a webserver called Deep-Kcr was established and is freely accessible at http://lin-group.cn/server/Deep-Kcr.

nhKcr: a new bioinformatics tool for predicting crotonylation sites on human nonhistone proteins based on deep learning.

Lysine crotonylation (Kcr) is a newly discovered type of protein post-translational modification and has been reported to be involved in various pathophysiological processes. High-resolution mass spectrometry is the primary approach for identification of Kcr sites. However, experimental approaches for identifying Kcr sites are often time-consuming and expensive when compared with computational approaches. To date, several predictors for Kcr site prediction have been developed, most of which are capable of predicting crotonylation sites on either histones alone or mixed histone and nonhistone proteins together. These methods exhibit high diversity in their algorithms, encoding schemes, feature selection techniques and performance assessment strategies. However, none of them were designed for predicting Kcr sites on nonhistone proteins. Therefore, it is desirable to develop an effective predictor for identifying Kcr sites from the large amount of nonhistone sequence data. For this purpose, we first provide a comprehensive review on six methods for predicting crotonylation sites. Second, we develop a novel deep learning-based computational framework termed as CNNrgb for Kcr site prediction on nonhistone proteins by integrating different types of features. We benchmark its performance against multiple commonly used machine learning classifiers (including random forest, logitboost, naïve Bayes and logistic regression) by performing both 10-fold cross-validation and independent test. The results show that the proposed CNNrgb framework achieves the best performance with high computational efficiency on large datasets. Moreover, to facilitate users' efforts to investigate Kcr sites on human nonhistone proteins, we implement an online server called nhKcr and compare it with other existing tools to illustrate the utility and robustness of our method. The nhKcr web server and all the datasets utilized in this study are freely accessible at http://nhKcr.erc.monash.edu/.

Development of Lysine Crotonyl-Mimic Probe to Covalently Identify H3K27Cr Interacting Proteins.

Histone lysine crotonylation (Kcr) is one newly discovered acylation modification and regulates numerous pathophysiological processes. The binding affinity between Kcr and its interacting proteins is generally weak, which makes it difficult to effectively identify Kcr-interacting partners. Changing the amide of crotonyl to an ester increased reactivity with proximal cysteines and retained specificity for Kcr antibody. The probe "H3g27Cr" was designed by incorporating the ester functionality into a H3K27 peptide. Using this probe, multiple Kcr-interacting partners including STAT3 were successfully identified, and this has not been reported previously. Further experiments suggested that STAT3 possibly could form complexes with Histone deacetylase HDACs to downregulate the acetylation and crotonylation of Histone H3K27. Our unique design provided intriguing tools to further explore Kcr-interacting proteins and elucidate their working mechanisms.

Systematic proteomics profiling of lysine crotonylation of the lung at Pseudoglandular and Canalicular phases in human fetus.

Lung tissue is an important organ of the fetus, and genomic research on its development has improved our understanding of the biology of this tissue. However, the proteomic research of developing fetal lung tissue is still very scarce. We conducted comprehensive analysis of two developmental stages of fetal lung tissue of proteomics. It showed the developmental characteristics of lung tissue, such as the down-regulation of metabolism-related protein expression, the up-regulation of cell cycle-related proteins, and the regulation in proteins and pathways related to lung development. In addition, we also discovered some key core proteins related to lung development, and provided some key crotonylation modification sites that regulation during lung tissue development. Our comprehensive analysis of lung proteomics can provide a more comprehensive understanding of the developmental status of lung tissue, and provide a certain reference for future research and epigenetics of lung tissue.