Frequent mutation of histone-modifying genes in non-Hodgkin lymphoma.

Follicular lymphoma (FL) and diffuse large B-cell lymphoma (DLBCL) are the two most common non-Hodgkin lymphomas (NHLs). Here we sequenced tumour and matched normal DNA from 13 DLBCL cases and one FL case to identify genes with mutations in B-cell NHL. We analysed RNA-seq data from these and another 113 NHLs to identify genes with candidate mutations, and then re-sequenced tumour and matched normal DNA from these cases to confirm 109 genes with multiple somatic mutations. Genes with roles in histone modification were frequent targets of somatic mutation. For example, 32% of DLBCL and 89% of FL cases had somatic mutations in MLL2, which encodes a histone methyltransferase, and 11.4% and 13.4% of DLBCL and FL cases, respectively, had mutations in MEF2B, a calcium-regulated gene that cooperates with CREBBP and EP300 in acetylating histones. Our analysis suggests a previously unappreciated disruption of chromatin biology in lymphomagenesis.

Conserved role of intragenic DNA methylation in regulating alternative promoters.

Although it is known that the methylation of DNA in 5' promoters suppresses gene expression, the role of DNA methylation in gene bodies is unclear. In mammals, tissue- and cell type-specific methylation is present in a small percentage of 5' CpG island (CGI) promoters, whereas a far greater proportion occurs across gene bodies, coinciding with highly conserved sequences. Tissue-specific intragenic methylation might reduce, or, paradoxically, enhance transcription elongation efficiency. Capped analysis of gene expression (CAGE) experiments also indicate that transcription commonly initiates within and between genes. To investigate the role of intragenic methylation, we generated a map of DNA methylation from the human brain encompassing 24.7 million of the 28 million CpG sites. From the dense, high-resolution coverage of CpG islands, the majority of methylated CpG islands were shown to be in intragenic and intergenic regions, whereas less than 3% of CpG islands in 5' promoters were methylated. The CpG islands in all three locations overlapped with RNA markers of transcription initiation, and unmethylated CpG islands also overlapped significantly with trimethylation of H3K4, a histone modification enriched at promoters. The general and CpG-island-specific patterns of methylation are conserved in mouse tissues. An in-depth investigation of the human SHANK3 locus and its mouse homologue demonstrated that this tissue-specific DNA methylation regulates intragenic promoter activity in vitro and in vivo. These methylation-regulated, alternative transcripts are expressed in a tissue- and cell type-specific manner, and are expressed differentially within a single cell type from distinct brain regions. These results support a major role for intragenic methylation in regulating cell context-specific alternative promoters in gene bodies.

BMPER variants associated with a novel, attenuated subtype of diaphanospondylodysostosis.

Diaphanospondylodysostosis (DSD), caused by loss of bone morphogenetic protein-binding endothelial regulator (BMPER), has been considered a lethal skeletal dysplasia characterized by severe deficiency of vertebral body and sacral ossification, reduced rib number and cystic kidneys. In this study, however, we have demonstrated that variants in BMPER may cause a milder disorder, without renal anomalies, that is compatible with long-term survival. Four siblings, three males and one female, presented with severe congenital scoliosis associated with rib and vertebral malformations as well as strikingly delayed ossification of the pedicles. The female was stillborn from an unrelated cause. Stabilization of the scoliosis with expandable titanium rods was successful in the three boys, all of whom have short stature. An autosomal recessive mode of inheritance was hypothesized. Single nucleotide polymorphism microarray analysis was performed for three of the siblings to identify autosomal genes with shared allele patterns, suggesting possible linkage. Exome sequencing of one sibling was then performed. Rare variants were identified in 347 genes with shared alleles. Only one of these genes had bi-allelic variants in a gene strongly expressed in paraxial mesenchyme: BMPER, which is the cause of DSD, an autosomal recessive disorder. The disorder described herein could represent an attenuated form of DSD or could be designated a separate entity such as spondylopedicular dysplasia.

Mutational evolution in a lobular breast tumour profiled at single nucleotide resolution.

