Interferon-Induced Gene Expression Is a Stronger Predictor of Treatment Response Than IL28B Genotype in Patients With Hepatitis C.

BACKGROUND & AIMS: The host immune response during the chronic phase of hepatitis C virus infection varies among individuals; some patients have a no interferon (IFN) response in the liver, whereas others have full activation IFN-stimulated genes (ISGs). Preactivation of this endogenous IFN system is associated with nonresponse to pegylated IFN-α (pegIFN-α) and ribavirin. Genome-wide association studies have associated allelic variants near the IL28B (IFNλ3) gene with treatment response. We investigated whether IL28B genotype determines the constitutive expression of ISGs in the liver and compared the abilities of ISG levels and IL28B genotype to predict treatment outcome. METHODS: We genotyped 109 patients with chronic hepatitis C for IL28B allelic variants and quantified the hepatic expression of ISGs and of IL28B. Decision tree ensembles, in the form of a random forest classifier, were used to calculate the relative predictive power of these different variables in a multivariate analysis. RESULTS: The minor IL28B allele was significantly associated with increased expression of ISG. However, stratification of the patients according to treatment response revealed increased ISG expression in nonresponders, irrespective of IL28B genotype. Multivariate analysis of ISG expression, IL28B genotype, and several other factors associated with response to therapy identified ISG expression as the best predictor of treatment response. CONCLUSIONS: IL28B genotype and hepatic expression of ISGs are independent predictors of response to treatment with pegIFN-α and ribavirin in patients with chronic hepatitis C. The most accurate prediction of response was obtained with a 4-gene classifier comprising IFI27, ISG15, RSAD2, and HTATIP2.

Cross-species and cross-platform gene expression studies with the Bioconductor-compliant R package 'annotationTools'.

Background: The variety of DNA microarray formats and datasets presently available offers an unprecedented opportunity to perform insightful comparisons of heterogeneous data. Cross-species studies, in particular, have the power of identifying conserved, functionally important molecular processes. Validation of discoveries can now often be performed in readily available public data which frequently requires cross-platform studies.Cross-platform and cross-species analyses require matching probes on different microarray formats. This can be achieved using the information in microarray annotations and additional molecular biology databases, such as orthology databases. Although annotations and other biological information are stored using modern database models ( e. g. relational), they are very often distributed and shared as tables in text files, i.e. flat file databases. This common flat database format thus provides a simple and robust solution to flexibly integrate various sources of information and a basis for the combined analysis of heterogeneous gene expression profiles.Results: We provide annotationTools, a Bioconductor-compliant R package to annotate microarray experiments and integrate heterogeneous gene expression profiles using annotation and other molecular biology information available as flat file databases. First, annotationTools contains a specialized set of functions for mining this widely used database format in a systematic manner. It thus offers a straightforward solution for annotating microarray experiments. Second, building on these basic functions and relying on the combination of information from several databases, it provides tools to easily perform cross-species analyses of gene expression data.Here, we present two example applications of annotationTools that are of direct relevance for the analysis of heterogeneous gene expression profiles, namely a cross-platform mapping of probes and a cross-species mapping of orthologous probes using different orthology databases. We also show how to perform an explorative comparison of disease-related transcriptional changes in human patients and in a genetic mouse model.Conclusion: The R package annotationTools provides a simple solution to handle microarray annotation and orthology tables, as well as other flat molecular biology databases. Thereby, it allows easy integration and analysis of heterogeneous microarray experiments across different technological platforms or species.

Genomic organization and gene expression in a chromosomal region of Leishmania major.

