Integrated bulk and single-cell RNA-sequencing reveals SPOCK2 as a novel biomarker gene in the development of congenital pulmonary airway malformation

Background Congenital pulmonary airway malformation (CPAM) is the most frequent pulmonary developmental malformation and the pathophysiology remains poorly understood. This study aimed to identify the characteristic gene expression patterns and the marker genes essential to CPAM. Methods Tissues from the cystic area displaying CPAM and the area of normal appearance were obtained during surgery. Bulk RNA sequencing (RNA-seq) and single-cell RNA sequencing (scRNA-seq) were performed for integrating analysis. Iterative weighted gene correlation network analysis (iWGCNA) was used to identify specifically expressed genes to CPAM. Results In total, 2074 genes were significantly differentially expressed between the CPAM and control areas. Of these differentially expressed genes (DEGs), 1675 genes were up-regulated and 399 genes were down-regulated. Gene ontology analysis revealed these DEGs were specifically enriched in ciliated epithelium and involved in immune response. We also identified several CPAM-related modules by iWGCNA, among them, P15_I4_M3 module was the most influential module for distinguishing CPAMs from controls. By combining the analysis of the expression dataset from RNA-seq and scRNA-seq, SPOCK2, STX11 , and ZNF331 were highlighted in CPAM. Conclusions Through our analysis of expression datasets from both scRNA-seq and bulk RNA-seq of tissues obtained from patients with CPAM, we identified the characteristic gene expression patterns associated with the condition. Our findings suggest that SPOCK2 could be a potential biomarker gene for the diagnosis and therapeutic target in the development of CPAM, whereas STX11 and ZNF331 might serve as prognostic markers for this condition. Further investigations with larger samples and function studies are necessary to confirm the involvement of these genes in CPAM.