Phosphorylation relays are a major mechanism by which bacteria alter transcription in response to environmental signals, but understanding of the functional consequences of bacterial response regulator phosphorylation is limited. We sought to characterize how phosphorylation of the control of virulence regulator (CovR) protein from the major human pathogen group A Streptococcus (GAS) influences GAS global gene expression and pathogenesis. CovR mainly serves to repress GAS virulence factor-encoding genes and has been shown to homodimerize following phosphorylation on aspartate-53 (D53) in vitro. We discovered that CovR is phosphorylated in vivo and that such phosphorylation is partially heat-stable, suggesting additional phosphorylation at non-aspartate residues. Using mass spectroscopy along with targeted mutagenesis, we identified threonine-65 (T65) as an additional CovR phosphorylation site under control of the serine/threonine kinase (Stk). Phosphorylation on T65, as mimicked by the recombinant CovR T65E variant, abolished in vitro CovR D53 phosphorylation. Similarly, isoallelic GAS strains that were either unable to be phosphorylated at D53 (CovR-D53A) or had functional constitutive phosphorylation at T65 (CovR-T65E) had essentially an identical gene repression profile to each other and to a CovR-inactivated strain. However, the CovR-D53A and CovR-T65E isoallelic strains retained the ability to positively influence gene expression that was abolished in the CovR-inactivated strain. Consistent with these observations, the CovR-D53A and CovR-T65E strains were hypervirulent compared to the CovR-inactivated strain in a mouse model of invasive GAS disease. Surprisingly, an isoalleic strain unable to be phosphorylated at CovR T65 (CovR-T65A) was hypervirulent compared to the wild-type strain, as auto-regulation of covR gene expression resulted in lower covR gene transcript and CovR protein levels in the CovR-T65A strain. Taken together, these data establish that CovR is phosphorylated in vivo and elucidate how the complex interplay between CovR D53 activating phosphorylation, T65 inhibiting phosphorylation, and auto-regulation impacts streptococcal host-pathogen interaction.

The genetically diverse viridans group streptococci (VGS) are increasingly recognized as the cause of a variety of human diseases. We used a recently developed multilocus sequence analysis scheme to define the species of 118 unique VGS strains causing bacteremia in patients with cancer; Streptococcus mitis (68 patients) and S. oralis (22 patients) were the most frequently identified strains. Compared with patients infected with non-S. mitis strains, patients infected with S. mitis strains were more likely to have moderate or severe clinical disease (e.g., VGS shock syndrome). Combined with the sequence data, whole-genome analyses showed that S. mitis strains may more precisely be considered as >2 species. Furthermore, we found that multiple S. mitis strains induced disease in neutropenic mice in a dose-dependent fashion. Our data define the prominent clinical effect of the group of organisms currently classified as S. mitis and lay the groundwork for increased understanding of this understudied pathogen.

To define the mutation spectrum in non-Down syndrome acute megakaryoblastic leukemia (non-DS-AMKL), we performed transcriptome sequencing on diagnostic blasts from 14 pediatric patients and validated our findings in a recurrency/validation cohort consisting of 34 pediatric and 28 adult AMKL samples. Our analysis identified a cryptic chromosome 16 inversion (inv(16)(p13.3q24.3)) in 27% of pediatric cases, which encodes a CBFA2T3-GLIS2 fusion protein. Expression of CBFA2T3-GLIS2 in Drosophila and murine hematopoietic cells induced bone morphogenic protein (BMP) signaling and resulted in a marked increase in the self-renewal capacity of hematopoietic progenitors. These data suggest that expression of CBFA2T3-GLIS2 directly contributes to leukemogenesis.

Cutaneous squamous cell carcinoma (cuSCC) comprises 15-20% of all skin cancers, accounting for over 700,000 cases in USA annually. Most cuSCC arise in association with a distinct precancerous lesion, the actinic keratosis (AK). To identify potential targets for molecularly targeted chemoprevention, here we perform integrated cross-species genomic analysis of cuSCC development through the preneoplastic AK stage using matched human samples and a solar ultraviolet radiation-driven Hairless mouse model. We identify the major transcriptional drivers of this progression sequence, showing that the key genomic changes in cuSCC development occur in the normal skin to AK transition. Our data validate the use of this ultraviolet radiation-driven mouse cuSCC model for cross-species analysis and demonstrate that cuSCC bears deep molecular similarities to multiple carcinogen-driven SCCs from diverse sites, suggesting that cuSCC may serve as an effective, accessible model for multiple SCC types and that common treatment and prevention strategies may be feasible.

