Evolution of gene expression in fire ants: the effects of developmental stage, caste, and species.

Ants provide remarkable examples of equivalent genotypes developing into divergent and discrete phenotypes. Diploid eggs can develop either into queens, which specialize in reproduction, or workers, which participate in cooperative tasks such as building the nest, collecting food, and rearing the young. In contrast, the differentiation between males and females generally depends upon whether eggs are fertilized, with fertilized (diploid) eggs giving rise to females and unfertilized (haploid) eggs giving rise to males. To obtain a comprehensive picture of the relative contributions of gender (sex), caste, developmental stage, and species divergence to gene expression evolution, we investigated gene expression patterns in pupal and adult queens, workers, and males of two species of fire ants, Solenopsis invicta and S. richteri. Microarray hybridizations revealed that variation in gene expression profiles is influenced more by developmental stage than by caste membership, sex, or species identity. The second major contributor to variation in gene expression was the combination of sex and caste. Although workers and queens share equivalent diploid nuclear genomes, they have highly distinctive patterns of gene expression in both the pupal and the adult stages, as might be expected given their extraordinary level of phenotypic differentiation. Overall, the difference in the proportion of differentially expressed genes was greater between workers and males than between workers and queens or queens and males, consistent with the fact that workers and males share neither gender nor reproductive capability. Moreover, between-species comparisons revealed that the greatest difference in gene expression patterns occurred in adult workers, a finding consistent with the fact that adult workers most directly experience the distinct external environments characterizing the different habitats occupied by the two species. Thus, much of the evolution of gene expression in ants may occur in the worker caste, despite the fact that these individuals are largely or completely sterile. Analyses of gene expression evolution revealed a combination of positive selection and relaxation of stabilizing selection as important factors driving the evolution of such genes.

Differentiation associated regulation of microRNA expression in vivo in human CD8+ T cell subsets.

BACKGROUND: The differentiation of CD8+ T lymphocytes following priming of naïve cells is central in the establishment of the adaptive immune response. Yet, the molecular events underlying this process are not fully understood. MicroRNAs have been recently shown to play a key role in the regulation of haematopoiesis in mouse, but their implication in peripheral lymphocyte differentiation in humans remains largely unknown. METHODS: In order to explore the potential implication of microRNAs in CD8+ T cell differentiation in humans, microRNA expression profiles were analysed using microarrays and quantitative PCR in several human CD8+ T cell subsets defining the major steps of the T cell differentiation pathway. RESULTS: We found expression of a limited set of microRNAs, including the miR-17~92 cluster. Moreover, we reveal the existence of differentiation-associated regulation of specific microRNAs. When compared to naive cells, miR-21 and miR-155 were indeed found upregulated upon differentiation to effector cells, while expression of the miR-17~92 cluster tended to concomitantly decrease. CONCLUSIONS: This study establishes for the first time in a large panel of individuals the existence of differentiation associated regulation of microRNA expression in human CD8+ T lymphocytes in vivo, which is likely to impact on specific cellular functions.

Gene expression patterns unveil a new level of molecular heterogeneity in colorectal cancer.

