Genomic insights into triple-negative and HER2-positive breast cancers using isogenic model systems

Introduction In general, genomic signatures of breast cancer subtypes have little or no overlap owing to the heterogeneous genetic backgrounds of study samples. Thus, obtaining a reliable signature in the context of isogenic nature of the cells has been challenging and the precise contribution of isogenic triple negative breast cancer (TNBC) versus non-TNBC remains poorly defined. Methods We established isogenic stable cell lines representing TNBC and Human Epidermal Growth Factor Receptor 2 positive (HER2+) breast cancers by introducing HER2 in TNBC cell lines MDA-MB-231 and MDA-MB-468. We examined protein level expression and functionality of the transfected receptor by treatment with an antagonist of HER2. Using microarray profiling, we obtained a comprehensive gene list of differentially expressed between TNBC and HER2+ clones. We identified and validated underlying isogenic components using qPCR and also compared results with expression data from patients with similar breast cancer subtypes. Results We identified 544 and 1087 statistically significant differentially expressed genes between isogenic TNBC and HER2+ samples in MDA-MB-231 and MDA-MB-468 backgrounds respectively and a shared signature of 49 genes. By comparing results from MDA-MB-231 and MDA-MB-468 backgrounds with two patient microarray datasets, we identified 17 and 22 common genes with same expression trend respectively. Additionally, we identified 56 and 78 genes from MDA-MB-231 and MDA-MB-468 comparisons respectively present in our published RNA-seq data. Conclusions Using our unique model system, we have identified an isogenic gene expression signature between TNBC and HER2+ breast cancer. A portion of our results was also verified in patient data samples, indicating an existence of isogenic element associated with HER2 status between genetically heterogeneous breast cancer samples. These findings may potentially contribute to the development of molecular platform that would be valuable for diagnostic and therapeutic decision for TNBC and in distinguishing it from HER2+ subtype

Sex-dependent transcription of cardiac electrophysiology and links to acetylation modifiers based on the GTEx database

Development of safer drugs based on epigenetic modifiers, e.g., histone deacetylase inhibitors (HDACi), requires better understanding of their effects on cardiac electrophysiology. Using RNAseq data from the genotype-tissue-expression database (GTEx), we created models that link the abundance of acetylation enzymes (HDAC/SIRT/HATs), and the gene expression of ion channels (IC) via select cardiac transcription factors (TFs) in male and female adult human hearts (left ventricle, LV). Gene expression data (transcripts per million, TPM) from GTEx donors (21–70 y.o.) were filtered, normalized and transformed to Euclidian space to allow quantitative comparisons in 84 female and 158 male LVs. Sex-specific partial least-square (PLS) regression models, linking gene expression data for HDAC/SIRT/HATs to TFs and to ICs gene expression, revealed tight co-regulation of cardiac ion channels by HDAC/SIRT/HATs, with stronger clustering in the male LV. Co-regulation of genes encoding excitatory and inhibitory processes in cardiac tissue by the acetylation modifiers may help explain their predominantly net-neutral effects on cardiac electrophysiology. ATP1A1, encoding for the Na/K pump, represented an outlier—with orthogonal regulation by the acetylation modifiers to most of the ICs. The HDAC/SIRT/HAT effects were mediated by strong (+) TF regulators of ICs, e.g., MEF2A and TBX5, in both sexes. Furthermore, for male hearts, PLS models revealed a stronger (+/-) mediatory role on ICs for NKX25 and TGF1B/KLF4, respectively, while RUNX1 exhibited larger (-) TF effects on ICs in females. Male-trained PLS models of HDAC/SIRT/HAT effects on ICs underestimated the effects on some ICs in females. Insights from the GTEx dataset about the co-expression and transcriptional co-regulation of acetylation-modifying enzymes, transcription factors and key cardiac ion channels in a sex-specific manner can help inform safer drug design

Analysis for co-occurring sequence features identifies link between common synonymous variant and an early-terminated NPC1 isoform