Tumors of the same histologic type often comprise clinically and molecularly distinct subgroups; however, the etiology of these subgroups is unknown. Here, we report that histologically identical, but genetically distinct, ependymomas exhibit patterns of gene expression that recapitulate those of radial glia cells in the corresponding region of the central nervous system. Cancer stem cells isolated from ependymomas displayed a radial glia phenotype and formed tumors when orthotopically transplanted in mice. These findings identify restricted populations of radial glia cells as candidate stem cells of the different subgroups of ependymoma, and they support a general hypothesis that subgroups of the same histologic tumor type are generated by different populations of progenitor cells in the tissues of origin.

Substantial heterogeneity exists within cervical cancer that is generally infected by human papillomavirus (HPV). However, the most common histological subtype of cervical cancer, cervical squamous cell carcinoma (CSCC), is poorly characterized regarding the association between its heterogeneity and HPV oncoprotein expression. We filtered out 138 CSCC samples with infection of HPV16 only as the first step; then we compressed HPV16 E6/E7 expression as HPVpca and correlated HPVpca with the immunological profiling of CSCC based on supervised clustering to discover subtypes and to characterize the differences between subgroups in terms of the HPVpca level, pathway activity, epigenetic dysregulation, somatic mutation frequencies, and likelihood of responding to chemo/immunotherapies. Supervised clustering of immune signatures revealed two HPV16 subtypes (namely, HPV16-IMM and HPV16-KRT) that correlated with HPVpca and clinical outcomes. HPV16-KRT is characterized by elevated expression of genes in keratinization, biological oxidation, and Wnt signaling, whereas HPV16-IMM has a strong immune response and mesenchymal features. HPV16-IMM exhibited much more epigenetic silencing and significant mutation at FBXW7, while MUC4 and PIK3CA were mutated frequently for HPV16-KRT. We also imputed that HPV16-IMM is much more sensitive to chemo/immunotherapy than is HPV16-KRT. Our characterization tightly links the expression of HPV16 E6/E7 with biological and clinical outcomes of CSCC, providing valuable molecular-level information that points to decoding heterogeneity. Together, these results shed light on stratifications of CSCC infected by HPV16 and shall help to guide personalized management and treatment of patients.

Tumor sequencing studies have emphasized the role of epigenetics and altered chromatin homeostasis in cancer. Mutations in DAXX, which encodes a chaperone for the histone 3.3 variant, occur in 25% of pancreatic neuroendocrine tumors (PanNETs). To advance our understanding of physiological functions of Daxx, we developed a conditional Daxx allele in mice. We demonstrate that Daxx loss is well tolerated in the pancreas but creates a permissive transcriptional state that cooperates with environmental stress (inflammation) and other genetic lesions (Men1 loss) to alter gene expression and cell state, impairing pancreas recovery from inflammatory stress in vivo. The transcriptional changes are associated with dysregulation of endogenous retroviral elements (ERVs), and dysregulation of endogenous genes near ERVs is also observed in human PanNETs with DAXX mutations. Our results reveal a physiologic function of DAXX, provide a mechanism associated with impaired tissue regeneration and tumorigenesis, and expand our understanding of ERV regulation in somatic cells.

Human endogenous retroviruses (HERVs), which make up approximately 8% of the human genome, are overexpressed in some breast cancer cells and tissues but without regard to cancer subtype. We, therefore, analyzed TCGA RNA-Seq data to evaluate differences in expression of the HERV-K family in breast cancers of the various subtypes. Four HERV-K loci on different chromosomes were analyzed in basal, Her2E, LumA, and LumB breast cancer subtypes of 512 breast cancer patients with invasive ductal carcinoma (IDC). The results for all four loci showed higher HERV-K expression in the basal subtype, suggesting similar mechanisms of regulation regardless of locus. Expression of the HERV-K envelope gene (env) was highly significantly increased in basal tumors in comparison with the also-upregulated expression of other HERV-K genes. Analysis of reverse-phase protein array data indicated that increased expression of HERV-K is associated with decreased mutation of H-Ras (wild-type). Our results show elevation of HERV-K expression exclusively in the basal subtype of IDC breast cancer (as opposed to the other subtypes) and suggest HERV-K as a possible target for cancer vaccines or immunotherapy against this highly aggressive form of breast cancer.

Crosstalk between genetic, epigenetic, and immune alterations in upper tract urothelial carcinomas and their role in shaping muscle invasiveness and patient outcome are poorly understood.

