Impaired expression of NADH dehydrogenase subunit 1 and PPARgamma coactivator-1 in skeletal muscle of ZDF rats: restoration by troglitazone.

Type 2 diabetes has been related to a decrease of mitochondrial DNA (mtDNA) content. In this study, we show increased expression of the peroxisome proliferator-activated receptor-alpha (PPARalpha) and its target genes involved in fatty acid metabolism in skeletal muscle of Zucker Diabetic Fatty (ZDF) (fa/fa) rats. In contrast, the mRNA levels of genes involved in glucose transport and utilization (GLUT4 and phosphofructokinase) were decreased, whereas the expression of pyruvate dehydrogenase kinase 4 (PDK-4), which suppresses glucose oxidation, was increased. The shift from glucose to fatty acids as the source of energy in skeletal muscle of ZDF rats was accompanied by a reduction of subunit 1 of complex I (NADH dehydrogenase subunit 1, ND1) and subunit II of complex IV (cytochrome c oxidase II, COII), two genes of the electronic transport chain encoded by mtDNA. The transcript levels of PPARgamma Coactivator 1 (PGC-1) showed a significant reduction. Treatment with troglitazone (30 mg/kg/day) for 15 days reduced insulin values and reversed the increase in PDK-4 mRNA levels, suggesting improved insulin sensitivity. In addition, troglitazone treatment restored ND1 and PGC-1 expression in skeletal muscle. These results suggest that troglitazone may avoid mitochondrial metabolic derangement during the development of diabetes mellitus 2 in skeletal muscle.

Comparative analysis of human and mouse expression data illuminates tissue-specific evolutionary patterns of miRNAs.

MicroRNAs (miRNAs) constitute an important class of gene regulators. While models have been proposed to explain their appearance and expansion, the validation of these models has been difficult due to the lack of comparative studies. Here, we analyze miRNA evolutionary patterns in two mammals, human and mouse, in relation to the age of miRNA families. In this comparative framework, we confirm some predictions of previously advanced models of miRNA evolution, e.g. that miRNAs arise more frequently de novo than by duplication, or that the number of protein-coding gene targeted by miRNAs decreases with evolutionary time. We also corroborate that miRNAs display an increase in expression level with evolutionary time, however we show that this relation is largely tissue-dependent, and especially low in embryonic or nervous tissues. We identify a bias of tag-sequencing techniques regarding the assessment of breadth of expression, leading us, contrary to predictions, to find more tissue-specific expression of older miRNAs. Together, our results refine the models used so far to depict the evolution of miRNA genes. They underline the role of tissue-specific selective forces on the evolution of miRNAs, as well as the potential co-evolution patterns between miRNAs and the protein-coding genes they target.

Vitellogenin Underwent Subfunctionalization to Acquire Caste and Behavioral Specific Expression in the Harvester Ant Pogonomyrmex barbatus.

The reproductive ground plan hypothesis (RGPH) proposes that the physiological pathways regulating reproduction were co-opted to regulate worker division of labor. Support for this hypothesis in honeybees is provided by studies demonstrating that the reproductive potential of workers, assessed by the levels of vitellogenin (Vg), is linked to task performance. Interestingly, contrary to honeybees that have a single Vg ortholog and potentially fertile nurses, the genome of the harvester ant Pogonomyrmex barbatus harbors two Vg genes (Pb_Vg1 and Pb_Vg2) and nurses produce infertile trophic eggs. P. barbatus, thus, provides a unique model to investigate whether Vg duplication in ants was followed by subfunctionalization to acquire reproductive and non-reproductive functions and whether Vg reproductive function was co-opted to regulate behavior in sterile workers. To investigate these questions, we compared the expression patterns of P. barbatus Vg genes and analyzed the phylogenetic relationships and molecular evolution of Vg genes in ants. qRT-PCRs revealed that Pb_Vg1 is more highly expressed in queens compared to workers and in nurses compared to foragers. By contrast, the level of expression of Pb_Vg2 was higher in foragers than in nurses and queens. Phylogenetic analyses show that a first duplication of the ancestral Vg gene occurred after the divergence between the poneroid and formicoid clades and subsequent duplications occurred in the lineages leading to Solenopsis invicta, Linepithema humile and Acromyrmex echinatior. The initial duplication resulted in two Vg gene subfamilies preferentially expressed in queens and nurses (subfamily A) or in foraging workers (subfamily B). Finally, molecular evolution analyses show that the subfamily A experienced positive selection, while the subfamily B showed overall relaxation of purifying selection. Our results suggest that in P. barbatus the Vg gene underwent subfunctionalization after duplication to acquire caste- and behavior- specific expression associated with reproductive and non-reproductive functions, supporting the validity of the RGPH in ants.

