An integrated genomic approach identifies follistatin as a target of the p63-epidermal growth factor receptor oncogenic network in head and neck squamous cell carcinoma.

Although numerous putative oncogenes have been associated with the etiology of head and neck squamous cell carcinoma (HNSCC), the mechanisms by which these oncogenes and their downstream targets mediate tumor progression have not been fully elucidated. We performed an integrative analysis to identify a crucial set of targets of the oncogenic transcription factor p63 that are common across multiple transcriptomic datasets obtained from HNSCC patients, and representative cell line models. Notably, our analysis revealed FST which encodes follistatin, a secreted glycoprotein that inhibits the transforming growth factor TGFß/activin signaling pathways, to be a direct transcriptional target of p63. In addition, we found that FST expression is also driven by epidermal growth factor receptor EGFR signaling, thus mediating a functional link between the TGF-ß and EGFR pathways. We show through loss- and gain-of-function studies that FST predominantly imparts a tumor-growth and migratory phenotype in HNSCC cells. Furthermore, analysis of single-cell RNA sequencing data from HNSCC patients unveiled cancer cells as the dominant source of FST within the tumor microenvironment and exposed a correlation between the expression of FST and its regulators with immune infiltrates. We propose FST as a prognostic biomarker for patient survival and a compelling candidate mediating the broad effects of p63 on the tumor and its associated microenvironment.

Robust Performance of SARS-CoV-2 Whole-Genome Sequencing from Wastewater with a Nonselective Virus Concentration Method.

The sequencing of human virus genomes from wastewater samples is an efficient method for tracking viral transmission and evolution at the community level. However, this requires the recovery of viral nucleic acids of high quality. We developed a reusable tangential-flow filtration system to concentrate and purify viruses from wastewater for genome sequencing. A pilot study was conducted with 94 wastewater samples from four local sewersheds, from which viral nucleic acids were extracted, and the whole genome of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) was sequenced using the ARTIC V4.0 primers. Our method yielded a high probability (0.9) of recovering complete or near-complete SARS-CoV-2 genomes (>90% coverage at 10× depth) from wastewater when the COVID-19 incidence rate exceeded 33 cases per 100 000 people. The relative abundances of sequenced SARS-CoV-2 variants followed the trends observed from patient-derived samples. We also identified SARS-CoV-2 lineages in wastewater that were underrepresented or not present in the clinical whole-genome sequencing data. The developed tangential-flow filtration system can be easily adopted for the sequencing of other viruses in wastewater, particularly those at low concentrations.

Multimodal Dimension Reduction and Subtype Classification of Head and Neck Squamous Cell Tumors.

Traditional analysis of genomic data from bulk sequencing experiments seek to group and compare sample cohorts into biologically meaningful groups. To accomplish this task, large scale databases of patient-derived samples, like that of TCGA, have been established, giving the ability to interrogate multiple data modalities per tumor. We have developed a computational strategy employing multimodal integration paired with spectral clustering and modern dimension reduction techniques such as PHATE to provide a more robust method for cancer sub-type classification. Using this integrated approach, we have examined 514 Head and Neck Squamous Carcinoma (HNSC) tumor samples from TCGA across gene-expression, DNA-methylation, and microbiome data modalities. We show that these approaches, primarily developed for single-cell sequencing can be efficiently applied to bulk tumor sequencing data. Our multimodal analysis captures the dynamic heterogeneity, identifies new and refines subtypes of HNSC, and orders tumor samples along well-defined cellular trajectories. Collectively, these results showcase the inherent molecular complexity of tumors and offer insights into carcinogenesis and importance of targeted therapy. Computational techniques as highlighted in our study provide an organic and powerful approach to identify granular patterns in large and noisy datasets that may otherwise be overlooked.

Reactivation of super-enhancers by KLF4 in human Head and Neck Squamous Cell Carcinoma.

