MicroRNA Profile of Circulating CD4-positive Regulatory T Cells in Human Adults and Impact of Differentially Expressed MicroRNAs on Expression of Two Genes Essential to Their Function.

Regulatory T cells (Tregs) are characterized by a high expression of IL-2 receptor α chain (CD25) and of forkhead box P3 (FOXP3), the latter being essential for their development and function. Another major player in the regulatory function is the cytotoxic T-lymphocyte associated molecule-4 (CTLA-4) that inhibits cytotoxic responses. However, the regulation of CTLA-4 expression remains less well explored. We therefore studied the microRNA signature of circulating CD4(+) Tregs isolated from adult healthy donors and identified a signature composed of 15 differentially expressed microRNAs. Among those, miR-24, miR-145, and miR-210 were down-regulated in Tregs compared with controls and were found to have potential target sites in the 3'-UTR of FOXP3 and CTLA-4; miR-24 and miR-210 negatively regulated FOXP3 expression by directly binding to their two target sites in its 3'-UTR. On the other hand, miR-95, which is highly expressed in adult peripheral blood Tregs, positively regulated FOXP3 expression via an indirect mechanism yet to be identified. Finally, we showed that miR-145 negatively regulated CTLA-4 expression in human CD4(+) adult peripheral blood Tregs by binding to its target site in CTLA-4 transcript 3'-UTR. To our knowledge, this is the first identification of a human adult peripheral blood CD4(+) Treg microRNA signature. Moreover, unveiling one mechanism regulating CTLA-4 expression is novel and may lead to a better understanding of the regulation of this crucial gene.

Differentially expressed gene profile in the 6-hydroxy-dopamine-induced cell culture model of Parkinson's disease.

Parkinson's disease (PD) is a chronic neurodegenerative disorder characterized by progressive loss of dopaminergic (DA) neurons of the substantia nigra pars compacta with unknown aetiology. 6-Hydroxydopamine (6-OHDA) treatment of neuronal cells is an established in vivo model for mimicking the effect of oxidative stress found in PD brains. We examined the effects of 6-OHDA treatment on human neuroblastoma cells (SH-SY5Y) and primary mesencephalic cultures. Using a reverse arbitrarily primed polymerase chain reaction (RAP-PCR) approach we generated reproducible genetic fingerprints of differential expression levels in cell cultures treated with 6-OHDA. Of the resulting sequences, 23 showed considerable homology to known human coding sequences. The results of the RAP-PCR were validated by reverse transcription PCR, real-time PCR and, for selected genes, by Western blot analysis and immunofluorescence. In four cases, [tomoregulin-1 (TMEFF-1), collapsin response mediator protein 1 (CRMP-1), neurexin-1, and phosphoribosylaminoimidazole synthetase (GART)], a down-regulation of mRNA and protein levels was detected. Further studies will be necessary on the physiological role of the identified proteins and their impact on pathways leading to neurodegeneration in PD.

Gene expression variation and expression quantitative trait mapping of human chromosome 21 genes.

Inter-individual differences in gene expression are likely to account for an important fraction of phenotypic differences, including susceptibility to common disorders. Recent studies have shown extensive variation in gene expression levels in humans and other organisms, and that a fraction of this variation is under genetic control. We investigated the patterns of gene expression variation in a 25 Mb region of human chromosome 21, which has been associated with many Down syndrome (DS) phenotypes. Taqman real-time PCR was used to measure expression variation of 41 genes in lymphoblastoid cells of 40 unrelated individuals. For 25 genes found to be differentially expressed, additional analysis was performed in 10 CEPH families to determine heritabilities and map loci harboring regulatory variation. Seventy-six percent of the differentially expressed genes had significant heritabilities, and genomewide linkage analysis led to the identification of significant eQTLs for nine genes. Most eQTLs were in trans, with the best result (P=7.46 x 10(-8)) obtained for TMEM1 on chromosome 12q24.33. A cis-eQTL identified for CCT8 was validated by performing an association study in 60 individuals from the HapMap project. SNP rs965951 located within CCT8 was found to be significantly associated with its expression levels (P=2.5 x 10(-5)) confirming cis-regulatory variation. The results of our study provide a representative view of expression variation of chromosome 21 genes, identify loci involved in their regulation and suggest that genes, for which expression differences are significantly larger than 1.5-fold in control samples, are unlikely to be involved in DS-phenotypes present in all affected individuals.

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

