Arindam Bhattacharjee

Classification of human lung carcinomas by mRNA expression profiling reveals distinct adenocarcinoma subclasses

We have generated a molecular taxonomy of lung carcinoma, the leading cause of cancer death in the United States and worldwide. Using oligonucleotide microarrays, we analyzed mRNA expression levels corresponding to 12,600 transcript sequences in 186 lung tumor samples, including 139 adenocarcinomas resected from the lung. Hierarchical and probabilistic clustering of expression data defined distinct subclasses of lung adenocarcinoma. Among these were tumors with high relative expression of neuroendocrine genes and of type II pneumocyte genes, respectively. Retrospective analysis revealed a less favorable outcome for the adenocarcinomas with neuroendocrine gene expression. The diagnostic potential of expression profiling is emphasized by its ability to discriminate primary lung adenocarcinomas from metastases of extra-pulmonary origin. These results suggest that integration of expression profile data with clinical parameters could aid in diagnosis of lung cancer patients.

Targeted capture and massively parallel sequencing of 12 human exomes

Genome-wide association studies suggest that common genetic variants explain only a modest fraction of heritable risk for common diseases, raising the question of whether rare variants account for a significant fraction of unexplained heritability. Although DNA sequencing costs have fallen markedly, they remain far from what is necessary for rare and novel variants to be routinely identified at a genome-wide scale in large cohorts. We have therefore sought to develop second-generation methods for targeted sequencing of all protein-coding regions (‘exomes’), to reduce costs while enriching for discovery of highly penetrant variants. Here we report on the targeted capture and massively parallel sequencing of the exomes of 12 humans. These include eight HapMap individuals representing three populations, and four unrelated individuals with a rare dominantly inherited disorder, Freeman–Sheldon syndrome (FSS). We demonstrate the sensitive and specific identification of rare and common variants in over 300 megabases of coding sequence. Using FSS as a proof-of-concept, we show that candidate genes for Mendelian disorders can be identified by exome sequencing of a small number of unrelated, affected individuals. This strategy may be extendable to diseases with more complex genetics through larger sample sizes and appropriate weighting of non-synonymous variants by predicted functional impact.

Short interfering RNAs can induce unexpected and divergent changes in the levels of untargeted proteins in mammalian cells

RNA interference (RNAi) mediated by short interfering RNAs (siRNAs) is a widely used method to analyze gene function. To use RNAi knockdown accurately to infer gene function, it is essential to determine the specificity of siRNA-mediated RNAi. We have assessed the specificity of 10 different siRNAs corresponding to the MEN1 gene by examining the expression of two additional genes, TP53 (p53) and CDKN1A (p21), which are considered functionally unrelated to menin but are sensitive markers of cell state. MEN1 RNA and corresponding protein levels were all reduced after siRNA transfection of HeLa cells, although the degree of inhibition mediated by individual siRNAs varied. Unexpectedly, we observed dramatic and significant changes in protein levels of p53 and p21 that were unrelated to silencing of the target gene. The modulations in p53 and p21 levels were not abolished on titration of the siRNAs, and similar results were obtained in three other cell lines; in none of the cell lines tested did we see an effect on the protein levels of actin. These data suggest that siRNAs can induce nonspecific effects on protein levels that are siRNA sequence dependent but that these effects may be difficult to detect until genes central to a pivotal cellular response, such as p53 and p21, are studied. We find no evidence that activation of the double-stranded RNA-triggered IFN-associated antiviral pathways accounts for these effects, but we speculate that partial complementary sequence matches to off-target genes may result in a micro-RNA-like inhibition of translation.

Menin associates with a trithorax family histone methyltransferase complex and with the Hoxc8 locus

The cellular function of the menin tumor suppressor protein, product of the MEN1 gene mutated in familial multiple endocrine neoplasia type 1, has not been defined. We now show that menin is associated with a histone methyltransferase complex containing two trithorax family proteins, MLL2 and Ash2L, and other homologs of the yeast Set1 assembly. This menin-associated complex methylates histone H3 on lysine 4. A subset of tumor-derived menin mutants lacks the associated histone methyltransferase activity. In addition, menin is associated with RNA polymerase II whose large subunit carboxyl-terminal domain is phosphorylated on Ser 5. Men1 knockout embryos and cells show decreased expression of the homeobox genes Hoxc6 and Hoxc8. Chromatin immunoprecipitation experiments reveal that menin is bound to the Hoxc8 locus. These results suggest that menin activates the transcription of differentiation-regulating genes by covalent histone modification, and that this activity is related to tumor suppression by MEN1.