Little is known about the relation between the genome organization and gene expression in Leishmania. Bioinformatic analysis can be used to predict genes and find homologies with known proteins. A model was proposed, in which genes are organized into large clusters and transcribed from only one strand, in the form of large polycistronic primary transcripts. To verify the validity of this model, we studied gene expression at the transcriptional, post-transcriptional and translational levels in a unique locus of 34kb located on chr27 and represented by cosmid L979. Sequence analysis revealed 115 ORFs on either DNA strand. Using computer programs developed for Leishmania genes, only nine of these ORFs, localized on the same strand, were predicted to code for proteins, some of which show homologies with known proteins. Additionally, one pseudogene, was identified. We verified the biological relevance of these predictions. mRNAs from nine predicted genes and proteins from seven were detected. Nuclear run-on analyses confirmed that the top strand is transcribed by RNA polymerase II and suggested that there is no polymerase entry site. Low levels of transcription were detected in regions of the bottom strand and stable transcripts were identified for four ORFs on this strand not predicted to be protein-coding. In conclusion, the transcriptional organization of the Leishmania genome is complex, raising the possibility that computer predictions may not be comprehensive.

Impact of salt on cardiac differential gene expression and coronary lesion in normotensive mineralocorticoid-treated mice.

We previously reported that excess of deoxycorticosterone-acetate (DOCA)/salt-induced cardiac hypertrophy in the absence of hypertension in one-renin gene mice. This model allows us to study molecular mechanisms of high-salt intake in the development of cardiovascular remodeling, independently of blood pressure in a high mineralocorticoid state. In this study, we compared the effect of 5-wk low- and high-salt intake on cardiovascular remodeling and cardiac differential gene expression in mice receiving the same amount of DOCA. Differential gene and protein expression was measured by high-density cDNA microarray assays, real-time PCR and Western blot analysis in DOCA-high salt (HS) vs. DOCA-low salt (LS) mice. DOCA-HS mice developed cardiac hypertrophy, coronary perivascular fibrosis, and left ventricular dysfunction. Differential gene and protein expression demonstrated that high-salt intake upregulated a subset of genes encoding for proteins involved in inflammation and extracellular matrix remodeling (e.g., Col3a1, Col1a2, Hmox1, and Lcn2). A major subset of downregulated genes encoded for transcription factors, including myeloid differentiation primary response (MyD) genes. Our data provide some evidence that vascular remodeling, fibrosis, and inflammation are important consequences of a high-salt intake in DOCA mice. Our study suggests that among the different pathogenic factors of cardiac and vascular remodeling, such as hypertension and mineralocorticoid excess and sodium intake, the latter is critical for the development of the profibrotic and proinflammatory phenotype observed in the heart of normotensive DOCA-treated mice.

Genome-wide analysis of cancer/testis gene expression.

Cancer/Testis (CT) genes, normally expressed in germ line cells but also activated in a wide range of cancer types, often encode antigens that are immunogenic in cancer patients, and present potential for use as biomarkers and targets for immunotherapy. Using multiple in silico gene expression analysis technologies, including twice the number of expressed sequence tags used in previous studies, we have performed a comprehensive genome-wide survey of expression for a set of 153 previously described CT genes in normal and cancer expression libraries. We find that although they are generally highly expressed in testis, these genes exhibit heterogeneous gene expression profiles, allowing their classification into testis-restricted (39), testis/brain-restricted (14), and a testis-selective (85) group of genes that show additional expression in somatic tissues. The chromosomal distribution of these genes confirmed the previously observed dominance of X chromosome location, with CT-X genes being significantly more testis-restricted than non-X CT. Applying this core classification in a genome-wide survey we identified >30 CT candidate genes; 3 of them, PEPP-2, OTOA, and AKAP4, were confirmed as testis-restricted or testis-selective using RT-PCR, with variable expression frequencies observed in a panel of cancer cell lines. Our classification provides an objective ranking for potential CT genes, which is useful in guiding further identification and characterization of these potentially important diagnostic and therapeutic targets.

Jasmonate and salicylate as global signals for defense gene expression.

Remarkably, only a few low molecular mass signals, including jasmonic acid, ethylene and salicylic acid, upregulate the expression of scores of defense-related genes. Using these regulators, the plant fine-tunes its defense gene expression against aggressors which, in some cases, may be able to disrupt or amplify plant defense signal pathways to their own ends.

