Celine, a long interspersed nuclear element retrotransposon, colonizes in the centromeres of poplar chromosomes.

Centromeres in most multicellular eukaryotes are composed of long arrays of repetitive DNA sequences. Interestingly, several transposable elements, including the well-known long terminal repeat (LTR) retrotransposon CRM (centromeric retrotransposon of maize), were found to be enriched in functional centromeres marked by the centromeric histone H3 (CENH3). Here we report a centromeric long interspersed nuclear element (LINE), Celine, in Populus species. Celine has colonized preferentially in the CENH3-associated chromatin of every poplar chromosome, with 84% of the Celine elements localized in the CENH3-binding domains. By contrast, only 51% of the CRM elements were bound to CENH3 domains in Populus trichocarpa. These results suggest different centromere targeting mechanisms employed by Celine and CRM elements. Nevertheless, the high target specificity seems to be detrimental to further amplification of the Celine elements, leading to a shorter life span and patchy distribution among plant species compared to the CRM elements. Using a phylogenetically guided approach we were able to identify Celine-like LINE elements in tea plant (Camellia sinensis) and green ash tree (Fraxinus pennsylvanica). The centromeric localization of these Celine-like LINEs was confirmed in both species. We demonstrate that the centromere targeting property of Celine-like LINEs is of primitive origin and has been conserved among distantly related plant species.

Tumor-associated macrophage-induced circMRCKα encodes a peptide to promote glycolysis and progression in hepatocellular carcinoma.

The tumor-associated macrophages (TAMs) play multifaceted roles in the progression of hepatocellular carcinoma (HCC). However, the involvement of circular RNAs in the interplay between TAMs and HCC remains unclear. Based on Transwell co-culturing and circular RNA sequencing, this study revealed that TAMs enhanced tumor glycolysis and progression by upregulating circMRCKα in HCC cells. Patients with HCC who exhibited elevated circMRCKα levels presented significantly reduced overall survival and greater cumulative recurrence. Notably, we identified a novel functional peptide of 227 amino acids named circMRCKα-227aa, encoded by circMRCKα. Mechanistically, circMRCKα-227aa bound to USP22 and enhanced its protein level to obstruct HIF-1α degradation via the ubiquitin-proteasome pathway, thereby augmenting HCC glycolysis and progression. In clinical HCC samples, a positive correlation was observed between the expression of circMRCKα and the number of infiltrating CD68+ TAMs and expression of USP22. Furthermore, circMRCKα emerged as an independent prognostic risk factor both individually and in conjunction with CD68+ TAMs and USP22. This study illustrated that circMRCKα-227aa, a novel TAM-induced peptide, promotes tumor glycolysis and progression via USP22 binding and HIF-1α upregulation, suggesting that circMRCKα and TAMs could be combined as therapeutic targets in HCC.

Characterization of tumor microbiome and associations with prognosis in intrahepatic cholangiocarcinoma.

BACKGROUND: The tumor microbiome has been characterized in several malignancies; however, no previous studies have investigated its role in intrahepatic cholangiocarcinoma (ICC). Hence, we explored the tumor microbiome and its association with prognosis in ICC. METHODS: One hundred and twenty-one ICC tumor samples and 89 adjacent normal tissues were profiled by 16S rRNA sequencing. Microbial differences between tumor and adjacent nontumoral liver tissues were assessed. Tumor microbial composition was then evaluated to detect its association with prognosis. Finally, a risk score calculated by the tumor microbiota was accessed by the least absolute shrinkage and selector operator method (Lasso) to predict prognosis of ICC. RESULTS: The tumor microbiome displayed a greater diversity than that in adjacent nontumoral liver tissues. Tumor samples were characterized by a higher abundance of Firmicutes, Actinobacteria, Bacteroidetes, and Acidobacteriota. Higher tumor microbial α diversity was associated with lymph node metastasis and predicted shortened overall survival (OS) and recurrence-free survival (RFS). A total of 11 bacteria were selected to generate the risk score by Lasso. This score showed potential in predicting OS, and was an independent risk factor for OS. CONCLUSION: In conclusion, our study characterized the tumor microbiome and revealed its role in predicting prognosis in ICC.

