Skip to main content

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

Abstract

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.

Introduction

Fetal development is crucial to the survival of newborns, among which lung hypoplasia is one of the main findings of autopsy of children who have died during the neonatal period [1]. Lung development had divided into embryonic stage, fetal stage, and postpartum lung development. The fetal stage includes Pseudoglandular, Canalicular, Saccular, and Alveolarization. The stages of lung development are mainly based on morphological criteria [2]. Nevertheless, lung development is very complex, and all stages overlap. In addition, cell metabolism, growth, and differentiation are involved in the process of lung development. However, research on fetal lung development is still limited.

Crotonylation is a new type of acylation that first discovered by the research team of Professor Yingming Zhao of the University of Chicago in 2011, and selected as the annual research highlight in the field of epigenetics by Cell magazine. It refers to the introduction of crotonyl groups on the amino acid residues of histones, mainly in the lysine residues of histones, and participates in important biological processes such as gene expression regulation. Crotonylation is evolutionarily conserved, usually combined with transcriptionally active chromatin regions, and plays an important role in regulating gene expression. Crotonylation closely related to reproductive regulation, and had shown to maintain the activity of sex chromosome-related genes in meiotic metaphase.

Now, proteomics technology has been widely used in the research and exploration of biological sciences [3]. Although many analyses describing gene expression in normal lung development have been reported [4, 5], there are few protein aspects. Meanwhile, the research on lung development has used mouse models for proteomic analysis, which provides guidance for human lung development to a certain extent [6], there is no report on human fetal lung development. Secondly, protein post-translational modification (PTM) is a vital regulation of protein function, which is essential to the structure and function of protein under physiological and pathological conditions, such as gene expression and regulation, cell growth, embryonic development, metabolism, and disease treatment [7]. Among them, the crotonylation modification is genetically related, and it had reported to have important significance in cell differentiation and metabolism [8, 9]. However, it had not reported in lung development.

In this study, we combined proteomics and crotonylation modification omics to analyze the two stages of lung development (Pseudoglandular and Canalicular), revealing regulation in proteins and modifications during lung development, especially metabolism, Cell cycle, and differentiation. We provided comprehensive insights into the biological functions of proteins and crotonylation modifications in lung development, which may pave the way for the identification of the stages of lung development and the mechanisms of lung hypoplasia.

Results

Protein difference between Pseudoglandular and Canalicular in the human fetal lung

An overview of the experimental measure had shown in Fig. 1A. We used Label Free technology to identified 6,943 proteins (5,929 quantitative proteins) for protein expression to Pseudoglandular and Canalicular in LC–MS/MS. Compared with the fetal lung of 17 weeks, there were 2,645 DEPs (different expression proteins) (Fig. 1B), including 1,402 proteins up-regulated (≥ 1.5-fold) and 1243 down-regulated (≤ 0.67-fold) in the fetal lung of 12 weeks. For protein crotonylation, we identified 15,427 crotonylation sites in 3,802 proteins (6,960 quantitative crotonylation sites in 2,217 proteins). Among then, 2,196 crotonylation sites were identified in 1,156 DCPs(different crotonylation proteins) between the fetal lung of 17 weeks and 12 weeks (Fig. 1C). One thousand three hundred fifty-nine crotonylation sites were up-regulated (≥ 1.5-fold) in 793 proteins, and 837 crotonylation sites were down-regulated (≤ 0.67-fold) in 586 proteins. The hierarchical clustering analysis of DEPs had shown in Fig. 1D. We overlapped to 514 proteins in all DEPs and DCPs (Fig. 1E). We analyzed the motifs of lysine crotonylation through Motif analysis, and there was no obvious amino acid enrichment in Top5 (Fig. 1F).

Fig. 1
figure 1

Numerous Proteins Differentially Expressed in the Lung of Fetal Were Identified. A The schematic flow to study the proteome and crotonylome of fetal lung. B, C DEPs and DCPs between 12 and 17 weeks of fetal lung, respectively. D Hierarchical clustering analysis of DEPs. E Venn diagram between DEPs and DCPs. F Significantly enriched crotonylation motifs of top 5

Differentially quantified proteins to subcellular, KOG, and functional category characterization

To clarify the functions of these DEPs and DCPs, we analyzed the DEPs and DCPs using subcellular location (Fig. 2A, B). We found the top three: most DEPs are nucleus (n = 912, 34.49%), cytoplasm (n = 734, 27.76%), and mitochondria (n = 311, 11.76%); most DCPs are cytoplasm (n = 469, 40.57%), nucleus (n = 335, 28.98%), and mitochondria (n = 107, 9.26%). We also investigated the KOG categories of the DEPs and DCPs (Fig. 2C, D). We found the top three: most DEPs are [K] Transcription (n = 303), [T] Signal transduction mechanisms (n = 269), and [O] Posttranslational modification, protein turnover, chaperones (n = 217); most DCPs are [Z] Cytoskeleton (n = 334), [O] Posttranslational modification, protein turnover, chaperones (n = 234), and [J] Translation, ribosomal structure and biogenesis (n = 192).

