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.

The impact of nucleosome structure on CRISPR/Cas9 fidelity

The clustered regularly interspaced short palindromic repeats (CRISPR) Cas system is a powerful tool that has the potential to become a therapeutic gene editor in the near future. Cas9 is the best studied CRISPR system and has been shown to have problems that restrict its use in therapeutic applications. Chromatin structure is a known impactor of Cas9 targeting and there is a gap in knowledge on Cas9's efficacy when targeting such locations. To quantify at a single base pair resolution how chromatin inhibits on-target gene editing relative to off-target editing of exposed mismatching targets, we developed the gene editor mismatch nucleosome inhibition assay (GEMiNI-seq). GEMiNI-seq utilizes a library of nucleosome sequences to examine all target locations throughout nucleosomes in a single assay. The results from GEMiNI-seq revealed that the location of the protospacer-adjacent motif (PAM) sequence on the nucleosome edge drives the ability for Cas9 to access its target sequence. In addition, Cas9 had a higher affinity for exposed mismatched targets than on-target sequences within a nucleosome. Overall, our results show how chromatin structure impacts the fidelity of Cas9 to potential targets and highlight how targeting sequences with exposed PAMs could limit off-target gene editing, with such considerations improving Cas9 efficacy and resolving current limitations.

Dietary carbohydrate intake is associated with the subgingival plaque oral microbiome abundance and diversity in a cohort of postmenopausal women

Limited research exists on carbohydrate intake and oral microbiome diversity and composition assessed with next-generation sequencing. We aimed to better understand the association between habitual carbohydrate intake and the oral microbiome, as the oral microbiome has been associated with caries, periodontal disease, and systemic diseases. We investigated if total carbohydrates, starch, monosaccharides, disaccharides, fiber, or glycemic load (GL) were associated with the diversity and composition of oral bacteria in subgingival plaque samples of 1204 post-menopausal women. Carbohydrate intake and GL were assessed from a food frequency questionnaire, and adjusted for energy intake. The V3-V4 region of the 16S rRNA gene from subgingival plaque samples were sequenced to identify the relative abundance of microbiome compositional data expressed as operational taxonomic units (OTUs). The abundance of OTUs were centered log(2)-ratio transformed to account for the compositional data structure. Associations between carbohydrate/GL intake and microbiome alpha-diversity measures were examined using linear regression. PERMANOVA analyses were conducted to examine microbiome beta-diversity measures across quartiles of carbohydrate/GL intake. Associations between intake of carbohydrates and GL and the abundance of the 245 identified OTUs were examined by using linear regression. Total carbohydrates, GL, starch, lactose, and sucrose intake were inversely associated with alpha-diversity measures. Beta-diversity across quartiles of total carbohydrates, fiber, GL, sucrose, and galactose, were all statistically significant (p for PERMANOVA p < 0.05). Positive associations were observed between total carbohydrates, GL, sucrose and Streptococcus mutans; GL and both Sphingomonas HOT 006 and Scardovia wiggsiae; and sucrose and Streptococcus lactarius. A negative association was observed between lactose and Aggregatibacter segnis, and between sucrose and both TM7_[G-1] HOT 346 and Leptotrichia HOT 223. Intake of total carbohydrate, GL, and sucrose were inversely associated with subgingival bacteria alpha-diversity, the microbial beta-diversity varied by their intake, and they were associated with the relative abundance of specific OTUs. Higher intake of sucrose, or high GL foods, may influence poor oral health outcomes (and perhaps systemic health outcomes) in older women via their influence on the oral microbiome.

The Subgingival Microbiome Relationship to Periodontal Disease in Older Women

