Therapeutic implications of transcriptomics in head and neck cancer patient-derived xenografts

There are currently no clinical strategies utilizing tumor gene expression to inform therapeutic selection for patients with head and neck squamous cell carcinoma (HNSCC). One of the challenges in developing predictive biomarkers is the limited characterization of preclinical HNSCC models. Patient-derived xenografts (PDXs) are increasingly recognized as translationally relevant preclinical avatars for human tumors; however, the overall transcriptomic concordance of HNSCC PDXs with primary human HNSCC is understudied, especially in human papillomavirus-associated (HPV+) disease. Here, we characterized 64 HNSCC PDXs (16 HPV+ and 48 HPV-) at the transcriptomic level using RNA-sequencing. The range of human-specific reads per PDX varied from 64.6%-96.5%, with a comparison of the most differentially expressed genes before and after removal of mouse transcripts revealing no significant benefit to filtering out mouse mRNA reads in this cohort. We demonstrate that four previously established HNSCC molecular subtypes found in The Cancer Genome Atlas (TCGA) are also clearly recapitulated in HNSCC PDXs. Unsupervised hierarchical clustering yielded a striking natural division of HNSCC PDXs by HPV status, with C19orf57 (BRME1), a gene previously correlated with positive response to cisplatin in cervical cancer, among the most significantly differentially expressed genes between HPV+ and HPV- PDXs. In vivo experiments demonstrated a possible relationship between increased C19orf57 expression and superior anti-tumor responses of PDXs to cisplatin, which should be investigated further. These findings highlight the value of PDXs as models for HPV+ and HPV- HNSCC, providing a resource for future discovery of predictive biomarkers to guide treatment selection in HNSCC.

Association between Ancestry-Specific 6q25 Variants and Breast Cancer Subtypes in Peruvian Women

BACKGROUND

Breast cancer incidence in the United States is lower in Hispanic/Latina (H/L) compared with African American/Black or Non-Hispanic White women. An Indigenous American breast cancer-protective germline variant (rs140068132) has been reported near the estrogen receptor 1 gene. This study tests the association of rs140068132 and other polymorphisms in the 6q25 region with subtype-specific breast cancer risk in H/Ls of high Indigenous American ancestry.

METHODS

Genotypes were obtained for 5,094 Peruvian women with (1,755) and without (3,337) breast cancer. Associations between genotype and overall and subtype-specific risk for the protective variant were tested using logistic regression models and conditional analyses, including other risk-associated polymorphisms in the region.

RESULTS

We replicated the reported association between rs140068132 and breast cancer risk overall [odds ratio (OR), 0.53; 95% confidence interval (CI), 0.47-0.59], as well as the lower odds of developing hormone receptor negative (HR-) versus HR+ disease (OR, 0.77; 95% CI, 0.61-0.97). Models, including HER2, showed further heterogeneity with reduced odds for HR+HER2+ (OR, 0.68; 95% CI, 0.51-0.92), HR-HER2+ (OR, 0.63; 95% CI, 0.44-0.90) and HR-HER2- (OR, 0.77; 95% CI, 0.56-1.05) compared with HR+HER2-. Inclusion of other risk-associated variants did not change these observations.

CONCLUSIONS

The rs140068132 polymorphism is associated with decreased risk of breast cancer in Peruvians and is more protective against HR- and HER2+ diseases independently of other breast cancer-associated variants in the 6q25 region.

IMPACT

These results could inform functional analyses to understand the mechanism by which rs140068132-G reduces risk of breast cancer development in a subtype-specific manner. They also illustrate the importance of including diverse individuals in genetic studies.

Caveolin-1 and Sox-2 are predictive biomarkers of cetuximab response in head and neck cancer

