Roadmap to DILI research in Europe. A proposal from COST action ProEuroDILINet

In the current article the aims for a constructive way forward in Drug-Induced Liver Injury (DILI) are to highlight the most important priorities in research and clinical science, therefore supporting a more informed, focused, and better funded future for European DILI research. This Roadmap aims to identify key challenges, define a shared vision across all stakeholders for the opportunities to overcome these challenges and propose a high-quality research program to achieve progress on the prediction, prevention, diagnosis and management of this condition and impact on healthcare practice in the field of DILI. This will involve 1. Creation of a database encompassing optimised case report form for prospectively identified DILI cases with well-characterised controls with competing diagnoses, biological samples, and imaging data; 2. Establishing of preclinical models to improve the assessment and prediction of hepatotoxicity in humans to guide future drug safety testing; 3. Emphasis on implementation science and 4. Enhanced collaboration between drug-developers, clinicians and regulatory scientists. This proposed operational framework will advance DILI research and may bring together basic, applied, translational and clinical research in DILI.

CAP-RNAseq: an integrated pipeline for functional annotation and prioritization of co-expression clusters

Cluster analysis is one of the most widely used exploratory methods for visualization and grouping of gene expression patterns across multiple samples or treatment groups. Although several existing online tools can annotate clusters with functional terms, there is no all-in-one webserver to effectively prioritize genes/clusters using gene essentiality as well as congruency of mRNA-protein expression. Hence, we developed CAP-RNAseq that makes possible (1) upload and clustering of bulk RNA-seq data followed by identification, annotation and network visualization of all or selected clusters; and (2) prioritization using DepMap gene essentiality and/or dependency scores as well as the degree of correlation between mRNA and protein levels of genes within an expression cluster. In addition, CAP-RNAseq has an integrated primer design tool for the prioritized genes. Herein, we showed using comparisons with the existing tools and multiple case studies that CAP-RNAseq can uniquely aid in the discovery of co-expression clusters enriched with essential genes and prioritization of novel biomarker genes that exhibit high correlations between their mRNA and protein expression levels. CAP-RNAseq is applicable to RNA-seq data from different contexts including cancer and available at http://konulabapps.bilkent.edu.tr:3838/CAPRNAseq/ and the docker image is downloadable from https://hub.docker.com/r/konulab/caprnaseq.

Long-Term Acetylcholinesterase Depletion Alters the Levels of Key Synaptic Proteins while Maintaining Neuronal Markers in the Aging Zebrafish (Danio rerio) Brain

INTRODUCTION

Interventions targeting cholinergic neurotransmission like acetylcholinesterase (AChE) inhibition distinguish potential mechanisms to delay age-related impairments and attenuate deficits related to neurodegenerative diseases. However, the chronic effects of these interventions are not well described.

METHODS

In the current study, global levels of cholinergic, cellular, synaptic, and inflammation-mediating proteins were assessed within the context of aging and chronic reduction of AChE activity. Long-term depletion of AChE activity was induced by using a mutant zebrafish line, and they were compared with the wildtype group at young and old ages.

RESULTS

Results demonstrated that AChE activity was lower in both young and old mutants, and this decrease coincided with a reduction in ACh content. Additionally, an overall age-related reduction in AChE activity and the AChE/ACh ratio was observed, and this decline was more prominent in wildtype groups. The levels of an immature neuronal marker were upregulated in mutants, while a glial marker showed an overall reduction. Mutants had preserved levels of inhibitory and presynaptic elements with aging, whereas glutamate receptor subunit levels declined.

CONCLUSION

Long-term AChE activity depletion induces synaptic and cellular alterations. These data provide further insights into molecular targets and adaptive responses following the long-term reduction of AChE activity that was also targeted pharmacologically to treat neurodegenerative diseases in human subjects.

Zebrafish Xenotransplantation Models for Studying Gene Function and Drug Treatment in Hepatocellular Carcinoma

INTRODUCTION

Zebrafish is a promising model organism for human disease including hepatocellular cancer (HCC). Recently, zebrafish has emerged also as a host for xenograft studies of liver cancer cell lines and patient derived tumors of HCC. Zebrafish embryos enable drug screening and gene function studies of xenografted cells via ease of microinjection and visualization of tumor growth and metastasis.

OBJECTIVES

In this review, we aimed to overview zebrafish HCC and liver cancer xenotransplantation studies focusing on 'gene functional analysis' and 'drug/chemical screening'.

METHODS

Herein, a comprehensive literature search was performed for liver and HCC xenografts in zebrafish on PubMed using different key words and filters for molecular modifications or drug exposure.

RESULTS

Our literature search revealed around 250 studies which were filtered and summarized in a table (Table 1) revealing comprehensive collection of experimental and technical details on microinjection, injected cell lines, molecular modifications of injected cells, types and doses of drug treatments as well as biological assessments.

CONCLUSION

This review provides a platform for HCC and liver xenografts and highlights studies performed to understand gene functionality and drug efficacy in vivo in zebrafish.

Global miRNA expression of bone marrow mesenchymal stem/stromal cells derived from Fanconi anemia patients

Fanconi anemia (FA) is a rare genetic disorder characterized by genomic instability, developmental defects, and bone marrow (BM) failure. Hematopoietic stem cells (HSCs) in BM interact with the mesenchymal stem/stromal cells (MSCs); and this partly sustains the tissue homeostasis. MicroRNAs (miRNAs) can play a critical role during these interactions possibly via paracrine mechanisms. This is the first study addressing the miRNA profile of FA BM-MSCs obtained before and after BM transplantation (preBMT and postBMT, respectively). Non-coding RNA expression profiling and quality control analyses were performed in Donors (n = 13), FA preBMT (n = 11), and FA postBMT (n = 6) BM-MSCs using GeneChip miRNA 2.0 Array. Six Donor-FA preBMT pairs were used to identify a differentially expressed miRNA expression signature containing 50 miRNAs, which exhibited a strong correlation with the signature obtained from unpaired samples. Five miRNAs (hsa-miR-146a-5p, hsa-miR-148b-3p, hsa-miR-187-3p, hsa-miR-196b-5p, and hsa-miR-25-3p) significantly downregulated in both the paired and unpaired analyses were used to generate the BM-MSCs' miRNA-BM mononuclear mRNA networks upon integration of a public dataset (GSE16334; studying Donor versus FA samples). Functionally enriched KEGG pathways included cellular senescence, miRNAs, and pathways in cancer. Here, we showed that hsa-miR-146a-5p and hsa-miR-874-3p were rescued upon BMT (n = 3 triplets). The decrease in miR-146a-5p was also validated using RT-qPCR and emerged as a strong candidate as a modulator of BM mRNAs in FA patients.