Head and neck squamous cell carcinoma (HNSCC) is a disease of significant morbidity and mortality and rarely diagnosed in early stages. Despite extensive genetic and genomic characterization, targeted therapeutics and diagnostic markers of HNSCC are lacking due to the inherent heterogeneity and complexity of the disease. Herein, we have generated the global histone mark based epigenomic and transcriptomic cartogram of SCC25, a representative cell type of mesenchymal HNSCC and its normal oral keratinocyte counterpart. Examination of genomic regions marked by differential chromatin states and associated with misregulated gene expression led us to identify SCC25 enriched regulatory sequences and transcription factors (TF) motifs. These findings were further strengthened by ATAC-seq based open chromatin and TF footprint analysis which unearthed Krüppel-like Factor 4 (KLF4) as a potential key regulator of the SCC25 cistrome. We reaffirm the results obtained from in silico and chromatin studies in SCC25 by ChIP-seq of KLF4 and identify ΔNp63 as a co-oncogenic driver of the cancer-specific gene expression milieu. Taken together, our results lead us to propose a model where elevated KLF4 levels sustains the oncogenic state of HNSCC by reactivating repressed chromatin domains at key downstream genes, often by targeting super-enhancers.

Molecular dissection of the oncogenic role of ETS1 in the mesenchymal subtypes of head and neck squamous cell carcinoma.

Head and Neck Squamous Cell Carcinoma (HNSCC) is a heterogeneous disease of significant mortality and with limited treatment options. Recent genomic analysis of HNSCC tumors has identified several distinct molecular classes, of which the mesenchymal subtype is associated with Epithelial to Mesenchymal Transition (EMT) and shown to correlate with poor survival and drug resistance. Here, we utilize an integrated approach to characterize the molecular function of ETS1, an oncogenic transcription factor specifically enriched in Mesenchymal tumors. To identify the global ETS1 cistrome, we have performed integrated analysis of RNA-Seq, ChIP-Seq and epigenomic datasets in SCC25, a representative ETS1high mesenchymal HNSCC cell line. Our studies implicate ETS1 as a crucial regulator of broader oncogenic processes and specifically Mesenchymal phenotypes, such as EMT and cellular invasion. We found that ETS1 preferentially binds cancer specific regulator elements, in particular Super-Enhancers of key EMT genes, highlighting its role as a master regulator. Finally, we show evidence that ETS1 plays a crucial role in regulating the TGF-ß pathway in Mesenchymal cell lines and in leading-edge cells in primary HNSCC tumors that are endowed with partial-EMT features. Collectively our study highlights ETS1 as a key regulator of TGF-ß associated EMT and reveals new avenues for sub-type specific therapeutic intervention.

Targeting TAZ-Driven Human Breast Cancer by Inhibiting a SKP2-p27 Signaling Axis.

Deregulated expression of the transcriptional coactivator with PDZ-binding motif (WWTR1/TAZ) is a common feature of basal-like breast cancer (BLBC). Yet, how oncogenic TAZ regulates cell-cycle progression and proliferation in breast cancer remains poorly understood, and whether TAZ is required for tumor maintenance has not been established. Here, using an integrative oncogenomic approach, TAZ-dependent cellular programs essential for tumor growth and progression were identified. Significantly, TAZ-driven tumor cells required sustained TAZ expression, given that its withdrawal impaired both genesis and maintenance of solid tumors. Moreover, temporal inhibition of TAZ diminished the metastatic burden in established macroscopic pulmonary metastases. Mechanistic investigation revealed that TAZ controls distinct gene profiles that determine cancer cell fate through cell-cycle networks, including a specific, causal role for S-phase kinase-associated protein 2 (SKP2) in mediating the neoplastic state. Together, this study elucidates the molecular events that underpin the role of TAZ in BLBC and link to SKP2, a convergent communication node for multiple cancer signaling pathways, as a key downstream effector molecule. IMPLICATIONS: Understanding the molecular role of TAZ and its link to SKP2, a signaling convergent point and key regulator in BLBC, represents an important step toward the identification of novel therapeutic targets for TAZ-dependent breast cancer.

RYR1 and CACNA1S genetic variants identified with statin-associated muscle symptoms.

