Genomic evolution and diverse models of systemic metastases in colorectal cancer.

OBJECTIVE: The systemic spread of colorectal cancer (CRC) is dominated by the portal system and exhibits diverse patterns of metastasis without systematical genomic investigation. Here, we evaluated the genomic evolution of CRC with multiorgan metastases using multiregion sequencing. DESIGN: Whole-exome sequencing was performed on multiple regions (n=74) of matched primary tumour, adjacent non-cancerous mucosa, liver metastasis and lung metastasis from six patients with CRC. Phylogenetic reconstruction and evolutionary analyses were used to investigate the metastatic seeding pattern and clonal origin. Recurrent driver gene mutations were analysed across patients and validated in two independent cohorts. Metastatic assays were performed to examine the effect of the novel driver gene on the malignant behaviour of CRC cells. RESULTS: Based on the migration patterns and clonal origins, three models were revealed (sequential, branch-off and diaspora), which not only supported the anatomic assumption that CRC cells spread to lung after clonally expanding in the liver, but also illustrated the direct seeding of extrahepatic metastases from primary tumours independently. Unlike other cancer types, polyphyletic seeding occurs in CRC, which may result in late metastases with intermetastatic driver gene heterogeneity. In cases with rapid dissemination, we found recurrent trunk loss-of-function mutations in ZFP36L2, which is enriched in metastatic CRC and associated with poor overall survival. CRISPR/Cas9-mediated knockout of ZFP36L2 enhances the metastatic potential of CRC cells. CONCLUSION: Our results provide genomic evidence for metastatic evolution and indicate that biopsy/sequencing of metastases may be considered for patients with CRC with multiorgan or late postoperative metastasis.

Association of recurrent APOBEC3B alterations with the prognosis of gastric-type cervical adenocarcinoma.

OBJECTIVE: Gastric-type cervical adenocarcinoma (GCA) is a rare and aggressive type of endocervical adenocarcinoma (ECA) with distinct histopathologic features and unfavorable treatment outcomes, but no genomic prognostic factor has been revealed. We aimed to systematically investigate the somatic alterations of GCA at genome-wide level and evaluate their prognostic value. METHODS: We performed whole-exome sequencing (WES) on 25 pairs of tumor and matched normal samples to characterize the genomic features of Chinese patients with GCA and investigated their relations to histopathological characterizations and prognosis. The prognostic value of the genomic alterations was evaluated in a total of 58 GCA patients. RESULTS: Mutations were commonly observed in reported GCA-related driver genes, including TP53 (32%), CDKN2A (20%), SKT11 (20%), BRCA2 (12%), SMAD4 (12%), and ERBB2 (12%). Recurrent novel trunk mutations were also observed in PBRM1 (12%), FRMPD4 (12%), and NOP2 (8%) with high variant allele frequency. Moreover, enrichment of the APOBEC signature was attributed to frequent gain of somatic copy number alteration (SCNA) of APOBEC3B (20%), which perfectly matched the nuclear-positive staining of APOBEC3B through immunohistochemistry. In contrast, APOBEC3B alteration was absent in patients with conventional type of ECA (N = 52). Notably, positive APOBEC3B was consistently enriched in patients with favorable prognosis in both the discovery cohort and an additional 33 GCA patients, thus indicating a significant association with lower relapse risk of GCA independent of cancer stage (P = 0.02). CONCLUSION: Our results can aid understanding of the molecular basis of GCA in the Chinese population by providing genomic profiles and highlighting the potential prognostic value of APOBEC3B for GCA through routine clinical IHC.

Characterization of immature ovarian teratomas through single-cell transcriptome.