Discovery of Cancer-Specific and Independent Prognostic Gene Subsets of the Slit-Robo Family Using TCGA-PANCAN Datasets

The Slit-Robo family of axon guidance molecules works in concert, playing important roles in organ development and cancer. Expressions of individual Slit-Robo genes have been used in calculating univariable hazard ratios (HR_uni) for predicting cancer prognosis in the literature. However, Slit-Robo members do not act independently; hence, hazard ratios from multivariable Cox regression (HR_multi) on the whole gene set can further lead to identification of cancer-specific, novel, and independent prognostic gene pairs or modules. Herein, we obtained mRNA expressions of the Slit-Robo family consisting of four Robos (ROBO1/2/3/4) and three Slits (SLIT1/2/3), along with four types of survival outcome across cancers found in the Cancer Genome Atlas (TCGA). We used cluster heat maps to visualize closely associated pairs/modules of prognostic genes across 33 different cancers. We found a smaller number of significant genes in HR_multi than in HR_uni, suggesting that the former analysis was less redundant. High ROBO4 expression emerged as relatively protective within the family, in both types of HR analyses. Multivariable Cox regression, on the other hand, revealed significantly more HR signatures containing Slit-Robo pairs acting in opposing directions than those containing Slit-Slit or Robo-Robo pairs for disease-specific survival. Furthermore, we discovered, through the online app SmulTCan's lasso regression, Slit-Robo gene subsets that significantly differentiated between high- versus low-risk prognosis patient groups, particularly for renal cancers and low-grade glioma. The statistical pipeline reported herein can help test independent and significant pairs/modules within a codependent gene family for cancer prognostication, and thus should also prove useful in personalized/precision medicine research.

Functional analysis of co-expression networks of zebrafish ace2 reveals enrichment of pathways associated with development and disease

Human Angiotensin I Converting Enzyme 2 (ACE2) plays an essential role in blood pressure regulation and SARS-CoV-2 entry. ACE2 has a highly conserved, one-to-one ortholog (ace2) in zebrafish, which is an important model for human diseases. However, the zebrafish ace2 expression profile has not yet been studied during early development, between genders, across different genotypes, or in disease. Moreover, a network-based meta-analysis for the extraction of functionally enriched pathways associated with differential ace2 expression is lacking in the literature. Herein, we first identified significant development-, tissue-, genotype-, and gender-specific modulations in ace2 expression via meta-analysis of zebrafish Affymetrix transcriptomics datasets (n_datasets = 107); and the correlation analysis of ace2 meta-differential expression profile revealed distinct positively and negatively correlated local functionally enriched gene networks. Moreover, we demonstrated that ace2 expression was significantly modulated under different physiological and pathological conditions related to development, tissue, gender, diet, infection, and inflammation using additional RNA-seq datasets. Our findings implicate a novel translational role for zebrafish ace2 in organ differentiation and pathologies observed in the intestines and liver.

Advanced preclinical models for evaluation of drug-induced liver injury - consensus statement by the European Drug-Induced Liver Injury Network [PRO-EURO-DILI-NET]

Drug-induced liver injury (DILI) is a major cause of acute liver failure (ALF) and one of the leading indications for liver transplantation in Western societies. Given the wide use of both prescribed and over the counter drugs, DILI has become a major health issue for which there is a pressing need to find novel and effective therapies. Although significant progress has been made in understanding the molecular mechanisms underlying DILI, our incomplete knowledge of its pathogenesis and inability to predict DILI is largely due to both discordance between human and animal DILI in preclinical drug development and a lack of models that faithfully recapitulate complex pathophysiological features of human DILI. This is exemplified by the hepatotoxicity of acetaminophen (APAP) overdose, a major cause of ALF because of its extensive worldwide use as an analgesic. Despite intensive efforts utilising current animal and in vitro models, the mechanisms involved in the hepatotoxicity of APAP are still not fully understood. In this expert Consensus Statement, which is endorsed by the European Drug-Induced Liver Injury Network, we aim to facilitate and outline clinically impactful discoveries by detailing the requirements for more realistic human-based systems to assess hepatotoxicity and guide future drug safety testing. We present novel insights and discuss major players in APAP pathophysiology, and describe emerging in vitro and in vivo pre-clinical models, as well as advanced imaging and in silico technologies, which may improve prediction of clinical outcomes of DILI.

Increased SGK1 activity potentiates mineralocorticoid/NaCl-induced kidney injury

Serum and glucocorticoid-regulated kinase 1 (SGK1) stimulates aldosterone-dependent renal Na⁺ reabsorption and modulates blood pressure. In addition, genetic ablation or pharmacological inhibition of SGK1 limits the development of kidney inflammation and fibrosis in response to excess mineralocorticoid signaling. In this work, we tested the hypothesis that a systemic increase in SGK1 activity would potentiate mineralocorticoid/salt-induced hypertension and kidney injury. To that end, we used a transgenic mouse model with increased SGK1 activity. Mineralocorticoid/salt-induced hypertension and kidney damage was induced by unilateral nephrectomy and treatment with deoxycorticosterone acetate and NaCl in the drinking water for 6 wk. Our results show that although SGK1 activation did not induce significantly higher blood pressure, it produced a mild increase in glomerular filtration rate, increased albuminuria, and exacerbated glomerular hypertrophy and fibrosis. Transcriptomic analysis showed that extracellular matrix- and immune response-related terms were enriched in the downregulated and upregulated genes, respectively, in transgenic mice. In conclusion, we propose that systemically increased SGK1 activity is a risk factor for the development of mineralocorticoid-dependent kidney injury in the context of low renal mass and independently of blood pressure.NEW & NOTEWORTHY Increased activity of the protein kinase serum and glucocorticoid-regulated kinase 1 may be a risk factor for accelerated renal damage. Serum and glucocorticoid-regulated kinase 1 expression could be a marker for the rapid progression toward chronic kidney disease and a potential therapeutic target to slow down the process.

Design, synthesis, anticancer activity, molecular docking and ADME studies of novel methylsulfonyl indole-benzimidazoles in comparison with ethylsulfonyl counterparts

Specific sidechain modifications on the indole-benzimidazole scaffold play fundamental roles for determining molecule's affinity against ERα and its anti-cancer activity.

ZenoFishDb v1.1: A Database for Xenotransplantation Studies in Zebrafish

