Detection of redundant fusion transcripts as biomarkers or disease-specific therapeutic targets in breast cancer

Genome-Wide Analysis of Rare Haplotypes Associated with Breast Cancer Risk

Numerous common genetic variants have been linked to breast cancer risk, but they only partially explain the total breast cancer heritability. Inference from Nordic population-based twin data indicates rare high-risk loci as the chief determinant of breast cancer risk. Here, we use haplotypes, rather than single variants, to identify rare high-risk loci for breast cancer. With computationally phased genotypes from 181,034 white British women in the UK Biobank, a genome-wide haplotype-breast cancer association analysis was conducted using sliding windows of 5 to 500 consecutive array-genotyped variants. In the discovery stage, haplotype-breast cancer associations were evaluated retrospectively in the prestudy-enrollment data including 5,487 breast cancer cases. Breast cancer hazard ratios (HR) for additive haplotypic effects were estimated using Cox regression. The replication analysis included a prospective cohort of women free of breast cancer at enrollment, of whom 3,524 later developed breast cancer. This two-stage analysis detected 13 rare loci (frequency <1%), each associated with an appreciable breast cancer-risk increase (discovery: HRs = 2.84-6.10, P < 5 × 10-8; replication: HRs = 2.08-5.61, P < 0.01). In contrast, the variants that formed these rare haplotypes individually exhibited much smaller effects. Functional annotation revealed extensive cis-regulatory DNA elements in breast cancer-related cells underlying the replicated rare haplotypes. Using phased, imputed genotypes from 30,064 cases and 25,282 controls in the DRIVE OncoArray case-control study, 6 of the 13 rare-loci associations were found generalizable (odds ratio estimates: 1.48-7.67, P < 0.05). This study demonstrates the complementary advantage of utilizing rare haplotypes to capture novel risk loci and suggests the potential for the discovery of more genetic elements contributing to cancer heritability as large data sets of germline whole-genome sequencing become available.

SIGNIFICANCE

A genome-wide two-stage haplotype analysis identifies rare haplotypes associated with breast cancer risk and suggests that the rare risk haplotypes represent long-range interactions with regulatory consequences influencing cancer risk.

Characterization of circular RNAs with advanced sequencing technologies in human complex diseases

Circular RNAs (circRNAs) are one category of non-coding RNAs that do not possess 5' caps and 3' free ends. Instead, they are derived in closed circle forms from pre-mRNAs by a non-canonical splicing mechanism named "back-splicing." CircRNAs were discovered four decades ago, initially called "scrambled exons." Compared to linear RNAs, the expression levels of circRNAs are considerably lower, and it is challenging to identify circRNAs specifically. Thus, the biological relevance of circRNAs has been underappreciated until the advancement of next generation sequencing (NGS) technology. The biological insights of circRNAs, such as their tissue-specific expression patterns, biogenesis factors, and functional effects in complex diseases, namely human cancers, have been extensively explored in the last decade. With the invention of the third generation sequencing (TGS) with longer sequencing reads and newly designed strategies to characterize full-length circRNAs, the panorama of circRNAs in human complex diseases could be further unveiled. In this review, we first introduce the history of circular RNA detection. Next, we describe widely adopted NGS-based methods and the recently established TGS-based approaches capable of characterizing circRNAs in full-length. We then summarize data resources and representative circRNA functional studies related to human complex diseases. In the last section, we reviewed computational tools and discuss the potential advantages of utilizing advanced sequencing approaches to a functional interpretation of full-length circRNAs in complex diseases. This article is categorized under: RNA Evolution and Genomics > Computational Analyses of RNA RNA in Disease and Development > RNA in Disease.

A Novel Bead-Capture Nanopore Sequencing Method for Large Structural Rearrangement Detection in Cancer

Rapid, cost-effective genomic stratification of structural rearrangements in cancer is often of vital importance when determining treatment; however, existing diagnostic cytogenetic and molecular testing fails to deliver the required speed when deployed at scale. Next-generation sequencing-based methods are widely used, but these can lack sensitivity and require batching of samples to be cost-effective, with long turnaround times. Here we present a novel method for rearrangement detection from genomic DNA based on third-generation long-read sequencing that overcomes these time and cost issues. The utility of this approach for the genomic stratification of patients with acute myeloid leukemia is shown based on detection of four of the most prevalent structural rearrangements. The method not only determines the precise genomic breakpoint for each expected rearrangement but also discovers and validates novel translocations in one-third of the tested samples, 80% of which involve known oncogenes. This method may prove to be a powerful tool for the diagnosis, genomic stratification, and characterization of cancers.

Novel read-through fusion transcript Bcl2l2-Pabpn1 in glioblastoma cells

Read‐through fusion transcripts have recently been identified as chimeric RNAs and have since been linked to tumour growth in some cases. Many fusion genes generated by chromosomal rearrangements have been described in glioblastoma. However, read‐through fusion transcripts between neighbouring genes in glioblastoma remain unexplored. We performed paired‐end RNA‐seq of rat C6 glioma cells and normal cells and discovered a read‐through fusion transcript Bcl2l2‐Pabpn1 in which exon 3 of Bcl‐2‐like protein 2 (Bcl2l2) fused to exon 2 of Polyadenylate‐binding protein 1 (Pabpn1). This fusion transcript was found in both human glioblastoma and normal cells. Unlike other fusions reported in glioblastoma, Bcl2l2‐Pabpn1 appeared to result from RNA processing rather than genomic rearrangement. Bcl2l2‐Pabpn1 fusion transcript encoded a fusion protein with BH4, BCL and RRM domains. Functionally, Bcl2l2‐Pabpn1 knockdown by targeting its fusion junction decreased its expression, and suppressed cell proliferation, migration and invasion in vitro. Mechanistically, Bcl2l2‐Pabpn1 blocked Bax activity and activated PI3K/AKT pathway to promote glioblastoma progression. Together, our work characterized a glioblastoma‐associated Bcl2l2‐Pabpn1 fusion transcript shared by humans and rats.

Resolving missing protein problems using functional class scoring

Despite technological advances in proteomics, incomplete coverage and inconsistency issues persist, resulting in “data holes”. These data holes cause the missing protein problem (MPP), where relevant proteins are persistently unobserved, or sporadically observed across samples, hindering biomarker discovery and proper functional characterization. Network-based approaches can provide powerful solutions for resolving these issues. Functional Class Scoring (FCS) is one such method that uses protein complex information to recover missing proteins with weak support. However, FCS has not been evaluated on more recent proteomic technologies with higher coverage, and there is no clear way to evaluate its performance. To address these issues, we devised a more rigorous evaluation schema based on cross-verification between technical replicates and evaluated its performance on data acquired under recent Data-Independent Acquisition (DIA) technologies (viz. SWATH). Although cross-replicate examination reveals some inconsistencies amongst same-class samples, tissue-differentiating signal is nonetheless strongly conserved, confirming that FCS selects for biologically meaningful networks. We also report that predicted missing proteins are statistically significant based on FCS p values. Despite limited cross-replicate verification rates, the predicted missing proteins as a whole have higher peptide support than non-predicted proteins. FCS also predicts missing proteins that are often lost due to weak specific peptide support.

The Landscape of Novel Expressed Chimeric RNAs in Rheumatoid Arthritis

In cancers and other complex diseases, the fusion of two genes can lead to the production of chimeric RNAs, which are associated with disease development. Several recurrent chimeric RNAs are expressed in different cancers and are thus used for clinical cancer diagnosis. Rheumatoid arthritis (RA) is an immune-mediated joint disorder resulting in synovial inflammation and joint destruction. Despite advances in therapy, many patients do not respond to treatment and present persistent inflammation. Understanding the landscape of chimeric RNA expression in RA patients could provide a better insight into RA pathogenesis, which might provide better treatment strategies and tailored therapies. Accordingly, we analyzed the publicly available RNA-seq data of synovium tissue from 151 RA patients and 28 healthy controls and were able to identify 37 recurrent chimeric RNAs found to be expressed in at least 3 RA samples. Furthermore, the parental genes of these 37 recurrent chimeric RNAs were found to be differentially expressed and enriched in immune-related processes, such as adaptive immune response and the positive regulation of B-cell activation. Interestingly, the appearance of 5 coding and 23 non-coding chimeric RNAs might be associated with regulating their parental gene expression, leading to the generation of dysfunctional immune responses, such as inflammation and bone destruction. Therefore, in this paper, we present the first study to demonstrate the novel chimeric RNAs that are highly expressed and functional in RA.

