Spotlight on differentially expressed genes in urinary bladder cancer

Introduction

We previously identified common differentially expressed (DE) genes in bladder cancer (BC). In the present study we analyzed in depth, the expression of several groups of these DE genes.

Materials and Methods

Samples from 30 human BCs and their adjacent normal tissues were analyzed by whole genome cDNA microarrays, qRT-PCR and Western blotting. Our attention was focused on cell-cycle control and DNA damage repair genes, genes related to apoptosis, signal transduction, angiogenesis, as well as cellular proliferation, invasion and metastasis. Four publicly available GEO Datasets were further analyzed, and the expression data of the genes of interest (GOIs) were compared to those of the present study. The relationship among the GOI was also investigated. GO and KEGG molecular pathway analysis was performed to identify possible enrichment of genes with specific biological themes.

Results

Unsupervised cluster analysis of DNA microarray data revealed a clear distinction in BC vs. control samples and low vs. high grade tumors. Genes with at least 2-fold differential expression in BC vs. controls, as well as in non-muscle invasive vs. muscle invasive tumors and in low vs. high grade tumors, were identified and ranked. Specific attention was paid to the changes in osteopontin (OPN, SPP1) expression, due to its multiple biological functions. Similarly, genes exhibiting equal or low expression in BC vs. the controls were scored. Significant pair-wise correlations in gene expression were scored. GO analysis revealed the multi-facet character of the GOIs, since they participate in a variety of mechanisms, including cell proliferation, cell death, metabolism, cell shape, and cytoskeletal re-organization. KEGG analysis revealed that the most significant pathway was that of Bladder Cancer (p = 1.5×10⁻³¹).

Conclusions

The present work adds to the current knowledge on molecular signature identification of BC. Such works should progress in order to gain more insight into disease molecular mechanisms.     

Clustered gene expression changes flank targeted gene loci in knockout mice

Background

Gene expression profiling using microarrays is a powerful technology widely used to study regulatory networks. Profiling of mRNA levels in mutant organisms has the potential to identify genes regulated by the mutated protein.

Methodology/Principle Findings

Using tissues from multiple lines of knockout mice we have examined genome-wide changes in gene expression. We report that a significant proportion of changed genes were found near the targeted gene.

Conclusions/Significance

The apparent clustering of these genes was explained by the presence of flanking DNA from the parental ES cell. We provide recommendations for the analysis and reporting of microarray data from knockout mice     

Methylcap-seq reveals novel DNA methylation markers for the diagnosis and recurrence prediction of bladder cancer in a Chinese population

Purpose

There is a need to supplement or supplant the conventional diagnostic tools, namely, cystoscopy and B-type ultrasound, for bladder cancer (BC). We aimed to identify novel DNA methylation markers for BC through genome-wide profiling of BC cell lines and subsequent methylation-specific PCR (MSP) screening of clinical urine samples.

Experimental Design

The methyl-DNA binding domain (MBD) capture technique, methylCap/seq, was performed to screen for specific hypermethylated CpG islands in two BC cell lines (5637 and T24). The top one hundred hypermethylated targets were sequentially screened by MSP in urine samples to gradually narrow the target number and optimize the composition of the diagnostic panel. The diagnostic performance of the obtained panel was evaluated in different clinical scenarios.

Results

A total of 1,627 hypermethylated promoter targets in the BC cell lines was identified by Illumina sequencing. The top 104 hypermethylated targets were reduced to eight genes (VAX1, KCNV1, ECEL1, TMEM26, TAL1, PROX1, SLC6A20, and LMX1A) after the urine DNA screening in a small sample size of 8 normal control and 18 BC subjects. Validation in an independent sample of 212 BC patients enabled the optimization of five methylation targets, including VAX1, KCNV1, TAL1, PPOX1, and CFTR, which was obtained in our previous study, for BC diagnosis with a sensitivity and specificity of 88.68% and 87.25%, respectively. In addition, the methylation of VAX1 and LMX1A was found to be associated with BC recurrence.

Conclusions

We identified a promising diagnostic marker panel for early non-invasive detection and subsequent BC surveillance.     

Genomic Responses during Acute Human Anaphylaxis Are Characterized by Upregulation of Innate Inflammatory Gene Networks

Background

Systemic spread of immune activation and mediator release is required for the development of anaphylaxis in humans. We hypothesized that peripheral blood leukocyte (PBL) activation plays a key role.

Objective

To characterize PBL genomic responses during acute anaphylaxis.

Methods

PBL samples were collected at three timepoints from six patients presenting to the Emergency Department (ED) with acute anaphylaxis and six healthy controls. Gene expression patterns were profiled on microarrays, differentially expressed genes were identified, and network analysis was employed to explore underlying mechanisms.

Results

Patients presented with moderately severe anaphylaxis after oral aspirin (2), peanut (2), bee sting (1) and unknown cause (1). Two genes were differentially expressed in patients compared to controls at ED arrival, 67 genes at 1 hour post-arrival and 2,801 genes at 3 hours post-arrival. Network analysis demonstrated that three inflammatory modules were upregulated during anaphylaxis. Notably, these modules contained multiple hub genes, which are known to play a central role in the regulation of innate inflammatory responses. Bioinformatics analyses showed that the data were enriched for LPS-like and TNF activation signatures.