Network-guided analysis of genes with altered somatic copy number and gene expression reveals pathways commonly perturbed in metastatic melanoma.

Cancer genomes frequently contain somatic copy number alterations (SCNA) that can significantly perturb the expression level of affected genes and thus disrupt pathways controlling normal growth. In melanoma, many studies have focussed on the copy number and gene expression levels of the BRAF, PTEN and MITF genes, but little has been done to identify new genes using these parameters at the genome-wide scale. Using karyotyping, SNP and CGH arrays, and RNA-seq, we have identified SCNA affecting gene expression ('SCNA-genes') in seven human metastatic melanoma cell lines. We showed that the combination of these techniques is useful to identify candidate genes potentially involved in tumorigenesis. Since few of these alterations were recurrent across our samples, we used a protein network-guided approach to determine whether any pathways were enriched in SCNA-genes in one or more samples. From this unbiased genome-wide analysis, we identified 28 significantly enriched pathway modules. Comparison with two large, independent melanoma SCNA datasets showed less than 10% overlap at the individual gene level, but network-guided analysis revealed 66% shared pathways, including all but three of the pathways identified in our data. Frequently altered pathways included WNT, cadherin signalling, angiogenesis and melanogenesis. Additionally, our results emphasize the potential of the EPHA3 and FRS2 gene products, involved in angiogenesis and migration, as possible therapeutic targets in melanoma. Our study demonstrates the utility of network-guided approaches, for both large and small datasets, to identify pathways recurrently perturbed in cancer.

Whole genome analysis of p38 SAPK-mediated gene expression upon stress

Background: Cells have the ability to respond and adapt to environmental changes through activation of stress-activated protein kinases (SAPKs). Although p38 SAPK signalling is known to participate in the regulation of gene expression little is known on the molecular mechanisms used by this SAPK to regulate stress-responsive genes and the overall set of genes regulated by p38 in response to different stimuli.Results: Here, we report a whole genome expression analyses on mouse embryonic fibroblasts (MEFs) treated with three different p38 SAPK activating-stimuli, namely osmostress, the cytokine TNFα and the protein synthesis inhibitor anisomycin. We have found that the activation kinetics of p38α SAPK in response to these insults is different and also leads to a complex gene pattern response specific for a given stress with a restricted set of overlapping genes. In addition, we have analysed the contribution of p38α the major p38 family member present in MEFs, to the overall stress-induced transcriptional response by using both a chemical inhibitor (SB203580) and p38α deficient (p38α-/-) MEFs. We show here that p38 SAPK dependency ranged between 60% and 88% depending on the treatments and that there is a very good overlap between the inhibitor treatment and the ko cells. Furthermore, we have found that the dependency of SAPK varies depending on the time the cells are subjected to osmostress. Conclusions: Our genome-wide transcriptional analyses shows a selective response to specific stimuli and a restricted common response of up to 20% of the stress up-regulated early genes that involves an important set of transcription factors, which might be critical for either cell adaptation or preparation for continuous extra-cellular changes. Interestingly, up to 85% of the up-regulated genes are under the transcriptional control of p38 SAPK. Thus, activation of p38 SAPK is critical to elicit the early gene expression program required for cell adaptation to stress.

Sleep deprivation effects on circadian clock gene expression in the cerebral cortex parallel electroencephalographic differences among mouse strains.