Recent advances in next generation sequencing have made it possible to precisely characterize all somatic coding mutations that occur during the development and progression of individual cancers. Here we used these approaches to sequence the genomes (>43-fold coverage) and transcriptomes of an oestrogen-receptor-alpha-positive metastatic lobular breast cancer at depth. We found 32 somatic non-synonymous coding mutations present in the metastasis, and measured the frequency of these somatic mutations in DNA from the primary tumour of the same patient, which arose 9 years earlier. Five of the 32 mutations (in ABCB11, HAUS3, SLC24A4, SNX4 and PALB2) were prevalent in the DNA of the primary tumour removed at diagnosis 9 years earlier, six (in KIF1C, USP28, MYH8, MORC1, KIAA1468 and RNASEH2A) were present at lower frequencies (1-13%), 19 were not detected in the primary tumour, and two were undetermined. The combined analysis of genome and transcriptome data revealed two new RNA-editing events that recode the amino acid sequence of SRP9 and COG3. Taken together, our data show that single nucleotide mutational heterogeneity can be a property of low or intermediate grade primary breast cancers and that significant evolution can occur with disease progression.

Tumor-associated macrophages and survival in classic Hodgkin's lymphoma.

BACKGROUND: Despite advances in treatments for Hodgkin's lymphoma, about 20% of patients still die from progressive disease. Current prognostic models predict the outcome of treatment with imperfect accuracy, and clinically relevant biomarkers have not been established to improve on the International Prognostic Score. METHODS: Using gene-expression profiling, we analyzed 130 frozen samples obtained from patients with classic Hodgkin's lymphoma during diagnostic lymph-node biopsy to determine which cellular signatures were correlated with treatment outcome. We confirmed our findings in an independent cohort of 166 patients, using immunohistochemical analysis. RESULTS: Gene-expression profiling identified a gene signature of tumor-associated macrophages that was significantly associated with primary treatment failure (P=0.02). In an independent cohort of patients, we found that an increased number of CD68+ macrophages was correlated with a shortened progression-free survival (P=0.03) and with an increased likelihood of relapse after autologous hematopoietic stem-cell transplantation (P=0.008), resulting in shortened disease-specific survival (P=0.003). In multivariate analysis, this adverse prognostic factor outperformed the International Prognostic Score for disease-specific survival (P=0.003 vs. P=0.03). The absence of an elevated number of CD68+ cells in patients with limited-stage disease defined a subgroup of patients with a long-term disease-specific survival of 100% with the use of current treatment strategies. CONCLUSIONS: An increased number of tumor-associated macrophages was strongly associated with shortened survival in patients with classic Hodgkin's lymphoma and provides a new biomarker for risk stratification.

ARID1A mutations in endometriosis-associated ovarian carcinomas.

BACKGROUND: Ovarian clear-cell and endometrioid carcinomas may arise from endometriosis, but the molecular events involved in this transformation have not been described. METHODS: We sequenced the whole transcriptomes of 18 ovarian clear-cell carcinomas and 1 ovarian clear-cell carcinoma cell line and found somatic mutations in ARID1A (the AT-rich interactive domain 1A [SWI-like] gene) in 6 of the samples. ARID1A encodes BAF250a, a key component of the SWI­SNF chromatin remodeling complex. We sequenced ARID1A in an additional 210 ovarian carcinomas and a second ovarian clear-cell carcinoma cell line and measured BAF250a expression by means of immunohistochemical analysis in an additional 455 ovarian carcinomas. RESULTS: ARID1A mutations were seen in 55 of 119 ovarian clear-cell carcinomas (46%), 10 of 33 endometrioid carcinomas (30%), and none of the 76 high-grade serous ovarian carcinomas. Seventeen carcinomas had two somatic mutations each. Loss of the BAF250a protein correlated strongly with the ovarian clear-cell carcinoma and endometrioid carcinoma subtypes and the presence of ARID1A mutations. In two patients, ARID1A mutations and loss of BAF250a expression were evident in the tumor and contiguous atypical endometriosis but not in distant endometriotic lesions. CONCLUSIONS: These data implicate ARID1A as a tumor-suppressor gene frequently disrupted in ovarian clear-cell and endometrioid carcinomas. Since ARID1A mutation and loss of BAF250a can be seen in the preneoplastic lesions, we speculate that this is an early event in the transformation of endometriosis into cancer. (Funded by the British Columbia Cancer Foundation and the Vancouver General Hospital­University of British Columbia Hospital Foundation.).

Alternative expression analysis by RNA sequencing.

In alternative expression analysis by sequencing (ALEXA-seq), we developed a method to analyze massively parallel RNA sequence data to catalog transcripts and assess differential and alternative expression of known and predicted mRNA isoforms in cells and tissues. As proof of principle, we used the approach to compare fluorouracil-resistant and -nonresistant human colorectal cancer cell lines. We assessed the sensitivity and specificity of the approach by comparison to exon tiling and splicing microarrays and validated the results with reverse transcription-PCR, quantitative PCR and Sanger sequencing. We observed global disruption of splicing in fluorouracil-resistant cells characterized by expression of new mRNA isoforms resulting from exon skipping, alternative splice site usage and intron retention. Alternative expression annotation databases, source code, a data viewer and other resources to facilitate analysis are available at http://www.alexaplatform.org/alexa_seq/.