Identifying the oncogenic potential of gene fusions exploiting miRNAs

It is estimated that oncogenic gene fusions cause about 20% of human cancer morbidity. Identifying potentially oncogenic gene fusions may improve affected patients' diagnosis and treatment. Previous approaches to this issue included exploiting specific gene-related information, such as gene function and regulation. Here we propose a model that profits from the previous findings and includes the microRNAs in the oncogenic assessment. We present ChimerDriver, a tool to classify gene fusions as oncogenic or not oncogenic. ChimerDriver is based on a specifically designed neural network and trained on genetic and post-transcriptional information to obtain a reliable classification. The designed neural network integrates information related to transcription factors, gene ontologies, microRNAs and other detailed information related to the functions of the genes involved in the fusion and the gene fusion structure. As a result, the performances on the test set reached 0.83 f1-score and 96% recall. The comparison with state-of-the-art tools returned comparable or higher results. Moreover, ChimerDriver performed well in a real-world case where 21 out of 24 validated gene fusion samples were detected by the gene fusion detection tool Starfusion. ChimerDriver integrates transcriptional and post-transcriptional information in an ad-hoc designed neural network to effectively discriminate oncogenic gene fusions from passenger ones. ChimerDriver source code is freely available at https://github.com/martalovino/ChimerDriver.

Computational analysis of sense-antisense chimeric transcripts reveals their potential regulatory features and the landscape of expression in human cells

Many human genes are transcribed from both strands and produce sense-antisense gene pairs. Sense-antisense (SAS) chimeric transcripts are produced upon the coalescing of exons/introns from both sense and antisense transcripts of the same gene. SAS chimera was first reported in prostate cancer cells. Subsequently, numerous SAS chimeras have been reported in the ChiTaRS-2.1 database. However, the landscape of their expression in human cells and functional aspects are still unknown. We found that longer palindromic sequences are a unique feature of SAS chimeras. Structural analysis indicates that a long hairpin-like structure formed by many consecutive Watson-Crick base pairs appears because of these long palindromic sequences, which possibly play a similar role as double-stranded RNA (dsRNA), interfering with gene expression. RNA-RNA interaction analysis suggested that SAS chimeras could significantly interact with their parental mRNAs, indicating their potential regulatory features. Here, 267 SAS chimeras were mapped in RNA-seq data from 16 healthy human tissues, revealing their expression in normal cells. Evolutionary analysis suggested the positive selection favoring sense-antisense fusions that significantly impacted the evolution of their function and structure. Overall, our study provides detailed insight into the expression landscape of SAS chimeras in human cells and identifies potential regulatory features.

Comparative study of bioinformatic tools for the identification of chimeric RNAs from RNA Sequencing

Chimeric RNAs are gaining more and more attention as they have broad implications in both cancer and normal physiology. To date, over 40 chimeric RNA prediction methods have been developed to facilitate their identification from RNA sequencing data. However, a limited number of studies have been conducted to compare the performance of these tools; additionally, previous studies have become outdated as more software tools have been developed within the last three years. In this study, we benchmarked 16 chimeric RNA prediction software, including seven top performers in previous benchmarking studies, and nine that were recently developed. We used two simulated and two real RNA-Seq datasets, compared the 16 tools for their sensitivity, positive prediction value (PPV), F-measure, and also documented the computational requirements (time and memory). We noticed that none of the tools are inclusive, and their performance varies depending on the dataset and objects. To increase the detection of true positive events, we also evaluated the pair-wise combination of these methods to suggest the best combination for sensitivity and F-measure. In addition, we compared the performance of the tools for the identification of three classes (read-through, inter-chromosomal and intra-others) of chimeric RNAs. Finally, we performed TOPSIS analyses and ranked the weighted performance of the 16 tools.

Transcriptome Analysis Reveals MFGE8-HAPLN3 Fusion as a Novel Biomarker in Triple-Negative Breast Cancer

Background

Triple-negative breast cancer (TNBC) is a highly aggressive cancer with poor prognosis. The lack of effective targeted therapies for TNBC remains a profound clinical challenge. Fusion transcripts play critical roles in carcinogenesis and serve as valuable diagnostic and therapeutic targets in cancer. The present study aimed to identify novel fusion transcripts in TNBC.

Methods

We analyzed the RNA sequencing data of 360 TNBC samples to identify and filter fusion candidates through SOAPfuse and ChimeraScan analysis. The characteristics, including recurrence, fusion type, chromosomal localization, TNBC subgroup distribution, and clinicopathological correlations, were analyzed in all candidates. Furthermore, we selected the promising fusion transcript and predicted its fusion type and protein coding capacity.

Results

Using the RNA sequencing data, we identified 189 fusion transcripts in TNBC, among which 22 were recurrent fusions. Compared to para-tumor tissues, TNBC tumor tissues accumulated more fusion events, especially in high-grade tumors. Interestingly, these events were enriched at specific chromosomal loci, and the distribution pattern varied in different TNBC subtypes. The vast majority of fusion partners were discovered on chromosomes 1p, 11q, 19p, and 19q. Besides, fusion events mainly clustered on chromosome 11 in the immunomodulatory subtype and chromosome 19 in the luminal androgen receptor subtype of TNBC. Considering the tumor specificity and frameshift mutation, we selected MFGE8-HAPLN3 as a novel biomarker and further validated it in TNBC samples using PCR and Sanger sequencing. Further, we successfully identified three types of MFGE8-HAPLN3 (E6-E2, E5-E3, and E6-E3) and predicted the ORF of E6-E2, which could encode a protein of 712 amino acids, suggesting its critical role in TNBC.

Conclusions

Improved bioinformatic stratification and comprehensive analysis identified the fusion transcript MFGE8-HAPLN3 as a novel biomarker with promising clinical application in the future.

Human testis-expressed (TEX) genes: a review focused on spermatogenesis and male fertility

Spermatogenesis is a complex process regulated by a multitude of genes. The identification and characterization of male-germ-cell-specific genes is crucial to understanding the mechanisms through which the cells develop. The term “TEX gene” was coined by Wang et al. (Nat Genet. 2001; 27: 422–6) after they used cDNA suppression subtractive hybridization (SSH) to identify new transcripts that were present only in purified mouse spermatogonia. TEX (Testis expressed) orthologues have been found in other vertebrates (mammals, birds, and reptiles), invertebrates, and yeasts. To date, 69 TEX genes have been described in different species and different tissues. To evaluate the expression of each TEX/tex gene, we compiled data from 7 different RNA-Seq mRNA databases in humans, and 4 in the mouse according to the expression atlas database.

Various studies have highlighted a role for many of these genes in spermatogenesis. Here, we review current knowledge on the TEX genes and their roles in spermatogenesis and fertilization in humans and, comparatively, in other species (notably the mouse). As expected, TEX genes appear to have a major role in reproduction in general and in spermatogenesis in humans but also in all mammals such as the mouse. Most of them are expressed specifically or predominantly in the testis. As most of the TEX genes are highly conserved in mammals, defects in the male (gene mutations in humans and gene-null mice) lead to infertility. In the future, cumulative data on the human TEX genes’ physiological functions and pathophysiological dysfunctions should become available and is likely to confirm the essential role of this family in the reproductive process. Thirteen TEX genes are now referenced in the OMIM database, and 3 have been linked to a specific phenotype. TEX11 (on Xq13.1) is currently the gene most frequently reported as being associated with azoospermia.

Landscape of Chimeric RNAs in Non-Cancerous Cells

Gene fusions and their products (RNA and protein) have been traditionally recognized as unique features of cancer cells and are used as ideal biomarkers and drug targets for multiple cancer types. However, recent studies have demonstrated that chimeric RNAs generated by intergenic alternative splicing can also be found in normal cells and tissues. In this study, we aim to identify chimeric RNAs in different non-neoplastic cell lines and investigate the landscape and expression of these novel candidate chimeric RNAs. To do so, we used HEK-293T, HUVEC, and LO2 cell lines as models, performed paired-end RNA sequencing, and conducted analyses for chimeric RNA profiles. Several filtering criteria were applied, and the landscape of chimeric RNAs was characterized at multiple levels and from various angles. Further, we experimentally validated 17 chimeric RNAs from different classifications. Finally, we examined a number of validated chimeric RNAs in different cancer and non-cancer cells, including blood from healthy donors, and demonstrated their ubiquitous expression pattern.

