Bimodal gene expression patterns in breast cancer

Subgroup identification of targeted therapy effects on biomarker for time to event data

BACKGROUND

The initiation biomarker-driven trials have revolutionized oncology drug development by challenging the traditional phased approach and introducing basket studies. Notable successes in non-small cell lung cancer (NSCLC) with ALK, ALK/ROS1, and EGFR inhibitors have prompted the need to expand this approach to other cancer sites.

OBJECTIVES

This study explores the use of dose response modeling and time-to-event algorithms on the biomarker molecular targeted agent (MTA). By simulating subgroup identification in MTA-related time-to-event data, the study aims to develop statistical methodology supporting biomarker-driven trials in oncology.

METHODS

A total of n patients are selected assigned for different doses. A dataset is prepared to mimic the situation on Subgroup Identification of MTA for time to event data analysis. The response is measured through MTA. The MTA value is also measured through ROC. The Markov Chain Monte Carlo (MCMC) techniques are prepared to perform the proposed algorithm. The analysis is carried out with a simulation study. The subset selection is performed through the Threshold Limit Value (TLV) by the Bayesian approach.

RESULTS

The MTA is observed with range 12-16. It is expected that there is a marginal level shift of the MTA from pre to post-treatment. The Cox time-varying model can be adopted further as causal-effect relation to establishing the MTA on prolonging the survival duration. The proposed work in the statistical methodology to support the biomarker-driven trial for oncology research.

CONCLUSION

This study extends the application of biomarker-driven trials beyond NSCLC, opening possibilities for implementation in other cancer sites. By demonstrating the feasibility and efficacy of utilizing MTA as a biomarker, the research lays the foundation for refining and validating biomarker use in clinical trials. These advancements aim to enhance the precision and effectiveness of cancer treatments, ultimately benefiting patients.

Bimodal gene expression in cancer patients provides interpretable biomarkers for drug sensitivity

Identifying biomarkers predictive of cancer cell response to drug treatment constitutes one of the main challenges in precision oncology. Recent large-scale cancer pharmacogenomic studies have opened new avenues of research to develop predictive biomarkers by profiling thousands of human cancer cell lines at the molecular level and screening them with hundreds of approved drugs and experimental chemical compounds. Many studies have leveraged these data to build predictive models of response using various statistical and machine learning methods. However, a common pitfall to these methods is the lack of interpretability as to how they make predictions, hindering the clinical translation of these models. To alleviate this issue, we used the recent logic modeling approach to develop a new machine learning pipeline that explores the space of bimodally expressed genes in multiple large in vitro pharmacogenomic studies and builds multivariate, nonlinear, yet interpretable logic-based models predictive of drug response. The performance of this approach was showcased in a compendium of the three largest in vitro pharmacogenomic data sets to build robust and interpretable models for 101 drugs that span 17 drug classes with high validation rates in independent datasets. These results along with in vivo and clinical validation, support a better translation of gene expression biomarkers between model systems using bimodal gene expression.

Silent transcription intervals and translational bursting lead to diverse phenotypic switching

For complex process of gene expression, we use theoretical analysis and stochastic simulations to study the phenotypic diversity induced by silent transcription intervals and translational bursting.

An integrated approach to identify bimodal genes associated with prognosis in câncer

Bimodal gene expression (where a gene expression distribution has two maxima) is associated with phenotypic diversity in different biological systems. A critical issue, thus, is the integration of expression and phenotype data to identify genuine associations. Here, we developed tools that allow both: i) the identification of genes with bimodal gene expression and ii) their association with prognosis in cancer patients from The Cancer Genome Atlas (TCGA). Bimodality was observed for 554 genes in expression data from 25 tumor types. Furthermore, 96 of these genes presented different prognosis when patients belonging to the two expression peaks were compared. The software to execute the method and the corresponding documentation are available at the Data access section.

Stochastic sequestration dynamics: a minimal model with extrinsic noise for bimodal distributions and competitors correlation