Next-generation tag sequencing for cancer gene expression profiling.

We describe a new method, Tag-seq, which employs ultra high-throughput sequencing of 21 base pair cDNA tags for sensitive and cost-effective gene expression profiling. We compared Tag-seq data to LongSAGE data and observed improved representation of several classes of rare transcripts, including transcription factors, antisense transcripts, and intronic sequences, the latter possibly representing novel exons or genes. We observed increases in the diversity, abundance, and dynamic range of such rare transcripts and took advantage of the greater dynamic range of expression to identify, in cancers and normal libraries, altered expression ratios of alternative transcript isoforms. The strand-specific information of Tag-seq reads further allowed us to detect altered expression ratios of sense and antisense (S-AS) transcripts between cancer and normal libraries. S-AS transcripts were enriched in known cancer genes, while transcript isoforms were enriched in miRNA targeting sites. We found that transcript abundance had a stronger GC-bias in LongSAGE than Tag-seq, such that AT-rich tags were less abundant than GC-rich tags in LongSAGE. Tag-seq also performed better in gene discovery, identifying >98% of genes detected by LongSAGE and profiling a distinct subset of the transcriptome characterized by AT-rich genes, which was expressed at levels below those detectable by LongSAGE. Overall, Tag-seq is sensitive to rare transcripts, has less sequence composition bias relative to LongSAGE, and allows differential expression analysis for a greater range of transcripts, including transcripts encoding important regulatory molecules.

Prospective isolation and molecular characterization of hematopoietic stem cells with durable self-renewal potential.

Hematopoietic stem cells (HSCs) are generally defined by their dual properties of pluripotency and extensive self-renewal capacity. However, a lack of experimental clarity as to what constitutes extensive self-renewal capacity coupled with an absence of methods to prospectively isolate long-term repopulating cells with defined self-renewal activities has made it difficult to identify the essential components of the self-renewal machinery and investigate their regulation. We now show that cells capable of repopulating irradiated congenic hosts for 4 months and producing clones of cells that can be serially transplanted are selectively and highly enriched in the CD150(+) subset of the EPCR(+)CD48(-)CD45(+) fraction of mouse fetal liver and adult bone marrow cells. In contrast, cells that repopulate primary hosts for the same period but show more limited self-renewal activity are enriched in the CD150(-) subset. Comparative transcriptome analyses of these 2 subsets with each other and with HSCs whose self-renewal activity has been rapidly extinguished in vitro revealed 3 new genes (VWF, Rhob, Pld3) whose elevated expression is a consistent and selective feature of the long-term repopulating cells with durable self-renewal capacity. These findings establish the identity of a phenotypically and molecularly distinct class of pluripotent hematopoietic cells with lifelong self-renewal capacity.

High resolution analysis of follicular lymphoma genomes reveals somatic recurrent sites of copy-neutral loss of heterozygosity and copy number alterations that target single genes.

A multiplatform approach, including conventional cytogenetic techniques, BAC array comparative genomic hybridization, and Affymetrix 500K SNP arrays, was applied to the study of the tumor genomes of 25 follicular lymphoma biopsy samples with paired normal DNA samples to characterize balanced translocations, copy number imbalances, and copy-neutral loss of heterozygosity (cnLOH). In addition to the t(14;18), eight unique balanced translocations were found. Commonly reported FL-associated copy number regions were revealed including losses of 1p32-36, 6q, and 10q, and gains of 1q, 6p, 7, 12, 18, and X. The most frequent regions affected by copy-neutral loss of heterozygosity were 1p36.33 (28%), 6p21.3 (20%), 12q21.2-q24.33 (16%), and 16p13.3 (24%). We also identified by SNP analysis, 45 aberrant regions that each affected one gene, including CDKN2A, CDKN2B, FHIT, KIT, PEX14, and PTPRD, which were associated with canonical pathways involved in tumor development. This study illustrates the power of using complementary high-resolution platforms on paired tumor/normal specimens and computational analysis to provide potential insights into the significance of single-gene somatic aberrations in FL tumorigenesis.

A maternally inherited chromosome 18q22.1 deletion in a male with late-presenting diaphragmatic hernia and microphthalmia-evaluation of DSEL as a candidate gene for the diaphragmatic defect.

