DNA microarray data integration by ortholog gene analysis reveals potential molecular mechanisms of estrogen-dependent growth of human uterine fibroids

BACKGROUND: Uterine fibroids or leiomyoma are a common benign smooth muscle tumor. The tumor growth is well known to be estrogen-dependent. However, the molecular mechanisms of its estrogen-dependency is not well understood. METHODS: Differentially expressed genes in human uterine fibroids were either retrieved from published papers or from our own statistical analysis of downloaded array data. Probes for the same genes on different Affymetrix chips were mapped based on probe comparison information provided by Affymetrix. Genes identified by two or three array studies were submitted for ortholog analysis. Human and rat ortholog genes were identified by using ortholog gene databases, HomoloGene and TOGA and were confirmed by synteny analysis with MultiContigView tool in the Ensembl genome browser. RESULTS: By integrated analysis of three recently published DNA microarray studies with human tissue, thirty-eight genes were found to be differentially expressed in the same direction in fibroid compared to adjacent uterine myometrium by at least two research groups. Among these genes, twelve with rat orthologs were identified as estrogen-regulated from our array study investigating uterine expression in ovariectomized rats treated with estrogen. Functional and pathway analyses of the twelve genes suggested multiple molecular mechanisms for estrogen-dependent cell survival and tumor growth. Firstly, estrogen increased expression of the anti-apoptotic PCP4 gene and suppressed the expression of growth inhibitory receptors PTGER3 and TGFBR2. Secondly, estrogen may antagonize PPARγ signaling, thought to inhibit fibroid growth and survival, at two points in the PPAR pathway: 1) through increased ANXA1 gene expression which can inhibit phospholipase A2 activity and in turn decrease arachidonic acid synthesis, and 2) by decreasing L-PGDS expression which would reduce synthesis of PGJ2, an endogenous ligand for PPARγ. Lastly, estrogen affects retinoic acid (RA) synthesis and mobilization by regulating expression of CRABP2 and ALDH1A1. RA has been shown to play a significant role in the development of uterine fibroids in an animal model. CONCLUSION: Integrated analysis of multiple array datasets revealed twelve human and rat ortholog genes that were differentially expressed in human uterine fibroids and transcriptionally responsive to estrogen in the rat uterus. Functional and pathway analysis of these genes suggest multiple potential molecular mechanisms for the poorly understood estrogen-dependent growth of uterine fibroids. Fully understanding the exact molecular interactions among these gene products requires further study to validate their roles in uterine fibroids. This work provides new avenues of study which could influence the future direction of therapeutic intervention for the disease

Reconstruction of metabolic pathways for the cattle genome

<p>Abstract</p> <p>Background</p> <p>Metabolic reconstruction of microbial, plant and animal genomes is a necessary step toward understanding the evolutionary origins of metabolism and species-specific adaptive traits. The aims of this study were to reconstruct conserved metabolic pathways in the cattle genome and to identify metabolic pathways with missing genes and proteins. The MetaCyc database and PathwayTools software suite were chosen for this work because they are widely used and easy to implement.</p> <p>Results</p> <p>An amalgamated cattle genome database was created using the NCBI and Ensembl cattle genome databases (based on build 3.1) as data sources. PathwayTools was used to create a cattle-specific pathway genome database, which was followed by comprehensive manual curation for the reconstruction of metabolic pathways. The curated database, CattleCyc 1.0, consists of 217 metabolic pathways. A total of 64 mammalian-specific metabolic pathways were modified from the reference pathways in MetaCyc, and two pathways previously identified but missing from MetaCyc were added. Comparative analysis of metabolic pathways revealed the absence of mammalian genes for 22 metabolic enzymes whose activity was reported in the literature. We also identified six human metabolic protein-coding genes for which the cattle ortholog is missing from the sequence assembly.</p> <p>Conclusion</p> <p>CattleCyc is a powerful tool for understanding the biology of ruminants and other cetartiodactyl species. In addition, the approach used to develop CattleCyc provides a framework for the metabolic reconstruction of other newly sequenced mammalian genomes. It is clear that metabolic pathway analysis strongly reflects the quality of the underlying genome annotations. Thus, having well-annotated genomes from many mammalian species hosted in BioCyc will facilitate the comparative analysis of metabolic pathways among different species and a systems approach to comparative physiology.</p

Combination of novel and public RNA-seq datasets to generate an mRNA expression atlas for the domestic chicken

Background: The domestic chicken (Gallus gallus) is widely used as a model in developmental biology and is also an important livestock species. We describe a novel approach to data integration to generate an mRNA expression atlas for the chicken spanning major tissue types and developmental stages, using a diverse range of publicly-archived RNA-seq datasets and new data derived from immune cells and tissues. Results: Randomly down-sampling RNA-seq datasets to a common depth and quantifying expression against a reference transcriptome using the mRNA quantitation tool Kallisto ensured that disparate datasets explored comparable transcriptomic space. The network analysis tool Graphia was used to extract clusters of co-expressed genes from the resulting expression atlas, many of which were tissue or cell-type restricted, contained transcription factors that have previously been implicated in their regulation, or were otherwise associated with biological processes, such as the cell cycle. The atlas provides a resource for the functional annotation of genes that currently have only a locus ID. We cross-referenced the RNA-seq atlas to a publicly available embryonic Cap Analysis of Gene Expression (CAGE) dataset to infer the developmental time course of organ systems, and to identify a signature of the expansion of tissue macrophage populations during development. Conclusion: Expression profiles obtained from public RNA-seq datasets - despite being generated by different laboratories using different methodologies - can be made comparable to each other. This meta-analytic approach to RNA-seq can be extended with new datasets from novel tissues, and is applicable to any species