Many biological processes are known to be based on molecular sequestration. This kind of dynamics involves two types of molecular species, namely targets and sequestrants, that bind to form a complex. In the simple framework of mass-action law, key features of these systems appear to be threshold-like profiles of the amounts of free molecules as a function of the parameters determining their possible maximum abundance. However, biochemical processes are probabilistic and take place in stochastically fluctuating environments. How these different sources of noise affect the final outcome of the network is not completely characterised yet. In this paper we specifically investigate the effects induced by a source of extrinsic noise onto a minimal stochastic model of molecular sequestration. We analytically show how bimodal distributions of the targets can appear and characterise them as a result of noise filtering mediated by the threshold response. We then address the correlations between target species induced by the sequestrant and discuss how extrinsic noise can turn the negative correlation caused by competition into a positive one. Finally, we consider the more complex scenario of competitive inhibition for enzymatic kinetics and discuss the relevance of our findings with respect to applications.

On the role of extrinsic noise in microRNA-mediated bimodal gene expression

Several studies highlighted the relevance of extrinsic noise in shaping cell decision making and differentiation in molecular networks. Bimodal distributions of gene expression levels provide experimental evidence of phenotypic differentiation, where the modes of the distribution often correspond to different physiological states of the system. We theoretically address the presence of bimodal phenotypes in the context of microRNA (miRNA)-mediated regulation. MiRNAs are small noncoding RNA molecules that downregulate the expression of their target mRNAs. The nature of this interaction is titrative and induces a threshold effect: below a given target transcription rate almost no mRNAs are free and available for translation. We investigate the effect of extrinsic noise on the system by introducing a fluctuating miRNA-transcription rate. We find that the presence of extrinsic noise favours the presence of bimodal target distributions which can be observed for a wider range of parameters compared to the case with intrinsic noise only and for lower miRNA-target interaction strength. Our results suggest that combining threshold-inducing interactions with extrinsic noise provides a simple and robust mechanism for obtaining bimodal populations without requiring fine tuning. Furthermore, we characterise the protein distribution's dependence on protein half-life.

MCM2: An alternative to Ki-67 for measuring breast cancer cell proliferation

Breast cancer is a heterogeneous disease comprising a diversity of tumor subtypes that manifest themselves in a wide variety of clinical, pathological, and molecular features. One important subset, luminal breast cancers, comprises two clinically distinct subtypes luminal A and B each of them endowed with its own genetic program of differentiation and proliferation. Luminal breast cancers were operationally defined as follows: Luminal A: ER+, PR+, HER2-, Ki-67<14% and Luminal B: ER+ and/or PR+, HER2-,Ki-67≥14% or, alternatively ER+ and/or PR+, HER2+, any Ki-67. There is currently a need for a clinically robust and validated immunohistochemical assay that can help distinguish between luminal A and B breast cancer. MCM2 is a family member of the minichromosome maintenance protein complex whose role in DNA replication and cell proliferation is firmly established. As MCM2 appears to be an attractive alternative to Ki-67, we sought to study the expression of MCM2 and Ki-67 in different histological grades and molecular subtypes of breast cancer focusing primarily on ER-positive tumors. MCM2 and Ki-67 mRNA expression were studied using in silico analysis of available DNA microarray and RNA-sequencing data of human breast cancer. We next used immunohistochemistry to evaluate protein expression of MCM2 and Ki-67 on tissue microarrays of invasive breast carcinoma. We found that MCM2 and Ki-67 are highly expressed in breast tumors of high histological grades, comprising clinically aggressive tumors such as triple-negative, HER2-positive and luminal B subtypes. MCM2 expression was detected at higher levels than that of Ki-67 in normal breast tissues and in breast cancers. The bimodal distribution of MCM2 scores in ER+/HER2- breast tumors led to the identification of two distinct subgroups with different relapse-free survival rates. In conclusion, MCM2 expression can help sorting out two clinically important subsets of luminal breast cancer whose treatment and clinical outcomes are likely to diverge.

RNAs competing for microRNAs mutually influence their fluctuations in a highly non-linear microRNA-dependent manner in single cells

Background

Distinct RNA species may compete for binding to microRNAs (miRNAs). This competition creates an indirect interaction between miRNA targets, which behave as miRNA sponges and eventually influence each other’s expression levels. Theoretical predictions suggest that not only the mean expression levels of targets but also the fluctuations around the means are coupled through miRNAs. This may result in striking effects on a broad range of cellular processes, such as cell differentiation and proliferation. Although several studies have reported the functional relevance of this mechanism of interaction, detailed experiments are lacking that study this phenomenon in controlled conditions by mimicking a physiological range.

Results