FTZ-F1-related orphan receptors in Xenopus laevis: transcriptional regulators differentially expressed during early embryogenesis.

Orphan receptors of the FTZ-F1-related group of nuclear receptors (xFF1r) were identified in Xenopus laevis by isolation of cDNAs from a neurula stage library. Two cDNAs were found, which encode full length, highly related receptor proteins, xFF1rA and B, whose closet relative known so far is the murine LRH-1 orphan receptor. xFF1rA protein expressed by a recombinant vaccinia virus system specifically binds to FTZ-F1 response elements (FRE; PyCAAGGPyCPu). In cotransfection studies, xFF1rA constitutively activates transcription, in a manner dependent on the number of FREs. The amounts of at least four mRNAs encoding full-length receptors greatly increase between gastrula and early tailbud stages and decrease at later stages. At early tailbud stages, xFTZ-F1-related antigens are found in all nuclei of the embryo.

Meta-analysis of estrogen response in MCF-7 distinguishes early target genes involved in signaling and cell proliferation from later target genes involved in cell cycle and DNA repair.

ABSTRACT: BACKGROUND: Many studies have been published outlining the global effects of 17 beta-estradiol (E2) on gene expression in human epithelial breast cancer derived MCF-7 cells. These studies show large variation in results, reporting between ~100 and ~1500 genes regulated by E2, with poor overlap. RESULTS: We performed a meta-analysis of these expression studies, using the Rank product method to obtain a more accurate and stable list of the differentially expressed genes, and of pathways regulated by E2. We analyzed 9 time-series data sets, concentrating on response at 3-4 hrs (early) and at 24 hrs (late). We found >1000 statistically significant probe sets after correction for multiple testing at 3-4 hrs, and >2000 significant probe sets at 24 hrs. Differentially expressed genes were examined by pathway analysis. This revealed 15 early response pathways, mostly related to cell signaling and proliferation, and 20 late response pathways, mostly related to breast cancer, cell division, DNA repair and recombination. CONCLUSIONS: Our results show that meta-analysis identified more differentially expressed genes than the individual studies, and that these genes act together in networks. These results provide new insight into E2 regulated mechanisms, especially in the context of breast cancer.

Identification of HMX1 target genes: a predictive promoter model approach.

PURPOSE: A homozygous mutation in the H6 family homeobox 1 (HMX1) gene is responsible for a new oculoauricular defect leading to eye and auricular developmental abnormalities as well as early retinal degeneration (MIM 612109). However, the HMX1 pathway remains poorly understood, and in the first approach to better understand the pathway's function, we sought to identify the target genes. METHODS: We developed a predictive promoter model (PPM) approach using a comparative transcriptomic analysis in the retina at P15 of a mouse model lacking functional Hmx1 (dmbo mouse) and its respective wild-type. This PPM was based on the hypothesis that HMX1 binding site (HMX1-BS) clusters should be more represented in promoters of HMX1 target genes. The most differentially expressed genes in the microarray experiment that contained HMX1-BS clusters were used to generate the PPM, which was then statistically validated. Finally, we developed two genome-wide target prediction methods: one that focused on conserving PPM features in human and mouse and one that was based on the co-occurrence of HMX1-BS pairs fitting the PPM, in human or in mouse, independently. RESULTS: The PPM construction revealed that sarcoglycan, gamma (35kDa dystrophin-associated glycoprotein) (Sgcg), teashirt zinc finger homeobox 2 (Tshz2), and solute carrier family 6 (neurotransmitter transporter, glycine) (Slc6a9) genes represented Hmx1 targets in the mouse retina at P15. Moreover, the genome-wide target prediction revealed that mouse genes belonging to the retinal axon guidance pathway were targeted by Hmx1. Expression of these three genes was experimentally validated using a quantitative reverse transcription PCR approach. The inhibitory activity of Hmx1 on Sgcg, as well as protein tyrosine phosphatase, receptor type, O (Ptpro) and Sema3f, two targets identified by the PPM, were validated with luciferase assay. CONCLUSIONS: Gene expression analysis between wild-type and dmbo mice allowed us to develop a PPM that identified the first target genes of Hmx1.