Fig. 2
figure 2

Subcellular, KOG, and Functional Category Characterization of Differentially Quantified Proteins. A, B Subcellular location of DEPs and DCPs, respectively. C, D KOG categories of DEPs and DCPs, respectively. E, F Functional category DEPs and DCPs in GO terms.

We adopted GO analysis tool to analyze DEPs and DCPs for functional category distribution (Fig. 2E, F). Three categories of GO of DEPs displayed: cellular component (92%), molecular function (81%), biological process (91%), and 24 biological groups; DCPs: functions were established for the DEPs and further divided into three classifications: cellular component (96%), molecular function (87%), biological process (94%), and 25 biological groups. Among them, most biological process of DEPs and DCPs: the cellular process (2,185 and 1,014, respectively), the biological regulation (1,726 and 823, respectively), and the metabolic process (1584 and 694, respectively). At the same, most cellular component of DEPs and DCPs: the cell (2,401 and 1,098, respectively) and the organelle (2,032 and 972, respectively). Most molecular function of DEPs and DCPs: the binding (1,667 and859, respectively) and the catalytic activity (1000 and 376, respectively).

Functional enrichment analysis of differentially quantified proteins

For the pathways analysis of GO in DEPs, we had been enriched: cellular component (n = 8), molecular function (n = 8), and biological process (n = 14) (Fig. 3A). The enrichment of biological process shows that protein targeting has regulated significantly, molecular function shows that histone binding has regulated significantly, and the cellular components show significant regulation in ribosomes. For the pathways analysis of GO in DCPs, we had been enriched: cellular component (n = 8), molecular function (n = 8), and biological process (n = 14) (Fig. 3B). The enrichment of the biological process shows that protein targeting has regulated significantly, the molecular function shows that the chromatin binding has regulated significantly, and the cellular components show that the vesicle has regulated significantly.

Fig. 3
figure 3

Functional Enrichment Analysis of Differentially Quantified Proteins. A, B GO enrichment analysis of DEPs and DCPs, respectively. C, D Kyoto Encyclopedia of Genes and Genomes (KEGG) functional enrichment analysis of DEPs and DCPs, respectively

In KEGG pathway of DEPs, the significant ones were in Carbon metabolism, Starch and sucrose metabolism, Nitrogen metabolism, Tyrosine metabolism, DNA replication, and so on (Fig. 3C). Regarding crotonylation, the significantones were in Carbon metabolism, Starch and sucrose metabolism, Cell cycle, Spliceosome, and so on (Fig. 3D). Furthermore, we found three same pathways (Carbon metabolism, Starch and sucrose metabolism, and Ribosome) In DEPs and DCPs.

Metabolism regulation in developing fetal lung

As an important part of the living organism, metabolism is especially vital in fetal development. Cell adhesion, extracellular matrix tissue, vasculature development, and lipid metabolism are the main biological processes induced in early developmental stages. Defense/immune response, ion transport, and signal pathway had identified as the main biological processes induced in the late pregnancy [5]. We found that many metabolic pathways significantly down-regulate at the protein expression level (Fig. 4A). Such as carbon metabolism, glycolysis, TCA cycle, amino acid metabolism, and lipid metabolism, and so on. The level of crotonylation modification showed that it was significantly up-regulate in carbon metabolism and glycolysis (Figure S1A). Furthermore, we analyzed the relationship between proteins regulated in metabolism and MCODE through PPI (Fig. 4D, E, F) (Table 1). Hubba analysis can look the key protein. And by selected Degree and MCC Top10 through Hubba in Cytoscape (Fig. 4B, C), we found a total of 14 core proteins, of which 12 proteins have undergone regulation of crotonylation modification site (Table 2). These results indicate that with the development of the lung, the metabolism gradually increases, and the crotonylation modification may regulate the metabolism of the fetal lung.

Fig. 4
figure 4

Metabolism regulation in Developing Fetal Lung. A KEGG functional enrichment analysis of down-regulated DEPs. B The Top 10 of Degree in PPI of metabolism of DEPs. (C) The Top 10 of MCC in PPI of metabolism of DEPs. D Metabolism of DEPs of Enrichment_GO_Color by cluster. E Metabolism of DEPs of Enrichment_GO_Color by PValue. F Metabolism of DEPs of PPI and MCODE by cluster

Table 1 Metabolism of DEPs of PPI and MCODE by cluster
Table 2 Metabolism of DEPs of 12 proteins have undergone regulation of crotonylation modification site

Gene expression regulation in developing fetal lung