Coat color is of great importance in animal breed characteristics; it is not only a significant productive trait but also an indispensable economic trait, especially in the rabbit industry. In the present study, the relationship between melanocortin 1 receptor (MC1R) genotypes and coat color phenotypes was observed in five rabbit breeds with popular coat colors that are present in China. These breeds comprised the Lianshan black rabbit (BR), Fujian yellow rabbit (YR), New Zealand white rabbit (WR), Gray Giant rabbit (GR), and Checkered Giant rabbit (CR), which were firstly determined, and the results showed that GR had an E allele; WR, CR, and BR had a 6-bp in-frame deletion (c.281_286del6, ED allele); and YR had a 30-bp deletion (c.304_333del30 E allele). To explore the feasibility of obtaining a novel rabbit coat color through the mutation of MC1R with the CRISPR/Cas9 system, two single-guide RNAs (sgRNAs) were designed for the MC1R gene, and the editing efficiency was confirmed by injection of rabbits' zygotes. Unlike the donor rabbits whose coat color was originally black, two novel pale-yellow-coated rabbits were generated in the founders. A total of six novel MC1R gene deletions were identified in the two founder rabbits, in which the longest deletion was more than 700 bp. The histological hematoxylin-and-eosin (H&E) staining results indicated that eumelanin amounts were absent in hair follicles of MC1R-knockout (KO) rabbits, when compared with that of donor BR. In addition, the messenger RNA (mRNA) levels of some key downstream genes in the MC1R pathway were all downregulated in MC1R-KO rabbits compared with BR and YR. These results further indicate that loss-of-function MC1R contributed to blocking the synthesis of eumelanin and created a novel pale-yellow coat color in the MC1R-KO rabbits, and gene editing technology may be a useful tool to generate novel phenotypes in rabbit breeding.

Renal medullary carcinoma (RMC) is a highly lethal malignancy that mainly afflicts young individuals of African descent and is resistant to all targeted agents used to treat other renal cell carcinomas. Comprehensive genomic and transcriptomic profiling of untreated primary RMC tissues was performed to elucidate the molecular landscape of these tumors. We found that RMC was characterized by high replication stress and an abundance of focal copy-number alterations associated with activation of the stimulator of the cyclic GMP-AMP synthase interferon genes (cGAS-STING) innate immune pathway. Replication stress conferred a therapeutic vulnerability to drugs targeting DNA-damage repair pathways. Elucidation of these previously unknown RMC hallmarks paves the way to new clinical trials for this rare but highly lethal malignancy.

Esthesioneuroblastoma (ENB) is a rare cancer of the olfactory mucosa, with no established molecular stratification to date. We report similarities of ENB with tumors arising in the neural crest and perform integrative analysis of these tumors. We propose a molecular-based subtype classification of ENB as basal or neural, both of which have distinct pathological, transcriptomic, proteomic, and immune features. Among the basal subtype, we uncovered an IDH2 R172 mutant-enriched subgroup (∼35%) harboring a CpG island methylator phenotype reminiscent of IDH2 mutant gliomas. Compared with the basal ENB methylome, the neural ENB methylome shows genome-wide reprogramming with loss of DNA methylation at the enhancers of axonal guidance genes. Our study reveals insights into the molecular pathogenesis of ENB and provides classification information of potential therapeutic relevance.

Renal medullary carcinoma (RMC) is an aggressive tumour driven by bi-allelic loss of SMARCB1 and tightly associated with sickle cell trait. However, the cell-of-origin and oncogenic mechanism remain poorly understood. Using single-cell sequencing of human RMC, we defined transformation of thick ascending limb (TAL) cells into an epithelial-mesenchymal gradient of RMC cells associated with loss of renal epithelial transcription factors TFCP2L1, HOXB9 and MITF and gain of MYC and NFE2L2-associated oncogenic and ferroptosis resistance programs. We describe the molecular basis for this transcriptional switch that is reversed by SMARCB1 re-expression repressing the oncogenic and ferroptosis resistance programs leading to ferroptotic cell death. Ferroptosis resistance links TAL cell survival with the high extracellular medullar iron concentrations associated with sickle cell trait, an environment propitious to the mutagenic events associated with RMC development. This unique environment may explain why RMC is the only SMARCB1-deficient tumour arising from epithelial cells, differentiating RMC from rhabdoid tumours arising from neural crest cells.