Features of Mammalian microRNA Promoters Emerge from Polymerase II Chromatin Immunoprecipitation Data

Background MicroRNAs (miRNAs) are short, non-coding RNA regulators of protein coding genes. miRNAs play a very important role in diverse biological processes and various diseases. Many algorithms are able to predict miRNA genes and their targets, but their transcription regulation is still under investigation. It is generally believed that intragenic miRNAs (located in introns or exons of protein coding genes) are co-transcribed with their host genes and most intergenic miRNAs transcribed from their own RNA polymerase II (Pol II) promoter. However, the length of the primary transcripts and promoter organization is currently unknown. Methodology We performed Pol II chromatin immunoprecipitation (ChIP)-chip using a custom array surrounding regions of known miRNA genes. To identify the true core transcription start sites of the miRNA genes we developed a new tool (CPPP). We showed that miRNA genes can be transcribed from promoters located several kilobases away and that their promoters share the same general features as those of protein coding genes. Finally, we found evidence that as many as 26% of the intragenic miRNAs may be transcribed from their own unique promoters. Conclusion miRNA promoters have similar features to those of protein coding genes, but miRNA transcript organization is more complex.

Targeted Investigation of the Neandertal Genome by Array-Based Sequence Capture.

Kissing Cousins Neandertals, our closest relatives, ranged across Europe and Southwest Asia before their extinction approximately 30,000 years ago. Green et al. (p. 710) report a draft sequence of the Neandertal genome, created from three individuals, and compare it with genomes of five modern humans. The results suggest that ancient genomes of human relatives can be recovered with acceptably low contamination from modern human DNA. Because ancient DNA can be contaminated with microbial DNA, Burbano et al. (p. 723) developed a target sequence capture approach to obtain 14 kilobases of Neandertal DNA from a fairly poorly preserved sample with a high microbial load. A number of genomic regions and genes were revealed as candidates for positive selection early in modern human history. The genomic data suggest that Neandertals mixed with modern human ancestors some 120,000 years ago, leaving traces of Neandertal DNA in contemporary humans. Array capture of Neandertal DNA identifies amino acid substitutions that occurred after the split between humans and Neandertals. It is now possible to perform whole-genome shotgun sequencing as well as capture of specific genomic regions for extinct organisms. However, targeted resequencing of large parts of nuclear genomes has yet to be demonstrated for ancient DNA. Here we show that hybridization capture on microarrays can successfully recover more than a megabase of target regions from Neandertal DNA even in the presence of ~99.8% microbial DNA. Using this approach, we have sequenced ~14,000 protein-coding positions inferred to have changed on the human lineage since the last common ancestor shared with chimpanzees. By generating the sequence of one Neandertal and 50 present-day humans at these positions, we have identified 88 amino acid substitutions that have become fixed in humans since our divergence from the Neandertals.

Histone H4K20/H3K9 demethylase PHF8 regulates zebrafish brain and craniofacial development

X-linked mental retardation (XLMR) is a complex human disease that causes intellectual disability. Causal mutations have been found in approximately 90 X-linked genes; however, molecular and biological functions of many of these genetically defined XLMR genes remain unknown. PHF8 (PHD (plant homeo domain) finger protein 8) is a JmjC domain-containing protein and its mutations have been found in patients with XLMR and craniofacial deformities. Here we provide multiple lines of evidence establishing PHF8 as the first mono-methyl histone H4 lysine 20 (H4K20me1) demethylase, with additional activities towards histone H3K9me1 and me2. PHF8 is located around the transcription start sites (TSS) of ∼7,000 RefSeq genes and in gene bodies and intergenic regions (non-TSS). PHF8 depletion resulted in upregulation of H4K20me1 and H3K9me1 at the TSS and H3K9me2 in the non-TSS sites, respectively, demonstrating differential substrate specificities at different target locations. PHF8 positively regulates gene expression, which is dependent on its H3K4me3-binding PHD and catalytic domains. Importantly, patient mutations significantly compromised PHF8 catalytic function. PHF8 regulates cell survival in the zebrafish brain and jaw development, thus providing a potentially relevant biological context for understanding the clinical symptoms associated with PHF8 patients. Lastly, genetic and molecular evidence supports a model whereby PHF8 regulates zebrafish neuronal cell survival and jaw development in part by directly regulating the expression of the homeodomain transcription factor MSX1/MSXB, which functions downstream of multiple signalling and developmental pathways. Our findings indicate that an imbalance of histone methylation dynamics has a critical role in XLMR.