Cell proliferation requires storage materials, and gene expression is an important part of cell proliferation and plays a decisive role in development. It had been reported that mRNA related to cell cycle, DNA repair/replication, RNA processing, and translation is the main biological process that decreases as the lung matures [5]. Our data show that Cell cycle, DNA replication, Homologous recombination, Mismatch repair, Base excision repair, Nucleotide excision repair, Basal transcription factors, RNA polymerase, Ribosome biogenesis in eukaryotes, and Ribosome are all significantly up-regulated (Fig. 5A). Meanwhile the level of crotonylation modification shows that the Cell cycle and Ribosome have enriched in DCPs (Figure S1B). Some protein sites was up-regulated, and some protein sites are down-regulated. The biological significance of crotonylation regulation needs to explore. Meanwhile we also found that the protein levels of 8 histones were up-regulated, and the crotonylation sites of 5 histones regulated to varying degrees. K55, K59, and K82 of H1.0 are down. K120 of H1.1 is up. K49, K78, and K100 of H1.5 is down, K55 is down. K47 of H1.10 is up, K75 reversly. K92 and K78 of H4 are down, K80 reversely (Fig. 5d). We further analyzed the proteins in the Cell cycle, and selected Degree and MCC Top5 through Hubba (Fig. 5B, C), and found a total of 8 core proteins, of which Cdk family proteins and Mcm proteins are the most. We found that there are more proteins of crotonylation in nuclear through subcellular localization, especially proteins related to chromatin, such as Crebbp and Mcm family proteins (Fig. 5E). However, crotonylation modification related to chromosome recombination and gene expression. Therefore, we speculate that protein crotonylation modification may affect lung development through cell proliferation.

Fig. 5
figure 5

Gene Expression regulation in Developing Fetal Lung. A KEGG functional enrichment analysis of up-regulated DEPs. B The Top 5 of Degree in PPI of Cell cycle of DEPs. C The Top 5 of MCC in PPI of Cell cycle of DEPs. D Schematic diagram of gene expression of DEPs involved. E Schematic diagram of Cell cycle of DEPs and DCPs involved

Lung development regulation in developing fetal lung

As the fetus develops and grows, the lungs continue to develop and improve. We found 42 DEPs and 13 DCPs protein of GO Terms regulation in various parts of the lung through protein annotation (Tables 3 and 4). Lung epithelium development, lung cell differentiation, lung epithelial cell differentiation, lung alveolus development, lung secretory cell differentiation, lung morphogenesis, and epithelial tube branching involved in lung morphogenesis are more proteins regulation. There are also modifications in these biological processes. We found 24 and 14 terms of biological process related to lungs in DEPs and DCPs (Fig. 6A, B). Among them, the 41 DEPs and 13 DCPs are in the term of lung development. As a result, the four proteins (Hmgb1, Cux1, Itga3, and Asah1) have dysregulated in expression and crotonylated of protein. For example, in the nucleus, Hmgb1 is one of the major chromatin-associated non-histone proteins. It acts as a DNA chaperone involved in replication, transcription, chromatin remodeling, V(D)J recombination, DNA repair, and genome stability [10]. At the same time, we also found the regulation in transcription factor–for instance, Foxp1, Foxp4, Sox9, Nfib, Gata6, and Nkx2-1. Nfib, Gata6, and Nkx2-1 are participating in many lungs of GO terms. Regarding protein of crotonylation of the lung, Ctnnb1 and Yap1 of crotonylation protein are participating in many GO terms of lungs.

Table 3 Lung development regulation of 42 DEPs in developing fetal lung
Table 4 Lung development regulation of 13 DCPs in developing fetal lung
Fig. 6
figure 6

Lung Development regulation in Developing Fetal Lung. A, B GO enrichment analysis of lung related DEPs and DCPs, respectively. C KEGG functional enrichment analysis of (Notch, TGF-beta, and Wnt) signaling pathway DEPs and DCPs. D Schematic diagram of Notch signaling pathway and proteins up-regulated (red)

Six signaling pathways are involved in lung organogenesis, such as Notch, Tgfβ / Bmp, Sonic hedgehog (Shh), Fgf, Egf, and Wnt [11]. Studies have reported that the predicted inhibitory role of Notch signaling during lung maturation [5]. Our data is also significantly enriched in the Notch signaling pathway and up-regulate (Fig. 6D). In KEGG analysis, the 14 DEPs and 9 DCPs belong to the Wnt signaling pathway (Fig. 6C). Crebbp is only both dysregulated. The protein is up to expression and the position of K1524 occurred down in modification. Meanwhile, we found the 9 DEPs and 2 DCPs belong to Notch signaling pathway. Wherein, Crebbp is only both dysregulated. We also found the 11 DEPs and 7 DCPs belong to the Tgf-beta signaling pathway. The Crebbp and Dcn are both dysregulated. Therefore, Crebbp is the role of importance in these pathways; it maybe regulates lung development. In the PPI network of all the dysregulated proteins of expression, Crebbp interacts with many proteins. Maybe it can influence all the dysregulated proteins of expression. Hence, we suggest that Crebbp is a key protein of lung organogenesis; it can be a biomarker to detect lung development.

