Using GOstats to test gene lists for GO term association

Antiviral response dictated by choreographed cascade of transcription factors

The dendritic cell (DC) is a master regulator of immune responses. Pathogenic viruses subvert normal immune function in DCs through the expression of immune antagonists. Understanding how these antagonists interact with the host immune system requires knowledge of the underlying genetic regulatory network that operates during an uninhibited antiviral response. To isolate and identify this network, we studied DCs infected with Newcastle disease virus, which is able to stimulate innate immunity and DC maturation through activation of RIG-I signaling, but lacks the ability to evade the human IFN response. To analyze this experimental model, we developed a new approach integrating genome-wide expression kinetics and time-dependent promoter analysis. We found that the genetic program underlying the antiviral cell-state transition during the first 18 h postinfection could be explained by a single convergent regulatory network. Gene expression changes were driven by a stepwise multifactor cascading control mechanism, where the specific transcription factors controlling expression changed over time. Within this network, most individual genes were regulated by multiple factors, indicating robustness against virus-encoded immune evasion genes. In addition to effectively recapitulating current biological knowledge, we predicted, and validated experimentally, antiviral roles for several novel transcription factors. More generally, our results show how a genetic program can be temporally controlled through a single regulatory network to achieve the large-scale genetic reprogramming characteristic of cell-state transitions.

Aberrant gene expression profile in a mouse model of endometriosis mirrors that observed in women

OBJECTIVE

To define the altered gene expression profile of endometriotic lesions in a mouse model of surgically induced endometriosis.

DESIGN

Autologous experimental mouse model.

SETTING

Medical school department.

ANIMAL(S)

Adult C57Bl6 mice.

INTERVENTION(S)

Endometriosis was surgically induced by autotransplantation of uterine tissue to the intestinal mesentery. Endometriotic lesions and eutopic uteri were recovered at 3 or 29 days after induction.

MAIN OUTCOME MEASURE(S)

Altered gene expression was measured in the endometriotic lesion relative to the eutopic uterus by genome-wide complementary DNA microarray analysis and was confirmed by real-time reverse transcriptase-polymerase chain reaction for six genes. Relevant categories of altered genes were identified using gene ontology analysis to determine groups of genes enriched for altered expression.

RESULT(S)

The expression of 479 and 114 genes was altered in the endometriotic lesion compared with the eutopic uterus at 3 or 29 days after induction, respectively. Gene ontology enrichment analysis revealed that genes associated with the extracellular matrix, cell adhesions, immune function, cell growth, and angiogenesis were altered in the endometriotic lesion compared with the eutopic uterus.

CONCLUSION(S)

According to gene expression analysis, the mouse model of surgically induced endometriosis is a good model for studying the pathophysiology and treatment of endometriosis.

High-fat diet leads to tissue-specific changes reflecting risk factors for diseases in DBA/2J mice

The aim of this study was to characterize the responses of individual tissues to high-fat feeding as a function of mass, fat composition, and transcript abundance. We examined a panel of eight tissues [5 white adipose tissues (WAT), brown adipose tissue (BAT), liver, muscle] obtained from DBA/2J mice on either a standard breeding diet (SBD) or a high-fat diet (HFD). HFD led to weight gain, decreased insulin sensitivity, and tissue-specific responses, including inflammation, in these mice. The dietary fatty acids were partially metabolized and converted in both liver and fat tissues. Saturated fatty acids (SFA) were converted in the liver to monounsaturated fatty acids (MUFA) and polyunsaturated fatty acids (PUFA), and oleic acid (C18:1) was the preferred MUFA for storage of excess energy in all tissues of HFD-fed mice. Transcriptional changes largely reflected the tissue-specific fat deposition. SFA were negatively correlated with genes in the collagen family and processes involving the extracellular matrix. We propose a novel role of the tryptophan hydroxylase 2 (Tph2) gene in adipose tissues of diet-induced obesity. Tissue-specific responses to HFD were identified. Liver steatosis was evident in HFD-fed mice. Gonadal, retroperitoneal and subcutaneous adipose tissue and BAT exhibited severe inflammatory and immune responses. Mesenteric adipose tissue was the most metabolically active adipose tissue. Gluteal adipose tissue had the highest mass gain but was sluggish in its metabolism. In HFD conditions, BAT functioned largely like WAT in its role as a depot for excess energy, whereas WAT played a role in thermogenesis.

Promiscuity of enhancer, coding and non-coding transcription functions in ultraconserved elements

BACKGROUND

Ultraconserved elements (UCEs) are highly constrained elements of mammalian genomes, whose functional role has not been completely elucidated yet. Previous studies have shown that some of them act as enhancers in mouse, while some others are expressed in both normal and cancer-derived human tissues. Only one UCE element so far was shown to present these two functions concomitantly, as had been observed in other isolated instances of single, non ultraconserved enhancer elements.

RESULTS

We used a custom microarray to assess the levels of UCE transcription during mouse development and integrated these data with published microarray and next-generation sequencing datasets as well as with newly produced PCR validation experiments. We show that a large fraction of non-exonic UCEs is transcribed across all developmental stages examined from only one DNA strand. Although the nature of these transcripts remains a mistery, our meta-analysis of RNA-Seq datasets indicates that they are unlikely to be short RNAs and that some of them might encode nuclear transcripts. In the majority of cases this function overlaps with the already established enhancer function of these elements during mouse development. Utilizing several next-generation sequencing datasets, we were further able to show that the level of expression observed in non-exonic UCEs is significantly higher than in random regions of the genome and that this is also seen in other regions which act as enhancers.

CONCLUSION

Our data shows that the concurrent presence of enhancer and transcript function in non-exonic UCE elements is more widespread than previously shown. Moreover through our own experiments as well as the use of next-generation sequencing datasets, we were able to show that the RNAs encoded by non-exonic UCEs are likely to be long RNAs transcribed from only one DNA strand.

The submergence tolerance regulator Sub1A mediates stress-responsive expression of AP2/ERF transcription factors

We previously characterized the rice (Oryza sativa) Submergence1 (Sub1) locus encoding three ethylene-responsive factor (ERF) transcriptional regulators. Genotypes carrying the Sub1A-1 allele are tolerant of prolonged submergence. To elucidate the mechanism of Sub1A-1-mediated tolerance, we performed transcriptome analyses comparing the temporal submergence response of Sub1A-1-containing tolerant M202(Sub1) with the intolerant isoline M202 lacking this gene. We identified 898 genes displaying Sub1A-1-dependent regulation. Integration of the expression data with publicly available metabolic pathway data identified submergence tolerance-associated pathways governing anaerobic respiration, hormone responses, and antioxidant systems. Of particular interest were a set of APETALA2 (AP2)/ERF family transcriptional regulators that are associated with the Sub1A-1-mediated response upon submergence. Visualization of expression patterns of the AP2/ERF superfamily members in a phylogenetic context resolved 12 submergence-regulated AP2/ERFs into three putative functional groups: (1) anaerobic respiration and cytokinin-mediated delay in senescence via ethylene accumulation during submergence (three ERFs); (2) negative regulation of ethylene-dependent gene expression (five ERFs); and (3) negative regulation of gibberellin-mediated shoot elongation (four ERFs). These results confirm that the presence of Sub1A-1 impacts multiple pathways of response to submergence.

Transcription factor Stat3 stimulates metastatic behavior of human prostate cancer cells in vivo, whereas Stat5b has a preferential role in the promotion of prostate cancer cell viability and tumor growth

Identification of the molecular changes that promote viability and metastatic behavior of prostate cancer is critical for the development of improved therapeutic interventions. Stat5a/b and Stat3 are both constitutively active in locally-confined and advanced prostate cancer, and both transcription factors have been reported to be critical for the viability of prostate cancer cells. We recently showed that Stat3 promotes metastatic behavior of human prostate cancer cells not only in vitro but also in an in vivo experimental metastases model. In this work, we compare side-by-side Stat5a/b versus Stat3 in the promotion of prostate cancer cell viability, tumor growth, and induction of metastatic colonization in vivo. Inhibition of Stat5a/b induced massive death of prostate cancer cells in culture and reduced both subcutaneous and orthotopic prostate tumor growth, whereas Stat3 had a predominant role over Stat5a/b in promoting metastases formation of prostate cancer cells in vivo in nude mice. The molecular mechanisms underlying the differential biological effects induced by these two transcription factors involve largely different sets of genes regulated by Stat5a/b versus Stat3 in human prostate cancer model systems. Of the two Stat5 homologs, Stat5b was more important for supporting growth of prostate cancer cells than Stat5a. This work provides the first mechanistic comparison of the biological effects induced by transcription factors Stat5a/b versus Stat3 in prostate cancer.