Understanding of the oral microbiome in relation to periodontal disease in older adults is limited. The composition and diversity of the subgingival microflora and their oligotypes in health and levels of periodontal disease were investigated in this study on older postmenopausal women. The 16S rRNA gene was sequenced using the Illumina MiSeq platform in 1,206 women aged 53 to 81 y. Presence and severity of periodontal disease were defined by Centers for Disease Control and Prevention/American Academy of Periodontology criteria. Composition of the microbiome was determined by 16S rRNA amplicon sequencing and the abundance of taxa described by the centered log₂-ratio (CLR) transformed operational taxonomic unit (OTU) values. Differences according to periodontal disease status were determined by analysis of variance with Bonferroni correction. Bacteria oligotypes associated with periodontal disease and health were determined by minimum entropy decomposition and their functions estimated in silico using PICRUSt. Prevalence of none/mild, moderate, and severe periodontal disease was 25.1%, 58.3%, and 16.6%, respectively. Alpha diversity of the microbiome differed significantly across the 3 periodontal disease categories. β-Diversity differed between no/mild and severe periodontal disease, although considerable overlap was noted. Of the 267 bacterial species identified at ≥0.02% abundance, 56 (20.9%) differed significantly in abundance according to periodontal disease status. Significant linear correlations for pocket depth and clinical attachment level with bacterial amounts were observed for several taxa. Of the taxa differing in abundance according to periodontal disease status, 53% had multiple oligotypes appearing to differ between none/mild and severe periodontal disease. Among older women, taxonomic differences in subgingival microbiome composition and diversity were observed in relation to clinical periodontal disease measures. Potential differences in bacterial subspecies (oligotypes) and their function were also identified in periodontal disease compared with health.

Composition and diversity of the subgingival microbiome and its relationship with age in postmenopausal women: an epidemiologic investigation

BACKGROUND

The extent to which the composition and diversity of the oral microbiome varies with age is not clearly understood.

METHODS

The 16S rRNA gene of subgingival plaque in 1219 women, aged 53-81 years, was sequenced and its taxonomy annotated against the Human Oral Microbiome Database (v.14.5). Composition of the subgingival microbiome was described in terms of centered log(2)-ratio (CLR) transformed OTU values, relative abundance, and prevalence. Correlations between microbiota abundance and age were evelauted using Pearson Product Moment correlations. P-values were corrected for multiple testing using the Bonferroni method.

RESULTS

Of the 267 species identified overall, Veillonella dispar was the most abundant bacteria when described by CLR OTU (mean 8.3) or relative abundance (mean 8.9%); whereas Streptococcus oralis, Veillonella dispar and Veillonella parvula were most prevalent (100%, all) when described as being present at any amount. Linear correlations between age and several CLR OTUs (Pearson r = - 0.18 to 0.18), of which 82 (31%) achieved statistical significance (P < 0.05). The correlations lost significance following Bonferroni correction. Twelve species that differed across age groups (each corrected P < 0.05); 5 (42%) were higher in women ages 50-59 compared to ≥70 (corrected P < 0.05), and 7 (48%) were higher in women 70 years and older.

CONCLUSIONS

We identified associations between several bacterial species and age across the age range of postmenopausal women studied. Understanding the functions of these bacteria could identify intervention targets to enhance oral health in later life.

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 ETS1^high 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.

The expression of the transcriptional regulator, E26 transformation-specific 1 (ETS1), is elevated in many epithelial cancers and portends aggressive tumor behavior and poor survival. Within these carcinomas, ETS1 function has been shown to be associated with a wide range of cellular responses that include increased proliferation, angiogenesis, metastasis and drug resistance. Here we focus on Head and Neck Squamous Cell Carcinoma (HNSCC) and discover that higher expression of ETS1 is specifically more pronounced in the mesenchymal subtypes of HNSCC, which represent tumors with enriched expression of Epithelial to Mesenchymal Transition (EMT) markers and inflammation. By using genomic and epigenomic strategies, we have identified the global targets of ETS1 in a preclinical Mesenchymal HNSCC cell model and determined the crucial gene network that is most dependent upon its function. We further validate this ETS1-driven gene expression signature within several HNSCC patient derived datasets and conclude that ETS1 acts as a crucial regulator of the TGFβ signaling cascade to drive EMT. Our findings reinforce the challenges of epithelial tumor heterogeneity and offer insights into molecular underpinning of a specific subtype that can be mined for cancer vulnerability.

Is the Oral Microbiome Associated with Blood Pressure in Older Women?

INTRODUCTION

A possible role of the oral microbiome, specifically oral nitrate reducing flora, in blood pressure (BP) homeostasis, if proven etiologic in nature, could lead to novel mechanism-based therapy to improve hypertension prevention and control.

AIM

This cross-sectional study characterized and compared the oral microbiome between four study groups based on BP status among 446 postmenopausal women aged 53-82 years.