CB1 blockade potentiates down-regulation of lipogenic gene expression in perirenal adipose tissue in high carbohydrate diet-induced obesity.

De novo lipogenesis and hypercaloric diets are thought to contribute to increased fat mass, particularly in abdominal fat depots. CB1 is highly expressed in adipose tissue, and CB1-mediated signalling is associated with stimulation of lipogenesis and diet-induced obesity, though its contribution to increasing fat deposition in adipose tissue is controversial. Lipogenesis is regulated by transcription factors such as liver X receptor (LXR), sterol-response element binding protein (SREBP) and carbohydrate-responsive-element-binding protein (ChREBP). We evaluated the role of CB1 in the gene expression of these factors and their target genes in relation to lipogenesis in the perirenal adipose tissue (PrAT) of rats fed a high-carbohydrate diet (HCHD) or a high-fat diet (HFD). Both obesity models showed an up-regulated gene expression of CB1 and Lxrα in this adipose pad. The Srebf-1 and ChREBP gene expressions were down-regulated in HFD but not in HCHD. The expression of their target genes encoding for lipogenic enzymes showed a decrease in diet-induced obesity and was particularly dramatic in HFD. In HCHD, CB1 blockade by AM251 reduced the Srebf-1 and ChREBP expression and totally abrogated the remnant gene expression of their target lipogenic enzymes. The phosphorylated form of the extracellular signal-regulated kinase (ERK-p), which participates in the CB1-mediated signalling pathway, was markedly present in the PrAT of obese rats. ERK-p was drastically repressed by AM251 indicating that CB1 is actually functional in PrAT of obese animals, though its activation loses the ability to stimulate lipogenesis in PrAT of obese rats. Even so, the remnant expression levels of lipogenic transcription factors found in HCHD-fed rats are still dependent on CB1 activity. Hence, in HCHD-induced obesity, CB1 blockade may help to further potentiate the reduction of lipogenesis in PrAT by means of inducing down-regulation of the ChREBP and Srebf-1 gene expression, and consequently in the expression of lipogenic enzymes.

Response to Infliximab in Crohn's Disease: Genetic Analysis Supporting Expression Profile.

Substantial proportion of Crohn's disease (CD) patients shows no response or a limited response to treatment with infliximab (IFX) and to identify biomarkers of response would be of great clinical and economic benefit. The expression profile of five genes (S100A8-S100A9, G0S2, TNFAIP6, and IL11) reportedly predicted response to IFX and we aimed at investigating their etiologic role through genetic association analysis. Patients with active CD (350) who received at least three induction doses of IFX were included and classified according to IFX response. A tagging strategy was used to select genetic polymorphisms that cover the variability present in the chromosomal regions encoding the identified genes with altered expression. Following genotyping, differences between responders and nonresponders to IFX were observed in haplotypes of the studied regions: S100A8-S100A9 (rs11205276* G/rs3014866* C/rs724781* C/rs3006488* A; P = 0.05); G0S2 (rs4844486* A/rs1473683* T; P = 0.15); TNFAIP6 (rs11677200* C/rs2342910* A/rs3755480* G/rs10432475* A; P = 0.10); and IL11 (rs1126760* C/rs1042506* G; P = 0.07). These differences were amplified in patients with colonic and ileocolonic location for all but the TNFAIP6 haplotype, which evidenced significant difference in ileal CD patients. Our results support the role of the reported expression signature as predictive of anti-TNF outcome in CD patients and suggest an etiological role of those top-five genes in the IFX response pathway.

Two waves of recombinase gene expression in developing thymocytes.