The epidermal growth factor receptor (EGFR) inhibitor cetuximab is the only FDA-approved oncogene-targeting therapy for head and neck squamous cell carcinoma (HNSCC). Despite variable treatment response, no biomarkers exist to stratify patients for cetuximab therapy in HNSCC. Here, we applied unbiased hierarchical clustering to reverse-phase protein array molecular profiles from patient-derived xenograft (PDX) tumors and revealed 2 PDX clusters defined by protein networks associated with EGFR inhibitor resistance. In vivo validation revealed unbiased clustering to classify PDX tumors according to cetuximab response with 88% accuracy. Next, a support vector machine classifier algorithm identified a minimalist biomarker signature consisting of 8 proteins — caveolin-1, Sox-2, AXL, STING, Brd4, claudin-7, connexin-43, and fibronectin — with expression that strongly predicted cetuximab response in PDXs using either protein or mRNA. A combination of caveolin-1 and Sox-2 protein levels was sufficient to maintain high predictive accuracy, which we validated in tumor samples from patients with HNSCC with known clinical response to cetuximab. These results support further investigation into the combined use of caveolin-1 and Sox-2 as predictive biomarkers for cetuximab response in the clinic.

In utero and early-life exposure to thirdhand smoke causes profound changes to the immune system

Acute lymphoblastic leukemia (ALL) is the most common cancer in children. Thirdhand smoke (THS) is the residual tobacco contamination that remains after the smoke clears. We investigated the effects of THS exposure in utero and during early life in a transgenic Cdkn2a knockout mouse model that is vulnerable to the development of leukemia/lymphoma. Female mice, and their offspring, were exposed from the first day of pregnancy to weaning. Plasma cytokines, body weight and hematologic parameters were measured in the offspring. To investigate THS exposure effects on the development of leukemia/lymphoma, bone marrow (BM) was collected from control and THS-exposed mice and transplanted into BM-ablated recipient mice, which were followed for tumor development for 1 year. We found that in utero and early-life THS exposure caused significant changes in plasma cytokine concentrations and in immune cell populations; changes appeared more pronounced in male mice. Spleen (SP) and BM B-cell populations were significantly lower in THS-exposed mice. We furthermore observed that THS exposure increased the leukemia/lymphoma-free survival in BM transplantation recipient mice, potentially caused by THS-induced B-cell toxicity. A trend towards increased solid tumors in irradiated mice reconstituted with THS-exposed BM stimulates the hypothesis that the immunosuppressive effects of in utero and early-life THS exposure might contribute to carcinogenesis by lowering the host defense to other toxic exposures. Our study adds to expanding evidence that THS exposure alters the immune system and that in utero and early-life developmental periods represent vulnerable windows of susceptibility for these effects.

Molecular and phenotypic diversity of CBL-mutated juvenile myelomonocytic leukemia

Mutations in the CBL gene were first identified in adults with various myeloid malignancies. Some patients with juvenile myelomonocytic leukemia (JMML) were also noted to harbor mutations in CBL, but were found to have generally less aggressive disease courses compared to patients with other forms of Ras pathway-mutant JMML. Importantly, and in contrast to most reports in adults, the majority of CBL mutations in JMML patients are germline with acquired uniparental disomy occurring in affected marrow cells. Here, we systematically studied a large cohort of 33 JMML patients with CBL mutations and found that this disease is highly diverse in presentation and overall outcome. Moreover, we discovered somatically acquired CBL mutations in 15% of pediatric patients who presented with more aggressive disease. Neither clinical features nor methylation profiling were able to distinguish patients with somatic CBL mutations from those with germline CBL mutations, highlighting the need for germline testing. Overall, we demonstrate that disease courses are quite heterogeneous even among patients with germline CBL mutations. Prospective clinical trials are warranted to find ideal treatment strategies for this diverse cohort of patients.

Glucocorticoids paradoxically facilitate steroid resistance in T cell acute lymphoblastic leukemias and thymocytes