METHODS

Three study groups were not taking hypertension medication and were separated based on BP, as follows: normal BP (systolic < 120 and diastolic < 80; N = 179), elevated BP/Stage I hypertension (systolic 120-139 or diastolic 80-90; N = 106), Stage II hypertension (systolic > 140 or diastolic > 90; N = 42). The forth group consisted of anyone taking hypertension medications, regardless of BP (N = 119). Subgingival microbiome composition was determined using 16S rRNA sequencing with the Illumina MiSeq platform. Kruskal-Wallis tests were used to compare species-level relative abundance of bacterial operational taxonomic units across the four groups.

RESULTS

Sixty-five bacterial species demonstrated significant differences in relative abundance in women with elevated BP or using hypertension medication as compared to those with normal BP. After correction for multiple testing, two species, Prevotella oral (species 317) and Streptococcus oralis, remained significant and were lower in abundance among women taking antihypertension medications compared to those with normal BP (corrected P < 0.05).

CONCLUSIONS

These data provide novel description of oral subgingival bacteria grouped according to BP status. Additional larger studies including functional analysis and prospective designs will help further assess the potential role of the oral microbiome in BP regulation and hypertension.

Cohort profile: the Buffalo OsteoPerio microbiome prospective cohort study

PURPOSE

The Buffalo Osteoporosis and Periodontal Disease (OsteoPerio) study is a prospective cohort study focused on the relationship between the microbiome and oral and systemic health outcomes in postmenopausal women. The cohort was established to examine how the oral microbiome is affected by (and how it affects) periodontal disease presence, severity and progression and to characterise the relationship between the microbiome, lifestyle habits and systemic disease outcomes.

PARTICIPANTS

Participants (n=1342) were postmenopausal women who were participating in the Women's Health Initiative observational study at the Buffalo, New York clinical centre. There were 1026 participants at the 5-year follow-up visit and 518 at the 15-year visit.

FINDINGS TO DATE

Data collected include questionnaires, anthropometric measures, serum blood and saliva samples. At each clinic visit, participants completed a comprehensive oral examination to measure oral health and the oral microbiome. Preliminary findings have contributed to our understanding of risk factors for periodontal disease and the relationship between the oral microbiome and periodontal disease.

FUTURE PLANS

The novel microbiome data collected on a large sample of participants at three time points will be used to answer a variety of research questions focused on temporal changes in the microbiome and the relationship between the oral microbiome and oral and systemic disease outcomes. Little is currently known about the relationship between the oral microbiome and health outcomes in older adults; data from the OsteoPerio cohort will fill this gap. Microbiome samples are currently being analysed using next-generation sequencing technology with an anticipated completion date of late 2018.

Suppressed hepatic bile acid signalling despite elevated production of primary and secondary bile acids in NAFLD

OBJECTIVE

Bile acids are regulators of lipid and glucose metabolism, and modulate inflammation in the liver and other tissues. Primary bile acids such as cholic acid and chenodeoxycholic acid (CDCA) are produced in the liver, and converted into secondary bile acids such as deoxycholic acid (DCA) and lithocholic acid by gut microbiota. Here we investigated the possible roles of bile acids in non-alcoholic fatty liver disease (NAFLD) pathogenesis and the impact of the gut microbiome on bile acid signalling in NAFLD.

DESIGN

Serum bile acid levels and fibroblast growth factor 19 (FGF19), liver gene expression profiles and gut microbiome compositions were determined in patients with NAFLD, high-fat diet-fed rats and their controls.

RESULTS

Serum concentrations of primary and secondary bile acids were increased in patients with NAFLD. In per cent, the farnesoid X receptor (FXR) antagonistic DCA was increased, while the agonistic CDCA was decreased in NAFLD. Increased mRNA expression for cytochrome P450 7A1, Na+-taurocholate cotransporting polypeptide and paraoxonase 1, no change in mRNA expression for small heterodimer partner and bile salt export pump, and reduced serum FGF19 were evidence of impaired FXR and fibroblast growth factor receptor 4 (FGFR4)-mediated signalling in NAFLD. Taurine and glycine metabolising bacteria were increased in the gut of patients with NAFLD, reflecting increased secondary bile acid production. Similar changes in liver gene expression and the gut microbiome were observed in high-fat diet-fed rats.

