Hyperactivation of ß-catenin signal in hepatocellular carcinoma recruits myeloid-derived suppressor cells through PF4-CXCR3 axis.

The high mutation rate of CTNNB1 (37 %) and Wnt-ß-catenin signal-associated genes (54 %) has been notified in hepatocellular carcinoma (HCC). The activation of Wnt-ß-catenin signal pathway was reported to be associated with an immune "desert" phenotype, but the underlying mechanism remains unclear. Here we mainly employed orthotopic HCC models to explore on it. Mass cytometry depicted the immune contexture of orthotopic HCC syngeneic grafts, unveiling that the exogenous expression of ß-catenin significantly increased the percentage of myeloid-derived suppressor cells (MDSCs) and decreased the percentage of CD8+ T-cells. Flow cytometry and immunohistochemistry further confirmed the findings. The protein microarray analysis, Western blot and PCR identified PF4 as its downstream regulating cytokine. Intratumorally injection of cytokine PF4 enhanced the accumulation of MDSCs. Knockout of PF4 abolished the effect of ß-catenin on recruiting MDSCs. Chromatin immunoprecipitation and luciferase reporter assay demonstrated that ß-catenin increases the mRNA level of PF4 via binding to PF4's promoter region. In vitro chemotaxis assay and in vivo administration of specific inhibitors identified CXCR3 on MDSCs as receptor for recruiting PF4. Lastly, the significant correlations across ß-catenin, PF4 and MDSCs and CD8+ T-cells infiltration were verified in HCC clinical samples. Our results unveiled HCC tumor cell intrinsic hyperactivation of ß-catenin can recruit MDSC through PF4-CXCR3, which contributes to the formation of immune "desert" phenotype. Our study provided new insights into the development of immunotherapeutic strategy of HCC with CTNNB1 mutation. SIGNIFICANCE: This study identifies PF4-CXCR3-MDSCs as a downstream mechanism underlying CTNNB1 mutation associated immune "desert" phenotype.

Mutant Allele Frequency-Based Intra-Tumoral Genetic Heterogeneity Related to the Tumor Shrinkage Mode After Neoadjuvant Chemotherapy in Breast Cancer Patients.

The shrinkage mode of tumor extent after neoadjuvant chemotherapy (NAC) is an important index to evaluate the odds of breast-conserving surgery. However, there is no sufficient measurement to predict the shrinkage mode after NAC. In this study, we analyzed 24 patients' formalin-fixed, paraffin-embedded samples before and after treatment and analyzed 456 cancer-related genes panel by using target next-generation sequencing. Meanwhile, the pathological shrinkage mode was reconstructed in three dimensions after surgery, and the genetic heterogeneity level was estimated by mutant-allele tumor heterogeneity (MATH). We measured the genetic intra-tumor heterogeneity and explored its correlation with the shrinkage mode after NAC. A total of 17 matched pair samples of primary tumor tissue and residual tumor tissue were successfully accessed. It was found that the most common mutated genes were TP53 and PIK3CA in both samples before and after NAC, and no recurrent mutations were significantly associated with the shrinkage mode. Besides, the MATH value of formalin-fixed, paraffin-embedded samples before and after NAC was analyzed by the area under the curve of the receiver operating characteristic, and it is feasible to classify patients into concentric shrinkage mode and non-concentric shrinkage mode in NAC based on the MATH threshold of 58. Our findings indicate that the MATH value was associated with the shrinkage mode of breast cancer in a non-linear model. Patients with the MATH value below the threshold of 58 before and after NAC displayed a concentric shrinkage mode. The area under the curve was 0.89, with a sensitivity of 0.69 and specificity of 1. Our study might provide a promising application of intra-tumor heterogeneity that is measured by MATH to make a choice of surgery.

Genetic profiling of primary and secondary tumors from patients with lung adenocarcinoma and bone metastases reveals targeted therapy options.