AIM: To examine the genetic differences between subjects with statin-associated muscle symptoms and statin-tolerant controls. MATERIALS & METHODS: Next-generation sequencing was used to characterize the exomes of 76 subjects with severe statin-associated muscle symptoms and 50 statin-tolerant controls. RESULTS: 12 probably pathogenic variants were found within the RYR1 and CACNA1S genes in 16% of cases with severe statin-induced myopathy representing a fourfold increase over variants found in statin-tolerant controls. Subjects with probably pathogenic RYR1 or CACNA1S variants had plasma CK 5X to more than 400X the upper limit of normal in addition to having muscle symptoms. CONCLUSIONS: Genetic variants within the RYR1 and CACNA1S genes are likely to be a major contributor to the susceptibility to statin-associated muscle symptoms.

Chronic vitamin D insufficiency impairs physical performance in C57BL/6J mice.

Vitamin D insufficiency (serum 25-OH vitamin D < 30 ng/ml) affects 70-80% of the general population, yet the long-term impacts on physical performance and the progression of sarcopenia are poorly understood. We therefore followed 6-month-old male C57BL/6J mice (n=6) consuming either sufficient (STD, 1000 IU) or insufficient (LOW, 125 IU) vitamin D3/kg chow for 12 months (equivalent to 20-30 human years). LOW supplemented mice exhibited a rapid decline of serum 25-OH vitamin D levels by two weeks that remained between 11-15 ng/mL for all time points thereafter. After 12 months LOW mice displayed worse grip endurance (34.6 ± 14.1 versus 147.5 ± 50.6 seconds, p=0.001), uphill sprint speed (16.0 ± 1.0 versus 21.8 ± 2.4 meters/min, p=0.0007), and stride length (4.4 ± 0.3 versus 5.1 ± 0.3, p=0.002). LOW mice also showed less lean body mass after 8 months (57.5% ± 5.1% versus 64.5% ± 4.0%, p=0.023), but not after 12 months of supplementation, as well as greater protein expression of atrophy pathway gene atrogin­1. Additionally, microRNA sequencing revealed differential expression of mIR­26a in muscle tissue of LOW mice. These data suggest chronic vitamin D insufficiency may be an important factor contributing to functional decline and sarcopenia.

Thymosin beta 4 induces significant changes in the plasma miRNA profile following severe traumatic brain injury in the rat lateral fluid percussion injury model.

OBJECTIVES: Thymosin beta 4 (Tß4) has demonstrated neuroprotective potential in models of neurlogical injury. The neuroprotective potential of Tß4 has been associated with increased miR-200a and miR-200b within the brain following stroke. Here we tested the hypothesis that Tß4 treatment could also alter miRNA profiles within the plasma following severe traumatic brain injury (TBI). METHODS: We used the rat lateral fluid percusion injury model of severe TBI to test this hypothesis. Highly sensitive and quantitative droplet digital polymerase chain reaction (ddPCR) was used to measure the plasma concentrations of miR-200 family members. In addition, we conducted RNAseq analysis of plasma miRNA to further identify changes associated with TBI and treatment with Tß4. RESULTS: ddPCR demonstrated that miR-200a-3p andmiR-200b-3p were both significantly increased in plasma following treatment with Tß4 after severe TBI. RNAseq analysis suggested that miR-300-3p and miR-598-3p increased while miR-450-3p and miR-194-5p significantly decreased following TBI. In contrast, miR-194-5p significantly increased in Tß4 treated rats following TBI. In addition, we identified nine plasma miRNAs whose expression significantly changed following treatment with Tß4. CONCLUSIONS: Tß4 treatment significantly increased plasma levels of miR-200a-3p and miR-200b-3p, while RNAseq analysis identified miR-194-5p as a candidate miRNA that may be critical for neuroprotection.

Computational approach for deriving cancer progression roadmaps from static sample data.