CONCLUSIONS

The serum bile acid profile, the hepatic gene expression pattern and the gut microbiome composition consistently support an elevated bile acid production in NAFLD. The increased proportion of FXR antagonistic bile acid explains, at least in part, the suppression of hepatic FXR-mediated and FGFR4-mediated signalling. Our study suggests that future NAFLD intervention may target the components of FXR signalling, including the bile acid converting gut microbiome.

Multi-targeting therapeutic mechanisms of the Chinese herbal medicine QHD in the treatment of non-alcoholic fatty liver disease

Beneficial effects of the Chinese herbal medicine Qushi Huayu Decoction (QHD) were observed with non-alcoholic fatty liver disease (NAFLD) patients and animal models. The impact of QHD or its active components (geniposide and chlorogenic acid, GC) on NAFLD liver transcriptome and gut microbiota was examined with NAFLD rats. Increased expression for genes required for glutathione production and decreased expression for genes required for lipid synthesis was observed in NAFLD livers treated with QHD and GC. GC treatment decreased serum LPS, which could be explained by reduced mucosal damage in the colon of GC-treated rats. Further, our data suggest an increased abundance of Treg-inducing bacteria that stimulated the Treg activity in GC treated colon, which in turn down-regulated inflammatory signals, improved gut barrier function and consequently reduced hepatic exposure to microbial products. Our study suggests that QHD simultaneously enhanced the hepatic anti-oxidative mechanism, decreased hepatic lipid synthesis, and promoted the regulatory T cell inducing microbiota in the gut.

Gut microbiome may contribute to insulin resistance and systemic inflammation in obese rodents: a meta-analysis

A number of studies have associated obesity with altered gut microbiota, although results are discordant regarding compositional changes in the gut microbiota of obese animals. Herein we used a meta-analysis to obtain an unbiased evaluation of structural and functional changes of the gut microbiota in diet-induced obese rodents. The raw sequencing data of nine studies generated from high-fat diet (HFD)-induced obese rodent models were processed with QIIME to obtain gut microbiota compositions. Biological functions were predicted and annotated with KEGG pathways with PICRUSt. No significant difference was observed for alpha diversity and Bacteroidetes-to-Firmicutes ratio between obese and lean rodents. Bacteroidia, Clostridia, Bacilli, and Erysipelotrichi were dominant classes, but gut microbiota compositions varied among studies. Meta-analysis of the nine microbiome data sets identified 15 differential taxa and 57 differential pathways between obese and lean rodents. In obese rodents, increased abundance was observed for Dorea, Oscillospira, and Ruminococcus, known for fermenting polysaccharide into short chain fatty acids (SCFAs). Decreased Turicibacter and increased Lactococcus are consistent with elevated inflammation in the obese status. Differential functional pathways of the gut microbiome in obese rodents included enriched pyruvate metabolism, butanoate metabolism, propanoate metabolism, pentose phosphate pathway, fatty acid biosynthesis, and glycerolipid metabolism pathways. These pathways converge in the function of carbohydrate metabolism, SCFA metabolism, and biosynthesis of lipid. HFD-induced obesity results in structural and functional dysbiosis of gut microbiota. The altered gut microbiome may contribute to obesity development by promoting insulin resistance and systemic inflammation.

Abstract 07: An oncogenic transcriptional network anchored by ETS1, p63 and subtype specific drivers of HNSCC revealed by epigenomic and genomic interrogation

Head and neck squamous cell carcinoma (HNSCC) is a heterogeneous disease of significant mortality and morbidity, with a 5-year survival rate that has shown only modest improvement. This persistent patient mortality and morbidity can be attributed, in part, to an overall lack of novel targeted treatment options. Indeed, recent genome sequencing studies have highlighted the striking heterogeneity of HNSCC tumors and have failed to identify pervasive and broadly actionable driver mutations and/or other genetic alterations. Despite this genetic heterogeneity, HNSCC tumors can be classified into gene-expression defined subtypes, raising the possibility that such distinct gene expression programs and regulatory networks are closely interwoven with an intrinsic heritable and signature epigenetic state. We hypothesized that there exist key oncogenic transcriptional regulators of subtype dependency and that defining the epigenetic landscape of HNSCC tumors will reveal the identity of these crucial factors. This will be impactful since targeting of tumor addiction to key oncogenic Transcription Factors (TFs), although successful in a number of other tumor types, remains largely unexplored in HNSCC due to limited knowledge of the identity of these regulators in this cancer.