During T cell development in the thymus, T cell receptor (TCR) alpha, beta, gamma, and delta genes are rearranged and expressed. TCR rearrangement strictly depends upon the coordinate activity of two recombinase activating genes, Rag-1 and Rag-2. In this study we have followed the expression of these genes at different stages of intrathymic development. The results indicate that there are two periods of high Rag-1 and Rag-2 mRNA expression. The first wave peaks early at the CD25+CD4-CD8-CD3- stage of development and coincides with the initial appearance of transcripts derived from fully rearranged TCR beta, gamma, and delta genes, whereas the second wave occurs later at the CD4+CD8+ stage coincident with full-length TCR alpha mRNA expression. Active downregulation of Rag-1 and Rag-2 mRNA expression appears to occur in vivo between the two peaks of recombinase activity. This phenomenon can be mimicked in vitro in response to artificial stimuli such as phorbol myristate acetate and calcium ionophore. Collectively our data suggest that recombinase expression is actively regulated during early thymus development independently of cell surface expression of a mature heterodimeric TCR protein complex.

Mitofusins 1/2 and ERRalpha expression are increased in human skeletal muscle after physical exercise

Mitochondrial impairment is hypothesized to contribute to the pathogenesis of insulin resistance. Mitofusin (Mfn) proteins regulate the biogenesis and maintenance of the mitochondrial network, and when inactivated, cause a failure in the mitochondrial architecture and decreases in oxidative capacity and glucose oxidation. Exercise increases muscle mitochondrial content, size, oxidative capacity and aerobic glucose oxidation. To address if Mfn proteins are implicated in these exercise-induced responses, we measured Mfn1 and Mfn2 mRNA levels, pre-, post-, 2 and 24 h post-exercise. Additionally, we measured the expression levels of transcriptional regulators that control mitochondrial biogenesis and functions, including PGC-1alpha, NRF-1, NRF-2 and the recently implicated ERRalpha. We show that Mfn1, Mfn2, NRF-2 and COX IV mRNA were increased 24 h post-exercise, while PGC-1alpha and ERRalpha mRNA increased 2 h post-exercise. Finally, using in vitro cellular assays, we demonstrate that Mfn2 gene expression is driven by a PGC-1alpha programme dependent on ERRalpha. The PGC-1alpha/ERRalpha-mediated induction of Mfn2 suggests a role of these two factors in mitochondrial fusion. Our results provide evidence that PGC-1alpha not only mediates the increased expression of oxidative phosphorylation genes but also mediates alterations in mitochondrial architecture in response to aerobic exercise in humans

Evolution of the correlation between expression divergence and protein divergence in mammals.

Divergence of protein sequences and gene expression patterns are two fundamental mechanisms that generate organismal diversity. Here, we have used genome and transcriptome data from eight mammals and one bird to study the positive correlation of these two processes throughout mammalian evolution. We demonstrate that the correlation is stable over time and most pronounced in neural tissues, which indicates that it is the result of strong negative selection. The correlation is not driven by genes with specific functions and may instead best be viewed as an evolutionary default state, which can nevertheless be evaded by certain gene types. In particular, genes with developmental and neural functions are skewed toward changes in gene expression, consistent with selection against pleiotropic effects associated with changes in protein sequences. Surprisingly, we find that the correlation between expression divergence and protein divergence is not explained by between-gene variation in expression level, tissue specificity, protein connectivity, or other investigated gene characteristics, suggesting that it arises independently of these gene traits. The selective constraints on protein sequences and gene expression patterns also fluctuate in a coordinate manner across phylogenetic branches: We find that gene-specific changes in the rate of protein evolution in a specific mammalian lineage tend to be accompanied by similar changes in the rate of expression evolution. Taken together, our findings highlight many new aspects of the correlation between protein divergence and expression divergence, and attest to its role as a fundamental property of mammalian genome evolution.

Integrative analysis of RUNX1 downstream pathways and target genes.