Discussion

In the present study, we performed LC–MS/MS-based proteomic analysis. We identified 2,645 DEPs and 1,156 DCPs between the fetal lung of 17 weeks and 12 weeks and 514 proteins with both expression and crotonylation difference. The lungs of the fetus at 12 weeks are in Pseudoglandular and at 17 weeks are in Canalicular. We revealed the regulation of protein related to lung development from the differences between the two phases.

Energy metabolism not only guarantees the metabolism of the fetus, but it is also the basis for the growth and development of fetal tissues and organs. Several new protein post-translational modifications that use intermediates in metabolism have been discovered, and interestingly, several of these modifications, in turn, regulate the activity of metabolic enzymes [8]. We analyzed the metabolic regulation of different proteins, in which lipid metabolism is consistent with what has been reported [5]. Our data shows that multiple metabolic pathways enhanced with the development of the lung. Secondly, we also found 11 core proteins from various metabolically regulated proteins, and most of these core proteins have undergone crotonylation modification [9]. Here, we suggest that if you need to study the metabolism of the lungs, you can start from these aspects, and you can carry out in-depth research from the perspective and direction of crotonylation modification.

Cell proliferation and differentiation are the conditions for the growth of fetal tissues and organs, and the basis for their functional development [12]. In the process of cell proliferation, the most involved is the gene expression in the cell. The results revealed by our data are consistent with previous reported. As the lungs develop, cell proliferation related protein express down, and we speculate it maybe decrease the rate of cell proliferation. Among them, Cdk2 and Cdk1 had reported as important genes in cell cycle control, and they are the core proteins in our cell cycle as same [4]. Secondly, our mini-chromosome maintenance protein Mcm family is at the core position [5]. It is reported that it decreases with development, and most of the proteins show regulation in crotonylation modification. We speculate that the crotonylation modification may affect the epigenetics of lung development through the Mcm family. The crotonylation firstly be found on histones and was related to genetic regulation [13]. Our data also detected regulation in the level of crotonylation modification in H1 and H4. These proteins may regulate some epigenetic change–specific needs further research.

Lung development has regulated by transcription factor [12]. Transcription factors determine the correct fate and differentiation of lung epithelial cells. We analyzed transcription factors related to lung development. Nkx2-1 is located at the center of the transcription network, which is essential for imposing the fate of lung epithelial cells through activation and inhibition mechanisms. Our data showed that Nkx2-1 expressed at 12 and 17 weeks, and it down regulated in lung development. A member of the Gata family of zinc finger proteins, highly expressed in the lung’s endoderm and vascular smooth muscle, and the loss of lung epithelial-specific Gata6 during development leads to defective epithelial differentiation and increased proliferation [14]. Our data show that Gata6 expression is present at 12 weeks and 17 weeks, and it down regulated with the development of lungs like Nkx2-1 [15]. We speculate that they have the same interaction relationship as reported in the literature. The Fox protein family plays an important role in many aspects of lung development [16]. Foxp1, Foxp2, and Foxp4 expressed in overlapping patterns in the developing lung and the postnatal lung. Our data show that Foxp1 and Foxp4 down regulated with lung development. There are also reports that Sox9 plays multiple roles in lung epithelial cells during branch morphogenesis, and Sox9 is essential for correct branch morphogenesis [17]. Our data show that Sox9 expression is at 12 and 17 weeks, and it down regulated with lung development. Hmgb protein is a non-histone chromatin-related molecule that affects transcription and cell differentiation through histone binding and RNA polymerase II activity [18]. Our data show that Hmgb1 has undergone crotonylation modification. Combined with the role of this protein, we believe that can study the epigenetic mechanism of lung development through Hmgb1.

Non-neurologically regulated body fluids or secreted factors, including growth factors or hormones, are one of the main factors that regulate the growth process of the fetus in the uterus. The molecular signaling pathways that regulate the different stages of lung development studied. Tgf-β, Wnt, hedgehog, Notch, and fibroblast growth factor (Fgf) signaling pathways are involved in regulating lung size, branching, and pattern during lung development. Our KEGG analysis enriched the Tgf-β, Wnt, and Notch signaling pathways. The expression of the three Tgf-β gene isoforms, receptors (I, II, and III), and signal transduction mediators (Smad-2, -3, -4, -6, and -7) are present throughout lung development [19]. Our data show that Smad3 and Smad4 in this signaling pathway different expressed as the lungs develop. In mice, Tgf-β1 showed to inhibit branch morphogenesis. Up-regulation of Tgf-β receptor II or Smad-2, -3, -4, and -6 also showed the effect of inhibiting branch morphogenesis. Our data showed that Smad3 and Smad4 had down-regulate with lung development, which can promote branch morphology Formation to meet the normal developmental needs of the lungs. The expected inhibitory effect is Notch signaling during lung maturation. Our data also shows that most of the proteins in the Notch signaling pathway are down-regulate with development. Crebbp gene encodes Creb-binding protein, Acetylates histones, giving a specific tag for transcriptional activation [20]. It also acetylates non-histone proteins, like Ddx21, Fbl, Irf2, Mafg, Ncoa3, Polr1e, Paf53, and Foxo1 [21], and Functions as a transcriptional coactivator for Smad4 in the Tgf-beta signaling pathway [22]. Because of a common protein has three important signal pathways, it also exists in the cell cycle and linked to various proteins in PPI. We suggest that it may be used to an important protein for lung research in the future.