We used an experimental design based on two bidirectional plasmids and flow cytometry measurements of cotransfected mammalian cells. We validated a stochastic gene interaction model that describes how mRNAs can influence each other’s fluctuations in a miRNA-dependent manner in single cells. We show that miRNA–target correlations eventually lead to either bimodal cell population distributions with high and low target expression states, or correlated fluctuations across targets when the pool of unbound targets and miRNAs are in near-equimolar concentration. We found that there is an optimal range of conditions for the onset of cross-regulation, which is compatible with 10–1000 copies of targets per cell.

Conclusions

Our results are summarized in a phase diagram for miRNA-mediated cross-regulation that links experimentally measured quantities and effective model parameters. This phase diagram can be applied to in vivo studies of RNAs that are in competition for miRNA binding.

Electronic supplementary material

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

Characteristic patterns of cancer incidence: Epidemiological data, biological theories, and multistage models

We investigate and classify several patterns in cancer incidence and relative risk data which persist across different countries and multiple published studies. We then explore biological hypotheses as well as many mathematical models in the literature that attempt to explain these patterns. A general modeling framework is presented which is general enough to model most of observed behaviors. It is our belief that this model has sufficient flexibility to be adapted to new information as it is discovered. As one application of this framework, we give a model for the effect of aging on the process of carcinogenesis.

Transcriptomic Biomarkers for Tuberculosis: Evaluation of DOCK9. EPHA4, and NPC2 mRNA Expression in Peripheral Blood

Lately, much effort has been made to find mRNA biomarkers for tuberculosis (TB) disease/infection with microarray-based approaches. In a pilot investigation, through RNA sequencing technology, we observed a prominent modulation of DOCK9, EPHA4, and NPC2 mRNA abundance in the blood of TB patients. To corroborate these findings, independent validations were performed in cohorts from different areas. Gene expression levels in blood were evaluated by quantitative real-time PCR (Brazil, n = 129) or reanalysis of public microarray data (UK: n = 96; South Africa: n = 51; Germany: n = 26; and UK/France: n = 63). In the Brazilian cohort, significant modulation of all target-genes was observed comparing TB vs. healthy recent close TB contacts (rCt). With a 92% specificity, NPC2 mRNA high expression (NPC2^high) showed the highest sensitivity (85%, 95% CI 65%–96%; area under the ROC curve [AUROC] = 0.88), followed by EPHA4 (53%, 95% CI 33%–73%, AUROC = 0.73) and DOCK9 (19%, 95% CI 7%–40%; AUROC = 0.66). All the other reanalyzed cohorts corroborated the potential of NPC2^high as a biomarker for TB (sensitivity: 82–100%; specificity: 94–97%). An NPC2^high profile was also observed in 60% (29/48) of the tuberculin skin test positive rCt, and additional follow-up evaluation revealed changes in the expression levels of NPC2 during the different stages of Mycobacterium tuberculosis infection, suggesting that further studies are needed to evaluate modulation of this gene during latent TB and/or progression to active disease. Considering its high specificity, our data indicate, for the first time, that NPC2^high might serve as an accurate single-gene biomarker for TB.

Outlier analysis of functional genomic profiles enriches for oncology targets and enables precision medicine

Background

Genome-scale functional genomic screens across large cell line panels provide a rich resource for discovering tumor vulnerabilities that can lead to the next generation of targeted therapies. Their data analysis typically has focused on identifying genes whose knockdown enhances response in various pre-defined genetic contexts, which are limited by biological complexities as well as the incompleteness of our knowledge. We thus introduce a complementary data mining strategy to identify genes with exceptional sensitivity in subsets, or outlier groups, of cell lines, allowing an unbiased analysis without any a priori assumption about the underlying biology of dependency.

Results

Genes with outlier features are strongly and specifically enriched with those known to be associated with cancer and relevant biological processes, despite no a priori knowledge being used to drive the analysis. Identification of exceptional responders (outliers) may not lead only to new candidates for therapeutic intervention, but also tumor indications and response biomarkers for companion precision medicine strategies. Several tumor suppressors have an outlier sensitivity pattern, supporting and generalizing the notion that tumor suppressors can play context-dependent oncogenic roles.

Conclusions

The novel application of outlier analysis described here demonstrates a systematic and data-driven analytical strategy to decipher large-scale functional genomic data for oncology target and precision medicine discoveries.