BACKGROUND: The RUNX1 transcription factor gene is frequently mutated in sporadic myeloid and lymphoid leukemia through translocation, point mutation or amplification. It is also responsible for a familial platelet disorder with predisposition to acute myeloid leukemia (FPD-AML). The disruption of the largely unknown biological pathways controlled by RUNX1 is likely to be responsible for the development of leukemia. We have used multiple microarray platforms and bioinformatic techniques to help identify these biological pathways to aid in the understanding of why RUNX1 mutations lead to leukemia. RESULTS: Here we report genes regulated either directly or indirectly by RUNX1 based on the study of gene expression profiles generated from 3 different human and mouse platforms. The platforms used were global gene expression profiling of: 1) cell lines with RUNX1 mutations from FPD-AML patients, 2) over-expression of RUNX1 and CBFbeta, and 3) Runx1 knockout mouse embryos using either cDNA or Affymetrix microarrays. We observe that our datasets (lists of differentially expressed genes) significantly correlate with published microarray data from sporadic AML patients with mutations in either RUNX1 or its cofactor, CBFbeta. A number of biological processes were identified among the differentially expressed genes and functional assays suggest that heterozygous RUNX1 point mutations in patients with FPD-AML impair cell proliferation, microtubule dynamics and possibly genetic stability. In addition, analysis of the regulatory regions of the differentially expressed genes has for the first time systematically identified numerous potential novel RUNX1 target genes. CONCLUSION: This work is the first large-scale study attempting to identify the genetic networks regulated by RUNX1, a master regulator in the development of the hematopoietic system and leukemia. The biological pathways and target genes controlled by RUNX1 will have considerable importance in disease progression in both familial and sporadic leukemia as well as therapeutic implications

Inducible gene and shRNA expression in resident hematopoietic stem cells in vivo.

Hematopoietic stem cells (HSC) are probably the best understood somatic stem cells and often serve as a paradigm for other stem cells. Nevertheless, most current techniques to genetically manipulate them in vivo are either constitutive and/or induced in settings of hematopoietic stress such as after irradiation. Here, we present a conditional expression system that allows for externally controllable transgenesis and knockdown in resident HSCs, based on a lentiviral vector containing a tet-O sequence and a transgenic mouse line expressing a doxycyclin-regulated tTR-KRAB repressor protein. HSCs harvested from tTR-KRAB mice are transduced with the lentiviral vector containing a cDNA (i.e., Green Fluorescent Protein (GFP)) and/or shRNA (i.e., p53) of interest and then transplanted into lethally irradiated recipients. While the vector is effectively repressed by tTR-KRAB during homing and engraftment, robust GFP/shp53 expression is induced on doxycyclin treatment in HSCs and their progeny. Doxycylin-controllable transcription is maintained on serial transplantation, indicating that repopulating HSCs are stably modified by this approach. In summary, this easy to implement conditional system provides inducible and reversible overexpression or knock down of genes in resident HSCs in vivo using a drug devoid of toxic or activating effects.

Functional diversification of duplicate genes through subcellular adaptation of encoded proteins.

BACKGROUND: Gene duplication is the primary source of new genes with novel or altered functions. It is known that duplicates may obtain these new functional roles by evolving divergent expression patterns and/or protein functions after the duplication event. Here, using yeast (Saccharomyces cerevisiae) as a model organism, we investigate a previously little considered mode for the functional diversification of duplicate genes: subcellular adaptation of encoded proteins. RESULTS: We show that for 24-37% of duplicate gene pairs derived from the S. cerevisiae whole-genome duplication event, the two members of the pair encode proteins that localize to distinct subcellular compartments. The propensity of yeast duplicate genes to evolve new localization patterns depends to a large extent on the biological function of their progenitor genes. Proteins involved in processes with a wider subcellular distribution (for example, catabolism) frequently evolved new protein localization patterns after duplication, whereas duplicate proteins limited to a smaller number of organelles (for example, highly expressed biosynthesis/housekeeping proteins with a slow rate of evolution) rarely relocate within the cell. Paralogous proteins evolved divergent localization patterns by partitioning of ancestral localizations ('sublocalization'), but probably more frequently by relocalization to new compartments ('neolocalization'). We show that such subcellular reprogramming may occur through selectively driven substitutions in protein targeting sequences. Notably, our data also reveal that relocated proteins functionally adapted to their new subcellular environments and evolved new functional roles through changes of their physico-chemical properties, expression levels, and interaction partners. CONCLUSION: We conclude that protein subcellular adaptation represents a common mechanism for the functional diversification of duplicate genes.

Comprehensive analysis of Arabidopsis expression level polymorphisms with simple inheritance.