Histone H3 lysine 27 methylation asymmetry on developmentally-regulated promoters distinguish the first two lineages in mouse preimplantation embryos

First lineage specification in the mammalian embryo leads to formation of the inner cell mass (ICM) and trophectoderm (TE), which respectively give rise to embryonic and extraembryonic tissues. We show here that this first differentiation event is accompanied by asymmetric distribution of trimethylated histone H3 lysine 27 (H3K27me3) on promoters of signaling and developmentally-regulated genes in the mouse ICM and TE. A genome-wide survey of promoter occupancy by H3K4me3 and H3K27me3 indicates that both compartments harbor promoters enriched in either modification, and promoters co-enriched in trimethylated H3K4 and H3K27 linked to developmental and signaling functions. The majority of H3K4/K27me3 co-enriched promoters are distinct between the two lineages, primarily due to differences in the distribution of H3K27me3. Derivation of embryonic stem cells leads to significant losses and gains of H3K4/K27me3 co-enriched promoters relative to the ICM, with distinct contributions of (de)methylation events on K4 and K27. Our results show histone trimethylation asymmetry on promoters in the first two developmental lineages, and highlight an epigenetic skewing associated with embryonic stem cell derivation.

The transcriptome of cervical ripening in human pregnancy before the onset of labor at term: identification of novel molecular functions involved in this process

OBJECTIVE

The aim of this study was to identify changes in the cervical transcriptome in the human uterine cervix as a function of ripening before the onset of labor.

STUDY DESIGN

Human cervical tissue was obtained from women at term not in labor with ripe (n = 11) and unripe (n = 11) cervices and profiled using Affymetrix GeneChip HGU133Plus2.0 arrays. Gene expression was analyzed using a moderated t-test (False Discovery Rate 5%). Gene ontology and pathway analysis were performed. Quantitative reverse transcription-polymerase chain reaction (qRT-PCR) was used for confirmation of selected differentially expressed genes.

RESULTS

(1) Ninety-one genes were differentially expressed between ripe and unripe groups. (2) Cervical ripening was associated with enrichment of specific biological processes (e.g. cell adhesion, regulation of anatomical structure), pathways and 11 molecular functions (e.g. extracelluar matrix (ECM)-structural constituent, protein binding, glycosaminoglycan binding). (3) qRT-PCR confirmed that 9 of 11 tested differentially expressed genes (determined by microarray) were upregulated in a ripe cervix (e.g. MYOCD, VCAN, THBS1, COL5A1). (4) Twenty-three additional genes related to ECM metabolism and adhesion molecules were differentially regulated (by qRT-PCR) in ripe cervices.

CONCLUSION

(1) This is the first description of the changes in the human cervical transcriptome with ripening before the onset of labor. (2) Biological processes, pathways and molecular functions were identified with the use of this unbiased approach. (3) In contrast to cervical dilation after term labor, inflammation-related genes did not emerge as differentially regulated with cervical ripening. (4) Myocardin was identified as a novel gene upregulated in human cervical ripening.

Low malignant potential tumors with micropapillary features are molecularly similar to low-grade serous carcinoma of the ovary

OBJECTIVE

Low-grade serous carcinoma (LGSC) is a chemoresistant ovarian neoplasm thought to potentially arise in a background of low malignant potential tumors (LMP), which are typically non-aggressive. However, LMP with micropapillary features (LMP-MP) have more aggressive clinical behavior and may represent an intermediate in progression to LGSC. The objective of this study was to obtain and compare gene expression profiles of LMP, LMP-MP and LGSC to determine if LMP-MP more closely resembles LGSC, and to identify genes involved in LGSC carcinogenesis.

METHODS

Epithelial cells from LMP (n=17), LMP-MP (n=9) and LGSC (n=11) were isolated by laser capture microdissection. RNA was extracted, reverse transcribed to cDNA, amplified and hybridized to Affymetrix U133 Plus2 genechip arrays. Gene expression data were checked for quality, filtered and significantly altered genes between subgroups were identified. Differential expression of selected genes was verified by RT-qPCR and immunohistochemistry.

RESULTS

Gene expression analysis identified differential expression between LMP and LMP-MP, LMP and LGSC but not LMP-MP and LGSC. Integration of differentially expressed genes into the protein interaction database CytoScape highlighted gene products in the MAPK pathway as differentially regulated between LMP and LGSC. Four genes were selected and validated by RT-qPCR performed on microarray samples (n=15) and immunohistochemistry on a representative microarray (n=57).

CONCLUSION

The gene expression profile of LMP-MP is similar to LGSC and distinct from LMP, reflecting their more aggressive clinical behavior. Candidate genes in the MAPK pathway were highlighted which may play a role in LGSC carcinogenesis and indicate potential therapeutic targets.

Long-term genome-wide blood RNA expression profiles yield novel molecular response candidates for IFN-β-1b treatment in relapsing remitting MS

Aims: In multiple sclerosis patients, treatment with recombinant IFN-β (rIFN-β) is partially efficient in reducing clinical exacerbations. However, its molecular mechanism of action is still under scrutiny. Materials & methods: We used DNA microarrays (Affymetrix, CA, USA) and peripheral mononuclear blood cells from 25 relapsing remitting multiple sclerosis patients to analyze the longitudinal transcriptional profile within 2 years of rIFN-β administration. Sets of differentially expressed genes were attained by applying a combination of independent criteria, thereby providing efficient data curation and gene filtering that accounted for technical and biological noise. Gene ontology term-association analysis and scientific literature text mining were used to explore evidence of gene interaction. Results: Post-therapy initiation, we identified 42 (day 2), 175 (month 1), 103 (month 12) and 108 (month 24) differentially expressed genes. Increased expression of established IFN-β marker genes, as well as differential expression of circulating IFN-β-responsive candidate genes, were observed. MS4A1 (CD20), a known target of B-cell depletion therapy, was significantly downregulated after one month. CMPK2, FCER1A, and FFAR2 appeared as hitherto unrecognized multiple sclerosis treatment-related differentially expressed genes that were consistently modulated over time. Overall, 84 interactions between 54 genes were attained, of which two major gene networks were identified at an earlier stage of therapy: the first (n = 15 genes) consisted of mostly known IFN-β-activated genes, whereas the second (n = 12) mainly contained downregulated genes that to date have not been associated with IFN-β effects in multiple sclerosis array research. Conclusion: We achieved both a broadening of the knowledge of IFN-β mechanism-of-action-related constituents and the identification of time-dependent interactions between IFN-β regulated genes.

Gene expression changes in the host response between resistant and susceptible inbred mouse strains after influenza A infection

Inbred mouse strains exhibit differences in susceptibility to influenza A infections. However, the molecular mechanisms underlying these differences are unknown. Therefore, we infected a highly susceptible mouse strain (DBA/2J) and a resistant strain (C57BL/6J) with influenza A H1N1 (PR8) and performed genome-wide expression analysis. We found genes expressed in lung epithelium that were specifically down-regulated in DBA/2J mice, whereas a cluster of genes on chromosome 3 was only down-regulated in C57BL/6J. In both mouse strains, chemokines, cytokines and interferon-response genes were up-regulated, indicating that the main innate immune defense pathways were activated. However, many immune response genes were up-regulated in DBA/2J much stronger than in C57BL/6J, and several immune response genes were exclusively regulated in DBA/2J. Thus, susceptible DBA/2J mice showed a hyper-inflammatory response. This response is similar to infections with highly pathogenic influenza virus and may serve as a paradigm for a hyper-inflammatory host response to influenza A virus.

Heteronemin, a spongean sesterterpene, inhibits TNF alpha-induced NF-kappa B activation through proteasome inhibition and induces apoptotic cell death