Introduction: Immature ovarian teratomas are a type of malignant germ cell tumor composed of complicated cell types and are characterized by pathological features of immature neuroectodermal tubules/rosettes. However, there is a lack of understanding of patient-derived immature ovarian teratomas (PDT) at the single cell level. Moreover, whether stem cell lines derived from immature teratomas (CDT) can be used as models for research on PDT remains to be elucidated. Methods: Single-cell RNA sequencing (scRNA-seq) and subsequent bioinformatic analysis was performed on three patient-derived immature ovarian teratomas (PDT) samples to reveal the heterogeneity, evolution trajectory, and cell communication within the tumor microenvironment of PDT. Validations were conducted in additional seven samples through multiplex immunofluorescence. Result: A total of qualified 22,153 cells were obtained and divided into 28 clusters, which can match to the scRNA-seq annotation of CDT as well as human fetal Cell Atlas, but with higher heterogeneity and more prolific cell-cell crosstalk. Radial glia cells (tagged by SOX2) and immature neuron (tagged by DCX) exhibited mutually exclusive expression and differentiated along distinct evolutionary trajectory from cycling neural progenitors. Proportions of these neuroectodermal cell subtypes may play important roles in PDT through contributing to the internal heterogeneity of PDTs. Moreover, the immune cells in PDTs were infiltrated rather than teratoma-derived, with more abundant macrophage in immature neuron than those in radial glia cells, and the infiltrated macrophage subtypes (i.e., M1 and M2) were significantly correlated to clinical grade. Overall, suppressed evolution process and transcriptome regulation in neuroectodermal cells, reduced cell-cell crosstalk, higher M1/M2 proportion ratio, and enhanced T cell effects in tumor microenvironment are enriched in patients with favorable prognosis. Discussion: This study provides a comprehensive profile of PDT at the single cell level, shedding light on the heterogeneity and evolution of neuroectodermal cells within PDTs and the role of immune cells within the tumor microenvironment. Also, our findings highlight the potential usage of CDTs as a model for research on PDT.

Novel genetic associations with five aesthetic facial traits: A genome-wide association study in the Chinese population.

Background: The aesthetic facial traits are closely related to life quality and strongly influenced by genetic factors, but the genetic predispositions in the Chinese population remain poorly understood. Methods: A genome-wide association studies (GWAS) and subsequent validations were performed in 26,806 Chinese on five facial traits: widow's peak, unibrow, double eyelid, earlobe attachment, and freckles. Functional annotation was performed based on the expression quantitative trait loci (eQTL) variants, genome-wide polygenic scores (GPSs) were developed to represent the combined polygenic effects, and single nucleotide polymorphism (SNP) heritability was presented to evaluate the contributions of the variants. Results: In total, 21 genetic associations were identified, of which ten were novel: GMDS-AS1 (rs4959669, p = 1.29 × 10-49) and SPRED2 (rs13423753, p = 2.99 × 10-14) for widow's peak, a previously unreported trait; FARSB (rs36015125, p = 1.96 × 10-21) for unibrow; KIF26B (rs7549180, p = 2.41 × 10-15), CASC2 (rs79852633, p = 4.78 × 10-11), RPGRIP1L (rs6499632, p = 9.15 × 10-11), and PAX1 (rs147581439, p = 3.07 × 10-8) for double eyelid; ZFHX3 (rs74030209, p = 9.77 × 10-14) and LINC01107 (rs10211400, p = 6.25 × 10-10) for earlobe attachment; and SPATA33 (rs35415928, p = 1.08 × 10-8) for freckles. Functionally, seven identified SNPs tag the missense variants and six may function as eQTLs. The combined polygenic effect of the associations was represented by GPSs and contributions of the variants were evaluated using SNP heritability. Conclusion: These identifications may facilitate a better understanding of the genetic basis of features in the Chinese population and hopefully inspire further genetic research on facial development.

Pan-cancer single-cell analysis reveals the heterogeneity and plasticity of cancer-associated fibroblasts in the tumor microenvironment.

Cancer-associated fibroblasts (CAFs) are the predominant components of the tumor microenvironment (TME) and influence cancer hallmarks, but without systematic investigation on their ubiquitous characteristics across different cancer types. Here, we perform pan-cancer analysis on 226 samples across 10 solid cancer types to profile the TME at single-cell resolution, illustrating the commonalities/plasticity of heterogenous CAFs. Activation trajectory of the major CAF types is divided into three states, exhibiting distinct interactions with other cell components, and relating to prognosis of immunotherapy. Moreover, minor CAF components represent the alternative origin from other TME components (e.g., endothelia and macrophages). Particularly, the ubiquitous presentation of endothelial-to-mesenchymal transition CAF, which may interact with proximal SPP1+ tumor-associated macrophages, is implicated in endothelial-to-mesenchymal transition and survival stratifications. Our study comprehensively profiles the shared characteristics and dynamics of CAFs, and highlight their heterogeneity and plasticity across different cancer types. Browser of integrated pan-cancer single-cell information is available at https://gist-fgl.github.io/sc-caf-atlas/ .