t(15;21) translocations leading to the concurrent downregulation of RUNX1 and its transcription factor partner genes SIN3A and TCF12 in myeloid disorders.

Through a combined approach integrating RNA-Seq, SNP-array, FISH and PCR techniques, we identified two novel t(15;21) translocations leading to the inactivation of RUNX1 and its partners SIN3A and TCF12. One is a complex t(15;21)(q24;q22), with both breakpoints mapped at the nucleotide level, joining RUNX1 to SIN3A and UBL7-AS1 in a patient with myelodysplasia. The other is a recurrent t(15;21)(q21;q22), juxtaposing RUNX1 and TCF12, with an opposite transcriptional orientation, in three myeloid leukemia cases. Since our transcriptome analysis indicated a significant number of differentially expressed genes associated with both translocations, we speculate an important pathogenetic role for these alterations involving RUNX1.

A simple genetic basis for complex social behaviour mediates widespread gene expression differences.

A remarkable social polymorphism is controlled by a single Mendelian factor in the fire ant Solenopsis invicta. A genomic element marked by the gene Gp-9 determines whether workers tolerate one or many fertile queens in their colony. Gp-9 was recently shown to be part of a supergene with two nonrecombining variants, SB and Sb. SB/SB and SB/Sb queens differ in how they initiate new colonies, and in many physiological traits, for example odour and maturation rate. To understand how a single genetic element can affect all these traits, we used a microarray to compare gene expression patterns between SB/SB and SB/Sb queens of three different age classes: 1-day-old unmated queens, 11-day-old unmated queens and mated, fully reproductive queens collected from mature field colonies. The number of genes that were differentially expressed between SB/SB and SB/Sb queens of the same age class was smallest in 1-day-old queens, maximal in 11-day-old queens and intermediate in reproductive queens. Gene ontology analysis showed that SB/SB queens upregulate reproductive genes faster than SB/Sb queens. For all age classes, genes inside the supergene were overrepresented among the differentially expressed genes. Consistent with the hypothesized greater number of transposons in the Sb supergene, 13 transposon genes were upregulated in SB/Sb queens. Viral genes were also upregulated in SB/Sb mature queens, consistent with the known greater parasite load in colonies headed by SB/Sb queens compared with colonies headed by SB/SB queens. Eighteen differentially expressed genes between reproductive queens were involved in chemical signalling. Our results suggest that many genes in the supergene are involved in regulating social organization and queen phenotypes in fire ants.

Phenotypic consequences of copy number variation: insights from Smith-Magenis and Potocki-Lupski syndrome mouse models.

A large fraction of genome variation between individuals is comprised of submicroscopic copy number variation of genomic DNA segments. We assessed the relative contribution of structural changes and gene dosage alterations on phenotypic outcomes with mouse models of Smith-Magenis and Potocki-Lupski syndromes. We phenotyped mice with 1n (Deletion/+), 2n (+/+), 3n (Duplication/+), and balanced 2n compound heterozygous (Deletion/Duplication) copies of the same region. Parallel to the observations made in humans, such variation in gene copy number was sufficient to generate phenotypic consequences: in a number of cases diametrically opposing phenotypes were associated with gain versus loss of gene content. Surprisingly, some neurobehavioral traits were not rescued by restoration of the normal gene copy number. Transcriptome profiling showed that a highly significant propensity of transcriptional changes map to the engineered interval in the five assessed tissues. A statistically significant overrepresentation of the genes mapping to the entire length of the engineered chromosome was also found in the top-ranked differentially expressed genes in the mice containing rearranged chromosomes, regardless of the nature of the rearrangement, an observation robust across different cell lineages of the central nervous system. Our data indicate that a structural change at a given position of the human genome may affect not only locus and adjacent gene expression but also "genome regulation." Furthermore, structural change can cause the same perturbation in particular pathways regardless of gene dosage. Thus, the presence of a genomic structural change, as well as gene dosage imbalance, contributes to the ultimate phenotype.