As with any biological process, cancer development is inherently dynamic. While major efforts continue to catalog the genomic events associated with human cancer, it remains difficult to interpret and extrapolate the accumulating data to provide insights into the dynamic aspects of the disease. Here, we present a computational strategy that enables the construction of a cancer progression model using static tumor sample data. The developed approach overcame many technical limitations of existing methods. Application of the approach to breast cancer data revealed a linear, branching model with two distinct trajectories for malignant progression. The validity of the constructed model was demonstrated in 27 independent breast cancer data sets, and through visualization of the data in the context of disease progression we were able to identify a number of potentially key molecular events in the advance of breast cancer to malignancy.

An automated method for efficient, accurate and reproducible construction of RNA-seq libraries.

BACKGROUND: Integration of RNA-seq expression data with knowledge on chromatin accessibility, histone modifications, DNA methylation, and transcription factor binding has been instrumental for the unveiling of cell-specific local and long-range regulatory patterns, facilitating further investigation on the underlying rules of transcription regulation at an individual and allele-specific level. However, full genome transcriptome characterization has been partially limited by the complexity and increased time-requirements of available RNA-seq library construction protocols. FINDINGS: Use of the SX-8G IP-Star® Compact System significantly reduces the hands-on time for RNA-seq library synthesis, adenylation, and adaptor ligation providing with high quality RNA-seq libraries tailored for Illumina high-throughput next-generation sequencing. Generated data exhibits high technical reproducibility compared to data from RNA-seq libraries synthesized manually for the same samples. Obtained results are consistent regardless the researcher, day of the experiment, and experimental run. CONCLUSIONS: Overall, the SX-8G IP-Star® Compact System proves an efficient, fast and reliable tool for the construction of next-generation RNA-seq libraries especially for trancriptome-based annotation of larger genomes.

Molecular profiling and computational network analysis of TAZ-mediated mammary tumorigenesis identifies actionable therapeutic targets.

Triple-negative breast cancer (TNBC) accounts for approximately 15-20% of all breast cancer (BC) cases and contributes disproportionately to BC mortality. TAZ, a key transducer of the Hippo pathway, has recently been demonstrated to confer breast cancer stem cell (CSC) traits. However, TAZ target genes and the underlying transcriptional regulatory pathways responsible for the CSC phenomenon remain unknown. Here, we demonstrate that the oncogenic activity of TAZ is essential for propagation of the malignant phenotype. We further show that constitutively active TAZ tumor-derived cells exhibit unique tumor-initiating properties, including increased self-renewal and metastatic seeding potential, acquired chemotherapy resistance and the ability to efficiently regenerate tumor formation in vivo. Combined digital RNA expression analysis and computational network approaches identify several signaling pathways that distinguish breast cancer tumor-initiating cells (T-ICs) from bulk tumor cells. We demonstrate the utility of this approach by repositioning the small molecule tyrosine kinase inhibitor, Dasatinib, which selectively targets T-ICs and inhibits TNBC growth in vivo.

Cancer progression modeling using static sample data.

As molecular profiling data continues to accumulate, the design of integrative computational analyses that can provide insights into the dynamic aspects of cancer progression becomes feasible. Here, we present a novel computational method for the construction of cancer progression models based on the analysis of static tumor samples. We demonstrate the reliability of the method with simulated data, and describe the application to breast cancer data. Our findings support a linear, branching model for breast cancer progression. An interactive model facilitates the identification of key molecular events in the advance of disease to malignancy.

FOXO1 regulates expression of a microRNA cluster on X chromosome.