Electronic supplementary material

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

Multi-omic measurement of mutually exclusive loss-of-function enriches for candidate synthetic lethal gene pairs

Background

Identification of synthetic lethal interactions in cancer cells could offer promising new therapeutic targets. Large-scale functional genomic screening presents an opportunity to test large numbers of cancer synthetic lethal hypotheses. Methods enriching for candidate synthetic lethal targets in molecularly defined cancer cell lines can steer effective design of screening efforts. Loss of one partner of a synthetic lethal gene pair creates a dependency on the other, thus synthetic lethal gene pairs should never show simultaneous loss-of-function. We have developed a computational approach to mine large multi-omic cancer data sets and identify gene pairs with mutually exclusive loss-of-function. Since loss-of-function may not always be genetic, we look for deleterious mutations, gene deletion and/or loss of mRNA expression by bimodality defined with a novel algorithm BiSEp.

Results

Applying this toolkit to both tumour cell line and patient data, we achieve statistically significant enrichment for experimentally validated tumour suppressor genes and synthetic lethal gene pairings. Notably non-reliance on genetic loss reveals a number of known synthetic lethal relationships otherwise missed, resulting in marked improvement over genetic-only predictions. We go on to establish biological rationale surrounding a number of novel candidate synthetic lethal gene pairs with demonstrated dependencies in published cancer cell line shRNA screens.

Conclusions

This work introduces a multi-omic approach to define gene loss-of-function, and enrich for candidate synthetic lethal gene pairs in cell lines testable through functional screens. In doing so, we offer an additional resource to generate new cancer drug target and combination hypotheses. Algorithms discussed are freely available in the BiSEp CRAN package at http://cran.r-project.org/web/packages/BiSEp/index.html.

Electronic supplementary material

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

Bayesian nonparametric estimation of targeted agent effects on biomarker change to predict clinical outcome

The effect of a targeted agent on a cancer patient's clinical outcome putatively is mediated through the agent's effect on one or more early biological events. This is motivated by pre-clinical experiments with cells or animals that identify such events, represented by binary or quantitative biomarkers. When evaluating targeted agents in humans, central questions are whether the distribution of a targeted biomarker changes following treatment, the nature and magnitude of this change, and whether it is associated with clinical outcome. Major difficulties in estimating these effects are that a biomarker's distribution may be complex, vary substantially between patients, and have complicated relationships with clinical outcomes. We present a probabilistically coherent framework for modeling and estimation in this setting, including a hierarchical Bayesian nonparametric mixture model for biomarkers that we use to define a functional profile of pre-versus-post-treatment biomarker distribution change. The functional is similar to the receiver operating characteristic used in diagnostic testing. The hierarchical model yields clusters of individual patient biomarker profile functionals, and we use the profile as a covariate in a regression model for clinical outcome. The methodology is illustrated by analysis of a dataset from a clinical trial in prostate cancer using imatinib to target platelet-derived growth factor, with the clinical aim to improve progression-free survival time.

Gene co-expression modules as clinically relevant hallmarks of breast cancer diversity

Co-expression modules are groups of genes with highly correlated expression patterns. In cancer, differences in module activity potentially represent the heterogeneity of phenotypes important in carcinogenesis, progression, or treatment response. To find gene expression modules active in breast cancer subpopulations, we assembled 72 breast cancer-related gene expression datasets containing ∼5,700 samples altogether. Per dataset, we identified genes with bimodal expression and used mixture-model clustering to ultimately define 11 modules of genes that are consistently co-regulated across multiple datasets. Functionally, these modules reflected estrogen signaling, development/differentiation, immune signaling, histone modification, ERBB2 signaling, the extracellular matrix (ECM) and stroma, and cell proliferation. The Tcell/Bcell immune modules appeared tumor-extrinsic, with coherent expression in tumors but not cell lines; whereas most other modules, interferon and ECM included, appeared intrinsic. Only four of the eleven modules were represented in the PAM50 intrinsic subtype classifier and other well-established prognostic signatures; although the immune modules were highly correlated to previously published immune signatures. As expected, the proliferation module was highly associated with decreased recurrence-free survival (RFS). Interestingly, the immune modules appeared associated with RFS even after adjustment for receptor subtype and proliferation; and in a multivariate analysis, the combination of Tcell/Bcell immune module down-regulation and proliferation module upregulation strongly associated with decreased RFS. Immune modules are unusual in that their upregulation is associated with a good prognosis without chemotherapy and a good response to chemotherapy, suggesting the paradox of high immune patients who respond to chemotherapy but would do well without it. Other findings concern the ECM/stromal modules, which despite common themes were associated with different sites of metastasis, possibly relating to the “seed and soil” hypothesis of cancer dissemination. Overall, co-expression modules provide a high-level functional view of breast cancer that complements the “cancer hallmarks” and may form the basis for improved predictors and treatments.

