B Cells and T Follicular Helper Cells Mediate Response to Checkpoint Inhibitors in High Mutation Burden Mouse Models of Breast Cancer.

This study identifies mechanisms mediating responses to immune checkpoint inhibitors using mouse models of triple-negative breast cancer. By creating new mammary tumor models, we find that tumor mutation burden and specific immune cells are associated with response. Further, we developed a rich resource of single-cell RNA-seq and bulk mRNA-seq data of immunotherapy-treated and non-treated tumors from sensitive and resistant murine models. Using this, we uncover that immune checkpoint therapy induces T follicular helper cell activation of B cells to facilitate the anti-tumor response in these models. We also show that B cell activation of T cells and the generation of antibody are key to immunotherapy response and propose a new biomarker for immune checkpoint therapy. In total, this work presents resources of new preclinical models of breast cancer with large mRNA-seq and single-cell RNA-seq datasets annotated for sensitivity to therapy and uncovers new components of response to immune checkpoint inhibitors.

Subependymal giant cell astrocytomas are characterized by mTORC1 hyperactivation, a very low somatic mutation rate, and a unique gene expression profile.

Subependymal giant-cell astrocytomas (SEGAs) are slow-growing brain tumors that are a hallmark feature seen in 5-10% of patients with Tuberous Sclerosis Complex (TSC). Though histologically benign, they can cause serious neurologic symptoms, leading to death if untreated. SEGAs consistently show biallelic loss of TSC1 or TSC2. Herein, we aimed to define other somatic events beyond TSC1/TSC2 loss and identify potential transcriptional drivers that contribute to SEGA formation. Paired tumor-normal whole-exome sequencing was performed on 21 resected SEGAs from 20 TSC patients. Pathogenic variants in TSC1/TSC2 were identified in 19/21 (90%) SEGAs. Copy neutral loss of heterozygosity (size range: 2.2-46 Mb) was seen in 76% (16/21) of SEGAs (44% chr9q and 56% chr16p). An average of 1.4 other somatic variants (range 0-7) per tumor were identified, unlikely of pathogenic significance. Whole transcriptome RNA-sequencing analyses revealed 190 common differentially expressed genes in SEGA (n = 16, 13 from a prior study) in pairwise comparison to each of: low grade diffuse gliomas (n = 530) and glioblastoma (n = 171) from The Cancer Genome Atlas (TCGA) consortium, ganglioglioma (n = 10), TSC cortical tubers (n = 15), and multiple normal tissues. Among these, homeobox transcription factors (TFs) HMX3, HMX2, VAX1, SIX3; and TFs IRF6 and EOMES were all expressed >12-fold higher in SEGAs (FDR/q-value < 0.05). Immunohistochemistry supported the specificity of IRF6, VAX1, SIX3 for SEGAs in comparison to other tumor entities and normal brain. We conclude that SEGAs have an extremely low somatic mutation rate, suggesting that TSC1/TSC2 loss is sufficient to drive tumor growth. The unique and highly expressed SEGA-specific TFs likely reflect the neuroepithelial cell of origin, and may also contribute to the transcriptional and epigenetic state that enables SEGA growth following two-hit loss of TSC1 or TSC2 and mTORC1 activation.

Virus expression detection reveals RNA-sequencing contamination in TCGA.

BACKGROUND: Contamination of reagents and cross contamination across samples is a long-recognized issue in molecular biology laboratories. While often innocuous, contamination can lead to inaccurate results. Cantalupo et al., for example, found HeLa-derived human papillomavirus 18 (H-HPV18) in several of The Cancer Genome Atlas (TCGA) RNA-sequencing samples. This work motivated us to assess a greater number of samples and determine the origin of possible contaminations using viral sequences. To detect viruses with high specificity, we developed the publicly available workflow, VirDetect, that detects virus and laboratory vector sequences in RNA-seq samples. We applied VirDetect to 9143 RNA-seq samples sequenced at one TCGA sequencing center (28/33 cancer types) over 5 years. RESULTS: We confirmed that H-HPV18 was present in many samples and determined that viral transcripts from H-HPV18 significantly co-occurred with those from xenotropic mouse leukemia virus-related virus (XMRV). Using laboratory metadata and viral transcription, we determined that the likely contaminant was a pool of cell lines known as the "common reference", which was sequenced alongside TCGA RNA-seq samples as a control to monitor quality across technology transitions (i.e. microarray to GAII to HiSeq), and to link RNA-seq to previous generation microarrays that standardly used the "common reference". One of the cell lines in the pool was a laboratory isolate of MCF-7, which we discovered was infected with XMRV; another constituent of the pool was likely HeLa cells. CONCLUSIONS: Altogether, this indicates a multi-step contamination process. First, MCF-7 was infected with an XMRV. Second, this infected cell line was added to a pool of cell lines, which contained HeLa. Finally, RNA from this pool of cell lines contaminated several TCGA tumor samples most-likely during library construction. Thus, these human tumors with H-HPV or XMRV reads were likely not infected with H-HPV 18 or XMRV.