Molecular dissection of an intronic enhancer governing cold-induced expression of the vacuolar invertase gene in potato.

Potato (Solanum tuberosum) is the third most important food crop in the world. Potato tubers must be stored at cold temperatures to minimize sprouting and losses due to disease. However, cold temperatures strongly induce the expression of the potato vacuolar invertase gene (VInv) and cause reducing sugar accumulation. This process, referred to as "cold-induced sweetening," is a major postharvest problem for the potato industry. We discovered that the cold-induced expression of VInv is controlled by a 200 bp enhancer, VInvIn2En, located in its second intron. We identified several DNA motifs in VInvIn2En that bind transcription factors involved in the plant cold stress response. Mutation of these DNA motifs abolished VInvIn2En function as a transcriptional enhancer. We developed VInvIn2En deletion lines in both diploid and tetraploid potato using clustered regularly interspaced short palindromic repeat (CRISPR)/CRISPR-associated nuclease 9 (Cas9)-mediated gene editing. VInv transcription in cold-stored tubers was significantly reduced in the deletion lines. Interestingly, the VInvIn2En sequence is highly conserved among distantly related Solanum species, including tomato (Solanum lycopersicum) and other non-tuber-bearing species. We conclude that the VInv gene and the VInvIn2En enhancer have adopted distinct roles in the cold stress response in tubers of tuber-bearing Solanum species.

Association of BRAF Variants With Disease Characteristics, Prognosis, and Targeted Therapy Response in Intrahepatic Cholangiocarcinoma.

Importance: BRAF variants are associated with tumor progression; however, the prevalence of BRAF variant subtypes and their association with disease characteristics, prognosis, and targeted therapy response in patients with intrahepatic cholangiocarcinoma (ICC) are largely unknown. Objective: To explore the association of BRAF variant subtypes with disease characteristics, prognosis, and targeted therapy response in patients with ICC. Design, Setting, and Participants: In this cohort study, 1175 patients who underwent curative resection for ICC from January 1, 2009, through December 31, 2017, were evaluated at a single hospital in China. Whole-exome sequencing, targeted sequencing, and Sanger sequencing were performed to identify BRAF variants. The Kaplan-Meier method and log-rank test were used to compare overall survival (OS) and disease-free survival (DFS). Univariate and multivariate analyses were performed using Cox proportional hazards regression. Associations between BRAF variants and targeted therapy response were tested in 6 BRAF-variant, patient-derived organoid lines and in 3 of the patient donors of those lines. Data were analyzed from June 1, 2021, to March 15, 2022. Interventions: Hepatectomy in patients with ICC. Main Outcomes and Measures: The association of BRAF variant subtypes with OS and DFS. Results: Of 1175 patients with ICC, the mean (SD) age was 59.4 (10.4) years and 701 (59.7%) were men. A total of 20 different subtypes of BRAF somatic variance affecting 49 patients (4.2%) were identified; V600E was the most frequent allele in this cohort, accounting for 27% of the identified BRAF variants, followed by K601E (14%), D594G (12%), and N581S (6%). Compared with patients with non-V600E BRAF variants, patients with BRAF V600E variants were more likely to have large tumor size (10 of 13 [77%] vs 12 of 36 [33%]; P = .007), multiple tumors (7 of 13 [54%] vs 8 of 36 [22%]; P = .04), and more vascular/bile duct invasion (7 of 13 [54%] vs 8 of 36 [22%]; P = .04). Multivariate analysis revealed that BRAF V600E variants, but not overall BRAF variants or non-V600E BRAF variants, were associated with poor OS (hazard ratio [HR], 1.87; 95% CI, 1.05-3.33; P = .03) and DFS (HR, 1.66; 95% CI, 1.03-2.97; P = .04). There were also broad differences among organoids with different BRAF variant subtypes in sensitivity to BRAF or MEK inhibitors. Conclusions and Relevance: The findings of this cohort study suggest that there are broad differences among organoids with different BRAF variant subtypes in sensitivity to BRAF or MEK inhibitors. Identifying and classifying BRAF variants may be able to help guide precise treatment for patients with ICC.