Missense mutations in the DNA binding domain of p53 are characterized as structural or contact mutations based on their effect on the conformation of the protein. These mutations show gain-of-function (GOF) activities, such as promoting increased metastatic incidence compared with p53 loss, often mediated by the interaction of mutant p53 with a set of transcription factors. These interactions are largely context specific. To understand the mechanisms by which p53 DNA binding domain mutations drive osteosarcoma progression, we created mouse models, in which either the p53 structural mutant p53R172H or the contact mutant p53R245W are expressed specifically in osteoblasts, yielding osteosarcoma tumor development. Survival significantly decreased and metastatic incidence increased in mice expressing p53 mutants compared with p53-null mice, suggesting GOF. RNA sequencing of primary osteosarcomas revealed vastly different gene expression profiles between tumors expressing the missense mutants and p53-null tumors. Further, p53R172H and p53R245W each regulated unique transcriptomes and pathways through interactions with a distinct repertoire of transcription factors. Validation assays showed that p53R245W, but not p53R172H, interacts with KLF15 to drive migration and invasion in osteosarcoma cell lines and promotes metastasis in allogeneic transplantation models. In addition, analyses of p53R248W chromatin immunoprecipitation peaks showed enrichment of KLF15 motifs in human osteoblasts. Taken together, these data identify unique mechanisms of action of the structural and contact mutants of p53.

Protein ubiquitination regulated by ubiquitin ligases plays important roles in innate immunity. However, key regulators of ubiquitination during innate response and roles of new types of ubiquitination (apart from Lys48- and Lys63-linkage) in control of innate signaling have not been clearly understood. Here we report that F-box only protein Fbxo21, a functionally unknown component of SCF (Skp1-Cul1-F-box protein) complex, facilitates Lys29-linkage and activation of ASK1 (apoptosis signal-regulating kinase 1), and promotes type I interferon production upon viral infection. Fbxo21 deficiency in mice cells impairs virus-induced Lys29-linkage and activation of ASK1, attenuates c-Jun N-terminal kinase (JNK) and p38 signaling pathway, and decreases the production of proinflammatory cytokines and type I interferon, resulting in reduced antiviral innate response and enhanced virus replication. Therefore Fbxo21 is required for ASK1 activation via Lys29-linkage of ASK1 during antiviral innate response, providing mechanistic insights into non-proteolytic roles of SCF complex in innate immune response.

Kawasaki disease (KD) is an acute, self‑limited vasculitis that predominantly affects medium‑sized arteries, particularly the coronary arteries. Recent studies have indicated that microRNAs are involved in many diseases, including KD. However, the detailed mechanism remains unclear. The aim of the present study was to explore the role of miR‑186 in KD and potentially discover a new target for KD treatment. The results demonstrated that miR‑186 was upregulated in serum from patients with KD and KD serum could increase miR‑186 transcript levels in endothelial cells (HUVECs). Overexpression of miR‑186 mimic induced HUVEC apoptosis through mitogen‑activated protein kinase (MAPK) activation by targeting and inhibiting SMAD family member 6 (SMAD6). Furthermore, KD serum induced HUVEC apoptosis through miR‑186. In conclusion, the present results suggested that KD serum‑associated miR‑186 has an essential role in endothelial cell apoptosis by activating the MAPK pathway through targeting the SMAD6 gene.

Blastic plasmacytoid dendritic cell neoplasm (BPDCN) is a rare hematologic malignancy with poor outcomes with conventional therapy. Nearly 100% of BPDCNs overexpress interleukin 3 receptor subunit alpha (CD123). Given that CD123 is differentially expressed on the surface of BPDCN cells, it has emerged as an attractive therapeutic target. UCART123 is an investigational product consisting of allogeneic T cells expressing an anti-CD123 chimeric antigen receptor (CAR), edited with TALEN® nucleases. In this study, we examine the antitumor activity of UCART123 in preclinical models of BPDCN. We report that UCART123 have selective antitumor activity against CD123-positive primary BPDCN samples (while sparing normal hematopoietic progenitor cells) in the in vitro cytotoxicity and T cell degranulation assays; supported by the increased secretion of IFNγ by UCART123 cells when cultured in the presence of BPDCN cells. UCART123 eradicate BPDCN and result in long-term disease-free survival in a subset of primary patient-derived BPDCN xenograft mouse models. One potential challenge of CD123 targeting therapies is the loss of CD123 antigen through diverse genetic mechanisms, an event observed in one of three BPDCN PDX studied. In summary, these results provide a preclinical proof-of-principle that allogeneic UCART123 cells have potent anti-BPDCN activity.