Using an Affymetrix GeneChip(R) Human Mapping 100K Set to study a patient with a late-presenting, right-sided diaphragmatic hernia and microphthalmia, we found a maternally inherited deletion that was 2.7 Mb in size at chromosome 18q22.1. Mapping of this deletion using fluorescence in situ hybridization revealed three deleted genes-CDH19, DSEL, and TXNDC10, and one gene that contained the deletion breakpoint, CCDC102B. We selected DSEL for further study in 125 patients with diaphragmatic hernias, as it is involved in the synthesis of decorin, a protein that is required for normal collagen formation and that is upregulated during myogenesis. We found p.Met14Ile in an unrelated patient with a late-presenting, anterior diaphragmatic hernia. In the murine diaphragm, Dsel was only weakly expressed at the time of diaphragm closure and its expression in C2C12 myoblast cells did not change significantly during myoblast differentiation, thus reducing the likelihood that the gene is involved in myogenesis of the diaphragm. Although it is possible that the 18q22.1 deletion and haploinsufficiency for DSEL contributed to the diaphragmatic defect in the patient, a definite role for DSEL and decorin in the formation of the collagen-containing, central tendon of the diaphragm has not yet been established.

Global analysis of in vivo Foxa2-binding sites in mouse adult liver using massively parallel sequencing.

Foxa2 (HNF3 beta) is a one of three, closely related transcription factors that are critical to the development and function of the mouse liver. We have used chromatin immunoprecipitation and massively parallel Illumina 1G sequencing (ChIP-Seq) to create a genome-wide profile of in vivo Foxa2-binding sites in the adult liver. More than 65% of the approximately 11.5 k genomic sites associated with Foxa2 binding, mapped to extended gene regions of annotated genes, while more than 30% of intragenic sites were located within first introns. 20.5% of all sites were further than 50 kb from any annotated gene, suggesting an association with novel gene regions. QPCR analysis demonstrated a strong positive correlation between peak height and fold enrichment for Foxa2-binding sites. We measured the relationship between Foxa2 and liver gene expression by overlapping Foxa2-binding sites with a SAGE transcriptome profile, and found that 43.5% of genes expressed in the liver were also associated with Foxa2 binding. We also identified potential Foxa2-interacting transcription factors whose motifs were enriched near Foxa2-binding sites. Our comprehensive results for in vivo Foxa2-binding sites in the mouse liver will contribute to resolving transcriptional regulatory networks that are important for adult liver function.

A comparison of MS/MS-based, stable-isotope-labeled, quantitation performance on ESI-quadrupole TOF and MALDI-TOF/TOF mass spectrometers.

The peptide-based quantitation accuracy and precision of LC-ESI (QSTAR Elite) and LC-MALDI (4800 MALDI TOF/TOF) were compared by analyzing identical Escherichia coli tryptic digests containing iTRAQ-labeled peptides of defined abundances (1:1, 2.5:1, 5:1, and 10:1). Only 51.4% of QSTAR spectra were used for quantitation by ProteinPilot Software versus 66.7% of LC-MALDI spectra. The average protein sequence coverages for LC-ESI and LC-MALDI were 24.0 and 18.2% (14.9 and 8.4 peptides per protein), respectively. The iTRAQ-based expression ratios determined by ProteinPilot from the 57 467 ESI-MS/MS and 26 085 MALDI-MS/MS spectra were analyzed for measurement accuracy and reproducibility. When the relative abundances of peptides within a sample were increased from 1:1 to 10:1, the mean ratios calculated on both instruments differed by only 0.7-6.7% between platforms. In the 10:1 experiment, up to 64.7% of iTRAQ ratios from LC-ESI MS/MS spectra failed S/N thresholds and were excluded from quantitation, while only 0.1% of the equivalent LC-MALDI iTRAQ ratios were rejected. Re-analysis of an archived LC-MALDI sample set stored for 5 months generated 3715 MS/MS spectra for quantitation, compared with 3845 acquired originally, and the average ratios differed by only 3.1%. Overall, MS/MS-based peptide quantitation performance of offline LC-MALDI was comparable with on-line LC-ESI, which required threefold less time. However, offline LC-MALDI allows the re-analysis of archived HPLC-separated samples.

Detection of pathogenic copy number variants in children with idiopathic intellectual disability using 500 K SNP array genomic hybridization.