Assembly-based inference of B-cell receptor repertoires from short read RNA sequencing data with V'DJer.

MOTIVATION: B-cell receptor (BCR) repertoire profiling is an important tool for understanding the biology of diverse immunologic processes. Current methods for analyzing adaptive immune receptor repertoires depend upon PCR amplification of VDJ rearrangements followed by long read amplicon sequencing spanning the VDJ junctions. While this approach has proven to be effective, it is frequently not feasible due to cost or limited sample material. Additionally, there are many existing datasets where short-read RNA sequencing data are available but PCR amplified BCR data are not. RESULTS: We present here V'DJer, an assembly-based method that reconstructs adaptive immune receptor repertoires from short-read RNA sequencing data. This method captures expressed BCR loci from a standard RNA-seq assay. We applied this method to 473 Melanoma samples from The Cancer Genome Atlas and demonstrate V'DJer's ability to accurately reconstruct BCR repertoires from short read mRNA-seq data. AVAILABILITY AND IMPLEMENTATION: V'DJer is implemented in C/C ++, freely available for academic use and can be downloaded from Github: https://github.com/mozack/vdjer CONTACT: benjamin_vincent@med.unc.edu or parkerjs@email.unc.eduSupplementary information: Supplementary data are available at Bioinformatics online.

Post-irradiation intratumoral heterogeneity modulates response to immune checkpoint inhibition therapy in a murine melanoma model.

PURPOSE: The underlying mechanism for radiation as a potentiator of immune checkpoint inhibition (ICI) is unclear. We developed a novel murine model to investigate the effects of post-irradiation intratumoral heterogeneity (ITH) on response to ICI. EXPERIMENTAL DESIGN: Parental mouse melanoma B16F10 cells were irradiated in vitro (5Gy x 3 fractions), then an a priori determined number of resulting colonies were implanted in C57BL/6J immunocompetent mice creating syngeneic models of unirradiated (parental) and irradiated tumors with low (irradiated-L) and high (irradiated-H) ITH. Mice were treated with placebo, α-PD-L1, α-CTLA-4 or dual ICI. Murine tumors underwent whole exome sequencing (WES). Clinically correlated paired pre- and post-irradiation patient rectal adenocarcinoma samples underwent WES. RESULTS: Irradiated-L tumors showed increased tumor mutational burden (TMB) and a sustained decrease in ITH. Irradiated-L tumors were predicted to express five neoantigens with high variant allele frequency/clonal distribution. Mice with irradiated-L and irradiated-H versus parental B16F10 tumors demonstrated longer overall survival with dual ICI. Only mice with irradiated-L tumors experienced an overall survival benefit with single agent ICI. Clinically correlated rectal adenocarcinoma samples showed similarly increased TMB and decreased ITH following irradiation. CONCLUSIONS: Post-irradiation ITH modulates ICI response in a murine melanoma model. Irradiation may offer a mechanism to widen the therapeutic window of ICI.

Integrative Analysis of miRNAs Identifies Clinically Relevant Epithelial and Stromal Subtypes of Head and Neck Squamous Cell Carcinoma.