Identification and functional validation of super-enhancers in Arabidopsis thaliana.

Super-enhancers (SEs) are exceptionally large enhancers and are recognized to play prominent roles in cell identity in mammalian species. We surveyed the genomic regions containing large clusters of accessible chromatin regions (ACRs) marked by deoxyribonuclease (DNase) I hypersensitivity in Arabidopsis thaliana. We identified a set of 749 putative SEs, which have a minimum length of 1.5 kilobases and represent the top 2.5% of the largest ACR clusters. We demonstrate that the genomic regions associating with these SEs were more sensitive to DNase I than other nonpromoter ACRs. The SEs were preferentially associated with topologically associating domains. Furthermore, the SEs and their predicted cognate genes were frequently associated with organ development and tissue identity in A. thaliana. Therefore, the A. thaliana SEs and their cognate genes mirror the functional characteristics of those reported in mammalian species. We developed CRISPR/Cas-mediated deletion lines of a 3,578-bp SE associated with the thalianol biosynthetic gene cluster (BGC). Small deletions (131-157 bp) within the SE resulted in distinct phenotypic changes and transcriptional repression of all five thalianol genes. In addition, T-DNA insertions in the SE region resulted in transcriptional alteration of all five thalianol genes. Thus, this SE appears to play a central role in coordinating the operon-like expression pattern of the thalianol BGC.

The role of tumor-infiltrating B cells in the tumor microenvironment of hepatocellular carcinoma and its prognostic value: a bioinformatics analysis.

Background: Accumulating evidence indicates that tumor heterogeneity is characterized by distinct immunosubtypes. However, prior studies have mainly focused on the functions of T cells. The role of tumor-infiltrating B cells in the microenvironment of hepatocellular carcinoma (HCC) requires further investigation. Methods: We conducted an integrative analysis of single cell RNA sequencing (scRNA-seq) datasets in HCC tumor samples from Gene Expression Omnibus database. We analyzed the features of B cells in normal liver tissue and HCC. Additionally, we conducted a deconvolution analysis using the matrix of scRNA-seq datasets and the RNA-seq datasets in The Cancer Genome Atlas-Liver Hepatocellular Carcinoma (TCGA-LIHC) database. The survival analyses of the TCGA-LIHC cohort with different B cell infiltration rates and was further validated. Finally, we performed immunohistochemistry analysis of primary tumor tissue of HCC patients using antibodies against CD79A and validated the impact of tumor-infiltrating B cells in the prognosis of LIHC. Results: We identified several subtypes of B cells in the microenvironment of HCC, including the plasma cells and naïve B cells. The relative ratio of B cells, but not the plasma cells, was significantly decreased in HCC as compared to the normal liver tissue (P<0.05). In addition, genes related to antigen presentation and cell proliferation were decreased in tumor-infiltrating B cells (P<0.05). The observation of B cell infiltration was further validated with the TCGA-LIHC cohort. The overall survival and disease-free survival in HCC patients with higher B-cell infiltration rate were significantly longer than those in the lower infiltration group (P<0.05) in the TCGA-LIHC cohort. Moreover, we demonstrated higher infiltration rates of B cells were significantly associated with a better prognosis of HCC in our cohort. Conclusions: Tumor-infiltrating B cells potentially exert a tumor-suppressive function in the microenvironment of HCC and the higher levels of B cell infiltration are associated with a favorable outcome of HCC. Targeted activation of B cells may improve the tumor immune-targeted therapy.

Chromosome-scale genome assemblies and annotations for Poales species Carex cristatella, Carex scoparia, Juncus effusus, and Juncus inflexus.