Direct assessment of allelic phase for DNA and RNA features of diploid genomes has been challenging for Sanger sequencing, due to its allele-conflating base-calling signal. Massively parallel sequencing technologies are based on the generation of a continuous copy of a single strand sequence segments, thus preserving the allelic relation between the features of the original molecules. We have performed a transcriptome-wide search for co-occurrence of variant nucleotides and exon-intron boundaries positioned within the length of a single sequencing read. Analysis of 75 human transcriptomes from retinal pigment epithelia (RPE), glioblastoma, low-grade brain tumor, breast cancer and colon cancer, have identified an association between the synonymous variant rs1140458 and an early-terminated NPC1 isoform lacking exons 19–25. Higher proportion of molecules bearing the variant nucleotide (versus the reference) incorporates the intron (P \u3c0.0001), which turns the last codon of exon 18 into a stop codon. The significance is highest in RPE cells (P = 3.88 × 10−12). NPC1 protein is involved in the control of the cholesterol trafficking. NPC1 mutations lead, in an autosomal recessive manner, to the neurological disorder Niemann-Pick syndrome type C (NP-C), and, ablation of NPC1 causes age-progressive retinal degeneration in mice and drosophila. The vast majority of the NP-C causative variants consist of missense/nonsense substitutions, small indels, and, intronic splice variants. Rs1140458 is a common exonic synonymous substitution that has never been linked to alternative splicing or pathogenicity. Our analysis suggests that rs1140458 may affect the levels of the functional NPC1 protein, and to contribute to some of the cholesterol-implicated cellular phenotype

Exosomes containing HIV protein Nef reorganize lipid rafts potentiating inflammatory response in bystander cells.

HIV infection has a profound effect on "bystander" cells causing metabolic co-morbidities. This may be mediated by exosomes secreted by HIV-infected cells and containing viral factors. Here we show that exosomes containing HIV-1 protein Nef (exNef) are rapidly taken up by macrophages releasing Nef into the cell interior. This caused down-regulation of ABCA1, reduction of cholesterol efflux and sharp elevation of the abundance of lipid rafts through reduced activation of small GTPase Cdc42 and decreased actin polymerization. Changes in rafts led to re-localization of TLR4 and TREM-1 to rafts, phosphorylation of ERK1/2, activation of NLRP3 inflammasome, and increased secretion of pro-inflammatory cytokines. The effects of exNef on lipid rafts and on inflammation were reversed by overexpression of a constitutively active mutant of Cdc42. Similar effects were observed in macrophages treated with exosomes produced by HIV-infected cells or isolated from plasma of HIV-infected subjects, but not with exosomes from cells and subjects infected with ΔNef-HIV or uninfected subjects. Mice injected with exNef exhibited monocytosis, reduced ABCA1 in macrophages, increased raft abundance in monocytes and augmented inflammation. Thus, Nef-containing exosomes potentiated pro-inflammatory response by inducing changes in cholesterol metabolism and reorganizing lipid rafts. These mechanisms may contribute to HIV-associated metabolic co-morbidities

Systematic pan-cancer analysis of somatic allele frequency

Imbalanced expression of somatic alleles in cancer can suggest functional and selective features, and can therefore indicate possible driving potential of the underlying genetic variants. To explore the correlation between allele frequency of somatic variants and total gene expression of their harboring gene, we used the unique data set of matched tumor and normal RNA and DNA sequencing data of 5523 distinct single nucleotide variants in 381 individuals across 10 cancer types obtained from The Cancer Genome Atlas (TCGA). We analyzed the allele frequency in the context of the variant and gene functional features and linked it with changes in the total gene expression. We documented higher allele frequency of somatic variants in cancer-implicated genes (Cancer Gene Census, CGC). Furthermore, somatic alleles bearing premature terminating variants (PTVs), when positioned in CGC genes, appeared to be less frequently degraded via nonsense-mediated mRNA decay, indicating possible favoring of truncated proteins by the tumor transcriptome. Among the genes with multiple PTVs with high allele frequency, ARID1, TP53 and NSD1 were known key cancer genes. All together, our analyses suggest that high allele frequency of tumor somatic variants can indicate driving functionality and can serve to identify potential cancer-implicated genes

RNA sequencing of cancer reveals novel splicing alterations

Breast cancer transcriptome acquires a myriad of regulation changes, and splicing is critical for the cell to “tailor-make” specific functional transcripts. We systematically revealed splicing signatures of the three most common types of breast tumors using RNA sequencing: TNBC, non-TNBC and HER2-positive breast cancer. We discovered subtype specific differentially spliced genes and splice isoforms not previously recognized in human transcriptome. Further, we showed that exon skip and intron retention are predominant splice events in breast cancer. In addition, we found that differential expression of primary transcripts and promoter switching are significantly deregulated in breast cancer compared to normal breast. We validated the presence of novel hybrid isoforms of critical molecules like CDK4, LARP1, ADD3, and PHLPP2. Our study provides the first comprehensive portrait of transcriptional and splicing signatures specific to breast cancer sub-types, as well as previously unknown transcripts that prompt the need for complete annotation of tissue and disease specific transcriptome