Sleep deprivation (SD) results in increased electroencephalographic (EEG) delta power during subsequent non-rapid eye movement sleep (NREMS) and is associated with changes in the expression of circadian clock-related genes in the cerebral cortex. The increase of NREMS delta power as a function of previous wake duration varies among inbred mouse strains. We sought to determine whether SD-dependent changes in circadian clock gene expression parallel this strain difference described previously at the EEG level. The effects of enforced wakefulness of incremental durations of up to 6 h on the expression of circadian clock genes (bmal1, clock, cry1, cry2, csnk1epsilon, npas2, per1, and per2) were assessed in AKR/J, C57BL/6J, and DBA/2J mice, three strains that exhibit distinct EEG responses to SD. Cortical expression of clock genes subsequent to SD was proportional to the increase in delta power that occurs in inbred strains: the strain that exhibits the most robust EEG response to SD (AKR/J) exhibited dramatic increases in expression of bmal1, clock, cry2, csnkIepsilon, and npas2, whereas the strain with the least robust response to SD (DBA/2) exhibited either no change or a decrease in expression of these genes and cry1. The effect of SD on circadian clock gene expression was maintained in mice in which both of the cryptochrome genes were genetically inactivated. cry1 and cry2 appear to be redundant in sleep regulation as elimination of either of these genes did not result in a significant deficit in sleep homeostasis. These data demonstrate transcriptional regulatory correlates to previously described strain differences at the EEG level and raise the possibility that genetic differences underlying circadian clock gene expression may drive the EEG differences among these strains.

Renin and angiotensin II receptor gene expression in kidneys of renal hypertensive rats.

The aim of this investigation was to examine the interrelation between renal mRNA levels of renin and angiotensin II receptor type 1 (AT1) in a renin-dependent form of experimental hypertension. Rats were studied 4 weeks after unilateral renal artery clipping. Mean blood pressure and plasma renin activity were significantly higher in the hypertensive rats (n = 10 206 +/- mm Hg and 72.4 +/- 20.9 ng/mL-1/h-1, respectively) than in sham-operated controls (n = 10, 136 +/- 3 mm Hg and 3.3 +/- 0.5 ng/mL-1/h, respectively). Northern blot analysis of polyA+ RNA obtained from the kidneys of renal hypertensive rats showed increased levels of renin mRNA in the clipped kidney, whereas a decrease was observed in the unclipped kidney. Plasma renin activity was directly correlated with renin mRNA expression of the poststenotic kidney (r = .94, P < .01). AT1 mRNA expression was lower in both kidneys of the hypertensive rats. This downregulation was specific for the AT1A subtype since the renal expression of the AT1B subtype remained normal in hypertensive rats. The downregulation of the renal AT1A receptor may be due to high circulating angiotensin II levels. This is supported by the significant inverse correlation (r = .71, P < .01) between plasma renin activity and AT1A mRNA expression measured in the clipped kidney of the hypertensive rats.

Characterization of human gene expression changes after olive oil ingestion: an exploratory approach

Olive oil consumption is protective against risk factors for cardiovascular and cancer diseases. A nutrigenomic approach was performed to assess whether changes in gene expression could occur in human peripheral blood mononuclear cells after oli ve oil ingestion at postprandial state. Six healthy male volunteers ingested, at fasting state, 50 ml of olive oil. Prior to intervention a 1-week washout period with a controlled diet and sunflower oil as the only source of fat was followed. During the 3 days before and on the intervention day, a very low-phenolic compound diet was followed. At baseline (0 h) and at post-ingestion (6 h), total RNA was isolated and gene expression (29,082 genes) was evaluated by microarray. From microarray data, nutrient-gene interactions were observed in genes related to metabolism, cellular processes, cancer, and atherosclerosis (e.g. USP48 by 2.16; OGT by 1.68-fold change) and associated processes such as inflammation (e.g. AKAP13 by 2.30; IL-10 by 1.66-fold change) and DNA damage (e.g. DCLRE1C by 1.47; POLK by 1.44- fold change). When results obtained by microarray were verified by qRT-PCR in nine genes, full concordance was achieved only in the case of up-regulated genes. Changes were observed at a real-life dose of olive oil, as it is daily consumed in some Mediterranean areas. Our results support the hypothesis that postprandial protective changes related to olive oil consumption could be mediated through gene expression changes.

Increased endometrial placenta growth factor (PLGF) gene expression in women with successful implantation.