The majority of sequenced genomes in the monocots are from species belonging to Poaceae, which include many commercially important crops. Here, we expand the number of sequenced genomes from the monocots to include the genomes of 4 related cyperids: Carex cristatella and Carex scoparia from Cyperaceae and Juncus effusus and Juncus inflexus from Juncaceae. The high-quality, chromosome-scale genome sequences from these 4 cyperids were assembled by combining whole-genome shotgun sequencing of Nanopore long reads, Illumina short reads, and Hi-C sequencing data. Some members of the Cyperaceae and Juncaceae are known to possess holocentric chromosomes. We examined the repeat landscapes in our sequenced genomes to search for potential repeats associated with centromeres. Several large satellite repeat families, comprising 3.2-9.5% of our sequenced genomes, showed dispersed distribution of large satellite repeat clusters across all Carex chromosomes, with few instances of these repeats clustering in the same chromosomal regions. In contrast, most large Juncus satellite repeats were clustered in a single location on each chromosome, with sporadic instances of large satellite repeats throughout the Juncus genomes. Recognizable transposable elements account for about 20% of each of the 4 genome assemblies, with the Carex genomes containing more DNA transposons than retrotransposons while the converse is true for the Juncus genomes. These genome sequences and annotations will facilitate better comparative analysis within monocots.

Tumor-derived exosomes induce immunosuppressive macrophages to foster intrahepatic cholangiocarcinoma progression.

BACKGROUND AND AIMS: Macrophages are prominent components of solid tumors and exhibit distinct functions in different tumor microenvironments. Exosomes are emerging as necessary mediators of the cross-talk between tumor cells and the microenvironment. However, the underlying mechanisms of exosomes involving into crosstalk between tumor cells and macrophages during disease progression of intrahepatic cholangiocarcinoma (ICC) have not been yet fully realized. APPROACH AND RESULTS: We found that the macrophages of ICC tumor tissues up-regulated the expression levels of immunosuppressive molecule programmed death-ligand 1 (PD-L1). Increased PD-L1+ macrophages in tumor tissues effectively suppressed T-cell immunity and correlated with poor survival rates in patients with ICC. High-throughput RNA-sequencing analysis that was performed to identify differential levels of microRNAs (miRNAs) between exosomes derived from ICC cells and primary human intrahepatic biliary epithelial cells revealed that miR-183-5p was increased in ICC cell-derived exosomes. Exosomal miR-183-5p inhibited phosphatase and tensin homolog (PTEN) expression, to subsequently affect the elevations on both phosphorylated AKT and PD-L1 expression in macrophages. Furthermore, macrophages that treated with ICC cell-derived exosomes significantly suppressed T-cell immunity in vitro and contributed to the growth and progression of ICC in vivo, which were reversible through blockages on PD-L1 of these macrophages. Finally, clinical data showed that up-regulated levels of plasma exosomal miR-183-5p correlated with poor prognosis of patients with ICC after curative resection. CONCLUSIONS: Tumor-derived exosomal miR-183-5p up-regulates PD-L1-expressing macrophages to foster immune suppression and disease progression in ICC through the miR-183-5p/PTEN/AKT/PD-L1 pathway. Exosomal miR-183-5p is a potential predictive biomarker for ICC progression and a potential target for development of therapeutic strategies against immune tolerance feature of ICC.

Whole-genome sequencing reveals the evolutionary trajectory of HBV-related hepatocellular carcinoma early recurrence.

Patients with hepatocellular carcinoma (HCC) have poor long-term survival following curative resection because of the high rate of tumor early recurrence. Little is known about the trajectory of genomic evolution from primary to early-recurrent HCC. In this study, we performed whole-genome sequencing (WGS) on 40 pairs of primary and early-recurrent hepatitis B virus (HBV)-related HCC tumors from patients who received curative resection, and from four patients whose primary and recurrent tumor were extensively sampled. We identified two recurrence patterns: de novo recurrence (18/40), which developed genetically independently of the primary tumor and carried different HCC drivers, and ancestral recurrence (22/40), which was clonally related to the primary tumor and progressed more rapidly than de novo recurrence. We found that the recurrence location was predictive of the recurrence pattern: distant recurrence tended to display the de novo pattern, whereas local recurrence tended to display the ancestral pattern. We then uncovered the evolutionary trajectories based on the subclonal architecture, driver-gene mutations, and mutational processes observed in the primary and recurrent tumors. Multi-region WGS demonstrated spatiotemporal heterogeneity and polyclonal, monophyletic dissemination in HCC ancestral recurrence. In addition, we identified recurrence-specific mutations and copy-number gains in BCL9, leading to WNT/ß-catenin signaling activation and an immune-excluded tumor microenvironment, which suggests that BCL9 might serve as a new therapeutic target for recurrent HCC. Collectively, our results allow us to view with unprecedented clarity the genomic evolution during HBV-related HCC early recurrence, providing an important molecular foundation for enhanced understanding of HCC with implications for personalized therapy to improve patient survival.