High definition profiling of mammalian DNA methylation by array capture and single molecule bisulfite sequencing

DNA methylation stabilizes developmentally programmed gene expression states. Aberrant methylation is associated with disease progression and is a common feature of cancer genomes. Presently, few methods enable quantitative, large-scale, single-base resolution mapping of DNA methylation states in desired regions of a complex mammalian genome. Here, we present an approach that combines array-based hybrid selection and massively parallel bisulfite sequencing to profile DNA methylation in genomic regions spanning hundreds of thousands of bases. This single molecule strategy enables methylation variable positions to be quantitatively examined with high sampling precision. Using bisulfite capture, we assessed methylation patterns across 324 randomly selected CpG islands (CGI) representing more than 25,000 CpG sites. A single lane of Illumina sequencing permitted methylation states to be definitively called for >90% of target sties. The accuracy of the hybrid-selection approach was verified using conventional bisulfite capillary sequencing of cloned PCR products amplified from a subset of the selected regions. This confirmed that even partially methylated states could be successfully called. A comparison of human primary and cancer cells revealed multiple differentially methylated regions. More than 25% of islands showed complex methylation patterns either with partial methylation states defining the entire CGI or with contrasting methylation states appearing in specific regional blocks within the island. We observed that transitions in methylation state often correlate with genomic landmarks, including transcriptional start sites and intron-exon junctions. Methylation, along with specific histone marks, was enriched in exonic regions, suggesting that chromatin states can foreshadow the content of mature mRNAs.

Gene Expression Profiling Reveals Reproducible Human Lung Adenocarcinoma Subtypes in Multiple Independent Patient Cohorts

PURPOSE Published reports suggest that DNA microarrays identify clinically meaningful subtypes of lung adenocarcinomas not recognizable by other routine tests. This report is an investigation of the reproducibility of the reported tumor subtypes. METHODS Three independent cohorts of patients with lung cancer were evaluated using a variety of DNA microarray assays. Using the integrative correlations method, a subset of genes was selected, the reliability of which was acceptable across the different DNA microarray platforms. Tumor subtypes were selected using consensus clustering and genes distinguishing subtypes were identified using the weighted difference statistic. Gene lists were compared across cohorts using centroids and gene set enrichment analysis. RESULTS Cohorts of 31, 72, and 128 adenocarcinomas were generated for a total of 231 microarrays, each with 2,553 reliable genes. Three adenocarcinoma subtypes were identified in each cohort. These were named bronchioid, squamoid, and magnoid according to their respective correlations with gene expression patterns from histologically defined bronchioalveolar carcinoma, squamous cell carcinoma, and large-cell carcinoma. Tumor subtypes were distinguishable by many hundreds of genes, and lists generated in one cohort were predictive of tumor subtypes in the two other cohorts. Tumor subtypes correlated with clinically relevant covariates, including stage-specific survival and metastatic pattern. Most notably, bronchioid tumors were correlated with improved survival in early-stage disease, whereas squamoid tumors were associated with better survival in advanced disease. CONCLUSION DNA microarray analysis of lung adenocarcinomas identified reproducible tumor subtypes which differ significantly in clinically important behaviors such as stage-specific survival.

Hybrid selection of discrete genomic intervals on custom-designed microarrays for massively parallel sequencing

Complementary techniques that deepen information content and minimize reagent costs are required to realize the full potential of massively parallel sequencing. Here, we describe a resequencing approach that directs focus to genomic regions of high interest by combining hybridization-based purification of multi-megabase regions with sequencing on the Illumina Genome Analyzer (GA). The capture matrix is created by a microarray on which probes can be programmed as desired to target any non-repeat portion of the genome, while the method requires only a basic familiarity with microarray hybridization. We present a detailed protocol suitable for 1–2 μg of input genomic DNA and highlight key design tips in which high specificity (>65% of reads stem from enriched exons) and high sensitivity (98% targeted base pair coverage) can be achieved. We have successfully applied this to the enrichment of coding regions, in both human and mouse, ranging from 0.5 to 4 Mb in length. From genomic DNA library production to base-called sequences, this procedure takes approximately 9–10 d inclusive of array captures and one Illumina flow cell run.