Population genomics of human gene expression

Genetic variation influences gene expression, and this variation in gene expression can be efficiently mapped to specific genomic regions and variants. Here we have used gene expression profiling of Epstein-Barr virus-transformed lymphoblastoid cell lines of all 270 individuals genotyped in the HapMap Consortium to elucidate the detailed features of genetic variation underlying gene expression variation. We find that gene expression is heritable and that differentiation between populations is in agreement with earlier small-scale studies. A detailed association analysis of over 2.2 million common SNPs per population (5% frequency in HapMap) with gene expression identified at least 1,348 genes with association signals in cis and at least 180 in trans. Replication in at least one independent population was achieved for 37% of cis signals and 15% of trans signals, respectively. Our results strongly support an abundance of cis-regulatory variation in the human genome. Detection of trans effects is limited but suggests that regulatory variation may be the key primary effect contributing to phenotypic variation in humans. We also explore several methodologies that improve the current state of analysis of gene expression variation.     

Adenovirus-mediated in vivo gene transfer and expression in normal rat liver.

Replication deficient, recombinant adenovirus (Ad) vectors do not require target cell replication for transfer and expression of exogenous genes and thus may be useful for in vivo gene therapy in hepatocytes. In vitro, primary cultures of rat hepatocytes infected with a recombinant Ad containing a human alpha 1-antitrypsin cDNA (Ad-alpha 1AT) synthesized and secreted human alpha 1AT for 4 weeks. In rats, in vivo intraportal administration of a recombinant Ad containing the E. coli lacZ gene, was followed by expression of beta-galactosidase in hepatocytes 3 days after infection. Intraportal infusion of Ad-alpha 1AT produced detectable serum levels of human alpha 1AT for 4 weeks. Thus, targeted gene expression has been achieved in the liver, albeit at low levels, suggesting that adenovirus vectors may be a useful means for in vivo gene therapy in liver disorders.     

Systematic variation in gene expression patterns in human cancer cell lines

We used cDNA microarrays to explore the variation in expression of approximately 8,000 unique genes among the 60 cell lines used in the National Cancer Institute's screen for anti-cancer drugs. Classification of the cell lines based solely on the observed patterns of gene expression revealed a correspondence to the ostensible origins of the tumours from which the cell lines were derived. The consistent relationship between the gene expression patterns and the tissue of origin allowed us to recognize outliers whose previous classification appeared incorrect. Specific features of the gene expression patterns appeared to be related to physiological properties of the cell lines, such as their doubling time in culture, drug metabolism or the interferon response. Comparison of gene expression patterns in the cell lines to those observed in normal breast tissue or in breast tumour specimens revealed features of the expression patterns in the tumours that had recognizable counterparts in specific cell lines, reflecting the tumour, stromal and inflammatory components of the tumour tissue. These results provided a novel molecular characterization of this important group of human cell lines and their relationships to tumours in vivo.     

MLL translocations specify a distinct gene expression profile that distinguishes a unique leukemia.

Acute lymphoblastic leukemias carrying a chromosomal translocation involving the mixed-lineage leukemia gene (MLL, ALL1, HRX) have a particularly poor prognosis. Here we show that they have a characteristic, highly distinct gene expression profile that is consistent with an early hematopoietic progenitor expressing select multilineage markers and individual HOX genes. Clustering algorithms reveal that lymphoblastic leukemias with MLL translocations can clearly be separated from conventional acute lymphoblastic and acute myelogenous leukemias. We propose that they constitute a distinct disease, denoted here as MLL, and show that the differences in gene expression are robust enough to classify leukemias correctly as MLL, acute lymphoblastic leukemia or acute myelogenous leukemia. Establishing that MLL is a unique entity is critical, as it mandates the examination of selectively expressed genes for urgently needed molecular targets.     

A literature network of human genes for high-throughput analysis of gene expression

We have carried out automated extraction of explicit and implicit biomedical knowledge from publicly available gene and text databases to create a gene-to-gene co-citation network for 13,712 named human genes by automated analysis of titles and abstracts in over 10 million MEDLINE records. The associations between genes have been annotated by linking genes to terms from the medical subject heading (MeSH) index and terms from the gene ontology (GO) database. The extracted database and accompanying web tools for gene-expression analysis have collectively been named 'PubGene'. We validated the extracted networks by three large-scale experiments showing that co-occurrence reflects biologically meaningful relationships, thus providing an approach to extract and structure known biology. We validated the applicability of the tools by analyzing two publicly available microarray data sets.     

A signature of chromosomal instability inferred from gene expression profiles predicts clinical outcome in multiple human cancers