Objective: To analyze the vascularization of the endometrium via hysteroscopy and to assess its correlation with angiogenic factor gene expression and embryo implantation rate.Design: Cross-sectional study.Setting: Public university hospital.Patient(s): Patients undergoing hysteroscopy for supposed infertility.Intervention(s): Endometrial quality assessment according to Sakumoto-Masamoto, performed in the early secretory phase of the cycle. Collection of an endometrial tissue biopsy.Main Outcome Measure(s): RNA extraction, reverse transcription, and determination of gene expression of angiogenesis- and implantation-relevant factors using quantitative polymerase chain reaction. Retrieval of pregnancy information from the medical records.Result(s): Good quantity/quality RNA with infertility history was obtained from 63 participating women. Those with a "good" endometrium and subsequent pregnancy showed increased gene expression for placenta growth factor when compared with patients with a "bad" endometrium and who did not succeed with pregnancy to date. Nonpregnant women with a "good" endometrium presented an intermediate result. No significant differences were observed for several other genes tested, but trends in the same direction were observed.Conclusion(s): This study demonstrates for the first time that endometrial PLGF expression corresponds to the hysteroscopic appearance of the endometrium, and therefore has potential as a clinically relevant prognosticator for infertility treatment success. (Fertil Steril (R) 2011;96:663-8. (C)2011 by American Society for Reproductive Medicine.)

Differential Gene Expression between Fall- and Spring-Run Chinook Salmon Assessed by Long Serial Analysis of Gene Expression

Of all Pacific salmonids, Chinook salmon Oncorhynchus tshawytscha display the greatest variability in return times to freshwater. The molecular mechanisms of these differential return times have not been well described. Current methods, such as long serial analysis of gene expression (LongSAGE) and microarrays, allow gene expression to be analyzed for thousands of genes simultaneously. To investigate whether differential gene expression is observed between fall- and spring-run Chinook salmon from California's Central Valley, LongSAGE libraries were constructed. Three libraries containing between 25,512 and 29,372 sequenced tags (21 base pairs/tag) were generated using messenger RNA from the brains of adult Chinook salmon returning in fall and spring and from one ocean-caught Chinook salmon. Tags were annotated to genes using complementary DNA libraries from Atlantic salmon Salmo salar and rainbow trout O. mykiss. Differentially expressed genes, as estimated by differences in the number of sequence tags, were found in all pairwise comparisons of libraries (freshwater versus saltwater = 40 genes; fall versus spring = 11 genes: and spawning versus nonspawning = 51 genes). The gene for ependymin, an extracellular glycoprotein involved in behavioral plasticity in fish, exhibited the most differential expression among the three groupings. Reverse transcription polymerase chain reaction analysis verified the differential expression of ependymin between the fall- and spring-run samples. These LongSAGE libraries, the first reported for Chinook salmon, provide a window of the transcriptional changes during Chinook salmon return migration to freshwater and spawning and increase the amount of expressed sequence data.

Genetically programmed autoinducer destruction reduces virulence gene expression and swarming motility in Pseudomonas aeruginosa PAO1.

Virulence in the opportunistic human pathogen Pseudomonas aeruginosa is controlled by cell density via diffusible signalling molecules ('autoinducers') of the N-acylhomoserine lactone (AHL) type. Two Bacillus sp. isolates (A23 and A24) with AHL-degrading activity were identified among a large collection of rhizosphere bacteria. From isolate A24 a gene was cloned which was similar to the aiiA gene, encoding an AHL lactonase in another Bacillus strain. Expression of the aiiA homologue from isolate A24 in P. aeruginosa PAO1 reduced the amount of the quorum sensing signal N-oxododecanoyl-L-homoserine lactone and completely prevented the accumulation of the second AHL signal, N-butyryl-L-homoserine lactone. This strongly reduced AHL content correlated with a markedly decreased expression and production of several virulence factors and cytotoxic compounds such as elastase, rhamnolipids, hydrogen cyanide and pyocyanin, and strongly reduced swarming. However, no effect was observed on flagellar swimming or on twitching motility, and aiiA expression did not affect bacterial adhesion to a polyvinylchloride surface. In conclusion, introduction of an AHL degradation gene into P. aeruginosa could block cell-cell communication and exoproduct formation, but failed to interfere with surface colonization.