However, there were some limitations to consider in the current study. First, it is a small number size cross-sectional study. Second, the results of DEPs and DCPs were two different samples at different developmental periods. The LC–MS/MS measurement could have neglected individual information. Third, although thousands of proteins had detected, there were no more than one thousand proteins with expression and crotonylation quantitative data. Some key proteins could have neglected key information. The most studied Fgf signaling pathway and Fgf10 expression of the protein were not detected in the lung [23]. Future studies should explore these limiting factors in more detail.

Conclusion

The quantitative proteomic analysis revealed differentially expressed and crotonylated proteins between 12 and 17 weeks of the fetal lung. The key features of DEPs and DCPs indicate that lung development is a complex process. We were the first to report on the epigenetic characteristics of fetal lung development and found some key proteins from it, laying a certain foundation for follow-up related research. Simultaneously, it can also be used as a reference for some clinical development-related diagnoses or provide candidate markers for diagnosis and treatment targets for the treatment of pulmonary dysplasia.

Details of material usage and methods

Sample collection and processing

All the schemes in this study are in line with the Interim Measures for the Management of Human Genetic Resources of the Ministry of Science and Technology of China. The ethics committee approved the human embryo collection and research protocol of Shenzhen People's Hospital (LL-KY-2019591). This study strictly abides by Shenzhen People’s Hospital’s laws and institutional ethics regulations on selective termination of pregnancy specimens. After the donor patients signed the informed consent document, the lung tissue samples of the developing fetus had collected. Process the tissue sample immediately after separation. Transfer the tissue to a sterile 10 cm2 tissue culture dish, cut into pieces < 0.6 cm3, and then transfer to a 50 mL conical tube.

Protein extraction

Under the guidance of experimental method by Xu et al [24], we carried out protein extraction. In short, prepare liquid nitrogen, lysis buffer A (8 M urea, 1% Protease Inhibitor Cocktail) for Proteomics, lysis buffer B (8 M urea, 1% Protease Inhibitor Cocktail, 3 μM TSA, 50 mM NAM) for Crotonylation and a mortar (pre-cooled with liquid nitrogen in advance) and an ultrasound system. Remove the sample from -80 °C, weigh 5 mg of tissue sample into a mortar, and then add liquid nitrogen to grind until the sample becomes a powder. Next, divide the sample into two parts 1:4, add the smaller part with 4 times the powder lysis buffer A, and the larger part with 4 times the powder lysis buffer B, and lyse the tissue by ultrasound. Finally, the cell debris was removed (4 °C, 12,000 g centrifugation for 10 min), and the protein concentration of the supernatant was determined using the BCA kit.

BCA Reagent A, 500 ml, containing sodium carbonate, sodium bicarbonate, bicinchoninic acid and sodium tartrate in 0.1 M sodium hydroxide.

BCA Reagent B, 25 ml, containing 4% cupric sulfate Albumin Standard Ampules, 2 mg/ml, 10 × 1 ml ampules, containing bovine serum albumin (BSA) at 2.0 mg/ml in 0.9% saline and 0.05% sodium azide.

Trypsin digestion

We used our own experimental method to trypsin digestion by Lu et al [25]. Overall, prepare 5 mM dithiothreitol and 11 mM iodoacetamide for alkylation, 100 mM triethylammonium bicarbonate (TEAB) buffer (Sigma-Aldrich; Merck KGaA). First, the protein solution was reduced with 5 mM dithiothreitol for 30 min at 56 °C, and then diluted with TEAB to obtain a urea concentration of < 2 M. Finally, digest twice (the first digestion is overnight at 37 °C with a trypsin to protein mass ratio of 1:50, and the second digestion is 4 h at 1:100).

Affinity enrichment