Exploring Parallels Between Molecular Changes Induced in PNS by Aging and Demyelinating Neuropathies

The peripheral nervous system (PNS) is involved in many age-dependent neurological deficits, including numbness, pain, restless legs, trouble with walking and balance that are commonly found in the elderly. These symptoms generally result from demyelination and/or loss of axonal integrity. However, the precise identity of age-regulated molecular changes in either neuronal or glial compartments of the nerve is unclear. Interestingly, these deficiencies are also present in inherited neuropathies, where the expressivity of the rapid and early onset phenotypes is undeniably more severe than in normal aging. Nevertheless, especially the molecular changes underlying loss of axonal integrity in neuropathy condition are also poorly understood. To unravel molecular mechanisms affected by PNS aging, we used wildtype mice at 17 time-points from day of birth until senescence (28 months-old). For the neuropathy study, we focused on 56 day-old Schwann cell-specific neuropathy-inducing mutants, MPZCre/1/ LpinfE2-3/fE2-3 and MPZCre/1/ScapfE1/fE1 mice, that have, at this age, already developed neuropathic symptoms. Transcriptomes of dissected Schwann cell-containing endoneurium or sensory neuron-containing dorsal root ganglia have been analyzed throughout time or genotypes, using Illumina Bead Chips. Following data validation, we identified groups of differentially expressed genes in the development, aging and in the neuropathic mutants, in both glial and neuronal compartments. We detected substantial differences in the dynamics of changes in gene expression during development and aging between these two compartments. Furthermore, considering the above-mentioned phenotypic similarities, we integrated aging and mutant data. Interestingly, we observed that there are some parallels at the molecular level between processes involved in aging, which leads to less severe and more progressive PNS alterations, and in the rapid onset peripheral neuropathies. Apart from helping the understanding of molecular alterations underlying age-related PNS phenotypes, this data should also contribute to the identification of pathways that could be used as targets for therapeutical approaches to prevent complications associated with both aging and inherited forms of neuropathies.

EuroDia: a beta-cell gene expression resource.

Type 2 diabetes mellitus (T2DM) is a major disease affecting nearly 280 million people worldwide. Whilst the pathophysiological mechanisms leading to disease are poorly understood, dysfunction of the insulin-producing pancreatic beta-cells is key event for disease development. Monitoring the gene expression profiles of pancreatic beta-cells under several genetic or chemical perturbations has shed light on genes and pathways involved in T2DM. The EuroDia database has been established to build a unique collection of gene expression measurements performed on beta-cells of three organisms, namely human, mouse and rat. The Gene Expression Data Analysis Interface (GEDAI) has been developed to support this database. The quality of each dataset is assessed by a series of quality control procedures to detect putative hybridization outliers. The system integrates a web interface to several standard analysis functions from R/Bioconductor to identify differentially expressed genes and pathways. It also allows the combination of multiple experiments performed on different array platforms of the same technology. The design of this system enables each user to rapidly design a custom analysis pipeline and thus produce their own list of genes and pathways. Raw and normalized data can be downloaded for each experiment. The flexible engine of this database (GEDAI) is currently used to handle gene expression data from several laboratory-run projects dealing with different organisms and platforms. Database URL: http://eurodia.vital-it.ch.

Effects of mitotane on gene expression in the adrenocortical cell line NCI-H295R: a microarray study.