Identification of m6A-Associated RNA Binding Proteins Using an Integrative Computational Framework

N6-methyladenosine (m⁶A) is an abundant modification on mRNA that plays an important role in regulating essential RNA activities. Several wet lab studies have identified some RNA binding proteins (RBPs) that are related to m⁶A's regulation. The objective of this study was to identify potential m⁶A-associated RBPs using an integrative computational framework. The framework was composed of an enrichment analysis and a classification model. Utilizing RBPs' binding data, we analyzed reproducible m⁶A regions from independent studies using this framework. The enrichment analysis identified known m⁶A-associated RBPs including YTH domain-containing proteins; it also identified RBM3 as a potential m⁶A-associated RBP for mouse. Furthermore, a significant correlation for the identified m⁶A-associated RBPs is observed at the protein expression level rather than the gene expression level. On the other hand, a Random Forest classification model was built for the reproducible m⁶A regions using RBPs' binding data. The RBP-based predictor demonstrated not only competitive performance when compared with sequence-based predictions but also reflected m⁶A's action of repelling against RBPs, which suggested that our framework can infer interaction between m⁶A and m⁶A-associated RBPs beyond sequence level when utilizing RBPs' binding data. In conclusion, we designed an integrative computational framework for the identification of known and potential m⁶A-associated RBPs. We hope the analysis will provide more insights on the studies of m⁶A and RNA modifications.

ChiRA: an integrated framework for chimeric read analysis from RNA-RNA interactome and RNA structurome data

BACKGROUND

With the advances in next-generation sequencing technologies, it is possible to determine RNA-RNA interaction and RNA structure predictions on a genome-wide level. The reads from these experiments usually are chimeric, with each arm generated from one of the interaction partners. Owing to short read lengths, often these sequenced arms ambiguously map to multiple locations. Thus, inferring the origin of these can be quite complicated. Here we present ChiRA, a generic framework for sensitive annotation of these chimeric reads, which in turn can be used to predict the sequenced hybrids.

RESULTS

Grouping reference loci on the basis of aligned common reads and quantification improved the handling of the multi-mapped reads in contrast to common strategies such as the selection of the longest hit or a random choice among all hits. On benchmark data ChiRA improved the number of correct alignments to the reference up to 3-fold. It is shown that the genes that belong to the common read loci share the same protein families or similar pathways. In published data, ChiRA could detect 3 times more new interactions compared to existing approaches. In addition, ChiRAViz can be used to visualize and filter large chimeric datasets intuitively.

CONCLUSION

ChiRA tool suite provides a complete analysis and visualization framework along with ready-to-use Galaxy workflows and tutorials for RNA-RNA interactome and structurome datasets. Common read loci built by ChiRA can rescue multi-mapped reads on paralogous genes without requiring any information on gene relations. We showed that ChiRA is sensitive in detecting new RNA-RNA interactions from published RNA-RNA interactome datasets.

Five multicopy gene family genes expressed during the maternal-to-zygotic transition are not essential for mouse development

Upon fertilization, oocytes transform into totipotent and pluripotent cleavage stage cells through the maternal-to-zygotic transition (MZT), which is regulated by maternal factors and zygotic genome activation (ZGA). Here, we investigated the in vivo function of 16 genes expressed with strong biases in oocytes and cleavage stage embryos by generating knockout (KO) mice. These MZT-associated genes are conserved across many mammalian species and include five multicopy gene family genes: the Nlrp9, Khdc1, Rfpl4, Trim43, and Zscan5 genes. Intercrosses between female KO and male KO mice, including Nlrp9a/b/c triple KO (TKO), Khdc1a/b/c TKO, Rfpl4a/b double KO (DKO), Trim43a/b/c TKO, and Zscan5b KO mice led to the birth to healthy offspring that in turn produced healthy offspring. Our study not only demonstrated that these MZT-associated genes are not essential for mouse development, but also provides valuable resources for analyzing the functions of these genes in other genetic backgrounds, in the presence of stressors, and under pathogenic conditions.

Update on SLC6A14 in lung and gastrointestinal physiology and physiopathology: focus on cystic fibrosis

The solute carrier family 6 member 14 (SLC6A14) protein imports and concentrates all neutral amino acids as well as the two cationic acids lysine and arginine into the cytoplasm of different cell types. Primarily described as involved in several cancer and colonic diseases physiopathological mechanisms, the SLC6A14 gene has been more recently identified as a genetic modifier of cystic fibrosis (CF) disease severity. It was indeed shown to have a pleiotropic effect, modulating meconium ileus occurrence, lung disease severity, and precocity of P. aeruginosa airway infection. The biological mechanisms explaining the impact of SLC6A14 on intestinal and lung phenotypes of CF patients are starting to be elucidated. This review focuses on SLC6A14 in lung and gastrointestinal physiology and physiopathology, especially its involvement in the pathophysiology of CF disease.

The landscape of chimeric RNAs in non-diseased tissues and cells

Chimeric RNAs and their encoded proteins have been traditionally viewed as unique features of neoplasia, and have been used as biomarkers and therapeutic targets for multiple cancers. Recent studies have demonstrated that chimeric RNAs also exist in non-cancerous cells and tissues, although large-scale, genome-wide studies of chimeric RNAs in non-diseased tissues have been scarce. Here, we explored the landscape of chimeric RNAs in 9495 non-diseased human tissue samples of 53 different tissues from the GTEx project. Further, we established means for classifying chimeric RNAs, and observed enrichment for particular classifications as more stringent filters are applied. We experimentally validated a subset of chimeric RNAs from each classification and demonstrated functional relevance of two chimeric RNAs in non-cancerous cells. Importantly, our list of chimeric RNAs in non-diseased tissues overlaps with some entries in several cancer fusion databases, raising concerns for some annotations. The data from this study provides a large repository of chimeric RNAs present in non-diseased tissues, which can be used as a control dataset to facilitate the identification of true cancer-specific chimeras.

Peptidyl Arginine Deiminase 2 (PADI2)-Mediated Arginine Citrullination Modulates Transcription in Cancer

Protein arginine deimination leading to the non-coded amino acid citrulline remains a key question in the field of post-translational modifications ever since its discovery by Rogers and Simmonds in 1958. Citrullination is catalyzed by a family of enzymes called peptidyl arginine deiminases (PADIs). Initially, increased citrullination was associated with autoimmune diseases, including rheumatoid arthritis and multiple sclerosis, as well as other neurological disorders and multiple types of cancer. During the last decade, research efforts have focused on how citrullination contributes to disease pathogenesis by modulating epigenetic events, pluripotency, immunity and transcriptional regulation. However, our knowledge regarding the functional implications of citrullination remains quite limited, so we still do not completely understand its role in physiological and pathological conditions. Here, we review the recently discovered functions of PADI2-mediated citrullination of the C-terminal domain of RNA polymerase II in transcriptional regulation in breast cancer cells and the proposed mechanisms to reshape the transcription regulatory network that promotes cancer progression.

ChiTaRS 5.0: the comprehensive database of chimeric transcripts matched with druggable fusions and 3D chromatin maps