In this study, we investigated the biological effects of heteronemin, a marine sesterterpene isolated from the sponge Hyrtios sp. on chronic myelogenous leukemia cells. To gain further insight into the molecular mechanisms triggered by this compound, we initially performed DNA microarray profiling and determined which genes respond to heteronemin stimulation in TNFalpha-treated cells and which genes display an interaction effect between heteronemin and TNFalpha. Within the differentially regulated genes, we found that heteronemin was affecting cellular processes including cell cycle, apoptosis, mitogen-activated protein kinases (MAPKs) pathway and the nuclear factor kappaB (NF-kappaB) signaling cascade. We confirmed in silico experiments regarding NF-kappaB inhibition by reporter gene analysis, electrophoretic mobility shift analysis and I-kappaB degradation. In order to assess the underlying molecular mechanisms, we determined that heteronemin inhibits both trypsin and chymotrypsin-like proteasome activity at an IC(50) of 0.4 microM. Concomitant to the inhibition of the NF-kappaB pathway, we also observed a reduction in cellular viability. Heteronemin induces apoptosis as shown by annexin V-FITC/propidium iodide-staining, nuclear morphology analysis, pro-caspase-3, -8 and -9 and poly(ADP-ribose) polymerase (PARP) cleavage as well as truncation of Bid. Altogether, results show that this compound has potential as anti-inflammatory and anti-cancer agent.

Genome-wide RNAi screen identifies human host factors crucial for influenza virus replication

Influenza A virus, being responsible for seasonal epidemics and reoccurring pandemics, represents a worldwide threat to public health. High mutation rates facilitate the generation of viral escape mutants, rendering vaccines and drugs directed against virus-encoded targets potentially ineffective. In contrast, targeting host cell determinants temporarily dispensable for the host but crucial for virus replication could prevent viral escape. Here we report the discovery of 287 human host cell genes influencing influenza A virus replication in a genome-wide RNA interference (RNAi) screen. Using an independent assay we confirmed 168 hits (59%) inhibiting either the endemic H1N1 (119 hits) or the current pandemic swine-origin (121 hits) influenza A virus strains, with an overlap of 60%. Notably, a subset of these common hits was also essential for replication of a highly pathogenic avian H5N1 strain. In-depth analyses of several factors provided insights into their infection stage relevance. Notably, SON DNA binding protein (SON) was found to be important for normal trafficking of influenza virions to late endosomes early in infection. We also show that a small molecule inhibitor of CDC-like kinase 1 (CLK1) reduces influenza virus replication by more than two orders of magnitude, an effect connected with impaired splicing of the viral M2 messenger RNA. Furthermore, influenza-virus-infected p27(-/-) (cyclin-dependent kinase inhibitor 1B; Cdkn1b) mice accumulated significantly lower viral titres in the lung, providing in vivo evidence for the importance of this gene. Thus, our results highlight the potency of genome-wide RNAi screening for the dissection of virus-host interactions and the identification of drug targets for a broad range of influenza viruses.

Using transcriptomics to identify and validate novel biomarkers of human skeletal muscle cancer cachexia

Background

Cancer cachexia is a multi-organ tissue wasting syndrome that contributes to morbidity and mortality in many cancer patients. Skeletal muscle loss represents an established key feature yet there is no molecular understanding of the disease process. In fact, the postulated molecular regulators of cancer cachexia originate largely from pre-clinical models and it is unclear how these translate to the clinical environment.

Methods

Rectus abdominis muscle biopsies were obtained from 65 upper gastrointestinal (UGI) cancer patients during open surgery and RNA profiling was performed on a subset of this cohort (n = 21) using the Affymetrix U133+2 platform. Quantitative analysis revealed a gene signature, which underwent technical validation and independent confirmation in a separate clinical cohort.

Results

Quantitative significance analysis of microarrays produced an 83-gene signature that was able to identify patients with greater than 5% weight loss, while this molecular profile was unrelated to markers of systemic inflammation. Selected genes correlating with weight loss were validated using quantitative real-time PCR and independently studied as general cachexia biomarkers in diaphragm and vastus lateralis from a second cohort (n = 13; UGI cancer patients). CaMKIIβ correlated positively with weight loss in all muscle groups and CaMKII protein levels were elevated in rectus abdominis. TIE1 was also positively associated with weight loss in both rectus abdominis and vastus lateralis muscle groups while other biomarkers demonstrated tissue-specific expression patterns. Candidates selected from the pre-clinical literature, including FOXO protein and ubiquitin E3 ligases, were not related to weight loss in this human clinical study. Furthermore, promoter analysis identified that the 83 weight loss-associated genes had fewer FOXO binding sites than expected by chance.

Conclusion

We were able to discover and validate new molecular biomarkers of human cancer cachexia. The exercise activated genes CaMKIIβ and TIE1 related positively to weight-loss across muscle groups, indicating that this cachexia signature is not simply due to patient inactivity. Indeed, excessive CaMKIIβ activation is a potential mechanism for reduced muscle protein synthesis. Our genomics analysis also supports the view that the available preclinical models do not accurately reflect the molecular characteristics of human muscle from cancer cachexia patients.

Comparative gene expression analysis of zebrafish and mammals identifies common regulators in hematopoietic stem cells

Hematopoiesis in teleost fish is maintained in the kidney. We previously reported that Hoechst dye efflux activity of hematopoietic stem cells (HSCs) is highly conserved in vertebrates, and that Hoechst can be used to purify HSCs from teleost kidneys. Regulatory molecules that are strongly associated with HSC activity may also be conserved in vertebrates. In this study, we identified evolutionarily conserved molecular components in HSCs by comparing the gene expression profiles of zebrafish, murine, and human HSCs. Microarray data of zebrafish kidney side population cells (zSPs) showed that genes involved in cell junction and signal transduction tended to be up-regulated in zSPs, whereas genes involved in DNA replication tended to be down-regulated. These properties of zSPs were similar to those of mammalian HSCs. Overlapping gene expression analysis showed that 40 genes were commonly up-regulated in these 3 HSCs. Some of these genes, such as egr1, gata2, and id1, have been previously implicated in the regulation of HSCs. In situ hybridization in zebrafish kidney revealed that expression domains of egr1, gata2, and id1 overlapped with that of abcg2a, a marker for zSPs. These results suggest that the overlapping genes identified in this study are regulated in HSCs and play important roles in their functions.

Regional regulation of transcription in the chicken genome

Background

Over the past years, the relationship between gene transcription and chromosomal location has been studied in a number of different vertebrate genomes. Regional differences in gene expression have been found in several different species. The chicken genome, as the closest sequenced genome relative to mammals, is an important resource for investigating regional effects on transcription in birds and studying the regional dynamics of chromosome evolution by comparative analysis.

Results

We used gene expression data to survey eight chicken tissues and create transcriptome maps for all chicken chromosomes. The results reveal the presence of two distinct types of chromosomal regions characterized by clusters of highly or lowly expressed genes. Furthermore, these regions correlate highly with a number of genome characteristics. Regions with clusters of highly expressed genes have higher gene densities, shorter genes, shorter average intron and higher GC content compared to regions with clusters of lowly expressed genes. A comparative analysis between the chicken and human transcriptome maps constructed using similar panels of tissues suggests that the regions with clusters of highly expressed genes are relatively conserved between the two genomes.

Conclusions

Our results revealed the presence of a higher order organization of the chicken genome that affects gene expression, confirming similar observations in other species. These results will aid in the further understanding of the regional dynamics of chromosome evolution.

The microarray data used in this analysis have been submitted to NCBI GEO database under accession number GSE17108. The reviewer access link is: http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?token=tjwjpscyceqawjk&acc=GSE17108

Quantitative phosphoproteomics reveals widespread full phosphorylation site occupancy during mitosis

Eukaryotic cells replicate by a complex series of evolutionarily conserved events that are tightly regulated at defined stages of the cell division cycle. Progression through this cycle involves a large number of dedicated protein complexes and signaling pathways, and deregulation of this process is implicated in tumorigenesis. We applied high-resolution mass spectrometry-based proteomics to investigate the proteome and phosphoproteome of the human cell cycle on a global scale and quantified 6027 proteins and 20,443 unique phosphorylation sites and their dynamics. Co-regulated proteins and phosphorylation sites were grouped according to their cell cycle kinetics and compared to publicly available messenger RNA microarray data. Most detected phosphorylation sites and more than 20% of all quantified proteins showed substantial regulation, mainly in mitotic cells. Kinase-motif analysis revealed global activation during S phase of the DNA damage response network, which was mediated by phosphorylation by ATM or ATR or DNA-dependent protein kinases. We determined site-specific stoichiometry of more than 5000 sites and found that most of the up-regulated sites phosphorylated by cyclin-dependent kinase 1 (CDK1) or CDK2 were almost fully phosphorylated in mitotic cells. In particular, nuclear proteins and proteins involved in regulating metabolic processes have high phosphorylation site occupancy in mitosis. This suggests that these proteins may be inactivated by phosphorylation in mitotic cells.