Overexpressed Somatic Alleles are Enriched in Functional Elements in Breast Cancer.

Asymmetric allele content in the transcriptome can be indicative of functional and selective features of the underlying genetic variants. Yet, imbalanced alleles, especially from diploid genome regions, are poorly explored in cancer. Here we systematically quantify and integrate the variant allele fraction from corresponding RNA and DNA sequence data from patients with breast cancer acquired through The Cancer Genome Atlas (TCGA). We test for correlation between allele prevalence and functionality in known cancer-implicated genes from the Cancer Gene Census (CGC). We document significant allele-preferential expression of functional variants in CGC genes and across the entire dataset. Notably, we find frequent allele-specific overexpression of variants in tumor-suppressor genes. We also report a list of over-expressed variants from non-CGC genes. Overall, our analysis presents an integrated set of features of somatic allele expression and points to the vast information content of the asymmetric alleles in the cancer transcriptome

RNA2DNAlign: nucleotide resolution allele asymmetries through quantitative assessment of RNA and DNA paired sequencing data.

We introduce RNA2DNAlign, a computational framework for quantitative assessment of allele counts across paired RNA and DNA sequencing datasets. RNA2DNAlign is based on quantitation of the relative abundance of variant and reference read counts, followed by binomial tests for genotype and allelic status at SNV positions between compatible sequences. RNA2DNAlign detects positions with differential allele distribution, suggesting asymmetries due to regulatory/structural events. Based on the type of asymmetry, RNA2DNAlign outlines positions likely to be implicated in RNA editing, allele-specific expression or loss, somatic mutagenesis or loss-of-heterozygosity (the first three also in a tumor-specific setting). We applied RNA2DNAlign on 360 matching normal and tumor exomes and transcriptomes from 90 breast cancer patients from TCGA. Under high-confidence settings, RNA2DNAlign identified 2038 distinct SNV sites associated with one of the aforementioned asymetries, the majority of which have not been linked to functionality before. The performance assessment shows very high specificity and sensitivity, due to the corroboration of signals across multiple matching datasets. RNA2DNAlign is freely available from http://github.com/HorvathLab/NGS as a self-contained binary package for 64-bit Linux systems

Novel insights into breast cancer genetic variance through RNA sequencing

Using RNA sequencing of triple-negative breast cancer (TNBC), non-TBNC and HER2-positive breast cancer sub-types, here we report novel expressed variants, allelic prevalence and abundance, and coexpression with other variation, and splicing signatures. To reveal the most prevalent variant alleles, we overlaid our findings with cancer- and population-based datasets and validated a subset of novel variants of cancer-related genes: ESRP2, GBP1, TPP1, MAD2L1BP, GLUD2 and SLC30A8. As a proof-of-principle, we demonstrated that a rare substitution in the splicing coordinator ESRP2(R353Q) impairs its ability to bind to its substrate FGFR2 pre-mRNA. In addition, we describe novel SNPs and INDELs in cancer relevant genes with no prior reported association of point mutations with cancer, such as MTAP and MAGED1. For the first time, this study illustrates the power of RNA-sequencing in revealing the variation landscape of breast transcriptome and exemplifies analytical strategies to search regulatory interactions among cancer relevant molecules

Alternative Splicing Promotes Tumour Aggressiveness and Drug Resistance in African American Prostate Cancer.

linical challenges exist in reducing prostate cancer (PCa) disparities. The RNA splicing landscape of PCa across racial populations has not been fully explored as a potential molecular mechanism contributing to race-related tumour aggressiveness. Here, we identify novel genome-wide, race-specific RNA splicing events as critical drivers of PCa aggressiveness and therapeutic resistance in African American (AA) men. AA-enriched splice variants of PIK3CD, FGFR3, TSC2 and RASGRP2 contribute to greater oncogenic potential compared with corresponding European American (EA)-expressing variants. Ectopic overexpression of the newly cloned AA-enriched variant, PIK3CD-S, in EA PCa cell lines enhances AKT/mTOR signalling and increases proliferative and invasive capacity in vitro and confers resistance to selective PI3Kδ inhibitor, CAL-101 (idelalisib), in mouse xenograft models. High PIK3CD-S expression in PCa specimens associates with poor survival. These results highlight the potential of RNA splice variants to serve as novel biomarkers and molecular targets for developmental therapeutics in aggressive PCa