In Arabidopsis thaliana, gene expression level polymorphisms (ELPs) between natural accessions that exhibit simple, single locus inheritance are promising quantitative trait locus (QTL) candidates to explain phenotypic variability. It is assumed that such ELPs overwhelmingly represent regulatory element polymorphisms. However, comprehensive genome-wide analyses linking expression level, regulatory sequence and gene structure variation are missing, preventing definite verification of this assumption. Here, we analyzed ELPs observed between the Eil-0 and Lc-0 accessions. Compared with non-variable controls, 5' regulatory sequence variation in the corresponding genes is indeed increased. However, approximately 42% of all the ELP genes also carry major transcription unit deletions in one parent as revealed by genome tiling arrays, representing a >4-fold enrichment over controls. Within the subset of ELPs with simple inheritance, this proportion is even higher and deletions are generally more severe. Similar results were obtained from analyses of the Bay-0 and Sha accessions, using alternative technical approaches. Collectively, our results suggest that drastic structural changes are a major cause for ELPs with simple inheritance, corroborating experimentally observed indel preponderance in cloned Arabidopsis QTL.

The JNK binding domain of islet-brain 1 inhibits IL-1 induced JNK activity and apoptosis but not the transcription of key proapoptotic or protective genes in insulin-secreting cell lines.

The stress-activated protein kinase c-Jun NH2-terminal kinase (JNK) is a central signal for interleukin-1beta (IL-1beta)-induced apoptosis in insulin-producing beta-cells. The cell-permeable peptide inhibitor of JNK (JNKI1), that introduces the JNK binding domain (JBD) of the scaffold protein islet-brain 1 (IB1) inside cells, effectively prevents beta-cell death caused by this cytokine. To define the molecular targets of JNK involved in cytokine-induced beta-cell apoptosis we investigated whether JNKI1 or stable expression of JBD affected the expression of selected pro- and anti-apoptotic genes induced in rat (RIN-5AH-T2B) and mouse (betaTC3) insulinoma cells exposed to IL-1beta. Inhibition of JNK significantly reduced phosphorylation of the specific JNK substrate c-Jun (p<0.05), IL-1beta-induced apoptosis (p<0.001), and IL-1beta-mediated c-fos gene expression. However, neither JNKI1 nor JBD did influence IL-1beta-induced NO synthesis or iNOS expression or the transcription of the genes encoding mitochondrial manganese superoxide dismutase (MnSOD), catalase (CAT), glutathione peroxidase (GPx), glutathione-S-transferase rho (GSTrho), heat shock protein (HSP) 70, IL-1beta-converting enzyme (ICE), caspase-3, apoptosis-inducing factor (AIF), Bcl-2 or Bcl-xL. We suggest that the anti-apoptotic effect of JNK inhibition by JBD is independent of the transcription of major pro- and anti-apoptotic genes, but may be exerted at the translational or posttranslational level.

Functional analysis: evaluation of response intensities--tailoring ANOVA for lists of expression subsets.

Background: Microarray data is frequently used to characterize the expression profile of a whole genome and to compare the characteristics of that genome under several conditions. Geneset analysis methods have been described previously to analyze the expression values of several genes related by known biological criteria (metabolic pathway, pathology signature, co-regulation by a common factor, etc.) at the same time and the cost of these methods allows for the use of more values to help discover the underlying biological mechanisms. Results: As several methods assume different null hypotheses, we propose to reformulate the main question that biologists seek to answer. To determine which genesets are associated with expression values that differ between two experiments, we focused on three ad hoc criteria: expression levels, the direction of individual gene expression changes (up or down regulation), and correlations between genes. We introduce the FAERI methodology, tailored from a two-way ANOVA to examine these criteria. The significance of the results was evaluated according to the self-contained null hypothesis, using label sampling or by inferring the null distribution from normally distributed random data. Evaluations performed on simulated data revealed that FAERI outperforms currently available methods for each type of set tested. We then applied the FAERI method to analyze three real-world datasets on hypoxia response. FAERI was able to detect more genesets than other methodologies, and the genesets selected were coherent with current knowledge of cellular response to hypoxia. Moreover, the genesets selected by FAERI were confirmed when the analysis was repeated on two additional related datasets. Conclusions: The expression values of genesets are associated with several biological effects. The underlying mathematical structure of the genesets allows for analysis of data from several genes at the same time. Focusing on expression levels, the direction of the expression changes, and correlations, we showed that two-step data reduction allowed us to significantly improve the performance of geneset analysis using a modified two-way ANOVA procedure, and to detect genesets that current methods fail to detect.