Expression profiling in peripheral blood reveals signature for penetrance in DYT1 dystonia

DYT1 dystonia is an autosomal-dominantly inherited movement disorder, which is usually caused by a GAG deletion in the TOR1A gene. Due to the reduced penetrance of approximately 30-40%, the determination of the mutation in a subject is of limited use with regard to actual manifestation of symptoms. In the present study, we used Affymetrix oligonucleotide microarrays to analyze global gene expression in blood samples of 15 manifesting and 15 non-manifesting mutation carriers in order to identify a susceptibility profile beyond the GAG deletion which is associated with the manifestation of symptoms in DYT1 dystonia. We identified a genetic signature which distinguished between asymptomatic mutation carriers and symptomatic DYT1 patients with 86.7% sensitivity and 100% specificity. This genetic signature could correctly predict the disease state in an independent test set with a sensitivity of 87.5% and a specificity of 85.7%. Conclusively, this genetic signature might provide a possibility to distinguish DYT1 patients from asymptomatic mutation carriers.

Functional profiling methods in cancer

The introduction of new high-throughput methodologies such as DNA microarrays constitutes a major breakthrough in cancer research. The unprecedented amount of data produced by such technologies has opened new avenues for interrogating living systems although, at the same time, it has demanded of the development of new data analytical methods as well as new strategies for testing hypotheses. A history of early successful applications in cancer boosted the use of microarrays and fostered further applications in other fields. Keeping the pace with these technologies, bioinformatics offers new solutions for data analysis and, what is more important, permits the formulation of a new class of hypotheses inspired in systems biology, more oriented to pathways or, in general, to modules of functionally related genes. Although these analytical methodologies are new, some options are already available and are discussed in this chapter.

Altered gene expression in the Werner and Bloom syndromes is associated with sequences having G-quadruplex forming potential

The human Werner and Bloom syndromes (WS and BS) are caused by deficiencies in the WRN and BLM RecQ helicases, respectively. WRN, BLM and their Saccharomyces cerevisiae homologue Sgs1, are particularly active in vitro in unwinding G-quadruplex DNA (G4-DNA), a family of non-canonical nucleic acid structures formed by certain G-rich sequences. Recently, mRNA levels from loci containing potential G-quadruplex-forming sequences (PQS) were found to be preferentially altered in sgs1Δ mutants, suggesting that G4-DNA targeting by Sgs1 directly affects gene expression. Here, we extend these findings to human cells. Using microarrays to measure mRNAs obtained from human fibroblasts deficient for various RecQ family helicases, we observe significant associations between loci that are upregulated in WS or BS cells and loci that have PQS. No such PQS associations were observed for control expression datasets, however. Furthermore, upregulated genes in WS and BS showed no or dramatically reduced associations with sequences similar to PQS but that have considerably reduced potential to form intramolecular G4-DNA. These findings indicate that, like Sgs1, WRN and BLM can regulate transcription globally by targeting G4-DNA.

Extracellular growth factors and mitogens cooperate to drive mitochondrial biogenesis

Cells generate new organelles when stimulated by extracellular factors to grow and divide; however, little is known about how growth and mitogenic signalling pathways regulate organelle biogenesis. Using mitochondria as a model organelle, we have investigated this problem in primary Schwann cells, for which distinct factors act solely as mitogens (neuregulin) or as promoters of cell growth (insulin-like growth factor 1; IGF1). We find that neuregulin and IGF1 act synergistically to increase mitochondrial biogenesis and mitochondrial DNA replication, resulting in increased mitochondrial density in these cells. Moreover, constitutive oncogenic Ras signalling results in a further increase in mitochondrial density. This synergistic effect is seen at the global transcriptional level, requires both the ERK and phosphoinositide 3-kinase (PI3K) signalling pathways and is mediated by the transcription factor ERRalpha. Interestingly, the effect is independent of Akt-TOR signalling, a major regulator of cell growth in these cells. This separation of the pathways that drive mitochondrial biogenesis and cell growth provides a mechanism for the modulation of mitochondrial density according to the metabolic requirements of the cell.

Transcriptional profiling at different sites in lungs of pigs during acute bacterial respiratory infection

The local transcriptional response was studied in different locations of lungs from pigs experimentally infected with the respiratory pathogen Actinobacillus pleuropneumoniae serotype 5B, using porcine cDNA microarrays. This infection gives rise to well-demarcated infection loci in the lung, characterized by necrotic and haemorrhagic lesions. Lung tissue was sampled from necrotic areas, from visually unaffected areas and from areas bordering on necrotic areas. Expression pattern of these areas from infected pigs was compared to healthy lung tissue from un-infected pigs. Transcription of selected genes important in the innate defence response were further analysed by quantitative real-time reverse-transcriptase PCR. A clear correlation was observed between the number of differentially expressed genes as well as the magnitude of their induction and the sampling location in the infected lung, with the highest number of differentially expressed genes, and the most highly induced genes found in necrotic areas. Interestingly, a group of differentially regulated genes was represented in all three areas, comprising genes encoding cytokines, acute phase proteins, and factors related to regulation of apoptosis and the complement system. Interferon-γ was downregulated in both necrotic and bordering areas. Evidence of neutrophil recruitment was seen by the up-regulation of chemotactic factors for neutrophils. In conclusion, we found subsets of genes expressed at different levels in the three selected areas of the infected lung as compared to the control group. Thus it is demonstrated that an infection with clearly defined infected loci leads to a rapid disseminated intra-organ response in neighbouring seemingly unaffected tissue areas of the infected organ. Within the lung, we found a clear division of induced genes as, in unaffected areas a large part of differently expressed genes were involved in systemic reactions to infections, while differently expressed genes in necrotic areas were mainly concerned with homeostasis regulation.

Graph ranking for exploratory gene data analysis

Background

Microarray technology has made it possible to simultaneously monitor the expression levels of thousands of genes in a single experiment. However, the large number of genes greatly increases the challenges of analyzing, comprehending and interpreting the resulting mass of data. Selecting a subset of important genes is inevitable to address the challenge. Gene selection has been investigated extensively over the last decade. Most selection procedures, however, are not sufficient for accurate inference of underlying biology, because biological significance does not necessarily have to be statistically significant. Additional biological knowledge needs to be integrated into the gene selection procedure.

Results

We propose a general framework for gene ranking. We construct a bipartite graph from the Gene Ontology (GO) and gene expression data. The graph describes the relationship between genes and their associated molecular functions. Under a species condition, edge weights of the graph are assigned to be gene expression level. Such a graph provides a mathematical means to represent both species-independent and species-dependent biological information. We also develop a new ranking algorithm to analyze the weighted graph via a kernelized spatial depth (KSD) approach. Consequently, the importance of gene and molecular function can be simultaneously ranked by a real-valued measure, KSD, which incorporates the global and local structure of the graph. Over-expressed and under-regulated genes also can be separately ranked.

Conclusion

The gene-function bigraph integrates molecular function annotations into gene expression data. The relevance of genes is described in the graph (through a common function). The proposed method provides an exploratory framework for gene data analysis.

Low oxygen levels as a trigger for enhancement of respiratory metabolism in Saccharomyces cerevisiae

Background

The industrially important yeast Saccharomyces cerevisiae is able to grow both in the presence and absence of oxygen. However, the regulation of its metabolism in conditions of intermediate oxygen availability is not well characterised. We assessed the effect of oxygen provision on the transcriptome and proteome of S. cerevisiae in glucose-limited chemostat cultivations in anaerobic and aerobic conditions, and with three intermediate (0.5, 1.0 and 2.8% oxygen) levels of oxygen in the feed gas.

Results