Conclusion

PBL genomic responses during human anaphylaxis are characterized by dynamic expression of innate inflammatory modules. Upregulation of these modules was observed in patients with different reaction triggers. Our findings indicate a role for innate immune pathways in the pathogenesis of human anaphylaxis, and the hub genes identified in this study represent logical candidates for follow-up studies.     

Coordinated Expression Domains in Mammalian Genomes

Background

Gene order in eukaryotic genomes is not random. Genes showing similar expression (coexpression) patterns are often clustered along the genome. The goal of this study is to characterize coexpression clustering in mammalian genomes and to investigate the underlying mechanisms.

Methodology/Principal Findings

We detect clustering of coexpressed genes across multiple scales, from neighboring genes to chromosomal domains that span tens of megabases and, in some cases, entire chromosomes. Coexpression domains may be positively or negatively correlated with other domains, within and between chromosomes. We find that long-range expression domains are associated with gene density, which in turn is related to physical organization of the chromosomes within the nucleus. We show that gene expression changes between healthy and diseased tissue samples occur in a gene density-dependent manner.

Conclusions/Significance

We demonstrate that coexpression domains exist across multiple scales. We identify potential mechanisms for short-range as well as long-range coexpression domains. We provide evidence that the three-dimensional architecture of the chromosomes may underlie long-range coexpression domains. Chromosome territory reorganization may play a role in common human diseases such as Alzheimer''s disease and psoriasis.     

Effects of Warm Ischemic Time on Gene Expression Profiling in Colorectal Cancer Tissues and Normal Mucosa

Background

Genome-wide gene expression analyses of tumors are a powerful tool to identify gene signatures associated with biologically and clinically relevant characteristics and for several tumor types are under clinical validation by prospective trials. However, handling and processing of clinical specimens may significantly affect the molecular data obtained from their analysis. We studied the effects of tissue handling time on gene expression in human normal and tumor colon tissues undergoing routine surgical procedures.

Methods

RNA extracted from specimens of 15 patients at four time points (for a total of 180 samples) after surgery was analyzed for gene expression on high-density oligonucleotide microarrays. A mixed-effects model was used to identify probes with different expression means across the four different time points. The p-values of the model were adjusted with the Bonferroni method.

Results

Thirty-two probe sets associated with tissue handling time in the tumor specimens, and thirty-one in the normal tissues, were identified. Most genes exhibited moderate changes in expression over the time points analyzed; however four of them were oncogenes, and two confirmed the effect of tissue handling by independent validation.

Conclusions

Our results suggest that a critical time point for tissue handling in colon seems to be 60 minutes at room temperature. Although the number of time-dependent genes we identified was low, the three genes that already showed changes at this time point in tumor samples were all oncogenes, hence recommending standardization of tissue-handling protocols and effort to reduce the time from specimen removal to snap freezing accounting for warm ischemia in this tumor type.     

Deep Sequencing of Target Linkage Assay-Identified Regions in Familial Breast Cancer: Methods,Analysis Pipeline and Troubleshooting

Background

The classical candidate-gene approach has failed to identify novel breast cancer susceptibility genes. Nowadays, massive parallel sequencing technology allows the development of studies unaffordable a few years ago. However, analysis protocols are not yet sufficiently developed to extract all information from the huge amount of data obtained.

Methodology/Principal Findings

In this study, we performed high throughput sequencing in two regions located on chromosomes 3 and 6, recently identified by linkage studies by our group as candidate regions for harbouring breast cancer susceptibility genes. In order to enrich for the coding regions of all described genes located in both candidate regions, a hybrid-selection method on tiling microarrays was performed.

Conclusions/Significance

We developed an analysis pipeline based on SOAP aligner to identify candidate variants with a high real positive confirmation rate (0.89), with which we identified eight variants considered candidates for functional studies. The results suggest that the present strategy might be a valid second step for identifying high penetrance genes.     

Polymorphic Repeat Length in the AIB1 Gene and Breast Cancer Risk in BRCA1 and BRCA2 Mutation Carriers: A Meta-Analysis of Observational Studies

Objectives

We carried out a meta-analysis focusing on the relationship between length of AIB1 gene poly-Q repeat domain as a modifier of breast cancer (BC) susceptibility in patients with BRCA1 and BRCA2 mutation carriers.

Data sources

We searched MEDLINE and EMBASE for all medical literature published until February, 2012.

Study Eligibility criteria

Studies were included in the meta-analysis if they met all the predetermined criteria, such as: (a) case-control or cohort studies; (b) the primary outcome was clearly defined as BC; (c) the exposure of interest measured was AIB1 polyglutamine repeat length genotype; (d) provided relative risk (RR) or odds ratio (OR) estimates and their 95% confidence intervals (CIs).

Synthesis methods

Two of the authors independently evaluated the quality of the studies included and extracted the data. Meta-analyses were performed for case-control and cohort studies separately. Heterogeneity was examined and the publication bias was assessed with a funnel plot for asymmetry.