Association of KRAS Variant Subtypes With Survival and Recurrence in Patients With Surgically Treated Intrahepatic Cholangiocarcinoma.

Importance: KRAS variants are associated with tumor progression; however, the prevalence of KRAS variant subtypes and their association with survival and recurrence in patients with intrahepatic cholangiocarcinoma (ICC) after curative resection are largely unknown. Objective: To explore the prognostic association of KRAS variant subtypes with survival and recurrence in patients with ICC. Design, Setting, and Participants: In this cohort study, patients who underwent curative resection for ICC from January 2009 through December 2016 at a single hospital in China were recruited, and whole-exome sequencing, targeted sequencing, and Sanger sequencing were performed to identify KRAS variants. Kaplan-Meier and log-rank tests were used to compare overall survival (OS) and disease-free survival (DFS). Univariate and multivariate analyses were performed using the Cox proportional hazards regression model. Data were analyzed from April 2020 to January 2021. Interventions: Hepatectomy in patients with ICC. Main Outcomes and Measures: The association of KRAS variant subtypes with OS and DFS. Results: Of 1024 included patients with ICC, 621 (60.6%) were male, and the mean (SD) age was 59.2 (10.2) years. A total of 14 different subtypes of KRAS somatic variants affecting 127 patients (12.4%) were identified. G12D was the most frequent allele in this cohort, accounting for 55 of 127 identified KRAS variants (43.3%), followed by G12V (25 [19.7%]), G12C (9 [7.1%]), and G13D (8 [6.3%]). Compared with patients with wild-type KRAS, patients with variant KRAS were more likely to have high levels of carbohydrate antigen 19-9 (92 of 127 [72.4%] vs 546 of 897 [60.9%]; P = .01) and γ-glutamyltransferase (72 of 127 [56.7%] vs 420 of 897 [46.8%]; P = .04). Multivariable analysis revealed that G12 KRAS variants but not non-G12 KRAS variants were independently associated with worse OS (hazard ratio [HR], 1.69; 95% CI, 1.31-2.18; P < .001) and DFS (HR, 1.47; 95% CI, 1.16-1.88; P = .002). Among the patients with G12 KRAS variants, the G12V KRAS variant was the strongest prognostic determinant for the worst OS (HR, 3.05; 95% CI, 1.94-4.79; P < .001) and DFS (HR, 1.79; 95% CI, 1.13-2.85; P = .01). Conclusions and Relevance: In this cohort study, the distribution of KRAS variant subtypes was characterized in a large cohort of patients with ICC from China. The presence of G12 KRAS variants but not non-G12 KRAS variants was associated with worse survival and increased risk of recurrence. Patients with the G12V variant exhibited the worst outcomes in the whole cohort.

Genomic evolution and the impact of SLIT2 mutation in relapsed intrahepatic cholangiocarcinoma.