To decipher the genomic and epigenomic drivers of HNSCC tumor subtypes, we examined a large cohort of cell lines that represent the Atypical, Basal and Mesenchymal subtypes. By using global H3K27 acetylation levels as a surrogate marker, we first generated a consensus genome-wide regulatory enhancer and Super-Enhancer map. Interestingly, we observed a differential enrichment pattern of H3K27Ac signals that recapitulated the segregation of cell-lines according to subtype based on differential gene expression profiles. We leveraged the differential enhancer map to identify enriched DNA-Binding Motifs to generate a Transcription Factor Combinatorial Regulatory Network, based on the co-occurrence frequency of TFs relative to genomic background. Our analysis highlighted a Cis-Regulatory Module (CRM) Network, specific to each molecular subtype of HNSCC. Furthermore, we identified key TF motifs that are embedded within these CRMs and enriched for each subtype-specific network, which we posit represent the pertinent TF drivers of the distinct gene-expression network found in various HNSCC tumors. Here, as proof of principle, we have examined the potential interaction between TP63 and ETS1, two oncogenic TFs revealed by our epigenomic and genomic analysis to be crucial mediators of HPV-negative HNSCC tumors.

To better understand the TP63-ETS1 driven molecular processes in HPV-negative tumors, we performed transcriptomic analysis of the effects of knockdown of TP63 and ETS1 by RNA-Seq and identified their global targets by ChIP-Seq. Interestingly, our comprehensive analysis showed that TP63 and ETS1 share a significant number of genomic binding sites, exhibit a marked preference for binding to Super-Enhancers and together act in a concerted fashion to drive expression of unique set of pathways and processes. By integrating our findings with meta-analysis of HNSCC patient datasets, we have uncovered a core TP63-ETS1 driven signature network that serves as a regulatory anchor of common facets of cancer biology, such as proliferation and angiogenesis, as well as drug-targetable tumor-specific processes involving cancer stem cells. Our study has harnessed the power of an integrative genomic/epigenomic approach to characterize HNSCC subtypes and heterogeneity and to provide a robust framework for the discovery of novel tumor drivers that can be leveraged for targeted therapeutics.

Citation Format: Christian Gluck, Isha Sethi, Maria Tsompana, Satrajit Sinha. An oncogenic transcriptional network anchored by ETS1, p63 and subtype specific drivers of HNSCC revealed by epigenomic and genomic interrogation [abstract]. In: Proceedings of the AACR-AHNS Head and Neck Cancer Conference: Optimizing Survival and Quality of Life through Basic, Clinical, and Translational Research; April 23-25, 2017; San Diego, CA. Philadelphia (PA): AACR; Clin Cancer Res 2017;23(23_Suppl):Abstract nr 07.

A time-series analysis of altered histone H3 acetylation and gene expression during the course of MMAIII-induced malignant transformation of urinary bladder cells

Our previous studies have shown that chronic exposure to low doses of monomethylarsonous acid (MMAIII) causes global histone acetylation dysregulation in urothelial cells (UROtsa cells) during the course of malignant transformation. To reveal the relationship between altered histone acetylation patterns and aberrant gene expression, more specifically, the carcinogenic relevance of these alterations, we performed a time-course analysis of the binding patterns of histone 3 lysine 18 acetylation (H3K18ac) across the genome and generated global gene-expression profiles from this UROtsa cell malignant transformation model. We showed that H3K18ac, one of the most significantly upregulated histone acetylation sites following MMAIII exposure, was enriched at gene promoter-specific regions across the genome and that MMAIII-induced upregulation of H3K18ac led to an altered binding pattern in a large number of genes that was most significant during the critical window for MMAIII-induced UROtsa cells' malignant transformation. Some genes identified as having a differential binding pattern with H3K18ac, acted as upstream regulators of critical gene networks with known functions in tumor development and progression. The altered H3K18ac binding patterns not only led to changes in expression of these directly affected upstream regulators but also resulted in gene-expression changes in their regulated networks. Collectively, our data suggest that MMAIII-induced alteration of histone acetylation patterns in UROtsa cells led to a time- and malignant stage-dependent aberrant gene-expression pattern, and that some gene regulatory networks were altered in accordance with their roles in carcinogenesis, probably contributing to MMAIII-induced urothelial cell malignant transformation and carcinogenesis.