Result

7 studies met our predetermined inclusion criteria and were included in the meta-analysis. Overall quality ratings of the studies varied from 0.36 to 0.77, with a median of 0.5. The overall RR estimates of 29/29 poly-Q repeats on risk of BC in BRCA1/2, BRCA1, and BRCA2, were always greater than 1.00; however, this effect was not statistically significant. In the meta-analysis of studies reporting the effect of 28/28 poly-Q repeats on risk of BC in BRCA1/2, BRCA1, and BRCA2, the overall RR decreased below 1.00; however, this effect was not statistically significant. Similar estimates were shown for at least 1 allele of ≤26 repeats.

Conclusions

Genotypes of AIB1 polyglutamine polymorphism analyzed do not appear to be associated to a modified risk of BC development in BRCA1 and BRCA2 mutation carriers. Future research on length of poly-Q repeat domain and BC susceptibility should be discouraged and more promising potential sources of penetrance variation among BRCA1 and BRCA2 mutation carriers should be investigated.     

GWAS Meets Microarray: Are the Results of Genome-Wide Association Studies and Gene-Expression Profiling Consistent? Prostate Cancer as an Example

Background

Genome-wide association studies (GWASs) and global profiling of gene expression (microarrays) are two major technological breakthroughs that allow hypothesis-free identification of candidate genes associated with tumorigenesis. It is not obvious whether there is a consistency between the candidate genes identified by GWAS (GWAS genes) and those identified by profiling gene expression (microarray genes).

Methodology/Principal Findings

We used the Cancer Genetic Markers Susceptibility database to retrieve single nucleotide polymorphisms from candidate genes for prostate cancer. In addition, we conducted a large meta-analysis of gene expression data in normal prostate and prostate tumor tissue. We identified 13,905 genes that were interrogated by both GWASs and microarrays. On the basis of P values from GWASs, we selected 1,649 most significantly associated genes for functional annotation by the Database for Annotation, Visualization and Integrated Discovery. We also conducted functional annotation analysis using same number of the top genes identified in the meta-analysis of the gene expression data. We found that genes involved in cell adhesion were overrepresented among both the GWAS and microarray genes.

Conclusions/Significance

We conclude that the results of these analyses suggest that combining GWAS and microarray data would be a more effective approach than analyzing individual datasets and can help to refine the identification of candidate genes and functions associated with tumor development.     

Unacylated Ghrelin Rapidly Modulates Lipogenic and Insulin Signaling Pathway Gene Expression in Metabolically Active Tissues of GHSR Deleted Mice

Background

There is increasing evidence that unacylated ghrelin (UAG) improves insulin sensitivity and glucose homeostasis; however, the mechanism for this activity is not fully understood since a UAG receptor has not been discovered.

Methodology/Principal Findings

To assess potential mechanisms of UAG action in vivo, we examined rapid effects of UAG on genome-wide expression patterns in fat, muscle and liver of growth hormone secretagogue receptor (GHSR)-ablated mice using microarrays. Expression data were analyzed using Ingenuity Pathways Analysis and Gene Set Enrichment Analysis. Regulation of subsets of these genes was verified by quantitative PCR in an independent experiment. UAG acutely regulated clusters of genes involved in glucose and lipid metabolism in all three tissues, consistent with enhancement of insulin sensitivity.

Conclusions/Significance

Fat, muscle and liver are central to the control of lipid and glucose homeostasis. UAG rapidly modulates the expression of metabolically important genes in these tissues in GHSR-deleted mice indicating a direct, GHSR-independent, action of UAG to improve insulin sensitivity and metabolic profile.     

Comparison of Disease Activity in SPMS and PPMS in the Context of Multicenter Clinical Trials

Background

Retrospective single center natural history studies have shown that times to reach disability milestones and ages at which they are reached are similar in primary (PPMS) and secondary (SPMS) progressive multiple sclerosis suggesting that they may be phenotypic variations of the same disease.

Objective

Here we compared longitudinal disease activity in SPMS and PPMS in the context of international multicenter clinical trials.

Methods

We analyzed all objective outcome measures that were systematically collected over 2 years for all subjects randomized to placebo arms in one SPMS and one PPMS clinical trial over the last decade. Conventional and exploratory definitions of clinical disease activity were used. Disease activity was analyzed in 3 different categories intermittent activity, progression, and improvement. Conventional MRI measures and one patient reported outcome measure of quality of life were included when available for comparison. Heat maps were drawn for all results followed by hierarchical clustering.

Results

There were 101 outcome variables from 206 SPMS subjects and 79 outcome variables from 135 PPMS subjects. The comparison revealed that SPMS and PPMS subjects exhibited similar disease activity over 2 years in all but two of the variables in common worsening in the EDSS sensory system was more common in PPMS while worsening on the 9 hole PEG was more common in SPMS. Intermittent activity was the most common pattern of disease activity in SPMS and PPMS. Clinical worsening and improvement occurred at similar frequency in both.

Conclusion

Longitudinal disease activity was nearly identical in SPMS and PPMS subjects in the context of the two multicenter international clinical trials we examined.