BACKGROUND AND AIMS: Intrahepatic cholangiocarcinoma (ICC) is aggressive and has high rates of relapse, conferring poor long-term survival after curative resection. Little is known about the genomic evolution that occurs during ICC relapse. APPROACH AND RESULTS: We conducted whole-exome sequencing of 30 paired primary and relapsed tumors from 10 patients with ICC who received curative resection. We sought to identify frequently altered genes, infer tumor subclonal architectures, and track genomic evolution from primary to relapsed tumors. We examined functional effects and the mechanism of action of SLIT2, a gene specifically mutated in relapsed tumors, on tumor growth and metastasis and the tumor microenvironment. Our results indicated that relapsed ICCs were genetically derived from intrahepatic dissemination of primary tumors. However, they acquired additional mutations while maintaining most drivers, such as TP53 and IDH1. Multiregion sequencing suggested polyclonal seeding of ICC dissemination. Four of 10 relapsed ICCs acquired SLIT2 mutations that were not present in the corresponding primary tumors. Validation in an expanded sample revealed SLIT2 mutations in 2.3% (1/44) of primary ICCs and 29.5% (13/44) of relapsed ICCs. Biofunctional investigations revealed that inactivating mutation of SLIT2 resulted in activation of PI3K-Akt signaling in ICC cells, directly enhanced neutrophil chemotaxis, mediated tumor-associated neutrophil infiltration, and contributed to ICC growth and metastasis. CONCLUSIONS: We characterized genomic evolution during ICC relapse and identified SLIT2 as a driver of tumor dissemination and tumor-associated neutrophil infiltration.

Characterization of immune infiltration in sarcomatoid hepatocellular carcinoma.

Sarcomatoid hepatocellular carcinoma (sHCC) is a rare type of liver malignancy. Currently, the tumor immune features of sHCC are poorly understood. We recruited 31 patients with resected sHCC for whom tissue samples and complete clinicopathologic and follow-up data were available. To understand the immune infiltration of sHCC, immunohistochemical staining was performed on the resected sHCC samples to compare the expressions of programmed death-1 (PD-1), programmed death-ligand 1 (PD-L1), B7-H3, indoleamine 2,3-dioxygenase (IDO), lymphocyte-activation gene 3 (LAG-3), CD8, FOXP3, and CD68 in tumor and peritumoral tissues. Kaplan-Meier and Cox regression analyses were used to assess the predictive value of immune markers. Sarcomatoid components were characterized with significantly higher expression of PD-L1 and B7-H3 in tumor cells than in conventional HCC components, as well as in peritumoral tissue. Additionally, sarcomatoid components had a higher density of FOXP3+ and LAG-3+ cells and a lower density of CD8+ cells than conventional HCC components or peritumoral tissue. Higher expression of PD-L1 in tumor cells significantly correlated with higher densities of CD8+, PD-1+, and LAG-3+ cells. Increased tumor PD-L1 expression and decreased CD8+ T-cell density were associated with poor overall survival (OS) and disease-free survival (DFS) in patients of sHCC. These findings suggest further characterization on relative mechanism of sHCC immune infiltration may identify therapeutic targets for immunotherapy.

LINC01133 promotes hepatocellular carcinoma progression by sponging miR-199a-5p and activating annexin A2.

BACKGROUND: Long noncoding RNAs (lncRNAs) are functionally associated with cancer development and progression. Although gene copy number variation (CNV) is common in hepatocellular carcinoma (HCC), it is not known how CNV in lncRNAs affects HCC progression and recurrence. We aimed to identify a CNV-related lncRNA involved in HCC progression and recurrence and illustrate its underlying mechanisms and prognostic value. METHODS: We analyzed the whole genome sequencing (WGS) data of matched cancerous and noncancerous liver samples from 49 patients with HCC to identify lncRNAs with CNV. The results were validated in another cohort of 238 paired HCC and nontumor samples by TaqMan copy number assay. We preformed Kaplan-Meier analysis and log-rank test to identify lncRNA CNV with prognostic value. We conducted loss- and gain-of-function studies to explore the biological functions of LINC01133 in vitro and in vivo. The competing endogenous RNAs (ceRNAs) mechanism was clarified by microRNA sequencing (miR-seq), quantitative real-time PCR (qRT-PCR), western blot, and dual-luciferase reporter assays. We confirmed the binding mechanism between lncRNA and protein by RNA pull-down, RNA immunoprecipitation, qRT-PCR, and western blot analyses. RESULTS: Genomic copy numbers of LINC01133 were increased in HCC, which were positively related with the elevated expression of LINC01133. Increased copy number of LINC01133 predicted the poor prognosis in HCC patients. LINC01133 overexpression in HCC cells promoted proliferation and aggressive phenotypes in vitro, and facilitated tumor growth and lung metastasis in vivo, whereas LINC01133 knockdown had the opposite effects. LINC01133 sponged miR-199a-5p, resulting in enhanced expression of SNAI1, which induced epithelial-to-mesenchymal transition (EMT) in HCC cells. In addition, LINC01133 interacted with Annexin A2 (ANXA2) to activate the ANXA2/STAT3 signaling pathway. CONCLUSIONS: LINC01133 promotes HCC progression by sponging miR-199a-5p and interacting with ANXA2. LINC01133 CNV gain is predictive of poor prognosis in patients with HCC.