PURPOSE: The objective of this study is to characterize the role of miRNAs in the classification of head and neck squamous cell carcinoma (HNSCC). EXPERIMENTAL DESIGN: Here, we analyzed 562 HNSCC samples, 88 from a novel cohort and 474 from The Cancer Genome Atlas, using miRNA microarray and miRNA sequencing, respectively. Using an integrative correlations method followed by miRNA expression-based hierarchical clustering, we validated miRNA clusters across cohorts. Evaluation of clusters by logistic regression and gene ontology approaches revealed subtype-based clinical and biological characteristics. RESULTS: We identified two independently validated and statistically significant (P < 0.01) tumor subtypes and named them "epithelial" and "stromal" based on associations with functional target gene ontology relating to differing stages of epithelial cell differentiation. miRNA-based subtypes were correlated with individual gene expression targets based on miRNA seed sequences, as well as with miRNA families and clusters including the miR-17 and miR-200 families. These correlated genes defined pathways relevant to normal squamous cell function and pathophysiology. miRNA clusters statistically associated with differential mutation patterns including higher proportions of TP53 mutations in the stromal class and higher NSD1 and HRAS mutation frequencies in the epithelial class. miRNA classes correlated with previously reported gene expression subtypes, clinical characteristics, and clinical outcomes in a multivariate Cox proportional hazards model with stromal patients demonstrating worse prognoses (HR, 1.5646; P = 0.006). CONCLUSIONS: We report a reproducible classification of HNSCC based on miRNA that associates with known pathologically altered pathways and mutations of squamous tumors and is clinically relevant.

Composition-Dependent Passivation Efficiency at the CdS/CuIn1- x Gax Se2 Interface.

The bandgap of CuIn1- x Gax Se2 (CIGS) chalcopyrite semiconductors can be tuned between ≈1.0 and ≈1.7 eV for Ga contents ranging between x = 0 and x = 1. While an optimum bandgap of 1.34 eV is desirable for achieving maximum solar energy conversion in solar cells, state-of-the-art CIGS-based devices experience a drop in efficiency for Ga contents x > 0.3 (i.e., for bandgaps >1.2 eV), an aspect that is limiting the full potential of these devices. The mechanism underlying the limited performance as a function of CIGS composition has remained elusive: both surface and bulk recombination effects are proposed. Here, the disentanglement between surface and bulk effects in CIGS absorbers as a function of Ga content is achieved by comparing photogenerated charge carrier dynamics in air/CIGS and surface-passivated ZnO/CdS/CIGS samples. While surface passivation prevents surface recombination of charge carriers for low Ga content (x < 0.3; up to 1.2 eV bandgap), surface recombination dominates for higher-bandgap materials. The results thus demonstrate that surface, rather than bulk effects, is responsible for the drop in efficiency for Ga contents larger than x ≈ 0.3.

How machine learning can help select capping layers to suppress perovskite degradation.

Environmental stability of perovskite solar cells (PSCs) has been improved by trial-and-error exploration of thin low-dimensional (LD) perovskite deposited on top of the perovskite absorber, called the capping layer. In this study, a machine-learning framework is presented to optimize this layer. We featurize 21 organic halide salts, apply them as capping layers onto methylammonium lead iodide (MAPbI3) films, age them under accelerated conditions, and determine features governing stability using supervised machine learning and Shapley values. We find that organic molecules' low number of hydrogen-bonding donors and small topological polar surface area correlate with increased MAPbI3 film stability. The top performing organic halide, phenyltriethylammonium iodide (PTEAI), successfully extends the MAPbI3 stability lifetime by 4 ± 2 times over bare MAPbI3 and 1.3 ± 0.3 times over state-of-the-art octylammonium bromide (OABr). Through characterization, we find that this capping layer stabilizes the photoactive layer by changing the surface chemistry and suppressing methylammonium loss.

Author Correction: How machine learning can help select capping layers to suppress perovskite degradation.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

FGFR4 regulates tumor subtype differentiation in luminal breast cancer and metastatic disease.

Mechanisms driving tumor progression from less aggressive subtypes to more aggressive states represent key targets for therapy. We identified a subset of luminal A primary breast tumors that give rise to HER2-enriched (HER2E) subtype metastases, but remain clinically HER2 negative (cHER2-). By testing the unique genetic and transcriptomic features of these cases, we developed the hypothesis that FGFR4 likely participates in this subtype switching. To evaluate this, we developed 2 FGFR4 genomic signatures using a patient-derived xenograft (PDX) model treated with an FGFR4 inhibitor, which inhibited PDX growth in vivo. Bulk tumor gene expression analysis and single-cell RNA sequencing demonstrated that the inhibition of FGFR4 signaling caused molecular switching. In the Molecular Taxonomy of Breast Cancer International Consortium (METABRIC) breast cancer cohort, FGFR4-induced and FGFR4-repressed signatures each predicted overall survival. Additionally, the FGFR4-induced signature was an independent prognostic factor beyond subtype and stage. Supervised analysis of 77 primary tumors with paired metastases revealed that the FGFR4-induced signature was significantly higher in luminal/ER+ tumor metastases compared with their primaries. Finally, multivariate analysis demonstrated that the FGFR4-induced signature also predicted site-specific metastasis for lung, liver, and brain, but not for bone or lymph nodes. These data identify a link between FGFR4-regulated genes and metastasis, suggesting treatment options for FGFR4-positive patients, whose high expression is not caused by mutation or amplification.