Phosphoinositol-3-kinase (PI3K) pathway is a crucial modulator of many physiological and pathophysiological phenomena, including aging, diabetes and cancer. Protein kinase Akt, a downstream effector of PI3K, controls a plethora of cellular functions, including gene transcription. A key mechanism connecting Akt activity to changes in gene expression is inhibitory phosphorylation of FOXO family of transcription factors. Accordingly, altered expression of FOXO targets may account for many biological consequences of PI3K/Akt signaling. While the previous efforts focused on FOXO-dependent regulation of protein-coding genes, non-coding RNA genes have emerged as equally important targets of many transcription factors. Therefore, we utilized a regulated form of FOXO1 to profile FOXO1-dependent changes in miRNA expression in human cells. Both microarray hybridization and next-generation sequencing revealed changes in the products of a miRNA cluster on X chromosome. Rapid induction of these miRNAs occurred independently of de novo protein synthesis. Furthermore, inhibition of PI3K in cancer cell lines caused derepression of these miRNAs, as would be expected for FOXO-regulated genes. Members of the major oncogenic cascades are significantly overrepresented among the predicted targets of the miRNAs, consistent with tumor-suppressive role of FOXO1. The discovered miRNAs represent new candidate mediators of FOXO1 functions and possible biomarkers of its activity.

Development and characterization of a novel human Waldenström macroglobulinemia cell line: RPCI-WM1, Roswell Park Cancer Institute - Waldenström Macroglobulinemia 1.

Understanding the biology of Waldenström macroglobulinemia is hindered by a lack of preclinical models. We report a novel cell line, RPCI-WM1, from a patient treated for WM. The cell line secretes human immunoglobulin M (h-IgM) with κ-light chain restriction identical to the primary tumor. The cell line has a modal chromosomal number of 46 and harbors chromosomal changes such as deletion of 6q21, monoallelic deletion of 9p21 (CDKN2A), 13q14 (RB1) and 18q21 (BCL-2), with a consistent amplification of 14q32 (immunoglobulin heavy chain; IgH) identical to its founding tumor sample. The clonal relationship is confirmed by identical CDR3 length and single nucleotide polymorphisms as well as a matching IgH sequence of the cell line and founding tumor. Both also harbor a heterozygous, non-synonymous mutation at amino acid 265 in the MYD88 gene (L265P). The cell line expresses most of the cell surface markers present on the parent cells. Overall, RPCI-WM1 represents a valuable model to study Waldenström macroglobulinemia.

Differential copy number aberrations in novel candidate genes associated with progression from in situ to invasive ductal carcinoma of the breast.

Only a minority of intraductal carcinomas of the breast give rise to stromally invasive disease. We microdissected 206 paraffin blocks representing 116 different cases of low-grade ductal carcinoma in situ (DCIS). Fifty-five were pure DCIS (PD) cases without progression to invasive carcinoma. Sixty-one cases had a small invasive component. DNA was extracted from microdissected sections and hybridized to high-density bacterial artificial chromosome arrays. Array comparative genomic hybridization analysis of 118 hybridized DNA samples yielded data on 69 samples that were suitable for further statistical analysis. This cohort included 20 pure DCIS cases, 25 mixed DCIS (MD), and 24 mixed invasive carcinoma samples. PD cases had a higher frequency of DNA copy number changes than MD cases, and the latter had similar DNA profiles compared to paired invasive carcinomas. Copy number changes on 13 chromosomal arms occurred at different rates in PD versus MD lesions. Eight of 19 candidate genes residing at those loci were confirmed to have differential copy number changes by quantitative PCR. NCOR2/SMRT and NR4A1 (both on 12q), DYNLRB2 (16q), CELSR1, UPK3A, and ST13 (all on 22q) were more frequently amplified in PD. Moreover, NCOR2, NR4A1, and DYNLRB2 showed more frequent copy number losses in MD. GRAP2 (22q) was more often amplified in MD, whereas TAF1C (16q) was more commonly deleted in PD. A multigene model comprising these candidate genes discriminated between PD and MD lesions with high accuracy. These findings suggest that the propensity to invade the stroma may be encoded in the genome of intraductal carcinomas.

Common genetic variants are associated with accelerated bone mineral density loss after hematopoietic cell transplantation.

