ORIGINAL RESEARCH article
Whole-Genome Methylation Study of Congenital Lung Malformations in Children
- 1Medical, Surgical and Health Sciences Department, University of Trieste, Trieste, Italy
- 2Pediatric Surgery, Institute for Maternal and Child Health–IRCCS “Burlo Garofolo”, Trieste, Italy
- 3Laboratory of Medical Genetics, Institute for Maternal and Child Health–IRCCS “Burlo Garofolo”, Trieste, Italy
Background and Objectives: The treatment of asymptomatic patients with congenital pulmonary malformations (CPMs) remains controversial, partially because the relationship between congenital lung malformations and malignancy is still undefined. Change in methylation pattern is a crucial event in human cancer, including lung cancer. We therefore studied all differentially methylated regions (DMRs) in a series of CPMs in an attempt to find methylation anomalies in genes already described in association with malignancy.
Methods: The DNA extracted from resected congenital lung malformations and control lung tissue was screened using Illumina MethylationEPIC arrays. Comparisons between the group of malformed samples or the malformed samples of same histology or each malformed sample and the controls and between a pleuropulmonary blastoma (PPB) and controls were performed. Moreover, each malformed sample was pairwise compared with its respective control. All differentially methylated regions (DMRs) with an adjusted p-value <0,05 were studied.
Results: Every comparison highlighted a number of DMRs closed to genes involved either in cell proliferation or in embryonic development or included in the Cancer Gene Census. Their abnormal methylation had been already described in lung tumors.
Conclusions: Methylation anomalies already described in lung tumors and also shared by the PPB were found in congenital lung malformations, regardless the histology. The presence of methylation abnormalities is suggestive of a correlation between congenital lung malformations and some step of malignant transformation.
There is a general consensus that symptomatic congenital pulmonary malformations (CPMs) should be removed surgically. However, the treatment of asymptomatic patients remains controversial (1). Some authors recommend prophylactic pulmonary resection to avoid the long-term risk of pulmonary recurrent infections, pneumothorax, or development of lung malignancy (1, 2), while others suggest a conservative approach based on the observation of the patient (1, 3).
The relationship between congenital lung malformations and malignancy remains undefined and continues to be a critical consideration in surgical decision making. Hartman and Stochat reported that 4% of pulmonary malignant tumors were associated with congenital cystic malformations. Tumors developing within these malformations included rhabdomyosarcoma, pleuropulmonary blastoma (PPB), adenocarcinoma, squamous cell carcinoma, and mesenchymoma (4). Later, Ozcan et al. reported 29 cases of primary rhabdomyosarcoma, 15 of which arose in a preexisting congenital lung malformations (5). Nasr et al. (6) found 2% of association between PPB and congenital pulmonary malformations. Recently, in a systematic review (7) we highlighted 168 cases, 76 children and 92 adults, in whom a lung tumor was found in association with a CPM. We concluded that all histological types of CPMs could be associated with malignant lung lesions and that the malignant transformation could happen at any age.
One hallmark of cancer cells is their completely different methylation pattern. In many malignant tumors, the levels of methylation are decreased, while promoter regions of important regulatory and tumor suppressor genes are hypermethylated and therefore silenced. Hypermethylation associated with tumor suppressor genes is uncommon in normal cells. However, it is widely represented in cancer cells (8, 9). Abnormal DNA hypomethylation has been demonstrated to also play an important role in tumor development, both increasing genome instability (10) and activating the transcription of oncogenes that are normally silenced (11).
The aim of this study was to investigate the possible biologic relationship between congenital pulmonary malformations and lung tumors. Using Illumina MethylationEPIC array analysis that is easy to use, time efficient, and cost effective technique (12), we studied all differentially methylated regions (DMRs) in a series of congenital lung malformations in an attempt to find methylation anomalies in genes already described in association with malignancy.
Material and Methods
All congenital lung malformations resected at the Institute for Maternal and Child Health–IRCCS “Burlo Garofolo” (Trieste, Italy) from January 2010 to January 2019 were assessed. After approval by the Institutional Ethical Committee, the medical records of the patients were analyzed.
Lung biopsies from the resected malformed lobes were snap-frozen and stored at −80°C. Samples used as control included biopsies from macroscopically normal lung tissue adjacent to the malformation of seven patients, three with intralobar sequestration (ILS) and four with congenital pulmonary airway malformation type 2 (CPAM2) and from the lung of a patient who was thought to have a congenital malformation until histological analysis proved the tissue to be normal. Lung biopsies from a patient with PPB were also analyzed. We were not able to recruit other patients with the same tumor firstly because of its rarity and, secondly, because it is frequently mistaken with other diagnostic entities and formalin-fixed and paraffin-embedded. This technique of storage could lower the quality of the DNA, and therefore the comparison with higher quality DNA from fresh frozen tissues could generate artifacts. Moreover, pediatric patients with lung malignancies are also absent from public databases such as The Cancer Genome Atlas (TCGA) (13) or ENCODE (14). It is also not possible to compare methylation data from adult lung malignancies to our pediatric malformations because the methylation pattern changes drastically with age (15).
Genomic DNA was extracted from lung tissue samples (16), and its concentration was measured with Qubit dsDNA Broad Range Assay Kit (Thermo Fisher Scientific). Of each DNA sample, 1 μg was bisulfite-converted with EZ DNA Methylation Kit and screened using MethylationEPIC Beadchips according to the manufacturer’s instructions (Zymo Research and Illumina Inc. respectively). Raw methylation data were analyzed with R version 3.6.2 (2019-12-12), using package ChAMP (17) (Chip Analysis Methylation Pipeline) version 2.16.1 in Rstudio.
Standard parameters of function champ.load were used to load and filter the dataset. After normalization with BMIQ method, the effect of two confounding variables (age of the samples and beadchip of origin) was removed using function removeBatchEffect from package Limma. All analyses were performed on the corrected dataset.
Significant differentially methylated regions (DMRs) were then calculated using the Bumphunter algorithm applying 1,000 permutations to each comparison, and the p-value was adjusted for multiple testing.
In our first analysis we compared the cohort of malformed samples with controls, looking for changes to the methylation pattern common to all malformations, regardless of histology. Then, we compared groups of the same histology with controls to identify any methylation signatures shared by cases of the same histological profile. Finally, we compared each individual malformed sample with the group of controls for individual differences. For eight samples that had a control tissue from the same patient available, we also performed a pairwise DMR analysis, using the standard parameters of R package DMRforPairs (18). Furthermore, the methylation profile of the malformed samples was compared with that of the single tumor sample to highlight any similarities.
At each comparison, the identified DMRs were classified according to their location, mapping to either the gene body, the TSS (transcription start site), the 5′ UTR (untranslated region) or the 3′ UTR of protein-coding genes of interest (that we refer to in the manuscript as GOI); or mapping to an intergenic region that according to the ENCODE database contains a candidate cis-regulatory element (ccRE) close to a GOI. These genes were either included in the Cancer Gene Census (CGC, Cosmic (19) or implicated in embryonic development or involved in cell proliferation (both according to the AmiGO 2 database (20), and the overlap between the three groups of genes was also taken into consideration.
Raw data was uploaded to Gene Expression Omnibus as dataset GSE174625.
Eighteen patients, nine girls and nine boys, all Caucasian except for one of African ethnicity, who underwent lung resection for congenital lung malformations, were considered. Thirteen of them had prenatal diagnosis. Ten patients remained asymptomatic, while eight had different degrees of respiratory infection. Six patients had intralobar sequestration (ILS), nine had congenital pulmonary airway malformations (CPAM) associated with extralobar sequestration (ELS) in two cases. Clinical data and histology are summarized in Table 1.
No clear separation between the malformed and the control samples was found on a multidimensional scaling (MDS) plot representing the similarity of all samples based on the signal of the top 1,000 most variable probes. They form a homogeneous group (Figure 1A) completely separated from the PPB sample (Figure 1B).
Figure 1 Multidimensional scaling (MDS) plot of (A) malformed and control samples; (B) malformed, control and pleuropulmonary blastoma samples. In each plot, the X axis represents Dimension 1 and the Y axis Dimension 2 of the MDS statistical analysis, which better expresses the mathematical distance between the samples.
The comparison between all the malformed samples and the controls identified 10 statistically significant DMRs (with and adjusted p-value < 0,05). Among them, one involved a ccRE located in an intergenic region near exon 1 of the gene ZFP57. In this region the mean methylation beta value was slightly higher in the malformed group than in the control group (Figure 2).
Figure 2 Plot of the beta methylation values (Y axis) of each probe (X axis) of the DMR located near the gene ZFP57. The beta values range from 1 (completely methylated) to 0 (completely unmethylated). The position of the probes inside the DMR is ordered according to their genomic position.
We also checked whether samples with the same histology have DMRs in common and that differ from controls. We found 10 significant DMRs in the ILS samples, none of which was close to a GOI. In the other histological type there were no DMRs near genes involved in cancer, but some were close to other GOIs, related to cell proliferation or embryonic development: in particular, seven out of 31 in the ELS samples, four out of 19 in the CPAM1 samples, two out of 14 in the CPAM2 samples, and two out of 20 in the CPAM3 samples (Figure 3 and Table 2).
Figure 3 Number of significant DMRs (adjusted p-value < 0,05), close to GOI (black) and to different genes (gray), identified when the cases of same histology were compared with controls.
Table 2 GOI close to significant DMRs elicited by the comparison between cases of same histology with controls.
When each malformed sample was compared with the controls, the number of DMRs ranged from eight, in samples 4 and 18B, to 100 in sample 10B (Figure 4). Nine of them were repeated at least three times and were localized close to genes which act in cell proliferation or in embryonic development or are included in the CGC (Table 3). The percentage of DMRs that each sample had in common with the PPB sample ranges from 30% in sample 1 to 84% in sample 2 (Figure 5).
Figure 4 Number of significant DMRs (adjusted p-value < 0,05), close to GOI (black) and to different genes (gray), identified comparing each malformed case with controls.
Table 3 GOI close to significant DMRs repeated at least three times when each malformed sample was compared to the controls.
Figure 5 Bar plot representing the percentage of significant DMRs (adjusted p-value < 0,05) shared by each sample with the PPB.
The pairwise DMR analysis, limited to the samples that had a control tissue from the same patient available, identified statistically significant DMRs in three out of eight samples: 23 in sample 10B, four in sample 12, and two in sample 17. The number of DMRs near the genes of interest was respectively thirteen, two, and one (Table 4).
Table 4 GOI close to significant DMRs elicited by the pairwise comparison between each malformed sample and the corresponding control tissue.
Congenital pulmonary malformations are developmental abnormalities of the lung occurring in approximately one to 4.2 per 10,000 births (21). Postnatal presentation varies from severe respiratory distress to complete lack of symptoms (22). Although all authors agree that surgical resection is the standard of care for symptomatic cases, the management of asymptomatic lesions remains controversial. The main reason for recommending elective surgery is late development of complications, including the highly debated malignancy. The relationship between congenital lung malformations and malignancy remains unknown and continues to be a critical consideration in surgical decision making. Some studies have addressed this issue in an attempt to find premalignant characteristics in congenital lung malformations.
Vargas and colleagues analyzed the karyotype and looked for p53 mutations in congenital cystic adenomatoid malformations. As they did not find cytogenetic abnormalities in congenital cystic adenomatoid malformations by conventional karyotype analysis, they concluded that congenital lung anomalies were non-neoplastic (23). Later, Rossi and colleagues found expression of mucins and K-RAS mutations in mucinogenic proliferations of congenital pulmonary airway malformations and concluded that these findings supported the neoplastic nature of the malformations (24). Recently, Hsu et al. (25) described in the blood of 19 cases of CPAM several damaging variants in genes involved in lung carcinoma.
The aim of this study was to look for further insight into the biologic relationship between congenital pulmonary malformations and lung tumors. As the change in methylation pattern is a pivotal event in human cancer including lung cancer, in this study we focused on DNA methylation in a series of congenital lung malformations.
Firstly, we compared all the malformed samples with the controls identifying, among ten, only one DMR of interest, involving a candidate cis-regulatory element (ccRE) and located near exon 1 of the gene ZFP57, which acts in genome imprinting, regulation of gene expression, and cell signaling. It has recently been hypothesized that ZFP57 is a potential susceptibility gene for lung cancer development through the increase of IGF2 expression (26).
The analysis did not uncover a common methylation pattern since the significant DMRs were very few. This was not surprising, as the malformed samples were a very heterogeneous group, and the individual differences in the methylation patterns of each malformed sample, combined with the likely cellular heterogeneity of the samples (that was beyond the control of the investigators), masked the differences between the two groups.
Looking for a common methylation signature among the cases of CPMs, we then compared the cases of same histology with controls. The results highlighted seven DMRs in the ELS close to genes mostly involved in embryonic development (at the TSS of HOXB1, HOXD4) and cell proliferation (at the TSS of CTNNA1, NR2F2, HSF4, MEIS1, and at the gene body of PTN). Both HOXB1 and HOXD4 were abnormally expressed in lung tumors. HOXB1 was described down-regulated in lung cancer tissue (27) and, similarly, we found the DMR close to this gene was hypermethylated. HOXB1 is an anti-tumor gene acting by inhibiting the expression of survival oncogenic genes (28); its decreased expression had been correlated with cell invasion and proliferation and inhibited apoptosis in glioma (29). HOXD4 was reported up-regulated in lung squamous cell carcinoma resulting in more aggressive invasiveness of lung cancer cells (30). In our samples of ELS HOXD4 was hypermethylated in an opposite fashion compared to those in the literature. Either CTNNA1 (31) or NF2F2 (32) has been found significantly down-regulated in lung cancer. The DMR close to CTNNA1 was hypermethylated in our samples in line with previous report of the literature. NF2F2 was also known to play a role in angiogenesis and development and to contribute to transform a non-invasive lung cancer in an invasive one through its expression (33). MEIS1 might limit the proliferation of non-small-cell lung adenocarcinoma (34) and was described methylated in squamous cell carcinomas (35). PTN is highly expressed in embryonic and postnatal development, while it is quite down-regulated in adult life; however, it is strongly expressed in lung tumor and other types of cancer (36). Hypermethylation was found in our samples for the DMR close to MEIS1 and PTN.
Four DMRs were found in CPAM1 cases close to genes related to cell proliferation or embryonic development (at the TSS of PLD6, S100A13; MXS2, and at the 5′ UTR of TXNRD1). S100A13’s overexpression was associated with intratumoral angiogenesis and more aggressive invasive phenotype in non-small-cell lung cancer (37). In our samples the DMR was hypermethylated differently from those in the literature. MSX2 is a regulator of embryonic development, and it is involved in pancreatic and breast cancer, but recently it has been found that it plays a role also in lung adenocarcinoma (38). The DMR was hypermethylated in our samples.
We found two DMRs in CPAM2 cases, close to the already discussed intergenic region near ZFP57, acting in embryonic development, and the gene body of MEIS1, whose closed DMR was hypomethylated differently from the mechanism in literature.
Finally, CPAM3 cases presented two DMRs, one close to the TSS of MSX2 which was this time hypomethylated, and the other in an intergenic region involving a ccRE near two genes involved in proliferation that was hypomethylated as well. The first gene near the second DMR was PITX2, whose low methylation correlated with higher risk of lung cancer progression (39), and the other was ENPEP that was under-expressed in lung adenocarcinoma (40).
One more time, the heterogenicity of the malformed samples likely prevented to highlight a common methylation signature. However, the presence of methylation abnormalities in genes already described in association with malignancy is suggestive, especially considering their role in lung tumors.
Due to the heterogeneity in our malformed samples, we also checked for individual differences by comparing each sample with the control group. It yielded a wide range of DMRs, but nine of them recurring at least three times: three including one ccRE each, near the genes ZFP57, SIX1-SIX4-MNAT1, and MMP2-IRX5; four at the TSS of genes HOXB1, HOXD4, NR2F2 and MSX2; one at the 5′ UTR of IRAK4, and one at the gene body of MEIS1. ZFP57, HOXB1, HOXD4, NR2F2, MEIS1, and MSX2 were already described when the cases of same histology were compared with controls. Interestingly, the DMR close to HOXB1 was found hypomethylated and the one close to HOXD4 hypermethylated, in the opposite way to the literature. SIX1 and SIX4 were associated with increased risk of tumorigenesis when their expression was increased. SIX4 controlled the expression of oncogenes, and it correlated with higher stages of the tumor, poor survival in NSCLC, and worse rate of relapses in lung adenocarcinoma (41). MNAT1 might have a role in promoting the development of NSCLC (42). MMP2 mRNA and protein levels were found increased in NSCLC (43).
The pairwise DMR analysis found statistically significant DMRs in three out of eight samples, when comparing each one of them with the corresponding control tissue from the same patient. The number of DMRs near to genes of interest was 12 out of 23 in one of the ELS (sample 10B), two out of four in an ILS (sample 12), and one out of two in another ILS (sample 17). In the ELS, the DMRs of interest were located at the gene body or the TSS of genes involved in embryo development (HOXA3, HOXC6, HOXB1, TGFB1I1, BRD2), cell proliferation (CTNNA1, CTSZ, GPR37L1, S100A13, TSPAN32, FOXP2), or both processes (IGF2BP1, NR2F2). In one of the ILS (sample 12), one DMR of interest was hypomethylated with respect to the control tissue, and located at the gene body of SMAD6 that has been reported as a “master regulator” of lung adenocarcinoma, for which both oncogenic and tumor suppressing activities have been proposed (44). The second DMR of interest for the same sample was a hypermethylated region at the gene body of FOXP1, a transcription factor involved in cell proliferation and included in the CGC that according to Sheng et al (45), prevented the growth of lung adenocarcinoma cells by the suppression of chemokine signaling pathways. The DMR of interest in the other ILS (sample 17) was located at the TSS of HOXD4, a gene that has been discussed earlier, and it was less methylated in the malformed than in the control tissue, contrarily to our previous findings but in concordance with the literature.
Both types of single sample analyses have strengths and weaknesses. Comparison of one sample to the entire control group allowed the analysis of samples that do not have a control tissue from the same patient, but it could introduce artifacts due to inter-patient variability of the methylation profile. The pairwise analysis, on the other hand, was possible in this case for a smaller number of samples, but it excluded any underlying methylation alterations not strictly related to the malformation, thus producing more specific results.
Finally, we compared the PPB sample with the controls and considered to discuss, among all, the DMRs in common with the malformed samples, mapping to the 3′ UTR of CACNA1C, the 5′ UTR of HOXA5, the TSS of CTSZ, ESRP2, HAND2, HOXA2, HOXA3, MAGI2, TWIST1, ccREs near MMP2, MSX1-OTX1-RUNX1, TBX3–TBX5, and the gene body of SH3PXD2A and WT1).
CACNA1C was found down-regulated in lung cancer (46). CTSZ was involved in promoting NSCLC cell migration and invasion (47). ESRP2 maintained the epithelial phenotype, avoiding the epithelial to mesenchymal transition that contributed to metastases. In NSCLC it was inhibited (48). HAND2 was found either hypermethylated or hypomethylated in different stages of lung adenocarcinoma (49). HOXA2, HOXA3, and HOXA5 have been recognized as target for DNA methylation in lung cancer, and they promoted carcinogenesis, but also acted as tumor-suppressor factors (50). MAGI2 was reported to act as an anti-tumor in hepatocellular cancer and breast cancer; its down-regulation has been demonstrated in NSCLC (51). CpG islands associated with MSX1 and OTX1 were methylated in the majority of lung squamous cell carcinomas, while the ones associated with RUNX1 were methylated in more than 80% of lung adenocarcinomas, being well known that hypermethylation of CpG islands is a signature of malignant progression (35). SH3PXD2A’s increased expression in lung adenocarcinoma directly correlated with metastasis and worse prognosis for the patient (52). TBX5 and TBX3 were highly expressed in normal lungs, but significantly suppressed in lung adenocarcinoma (53). TWIST1 expression was found increased in lung cancer tissue (54). WT1 was described as oncogene in lung cancer, among other malignancies (55).
The DMRs of interest were spread throughout the cases, and we were not able to find a recurrent pattern of abnormalities in the different types of congenital lung malformations.
It has been demonstrated that lung cancer develops through the acquisition of alterations in oncogenes and tumor suppressor genes. Prolonged exposure to carcinogens that interact with various genetic susceptibility and/or resistance factors contributes to the accumulation of genomic alterations, including nucleotide substitutions, small insertions and deletions, and chromosomal rearrangements in human lung cancer (56). This mechanism of action may explain the cases of association between congenital pulmonary malformations, sometimes asymptomatic for many years, and lung tumors described in literature (7), in which the genetic susceptibility together with the progressive exposure to carcinogens might trigger the development of malignancy.
This study has highlighted some key points. First, methylation anomalies already described in lung tumors could be found in samples of congenital lung malformations, regardless the histology, both when compared to the control group and when compared to a specific control tissue. Second, some methylation anomalies of the congenital lung malformations were shared by the PPB. Third, it seems unlikely that the presence of methylation abnormalities, which have been reported in association with lung tumors, could be considered a casual event in congenital lung malformations, in which malignant transformation has been described (7). Pulmonary malformations are essentially due to a dysregulation of cell proliferation which, during organogenesis, creates an abnormal development of some areas of the lower respiratory tract. Our study indicates that, at least in part, this dysregulation is caused by some genes which, due to their role in the cell cycle, are also involved in some stages of tumor development.
Our study, however, had some limits that should be acknowledged, such as the small number of samples of congenital lung malformations included and their heterogenicity. Moreover, the difficulty in finding proper control samples has limited their number. Another limitation is the lack of material from pediatric adenocarcinoma of the lung, which is the second tumor for frequency that has been found in association with CPMs in children (7), and the impossibility to compare the methylation of pediatric CPMs to the available adult adenocarcinoma samples due to the age-related changes in the methylation pattern. To our knowledge, this study represents the first attempt to address the methylation anomalies in pediatric congenital lung malformations using a whole genome approach. We have also described some methylation changes which, in some genes, appeared to have the opposite sign to what is described in the literature regarding their expression. These results, however, should not be discarded, because the relationship between DNA methylation and gene expression could be more complicated than previously understood [for example, it can have different effects according to its position with respect to the gene (57)], and DNA methylation is only one of the factors that regulate the expression of a gene. Thus, since we did not study the gene expression in our samples, we cannot be certain that hypomethylation corresponds to hyperexpression or vice versa.
Methylation anomalies already described in lung tumors and also shared by the PPB were found in congenital lung malformations, regardless the histology. This is suggestive of a correlation between congenital lung malformations and some step of malignant transformation.
More detailed analysis of genetic and epigenetic interactions as well as functional interactions among genes altered in congenital pulmonary malformations will further provide insights into the molecular mechanism of lung carcinogenesis.
Data Availability Statement
The original contributions presented in the study are publicly available. This data can be found here https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE174625.
The studies involving human participants were reviewed and approved by IRCCS Burlo Garofolo. Written informed consent to participate in this study was provided by the participants’ legal guardian/next of kin.
FP and AD’A conceptualized and designed the study, drafted the initial manuscript, and reviewed and revised the manuscript. SP designed the data collection instruments, collected data, carried out the initial analyses, and reviewed and revised the manuscript. All authors contributed to the article and approved the submitted version.
This study was supported by the IRCCS “Burlo Garofolo” (RC 22/11). The funder/sponsor did not participate in the work.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
ccRE, cis-regulatory element; CGC, Cancer Gene Census; CPAM1, 2, 3, congenital pulmonary airway malformation types 1, 2, 3; CPM, congenital pulmonary malformation; DMR, differentially methylated region; DNA, deoxyribonucleic acid; ELS, extralobar sequestration; PPB, pleuropulmonary blastoma; GOI, genes of interest; ILS, intralobar sequestration; MDS, multidimensional scaling; TCGA, The Cancer Genome Atlas; TSS, transcription start site; UTR, untranslated region.
1. Lo AY, Jones S. Lack of Consensus Among Canadian Pediatric Surgeons Regarding the Management of Congenital Cystic Adenomatoid Malformation of the Lung. J Pediatr Surg (2008) 43:797–9. doi: 10.1016/j.jpedsurg.2007.12.016
3. Aziz D, Langer JC, Tuuha SE, Ryan G, Ein SH, Kim PC. Perinatally Diagnosed Asymptomatic Congenital Cystic Adenomatoid Malformation: to Resect or Not? J Pediatr Surg (2004) 39:329–34; discussion 329-34. doi: 10.1016/j.jpedsurg.2003.11.021
5. Ozcan C, Celik A, Ural Z, Veral A, Kandiloglu G, Balik E. Primary Pulmonary Rhabdomyosarcoma Arising Within Cystic Adenomatoid Malformation: A Case Report and Review of the Literature. J Pediatr Surg (2001) 36:1062–5. doi: 10.1053/jpsu.2001.24747
6. Nasr A, Himidan S, Pastor AC, Taylor G, Kim PC. Is Congenital Cystic Adenomatoid Malformation a Premalignant Lesion for Pleuropulmonary Blastoma? J Pediatr Surg (2010) 45:1086–9. doi: 10.1016/j.jpedsurg.2010.02.067
7. Casagrande A, Pederiva F. Association Between Congenital Lung Malformations and Lung Tumors in Children and Adults: A Systematic Review. J Thorac Oncol (2016) 11:1837–45. doi: 10.1016/j.jtho.2016.06.023
9. Liu XR, Zhang RY, Gong H, Rugo HS, Chen LB, Fu Y, et al. Methylome Variation Predicts Exemestane Resistance in Advanced Er(+) Breast Cancer. Technol Cancer Res Treat (2020) 19:1533033819896331. doi: 10.1177/1533033819896331
12. Pidsley R, Zotenko E, Peters TJ, Lawrence MG, Risbridger GP, Molloy P, et al. Critical Evaluation of the Illumina Methylationepic BeadChip Microarray for Whole-Genome DNA Methylation Profiling. Genome Biol (2016) 17:208. doi: 10.1186/s13059-016-1066-1
17. Tian Y, Morris TJ, Webster AP, Yang Z, Beck S, Feber A, et al. ChAMP: Updated Methylation Analysis Pipeline for Illumina Beadchips. Bioinformatics (2017) 33:3982–4. doi: 10.1093/bioinformatics/btx513
18. Rijlaarsdam MA, van der Zwan YG, Dorssers LC, Looijenga LH. DmrforPairs: Identifying Differentially Methylated Regions Between Unique Samples Using Array Based Methylation Profiles. BMC Bioinf (2014) 15:141. doi: 10.1186/1471-2105-15-141
19. Sondka Z, Bamford S, Cole CG, Ward SA, Dunham I, Forbes SA. The COSMIC Cancer Gene Census: Describing Genetic Dysfunction Across All Human Cancers. Nat Rev Cancer (2018) 18:696–705. doi: 10.1038/s41568-018-0060-1
21. Raychaudhuri P, Pasupati A, James A, Whitehead B, Kumar R. Prospective Study of Antenatally Diagnosed Congenital Cystic Adenomatoid Malformations. Pediatr Surg Int (2011) 27:1159–64. doi: 10.1007/s00383-011-2909-1
22. Hammond PJ, Devdas JM, Ray B, Ward-Platt M, Barrett AM, McKean M. The Outcome of Expectant Management of Congenital Cystic Adenomatoid Malformations (CCAM) of the Lung. Eur J Pediatr Surg (2010) 20:145–9. doi: 10.1055/s-0030-1249047
23. Vargas SO, Korpershoek E, Kozakewich HP, de Krijger RR, Fletcher JA, Perez-Atayde AR. Cytogenetic and p53 Profiles in Congenital Cystic Adenomatoid Malformation: Insights Into Its Relationship With Pleuropulmonary Blastoma. Pediatr Dev Pathol (2006) 9:190–5. doi: 10.2350/06-01-0025.1
24. Rossi G, Gasser B, Sartori G, Migaldi M, Costantini M, Mengoli MC, et al. MUC5AC, Cytokeratin 20 and HER2 Expression and K-RAS Mutations Within Mucinogenic Growth in Congenital Pulmonary Airway Malformations. Histopathology (2012) 60:1133–43. doi: 10.1111/j.1365-2559.2011.04170.x
25. Hsu JS, Zhang R, Yeung F, Tang CSM, Wong JKL, So MT, et al. Cancer Gene Mutations in Congenital Pulmonary Airway Malformation Patients. ERJ Open Res (2019) 5:00196–2018. doi: 10.1183/23120541.00196-2018
26. Lamontagne M, Joubert P, Timens W, Postma DS, Hao K, Nickle D, et al. Susceptibility Genes for Lung Diseases in the Major Histocompatibility Complex Revealed by Lung Expression Quantitative Trait Loci Analysis. Eur Respir J (2016) 48:573–6. doi: 10.1183/13993003.00114-2016
30. Bao L, Zhang Y, Wang J, Wang H, Dong N, Su X, et al. Variations of Chromosome 2 Gene Expressions Among Patients With Lung Cancer or Non-Cancer. Cell Biol Toxicol (2016) 32:419–35. doi: 10.1007/s10565-016-9343-z
31. Srivastava M, Khurana P, Sugadev R. Lung Cancer Signature Biomarkers: Tissue Specific Semantic Similarity Based Clustering of Digital Differential Display (DDD) Data. BMC Res Notes (2012) 5:617. doi: 10.1186/1756-0500-5-617
32. Ning J, Li P, Zhang B, Han B, Su X, Wang Q, et al. miRNAs Deregulation in Serum of Mice is Associated With Lung Cancer Related Pathway Deregulation Induced by PM2. 5 Environ Pollut (2019) 254:112875. doi: 10.1016/j.envpol.2019.07.043
33. Navab R, Gonzalez-Santos JM, Johnston MR, Liu J, Brodt P, Tsao MS, et al. Expression of Chicken Ovalbumin Upstream Promoter-Transcription Factor II Enhances Invasiveness of Human Lung Carcinoma Cells. Cancer Res (2004) 64:5097–105. doi: 10.1158/0008-5472.CAN-03-1185
36. Zha L, He L, Xie W, Cheng J, Li T, Mohsen MO, et al. Therapeutic Silence of Pleiotrophin by Targeted Delivery of siRNA and Its Effect on the Inhibition of Tumor Growth and Metastasis. PLoS One (2017) 12:e0177964. doi: 10.1371/journal.pone.0177964
37. Miao S, Qiu T, Zhao Y, Wang H, Sun X, Wang Y, et al. Overexpression of S100A13 Protein Is Associated With Tumor Angiogenesis and Poor Survival in Patients With Early-Stage Non-Small Cell Lung Cancer. Thorac Cancer (2018) 9:1136–44. doi: 10.1111/1759-7714.12797
38. Bidkhori G, Narimani Z, Hosseini Ashtiani S, Moeini A, Nowzari-Dalini A, Masoudi-Nejad A. Reconstruction of an Integrated Genome-Scale Co-Expression Network Reveals Key Modules Involved in Lung Adenocarcinoma. PLoS One (2013) 8:e67552. doi: 10.1371/journal.pone.0067552
39. Luo J, Yao Y, Ji S, Sun Q, Xu Y, Liu K, et al. PITX2 Enhances Progression of Lung Adenocarcinoma by Transcriptionally Regulating WNT3A and Activating Wnt/beta-Catenin Signaling Pathway. Cancer Cell Int (2019) 19:96. doi: 10.1186/s12935-019-0800-7
40. Goldstein B, Trivedi M, Speth RC. Alterations in Gene Expression of Components of the Renin-Angiotensin System and Its Related Enzymes in Lung Cancer. Lung Cancer Int (2017) 2017:6914976. doi: 10.1155/2017/6914976
41. Liu Q, Li A, Tian Y, Liu Y, Li T, Zhang C, et al. The Expression Profile and Clinic Significance of the SIX Family in Non-Small Cell Lung Cancer. J Hematol Oncol (2016) 9:119. doi: 10.1186/s13045-016-0339-1
42. Chen Y, Guan JX, Shen H. [MNAT1 Expression in Non-Small Cell Lung Cancer and Its Biological Cellular Impact]. Zhonghua Bing Li Xue Za Zhi (2019) 48:626–32. doi: 10.3760/cma.j.issn.0529-5807.2019.08.008
43. Zeng F, Yu N, Han Y, Ainiwaer J. The Long non-Coding RNA Miat/miR-139-5p/MMP2 Axis Regulates Cell Migration and Invasion in Non-Small-Cell Lung Cancer. J Biosci (2020) 45:51. doi: 10.1007/s12038-020-0019-8
44. De Bastiani MA, Klamt F. Integrated Transcriptomics Reveals Master Regulators of Lung Adenocarcinoma and Novel Repositioning of Drug Candidates. Cancer Med (2019) 8:6717–29. doi: 10.1002/cam4.2493
47. Li W, Yu X, Ma X, Xie L, Xia Z, Liu L, et al. Deguelin Attenuates Non-Small Cell Lung Cancer Cell Metastasis Through Inhibiting the CtsZ/FAK Signaling Pathway. Cell Signal (2018) 50:131–41. doi: 10.1016/j.cellsig.2018.07.001
50. Rauch T, Wang Z, Zhang X, Zhong X, Wu X, Lau SK, et al. Homeobox Gene Methylation in Lung Cancer Studied by Genome-Wide Analysis With a Microarray-Based Methylated CpG Island Recovery Assay. Proc Natl Acad Sci USA (2007) 104:5527–32. doi: 10.1073/pnas.0701059104
51. He J, Zhou X, Li L, Han Z. Long Noncoding Magi2-As3 Suppresses Several Cellular Processes of Lung Squamous Cell Carcinoma Cells by Regulating miR-374a/b-5p/CADM2 Axis. Cancer Manag Res (2020) 12:289–302. doi: 10.2147/CMAR.S232595
52. Li CM, Chen G, Dayton TL, Kim-Kiselak C, Hoersch S, Whittaker CA, et al. Differential Tks5 Isoform Expression Contributes to Metastatic Invasion of Lung Adenocarcinoma. Genes Dev (2013) 27:1557–67. doi: 10.1101/gad.222745.113
53. Khalil A, Dekmak B, Boulos F, Kantrowitz J, Spira A, Fujimoto J, et al. Transcriptomic Alterations in Lung Adenocarcinoma Unveil New Mechanisms Targeted by the TBX2 Subfamily of Tumor Suppressor Genes. Front Oncol (2018) 8:482. doi: 10.3389/fonc.2018.00482
54. Pan J, Fang S, Tian H, Zhou C, Zhao X, Tian H, et al. Lncrna JPX/miR-33a-5p/Twist1 Axis Regulates Tumorigenesis and Metastasis of Lung Cancer by Activating Wnt/beta-Catenin Signaling. Mol Cancer (2020) 19:9. doi: 10.1186/s12943-020-1133-9
56. Kohno T, Otsuka A, Girard L, Sato M, Iwakawa R, Ogiwara H, et al. A Catalog of Genes Homozygously Deleted in Human Lung Cancer and the Candidacy of PTPRD as a Tumor Suppressor Gene. Genes Chromosomes Cancer (2010) 49:342–52. doi: 10.1002/gcc.20746
Keywords: congenital lung malformation, lung tumor, methylation, whole genome, children
Citation: Patrizi S, Pederiva F and d’Adamo AP (2021) Whole-Genome Methylation Study of Congenital Lung Malformations in Children. Front. Oncol. 11:689833. doi: 10.3389/fonc.2021.689833
Received: 01 April 2021; Accepted: 07 June 2021;
Published: 28 June 2021.
Edited by:Hongbo Liu, University of Pennsylvania, United States
Reviewed by:Alexandra Avgustinova, Institute for Research in Biomedicine, Spain
Julia Krushkal, National Cancer Institute (NCI), United States
Copyright © 2021 Patrizi, Pederiva and d’Adamo. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Federica Pederiva, firstname.lastname@example.org
†These authors have contributed equally to this work and share senior authorship