Salivary inflammatory markers and microbiome in normoglycemic lean and obese children compared to obese children with type 2 diabetes

BACKGROUND

There is emerging evidence linking diabetes with periodontal disease. Diabetes is a well-recognized risk factor for periodontal disease. Conversely, pro-inflammatory molecules released by periodontally-diseased tissues may enter the circulation to induce insulin resistance. While this association has been demonstrated in adults, there is little information regarding periodontal status in obese children with and without type 2 diabetes (T2D). We hypothesized that children with T2D have higher rates of gingivitis, elevated salivary inflammatory markers, and an altered salivary microbiome compared to children without T2D.

METHODS

Three pediatric cohorts ages 10-19 years were studied: lean (normal weight-C), obese (Ob), and obese with T2D (T2D). Each subject completed an oral health survey, received a clinical oral examination, and provided unstimulated saliva for measurement of inflammatory markers and microbiome analysis.

RESULTS

The diabetes group was less likely to have had a dental visit within the last six months. Body mass index (BMI) Z-scores and waist circumference/height ratios were similar between Ob and T2D cohorts. The number of carious lesions and fillings were similar for all three groups. The gingival index was greater in the T2D group compared to the Ob and C groups. Although salivary microbial diversity was minimal between groups, a few differences in bacterial genus composition were noted.

CONCLUSIONS

Obese children with T2D show a trend toward poorer oral health compared to normal weight and obese children without T2D. This study characterizes the salivary microbiome of children with and without obesity and T2D. This study supports a modest link between T2D and periodontal inflammation in the pediatric population.

Connection Map for Compounds (CMC): A Server for Combinatorial Drug Toxicity and Efficacy Analysis

Drug discovery and development is a costly and time-consuming process with a high risk for failure resulting primarily from a drug's associated clinical safety and efficacy potential. Identifying and eliminating inapt candidate drugs as early as possible is an effective way for reducing unnecessary costs, but limited analytical tools are currently available for this purpose. Recent growth in the area of toxicogenomics and pharmacogenomics has provided with a vast amount of drug expression microarray data. Web servers such as CMap and LTMap have used this information to evaluate drug toxicity and mechanisms of action independently; however, their wider applicability has been limited by the lack of a combinatorial drug-safety type of analysis. Using available genome-wide drug transcriptional expression profiles, we developed the first web server for combinatorial evaluation of toxicity and efficacy of candidate drugs named "Connection Map for Compounds" (CMC). Using CMC, researchers can initially compare their query drug gene signatures with prebuilt gene profiles generated from two large-scale toxicogenomics databases, and subsequently perform a drug efficacy analysis for identification of known mechanisms of drug action or generation of new predictions. CMC provides a novel approach for drug repositioning and early evaluation in drug discovery with its unique combination of toxicity and efficacy analyses, expansibility of data and algorithms, and customization of reference gene profiles. CMC can be freely accessed at http://cadd.tongji.edu.cn/webserver/CMCbp.jsp .

Abstract B13: Genomic mistargeting of p63 drives the cancer phenotype in head and neck squamous cell carcinoma

TP63 is a member of the p53 family and is a master regulator of epithelial development and differentiation. Overexpression and/or genomic amplification of TP63 is commonly found in a large number of epithelial tumors, including the majority of head and neck squamous cell carcinomas (HNSCC). In many tumors, p63 appears to function as an oncogene, contributing to both increased tumor cell proliferation and survival. To address the mechanisms of p63 function, we studied the interplay of p63 with the chromatin microenvironment at its genomic targets, under distinct physiological and pathological states.