The main differences in the transcriptome were observed in the comparison of fully aerobic, intermediate oxygen and anaerobic conditions, while the transcriptome was generally unchanged in conditions receiving different intermediate levels (0.5, 1.0 or 2.8% O₂) of oxygen in the feed gas. Comparison of the transcriptome and proteome data suggested post-transcriptional regulation was important, especially in 0.5% oxygen. In the conditions of intermediate oxygen, the genes encoding enzymes of the respiratory pathway were more highly expressed than in either aerobic or anaerobic conditions. A similar trend was also seen in the proteome and in enzyme activities of the TCA cycle. Further, genes encoding proteins of the mitochondrial translation machinery were present at higher levels in all oxygen-limited and anaerobic conditions, compared to fully aerobic conditions.

Conclusion

Global upregulation of genes encoding components of the respiratory pathway under conditions of intermediate oxygen suggested a regulatory mechanism to control these genes as a response to the need of more efficient energy production. Further, cells grown in three different intermediate oxygen levels were highly similar at the level of transcription, while they differed at the proteome level, suggesting post-transcriptional mechanisms leading to distinct physiological modes of respiro-fermentative metabolism.

Some molecular effectors of antidepressant action of quetiapine revealed by DNA microarray in the frontal cortex of anhedonic rats

OBJECTIVES AND METHODS

We have previously demonstrated that quetiapine (QTP) had antidepressant-like action by using the chronic mild stress (CMS) paradigm, an animal model of human depression. The aim of this study was to investigate the molecular mechanism(s) of QTP antidepressant effect by coupling the CMS protocol with Affymetrix microarray technology to screen the entire rat genome for gene changes in the frontal cortex.

RESULTS

The genes regulated by the administration of CMS whose transcription was reversed by chronic QTP treatment (2 mg/kg/day) were 42 (23 upregulated and 19 downregulated). The transcripts that showed no significant altered expression levels in anhedonic rats but were regulated by the administration of QTP were 19 (nine upregulated and 10 downregulated). On the whole, the action of QTP prevented the stress-induced impairment of some processes involved in central nervous system development or having a crucial role for viability of neural cells and cell-cell communications, like regulation of signal transduction, inorganic cation transport, membrane organization, and neurite morphogenesis. For 11 genes (Ptgs2, Gad1, Plcb1, Camk2a, Homer1, Senp2, Junb, Nfib, Hes5, Capon, and Marcks), significant differential expressions were confirmed by real-time reverse-transcriptase polymerase chain reaction.

CONCLUSION

We have shown that chronic QTP treatment prevented anhedonia and reversed, at least in part, the changes of gene expression induced by CMS in the rat frontal cortex. We have also identified and confirmed by two different methods that 11 genes, representing molecular targets of QTP, are presumably the effectors of its clinical efficacy.

Reciprocal regulation of microRNA and mRNA profiles in neuronal development and synapse formation

Background

Synapse formation and the development of neural networks are known to be controlled by a coordinated program of mRNA synthesis. microRNAs are now recognized to be important regulators of mRNA translation and stability in a wide variety of organisms. While specific microRNAs are known to be involved in neural development, the extent to which global microRNA and mRNA profiles are coordinately regulated in neural development is unknown.

Results

We examined mouse primary neuronal cultures, analyzing microRNA and mRNA expression. Three main developmental patterns of microRNA expression were observed: steady-state levels, up-regulated and down-regulated. Co-expressed microRNAs were found to have related target recognition sites and to be encoded in distinct genomic locations. A number of 43 differentially expressed miRNAs were located in five genomic clusters. Their predicted mRNA targets show reciprocal levels of expression. We identified a set of reciprocally expressed microRNAs that target mRNAs encoding postsynaptic density proteins and high-level steady-state microRNAs that target non-neuronal low-level expressed mRNAs.

Conclusion

We characterized hundreds of miRNAs in neuronal culture development and identified three major modes of miRNA expression. We predict these miRNAs to regulate reciprocally expressed protein coding genes, including many genes involved in synaptogenesis. The identification of miRNAs that target mRNAs during synaptogenesis indicates a new level of regulation of the synapse.

Arsenicals produce stable progressive changes in DNA methylation patterns that are linked to malignant transformation of immortalized urothelial cells

Aberrant DNA methylation participates in carcinogenesis and is a molecular hallmark of a tumor cell. Tumor cells generally exhibit a redistribution of DNA methylation resulting in global hypomethylation with regional hypermethylation; however, the speed in which these changes emerge has not been fully elucidated and may depend on the temporal location of the cell in the path from normal, finite lifespan to malignant transformation. We used a model of arsenical-induced malignant transformation of immortalized human urothelial cells and DNA methylation microarrays to examine the extent and temporal nature of changes in DNA methylation that occur during the transition from immortal to malignantly transformed. Our data presented herein suggest that during arsenical-induced malignant transformation, aberrant DNA methylation occurs non-randomly, progresses gradually at hundreds of gene promoters, and alters expression of the associated gene, and these changes are coincident with the acquisition of malignant properties, such as anchorage independent growth and tumor formation in immunocompromised mice. The DNA methylation changes appear stable, since malignantly transformed cells removed from the transforming arsenical exhibited no reversion in DNA methylation levels, associated gene expression, or malignant phenotype. These data suggest that arsenicals act as epimutagens and directly link their ability to induce malignant transformation to their actions on the epigenome.

Statistical methods for gene set co-expression analysis

Motivation: The power of a microarray experiment derives from the identification of genes differentially regulated across biological conditions. To date, differential regulation is most often taken to mean differential expression, and a number of useful methods for identifying differentially expressed (DE) genes or gene sets are available. However, such methods are not able to identify many relevant classes of differentially regulated genes. One important example concerns differentially co-expressed (DC) genes.

Results: We propose an approach, gene set co-expression analysis (GSCA), to identify DC gene sets. The GSCA approach provides a false discovery rate controlled list of interesting gene sets, does not require that genes be highly correlated in at least one biological condition and is readily applied to data from individual or multiple experiments, as we demonstrate using data from studies of lung cancer and diabetes.

Availability: The GSCA approach is implemented in R and available at www.biostat.wisc.edu/∼kendzior/GSCA/.

Contact:kendzior@biostat.wisc.edu

Supplementary information:Supplementary data are available at Bioinformatics online.

Effects of atherogenic diet on hepatic gene expression across mouse strains