Tumor-associated neutrophils and macrophages interaction contributes to intrahepatic cholangiocarcinoma progression by activating STAT3.

BACKGROUND: Tumor-associated neutrophils (TANs) and macrophages (TAMs) can each influence cancer growth and metastasis, but their combined effects in intrahepatic cholangiocarcinoma (ICC) remain unclear. METHODS: We explored the distributions of TANs and TAMs in patient-derived ICC samples by multiplex immunofluorescent staining and tested their separate and combined effects on ICC in vitro and in vivo. We then investigated the mechanistic basis of the effects using PCR array, western blot analysis and ELISA experiments. Finally, we validated our results in a tissue microarray composed of primary tumor tissues from 359 patients with ICC. RESULTS: The spatial distributions of TANs and TAMs were correlated with each other in patient-derived ICC samples. Interaction between TANs and TAMs enhanced the proliferation and invasion abilities of ICC cells in vitro and tumor progression in a mouse xenograft model of ICC. TANs and TAMs produced higher levels of oncostatin M and interleukin-11, respectively, in co-culture than in monoculture. Both of those cytokines activated STAT3 signaling in ICC cells. Knockdown of STAT3 abolished the protumor effect of TANs and TAMs on ICC. In tumor samples from patients with ICC, increased TAN and TAM levels were correlated with elevated p-STAT3 expression. All three of those factors were independent predictors of patient outcomes. CONCLUSIONS: TANs and TAMs interact to promote ICC progression by activating STAT3.

Associations among the mutational landscape, immune microenvironment, and prognosis in Chinese patients with hepatocellular carcinoma.

Recent studies suggested that the immune microenvironment and mutational landscape are associated with the response to immune-based therapy in several types of cancer. The roles of those factors in Chinese HCC remain largely unknown. In this study, we obtained 182 FFPE samples of HCC cohort that were previously subjected to NGS (49 WGS, 18 WES, and 115 targeted sequencing). We performed immunohistochemistry to detect CD3, CD4, CD8, CD57, Foxp3, CD68, CD66b, and PD-L1 expression in the samples. We identified diverse associations between the mutational landscape and the immune microenvironment in the HCC samples. High mutational burden and an aristolochic acid-dominated mutational signature were both correlated with elevated tumoral PD-L1 expression and CD3+ T-cell infiltration and high numbers of CD68+ TAMs and CD66b+ TANs. CD4+ and CD8+ T cells exhibited lower infiltration levels in tumors with mutations in AXIN1/CTNNB1 and in tumors with aflatoxin-dominant mutational signatures. Moreover, tumors with TP53 mutations had less CD8+ T-cell infiltration and more Foxp3+ Treg-cell infiltration than those without TP53 mutations. Kaplan-Meier survival analysis revealed that the presence of CD8+, Foxp3+, CD66b+, or CD68+ immune cells; tumoral PD-L1 expression alone; or the presence of CD8+ or Foxp3+ cells combined with TP53 mutation were predictive of recurrence and poor overall survival after curative resection. In conclusion, the association between the mutational landscape and the immune microenvironment warrants further analysis to determine its impact on patient outcomes to guide personalized immune-based therapy for Chinese patients with HCC.