Rapid Clearance Profile of Plasma Circulating Tumor HPV Type 16 DNA during Chemoradiotherapy Correlates with Disease Control in HPV-Associated Oropharyngeal Cancer.

PURPOSE: To identify a profile of circulating tumor human papilloma virus (HPV) DNA (ctHPVDNA) clearance kinetics that is associated with disease control after chemoradiotherapy (CRT) for HPV-associated oropharyngeal squamous cell carcinoma (OPSCC). EXPERIMENTAL DESIGN: A multi-institutional prospective biomarker trial was conducted in 103 patients with (i) p16-positive OPSCC, (ii) M0 disease, and (iii) receipt of definitive CRT. Blood specimens were collected at baseline, weekly during CRT, and at follow-up visits. Optimized multianalyte digital PCR assays were used to quantify ctHPVDNA (types 16/18/31/33/35) in plasma. A control cohort of 55 healthy volunteers and 60 patients with non-HPV-associated malignancy was also analyzed. RESULTS: Baseline plasma ctHPVDNA had high specificity (97%) and high sensitivity (89%) for detecting newly diagnosed HPV-associated OPSCC. Pretreatment ctHPV16DNA copy number correlated with disease burden, tumor HPV copy number, and HPV integration status. We define a ctHPV16DNA favorable clearance profile as having high baseline copy number (>200 copies/mL) and >95% clearance of ctHPV16DNA by day 28 of CRT. Nineteen of 67 evaluable patients had a ctHPV16DNA favorable clearance profile, and none had persistent or recurrent regional disease after CRT. In contrast, patients with adverse clinical risk factors (T4 or >10 pack years) and an unfavorable ctHPV16DNA clearance profile had a 35% actuarial rate of persistent or recurrent regional disease after CRT (P = 0.0049). CONCLUSIONS: A rapid clearance profile of ctHPVDNA may predict likelihood of disease control in patients with HPV-associated OPSCC patients treated with definitive CRT and may be useful in selecting patients for deintensified therapy.

Epstein-Barr Virus-Positive Cancers Show Altered B-Cell Clonality.

Epstein-Barr virus (EBV) is convincingly associated with gastric cancer, nasopharyngeal carcinoma, and certain lymphomas, but its role in other cancer types remains controversial. To test the hypothesis that there are additional cancer types with high prevalence of EBV, we determined EBV viral expression in all the Cancer Genome Atlas Project (TCGA) mRNA sequencing (mRNA-seq) samples (n = 10,396) from 32 different tumor types. We found that EBV was present in gastric adenocarcinoma and lymphoma, as expected, and was also present in >5% of samples in 10 additional tumor types. For most samples, EBV transcript levels were low, which suggests that EBV was likely present due to infected infiltrating B cells. In order to determine if there was a difference in the B-cell populations, we assembled B-cell receptors for each sample and found B-cell receptor abundance (P ≤ 1.4 × 10-20) and diversity (P ≤ 8.3 × 10-27) were significantly higher in EBV-positive samples. Moreover, diversity was independent of B-cell abundance, suggesting that the presence of EBV was associated with an increased and altered B-cell population. IMPORTANCE Around 20% of human cancers are associated with viruses. Epstein-Barr virus (EBV) contributes to gastric cancer, nasopharyngeal carcinoma, and certain lymphomas, but its role in other cancer types remains controversial. We assessed the prevalence of EBV in RNA-seq from 32 tumor types in the Cancer Genome Atlas Project (TCGA) and found EBV to be present in >5% of samples in 12 tumor types. EBV infects epithelial cells and B cells and in B cells causes proliferation. We hypothesized that the low expression of EBV in most of the tumor types was due to infiltration of B cells into the tumor. The increase in B-cell abundance and diversity in subjects where EBV was detected in the tumors strengthens this hypothesis. Overall, we found that EBV was associated with an increased and altered immune response. This result is not evidence of causality, but a potential novel biomarker for tumor immune status.

Integrated RNA and DNA sequencing reveals early drivers of metastatic breast cancer.