Diets high in fat and cholesterol are associated with increased obesity and metabolic disease in mice and humans. To study the molecular basis of the metabolic response to dietary fat, 10 inbred strains of mice were fed atherogenic high-fat and control low-fat diets. Liver gene expression and whole animal phenotypes were measured and analyzed in both sexes. The effects of diet, strain, and sex on gene expression were determined irrespective of complex processes, such as feedback mechanisms, that could have mediated the genomic responses. Global gene expression analyses demonstrated that animals of the same strain and sex have similar transcriptional profiles on a low-fat diet, but strains may show considerable variability in response to high-fat diet. Functional profiling indicated that high-fat feeding induced genes in the immune response, indicating liver damage, and repressed cholesterol biosynthesis. The physiological significance of the transcriptional changes was confirmed by a correlation analysis of transcript levels with whole animal phenotypes. The results found here were used to confirm a previously identified quantitative trait locus on chromosome 17 identified in males fed a high-fat diet in two crosses, PERA x DBA/2 and PERA x I/Ln. The gene expression data and phenotype data have been made publicly available as an online tool for exploring the effects of atherogenic diet in inbred mouse strains (http://cgd-array.jax.org/DietStrainSurvey).

Computational identification of hepatitis C virus associated microRNA-mRNA regulatory modules in human livers

BACKGROUND

Hepatitis C virus (HCV) is a major cause of chronic liver disease by infecting over 170 million people worldwide. Recent studies have shown that microRNAs (miRNAs), a class of small non-coding regulatory RNAs, are involved in the regulation of HCV infection, but their functions have not been systematically studied. We propose an integrative strategy for identifying the miRNA-mRNA regulatory modules that are associated with HCV infection. This strategy combines paired expression profiles of miRNAs and mRNAs and computational target predictions. A miRNA-mRNA regulatory module consists of a set of miRNAs and their targets, in which the miRNAs are predicted to coordinately regulate the level of the target mRNA.

RESULTS

We simultaneously profiled the expression of cellular miRNAs and mRNAs across 30 HCV positive or negative human liver biopsy samples using microarray technology. We constructed a miRNA-mRNA regulatory network, and using a graph theoretical approach, identified 38 miRNA-mRNA regulatory modules in the network that were associated with HCV infection. We evaluated the direct miRNA regulation of the mRNA levels of targets in regulatory modules using previously published miRNA transfection data. We analyzed the functional roles of individual modules at the systems level by integrating a large-scale protein interaction network. We found that various biological processes, including some HCV infection related canonical pathways, were regulated at the miRNA level during HCV infection.

CONCLUSION

Our regulatory modules provide a framework for future experimental analyses. This report demonstrates the utility of our approach to obtain new insights into post-transcriptional gene regulation at the miRNA level in complex human diseases.

Salmonella Typhimurium type III secretion effectors stimulate innate immune responses in cultured epithelial cells

Recognition of conserved bacterial products by innate immune receptors leads to inflammatory responses that control pathogen spread but that can also result in pathology. Intestinal epithelial cells are exposed to bacterial products and therefore must prevent signaling through innate immune receptors to avoid pathology. However, enteric pathogens are able to stimulate intestinal inflammation. We show here that the enteric pathogen Salmonella Typhimurium can stimulate innate immune responses in cultured epithelial cells by mechanisms that do not involve receptors of the innate immune system. Instead, S. Typhimurium stimulates these responses by delivering through its type III secretion system the bacterial effector proteins SopE, SopE2, and SopB, which in a redundant fashion stimulate Rho-family GTPases leading to the activation of mitogen-activated protein (MAP) kinase and NF-κB signaling. These observations have implications for the understanding of the mechanisms by which Salmonella Typhimurium induces intestinal inflammation as well as other intestinal inflammatory pathologies.

Salmonella Typhimurium is one of the most common causes of food-borne illness in the United States and a major cause of diarrheal diseases in developing countries. This pathogen induces diarrhea by stimulating inflammation in the intestinal tract. This study shows that S. Typhimurium delivers molecules into epithelial cells with the capacity to stimulate the production of pro-inflammatory substances. This mechanism may help the pathogen to initiate the inflammatory response in the intestinal epithelium. This study provides insight into the mechanisms by which Salmonella Typhimurium causes diarrhea.

R/BHC: fast Bayesian hierarchical clustering for microarray data

BACKGROUND

Although the use of clustering methods has rapidly become one of the standard computational approaches in the literature of microarray gene expression data analysis, little attention has been paid to uncertainty in the results obtained.

RESULTS

We present an R/Bioconductor port of a fast novel algorithm for Bayesian agglomerative hierarchical clustering and demonstrate its use in clustering gene expression microarray data. The method performs bottom-up hierarchical clustering, using a Dirichlet Process (infinite mixture) to model uncertainty in the data and Bayesian model selection to decide at each step which clusters to merge.

CONCLUSION

Biologically plausible results are presented from a well studied data set: expression profiles of A. thaliana subjected to a variety of biotic and abiotic stresses. Our method avoids several limitations of traditional methods, for example how many clusters there should be and how to choose a principled distance metric.

From disease ontology to disease-ontology lite: statistical methods to adapt a general-purpose ontology for the test of gene-ontology associations

Subjective methods have been reported to adapt a general-purpose ontology for a specific application. For example, Gene Ontology (GO) Slim was created from GO to generate a highly aggregated report of the human-genome annotation. We propose statistical methods to adapt the general purpose, OBO Foundry Disease Ontology (DO) for the identification of gene-disease associations. Thus, we need a simplified definition of disease categories derived from implicated genes. On the basis of the assumption that the DO terms having similar associated genes are closely related, we group the DO terms based on the similarity of gene-to-DO mapping profiles. Two types of binary distance metrics are defined to measure the overall and subset similarity between DO terms. A compactness-scalable fuzzy clustering method is then applied to group similar DO terms. To reduce false clustering, the semantic similarities between DO terms are also used to constrain clustering results. As such, the DO terms are aggregated and the redundant DO terms are largely removed. Using these methods, we constructed a simplified vocabulary list from the DO called Disease Ontology Lite (DOLite). We demonstrated that DOLite results in more interpretable results than DO for gene-disease association tests. The resultant DOLite has been used in the Functional Disease Ontology (FunDO) Web application at http://www.projects.bioinformatics.northwestern.edu/fundo.

Contact:s-lin2@northwestern.edu

Meta-analysis of gene expression in human pancreatic islets after in vitro expansion

Pancreatic islet transplantation as a potential cure for type 1 diabetes (T1D) cannot be scaled up due to a scarcity of human pancreas donors. In vitro expansion of beta-cells from mature human pancreatic islets provides an alternative source of insulin-producing cells. The exact nature of the expanded cells produced by diverse expansion protocols and their potential for differentiation into functional beta-cells remain elusive. We performed a large-scale meta-analysis of gene expression in human pancreatic islet cells, which were processed using three different previously described protocols for expansion and for which redifferentiation was attempted. All three expansion protocols induced dramatic changes in the expression profiles of pancreatic islets; many of these changes are shared among the three protocols. Attempts at redifferentiation of expanded cells induce a limited number of gene expression changes. Nevertheless, these fail to restore a pancreatic islet-like gene expression pattern. Comparison with a collection of public microarray datasets confirmed that expanded cells are highly comparable to mesenchymal stem cells. Genes induced in expanded cells are also enriched for targets of transcription factors important for pluripotency induction. The present data increase our understanding of the active pathways in expanded and redifferentiated islets. Knowledge of the mesenchymal stem cell potential may help development of drug therapeutics to restore beta-cell mass in T1D patients.

Importance of randomization in microarray experimental designs with Illumina platforms

Measurements of gene expression from microarray experiments are highly dependent on experimental design. Systematic noise can be introduced into the data at numerous steps. On Illumina BeadChips, multiple samples are assayed in an ordered series of arrays. Two experiments were performed using the same samples but different hybridization designs. An experiment confounding genotype with BeadChip and treatment with array position was compared to another experiment in which these factors were randomized to BeadChip and array position. An ordinal effect of array position on intensity values was observed in both experiments. We demonstrate that there is increased rate of false-positive results in the confounded design and that attempts to correct for confounded effects by statistical modeling reduce power of detection for true differential expression. Simple analysis models without post hoc corrections provide the best results possible for a given experimental design. Normalization improved differential expression testing in both experiments but randomization was the most important factor for establishing accurate results. We conclude that lack of randomization cannot be corrected by normalization or by analytical methods. Proper randomization is essential for successful microarray experiments.

Methods for interpreting lists of affected genes obtained in a DNA microarray experiment

BACKGROUND

The aim of this paper was to describe and compare the methods used and the results obtained by the participants in a joint EADGENE (European Animal Disease Genomic Network of Excellence) and SABRE (Cutting Edge Genomics for Sustainable Animal Breeding) workshop focusing on post analysis of microarray data. The participating groups were provided with identical lists of microarray probes, including test statistics for three different contrasts, and the normalised log-ratios for each array, to be used as the starting point for interpreting the affected probes. The data originated from a microarray experiment conducted to study the host reactions in broilers occurring shortly after a secondary challenge with either a homologous or heterologous species of Eimeria.

RESULTS

Several conceptually different analytical approaches, using both commercial and public available software, were applied by the participating groups. The following tools were used: Ingenuity Pathway Analysis, MAPPFinder, LIMMA, GOstats, GOEAST, GOTM, Globaltest, TopGO, ArrayUnlock, Pathway Studio, GIST and AnnotationDbi. The main focus of the approaches was to utilise the relation between probes/genes and their gene ontology and pathways to interpret the affected probes/genes. The lack of a well-annotated chicken genome did though limit the possibilities to fully explore the tools. The main results from these analyses showed that the biological interpretation is highly dependent on the statistical method used but that some common biological conclusions could be reached.

CONCLUSION

It is highly recommended to test different analytical methods on the same data set and compare the results to obtain a reliable biological interpretation of the affected genes in a DNA microarray experiment.

Comparison of three microarray probe annotation pipelines: differences in strategies and their effect on downstream analysis

BACKGROUND

Reliable annotation linking oligonucleotide probes to target genes is essential for functional biological analysis of microarray experiments. We used the IMAD, OligoRAP and sigReannot pipelines to update the annotation for the ARK-Genomics Chicken 20 K array as part of a joined EADGENE/SABRE workshop. In this manuscript we compare their annotation strategies and results. Furthermore, we analyse the effect of differences in updated annotation on functional analysis for an experiment involving Eimeria infected chickens and finally we propose guidelines for optimal annotation strategies.

RESULTS

IMAD, OligoRAP and sigReannot update both annotation and estimated target specificity. The 3 pipelines can assign oligos to target specificity categories although with varying degrees of resolution. Target specificity is judged based on the amount and type of oligo versus target-gene alignments (hits), which are determined by filter thresholds that users can adjust based on their experimental conditions. Linking oligos to annotation on the other hand is based on rigid rules, which differ between pipelines.For 52.7% of the oligos from a subset selected for in depth comparison all pipelines linked to one or more Ensembl genes with consensus on 44.0%. In 31.0% of the cases none of the pipelines could assign an Ensembl gene to an oligo and for the remaining 16.3% the coverage differed between pipelines. Differences in updated annotation were mainly due to different thresholds for hybridisation potential filtering of oligo versus target-gene alignments and different policies for expanding annotation using indirect links. The differences in updated annotation packages had a significant effect on GO term enrichment analysis with consensus on only 67.2% of the enriched terms.

CONCLUSION

In addition to flexible thresholds to determine target specificity, annotation tools should provide metadata describing the relationships between oligos and the annotation assigned to them. These relationships can then be used to judge the varying degrees of reliability allowing users to fine-tune the balance between reliability and coverage. This is important as it can have a significant effect on functional microarray analysis as exemplified by the lack of consensus on almost one third of the terms found with GO term enrichment analysis based on updated IMAD, OligoRAP or sigReannot annotation.

Microarray data mining using Bioconductor packages

Background

This paper describes the results of a Gene Ontology (GO) term enrichment analysis of chicken microarray data using the Bioconductor packages. By checking the enriched GO terms in three contrasts, MM8-PM8, MM8-MA8, and MM8-MM24, of the provided microarray data during this workshop, this analysis aimed to investigate the host reactions in chickens occurring shortly after a secondary challenge with either a homologous or heterologous species of Eimeria. The results of GO enrichment analysis using GO terms annotated to chicken genes and GO terms annotated to chicken-human orthologous genes were also compared. Furthermore, a locally adaptive statistical procedure (LAP) was performed to test differentially expressed chromosomal regions, rather than individual genes, in the chicken genome after Eimeria challenge.

Results

GO enrichment analysis identified significant (raw p-value < 0.05) GO terms for all three contrasts included in the analysis. Some of the GO terms linked to, generally, primary immune responses or secondary immune responses indicating the GO enrichment analysis is a useful approach to analyze microarray data. The comparisons of GO enrichment results using chicken gene information and chicken-human orthologous gene information showed more refined GO terms related to immune responses when using chicken-human orthologous gene information, this suggests that using chicken-human orthologous gene information has higher power to detect significant GO terms with more refined functionality. Furthermore, three chromosome regions were identified to be significantly up-regulated in contrast MM8-PM8 (q-value < 0.01).

Conclusion

Overall, this paper describes a practical approach to analyze microarray data in farm animals where the genome information is still incomplete. For farm animals, such as chicken, with currently limited gene annotation, borrowing gene annotation information from orthologous genes in well-annotated species, such as human, will help improve the pathway analysis results substantially. Furthermore, LAP analysis approach is a relatively new and very useful way to be applied in microarray analysis.

Characterisation of a novel NR4A2 mutation in Parkinson's disease brain

OBJECTIVE

We performed a mutation screen of NR4A2 (also known as NURR1) in 409 Parkinson's disease (PD) patients. We identified a novel single base substitution in the 5'UTR of the NR4A2 (also known as NURR1) gene (c.-309C>T).

RESULTS

We have performed expression studies in neuronal cell lines showing that the c.-309C>T mutation reduces NR4A2 mRNA expression in vitro. We have confirmed this finding in vivo by performing allele specific real-time PCR from brain tissue harbouring the 309C>T mutation and show a 3.48+/-1.62 fold reduction in mRNA expression of the mutant allele compared to wild-type. In addition we have undertaken genome wide expression analysis of the mutant NR4A2 brain and shown underexpressed genes were significantly enriched for gene ontology categories in nervous system development and synaptic transmission and overexpressed genes were enriched for unfolded protein response and morphogenesis. Lastly we have shown that the c.-309C>T mutation abrogates the protective effect of wild-type NR4A2 against apoptopic stress.

CONCLUSIONS

Our findings indicate the c.-309C>T mutation reduces NR4A2 expression resulting in the downregulation of genes involved in the development and maintenance of the nervous system and synaptic transmission. These downregulated pathways contained genes known to be transactivated by NR4A2 and were not disrupted in idiopathic PD brain suggesting causality of the mutation.

Host cell transcriptional profiling during malaria liver stage infection reveals a coordinated and sequential set of biological events

Background

Plasmodium sporozoites migrate to the liver where they traverse several hepatocytes before invading the one inside which they will develop and multiply into thousands of merozoites. Although this constitutes an essential step of malaria infection, the requirements of Plasmodium parasites in liver cells and how they use the host cell for their own survival and development are poorly understood.

Results

To gain new insights into the molecular host-parasite interactions that take place during malaria liver infection, we have used high-throughput microarray technology to determine the transcriptional profile of P. berghei-infected hepatoma cells. The data analysis shows differential expression patterns for 1064 host genes starting at 6 h and up to 24 h post infection, with the largest proportion correlating specifically with the early stages of the infection process. A considerable proportion of those genes were also found to be modulated in liver cells collected from P. yoelii-infected mice 24 and 40 h after infection, strengthening the data obtained with the in vitro model and highlighting genes and pathways involved in the host response to rodent Plasmodium parasites.

Conclusion

Our data reveal that host cell infection by Plasmodium sporozoites leads to a coordinated and sequential set of biological events, ranging from the initial stage of stress response up to the engagement of host metabolic processes and the maintenance of cell viability throughout infection.

Stepwise DNA methylation changes are linked to escape from defined proliferation barriers and mammary epithelial cell immortalization

The timing and progression of DNA methylation changes during carcinogenesis are not completely understood. To develop a timeline of aberrant DNA methylation events during malignant transformation, we analyzed genome-wide DNA methylation patterns in an isogenic human mammary epithelial cell (HMEC) culture model of transformation. To acquire immortality and malignancy, the cultured finite lifespan HMEC must overcome two distinct proliferation barriers. The first barrier, stasis, is mediated by the retinoblastoma protein and can be overcome by loss of p16(INK4A) expression. HMEC that escape stasis and continue to proliferate become genomically unstable before encountering a second more stringent proliferation barrier, telomere dysfunction due to telomere attrition. Rare cells that acquire telomerase expression may escape this barrier, become immortal, and develop further malignant properties. Our analysis of HMEC transitioning from finite lifespan to malignantly transformed showed that aberrant DNA methylation changes occur in a stepwise fashion early in the transformation process. The first aberrant DNA methylation step coincides with overcoming stasis, and results in few to hundreds of changes, depending on how stasis was overcome. A second step coincides with immortalization and results in hundreds of additional DNA methylation changes regardless of the immortalization pathway. A majority of these DNA methylation changes are also found in malignant breast cancer cells. These results show that large-scale epigenetic remodeling occurs in the earliest steps of mammary carcinogenesis, temporally links DNA methylation changes and overcoming cellular proliferation barriers, and provides a bank of potential epigenetic biomarkers that may prove useful in breast cancer risk assessment.

PTEN deficiency in a luminal ErbB-2 mouse model results in dramatic acceleration of mammary tumorigenesis and metastasis

Overexpression and/or amplification of the ErbB-2 oncogene as well as inactivation of the PTEN tumor suppressor are two important genetic events in human breast carcinogenesis. To address the biological impact of conditional inactivation of PTEN on ErbB-2-induced mammary tumorigenesis, we generated a novel transgenic mouse model that utilizes the murine mammary tumor virus (MMTV) promoter to directly couple expression of activated ErbB-2 and Cre recombinase to the same mammary epithelial cell (MMTV-NIC). Disruption of PTEN in the mammary epithelium of the MMTV-NIC model system dramatically accelerated the formation of multifocal and highly metastatic mammary tumors, which exhibited homogenous pathology. PTEN-deficient/NIC-induced tumorigenesis was associated with an increase in angiogenesis. Moreover, inactivation of PTEN in the MMTV-NIC mouse model resulted in hyperactivation of the phosphatidylinositol 3'-kinase/Akt signaling pathway. However, like the parental strain, tumors obtained from PTEN-deficient/NIC mice displayed histopathological and molecular features of the luminal subtype of primary human breast cancer. Taken together, our findings provide important implications in understanding the molecular determinants of mammary tumorigenesis driven by PTEN deficiency and ErbB-2 activation and could provide a valuable tool for testing the efficacy of therapeutic strategies that target these critical signaling pathways.

Transcriptional profiling of cattle infected with Trypanosoma congolense highlights gene expression signatures underlying trypanotolerance and trypanosusceptibility

Background

African animal trypanosomiasis (AAT) caused by tsetse fly-transmitted protozoa of the genus Trypanosoma is a major constraint on livestock and agricultural production in Africa and is among the top ten global cattle diseases impacting on the poor. Here we show that a functional genomics approach can be used to identify temporal changes in host peripheral blood mononuclear cell (PBMC) gene expression due to disease progression. We also show that major gene expression differences exist between cattle from trypanotolerant and trypanosusceptible breeds. Using bovine long oligonucleotide microarrays and real time quantitative reverse transcription PCR (qRT-PCR) validation we analysed PBMC gene expression in naïve trypanotolerant and trypanosusceptible cattle experimentally challenged with Trypanosoma congolense across a 34-day infection time course.

Results

Trypanotolerant N'Dama cattle displayed a rapid and distinct transcriptional response to infection, with a ten-fold higher number of genes differentially expressed at day 14 post-infection compared to trypanosusceptible Boran cattle. These analyses identified coordinated temporal gene expression changes for both breeds in response to trypanosome infection. In addition, a panel of genes were identified that showed pronounced differences in gene expression between the two breeds, which may underlie the phenomena of trypanotolerance and trypanosusceptibility. Gene ontology (GO) analysis demonstrate that the products of these genes may contribute to increased mitochondrial mRNA translational efficiency, a more pronounced B cell response, an elevated activation status and a heightened response to stress in trypanotolerant cattle.

Conclusion

This study has revealed an extensive and diverse range of cellular processes that are altered temporally in response to trypanosome infection in African cattle. Results indicate that the trypanotolerant N'Dama cattle respond more rapidly and with a greater magnitude to infection compared to the trypanosusceptible Boran cattle. Specifically, a subset of the genes analyzed by real time qRT-PCR, which display significant breed differences, could collectively contribute to the trypanotolerance trait in N'Dama.

Chromosome Y variants from different inbred mouse strains are linked to differences in the morphologic and molecular responses of cardiac cells to postpubertal testosterone

BACKGROUND

We have reported previously that when chromosome Y (chrY) from the mouse strain C57BL/6J (ChrYC57) was substituted for that of A/J mice (ChrYA), cardiomyocytes from the resulting "chromosome substitution" C57BL/6J-chrYA strain were smaller than that of their C57BL/6J counterparts. In reverse, when chrYA from A/J mice was substituted for that of chrYC57, cardiomyocytes from the resulting A/J-chrYC57 strain were larger than in their A/J counterparts. We further used these strains to test whether: 1) the origin of chrY could also be linked to differences in the profile of gene expression in the hearts of adult male mice, and 2) post-pubertal testosterone could play a role in the differential morphologic and/or molecular effects of chrYC57 and chrYA.