Integrative DNA methylation and gene expression analysis in high-grade soft tissue sarcomas

Background

High-grade soft tissue sarcomas are a heterogeneous, complex group of aggressive malignant tumors showing mesenchymal differentiation. Recently, soft tissue sarcomas have increasingly been classified on the basis of underlying genetic alterations; however, the role of aberrant DNA methylation in these tumors is not well understood and, consequently, the usefulness of methylation-based classification is unclear.

Results

We used the Infinium HumanMethylation27 platform to profile DNA methylation in 80 primary, untreated high-grade soft tissue sarcomas, representing eight relevant subtypes, two non-neoplastic fat samples and 14 representative sarcoma cell lines. The primary samples were partitioned into seven stable clusters. A classification algorithm identified 216 CpG sites, mapping to 246 genes, showing different degrees of DNA methylation between these seven groups. The differences between the clusters were best represented by a set of eight CpG sites located in the genes SPEG, NNAT, FBLN2, PYROXD2, ZNF217, COL14A1, DMRT2 and CDKN2A. By integrating DNA methylation and mRNA expression data, we identified 27 genes showing negative and three genes showing positive correlation. Compared with non-neoplastic fat, NNAT showed DNA hypomethylation and inverse gene expression in myxoid liposarcomas, and DNA hypermethylation and inverse gene expression in dedifferentiated and pleomorphic liposarcomas. Recovery of NNAT in a hypermethylated myxoid liposarcoma cell line decreased cell migration and viability.

Conclusions

Our analysis represents the first comprehensive integration of DNA methylation and transcriptional data in primary high-grade soft tissue sarcomas. We propose novel biomarkers and genes relevant for pathogenesis, including NNAT as a potential tumor suppressor in myxoid liposarcomas.

A variational Bayes beta mixture model for feature selection in DNA methylation studies

An increasing number of studies are using beadarrays to measure DNA methylation on a genome-wide basis. The purpose is to identify novel biomarkers in a wide range of complex genetic diseases including cancer. A common difficulty encountered in these studies is distinguishing true biomarkers from false positives. While statistical methods aimed at improving the feature selection step have been developed for gene expression, relatively few methods have been adapted to DNA methylation data, which is naturally beta-distributed. Here we explore and propose an innovative application of a recently developed variational Bayesian beta-mixture model (VBBMM) to the feature selection problem in the context of DNA methylation data generated from a highly popular beadarray technology. We demonstrate that VBBMM offers significant improvements in inference and feature selection in this type of data compared to an Expectation-Maximization (EM) algorithm, at a significantly reduced computational cost. We further demonstrate the added value of VBBMM as a feature selection and prioritization step in the context of identifying prognostic markers in breast cancer. A variational Bayesian approach to feature selection of DNA methylation profiles should thus be of value to any study undergoing large-scale DNA methylation profiling in search of novel biomarkers.

Mendelian and non-Mendelian regulation of gene expression in maize

Transcriptome variation plays an important role in affecting the phenotype of an organism. However, an understanding of the underlying mechanisms regulating transcriptome variation in segregating populations is still largely unknown. We sought to assess and map variation in transcript abundance in maize shoot apices in the intermated B73×Mo17 recombinant inbred line population. RNA–based sequencing (RNA–seq) allowed for the detection and quantification of the transcript abundance derived from 28,603 genes. For a majority of these genes, the population mean, coefficient of variation, and segregation patterns could be predicted by the parental expression levels. Expression quantitative trait loci (eQTL) mapping identified 30,774 eQTL including 96 trans-eQTL “hotspots,” each of which regulates the expression of a large number of genes. Interestingly, genes regulated by a trans-eQTL hotspot tend to be enriched for a specific function or act in the same genetic pathway. Also, genomic structural variation appeared to contribute to cis-regulation of gene expression. Besides genes showing Mendelian inheritance in the RIL population, we also found genes whose expression level and variation in the progeny could not be predicted based on parental difference, indicating that non-Mendelian factors also contribute to expression variation. Specifically, we found 145 genes that show patterns of expression reminiscent of paramutation such that all the progeny had expression levels similar to one of the two parents. Furthermore, we identified another 210 genes that exhibited unexpected patterns of transcript presence/absence. Many of these genes are likely to be gene fragments resulting from transposition, and the presence/absence of their transcripts could influence expression levels of their ancestral syntenic genes. Overall, our results contribute to the identification of novel expression patterns and broaden the understanding of transcriptional variation in plants.