Peritumoral plasmacytoid dendritic cells predict a poor prognosis for intrahepatic cholangiocarcinoma after curative resection.

BACKGROUND: Plasmacytoid dendritic cells (pDCs) are present in various primary and metastatic human neoplasms; however, their clinical significance in intrahepatic cholangiocarcinoma is not clear. METHODS: To evaluate pDCs' distributions in and around tumors as well as their potential function and predictive value for prognosis in patients undergoing curative resection, we performed immunohistochemistry to examine the expression of pDC marker BDCA2, and CD3, CD4, CD8 and Foxp3 in intratumoral and peritumoral tissues from 359 patients with intrahepatic cholangiocarcinoma and compared with prognostic and clinicopathologic factors. RESULTS: Results showed that patients with high numbers of BDCA2+ pDCs in peritumoral tissues were more likely to have elevated levels of carbohydrate antigen 19-9 and gamma-glutamyl transferase, larger and more tumors, advanced tumor-node-metastasis staging, more vascular/bile duct invasion, and lymphatic metastasis in association with greater chance of recurrence and shorter overall survival. Peritumoral tissues with larger numbers of pDCs also showed increased Foxp3+ regulatory T cell infiltration, both of which were found to be independent factors for predicting time to recurrence and overall survival. By contrast, patient outcomes were not associated with the presence of intratumoral pDCs. CONCLUSIONS: Peritumoral pDC infiltration may indicate an immune tolerogenic peritumor microenvironment and can be used to predict a poor prognosis for patients undergoing curative resection for intrahepatic cholangiocarcinoma.

MACROD2 deficiency promotes hepatocellular carcinoma growth and metastasis by activating GSK-3ß/ß-catenin signaling.

Structural variations (SVs) influence the development and progression of multiple types of cancer. The genes affected by SVs in hepatocellular carcinoma (HCC) and their contribution to tumor growth and metastasis remain unknown. In this study, through whole-genome sequencing (WGS), we identified MACROD2 as the gene most frequently affected by SVs, which were associated with low MACROD2 expression levels. Low MACROD2 expression was predictive of tumor recurrence and poor overall survival. MACROD2 expression was decreased in HCC cell lines, especially those with high metastatic potential. MACROD2 knockdown in HCC cells markedly enhanced proliferation and invasiveness in vitro and tumor progression in vivo and promoted epithelial-mesenchymal transition (EMT). By contrast, MACROD2 overexpression reversed EMT and inhibited HCC growth and metastasis. Mechanistically, MACROD2 deficiency suppressed glycogen synthase kinase-3ß (GSK-3ß) activity and activated ß-catenin signaling, which mediated the effect of MACROD2 on HCC. In clinical HCC samples, decreased MACROD2 expression was correlated with the activation of GSK-3ß/ß-catenin signaling and the EMT phenotype. Overall, our results revealed that MACROD2 is frequently affected by SVs in HCC, and its deficiency promotes tumor growth and metastasis by activating GSK-3ß/ß-catenin signaling.

An extraordinarily stable karyotype of the woody Populus species revealed by chromosome painting.

The karyotype represents the basic genetic make-up of a eukaryotic species. Comparative cytogenetic analysis of related species based on individually identified chromosomes has been conducted in only a few plant groups and not yet in woody plants. We have developed a complete set of 19 chromosome painting probes based on the reference genome of the model woody plant Populus trichocarpa. Using sequential fluorescence in situ hybridization we were able to identify all poplar chromosomes in the same metaphase cells, which led to the development of poplar karyotypes based on individually identified chromosomes. We demonstrate that five Populus species, belonging to five different sections within Populus, have maintained a remarkably conserved karyotype. No inter-chromosomal structural rearrangements were observed on any of the 19 chromosomes among the five species. Thus, the chromosomal synteny in Populus has been remarkably maintained after nearly 14 million years of divergence. We propose that the karyotypes of woody species are more stable than those of herbaceous plants since it may take a longer period of time for woody plants to fix chromosome number or structural variants in natural populations.