RESULTS

The increased size of cardiomyocytes from adult male C57BL/6J mice compared to C57BL/6J-chrYA resulted from the absence of hypertrophic effects of post-pubertal testosterone on cells from the latter strain. However, gene profiling revealed that the latter effect could not be explained on the basis of an insensitivity of cells from C57BL/6J-chrYA to androgens, since even more cardiac genes were affected by post-pubertal testosterone in C57BL/6J-chrYA hearts than in C57BL/6J. By testing for interaction between the effects of surgery and strain, we identified 249 "interaction genes" whose expression was affected by post-pubertal testosterone differentially according to the genetic origin of chrY. These interaction genes were found to be enriched within a limited number of signaling pathways, including: 1) p53 signaling, which comprises the interacting genes Ccnd1, Pten and Cdkn1a that are also potential co-regulators of the androgen receptors, and 2) circadian rhythm, which comprises Arntl/Bmal1, which may in turn regulate cell growth via the control of Cdkn1a.

CONCLUSION

Although post-pubertal testosterone increased the size of cardiomyocytes from male C56BL/6J mice but not that from their C57BL/6J-chrYA counterparts, it affected gene expression in the hearts from both strains. However, several cardiac genes responded to post-pubertal testosterone in a strict strain-selective manner, which provides possible mechanisms explaining how chrY may, in part via interference with androgen regulatory events, be linked to morphologic differences of cardiac cells of adult male mice.