Chimeric RNA transcripts are formed when exons from two genes fuse together, often due to chromosomal translocations, transcriptional errors or trans-splicing effect. While these chimeric RNAs produce functional proteins only in certain cases, they play a significant role in disease phenotyping and progression. ChiTaRS 5.0 (http://chitars.md.biu.ac.il/) is the latest and most comprehensive chimeric transcript repository, with 111 582 annotated entries from eight species, including 23 167 known human cancer breakpoints. The database includes unique information correlating chimeric breakpoints with 3D chromatin contact maps, generated from public datasets of chromosome conformation capture techniques (Hi-C). In this update, we have added curated information on druggable fusion targets matched with chimeric breakpoints, which are applicable to precision medicine in cancers. The introduction of a new section that lists chimeric RNAs in various cell-lines is another salient feature. Finally, using text-mining techniques, novel chimeras in Alzheimer's disease, schizophrenia, dyslexia and other diseases were collected in ChiTaRS. Thus, this improved version is an extensive catalogue of chimeras from multiple species. It extends our understanding of the evolution of chimeric transcripts in eukaryotes and contributes to the analysis of 3D genome conformational changes and the functional role of chimeras in the etiopathogenesis of cancers and other complex diseases.

Multivariable regulation of gene expression plasticity in metazoans

An important capacity of genes is the rapid change of expression levels to cope with the environment, known as expression responsiveness or plasticity. Elucidating the genomic mechanisms determining expression plasticity is critical for understanding the molecular basis of phenotypic plasticity, fitness and adaptation. In this study, we systematically quantified gene expression plasticity in four metazoan species by integrating changes of expression levels under a large number of genetic and environmental conditions. From this, we demonstrated that expression plasticity measures a distinct feature of gene expression that is orthogonal to other well-studied features, including gene expression level and tissue specificity/broadness. Expression plasticity is conserved across species with important physiological implications. The magnitude of expression plasticity is highly correlated with gene function and genes with high plasticity are implicated in disease susceptibility. Genome-wide analysis identified many conserved promoter cis-elements, trans-acting factors (such as CTCF), and gene body histone modifications (H3K36me3, H3K79me2 and H4K20me1) that are significantly associated with expression plasticity. Analysis of expression changes in perturbation experiments further validated a causal role of specific transcription factors and histone modifications. Collectively, this work reveals the general properties, physiological implications and multivariable regulation of gene expression plasticity in metazoans, extending the mechanistic understanding of gene regulation.

Toward a holistic view of multiscale breast cancer molecular biomarkers

Powered by rapid technology developments, biomarkers become increasingly diverse, including those detected at genomic, transcriptomic, proteomic, metabolomic and cellular levels. While diverse sets of biomarkers have been utilized in breast cancer predisposition, diagnosis, prognosis, treatment and management, recent additions derived from lincRNA, circular RNA, circulating DNA together with its methylated and hydroxymethylated forms and immune signatures are likely to further transform clinical practice. Here, we take breast cancer as an example of heterogeneous diseases that require many informed decisions from treatment to care to review the huge variety of biomarkers. By assessing the advantages and limitations of modern biomarkers in diverse use scenarios, this article outlines the prospects and challenges of releasing complimentary advantages by augmentation of multiscale molecular biomarkers.

High expression of the vacuole membrane protein 1 (VMP1) is a potential marker of poor prognosis in HER2 positive breast cancer

Background

Fusion genes result from genomic structural changes, which can lead to alterations in gene expression that supports tumor development. The aim of the study was to use fusion genes as a tool to identify new breast cancer (BC) genes with a role in BC progression.

Methods

Fusion genes from breast tumors and BC cell lines were collected from publications. RNA-Seq data from tumors and cell lines were retrieved from databanks and analyzed for fusions with SOAPfuse or the analysis was purchased. Fusion genes identified in both tumors (n = 1724) and cell lines (n = 45) were confirmed by qRT-PCR and sequencing. Their individual genes were ranked by selection criteria that included correlation of their mRNA level with copy number. The expression of the top ranked gene was measured by qRT-PCR in normal tissue and in breast tumors from an exploratory cohort (n = 141) and a validation cohort (n = 277). Expression levels were correlated with clinical and pathological factors as well as the patients’ survival. The results were followed up in BC cohorts from TCGA (n = 818) and METABRIC (n = 2509).

Results

Vacuole membrane protein 1 (VMP1) was the most promising candidate based on specific selection criteria. Its expression was higher in breast tumor tissue than normal tissue (p = 1x10^-4), and its expression was significantly higher in HER2 positive than HER2 negative breast tumors in all four cohorts analyzed. High expression of VMP1 associated with breast cancer specific survival (BCSS) in cohort 1 (hazard ratio (HR) = 2.31, CI 1.27–4.18) and METABRIC (HR = 1.26, CI 1.02–1.57), and also after adjusting for HER2 expression in cohort 1 (HR = 2.03, CI 1.10–3.72). BCSS was not significant in cohort 2 or TCGA cohort, which may be due to differences in treatment regimens.

Conclusions

The results suggest that high VMP1 expression is a potential marker of poor prognosis in HER2 positive BC. Further studies are needed to elucidate how VMP1 could affect pathways supportive of tumorigenesis.

The pivotal role of sampling recurrent tumors in the precision care of patients with tumors of the central nervous system

Effective management of brain and spine tumors relies on a multidisciplinary approach encompassing surgery, radiation, and systemic therapy. In the era of personalized oncology, the latter is complemented by various molecularly targeting agents. Precise identification of cellular targets for these drugs requires comprehensive profiling of the cancer genome coupled with an efficient analytic pipeline, leading to an informed decision on drug selection, prognosis, and confirmation of the original pathological diagnosis. Acquisition of optimal tumor tissue for such analysis is paramount and often presents logistical challenges in neurosurgery. Here, we describe the experience and results of the Personalized OncoGenomics (POG) program with a focus on tumors of the central nervous system (CNS). Patients with recurrent CNS tumors were consented and enrolled into the POG program prior to accrual of tumor and matched blood followed by whole-genome and transcriptome sequencing and processing through the POG bioinformatic pipeline. Sixteen patients were enrolled into POG. In each case, POG analyses identified genomic drivers including novel oncogenic fusions, aberrant pathways, and putative therapeutic targets. POG has highlighted that personalized oncology is truly a multidisciplinary field, one in which neurosurgeons must play a vital role if these programs are to succeed and benefit our patients.

Salivary Extracellular Vesicle-Associated exRNA as Cancer Biomarker

Extracellular vesicles (EVs) secreted in biological fluids contain several transcripts of the cell of origin, which may modify the functions and phenotype of proximal and distant cells. Cancer-derived EVs may promote a favorable microenvironment for cancer growth and invasion by acting on stroma and endothelial cells and may favor metastasis formation. The transcripts contained in cancer EVs may be exploited as biomarkers. Protein and extracellular RNA (exRNA) profiling in patient bio-fluids, such as blood and urine, was performed to identify molecular features with potential diagnostic and prognostic values. EVs are concentrated in saliva, and salivary EVs are particularly enriched in exRNAs. Several studies were focused on salivary EVs for the detection of biomarkers either of non-oral or oral cancers. The present paper provides an overview of the available studies on the diagnostic potential of exRNA profiling in salivary EVs.

Detection of novel fusion-transcripts by RNA-Seq in T-cell lymphoblastic lymphoma

Fusions transcripts have been proven to be strong drivers for neoplasia-associated mutations, although their incidence in T-cell lymphoblastic lymphoma needs to be determined yet. Using RNA-Seq we have selected 55 fusion transcripts identified by at least two of three detection methods in the same tumour. We confirmed the existence of 24 predicted novel fusions that had not been described in cancer or normal tissues yet, indicating the accuracy of the prediction. Of note, one of them involves the proto oncogene TAL1. Other confirmed fusions could explain the overexpression of driver genes such as COMMD3-BMI1, LMO1 or JAK3. Five fusions found exclusively in tumour samples could be considered pathogenic (NFYG-TAL1, RIC3-TCRBC2, SLC35A3-HIAT1, PICALM MLLT10 and MLLT10-PICALM). However, other fusions detected simultaneously in normal and tumour samples (JAK3-INSL3, KANSL1-ARL17A/B and TFG-ADGRG7) could be germ-line fusions genes involved in tumour-maintaining tasks. Notably, some fusions were confirmed in more tumour samples than predicted, indicating that the detection methods underestimated the real number of existing fusions. Our results highlight the potential of RNA-Seq to identify new cryptic fusions, which could be drivers or tumour-maintaining passenger genes. Such novel findings shed light on the searching for new T-LBL biomarkers in these haematological disorders.

WWOX Tumor Suppressor Gene in Breast Cancer, a Historical Perspective and Future Directions

The WWOX tumor suppressor gene is located at 16q23. 1–23.2, which covers the region of FRA16D—a common fragile sites. Deletions within the WWOX coding sequence are observed in up to 80% of breast cancer cases, which makes it one of the most common genetic alterations in this tumor type. The WWOX gene is known to play a role in breast cancer: increased expression of WWOX inhibits cell proliferation in suspension, reduces tumor growth rates in xenographic transplants, but also enhances cell migration through the basal membrane and contributes to morphological changes in 3D matrix-based cell cultures. The WWOX protein may act in several ways, as it has three functional domains—two WW domains, responsible for protein-protein interactions and an SDR domain (short dehydrogenase/reductase domain) which catalyzes conversions of low molecular weight ligands, most likely steroids. In epithelial cells, WWOX modulates gene transcription through interaction with p73, AP-2γ, and ERBB4 proteins. In steroid hormone-regulated tissues like mammary gland epithelium, the WWOX SDR domain acts as a steroid dehydrogenase. The relationship between WWOX and hormone receptors was shown in an animal model, where WWOX(C3H)+/–mice exhibited loss of both ER and PR receptors. Moreover, in breast cancer specimens, a positive correlation was observed between WWOX expression and ER status. On the other hand, decreased WWOX expression was associated with worse prognosis, namely higher relapse and mortality rates in BC patients. Recently, it was shown that genomic instability might be driven by the loss of WWOX expression. It was reported that WWOX plays role in DNA damage response (DDR) and DNA repair by regulating ATM activation through physical interaction. A genome caretaker function has also been proposed for WWOX, as it was found that WWOX sufficiency decreases homology directed repair (HDR) and supports non-homologous end-joining (NHEJ) repair as the dominant DSB repair pathway by Brca1-Wwox interaction. In breast cancer cells, WWOX was also found to modulate the expression of glycolysis pathway genes, through hypoxia-inducible transcription factor 1α (HIF1α) regulation. The paper presents the current state of knowledge regarding the WWOX tumor suppressor gene in breast cancer, as well as future research perspectives.

Tracing the evolution of the heterotrimeric G protein α subunit in Metazoa

BACKGROUND

Heterotrimeric G proteins are fundamental signaling proteins composed of three subunits, Gα and a Gβγ dimer. The role of Gα as a molecular switch is critical for transmitting and amplifying intracellular signaling cascades initiated by an activated G protein Coupled Receptor (GPCR). Despite their biochemical and therapeutic importance, the study of G protein evolution has been limited to the scope of a few model organisms. Furthermore, of the five primary Gα subfamilies, the underlying gene structure of only two families has been thoroughly investigated outside of Mammalia evolution. Therefore our understanding of Gα emergence and evolution across phylogeny remains incomplete.

RESULTS

We have computationally identified the presence and absence of every Gα gene (GNA-) across all major branches of Deuterostomia and evaluated the conservation of the underlying exon-intron structures across these phylogenetic groups. We provide evidence of mutually exclusive exon inclusion through alternative splicing in specific lineages. Variations of splice site conservation and isoforms were found for several paralogs which coincide with conserved, putative motifs of DNA-/RNA-binding proteins. In addition to our curated gene annotations, within Primates, we identified 15 retrotranspositions, many of which have undergone pseudogenization. Most importantly, we find numerous deviations from previous findings regarding the presence and absence of individual GNA- genes, nuanced differences in phyla-specific gene copy numbers, novel paralog duplications and subsequent intron gain and loss events.

CONCLUSIONS

Our curated annotations allow us to draw more accurate inferences regarding the emergence of all Gα family members across Metazoa and to present a new, updated theory of Gα evolution. Leveraging this, our results are critical for gaining new insights into the co-evolution of the Gα subunit and its many protein binding partners, especially therapeutically relevant G protein - GPCR signaling pathways which radiated in Vertebrata evolution.

Genomic hotspots but few recurrent fusion genes in breast cancer

The advent of next generation sequencing technologies has boosted the interest in exploring the role of fusion genes in the development and progression of solid tumors. In breast cancer, most of the detected gene fusions seem to be "passenger" events while the presence of recurrent and driver fusions is still under study. We performed RNA sequencing in 55 well-characterized breast cancer samples and 10 adjacent normal breast tissues, complemented by an analysis of SNP array data. We explored the presence of fusion genes and defined their association with breast cancer subtypes, clinical-pathologic characteristics and copy number aberrations. Overall, 370 fusions were detected across the majority of the samples. HER2+ samples had significantly more fusions than triple negative and luminal subtypes. The number of fusions was correlated with histological grade, Ki67 and tumor size. Clusters of fusion genes were observed across the genome and a significant correlation of fusions with copy number aberrations and more specifically amplifications was also revealed. Despite the large number of fusion events, only a few were recurrent, while recurrent individual genes forming fusions with different partners were also detected including the estrogen receptor 1 gene in the previously detected ESR1-CCDC170 fusion. Overall we detected novel gene fusion events while we confirmed previously reported fusions. Genomic hotspots of fusion genes, differences between subtypes and small number of recurrent fusions are the most relevant characteristics of these events in breast cancer. Further investigation is necessary to comprehend the biological significance of these fusions.

Towards precision medicine: discovering novel gynecological cancer biomarkers and pathways using linked data

BACKGROUND

Next Generation Sequencing (NGS) is playing a key role in therapeutic decision making for the cancer prognosis and treatment. The NGS technologies are producing a massive amount of sequencing datasets. Often, these datasets are published from the isolated and different sequencing facilities. Consequently, the process of sharing and aggregating multisite sequencing datasets are thwarted by issues such as the need to discover relevant data from different sources, built scalable repositories, the automation of data linkage, the volume of the data, efficient querying mechanism, and information rich intuitive visualisation.

RESULTS

We present an approach to link and query different sequencing datasets (TCGA, COSMIC, REACTOME, KEGG and GO) to indicate risks for four cancer types - Ovarian Serous Cystadenocarcinoma (OV), Uterine Corpus Endometrial Carcinoma (UCEC), Uterine Carcinosarcoma (UCS), Cervical Squamous Cell Carcinoma and Endocervical Adenocarcinoma (CESC) - covering the 16 healthy tissue-specific genes from Illumina Human Body Map 2.0. The differentially expressed genes from Illumina Human Body Map 2.0 are analysed together with the gene expressions reported in COSMIC and TCGA repositories leading to the discover of potential biomarkers for a tissue-specific cancer.

CONCLUSION

We analyse the tissue expression of genes, copy number variation (CNV), somatic mutation, and promoter methylation to identify associated pathways and find novel biomarkers. We discovered twenty (20) mutated genes and three (3) potential pathways causing promoter changes in different gynaecological cancer types. We propose a data-interlinked platform called BIOOPENER that glues together heterogeneous cancer and biomedical repositories. The key approach is to find correspondences (or data links) among genetic, cellular and molecular features across isolated cancer datasets giving insight into cancer progression from normal to diseased tissues. The proposed BIOOPENER platform enriches mutations by filling in missing links from TCGA, COSMIC, REACTOME, KEGG and GO datasets and provides an interlinking mechanism to understand cancer progression from normal to diseased tissues with pathway components, which in turn helped to map mutations, associated phenotypes, pathways, and mechanism.

De novo assembly and characterization of breast cancer transcriptomes identifies large numbers of novel fusion-gene transcripts of potential functional significance

Background

Gene-fusion or chimeric transcripts have been implicated in the onset and progression of a variety of cancers. Massively parallel RNA sequencing (RNA-Seq) of the cellular transcriptome is a promising approach for the identification of chimeric transcripts of potential functional significance. We report here the development and use of an integrated computational pipeline for the de novo assembly and characterization of chimeric transcripts in 55 primary breast cancer and normal tissue samples.

Methods

An integrated computational pipeline was employed to screen the transcriptome of breast cancer and control tissues for high-quality RNA-sequencing reads. Reads were de novo assembled into contigs followed by reference genome mapping. Chimeric transcripts were detected, filtered and characterized using our R-SAP algorithm. The relative abundance of reads was used to estimate levels of gene expression.

Results

De novo assembly allowed for the accurate detection of 1959 chimeric transcripts to nucleotide level resolution and facilitated detailed molecular characterization and quantitative analysis. A number of the chimeric transcripts are of potential functional significance including 79 novel fusion-protein transcripts and many chimeric transcripts with alterations in their un-translated leader regions. A number of chimeric transcripts in the cancer samples mapped to genomic regions devoid of any known genes. Several ‘pro-neoplastic’ fusions comprised of genes previously implicated in cancer are expressed at low levels in normal tissues but at high levels in cancer tissues.

Conclusions

Collectively, our results underscore the utility of deep sequencing technologies and improved bioinformatics workflows to uncover novel and potentially significant chimeric transcripts in cancer and normal somatic tissues.

Electronic supplementary material

The online version of this article (doi:10.1186/s12920-017-0289-7) contains supplementary material, which is available to authorized users.

Linked read sequencing resolves complex genomic rearrangements in gastric cancer metastases

Background

Genome rearrangements are critical oncogenic driver events in many malignancies. However, the identification and resolution of the structure of cancer genomic rearrangements remain challenging even with whole genome sequencing.

Methods

To identify oncogenic genomic rearrangements and resolve their structure, we analyzed linked read sequencing. This approach relies on a microfluidic droplet technology to produce libraries derived from single, high molecular weight DNA molecules, 50 kb in size or greater. After sequencing, the barcoded sequence reads provide long range genomic information, identify individual high molecular weight DNA molecules, determine the haplotype context of genetic variants that occur across contiguous megabase-length segments of the genome and delineate the structure of complex rearrangements. We applied linked read sequencing of whole genomes to the analysis of a set of synchronous metastatic diffuse gastric cancers that occurred in the same individual.

Results

When comparing metastatic sites, our analysis implicated a complex somatic rearrangement that was present in the metastatic tumor. The oncogenic event associated with the identified complex rearrangement resulted in an amplification of the known cancer driver gene FGFR2. With further investigation using these linked read data, the FGFR2 copy number alteration was determined to be a deletion-inversion motif that underwent tandem duplication, with unique breakpoints in each metastasis. Using a three-dimensional organoid tissue model, we functionally validated the metastatic potential of an FGFR2 amplification in gastric cancer.

Conclusions

Our study demonstrates that linked read sequencing is useful in characterizing oncogenic rearrangements in cancer metastasis.

Electronic supplementary material

The online version of this article (doi:10.1186/s13073-017-0447-8) contains supplementary material, which is available to authorized users.

RNA Biomarkers: Frontier of Precision Medicine for Cancer

As an essential part of central dogma, RNA delivers genetic and regulatory information and reflects cellular states. Based on high-throughput sequencing technologies, cumulating data show that various RNA molecules are able to serve as biomarkers for the diagnosis and prognosis of various diseases, for instance, cancer. In particular, detectable in various bio-fluids, such as serum, saliva and urine, extracellular RNAs (exRNAs) are emerging as non-invasive biomarkers for earlier cancer diagnosis, tumor progression monitor, and prediction of therapy response. In this review, we summarize the latest studies on various types of RNA biomarkers, especially extracellular RNAs, in cancer diagnosis and prognosis, and illustrate several well-known RNA biomarkers of clinical utility. In addition, we describe and discuss general procedures and issues in investigating exRNA biomarkers, and perspectives on utility of exRNAs in precision medicine.

Functions of bromodomain-containing proteins and their roles in homeostasis and cancer

Bromodomains (BRDs) are evolutionarily conserved protein-protein interaction modules that are found in a wide range of proteins with diverse catalytic and scaffolding functions and are present in most tissues. BRDs selectively recognize and bind to acetylated Lys residues - particularly in histones - and thereby have important roles in the regulation of gene expression. BRD-containing proteins are frequently dysregulated in cancer, they participate in gene fusions that generate diverse, frequently oncogenic proteins, and many cancer-causing mutations have been mapped to the BRDs themselves. Importantly, BRDs can be targeted by small-molecule inhibitors, which has stimulated many translational research projects that seek to attenuate the aberrant functions of BRD-containing proteins in disease.

Integrated genomic analyses identify frequent gene fusion events and VHL inactivation in gastrointestinal stromal tumors

Gastrointestinal stromal tumors (GISTs) are the most common mesenchymal tumors of the gastrointestinal tract. We sequenced nine exomes and transcriptomes, and two genomes of GISTs for integrated analyses. We detected 306 somatic variants in nine GISTs and recurrent protein-altering mutations in 29 genes. Transcriptome sequencing revealed 328 gene fusions, and the most frequently involved fusion events were associated with IGF2 fused to several partner genes including CCND1, FUS, and LASP1. We additionally identified three recurrent read-through fusion transcripts: POLA2-CDC42EP2, C8orf42-FBXO25, and STX16-NPEPL1. Notably, we found intragenic deletions in one of three exons of the VHL gene and increased mRNAs of VEGF, PDGF-β, and IGF-1/2 in 56% of GISTs, suggesting a mechanistic link between VHL inactivation and overexpression of hypoxia-inducible factor target genes in the absence of hypoxia. We also identified copy number gain and increased mRNA expression of AMACR, CRIM1, SKP2, and CACNA1E. Mapping of copy number and gene expression results to the KEGG pathways revealed activation of the JAK-STAT pathway in small intestinal GISTs and the MAPK pathway in wild-type GISTs. These observations will allow us to determine the genetic basis of GISTs and will facilitate further investigation to develop new therapeutic options.

Next generation sequencing approach for detecting 491 fusion genes from human cancer

Next-generation DNA sequencing (NGS) of the genomes of cancer cells is contributing to new discoveries that illuminate the mechanisms of tumorigenesis. To this end, the International Cancer Genome Consortium and The Cancer Genome Atlas are investigating novel alterations of genes that will define the pathways and mechanisms of the development and growth of cancers. These efforts contribute to the development of innovative pharmaceuticals as well as to the introduction of genome sequencing as a component of personalized medicine. In particular, chromosomal translocations that fuse coding sequences serve as important pharmaceutical targets and diagnostic markers given their association with tumorigenesis. Although increasing numbers of fusion genes are being discovered using NGS, the methodology used to identify such fusion genes is complicated, expensive, and requires relatively large samples. Here, to address these problems, we describe the design and development of a panel of 491 fusion genes that performed well in the analysis of cultured human cancer cell lines and 600 clinical tumor specimens.

ChiTaRS-3.1-the enhanced chimeric transcripts and RNA-seq database matched with protein-protein interactions

Discovery of chimeric RNAs, which are produced by chromosomal translocations as well as the joining of exons from different genes by trans-splicing, has added a new level of complexity to our study and understanding of the transcriptome. The enhanced ChiTaRS-3.1 database (http://chitars.md.biu.ac.il) is designed to make widely accessible a wealth of mined data on chimeric RNAs, with easy-to-use analytical tools built-in. The database comprises 34 922 chimeric transcripts along with 11 714 cancer breakpoints. In this latest version, we have included multiple cross-references to GeneCards, iHop, PubMed, NCBI, Ensembl, OMIM, RefSeq and the Mitelman collection for every entry in the ‘Full Collection’. In addition, for every chimera, we have added a predicted chimeric protein–protein interaction (ChiPPI) network, which allows for easy visualization of protein partners of both parental and fusion proteins for all human chimeras. The database contains a comprehensive annotation for 34 922 chimeric transcripts from eight organisms, and includes the manual annotation of 200 sense-antiSense (SaS) chimeras. The current improvements in the content and functionality to the ChiTaRS database make it a central resource for the study of chimeric transcripts and fusion proteins.

Targeted Cancer Therapy: Vital Oncogenes and a New Molecular Genetic Paradigm for Cancer Initiation Progression and Treatment

It has been declared repeatedly that cancer is a result of molecular genetic abnormalities. However, there has been no working model describing the specific functional consequences of the deranged genomic processes that result in the initiation and propagation of the cancer process during carcinogenesis. We no longer need to question whether or not cancer arises as a result of a molecular genetic defect within the cancer cell. The legitimate questions are: how and why? This article reviews the preeminent data on cancer molecular genetics and subsequently proposes that the sentinel event in cancer initiation is the aberrant production of fused transcription activators with new molecular properties within normal tissue stem cells. This results in the production of vital oncogenes with dysfunctional gene activation transcription properties, which leads to dysfunctional gene regulation, the aberrant activation of transduction pathways, chromosomal breakage, activation of driver oncogenes, reactivation of stem cell transduction pathways and the activation of genes that result in the hallmarks of cancer. Furthermore, a novel holistic molecular genetic model of cancer initiation and progression is presented along with a new paradigm for the approach to personalized targeted cancer therapy, clinical monitoring and cancer diagnosis.

MXRA5 is a TGF-β1-regulated human protein with anti-inflammatory and anti-fibrotic properties

Current therapy for chronic kidney disease (CKD) is unsatisfactory because of an insufficient understanding of its pathogenesis. Matrix remodelling-associated protein 5 (MXRA5, adlican) is a human protein of unknown function with high kidney tissue expression, not present in rodents. Given the increased expression of MXRA5 in injured tissues, including the kidneys, we have suggested that MXRA5 may modulate kidney injury. MXRA5 immunoreactivity was observed in tubular cells in human renal biopsies and in urine from CKD patients. We then explored factors regulating MXRA5 expression and MXRA5 function in cultured human proximal tubular epithelial cells and explored MXRA5 expression in kidney cancer cells and kidney tissue. The fibrogenic cytokine transforming growth factor-β1 (TGFβ1) up-regulated MXRA5 mRNA and protein expression. TGFβ1-induced MXRA5 up-regulation was prevented by either interference with TGFβ1 activation of the TGFβ receptor 1 (TGFBR1, ALK5) or by the vitamin D receptor agonist paricalcitol. By contrast, the pro-inflammatory cytokine TWEAK did not modulate MXRA5 expression. MXRA5 siRNA-induced down-regulation of constitutive MXRA5 expression resulted in higher TWEAK-induced expression of chemokines. In addition, MXRA5 down-regulation resulted in a magnified expression of genes encoding extracellular matrix proteins in response to TGFβ1. Furthermore, in clear cell renal cancer, von Hippel-Lindau (VHL) regulated MXRA5 expression. In conclusion, MXRA5 is a TGFβ1- and VHL-regulated protein and, for the first time, we identify MXRA5 functions as an anti-inflammatory and anti-fibrotic molecule. This information may yield clues to design novel therapeutic strategies in diseases characterized by inflammation and fibrosis.

Deciphering and Targeting Oncogenic Mutations and Pathways in Breast Cancer

: Advances in DNA and RNA sequencing revealed substantially greater genomic complexity in breast cancer than simple models of a few driver mutations would suggest. Only very few, recurrent mutations or copy-number variations in cancer-causing genes have been identified. The two most common alterations in breast cancer are TP53 (affecting the majority of triple-negative breast cancers) and PIK3CA (affecting almost half of estrogen receptor-positive cancers) mutations, followed by a long tail of individually rare mutations affecting <1%-20% of cases. Each cancer harbors from a few dozen to a few hundred potentially high-functional impact somatic variants, along with a much larger number of potentially high-functional impact germline variants. It is likely that it is the combined effect of all genomic variations that drives the clinical behavior of a given cancer. Furthermore, entirely new classes of oncogenic events are being discovered in the noncoding areas of the genome and in noncoding RNA species driven by errors in RNA editing. In light of this complexity, it is not unexpected that, with the exception of HER2 amplification, no robust molecular predictors of benefit from targeted therapies have been identified. In this review, we summarize the current genomic portrait of breast cancer, focusing on genetic aberrations that are actively being targeted with investigational drugs.

IMPLICATIONS FOR PRACTICE

Next-generation sequencing is now widely available in the clinic, but interpretation of the results is challenging, and its impact on treatment selection is often limited. This work provides an overview of frequently encountered molecular abnormalities in breast cancer and discusses their potential therapeutic implications. This review emphasizes the importance of administering investigational targeted therapies, or off-label use of approved targeted drugs, in the context of a formal clinical trial or registry programs to facilitate learning about the clinical utility of tumor target profiling.

Diverse, Biologically Relevant, and Targetable Gene Rearrangements in Triple-Negative Breast Cancer and Other Malignancies

Triple-negative breast cancer (TNBC) and other molecularly heterogeneous malignancies present a significant clinical challenge due to a lack of high-frequency "driver" alterations amenable to therapeutic intervention. These cancers often exhibit genomic instability, resulting in chromosomal rearrangements that affect the structure and expression of protein-coding genes. However, identification of these rearrangements remains technically challenging. Using a newly developed approach that quantitatively predicts gene rearrangements in tumor-derived genetic material, we identified and characterized a novel oncogenic fusion involving the MER proto-oncogene tyrosine kinase (MERTK) and discovered a clinical occurrence and cell line model of the targetable FGFR3-TACC3 fusion in TNBC. Expanding our analysis to other malignancies, we identified a diverse array of novel and known hybrid transcripts, including rearrangements between noncoding regions and clinically relevant genes such as ALK, CSF1R, and CD274/PD-L1 The over 1,000 genetic alterations we identified highlight the importance of considering noncoding gene rearrangement partners, and the targetable gene fusions identified in TNBC demonstrate the need to advance gene fusion detection for molecularly heterogeneous cancers. Cancer Res; 76(16); 4850-60. ©2016 AACR.

Recurrent hormone-binding domain truncated ESR1 amplifications in primary endometrial cancers suggest their implication in hormone independent growth

The estrogen receptor alpha (ERα) is highly expressed in both endometrial and breast cancers, and represents the most prevalent therapeutic target in breast cancer. However, anti-estrogen therapy has not been shown to be effective in endometrial cancer. Recently it has been shown that hormone-binding domain alterations of ERα in breast cancer contribute to acquired resistance to anti-estrogen therapy. In analyses of genomic data from The Cancer Genome Atlas (TCGA), we observe that endometrial carcinomas manifest recurrent ESR1 gene amplifications that truncate the hormone-binding domain encoding region of ESR1 and are associated with reduced mRNA expression of exons encoding the hormone-binding domain. These findings support a role for hormone-binding alterations of ERα in primary endometrial cancer, with potentially important therapeutic implications.

Dissecting the expression relationships between RNA-binding proteins and their cognate targets in eukaryotic post-transcriptional regulatory networks

RNA-binding proteins (RBPs) are pivotal in orchestrating several steps in the metabolism of RNA in eukaryotes thereby controlling an extensive network of RBP-RNA interactions. Here, we employed CLIP (cross-linking immunoprecipitation)-seq datasets for 60 human RBPs and RIP-ChIP (RNP immunoprecipitation-microarray) data for 69 yeast RBPs to construct a network of genome-wide RBP- target RNA interactions for each RBP. We show in humans that majority (~78%) of the RBPs are strongly associated with their target transcripts at transcript level while ~95% of the studied RBPs were also found to be strongly associated with expression levels of target transcripts when protein expression levels of RBPs were employed. At transcript level, RBP - RNA interaction data for the yeast genome, exhibited a strong association for 63% of the RBPs, confirming the association to be conserved across large phylogenetic distances. Analysis to uncover the features contributing to these associations revealed the number of target transcripts and length of the selected protein-coding transcript of an RBP at the transcript level while intensity of the CLIP signal, number of RNA-Binding domains, location of the binding site on the transcript, to be significant at the protein level. Our analysis will contribute to improved modelling and prediction of post-transcriptional networks.

Comparative assessment of methods for the fusion transcripts detection from RNA-Seq data

RNA-Seq made possible the global identification of fusion transcripts, i.e. "chimeric RNAs". Even though various software packages have been developed to serve this purpose, they behave differently in different datasets provided by different developers. It is important for both users, and developers to have an unbiased assessment of the performance of existing fusion detection tools. Toward this goal, we compared the performance of 12 well-known fusion detection software packages. We evaluated the sensitivity, false discovery rate, computing time, and memory usage of these tools in four different datasets (positive, negative, mixed, and test). We conclude that some tools are better than others in terms of sensitivity, positive prediction value, time consumption and memory usage. We also observed small overlaps of the fusions detected by different tools in the real dataset (test dataset). This could be due to false discoveries by various tools, but could also be due to the reason that none of the tools are inclusive. We have found that the performance of the tools depends on the quality, read length, and number of reads of the RNA-Seq data. We recommend that users choose the proper tools for their purpose based on the properties of their RNA-Seq data.

INTEGRATE: gene fusion discovery using whole genome and transcriptome data

While next-generation sequencing (NGS) has become the primary technology for discovering gene fusions, we are still faced with the challenge of ensuring that causative mutations are not missed while minimizing false positives. Currently, there are many computational tools that predict structural variations (SV) and gene fusions using whole genome (WGS) and transcriptome sequencing (RNA-seq) data separately. However, as both WGS and RNA-seq have their limitations when used independently, we hypothesize that the orthogonal validation from integrating both data could generate a sensitive and specific approach for detecting high-confidence gene fusion predictions. Fortunately, decreasing NGS costs have resulted in a growing quantity of patients with both data available. Therefore, we developed a gene fusion discovery tool, INTEGRATE, that leverages both RNA-seq and WGS data to reconstruct gene fusion junctions and genomic breakpoints by split-read mapping. To evaluate INTEGRATE, we compared it with eight additional gene fusion discovery tools using the well-characterized breast cell line HCC1395 and peripheral blood lymphocytes derived from the same patient (HCC1395BL). The predictions subsequently underwent a targeted validation leading to the discovery of 131 novel fusions in addition to the seven previously reported fusions. Overall, INTEGRATE only missed six out of the 138 validated fusions and had the highest accuracy of the nine tools evaluated. Additionally, we applied INTEGRATE to 62 breast cancer patients from The Cancer Genome Atlas (TCGA) and found multiple recurrent gene fusions including a subset involving estrogen receptor. Taken together, INTEGRATE is a highly sensitive and accurate tool that is freely available for academic use.

The Impact of Normalization Methods on RNA-Seq Data Analysis

High-throughput sequencing technologies, such as the Illumina Hi-seq, are powerful new tools for investigating a wide range of biological and medical problems. Massive and complex data sets produced by the sequencers create a need for development of statistical and computational methods that can tackle the analysis and management of data. The data normalization is one of the most crucial steps of data processing and this process must be carefully considered as it has a profound effect on the results of the analysis. In this work, we focus on a comprehensive comparison of five normalization methods related to sequencing depth, widely used for transcriptome sequencing (RNA-seq) data, and their impact on the results of gene expression analysis. Based on this study, we suggest a universal workflow that can be applied for the selection of the optimal normalization procedure for any particular data set. The described workflow includes calculation of the bias and variance values for the control genes, sensitivity and specificity of the methods, and classification errors as well as generation of the diagnostic plots. Combining the above information facilitates the selection of the most appropriate normalization method for the studied data sets and determines which methods can be used interchangeably.

Common Oncogene Mutations and Novel SND1-BRAF Transcript Fusion in Lung Adenocarcinoma from Never Smokers

Lung adenocarcinomas from never smokers account for approximately 15 to 20% of all lung cancers and these tumors often carry genetic alterations that are responsive to targeted therapy. Here we examined mutation status in 10 oncogenes among 89 lung adenocarcinomas from never smokers. We also screened for oncogene fusion transcripts in 20 of the 89 tumors by RNA-Seq. In total, 62 tumors had mutations in at least one of the 10 oncogenes, including EGFR (49 cases, 55%), K-ras (5 cases, 6%), BRAF (4 cases, 5%), PIK3CA (3 cases, 3%), and ERBB2 (4 cases, 5%). In addition to ALK fusions identified by IHC/FISH in four cases, two previously known fusions involving EZR- ROS1 and KIF5B-RET were identified by RNA-Seq as well as a third novel fusion transcript that was formed between exons 1–9 of SND1 and exons 2 to 3′ end of BRAF. This in-frame fusion was observed in 3/89 tested tumors and 2/64 additional never smoker lung adenocarcinoma samples. Ectopic expression of SND1-BRAF in H1299 cells increased phosphorylation levels of MEK/ERK, cell proliferation, and spheroid formation compared to parental mock-transfected control. Jointly, our results suggest a potential role of the novel BRAF fusion in lung cancer development and therapy.

Polycomb repressive complex 2 epigenomic signature defines age-associated hypermethylation and gene expression changes

Although age-associated gene expression and methylation changes have been reported throughout the literature, the unifying epigenomic principles of aging remain poorly understood. Recent explosion in availability and resolution of functional/regulatory genome annotation data (epigenomic data), such as that provided by the ENCODE and Roadmap Epigenomics projects, provides an opportunity for the identification of epigenomic mechanisms potentially altered by age-associated differentially methylated regions (aDMRs) and regulatory signatures in the promoters of age-associated genes (aGENs). In this study we found that aDMRs and aGENs identified in multiple independent studies share a common Polycomb Repressive Complex 2 signature marked by EZH2, SUZ12, CTCF binding sites, repressive H3K27me3, and activating H3K4me1 histone modification marks, and a "poised promoter" chromatin state. This signature is depleted in RNA Polymerase II-associated transcription factor binding sites, activating H3K79me2, H3K36me3, H3K27ac marks, and an "active promoter" chromatin state. The PRC2 signature was shown to be generally stable across cell types. When considering the directionality of methylation changes, we found the PRC2 signature to be associated with aDMRs hypermethylated with age, while hypomethylated aDMRs were associated with enhancers. In contrast, aGENs were associated with the PRC2 signature independently of the directionality of gene expression changes. In this study we demonstrate that the PRC2 signature is the common epigenomic context of genomic regions associated with hypermethylation and gene expression changes in aging.

Assessing the consistency of public human tissue RNA-seq data sets

Sequencing-based gene expression methods like RNA-sequencing (RNA-seq) have become increasingly common, but it is often claimed that results obtained in different studies are not comparable owing to the influence of laboratory batch effects, differences in RNA extraction and sequencing library preparation methods and bioinformatics processing pipelines. It would be unfortunate if different experiments were in fact incomparable, as there is great promise in data fusion and meta-analysis applied to sequencing data sets. We therefore compared reported gene expression measurements for ostensibly similar samples (specifically, human brain, heart and kidney samples) in several different RNA-seq studies to assess their overall consistency and to examine the factors contributing most to systematic differences. The same comparisons were also performed after preprocessing all data in a consistent way, eliminating potential bias from bioinformatics pipelines. We conclude that published human tissue RNA-seq expression measurements appear relatively consistent in the sense that samples cluster by tissue rather than laboratory of origin given simple preprocessing transformations. The article is supplemented by a detailed walkthrough with embedded R code and figures.

Functional characterization of BC039389-GATM and KLK4-KRSP1 chimeric read-through transcripts which are up-regulated in renal cell cancer

Background

Chimeric read-through RNAs are transcripts originating from two directly adjacent genes (<10 kb) on the same DNA strand. Although they are found in next-generation whole transcriptome sequencing (RNA-Seq) data on a regular basis, investigating them further has usually been refrained from. Therefore, their expression patterns or functions in general, and in oncogenesis in particular, are poorly understood.

Results

We used paired-end RNA-Seq and a specifically designed computational data analysis pipeline (FusionSeq) to nominate read-through events in a small discovery set of renal cell carcinomas (RCC) and confirmed them in a larger validation cohort.

324 read-through events were called overall; 22/27 (81%) selected nominees passed validation with conventional PCR and were sequenced at the junction region. We frequently identified various isoforms of a given read-through event. 2/22 read-throughs were up-regulated: BC039389-GATM was higher expressed in RCC compared to benign adjacent kidney; KLK4-KRSP1 was expressed in 46/169 (27%) RCCs, but rarely in normal tissue. KLK4-KRSP1 expression was associated with worse clinical outcome in the patient cohort. In cell lines, both read-throughs influenced molecular mechanisms (i.e. target gene expression or migration/invasion) in a way that counteracted the effect of the respective parent transcript GATM or KLK4.

Conclusions

Our data suggests that the up-regulation of read-through RNA chimeras in tumors is not random but causes regulatory effects on cellular mechanisms and may impact patient survival.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1446-z) contains supplementary material, which is available to authorized users.

Genome-wide interrogation of longitudinal FEV1 in children with asthma

RATIONALE

Most genomic studies of lung function have used phenotypic data derived from a single time-point (e.g., presence/absence of disease) without considering the dynamic progression of a chronic disease.

OBJECTIVES

To characterize lung function change over time in subjects with asthma and identify genetic contributors to a longitudinal phenotype.

METHODS

We present a method that models longitudinal FEV1 data, collected from 1,041 children with asthma who participated in the Childhood Asthma Management Program. This longitudinal progression model was built using population-based nonlinear mixed-effects modeling with an exponential structure and the determinants of age and height.

MEASUREMENTS AND MAIN RESULTS

We found ethnicity was a key covariate for FEV1 level. Budesonide-treated children with asthma had a slight but significant effect on FEV1 when compared with those treated with placebo or nedocromil (P < 0.001). A genome-wide association study identified seven single-nucleotide polymorphisms nominally associated with longitudinal lung function phenotypes in 581 white Childhood Asthma Management Program subjects (P < 10(-4) in the placebo ["discovery"] and P < 0.05 in the nedocromil treatment ["replication"] group). Using ChIP-seq and RNA-seq data, we found that some of the associated variants were in strong enhancer regions in human lung fibroblasts and may affect gene expression in human lung tissue. Genetic mapping restricted to genome-wide enhancer single-nucleotide polymorphisms in lung fibroblasts revealed a highly significant variant (rs6763931; P = 4 × 10(-6); false discovery rate < 0.05).

CONCLUSIONS

This study offers a strategy to explore the genetic determinants of longitudinal phenotypes, provide a comprehensive picture of disease pathophysiology, and suggest potential treatment targets.