Huai Li

Genome wide profiling of human embryonic stem cells (hESCs), their derivatives and embryonal carcinoma cells to develop base profiles of U.S. Federal government approved hESC lines

BackgroundIn order to compare the gene expression profiles of human embryonic stem cell (hESC) lines and their differentiated progeny and to monitor feeder contaminations, we have examined gene expression in seven hESC lines and human fibroblast feeder cells using Illumina® bead arrays that contain probes for 24,131 transcript probes.ResultsA total of 48 different samples (including duplicates) grown in multiple laboratories under different conditions were analyzed and pairwise comparisons were performed in all groups. Hierarchical clustering showed that blinded duplicates were correctly identified as the closest related samples. hESC lines clustered together irrespective of the laboratory in which they were maintained. hESCs could be readily distinguished from embryoid bodies (EB) differentiated from them and the karyotypically abnormal hESC line BG01V. The embryonal carcinoma (EC) line NTera2 is a useful model for evaluating characteristics of hESCs. Expression of subsets of individual genes was validated by comparing with published databases, MPSS (Massively Parallel Signature Sequencing) libraries, and parallel analysis by microarray and RT-PCR.Conclusionwe show that Illuminas bead array platform is a reliable, reproducible and robust method for developing base global profiles of cells and identifying similarities and differences in large number of samples.

Gene expression atlas of the mouse central nervous system: impact and interactions of age, energy intake and gender.

BackgroundThe structural and functional complexity of the mammalian central nervous system (CNS) is organized and modified by complicated molecular signaling processes that are poorly understood.ResultsWe measured transcripts of 16,896 genes in 5 CNS regions from cohorts of young, middle-aged and old male and female mice that had been maintained on either a control diet or a low energy diet known to retard aging. Each CNS region (cerebral cortex, hippocampus, striatum, cerebellum and spinal cord) possessed its own unique transcriptome fingerprint that was independent of age, gender and energy intake. Less than 10% of genes were significantly affected by age, diet or gender, with most of these changes occurring between middle and old age. The transcriptome of the spinal cord was the most responsive to age, diet and gender, while the striatal transcriptome was the least responsive. Gender and energy restriction had particularly robust influences on the hippocampal transcriptome of middle-aged mice. Prominent functional groups of age- and energy-sensitive genes were those encoding proteins involved in DNA damage responses (Werner and telomere-associated proteins), mitochondrial and proteasome functions, cell fate determination (Wnt and Notch signaling) and synaptic vesicle trafficking.ConclusionMouse CNS transcriptomes responded to age, energy intake and gender in a regionally distinctive manner. The systematic transcriptome dataset also provides a window into mechanisms of age-, diet- and sex-related CNS plasticity and vulnerability.

Evolutionarily Conserved Transcriptional Co-Expression Guiding Embryonic Stem Cell Differentiation

Background Understanding the molecular mechanisms controlling pluripotency in embryonic stem cells (ESCs) is of central importance towards realizing their potentials in medicine and science. Cross-species examination of transcriptional co-expression allows elucidation of fundamental and species-specific mechanisms regulating ESC self-renewal or differentiation. Methodology/Principal Findings We examined transcriptional co-expression of ESCs from pathways to global networks under the framework of human-mouse comparisons. Using generalized singular value decomposition and comparative partition around medoids algorithms, evolutionarily conserved and divergent transcriptional co-expression regulating pluripotency were identified from ESC-critical pathways including ACTIVIN/NODAL, ATK/PTEN, BMP, CELL CYCLE, JAK/STAT, PI3K, TGFβ and WNT. A set of transcription factors, including FOX, GATA, MYB, NANOG, OCT, PAX, SOX and STAT, and the FGF response element were identified that represent key regulators underlying the transcriptional co-expression. By transcriptional intervention conducted in silico, dynamic behavior of pathways was examined, which demonstrate how much and in which specific ways each gene or gene combination effects the behavior transition of a pathway in response to ESC differentiation or pluripotency induction. The global co-expression networks of ESCs were dominated by highly connected hub genes such as IGF2, JARID2, LCK, MYCN, NASP, OCT4, ORC1L, PHC1 and RUVBL1, which are possibly critical in determining the fate of ESCs. Conclusions/Significance Through these studies, evolutionary conservation at genomic, transcriptomic, and network levels is shown to be an effective predictor of molecular factors and mechanisms controlling ESC development. Various hypotheses regarding mechanisms controlling ESC development were generated, which could be further validated by in vitro experiments. Our findings shed light on the systems-level understanding of how ESC differentiation or pluripotency arises from the connectivity or networks of genes, and provide a “road-map” for further experimental investigation.

Unraveling transcriptional regulatory programs by integrative analysis of microarray and transcription factor binding data

Motivation: Unraveling the transcriptional regulatory program mediated by transcription factors (TFs) is a fundamental objective of computational biology, yet still remains a challenge. Method: Here, we present a new methodology that integrates microarray and TF binding data for unraveling transcriptional regulatory networks. The algorithm is based on a two-stage constrained matrix decomposition model. The model takes into account the non-linear structure in gene expression data, particularly in the TF-target gene interactions and the combinatorial nature of gene regulation by TFs. The gene expression profile is modeled as a linear weighted combination of the activity profiles of a set of TFs. The TF activity profiles are deduced from the expression levels of TF target genes, instead directly from TFs themselves. The TF-target gene relationships are derived from ChIP-chip and other TF binding data. The proposed algorithm can not only identify transcriptional modules, but also reveal regulatory programs of which TFs control which target genes in which specific ways (either activating or inhibiting). Results: In comparison with other methods, our algorithm identifies biologically more meaningful transcriptional modules relating to specific TFs. We applied the new algorithm on yeast cell cycle and stress response data. While known transcriptional regulations were confirmed, novel TF-gene interactions were predicted and provide new insights into the regulatory mechanisms of the cell. Contact: zhanmi@mail.nih.gov Supplementary information: Supplementary data are available at Bioinformatics online.