We generated transcriptomic (RNA-Seq) and epigenomic (ChIP-Seq, ATAC-Seq, and DNase-Seq) datasets in representative normal oral and epidermal keratinocytes (HGEP, NHEK) and HNSCC (SCC25, SCC13) cell-lines, to examine whether there exists a HNSCC-specific molecular pattern of p63 targeting. We utilized machine learning methods including Random Forests for classification and regression, to determine the contribution of epigenetic states, transcriptional co-factors and the p63 DNA binding motif in differential p63 binding between normal and cancer states. We found that a core epigenomic signature (H3K4me1 and H3K27me3 binding status and chromatin accessibility) is highly predictive of p63 binding. Furthermore, our analysis revealed that SCC25-specific p63 targets (corresponding to 1239 genomic regions) are inaccessible in normal cells, but in SCC25 they are marked by H3K4me1 and H3K27ac marks, indicative of active regulatory regions. Our data suggests that the complex interplay between local chromatin microenvironment, presence and/or absence of distinct co-factors and potential cell-type specific p63 activity might be key drivers of the differential binding events between normal and cancer epithelial cells.

Our investigation of SCC-specific p63 binding events also led to the identification of several oncogenic genes (for e.g. SOX2, TRIM32, HAS2, PRKCE) that are direct targets of p63. By performing Gene Set Enrichment Analysis we found that these p63 misregulated genes encompass both known (G-protein coupled receptor signaling pathway) and novel cancer pathways thus providing a framework for the oncogenic function of p63. To further evaluate the importance of these findings in the tumor context, we derived a p63 target genes signature and used it to cluster RNA-Seq data of HNSCC patients (obtained from The Cancer Genome Atlas). The gene signature could both classify normal samples away from tumors and importantly distinguish between different molecular subtypes of HNSCCs.

Collectively our studies suggest that aberrant targeting of p63 leads to dynamic remodeling of the chromatin at its target sites, trigger the expression of a cascade of pro-tumorigenic genes, ultimately condemning the cell to an oncogenic fate. These novel insights potentially provide new avenues to be harnessed for tumor subtype classification, development of biomarkers and for targeted therapies.

Citation Format: Isha Sethi, Maria Tsompana, Christian Gluck, Michael J. Buck, Satrajit Sinha. Genomic mistargeting of p63 drives the cancer phenotype in head and neck squamous cell carcinoma. [abstract]. In: Proceedings of the AACR Special Conference on Chromatin and Epigenetics in Cancer; Sep 24-27, 2015; Atlanta, GA. Philadelphia (PA): AACR; Cancer Res 2016;76(2 Suppl):Abstract nr B13.

An accurate and efficient experimental approach for characterization of the complex oral microbiota

BACKGROUND

Currently, taxonomic interrogation of microbiota is based on amplification of 16S rRNA gene sequences in clinical and scientific settings. Accurate evaluation of the microbiota depends heavily on the primers used, and genus/species resolution bias can arise with amplification of non-representative genomic regions. The latest Illumina MiSeq sequencing chemistry has extended the read length to 300 bp, enabling deep profiling of large number of samples in a single paired-end reaction at a fraction of the cost. An increasingly large number of researchers have adopted this technology for various microbiome studies targeting the 16S rRNA V3-V4 hypervariable region.

RESULTS

To expand the applicability of this powerful platform for further descriptive and functional microbiome studies, we standardized and tested an efficient, reliable, and straightforward workflow for the amplification, library construction, and sequencing of the 16S V1-V3 hypervariable region using the new 2 × 300 MiSeq platform. Our analysis involved 11 subgingival plaque samples from diabetic and non-diabetic human subjects suffering from periodontitis. The efficiency and reliability of our experimental protocol was compared to 16S V3-V4 sequencing data from the same samples. Comparisons were based on measures of observed taxonomic richness and species evenness, along with Procrustes analyses using beta(β)-diversity distance metrics. As an experimental control, we also analyzed a total of eight technical replicates for the V1-V3 and V3-V4 regions from a synthetic community with known bacterial species operon counts. We show that our experimental protocol accurately measures true bacterial community composition. Procrustes analyses based on unweighted UniFrac β-diversity metrics depicted significant correlation between oral bacterial composition for the V1-V3 and V3-V4 regions. However, measures of phylotype richness were higher for the V1-V3 region, suggesting that V1-V3 offers a deeper assessment of population diversity and community ecology for the complex oral microbiota.