Aim: The adrenolytic agent mitotane is widely used in the treatment of adrenocortical cancer; however, its mechanism of action is poorly elucidated. We have studied mitotane-induced mRNA expression changes in the NCI-H295R adrenocortical cancer cell line. Materials & methods: Cell viability and hormone assays were used to select the optimal mitotane concentration effectively inhibiting hormone secretion without affecting cell viability. RNA isolated from cultures treated for 48 and 72 h was subjected to Agilent 4×44K microarray platforms. Microarray results were validated by quantitative reverse-transcription PCR. Results: Altogether, 117 significantly differentially expressed genes were detected at 48 h and 72 h (p < 0.05) in mitotane-treated samples relative to controls. Three significantly underexpressed genes involved in steroid hormone biosynthesis (HSD3B1, HSD3B2 and CYP21A2) and four significantly overexpressed genes (GDF15, ALDH1L2, TRIB3 and SERPINE2) have been validated. Conclusion: Gene-expression changes might be involved in the adrenal action of mitotane and in the inhibition of hormone secretion. Original submitted 20 January 2012; Revision submitted 17 May 2012.

Genome-Wide Analysis of Salicylate and Dibenzofuran Metabolism in Sphingomonas Wittichii RW1.

Sphingomonas wittichii RW1 is a bacterium isolated for its ability to degrade the xenobiotic compounds dibenzodioxin and dibenzofuran (DBF). A number of genes involved in DBF degradation have been previously characterized, such as the dxn cluster, dbfB, and the electron transfer components fdx1, fdx3, and redA2. Here we use a combination of whole genome transcriptome analysis and transposon library screening to characterize RW1 catabolic and other genes implicated in the reaction to or degradation of DBF. To detect differentially expressed genes upon exposure to DBF, we applied three different growth exposure experiments, using either short DBF exposures to actively growing cells or growing them with DBF as sole carbon and energy source. Genome-wide gene expression was examined using a custom-made microarray. In addition, proportional abundance determination of transposon insertions in RW1 libraries grown on salicylate or DBF by ultra-high throughput sequencing was used to infer genes whose interruption caused a fitness loss for growth on DBF. Expression patterns showed that batch and chemostat growth conditions, and short or long exposure of cells to DBF produced very different responses. Numerous other uncharacterized catabolic gene clusters putatively involved in aromatic compound metabolism increased expression in response to DBF. In addition, only very few transposon insertions completely abolished growth on DBF. Some of those (e.g., in dxnA1) were expected, whereas others (in a gene cluster for phenylacetate degradation) were not. Both transcriptomic data and transposon screening suggest operation of multiple redundant and parallel aromatic pathways, depending on DBF exposure. In addition, increased expression of other non-catabolic genes suggests that during initial exposure, S. wittichii RW1 perceives DBF as a stressor, whereas after longer exposure, the compound is recognized as a carbon source and metabolized using several pathways in parallel.

Gene expression profiling in breast cancer: understanding the molecular basis of histologic grade to improve prognosis.

BACKGROUND: Histologic grade in breast cancer provides clinically important prognostic information. However, 30%-60% of tumors are classified as histologic grade 2. This grade is associated with an intermediate risk of recurrence and is thus not informative for clinical decision making. We examined whether histologic grade was associated with gene expression profiles of breast cancers and whether such profiles could be used to improve histologic grading. METHODS: We analyzed microarray data from 189 invasive breast carcinomas and from three published gene expression datasets from breast carcinomas. We identified differentially expressed genes in a training set of 64 estrogen receptor (ER)-positive tumor samples by comparing expression profiles between histologic grade 3 tumors and histologic grade 1 tumors and used the expression of these genes to define the gene expression grade index. Data from 597 independent tumors were used to evaluate the association between relapse-free survival and the gene expression grade index in a Kaplan-Meier analysis. All statistical tests were two-sided. RESULTS: We identified 97 genes in our training set that were associated with histologic grade; most of these genes were involved in cell cycle regulation and proliferation. In validation datasets, the gene expression grade index was strongly associated with histologic grade 1 and 3 status; however, among histologic grade 2 tumors, the index spanned the values for histologic grade 1-3 tumors. Among patients with histologic grade 2 tumors, a high gene expression grade index was associated with a higher risk of recurrence than a low gene expression grade index (hazard ratio = 3.61, 95% confidence interval = 2.25 to 5.78; P < .001, log-rank test). CONCLUSIONS: Gene expression grade index appeared to reclassify patients with histologic grade 2 tumors into two groups with high versus low risks of recurrence. This approach may improve the accuracy of tumor grading and thus its prognostic value.