Phenotypes are determined by the expression of genes, the environment, and the interaction of gene expression and the environment. However, a complete understanding of the inheritance of and genome-wide regulation of gene expression is lacking. One approach, called expression quantitative trait locus (eQTL) mapping provides the opportunity to examine the genome-wide inheritance and regulation of gene expression. In this paper, we conducted high-throughput sequencing of gene transcripts to examine gene expression in the shoot apex of a maize biparental mapping population. We quantified expression levels from 28,603 genes in the population and showed that the vast majority of genes exhibited the expected pattern of Mendelian inheritance. We genetically mapped the expression patterns and identified genomic regions associated with gene expression. Notably, we detected gene expression patterns that exhibited non-Mendelian inheritance. These included 145 genes that exhibited expression patterns in the progeny that were similar to only one of the parents and 210 genes with unexpected presence/absence expression patterns. The findings of non-Mendelian inheritance underscore the complexity of gene expression and provide a framework for understanding these complexities.

BiC: a web server for calculating bimodality of coexpression between gene and protein networks

Bimodal patterns of expression have recently been shown to be useful not only in prioritizing genes that distinguish phenotypes, but also in prioritizing network models that correlate with proteomic evidence. In particular, subgroups of strongly coexpressed gene pairs result in an increased variance of the correlation distribution. This variance, a measure of association between sets of genes (or proteins), can be summarized as the bimodality of coexpression (BiC). We developed an online tool to calculate the BiC for user-defined gene lists and associated mRNA expression data. BiC is a comprehensive application that provides researchers with the ability to analyze both publicly available and user-collected array data.

AVAILABILITY

The freely available web service and the documentation can be accessed at http://gurkan.case.edu/software.

CONTACT

gurkan@case.edu.

Bimodal distribution of RNA expression levels in human skeletal muscle tissue

Background

Many human diseases and phenotypes are related to RNA expression, levels of which are influenced by a wide spectrum of genetic and exposure-related factors. In a large genome-wide study of muscle tissue expression, we found that some genes exhibited a bimodal distribution of RNA expression, in contrast to what is usually assumed in studies of a single healthy tissue. As bimodality has classically been considered a hallmark of genetic control, we assessed the genome-wide prevalence, cause, and association of this phenomenon with diabetes-related phenotypes in skeletal muscle tissue from 225 healthy Pima Indians using exon array expression chips.

Results

Two independent batches of microarrays were used for bimodal assessment and comparison. Of the 17,881 genes analyzed, eight (GSTM1, HLA-DRB1, ERAP2, HLA-DRB5, MAOA, ACTN3, NR4A2, and THNSL2) were found to have bimodal expression replicated in the separate batch groups, while 24 other genes had evidence of bimodality in only one group. Some bimodally expressed genes had modest associations with pre-diabetic phenotypes, of note ACTN3 with insulin resistance. Most of the other bimodal genes have been reported to be involved with various other diseases and characteristics. Association of expression with cis genetic variation in a subset of 149 individuals found all but one of the confirmed bimodal genes and nearly half of all potential ones to be highly significant expression quantitative trait loci (eQTL). The rare prevalence of these bimodally expressed genes found after controlling for batch effects was much lower than the prevalence reported in other studies. Additional validation in data from separate muscle expression studies confirmed the low prevalence of bimodality we observed.

Conclusions

We conclude that the prevalence of bimodal gene expression is quite rare in healthy muscle tissue (<0.2%), and is much lower than limited reports from other studies. The major cause of these clearly bimodal expression patterns in homogeneous tissue appears to be cis-polymorphisms, indicating that such bimodal genes are, for the most part, eQTL. The high frequency of disease associations reported with these genes gives hope that this unique feature may identify or actually be an underlying factor responsible for disease development.