Sarcomatoid clear-cell renal cell carcinomas (sRCC) are associated with dismal prognosis. Genomic alterations associated with sarcomatoid dedifferentiation are poorly characterized. We sought to define the genomic landscape of sRCC and uncover potentially actionable therapeutic targets. We assessed the genomic landscape of sRCC using targeted panel sequencing including patients with microdissected sarcomatoid and epithelial components. Along with common genomic alterations associated with clear-cell histology, we found that Hippo was one of the most frequently altered pathways in these tumours. Hippo alterations were differentially enriched in sRCC compared to non-sRCC. Functional analysis showed that Hippo members mutations were associated with higher nuclear accumulation of YAP/TAZ, core effectors of the Hippo pathway. In a NF2-mutant sRCC model, YAP1 knockdown and NF2 reconstitution suppressed cell proliferation, tumour growth and invasion, both in vitro and in vivo. Overall, we show that Hippo pathway alterations are a feature of sRCC, and enable the exploration of the Hippo pathway as a novel potential therapeutic target.

Objectives: This study was designed to identify a messenger RNA (mRNA) expression signature to predict survival in patients with oral squamous cell carcinoma (OSCC). Methods: mRNA expression profiles were integrated with clinical data from 280 samples, including 19 normal tissues and 261 OSCC tissues in The Cancer Genome Atlas. We identified differentially expressed mRNAs (DEmRNAs) between the OSCC and normal tissue samples and developed a novel mRNA-focused expression signature using a Cox regression analysis and other bioinformatic methods. The prognostic value of this signature was evaluated by Kaplan-Meier analysis, multivariable COX regression, and receiver operating characteristic (ROC) curve analysis. Protein-protein interaction (PPI) network, gene ontology, and Kyoto Encyclopedia of Genes and Genomes enrichment analysis were performed to predict the function of the DEmRNAs. Signature-related mRNAs were analyzed by gene set enrichment analyses (GSEA) and validated by quantitative real-time polymerase chain reaction (qRT-PCR) in 20 paired OSCC and adjacent healthy tissues. Results: We identified a novel 5-mRNA expression signature (HOXA1, CELSR3, HIST1H3J, ZFP42, and ASCL4) that could predict patient outcomes in OSCC. The risk score based on the signature was able to separate OSCC patients into high- and low-risk groups that showed significantly different overall survival (p < 0.001, log-rank test). The signature was further validated as an effective independent prognostic predictor of OSCC by multivariate Cox regression analysis (hazard ratio = 3.747, confidence interval: 2.279-5.677, p < 0.001) and ROC curve of the third year (area under the curve = 0.733). Functional analysis demonstrated that the key hub genes in the PPI network were mainly enriched in cell division, cell proliferation, and the p53 signaling pathway. GSEA results showed that the 5 mRNAs were significantly enriched in mismatch repair, DNA replication, and the NOTCH signaling pathway. Finally, qRT-PCR results showed that the 5 mRNAs were upregulated in OSCC tissue in agreement with the predictions from our bioinformatics analysis. Conclusions: We identified a novel 5-mRNA signature that could predict the survival of patients with OSCC and may be a promising biomarker for personalized cancer treatments.

Racial disparities have been documented in the biology and outcome of certain renal cell carcinomas (RCCs) among Black patients. However, little is known about racial differences in MiT family translocation RCC (TRCC). To investigate this issue, we performed a case-control study using data from The Cancer Genome Atlas (TCGA) and the Chinese OrigiMed2020 cohort. A total of 676 patients with RCC (14 Asian, 113 Black, and 525 White) were identified in TCGA, and TRCC was defined as RCC with TFE3/TFEB translocation or TFEB amplification, leading to 21 patients with TRCC (2 Asian, 8 Black, 10 White, and 1 unknown). Asian (2 of 14 [14.3%] vs 10 of 525 [1.9%]; P = .036) and Black (8 of 113 [7.1%] vs 1.9%; P = .007) patients with RCC showed significantly higher prevalence of TRCC compared with White patients with RCC. The overall mortality rate of TRCC was slightly higher in Asian and Black patients compared with White patients (HR: 6.05, P = .069). OrigiMed2020 Chinese patients with RCC had a significantly higher proportion of TRCC with TFE3 fusions than TCGA White patients with RCC (13 of 250 [5.2%] vs 7 of 525 [1.3%]; P = .003). Black patients with TRCC were more likely to exhibit the proliferative subtype than White patients (6 of 8 [75%] vs 2 of 9 [22.2%]; P = .057) for those who had RNA-seq profiles. We present evidence of higher prevalence of TRCC in Asian and Black patients with RCC compared with White patients and show that these tumors in Asian and Black patients have distinct transcriptional signatures and are associated with poor outcomes.