We developed a computational method to characterize aneuploidy in tumor samples based on coordinated aberrations in expression of genes localized to each chromosomal region. We summarized the total level of chromosomal aberration in a given tumor in a univariate measure termed total functional aneuploidy. We identified a signature of chromosomal instability from specific genes whose expression was consistently correlated with total functional aneuploidy in several cancer types. Net overexpression of this signature was predictive of poor clinical outcome in 12 cancer data sets representing six cancer types. Also, the signature of chromosomal instability was higher in metastasis samples than in primary tumors and was able to stratify grade 1 and grade 2 breast tumors according to clinical outcome. These results provide a means to assess the potential role of chromosomal instability in determining malignant potential over a broad range of tumors.     

Global assessment of promoter methylation in a mouse model of cancer identifies ID4 as a putative tumor-suppressor gene in human leukemia

DNA methylation is associated with malignant transformation, but limitations imposed by genetic variability, tumor heterogeneity, availability of paired normal tissues and methodologies for global assessment of DNA methylation have limited progress in understanding the extent of epigenetic events in the initiation and progression of human cancer and in identifying genes that undergo methylation during cancer. We developed a mouse model of T/natural killer acute lymphoblastic leukemia that is always preceded by polyclonal lymphocyte expansion to determine how aberrant promoter DNA methylation and consequent gene silencing might be contributing to leukemic transformation. We used restriction landmark genomic scanning with this mouse model of preleukemia reproducibly progressing to leukemia to show that specific genomic methylation is associated with only the leukemic phase and is not random. We also identified Idb4 as a putative tumor-suppressor gene that is methylated in most mouse and human leukemias but in only a minority of other human cancers.     

Small changes in expression affect predisposition to tumorigenesis.

We have used quantitative measures of gene expression to show that constitutional 50% decreases in expression of one adenomatous polyposis coli tumor suppressor gene (APC) allele can lead to the development of familial adenomatous polyposis.     

Bmi1 is expressed in vivo in intestinal stem cells

Bmi1 plays an essential part in the self-renewal of hematopoietic and neural stem cells. To investigate its role in other adult stem cell populations, we generated a mouse expressing a tamoxifen-inducible Cre from the Bmi1 locus. We found that Bmi1 is expressed in discrete cells located near the bottom of crypts in the small intestine, predominantly four cells above the base of the crypt (+4 position). Over time, these cells proliferate, expand, self-renew and give rise to all the differentiated cell lineages of the small intestine epithelium. The induction of a stable form of beta-catenin in these cells was sufficient to rapidly generate adenomas. Moreover, ablation of Bmi1(+) cells using a Rosa26 conditional allele, expressing diphtheria toxin, led to crypt loss. These experiments identify Bmi1 as an intestinal stem cell marker in vivo. Unexpectedly, the distribution of Bmi1-expressing stem cells along the length of the small intestine suggested that mammals use more than one molecularly distinguishable adult stem cell subpopulation to maintain organ homeostasis.     

Genetics of gene expression in primary immune cells identifies cell type-specific master regulators and roles of HLA alleles

Trans-acting genetic variants have a substantial, albeit poorly characterized, role in the heritable determination of gene expression. Using paired purified primary monocytes and B cells, we identify new predominantly cell type-specific cis and trans expression quantitative trait loci (eQTLs), including multi-locus trans associations to LYZ and KLF4 in monocytes and B cells, respectively. Additionally, we observe a B cell-specific trans association of rs11171739 at 12q13.2, a known autoimmune disease locus, with IP6K2 (P = 5.8 × 10(-15)), PRIC285 (P = 3.0 × 10(-10)) and an upstream region of CDKN1A (P = 2 × 10(-52)), suggesting roles for cell cycle regulation and peroxisome proliferator-activated receptor γ (PPARγ) signaling in autoimmune pathogenesis. We also find that specific human leukocyte antigen (HLA) alleles form trans associations with the expression of AOAH and ARHGAP24 in monocytes but not in B cells. In summary, we show that mapping gene expression in defined primary cell populations identifies new cell type-specific trans-regulated networks and provides insights into the genetic basis of disease susceptibility.     

MicroRNAs can generate thresholds in target gene expression

MicroRNAs (miRNAs) are short, highly conserved noncoding RNA molecules that repress gene expression in a sequence-dependent manner. We performed single-cell measurements using quantitative fluorescence microscopy and flow cytometry to monitor a target gene's protein expression in the presence and absence of regulation by miRNA. We find that although the average level of repression is modest, in agreement with previous population-based measurements, the repression among individual cells varies dramatically. In particular, we show that regulation by miRNAs establishes a threshold level of target mRNA below which protein production is highly repressed. Near this threshold, protein expression responds sensitively to target mRNA input, consistent with a mathematical model of molecular titration. These results show that miRNAs can act both as a switch and as a fine-tuner of gene expression.