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.

The adult pancreatic duct system accommodates endocrine cells that have the potential to produce insulin. Here we report the characterization and distribution of insulin-immunoreactive cells and endocrine cells within the ductal units of adult human pancreas. Sequential pancreas sections from 12 nondiabetic adults were stained with biomarkers of ductal epithelial cells (cytokeratin 19), acinar cells (amylase), endocrine cells (chromogranin A; neuron-specific enolase), islet hormones (insulin, glucagon, somatostatin, pancreatic polypeptide), cell proliferation (Ki-67), and neogenesis (CD29). The number of islet hormone-immunoreactive cells increased from large ducts to the terminal branches. The insulin-producing cells outnumbered endocrine cells reactive for glucagon, somatostatin, or pancreatic polypeptide. The proportions of insulin-immunoreactive count compared with local islets (100% as a baseline) were 1.5% for the main ducts, 7.2% for interlobular ducts, 24.8% for intralobular ducts, 67.9% for intercalated ducts, and 348.9% for centroacinar cells. Both Ki-67- and CD29-labeled cells were predominantly localized in the terminal branches around the islets. The terminal branches also showed cells coexpressing islet hormones and cytokeratin 19. The adult human pancreatic ducts showed islet hormone-producing cells. The insulin-reactive cells predominantly localized in terminal branches where they may retain potential capability for β-cell neogenesis.

Purpose: 46,XY disorders of sex development (46,XY DSD) are characterized by incomplete masculinization of genitalia with reduced androgenization. Accurate clinical management remains challenging, especially based solely on physical examination. Targeted next-generation sequencing (NGS) with known pathogenic genes provides a powerful tool for diagnosis efficiency. This study aims to identify the prevalent genetic variants by targeted NGS technology and investigate the diagnostic rate in a large cohort of 46,XY DSD patients, with most of them presenting atypical phenotypes. Methods: Two different DSD panels were developed for sequencing purposes, targeting a cohort of 402 patients diagnosed with 46,XY DSD, who were recruited from the Department of Urology at Children’s Hospital, Zhejiang University School of Medicine (Hangzhou, China). The detailed clinical characteristics were evaluated, and peripheral blood was collected for targeted panels to find the patients’ variants. The clinical significance of these variants was annotated according to American College of Medical Genetics and Genomics (ACMG) guidelines. Results: A total of 108 variants across 42 genes were found in 107 patients, including 46 pathogenic or likely pathogenic variants, with 45.7%(21/46) being novel. Among these genes, SRD5A2, AR, FGFR1, LHCGR, NR5A1, CHD7 were the most frequently observed. Besides, we also detected some uncommon causative genes like SOS1, and GNAS. Oligogenic variants were also identified in 9 patients, including several combinations PROKR2/FGFR1/CYP11B1, PROKR2/ATRX, PROKR2/AR, FGFR1/LHCGR/POR, FGFR1/NR5A1, GATA4/NR5A1, WNT4/AR, MAP3K1/FOXL2, WNT4/AR, and SOS1/FOXL2. Conclusion: The overall genetic diagnostic rate was 11.2%(45/402), with an additional 15.4% (62/402) having variants of uncertain significance. Additionally, trio/duo patients had a higher genetic diagnostic rate (13.4%) compared to singletons (8.6%), with a higher proportion of singletons (15.1%) presenting variants of uncertain significance. In conclusion, targeted gene panels identified pathogenic variants in a Chinese 46,XY DSD cohort, expanding the genetic understanding and providing evidence for known pathogenic genes’ involvement.