Rapidly accumulating literature has proven feasibility of the zebrafish xenograft models in cancer research. Nevertheless, online databases for searching the current zebrafish xenograft literature are in great demand. Herein, we have developed a manually curated database, called ZenoFishDb v1.1 (https://konulab.shinyapps.io/zenofishdb), based on R Shiny platform aiming to provide searchable information on ever increasing collection of zebrafish studies for cancer cell line transplantation and patient-derived xenografts (PDXs). ZenoFishDb v1.1 user interface contains four modules: DataTable, Visualization, PDX Details, and PDX Charts. The DataTable and Visualization pages represent xenograft study details, including injected cell lines, PDX injections, molecular modifications of cell lines, zebrafish strains, as well as technical aspects of the xenotransplantation procedures in table, bar, and/or pie chart formats. The PDX Details module provides comprehensive information on the patient details in table format and can be searched and visualized. Overall, ZenoFishDb v1.1 enables researchers to effectively search, list, and visualize different technical and biological attributes of zebrafish xenotransplantation studies particularly focusing on the new trends that make use of reporters, RNA interference, overexpression, or mutant gene constructs of transplanted cancer cells, stem cells, and PDXs, as well as distinguished host modifications.

Design, synthesis and anticancer/antiestrogenic activities of novel indole-benzimidazoles

Indole-benzimidazoles have recently gained attention due to their antiproliferative and antiestrogenic effects. However, their structural similarities and molecular mechanisms shared with selective estrogen receptor modulators (SERMs) have not yet been investigated. In this study, we synthesized novel ethylsulfonyl indole-benzimidazole derivatives by substituting the first (R₁) and fifth (R₂) positions of benzimidazole and indole groups, respectively. Subsequently, we performed ¹H NMR, ¹³C NMR, and Mass spectral and in silico docking analyses, and anticancer activity screening studies of these novel indole-benzimidazoles. The antiproliferative effects of indole-benzimidazoles were found to be more similar between the estrogen (E2) responsive cell lines MCF-7 and HEPG2 in comparison to the Estrogen Receptor negative (ER-) cell line MDA-MB-231. R₁:p-fluorobenzyl group members were selected as lead compounds for their potent anticancer effects and moderate structural affinity to ER. Microarray expression profiling and gene enrichment analyses (GSEA) of the selected compounds (R₁:p-fluorobenzyl: 48, 49, 50, 51; R₁:3,4-difluorobenzyl: 53) helped determine the similarly modulated cellular signaling pathways among derivatives. Moreover, we identified known compounds that have significantly similar gene signatures to that of 51 via queries performed in LINCS database; and further transcriptomics comparisons were made using public GEO datasets (GSE35428, GSE7765, GSE62673). Our results strongly demonstrate that these novel indole-benzimidazoles can modulate ER target gene expression as well as dioxin-mediated aryl hydrocarbon receptor and amino acid deprivation-mediated integrated stress response signaling in a dose-dependent manner.

Cholinergic Receptor Nicotinic Alpha 5 (CHRNA5) RNAi is associated with cell cycle inhibition, apoptosis, DNA damage response and drug sensitivity in breast cancer

Cholinergic Receptor Nicotinic Alpha 5 (CHRNA5) is an important susceptibility locus for nicotine addiction and lung cancer. Depletion of CHRNA5 has been associated with reduced cell viability, increased apoptosis and alterations in cellular motility in different cancers yet not in breast cancer. Herein we first showed the expression of CHRNA5 was variable and positively correlated with the fraction of total genomic alterations in breast cancer cell lines and tumors indicating its potential role in DNA damage response (DDR). Next, we demonstrated that silencing of CHRNA5 expression in MCF7 breast cancer cell line by RNAi affected expression of genes involved in cytoskeleton, TP53 signaling, DNA synthesis and repair, cell cycle, and apoptosis. The transcription profile of CHRNA5 depleted MCF7 cells showed a significant positive correlation with that of A549 lung cancer cell line while exhibiting a negative association with the CHRNA5 co-expression profile obtained from Cancer Cell Line Encylopedia (CCLE). Moreover, it exhibited high similarities with published MCF7 expression profiles obtained from exposure to TP53 inducer nutlin-3a and topoisomerase inhibitors. We then demonstrated that CHRNA5 siRNA treatment reduced cell viability and DNA synthesis indicating G1 arrest while it significantly increased apoptotic sub-G1 cell population. Accordingly, we observed lower levels of phosphorylated RB (Ser807/811) and an increased BAX/BCL2 ratio in RNAi treated MCF7 cells. We also showed that CHRNA5 RNAi transcriptome correlated negatively with DDR relevant gene expression profile in breast cancer gene expression datasets while the coexposure to topoisomerase inhibitors in the presence of CHRNA5 RNAi enhanced chemosensitivity potentially due to reduced DDR. CHRNA5 RNAi consistently lowered total CHEK1 mRNA and protein levels as well as phosphorylated CHEK1 (Ser345) in MCF7 cells. We also detected a significant positive correlation between the expression levels of CHRNA5 and CHEK1 in CCLE, TCGA and METABRIC breast cancer datasets. Our study suggests CHRNA5 RNAi is associated with cell cycle inhibition, apoptosis as well as reduced DDR and increased drug sensitivity in breast cancer yet future studies are warranted since dose- and cell line-specific differences exist in response to CHRNA5 depletion. Gene expression microarray data can be accessed from GEO database under the accession number GSE89333.

Development of a novel zebrafish xenograft model in ache mutants using liver cancer cell lines

Acetylcholinesterase (AChE), an enzyme responsible for degradation of acetylcholine, has been identified as a prognostic marker in liver cancer. Although in vivo Ache tumorigenicity assays in mouse are present, no established liver cancer xenograft model in zebrafish using an ache mutant background exists. Herein, we developed an embryonic zebrafish xenograft model using epithelial (Hep3B) and mesenchymal (SKHep1) liver cancer cell lines in wild-type and ache sb55 sibling mutant larvae after characterization of cholinesterase expression and activity in cell lines and zebrafish larvae. The comparison of fluorescent signal reflecting tumor size at 3-days post-injection (dpi) revealed an enhanced tumorigenic potential and a reduced migration capacity in cancer cells injected into homozygous ache sb55 mutants when compared with the wild-type. Increased tumor load was confirmed using an ALU based tumor DNA quantification method modified for use in genotyped xenotransplanted zebrafish embryos. Confocal microscopy using the Huh7 cells stably expressing GFP helped identify the distribution of tumor cells in larvae. Our results imply that acetylcholine accumulation in the microenvironment directly or indirectly supports tumor growth in liver cancer. Use of this model system for drug screening studies holds potential in discovering new cholinergic targets for treatment of liver cancers.

A Combined ULBP2 and SEMA5A Expression Signature as a Prognostic and Predictive Biomarker for Colon Cancer

Background: Prognostic biomarkers for cancer have the power to change the course of disease if they add value beyond known prognostic factors, if they can help shape treatment protocols, and if they are reliable. The aim of this study was to identify such biomarkers for colon cancer and to understand the molecular mechanisms leading to prognostic stratifications based on these biomarkers. Methods and Findings: We used an in house R based script (SSAT) for the in silico discovery of stage-independent prognostic biomarkers using two cohorts, GSE17536 and GSE17537, that include 177 and 55 colon cancer patients, respectively. This identified 2 genes, ULBP2 and SEMA5A, which when used jointly, could distinguish patients with distinct prognosis. We validated our findings using a third cohort of 48 patients ex vivo. We find that in all cohorts, a combined ULBP2/SEMA5A classification (SU-GIB) can stratify distinct prognostic sub-groups with hazard ratios that range from 2.4 to 4.5 (p≤0.01) when overall- or cancer-specific survival is used as an end-measure, independent of confounding prognostic parameters. In addition, our preliminary analyses suggest SU-GIB is comparable to Oncotype DX colon(®) in predicting recurrence in two different cohorts (HR: 1.5-2; p≤0.02). SU-GIB has potential as a companion diagnostic for several drugs including the PI3K/mTOR inhibitor BEZ235, which are suitable for the treatment of patients within the bad prognosis group. We show that tumors from patients with worse prognosis have low EGFR autophosphorylation rates, but high caspase 7 activity, and show upregulation of pro-inflammatory cytokines that relate to a relatively mesenchymal phenotype. Conclusions: We describe two novel genes that can be used to prognosticate colon cancer and suggest approaches by which such tumors can be treated. We also describe molecular characteristics of tumors stratified by the SU-GIB signature.

Functionally conserved effects of rapamycin exposure on zebrafish

Mechanistic target of rapamycin (mTOR) is a conserved serine/threonine kinase important in cell proliferation, growth and protein translation. Rapamycin, a well‑known anti‑cancer agent and immunosuppressant drug, inhibits mTOR activity in different taxa including zebrafish. In the present study, the effect of rapamycin exposure on the transcriptome of a zebrafish fibroblast cell line, ZF4, was investigated. Microarray analysis demonstrated that rapamycin treatment modulated a large set of genes with varying functions including protein synthesis, assembly of mitochondrial and proteasomal machinery, cell cycle, metabolism and oxidative phosphorylation in ZF4 cells. A mild however, coordinated reduction in the expression of proteasomal and mitochondrial ribosomal subunits was detected, while the expression of numerous ribosomal subunits increased. Meta‑analysis of heterogeneous mouse rapamycin microarray datasets enabled the comparison of zebrafish and mouse pathways modulated by rapamycin, using Kyoto Encyclopedia of Genes and Genomes and Gene Ontology pathway analysis. The analyses demonstrated a high degree of functional conservation between zebrafish and mice in response to rapamycin. In addition, rapamycin treatment resulted in a marked dose‑dependent reduction in body size and pigmentation in zebrafish embryos. The present study is the first, to the best of our knowledge, to evaluate the conservation of rapamycin‑modulated functional pathways between zebrafish and mice, in addition to the dose‑dependent growth curves of zebrafish embryos upon rapamycin exposure.

Characterization of a novel zebrafish (Danio rerio) gene, wdr81, associated with cerebellar ataxia, mental retardation and dysequilibrium syndrome (CAMRQ)

Background

WDR81 (WD repeat-containing protein 81) is associated with cerebellar ataxia, mental retardation and disequilibrium syndrome (CAMRQ2, [MIM 610185]). Human and mouse studies suggest that it might be a gene of importance during neurodevelopment. This study aimed at fully characterizing the structure of the wdr81 transcript, detecting the possible transcript variants and revealing its expression profile in zebrafish, a powerful model organism for studying development and disease.

Results

As expected in human and mouse orthologous proteins, zebrafish wdr81 is predicted to possess a BEACH (Beige and Chediak-Higashi) domain, a major facilitator superfamily domain and WD40-repeats, which indicates a conserved function in these species. We observed that zebrafish wdr81 encodes one open reading frame while the transcript has one 5′ untranslated region (UTR) and the prediction of the 3′ UTR was mainly confirmed along with a detected insertion site in the embryo and adult brain. This insertion site was also found in testis, heart, liver, eye, tail and muscle, however, there was no amplicon in kidney, intestine and gills, which might be the result of possible alternative polyadenylation processes among tissues. The 5 and 18 hpf were critical timepoints of development regarding wdr81 expression. Furthermore, the signal of the RNA probe was stronger in the eye and brain at 18 and 48 hpf, then decreased at 72 hpf. Finally, expression of wdr81 was detected in the adult brain and eye tissues, including but not restricted to photoreceptors of the retina, presumptive Purkinje cells and some neurogenic brains regions.

Conclusions

Taken together these data emphasize the importance of this gene during neurodevelopment and a possible role for neuronal proliferation. Our data provide a basis for further studies to fully understand the function of wdr81.

Electronic supplementary material

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

Transgelin gene is frequently downregulated by promoter DNA hypermethylation in breast cancer

Background

CpG hypermethylation in gene promoters is a frequent mechanism of tumor suppressor gene silencing in various types of cancers. It usually occurs at early steps of cancer progression and can be detected easily, giving rise to development of promising biomarkers for both detection and progression of cancer, including breast cancer. 5-aza-2′-deoxycytidine (AZA) is a DNA demethylating and anti-cancer agent resulting in induction of genes suppressed via DNA hypermethylation.

Results

Using microarray expression profiling of AZA- or DMSO-treated breast cancer and non-tumorigenic breast (NTB) cells, we identified for the first time TAGLN gene as a target of DNA hypermethylation in breast cancer. TAGLN expression was significantly and frequently downregulated via promoter DNA hypermethylation in breast cancer cells compared to NTB cells, and also in 13/21 (61.9 %) of breast tumors compared to matched normal tissues. Analyses of public microarray methylation data showed that TAGLN was also hypermethylated in 63.02 % of tumors compared to normal tissues; relapse-free survival of patients was worse with higher TAGLN methylation; and methylation levels could discriminate between tumors and healthy tissues with 83.14 % sensitivity and 100 % specificity. Additionally, qRT-PCR and immunohistochemistry experiments showed that TAGLN expression was significantly downregulated in two more independent sets of breast tumors compared to normal tissues and was lower in tumors with poor prognosis. Colony formation was increased in TAGLN silenced NTB cells, while decreased in overexpressing BC cells.

Conclusions

TAGLN gene is frequently downregulated by DNA hypermethylation, and TAGLN promoter methylation profiles could serve as a future diagnostic biomarker, with possible clinical impact regarding the prognosis in breast cancer.

Electronic supplementary material

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

Relative expression of rRNA transcripts and 45S rDNA promoter methylation status are dysregulated in tumors in comparison with matched-normal tissues in breast cancer

Ribosomal RNA (rRNA) expression, one of the most important factors regulating ribosome production, is primarily controlled by a CG-rich 45 S rDNA promoter. However, the DNA methylation state of the 45 S rDNA promoter, as well as its effect on rRNA gene expression in types of human cancers is controversial. In the present study we analyzed the methylation status of the rDNA promoter (-380 to +53 bp) as well as associated rRNA expression levels in breast cancer cell lines and breast tumor-normal tissue pairs. We found that the aforementioned regulatory region was extensively methylated (74-96%) in all cell lines and in 68% (13/19 tumor-normal pairs) of the tumors. Expression levels of rRNA transcripts 18 S, 28 S, 5.8 S and 45 S external transcribed spacer (45 S ETS) greatly varied in the breast cancer cell lines regardless of their methylation status. Analyses of rRNA transcript expression levels in the breast tumor and normal matched tissues showed no significant difference when normalized with TBP. On the other hand, using the geometric mean of the rRNA expression values (GM-rRNA) as reference enabled us to identify significant changes in the relative expression of rRNAs in the tissue samples. We propose GM-rRNA normalization as a novel strategy to analyze expression differences between rRNA transcripts. Accordingly, the 18S rRNA/GM-rRNA ratio was significantly higher whereas the 5.8S rRNA/GM-rRNA ratio was significantly lower in breast tumor samples than this ratio in the matched normal samples. Moreover, the 18S rRNA/GM-rRNA ratio was negatively correlated with the 45 S rDNA promoter methylation level in the normal breast tissue samples, yet not in the breast tumors. Significant correlations observed between the expression levels of rRNA transcripts in the normal samples were lost in the tumor samples. We showed that the expression of rRNA transcripts may not be based solely on promoter methylation. Carcinogenesis may cause dysregulation of the correlation between spliced rRNA expression levels, possibly due to changes in rRNA processing, which requires further investigation.

Identification of novel reference genes based on MeSH categories

Transcriptome experiments are performed to assess protein abundance through mRNA expression analysis. Expression levels of genes vary depending on the experimental conditions and the cell response. Transcriptome data must be diverse and yet comparable in reference to stably expressed genes, even if they are generated from different experiments on the same biological context from various laboratories. In this study, expression patterns of 9090 microarray samples grouped into 381 NCBI-GEO datasets were investigated to identify novel candidate reference genes using randomizations and Receiver Operating Characteristic (ROC) curves. The analysis demonstrated that cell type specific reference gene sets display less variability than a united set for all tissues. Therefore, constitutively and stably expressed, origin specific novel reference gene sets were identified based on their coefficient of variation and percentage of occurrence in all GEO datasets, which were classified using Medical Subject Headings (MeSH). A large number of MeSH grouped reference gene lists are presented as novel tissue specific reference gene lists. The most commonly observed 17 genes in these sets were compared for their expression in 8 hepatocellular, 5 breast and 3 colon carcinoma cells by RT-qPCR to verify tissue specificity. Indeed, commonly used housekeeping genes GAPDH, Actin and EEF2 had tissue specific variations, whereas several ribosomal genes were among the most stably expressed genes in vitro. Our results confirm that two or more reference genes should be used in combination for differential expression analysis of large-scale data obtained from microarray or next generation sequencing studies. Therefore context dependent reference gene sets, as presented in this study, are required for normalization of expression data from diverse technological backgrounds.

Cancer-testis gene expression is associated with the methylenetetrahydrofolate reductase 677 C>T polymorphism in non-small cell lung carcinoma

Background

Tumor-specific, coordinate expression of cancer-testis (CT) genes, mapping to the X chromosome, is observed in more than 60% of non-small cell lung cancer (NSCLC) patients. Although CT gene expression has been unequivocally related to DNA demethylation of promoter regions, the underlying mechanism leading to loss of promoter methylation remains elusive. Polymorphisms of enzymes within the 1-carbon pathway have been shown to affect S-adenosyl methionine (SAM) production, which is the sole methyl donor in the cell. Allelic variants of several enzymes within this pathway have been associated with altered SAM levels either directly, or indirectly as reflected by altered levels of SAH and Homocysteine levels, and altered levels of DNA methylation. We, therefore, asked whether the five most commonly occurring polymorphisms in four of the enzymes in the 1-carbon pathway associated with CT gene expression status in patients with NSCLC.

Methods

Fifty patients among a cohort of 763 with NSCLC were selected based on CT gene expression status and typed for five polymorphisms in four genes known to affect SAM generation by allele specific q-PCR and RFLP.

Results

We identified a significant association between CT gene expression and the MTHFR 677 CC genotype, as well as the C allele of the SNP, in this cohort of patients. Multivariate analysis revealed that the genotype and allele strongly associate with CT gene expression, independent of potential confounders.

Conclusions

Although CT gene expression is associated with DNA demethylation, in NSCLC, our data suggests this is unlikely to be the result of decreased MTHFR function.

Genome-wide transcriptional reorganization associated with senescence-to-immortality switch during human hepatocellular carcinogenesis

Senescence is a permanent proliferation arrest in response to cell stress such as DNA damage. It contributes strongly to tissue aging and serves as a major barrier against tumor development. Most tumor cells are believed to bypass the senescence barrier (become “immortal”) by inactivating growth control genes such as TP53 and CDKN2A. They also reactivate telomerase reverse transcriptase. Senescence-to-immortality transition is accompanied by major phenotypic and biochemical changes mediated by genome-wide transcriptional modifications. This appears to happen during hepatocellular carcinoma (HCC) development in patients with liver cirrhosis, however, the accompanying transcriptional changes are virtually unknown. We investigated genome-wide transcriptional changes related to the senescence-to-immortality switch during hepatocellular carcinogenesis. Initially, we performed transcriptome analysis of senescent and immortal clones of Huh7 HCC cell line, and identified genes with significant differential expression to establish a senescence-related gene list. Through the analysis of senescence-related gene expression in different liver tissues we showed that cirrhosis and HCC display expression patterns compatible with senescent and immortal phenotypes, respectively; dysplasia being a transitional state. Gene set enrichment analysis revealed that cirrhosis/senescence-associated genes were preferentially expressed in non-tumor tissues, less malignant tumors, and differentiated or senescent cells. In contrast, HCC/immortality genes were up-regulated in tumor tissues, or more malignant tumors and progenitor cells. In HCC tumors and immortal cells genes involved in DNA repair, cell cycle, telomere extension and branched chain amino acid metabolism were up-regulated, whereas genes involved in cell signaling, as well as in drug, lipid, retinoid and glycolytic metabolism were down-regulated. Based on these distinctive gene expression features we developed a 15-gene hepatocellular immortality signature test that discriminated HCC from cirrhosis with high accuracy. Our findings demonstrate that senescence bypass plays a central role in hepatocellular carcinogenesis engendering systematic changes in the transcription of genes regulating DNA repair, proliferation, differentiation and metabolism.

mESAdb: microRNA expression and sequence analysis database

microRNA expression and sequence analysis database (http://konulab.fen.bilkent.edu.tr/mirna/) (mESAdb) is a regularly updated database for the multivariate analysis of sequences and expression of microRNAs from multiple taxa. mESAdb is modular and has a user interface implemented in PHP and JavaScript and coupled with statistical analysis and visualization packages written for the R language. The database primarily comprises mature microRNA sequences and their target data, along with selected human, mouse and zebrafish expression data sets. mESAdb analysis modules allow (i) mining of microRNA expression data sets for subsets of microRNAs selected manually or by motif; (ii) pair-wise multivariate analysis of expression data sets within and between taxa; and (iii) association of microRNA subsets with annotation databases, HUGE Navigator, KEGG and GO. The use of existing and customized R packages facilitates future addition of data sets and analysis tools. Furthermore, the ability to upload and analyze user-specified data sets makes mESAdb an interactive and expandable analysis tool for microRNA sequence and expression data.

Dose- and time-dependent expression patterns of zebrafish orthologs of selected E2F target genes in response to serum starvation/replenishment

Targets of E2F transcription factors effectively regulate the cell cycle from worms to humans. Furthermore, the dysregulation of E2F transcription modules plays a highly conserved role in cancers of human and zebrafish. Studying E2F target expression under a given cellular state, such as quiescence, might lead to a better understanding of the conserved patterns of expression in different taxa. In the present study, we used literature searches and phylogeny to identify several targets of E2F transcription factors that are known to be serum-responsive; namely, PCNA, MYBL2, MCM7, TYMS, and CTGF. The transcriptional serum response of zebrafish orthologs of these genes were quantified under different doses (i.e., 0, 0.1, 1, 3, and 10% FBS) and time points (i.e., 6, 24 and 48 hours, h) using quantitative RT-PCR (qRT-PCR) in the zebrafish fibroblast cells (ZF4). Our results indicated that mRNA expression of zebrafish pcna, mybl2, mcm7 and tyms drastically decreased while that of ctgf increased with decreasing serum levels as observed in mammals. These genes responded to serum starvation at 24 and 48 h and to the mitogenic stimuli as early as 6 h except for ctgf whose expression was significantly altered at 24 h. The zebrafish Mcm7 protein levels also were modulated by serum starvation/replenishment. The present study provides a foundation for the comparative analysis of quantitative expression patterns for genes involved in regulation of cell cycle using a zebrafish serum response model.

Identification of endogenous reference genes for qRT-PCR analysis in normal matched breast tumor tissues

Quantitative gene expression measurements from tumor tissue are frequently compared with matched normal and/or adjacent tumor tissue expression for diagnostic marker gene selection as well as assessment of the degree of transcriptional deregulation in cancer. Selection of an appropriate reference gene (RG) or an RG panel, which varies depending on cancer type, molecular subtypes, and the normal tissues used for interindividual calibration, is crucial for the accurate quantification of gene expression. Several RG panels have been suggested in breast cancer for making comparisons among tumor subtypes, cell lines, and benign/malignant tumors. In this study, expression patterns of 15 widely used endogenous RGs (ACTB, TBP, GAPDH, SDHA, HPRT, HMBS, B2M, PPIA, GUSB, YWHAZ2, PGK1, RPLP0, PUM1, MRPL19, and RPL41), and three candidate genes that were selected through analysis of two independent microarray datasets (IL22RA1, TC22, ZNF224) were determined in 23 primary breast tumors and their matched normal tissues using qRT-PCR. Additionally, 18S rRNA, ACTB, and SDHA were tested using randomly primed cDNAs from 13 breast tumor pairs to assess the rRNA/mRNA ratio. The tumors exhibited significantly lower rRNA/mRNA ratio when compared to their normals, on average. The expression of the studied RGs in breast tumors did not exhibit differences in terms of grade, ER, or PR status. The stability of RGs was examined based on two different statistical models, namely GeNorm and NormFinder. Among the 18 tested endogenous reference genes, ACTB and SDHA were identified as the most suitable reference genes for the normalization of qRT-PCR data in the analysis of normal matched tumor breast tissue pairs by both programs. In addition, the expression of the gelsolin (GSN) gene, a well-known downregulated target in breast tumors, was analyzed using the two most suitable genes and different RG combinations to validate their effectiveness as a normalization factor (NF). The GSN expression of the tumors used in this study was significantly lower than that of normals showing the effectivity of using ACTB and SDHA as suitable RGs in this set of tumor-normal tissue panel. The combinational use of the best performing two RGs (ACTB and SDHA) as a normalization factor can be recommended to minimize sample variability and to increase the accuracy and resolution of gene expression normalization in tumor-normal paired breast cancer qRT-PCR studies.

A resampling-based meta-analysis for detection of differential gene expression in breast cancer

Background

Accuracy in the diagnosis of breast cancer and classification of cancer subtypes has improved over the years with the development of well-established immunohistopathological criteria. More recently, diagnostic gene-sets at the mRNA expression level have been tested as better predictors of disease state. However, breast cancer is heterogeneous in nature; thus extraction of differentially expressed gene-sets that stably distinguish normal tissue from various pathologies poses challenges. Meta-analysis of high-throughput expression data using a collection of statistical methodologies leads to the identification of robust tumor gene expression signatures.

Methods

A resampling-based meta-analysis strategy, which involves the use of resampling and application of distribution statistics in combination to assess the degree of significance in differential expression between sample classes, was developed. Two independent microarray datasets that contain normal breast, invasive ductal carcinoma (IDC), and invasive lobular carcinoma (ILC) samples were used for the meta-analysis. Expression of the genes, selected from the gene list for classification of normal breast samples and breast tumors encompassing both the ILC and IDC subtypes were tested on 10 independent primary IDC samples and matched non-tumor controls by real-time qRT-PCR. Other existing breast cancer microarray datasets were used in support of the resampling-based meta-analysis.

Results

The two independent microarray studies were found to be comparable, although differing in their experimental methodologies (Pearson correlation coefficient, R = 0.9389 and R = 0.8465 for ductal and lobular samples, respectively). The resampling-based meta-analysis has led to the identification of a highly stable set of genes for classification of normal breast samples and breast tumors encompassing both the ILC and IDC subtypes. The expression results of the selected genes obtained through real-time qRT-PCR supported the meta-analysis results.

Conclusion

The proposed meta-analysis approach has the ability to detect a set of differentially expressed genes with the least amount of within-group variability, thus providing highly stable gene lists for class prediction. Increased statistical power and stringent filtering criteria used in the present study also make identification of novel candidate genes possible and may provide further insight to improve our understanding of breast cancer development.

Quantification of SLIT-ROBO transcripts in hepatocellular carcinoma reveals two groups of genes with coordinate expression

BACKGROUND

SLIT-ROBO families of proteins mediate axon pathfinding and their expression is not solely confined to nervous system. Aberrant expression of SLIT-ROBO genes was repeatedly shown in a wide variety of cancers, yet data about their collective behavior in hepatocellular carcinoma (HCC) is missing. Hence, we quantified SLIT-ROBO transcripts in HCC cell lines, and in normal and tumor tissues from liver.

METHODS

Expression of SLIT-ROBO family members was quantified by real-time qRT-PCR in 14 HCC cell lines, 8 normal and 35 tumor tissues from the liver. ANOVA and Pearson's correlation analyses were performed in R environment, and different clinicopathological subgroups were pairwise compared in Minitab. Gene expression matrices of cell lines and tissues were analyzed by Mantel's association test.

RESULTS

Genewise hierarchical clustering revealed two subgroups with coordinate expression pattern in both the HCC cell lines and tissues: ROBO1, ROBO2, SLIT1 in one cluster, and ROBO4, SLIT2, SLIT3 in the other, respectively. Moreover, SLIT-ROBO expression predicted AFP-dependent subgrouping of HCC cell lines, but not that of liver tissues. ROBO1 and ROBO2 were significantly up-regulated, whereas SLIT3 was significantly down-regulated in cell lines with high-AFP background. When compared to normal liver tissue, ROBO1 was found to be significantly overexpressed, while ROBO4 was down-regulated in HCC. We also observed that ROBO1 and SLIT2 differentiated histopathological subgroups of liver tissues depending on both tumor staging and differentiation status. However, ROBO4 could discriminate poorly differentiated HCC from other subgroups.

CONCLUSION

The present study is the first in comprehensive and quantitative evaluation of SLIT-ROBO family gene expression in HCC, and suggests that the expression of SLIT-ROBO genes is regulated in hepatocarcinogenesis. Our results implicate that SLIT-ROBO transcription profile is bi-modular in nature, and that each module shows intrinsic variability. We also provide quantitative evidence for potential use of ROBO1, ROBO4 and SLIT2 for prediction of tumor stage and differentiation status.

Cloning and expression profile of FLT3 gene during progenitor cell-dependent liver regeneration

BACKGROUND AND AIM

The liver has a unique capacity to regenerate upon exposure to viral infections, toxic reactions and cancer formation. Liver regeneration is a complex phenomenon in which several factors participate during its onset. Cellular proliferation is an important component of this process and the factors that regulate this proliferation have a vital role. FLT3, a well-known hematopoietic stem cell and hepatic lineage surface marker, is involved in proliferative events of hematopoietic stem cells. However, its contribution to liver regeneration is not known. Therefore, the aim of this study was to clone and examine the role of FLT3 during liver regeneration in rats.

METHODS

Partial cDNA of rat homolog of FLT3 gene was cloned from thymus and the tissue specific expression of this gene at mRNA and protein levels was examined by RT-PCR and Western blot. After treating with 2-AAF and performing hepatectomy in rats to induce progenitor-dependent liver regeneration, the mRNA and protein expression profile of FLT3 was investigated by real-time PCR and Western blot during liver regeneration. In addition, cellular localization of FLT3 protein was determined by immunohistochemistry.

RESULTS

The results indicated that rat FLT3 cDNA has high homology with mouse and human FLT3 cDNA. It was also found that FLT3 is expressed in most of the rat tissues and during liver regeneration. In addition, its intracellular localization is altered during the late stages of liver regeneration.

CONCLUSION

The FLT3 receptor is activated at the late stages of liver regeneration and participates in the proliferation response that is observed during progenitor-dependent liver regeneration.

Strain- and region-specific gene expression profiles in mouse brain in response to chronic nicotine treatment

A pathway-focused complementary DNA microarray and gene ontology analysis were used to investigate gene expression profiles in the amygdala, hippocampus, nucleus accumbens, prefrontal cortex (PFC) and ventral tegmental area of C3H/HeJ and C57BL/6J mice receiving nicotine in drinking water (100 mug/ml in 2% saccharin for 2 weeks). A balanced experimental design and rigorous statistical analysis have led to the identification of 3.5-22.1% and 4.1-14.3% of the 638 sequence-verified genes as significantly modulated in the aforementioned brain regions of the C3H/HeJ and C57BL/6J strains, respectively. Comparisons of differential expression among brain tissues showed that only a small number of genes were altered in multiple brain regions, suggesting presence of a brain region-specific transcriptional response to nicotine. Subsequent principal component analysis and Expression Analysis Systematic Explorer analysis showed significant enrichment of biological processes both in C3H/HeJ and C57BL/6J mice, i.e. cell cycle/proliferation, organogenesis and transmission of nerve impulse. Finally, we verified the observed changes in expression using real-time reverse transcriptase polymerase chain reaction for six representative genes in the PFC region, providing an independent replication of our microarray results. Together, this report represents the first comprehensive gene expression profiling investigation of the changes caused by nicotine in brain tissues of the two mouse strains known to exhibit differential behavioral and physiological responses to nicotine.

Alternatively spliced Robo2 isoforms in zebrafish and rat

Robo2, a member of the robo gene family, functions as a repulsive axon guidance receptor as well as a regulator of cell migration and tissue morphogenesis in different taxa. In this study, a novel isoform of the zebrafish robo2 (robo2_tv2), which included an otherwise alternatively spliced exon (CAE), has been characterized. Robo2_tv2 is expressed differentially in most non-neuronal tissues of adult zebrafish whereas robo2_tv1 expression to a great extent is restricted to the brain and eye. In zebrafish, robo2_tv2 exhibits a very-low-level basal expression starting from 1 day post fertilization until the mid-larval stages, at which time its expression increases dramatically and could be detected throughout adulthood. Our findings demonstrate that the amino acid sequence coded by CAE of the robo2 gene is highly conserved between zebrafish and mammals, and also contains conserved motifs shared with robo1 and robo4 but not with robo3. Furthermore, we provide an account of differential transcription of the CAE homolog in various tissues of the adult rat. These results suggest that the alternatively spliced robo2 isoforms may exhibit tissue specificity.

Regulation of Homer and group I metabotropic glutamate receptors by nicotine

The present study focuses on the nicotine-induced modulation of mRNA and protein expression of a number of genes involved in glutamatergic synaptic transmission in rat brain over different time periods of exposure. A subchronic (3 days) but not the chronic (7 or 14 days) administration of nicotine resulted in the up-regulation of Homer2a/b mRNA in the amygdala while in the ventral tegmental area (VTA) no change in expression of either Homer2a/b or Homer1b/c was observed. Although the increase in Homer2a/b mRNA was not translated into the protein level in the amygdala, a slight but significant up-regulation of Homer1b/c protein was observed in the same region at day 3. Both Homer forms were up-regulated at the protein level in the VTA at day 3. In the nucleus accumbens, 14 days of nicotine treatment up-regulated mRNA of Homer2b/c by 68.2% (P < 0.05), while the short form Homer1a gene was down-regulated by 65.0% at day 3 (P < 0.05). In regard to other components of the glutamatergic signalling, we identified an acute and intermittent increase in the mRNA and protein levels of mGluR1 and mGluR5 in the amygdala. In the VTA, however, the effects of nicotine on mGluR mRNA expression were long-lasting but rather specific to mGluR1. Nevertheless, mGluR1 protein levels in the VTA area were up-regulated only at day 3, as in the amygdala. These data provide further evidence for the involvement of nicotine in the glutamatergic neuronal synaptic activity in vivo, suggesting a role for the newly identified Homer proteins in this paradigm.

Nicotine coregulates multiple pathways involved in protein modification/degradation in rat brain

Previously, we used cDNA microarrays to demonstrate that the phosphatidylinositol and MAP kinase signaling pathways are regulated by nicotine in different rat brain regions. In the present report, we show that, after exposure to nicotine for 14 days, ubiquitin, ubiquitin-conjugating enzymes, 20S and 19S proteasomal subunits, and chaperonin-containing TCP-1 protein (CCT) complex members are upregulated in rat prefrontal cortex (PFC) while being downregulated in the medial basal hypothalamus (MBH). In particular, relative to saline controls, ubiquitins B and C were upregulated by 33% and 47% (P<0.01), respectively, in the PFC. The proteasome beta subunit 1 (PSMB1) and 26S ATPase 3 (PSMC3) genes were upregulated in the PFC by 95% and 119% (P<0.001), respectively. In addition to the protein degradation pathway of the ubiquitin-proteasome complexes, we observed in the PFC an increase in the expression of small, ubiquitin-related modifiers (SUMO) 1 and 2 by 80% and 33%, respectively (P<0.001), and in 3 of 6 CCT subunits by up to 150% (P<0.0001). To a lesser extent, a change in the opposite direction was obtained in the expression of the same gene families in the MBH. Quantitative real-time RT-PCR was used to validate the microarray results obtained with some representative genes involved in these pathways. Taken together, our results suggest that, in response to systemic nicotine administration, the ubiquitin-proteasome, SUMO, and chaperonin complexes provide an intricate control mechanism to maintain cellular homeostasis, possibly by regulating the composition and signaling of target neurons in a region-specific manner.

Application of a customized pathway-focused microarray for gene expression profiling of cellular homeostasis upon exposure to nicotine in PC12 cells

Maintenance of cellular homeostasis is integral to appropriate regulation of cellular signaling and cell growth and division. In this study, we report the development and quality assessment of a pathway-focused microarray comprising genes involved in cellular homeostasis. Since nicotine is known to have highly modulatory effects on the intracellular calcium homeostasis, we therefore tested the applicability of the homeostatic pathway-focused microarray on the gene expression in PC-12 cells treated with 1 mM nicotine for 48 h relative to the untreated control cells. We first provided a detailed description of the focused array with respect to its gene and pathway content and then assessed the array quality using a robust regression procedure that allows for the exclusion of unreliable measurements while decreasing the number of false positives. As a result, the mean correlation coefficient between duplicate measurements of the arrays used in this study (control vs. nicotine treatment, three samples each) has increased from 0.974+/-0.017 to 0.995+/-0.002. Furthermore, we found that nicotine affected various structural and signaling components of the AKT/PKB signaling pathway and protein synthesis and degradation processes in PC-12 cells. Since modulation of intracellular calcium concentrations ([Ca(2+)](i)) and phosphatidylinositol signaling are important in various biological processes such as neurotransmitter release and tissue pathogenesis including tumor formation, we expect that the homeostatic pathway-focused microarray potentially can be used for the identification of unique gene expression profiles in comparative studies of drugs of abuse and diverse environmental stimuli, such as starvation and oxidative stress.

Nicotine, body weight and potential implications in the treatment of obesity

Obesity is an epidemic problem in the U. S. and many other industrialized nations. Historically, the drugs used for the treatment of obesity generally targeted small molecule neurotransmitters. As research grows to decipher the underlying molecular mechanisms behind energy homeostasis, it is becoming evident that the modulating effects of neuropeptides also are critical in the regulation of appetite and metabolism. The search for drugs to modify these monoaminergic and peptidergic pathways may eventually prove successful in the treatment of obesity. While tobacco smoking has long been used as one strategy to maintain a lower body weight, especially in female smokers, its adverse associations with addiction and disease overshadow its potential use as an antiobesity agent. Potential pharmacological effects of nicotine could be better understood as the intricacies of the nicotinic acetylcholine receptor are revealed. The objective of this review is threefold: first is to provide an overview of the physiological effects of nicotine on body weight while focusing on the drugs that are available as antiobesity and smoking cessation agents. Second is to provide the present status of the nicotinic acetylcholine receptor as it relates to energy homeostasis and its potential as an effective treatment modality for obesity. Third is to present the current knowledge with respect to nicotine's effects on energy homeostatic and reward related pathways at the molecular level. A better understanding of the regulatory mechanisms underlying the pharmacological effects of nicotine on body weight will provide insights in identification of potential targets for the development of appropriate medicines in the treatment of obesity.

Correlations between mRNA expression levels and GC contents of coding and untranslated regions of genes in rodents

Gene expression is regulated by a highly coordinated network of events whose efficiency may constrain the level of expression. Among other factors, natural selection for increased translational efficiency and/or fidelity may shape nucleotide composition and, hence, codon usage during evolution. Previous studies have shown that highly expressed genes in Saccharomyces cerevisiae, Caenorhabditis elegans, and Drosophila melanogaster have relatively higher codon usage biases. However, in the case of mammals, results have been equivocal. In this study, we assessed the correlation between nucleotide composition and mRNA expression levels of rodent genes measured by cDNA microarray and serial analysis of gene expression (SAGE) techniques. We found that mRNA expression levels were correlated with the third nucleotide position GC (GC3) content for both Rattus norvegicus (r = 0.246, p = 0.01; N = 110) and Mus musculus (r = 0.21, p = 0.0026; N = 203) genes. However, no significant correlation was evident between mRNA expression level and GC contents of 5'- and 3'-untranslated regions (UTRs) for either species. This suggests that, in rodents, nucleotide composition of coding sequences and UTRs might evolve differentially when considered along an expression gradient. Accordingly, it is possible that higher GC levels may present the rodent genes with a selective advantage for translational efficiency. However, the increase in GC3 content seems to level off above an expressional threshold (e.g., >or=threefold the median expression for R. norvegicus), suggesting that conflicting demands posed by different aspects of transcriptional and translational machineries (e.g., efficiency versus fidelity) may set an upper limit for GC3.