BACKGROUND: Array genomic hybridization is being used clinically to detect pathogenic copy number variants in children with intellectual disability and other birth defects. However, there is no agreement regarding the kind of array, the distribution of probes across the genome, or the resolution that is most appropriate for clinical use. RESULTS: We performed 500 K Affymetrix GeneChip array genomic hybridization in 100 idiopathic intellectual disability trios, each comprised of a child with intellectual disability of unknown cause and both unaffected parents. We found pathogenic genomic imbalance in 16 of these 100 individuals with idiopathic intellectual disability. In comparison, we had found pathogenic genomic imbalance in 11 of 100 children with idiopathic intellectual disability in a previous cohort who had been studied by 100 K GeneChip array genomic hybridization. Among 54 intellectual disability trios selected from the previous cohort who were re-tested with 500 K GeneChip array genomic hybridization, we identified all 10 previously-detected pathogenic genomic alterations and at least one additional pathogenic copy number variant that had not been detected with 100 K GeneChip array genomic hybridization. Many benign copy number variants, including one that was de novo, were also detected with 500 K array genomic hybridization, but it was possible to distinguish the benign and pathogenic copy number variants with confidence in all but 3 (1.9%) of the 154 intellectual disability trios studied. CONCLUSION: Affymetrix GeneChip 500 K array genomic hybridization detected pathogenic genomic imbalance in 10 of 10 patients with idiopathic developmental disability in whom 100 K GeneChip array genomic hybridization had found genomic imbalance, 1 of 44 patients in whom 100 K GeneChip array genomic hybridization had found no abnormality, and 16 of 100 patients who had not previously been tested. Effective clinical interpretation of these studies requires considerable skill and experience.

Novel deletions of 14q11.2 associated with developmental delay, cognitive impairment and similar minor anomalies in three children.

METHODS AND RESULTS: We identified de novo submicroscopic chromosome 14q11.2 deletions in two children with idiopathic developmental delay and cognitive impairment. Vancouver patient 5566 has a approximately 200 kb deletion and Vancouver patient 8326 has a approximately 1.6 Mb deletion. The Database of Chromosomal Imbalance and Phenotype in Humans using Ensembl Resources (DECIPHER) revealed a third patient with idiopathic developmental delay and cognitive impairment, DECIPHER patient 126, who has a approximately 1.1 Mb deletion of 14q11.2. The deletion of patient 5566 overlaps that of patient 126 and both of these deletions lie entirely within that of patient 8326. All three children have similar dysmorphic features, including widely-spaced eyes, short nose with flat nasal bridge, long philtrum, prominent Cupid's bow of the upper lip, full lower lip and similar auricular anomalies. CONCLUSION: The minimal common deletion region on chromosome 14q11.2 is only approximately 35 kb (from 20.897 to 20.932, University of California at Santa Cruz (UCSC) Genome Browser; build hg18, March 2006) and includes only two genes, SUPT16H and CHD8, which are good candidate genes for the phenotypes. The non-recurrent breakpoints of these patients, the presence of normal copy number variants in the region and the local genomic structure support the notion that this region has reduced stability.

Sequence biases in large scale gene expression profiling data.

We present the results of a simple, statistical assay that measures the G+C content sensitivity bias of gene expression experiments without the requirement of a duplicate experiment. We analyse five gene expression profiling methods: Affymetrix GeneChip, Long Serial Analysis of Gene Expression (LongSAGE), LongSAGELite, 'Classic' Massively Parallel Signature Sequencing (MPSS) and 'Signature' MPSS. We demonstrate the methods have systematic and random errors leading to a different G+C content sensitivity. The relationship between this experimental error and the G+C content of the probe set or tag that identifies each gene influences whether the gene is detected and, if detected, the level of gene expression measured. LongSAGE has the least bias, while Signature MPSS shows a strong bias to G+C rich tags and Affymetrix data show different bias depending on the data processing method (MAS 5.0, RMA or GC-RMA). The bias in the Affymetrix data primarily impacts genes expressed at lower levels. Despite the larger sampling of the MPSS library, SAGE identifies significantly more genes (60% more RefSeq genes in a single comparison).

Assessment of algorithms for high throughput detection of genomic copy number variation in oligonucleotide microarray data.