Breast cancer metastasis remains a clinical challenge, even within a single patient across multiple sites of the disease. Genome-wide comparisons of both the DNA and gene expression of primary tumors and metastases in multiple patients could help elucidate the underlying mechanisms that cause breast cancer metastasis. To address this issue, we performed DNA exome and RNA sequencing of matched primary tumors and multiple metastases from 16 patients, totaling 83 distinct specimens. We identified tumor-specific drivers by integrating known protein-protein network information with RNA expression and somatic DNA alterations and found that genetic drivers were predominantly established in the primary tumor and maintained through metastatic spreading. In addition, our analyses revealed that most genetic drivers were DNA copy number changes, the TP53 mutation was a recurrent founding mutation regardless of subtype, and that multiclonal seeding of metastases was frequent and occurred in multiple subtypes. Genetic drivers unique to metastasis were identified as somatic mutations in the estrogen and androgen receptor genes. These results highlight the complexity of metastatic spreading, be it monoclonal or multiclonal, and suggest that most metastatic drivers are established in the primary tumor, despite the substantial heterogeneity seen in the metastases.

The Prognostic Significance of Low-Frequency Somatic Mutations in Metastatic Cutaneous Melanoma.

Background: Little is known about the prognostic significance of somatically mutated genes in metastatic melanoma (MM). We have employed a combined clinical and bioinformatics approach on tumor samples from cutaneous melanoma (SKCM) as part of The Cancer Genome Atlas project (TCGA) to identify mutated genes with potential clinical relevance. Methods: After limiting our DNA sequencing analysis to MM samples (n = 356) and to the CANCER CENSUS gene list, we filtered out mutations with low functional significance (snpEFF). We performed Cox analysis on 53 genes that were mutated in ≥3% of samples, and had ≥50% difference in incidence of mutations in deceased subjects versus alive subjects. Results: Four genes were potentially prognostic [RAC1, FGFR1, CARD11, CIITA; false discovery rate (FDR) < 0.2]. We identified 18 additional genes (e.g., SPEN, PDGFRB, GNAS, MAP2K1, EGFR, TSC2) that were less likely to have prognostic value (FDR < 0.4). Most somatic mutations in these 22 genes were infrequent (< 10%), associated with high somatic mutation burden, and were evenly distributed across all exons, except for RAC1 and MAP2K1. Mutations in only 9 of these 22 genes were also identified by RNA sequencing in >75% of the samples that exhibited corresponding DNA mutations. The low frequency, UV signature type and RNA expression of the 22 genes in MM samples were confirmed in a separate multi-institution validation cohort (n = 413). An underpowered analysis within a subset of this validation cohort with available patient follow-up (n = 224) showed that somatic mutations in SPEN and RAC1 reached borderline prognostic significance [log-rank favorable (p = 0.09) and adverse (p = 0.07), respectively]. Somatic mutations in SPEN, and to a lesser extent RAC1, were not associated with definite gene copy number or RNA expression alterations. High (>2+) nuclear plus cytoplasmic expression intensity for SPEN was associated with longer melanoma-specific overall survival (OS) compared to lower (≤ 2+) nuclear intensity (p = 0.048). We conclude that expressed somatic mutations in infrequently mutated genes beyond the well-characterized ones (e.g., BRAF, RAS, CDKN2A, PTEN, TP53), such as RAC1 and SPEN, may have prognostic significance in MM.

Genomic analysis of oesophageal squamous-cell carcinoma identifies alcohol drinking-related mutation signature and genomic alterations.

Approximately half of the world's 500,000 new oesophageal squamous-cell carcinoma (ESCC) cases each year occur in China. Here, we show whole-genome sequencing of DNA and RNA in 94 Chinese individuals with ESCC. We identify six mutational signatures (E1-E6), and Signature E4 is unique in ESCC linked to alcohol intake and genetic variants in alcohol-metabolizing enzymes. We discover significantly recurrent mutations in 20 protein-coding genes, 4 long non-coding RNAs and 10 untranslational regions. Functional analyses show six genes that have recurrent copy-number variants in three squamous-cell carcinomas (oesophageal, head and neck and lung) significantly promote cancer cell proliferation, migration and invasion. The most frequently affected genes by structural variation are LRP1B and TTC28. The aberrant cell cycle and PI3K-AKT pathways seem critical in ESCC. These results establish a comprehensive genomic landscape of ESCC and provide potential targets for precision treatment and prevention of the cancer.