Transcriptome coexpression map of human embryonic stem cells

BackgroundHuman embryonic stem (ES) cells hold great promise for medicine and science. The transcriptome of human ES cells has been studied in detail in recent years. However, no systematic analysis has yet addressed whether gene expression in human ES cells may be regulated in chromosomal domains, and no chromosomal domains of coexpression have been identified.ResultsWe report the first transcriptome coexpression map of the human ES cell and the earliest stage of ES differentiation, the embryoid body (EB), for the analysis of how transcriptional regulation interacts with genomic structure during ES self-renewal and differentiation. We determined the gene expression profiles from multiple ES and EB samples and identified chromosomal domains showing coexpression of adjacent genes on the genome. The coexpression domains were not random, with significant enrichment in chromosomes 8, 11, 16, 17, 19, and Y in the ES state, and 6, 11, 17, 19 and 20 in the EB state. The domains were significantly associated with Giemsa-negative bands in EB, yet showed little correlation with known cytogenetic structures in ES cells. Different patterns of coexpression were revealed by comparative transcriptome mapping between ES and EB.ConclusionThe findings and methods reported in this investigation advance our understanding of how genome organization affects gene expression in human ES cells and help to identify new mechanisms and pathways controlling ES self-renewal or differentiation.

Unraveling Regulatory Programs for NF-kappaB, p53 and MicroRNAs in Head and Neck Squamous Cell Carcinoma

In head and neck squamous cell carcinoma (HNSCC), mutations of p53 usually coexist with aberrant activation of NF-kappaB (NF-κB), other transcription factors and microRNAs, which promote tumor pathogenesis. However, how these factors and microRNAs interact to globally modulate gene expression and mediate oncogenesis is not fully understood. We devised a novel bioinformatics method to uncover interactive relationships between transcription factors or microRNAs and genes. This approach is based on matrix decomposition modeling under the joint constraints of sparseness and regulator-target connectivity, and able to integrate gene expression profiling and binding data of regulators. We employed this method to infer the gene regulatory networks in HNSCC. We found that the majority of the predicted p53 targets overlapped with those for NF-κB, suggesting that the two transcription factors exert a concerted modulation on regulatory programs in tumor cells. We further investigated the interrelationships of p53 and NF-κB with five additional transcription factors, AP1, CEBPB, EGR1, SP1 and STAT3, and microRNAs mir21 and mir34ac. The resulting gene networks indicate that interactions among NF-κB, p53, and the two miRNAs likely regulate progression of HNSCC. We experimentally validated our findings by determining expression of the predicted NF-κB and p53 target genes by siRNA knock down, and by examining p53 binding activity on promoters of predicted target genes in the tumor cell lines. Our results elucidating the cross-regulations among NF-κB, p53, and microRNAs provide insights into the complex regulatory mechanisms underlying HNSCC, and shows an efficient approach to inferring gene regulatory programs in biological complex systems.

Analysis of gene coexpression by B-spline based CoD estimation

The gene coexpression study has emerged as a novel holistic approach for microarray data analysis. Different indices have been used in exploring coexpression relationship, but each is associated with certain pitfalls. The Pearsons correlation coefficient, for example, is not capable of uncovering nonlinear pattern and directionality of coexpression. Mutual information can detect nonlinearity but fails to show directionality. The coefficient of determination (CoD) is unique in exploring different patterns of gene coexpression, but so far only applied to discrete data and the conversion of continuous microarray data to the discrete format could lead to information loss. Here, we proposed an effective algorithm, CoexPro, for gene coexpression analysis. The new algorithm is based on B-spline approximation of coexpression between a pair of genes, followed by CoD estimation. The algorithm was justified by simulation studies and by functional semantic similarity analysis. The proposed algorithm is capable of uncovering both linear and a specific class of nonlinear relationships from continuous microarray data. It can also provide suggestions for possible directionality of coexpression to the researchers. The new algorithm presents a novel model for gene coexpression and will be a valuable tool for a variety of gene expression and network studies. The application of the algorithm was demonstrated by an analysis on ligand-receptor coexpression in cancerous and noncancerous cells. The software implementing the algorithm is available upon request to the authors.

Exploring pathways from gene co-expression to network dynamics.

One of the major challenges in post-genomic research is to understand how physiological and pathological phenotypes arise from the networks or connectivity of expressed genes. In addressing this issue, we have developed two computational algorithms, CoExMiner and PathwayPro, to explore static features of gene co-expression and dynamic behaviors of gene networks. CoExMiner is based on B-spline approximation followed by the coefficient of determination (CoD) estimation for modeling gene co-expression patterns. The algorithm allows the exploration of transcriptional responses that involve coordinated expression of genes encoding proteins which work in concert in the cell. PathwayPro is based on a finite-state Markov chain model for mimicking dynamic behaviors of a transcriptional network. The algorithm allows quantitative assessment of a wide range of network responses, including susceptibility to disease, potential usefulness of a given drug, and consequences of such external stimuli as pharmacological interventions or caloric restriction. We demonstrated the applications of CoExMiner and PathwayPro by examining gene expression profiles of ligands and receptors in cancerous and non-cancerous cells and network dynamics of the leukemia-associated BCR-ABL pathway. The examinations disclosed both linear and nonlinear relationships of ligand-receptor interactions associated with cancer development, identified disease and drug targets of leukemia, and provided new insights into biology of the diseases. The analysis using these newly developed algorithms show the great usefulness of computational systems biology approaches for biological and medical research.