BACKGROUND: Genomic deletions and duplications are important in the pathogenesis of diseases, such as cancer and mental retardation, and have recently been shown to occur frequently in unaffected individuals as polymorphisms. Affymetrix GeneChip whole genome sampling analysis (WGSA) combined with 100 K single nucleotide polymorphism (SNP) genotyping arrays is one of several microarray-based approaches that are now being used to detect such structural genomic changes. The popularity of this technology and its associated open source data format have resulted in the development of an increasing number of software packages for the analysis of copy number changes using these SNP arrays. RESULTS: We evaluated four publicly available software packages for high throughput copy number analysis using synthetic and empirical 100 K SNP array data sets, the latter obtained from 107 mental retardation (MR) patients and their unaffected parents and siblings. We evaluated the software with regards to overall suitability for high-throughput 100 K SNP array data analysis, as well as effectiveness of normalization, scaling with various reference sets and feature extraction, as well as true and false positive rates of genomic copy number variant (CNV) detection. CONCLUSION: We observed considerable variation among the numbers and types of candidate CNVs detected by different analysis approaches, and found that multiple programs were needed to find all real aberrations in our test set. The frequency of false positive deletions was substantial, but could be greatly reduced by using the SNP genotype information to confirm loss of heterozygosity.

Novel expressed sequences identified in a model of androgen independent prostate cancer.

BACKGROUND: Prostate cancer is the most frequently diagnosed cancer in American men, and few effective treatment options are available to patients who develop hormone-refractory prostate cancer. The molecular changes that occur to allow prostate cells to proliferate in the absence of androgens are not fully understood. RESULTS: Subtractive hybridization experiments performed with samples from an in vivo model of hormonal progression identified 25 expressed sequences representing novel human transcripts. Intriguingly, these 25 sequences have small open-reading frames and are not highly conserved through evolution, suggesting many of these novel expressed sequences may be derived from untranslated regions of novel transcripts or from non-coding transcripts. Examination of a large metalibrary of human Serial Analysis of Gene Expression (SAGE) tags demonstrated that only three of these novel sequences had been previously detected. RT-PCR experiments confirmed that the 6 sequences tested were expressed in specific human tissues, as well as in clinical samples of prostate cancer. Further RT-PCR experiments for five of these fragments indicated they originated from large untranslated regions of unannotated transcripts. CONCLUSION: This study underlines the value of using complementary techniques in the annotation of the human genome. The tissue-specific expression of 4 of the 6 clones tested indicates the expression of these novel transcripts is tightly regulated, and future work will determine the possible role(s) these novel transcripts may play in the progression of prostate cancer.

Analysis of the prostate cancer cell line LNCaP transcriptome using a sequencing-by-synthesis approach.

BACKGROUND: High throughput sequencing-by-synthesis is an emerging technology that allows the rapid production of millions of bases of data. Although the sequence reads are short, they can readily be used for re-sequencing. By re-sequencing the mRNA products of a cell, one may rapidly discover polymorphisms and splice variants particular to that cell. RESULTS: We present the utility of massively parallel sequencing by synthesis for profiling the transcriptome of a human prostate cancer cell-line, LNCaP, that has been treated with the synthetic androgen, R1881. Through the generation of approximately 20 megabases (MB) of EST data, we detect transcription from over 10,000 gene loci, 25 previously undescribed alternative splicing events involving known exons, and over 1,500 high quality single nucleotide discrepancies with the reference human sequence. Further, we map nearly 10,000 ESTs to positions on the genome where no transcription is currently predicted to occur. We also characterize various obstacles with using sequencing by synthesis for transcriptome analysis and propose solutions to these problems. CONCLUSION: The use of high-throughput sequencing-by-synthesis methods for transcript profiling allows the specific and sensitive detection of many of a cell's transcripts, and also allows the discovery of high quality base discrepancies, and alternative splice variants. Thus, this technology may provide an effective means of understanding various disease states, discovering novel targets for disease treatment, and discovery of novel transcripts.

Epigenetic and transcriptional determinants of the human breast.

While significant effort has been dedicated to the characterization of epigenetic changes associated with prenatal differentiation, relatively little is known about the epigenetic changes that accompany post-natal differentiation where fully functional differentiated cell types with limited lifespans arise. Here we sought to address this gap by generating epigenomic and transcriptional profiles from primary human breast cell types isolated from disease-free human subjects. From these data we define a comprehensive human breast transcriptional network, including a set of myoepithelial- and luminal epithelial-specific intronic retention events. Intersection of epigenetic states with RNA expression from distinct breast epithelium lineages demonstrates that mCpG provides a stable record of exonic and intronic usage, whereas H3K36me3 is dynamic. We find a striking asymmetry in epigenomic reprogramming between luminal and myoepithelial cell types, with the genomes of luminal cells harbouring more than twice the number of hypomethylated enhancer elements compared with myoepithelial cells.