We adopted our own experimental method to affinity enrichment by Lin et al [26]. Prepare NETN buffer (100 mm NaCl; 1 mm EDTA; 50 mm Tris–HCl; 0.5% NP-40; pH 8.0) to enrich the crotonylated peptides. Dissolve trypsin peptides and pre-washed antibody beads (PTM Biolabs, Inc.) in NETN buffer, incubate gently at 4 °C, and shake overnight. Wash the beads four times (with NETN buffer) and then twice (with H2O). Next is the classification, drying and desalting of the peptides. The peptides bound to the beads were eluted (using 0.1% trifluoroacetic acid), combined and vacuum dried, and the resulting peptides were desalted (using C18 ZipTips (EMD Millipore) for liquid chromatography-MS (LC–MS/MS) analysis).

LC–MS/MS analysis

We use the EASY-nLC 1000 ultra-high performance liquid chromatography system (Thermo Fisher Scientific, Inc.) and the Q Exactive Plus instrument of the ultra-high performance liquid chromatograph (Thermo Fisher Scientific, Inc.). To prepare solvent A (0.1% formic acid reagent), Solvent B (0.1% formic acid in 98% acetonitrile) and a self-made reversed-phase analytical column (length, 15 cm; ID 75 µm), used to dissolve tryptic peptides, and load the sample onto the reversed-phase analytical column. Settings: chromatographic system (at a constant flow rate of 400 nL/min), gradient: solvent B, increased from 6 to 23% in 26 min, increased from 23 to 35% in 8 min, increased to 80% in 3 min and then maintained at 80% in the last 3 min. Use nanospray positive ion source to process peptides, MS/MS analysis: connect Q Exactive Plus instrument online; and settings: 2.0 kV (electrospray voltage) 350–1 (m/z scan range). The first step is to check 800 complete scans and complete peptides (using Orbitrap with a resolution of 70,000), and the second step is to set 28 (standard collision energy) MS/MS peptides, and analyze fragments (using Orbitrap with a resolution of 17,500). Rate), the third step is to perform the data-related process (alternating between MS scan and 20 MS/MS scans and 15.0 s of dynamic exclusion), and finally, set the automatic gain control to 5E4 and the fixed first quality to 100 m/z.

Database search

