RESUMO
Charting microRNA (miRNA) regulation across pathways is key to characterizing their function. Yet, no method currently exists that can quantify how miRNAs regulate multiple interconnected pathways or prioritize them for their ability to regulate coordinate transcriptional programs. Existing methods primarily infer one-to-one relationships between miRNAs and pathways using differentially expressed genes. We introduce PanomiR, an in silico framework for studying the interplay of miRNAs and disease functions. PanomiR integrates gene expression, mRNA-miRNA interactions and known biological pathways to reveal coordinated multi-pathway targeting by miRNAs. PanomiR utilizes pathway-activity profiling approaches, a pathway co-expression network and network clustering algorithms to prioritize miRNAs that target broad-scale transcriptional disease phenotypes. It directly resolves differential regulation of pathways, irrespective of their differential gene expression, and captures co-activity to establish functional pathway groupings and the miRNAs that may regulate them. PanomiR uses a systems biology approach to provide broad but precise insights into miRNA-regulated functional programs. It is available at https://bioconductor.org/packages/PanomiR.
Assuntos
MicroRNAs , MicroRNAs/metabolismo , Biologia de Sistemas , Perfilação da Expressão Gênica/métodos , Biologia Computacional/métodos , Redes Reguladoras de GenesRESUMO
Attempts to develop drugs that address sepsis based on leads developed in animal models have failed. We sought to identify leads based on human data by exploiting a natural experiment: the relative resistance of children to mortality from severe infections and sepsis. Using public datasets, we identified key differences in pathway activity (Pathprint) in blood transcriptome profiles of septic adults and children. To find drugs that could promote beneficial (child) pathways or inhibit harmful (adult) ones, we built an in silico pathway drug network (PDN) using expression correlation between drug, disease, and pathway gene signatures across 58,475 microarrays. Specific pathway clusters from children or adults were assessed for correlation with drug-based signatures. Validation by literature curation and by direct testing in an endotoxemia model of murine sepsis of the most correlated drug candidates demonstrated that the Pathprint-PDN methodology is more effective at generating positive drug leads than gene-level methods (e.g., CMap). Pathway-centric Pathprint-PDN is a powerful new way to identify drug candidates for intervention against sepsis and provides direct insight into pathways that may determine survival.
Assuntos
Descoberta de Drogas , Avaliação Pré-Clínica de Medicamentos , Sepse/tratamento farmacológico , Sepse/mortalidade , Adolescente , Adulto , Idoso , Idoso de 80 Anos ou mais , Animais , Criança , Pré-Escolar , Análise por Conglomerados , Simulação por Computador , Modelos Animais de Doenças , Resistência à Doença , Feminino , Perfilação da Expressão Gênica , Humanos , Masculino , Camundongos , Camundongos Endogâmicos C57BL , Análise em Microsséries , Pessoa de Meia-Idade , Reprodutibilidade dos Testes , Sensibilidade e Especificidade , Transcriptoma , Adulto JovemRESUMO
Sonic hedgehog (Shh) signaling patterns the vertebrate spinal cord by activating a group of transcriptional repressors in distinct neural progenitors of somatic motor neuron and interneuron subtypes. To identify the action of this network, we performed a genome-wide analysis of the regulatory actions of three key ventral determinants in mammalian neural tube patterning: Nkx2.2, Nkx6.1 and Olig2. Previous studies have demonstrated that each factor acts predominantly as a transcriptional repressor, at least in part, to inhibit alternative progenitor fate choices. Here, we reveal broad and direct repression of multiple alternative fates as a general mechanism of repressor action. Additionally, the repressor network targets multiple Shh signaling components providing negative feedback to ongoing Shh signaling. Analysis of chromatin organization around Nkx2.2-, Nkx6.1- and Olig2-bound regions, together with co-analysis of engagement of the transcriptional activator Sox2, indicate that repressors bind to, and probably modulate the action of, neural enhancers. Together, the data suggest a model for neural progenitor specification downstream of Shh signaling, in which Nkx2.2 and Olig2 direct repression of alternative neural progenitor fate determinants, an action augmented by the overlapping activity of Nkx6.1 in each cell type. Integration of repressor and activator inputs, notably activator inputs mediated by Sox2, is probably a key mechanism in achieving cell type-specific transcriptional outcomes in mammalian neural progenitor fate specification.
Assuntos
Regulação da Expressão Gênica no Desenvolvimento/fisiologia , Proteínas Hedgehog/metabolismo , Modelos Biológicos , Células-Tronco Neurais/fisiologia , Tubo Neural/embriologia , Proteínas Repressoras/metabolismo , Transdução de Sinais/fisiologia , Fatores de Transcrição Hélice-Alça-Hélice Básicos/metabolismo , Cromatina/genética , Imunoprecipitação da Cromatina , Proteína Homeobox Nkx-2.2 , Proteínas de Homeodomínio/metabolismo , Humanos , Interneurônios/citologia , Neurônios Motores/citologia , Proteínas do Tecido Nervoso/metabolismo , Tubo Neural/metabolismo , Proteínas Nucleares , Fator de Transcrição 2 de Oligodendrócitos , Fases de Leitura Aberta/genética , Análise Serial de Proteínas , Ligação Proteica , Fatores de Transcrição SOXB1/metabolismo , Fatores de Transcrição/metabolismo , Proteínas de Peixe-ZebraRESUMO
Alzheimer's disease (AD) is a complex neurodegenerative disorder. The relative contribution of the numerous underlying functional mechanisms is poorly understood. To comprehensively understand the context and distribution of pathways that contribute to AD, we performed text-mining to generate an exhaustive, systematic assessment of the breadth and diversity of biological pathways within a corpus of 206,324 dementia publication abstracts. A total of 91% (325/335) of Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways have publications containing an association via at least 5 studies, while 63% of pathway terms have at least 50 studies providing a clear association with AD. Despite major technological advances, the same set of top-ranked pathways have been consistently related to AD for 30 years, including AD, immune system, metabolic pathways, cholinergic synapse, long-term depression, proteasome, diabetes, cancer, and chemokine signaling. AD pathways studied appear biased: animal model and human subject studies prioritize different AD pathways. Surprisingly, human genetic discoveries and drug targeting are not enriched in the most frequently studied pathways. Our findings suggest that not only is this disorder incredibly complex, but that its functional reach is also nearly global. As a consequence of our study, research results can now be assessed in the context of the wider AD literature, supporting the design of drug therapies that target a broader range of mechanisms. The results of this study can be explored at www.adpathways.org.
RESUMO
BACKGROUND: About 5% of western populations are afflicted by autoimmune diseases many of which are affected by sex hormones. Autoimmune diseases are complex and involve many genes. Identifying these disease-associated genes contributes to development of more effective therapies. Also, association studies frequently imply genomic regions that contain disease-associated genes but fall short of pinpointing these genes. The identification of disease-associated genes has always been challenging and to date there is no universal and effective method developed. RESULTS: We have developed a method to prioritize disease-associated genes for diseases affected strongly by sex hormones. Our method uses various types of information available for the genes, but no information that directly links genes with the disease. It generates a score for each of the considered genes and ranks genes based on that score. We illustrate our method on early-onset myasthenia gravis (MG) using genes potentially controlled by estrogen and localized in a genomic segment (which contains the MHC and surrounding region) strongly associated with MG. Based on the considered genomic segment 283 genes are ranked for their relevance to MG and responsiveness to estrogen. The top three ranked genes, HLA-G, TAP2 and HLA-DRB1, are implicated in autoimmune diseases, while TAP2 is associated with SNPs characteristic for MG. Within the top 35 prioritized genes our method identifies 90% of the 10 already known MG-associated genes from the considered region without using any information that directly links genes to MG. Among the top eight genes we identified HLA-G and TUBB as new candidates. We show that our ab-initio approach outperforms the other methods for prioritizing disease-associated genes. CONCLUSION: We have developed a method to prioritize disease-associated genes under the potential control of sex hormones. We demonstrate the success of this method by prioritizing the genes localized in the MHC and surrounding region and evaluating the role of these genes as potential candidates for estrogen control as well as MG. We show that our method outperforms the other methods. The method has a potential to be adapted to prioritize genes relevant to other diseases.
Assuntos
Hormônios Esteroides Gonadais/metabolismo , Miastenia Gravis/genética , Miastenia Gravis/metabolismo , Biologia de Sistemas/métodos , Estrogênios/metabolismo , Regulação da Expressão Gênica , Humanos , Timo/metabolismoRESUMO
Genome-wide techniques such as microarray analysis, Serial Analysis of Gene Expression (SAGE), Massively Parallel Signature Sequencing (MPSS), linkage analysis and association studies are used extensively in the search for genes that cause diseases, and often identify many hundreds of candidate disease genes. Selection of the most probable of these candidate disease genes for further empirical analysis is a significant challenge. Additionally, identifying the genes that cause complex diseases is problematic due to low penetrance of multiple contributing genes. Here, we describe a novel bioinformatic approach that selects candidate disease genes according to their expression profiles. We use the eVOC anatomical ontology to integrate text-mining of biomedical literature and data-mining of available human gene expression data. To demonstrate that our method is successful and widely applicable, we apply it to a database of 417 candidate genes containing 17 known disease genes. We successfully select the known disease gene for 15 out of 17 diseases and reduce the candidate gene set to 63.3% (+/-18.8%) of its original size. This approach facilitates direct association between genomic data describing gene expression and information from biomedical texts describing disease phenotype, and successfully prioritizes candidate genes according to their expression in disease-affected tissues.
Assuntos
Biologia Computacional/métodos , Predisposição Genética para Doença , Genômica/métodos , PubMed , Anatomia , Bases de Dados Genéticas , Perfilação da Expressão Gênica , Humanos , Integração de Sistemas , Vocabulário ControladoRESUMO
BACKGROUND: Rates of substitution in protein-coding sequences can provide important insights into evolutionary processes that are of biomedical and theoretical interest. Increased availability of coding sequence data has enabled researchers to estimate more accurately the coding sequence divergence of pairs of organisms. However the use of different data sources, alignment protocols and methods to estimate substitution rates leads to widely varying estimates of key parameters that define the coding sequence divergence of orthologous genes. Although complete genome sequence data are not available for all organisms, fragmentary sequence data can provide accurate estimates of substitution rates provided that an appropriate and consistent methodology is used and that differences in the estimates obtainable from different data sources are taken into account. RESULTS: We have developed FRAGS, an application framework that uses existing, freely available software components to construct in-frame alignments and estimate coding substitution rates from fragmentary sequence data. Coding sequence substitution estimates for human and chimpanzee sequences, generated by FRAGS, reveal that methodological differences can give rise to significantly different estimates of important substitution parameters. The estimated substitution rates were also used to infer upper-bounds on the amount of sequencing error in the datasets that we have analysed. CONCLUSION: We have developed a system that performs robust estimation of substitution rates for orthologous sequences from a pair of organisms. Our system can be used when fragmentary genomic or transcript data is available from one of the organisms and the other is a completely sequenced genome within the Ensembl database. As well as estimating substitution statistics our system enables the user to manage and query alignment and substitution data.
Assuntos
Mutagênese/genética , RNA Mensageiro/genética , Software , Animais , Biologia Computacional/métodos , Biologia Computacional/estatística & dados numéricos , Biologia Computacional/tendências , Bases de Dados Genéticas/estatística & dados numéricos , Humanos , Modelos Genéticos , Modelos Estatísticos , Pan troglodytes/genética , Software/estatística & dados numéricos , Software/tendênciasRESUMO
BACKGROUND: The intraerythrocytic development of Plasmodium falciparum, the most virulent human malaria parasite involves asexual and gametocyte stages. There has been a significant increase in disparate datasets derived from genomic and post-genomic analysis of the parasite that necessitates delivery of integrated analysis from which biological processes important to the survival of the parasite can be determined. METHODS: In order to resolve genes associated with stage differentially expressed transcripts, we have developed and implemented an integrative approach that combines evidence from P. falciparum expressed sequence tags (ESTs), genomic, microarray, proteomic and gene ontology data. RESULTS: A total of 143 gametocyte-overexpressed and 51 asexual-overexpressed transcripts were identified. A subset of 74 genes associated with these transcripts showed evidence of stage-correlated protein expression, of which 53 have not been experimentally characterised. Our study has revealed (1) possible regulatory mechanisms in malaria parasites' gametocyte maturation, (2) correlation between EST and microarray data for a P. falciparum gene family to present unique EST-derived information, (3) candidate drug and antigenic targets on which computational and experimental studies can be performed, and (4) the need for more empirical studies on gene and protein expression in malaria parasites. CONCLUSION: Applying different domains of data to the same underlying gene set has yielded novel insights into the biology of the parasite and presents an approach to appraise critically the data quality of post-genomic datasets from malaria parasites.
RESUMO
Parasites successfully exist within the host as a result of highly specific genetic adaptations. Therefore, detecting genes that contain relevant adaptive mutations can provide a guide to biological processes that are potentially essential to the parasite. Random genetic mutations that confer selective advantage can act to alter amino acids so as to confer gain of function that has positive impact on the survival of the parasite. Directional selection of advantageous mutations results from an elevated rate of nonsynonymous substitutions in rapidly evolving genes. Genes on which this positive selection operates are considered to have an evolutionary characteristic such that the normalized number of nonsynonymous (dn) substitutions is greater than that of synonymous (ds) substitutions. By searching in a statistically robust way for genes that contain this characteristic, it is possible to apply a stringent method to identify genes that may be under positive selection and thus to identify biological processes involving those genes that are essential to the survival of the parasite. Genes detected typically class into those under host immune surveillance and those intrinsic to pathways essential to survival of the parasite within the host. Depending on their function and location of protein expression, such genes have the potential to provide exceptional vaccine and drug candidates.
Assuntos
DNA/genética , Parasitos/genética , Sequência de Aminoácidos , Animais , Sequência de Bases , Dados de Sequência Molecular , Homologia de Sequência de AminoácidosRESUMO
Acute kidney injury (AKI) promotes an abrupt loss of kidney function that results in substantial morbidity and mortality. Considerable effort has gone toward identification of diagnostic biomarkers and analysis of AKI-associated molecular events; however, most studies have adopted organ-wide approaches and have not elucidated the interplay among different cell types involved in AKI pathophysiology. To better characterize AKI-associated molecular and cellular events, we developed a mouse line that enables the identification of translational profiles in specific cell types. This strategy relies on CRE recombinase-dependent activation of an EGFP-tagged L10a ribosomal protein subunit, which allows translating ribosome affinity purification (TRAP) of mRNA populations in CRE-expressing cells. Combining this mouse line with cell type-specific CRE-driver lines, we identified distinct cellular responses in an ischemia reperfusion injury (IRI) model of AKI. Twenty-four hours following IRI, distinct translational signatures were identified in the nephron, kidney interstitial cell populations, vascular endothelium, and macrophages/monocytes. Furthermore, TRAP captured known IRI-associated markers, validating this approach. Biological function annotation, canonical pathway analysis, and in situ analysis of identified response genes provided insight into cell-specific injury signatures. Our study provides a deep, cell-based view of early injury-associated molecular events in AKI and documents a versatile, genetic tool to monitor cell-specific and temporal-specific biological processes in disease modeling.
Assuntos
Injúria Renal Aguda/metabolismo , Traumatismo por Reperfusão/metabolismo , Transcriptoma , Injúria Renal Aguda/genética , Injúria Renal Aguda/patologia , Animais , Ontologia Genética , Proteínas de Fluorescência Verde/biossíntese , Proteínas de Fluorescência Verde/genética , Rim/irrigação sanguínea , Rim/metabolismo , Rim/patologia , Masculino , Camundongos , Camundongos da Linhagem 129 , Camundongos Endogâmicos C57BL , Camundongos Transgênicos , Especificidade de Órgãos , Biossíntese de Proteínas , Proteínas Recombinantes de Fusão/biossíntese , Proteínas Recombinantes de Fusão/genética , Traumatismo por Reperfusão/genética , Proteínas Ribossômicas/biossíntese , Proteínas Ribossômicas/genéticaRESUMO
New strategies to combat complex human disease require systems approaches to biology that integrate experiments from cell lines, primary tissues and model organisms. We have developed Pathprint, a functional approach that compares gene expression profiles in a set of pathways, networks and transcriptionally regulated targets. It can be applied universally to gene expression profiles across species. Integration of large-scale profiling methods and curation of the public repository overcomes platform, species and batch effects to yield a standard measure of functional distance between experiments. We show that pathprints combine mouse and human blood developmental lineage, and can be used to identify new prognostic indicators in acute myeloid leukemia. The code and resources are available at http://compbio.sph.harvard.edu/hidelab/pathprint.
RESUMO
The exponential increase in expressed sequence tag (EST) sequence data amplifies the computational cost of clustering sequences such that new algorithms are required to analyze data at a greater rate. We have parallelized d2_cluster on a SGI Origin 2000 multiprocessor and observed a speedup of approximately 100x on 126 processors when processing a 15,876 EST dataset. The parallelized d2_cluster code is obtainable from the SANBI website (http://www.sanbi.ac.za/CODES).