CONCLUSION

This study provides researchers with valuable experimental evidence for the selection of appropriate 16S amplicons for future human oral microbiome studies. We expect that the tested 16S V1-V3 framework will be widely applicable to other types of microbiota, allowing robust, time-efficient, and inexpensive examination of thousands of samples for population, phylogenetic, and functional crossectional and longitutidal studies.

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.

Sunitinib dose escalation overcomes transient resistance in clear cell renal cell carcinoma and is associated with epigenetic modifications

Sunitinib is considered a first-line therapeutic option for patients with advanced clear cell renal cell carcinoma (ccRCC). Despite sunitinib's clinical efficacy, patients eventually develop drug resistance and disease progression. Herein, we tested the hypothesis whether initial sunitinib resistance may be transient and could be overcome by dose increase. In selected patients initially treated with 50 mg sunitinib and presenting with minimal toxicities, sunitinib dose was escalated to 62.5 mg and/or 75 mg at the time of tumor progression. Mice bearing two different patient-derived ccRCC xenografts (PDX) were treated 5 days per week with a dose-escalation schema (40-60-80 mg/kg sunitinib). Tumor tissues were collected before dose increments for immunohistochemistry analyses and drug levels. Selected intrapatient sunitinib dose escalation was safe and several patients had added progression-free survival. In parallel, our preclinical results showed that PDXs, although initially responsive to sunitinib at 40 mg/kg, eventually developed resistance. When the dose was incrementally increased, again we observed tumor response to sunitinib. A resistant phenotype was associated with transient increase of tumor vasculature despite intratumor sunitinib accumulation at higher dose. In addition, we observed associated changes in the expression of the methyltransferase EZH2 and histone marks at the time of resistance. Furthermore, specific EZH2 inhibition resulted in increased in vitro antitumor effect of sunitinib. Overall, our results suggest that initial sunitinib-induced resistance may be overcome, in part, by increasing the dose, and highlight the potential role of epigenetic changes associated with sunitinib resistance that can represent new targets for therapeutic intervention.

Chromatin accessibility: a window into the genome

Transcriptional activation throughout the eukaryotic lineage has been tightly linked with disruption of nucleosome organization at promoters, enhancers, silencers, insulators and locus control regions due to transcription factor binding. Regulatory DNA thus coincides with open or accessible genomic sites of remodeled chromatin. Current chromatin accessibility assays are used to separate the genome by enzymatic or chemical means and isolate either the accessible or protected locations. The isolated DNA is then quantified using a next-generation sequencing platform. Wide application of these assays has recently focused on the identification of the instrumental epigenetic changes responsible for differential gene expression, cell proliferation, functional diversification and disease development. Here we discuss the limitations and advantages of current genome-wide chromatin accessibility assays with especial attention on experimental precautions and sequence data analysis. We conclude with our perspective on future improvements necessary for moving the field of chromatin profiling forward.

The molecular population genetics of the Tomato spotted wilt virus (TSWV) genome

RNA viruses are characterized by high genetic variability resulting in rapid adaptation to new or resistant hosts. Research for plant RNA virus genetic structure and its variability has been relatively scarce compared to abundant research done for human and animal RNA viruses. Here, we utilized a molecular population genetic framework to characterize the evolution of a highly pathogenic plant RNA virus [Tomato spotted wilt virus (TSWV), Tospovirus, Bunyaviridae]. Data from genes encoding five viral proteins were used for phylogenetic analysis, and for estimation of population parameters, subpopulation differentiation, recombination, divergence between Tospovirus species, and selective constraints on the TSWV genome. Our analysis has defined the geographical structure of TSWV, attributed possibly to founder effects. Also, we identify positive selection favouring divergence between Tospovirus species. At the species level, purifying selection has acted to preserve protein function, although certain amino acids appear to be under positive selection. This analysis provides demonstration of population structuring and species-wide population expansions in a multisegmented plant RNA virus, using sequence-based molecular population genetic analyses. It also identifies specific amino acid sites subject to selection within Bunyaviridae and estimates the level of genetic heterogeneity of a highly pathogenic plant RNA virus. The study of the variability of TSWV populations lays the foundation in the development of strategies for the control of other viral diseases in floral crops.