Gene expression profiling differentiates autism case-controls and phenotypic variants of autism spectrum disorders: evidence for circadian rhythm dysfunction in severe autism

Autism spectrum disorders (ASD) are neurodevelopmental disorders characterized by delayed/abnormal language development, deficits in social interaction, repetitive behaviors and restricted interests. The heterogeneity in clinical presentation of ASD, likely due to different etiologies, complicates genetic/biological analyses of these disorders. DNA microarray analyses were conducted on 116 lymphoblastoid cell lines (LCL) from individuals with idiopathic autism who are divided into three phenotypic subgroups according to severity scores from the commonly used Autism Diagnostic Interview-Revised questionnaire and age-matched, nonautistic controls. Statistical analyses of gene expression data from control LCL against that of LCL from ASD probands identify genes for which expression levels are either quantitatively or qualitatively associated with phenotypic severity. Comparison of the significant differentially expressed genes from each subgroup relative to the control group reveals differentially expressed genes unique to each subgroup as well as genes in common across subgroups. Among the findings unique to the most severely affected ASD group are 15 genes that regulate circadian rhythm, which has been shown to have multiple effects on neurological as well as metabolic functions commonly dysregulated in autism. Among the genes common to all three subgroups of ASD are 20 novel genes mostly in putative noncoding regions, which appear to associate with androgen sensitivity and which may underlie the strong 4:1 bias toward affected males.

CORNA: testing gene lists for regulation by microRNAs

MOTIVATION

With the increasing use of post-genomics techniques to examine a wide variety of biological systems in laboratories throughout the world, scientists are often presented with lists of genes that they must make sense of. A consistently challenging problem is that of defining co-regulated genes within those gene lists. In recent years, microRNAs have emerged as a mechanism for regulating several cellular processes. In this article, we report on how gene lists and microRNA targets data may be integrated to test for significant associations between gene lists and microRNAs.

RESULTS

We discuss CORNA, a package written in R and released under the GNU GPL, which allows users to test gene lists for significant microRNA-target associations using one of three separate statistical tests, to link microRNA targets to functional annotation and to visualize quantitative data associated with those data.

AVAILABILITY

CORNA is available as an R package from http://corna.sf.net.

Comparative proteomic phenotyping of cell lines and primary cells to assess preservation of cell type-specific functions

Biological experiments are most often performed with immortalized cell lines because they are readily available and can be expanded without limitation. However, cell lines may differ from the in vivo situation in important aspects. Here we introduce a straightforward methodology to compare cell lines to their cognate primary cells and to derive a comparative functional phenotype. We used SILAC (stable isotope labeling by amino acids in cell culture) for quantitative, mass spectrometry-based comparison of the hepatoma cell line Hepa1-6 with primary hepatocytes. The resulting quantitative proteome of 4,063 proteins had an asymmetric distribution, with many proteins down-regulated in the cell line. Bioinformatic analysis of the quantitative proteomics phenotypes revealed that Hepa1-6 cells were deficient in mitochondria, reflecting re-arrangement of metabolic pathways, drastically up-regulate cell cycle-associated functions and largely shut down drug metabolizing enzymes characteristic for the liver. This quantitative knowledge of changes provides an important basis to adapt cell lines to more closely resemble physiological conditions.