BACKGROUND: Bone mineral density (BMD) loss commonly occurs after hematopoietic cell transplantation (HCT). Hypothesizing that genetic variants may influence post-HCT BMD loss, we conducted a prospective study to examine the associations of single nucleotide polymorphisms (SNP) in bone metabolism pathways and acute BMD loss after HCT. METHODS AND FINDINGS: We genotyped 122 SNPs in 45 genes in bone metabolism pathways among 121 autologous and allogeneic HCT patients. BMD changes from pre-HCT to day +100 post-HCT were analyzed in relation to these SNPs in linear regression models. After controlling for clinical risk factors, we identified 16 SNPs associated with spinal or femoral BMD loss following HCT, three of which have been previously implicated in genome-wide association studies of bone phenotypes, including rs2075555 in COL1A1, rs9594738 in RANKL, and rs4870044 in ESR1. When multiple SNPs were considered simultaneously, they explained 5-35% of the variance in post-HCT BMD loss. There was a significant trend between the number of risk alleles and the magnitude of BMD loss, with patients carrying the most risk alleles having the greatest loss. CONCLUSION: Our data provide the first evidence that common genetic variants play an important role in BMD loss among HCT patients similar to age-related BMD loss in the general population. This infers that the mechanism for post-HCT bone loss is a normal aging process that is accelerated during HCT. A limitation of our study comes from its small patient population; hence future larger studies are warranted to validate our findings.

Upregulation of hemoglobin expression by oxidative stress in hepatocytes and its implication in nonalcoholic steatohepatitis.

Recent studies revealed that hemoglobin is expressed in some non-erythrocytes and it suppresses oxidative stress when overexpressed. Oxidative stress plays a critical role in the pathogenesis of non-alcoholic steatohepatitis (NASH). This study was designed to investigate whether hemoglobin is expressed in hepatocytes and how it is related to oxidative stress in NASH patients. Analysis of microarray gene expression data revealed a significant increase in the expression of hemoglobin alpha (HBA1) and beta (HBB) in liver biopsies from NASH patients. Increased hemoglobin expression in NASH was validated by quantitative real time PCR. However, the expression of hematopoietic transcriptional factors and erythrocyte specific marker genes were not increased, indicating that increased hemoglobin expression in NASH was not from erythropoiesis, but could result from increased expression in hepatocytes. Immunofluorescence staining demonstrated positive HBA1 and HBB expression in the hepatocytes of NASH livers. Hemoglobin expression was also observed in human hepatocellular carcinoma HepG2 cell line. Furthermore, treatment with hydrogen peroxide, a known oxidative stress inducer, increased HBA1 and HBB expression in HepG2 and HEK293 cells. Importantly, hemoglobin overexpression suppressed oxidative stress in HepG2 cells. We concluded that hemoglobin is expressed by hepatocytes and oxidative stress upregulates its expression. Suppression of oxidative stress by hemoglobin could be a mechanism to protect hepatocytes from oxidative damage in NASH.

Phase 1 study of arsenic trioxide, high-dose cytarabine, and idarubicin to down-regulate constitutive signal transducer and activator of transcription 3 activity in patients aged <60 years with acute myeloid leukemia.

BACKGROUND: Constitutive activation of signal transducer and activator of transcription-3 (STAT3) was detected in blasts from approximately 50% of patients with acute myeloid leukemia (AML) and was correlated with an adverse outcome. In vitro treatment of AML blasts with arsenic trioxide (ATO) down-regulated STAT3 activity within 6 hours associated with a reduced viability within 48 hours. METHODS: A phase 1 clinical trial to evaluate the biologically effective dose and/or the maximally tolerated dose (MTD) of ATO in vivo in conjunction with high-dose cytarabine (Hidac) and idarubicin (Ida) in patients with AML aged <60 years was conducted. Data were compared with 117 historic AML patients who had received treatment with Hidac/Ida. RESULTS: In total, 61 patients were enrolled onto 11 different dose levels (from 0.01 to 0.65 mg/kg ideal body weight). The MTD was 0.5 mg/kg. Compared with historic controls, patients who received ATO/Hidac/Ida, although they had similar pretreatment characteristics, had better overall survival (P = .039). CONCLUSIONS: ATO priming may have improved the outcome of patients aged <60 years with AML who received Hidac/Ida. The current data suggested that ATO may enhance the effect of chemotherapy. The authors concluded that further studies of this novel combination are warranted.