Glucocorticoids (GCs) are a central component of therapy for patients with T cell acute lymphoblastic leukemia (T-ALL), and although resistance to GCs is a strong negative prognostic indicator in T-ALL, the mechanisms of GC resistance remain poorly understood. Using diagnostic samples from patients enrolled in the frontline Children's Oncology Group (COG) T-ALL clinical trial AALL1231, we demonstrated that one-third of primary T-ALLs were resistant to GCs when cells were cultured in the presence of IL-7, a cytokine that is critical for normal T cell function and that plays a well-established role in leukemogenesis. We demonstrated that in these T-ALLs and in distinct populations of normal developing thymocytes, GCs paradoxically induced their own resistance by promoting upregulation of IL-7 receptor (IL-7R) expression. In the presence of IL-7, this augmented downstream signal transduction, resulting in increased STAT5 transcriptional output and upregulation of the prosurvival protein BCL-2. Taken together, we showed that IL-7 mediates an intrinsic and physiologic mechanism of GC resistance in normal thymocyte development that is retained during leukemogenesis in a subset of T-ALLs and is reversible with targeted inhibition of the IL-7R/JAK/STAT5/BCL-2 axis.

Developmental roles of 21 Drosophila transcription factors are determined by quantitative differences in binding to an overlapping set of thousands of genomic regions

Distinct developmental fates in Drosophila melanogaster are specified by quantitative differences in transcription factor occupancy on a common set of bound regions.

Background

We previously established that six sequence-specific transcription factors that initiate anterior/posterior patterning in Drosophila bind to overlapping sets of thousands of genomic regions in blastoderm embryos. While regions bound at high levels include known and probable functional targets, more poorly bound regions are preferentially associated with housekeeping genes and/or genes not transcribed in the blastoderm, and are frequently found in protein coding sequences or in less conserved non-coding DNA, suggesting that many are likely non-functional.

Results

Here we show that an additional 15 transcription factors that regulate other aspects of embryo patterning show a similar quantitative continuum of function and binding to thousands of genomic regions in vivo. Collectively, the 21 regulators show a surprisingly high overlap in the regions they bind given that they belong to 11 DNA binding domain families, specify distinct developmental fates, and can act via different cis-regulatory modules. We demonstrate, however, that quantitative differences in relative levels of binding to shared targets correlate with the known biological and transcriptional regulatory specificities of these factors.

Conclusions

It is likely that the overlap in binding of biochemically and functionally unrelated transcription factors arises from the high concentrations of these proteins in nuclei, which, coupled with their broad DNA binding specificities, directs them to regions of open chromatin. We suggest that most animal transcription factors will be found to show a similar broad overlapping pattern of binding in vivo, with specificity achieved by modulating the amount, rather than the identity, of bound factor.

Transcription factors bind thousands of active and inactive regions in the Drosophila blastoderm

Identifying the genomic regions bound by sequence-specific regulatory factors is central both to deciphering the complex DNA cis-regulatory code that controls transcription in metazoans and to determining the range of genes that shape animal morphogenesis. We used whole-genome tiling arrays to map sequences bound in Drosophila melanogaster embryos by the six maternal and gap transcription factors that initiate anterior-posterior patterning. We find that these sequence-specific DNA binding proteins bind with quantitatively different specificities to highly overlapping sets of several thousand genomic regions in blastoderm embryos. Specific high- and moderate-affinity in vitro recognition sequences for each factor are enriched in bound regions. This enrichment, however, is not sufficient to explain the pattern of binding in vivo and varies in a context-dependent manner, demonstrating that higher-order rules must govern targeting of transcription factors. The more highly bound regions include all of the over 40 well-characterized enhancers known to respond to these factors as well as several hundred putative new cis-regulatory modules clustered near developmental regulators and other genes with patterned expression at this stage of embryogenesis. The new targets include most of the microRNAs (miRNAs) transcribed in the blastoderm, as well as all major zygotically transcribed dorsal-ventral patterning genes, whose expression we show to be quantitatively modulated by anterior-posterior factors. In addition to these highly bound regions, there are several thousand regions that are reproducibly bound at lower levels. However, these poorly bound regions are, collectively, far more distant from genes transcribed in the blastoderm than highly bound regions; are preferentially found in protein-coding sequences; and are less conserved than highly bound regions. Together these observations suggest that many of these poorly bound regions are not involved in early-embryonic transcriptional regulation, and a significant proportion may be nonfunctional. Surprisingly, for five of the six factors, their recognition sites are not unambiguously more constrained evolutionarily than the immediate flanking DNA, even in more highly bound and presumably functional regions, indicating that comparative DNA sequence analysis is limited in its ability to identify functional transcription factor targets.