Ran Yang1†
Yong Zhan1†
Yi Li1†
Shu-Yang Dai1
Shi-Wei He1
Chun-Jing Ye1
Ling-Du Meng1
De-Qian Chen1Chen-Bin Dong1
Lian Chen2
Gong Chen1Kui-Ran Dong1
Kai Li1
Shan Zheng1Jun Li1*Wei Yao1*
Rui Dong1*- 1Department of Pediatric Surgery, Children’s Hospital of Fudan University, Shanghai Key Laboratory of Birth Defect, Shanghai, China
- 2Department of Pathology, Children’s Hospital of Fudan University, Shanghai, China
Sialoblastoma (SBL) is an infrequent embryonal malignant tumor originating from the salivary gland, resembling primitive salivary gland anlage, whereas hepatoblastoma (HB) is the most common pediatric liver malignancy. The simultaneous occurrence of both tumors is extremely rare. Here we reported a case of a 6-month-old infant diagnosed with synchronous SBL and HB. The patient received neoadjuvant chemotherapy followed by surgical resection. Fresh tissues of both tumors were collected before and after chemotherapy, which were further profiled by whole exome sequencing (WES) and single-cell RNA sequencing (scRNA-seq). WES analysis revealed potential somatic driver mutation PIK3CA p.Glu454Lys for SBL and canonical mutation CTNNB1 p.Ser45Pro for HB. No shared somatic variants or common copy number alterations were found between SBL and HB primary tumor samples. Though scRNA-seq, single-cell atlases were constructed for both tumors. SBL may recapitulate a pre-acinar stage in the development of salivary gland, including basaloid, duct-like, myoepithelial-like, and cycling phenotypes. In the meantime, HB was composed of tumor cells resembling different stages of the liver, including hepatocyte-like, hepatic progenitor-like, and hepatoblast-like cells. After chemotherapy, both tumors were induced into a more mature phenotype. In terms of transcriptional signatures, SBL and HB showed enhanced expression of epithelial markers KRT8, KRT18, and essential embryo development genes SDC1, MDK, indicating the disruption of normal embryo epithelium development. Finally, heterozygous deleterious germline mutation BLM and FANCI were identified which could predispose the patient to higher cancer risk. It partially explained the reason for the co-occurrence of SBL and HB. Taken together, we provided valuable resources for deciphering cellular heterogeneity and adaptive change of tumor cells after chemotherapy for synchronous SBL and HB, providing insights into the mechanisms leading to synchronous pediatric tumors.
Introduction
Sialoblastoma (SBL) is an extremely rare embryonal epithelial malignancy derived from the salivary gland, which recapitulates the structure of primitive salivary gland. Only about 80 cases have been reported in the literature worldwide. SBL is usually diagnosed at birth or infancy with identifiable cheek mass enlargement. It predominantly occurs in the major salivary gland, with parotid gland more common than submandibular gland (1). As for histological features, the tumor is composed of solid nests or islands of basaloid cells with scant cytoplasm at the periphery and central ductal-like cells. It is separated by fibrous tissue and myoepithelial-like cells (2). Immunohistochemistry (IHC) results showed tumor cells were positive for epithelial markers, e.g., EPCAM, CK7, and CK8. Vimentins and SMA were positive in fibrous stroma with myoepithelial cells. Ki-67 proliferation index varies between patients, with a high index indicating unfavorable outcomes (3, 4). Out of its rarity, few studies have revealed transcriptional features and pathogenesis of SBL. Consequently, there is no agreement on the management of this disease. Surgical resection with negative margins plays a core role, while radiotherapy or neoadjuvant chemotherapies protocols from other solid tumors have also been applied (5).
Unlike the rarity of SBL, hepatoblastoma (HB) is the most common primary liver tumor of childhood. The incidence of HB is 1.5 cases per million each year (6). HB is proposed to be derived from hepatic precursor cells and is morphologically similar to immature fetal hepatocytes. The mutation or deletion of ß-catenin encoding gene (CTNNB1) exon is the most frequently detected gene in HB, resulting in the activation of WNT signally pathways (7). Cairo et al. classified HB as C1 and C2 subtype at the transcriptome level, with C2 presenting a highly proliferative phenotype and signifying poor prognosis (8). This classification was further improved that C2a and C2b subtypes were defined based on VIM expression (9). The current standard of care is comprised of neoadjuvant or adjuvant chemotherapy regimens based on PRETEXT staging system, combined with surgical resection or liver transplantation. Though many studies have investigated genomic and transcriptomic features of HB, few have elucidated intratumoral heterogeneity and tumor evolution after chemotherapy.
Synchronous tumors are defined as cases in which more than one primary tumor in a single patient is identified at initial presentation, which is extremely rare, especially for children. The most common reported pediatric synchronous tumor was Wilms tumor and neuroblastoma. The majority of these patients were also diagnosed with Fanconi Anemia (FA), an autosomal recessive inherited syndrome characterized by congenital defects, aplastic anemia, and a high likelihood for cancer (10, 11). As for simultaneous SBL and HB, only 5 cases have been previously reported without in-depth investigation of potential mechanisms (12–16). Underlying molecular and cellular features of simultaneous SBL and HB were poorly understood.
In the current study, we used whole exome sequencing (WES) and single-cell RNA sequencing (scRNA-seq) to profile the genome and transcriptome of SBL and HB tumor tissues from the same patient before and after chemotherapy (Figure 1A). We have built a cellular atlas and characterized the genomic as well as transcriptomic profiles of synchronous SBL and HB, elucidating how they evolve after treatment. Our study highlighted the intratumoral heterogeneity of both tumors and each subtype displayed distinct signatures reminiscent of different developmental stages of normal cells. Furthermore, we have found heterozygous deleterious germline mutations that might contribute to the co-occurrence of the two embryonal tumors. These results shed light on the cellular and molecular features of synchronous tumors and their potential mechanisms.
Figure 1 The occurrence of two synchronous pediatric tumors. (A) Diagrams of the study design and sample workflow (Created with BioRender.com). Fresh tumor specimens were collected at the timepoint of surgical biopsy and tumor resection. (B) CT images and H&E staining of SBL before and after neoadjuvant chemotherapy. (C) CT images and H&E staining of HB before and after neoadjuvant chemotherapy. CT, Computed tomography; H&E, hematoxylin and eosin; SBL, sialoblastoma; HB, hepatoblastoma; PT, primary tumor; AT, tumor after treatment.
Materials and Methods
Patient and Tumor Samples
This study was approved by the Institutional Review Board of Children’s Hospital of Fudan University. Informed consent was obtained from the legal guardians of the patients. The diagnosis of sialoblastoma and hepatoblastoma was confirmed by two pathologists according to the hematoxylin and eosin (H&E)-stained sections. Samples of tumor tissue were collected during the biopsy and resection surgery.
Whole-Exome Sequencing (WES) DNA Isolation, Library Construction and Sequencing
Tumor tissues collected at diagnosis and surgical resection were stored as fresh frozen tissue. Peripheral blood collected from the patient was used as a matched control for tumor samples. DNA was prepared using the HiPure Universal DNA Kit (Angen Biotech) according to the manufacturer’s instructions. The library was constructed using the Vazyme TruePrep™ DNA Library Prep Kit V2 for Illumina (Vazyme Biotech). The SeqCap EZ Exome Probes v3.0 (Roche) was used to capture exome regions. The libraries were further sequenced on NovaSeq 6000 (Illumina).
SNV and Indel Mutation Identification
Reads were aligned to the human reference genome (GRCh38) using Burrows–Wheeler Aligner (0.7.15). We used Genome Analysis Toolkit (GATK version 4.2.0) to process BAM files, including sorting bam files, marking duplicates, and local realignment around high confidence insertion and deletions. Single nucleotide variants (SNVs) and short insertions/deletions (indels) for each sample were called using GATK HaplotypeCaller (https://gatk.broadinstitute.org/hc/en-us/articles/360035535932-Germline-short-variant-discovery-SNPs-Indels-), followed by joint genotyping using GATK GenotypeGVCFs. Hard-filtering with VariantFiltration was undertaken independently for SNPs and indels. Candidate germline variants were searched in the COSMIC database (the Catalogue Of Somatic Mutations In Cancer; http://cancer.sanger.ac.uk/cosmic) and ClinVar (https://www.ncbi.nlm.nih.gov/clinvar), and their potential effect on protein function was predicted with CADD, SIFT and PolyPhen. Germline variants with CADD PHRED> 20, SIFT deleterious, PolyPhen probably/possibly damaging and gnomAD exome frequency <0.002 were considered as potential pathogenic germline variants. Somatic SNVs and indels in tumor samples were called using the MuTect2 workflow (https://gatk.broadinstitute.org/hc/en-us/articles/360035894731-Somatic-short-variant-discovery-SNVs-Indels-). Variants from germline resources of population data, e.g., gnomAD, 1000 Genomes were also removed. All variants were annotated using ANNOVAR (http://www.openbioinformatics.org/annovar/. Variants with depth were removed from subsequent analyses if the position was <10× in normal and tumor samples. The R package “maftools” (http://bioconductor.org/packages/maftools/) was used for the visualization of somatic mutation data.
Copy Number Variation Analysis
Copy Number Variation (CNV) calling was done using the GATK (v 4.2.0) somatic CNV calling pipeline (https://gatk.broadinstitute.org/hc/en-us/articles/360035531092–How-to-part-I-Sensitively-detect-copy-ratio-alterations-and-allelic-segments#ref2). The peripheral blood sample from the patient was used as a control to identify tumor-specific genomic alterations. Called segments were further annotated with gene-based and region-based information using Annovar.
Validation of Variant Using Sanger Sequencing
Candidate germline variants identified through WES were validated using Sanger sequencing. Whole blood DNA was isolated according to the manufacturer’s protocol with QIAamp DNA Mini Kit (QIAGEN, Cat.51306) and quantified by NnodropOne(Thermo fisher). Premiers were designed as follows.
Sequencing was performed with BigDye™ Terminator v3.1 cycle sequencing kit (Applied Biosystems. cat.4337457) on ABI 3730XL genetic analyzer (Applied Biosystems). Sequencing chromatograms were analyzed by Chromas (http://technelysium.com.au/wp/chromas/).
scRNA-Seq Sample Preparation
Samples were processed within an hour after surgery. Each tissue was cut into 1×1×1 mm fragments and minced on a plate. Then the specimens were enzymatically digested with collagenase IV (Gibco) and DNase I (Sigma) for 30 min at 37°C with agitation. After digestion, specimens were sieved through a 70-μm cell strainer, centrifuged for 5 minutes at 400 g, suspended in Dulbecco’s modified Eagle’s medium (DMEM, Gibco) with 10% fetal bovine serum albumin (FBS), and centrifuged through Lympholyte-H separation (CL5020; Cedarlane) to remove RBCs for 20min, 800 g. Pelleted cells were then re-suspended in DMEM with 10% FBS and assessed for scRNA-seq.
Generation of scRNA-Seq Data
An estimated 8,000 to 10,000 cells were targeted for capture per sample. Then, cell suspensions of each sample were run in the Chromium Controller with appropriate reagents to generate single-cell gel bead-in-emulsions (GEMs) for cell barcoding. The libraries were then pooled and sequenced on NovaSeq 6000 (Illumina) at a depth of approximately 400M reads per sample. Raw sequencing data were converted to FASTQ files with Illumina bcl2fastq, version 2.19.1, and aligned to the human genome reference sequence (GRCH38). The CellRanger (10X Genomics, 3.0.1 version) analysis pipeline was used to sample demultiplexing, barcode processing and single-cell 3′ gene counting to generate a digital gene-cell matrix from this data.
Quality Control and Unsupervised Clustering
Scrublet Package was used to remove cell doublets (17). The Pegasus package (18)was used to process and analyze the gene expression matrix. To filter out low-quality cells, we removed cells with less than 200 genes or more than 6500 genes detected. Cells with over 20% genes derived from the mitochondrial genome were also removed. We excluded nuclear mitochondrial genes and ribosomal genes from the following analysis. In addition, we also removed cells expressing two cell-type markers. In SBL data, one cluster of cells highly expressed T/NK cell markers (NKG7, GNLY, GMZA) and keratins (KRT19, KRT14), which might be T/NK cells contaminated by tumor RNA during the process of sample preparation and library construction. These cells were excluded from the following analysis.
After quality control, the robust genes were identified by pg.identify_ robust_genes function. Then data were log normalized with the function pg.log_norm. Highly variable genes were identified and analyzed by pg.highly_variable_features function. Principal component analysis was performed with the pg.pca function. The batch correction was achieved by pg.run_harmony function. Clusters were identified by pg.neighbors, pg.louvain function and visualized by t-distributed stochastic neighbor embedding (t-SNE) and uniform manifold approximation and projection (UMAP) using pg.tsne and pg.umap function, respectively. DEGs of a given cell type compared with all other cell types were determined with the pg.de_analysis and pg.markers function from the Pegasus package by Mann-Whitney U (MWU) test.
InferCNV Analysis
We identified malignant tumor cells by inferring large-scale chromosomal copy-number variations (CNVs) in each cell (19). We used endothelial cells and fibroblasts in each dataset as reference cells. Since whole exome sequencing showed chr5q loss before and after chemotherapy in SBL, we extracted chr5q loss status and calculated chr5q loss score from ”map_metadata_from_infercnv. txt” in inferCNV results.
Differentially Expressed Genes Analysis
The Deseq2 package (20) was used to perform the differentially expressed genes (DEGs) analysis between primary (PT) and post-chemotherapy (AT) tumor samples. The gene expression matrix was extracted from the Pegasus object. And the DEGs were screened by satisfying p-value <0.05, abs (fold change (FC))>1.5.
Functional Enrichment Analysis
The Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis were performed to investigate the cell function status by using the R package, clusterProfiler (21). FDR p-value < 0.05 was used to distinguish significantly enriched terms. “enrichGO” and “enrichKEGG” are used. Gene set enrichment analysis (GSEA) was performed using MSigDB collections(https://www.gsea-msigdb.org/gsea/msigdb/index.jsp) through the function “GSEA” in clusterProfiler.
Assessing Similarity Between Tumor Cell Subpopulations and Developmental Cell Types
We utilized murine salivary gland datasets GSE150327 as the reference dataset to evaluate the similarity between tumor cell subpopulations and developmental cell types (22). To increase the accuracy of cell classification, we simplified original cell definitions into embryo end bud, embryo duct, embryo basal duct, embryo myoepithelial, postnatal basal duct, postnatal duct, postnatal proacinar, postnatal acinar, and postnatal myoepithelial cells. R package biomaRt (v 2.46.3) was used to convert mouse to human gene names. R package Garnett for Monocle 3 (v 0.2.17) (23)was used to generate cell type classifiers in the maker-free mode with parameters: max_training_samples = 2000, num_unknown = 5. We applied the function classifiy_cells to sialoblastoma malignant cell datasets using the trained classifiers with the following non-default parameters: cluster_extend = TRUE, rank_prob_ratio=1.2, cluster_extend_max_frac_unknown=0.2.
SCENIC Analysis
pySCENIC was applied to run the SCENIC analysis (24). The analysis is composed of pre-processing data, network inference, module generation, motif enrichment, TF-regulon prediction and cellular enrichment. We calculated Regulon specificity score (RSS) for each metaprogram. The top 10 regulons with the highest RSS were used for plotting heatmap of regulon activities. All tumor cells were further clustered by regulon activities.
Determination of Cell Developmental Potential
The R package CytoTRACE (0.3.3) was used for investigating the developmental potential of hepatoblastoma tumor subpopulations. The top 200 genes correlated with CytoTRACE values were chosen as markers for immaturity, respectively. Gene set enrichment analysis was further performed using MSigDB (https://www.gsea-msigdb.org/gsea/msigdb/index.jsp).
Results
Clinical Characteristics of Patient With Synchronous SBL and HB
A 6-month-old infant with synchronous growing left-sided firm cheek mass and abdominal mass was referred to our institution. Computed tomography (CT) scan showed the existence of a soft tissue tumor in the right parotid gland (about 66×61×43 mm3) as well as a heterogeneous lesion affecting the left and right lobes of the liver (about 91×85×85 mm3). Whole-body bone scan showed no systemic metastasis. Lab tests demonstrated elevated levels of blood tumor markers, including alpha-fetoprotein (AFP), carcinoembryonic antigen (CEA) and neuron-specific enolase (NSE) (Table 1). The patient was then subjected to excisional biopsy and diagnosed as SBL with HB based on pathological examination results (Figures 1B, C and Supplementary Figure S1). No other congenital abnormalities were identified. After receiving four cycles of neoadjuvant chemotherapy adopting an intermediate-risk neuroblastoma regimen (carboplatin, cyclophosphamide/etoposide, and doxorubicin), the volume of both tumors showed an obvious decrease (post-treatment size: 60×59×42 mm3 for SBL, 50×30×15 mm3 for HB). Though still higher than normal, the level of blood tumor markers decreased. Segmental liver resection (4b, 5, 6) was performed at the age of 10 months, achieving complete resection of the tumor. Surgical resection of the whole facial mass and partial parotid sialoadenectomy was performed simultaneously. The buccal branch of the facial nerve was infiltrated and thus not preserved. Then the patient received another four courses of chemotherapy. Clinical characteristics and immunohistochemical (IHC) features of SBL and HB at the timepoint of biopsy and surgical resection were summarized in Table 1. No signs of tumor recurrence or metastasis were observed during the 32-month follow-up after surgery.
Table 1 Clinical information and immunohistochemical results.
Uncovering Genomic and Cellular Landscape of SBL
In order to characterize the cellular landscape of SBL and HB, and tumor remodeling after chemotherapy, we collected fresh specimens before and after chemotherapy, encoded as primary tumor (PT) and tumor after treatment (AT), respectively, which were further subjected to scRNA-seq and WES (Figure 1A; Supplementary Tables 1, 2).
For SBL, WES analysis revealed that PT and AT harbored 30 and 24 non-silent coding mutations, respectively (Figure 2A, Supplementary Table 3). No common variants were found between two samples, possibly due to tumor heterogeneity or tumor evolution after chemotherapy. A hotspot gain-of-function mutation PIK3CA E542K was identified in the PT sample, located in the helical domain of PIK3CA (Figure 2C). This variant has been proven to be oncogenic and associated with tumor proliferation and progression in several other tumors, including salivary duct carcinoma, breast and cervical cancer (25, 26). It may also play an essential role in the carcinogenesis of SBL. Mutated genes in PT and AT samples were enriched in similar pathways such as PI3K-Akt and Wnt signaling pathways (Supplementary Figure S2A, B). Although there were no common variants shared between T1 and T2, they showed consistent chromosome 5q34-5q35.3 loss (Figures 2B, D; Supplementary Table S4).
Figure 2 The genomic and transcriptomic landscape of SBL. (A) Oncoplot outlining the distribution of mutations across 4 tumor samples. The color indicated different variant types. (B) Phylogenetic tree diagram showing key somatic mutations as well as shared mutated genes and copy number alterations. (C) Lollipop plot showing the hot spot mutation of PI3KCA p.E542K. Boxes represent functional domains. p85B, p85 binding domain; RBD, Ras binding domain. (D) Scatter plots showing copy number alterations in SBL PT and AT samples. (E) tSNE visualization of 22,107 cells from SBL colored by cell types (left) and density of sample distribution(middle). The right panel shows the proportions of different cell types in SBL PT and AT samples. tSNE, t-distributed stochastic neighbor embedding. (F) Dot plot for expression of cell-type-specific markers in SBL. The color indicated the scaled mean expression of marker genes, and the size indicated the proportion of cells expressing marker genes. SBL, sialoblastoma; HB, hepatoblastoma; PT, primary tumor; AT, tumor after treatment.
After analyzing its genomic landscape, we next turned to explore its transcriptional features at the single-cell level. After removing low-quality cells and potential doublets, 22,107 high-quality cells were proceeded to subsequent analysis. Following normalization of gene expression profiles and graph-based clustering on principal-component analysis (PCA) analysis, t-distributed stochastic neighbor embedding (tSNE) plots were used for visualization (Figure 2E). We identified 5 major cell types based on canonical markers, including epithelial cells (KRT5, KRT14, KRT18, KRT19, S100A2) (27), myeloid cells (CD14, CD68, C1QA, C1QB), fibroblasts (COL1A1, COL3A1), endothelial cells (CLDN5, ENG, PECAM1) and a cluster of sebocytes which highly expressed fatty-acid and sterol synthesis genes such as SCD, SOAT1 and DHCR24, as well as high quantity cytokeratin gene KRT79 (28) (Figure 2F).
In order to identify malignant cells from non-malignant cells, we first calculated the average score of sialoblastoma markers reported in previous studies (KRT5, KRT14, KRT23, TP63, EPCAM, KRT7, KRT19, KRT18, ACTA2, CNN2, S100A2), showing that epithelial cells have the highest tumor score (Supplementary Figure S2C). Next, we utilized inferred copy-number variation (CNV) to differentiate malignant cells from non-malignant ones (19). Consistent with WES results showing tumor cells had chromosome 5q loss before and after chemotherapy, inferred CNV analysis revealed that epithelial cells displayed significantly higher chromosome 5q loss scores than other cells (two-sided Wilcoxon test, pvalue<2.2e-16) (Supplementary Figure S2D, E). Therefore, epithelial cells were identified as malignant cells, constituting the major components of SBL tumor ecosystems.
SBL May Resemble a Pre-Acinar Stage in the Development of Salivary Gland
Following identifying malignant populations, we further explored the transcriptional diversity of malignant cells. Following extracting and re-clustering all SBL tumor cells, we obtained a total of 12 clusters (Figure 3A). Based on gene expression pattern and transcriptional similarity, 12 clusters could be further converged into 4 different programs (Figure 3B; Supplementary Figure S2F, Supplementary Table 5).
Figure 3 The intratumoral heterogeneity within SBL tumor cells. (A) UMAP visualization of 20,111 tumor cells from SBL colored by clusters (left) and density of sample distribution(right). UMAP, uniform manifold approximation and projection. (B) Heatmap displaying the Pearson correlation coefficients calculated between average gene expressions of tumor clusters. (C) Violin plots showing the expression distribution of selected genes involved in cell cycle genes, Wnt signaling pathway, and markers specific to duct, myoepithelial, basal and acinar cells in SBL tumor expression programs. (D) UMAP visualization of 20,111 tumor cells from SBL based on Scenic regulon AUC scores. AUC, area under the ROC Curve. (E) UMAP visualization of specific transcription factor expression and regulon activity. For four transcription factors: (left) histogram of AUC values, together with the chosen cutoff (red dashed line); (middle)the cells with AUC value over the cutoff value are shown in purple, where the regulon is considered active; (right) the expression of the transcription factor is shown. (F) UMAP visualization of murine embryo and postnatal salivary gland datasets curated from GSE150327.Colors indicated different cell types used for Garnett classifier training. (G) Bar plot indicating the percentage of tumor subtypes corresponding to 9 developing cell types. Tumor cells were assigned to different identities predicted by Garnett. (H) Bar plot indicating the percentage of tumor subtypes in SBL PT and AT samples. (I) Volcano plot showing differentially expressed genes between PT (green dots) and RT tumor cells (blue dots). The names of selected important genes are indicated in the plots. (J) Bar chart showing the enrichment of GO terms, based on top50 upregulated genes in PT or AT malignant cells. GO, gene ontology. (K) GSEA enrichment plot of expression signatures of Hallmark_G2M checkpoint in SBL PT malignant cells. (L) GSEA enrichment plot of expression signatures of GO_Morphogenesis of embryonic epithelium in SBL AT malignant cells. SBL, sialoblastoma; PT, primary tumor; AT, tumor after treatment.
P1 was enriched for cell cycle genes (e.g., TOP2A, UBE2C, CKS2), thus defined as cycling tumor cells. P2 was identified as myoepithelial-like cells by expressing muscle contraction and myofibril genes (e.g., ACTA2, MYL, FHL2), resembling myoepithelial cells in the salivary gland. P3 was potentially basaloid cells and expressed significantly higher canonical basal and progenitor cell markers, e.g., KRT5, KRT14, KRT23 (29), suggesting an immature cell state (Figure 3C). Wnt signaling pathway inhibitors (e.g., DKK4, APCDD1, NKD1) and fibroblast growth factors (e.g., FGF4, FGF19, FGF20) were also tightly involved in its signatures. It has been revealed that the FGF signaling pathway could inhibit the Wnt pathway to repress ductal fate and help maintain the undifferentiated state of bud cells (30), which we supposed might be one of the mechanisms leading to the immature state of basaloid cells. P4 (duct-like cells) exhibited higher expression of common epithelial differentiation markers (e.g., KRT6B, KRT7, KRT8, KRT18, KRT19), and lower expression of basal or myoepithelial genes, demonstrating a more mature phenotype. Identification of divergent tumor phenotype was consistent with histology results, which revealed SBL was composed of solid nests of basaloid epithelium with central ductal-like cells and fibroconnective stroma. KI67, SMA, CK5, CK7 specifically marked P1, P2, P3 and P4, respectively. Acinar genes were barely expressed by SBL tumor cells (Supplementary Figure S1A, B).
Since transcriptional programs are determined by underlying gene regulatory network (GRN), we next applied single-cell regulatory network inference and clustering (SCENIC) analysis to interrogate stable GRNs in each tumor subpopulation (24). In line with expectations, this analysis further confirmed a clear distinction in GRNs between different states (Figure 3D). Corresponding to its epithelium morphogenesis phenotype, activated regulons such as PRRX1 and MSX2 in basaloid cells were all associated with epithelial cell programming and differentiation (31, 32). As for myoepithelial-like cells, dominant regulons were transcriptional factors (TFs) that collectively control myoepithelial cell state (e.g., MYLK, SNAI2, TP73), crucial for maintaining mesenchymal cell state and inhibiting transition into mature luminal cells (33–35). As duct-like cells expressed ductal markers, e.g., KRT8, KRT18, TFs that regulated its phenotype equally promotes ductal cell differentiation and tube formation (EHF, TFCP2L1) (36, 37)(Figure 3E). Cellular stress regulators (e.g., MAFF, KLF9) were also enriched in duct-like cells (Supplementary Figure S2G). The TFs for cycling tumor cells mainly regulate the expression of genes required for progression through the cell cycle (e.g., BRCA1, E2F8) (38), consistent with its proliferation phenotype.
SBL is an embryonal tumor transcriptionally similar to primitive salivary gland anlage.We next explored what developmental stage it resembles. We utilized murine embryo and postnatal salivary gland datasets comprising comprehensive cell types, including end bud, basal duct, duct, and acinar cells (22)(Figure 3F; Supplementary Figures S2H, I). Without manually and selectively refining gene signature for each cell type, we used a marker-free method of Garnett to generate an automated classifier of salivary gland (23). Classifiers by Garnett were concordant with manual classifications during self-validation (Supplementary Figure S2J). Using this Garnett model to train SBL tumor cells, we found most tumor cells were classified as duct or basal duct cells (Figure 3G). Almost no tumor cell was identified as acinar cells, in accordance with the scarce expression of acinar genes (e.g., AQP5, NKCC1, AMY1A) (Figure 3C). We further validated this finding on IHC images, which showed expression of basal and ductal genes exclusively without acinar cell differentiation. These characteristics reflect that SBL may recapitulate the pre-acinar stage of the salivary gland (Supplementary Figures S1A, B). Consistent with its phenotype, more myoepithelial-like tumor cells were classified as myoepithelial cells, while more duct-like tumor cells were identified as duct cells.
In order to investigate therapy-induced dynamic change of SBL tumor cells, we first compared the shift of tumor cell composition, which showed a decrease in the proportion of cycling tumor cells, demonstrating that chemotherapy significantly inhibited proliferation of tumor cells (Pearson’s chi-squared test, pvalue< 2.2e-16) (Figure 3H). We also examined gene expression patterns before and after chemotherapy, and identified 2,932 differentially expressed genes (DEGs) (pvalue<0.05, |logFC|>1.5) (Figure 3I, Supplementary Table 6). Functional gene ontology (GO) enrichment analysis indicated that genes with higher expression in SBL PT tumor cells were mainly enriched in cell division (CENPF, TOP2A), cell migration (BST2, NR4A1, KIF2A, TJP1) and negative regulation of cell differentiation (YBX1, EFNB2, DNMT1), supporting that T1 possessed a more malignant and invasive phenotype, which is consistent with gene set enrichment analysis (GSEA) (Figures 3J, K). In contrast, genes up-regulated in RT tumor cells were associated with negative regulation of cell population proliferation (RARRES1, NDRG2, SPINT2), epithelial cell differentiation (KRT7, KRT8, KRT14, KRT15, KRT19), muscle structure (MYLK, TAGLN), leukocyte mediated immunity (GRN, CD74, SLPI) and cell apoptosis (HSPA1B, HSPA1A, LAMP1) (Figures 3J, L), demonstrating that tumor cells turned into a less proliferated state and displaying a potentially more mature structure after treatment, consistent with the shift of tumor subpopulations.
Uncovering Genomic and Cellular Landscape of HB
After revealing genomic and transcriptional features of SBL, we next sought to uncover characteristics of HB. For HB, WES analysis revealed that PT and AT samples both harbored 41 non-silent coding mutations, with four shared mutated genes including CTNNB1, ANP32C, GOLGA8O and LAMB4 (Figures 2A, B; Supplementary Figure S3A, B, Supplementary Table 3). As the most frequent genes, CTNNB1 p.Ser45Pro missense mutation in exon 3 affects phosphorylation sites of casein kinase I epsilon, leading to the downstream activation of Wnt pathway (Figure 4A) (39). Apart from CTNNB1, we’ve identified several cancer-related gene mutations, including PML, HSP90AA1, BAZ1A, ARHGAP5, LHFPL6 in the primary tumor of HB, which may also contribute to the carcinogenesis of HB (40). Mutated genes in HB were enriched in response to DNA damage, epithelium development, Wnt and PI3K-Akt signaling pathways (Supplementary Figures S3C, D).
Figure 4 The genomic and transcriptomic landscape of HB (A) Lollipop plot showing CTNNB1 p.S45P mutation. Boxes represented functional domains. VCL, interaction with VCL; ARM, Armadillo/beta-catenin-like repeats; SCRIB, interaction with SCRIB. (B) Lollipop plot showing ZMIZ2 in-frame deletion. Boxes represented functional domains. Zf-MIZ, Zinc finger MIZ domain. (C) Scatter plots showing copy number alterations in HB PT and AT samples. (D) UMAP visualization of 21,450 cells from HB colored by cell types (left) and density of sample distribution(middle). Right panel showing the percentage of different cell types in HB PT and AT samples. (E) Dot plot for expression of cell-type-specific markers in HB. The color indicated scaled mean expression of marker genes, and the size indicated the proportion of cells expressing marker genes. HB, hepatoblastoma; PT, primary tumor; AT, tumor after treatment.
Though no common somatic variants were found between primary tumors of SBL and HB, in-frame deletion of PAN3 and ZMIZ2 were identified in both tumor samples after treatment (Figures 2A, B). PAN3, a regulatory subunit of the PAN2-PAN3 complex, was involved in the metabolism of messenger RNA (41). As one of PIAS-like proteins, ZMIZ2 has been revealed to interact with β-catenin physically, augmenting β-catenin-mediated transcription and tumor cell growth (42, 43). The disruption of ZMIZ2 in tumor tissues of both SBL and HB may lead to down-regulation of β-catenin downstream targets and inhibition of the Wnt signaling pathway after treatment (Figure 4B).
In terms of somatic copy number alterations, whole chromosome 2 gain and segmental 9p12-q21 loss were shared in both PT and AT samples of HB. Chromosome 1p12-1q21 deletion was exclusively detected in the PT sample while 1p gain was identified in AT sample of HB (Figure 4C; Supplementary Table 4), suggesting tumor evolution after chemotherapy.
Following revealing its genomic features, we next focused on the transcriptional characteristics of HB. After quality control, we obtained 21,450 single-cell transcriptomes from PT and AT samples of HB. Based on the canonical marker gene expression, 3,085 cells were annotated as tumor cells (GPC3, EPCAM, AFP, ALB, DLK1), while 18,365 cells were designed as non-tumor cells (Figures 4D, E) (44, 45). We further confirmed the identity of tumor cells by evaluation of inferred segmental chromosomal alterations. Consistent with WES findings, PT tumor cells showed chromosome 1q loss and 2 gain, while AT tumor cells exhibited chromosome 1q and 2 gain (Supplementary Figure S3E).
To elucidate various expression states of malignant cells within HB, we performed unsupervised clustering of all tumor cells (Figure 5A). A total of 10 subclusters were generated, which can be aggregated into 4 expression programs using hierarchical clustering based on expression profiles (Figures 5B, C; Supplementary Table 7). The 16-gene signature developed by Cairo et al. and the 4-gene signature developed by Hooks et al., which have been used to discriminate mature and invasive HB, was also evaluated for each program (8, 9) (Figures 5D, E).
Figure 5 The intratumoral heterogeneity within HB malignant cells. (A) UMAP visualization of 3,085 tumor cells from HB colored by clusters (left) and density of sample distribution(right). (B) Heatmap displaying the Pearson correlation coefficients calculated between average gene expressions of tumor clusters. (C) Violin plots showing the expression distribution of selected genes involved in cell cycle genes, and markers specific to mature hepatocytes, hepatoblasts, and mesenchymal cells in HB tumor subtypes. (D) Boxplots showing the expression score of HB C1 and C2 group gene sets curated from Cairo et al. in HB tumor subtypes. *P < 0.05, **P < 0.01, ***P < 0.001, ****P < 0.0001, two-sided Wilcoxon test. (E) Violin plots showing the expression of four-gene signature curated from Hooks et al. in HB tumor subtypes. (F) UMAP visualization of CytoTRACE values of HB malignant cells. (G) Boxplots showing CytoTRACE values in HB tumor subtypes. *P < 0.05, **P < 0.01, ***P < 0.001, ****P < 0.0001, two-sided Wilcoxon test. (H) GSEA analysis of genes correlated with high CytoTRACE values (associated with immaturity). (I) Volcano plot showing differentially expressed genes between HB PT (green dots) and RT tumor cells (blue dots). The names of selected important genes are indicated in the plots. (J) Bar chart showing the enrichment of GO terms, based on top50 upregulated genes in PT or AT malignant cells. (K) Bar plot indicating the percentage of tumor subtypes in HB PT and AT samples. HB, hepatoblastoma; PT, primary tumor; AT, tumor after treatment.
P1, which constituted a major fraction of HB tumor cells and had the highest C1 group score, exhibited a more mature hepatocyte-like phenotype, expressing the highest level of ALB as well as liver metabolism genes (e.g., CYP2E1, GSTA1, G6PC) (Figures 5C, D; Supplementary Figure S3F) (46). Therefore, P1 corresponded to liver features at later stages of fetal development and the fetal histology pattern of HB. Unlike P1, P2 expressed higher expression of liver progenitor markers AFP, EPCAM, KRT19, LGR5, and abundant immune-related genes, including MHC class II antigen presentation genes (e.g., HLA-DRA, HLA-DRB1). P3 reflected distinguished cycling cells, enriching cell-cycle genes (e.g., MKI67, TOP2A, HMGB2). P2 and P3 had higher C2 group scores and relatively lower C1 group scores. And P3 may resemble the C2a group with enhanced expression of TOP2A (Figures 5C–E; Supplementary Figure S3F). The features suggested that P2 and P3 obtained more invasive signatures, recapitulating earlier stages of liver development and embryonal histology pattern of HB. Immunohistology analysis also confirmed the existence of P1 and P2/P3 programs, reflecting the co-occurrence of fetal and embryonal patterns (Supplementary Figure S1C).
In contrast with the other three programs, a fourth program, P4, which was only comprised of a minor fraction of malignant cells, presented moderate expression of both C1 and C2 group genes. Nevertheless, P4 up-regulated mesenchymal features (e.g., SPARC, COL4A2) as well as C2b signature gene VIM, hence annotated as a group of mesenchymal-like cells (Figures 5C–E; Supplementary Figure S3F). In the meantime, enhanced expression of early hepatoblast markers (e.g., KDR, CAV1, ID3) was also found in P4, reminiscent of a subpopulation of ID3+ hepatoblasts restricted to the development stages of liver bud (47–49).
We validated the classification of different transcriptional programs by SCENIC regulon analysis. We selected the regulons with the highest regulon specificity score for each program (Supplementary Figures S3G, H). Our analysis identified HNF4A and NR1I3 as specific regulons associated with P1, both of which were well-characterized regulators for core hepatic genes (50, 51). For P2, the most specific regulons included HMGA2 and MYCN. The deregulated MYCN-HMGA2-LIN28B pathway, which contributed to tumor development and progression (52, 53), was activated in P2, further strengthening its invasive phenotype. P3 was associated with cell-cycle regulators, e.g., FOXM1, BRCA1. Consistent with its mesenchymal state, P4 was enriched for tumor epithelial-mesenchymal transition regulons, e.g., SOX18, ZEB1 and NFIB (54–56).
To further define divergent differential status in HB, we used CytoTRACE (57) to predict the developmental potential for each tumor program. In accordance with immature transcriptional phenotype, cells in P2 and P3 are with the highest CytoTRACE score (Figures 5F, G). Gene sets positively correlated with HB immaturity included MYC targets, kinase activity, and stem cells (Figure 5H).
Finally, to gain insights into the transcriptional changes of HB tumor cells after chemotherapy, we performed a DEG analysis between PT and AT malignant cells (Figure 5I; Supplementary Table 8). We observed an enrichment of genes involved in cell adhesion and immune-response pathways (e.g., antigen processing and presentation, cell chemotaxis) in PT which were signatures of P2, whereas the genes up-regulated in AT mainly belonged to normal liver metabolism pathways (e.g., blood microparticle, steroid metabolic process), the signatures for P1 (Figure 5J). This led to the hypothesis that the composition of tumor cells could change after treatment. We next compared the proportion of each tumor group in HB PT and AT. Consistent with DEG analysis, the proportion of P1 was higher in AT than in PT, while the proportion of P2 was higher in PT. These results suggested a more mature phenotype was induced by chemotherapy. (Figure 5K). H&E further confirmed this finding, exhibiting mixed fetal and embryonal patterns in PT were transformed into fetal-dominant histology patterns in AT (Figure 1C; Supplementary Figures S1C, D).
Investigation of Similarities between Synchronous SBL and HB
The salivary gland and liver have different embryonic origins. It’s been suggested salivary gland epithelium arises from the oral ectoderm, whereas liver parenchyma comes from the foregut endoderm (58, 59). Despite different origins, they share functional similarities as accessory organs of the digestive system, playing a major function in digestion. Recent studies suggested that adult salivary gland progenitors can be induced into functional hepatocytes and pancreatic cells both in vitro and in vivo culture conditions (60, 61). These prompted us to investigate the similarities of SBL and HB. By calculating the DEGs of each PT versus non-tumor cells (stromal and immune cells) respectively, we searched for overlapping genes between the signature genes of SBL and HB. This analysis found 67 DEGs were shared between SBL and HB malignant cells (Figures 6A, B). As one of the hallmark traits of the tumor, deregulating cellular metabolism (ATP metabolic, amino acid, and carbohydrate metabolic process) pathways were up-regulated (Figures 6C, D). Apart from these, SBL and HB both showed enhanced expression of epithelial markers (e.g., KRT8, KRT18, PERP, DSP) and essential embryo development genes (e.g., SDC1, MDK, ERBB3) (62–64), demonstrating that disruption of normal embryo epithelial development gave rise to the tumorigenesis of both tumors. In spite of transcriptional similarities, there still existed apparent differences between SBL and HB. SBL uniquely expressed cytokeratins specific to salivary duct and basal cells (e.g., KRT5, KRT15, KRT23, KRT17, KRT19) while HB showed significant expression of liver metabolism genes (e.g., APOE, APOC3, ALB, AHSG), reflecting signatures of original organs (Figure 6B; Supplementary Table 9).
Figure 6 Transcriptional similarities between SBL and HB malignant cells. (A) Venn diagram displaying shared signature genes between SBL and HB malignant cells. (B) Heatmap showing the expression of shared genes and selected top20 unique signature genes in SBL and HB malignant cells. (C) Bar chart showing the enrichment of GO terms, based on shared signature genes between SBL and HB malignant cells. (D) Bar chart showing the enrichment of KEGG terms, based on shared signature genes between SBL and HB malignant cells. (E) IGV screenshot (upper) and DNA sequence chromatogram of Sanger sequencing (bottom) demonstrating BLM c.G2293A, pV765I germline mutation. (F) IGV screenshot (upper) and DNA sequence chromatogram of Sanger sequencing (bottom) demonstrating FANCI c.G2113A, pE96K germline mutation. SBL, sialoblastoma; HB, hepatoblastoma; PT, primary tumor; DEG, differentially expressed gene.
Out of its rarity of synchronous multiple primary tumors in a single patient, especially for children, we proposed that this patient may carry pathogenic germline mutations that predispose him to higher cancer risk. Based on known cancer predisposing genes curated from COSMIC database (40) and other pan-cancer germline genomic studies (65–67), we found a total of 6 damaging or potentially damaging heterozygous germline variants including missense mutation of BLM and FANCI (Table 2). Sanger sequencing was performed for variant validation (Figures 6E, F; Supplementary Figure S4A–D). BLM and FANCI play an essential role in homologous recombination (HR)-mediated repair of DNA interstrand cross-links, the mutation of which resulted in the deficiency of DNA repair and lead to genome instability (68). The biallelic germline mutations of BLM and FANCI resulted in Bloom syndrome and Fanconi anemia (FA), respectively, characterized by congenital disabilities and cancer predisposition (69, 70). Recently, several studies have also revealed the gene-dose effect (haploinsufficiency) of BLM and FA family genes may be associated with cancer predisposition (71–73). BLM haploinsufficiency led to an early onset of murine T cell lymphoma in response to challenge with leukemia virus (72). Carriers of deleterious BLM mutations are at increased risk to develop colorectal cancer and breast cancer (71, 73). Heterozygous mutations in FA-family genes are associated with hereditary breast cancer and familial ovarian cancer (73, 74). Since BLM helicase and FA proteins work concurrently in the DNA repair and maintenance of genome stability, we suppose the patient with predicted deleterious heterozygous germline mutations of two genes may be sensitive to DNA-damaging agents and thus have an increased cancer risk. In addition, we sought to identify potential pathogenic germline CNVs (Supplementary Table 10). Two CNVs (chromosome 8q24.22-q24.3 and 11p11.2-p12 deletion) encompassed known germline cancer-predisposing genes (RECQL4 and EXT2, respectively) (65–67) and were not appeared in the DGV database, indicating both CNVs were probably damaging (Supplementary Figure S4E, F). However, since the infant was an abandoned child, we could not collect his family history. More cases are needed to validate these findings.
Table 2 Potential pathogenic germline variants.
Discussion
In the current study, we presented a rare case of synchronous SBL and HB, and constructed a cellular and molecular landscape of both tumors from a single-cell level, providing a valuable resource for synchronous tumor research.
In the first part, we focused on the genome and transcriptomic features of SBL. SBL is a rare embryonal salivary gland tumor, first reported by Vawter and Tefft in 1966 (75). Only around 80 cases of SBL have been reported around the world, with few studies elucidating its molecular characteristics (76). We identified a hotspot mutation PIK3CA E542K and chromosome 5q34-5q35.3 loss in PT, which might play an important role in tumorigenesis. Previous studies of SBL have identified 47, XX, del(3)(q13), inv(9) (p11q12)c, +?r[5]/46,XX,inv(9) (p11q12)c[41] in a female patient with submandibular SBL (77), and FGFR2 p.Cys382Arg activated somatic mutation in a child with unresectable SBL (78). Our study did not detect recurrent patterns, suggesting potential genetic heterogeneity of this disease. Nevertheless, a gain of function on the FGFR2 protein results in activation of downstream pathways, including PI3K/AKT, strengthening the importance of PI3K signaling pathway in SBL (79). Through scRNA-seq, we uncovered intra-tumoral heterogeneity of SBL. Tumor cells can be divided into 4 different subtypes, including basaloid cells (KRT5, KRT14), myoepithelial cells (ACTA2, MYL), cycling cells (MKI67, TOP2A) and duct-like cells (KRT7, KRT19). After chemotherapy, tumor cells were induced into a more mature phenotype, indicating the effectiveness of the current regimen. By comparing SBL tumor cells with its normal counterparts, we found SBL may resemble a pre-acinar stage of the salivary gland, while few cells are classified as acinar cells. Immunostaining results also confirmed this finding.
In the second part, we turned to explore the characteristics of HB. Canonical CTNNB1 p.Ser45Pro missense mutation, as well as several other cancer-related genes (e.g., PML, HSP90AA1, BAZ1A, ARHGAP5, LHFPL6), were detected in the PT sample. No common somatic variant was found between SBL and HB PT samples. Previous transcriptomic analysis identified 3 major subtypes of HB samples, including C1, C2a, and C2b, revealing inter-tumoral heterogeneity(9). Our results supported this notion and expanded it to intra-tumoral heterogeneity. Our analysis indicated four subtypes of HB tumor cells, with P1, P3, P4 corresponding to C1, C2a, and C2b, respectively. P1 represented fetal-hepatocyte-like tumor cells in this tumor sample, while P2/P3 represented embryonal-hepatocyte-like tumor cells.
Finally, we investigated the similarities of both tumors and detected hazardous germline mutations of this patient. SBL and HB both showed enhanced expression of epithelial markers (e.g., KRT8, KRT18 and embryo development genes (e.g., SDC1, MDK), suggesting the interruption of normal epithelial development. In addition, we found 6 damaging or potentially damaging heterozygous germline variants. Deleterious germline mutation BLM and FANCI play an essential role in DNA repair. Heterozygous mutation of both genes may predispose the patient to higher cancer risk, contributing to the tumorigenesis of SBL and HB. However, despite several studies suggesting that heterozygous mutations in DNA repair genes may predispose carriers to increased cancer susceptibility (71–73), this concept still remains controversial and association between genes and specific cancer type warranted further validation. Though predicted to be deleterious and pathogenic by in-silico algorithms, the germline variants of BLM and FANCI identified in this study were missense mutations with unclear impact on protein function, which requires further validation. In addition, we identified potentially damaging chr8q24.22-q24.3 and chr11p11.2-p12 deletion affecting tumor suppressor gene RECQL4 and EXT2, respectively (80, 81). In 1999, Choi HJ et al. reported a case of chr1q deletion with synchronous SBL and HB. We’ve also detected germline chr1q12-q21.2 deletion (Supplementary Table 10) and somatic chr1p12-1q21 deletion in the primary HB tumor tissue, indicating that chr1q deletion may participate in the pathogenesis of this disease. Rigorous functional experiments and more cases are needed to verify our findings.
As tumors from different organs, the synchronous occurrence of SBL and HB has rarely been revealed, with only 5 previously reported cases (Table 3) (12–16). After summarizing all cases, we found that almost all patients were prenatally diagnosed with well-differentiated histology, further revealing the embryonal origin of both tumors and providing clues into the development similarities of salivary gland and liver. The salivary gland and liver are accessory organs of the digestive system, sharing functional similarities in food digestion. The development of both organs during embryo development consisted of primitive bud formation followed by branching morphogenesis. Interaction of epithelium and mesenchyme plays a vital role in the maturation of both organs. More interestingly, epithelial cells from the salivary gland can be induced into functional hepatocytes in vitro, which can alleviate severe acute liver damage caused by carbon tetrachloride in SCID mouse (82). Several studies have proposed that both organs might share the same embryonal origin. However, recent genetic tracing results indicated the ectodermal origin of the salivary gland and the endodermal origin of the liver (58, 59). More studies are needed to explain the potential developmental similarities of both organs. Though with limited cases, there might be a slight male preponderance for synchronous SBL and HB, with four male patients and only one female patient. This could be partially explained by male predisposition in both SBL (Male : Female=~1.3:1) (1, 4, 83, 84) and HB (Male : Female=1.4-2:1) (85–87). In addition, two reported cases adopted neoadjuvant chemotherapy (cisplatin, vincristine, 5-fluorouracil or vincristine, actinomycin, cyclophosphamide) before surgical resection. Though the volume of HB decreased, SBL showed no response to treatment with a static size. In our case, we utilized an intermediate neuroblastoma regimen (carboplatin, cyclophosphamide/etoposide, and doxorubicin), which were effective in decreasing tumor volume and inducing a more mature phenotype for both tumors. This therapy may provide a reference for treating synchronous SBL and HB.
Table 3 Summary of synchronous pediatric sialoblastoma and hepatoblastoma cases.
Our study also has its limitations. First, out of the rarity of synchronous SBL and HB, our analysis was confined to only one case and did not have enough sample size. Results may not be comprehensive and offer abundant information. Much remains to be gained by enlarging the cohort. Second, because of the specificity of the case, we were not able to obtain a family history of the patient and further investigate the inheritance pattern of germline mutations. Also due to limitations of WES technology (e.g., non-uniform read-depth distribution and restricted to exonic regions), it’s difficult to accurately detect large CNVs or structural variants (SVs) in genomic areas. Therefore, identified potential deleterious CNVs could be confirmed with orthogonal methods. The possible role of pathogenic SVs should be further explored if whole genome sequencing data was available. Future efforts in collecting more cases and entire family history will facilitate the identification of genetic abnormalities of the disease.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: Gene Expression Omnibus (GEO), accession number GSE166345, and Genome Sequence Archive (GSA) for Human in National Genomics Data Center (NGDC), accession number HRA000610.
Ethics Statement
The studies involving human participants were reviewed and approved by Children’s Hospital of Fudan University. Written informed consent to participate in this study was provided by the participants’ legal guardian/next of kin.
Author Contributions
RY, YZ, and YL contributed to conception and design, data acquisition, analysis and interpretation, drafted and critically revised the manuscript. C-JY, L-DM, and D-QC contributed to data curation and investigation. S-WH and S-YD contributed to materials collection and literature review. C-BD, LC, GC, K-RD, and SZ contributed to manuscript review and editing. JL, WY, and RD contributed to the experiment design and critically revised the manuscript. All authors contributed to the article and approved the submitted version.
Funding
This work was supported by the Cyrus Tang Foundation (No. ZSBK0070), Shanghai Municipal Key Clinical Specialty (No. shslczdzk05703), National Natural Science Foundation of China (No. 82072782) and Shanghai Municipal Health Commission (No. EK112520180301).
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.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fonc.2022.893206/full#supplementary-material
References
1. Di Micco R, Prufer F, Bruder E, Schifferli A, Gurtler N. Sialoblastoma of the Submandibular Gland: A Distinct Entity? Eur J Pediatr (2019) 178(8):1301–4. doi: 10.1007/s00431-019-03411-x
2. Worasakwutiphong S. Sialoblastoma of the Cheek: A Case Report and Review of the Literature. Mol Clin Oncol (2016) 4:925. doi: 10.3892/mco.2016.1060
3. Saffari Y, Blei F, Warren SM, Milla S, Greco MA. Congenital Minor Salivary Gland Sialoblastoma: A Case Report and Review of the Literature. Fetal Pediatr Pathol (2011) 30(1):32–9. doi: 10.3109/15513815.2010.502961
4. Irace AL, Adil EA, Archer NM, Silvera VM, Perez-Atayde A, Rahbar R. Pediatric Sialoblastoma: Evaluation and Management. Int J Pediatr Otorhinolaryngol (2016) 87:44–9. doi: 10.1016/j.ijporl.2016.04.037
5. Saravakos P, Hartwein J, Fayyazi A. Sialoblastoma of the Parotid Gland in a 13-Year-Old Girl With Multiple Recurrences and Long-Term Follow-Up. Head Neck-J Sci Specialties Head Neck (2016) 38(1):E13–E5. doi: 10.1002/hed.24084
6. Ramos-Gonzalez G, LaQuaglia M, O'Neill AF, Elisofon S, Zurakowski D, Kim HB, et al. Long-Term Outcomes of Liver Transplantation for Hepatoblastoma: A Single-Center 14-Year Experience. Pediatric Transplantation (2018) 22(6):e13250. doi: 10.1111/petr.13250
7. Koch A, Denkhaus D, Albrecht S, Leuschner I, von Schweinitz D, Pietsch T. Childhood Hepatoblastomas Frequently Carry a Mutated Degradation Targeting Box of the Beta-Catenin Gene. Cancer Res (1999) 59(2):269–73.
8. Cairo S, Armengol C, De Reynies A, Wei Y, Thomas E, Renard CA, et al. Hepatic Stem-Like Phenotype and Interplay of Wnt/Beta-Catenin and Myc Signaling in Aggressive Childhood Liver Cancer. Cancer Cell (2008) 14(6):471–84. doi: 10.1016/j.ccr.2008.11.002
9. Hooks KB, Audoux J, Fazli H, Lesjean S, Ernault T, Dugot-Senant N, et al. New Insights Into Diagnosis and Therapeutic Options for Proliferative Hepatoblastoma. Hepatology (2018) 68(1):89–102. doi: 10.1002/hep.29672
10. Compostella A, Toffolutti T, Soloni P, Dall'Igna P, Carli M, Bisogno G. Multiple Synchronous Tumors in a Child With Fanconi Anemia. J Pediatr Surg (2010) 45(2):e5–8. doi: 10.1016/j.jpedsurg.2009.11.015
11. Serra A, Eirich K, Winkler AK, Mrasek K, Göhring G, Barbi G, et al. Shared Copy Number Variation in Simultaneous Nephroblastoma and Neuroblastoma Due to Fanconi Anemia. Mol Syndromol (2012) 3(3):120–30. doi: 10.1159/000341935
12. Huh CY, Choi HJ, Kim SB, Lee S, Lim SJ, Yang MH. Case Reports : 1 Case of Chromosome 1q Deletion With Sialoblastoma and Hepatoblastoma in Neonate. Obstet Gynecol Sci (1999) 42(1):175–8.
13. Siddiqi SH, Solomon MP, Haller JO. Sialoblastoma and Hepatoblastoma in a Neonate. Pediatr Radiol (2000) 30(5):349–51. doi: 10.1007/s002470050758
14. Stones DK, Jansen JC, Griessel D. Sialoblastoma and Hepatoblastoma in a Newborn Infant. Pediatr Blood Cancer (2009) 52(7):883–5. doi: 10.1002/pbc.21945
15. Cheng YK, Chu WC, Law LW, Ting YH, Leung TY. A Fetus With a Huge Neck Mass and a Large Abdominal Circumference–a Rare Case of Sialoblastoma and Hepatoblastoma. Prenat Diagn (2012) 32(9):915–7. doi: 10.1002/pd.3927
16. Rodríguez-Zubieta M, Albarenque K, Lagues C, San Roman A, Varela M, Russo D, et al. Two Synchronous Neonatal Tumors: An Extremely Rare Case. Case Rep Pathol (2021) 2021:6674372. doi: 10.1155/2021/6674372
17. Wolock SL, Lopez R, Klein AM. Computational Identification of Cell Doublets in Single-Cell Transcriptomic Data. Cell Systems (2019) 8(4):281–91.e9. doi: 10.1016/j.cels.2018.11.005
18. Li B, Gould J, Yang Y, Sarkizova S, Tabaka M, Ashenberg O, et al. Cumulus Provides Cloud-Based Data Analysis for Large-Scale Single-Cell and Single-Nucleus Rna-Seq. Nat Methods (2020) 17(8):793–8. doi: 10.1038/s41592-020-0905-x
19. Venteicher AS, Tirosh I, Hebert C, Yizhak K, Neftel C, Filbin MG, et al. Decoupling Genetics, Lineages, and Microenvironment in Idh-Mutant Gliomas by Single-Cell Rna-Seq. Science (2017) 355(6332):eaai8478. doi: 10.1126/science.aai8478
20. Love MI, Huber W, Anders S. Moderated Estimation of Fold Change and Dispersion for Rna-Seq Data With Deseq2. Genome Biol (2014) 15(12):550. doi: 10.1186/s13059-014-0550-8
21. Yu G, Wang LG, Han Y, He QY. Clusterprofiler: An R Package for Comparing Biological Themes Among Gene Clusters. OMICS (2012) 16(5):284–7. doi: 10.1089/omi.2011.0118
22. Hauser BR, Aure MH, Kelly MC, Hoffman MP, Chibly AM. Generation of a Single-Cell Rnaseq Atlas of Murine Salivary Gland Development. iScience (2020) 23(12):101838. doi: 10.1016/j.isci.2020.101838
23. Pliner HA, Shendure J, Trapnell C. Supervised Classification Enables Rapid Annotation of Cell Atlases. Nat Methods (2019) 16(10):983–6. doi: 10.1038/s41592-019-0535-3
24. Van de Sande B, Flerin C, Davie K, De Waegeneer M, Hulselmans G, Aibar S, et al. A Scalable Scenic Workflow for Single-Cell Gene Regulatory Network Analysis. Nat Protoc (2020) 15(7):2247–76. doi: 10.1038/s41596-020-0336-2
25. Griffith CC, Seethala RR, Luvison A, Miller M, Chiosea SI. Pik3ca Mutations and Pten Loss in Salivary Duct Carcinomas. Am J Surg Pathol (2013) 37(8):1201–7. doi: 10.1097/PAS.0b013e3182880d5a
26. Jiang W, He T, Liu S, Zheng Y, Xiang L, Pei X, et al. The Pik3ca E542k and E545k Mutations Promote Glycolysis and Proliferation Via Induction of the B-Catenin/Sirt3 Signaling Pathway in Cervical Cancer. J Hematol Oncol (2018) 11(1):139. doi: 10.1186/s13045-018-0674-5
27. Williams SB, Ellis GL, Warnock GR. Sialoblastoma: A Clinicopathologic and Immunohistochemical Study of 7 Cases. Ann Diagn Pathol (2006) 10(6):320–6. doi: 10.1016/j.anndiagpath.2006.02.008
28. Veniaminova NA, Grachtchouk M, Doane OJ, Peterson JK, Quigley DA, Lull MV, et al. Niche-Specific Factors Dynamically Regulate Sebaceous Gland Stem Cells in the Skin. Dev Cell (2019) 51(3):326–40.e4. doi: 10.1016/j.devcel.2019.08.015
29. Su X, Long Q, Bo J, Shi Y, Zhao LN, Lin Y, et al. Mutational and Transcriptomic Landscapes of a Rare Human Prostate Basal Cell Carcinoma. Prostate (2020) 80(6):508–17. doi: 10.1002/pros.23965
30. Patel N, Sharpe PT, Miletich I. Coordination of Epithelial Branching and Salivary Gland Lumen Formation by Wnt and Fgf Signals. Dev Biol (2011) 358(1):156–67. doi: 10.1016/j.ydbio.2011.07.023
31. Zhang WL, Wang SS, Jiang YP, Liu Y, Yu XH, Wu JB, et al. Fatty Acid Synthase Contributes to Epithelial-Mesenchymal Transition and Invasion of Salivary Adenoid Cystic Carcinoma Through Prrx1/Wnt/B-Catenin Pathway. J Cell Mol Med (2020) 24(19):11465–76. doi: 10.1111/jcmm.15760
32. Jaskoll T, Luo W, Snead ML. Msx-2 Expression and Glucocorticoid-Induced Overexpression in Embryonic Mouse Submandibular Glands. J Craniofac Genet Dev Biol (1998) 18(2):79–87.
33. He WQ, Wang J, Sheng JY, Zha JM, Graham WV, Turner JR. Contributions of Myosin Light Chain Kinase to Regulation of Epithelial Paracellular Permeability and Mucosal Homeostasis. Int J Mol Sci (2020) 21(3):993. doi: 10.3390/ijms21030993
34. Idoux-Gillet Y, Nassour M, Lakis E, Bonini F, Theillet C, Du Manoir S, et al. Slug/Pcad Pathway Controls Epithelial Cell Dynamics in Mammary Gland and Breast Carcinoma. Oncogene (2018) 37(5):578–88. doi: 10.1038/onc.2017.355
35. Weber A, Langhanki L, Schütz A, Gerstner A, Bootz F, Wittekind C, et al. Expression Profiles of P53, P63, and P73 in Benign Salivary Gland Tumors. Virch Arch (2002) 441(5):428–36. doi: 10.1007/s00428-002-0705-y
36. Werth M, Schmidt-Ott KM, Leete T, Qiu A, Hinze C, Viltard M, et al. Transcription Factor Tfcp2l1 Patterns Cells in the Mouse Kidney Collecting Ducts. Elife (2017) 6:e24265. doi: 10.7554/eLife.24265
37. Tanaka J, Ogawa M, Hojo H, Kawashima Y, Mabuchi Y, Hata K, et al. Generation of Orthotopically Functional Salivary Gland From Embryonic Stem Cells. Nat Commun (2018) 9(1):4216. doi: 10.1038/s41467-018-06469-7
38. Liao GB, Li XZ, Zeng S, Liu C, Yang SM, Yang L, et al. Regulation of the Master Regulator Foxm1 in Cancer. Cell Commun Signal (2018) 16(1):57. doi: 10.1186/s12964-018-0266-6
39. Sakanaka C. Phosphorylation and Regulation of Beta-Catenin by Casein Kinase I Epsilon. J Biochem (2002) 132(5):697–703. doi: 10.1093/oxfordjournals.jbchem.a003276
40. 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(11):696–705. doi: 10.1038/s41568-018-0060-1
41. Wolf J, Passmore LA. Mrna Deadenylation by Pan2-Pan3. Biochem Soc Trans (2014) 42(1):184–7. doi: 10.1042/bst20130211
42. Zhu Y, Gu L, Lin X, Cui K, Liu C, Lu B, et al. Linc00265 Promotes Colorectal Tumorigenesis Via Zmiz2 and Usp7-Mediated Stabilization of B-Catenin. Cell Death Differ (2020) 27(4):1316–27. doi: 10.1038/s41418-019-0417-3
43. Lee SH, Zhu C, Peng Y, Johnson DT, Lehmann L, Sun Z. Identification of a Novel Role of Zmiz2 Protein in Regulating the Activity of the Wnt/B-Catenin Signaling Pathway. J Biol Chem (2013) 288(50):35913–24. doi: 10.1074/jbc.M113.529727
44. Bondoc A, Glaser K, Jin K, Lake C, Cairo S, Geller J, et al. Identification of Distinct Tumor Cell Populations and Key Genetic Mechanisms Through Single Cell Sequencing in Hepatoblastoma. Commun Biol (2021) 4(1):1049. doi: 10.1038/s42003-021-02562-8
45. Nagae G, Yamamoto S, Fujita M, Fujita T, Nonaka A, Umeda T, et al. Genetic and Epigenetic Basis of Hepatoblastoma Diversity. Nat Commun (2021) 12(1):5423. doi: 10.1038/s41467-021-25430-9
46. Xu J, Ma HY, Liang S, Sun M, Karin G, Koyama Y, et al. The Role of Human Cytochrome P450 2e1 in Liver Inflammation and Fibrosis. Hepatol Commun (2017) 1(10):1043–57. doi: 10.1002/hep4.1115
47. Goldman O, Han S, Sourisseau M, Dziedzic N, Hamou W, Corneo B, et al. Kdr Identifies a Conserved Human and Murine Hepatic Progenitor and Instructs Early Liver Development. Cell Stem Cell (2013) 12(6):748–60. doi: 10.1016/j.stem.2013.04.026
48. Wang X, Yang L, Wang YC, Xu ZR, Feng Y, Zhang J, et al. Comparative Analysis of Cell Lineage Differentiation During Hepatogenesis in Humans and Mice at the Single-Cell Transcriptome Level. Cell Res (2020) 30(12):1109–26. doi: 10.1038/s41422-020-0378-6
49. Tse EY, Ko FC, Tung EK, Chan LK, Lee TK, Ngan ES, et al. Caveolin-1 Overexpression Is Associated With Hepatocellular Carcinoma Tumourigenesis and Metastasis. J Pathol (2012) 226(4):645–53. doi: 10.1002/path.3957
50. Huang KW, Reebye V, Czysz K, Ciriello S, Dorman S, Reccia I, et al. Liver Activation of Hepatocellular Nuclear Factor-4α by Small Activating Rna Rescues Dyslipidemia and Improves Metabolic Profile. Mol Ther Nucleic Acids (2020) 19:361–70. doi: 10.1016/j.omtn.2019.10.044
51. Yamazaki Y, Moore R, Negishi M. Nuclear Receptor Car (Nr1i3) Is Essential for Ddc-Induced Liver Injury and Oval Cell Proliferation in Mouse Liver. Lab Invest (2011) 91(11):1624–33. doi: 10.1038/labinvest.2011.115
52. Eide HA, Halvorsen AR, Bjaanaes MM, Piri H, Holm R, Solberg S, et al. The Mycn-Hmga2-Cdkn2a Pathway in Non-Small Cell Lung Carcinoma–Differences in Histological Subtypes. BMC Cancer (2016) 16:71. doi: 10.1186/s12885-016-2104-9
53. Morishita A, Zaidi MR, Mitoro A, Sankarasharma D, Szabolcs M, Okada Y, et al. Hmga2 Is a Driver of Tumor Metastasis. Cancer Res (2013) 73(14):4289–99. doi: 10.1158/0008-5472.Can-12-3848
54. Huaqi Y, Caipeng Q, Qiang W, Yiqing D, Xiang D, Xu T, et al. Transcription Factor Sox18 Promotes Clear Cell Renal Cell Carcinoma Progression and Alleviates Cabozantinib-Mediated Inhibitory Effects. Mol Cancer Ther (2019) 18(12):2433–45. doi: 10.1158/1535-7163.Mct-19-0043
55. Liu Z, Chen J, Yuan W, Ruan H, Shu Y, Ji J, et al. Nuclear Factor I/B Promotes Colorectal Cancer Cell Proliferation, Epithelial-Mesenchymal Transition and 5-Fluorouracil Resistance. Cancer Sci (2019) 110(1):86–98. doi: 10.1111/cas.13833
56. Zhang P, Sun Y, Ma L. Zeb1: At the Crossroads of Epithelial-Mesenchymal Transition, Metastasis and Therapy Resistance. Cell Cycle (2015) 14(4):481–7. doi: 10.1080/15384101.2015.1006048
57. Gulati GS, Sikandar SS, Wesche DJ, Manjunath A, Bharadwaj A, Berger MJ, et al. Single-Cell Transcriptional Diversity Is a Hallmark of Developmental Potential. Science (2021) 367(6476):4057–11. doi: 10.1126/science.aax0249
58. Rothova M, Thompson H, Lickert H, Tucker AS. Lineage Tracing of the Endoderm During Oral Development. Dev Dyn (2012) 241(7):1183–91. doi: 10.1002/dvdy.23804
59. Ober EA, Lemaigre FP. Development of the Liver: Insights Into Organ and Tissue Morphogenesis. J Hepatol (2018) 68(5):1049–62. doi: 10.1016/j.jhep.2018.01.005
60. Petrakova OS, Terskikh VV, Chernioglo ES, Ashapkin VV, Bragin EY, Shtratnikova VY, et al. Comparative Analysis Reveals Similarities between Cultured Submandibular Salivary Gland Cells and Liver Progenitor Cells. Springerplus (2014) 3:183. doi: 10.1186/2193-1801-3-183
61. Hisatomi Y, Okumura K, Nakamura K, Matsumoto S, Satoh A, Nagano K, et al. Flow Cytometric Isolation of Endodermal Progenitors From Mouse Salivary Gland Differentiate Into Hepatic and Pancreatic Lineages. Hepatology (2004) 39(3):667–75. doi: 10.1002/hep.20063
62. Garcha D, Walker SP, MacDonald TM, Hyett J, Jellins J, Myers J, et al. Circulating Syndecan-1 Is Reduced in Pregnancies With Poor Fetal Growth and Its Secretion Regulated by Matrix Metalloproteinases and the Mitochondria. Sci Rep (2021) 11(1):16595. doi: 10.1038/s41598-021-96077-1
63. Espinosa-Medina I, Jevans B, Boismoreau F, Chettouh Z, Enomoto H, Müller T, et al. Dual Origin of Enteric Neurons in Vagal Schwann Cell Precursors and the Sympathetic Neural Crest. Proc Natl Acad Sci USA (2017) 114(45):11980–5. doi: 10.1073/pnas.1710308114
64. Filippou PS, Karagiannis GS, Constantinidou A. Midkine (Mdk) Growth Factor: A Key Player in Cancer Progression and a Promising Therapeutic Target. Oncogene (2020) 39(10):2040–54. doi: 10.1038/s41388-019-1124-8
65. Akhavanfard S, Padmanabhan R, Yehia L, Cheng F, Eng C. Comprehensive Germline Genomic Profiles of Children, Adolescents and Young Adults With Solid Tumors. Nat Commun (2020) 11(1):2206. doi: 10.1038/s41467-020-16067-1
66. Fiala EM, Jayakumaran G, Mauguen A, Kennedy JA, Bouvier N, Kemel Y, et al. Prospective Pan-Cancer Germline Testing Using Msk-Impact Informs Clinical Translation in 751 Patients With Pediatric Solid Tumors. Nat Cancer (2021) 2:357–65. doi: 10.1038/s43018-021-00172-1
67. Toss A, Quarello P, Mascarin M, Banna GL, Zecca M, Cinieri S, et al. Cancer Predisposition Genes in Adolescents and Young Adults (Ayas): A Review Paper From the Italian Aya Working Group. Curr Oncol Rep (2022) 24(7):843–60. doi: 10.1007/s11912-022-01213-3
68. Niraj J, Farkkila A, D'Andrea AD. The Fanconi Anemia Pathway in Cancer. Annu Rev Cancer Biol (2019) 3:457–78. doi: 10.1146/annurev-cancerbio-030617-050422
69. Ababou M. Bloom Syndrome and the Underlying Causes of Genetic Instability. Mol Genet Metab (2021) 133(1):35–48. doi: 10.1016/j.ymgme.2021.03.003
70. Nepal M, Che R, Zhang J, Ma C, Fei P. Fanconi Anemia Signaling and Cancer. Trends Cancer (2017) 3(12):840–56. doi: 10.1016/j.trecan.2017.10.005
71. de Voer RM, Hahn MM, Mensenkamp AR, Hoischen A, Gilissen C, Henkes A, et al. Deleterious Germline Blm Mutations and the Risk for Early-Onset Colorectal Cancer. Sci Rep (2015) 5:14060. doi: 10.1038/srep14060
72. Goss KH, Risinger MA, Kordich JJ, Sanz MM, Straughen JE, Slovek LE, et al. Enhanced Tumor Formation in Mice Heterozygous for Blm Mutation. Science (2002) 297(5589):2051–3. doi: 10.1126/science.1074340
73. Thompson ER, Doyle MA, Ryland GL, Rowley SM, Choong DY, Tothill RW, et al. Exome Sequencing Identifies Rare Deleterious Mutations in DNA Repair Genes Fancc and Blm as Potential Breast Cancer Susceptibility Alleles. PloS Genet (2012) 8(9):e1002894. doi: 10.1371/journal.pgen.1002894
74. Gordiev M, Brovkina O, Shigapova L, Shagimardanova E, Enikeev R, Nikitin A, et al. Heterozygous Mutation in Fanconi Anemia Genes Associated With Hereditary Breast Cancer. Ann Oncol (2019) 30:iii10. doi: 10.1093/annonc/mdz095.028
76. Choudhary K, Panda S, Beena VT, Rajeev R, Sivakumar R, Krishanan S. Sialoblastoma: A Literature Review From 1966-2011. Natl J Maxillofac Surg (2013) 4(1):13–8. doi: 10.4103/0975-5950.117821
77. Mertens F, Wahlberg P, Domanski HA. Clonal Chromosome Aberrations in a Sialoblastoma. Cancer Genet Cytogenet (2009) 189(1):68–9. doi: 10.1016/j.cancergencyto.2008.10.005
78. Harris MH, DuBois SG, Glade Bender JL, Kim A, Crompton BD, Parker E, et al. Multicenter Feasibility Study of Tumor Molecular Profiling to Inform Therapeutic Decisions in Advanced Pediatric Solid Tumors: The Individualized Cancer Therapy (Icat) Study. JAMA Oncol (2016) 2(5):608–15. doi: 10.1001/jamaoncol.2015.5689
79. Xie Y, Su N, Yang J, Tan Q, Huang S, Jin M, et al. Fgf/Fgfr Signaling in Health and Disease. Signal Transduction Targeted Ther (2020) 5(1):181. doi: 10.1038/s41392-020-00222-7
80. Croteau DL, Singh DK, Hoh Ferrarelli L, Lu H, Bohr VA. Recql4 in Genomic Instability and Aging. Trends Genet (2012) 28(12):624–31. doi: 10.1016/j.tig.2012.08.003
81. Heddar A, Fermey P, Coutant S, Angot E, Sabourin JC, Michelin P, et al. Familial Solitary Chondrosarcoma Resulting From Germline Ext2 Mutation. Genes Chromosomes Cancer (2017) 56(2):128–34. doi: 10.1002/gcc.22419
82. Zhang C, Li Y, Zhang XY, Liu L, Tong HZ, Han TL, et al. Therapeutic Potential of Human Minor Salivary Gland Epithelial Progenitor Cells in Liver Regeneration. Sci Rep (2017) 7(1):12707. doi: 10.1038/s41598-017-11880-z
83. Wang Q, Ma S, Chen H, Yang M, Cai W. Sialoblastoma in Chin and Management of Treatment. Int J Pediatr Otorhinolaryngol (2018) 109:168–73. doi: 10.1016/j.ijporl.2018.03.023
84. Ronchi A, Lucà S, Errico ME, D'Onofrio V, Mariani P, Laino L, et al. Congenital Sialoblastoma in a Newborn: Diagnostic Challenge of a Rare Entity. Minerva Dent Oral Sci (2021). doi: 10.23736/s2724-6329.21.04440-x
85. von Schweinitz D. Management of Liver Tumors in Childhood. Semin Pediatr Surg (2006) 15(1):17–24. doi: 10.1053/j.sempedsurg.2005.11.004
86. Devi LP, Kumar R, Handique A, Kumar M. Hepatoblastoma—a Rare Liver Tumor With Review of Literature. J Gastrointestinal Cancer (2014) 45(1):261–4. doi: 10.1007/s12029-014-9659-y
Keywords: hepatoblastoma, sialoblastoma, single cell analysis, embryonal neoplasm, genetic predisposition to cancer
Citation: Yang R, Zhan Y, Li Y, Dai S-Y, He S-W, Ye C-J, Meng L-D, Chen D-Q, Dong C-B, Chen L, Chen G, Dong K-R, Li K, Zheng S, Li J, Yao W and Dong R (2022) The Cellular and Molecular Landscape of Synchronous Pediatric Sialoblastoma and Hepatoblastoma. Front. Oncol. 12:893206. doi: 10.3389/fonc.2022.893206
Received: 10 March 2022; Accepted: 31 May 2022;
Published: 04 July 2022.
Edited by:
Jaume Mora, Hospital Sant Joan de Déu Barcelona, SpainReviewed by:
Monica Sanchez Guixe, Institute for Research in Biomedicine, SpainPawel Buczkowicz, PhenoTips, Canada
Copyright © 2022 Yang, Zhan, Li, Dai, He, Ye, Meng, Chen, Dong, Chen, Chen, Dong, Li, Zheng, Li, Yao and Dong. 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: Rui Dong, rdong@fudan.edu.cn; Wei Yao, 2842094044@qq.com; Jun Li, lovenvoyy@126.com
†These authors have contributed equally to this work and share first authorship