Minimising immunohistochemical false negative ER classification using a complementary 23 gene expression signature of ER status

Background

Expression of the oestrogen receptor (ER) in breast cancer predicts benefit from endocrine therapy. Minimising the frequency of false negative ER status classification is essential to identify all patients with ER positive breast cancers who should be offered endocrine therapies in order to improve clinical outcome. In routine oncological practice ER status is determined by semi-quantitative methods such as immunohistochemistry (IHC) or other immunoassays in which the ER expression level is compared to an empirical threshold[1], [2]. The clinical relevance of gene expression-based ER subtypes as compared to IHC-based determination has not been systematically evaluated. Here we attempt to reduce the frequency of false negative ER status classification using two gene expression approaches and compare these methods to IHC based ER status in terms of predictive and prognostic concordance with clinical outcome.

Methodology/Principal Findings

Firstly, ER status was discriminated by fitting the bimodal expression of ESR1 to a mixed Gaussian model. The discriminative power of ESR1 suggested bimodal expression as an efficient way to stratify breast cancer; therefore we identified a set of genes whose expression was both strongly bimodal, mimicking ESR expression status, and highly expressed in breast epithelial cell lines, to derive a 23-gene ER expression signature-based classifier. We assessed our classifiers in seven published breast cancer cohorts by comparing the gene expression-based ER status to IHC-based ER status as a predictor of clinical outcome in both untreated and tamoxifen treated cohorts. In untreated breast cancer cohorts, the 23 gene signature-based ER status provided significantly improved prognostic power compared to IHC-based ER status (P = 0.006). In tamoxifen-treated cohorts, the 23 gene ER expression signature predicted clinical outcome (HR = 2.20, P = 0.00035). These complementary ER signature-based strategies estimated that between 15.1% and 21.8% patients of IHC-based negative ER status would be classified with ER positive breast cancer.

Conclusion/Significance

Expression-based ER status classification may complement IHC to minimise false negative ER status classification and optimise patient stratification for endocrine therapies.

A global analysis of QTLs for expression variations in rice shoots at the early seedling stage

Analyses of quantitative trait loci (QTLs) for expression levels (eQTLs) of genes reveal a genetic relationship between expression variation and the regulator, thus unlocking information for identifying the regulatory network. Oligo-nucleotide expression microarrays hybridized with RNA can simultaneously provide data for molecular markers and transcript abundance. In this study, we used an Affymetrix GeneChip Rice Genome Array to analyze eQTLs in rice shoots at 72 h after germination from 110 recombinant inbred lines (RILs) derived from a cross between Zhenshan 97 and Minghui 63. In total, 1632 single-feature polymorphisms (SFPs) plus 23 PCR markers were identified and placed into 601 recombinant bins, spanning 1459 cM in length, which were used as markers to genotype the RILs. We obtained 16,372 expression traits (e-traits) each with at least one eQTL, resulting in 26,051 eQTLs in total, including both cis- and trans-eQTLs. We also identified 171 eQTL hot spots in the rice genome, each of which controls transcript variations of many e-traits. Gene ontology analysis revealed an enrichment of certain functional categories of genes in some of the eQTL hot spots. In particular, eQTLs for e-traits involving the DNA metabolic process was significantly enriched in several eQTL hot spots on chromosomes 3, 5 and 10. Several e-traits co-localizing with cis-eQTLs showed significant correlations with hundreds of e-traits, indicating possible co-regulation. We also detected correlations between QTLs for shoot dry weight and eQTLs, revealing possible candidate genes for the trait. These results provided clues for the identification and characterization of the regulatory network in the whole genome at the transcriptional level.

Rough set soft computing cancer classification and network: one stone, two birds

Gene expression profiling provides tremendous information to help unravel the complexity of cancer. The selection of the most informative genes from huge noise for cancer classification has taken centre stage, along with predicting the function of such identified genes and the construction of direct gene regulatory networks at different system levels with a tuneable parameter. A new study by Wang and Gotoh described a novel Variable Precision Rough Sets-rooted robust soft computing method to successfully address these problems and has yielded some new insights. The significance of this progress and its perspectives will be discussed in this article.