Novel Recurrent Altered Genes in Chinese Patients With Anaplastic Thyroid Cancer.

BACKGROUND: Anaplastic thyroid cancer (ATC) is a rare but lethal malignancy, and few systematic investigations on genomic profiles of ATC have been performed in Chinese patients. METHODS: Fifty-four ATC patients in West China Hospital between 2010 to 2020 were retrospectively analyzed, while 29 patients with available samples were sequenced by whole-exome sequencing (WES). The associations between genomic alterations and clinical characteristics were statistically evaluated. RESULTS: The median overall survival was 3.0 months in the entire cohort, which was impacted by multiple clinical features, including age, tumor size, and different treatment strategies. In the WES cohort, totally 797 nonsilent mutations were detected; the most frequently altered genes were TP53 (48%), BRAF (24%), PIK3CA (24%), and TERT promoter (21%). Although these mutations have been well-reported in previous studies, ethnic specificity was exhibited in terms of mutation frequency. Moreover, several novel significantly mutated genes were identified including RBM15 (17%), NOTCH2NL (14%), CTNNA3 (10%), and KATNAL2 (10%). WES-based copy number alteration analysis also revealed a high frequent gain of NOTCH2NL (41%), which induced its increased expression. Gene mutations and copy number alterations were enriched in phosphatidylinositol 3-kinase/AKT/mechanistic target of rapamycin (mTOR), NOTCH, and WNT pathways. CONCLUSIONS: This study reveals shared and ethnicity-specific genomic profiles of ATC in Chinese patients and suggests NOTCH2NL may act as a novel candidate driver gene for ATC tumorigenesis.

Age-related differences of genetic susceptibility to patients with acute lymphoblastic leukemia.

Inherited predispositions to acute lymphoblastic leukemia have been well investigated in pediatric patients, but studies on adults, particularly Chinese patients, are limited. In this study, we conducted a genome-wide association study in 466 all-age Chinese patients with Acute lymphoblastic leukemia (ALL) and 1,466 non-ALL controls to estimate the impact of age on ALL susceptibility in the Chinese population. Among the 17 reported loci, 8 have been validated in pediatric and 1 in adult patients. The strongest association signal was identified at ARID5B locus and gradually decreased with age, while the signal at GATA3 exhibited the opposite trend and significantly impact on adult patients. With genome-wide approaches, germline variants at 2q14.3 rank as the top inherited predisposition to adult patients (e.g., rs73956024, P = 4.3 × 10-5) and separate the genetic risk of pediatric vs. adult patients (P = 3.6 × 10-6), whereas variants at 15q25.3 (e.g., rs11638062) have a similar impact on patients in different age groups (overall P = 2.9 × 10-7). Our analysis highlights the impact of age on genetic susceptibility to ALL in Chinese patients.

Screening of Novel Pharmacogenetic Candidates for Mercaptopurine-Induced Toxicity in Patients With Acute Lymphoblastic Leukemia.

A small proportion of patients with acute lymphoblastic leukemia (ALL) may experience severe leukopenia after treating with 6-mercaptopurine (6MP), which can be largely explained by germline variants in TPMT and NUDT15. However, a minority of patients who suffered such adverse drug reaction have NUDT15 wt/wt TPMT wt/wt genotype, indicating that other genetic factors may take part in. In this study, we genotyped 539 exon-located nonsilent pharmacogenetic variants in genes involved in phase I/II of drug metabolism in 173 pediatric patients with ALL and conducted association screening for 6MP-induced leukopenia. Besides NUDT15 (rs116855232, P = 6.4 × 10-11) and TPMT (rs1142345, P = 0.003), a novel variant was identified in CYP2A7 gene (i.e., rs73032311, P = 0.0007), which is independent of NUDT15/TPMT variant. In addition, a variant (i.e., rs4680) in COMT is significantly associated with 6MP-induced hepatotoxicity (P = 0.007). In conclusion, variants in CYP2A7 and COMT may be considered as novel potential pharmacogenetic markers for 6MP-induced toxicities, but additional independent validations with large sample size and investigations on related mechanisms are further needed.