The recognition that colorectal cancer (CRC) is a heterogeneous disease in terms of clinical behaviour and response to therapy translates into an urgent need for robust molecular disease subclassifiers that can explain this heterogeneity beyond current parameters (MSI, KRAS, BRAF). Attempts to fill this gap are emerging. The Cancer Genome Atlas (TGCA) reported two main CRC groups, based on the incidence and spectrum of mutated genes, and another paper reported an EMT expression signature defined subgroup. We performed a prior free analysis of CRC heterogeneity on 1113 CRC gene expression profiles and confronted our findings to established molecular determinants and clinical, histopathological and survival data. Unsupervised clustering based on gene modules allowed us to distinguish at least five different gene expression CRC subtypes, which we call surface crypt-like, lower crypt-like, CIMP-H-like, mesenchymal and mixed. A gene set enrichment analysis combined with literature search of gene module members identified distinct biological motifs in different subtypes. The subtypes, which were not derived based on outcome, nonetheless showed differences in prognosis. Known gene copy number variations and mutations in key cancer-associated genes differed between subtypes, but the subtypes provided molecular information beyond that contained in these variables. Morphological features significantly differed between subtypes. The objective existence of the subtypes and their clinical and molecular characteristics were validated in an independent set of 720 CRC expression profiles. Our subtypes provide a novel perspective on the heterogeneity of CRC. The proposed subtypes should be further explored retrospectively on existing clinical trial datasets and, when sufficiently robust, be prospectively assessed for clinical relevance in terms of prognosis and treatment response predictive capacity. Original microarray data were uploaded to the ArrayExpress database (http://www.ebi.ac.uk/arrayexpress/) under Accession Nos E-MTAB-990 and E-MTAB-1026. © 2013 Swiss Institute of Bioinformatics. Journal of Pathology published by John Wiley & Sons Ltd on behalf of Pathological Society of Great Britain and Ireland.

Gene expression profiling reveals consistent differences between clinical samples of human leukaemias and their model cell lines.

Microarray gene expression profiles of fresh clinical samples of chronic myeloid leukaemia in chronic phase, acute promyelocytic leukaemia and acute monocytic leukaemia were compared with profiles from cell lines representing the corresponding types of leukaemia (K562, NB4, HL60). In a hierarchical clustering analysis, all clinical samples clustered separately from the cell lines, regardless of leukaemic subtype. Gene ontology analysis showed that cell lines chiefly overexpressed genes related to macromolecular metabolism, whereas in clinical samples genes related to the immune response were abundantly expressed. These findings must be taken into consideration when conclusions from cell line-based studies are extrapolated to patients.

Correcting for the bias due to expression specificity improves the estimation of constrained evolution of expression between mouse and human.

MOTIVATION: Comparative analyses of gene expression data from different species have become an important component of the study of molecular evolution. Thus methods are needed to estimate evolutionary distances between expression profiles, as well as a neutral reference to estimate selective pressure. Divergence between expression profiles of homologous genes is often calculated with Pearson's or Euclidean distance. Neutral divergence is usually inferred from randomized data. Despite being widely used, neither of these two steps has been well studied. Here, we analyze these methods formally and on real data, highlight their limitations and propose improvements. RESULTS: It has been demonstrated that Pearson's distance, in contrast to Euclidean distance, leads to underestimation of the expression similarity between homologous genes with a conserved uniform pattern of expression. Here, we first extend this study to genes with conserved, but specific pattern of expression. Surprisingly, we find that both Pearson's and Euclidean distances used as a measure of expression similarity between genes depend on the expression specificity of those genes. We also show that the Euclidean distance depends strongly on data normalization. Next, we show that the randomization procedure that is widely used to estimate the rate of neutral evolution is biased when broadly expressed genes are abundant in the data. To overcome this problem, we propose a novel randomization procedure that is unbiased with respect to expression profiles present in the datasets. Applying our method to the mouse and human gene expression data suggests significant gene expression conservation between these species. CONTACT: marc.robinson-rechavi@unil.ch; sven.bergmann@unil.ch SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Gene Expression Profiling of Desmoid Tumors by cDNA Microarrays and Correlation with Progression-Free Survival.

PURPOSE: Because desmoid tumors exhibit an unpredictable clinical course, translational research is crucial to identify the predictive factors of progression in addition to the clinical parameters. The main issue is to detect patients who are at a higher risk of progression. The aim of this work was to identify molecular markers that can predict progression-free survival (PFS). EXPERIMENTAL DESIGN: Gene-expression screening was conducted on 115 available independent untreated primary desmoid tumors using cDNA microarray. We established a prognostic gene-expression signature composed of 36 genes. To test robustness, we randomly generated 1,000 36-gene signatures and compared their outcome association to our define 36-genes molecular signature and we calculated positive predictive value (PPV) and negative predictive value (NPV). RESULTS: Multivariate analysis showed that our molecular signature had a significant impact on PFS while no clinical factor had any prognostic value. Among the 1,000 random signatures generated, 56.7% were significant and none was more significant than our 36-gene molecular signature. PPV and NPV were high (75.58% and 81.82%, respectively). Finally, the top two genes downregulated in no-recurrence were FECH and STOML2 and the top gene upregulated in no-recurrence was TRIP6. CONCLUSIONS: By analyzing expression profiles, we have identified a gene-expression signature that is able to predict PFS. This tool may be useful for prospective clinical studies. Clin Cancer Res; 21(18); 4194-200. ©2015 AACR.

Genomic, Epigenomic and Transcriptomic Molecular Characterization of Splenic Marginal Zone Lymphoma Reveals Recurrent Abnormalities in miR-182

Splenic marginal zone lymphoma (SMZL) is a low grade B-cell non-Hodgkin's lymphoma. The molecular pathology of this entity remains poorly understood. To characterise this lymphoma at the molecular level, we performed an integrated analysis of 1) genome wide genetic copy number alterations 2) gene expression profiles and 3) epigenetic DNA methylation profiles.We have previously shown that SMZL is characterised by recurrent alterations of chromosomes 7q, 6q, 3q, 9q and 18; however, gene resolution oligonucleotide array comparative genomic hybridisation did not reveal evidence of cryptic amplification or deletion in these regions. The most frequently lost 7q32 region contains a cluster of miRNAs. qRT-PCR revealed that three of these (miR-182/96/183) show underexpression in SMZL, and miR-182 is somatically mutated in >20% of cases of SMZL, as well as in >20% of cases of follicular lymphoma, and between 5-15% of cases of chronic lymphocytic leukaemia, MALT-lymphoma and hairy cell leukaemia. We conclude that miR-182 is a strong candidate novel tumour suppressor miRNA in lymphoma.The overall gene expression signature of SMZL was found to be strongly distinct fromthose of other lymphomas. Functional analysis of gene expression data revealed SMZL to be characterised by abnormalities in B-cell receptor signalling (especially through the CD19/21-PI3K/AKT pathway) and apoptotic pathways. In addition, genes involved in the response to viral infection appeared upregulated. SMZL shows a unique epigenetic profile, but analysis of differentially methylated genes showed few with methylation related transcriptional deregulation, suggesting that DNA methylation abnormalities are not a critical component of the SMZL malignant phenotype.

Genetics of Parkinson disease and essential tremor.

Purpose of review: Elucidating the genetic background of Parkinson disease and essential tremor is crucial to understand the pathogenesis and improve diagnostic and therapeutic strategies. Recent findings: A number of approaches have been applied including familial and association studies, and studies of gene expression profiles to identify genes involved in susceptibility to Parkinson disease. These studies have nominated a number of candidate Parkinson disease genes and novel loci including Omi/HtrA2, GIGYF2, FGF20, PDXK, EIF4G1 and PARK16. A recent notable finding has been the confirmation for the role of heterozygous mutations in glucocerebrosidase (GBA) as risk factors for Parkinson disease. Finally, association studies have nominated genetic variation in the leucine-rich repeat and Ig containing 1 gene (LINGO1) as a risk for both Parkinson disease and essential tremor, providing the first genetic evidence of a link between the two conditions. Summary: Although undoubtedly genes remain to be identified, considerable progress has been achieved in the understanding of the genetic basis of Parkinson disease. This same effort is now required for essential tremor. The use of next-generation high-throughput sequencing and genotyping technologies will help pave the way for future insight leading to advances in diagnosis, prevention and cure.

In vivo conditional Pax4 overexpression in mature islet β-cells prevents stress-induced hyperglycemia in mice.

OBJECTIVE To establish the role of the transcription factor Pax4 in pancreatic islet expansion and survival in response to physiological stress and its impact on glucose metabolism, we generated transgenic mice conditionally and selectively overexpressing Pax4 or a diabetes-linked mutant variant (Pax4R129 W) in β-cells. RESEARCH DESIGN AND METHODS Glucose homeostasis and β-cell death and proliferation were assessed in Pax4- or Pax4R129 W-overexpressing transgenic animals challenged with or without streptozotocin. Isolated transgenic islets were also exposed to cytokines, and apoptosis was evaluated by DNA fragmentation or cytochrome C release. The expression profiles of proliferation and apoptotic genes and β-cell markers were studied by immunohistochemistry and quantitative RT-PCR. RESULTS Pax4 but not Pax4R129 W protected animals against streptozotocin-induced hyperglycemia and isolated islets from cytokine-mediated β-cell apoptosis. Cytochrome C release was abrogated in Pax4 islets treated with cytokines. Interleukin-1β transcript levels were suppressed in Pax4 islets, whereas they were increased along with NOS2 in Pax4R129 W islets. Bcl-2, Cdk4, and c-myc expression levels were increased in Pax4 islets while MafA, insulin, and GLUT2 transcript levels were suppressed in both animal models. Long-term Pax4 expression promoted proliferation of a Pdx1-positive cell subpopulation while impeding insulin secretion. Suppression of Pax4 rescued this defect with a concomitant increase in pancreatic insulin content. CONCLUSIONS Pax4 protects adult islets from stress-induced apoptosis by suppressing selective nuclear factor-κB target genes while increasing Bcl-2 levels. Furthermore, it promotes dedifferentiation and proliferation of β-cells through MafA repression, with a concomitant increase in Cdk4 and c-myc expression.

IQRray, a new method for Affymetrix microarray quality control, and the homologous organ conservation score, a new benchmark method for quality control metrics.

MOTIVATION: Microarray results accumulated in public repositories are widely reused in meta-analytical studies and secondary databases. The quality of the data obtained with this technology varies from experiment to experiment, and an efficient method for quality assessment is necessary to ensure their reliability. RESULTS: The lack of a good benchmark has hampered evaluation of existing methods for quality control. In this study, we propose a new independent quality metric that is based on evolutionary conservation of expression profiles. We show, using 11 large organ-specific datasets, that IQRray, a new quality metrics developed by us, exhibits the highest correlation with this reference metric, among 14 metrics tested. IQRray outperforms other methods in identification of poor quality arrays in datasets composed of arrays from many independent experiments. In contrast, the performance of methods designed for detecting outliers in a single experiment like Normalized Unscaled Standard Error and Relative Log Expression was low because of the inability of these methods to detect datasets containing only low-quality arrays and because the scores cannot be directly compared between experiments. AVAILABILITY AND IMPLEMENTATION: The R implementation of IQRray is available at: ftp://lausanne.isb-sib.ch/pub/databases/Bgee/general/IQRray.R. CONTACT: Marta.Rosikiewicz@unil.ch SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Aging of myelinating glial cells predominantly affects lipid metabolism and immune response pathways.

Both the central and the peripheral nervous systems are prone to multiple age-dependent neurological deficits, often attributed to still unknown alterations in the function of myelinating glia. To uncover the biological processes affected in glial cells by aging, we analyzed gene expression of the Schwann cell-rich mouse sciatic nerve at 17 time points throughout life, from day of birth until senescence. By combining these data with the gene expression data of myelin mouse mutants carrying deletions of either Pmp22, SCAP, or Lpin1, we found that the majority of age-related transcripts were also affected in myelin mutants (54.4%) and were regulated during PNS development (59.5%), indicating a high level of overlap in implicated molecular pathways. The expression profiles in aging copied the direction of transcriptional changes observed in neuropathy models; however, they had the opposite direction when compared with PNS development. The most significantly altered biological processes in aging involved the inflammatory/immune response and lipid metabolism. Interestingly, both these pathways were comparably changed in the aging optic nerve, suggesting that similar biological processes are affected in aging of glia-rich parts of the central and peripheral nervous systems. Our comprehensive comparison of gene expression in three distinct biological conditions including development, aging, and myelin disease thus revealed a previously unanticipated relationship among themselves and identified lipid metabolism and inflammatory/immune response pathways as potential therapeutical targets to prevent or delay so far incurable age-related and inherited forms of neuropathies.

Bone anabolic potential of soy isoflavones and plant extracts and genomic changes during differentiation of hPOB-tert osteoblast-like cells

Summary For the nutritional management of bone health and the prevention of osteoporosis it is important to identify nutrients that positively influence the bone remodeling process at the cellular level. Soy isoflavones show promising osteoprotective effects in animals and humans but their mechanism of action in bone cells is yet poorly understood. Firstly, soy tissue cultures were characterized as a new and optimized source of isoflavones. A large variability in the isoflavone content was observed and high-producing strains (46.3 mg/g dry wt isoflavones) were identified. In the Ishikawa cells bioassay, the estrogenicity of isoflavones was confirmed to be 1000 to 10000 less than 17Mestradiol and that of the malonyl forms was shown for the first time (EC50 of 350 nM and 1880 nM for malonylgenistin and malonyldaidzin, respectively). The estrogenic activity of soya tissue culture extracts correlated to their isoflavone content. Secondly, the effects of phytonutrients on BMP-2 gene expression and on the mevalonate synthesis pathway, as key mediators of bone formation, were investigated. Dietary achievable concentrations of genistein and daidzein (10vM), and statins (4xM) but not 17M estradiol (10nM), induced BMP-2 gene expression (by up to 3-fold) and inhibited the cholesterol biosynthetic pathway (by up to 50%) in the human osteoblastic cell line hP0B¬tert. In addition, several plant extracts (Cyperus rotundus, Lindera benzoin and Cnidium monnieri) induced BMP-2 gene expression but this induction was not restricted to the inhibition of the cholesterol synthesis pathway neither to the estrogenicity. Finally, the gene expression profiles during hP0B-tert differentiation induced by vitamin D and dexamethasone were analyzed with the Affymetrix human GeneChip. 1665 different genes and 98 ESTs were significantly regulated. The expression profiles of bone-related genes was largely in agreement with previously documented patterns, supporting the physiological relevance of the genomic results and the hP0B-tert cell line as a valid model of human osteoblast differentiation. The expression of alternative differentiation markers during the osteogenic treatment of hP0B-tert cells indicated that the adipocyte and myoblast differentiation pathways were repressed, confirming that these culture conditions allowed only osteoblast differentiation. The gene ontology analysis identified further sub-groups of genes that may be involved in the bone formation process. Aims of the thesis In order to define new strategies for the nutritional management of bone health and for the prevention of osteoporosis the major goal of the present work was to investigate the potential of phytonutrients to positively modulate the bone formation process at the cellular level and, in particular: 1.To select and optimise alternative plant sources containing high levels of isoflavones with estrogenic activity (Chapter 3). 2.To compare the effects of statins and phytonutrients on BMP-2 gene expression and on the mevalonate synthesis pathway and to select new plant extracts with a bone anabolic potential (Chapter 4). 3.To further characterize the new human periosteal cell line, hP0B-tert, as a bone- formation model, by elucidating its gene expression profile during differentiation induced by vitamin D and dexamethasone (Chapter 5).

Cellular origin of T-cell lymphomas.

In this issue of Blood, Iqbal et al, having compiled gene expression profiles from >300 peripheral T-cell lymphomas, expand previous findings on the diagnostic value of molecular signatures that correlate with different histological types of T-cell lymphomas. They report the discovery of 2 molecular subgroups of peripheral T-cell lymphomas, not otherwise specified (PTCL, NOS), characterized by high expression of either GATA-binding protein 3 (GATA-3) or t-box 21 (TBX21) transcription factors and corresponding target genes, with the GATA3 subgroup being associated with distinctly worse prognosis. In an independent study, Wang et al(2) also show that GATA3 expression in a subset of PTCL, NOS identifies a subgroup of patients with inferior survival.

Changes in the transcriptional profile of transporters in the intestine along the anterior-posterior and crypt-villus axes.

BACKGROUND: The purpose of this work was to characterize the expression of drug and nutrient carriers along the anterior-posterior and crypt-villus axes of the intestinal epithelium and to study the validity of utilizing whole gut tissue rather than purified epithelial cells to examine regional variations in gene expression. RESULTS: We have characterized the mRNA expression profiles of 76 % of all currently known transporters along the anterior-posterior axis of the gut. This is the first study to describe the expression profiles of the majority of all known transporters in the intestine. The expression profiles of transporters, as defined according to the Gene Ontology consortium, were measured in whole tissue of the murine duodenum, jejunum, ileum and colon using high-density microarrays. For nine transporters (Abca1, Abcc1, Abcc3, Abcg8, Slc10a2, Slc28a2, Slc2a1, Slc34a2 and Slc5a8), the mRNA profiles were further measured by RT-PCR in laser micro-dissected crypt and villus epithelial cells corresponding to the aforementioned intestinal regions. With respect to differentially regulated transporters, the colon had a distinct expression profile from small intestinal segments. The majority (59 % for p cutoff < or = 0.05) of transporter mRNA levels were constant across the intestinal sections studied. For the transporter subclass "carrier activity", which contains the majority of known carriers for biologically active compounds, a significant change (p < or = 0.05) along the anterior-posterior axis was observed. CONCLUSION: All nine transporters examined in laser-dissected material demonstrated good replication of the region-specific profiles revealed by microarray. Furthermore, we suggest that the distribution characteristics of Slc5a8 along the intestinal tract render it a suitable candidate carrier for monocarboxylate drugs in the posterior portion of the intestine. Our findings also predict that there is a significant difference in the absorption of carrier-mediated compounds in the different intestinal segments. The most pronounced differences can be expected between the adjoining segments ileum and colon, but the differences between the other adjoining segments are not negligible. Finally, for the examined genes, profiles measured in whole intestinal tissue extracts are representative of epithelial cell-only gene expression.

Species differences in the response of liver drug-metabolizing enzymes to (S)-4-O-tolylsulfanyl-2-(4-trifluormethyl-phenoxy)-butyric acid (EMD 392949) in vivo and in vitro.

Induction of drug-metabolizing enzymes (DMEs) is highly species-specific and can lead to drug-drug interaction and toxicities. In this series of studies we tested the species specificity of the antidiabetic drug development candidate and mixed peroxisome proliferator-activated receptor (PPAR) alpha/gamma agonist (S)-4-O-tolylsulfanyl-2-(4-trifluormethyl-phenoxy)-butyric acid (EMD 392949, EMD) with regard to the induction of gene expression and activities of DMEs, their regulators, and typical PPAR target genes. EMD clearly induced PPARalpha target genes in rats in vivo and in rat hepatocytes but lacked significant induction of DMEs, except for cytochrome P450 (P450) 4A. CYP2C and CYP3A were consistently induced in livers of EMD-treated monkeys. Interestingly, classic rodent peroxisomal proliferation markers were induced in monkeys after 17 weeks but not after a 4-week treatment, a fact also observed in human hepatocytes after 72 h but not 24 h of EMD treatment. In human hepatocyte cultures, EMD showed similar gene expression profiles and induction of P450 activities as in monkeys, indicating that the monkey is predictive for human P450 induction by EMD. In addition, EMD induced a similar gene expression pattern as the PPARalpha agonist fenofibrate in primary rat and human hepatocyte cultures. In conclusion, these data showed an excellent correlation of in vivo data on DME gene expression and activity levels with results generated in hepatocyte monolayer cultures, enabling a solid estimation of human P450 induction. This study also clearly highlighted major differences between primates and rodents in the regulation of major inducible P450s, with evidence of CYP3A and CYP2C inducibility by PPARalpha agonists in monkeys and humans.