BACKGROUND: Patients newly diagnosed with lung adenocarcinoma with bone metastases (LABM) have poor survival rates after treatment with conventional therapies. To improve outcomes, we retrospectively investigated whether the application of a more comprehensive genetic test of tumor biopsies samples from LABM patients could provide the basis for treatment with more effective tyrosine kinase inhibitors (TKIs) regimens. METHODS: Fine needle biopsies were taken from the primary tumor (PT) and a secondary bone metastasis (BM) of 17 LABM patients before treatment. Simple genetic profiles for selecting therapies were initially obtained using an ARMS-PCR test for EGFR and ALK fusion mutations. More detailed genetic profiles of somatic exon SNVs and CNVs in 457 cancer-related genes were retrospectively derived using capture single molecule amplification and resequencing technology (capSMART). RESULTS: ARMS-PCR identified 14 EGFR positive, 3 EGFR negative and 1 ALK fusion positive patient. A therapy regimen incorporating TKIs Gefitinib and Crizotinib was offered to the EGFR and ALK fusion positive patients, respectively. With the exception of two patients, molecular profiling of matching PT and BM biopsies identified a highly shared somatic variant fingerprint, although the BMs exhibited additional genomic instability. In six of 13 EGFR positive patients and in all three EGFR negative patients, examination of the genetic profiles identified additional clinically significant mutations that are known or experimental drug targets for treatment of lung cancer. CONCLUSION: Our findings firstly suggest that treatment regimens based on comprehensive genetic assessment of newly diagnosed LABM patients should target both the PT and secondary BMs, including rogue clones with potential to form new BMs. Second, the additional information gained should allow clinicians to design and implement more personalized treatment regimens and potentially improve outcomes for LABM patients.

Contributions of Zea mays subspecies mexicana haplotypes to modern maize.

Maize was domesticated from lowland teosinte (Zea mays ssp. parviglumis), but the contribution of highland teosinte (Zea mays ssp. mexicana, hereafter mexicana) to modern maize is not clear. Here, two genomes for Mo17 (a modern maize inbred) and mexicana are assembled using a meta-assembly strategy after sequencing of 10 lines derived from a maize-teosinte cross. Comparative analyses reveal a high level of diversity between Mo17, B73, and mexicana, including three Mb-size structural rearrangements. The maize spontaneous mutation rate is estimated to be 2.17 × 10-8 ~3.87 × 10-8 per site per generation with a nonrandom distribution across the genome. A higher deleterious mutation rate is observed in the pericentromeric regions, and might be caused by differences in recombination frequency. Over 10% of the maize genome shows evidence of introgression from the mexicana genome, suggesting that mexicana contributed to maize adaptation and improvement. Our data offer a rich resource for constructing the pan-genome of Zea mays and genetic improvement of modern maize varieties.

Comparative genomic and functional analyses: unearthing the diversity and specificity of nematicidal factors in Pseudomonas putida strain 1A00316.

We isolated Pseudomonas putida (P. putida) strain 1A00316 from Antarctica. This bacterium has a high efficiency against Meloidogyne incognita (M. incognita) in vitro and under greenhouse conditions. The complete genome of P. putida 1A00316 was sequenced using PacBio single molecule real-time (SMRT) technology. A comparative genomic analysis of 16 Pseudomonas strains revealed that although P. putida 1A00316 belonged to P. putida, it was phenotypically more similar to nematicidal Pseudomonas fluorescens (P. fluorescens) strains. We characterized the diversity and specificity of nematicidal factors in P. putida 1A00316 with comparative genomics and functional analysis, and found that P. putida 1A00316 has diverse nematicidal factors including protein alkaline metalloproteinase AprA and two secondary metabolites, hydrogen cyanide and cyclo-(l-isoleucyl-l-proline). We show for the first time that cyclo-(l-isoleucyl-l-proline) exhibit nematicidal activity in P. putida. Interestingly, our study had not detected common nematicidal factors such as 2,4-diacetylphloroglucinol (2,4-DAPG) and pyrrolnitrin in P. putida 1A00316. The results of the present study reveal the diversity and specificity of nematicidal factors in P. putida strain 1A00316.

Functional analysis of long intergenic non-coding RNAs in phosphate-starved rice using competing endogenous RNA network.

Long intergenic non-coding RNAs (lincRNAs) may play widespread roles in gene regulation and other biological processes, however, a systematic examination of the functions of lincRNAs in the biological responses of rice to phosphate (Pi) starvation has not been performed. Here, we used a computational method to predict the functions of lincRNAs in Pi-starved rice. Overall, 3,170 lincRNA loci were identified using RNA sequencing data from the roots and shoots of control and Pi-starved rice. A competing endogenous RNA (ceRNA) network was constructed for each tissue by considering the competing relationships between lincRNAs and genes, and the correlations between the expression levels of RNAs in ceRNA pairs. Enrichment analyses showed that most of the communities in the networks were related to the biological processes of Pi starvation. The lincRNAs in the two tissues were individually functionally annotated based on the ceRNA networks, and the differentially expressed lincRNAs were biologically meaningful. For example, XLOC_026030 was upregulated from 3 days after Pi starvation, and its functional annotation was 'cellular response to Pi starvation'. In conclusion, we systematically annotated lincRNAs in rice and identified those involved in the biological response to Pi starvation.