Trypsin/P is designated as a lyase and allows ≤ 4 deletions to be cleavage. Search engine: Maxquant search engine (v.1.5.2.8), crotonylation database (https://www.uniprot.org/), analyze MS/MS data, and connect with reverse decoy database to search serial MS data. Search settings: the first time 20 ppm (mass tolerance of precursor ions), 5 ppm (main search), 0.02 Da (mass tolerance of fragment ions). It based on fixed modification (carbamylmethylation of cystine) and variable modification (crotonylation and oxidation of methionine). For the false discovery rate, we adjusted it to < 1% and set the minimum score to > 40.

Bioinformatics analysis to proteins

We refer to the analysis method of Chen et al [27]. First, perform gene ontology (GO) annotation by UniProt-GOA database (www. http://www.ebi.ac.uk/GOA/, gaf-version, 2.0) on all identified proteins, subcellular localization by WoLF PSORT (https://www.genscript.com/wolf-psort.html), and annotation of kyoto encyclopedia of genes and genomes (KEGG) pathway. Then, we performed GO and KEGG functional enrichment of the protein (statistical method: Two-tailed Fisher's exact test). Finally, we use the one-way hierarchical clustering in Genesis to visualize the cluster members through the heatmap using the "heatmap.2" function in the "gplots" R-package.

Motif analysis

We adopted the motif analysis method of Huang et al [28]. We used Motif-x software (http://motif-x.med.harvard.edu/) to model 21-mers (10 amino acids upstream and downstream of the site), and all identified crotonylated modified peptides Sequence analysis. All protein sequences in the database used to background database parameters, and the default values used for other parameters.

Protein–protein interaction analysis

Protein–protein interaction (PPI) had conducted using a procedure described by Huang et al [28]. Identified proteins had performed by PPI analysis using Cytoscape software (Version 3.3.0). The search tool for retrieving mutual genes/proteins (STRING) database (the confidence score by using a metric) to define interaction confidence can obtain PPI networks. We fetched high confidence (confidence score ≥ 0.9) of all mutual. Among them, MCC and Degree are good tools to analyze core proteins.

Availability of data and materials

Uploading.

References

  1. Cotten CM. Pulmonary hypoplasia. Semin Fetal Neonatal Med. 2017;22(4):250–5.

    Article  PubMed  Google Scholar 

  2. Schittny JC. Development of the lung. Cell Tissue Res. 2017;367(3):427–44.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Du Y, Clair GC, Al Alam D, Danopoulos S, Schnell D, Kitzmiller JA, Misra RS, Bhattacharya S, Warburton D, Mariani TJ, Pryhuber GS, Whitsett JA, Ansong C, Xu Y. Integration of transcriptomic and proteomic data identifies biological functions in cell populations from human infant lung. Am J Physiol Lung Cell Molecular Physiol. 2019;317(3):L347-l360.

    Article  CAS  Google Scholar 

  4. Bhattacharya S, Mariani TJ. Systems biology approaches to identify developmental bases for lung diseases. Pediatr Res. 2013;73(4 Pt 2):514–22.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  5. Xu Y, Wang Y, Besnard V, Ikegami M, Wert SE, Heffner C, Murray SA, Donahue LR, Whitsett JA. Transcriptional programs controlling perinatal lung maturation. PLoS One. 2012;7(8):e37046.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Cox B, Kislinger T, Wigle DA, Kannan A, Brown K, Okubo T, Hogan B, Jurisica I, Frey B, Rossant J, Emili A. Integrated proteomic and transcriptomic profiling of mouse lung development and Nmyc target genes. Mol Syst Biol. 2007;3:109.

    Article  PubMed  PubMed Central  Google Scholar 

  7. Swaffer MP, Jones AW, Flynn HR, Snijders AP, Nurse P. CDK Substrate phosphorylation and ordering the cell cycle. Cell. 2016;167(7):1750-1761.e16.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Lin H, Su X, He B. Protein lysine acylation and cysteine succination by intermediates of energy metabolism. ACS Chem Biol. 2012;7(6):947–60.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Fang Y, Xu X, Ding J, Yang L, Doan MT, Karmaus PWF, Snyder NW, Zhao Y, Li JL, Li X. Histone crotonylation promotes mesoendodermal commitment of human embryonic stem cells. Cell Stem Cell. 2021;28(4):748-763.e7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Wei J, Alfajaro MM, DeWeirdt PC, Hanna RE, Lu-Culligan WJ, Cai WL, Strine MS, Zhang SM, Graziano VR, Schmitz CO, Chen JS, Mankowski MC, Filler RB, Ravindra NG, Gasque V, de Miguel FJ, Patil A, Chen H, Oguntuyo KY, Abriola L, Surovtseva YV, Orchard RC, Lee B, Lindenbach BD, Politi K, van Dijk D, Kadoch C, Simon MD, Yan Q, Doench JG, Wilen CB. Genome-wide CRISPR Screens Reveal Host Factors Critical for SARS-CoV-2 Infection. Cell. 2021;184(1):76-91.e13.

    Article  CAS  PubMed  Google Scholar 

  11. Herriges M, Morrisey EE. Lung development: orchestrating the generation and regeneration of a complex organ. Development (Cambridge, England). 2014;141(3):502–13.

    Article  CAS  PubMed  Google Scholar 

  12. Warburton D, El-Hashash A, Carraro G, Tiozzo C, Sala F, Rogers O, De Langhe S, Kemp PJ, Riccardi D, Torday J, Bellusci S, Shi W, Lubkin SR, Jesudason E. Lung organogenesis. Curr Top Dev Biol. 2010;90:73–158.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Tan M, Luo H, Lee S, Jin F, Yang JS, Montellier E, Buchou T, Cheng Z, Rousseaux S, Rajagopal N, Lu Z, Ye Z, Zhu Q, Wysocka J, Ye Y, Khochbin S, Ren B, Zhao Y. Identification of 67 histone marks and histone lysine crotonylation as a new type of histone modification. Cell. 2011;146(6):1016–28.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Morrisey EE, Ip HS, Lu MM, Parmacek MS. GATA-6: a zinc finger transcription factor that is expressed in multiple cell lineages derived from lateral mesoderm. Dev Biol. 1996;177(1):309–22.

    Article  CAS  PubMed  Google Scholar 

  15. Zhang Y, Rath N, Hannenhalli S, Wang Z, Cappola T, Kimura S, Atochina-Vasserman E, Lu MM, Beers MF, Morrisey EE. GATA and Nkx factors synergistically regulate tissue-specific gene expression and development in vivo. Development (Cambridge, England). 2007;134(1):189–98.

    Article  CAS  PubMed  Google Scholar 

  16. Shu W, Lu MM, Zhang Y, Tucker PW, Zhou D, Morrisey EE. Foxp2 and Foxp1 cooperatively regulate lung and esophagus development. Development (Cambridge, England). 2007;134(10):1991–2000.

    Article  CAS  PubMed  Google Scholar 

  17. Rockich BE, Hrycaj SM, Shih HP, Nagy MS, Ferguson MA, Kopp JL, Sander M, Wellik DM, Spence JR. Sox9 plays multiple roles in the lung epithelium during branching morphogenesis. Proc Natl Acad Sci USA. 2013;110(47):E4456–64.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Furst DE, Clements PJ, Hillis S, Lachenbruch PA, Miller BL, Sterz MG, Paulus HE. Immunosuppression with chlorambucil, versus placebo, for scleroderma. Results of a three-year, parallel, randomized, double-blind study. Arthritis Rheum. 1989;32(5):584–93.

    Article  CAS  PubMed  Google Scholar 

  19. Aschner Y, Downey GP. Transforming growth factor-β: master regulator of the respiratory system in health and disease. Am J Respir Cell Mol Biol. 2016;54(5):647–55.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Das C, Roy S, Namjoshi S, Malarkey CS, Jones DN, Kutateladze TG, Churchill ME, Tyler JK. Binding of the histone chaperone ASF1 to the CBP bromodomain promotes histone acetylation. Proc Natl Acad Sci USA. 2014;111(12):E1072–81.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Iyer-Bierhoff A, Krogh N, Tessarz P, Ruppert T, Nielsen H, Grummt I. SIRT7-dependent deacetylation of fibrillarin controls histone H2A methylation and rRNA synthesis during the cell cycle. Cell Rep. 2018;25(11):2946-2954.e5.

    Article  CAS  PubMed  Google Scholar 

  22. Yang Y, Cui J, Xue F, Zhang C, Mei Z, Wang Y, Bi M, Shan D, Meredith A, Li H, Xu ZQ. Pokemon (FBI-1) interacts with Smad4 to repress TGF-β-induced transcriptional responses. Biochem Biophys Acta. 2015;1849(3):270–81.

    CAS  PubMed  Google Scholar 

  23. Yin Y, Ornitz DM. FGF9 and FGF10 activate distinct signaling pathways to direct lung epithelial specification and branching. Sci Signal. 2020;13(621):eaay4353.

  24. Xu JY, Zhang C, Wang X, Zhai L, Ma Y, Mao Y, Qian K, Sun C, Liu Z, Jiang S, Wang M, Feng L, Zhao L, Liu P, Wang B, Zhao X, Xie H, Yang X, Zhao L, Chang Y, Jia J, Wang X, Zhang Y, Wang Y, Yang Y, Wu Z, Yang L, Liu B, Zhao T, Ren S, Sun A, Zhao Y, Ying W, Wang F, Wang G, Zhang Y, Cheng S, Qin J, Qian X, Wang Y, Li J, He F, Xiao T, Tan M. Integrative Proteomic Characterization of Human Lung Adenocarcinoma. Cell. 2020;182(1):245-261.e17.

    Article  CAS  PubMed  Google Scholar 

  25. Lu YP, Zhang XL, Zheng F, Yun C, Zhu C, Cai W, Liu D, Hong X, Li Q, Hu B, Tang D, Yin LH, Dai Y. Quantitative proteomic analyses to reveal the key features of proteins in new onset ankylosing spondylitis patients. ACS Omega. 2020;5(32):20153–61.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Lin H, Tang D, Xu Y, Zhang R, Ou M, Zheng F, Chen J, Zhang Y, Zou G, Xue W, Zou Y, Dai W, Sui W, Dai Y. Quantitative analysis of protein crotonylation identifies its association with immunoglobulin A nephropathy. Mol Med Rep. 2020;21(3):1242–50.

    CAS  PubMed  PubMed Central  Google Scholar 

  27. Chen W, Tang D, Xu Y, Zou Y, Sui W, Dai Y, Diao H. Comprehensive analysis of lysine crotonylation in proteome of maintenance hemodialysis patients. Medicine. 2018;97(37):e12035.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Huang S, Zheng F, Liu L, Meng S, Cai W, Zhang C, Dai W, Liu D, Hong X, Tang D, Dai Y. Integrated proteome and phosphoproteome analyses of peripheral blood mononuclear cells in primary Sjögren syndrome patients. Aging. 2020;13(1):1071–95.

    Article  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

We thank Li Chen from Hangzhou LC-BIO Co., Ltd for help with the data analysis and programming in R.

Funding

This work was Supported by Guangdong Provincial Engineering Research Center of Autoimmune Disease Precision Medicine, Science and Technology Planning Project of Guangdong Province, China (No.2017B020209001), grants Sanming project of medicine in Shenzhen (SZSM201812078) and the science and technology plan of Shenzhen(NO.JCYJ20220530152015035).

Author information

Authors and Affiliations

Authors

Contributions

W.W. interpreted the data and wrote the main manuscript text. W.H. prepared all the figures. W.S. and Y.W. performed the experiments. Z.Z. contributed to the interpretation of data. Y.Y. and P.L. reviewed the manuscript. Y.C. performed the data analysis. D.T. and Y.D. supervised the experiments and contributed to the interpretation of data.

Corresponding authors

Correspondence to Donge Tang or Yong Dai.

Ethics declarations

Competing interests

The authors declare no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Figure S1.

 (A) KEGG functional enrichment analysis of up-regulated DCPs. (B) KEGG functional enrichment analysis of DCPs.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Wang, W., Shi, W., Wang, Y. et al. Systematic proteomics profiling of lysine crotonylation of the lung at Pseudoglandular and Canalicular phases in human fetus. Proteome Sci 21, 22 (2023). https://doi.org/10.1186/s12953-023-00215-8

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12953-023-00215-8

Keywords